Background & aims

Genome-wide association studies have identified 200 inflammatory bowel disease (IBD) loci, but the genetic architecture of Crohn's disease (CD) and ulcerative colitis remain incompletely defined. Here, we aimed to identify novel associations between IBD and functional genetic variants using the Illumina ExomeChip (San Diego, CA).

Methods

Genotyping was performed in 10,523 IBD cases and 5726 non-IBD controls. There were 91,713 functional single-nucleotide polymorphism loci in coding regions analyzed. A novel identified association was replicated further in 2 independent cohorts. We further examined the association of the identified single-nucleotide polymorphism with microbiota from 338 mucosal lavage samples in the Mucosal Luminal Interface cohort measured using 16S sequencing.

Results

We identified an association between CD and a missense variant encoding alanine or threonine at position 391 in the zinc transporter solute carrier family 39, member 8 protein (SLC39A8 alanine 391 threonine, rs13107325) and replicated the association with CD in 2 replication cohorts (combined meta-analysis P = 5.55 × 10(-13)). This variant has been associated previously with distinct phenotypes including obesity, lipid levels, blood pressure, and schizophrenia. We subsequently determined that the CD risk allele was associated with altered colonic mucosal microbiome composition in both healthy controls (P = .009) and CD cases (P = .0009). Moreover, microbes depleted in healthy carriers strongly overlap with those reduced in CD patients (P = 9.24 × 10(-16)) and overweight individuals (P = 6.73 × 10(-16)).

Conclusions

Our results suggest that an SLC39A8-dependent shift in the gut microbiome could explain its pleiotropic effects on multiple complex diseases including CD.

Ulcerative colitis is a chronic inflammatory disease of the colon that presents as diarrhea and gastrointestinal bleeding. We performed a genome-wide association study using DNA samples from 1,052 individuals with ulcerative colitis and preexisting data from 2,571 controls, all of European ancestry. In an analysis that controlled for gender and population structure, ulcerative colitis loci attaining genome-wide significance and subsequent replication in two independent populations were identified on chromosomes 1p36 (rs6426833, combined P = 5.1 x 10(-13), combined odds ratio OR = 0.73) and 12q15 (rs1558744, combined P = 2.5 x 10(-12), combined OR = 1.35). In addition, combined genome-wide significant evidence for association was found in a region spanning BTNL2 to HLA-DQB1 on chromosome 6p21 (rs2395185, combined P = 1.0 x 10(-16), combined OR = 0.66) and at the IL23R locus on chromosome 1p31 (rs11209026, combined P = 1.3 x 10(-8), combined OR = 0.56; rs10889677, combined P = 1.3 x 10(-8), combined OR = 1.29).

Crohn's disease and ulcerative colitis, the two common forms of inflammatory bowel disease (IBD), affect over 2.5 million people of European ancestry, with rising prevalence in other populations. Genome-wide association studies and subsequent meta-analyses of these two diseases as separate phenotypes have implicated previously unsuspected mechanisms, such as autophagy, in their pathogenesis and showed that some IBD loci are shared with other inflammatory diseases. Here we expand on the knowledge of relevant pathways by undertaking a meta-analysis of Crohn's disease and ulcerative colitis genome-wide association scans, followed by extensive validation of significant findings, with a combined total of more than 75,000 cases and controls. We identify 71 new associations, for a total of 163 IBD loci, that meet genome-wide significance thresholds. Most loci contribute to both phenotypes, and both directional (consistently favouring one allele over the course of human history) and balancing (favouring the retention of both alleles within populations) selection effects are evident. Many IBD loci are also implicated in other immune-mediated disorders, most notably with ankylosing spondylitis and psoriasis. We also observe considerable overlap between susceptibility loci for IBD and mycobacterial infection. Gene co-expression network analysis emphasizes this relationship, with pathways shared between host responses to mycobacteria and those predisposing to IBD.

Primary sclerosing cholangitis (PSC) is a severe liver disease of unknown etiology leading to fibrotic destruction of the bile ducts and ultimately to the need for liver transplantation. We compared 3,789 PSC cases of European ancestry to 25,079 population controls across 130,422 SNPs genotyped using the Immunochip. We identified 12 genome-wide significant associations outside the human leukocyte antigen (HLA) complex, 9 of which were new, increasing the number of known PSC risk loci to 16. Despite comorbidity with inflammatory bowel disease (IBD) in 72% of the cases, 6 of the 12 loci showed significantly stronger association with PSC than with IBD, suggesting overlapping yet distinct genetic architectures for these two diseases. We incorporated association statistics from 7 diseases clinically occurring with PSC in the analysis and found suggestive evidence for 33 additional pleiotropic PSC risk loci. Together with network analyses, these findings add to the genetic risk map of PSC and expand on the relationship between PSC and other immune-mediated diseases.