The genetic basis of antineutrophil cytoplasmic antibody, an important biomarker of inflammatory bowel disease (IBD), has never been thoroughly examined on a genome-wide scale. In this study, we performed a 2-stage genome-wide association study (GWAS) on antineutrophil cytoplasmic antibody in IBD cases. In the 2959 IBD cases in the discovery stage, we observed an association between a variant in the gene TNFRSF1B with antineutrophil cytoplasmic antibody level (rs5745994, minor allele frequency = 0.028, beta = 18.12, 95% CI, 11.82-24.22, P = 1.89 × 10). This association was replicated in an independent cohort of 419 IBD cases (beta = 16.91, 95% CI, 6.13-27.69, P = 2.38 × 10). With a Q-value of 0.036, we performed a fixed-effect meta-analysis for the association of rs5745994 in both cohorts and observed a stronger association signal (beta = 17.81, 95% CI, 12.36-23.25, P = 8.97 × 10). TNFRSF1B gene codes for tumor necrosis factor (TNF) receptor 2 (TNFR2), thereby we examined the reported TNFRSF1B variant with serum TNFR2 level. We observed a negative association with serum TNFR2 level being 8.23 EU/mL in carriers and 9.12 EU/mL in noncarriers (P = 0.033). This finding indicates the functional role of identified TNFRSF1B variant in IBD serology and may be reflective of the underlying biological mechanisms that determine clinical expression and/or response to certain therapies.

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.