Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(1): e54338.
Published online 2013 January 24. doi:  10.1371/journal.pone.0054338
PMCID: PMC3554777

IRGM Variants and Susceptibility to Inflammatory Bowel Disease in the German Population

David L. Boone, Editor


Background & Aims

Genome-wide association studies identified the autophagy gene IRGM to be strongly associated with Crohn's disease (CD) but its impact in ulcerative colitis (UC), its phenotypic effects and potential epistatic interactions with other IBD susceptibility genes are less clear which we therefore analyzed in this study.

Methodology/Principal Findings

Genomic DNA from 2060 individuals including 817 CD patients, 283 UC patients, and 961 healthy, unrelated controls (all of Caucasian origin) was analyzed for six IRGM single nucleotide polymorphisms (SNPs) (rs13371189, rs10065172 = p.Leu105Leu, rs4958847, rs1000113, rs11747270, rs931058). In all patients, a detailed genotype-phenotype analysis and testing for epistasis with the three major CD susceptibility genes NOD2, IL23R and ATG16L1 were performed.

Our analysis revealed an association of the IRGM SNPs rs13371189 (p = 0.02, OR 1.31 [95% CI 1.05–1.65]), rs10065172 = p.Leu105Leu (p = 0.016, OR 1.33 [95% CI 1.06–1.66]) and rs1000113 (p = 0.047, OR 1.27 [95% CI 1.01–1.61]) with CD susceptibility. There was linkage disequilibrium between these three IRGM SNPs. In UC, several IRGM haplotypes were weakly associated with UC susceptibility (p<0.05). Genotype-phenotype analysis revealed no significant associations with a specific IBD phenotype or ileal CD involvement. There was evidence for weak gene-gene-interaction between several SNPs of the autophagy genes IRGM and ATG16L1 (p<0.05), which, however, did not remain significant after Bonferroni correction.


Our results confirm IRGM as susceptibility gene for CD in the German population, supporting a role for the autophagy genes IRGM and ATG16L1 in the pathogenesis of CD.


Crohn's disease (CD) and ulcerative colitis (UC) are chronic inflammatory bowel diseases (IBD) resulting from an inappropriate immune response to microbial antigens in genetically susceptible individuals [1], [2], [3], [4]. Recent genome-wide association studies (GWAS) have provided valuable insights into the genetic architecture particularly of CD, identifying more than 70 CD susceptibility variants with the most significant findings in the gene regions of NOD2, IL23R and ATG16L1 [5], [6], [7], [8], [9]. These genetic findings confirm an important role for innate immunity, proinflammatory IL-23/Th17 immune responses as well as autophagy for both gut homeostasis and the development of chronic inflammation in IBD.

In addition to the CD susceptibility gene ATG16L1 involved in autophagy [9], [10], [11], recent GWAS identified the single nucleotide polymorphism (SNP) rs13361189 – a SNP lying immediately upstream of the autophagy gene IRGM (immunity-related GTPase family M) – and other IRGM SNPs to be strongly associated with CD [12], [13]. Since the discovery of the IRGM as a CD susceptibility gene, further studies have investigated IRGM gene variants in both adult and pediatric CD [14], [15], [16], [17], [18], [19], [20], [21], [22] as well as in UC [14], confirming its role in the IBD pathogenesis.

A functional study suggested a common, 20-kb deletion polymorphism upstream of IRGM, which is in perfect linkage disequilibrium (LD) with rs13361189, as a likely causal variant, since the deletion allele modulated the expression of IRGM in transformed cells [22]. Another recent study implicated a variant in the 5′-untranslated region (−308(GTTT)5) to be independently associated with CD [20], while the very recent study by Brest et al. [21] demonstrated functional effects of the synonymous SNP rs10065172 (c.313C>T), which is also in linkage disequilibrium with rs13361189 and the deletion polymorphism. This exonic, synonymous variant rs10065172 in IRGM alters a binding site for certain microRNAs (miR-196) and causes deregulation of IRGM-dependent xenophagy of bacteria in patients with CD [21], therefore suggesting rs10065172 as disease-causing variant.

These studies implicate that autophagy plays an important role in human inflammatory disorders by direct elimination of intracellular bacteria and activation of pattern recognition receptor (PRR) signaling which is involved in gut homeostasis and CD pathogenesis [10], [23]. The IRGM gene belongs to immunity-related GTPases (IRG), a family of genes in mammalian species induced by interferons (IFNs) and functioning as key mediators of IFN-regulated resistance to intracellular bacteria and protozoa [23]. IRGM has been shown to play a role in the autophagy-targeted destruction of Mycobacterium bovis BCG [23] and the IFN-γ-induced host defense against Salmonella typhimurium infection [10]. Interestingly, a recent study in CD patients demonstrated that autophagy limits the replication of intracellular adherent-invasive Escherichia coli (AIEC) associated with ileal CD and that IRGM- and ATG16L1-deficient cells had enhanced intracellular AIEC bacteria replication, suggesting a significant impact on the outcome of intestinal inflammation [24].

While several GWAS and replication studies established IRGM as a CD susceptibility gene, its effects on the IBD phenotype are less clear. In addition, epistatic interactions with other IBD susceptibility genes, in particular the second autophagy gene ATG16L1, have not been studied in detail. Therefore, in this study, we aimed to analyze the role of IRGM on CD and UC susceptibility as well as its effect on the IBD phenotype in a large patient-control cohort. In addition, we performed a detailed epistasis analysis of IRGM with the three major CD susceptibility genes NOD2, ATG16L1 and IL23R. In total, six major IRGM SNPs, for which associations with CD were shown in previous studies (see details in Methods), were genotyped in more than 2000 German IBD patients and controls.

Patients and Methods

Ethics statement

The study was approved by the Ethics committee of the Medical Faculty of the Ludwig-Maximilians-University Munich. Written, informed consent was obtained from all patients prior to the study. Study protocols were based on the ethical principles for medical research involving human subjects of the Helsinki Declaration (

Study population and definition of IBD phenotype

The study population (n = 2060) consisted of 1099 IBD patients including 817 patients with CD, 283 patients with UC, and 961 healthy, unrelated controls, all of Caucasian origin. Patient charts were analyzed for demographic and clinical parameters (disease behaviour and anatomic location of IBD, disease-related complications, history of surgery or immunosuppressive therapy) and all patients participated in a detailed questionnaire including an interview at time of enrolment. The diagnosis of CD or UC was determined according to endoscopic, histopathologic and radiological criteria of current international guidelines [25]. Patients with clinical features of both CD and UC (and therefore classified as “indeterminate colitis”) were excluded from this study. Patients with CD were assessed based on the Montreal classification including age at diagnosis (A), location (L), and behaviour (B) of disease [26]. In patients with UC, anatomic location was also assessed following the Montreal classification analyzing the criteria ulcerative proctitis (E1), left-sided UC (distal UC; E2), and extensive UC (pancolitis; E3) [26]. The demographic baseline characteristics of the study population were collected blind to the results of the genotype analyses and are summarized in Table 1.

Table 1
Demographic characteristics of the IBD study population.

DNA extraction

From all study participants, blood samples were taken and genomic DNA was isolated from peripheral blood leukocytes using the DNA blood mini kit from Qiagen (Hilden, Germany) according to the manufacturer's guidelines.

Genotyping of the IRGM variants

Six IRGM SNPs (rs13361189, rs10065172 = pLeu105Leu, rs4958847, rs1000113, rs931058, rs11747270) were genotyped. The selection of these six IRGM SNPs was based on previous studies showing associations for these SNPs in large case-control cohorts. The SNPs rs13361189, rs10065172 = p.Leu105Leu and rs4958847 were selected from the study of Parkes et al. [12], while the SNPs rs1000113 and rs931058 were tested in the study of the Wellcome Trust Case Control Consortium (WTCCC) [13]. The SNPs rs13361189 and rs10065172 = p.Leu105Leu served also as proxies for a common, 20-kb deletion polymorphism upstream of IRGM, since they are in perfect linkage disequilibrium (r2 = 1.0) with this deletion polymorphism [22]. Additionally rs11747270, which was the most strongly CD-associated SNP within the IRGM region in the meta-analysis of Barrett et al., was included.

IRGM genotyping was performed by PCR and melting curve analysis using a pair of fluorescence resonance energy transfer (FRET) probes in a LightCycler® 480 Instrument (Roche Diagnostics, Mannheim, Germany) as previously described in detail [27], [28], [29], [30], [31]. The donor fluorescent molecule (fluorescein) at the 3′-end of the sensor probe is excited at its specific fluorescence excitation wavelength (533 nm) and the energy is transferred to the acceptor fluorescent molecule at the 5′-end (LightCycler Red 610, 640 or 670) of the anchor probe. The specific fluorescence signal emitted by the acceptor molecule is detected by the optical unit of the LightCycler 480 instrument. The sensor probe is exactly matching to one allele of each SNP, preferentially to the rarer allele, whereas in the case of the other allele there is a mismatch resulting in a lower melting temperature. The total volume of the PCR was 5 µl containing 25 ng of genomic DNA, 1× Light Cycler 480 Genotyping Master (Roche Diagnostics), 2.5 pmol of each primer and 0.75 pmol of each FRET probe (TIB MOLBIOL, Berlin, Germany). In the case of rs11747270, the amount of the forward primer was reduced to one fifth and in the case of rs4958847 the reverse primer was reduced to one half. In the case of rs10065172 and rs931058, the reverse primers were reduced to one third, respectively. The PCR comprised an initial denaturation step (95°C for 10 min) and 45 cycles (50 cycles in the case of rs10065172) [95°C for 10°C sec, 60°C (55°C in the case of rs10065172) for 10 sec, 72°C for 15 sec]. The melting curve analysis comprised an initial denaturation step (95°C for 1 min), a step rapidly lowering the temperature to 40°C and holding for 2 min, and a heating step slowly (1 acquisition/°C) increasing the temperature up to 95°C and continuously measuring the fluorescence intensity. The results of melting curve analysis have been confirmed by analyzing two patient samples for each possible genotype using sequence analysis. For sequencing, the total volume of the PCR was 100 µl containing 250 ng of genomic DNA, 1× PCR buffer (Qiagen, Hilden, Germany), a final MgCl2 concentration of 2 mM, 0.5 mM of a dNTP mix (Sigma, Steinheim, Germany), 2.5 units of HotStar Plus Taq™ DNA polymerase (Qiagen) and 10 pmol of each primer (TIB MOLBIOL). The PCR comprised an initial denaturation step (95°C for 5 min), 35 cycles (denaturation at 94°C for 30 sec, primer annealing at 60°C for 30 sec, extension at 72°C for 30 sec) and a final extension step (72°C for 10 min). The PCR products were purified using the QIAquick PCR Purification Kit (Qiagen) and sequenced by a commercial sequencing company (Sequiserve, Vaterstetten, Germany). All sequences of primers and FRET probes used for genotyping and for sequence analysis are given in Tables S1 and S2.

Genotyping of NOD2, IL23R and ATG16L1 variants

For analysis of potential epistatic interactions, the genotypes of gene variants in NOD2, IL23R and ATG16L1 were available from previous studies for all patients and controls analyzed [11], [27], [32], [33], [34], [35], [36]. Genotyping of the NOD2 variants p.Arg702Trp (rs2066844), p.Gly908Arg (rs2066847), and p.Leu1007fsX1008 (rs2066847) was performed as described previously (primer sequences available on request) [34]. The 10 IL23R SNPs (rs1004819, rs7517847, rs10489629, rs2201841, rs11465804, rs11209026 = p.Arg381Gln, rs1343151, rs10889677, rs11209032, rs1495965) and nine ATG16L1 SNPs (rs13412102, rs12471449, rs6431660, rs1441090, rs2289472, rs2241880 ( = T300A), rs2241879, rs3792106, rs4663396) were genotyped by PCR and melting curve analysis as described previously [11], [35].

Statistical analyses

For evaluation of data, the SPSS 13.0 software (SPSS Inc., Chicago, IL, USA) and R-2.13.1 ( were used. Each genetic marker was tested for Hardy-Weinberg equilibrium in the control group. Fisher's exact test was used for comparison between categorical variables and Student's t test was applied for quantitative variables. All tests were two-tailed and p-values<0.05 were considered as significant. Odds ratios were calculated for the minor allele of each SNP. Correction for multiple testing was performed by Bonferroni correction where indicated. Haplotype analysis was calculated using the –hap-logistic command in PLINK (, epistasis analysis was performed with the –epistasis option. LD between SNPs was evaluated using the R-library genetics. Genotype-phenotype associations were also tested in R using logistic regression.


The IRGM gene variants are associated with susceptibility to CD

The allele frequencies of the SNPs rs13371189, rs10065172 = p.Leu105Leu, rs4958847, rs11747270, rs931058 and rs1000113 of all three subgroups (CD, UC, and controls) were in accordance with the predicted Hardy-Weinberg equilibrium and are summarized in Table 2. Overall, our analysis revealed an association of the IRGM variants rs13371189 (p = 0.02, OR 1.31 [95% CI 1.05–1.65]), rs10065172 = p.Leu105Leu (p = 0.016, OR 1.33 [95% CI 1.06–1.66]) and rs1000113 (p = 0.047, OR 1.27 [95% CI 1.01–1.61]) with the susceptibility to CD. Similar to previous studies, rs13371189 and rs10065172 = p.Leu105Leu were in perfect linkage disequilibrium (r2≈1.0) in all three subgroups (CD, UC, controls; Tables S3, S4, S5). Strong linkage disequilibrium was also shown for these two SNPs with the third CD-associated IRGM SNP rs1000113 (Tables S3, S4, S5). With exception of rs11747270, none of the genotyped IRGM SNPs was associated with UC susceptibility (Table 2).

Table 2
Associations of IRGM gene markers in CD and UC case-control association studies.

IRGM haplotype analysis

Next, we performed a detailed haplotype analysis investigating the role of IRGM haplotypes on CD and UC susceptibility. As demonstrated in tables 3 and and4,4, several IRGM haplotypes demonstrated an association with CD and UC susceptibility. In CD, the strongest associations were found for haplotypes containing at least one of the most strongly CD-associated SNP rs13361189 or rs10065172 (Table 3), while in UC, the strongest association was found for rs11747270-rs931058 (omnibus p-value 1.57×10−2) (Table 4). However, given the large number of haplotypes analyzed, none of these associations withstood Bonferroni correction for multiple testing.

Table 3
Haplotypes of IRGM SNPs in CD case-control sample and omnibus p-values for association with CD susceptibility.
Table 4
Haplotypes of IRGM SNPs in UC case-control sample and omnibus p-values for association with UC susceptibility.

Genotype-phenotype analysis

We further investigated whether IRGM SNPs are associated with certain phenotypic characteristics in IBD patients. Based on the Montreal classification of IBD, the phenotypic data of IBD patients were analyzed for anatomic localization. However, none of the IRGM SNPs investigated were associated with specific disease localization in CD (Table S6) or UC (Table S7). Moreover, a detailed genotype-phenotype analysis in CD patients of the exonic synonymous SNP rs10065172 = p.Leu105Leu, which was in linkage disequilibrium with rs13361189 and with the previously identified 20-kb deletion polymorphism immediately upstream of IRGM (r2 = 1.0), did not reveal any significant associations with the CD phenotype (Table S8).

Analysis for epistasis of IRGM with other major CD susceptibility genes

Finally, we analyzed potential evidence for gene-gene interactions of IRGM variants with other CD susceptibility genes such as variants in the NOD2, IL23R and ATG16L1 gene including their effect on CD susceptibility. Interestingly, there was evidence for weak gene-gene-interaction between several SNPs of the two autophagy genes IRGM and ATG16L1 (ATG16L1 rs12471449, ATG16L1 rs1441090, ATG16L1 rs4663396), which, however, did not remain significant after Bonferroni correction (Table 5). The odds ratios of gene-gene interactions, which were significant before Bonferroni correction, are given in Table 6. There was no epistasis between IRGM and the other two major CD susceptibility genes NOD2 and IL23R.

Table 5
Analysis for gene-gene interaction (epistasis) of IRGM SNPs with NOD2, ATGT16L1, and IL23R gene variants regarding CD susceptibility.
Table 6
Odds ratios (ORs) and 95% CI for the gene-gene interaction (epistasis) found to be significant (before Bonferroni correction) for CD susceptibility (shown in Table 5).


This study represents a detailed analysis of IRGM gene variants regarding their role in the susceptibility and phenotype of IBD in a large cohort of more than 2000 Caucasian individuals. In line with previous GWAS and replication studies [12], [13], [14], [15], [16], [17], [18], [19], our results confirm an association of the IRGM variant rs13371189 with CD susceptibility. A detailed functional study identified a deletion polymorphism directly upstream of the IRGM locus as a candidate SNP to explain the CD association at this locus [22], affecting the tissue-specific expression level of IRGM [37]. This 20-kb deletion polymorphism is in perfect linkage disequilibrium (r2 = 1.0) with SNP rs13361189, therefore implicating that rs13361189 is a proxy for this deletion polymorphism.

Similar to the study by McCarroll et al. [22], we demonstrate that the common exonic synonymous SNP rs10065172 = p.Leu105Leu is in linkage disequilibrium with rs13361189 and therefore also with the previously identified 20-kb deletion polymorphism. The exonic SNP rs10065172 (c.313C>T) has been previously classified as non-causative given the absence of an alteration in the IRGM protein sequence or splice sites, although this view is challenged by the results of a very recent study [21]. The study by Brest et al. demonstrated that a family of microRNAs (miRNAs), miR-196, is overexpressed in the inflamed intestinal epithelium of CD patients and downregulates the IRGM protective variant (c.313C) but not the CD-associated allele (c.313T) [21]. The same study demonstrated that the resulting loss of tight regulation of IRGM expression compromises the control of the intracellular replication of CD-associated adherent invasive Escherichia coli (AIEC) by affecting the efficacy of bacterial phagocytosis (xenophagy) [21]. Therefore, Brest et al. [21] suggest the synonymous SNP rs10065172 (c.313C>T) as a likely causal variant. rs10065172 has been also shown to be associated with susceptibility to tuberculosis [38] which is of interest, given evidence that certain mycobacteria may play a role in the pathogenesis of CD. Moreover, a functional study demonstrated that IRGM induces autophagy to eliminate intracellular mycobacteria [23].

Overall, the association signal of IRGM with CD found in our study was considerably weaker than that shown by us for the other autophagy gene ATG16L1 in a similar sized cohort [11]. Similarly, the recent CD meta-analyses showed a stronger association signal for ATG16L1 than for IRGM [5]. Anderson et al. performed a very large meta-analysis of CD and UC associated susceptibility loci [6]. In this analysis, the CD case-control cohort included n = 6,333 CD patients and n = 15,056 controls, while the UC case-control cohort consisted of n = 6,687 UC patients and 19,718 controls [6]. In this large meta-analysis, they demonstrated that the IRGM SNP rs7714584 is associated with both CD (p = 7.76×10−19; OR 1.37, 95% CI 1.28–1.47) and UC (p = 3.95×10−4; OR 1.14, 95% CI 1.06–1.22) [6]. However, the authors defined only SNPs with p-values of <1×10−4 to be significantly associated with both CD and UC and therefore included IRGM not in the list of susceptibility loci shared between CD and UC [6]. In addition, the meta-analysis of Palomino-Morales et al. demonstrated also an association of two IRGM SNPs with UC (rs13361189 p = 0.0069, pooled OR = 1.16; rs4958847 p = 0.014, pooled OR = 1.13) [14]. Most importantly, a recent very large IBD meta-analysis comprising 75,000 IBD patients and controls confirmed the IRGM gene region to be strongly associated with CD (p = 2.94×10−37) and to a much lesser degree with UC (p = 0.0025 for the UC GWAS cohort and p = 1.05×10−7 for the combined UC Immunochip and GWAS cohort) [39]. Therefore, based on these meta-analyses as well as on our results showing a trend for association with UC for the IRGM SNP rs11747270 (p = 0.084, OR 1.42 [0.96–2.12] for comparing minor allele frequencies; p = 0.033, OR 5.34 [1.03–34.64] for comparing genotype frequencies; Table 2) and associations of several IRGM haplotypes with UC (p<0.05), IRGM can be regarded to be weakly associated with UC and as a shared susceptibility gene of both UC and CD, although it has a much more prominent role in the pathogenesis of CD.

In addition, we also performed a detailed genotype-phenotype analysis of IRGM variants in CD and UC patients. In contrast to a recent study of Latiano et al. [40] demonstrating an association of IRGM variants with fistulizing CD, our genotype-phenotype analysis did not reveal any significant association of IRGM variants with the CD phenotype. We were also unable to confirm an association with ileal CD found in a previous study of a smaller CD cohort from New Zealand [18]. Our findings may be related to the rather weak association signal found for IRGM in the German CD cohort, although the results of this genotype-phenotype analysis are consistent with the lack of a well-defined phenotype in CD patients carrying risk alleles of the other autophagy gene ATG16L1 [11].

The identification of the two major CD susceptibility genes ATG16L1 and IRGM involved in autophagy has significantly strengthened the importance of autophagy and bacterial xenophagy in the complex and multifactorial etiology of IBD. However, potential epistatic interactions between ATG16L1 and IRGM have not been investigated in detail so far. We therefore analyzed epistasis between these two genes, demonstrating a weak gene-gene-interaction between several SNPs of the two autophagy genes IRGM and ATG16L1 which, however, did not remain significant after Bonferroni correction. Given their close functional relationship, this potential epistasis signal is highly interesting. Very recently, the largest IBD meta-analysis published so far (including 75,000 IBD patients and controls) was made publicly available [39]. Part of this meta-analysis was an epistasis analysis in IBD, UC and CD datasets of the Immunochip study. While the analyses of the CD and UC subsets were inconclusive, the results for the analysis with IBD showed only one suggestive result between SNPs near SLC7A10 (rs17694108) and IL2RA (rs12722515) with a p-value of 3.26×10−5 [39]. Therefore, the weak gene-gene interaction found between IRGM and ATG16L1 regarding CD susceptibility in our study, which was not significant after Bonferroni correction, could not been replicated on a significant level in this very large CD cohort. Thus, it is unlikely that epistasis between the two major autophagy genes contributes significantly to CD susceptibility.

There is increasing evidence for important intersections of autophagy and intracellular bacterial sensing (demonstrated by the importance of NOD2 in autophagy induction [41], [42]) in the pathogenesis of IBD. Moreover, recent studies identified a new pathway closely linked to autophagy and innate immunity, which is characterized by an unfolded protein response, stimulated by endoplasmic reticulum (ER) stress due to the accumulation of misfolded proteins. Several genes involved in ER stress, including XBP1 and ORMDL3 have been linked to the IBD pathogenesis on a genetic level [43], [44]. Interestingly, ATG16L1, NOD2, and XBP1 have been also demonstrated to affect the function of Paneth cells [43], [45], [46], suggesting a central role for this cell type in the development of IBD.

These recent findings are in line with raising evidence that NOD2 is involved in regulation of autophagy. Dendritic cells from CD patients expressing CD-associated NOD2 or ATG16L1 variants were shown to be defective in autophagy induction, bacterial trafficking and antigen presentation [41]. Most interestingly, a recent study demonstrated that the intracellular sensors NOD1 and NOD2 are critical for the autophagic response to invasive bacteria by recruiting the autophagy protein ATG16L1 to the plasma membrane at the bacterial entry site [42]. In cells homozygous for the CD-associated NOD2 frameshift mutation (p.Leu1007fsX1008), mutant NOD2 failed to recruit ATG16L1 to the plasma membrane and wrapping of invading bacteria by autophagosomes was impaired [42]. This is of particular interest, since we previously demonstrated a very severe stricturing phenotype in CD patients homozygous for the NOD2 p.Leu1007fsX1008 mutation associated with early disease onset, ileal stenosis, recurrent need for surgery and increased prevalence of entero-enteral fistulae [32], [33]. However, despite the central functional role of NOD2 in the induction of autophagic processes, our study could not demonstrate gene-gene-interactions between NOD2 and IRGM regarding CD susceptibility. Moreover, we could not identify significant epistatic interactions between IRGM and IL23R, the main IBD susceptibility gene involved in Th17 responses. Of interest, a very recent study demonstrated IL23R variants as susceptibility variants for leprosy and suggested a potential involvement of IL23R in the autophagocytosis of mycobacteria involved in the pathogenesis of leprosy [47].

In conclusion, our results confirm IRGM as susceptibility gene for CD in the German population, while we did not show an association with a specific IBD subphenotype. The strongest association signals for CD susceptibility were found for rs13361189 (proxy for the common, 20-kb deletion polymorphism upstream of IRGM) and the exonic synonymous SNP rs10065172 = p.Leu105Leu, supporting previous functional studies that these two SNPs may be the causal variants. However, the strength of the association signal with CD found here was several log-fold weaker than that demonstrated by us for the second autophagy gene ATG16L1 [11], suggesting a more important role for ATG16L1 in the CD pathogenesis. In UC, several IRGM haplotypes were weakly associated with UC susceptibility. This is consistent with recent meta-analyses which found weak associations with UC but very strong disease associations with CD. One might therefore hypothesize that autophagy genes such as IRGM and ATG16L1 play a more important role in the susceptibility to CD than UC. The potential epistasis signal between IRGM and ATG16L1 regarding CD susceptibility found in this study is highly interesting but could not been confirmed in a very large recent IBD meta-analysis [39] arguing against a major role of epistasis between IRGM and ATG16L1 regarding CD susceptibility.

Supporting Information

Table S1

Primer sequences and FRET probe sequences used for genotyping IRGM variants.


Table S2

Primer sequences used for the sequence analysis of IRGM variants.


Table S3

Analysis for linkage disequilibrium in CD patients. Values are given as r2/D′-measurements.


Table S4

Analysis for linkage disequilibrium in UC patients. Values are given as r2/D′-measurements.


Table S5

Analysis for linkage disequilibrium in controls. Values are given as r2/D′-measurements.


Table S6

P-values for allelic association of IRGM gene markers with the anatomic location of Crohn's disease (CD) according to the Montreal classification.


Table S7

P-values for allelic associations of IRGM gene markers with the anatomic location of ulcerative colitis (UC) according to the Montreal classification.


Table S8

Genotype-phenotype-analysis of the exonic synonymous SNP rs10065172 = p.Leu105Leu.


Funding Statement

JG was supported by a grant from the Broad Medical Foundation (IBD-0126R2). J. Seiderer and JD were supported by grants from the Ludwig-Maximilians-University Munich (FöFoLe Nr. 422; Habilitationsstipendium, LMUExcellent to J. Seiderer and Promotionsstipendium to JD); J. Seiderer was also supported by the Else Kröner-Fresenius-Stiftung (81/08//EKMS08/01). S. Brand was supported by grants from the German Research Foundation (BR 1912/6-1), the Else Kröner-Fresenius-Stiftung (Else Kröner Exzellenzstipendium 2010; 2010_EKES.32), and by grants of Ludwig-Maximilians-University Munich (Excellence Initiative, Investment Funds 2008 and FöFoLe program). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Khor B, Gardet A, Xavier RJ (2011) Genetics and pathogenesis of inflammatory bowel disease. Nature 474: 307–317 [PMC free article] [PubMed]
2. Brand S (2009) Crohn's disease: Th1, Th17 or both? The change of a paradigm: new immunological and genetic insights implicate Th17 cells in the pathogenesis of Crohn's disease. Gut 58: 1152–1167 [PubMed]
3. Danese S, Fiocchi C (2011) Ulcerative colitis. N Engl J Med 365: 1713–1725 [PubMed]
4. Abraham C, Cho JH (2009) Inflammatory bowel disease. N Engl J Med 361: 2066–2078 [PMC free article] [PubMed]
5. Franke A, McGovern DP, Barrett JC, Wang K, Radford-Smith GL, et al. (2010) Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci. Nat Genet 42: 1118–1125 [PMC free article] [PubMed]
6. Anderson CA, Boucher G, Lees CW, Franke A, D'Amato M, et al. (2011) Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47. Nat Genet 43: 246–252 [PMC free article] [PubMed]
7. Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr RH, et al. (2008) Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease. Nat Genet 40: 955–962 [PMC free article] [PubMed]
8. Duerr RH, Taylor KD, Brant SR, Rioux JD, Silverberg MS, et al. (2006) A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science 314: 1461–1463 [PubMed]
9. Hampe J, Franke A, Rosenstiel P, Till A, Teuber M, et al. (2007) A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1. Nat Genet 39: 207–211 [PubMed]
10. Rioux JD, Xavier RJ, Taylor KD, Silverberg MS, Goyette P, et al. (2007) Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis. Nat Genet 39: 596–604 [PMC free article] [PubMed]
11. Glas J, Konrad A, Schmechel S, Dambacher J, Seiderer J, et al. (2008) The ATG16L1 gene variants rs2241879 and rs2241880 (T300A) are strongly associated with susceptibility to Crohn's disease in the German population. Am J Gastroenterol 103: 682–691 [PubMed]
12. Parkes M, Barrett JC, Prescott NJ, Tremelling M, Anderson CA, et al. (2007) Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility. Nat Genet 39: 830–832 [PMC free article] [PubMed]
13. WTCC Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678 [PMC free article] [PubMed]
14. Palomino-Morales RJ, Oliver J, Gomez-Garcia M, Lopez-Nevot MA, Rodrigo L, et al. (2009) Association of ATG16L1 and IRGM genes polymorphisms with inflammatory bowel disease: a meta-analysis approach. Genes Immun 10: 356–364 [PubMed]
15. Weersma RK, Stokkers PC, Cleynen I, Wolfkamp SC, Henckaerts L, et al. (2009) Confirmation of multiple Crohn's disease susceptibility loci in a large Dutch-Belgian cohort. Am J Gastroenterol 104: 630–638 [PubMed]
16. Van Limbergen J, Russell RK, Nimmo ER, Drummond HE, Davies G, et al. (2009) Germline variants of IRGM in childhood-onset Crohn's disease. Gut 58: 610–611 [PubMed]
17. Amre DK, Mack DR, Morgan K, Krupoves A, Costea I, et al. (2009) Autophagy gene ATG16L1 but not IRGM is associated with Crohn's disease in Canadian children. Inflamm Bowel Dis 15: 501–507 [PubMed]
18. Roberts RL, Hollis-Moffatt JE, Gearry RB, Kennedy MA, Barclay ML, et al. (2008) Confirmation of association of IRGM and NCF4 with ileal Crohn's disease in a population-based cohort. Genes Immun 9: 561–565 [PubMed]
19. Fisher SA, Tremelling M, Anderson CA, Gwilliam R, Bumpstead S, et al. (2008) Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease. Nat Genet 40: 710–712 [PMC free article] [PubMed]
20. Prescott NJ, Dominy KM, Kubo M, Lewis CM, Fisher SA, et al. (2010) Independent and population-specific association of risk variants at the IRGM locus with Crohn's disease. Hum Mol Genet 19: 1828–1839 [PMC free article] [PubMed]
21. Brest P, Lapaquette P, Souidi M, Lebrigand K, Cesaro A, et al. (2011) A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn's disease. Nat Genet 43: 242–245 [PubMed]
22. McCarroll SA, Huett A, Kuballa P, Chilewski SD, Landry A, et al. (2008) Deletion polymorphism upstream of IRGM associated with altered IRGM expression and Crohn's disease. Nat Genet 40: 1107–1112 [PMC free article] [PubMed]
23. Singh SB, Davis AS, Taylor GA, Deretic V (2006) Human IRGM induces autophagy to eliminate intracellular mycobacteria. Science 313: 1438–1441 [PubMed]
24. Lapaquette P, Glasser AL, Huett A, Xavier RJ, Darfeuille-Michaud A (2010) Crohn's disease-associated adherent-invasive E. coli are selectively favoured by impaired autophagy to replicate intracellularly. Cell Microbiol 12: 99–113 [PubMed]
25. Lennard-Jones JE (1989) Classification of inflammatory bowel disease. Scand J Gastroenterol Suppl 170: 2–6 discussion 16–19. [PubMed]
26. Silverberg MS, Satsangi J, Ahmad T, Arnott ID, Bernstein CN, et al. (2005) 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 19 Suppl A: 5–36 [PubMed]
27. Glas J, Seiderer J, Fries C, Tillack C, Pfennig S, et al. (2011) CEACAM6 gene variants in inflammatory bowel disease. PLoS One 6: e19319. [PMC free article] [PubMed]
28. Glas J, Seiderer J, Nagy M, Fries C, Beigel F, et al. (2010) Evidence for STAT4 as a common autoimmune gene: rs7574865 is associated with colonic Crohn's disease and early disease onset. PLoS One 5: e10373. [PMC free article] [PubMed]
29. Glas J, Seiderer J, Pasciuto G, Tillack C, Diegelmann J, et al. (2009) rs224136 on chromosome 10q21.1 and variants in PHOX2B, NCF4, and FAM92B are not major genetic risk factors for susceptibility to Crohn's disease in the German population. Am J Gastroenterol 104: 665–672 [PubMed]
30. Glas J, Seiderer J, Fischer D, Tengler B, Pfennig S, et al. (2011) Pregnane X receptor (PXR/NR1I2) gene haplotypes modulate susceptibility to inflammatory bowel disease. Inflamm Bowel Dis 17: 1917–1924 [PubMed]
31. Glas J, Stallhofer J, Ripke S, Wetzke M, Pfennig S, et al. (2009) Novel genetic risk markers for ulcerative colitis in the IL2/IL21 region are in epistasis with IL23R and suggest a common genetic background for ulcerative colitis and celiac disease. Am J Gastroenterol 104: 1737–1744 [PubMed]
32. Seiderer J, Brand S, Herrmann KA, Schnitzler F, Hatz R, et al. (2006) Predictive value of the CARD15 variant 1007fs for the diagnosis of intestinal stenoses and the need for surgery in Crohn's disease in clinical practice: results of a prospective study. Inflamm Bowel Dis 12: 1114–1121 [PubMed]
33. Seiderer J, Schnitzler F, Brand S, Staudinger T, Pfennig S, et al. (2006) Homozygosity for the CARD15 frameshift mutation 1007fs is predictive of early onset of Crohn's disease with ileal stenosis, entero-enteral fistulas, and frequent need for surgical intervention with high risk of re-stenosis. Scand J Gastroenterol 41: 1421–1432 [PubMed]
34. Glas J, Seiderer J, Tillack C, Pfennig S, Beigel F, et al. (2010) The NOD2 single nucleotide polymorphisms rs2066843 and rs2076756 are novel and common Crohn's disease susceptibility gene variants. PLoS One 5: e14466. [PMC free article] [PubMed]
35. Glas J, Seiderer J, Wetzke M, Konrad A, Torok HP, et al. (2007) rs1004819 is the main disease-associated IL23R variant in German Crohn's disease patients: combined analysis of IL23R, CARD15, and OCTN1/2 variants. PLoS One 2: e819. [PMC free article] [PubMed]
36. Jurgens M, Laubender RP, Hartl F, Weidinger M, Seiderer J, et al. (2010) Disease activity, ANCA, and IL23R genotype status determine early response to infliximab in patients with ulcerative colitis. Am J Gastroenterol 105: 1811–1819 [PubMed]
37. Huett A, McCarroll SA, Daly MJ, Xavier RJ (2009) On the level: IRGM gene function is all about expression. Autophagy 5: 96–99 [PMC free article] [PubMed]
38. King KY, Lew JD, Ha NP, Lin JS, Ma X, et al. (2011) Polymorphic allele of human IRGM1 is associated with susceptibility to tuberculosis in African Americans. PLoS One 6: e16317. [PMC free article] [PubMed]
39. Jostins L, Ripke S, Weersma RK, Duerr RH, McGovern DP, et al. (2012) Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 491: 119–124 [PMC free article] [PubMed]
40. Latiano A, Palmieri O, Cucchiara S, Castro M, D'Inca R, et al. (2009) Polymorphism of the IRGM gene might predispose to fistulizing behavior in Crohn's disease. Am J Gastroenterol 104: 110–116 [PubMed]
41. Cooney R, Baker J, Brain O, Danis B, Pichulik T, et al. (2010) NOD2 stimulation induces autophagy in dendritic cells influencing bacterial handling and antigen presentation. Nat Med 16: 90–97 [PubMed]
42. Travassos LH, Carneiro LA, Ramjeet M, Hussey S, Kim YG, et al. (2010) Nod1 and Nod2 direct autophagy by recruiting ATG16L1 to the plasma membrane at the site of bacterial entry. Nat Immunol 11: 55–62 [PubMed]
43. Kaser A, Lee AH, Franke A, Glickman JN, Zeissig S, et al. (2008) XBP1 links ER stress to intestinal inflammation and confers genetic risk for human inflammatory bowel disease. Cell 134: 743–756 [PMC free article] [PubMed]
44. McGovern DP, Gardet A, Torkvist L, Goyette P, Essers J, et al. (2010) Genome-wide association identifies multiple ulcerative colitis susceptibility loci. Nat Genet 42: 332–337 [PMC free article] [PubMed]
45. Cadwell K, Liu JY, Brown SL, Miyoshi H, Loh J, et al. (2008) A key role for autophagy and the autophagy gene Atg16l1 in mouse and human intestinal Paneth cells. Nature 456: 259–263 [PMC free article] [PubMed]
46. Wehkamp J, Harder J, Weichenthal M, Schwab M, Schaffeler E, et al. (2004) NOD2 (CARD15) mutations in Crohn's disease are associated with diminished mucosal alpha-defensin expression. Gut 53: 1658–1664 [PMC free article] [PubMed]
47. Zhang F, Liu H, Chen S, Low H, Sun L, et al. (2011) Identification of two new loci at IL23R and RAB32 that influence susceptibility to leprosy. Nat Genet 43: 1247–1251 [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science