A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1.

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.

Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility.

A genome-wide association scan in individuals with Crohn's disease by the Wellcome Trust Case Control Consortium detected strong association at four novel loci. We tested 37 SNPs from these and other loci for association in an independent case-control sample. We obtained replication for the autophagy-inducing IRGM gene on chromosome 5q33.1 (replication P = 6.6 x 10(-4), combined P = 2.1 x 10(-10)) and for nine other loci, including NKX2-3, PTPN2 and gene deserts on chromosomes 1q and 5p13.

Independent and population-specific association of risk variants at the IRGM locus with Crohn's disease.

DNA polymorphisms in a region on chromosome 5q33.1 which contains two genes, immunity related GTPase related family, M (IRGM) and zinc finger protein 300 (ZNF300), are associated with Crohn's disease (CD). The deleted allele of a 20 kb copy number variation (CNV) upstream of IRGM was recently shown to be in strong linkage disequilibrium (LD) with the CD-associated single nucleotide polymorphisms and is itself associated with CD (P < 0.01). The deletion was correlated with increased or reduced expression of IRGM in transformed cells in a cell line-dependent manner, and has been proposed as a likely causal variant. We report here that small insertion/deletion polymorphisms in the promoter and 5' untranslated region of IRGM are, together with the CNV, strongly associated with CD (P = 1.37 x 10(-5) to 1.40 x 10(-9)), and that the CNV and the 5'-untranslated region variant -308(GTTT)(5) contribute independently to CD susceptibility (P = 2.6 x 10(-7) and P = 2 x 10(-5), respectively). We also show that the CD risk haplotype is associated with a significant decrease in IRGM expression (P < 10(-12)) in untransformed lymphocytes from CD patients. Further analysis of these variants in a Japanese CD case-control sample and of IRGM expression in HapMap populations revealed that neither the IRGM insertion/deletion polymorphisms nor the CNV was associated with CD or with altered IRGM expression in the Asian population. This suggests that the involvement of the IRGM risk haplotype in the pathogenesis of CD requires gene-gene or gene-environment interactions which are absent in Asian populations, or that none of the variants analysed are causal, and that the true causal variants arose after the European-Asian split.

Searching for genotype-phenotype structure: using hierarchical log-linear models in Crohn disease.

There has been considerable recent success in the detection of gene-disease associations. We consider here the development of tools that facilitate the more detailed characterization of the effect of a genetic variant on disease. We replace the simplistic classification of individuals according to a single binary disease indicator with classification according to a number of subphenotypes. This more accurately reflects the underlying biological complexity of the disease process, but it poses additional analytical difficulties. Notably, the subphenotypes that make up a particular disease are typically highly associated, and it becomes difficult to distinguish which genes might be causing which subphenotypes. Such problems arise in many complex diseases. Here, we concentrate on an application to Crohn disease (CD). We consider this problem as one of model selection based upon log-linear models, fitted in a Bayesian framework via reversible-jump Metropolis-Hastings approach. We evaluate the performance of our suggested approach with a simple simulation study and then apply the method to a real data example in CD, revealing a sparse disease structure. Most notably, the associated NOD2.908G-->R mutation appears to be directly related to more severe disease behaviors, whereas the other two associated NOD2 variants, 1007L-->FS and 702R-->W, are more generally related to disease in the small bowel (ileum and jejenum). The ATG16L1.300T-->A variant appears to be directly associated with only disease of the small bowel.

Investigation of Crohn's disease risk loci in ulcerative colitis further defines their molecular relationship.

BACKGROUND & AIMS: Identifying shared and disease-specific susceptibility loci for Crohn's disease (CD) and ulcerative colitis (UC) would help define the biologic relationship between the inflammatory bowel diseases. More than 30 CD susceptibility loci have been identified. These represent important candidate susceptibility loci for UC. Loci discovered by the index genome scans in CD have previously been tested for association with UC, but those identified in the recent meta-analysis await such investigation. Furthermore, the recently identified UC locus at ECM1 requires formal testing for association with CD. METHODS: We analyzed 45 single nucleotide polymorphisms, tagging 29 of the loci recently associated with CD in 2527 UC cases and 4070 population controls. We also genotyped the UC-associated ECM1 variant rs11205387 in 1560 CD patients and 3028 controls. RESULTS: Nine regions showed association with UC at a threshold corrected for the 29 loci tested (P < .0017). The strongest association (P = 4.13 x 10(-8); odds ratio = 1.27) was identified with a 170-kilobase region on chromosome 1q32 that contains 3 genes. We also found association with JAK2 and replicated a recently reported association with STAT3, further implicating the role of this signaling pathway in inflammatory bowel disease. Additional novel UC susceptibility genes were LYRM4 and CDKAL1. Twenty of the loci were not associated with UC, and several appear to be specific to CD. ECM1 variation was not associated with CD. CONCLUSIONS: Collectively, these data help define the genetic relationship between CD and UC and characterize common, as well as disease-specific mechanisms of pathogenesis.

Diverse effects of the CARD15 and IBD5 loci on clinical phenotype in 630 patients with Crohn's disease.

OBJECTIVES: Genetic variants at the CARD15 and IBD5 loci are strongly associated with Crohn's disease (CD), but evidence of the effect of these variants on the clinical expression of CD is conflicting and has often been hampered by small sample sizes. We studied 630 well-characterized patients to clarify the genotype/phenotype relationship in CD. METHODS: Patients and healthy controls were genotyped for three common mutations in CARD15 and a marker of the IBD5 risk haplotype. Allele frequencies were compared between phenotypic subgroups using chi2 or Fisher's exact tests. Genotype/phenotype analysis was carried out using multinomial logistic regression modelling allowing for adjustment for correlated or confounding factors. RESULTS: The mean age at diagnosis was significantly lower in carriers of the CARD15 or IBD5 risk alleles. After correction for age and smoking, CARD15 mutations were strongly associated with both ileal disease (P=8.8 x 10(-6)) and stenotic disease (P=0.003), but the association with stenotic disease appeared to be due to a confounding effect with ileal disease. CARD15 mutations were also associated with the presence of granulomas (P=5.7 x 10(-5)), which remained significant after adjustment for age at diagnosis and disease location (P=0.0047). The IBD5 risk haplotype frequency was significantly elevated in cases with perianal disease (P=0.028) and axial spondyloarthropathy (P=0.012). CONCLUSION: Genetic variants at the CARD15 and IBD5 loci have diverse effects on clinical expression in CD. CARD15 mutations are significantly correlated with the presence of granulomas.

Mutation, selection, and evolution of the Crohn disease susceptibility gene CARD15.

Three common mutations in the CARD15 (NOD2) gene are known to be associated with susceptibility to Crohn disease (CD), and genetic data suggest a gene dosage model with an increased risk of 2-4-fold in heterozygotes and 20-40-fold in homozygotes. However, the discovery of numerous rare variants of CARD15 indicates that some heterozygotes for the common mutations have a rare mutation on the other CARD15 allele, which would support a recessive model for CD. We addressed this issue by screening CARD15 for mutations in 100 CD patients who were heterozygous for one of the three common mutations. We also developed a strategy for evaluating potential disease susceptibility alleles (DSAs) that involves assessing the degree of evolutionary conservation of involved residues, predicted effects on protein structure and function, and genotyping in a large sample of cases and controls. The evolutionary analysis was aided by sequencing the entire coding region of CARD15 in three primates (chimp, gibbon, and tamarin) and aligning the human sequence with these and orthologs from other species. We found that 11 of the 100 CD patients screened had a second potential pathogenic mutation within the exonic and periexonic sequences examined. Assuming that there are no additional pathogenic mutations in noncoding regions, our study suggests that most carriers of the common DSAs are true heterozygotes, and supports previous evidence for a gene dosage model. Four novel nonsynonymous mutations were detected, one of which would produce premature termination of translation c.2686C>T (p.Arg896X). Two potential DSAs--c.2107C>T (p.Arg703Cys) and g.2238T>A (c.74-7T>A)--were significantly associated with CD in the case control sample. Analysis of the evolution of CARD15 revealed strong conservation of the encoded protein, with identity to the human sequence ranging from 99.1% in the chimp to 44.5% in fugu. Higher primates possess an open reading frame (ORF) upstream of the putative initiation site in other species that encodes a further 27 N-terminal amino acids, while four regions of high conservation are observed outside of the known domains of CARD15, indicative of additional residues of functional importance. The strategy developed here may have general application to the assessment of mutation pathogenicity and genetic models in other complex disorders.

Direct or indirect association in a complex disease: the role of SLC22A4 and SLC22A5 functional variants in Crohn disease.

A common haplotype spanning 250 kb on chromosome 5q31 is strongly associated with Crohn disease (CD). Recently, two functional variants within the SLC22A4 and SLC22A5 genes at this locus (IBD5), L503F (c.1507C > T) and G-207C (c.-207G > C), have been proposed to contribute directly to susceptibility to CD. However, extensive linkage disequilibrium at the IBD5 locus has complicated efforts to distinguish causal variants from association of the general risk haplotype. We genotyped the SLC22A4 and SLC22A5 variants and other polymorphisms across the risk haplotype in four populations of European origin, and applied regression-based haplotype analysis to over 1,200 fully genotyped case-control pairs, modeling case/control status on the presence of one or more SNPs to test for conditional association and to identify risk haplotypes. We found highly significant association of SNPs at the IBD5 locus with Crohn disease in all populations tested. However, the frequencies of L503F and G-207C in individuals who did not carry the general IBD5 risk haplotype were not significantly different in cases and controls, with associated disease odds ratios (ORs) of 0.90 (95% CI, 0.57-1.40) and 0.90 (95% CI, 0.65-1.23), respectively. Haplotype analysis showed that addition of the SLC22A4 and SLC22A5 variants to a null model that included the background risk haplotype did not significantly improve the model fit. In addition to the common risk haplotype, several rare haplotypes had an increased frequency in cases compared to controls. This study suggests that the molecular basis for Crohn disease susceptibility at the IBD5 locus remains to be defined, and highlights the challenge of the identification of causal variants in a complex disease in regions of extensive linkage disequilibrium.

Associations of allelic variants of the multidrug resistance gene (ABCB1 or MDR1) and inflammatory bowel disease and their effects on disease behavior: a case-control and meta-analysis study.

BACKGROUND: Allelic variants of the ATP-binding cassette, subfamily B member 1 (ABCB1), also known as the multidrug resistance gene (MDR1) that encodes the membrane-bound efflux transporter P-glycoprotein 170 (PGP-170), have been associated with inflammatory bowel disease but with conflicting results. METHODS: The present study examined the association of ABCB1 C3435T and G2677T/A in a large British case-control cohort of 828 Crohn's disease, 580 ulcerative colitis (UC) cases, and 285 healthy controls. The effect of these variants was further examined with respect to phenotypic and epidemiological characteristics. A meta-analysis was carried out of our results and those from 8 previously published association studies of the C3435T variant in inflammatory bowel disease. RESULTS: The 2677T allele was significantly increased in British UC cases compared with controls (45.2% vs. 39.6%; P = 0.034). In particular, the TT genotype was significantly associated with severe UC (odds ratio [OR] 1.90; 95% CI 1.01-3.55) and the use of steroids in UC (OR 1.77; 95% CI 1.08-2.88). No significant association was seen with C3435T and UC, Crohn's disease, or any clinical subgroup. A meta-analysis of 9 association studies of C3435T showed a significant association of the 3435T allele with UC (OR 1.12; 95% CI 1.02-1.23; P = 0.013) but not with CD. CONCLUSIONS: These results indicate that ABCB1 sequence variants are associated with a small increase in the risk of developing UC and may influence disease behavior.

Genome-wide association study for ulcerative colitis identifies risk loci at 7q22 and 22q13 (IL17REL).

We performed a genome-wide association analysis of 1,897,764 SNPs in 1,043 German ulcerative colitis (UC) cases and 1,703 controls. We discovered new associations at chromosome 7q22 (rs7809799) and at chromosome 22q13 in IL17REL (rs5771069) and confirmed these associations in six replication panels (2,539 UC cases and 5,428 controls) from different regions of Europe (overall study sample P(rs7809799) = 8.81 x 10(-11) and P(rs5771069) = 4.21 x 10(-8), respectively).

Mutational characterization of the bile acid receptor TGR5 in primary sclerosing cholangitis.

BACKGROUND: TGR5, the G protein-coupled bile acid receptor 1 (GPBAR1), has been linked to inflammatory pathways as well as bile homeostasis, and could therefore be involved in primary sclerosing cholangitis (PSC) a chronic inflammatory bile duct disease. We aimed to extensively investigate TGR5 sequence variation in PSC, as well as functionally characterize detected variants. METHODOLOGY/PRINCIPAL FINDINGS: Complete resequencing of TGR5 was performed in 267 PSC patients and 274 healthy controls. Six nonsynonymous mutations were identified in addition to 16 other novel single-nucleotide polymorphisms. To investigate the impact from the nonsynonymous variants on TGR5, we created a receptor model, and introduced mutated TGR5 constructs into human epithelial cell lines. By using confocal microscopy, flow cytometry and a cAMP-sensitive luciferase assay, five of the nonsynonymous mutations (W83R, V178M, A217P, S272G and Q296X) were found to reduce or abolish TGR5 function. Fine-mapping of the previously reported PSC and UC associated locus at chromosome 2q35 in large patient panels revealed an overall association between the TGR5 single-nucleotide polymorphism rs11554825 and PSC (odds ratio = 1.14, 95% confidence interval: 1.03-1.26, p = 0.010) and UC (odds ratio = 1.19, 95% confidence interval 1.11-1.27, p = 8.5 x 10(-7)), but strong linkage disequilibrium precluded demarcation of TGR5 from neighboring genes. CONCLUSIONS/SIGNIFICANCE: Resequencing of TGR5 along with functional investigations of novel variants provided unique insight into an important candidate gene for several inflammatory and metabolic conditions. While significant TGR5 associations were detected in both UC and PSC, further studies are needed to conclusively define the role of TGR5 variation in these diseases.

Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility.

Inflammatory bowel disease (IBD) typically manifests as either ulcerative colitis (UC) or Crohn's disease (CD). Systematic identification of susceptibility genes for IBD has thus far focused mainly on CD, and little is known about the genetic architecture of UC. Here we report a genome-wide association study with 440,794 SNPs genotyped in 1,167 individuals with UC and 777 healthy controls. Twenty of the most significantly associated SNPs were tested for replication in three independent European case-control panels comprising a total of 1,855 individuals with UC and 3,091 controls. Among the four consistently replicated markers, SNP rs3024505 immediately flanking the IL10 (interleukin 10) gene on chromosome 1q32.1 showed the most significant association in the combined verification samples (P = 1.35 x 10(-12); OR = 1.46 (1.31-1.62)). The other markers were located in ARPC2 and in the HLA-BTNL2 region. Association between rs3024505 and CD (1,848 cases, 1,804 controls) was weak (P = 0.013; OR = 1.17 (1.01-1.34)). IL10 is an immunosuppressive cytokine that has long been proposed to influence IBD pathophysiology. Our findings strongly suggest that defective IL10 function is central to the pathogenesis of the UC subtype of IBD.

Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease.

We report results of a nonsynonymous SNP scan for ulcerative colitis and identify a previously unknown susceptibility locus at ECM1. We also show that several risk loci are common to ulcerative colitis and Crohn's disease (IL23R, IL12B, HLA, NKX2-3 and MST1), whereas autophagy genes ATG16L1 and IRGM, along with NOD2 (also known as CARD15), are specific for Crohn's disease. These data provide the first detailed illustration of the genetic relationship between these common inflammatory bowel diseases.

Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease.

Several risk factors for Crohn's disease have been identified in recent genome-wide association studies. To advance gene discovery further, we combined data from three studies on Crohn's disease (a total of 3,230 cases and 4,829 controls) and carried out replication in 3,664 independent cases with a mixture of population-based and family-based controls. The results strongly confirm 11 previously reported loci and provide genome-wide significant evidence for 21 additional loci, including the regions containing STAT3, JAK2, ICOSLG, CDKAL1 and ITLN1. The expanded molecular understanding of the basis of this disease offers promise for informed therapeutic development.

Investigation of association of the DLG5 gene with phenotypes of inflammatory bowel disease in the British population.

BACKGROUND AND AIMS: Mutations in the DLG5 gene are associated with an increased risk of inflammatory bowel disease (IBD) in some European populations. Initial investigation of a British IBD population showed evidence of association of one of three DLG5 variants, R30Q, in a family-based collection but not in a case-control cohort. We have now examined the association of the R30Q polymorphism in a large cohort of British IBD cases, tested for interaction between the DLG5 and CARD15 genes and assessed possible association of DLG5 with clinical features of Crohn's disease (CD) and ulcerative colitis (UC). MATERIALS AND METHODS: DLG5 R30Q and the CARD15 polymorphisms, Arg702Trp, Gly908Arg and Leu1007fs were genotyped in 1,148 IBD cases and 749 controls. DLG5 R30Q was also analysed in cases stratified by CARD15 genotype, disease subtype and smoking history. RESULTS/FINDINGS: No significant difference in frequencies of the R30Q variant was observed between IBD cases (9.9%) and controls (10.1%) or in cases analysed separately as CD and UC. There was also no significant difference in the frequency of R30Q between CD cases carrying risk-associated CARD15 alleles and those that did not. The frequency of R30Q was higher in CD cases with ileal disease than cases without (p=0.042) and higher in CD cases who had smoked than in nonsmokers (p=0.009). INTERPRETATION/CONCLUSION: The R30Q variant in the DLG5 gene does not appear to be associated with an overall increase in the risk of disease in a British IBD cohort, but differences in its frequency in subgroups of CD patients warrant further investigation.

A nonsynonymous SNP in ATG16L1 predisposes to ileal Crohn's disease and is independent of CARD15 and IBD5.

BACKGROUND & AIMS: A genome-wide association scan of nonsynonymous DNA polymorphisms identified association of a threonine-to-alanine substitution (T300A) in the autophagy-related 16-like gene ATG16L1 with Crohn's disease. We investigated this association in independent U.K. cohorts of Crohn's disease and ulcerative colitis. METHODS: The T300A variant (rs2241880) was genotyped in an independent sample of 727 Crohn's disease and 877 ulcerative colitis cases, and in 579 controls. We then performed an extension analysis combining these data with the U.K. data from the initial study to give a total of 1236 U.K. Crohn's disease cases and 1235 controls to estimate disease risk and test for interaction with the CARD15 and IBD5 risk loci and for association with disease subtypes. RESULTS: The association of T300A was replicated in the independent sample of 727 Crohn's disease cases (P = .001), and was strongly associated in the extended analysis of 1236 Crohn's cases (P = 2.4 x 10(-6)). The 300A/A genotype conferred a 1.65-fold risk of Crohn's disease, with a 2.2-fold risk of ileal disease. Analysis of the interaction of ATG16L1 with CARD15 and IBD5 indicated that all 3 loci contribute independently to disease risk. Homozygosity for the risk allele at all 3 loci conferred a combined risk of 20.4 (95% confidence interval: 8.71, 47.7) for Crohn's disease. The ATG16L1 risk genotype showed a modest but significant association with ulcerative colitis (P = .026). CONCLUSIONS: The association of ATG16L1 with Crohn's disease and possibly with ulcerative colitis supports a role for autophagy in the pathogenesis of inflammatory bowel disease.

IL23R variation determines susceptibility but not disease phenotype in inflammatory bowel disease.

BACKGROUND & AIMS: Identification of inflammatory bowel disease (IBD) susceptibility genes is key to understanding pathogenic mechanisms. Recently, the North American IBD Genetics Consortium provided compelling evidence for an association between ileal Crohn's disease (CD) and the IL23R gene using genome-wide association scanning. External replication is a priority, both to confirm this finding in other populations and to validate this new technique. We tested for association between IL23R and IBD in a large independent UK panel to determine the size of the effect and explore subphenotype correlation and interaction with CARD15. METHODS: Eight single nucleotide polymorphism markers in IL23R tested in the North American study were genotyped in 1902 cases of Crohn's disease (CD), 975 cases of ulcerative colitis (UC), and 1345 controls using MassARRAY. Data were analyzed using chi(2) statistics, and subgroup association was sought. RESULTS: A highly significant association with CD was observed, with the strongest signal at coding variant Arg381Gln (allele frequency, 2.5% in CD vs 6.2% in controls [P = 1.1 x 10(-12)]; odds ratio, 0.38; 95% confidence interval, 0.29-0.50). A weaker effect was seen in UC (allele frequency, 4.6%; odds ratio, 0.73; 95% confidence interval, 0.55-0.96). Analysis accounting for Arg381Gln suggested that other loci within IL23R also influence IBD susceptibility. Within CD, there were no subphenotype associations or evidence of interaction with CARD15. CONCLUSIONS: This study shows an association between IL23R and all subphenotypes of CD with a smaller effect on UC. This extends the findings of the North American study, providing clear evidence that genome-wide association scanning can successfully identify true complex disease genes.

Systematic association mapping identifies NELL1 as a novel IBD disease gene.

Crohn disease (CD), a sub-entity of inflammatory bowel disease (IBD), is a complex polygenic disorder. Although recent studies have successfully identified CD-associated genetic variants, these susceptibility loci explain only a fraction of the heritability of the disease. Here, we report on a multi-stage genome-wide scan of 393 German CD cases and 399 controls. Among the 116,161 single-nucleotide polymorphisms tested, an association with the known CD susceptibility gene NOD2, the 5q31 haplotype, and the recently reported CD locus at 5p13.1 was confirmed. In addition, SNP rs1793004 in the gene encoding nel-like 1 precursor (NELL1, chromosome 11p15.1) showed a consistent disease-association in independent German population- and family-based samples (942 cases, 1082 controls, 375 trios). Subsequent fine mapping and replication in an independent sample of 454 French/Canadian CD trios supported the authenticity of the NELL1 association. Further confirmation in a large German ulcerative colitis (UC) sample indicated that NELL1 is a ubiquitous IBD susceptibility locus (combined p<10(-6); OR = 1.66, 95% CI: 1.30-2.11). The novel 5p13.1 locus was also replicated in the French/Canadian sample and in an independent UK CD patient panel (453 cases, 521 controls, combined p<10(-6) for SNP rs1992660). Several associations were replicated in at least one independent sample, point to an involvement of ITGB6 (upstream), GRM8 (downstream), OR5V1 (downstream), PPP3R2 (downstream), NM_152575 (upstream) and HNF4G (intron).

Genetic variation in myosin IXB is associated with ulcerative colitis.

BACKGROUND & AIMS: Common germline genetic variation in the 3' region of myosin IXB (MYO9B) has been associated recently with susceptibility to celiac disease, with a hypothesis that MYO9B variants might influence intestinal permeability. These findings suggested the current study investigating a possible further role for MYO9B variation in inflammatory bowel disease. METHODS: Eight single-nucleotide polymorphisms (SNPs) were selected to tag common haplotypes from the 35-kb 3' region of MYO9B. These included the strongest celiac disease-associated variants reported in a Dutch cohort. These SNPs were studied in 3 independently collected and genotyped case-control cohorts of European descent (UK, Dutch, and Canadian/Italian), comprising in total 2717 inflammatory bowel disease patients (1197 with Crohn's disease, 1520 with ulcerative colitis) and 4440 controls. RESULTS: Common variation in MYO9B was associated with susceptibility to inflammatory bowel disease in all 3 cohorts examined (most associated SNP, rs1545620; meta-analysis P = 1.9 x 10(-6); odds ratio, 1.2), with the same alleles showing association as reported for celiac disease. CONCLUSIONS: MYO9B genetic variants predispose to inflammatory bowel disease. Interestingly, rs1545620 is a nonsynonymous variant leading to an amino acid change (Ala1011Ser) in the third calmodulin binding IQ domain of MYO9B. Unlike previous variants (in other genes) reported to predispose to inflammatory bowel disease, the association at MYO9B was considerably stronger with ulcerative colitis, although weaker association with Crohn's disease also was observed. These data imply shared causal mechanisms underlying intestinal inflammatory diseases.