- Open Access
- Open Peer Review
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
Journal of Negative Results in BioMedicine volume 4, Article number: 8 (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.
Intragenic microsatellite or markers located in the immediate vicinity (<50 kb) of the specific gene were included. Information on the oligonucleotide sequences and location of markers are given at the website (Additional file 1; see also Tab. 1). As reported before, only markers with equal "intra-subgroup" allele distributions with ≥ 2 alleles were considered in subsequent analyses . Significantly associated markers were genotyped individually in order to exclude false-positive results due to possible pooling artefacts. All in all, 245 microsatellite markers representing distinct genes were analysed on an ABI377 slab-gel system (Applied Biosystems, Darmstadt, Germany).
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.
Markers with significantly different allele distributions between patients and controls were controlled by genotyping individual DNA samples of patients and controls in order to exclude false-positive results due to pooling artefacts. Individual genotyping was performed by capillary gel electrophoresis by using the BeckmanCoulter CEQ8000 genetic analysis system (Beckman Coulter, Germany). Results were analysed by comparing each microsatellite allele frequency from the CD cohort with the corresponding allele frequency of the control group by χ2 testing and corrected by the number of marker specific alleles according to Bonferroni (see Tab. 2 and URL: http://www.ruhr-uni-bochum.de/mhg/marker_information_SENW.pdf). Hardy-Weinberg equilibrium (HWE) was tested using the Genepop program http://wbiomed.curtin.edu.au/genepop.
SNPs in genes as represented by significantly associated markers after individual genotyping were investigated by analysis of restriction fragment length polymorphisms (RFLP; see Tab. 3). As the marker representing the TNFRSF17 gene is located in ~1 MBp distance to the MHC class II transactivator (MHC2TA) gene, a functional variation (rs3087456, ) of MHC2TA was genotyped by RFLP analyses in 147 CD patients and 463 healthy controls from the abovementioned control populations (see Tab. 3). The results were evaluated by means of χ2 -and HWE testing. Linkage disequilibrium (LD) between the marker alleles and the polymorphism was calculated by the Genepop program.
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.
Genotypes including allele 2 are overrepresented in the control cohort, whereas those with the apparently predisposing allele 3 are more frequent in the CD cohort, thus confirming the results of individual genotyping (see Fig. 1).
The analyses of the functionally significant polymorphism rs3087456 revealed a marginal association in our CD patients when allele or genotype frequencies were compared between the combined control (WeG and NoG did not differ in allele frequencies) and the patient cohorts (see Tab. 4). Analyses for LD between TNFRSF17 and MHC2TA alleles, however, did not reveal any significant deviations from equilibrium.
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.
Hugot JP, Thomas G: Genome-wide scanning in inflammatory bowel diseases. Dig Dis. 1998, 16 (6): 364-369. 10.1159/000016893.
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.
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.
Rogler G: Update in inflammatory bowel disease pathogenesis. Curr Opin Gastroenterol. 2004, 20 (4): 311-317. 10.1097/00001574-200407000-00003.
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.
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.
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.
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.
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.
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.
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.
Peppelenbosch MP, van Deventer SJ: T cell apoptosis and inflammatory bowel disease. Gut. 2004, 53 (11): 1556-1558. 10.1136/gut.2004.040824.
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.
Seidelin JB, Nielsen OH: [Apoptosis in chronic inflammatory bowel disease. The importance for pathogenesis and treatment]. Ugeskr Laeger. 2003, 165 (8): 790-792.
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.
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.
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.
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.
Sands BE: Why Do Anti-Tumor Necrosis Factor Antibodies Work in Crohn's Disease?. Rev Gastroenterol Disord. 2004, 4 (Suppl. 3): S10-S17.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Electronic supplementary material
Additional File 1: This file provides detailed information on the sequence of used oligonucleotides, represented gene, marker distance to gene and kind of nucleotide repeat (di, tri, etc.). Furthermore, the file includes graphical information on individually genotyped microsatellites markers with significant differences in allele distributions. (PDF 39 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.