On the genetic involvement of apoptosis-related genes in Crohn's disease as revealed by an extended association screen using 245 markers: no evidence for new predisposing factors
© Wagenleiter et al; licensee BioMed Central Ltd. 2005
Received: 07 July 2005
Accepted: 30 November 2005
Published: 30 November 2005
Crohn's disease (CD) presents as an inflammatory barrier disease with characteristic destructive processes in the intestinal wall. Although the pathomechanisms of CD are still not exactly understood, there is evidence that, in addition to e.g. bacterial colonisation, genetic predisposition contributes to the development of CD. In order to search for predisposing genetic factors we scrutinised 245 microsatellite markers in a population-based linkage mapping study. These microsatellites cover gene loci the encoded protein of which take part in the regulation of apoptosis and (innate) immune processes. Respective loci contribute to the activation/suppression of apoptosis, are involved in signal transduction and cell cycle regulators or they belong to the tumor necrosis factor superfamily, caspase related genes or the BCL2 family. Furthermore, several cytokines as well as chemokines were included. The approach is based on three steps: analyzing pooled DNAs of patients and controls, verification of significantly differing microsatellite markers by genotyping individual DNA samples and, finally, additional reinvestigation of the respective gene in the region covered by the associated microsatellite by analysing single-nucleotide polymorphisms (SNPs). Using this step-wise process we were unable to demonstrate evidence for genetic predisposition of the chosen apoptosis- and immunity-related genes with respect to susceptibility for CD.
Crohn's disease (CD) is a chronic inflammatory disorder characterized by destructive processes in the intestinal wall. Interactions between genetic and environmental factors potentially lead to an imbalance between the luminal bacterial flora, and the innate as well as the adaptive immune systems [1, 2]. Epidemiological and genome wide studies have lead to the identification of factors establishing genetic involvement in CD [1, 3, 4]. Despite of fundamental findings, namely the variation in the CARD15 receptor and their association with CD, the causative instances regulating the exaggerated mucosal response remained elusive. The proposed pathomechanisms of CD are manifold. The dysregulated response of the innate immune system is supposed to present a crucial step in the pathogenesis of CD . This fact has been confirmed genetically by several CD associations of genes such as CD14, TLR4 and in some instances the interaction of their variations with CARD15 [6, 7]. In regard to the polarized T helper (Th) response, the adaptive immune system appears affected in CD as well [8–10]. Moreover, several studies implicated a role of programmed cell death in CD [11–15]. Apoptosis mediates 'self-tolerance', the elimination of autoreactive immune compartments. In addition, the thoroughly controlled termination of a physiological immune response is due to the process of programmed cell death. In CD mucosal T cells show less susceptibility to apoptosis . In this context TNFα protein exerts multiple physiological effects, and anti-TNFα therapeutic strategies (e.g. infliximab) are effective in (maintaining) remission of CD . In several studies it has been revealed that treatment of CD patients with infliximab leads to an activation of T cells rendering them susceptible for apoptosis [18, 19]. Interestingly, the effect of this treatment may not be due to neutralisation of soluble TNFα (and its binding to the TNFRs), but rather it may be caused by its affinity to membrane-bound TNFα putatively changing the ratio of anti- and pro-apoptotic mediators towards induction of apoptosis [18, 20]. Although the mechanisms of the causal role of T cells responses in CD remain to be determined in detail, there is substantial clinical evidence that suggests a role for uncontrolled activated T lymphocytes in the pathogenic process of CD [21–24]. Nevertheless, it is uncertain, whether a genetic basis for a decreased activation/apoptosis of T lymphocytes in CD patients exists, and whether increased anti-apoptotic markers, found in T cells of these patients are due to the mucosal inflammation, secondarily .
In such a complex situation we used extended association screening (EAS) with markers representing 245 apoptosis- and (innate) immunity-related genes. The majority of the investigated markers have been successfully utilized in respective studies before [25, 26]. Our population based linkage mapping comprises a 3-stage analysis with pooled DNA in the initial phase and subsequently individual genotyping. In order to confirm such results, several tagging SNPs of the adjacent gene represented by the marker were analysed. Here, we investigated the role of distinct biological pathways for the susceptibility of CD.
Materials and methods
One hundred and fifty eight well-characterized patients with a clinical, endoscopical and histological diagnosis of CD were included. This patient cohort has been reported before [27, 28]. All patients were of German origin and the diagnosis of CD was adjusted according to the diagnostic criteria of the European Community Workshop on Inflammatory Bowel Diseases (IBD). As controls a group of healthy northern German (NoG) and western German (WeG) origin were analysed. In the initial step a group of ~100 NoG individuals were used. In order to exclude population stratification, genotyping of chosen SNPs was performed in 180–460 NoG and WeG individuals.
Pooling of DNA
The DNA concentration from each individual of the patient and control cohorts was quantified by spectrophotometry, carried out four times, and then diluted accordingly to 100 ng/μl. In a second step the DNA was diluted to a concentration of 65 ng/μl and once more measured by spectrophotometry. Finally, DNA diluted to 50 ng/μl was adjusted to a final amount of 1000 ng for each individual in a pool of 50 persons. In the initial stage, marker analyses were performed with two patient and two control subpools, respectively.
Tailed primer PCR
Tailed primer PCR was performed as described before : An 18 bp-tail was added to each sense oligonucleotide. PCR reaction included three oligonulceotides, two of which were target specific. The third one consists of the same sequence as the abovementioned tail that was additionally fluorescence-labelled.
Genes investigated for CD association as represented by an intra- or juxtagenic microsatellite marker (for additional information see URL: http://www.ruhr-uni-bochum.de/mhg/marker_information_SENW.pdf)
cell cycle regulators
Statistics for initial comparisons of allele frequencies
Raw data from ABI377 profiles were analysed by the Genotyper software (ABI) producing a marker specific allele image profile (AIP) which includes different heights of peaks reflecting the allele frequencies. In order to test differences of the AIPs between CD patients and the controls, all peak heights were summarized for each pool and set to 100 %. The total allele count for each distinct allele was then estimated. Thereupon, the AIPs of the case and control pools were compared statistically by means of contingency tables. Hence, P values are nominal and approximate, because estimated rather than observed counts were used for allele frequencies. The significance level was set at p = 0.05. In order to focus the statistics on major alleles, all minor alleles with a frequency of less than 0.05 were summarized to a virtual allele. Subsequently, a second statistical analysis by means of contingency tables was undertaken. A third step for statistical testing each allele individually was accomplished (and the summation of all other marker alleles), whereby the respective value of the patient group was compared with those of the controls and subsequent χ2 analyses. Despite of evidence that correction for multiple comparisons might eliminate 'real positive' results , Q value correction was performed with a cut off of 5% for the initial screening procedure .
Nevertheless, for selecting markers for further investigations, non-corrected P values were simply ranked according to their evidence for association including all performed statistical procedures.
P values for microsatellite markers located intragenically or in the immediate vicinity of represented genes after the initial step and individual genotyping.
gene (as represented by the respective marker)
after analysis with pooled DNA
after summation of alleles beneath 5%
after analyses of each single allele (most significant allele)
after individual genotyping 1 (p c value)
after correction by Q-value of pooled data
0.0044 (pc > 0.05; c = 9)
0.0012 (pc < 0.01; c = 6)
0.0050 (pc < 0.04; c = 7)
Investigated SNPs in genes as represented by significantly differing microsatellites of the individual genotyping step.
Allele: fragment length (bp)
* 1 GTGAAATTAATTTCAGAGC TGT/CTCAGCTTCCCCAAGGAT
Microsatellites representing 245 genes involved in apoptosis regulation (see Tab. 1) were investigated by using EAS. None of the markers presented with significant intra-subgroup differences confirming the homogeneity of the pools. The statistical evaluation of the microsatellite frequencies in the CD patient and the control cohorts revealed 9 significantly different allele distributions of intra- or juxtagenic markers for FLIP, BCL2A1, BAG1, BPI, erbB3, TP73, TLR9, TNFRSF17 and CARD15 (summarized in Tab. 2).
Individual genotyping confirmed significant P values only for the 3 markers FLIP (p = 0.0044, pc > 0.05, in HWE), TNFRSF17 (p = 0.0012, pc < 0.01, in HWE) and the positive control CARD15 (p = 0.0050, pc < 0.04, in HWE). The additional associations for the other markers were rejected (see Tab. 2 and Additional file 1). There were no differences analysing CARD15+ and CARD15- patients.
SNP markers (Tab. 3) were genotyped located in the respective genes in the vicinity of the microsatellites representing TNFRSF17 and FLIP. Thus, SNPs were analyzed spread across the genes representing haplotypes as predisposed by the 'LD Select' method reported before . RFLP analyses did not reveal any association of the selected SNPs, neither by comparing the CARD15+ nor the CARD15- patients with the control group.
Comparison of TNFRSF17 microsatellite alleles
The genotypes of the TNFRSF17 microsatellite alleles were compared between the patient and control cohorts. Analyses revealed evidence either for a predisposing (allele 3) and a protective allele (2) or linkage between these alleles and the marker alleles, respectively.
Allele and genotype frequencies of the functional MHC2TA polymorphism (rs3087456).
CD (n = 147)
controls (n = 463)
The pathomechanisms of CD are still not exactly understood, albeit certain CARD15 variations appear especially frequent in CD patients; thus genetic involvement is proven. These genetic predisposition factors, however, are neither sufficient nor explain they the pathogenesis in all CD patients. In this study we present an association screen mainly for apoptosis and immunity related genes by microsatellite markers as investigated in a 3-step approach.
Our initial analyses revealed 9 significantly different allele distributions of intra- or juxtagenic markers for FLIP, BCL2A1, BAG1, BPI, erbB3, TP73, TLR9, TNFRSF17 and CARD15 (see Tab. 2). Yet, after correction by Q-value, none of those markers remained significant. On the other hand, a recent study raised the question, whether the correction for multiple comparisons should be applied at all in EAS . For example, in these analyses a previously significantly associated microsatellite (representing the TNFα gene), which has been used as a positive control such as CARD15, would have been rejected by the correction procedure. Therefore, it remains conceivable that the abovementioned markers represent rather hints for additional predisposing factors/loci with low effect size.
The most promising markers (reflected by a significant p-value) were included in further analyses regardless of the correction procedure. Individual genotyping rejected most markers found to be significantly different in the initial step of our approach and only three markers remained significant representing the TNFRSF17, FLIP, CARD15 genes (Tab. 2). Obviously, pooled and individual genotyping yield somewhat contradictory results. Eight microsatellites revealed significantly differences between the patient and control cohorts after the pooling procedure, whereas individual genotyping results in the confirmation of 'only' 2 markers. These conspicuous differences might be due to several artefacts caused by analyses with pooled DNA. For example, a typical artefact is the length-dependent amplification of short alleles or the presence of null-alleles. Additionally, consistency of the analyses by a slab-gel system might reflect a further hindrance in this subtle procedure. Nevertheless, individual genotyping eliminates false positive results due to pooling artefacts and, in case of significant results, enables the thorough analyses of the marker alleles in detail (see Fig.1). In order to confirm the aforementioned positive results further markers (SNPs, Tab. 3) were genotyped located in the respective genes in the vicinity of the microsatellites representing TNFRSF17 and FLIP. Yet, RFLP analyses did not reveal any association of the selected SNPs and, therefore, the microsatellite data were not confirmed. On the other hand, these SNPs might not represent regions properly that encompass regulatory elements.
In some instances, the LD of distinct microsatellite alleles covers long genetic distances, thus further gene variations might be in linkage with these alleles. Since the significantly associated 'TNFRSF17' marker is located at the IBD8 region with 1MBp distance to the major histocompatibility class (MHC) II transactivator (MHCIITA), a previously reported functional variation of the MHC2TA gene was analysed (see Tab. 4; ). MHC2TA regulates the expression of human leukocyte antigen (HLA) genes regulating the adaptive immune system by presenting antigens to CD4+ T cells, thereby re-activating these cells. The HLA region has been implicated in IBD . In addition to the localisation of MHC2TA at IBD8 and the associated marker in the adjacent region, the putative biological relevance of the functional rs3087456 polymorphism for CD motivated us to genotype this variation. The analyses did reveal a marginal association in our CD patients when allele or genotype frequencies were compared between the combined control and patient cohorts (see Tab. 4). Yet there was no evidence for LD between TNFRSF17 and MHC2TA alleles. In order to validate these data further patient cohorts comprising more individuals must be scrutinised. In addition, other genes that might be linked with the 'TNFRSF17' marker must be analysed (at least 15 RefSeq genes in the region are encompassed by the microsatellite marker and MHC2TA).
In conclusion, this study did not reveal overt evidence for CD predisposition factors in apoptotic (and immune) pathways. Certainly, our approach depends on the LD between the investigated microsatellites and putative predisposing or protective alleles, depending on functional relevance to the disease. Thus, in some instances microsatellites might not be entirely representative for the adjacent genes. Furthermore, the investigated genes only cover part of the factors which coordinate programmed cell death. Yet, future information about haplotype blocks may facilitate more far-fetched interpretations of our analyses.
- Russel MG, Pastoor CJ, Janssen KM, van Deursen CT, Muris JW, van Wijlick EH, Stockbrugger RW: Familial aggregation of inflammatory bowel disease: a population-based study in South Limburg, The Netherlands. The South Limburg IBD Study Group. Scand J Gastroenterol Suppl. 1997, 223: 88-91.PubMedGoogle Scholar
- Hugot JP, Thomas G: Genome-wide scanning in inflammatory bowel diseases. Dig Dis. 1998, 16 (6): 364-369. 10.1159/000016893.View ArticlePubMedGoogle Scholar
- Hugot JP, Chamaillard M, Zouali H, Lesage S, Cezard JP, Belaiche J, Almer S, Tysk C, O'Morain CA, Gassull M: Association of NOD2 leucine-rich repeat variants with susceptibility to Crohn's disease. Nature. 2001, 411 (6837): 599-603. 10.1038/35079107.View ArticlePubMedGoogle Scholar
- Ogura Y, Bonen DK, Inohara N, Nicolae DL, Chen FF, Ramos R, Britton H, Moran T, Karaliuskas R, Duerr RH: A frameshift mutation in NOD2 associated with susceptibility to Crohn's disease. Nature. 2001, 411 (6837): 603-606. 10.1038/35079114.View ArticlePubMedGoogle Scholar
- Rogler G: Update in inflammatory bowel disease pathogenesis. Curr Opin Gastroenterol. 2004, 20 (4): 311-317. 10.1097/00001574-200407000-00003.View ArticlePubMedGoogle Scholar
- Franchimont D, Vermeire S, El Housni H, Pierik M, Van Steen K, Gustot T, Quertinmont E, Abramowicz M, Van Gossum A, Deviere J: Deficient host-bacteria interactions in inflammatory bowel disease? The toll-like receptor (TLR)-4 Asp299gly polymorphism is associated with Crohn's disease and ulcerative colitis. Gut. 2004, 53 (7): 987-992. 10.1136/gut.2003.030205.View ArticlePubMedPubMed CentralGoogle Scholar
- Klein W, Tromm A, Griga T, Folwaczny C, Hocke M, Eitner K, Marx M, Duerig N, Epplen JT: Interaction of polymorphisms in the CARD15 and CD14 genes in patients with Crohn disease. Scand J Gastroenterol. 2003, 38 (8): 834-836. 10.1080/00365520310003147.View ArticlePubMedGoogle Scholar
- Pallone F, Blanco Gdel V, Vavassori P, Monteleone I, Fina D, Monteleone G: Genetic and pathogenetic insights into inflammatory bowel disease. Curr Gastroenterol Rep. 2003, 5 (6): 487-492.View ArticlePubMedGoogle Scholar
- Okazawa A, Kanai T, Watanabe M, Yamazaki M, Inoue N, Ikeda M, Kurimoto M, Ishii H, Hibi T: Th1-mediated intestinal inflammation in Crohn's disease may be induced by activation of lamina propria lymphocytes through synergistic stimulation of interleukin-12 and interleukin-18 without T cell receptor engagement. Am J Gastroenterol. 2002, 97 (12): 3108-3117. 10.1111/j.1572-0241.2002.07107.x.View ArticlePubMedGoogle Scholar
- Singh UP, Singh S, Iqbal N, Weaver CT, McGhee JR, Lillard JW: IFN-gamma-inducible chemokines enhance adaptive immunity and colitis. J Interferon Cytokine Res. 2003, 23 (10): 591-600. 10.1089/107999003322485099.View ArticlePubMedGoogle Scholar
- Souza HS, Tortori CJ, Castelo-Branco MT, Carvalho AT, Margallo VS, Delgado CF, Dines I, Elia CC: Apoptosis in the intestinal mucosa of patients with inflammatory bowel disease: evidence of altered expression of FasL and perforin cytotoxic pathways. Int J Colorectal Dis. 2005, 20 (3): 277-286. 10.1007/s00384-004-0639-8.View ArticlePubMedGoogle Scholar
- Xia B, Yu YH, Guo QS, Li XY, Jiang L, Li J: Association of Fas-670 gene polymorphism with inflammatory bowel disease in Chinese patients. World J Gastroenterol. 2005, 11 (3): 415-417.View ArticlePubMedPubMed CentralGoogle Scholar
- Peppelenbosch MP, van Deventer SJ: T cell apoptosis and inflammatory bowel disease. Gut. 2004, 53 (11): 1556-1558. 10.1136/gut.2004.040824.View ArticlePubMedPubMed CentralGoogle Scholar
- Brannigan AE, O'Connell PR, Hurley H, O'Neill A, Brady HR, Fitzpatrick JM, Watson RW: Neutrophil apoptosis is delayed in patients with inflammatory bowel disease. Shock. 2000, 13 (5): 361-366.View ArticlePubMedGoogle Scholar
- Seidelin JB, Nielsen OH: [Apoptosis in chronic inflammatory bowel disease. The importance for pathogenesis and treatment]. Ugeskr Laeger. 2003, 165 (8): 790-792.PubMedGoogle Scholar
- Ina K, Itoh J, Fukushima K, Kusugami K, Yamaguchi T, Kyokane K, Imada A, Binion DG, Musso A, West GA: Resistance of Crohn's disease T cells to multiple apoptotic signals is associated with a Bcl-2/Bax mucosal imbalance. J Immunol. 1999, 163 (2): 1081-1090.PubMedGoogle Scholar
- Hommes DW, van Deventer SJ: Targeting tumor necrosis factor-alpha in inflammatory bowel disease: why, how, and when?. Curr Opin Gastroenterol. 2003, 19 (4): 350-357. 10.1097/00001574-200307000-00005.View ArticlePubMedGoogle Scholar
- ten Hove T, van Montfrans C, Peppelenbosch MP, van Deventer SJ: Infliximab treatment induces apoptosis of lamina propria T lymphocytes in Crohn's disease. Gut. 2002, 50 (2): 206-211. 10.1136/gut.50.2.206.View ArticlePubMedPubMed CentralGoogle Scholar
- Van den Brande JM, Braat H, van den Brink GR, Versteeg HH, Bauer CA, Hoedemaeker I, van Montfrans C, Hommes DW, Peppelenbosch MP, van Deventer SJ: Infliximab but not etanercept induces apoptosis in lamina propria T-lymphocytes from patients with Crohn's disease. Gastroenterology. 2003, 124 (7): 1774-1785. 10.1016/S0016-5085(03)00382-2.View ArticlePubMedGoogle Scholar
- Sands BE: Why Do Anti-Tumor Necrosis Factor Antibodies Work in Crohn's Disease?. Rev Gastroenterol Disord. 2004, 4 (Suppl. 3): S10-S17.PubMedGoogle Scholar
- Sturm A, Leite AZ, Danese S, Krivacic KA, West GA, Mohr S, Jacobberger JW, Fiocchi C: Divergent cell cycle kinetics underlie the distinct functional capacity of mucosal T cells in Crohn's disease and ulcerative colitis. Gut. 2004, 53 (11): 1624-1631. 10.1136/gut.2003.033613.View ArticlePubMedPubMed CentralGoogle Scholar
- Sartor RB: Current concepts of the etiology and pathogenesis of ulcerative colitis and Crohn's disease. Gastroenterol Clin North Am. 1995, 24 (3): 475-507.PubMedGoogle Scholar
- Doering J, Begue B, Lentze MJ, Rieux-Laucat F, Goulet O, Schmitz J, Cerf-Bensussan N, Ruemmele FM: Induction of T lymphocyte apoptosis by sulphasalazine in patients with Crohn's disease. Gut. 2004, 53 (11): 1632-1638. 10.1136/gut.2003.037911.View ArticlePubMedPubMed CentralGoogle Scholar
- Elson CO, Sartor RB, Tennyson GS, Riddell RH: Experimental models of inflammatory bowel disease. Gastroenterology. 1995, 109 (4): 1344-1367. 10.1016/0016-5085(95)90599-5.View ArticlePubMedGoogle Scholar
- Jagiello P, Gencik M, Arning L, Wieczorek S, Kunstmann E, Csernok E, Gross WL, Epplen JT: New genomic region for Wegener's granulomatosis as revealed by an extended association screen with 202 apoptosis-related genes. Hum Genet. 2004, 114 (5): 468-477. 10.1007/s00439-004-1092-z.View ArticlePubMedGoogle Scholar
- Wieczorek S, Jagiello P, Arning L, Dahmen N, Epplen JT: Screening for candidate gene regions in narcolepsy using a microsatellite based approach and pooled DNA. J Mol Med. 2004, 82 (10): 696-705. 10.1007/s00109-004-0569-5.View ArticlePubMedGoogle Scholar
- Klein W, Tromm A, Griga T, Fricke H, Folwaczny C, Hocke M, Eitner K, Marx M, Duerig N, Epplen JT: A polymorphism in the CD14 gene is associated with Crohn disease. Scand J Gastroenterol. 2002, 37 (2): 189-191. 10.1080/003655202753416867.View ArticlePubMedGoogle Scholar
- Klein W, Tromm A, Griga T, Fricke H, Folwaczny C, Hocke M, Eitner K, Marx M, Duerig N, Epplen JT: Interleukin-4 and interleukin-4 receptor gene polymorphisms in inflammatory bowel diseases. Genes Immun. 2001, 2 (5): 287-289. 10.1038/sj.gene.6363779.View ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.View ArticlePubMedPubMed CentralGoogle Scholar
- Swanberg M, Lidman O, Padyukov L, Eriksson P, Akesson E, Jagodic M, Lobell A, Khademi M, Borjesson O, Lindgren CM: MHC2TA is associated with differential MHC molecule expression and susceptibility to rheumatoid arthritis, multiple sclerosis and myocardial infarction. Nat Genet. 2005, 37 (5): 486-494. 10.1038/ng1544.View ArticlePubMedGoogle Scholar
- Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004, 74 (1): 106-120. 10.1086/381000.View ArticlePubMedGoogle Scholar
- Yap LM, Ahmad T, Jewell DP: The contribution of HLA genes to IBD susceptibility and phenotype. Best Pract Res Clin Gastroenterol. 2004, 18 (3): 577-596. 10.1016/j.bpg.2004.01.003.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.