- Research article
- Open Access
- Open Peer Review
The CTLA4 variants may interact with the IL23R- and NOD2-conferred risk in development of Crohn's disease
BMC Medical Geneticsvolume 11, Article number: 91 (2010)
The CTLA4 (cytotoxic T-lymphocyte antigen 4) gene is associated with several immunopathologic diseases and because of its important immuno-regulatory role it may be considered also a plausible candidate for a genetic association with inflammatory bowel diseases. Previously published studies found no association of CTLA4 with Crohn's disease itself, but some indicated an association with its subphenotypes. The aim of this study was to assess the association in the Czech population, using a set of markers shown to associate with other diseases.
Six polymorphisms within the CTLA4 region were investigated in 333 patients with Crohn's disease and 482 unrelated healthy controls, all Caucasians of Czech origin. The genotypes of the SNPs were determined using the TaqMan SNP genotyping assays. Haplotypes were reconstructed using an expectation-maximization algorithm, and their association with the condition was assessed using log-linear modeling. Then, potential interactions were tested between the CTLA4 variants and other genetic factors known to confer the disease susceptibility.
No crude associations with Crohn's disease were found for the tested CTLA4 variants under the log-additive or dominant models. However, when stratified for the genetic risk conferred by the variants in the NOD2 (the p.Leu1007fsX1008, rs5743293) or the IL23R (p.R381Q, rs11209026), a significant negative association emerged for the minor alleles of CTLA4 CT60 (rs3087243), JO31 (rs11571302), JO27-1 (rs11571297) polymorphisms. This negative association with CTLA4 was apparent only in the strata defined by presence minor alleles at the NOD2 rs5743293 (here the CTLA4 CT60 A coffered an OR = 0.43, 95%CI 0.19 - 0.95 for the presence of CT60 A), or IL23R rs11209026 (here the OR for presence of CT60 A was 0.23, 95%CI 0.07 - 0.71). We observed this effect also for the haplotype consisting of minor alleles of the three tightly linked CTLA4 markers. Furthermore, this haplotype was associated with the younger age at diagnosis (OR 1.52, 95%CI 1.09 - 2.11, p = 0.014).
A protective effect of a CTLA4 haplotype was unmasked after stratification for the risk variants in the NOD2 and IL23R genes, and may point towards the biological relevance of the molecule in the pathogenesis of the disease.
Crohn's disease (CD) belongs to inflammatory bowel diseases (IBD) that are characterized by chronic, relapsing and recurrent inflammation of intestinal mucosa. The disease is thought to result from the action of environmental factors in genetically susceptible individuals. Three variants in the NOD2 [1, 2], IBD5 locus  and one variant in the IL23R  and in the ATG16L1  have been independently confirmed to be associated with CD, including associations found previously in the Czech population [6, 7]. Recent studies, however, show that this list is far from being complete [5, 8–12].
The CTLA4 gene may also be considered as a plausible candidate for a genetic association with IBD. Its product, the cytotoxic T-lymphocyte-associated protein 4 (CTLA4) is a T-cell suppressor which plays an essential role in the function of the CD25(+)CD4(+) regulatory cells that control the process of intestinal inflammation [13, 14]. The CTLA4 gene maps within the 2q33 region that has been found to carry suggestive linkage significance for IBD . The CTLA4 gene is associated with other immunopathologic diseases (type 1 diabetes, Graves' disease, Addison's disease, celiac disease, systemic lupus erythematosus, rheumatoid arthritis, vitiligo) . Among the studied single nucleotide polymorphisms (SNPs), the CT60 (rs3087243) shows the most prominent associations, being followed by other three SNPs: JO31 (rs11571302), JO30 (rs7565213) and JO27-1 (rs11571297) . A recent publication has shown evidence for association of another SNP within the CTLA4 with type 1 diabetes, the rs1427676 . Previously published papers about genetic association with CD tested three variants in the CTLA4 gene: g.49A > G (rs231775), g.-318C > T (rs5742909) and the previously mentioned CT60, having found no association [18–21]. However, several works suggested that CTLA4 variants may influence the phenotype of CD [18, 19].
The aim of this study was to assess the association in the Czech population, using a set of markers previously shown to associate with other diseases.
In a case-control design, 333 Czech patients were compared to 482 unrelated healthy Czech controls representing a general population sample from the same geographical region. We tested 137 pediatric-onset patients (71 boys, 66 girls) who developed CD under or at the age of 18 years and were diagnosed according to the Porto criteria , and 196 adult onset patients (77 males, 119 females) diagnosed according to endoscopic, radiological, histological and clinical criteria. Phenotypic classification was done according to the Montreal Classification . The demographic and clinical characteristics of the patients are listed in Table 1 and Table 2. The control group included 482 individuals: 295 children, 187 adult; 311 males, 171 females; median age 12 years, interquartile range 7-34 years. The study was approved by the Ethics Committees of the authors' institutions, and a written informed consent was obtained from all participants or their guardians.
Genomic DNA was extracted from peripheral blood leukocytes with a routine salting out method, or from salivary samples using Oragene DNA Self-Collection Kit according the manufacturer's protocol (DNA Genotek Inc., Ottawa, Ontario, Canada). One SNP proximal to CTLA4 (rs736611), one from within the gene (g.49G > A, rs231775), and four SNPs located distally from the coding part of CTLA4 gene (rs3087243 also called CT60; rs11571302 called JO31; rs11571297 called JO27-1; and rs1427676) were selected based on available literature and genotyped using the TaqMan SNP genotyping assays (TaqMan SNP Genotyping Assay by Applied Biosystems, Foster City, CA, USA). The assays were run on an ABI 7300 machine (Applied Biosystems, Foster City, CA, USA) and evaluated according to manufacturer's instructions. To ensure consistency between runs, samples of known genotypes were repeated in every analysis. For testing interactions with other associated genes, we used genotypes generated in previously published studies on this sample set [6, 7].
The Hardy-Weinberg equilibrium was checked by comparing observed to expected genotype frequencies in the control subjects, and tested using exact tests. Associations of particular SNPs with CD were evaluated using odds ratios (OR) with 95% confidence intervals (CI). Haplotype analysis was performed by estimating the haplotype frequencies by the expectation-maximization algorithm implemented in the R-project package 'haplo.stats' version 1.3.1. Association of haplotypes with the conditions was tested using log-linear modeling. Then, a potential interaction between the CTLA4 variants and other genetic factors associated with the autoimmune conditions were tested. The statistical analysis was performed using the R-project package 'SNPassoc' version 1.5-2 .
The frequencies of the variants and respective OR are listed in Table 3. No crude associations with CD were found for the tested SNPs under the log-additive or dominant models. The genotype distributions in control subjects conformed to Hardy-Weinberg equilibrium in all SNPs (p > 0.20) except the rs1427676 (p = 0.014 in exact tests) which was therefore excluded from all further analyses.
As the part of chromosome under the CTLA4 gene is divided into the several blocks  we performed a haplotype analysis using the five SNPs; no crude association with CD was observed (data not shown).
Interaction of the CTLA4 SNPs with variants in IL23R andNOD2
We then tested possible interactions between variants in the CTLA4 and polymorphisms in other genes previously associated with CD: NOD2 gene p.Leu1007fsX1008 (c.3020insC), IL23R gene rs11209026 (c.1142G > A) [6, 7], see Figure 1. This was done in dominant models using an R-project package 'SNPassoc' version 1.5-2 . Significant interactions were observed between the three CTLA4 variants (CT60, JO31, JO27-1) and NOD2 p.Leu1007fsX1008, and the same variants in the CTLA4 and IL23R rs11209026.
For a quantification of the CTLA4 association stratified by the above NOD2 and IL23R polymorphisms see Table 4 a Table 5: the minor alleles of the CT60, JO31 and JO27-1 within the CTLA4 modified the risk of Crohn's disease in the stratum of subjects carrying the frameshift insertion p.Leu1007fsX1008 in NOD2, while no perceivable effect of CTLA4 was found in the stratum of p.Leu1007fsX1008 wild-type homozygotes. Similarly, the three CTLA4 variants clearly, albeit moderately, decreased the risk of CD in the stratum of subjects carrying minor alleles of rs11209026 within the IL23R (G/A and A/A), while no effect was observed in the IL23R wild-type homozygotes. The effect was observable also for the haplotype consisting of the three minor alleles of the variants in tight linkage disequilibrium, the CT60 "A", JO31 "T", JO27-1 "G" haplotype.
Genotype - phenotype analysis
Using a case-only design, we tested whether the phenotypic characteristics of the patients are dependent on carriage status of the minor alleles at CTLA4 variants.
Table 6 shows logistic regression with the outcomes of clinical characteristics and the three CTLA4 variants and the haplotype as the predictors. The pediatric-onset patients differed to the adult-onset patients in their frequencies of the minor allele at the CT60 (74% versus 63%, p = 0.03). The CT60, JO31 and JO27-1 SNPs, as well as their "ATG" haplotype retained their associations with the age at diagnosis after adjustment to the effect of the NOD2 variant (p.Leu1007fsX1008). No difference between pediatric and adult-onset group was found for the other two variants in the CTLA4 gene, data not shown.
The CTLA4 variants were weakly associated with the ileal-only (L1) and ileocolonic involvement (L3) in a dominant manner, while no association was observed with any of the remaining clinical characteristics: localization in the upper gastrointestinal tract, the stricturing or penetrating behavior of the disease, perianal disease, extraintestinal manifestation, or the need for abdominal surgery (data not shown).
The immunologic importance of the CTLA4 gene is in striking contrast to the lack of knowledge on the functional relevance of its numerous polymorphisms. Consequently, many groups have investigated various polymorphisms located within various regions of the gene. The first published study on CTLA4 variants in CD investigated the g.49A > G (rs231775) and g.-318C > T (rs5742909) in the Dutch and the Chinese populations, finding only an association with the age of onset . Similarly, in a Hungarian work, no association of g.49A > G with CD was detected . Since the work by Ueda et al  had been published on dissecting the association of CTLA4 with immunopathological diseases, further investigators focused on the CT60 polymorphism. This variant was studied in the Japanese  and the Spanish  populations, however no crude association with CD was detected. The G/G genotype of g.49A > G was associated with penetrating form of CD in the Japanese dataset . No association within 2q33 chromosomal region has been found by genome-wide studies [4, 5, 8–12].
Thus, compelling evidence has been gathered against simple association of the disease itself with the polymorphisms of CTLA4. In line with these findings, we observed no crude association unless further genetic factors were taken into account. However, when CTLA4 was considered as a modifier of the effects conferred by the NOD2 and IL23R genes, possible interactions substantiated. Interactions in multifactorial immunopathological diseases are not infrequent: in CD, the interactions with the NOD2 gene were detected in the IBD5 locus , IBD6 locus , TNFA , DLG5 , ATG16L1, IL23R , TLR4  and in CD14 . The interaction was also found between IBD5 locus and IL23R  and between Toll-like receptor-9 polymorphisms and variants in NOD2 and IL23R .
The interactions we found for the CT60, JO31, JO27-1 variants of CTLA4 (or their haplotype) with the p.Leu1007fsX1008 variant of NOD2 may imply that the effect of the strongest risk variant within the NOD2 (p.Leu1007fsX1008) can be expressed better on the background of the common CTLA4 haplotype. This is suggestive of a complex pattern of gene-gene interaction that may merit pursuing further functional studies. Similarly, the risk haplotype of CTLA4 also interacts with the IL23R protective variants. This rather weak interaction can be also due to the relatively limited size of the dataset.
In addition to the modifying effect on the risk of the disease itself, we observed an association with the age at onset and the disease subphenotypes. Indeed, the impact of genetic factors in early-onset patients with CD seems to be stronger than in adult-onset patients (reviewed by de Ridder L et al. ). In our dataset, the age at diagnosis was associated with CT60, JO31 and JO27-1. Influence of CTLA4 variants on the age at diagnosis has been previously described by Xia et al, although with a different SNP (g.-318C > T) . Moreover, their patients were divided into groups where 40 years of age was the cut off, not 18 years as in our study.
The CTLA4 was associated with the ileal and ileocolonic involvement in our case set: up to our knowledge, this is the first time when localization of CD is influenced by any variant within CTLA4. It should be however noted that these associations are weak, merit further investigation in other populations, and their clinical relevance can be only hardly envisaged. The ileal form of disease has been shown more common in adult-onset patients and more common in patients carrying minor variants of the NOD2 gene. A possible explanation of the association of CTLA4 with localization of disease might be found in the interaction between CTLA4 and NOD2 gene.
Similarly to Machida et al  we also tested whether the g.49A > G variant influences the occurrence of penetrating disease, but we were not able to confirm this association. However, the genetic background between Japanese and Czech populations differs markedly.
We present a study of genetic association of polymorphisms within the CTLA4 gene with CD and its subphenotypes, using a representative set of markers previously reported from other studies. We observed interactions of the CTLA4 haplotype with variants in NOD2 and IL23R genes, and detected an effect of three variants of the CTLA4 on the age at diagnosis and localization of the disease.
Hugot JP, Chamaillard M, Zouali H, Lesage S, Cézard JP, Belaiche J, Almer S, Tysk C, O'Morain CA, Gassull M, et al: Association of NOD2 leucine-rich repeat variants with susceptibility to Crohn's disease. Nature. 2001, 411 (6837): 599-603. 10.1038/35079107.
Ogura Y, Bonen DK, Inohara N, Nicolae DL, Chen FF, Ramos R, Britton H, Moran T, Karaliuskas R, Duerr RH, et al: A frameshift mutation in NOD2 associated with susceptibility to Crohn's disease. Nature. 2001, 411 (6837): 603-606. 10.1038/35079114.
Rioux JD, Daly MJ, Silverberg MS, Lindblad K, Steinhart H, Cohen Z, Delmonte T, Kocher K, Miller K, Guschwan S, et al: Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease. Nat Genet. 2001, 29 (2): 223-228. 10.1038/ng1001-223.
Duerr RH, Taylor KD, Brant SR, Rioux JD, Silverberg MS, Daly MJ, Steinhart AH, Abraham C, Regueiro M, Griffiths A, et al: A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science. 2006, 314 (5804): 1461-1463. 10.1126/science.1135245.
Rioux JD, Xavier RJ, Taylor KD, Silverberg MS, Goyette P, Huett A, Green T, Kuballa P, Barmada MM, Datta LW, et al: Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis. Nat Genet. 2007, 39 (5): 596-604. 10.1038/ng2032.
Hradsky O, Lenicek M, Dusatkova P, Bronsky J, Nevoral J, Valtrova V, Kotalova R, Szitanyi P, Petro R, Starzykova V, et al: Variants of CARD15, TNFA and PTPN22 and susceptibility to Crohn's disease in the Czech population: high frequency of the CARD15 1007fs. Tissue Antigens. 2008, 71 (6): 538-547. 10.1111/j.1399-0039.2008.01047.x.
Dusatkova P, Hradsky O, Lenicek M, Bronsky J, Nevoral J, Kotalova R, Bajerova K, Vitek L, Lukas M, Cinek O: Association of IL23R p.381Gln and ATG16L1 p.197Ala With Crohn Disease in the Czech Population. J Pediatr Gastroenterol Nutr. 2009, 49 (4): 405-10. 10.1097/MPG.0b013e31819344ee.
The Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678. 10.1038/nature05911.
Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, Rioux JD, Brant SR, Silverberg MS, Taylor KD, Barmada MM, et al: Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease. Nat Genet. 2008, 40 (8): 955-962. 10.1038/ng.175.
Tremelling M, Parkes M: Genome-wide association scans identify multiple confirmed susceptibility loci for Crohn's disease: lessons for study design. Inflamm Bowel Dis. 2007, 13 (12): 1554-1560. 10.1002/ibd.20239.
Parkes M, Barrett JC, Prescott NJ, Tremelling M, Anderson CA, Fisher SA, Roberts RG, Nimmo ER, Cummings FR, Soars D, et al: Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility. Nat Genet. 2007, 39 (7): 830-832. 10.1038/ng2061.
Libioulle C, Louis E, Hansoul S, Sandor C, Farnir F, Franchimont D, Vermeire S, Dewit O, de Vos M, Dixon A, et al: Novel Crohn disease locus identified by genome-wide association maps to a gene desert on 5p13.1 and modulates expression of PTGER4. PLoS Genet. 2007, 3 (4): e58-10.1371/journal.pgen.0030058.
Read S, Malmström V, Powrie F: Cytotoxic T lymphocyte-associated antigen 4 plays an essential role in the function of CD25(+)CD4(+) regulatory cells that control intestinal inflammation. J Exp Med. 2000, 192 (2): 295-302. 10.1084/jem.192.2.295.
Liu Z, Geboes K, Hellings P, Maerten P, Heremans H, Vandenberghe P, Boon L, van Kooten P, Rutgeerts P, Ceuppens JL: B7 interactions with CD28 and CTLA-4 control tolerance or induction of mucosal inflammation in chronic experimental colitis. J Immunol. 2001, 167 (3): 1830-1838.
van Heel DA, Fisher SA, Kirby A, Daly MJ, Rioux JD, Lewis CM: Inflammatory bowel disease susceptibility loci defined by genome scan meta-analysis of 1952 affected relative pairs. Hum Mol Genet. 2004, 13 (7): 763-770. 10.1093/hmg/ddh090.
Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, et al: Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature. 2003, 423 (6939): 506-511. 10.1038/nature01621.
Howson JM, Walker NM, Smyth DJ, Todd JA: Analysis of 19 genes for association with type I diabetes in the Type I Diabetes Genetics Consortium families. Genes Immun. 2009, 10 (Suppl 1): S74-84. 10.1038/gene.2009.96.
Xia B, Crusius JB, Wu J, Zwiers A, van Bodegraven AA, Pena AS: CTLA4 gene polymorphisms in Dutch and Chinese patients with inflammatory bowel disease. Scand J Gastroenterol. 2002, 37 (11): 1296-1300. 10.1080/003655202761020579.
Machida H, Tsukamoto K, Wen CY, Narumi Y, Shikuwa S, Isomoto H, Takeshima F, Mizuta Y, Niikawa N, Murata I, et al: Association of polymorphic alleles of CTLA4 with inflammatory bowel disease in the Japanese. World J Gastroenterol. 2005, 11 (27): 4188-4193.
Rueda B, Zhernakova A, Lopez-Nevot MA, Gomez-Garcia M, Ortega E, Pinero A, Correro F, Brieva JA, Nieto A, Koeleman BP, et al: CTLA4/CT60 polymorphism is not relevant in susceptibility to autoimmune inflammatory intestinal disorders. Hum Immunol. 2005, 66 (3): 321-325. 10.1016/j.humimm.2004.11.005.
Magyari L, Faragó B, Bene J, Horvatovich K, Lakner L, Varga M, Figler M, Gasztonyi B, Mózsik G, Melegh B: No association of the cytotoxic T-lymphocyte associated gene CTLA4 +49A/G polymorphisms with Crohn's disease and ulcerative colitis in Hungarian population samples. World J Gastroenterol. 2007, 13 (15): 2205-2208.
IBD Working Group of the European Society for Paediatric Gastroenterology Hepatology and Nutrition.: Inflammatory bowel disease in children and adolescents: recommendations for diagnosis--the Porto criteria. J Pediatr Gastroenterol Nutr. 2005, 41 (1): 1-7. 10.1097/01.MPG.0000163736.30261.82.
Silverberg MS, Satsangi J, Ahmad T, Arnott ID, Bernstein CN, Brant SR, Caprilli R, Colombel JF, Gasche C, Geboes K, et al: Toward an integrated clinical, molecular and serological classification of inflammatory bowel disease: Report of a Working Party of the 2005 Montreal World Congress of Gastroenterology. Can J Gastroenterol. 2005, 19 (Suppl A): 5-36.
González JR, Armengol L, Solé X, Guinó E, Mercader JM, Estivill X, Moreno V: SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007, 23 (5): 644-645. 10.1093/bioinformatics/btm025.
Armuzzi A, Ahmad T, Ling KL, de Silva A, Cullen S, van Heel D, Orchard TR, Welsh KI, Marshall SE, Jewell DP: Genotype-phenotype analysis of the Crohn's disease susceptibility haplotype on chromosome 5q31. Gut. 2003, 52 (8): 1133-1139. 10.1136/gut.52.8.1133.
Linderson Y, Bresso F, Buentke E, Pettersson S, D'Amato M: Functional interaction of CARD15/NOD2 and Crohn's disease-associated TNFalpha polymorphisms. Int J Colorectal Dis. 2005, 20 (4): 305-311. 10.1007/s00384-004-0732-z.
Stoll M, Corneliussen B, Costello CM, Waetzig GH, Mellgard B, Koch WA, Rosenstiel P, Albrecht M, Croucher PJ, Seegert D, et al: Genetic variation in DLG5 is associated with inflammatory bowel disease. Nat Genet. 2004, 36 (5): 476-480. 10.1038/ng1345.
Hampe J, Franke A, Rosenstiel P, Till A, Teuber M, Huse K, Albrecht M, Mayr G, De La Vega FM, Briggs J, et al: A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1. Nat Genet. 2007, 39 (2): 207-211. 10.1038/ng1954.
Roberts RL, Gearry RB, Hollis-Moffatt JE, Miller AL, Reid J, Abkevich V, Timms KM, Gutin A, Lanchbury JS, Merriman TR, et al: IL23R R381Q and ATG16L1 T300A are strongly associated with Crohn's disease in a study of New Zealand Caucasians with inflammatory bowel disease. Am J Gastroenterol. 2007, 102 (12): 2754-2761. 10.1111/j.1572-0241.2007.01525.x.
Browning BL, Annese V, Barclay ML, Bingham SA, Brand S, Büning C, Castro M, Cucchiara S, Dallapiccola B, Drummond H, et al: Gender-stratified analysis of DLG5 R30Q in 4707 patients with Crohn disease and 4973 controls from 12 Caucasian cohorts. J Med Genet. 2008, 45 (1): 36-42. 10.1136/jmg.2007.050773.
Cummings JR, Ahmad T, Geremia A, Beckly J, Cooney R, Hancock L, Pathan S, Guo C, Cardon LR, Jewell DP: Contribution of the novel inflammatory bowel disease gene IL23R to disease susceptibility and phenotype. Inflamm Bowel Dis. 2007, 13 (9): 1063-1068. 10.1002/ibd.20180.
Török HP, Glas J, Endres I, Tonenchi L, Teshome MY, Wetzke M, Klein W, Lohse P, Ochsenkühn T, Folwaczny M, et al: Epistasis between Toll-like receptor-9 polymorphisms and variants in NOD2 and IL23R modulates susceptibility to Crohn's disease. Am J Gastroenterol. 2009, 104 (7): 1723-1733. 10.1038/ajg.2009.184.
de Ridder L, Weersma RK, Dijkstra G, van der Steege G, Benninga MA, Nolte IM, Taminiau JA, Hommes DW, Stokkers PC: Genetic susceptibility has a more important role in pediatric-onset Crohn's disease than in adult-onset Crohn's disease. Inflamm Bowel Dis. 2007, 13 (9): 1083-1092. 10.1002/ibd.20171.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/91/prepub
This study was supported by grants from the Czech Ministry of Health (MZOFNM2005), the University Hospital Motol and IGA - NR/9219-3/2007. Drs V. Valtrova, P.Szitanyi, R. Petro, M. Bortlik and V. Starzykova are gratefully acknowledged for providing us with samples and clinical data from their patients.
The authors declare that they have no competing interests.
OH, PD and MLe performed the experiments; JB, JN, LV, MLu coordinated and performed the collection of the samples and were also involved in editing the manuscript; OC and OH designed the study and wrote the manuscript. All authors read and approved the final manuscript.