Association of ATG16L1 gene haplotype with inflammatory bowel disease in Indians

Abstract

Inflammatory bowel disease (IBD) is characterized by multigenic inheritance. Defects in autophagy related genes are considered to show genetic heterogeneity between populations. We evaluated the association of several single nucleotide polymorphisms (SNPs) in the autophagy related 16 like 1 (ATG16L1) gene with IBD in Indians. The ATG16L1 gene was genotyped for ten different SNPs using DNA extracted from peripheral blood of 234 patients with Crohn’s disease (CD), 249 patients with ulcerative colitis (UC) and 393 healthy controls The SNPs rs2241880, rs4663396, rs3792106, rs10210302, rs3792109, rs2241877, rs6737398, rs11682898, rs4663402 and rs4663421 were genotyped using the Sequenom MassArray platform. PLINK was used for the association analysis and pairwise linkage disequilibrium (LD) values. Haplotype analysis was done using Haploview. All SNPs were in Hardy Weinberg equilibrium in cases and controls. The G allele at rs6737398 exhibited a protective association with both CD and UC. The T allele at rs4663402 and C allele at rs4663421 were positively associated with CD and UC. The T allele at rs2241877 exhibited protective association with UC only. The AA genotype at rs4663402 and the GG genotype at rs4663421 were protectively associated with both CD and UC. Haplotype analysis revealed that all the SNPs in tight LD (D’ = 0.76–1.0) and organized in a single haplotype block. Haplotype D was positively associated with IBD (P = 5.8 x 10⁻⁶ for CD and 0.002 for UC). SNPs in ATG16L1 were associated with IBD in Indian patients. The relevance to management of individual patients requires further study.

About the Authors:

Srinivasan Pugazhendhi

Current address: University of Kansas Medical Center, Kansas City, Kansas, United States of America

Affiliation: Wellcome Trust Research Laboratory, Christian Medical College, Vellore, India

Kirankumar Baskaran

Affiliation: Wellcome Trust Research Laboratory, Christian Medical College, Vellore, India

Srikanth Santhanam

Affiliation: Wellcome Trust Research Laboratory, Christian Medical College, Vellore, India

Balakrishnan S. Ramakrishna

* E-mail: [email protected]

Current address: SRM Institutes for Medical Science, Vadapalani, Chennai, India

Affiliation: Wellcome Trust Research Laboratory, Christian Medical College, Vellore, India

ORCID http://orcid.org/0000-0001-5090-9501

Introduction

Inflammatory bowel disease (IBD), comprising Crohn’s disease (CD) and ulcerative colitis (UC), is considered to result from abnormal immune reactions to gut luminal microbiota in genetically predisposed individuals [1–3]. The genetic predisposition is contributed by multiple genes through common variants (single nucleotide polymorphisms or SNPs) that each appear to exert a small influence. The discovery that CD was associated with mutations in the NOD2 gene [4,5] provided the first indication of the important role played by innate immunity in inflammation and disease progression. Subsequent genome wide association studies (GWAS) provided further insights into the genetic architecture of IBD, unravelling 163 susceptibility loci by the end of 2012 [3]. Disturbances in autophagy have been implicated as a potential pathogenetic pathway in IBD following the discovery that SNPs of genes involved in this pathway were associated with IBD [6].

Autophagy is a cellular process that involves the sequestration and eventual destruction of aging proteins, damaged cell organelles, apoptotic bodies and intracellular bacterial components by a specialized double membrane vesicle called the autophagosome. Autophagy plays a pivotal role in maintenance of immune homeostasis in the gut, contributing to both innate and adaptive immunity [7,8].

The link between autophagy and CD pathogenesis became first apparent when an association was identified between autophagy related 16-like 1 gene (ATG16L1) and CD [6]. ATG16L1 encodes a small coiled coil protein which interacts with ATG5 and ATG12 to form a 350 kDa multimeric complex that plays a crucial role in the bulk degradation or autophagy of cytoplasmic proteins and organelles. ATG16L1 protein is expressed in the colon, small bowel, intestinal epithelial cells, leukocytes and spleen. A coding SNP, named rs2241880, in the ATG16L1 was found to have a disease association with CD and this was responsible for threonine to alanine substitution (T300A) at amino acid 300 of protein. This SNP appeared to account for all of the disease risk exerted by the ATG16L1 locus. Further replication studies in an independent UK cohort confirmed the association existing between autophagy and IBD, particularly CD [9,10]. However, the literature shows both presence of and lack of association of IBD with SNPs in the ATG16L1 locus [11–16]. All of 9 ATG16L1 SNPs that were genotyped in a German population displayed significant protective association with CD, the strongest association being with rs2241879 and rs2241880 [17]. Other ATG16L1 gene variants independent of rs2241880 also appear to contribute to CD susceptibility [18]. Studies in several Asian countries including Japan, Korea and China failed to show an association between ATG16L1 gene variants and CD [19–21].

In light of the conflicting findings and ethnic differences in aforementioned studies we therefore aimed at analyzing various polymorphisms in ATG16L1 gene for its association with CD and UC in Indian population.

Materials and methods

876 participants were recruited from patients attending the outpatient and inpatient services of the Department of Gastrointestinal Sciences at the Christian Medical College. The cohort comprised of 234 CD patients, 249 UC patients and 393 healthy controls (HC). The diagnosis of CD and UC was based on a composite of clinical, radiological, endoscopic, and histopathological findings according to the consensus criteria of the Indian Society of Gastroenterology [22,23]. Patients with proven intestinal or extra-intestinal tuberculosis were excluded. Participants who refused the consent to participate were also excluded. Healthy adults accompanying non-IBD patients to the Gastroenterology clinic were recruited as controls. The history and clinical details were recorded and samples of venous blood were obtained in EDTA-coated Vacutainer tubes.

SNP genotyping

Genomic DNA was isolated from 8ml of EDTA–anticoagulated venous blood by salting out method. Isolated DNA was checked for quality and concentration, and stored at -80°C until analysis. Ten SNPs—rs2241880, rs4663396, rs3792106, rs10210302, rs3792109, rs2241877, rs6737398, rs11682898, rs4663402 and rs4663421 in ATG16L1 gene were selected for genotyping. Our choice of SNPs for genotyping was based on the following considerations. The first four have been previously found associated with CD in German [6,17] and UK [24] patients, while the remaining six were selected from the Hapmap data of Gujarati Indians [25]. Of these SNPs, rs2241880 is a coding SNP, rs10210302 is located in the 5’ UTR region, and the remainder are located in intronic regions. In German patients, rs2241880 and rs2241879 have both been shown to be associated with CD, with the latter being in close LD with the former. In other studies, rs2241880 was the most replicated SNP, and we therefore chose to study rs2241880 but not rs2241879.

Genotyping of ATG16L1 polymorphisms was performed using the Sequenom-MassArray platform at NxGenBio Life Sciences, New Delhi, India. As described earlier, in these assays a locus specific PCR reaction was carried out initially, followed by a locus specific primer extension (iPLEX) reaction in which an oligonucleotide primer annealed immediately upstream of the polymorphic site being genotyped. SNPs and small insertion/deletion polymorphisms were detected, by matrix-assisted laser desorption ionization–time-of-flight mass spectrometry after incubating primer and amplified DNA with mass-modified dideoxynucleotide terminators [26]. The primer extension was made according to the sequence of the variant site, and was a single complementary mass-modified base. Genotypes were assigned at each SNP locus using SpectroTYPER software (Sequenom).

Statistical analysis

To determine association of ATG16L1 SNPs with IBD susceptibility, comparison of allele and genotype distributions among cases and controls was done using PLINK v. 1.07 (website: http://zzz.bwh.harvard.edu/plink/) [27]. Odds ratios (ORs), 95% confidence intervals (CIs) and P-values were calculated. LD pairwise values, haplotype structure, and haplotype frequencies were determined using the Haploview software v. 4.2 [28]. Significance of difference between groups was analysed using Chi square test.

Ethical considerations

The Institutional Review Board of the Christian Medical College approved the study protocol and consent forms. Informed written consent for genetic testing was obtained from all study participants and from parents of minors participating in the study.

Results

Of the CD patients 152 were male and 82 female and their ages ranged from 10–78 years (median 44). Eighty eight had ileocolonic disease, 78 had ileal disease, 56 had colonic disease, 7 had upper GI involvement alone, and 5 had ileocolonic disease with involvement of upper GI tract. Forty three had stricturing disease, 22 had penetrating disease and the remaining had non-stricturing non-penetrating disease. Of the UC patients 145 were male and 104 female, and their ages ranged from 15 to 76 years (median 39). Thirty four had proctitis, 63 had left sided colitis and 152 had extensive colitis or pancolitis. The healthy controls, 256 male and 143 female, ranged in age from 15–73 years (median 30).

SNP analysis

All investigated polymorphisms were in Hardy Weinberg equilibrium in cases and controls. The allele and genotype frequencies of ATG16L1 polymorphisms are presented in Tables 1 and 2 respectively (see S1 Table for individual level data). Association analysis was carried out by comparing the allele and genotype frequencies of ATG16L1 markers amidst cases and controls.

[Figure omitted. See PDF.]

Table 1. Allele frequencies of ATG16L1 SNPs in cases and controls in Indian population.

https://doi.org/10.1371/journal.pone.0178291.t001

[Figure omitted. See PDF.]

Table 2. Genotype distributions of ATG16L1 SNPs in cases and controls.

https://doi.org/10.1371/journal.pone.0178291.t002

rs4663402 and rs4663421.

rs4663402 and rs4663421 exhibited unequivocal association with CD and UC in terms of both allele and genotype frequencies. The T allele of rs4663402 and the C allele of rs4663421 were associated with both CD (OR 2.65, 95% CI 1.62–4.35 and OR 3.17, 95% 1.91–5.25 respectively) and UC (OR 2.19, 95% CI 1.32–3.62 and OR 2.55, 95% CI 1.52–4.28 respectively) (Table 1). The genotype frequencies at rs4663402 and rs4663421 were also significantly associated with IBD (Table 2).

rs6737398.

The minor G allele at this locus was protectively associated with IBD, its frequency being higher in controls compared to cases. Corresponding OR (95% CI) values were 0.74 (0.58–0.94) for CD and 0.78 (0.62–0.99) for UC (Table 1). The GG genotype at this locus was noted in 18% of HC compared to 10.4% of CD patients (Table 2).

rs10210302.

Although the minor C allele of rs10210302 did not reveal significant association with IBD, the CC genotype of rs10210302 was protectively associated with IBD being present in 26% of HC compared to 18% of CD patients and 17% of UC patients (Table 2).

rs2241877.

The minor T allele at the rs2241877 locus demonstrated a protective association with UC, but not CD (Table 1). Correspondingly, the TT genotype was protectively associated with UC (Table 2).

The allele and genotype frequencies of rest of the ATG16L1 SNPs did not show a significant difference between patients and healthy individuals.

Haplotype analysis

Nine haplotypes (A-I) were constructed from the ten ATG16L1 SNPs (Table 3). Despite the existence of weak LD between SNPs rs3792106 and rs4663396 (D’ = 0.13), rs6737398 and rs4663396 (D’ = 0.31), rs6737398 and rs3792106 (D’ = 0.47) and rs11682898 & rs2241877 (D’ = 0.64), all the ten SNPs were organized in a single haplotype block (D’ ranging from 0.76–1.0) (Fig 1). On assessing the differences in haplotype frequencies amidst cases and controls it was found that haplotype D was positively associated with IBD being found in 8% of CD patients and 6% of UC patients compared to 2% of controls (P = 5.8x10^-6 and 0.002 respectively for CD and UC). Haplotype B exhibited a protective association with UC being present in 16% of UC patients and 22% of HC (Table 4).

[Figure omitted. See PDF.]

Fig 1. Linkage disequilibrium (LD) pattern between the ten ATG16L1 single nucleotide polymorphisms (SNPs) genotyped in this study.

All SNPs were in tight LD (D’ = 0.76–1.0) and were organized in a single haplotype block.

https://doi.org/10.1371/journal.pone.0178291.g001

[Figure omitted. See PDF.]

Table 3. Haplotype structure of ten ATG16L1 SNPs.

https://doi.org/10.1371/journal.pone.0178291.t003

[Figure omitted. See PDF.]

Table 4. Haplotype frequency in healthy controls and IBD patients.

https://doi.org/10.1371/journal.pone.0178291.t004

Discussion

In the present study we explored the association of IBD with a number of genetic variants in the ATG16L1 gene. The SNPs rs4663421 and rs4663402 revealed a positive association with IBD. The SNPs rs6737398 and rs10210302 had significant protective associations, at the genotype level, with CD and UC. The protective association of rs2241877 was noted with UC alone. While specific alleles were associated, it became clear that a single haplotype was responsible for the strong association with both CD and UC.

Following the identification of genetic alterations within ATG16L1 to be associated with CD and UC development by two genome wide association studies [6,24], a significant number of studies have replicated this association in populations of Caucasian origin [11,17], corroborating the relevance of ATG16L1 to IBD genesis. Earlier data suggested that all the CD risk associated with ATG16L1 could be fully attributed to the T300A variant [6]. Hence most studies confined themselves to analyzing the association of this marker with IBD. Notable exceptions were studies in the German [17] and Japanese [19] populations in which 9 and 12 SNPs respectively were typed. Strong association of T300A (rs2241880) with CD has been established in populations in New Zealand [11] [P = 0.0001, OR = 1.41, 95% CI = 1.18–1.67], Australia [16] [P< 0.001, OR = 1.96, 95% CI = 1.49–2.58], Netherlands [12] [OR = 1.36, 95% CI = 1.21–1.66, P = 0.0017], and Hungary [29] [P = 0.004, OR = 1.69, 95% CI = 1.19–2.41]. On the other hand, a lack of association between rs2241880 and CD has been reported from Italy [30], Lithuania [31] and Brazil [32]. Protective and positive associations of the T300A variant with UC have been reported in Australian [16] and UK [9] studies respectively. However, several other studies failed to find a significant association of T300A with UC [6,10]. In the present study rs2241880 did not have any significant association with either CD or UC. This finding is consistent with a study reported in the abstract in which 216 Indian Asian IBD patients (60 CD: 156 UC) were compared with 203 healthy Indian Asian controls and no association was found between rs2241880 and IBD [33].

ATG16L1 is a member of a large family of genes involved in autophagy, a housekeeping cellular process by which degraded intracellular components, and excess and damaged organelles are digested and removed [34]. Autophagy is believed to have a role in immune defense against intracellular pathogens [35], which may explain why defects in this pathway are important in IBD pathogenesis. For instance, functional knockdown of the ATG16L1 gene abrogates the autophagy of Salmonella Typhimurium. The T300A variant at rs2241880 in ATG16L1 gene has been considered to be the key variant associated with IBD [6]. In a human intestinal epithelial cell line, defective antibacterial autophagy was noted in cells expressing the CD associated variant protein (Ala) compared to cells expressing the wild type protein (Thr) [36], demonstrating that the non-synonymous SNP is associated with altered responses to bacterial components. In a study of ileal tissue associated microbiota from patients with ileal CD, patients homozygous for the ATG16L1 risk allele had abundant Fusobacteriaceae, whereas patients homozygous for the ATG16L1 protective allele showed reduced abundance of Bacteroidaceae and Enterobacteriaceae [37]. Monocytes from these patients showed impaired killing of adherent invasive Escherichia coli in patients homozygous for the ATG16L1 risk allele compared with those homozygous for the protective allele [37].

The SNPs rs3792106 and rs4663396 were associated protectively with CD (OR 0.77, 95% CI 0.67–0.87 and OR 0.81, 95% CI 0.68–0.95 respectively) in a German population [17]. On the other hand, these SNPs did not show any association with CD in Japanese [19]. In the present study, no significant association was detected between these SNPs and CD.

The SNP rs10210302 was very significantly associated with CD in the Wellcome Trust Case Control Consortium study [34]. However, a study in UC patients failed to identify any significant association between rs10210302 and UC [38]. In the present study the homozygous CC genotype of rs10210302 was significantly protectively associated with both CD and UC. It should be noted that the rs10210302 polymorphism was not associated with CD in an Ashkenazi Jewish cohort [39]. A meta-analysis (including 15 studies with 9211 UC and 10899 controls) of the association between UC and rs2241880 indicated a significant association between the G allele and increased risk for UC (odds ratio 1.08, p = 0.0003) [40]. The above comparisons indicate that there is significant variation with regard to the frequencies of ATG16L1 variants, and their association with IBD, in ethnically divergent populations. A recent meta-analysis concluded that ATG16L1 was not associated with IBD in East Asians [41]. However an analysis of IBD genetics data from 86,640 European individuals and 9,846 individuals of East Asian, Indian or Iranian descent concluded that heterogeneity between populations in IBD genetics was driven, in the case of ATG16L1, by effect size rather than by allele frequency [42].

Of interest, the present study also revealed significant associations between SNPs rs4663421, rs4663402, rs6737398 and rs2241877 in the ATG16L1 gene and IBD. There is as yet no data on the association of these variants with IBD in other populations.

It is also to be noted that, with the exception of a single Japanese study [19], haplotype analysis has not been done with the SNPs of ATG16L1 in any of the association studies in various populations. Of the nine haplotypes (A-I) that we identified in our population, only one was majorly associated with both forms of IBD.

The functional consequence of these mutations can only be indirectly inferred as there are no studies which directly examine the influence of these SNPs on ATG16L1 expression. In the present study the coding SNP rs2241880 was not associated with IBD in the population studied. However SNPs rs4663421 and rs4663402 were significantly associated with IBD and were included in the haplotype that was associated with IBD. Using HaploReg (http://archive.broadinstitute.org/mammals/haploreg/haploreg_v3.php) we identified that these two SNPs are in close linkage disequilibrium (0.97 for rs4663421 and 0.91 for rs4663402) with rs140725980 (S2 Table). The latter SNP (rs140725980) is an enhancer histone mark for ATG16L1 gene in gastrointestinal tissue. Enhancer histone marks flank sites of transcription factor binding and enhancer activity and histone marked active enhancers are associated with level of gene expression. Methylation of histone H3K4 is tightly associated with the promoters of active genes. Thus it is possible that the haplotype that we identified is associated with altered expression of ATG16L1 in gastrointestinal tissues, explaining a role for the haplotype in the genesis of IBD.

This is the first study showing significant association of ATG16L1 variants rs4663421, rs4663402, rs6737398 and rs2241877 with CD and UC in Indian population, identifying autophagy as a host pathogen defense pathway crucial to IBD biology. Further elucidation of role of this process in IBD pathogenesis will be facilitated by detailed investigation into the functional consequences of these variants, besides exploring other proteins involved in autophagy pathway.

Supporting information

[Figure omitted. See PDF.]

S1 Table. This table shows the individual level data on each allele in each participant.

Following are the expansions and codes: PID is the participant ID. Sex, 1 = male, 2 = female. Diagnosis, 1 = control, 2 = Crohn’s disease, 3 = ulcerative colitis. For individual SNPs each allele is shown separately, 1 = A, 2 = C, 3 = G, 4 = T, and 0 = untypable.

https://doi.org/10.1371/journal.pone.0178291.s001

(XLSX)

S2 Table. This table shows the HaploReg output when queried for the associations of rs4663402 and rs4663421.

https://doi.org/10.1371/journal.pone.0178291.s002

(XLSX)

Author Contributions

1. Conceptualization: SP BSR.

2. Data curation: SP.

3. Formal analysis: SP.

4. Funding acquisition: BSR.

5. Investigation: SS KB BSR.

6. Methodology: SP BSR.

7. Project administration: BSR.

8. Resources: BSR.

9. Supervision: SP BSR.

10. Writing – original draft: KB.

11. Writing – review & editing: BSR.

Citation: Pugazhendhi S, Baskaran K, Santhanam S, Ramakrishna BS (2017) Association of ATG16L1 gene haplotype with inflammatory bowel disease in Indians. PLoS ONE 12(5): e0178291. https://doi.org/10.1371/journal.pone.0178291

References

1. McGovern DP, Kugathasan S, Cho JH. Genetics of inflammatory bowel diseases. Gastroenterology. 2015; 149: 1163–1176. pmid:26255561

2. de Souza HS, Fiocchi C. Immunopathogenesis of IBD: current state of the art. Nat Rev Gastroenterol Hepatol. 2016; 13: 13–27. pmid:26627550

3. Jostins L, Ripke S, Weersma RK, Duerr RH, McGovern DP, Hui KY et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature. 2012; 491: 119–24. pmid:23128233

4. Ogura Y, Bonen DK, Inohara N, Nicolae DL, Chen FF, Ramos R et al. A frameshift mutation in NOD2 associated with susceptibility to Crohn’s disease. Nature. 2001; 411: 603–6. pmid:11385577

5. Hampe J, Cuthbert A, Croucher PJ, Mirza MM, Mascheretti S, Fisher S et al. Association between insertion mutation in NOD2 gene and Crohn’s disease in German and British populations. Lancet. 2001; 357: 1925–8. pmid:11425413

6. Hampe J, Franke A, Rosenstiel P, Till A, Teuber M, Huse K, et al. A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1. Nat. Genet. 2007; 39: 207–211. pmid:17200669

7. Cooney R, Baker J, Brain O, Danis B, Pichulik T, Allan P et al. NOD2 stimulation induces autophagy in dendritic cells influencing bacterial handling and antigen presentation. Nat. Med. 2010; 16: 90–97. pmid:19966812

8. Munz C. Antigen processing for MHC class II presentation via autophagy. Front. Immunol. 2012; 3: 9. pmid:22566895

9. Prescott NJ, Fisher SA, Franke A, Hampe J, Onnie CM, Soars D et al. A nonsynonymous SNP in ATG16L1 predisposes to ileal Crohn’s disease and is independent of CARD15 and IBD5. Gastroenterology. 2007; 132: 1665–71. pmid:17484864

10. Cummings JR, Cooney R, Pathan S, Anderson CA, Barrett JC, Beckly J et al. Confirmation of the role of ATG16L1 as a Crohn’s disease susceptibility gene. Inflamm Bowel Dis. 2007; 13: 941–6. pmid:17455206

11. Roberts RL, Gearry RB, Hollis-Moffatt JE, Miller AL, Reid J, Abkevich V 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: 2754–2761. pmid:17894849

12. Weersma RK, Zhernakova A, Nolte IM, Lefebvre C, Rioux JD, Mulder F et al. ATG16L1 and IL23R are associated with inflammatory bowel disease but not with celiac disease in the Netherlands. Am J Gastroenterol. 2008; 103: 621–627. pmid:18047540

13. Okazaki T, Wang M-H, Rawsthorne P, Sargent M, Datta LW, Shugart YY et al. Contributions of IBD5, IL23R, ATG16L1 and NOD2 to Crohn’s disease risk in a population based case-control study: evidence of gene-gene interaction. Inflamm Bowel Dis. 2008; 14: 1528–41. pmid:18521914

14. Lappalainen M, Halme L, Turunen U, Saavalainen P, Einarsdottir E, Färkkilä M et al. Association of IL23R, TNFRSF1A, and HLA-DRB1*0103 allele variants with inflammatory bowel disease phenotypes in the Finnish population. Inflamm Bowel Dis. 2008; 14: 1118–1124. pmid:18338763

15. Lakatos PL, Szamosi T, Szilvasi A, Molnar E, Lakatos L, Kovacs A, et al. ATG16L1 and IL23 receptor (IL23R) genes are associated with disease susceptibility in Hungarian CD patients. Dig Liv Dis. 2008; 40: 867–873.

16. Fowler EV, Doecke J, Simms LA, Zhao ZZ, Webb PM, Hayward NK et al. ATG16L1 T300A shows strong associations with disease subgroups in a large Australian IBD population: further support for significant disease heterogeneity. Am J Gastroenterol. 2008; 103: 1–8.

17. Glas J, Konrad A, Schmechel S, Dambacher J, Seiderer J, Schroff F et al. The ATG16L1 gene variants rs2241879 and rs2241880 (T300A) are strongly associated with susceptibility to Crohn’s disease in the German population. Am J Gastroenterol. 2008; 103: 682–691. pmid:18162085

18. Van Limbergen J, Kabakchiev B, Stempak J, Schumm P, Xu W, Henderson P, et al. Hypothesis-free analysis of ATG16L1 demonstrates gene-wide extent of association with Crohn’s disease susceptibility. Gut. 2013; 62: 331–333. pmid:22833393

19. Yamazaki K, Onouchi Y, Takazoe M, Kubo M, Nakamura Y, Hata A. Association analysis of genetic variants in IL23R, ATG16L1 and 5p13.1 loci with Crohn’s disease in Japanese patients. J Hum Genet. 2007; 52: 575–583. pmid:17534574

20. Yang SK, Park M, Lim J, Park SH, Ye BD, Lee I et al. Contribution of IL23R but not ATG16L1 to Crohn’s disease susceptibility in Koreans. Inflamm Bowel Dis. 2009; 15: 1385–1390. pmid:19334001

21. Zhi J, Zhi FC, Chen ZY, Yao GP, Guan J, Lin Y et al. [Correlation of the autophagosome gene ATG16L1 polymorphism and inflammatory bowel disease.] Nan Fang Yi Ke Da Xue Xue Bao.2008;28:649–651 (in Chinese). pmid:18495612

22. Ramakrishna BS, Makharia GK, Ahuja V, Ghoshal UC, Jayanthi V, Perakath B et al. Indian Society of Gastroenterology consensus statements on Crohn's disease in India. Indian J Gastroenterol. 2015; 34: 3–22. pmid:25772856

23. Ramakrishna BS, Makharia GK, Abraham P, Ghoshal UC, Jayanthi V, Agarwal BK et al. Indian Society of Gastroenterology consensus on ulcerative colitis. Indian J Gastroenterol. 2012; 31: 307–23. pmid:23096266

24. 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: 661–78. pmid:17554300

25. http://www.ncbi.nlm.nih.gov/projects/SNP/snp_viewTable.cgi?pop=12159. Accessed on 03.05.2016.

26. Baskaran K, Pugazhendhi S, Ramakrishna BS. Protective association of tumor necrosis factor superfamily 15 (TNFSF15) polymorphic haplotype with ulcerative colitis and Crohn's disease in an Indian population. PLoS One. 2014; 9: e114665. pmid:25501099

27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D et al. PLINK: a toolset for whole-genome association and population-based linkage analysis. American Journal of Human Genetics. 2007; 81: 559–75. pmid:17701901

28. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21: 263–265. pmid:15297300

29. Csongei V, Jaromi L, Safrany E, Sipeky C, Magyari L, Farago B et al. Interaction of the major inflammatory bowel disease susceptibility alleles in Crohn’s disease patients. World J Gastroenterol. 2010; 16: 176–183. pmid:20066736

30. Perricone C, Borgiani P, Romano S, Ciccacci C, Fusco G, Novelli G et al. ATG16L1 Ala197Thr is not associated with susceptibility to Crohn’s disease or with phenotype in an Italian population. Gastroenterology. 2008; 134: 368–370. pmid:18166373

31. Sventoraityte J, Zvirbliene A, Franke A, Kwiatkowski R, Kiudelis G, Kupcinskas L et al. NOD2, IL23R and ATG16L1 polymorphisms in Lithuanian patients with inflammatory bowel disease. World J Gastroenterol. 2010; 16: 359–364. pmid:20082483

32. Baptista ML, Amarante H, Picheth G, Sdepanian VL, Peterson N, Babasukumar U et al. CARD15 and IL23R influences Crohn’s disease susceptibility but not disease phenotype in a Brazilian population. Inflamm Bowel Dis. 2008; 14: 674–679. pmid:18200510

33. Walker DG, Bancil AS, Rai PS, Williams HR, Chambers J, Kooner JS et al. Ethnic variation in the frequency of IBD related polymorphisms in IRGM, ATG16L1 and IL23R. Gut 2011; 60: A223.

34. Deretic V. Autophagy in innate and adaptive immunity. Trends Immunol. 2005; 26: 523–8. pmid:16099218

35. Schmid D, Dengiel J, Schoor O, Stevanovic S, Munz C. Autophagy in innate and adaptive immunity against intracellular pathogens. J Mol Med. 2006; 84: 194–202. pmid:16501849

36. Kuballa P, Huett A, Rioux JD, Daly MJ, Xavier RJ. Impaired autophagy of an intracellular pathogen induced by a Crohn's disease associated ATG16L1 variant. PLoS One 3; 2008: e33913227.

37. Sadaghian Sadabad M, Regeling A, de Goffau MC, Blokzijl T, Weersma RK, Penders J et al. The ATG16L1-T300A allele impairs clearance of pathosymbionts in the inflamed ileal mucosa of Crohn's disease patients. Gut. 2015; 64: 1546–52. pmid:25253126

38. Fisher SA, Tremelling M, Anderson CA, Gwilliam R, Bumpstead S, Prescott NJ et al. Genetic determinants of ulcerative colitis includes the ECM1 locus and five loci implicated in Crohn’s disease. Nat Genet. 2008; 40: 710–712. pmid:18438406

39. Peter I, Mitchell AA, Ozelius L, Erazo M, Hu J et al. Evaluation of 22 genetic variants with Crohn's disease risk in the Ashkenazi Jewish population: a case-control study. BMC Med Genet. 2011; 12: 63. pmid:21548950

40. Palomino-Morales RJ, Oliver J, Gómez-García M, López-Nevot MA, Rodrigo L, Nieto A et al. Association of ATG16L1 and IRGM genes polymorphisms with inflammatory bowel disease: a meta-analysis approach. Genes Immun. 2009; 10: 356–64. pmid:19491842

41. Zhang J, Chen J, Gu J, Guo H, Chen W. Association of IL23R and ATG16L1 with susceptibility of Crohn's disease in Chinese population. Scand J Gastroenterol. 2014; 49: 1201–6. pmid:25048429

42. Liu JZ, van Sommeren S, Huang H, Ng SC, Alberts R, Takahashi A et al. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat Genet. 2015; 47: 979–86. pmid:26192919