Comprehensive Association Analyses of Extraintestinal Manifestations in Inflammatory Bowel Disease

BACKGROUND & AIMS

Patients with inflammatory bowel disease (IBD) frequently develop extraintestinal manifestations (EIMs) that contribute substantially to morbidity. We assembled the largest multicohort data set to date to investigate the clinical, serologic, and genetic factors associated with EIM complications in IBD.

METHODS

Data were available in 12,083 unrelated European ancestry IBD cases with presence or absence of EIMs (eg, ankylosing spondylitis [ankylosing spondylitis and sacroiliitis], primary sclerosing cholangitis [PSC], peripheral arthritis, and skin and ocular manifestations) across 4 cohorts (Cedars-Sinai Medical Center, National Institute for Diabetes and Digestive and Kidney Diseases IBD Genetics Consortium, Sinai Helmsley Alliance for Research Excellence Consortium, and Risk Stratification and Identification of Immunogenetic and Microbial Markers of Rapid Disease Progression in Children with Crohn's Disease cohort). Clinical and serologic parameters were analyzed by means of univariable and multivariable regression analyses using a mixed-effects model. Within-case logistic regression was performed to assess genetic associations.

RESULTS

Most EIMs occurred more commonly in female subjects (overall EIM: P = 9.0E-05, odds ratio [OR], 1.2; 95% CI, 1.1-1.4), with CD (especially colonic disease location; P = 9.8E-09, OR, 1.7; 95% CI, 1.4-2.0), and in subjects who required surgery (both CD and UC; P = 3.6E-19, OR, 1.7; 95% CI, 1.5-1.9). Smoking increased risk of EIMs except for PSC, where there was a "protective" effect. Multiple serologic associations were observed, including with PSC (IgG and IgA, perinuclear anti-nuclear cytoplasmic antibody; anti-Saccharomyces cerevisiae antibodies; and anti-flagellin) and any EIM (IgG and IgA, perinuclear anti-nuclear cytoplasmic antibody; anti-Saccharomyces cerevisiae antibodies; and anti-Pseudomonas fluorescens-associated sequence). We identified genome-wide significant associations within major histocompatibility complex (ankylosing spondylitis and sacroiliitis, P = 1.4E-15; OR, 2.5; 95% CI, 2.0-3.1; PSC, P = 2.7E-10; OR, 2.8; 95% CI, 2.0-3.8; ocular, P = 2E-08, OR, 3.6; 95% CI, 2.3-5.6; and overall EIM, P = 8.4E-09; OR, 2.2; 95% CI, 1.7-2.9) and CPEB4 (skin, P = 2.7E-08; OR, 1.5; 95% CI, 1.3-1.8). Genetic associations implicated tumor necrosis factor, JAK-STAT, and IL6 as potential targets for EIMs. Contrary to previous reports, only 2% of our subjects had multiple EIMs and most co-occurrences were negatively correlated.

CONCLUSIONS

We have identified demographic, clinical, and genetic associations with EIMs that revealed underlying mechanisms and implicated novel and existing drug targets-important steps toward a more personalized approach to IBD management.

Clinical predictors of early and late endoscopic recurrence following ileocolonic resection in Crohn's disease

BACKGROUND AND AIM

Multiple factors are suggested to place Crohn's disease patients at risk of recurrence after ileocolic resection with conflicting associations. We aimed to identify clinical predictors of recurrence at first (early) and further (late) postoperative colonoscopy.

METHODS

Crohn's disease patients undergoing ileocolic resection were prospectively recruited at six North American centers. Clinical data was collected and endoscopic recurrence was defined as Rutgeerts score ≥i2. A multivariable model was fitted to analyze variables independently associated with recurrence.

RESULTS

A total of 365 patients undergoing 674 postoperative colonoscopies were included with a median age of 32 years, 189 (51.8%) were male and 37 (10.1%) non-Whites. Postoperatively, 133 (36.4%) used anti-TNF and 30 (8.2%) were smokers. At first colonoscopy, 109 (29.9%) had recurrence. Male gender (OR = 1.95, 95% CI 1.12 - 3.40), non-White ethnicity (OR = 2.48, 95% CI 1.09 - 5.63), longer interval between surgery and colonoscopy (OR = 1.09, 95% CI 1.002 - 1.18), and postoperative smoking (OR = 2.78, 95% CI 1.16 - 6.67) were associated with recurrence, while prophylactic anti-TNF reduced the risk (OR = 0.28, 95% CI 0.14 - 0.55). Postoperative anti-TNF prophylaxis had a protective effect on anti-TNF experienced patients but not on anti-TNF naïve patients. Among patients without recurrence at first colonoscopy, Rutgeerts score i1 was associated with subsequent recurrence (OR = 4.43, 95% CI 1.73 - 11.35).

CONCLUSIONS

We identified independent clinical predictors of early and late Crohn's disease postoperative endoscopic recurrence. Clinical factors traditionally used for risk stratification failed to predict recurrence and need to be revised.

Genetic coding variant in complement factor B (CFB) is associated with increased risk for perianal Crohn's disease and leads to impaired CFB cleavage and phagocytosis

OBJECTIVE

Perianal Crohn's disease (pCD) occurs in up to 40% of patients with CD and is associated with poor quality of life, limited treatment responses and poorly understood aetiology. We performed a genetic association study comparing CD subjects with and without perianal disease and subsequently performed functional follow-up studies for a pCD associated SNP in Complement Factor B (CFB).

DESIGN

Immunochip-based meta-analysis on 4056 pCD and 11 088 patients with CD from three independent cohorts was performed. Serological and clinical variables were analysed by regression analyses. Risk allele of rs4151651 was introduced into human CFB plasmid by site-directed mutagenesis. Binding of recombinant G252 or S252 CFB to C3b and its cleavage was determined in cell-free assays. Macrophage phagocytosis in presence of recombinant CFB or serum from CFB risk, or protective CD or healthy subjects was assessed by flow cytometry.

RESULTS

Perianal complications were associated with colonic involvement, OmpC and ASCA serology, and serology quartile sum score. We identified a genetic association for pCD (rs4151651), a non-synonymous SNP (G252S) in CFB, in all three cohorts. Recombinant S252 CFB had reduced binding to C3b, its cleavage was impaired, and complement-driven phagocytosis and cytokine secretion were reduced compared with G252 CFB. Serine 252 generates a de novo glycosylation site in CFB. Serum from homozygous risk patients displayed significantly decreased macrophage phagocytosis compared with non-risk serum.

CONCLUSION

pCD-associated rs4151651 in CFB is a loss-of-function mutation that impairs its cleavage, activation of alternative complement pathway, and pathogen phagocytosis thus implicating the alternative complement pathway and CFB in pCD aetiology.

Serum Lipidomic Screen Identifies Key Metabolites, Pathways, and Disease Classifiers in Crohn's Disease

BACKGROUND

There is an unmet medical need for biomarkers that capture host and environmental contributions in inflammatory bowel diseases (IBDs). This study aimed at testing the potential of circulating lipids as disease classifiers given their major roles in inflammation.

METHODS

We applied a previously validated comprehensive high-resolution liquid chromatography-mass spectrometry-based untargeted lipidomic workflow covering 25 lipid subclasses to serum samples from 100 Crohn's disease (CD) patients and 100 matched control subjects. Findings were replicated and expanded in another 200 CD patients and 200 control subjects. Key metabolites were tested for associations with disease behavior and location, and classification models were built and validated. Their association with disease activity was tested using an independent cohort of 42 CD patients.

RESULTS

We identified >70 metabolites with strong association (P < 1 × 10-4, q < 5 × 10-4) to CD. Highly performing classification models (area under the curve > 0.84-0.97) could be built with as few as 5 to 9 different metabolites, representing 6 major correlated lipid clusters. These classifiers included a phosphatidylethanolamine ether (O-16:0/20:4), a sphingomyelin (d18:1/21:0) and a cholesterol ester (14:1), a very long-chain dicarboxylic acid [28:1(OH)] and sitosterol sulfate. These classifiers and correlated lipids indicate a dysregulated metabolism in host cells, notably in peroxisomes, as well as dysbiosis, oxidative stress, compromised inflammation resolution, or intestinal membrane integrity. A subset of these were associated with disease behavior or location.

CONCLUSIONS

Untargeted lipidomic analyses uncovered perturbations in the circulating human CD lipidome, likely resulting from multiple pathogenic mechanisms. Models using as few as 5 biomarkers had strong disease classifier characteristics, supporting their potential use in diagnosis or prognosis.

Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

SECANT: a biology-guided semi-supervised method for clustering, classification, and annotation of single-cell multi-omics

The recent advance of single cell sequencing (scRNA-seq) technology such as Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq) allows researchers to quantify cell surface protein abundance and RNA expression simultaneously at single cell resolution. Although CITE-seq and other similar technologies have gained enormous popularity, novel methods for analyzing this type of single cell multi-omics data are in urgent need. A limited number of available tools utilize data-driven approach, which may undermine the biological importance of surface protein data. In this study, we developed SECANT, a biology-guided SEmi-supervised method for Clustering, classification, and ANnoTation of single-cell multi-omics. SECANT is used to analyze CITE-seq data, or jointly analyze CITE-seq and scRNA-seq data. The novelties of SECANT include (1) using confident cell type label identified from surface protein data as guidance for cell clustering, (2) providing general annotation of confident cell types for each cell cluster, (3) utilizing cells with uncertain or missing cell type label to increase performance, and (4) accurate prediction of confident cell types for scRNA-seq data. Besides, as a model-based approach, SECANT can quantify the uncertainty of the results through easily interpretable posterior probability, and our framework can be potentially extended to handle other types of multi-omics data. We successfully demonstrated the validity and advantages of SECANT via simulation studies and analysis of public and in-house datasets from multiple tissues. We believe this new method will be complementary to existing tools for characterizing novel cell types and make new biological discoveries using single-cell multi-omics data.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Comprehensive benchmarking of CITE-seq versus DOGMA-seq single cell multimodal omics

The recently developed method TEA-seq and similar DOGMA-seq single cell trimodal omics assays provide unprecedented opportunities for understanding cell biology, but independent evaluation is lacking. We explore the utility of DOGMA-seq compared to the bimodal CITE-seq assay in activated and stimulated human peripheral blood T cells. We find that single cell trimodal omics measurements after digitonin (DIG) permeabilization were generally better than after an alternative “low-loss lysis” (LLL) permeabilization condition. Next, we find that DOGMA-seq with optimized DIG permeabilization and its ATAC library provides more information, although its mRNA and cell surface protein libraries have slightly inferior quality, compared to CITE-seq.

A Role for CXCR3 Ligands as Biomarkers of Post-Operative Crohn's Disease Recurrence

BACKGROUND AND AIMS

Crohn's disease [CD] recurrence following ileocolic resection [ICR] is common. We sought to identify blood-based biomarkers associated with CD recurrence.

METHODS

CD patients undergoing ICR were recruited across six centres. Serum samples were obtained at post-operative colonoscopy. A multiplex immunoassay was used to analyse 92 inflammation-related proteins [Olink Proteomics]. Bayesian analysis was used to identify proteins associated with increasing Rutgeerts score. Identified proteins were used in receiver operating characteristic [ROC] analysis to examine the ability to identify CD recurrence [Rutgeerts score ≥i2]. Existing single cell data were interrogated to further elucidate the role of the identified proteins.

RESULTS

Data from 276 colonoscopies in 213 patients were available. Median time from surgery to first and second colonoscopy was 7 (interquartile range [IQR] 6-9) and 19 [IQR 16-23] months, respectively. Disease recurrence was evident at 60 [30%] first and 36 [49%] second colonoscopies. Of 14 proteins significantly associated with Rutgeerts score, the strongest signal was seen for CXCL9 and MMP1. Among patients on anti-tumour necrosis factor drugs, CXCL9 and CXCL11 were most strongly associated with Rutgeerts score. Both are CXCR3 ligands. Incorporation of identified proteins into ROC analysis improved the ability to identify disease recurrence as compared to C-reactive protein alone: area under the curve [AUC] 0.75 (95% confidence interval [CI]: 0.66-0.82] vs 0.64 [95% CI 0.56-0.72], p = 0.012. Single cell transcriptomic data provide evidence that innate immune cells are the primary source of the identified proteins.

CONCLUSIONS

CXCR3 ligands are associated with CD recurrence following ICR. Incorporation of novel blood-based candidate biomarkers may aid in identification of CD recurrence.

Inflamed Ulcerative Colitis Regions Associated With MRGPRX2-Mediated Mast Cell Degranulation and Cell Activation Modules, Defining a New Therapeutic Target

BACKGROUND & AIMS

Recent literature has implicated a key role for mast cells in murine models of colonic inflammation, but their role in human ulcerative colitis (UC) is not well established. A major advance has been the identification of mrgprb2 (human orthologue, MRGPX2) as mediating IgE-independent mast cell activation. We sought to define mechanisms of mast cell activation and MRGPRX2 in human UC.

METHODS

Colon tissues were collected from patients with UC for bulk RNA sequencing and lamina propria cells were isolated for MRGPRX2 activation studies and single-cell RNA sequencing. Genetic association of all protein-altering G-protein coupled receptor single-nucleotide polymorphism was performed in an Ashkenazi Jewish UC case-control cohort. Variants of MRGPRX2 were transfected into Chinese hamster ovary (CHO) and human mast cell (HMC) 1.1 cells to detect genotype-dependent effects on β-arrestin recruitment, IP-1 accumulation, and phosphorylated extracellular signal-regulated kinase.

RESULTS

Mast cell-specific mediators and adrenomedullin (proteolytic precursor of PAMP-12, an MRGPRX2 agonist) are up-regulated in inflamed compared to uninflamed UC. MRGPRX2 stimulation induces carboxypeptidase secretion from inflamed UC. Of all protein-altering GPCR alleles, a unique variant of MRGPRX2, Asn62Ser, was most associated with and was bioinformatically predicted to alter arrestin recruitment. We validated that the UC protective serine allele enhances β-arrestin recruitment, decreases IP-1, and increases phosphorylated extracellular signal-regulated kinase with MRGPRX2 agonists. Single-cell RNA sequencing defines that adrenomedullin is expressed by activated fibroblasts and epithelial cells and that interferon gamma is a key upstream regulator of mast cell gene expression.

CONCLUSION

Inflamed UC regions are distinguished by MRGPRX2-mediated activation of mast cells, with decreased activation observed with a UC-protective genetic variant. These results define cell modules of UC activation and a new therapeutic target.

Common and Rare Variant Prediction and Penetrance of IBD in a Large, Multi-ethnic, Health System-based Biobank Cohort

BACKGROUND AND AIMS

Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology.

METHODS

PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population.

RESULTS

Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation.

CONCLUSIONS

Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.

Whole-genome sequencing of African Americans implicates differential genetic architecture in inflammatory bowel disease

Whether or not populations diverge with respect to the genetic contribution to risk of specific complex diseases is relevant to understanding the evolution of susceptibility and origins of health disparities. Here, we describe a large-scale whole-genome sequencing study of inflammatory bowel disease encompassing 1,774 affected individuals and 1,644 healthy control Americans with African ancestry (African Americans). Although no new loci for inflammatory bowel disease are discovered at genome-wide significance levels, we identify numerous instances of differential effect sizes in combination with divergent allele frequencies. For example, the major effect at PTGER4 fine maps to a single credible interval of 22 SNPs corresponding to one of four independent associations at the locus in European ancestry individuals but with an elevated odds ratio for Crohn disease in African Americans. A rare variant aggregate analysis implicates Ca²⁺-binding neuro-immunomodulator CALB2 in ulcerative colitis. Highly significant overall overlap of common variant risk for inflammatory bowel disease susceptibility between individuals with African and European ancestries was observed, with 41 of 241 previously known lead variants replicated and overall correlations in effect sizes of 0.68 for combined inflammatory bowel disease. Nevertheless, subtle differences influence the performance of polygenic risk scores, and we show that ancestry-appropriate weights significantly improve polygenic prediction in the highest percentiles of risk. The median amount of variance explained per locus remains the same in African and European cohorts, providing evidence for compensation of effect sizes as allele frequencies diverge, as expected under a highly polygenic model of disease.

BREM-SC: a bayesian random effects mixture model for joint clustering single cell multi-omics data

Droplet-based single cell transcriptome sequencing (scRNA-seq) technology, largely represented by the 10× Genomics Chromium system, is able to measure the gene expression from tens of thousands of single cells simultaneously. More recently, coupled with the cutting-edge Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq), the droplet-based system has allowed for immunophenotyping of single cells based on cell surface expression of specific proteins together with simultaneous transcriptome profiling in the same cell. Despite the rapid advances in technologies, novel statistical methods and computational tools for analyzing multi-modal CITE-Seq data are lacking. In this study, we developed BREM-SC, a novel Bayesian Random Effects Mixture model that jointly clusters paired single cell transcriptomic and proteomic data. Through simulation studies and analysis of public and in-house real data sets, we successfully demonstrated the validity and advantages of this method in fully utilizing both types of data to accurately identify cell clusters. In addition, as a probabilistic model-based approach, BREM-SC is able to quantify the clustering uncertainty for each single cell. This new method will greatly facilitate researchers to jointly study transcriptome and surface proteins at the single cell level to make new biological discoveries, particularly in the area of immunology.

GMM-Demux: sample demultiplexing, multiplet detection, experiment planning, and novel cell-type verification in single cell sequencing

Identifying and removing multiplets are essential to improving the scalability and the reliability of single cell RNA sequencing (scRNA-seq). Multiplets create artificial cell types in the dataset. We propose a Gaussian mixture model-based multiplet identification method, GMM-Demux. GMM-Demux accurately identifies and removes multiplets through sample barcoding, including cell hashing and MULTI-seq. GMM-Demux uses a droplet formation model to authenticate putative cell types discovered from a scRNA-seq dataset. We generate two in-house cell-hashing datasets and compared GMM-Demux against three state-of-the-art sample barcoding classifiers. We show that GMM-Demux is stable and highly accurate and recognizes 9 multiplet-induced fake cell types in a PBMC dataset.

Challenges in IBD Research: Precision Medicine

Precision medicine is part of five focus areas of the Challenges in IBD research document, which also includes preclinical human IBD mechanisms, environmental triggers, novel technologies, and pragmatic clinical research. The Challenges in IBD Research document provides a comprehensive overview of current gaps in inflammatory bowel diseases (IBD) research and delivers actionable approaches to address them. It is the result of a multidisciplinary input from scientists, clinicians, patients, and funders, and represents a valuable resource for patient centric research prioritization. In particular, the precision medicine section is focused on highlighting the main gap areas that must be addressed to get closer to treatments tailored to the biological and clinical characteristics of each patient, which is the aim of precision medicine. The main gaps were identified in: 1) understanding and predicting the natural history of IBD: disease susceptibility, activity, and behavior; 2) predicting disease course and treatment response; and 3) optimizing current and developing new molecular technologies. Suggested approaches to bridge these gaps include prospective longitudinal cohort studies to identify and validate precision biomarkers for prognostication of disease course, and prediction and monitoring of treatment response. To achieve this, harmonization across studies is key as well as development of standardized methods and infrastructure. The implementation of state-of-the-art molecular technologies, systems biology and machine learning approaches for multi-omics and clinical data integration and analysis will be also fundamental. Finally, randomized biomarker-stratified trials will be critical to evaluate the clinical utility of validated signatures and biomarkers in improving patient outcomes and cost-effective care.

Functional variants in the LRRK2 gene confer shared effects on risk for Crohn's disease and Parkinson's disease

Crohn's disease (CD), a form of inflammatory bowel disease, has a higher prevalence in Ashkenazi Jewish than in non-Jewish European populations. To define the role of nonsynonymous mutations, we performed exome sequencing of Ashkenazi Jewish patients with CD, followed by array-based genotyping and association analysis in 2066 CD cases and 3633 healthy controls. We detected association signals in the LRRK2 gene that conferred risk for CD (N2081D variant, P = 9.5 × 10^-10) or protection from CD (N551K variant, tagging R1398H-associated haplotype, P = 3.3 × 10^-8). These variants affected CD age of onset, disease location, LRRK2 activity, and autophagy. Bayesian network analysis of CD patient intestinal tissue further implicated LRRK2 in CD pathogenesis. Analysis of the extended LRRK2 locus in 24,570 CD cases, patients with Parkinson's disease (PD), and healthy controls revealed extensive pleiotropy, with shared genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The LRRK2 N2081D CD risk allele is located in the same kinase domain as G2019S, a mutation that is the major genetic cause of familial and sporadic PD. Like the G2019S mutation, the N2081D variant was associated with increased kinase activity, whereas neither N551K nor R1398H variants on the protective haplotype altered kinase activity. We also confirmed that R1398H, but not N551K, increased guanosine triphosphate binding and hydrolyzing enzyme (GTPase) activity, thereby deactivating LRRK2. The presence of shared LRRK2 alleles in CD and PD provides refined insight into disease mechanisms and may have major implications for the treatment of these two seemingly unrelated diseases.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Appendectomy does not decrease the risk of future colectomy in UC: results from a large cohort and meta-analysis

OBJECTIVES

Early appendectomy is inversely associated with the development of UC. However, the impact of appendectomy on the clinical course of UC is controversial, generally favouring a milder disease course. We aim to describe the effect appendectomy has on the disease course of UC with focus on the timing of appendectomy in relation to UC diagnosis.

DESIGN

Using the National Institute of Diabetes and Digestive and Kidney Diseases Inflammatory Bowel Disease Genetics Consortium database of patients with UC, the risk of colectomy was compared between patients who did and did not undergo appendectomy. In addition, we performed a meta-analysis of studies that examined the association between appendectomy and colectomy.

RESULTS

2980 patients with UC were initially included. 111 (4.4%) patients with UC had an appendectomy; of which 63 were performed prior to UC diagnosis and 48 after diagnosis. In multivariable analysis, appendectomy performed at any time was an independent risk factor for colectomy (OR 1.9, 95% CI 1.1 to 3.1), with appendectomy performed after UC diagnosis most strongly associated with colectomy (OR 2.2, 95% CI 1.1 to 4.5). An updated meta-analysis showed appendectomy performed either prior to or after UC diagnosis had no effect on colectomy rates.

CONCLUSIONS

Appendectomy performed at any time in relation to UC diagnosis was not associated with a decrease in severity of disease. In fact, appendectomy after UC diagnosis may be associated with a higher risk of colectomy. These findings question the proposed use of appendectomy as treatment for UC.

Fine-mapping inflammatory bowel disease loci to single-variant resolution

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Re: Genome-Wide Association Study Identifies African-Specific Susceptibility Loci in African Americans With Inflammatory Bowel Disease

BACKGROUND & AIMS

The inflammatory bowel diseases (IBD) ulcerative colitis (UC) and Crohn’s disease (CD) cause significant morbidity and are increasing in prevalence among all populations, including African Americans. More than 200 susceptibility loci have been identified in populations of predominantly European ancestry, but few loci have been associated with IBD in other ethnicities.

METHODS

We performed 2 high-density, genome-wide scans comprising 2345 cases of African Americans with IBD (1646 with CD, 583 with UC, and 116 inflammatory bowel disease unclassified) and 5002 individuals without IBD (controls, identified from the Health Retirement Study and Kaiser Permanente database). Single-nucleotide polymorphisms (SNPs) associated at P < 5.0 × 10⁻⁸ in meta-analysis with a nominal evidence (P < .05) in each scan were considered to have genome-wide significance.

RESULTS

We detected SNPs at HLA-DRB1, and African-specific SNPs at ZNF649 and LSAMP, with associations of genome-wide significance for UC. We detected SNPs at USP25 with associations of genome-wide significance for IBD. No associations of genome-wide significance were detected for CD. In addition, 9 genes previously associated with IBD contained SNPs with significant evidence for replication (P < 1.6 × 10⁻⁶): ADCY3, CXCR6, HLA-DRB1 to HLA-DQA1 (genome-wide signifi-cance on conditioning), IL12B, PTGER4, and TNC for IBD; IL23R, PTGER4, and SNX20 (in strong linkage disequilibrium with NOD2) for CD; and KCNQ2 (near TNFRSF6B) for UC. Several of these genes, such as TNC (near TNFSF15), CXCR6, and genes associated with IBD at the HLA locus, contained SNPs with unique association patterns with African-specific alleles.

CONCLUSIONS

We performed a genome-wide association study of African Americans with IBD and identified loci associated with UC in only this population; we also replicated IBD, CD, and UC loci identified in European populations. The detection of variants associated with IBD risk in only people of African descent demonstrates the importance of studying the genetics of IBD and other complex diseases in populations beyond those of European ancestry.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications

OBJECTIVES

Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA.

METHODS

We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes.

RESULTS

The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes.

CONCLUSIONS

The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

Genome-Wide Association Study Identifies African-Specific Susceptibility Loci in African Americans With Inflammatory Bowel Disease

BACKGROUND & AIMS

The inflammatory bowel diseases (IBD) ulcerative colitis (UC) and Crohn's disease (CD) cause significant morbidity and are increasing in prevalence among all populations, including African Americans. More than 200 susceptibility loci have been identified in populations of predominantly European ancestry, but few loci have been associated with IBD in other ethnicities.

METHODS

We performed 2 high-density, genome-wide scans comprising 2345 cases of African Americans with IBD (1646 with CD, 583 with UC, and 116 inflammatory bowel disease unclassified) and 5002 individuals without IBD (controls, identified from the Health Retirement Study and Kaiser Permanente database). Single-nucleotide polymorphisms (SNPs) associated at P < 5.0 × 10^-8 in meta-analysis with a nominal evidence (P < .05) in each scan were considered to have genome-wide significance.

RESULTS

We detected SNPs at HLA-DRB1, and African-specific SNPs at ZNF649 and LSAMP, with associations of genome-wide significance for UC. We detected SNPs at USP25 with associations of genome-wide significance for IBD. No associations of genome-wide significance were detected for CD. In addition, 9 genes previously associated with IBD contained SNPs with significant evidence for replication (P < 1.6 × 10^-6): ADCY3, CXCR6, HLA-DRB1 to HLA-DQA1 (genome-wide significance on conditioning), IL12B,PTGER4, and TNC for IBD; IL23R, PTGER4, and SNX20 (in strong linkage disequilibrium with NOD2) for CD; and KCNQ2 (near TNFRSF6B) for UC. Several of these genes, such as TNC (near TNFSF15), CXCR6, and genes associated with IBD at the HLA locus, contained SNPs with unique association patterns with African-specific alleles.

CONCLUSIONS

We performed a genome-wide association study of African Americans with IBD and identified loci associated with UC in only this population; we also replicated IBD, CD, and UC loci identified in European populations. The detection of variants associated with IBD risk in only people of African descent demonstrates the importance of studying the genetics of IBD and other complex diseases in populations beyond those of European ancestry.

A Frameshift in CSF2RB Predominant Among Ashkenazi Jews Increases Risk for Crohn's Disease and Reduces Monocyte Signaling via GM-CSF

BACKGROUND & AIMS

Crohn's disease (CD) has the highest prevalence in Ashkenazi Jewish populations. We sought to identify rare, CD-associated frameshift variants of high functional and statistical effects.

METHODS

We performed exome sequencing and array-based genotype analyses of 1477 Ashkenazi Jewish individuals with CD and 2614 Ashkenazi Jewish individuals without CD (controls). To validate our findings, we performed genotype analyses of an additional 1515 CD cases and 7052 controls for frameshift mutations in the colony-stimulating factor 2-receptor β common subunit gene (CSF2RB). Intestinal tissues and blood samples were collected from patients with CD; lamina propria leukocytes were isolated and expression of CSF2RB and granulocyte-macrophage colony-stimulating factor-responsive cells were defined by adenomatous polyposis coli (APC) time-of-flight mass cytometry (CyTOF analysis). Variants of CSF2RB were transfected into HEK293 cells and the expression and functions of gene products were compared.

RESULTS

In the discovery cohort, we associated CD with a frameshift mutation in CSF2RB (P = 8.52 × 10(-4)); the finding was validated in the replication cohort (combined P = 3.42 × 10(-6)). Incubation of intestinal lamina propria leukocytes with granulocyte-macrophage colony-stimulating factor resulted in high levels of phosphorylation of signal transducer and activator of transcription (STAT5) and lesser increases in phosphorylation of extracellular signal-regulated kinase and AK straining transforming (AKT). Cells co-transfected with full-length and mutant forms of CSF2RB had reduced pSTAT5 after stimulation with granulocyte-macrophage colony-stimulating factor, compared with cells transfected with control CSF2RB, indicating a dominant-negative effect of the mutant gene. Monocytes from patients with CD who were heterozygous for the frameshift mutation (6% of CD cases analyzed) had reduced responses to granulocyte-macrophage colony-stimulating factor and markedly decreased activity of aldehyde dehydrogenase; activity of this enzyme has been associated with immune tolerance.

CONCLUSIONS

In a genetic analysis of Ashkenazi Jewish individuals, we associated CD with a frameshift mutation in CSF2RB. Intestinal monocytes from carriers of this mutation had reduced responses to granulocyte-macrophage colony-stimulating factor, providing an additional mechanism for alterations to the innate immune response in individuals with CD.

A Pleiotropic Missense Variant in SLC39A8 Is Associated With Crohn's Disease and Human Gut Microbiome Composition

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.

Genetic Predictors of Benign Course of Ulcerative Colitis-A North American Inflammatory Bowel Disease Genetics Consortium Study

BACKGROUND

A subset of patients with ulcerative colitis (UC) have a benign course and an overall favorable prognosis. Early identification of these low-risk patients may allow for a less aggressive therapeutic approach and possible reduction of therapy-associated risks. The aim of this project was to identify the genetic predictors of benign UC phenotype.

METHODS

UC patients were selected from the National Institute of Diabetes and Digestive and Kidney Diseases Inflammatory Bowel Disease Genetics Consortium. Benign phenotype was defined as no need for immunomodulatory or biological therapy, hospitalizations, or colectomy. The association between benign UC phenotype and known loci linked to the risk of inflammatory bowel disease (IBD) was evaluated. The results for 156 index single-nucleotide polymorphisms (SNPs) from the known IBD loci were extracted for the main analysis. The association of the benign phenotype to a genetic burden score was also evaluated.

RESULTS

None of the index SNPs from the IBD loci reached the predefined threshold of 1 × 10. In the exploratory analysis of the remaining Immunochip SNPs and imputed major histocompatibility complex data, 5 distinct suggestive association signals are identified (rs1697950, rs2523639, rs17836409, rs11742854, and rs75001121).

CONCLUSIONS

No SNPs from IBD susceptibility loci were found to be associated (at our predefined threshold of 1 × 10) with a benign UC disease course. The rs11742570 variant on chromosome 5 was the one with the greatest association to benign disease although the association did not reach the predefined significant threshold. Given the modest power of our study, the findings suggested on the exploratory analysis merit extension to larger discovery cohorts.

A protein-truncating R179X variant in RNF186 confers protection against ulcerative colitis

Protein-truncating variants protective against human disease provide in vivo validation of therapeutic targets. Here we used targeted sequencing to conduct a search for protein-truncating variants conferring protection against inflammatory bowel disease exploiting knowledge of common variants associated with the same disease. Through replication genotyping and imputation we found that a predicted protein-truncating variant (rs36095412, p.R179X, genotyped in 11,148 ulcerative colitis patients and 295,446 controls, MAF=up to 0.78%) in RNF186, a single-exon ring finger E3 ligase with strong colonic expression, protects against ulcerative colitis (overall P=6.89 × 10(-7), odds ratio=0.30). We further demonstrate that the truncated protein exhibits reduced expression and altered subcellular localization, suggesting the protective mechanism may reside in the loss of an interaction or function via mislocalization and/or loss of an essential transmembrane domain.

Inherited determinants of Crohn's disease and ulcerative colitis phenotypes: a genetic association study

BACKGROUND

Crohn's disease and ulcerative colitis are the two major forms of inflammatory bowel disease; treatment strategies have historically been determined by this binary categorisation. Genetic studies have identified 163 susceptibility loci for inflammatory bowel disease, mostly shared between Crohn's disease and ulcerative colitis. We undertook the largest genotype association study, to date, in widely used clinical subphenotypes of inflammatory bowel disease with the goal of further understanding the biological relations between diseases.

METHODS

This study included patients from 49 centres in 16 countries in Europe, North America, and Australasia. We applied the Montreal classification system of inflammatory bowel disease subphenotypes to 34,819 patients (19,713 with Crohn's disease, 14,683 with ulcerative colitis) genotyped on the Immunochip array. We tested for genotype-phenotype associations across 156,154 genetic variants. We generated genetic risk scores by combining information from all known inflammatory bowel disease associations to summarise the total load of genetic risk for a particular phenotype. We used these risk scores to test the hypothesis that colonic Crohn's disease, ileal Crohn's disease, and ulcerative colitis are all genetically distinct from each other, and to attempt to identify patients with a mismatch between clinical diagnosis and genetic risk profile.

FINDINGS

After quality control, the primary analysis included 29,838 patients (16,902 with Crohn's disease, 12,597 with ulcerative colitis). Three loci (NOD2, MHC, and MST1 3p21) were associated with subphenotypes of inflammatory bowel disease, mainly disease location (essentially fixed over time; median follow-up of 10·5 years). Little or no genetic association with disease behaviour (which changed dramatically over time) remained after conditioning on disease location and age at onset. The genetic risk score representing all known risk alleles for inflammatory bowel disease showed strong association with disease subphenotype (p=1·65 × 10(-78)), even after exclusion of NOD2, MHC, and 3p21 (p=9·23 × 10(-18)). Predictive models based on the genetic risk score strongly distinguished colonic from ileal Crohn's disease. Our genetic risk score could also identify a small number of patients with discrepant genetic risk profiles who were significantly more likely to have a revised diagnosis after follow-up (p=6·8 × 10(-4)).

INTERPRETATION

Our data support a continuum of disorders within inflammatory bowel disease, much better explained by three groups (ileal Crohn's disease, colonic Crohn's disease, and ulcerative colitis) than by Crohn's disease and ulcerative colitis as currently defined. Disease location is an intrinsic aspect of a patient's disease, in part genetically determined, and the major driver to changes in disease behaviour over time.

FUNDING

International Inflammatory Bowel Disease Genetics Consortium members funding sources (see Acknowledgments for full list).

Genomic architecture of inflammatory bowel disease in five families with multiple affected individuals

Currently, the best clinical predictor for inflammatory bowel disease (IBD) is family history. Over 163 sequence variants have been associated with IBD in genome-wide association studies, but they have weak effects and explain only a fraction of the observed heritability. It is expected that additional variants contribute to the genomic architecture of IBD, possibly including rare variants with effect sizes larger than the identified common variants. Here we applied a family study design and sequenced 38 individuals from five families, under the hypothesis that families with multiple IBD-affected individuals harbor one or more risk variants that (i) are shared among affected family members, (ii) are rare and (iii) have substantial effect on disease development. Our analysis revealed not only novel candidate risk variants but also high polygenic risk scores for common known risk variants in four out of the five families. Functional analysis of our top novel variant in the remaining family, a rare missense mutation in the ubiquitin ligase TRIM11, suggests that it leads to increased nuclear factor of kappa light chain enhancer in B-cells (NF-κB) signaling. We conclude that an accumulation of common weak-effect variants accounts for the high incidence of IBD in most, but not all families we analyzed and that a family study design can identify novel rare variants conferring risk for IBD with potentially large effect size, such as the TRIM11 p.H414Y mutation.

Ocular Manifestations in Inflammatory Bowel Disease Are Associated with Other Extra-intestinal Manifestations, Gender, and Genes Implicated in Other Immune-related Traits

BACKGROUND

There has been considerable progress in identifying inflammatory bowel disease [IBD] susceptibility genes but little progress in examining the role of genetic variation in the development of the extra-intestinal manifestations [EIMs] of IBD. This study identified clinical, serological, and genetic factors associated with ocular EIMs [O-EIMs] in IBD.

METHODS

We performed a retrospective case-control study of IBD patients, comparing those with and without O-EIMs using the Cedars-Sinai IBD Research Repository and the NIDDK IBD Genetics Consortium Repository. Genotyping was performed using Illumina whole genome platforms.

RESULTS

In all, 124 cases and 3328 controls with available clinical data were identified; 103 cases and 2808 controls had genetic data available. Erythema nodosum and peripheral arthritis particularly were common in patients with O-EIMs [p = 2.77 x 10(-13) and p = 2.58 x 10(-13), respectively] with increasing odds ratios for O-EIMs with each additional non-ocular-EIM [for ≥ 2 EIMs, odds ratio 14.72]. Nominal association with O-EIMs was observed at several known IBD susceptibility single nuclear polymorphisms. One locus, containing RBM19, achieved genome-wide level of significance for association with O-EIMs.

CONCLUSIONS

In IBD, O-EIMs co-occur with musculoskeletal and skin manifestations and, in this study, are nominally associated with known IBD loci. Additional cohorts are needed to verify these results and identify additional genes.

Characterization of genetic loci that affect susceptibility to inflammatory bowel diseases in African Americans

BACKGROUND & AIMS

Inflammatory bowel disease (IBD) has familial aggregation in African Americans (AAs), but little is known about the molecular genetic susceptibility. Mapping studies using the Immunochip genotyping array expand the number of susceptibility loci for IBD in Caucasians to 163, but the contribution of the 163 loci and European admixture to IBD risk in AAs is unclear. We performed a genetic mapping study using the Immunochip to determine whether IBD susceptibility loci in Caucasians also affect risk in AAs and identify new associated loci.

METHODS

We recruited AAs with IBD and without IBD (controls) from 34 IBD centers in the United States; additional controls were collected from 4 other Immunochip studies. Association and admixture loci were mapped for 1088 patients with Crohn's disease, 361 with ulcerative colitis, 62 with IBD type unknown, and 1797 controls; 130,241 autosomal single-nucleotide polymorphisms (SNPs) were analyzed.

RESULTS

The strongest associations were observed between ulcerative colitis and HLA rs9271366 (P = 7.5 × 10(-6)), Crohn's disease and 5p13.1 rs4286721 (P = 3.5 × 10(-6)), and IBD and KAT2A rs730086 (P = 2.3 × 10(-6)). Additional suggestive associations (P < 4.2 × 10(-5)) were observed between Crohn's disease and IBD and African-specific SNPs in STAT5A and STAT3; between IBD and SNPs in IL23R, IL12B, and C2orf43; and between ulcerative colitis and SNPs near HDAC11 and near LINC00994. The latter 3 loci have not been previously associated with IBD, but require replication. Established Caucasian associations were replicated in AAs (P < 3.1 × 10(-4)) at NOD2, IL23R, 5p15.3, and IKZF3. Significant admixture (P < 3.9 × 10(-4)) was observed for 17q12-17q21.31 (IZKF3 through STAT3), 10q11.23-10q21.2, 15q22.2-15q23, and 16p12.2-16p12.1. Network analyses showed significant enrichment (false discovery rate <1 × 10(-5)) in genes that encode members of the JAK-STAT, cytokine, and chemokine signaling pathways, as well those involved in pathogenesis of measles.

CONCLUSIONS

In a genetic analysis of 3308 AA IBD cases and controls, we found that many variants associated with IBD in Caucasians also showed association evidence with these diseases in AAs; we also found evidence for variants and loci not previously associated with IBD. The complex genetic factors that determine risk for or protection against IBD in different populations require further study.

High-density mapping of the MHC identifies a shared role for HLA-DRB1*01:03 in inflammatory bowel diseases and heterozygous advantage in ulcerative colitis

Genome-wide association studies of the related chronic inflammatory bowel diseases (IBD) known as Crohn’s disease and ulcerative colitis have shown strong evidence of association to the major histocompatibility complex (MHC). This region encodes a large number of immunological candidates, including the antigen-presenting classical HLA molecules¹. Studies in IBD have indicated that multiple independent associations exist at HLA and non-HLA genes, but lacked the statistical power to define the architecture of association and causal alleles^2,3. To address this, we performed high-density SNP typing of the MHC in >32,000 patients with IBD, implicating multiple HLA alleles, with a primary role for HLA-DRB1*01:03 in both Crohn’s disease and ulcerative colitis. Significant differences were observed between these diseases, including a predominant role of class II HLA variants and heterozygous advantage observed in ulcerative colitis, suggesting an important role of the adaptive immune response to the colonic environment in the pathogenesis of IBD.

Gene-gene and gene-environment interactions in ulcerative colitis

Genome-wide association studies (GWAS) have identified at least 133 ulcerative colitis (UC) associated loci. The role of genetic factors in clinical practice is not clearly defined. The relevance of genetic variants to disease pathogenesis is still uncertain because of not characterized gene-gene and gene-environment interactions. We examined the predictive value of combining the 133 UC risk loci with genetic interactions in an ongoing inflammatory bowel disease (IBD) GWAS. The Wellcome Trust Case-Control Consortium (WTCCC) IBD GWAS was used as a replication cohort. We applied logic regression (LR), a novel adaptive regression methodology, to search for high-order interactions. Exploratory genotype correlations with UC sub-phenotypes [extent of disease, need of surgery, age of onset, extra-intestinal manifestations and primary sclerosing cholangitis (PSC)] were conducted. The combination of 133 UC loci yielded good UC risk predictability [area under the curve (AUC) of 0.86]. A higher cumulative allele score predicted higher UC risk. Through LR, several lines of evidence for genetic interactions were identified and successfully replicated in the WTCCC cohort. The genetic interactions combined with the gene-smoking interaction significantly improved predictability in the model (AUC, from 0.86 to 0.89, P = 3.26E-05). Explained UC variance increased from 37 to 42 % after adding the interaction terms. A within case analysis found suggested genetic association with PSC. Our study demonstrates that the LR methodology allows the identification and replication of high-order genetic interactions in UC GWAS datasets. UC risk can be predicted by a 133 loci and improved by adding gene-gene and gene-environment interactions.

Deep resequencing of GWAS loci identifies rare variants in CARD9, IL23R and RNF186 that are associated with ulcerative colitis

Genome-wide association studies and follow-up meta-analyses in Crohn's disease (CD) and ulcerative colitis (UC) have recently identified 163 disease-associated loci that meet genome-wide significance for these two inflammatory bowel diseases (IBD). These discoveries have already had a tremendous impact on our understanding of the genetic architecture of these diseases and have directed functional studies that have revealed some of the biological functions that are important to IBD (e.g. autophagy). Nonetheless, these loci can only explain a small proportion of disease variance (~14% in CD and 7.5% in UC), suggesting that not only are additional loci to be found but that the known loci may contain high effect rare risk variants that have gone undetected by GWAS. To test this, we have used a targeted sequencing approach in 200 UC cases and 150 healthy controls (HC), all of French Canadian descent, to study 55 genes in regions associated with UC. We performed follow-up genotyping of 42 rare non-synonymous variants in independent case-control cohorts (totaling 14,435 UC cases and 20,204 HC). Our results confirmed significant association to rare non-synonymous coding variants in both IL23R and CARD9, previously identified from sequencing of CD loci, as well as identified a novel association in RNF186. With the exception of CARD9 (OR = 0.39), the rare non-synonymous variants identified were of moderate effect (OR = 1.49 for RNF186 and OR = 0.79 for IL23R). RNF186 encodes a protein with a RING domain having predicted E3 ubiquitin-protein ligase activity and two transmembrane domains. Importantly, the disease-coding variant is located in the ubiquitin ligase domain. Finally, our results suggest that rare variants in genes identified by genome-wide association in UC are unlikely to contribute significantly to the overall variance for the disease. Rather, these are expected to help focus functional studies of the corresponding disease loci.

Association between variants of PRDM1 and NDP52 and Crohn's disease, based on exome sequencing and functional studies

BACKGROUND & AIMS

Genome-wide association studies (GWAS) have identified 140 Crohn's disease (CD) susceptibility loci. For most loci, the variants that cause disease are not known and the genes affected by these variants have not been identified. We aimed to identify variants that cause CD through detailed sequencing, genetic association, expression, and functional studies.

METHODS

We sequenced whole exomes of 42 unrelated subjects with CD and 5 healthy subjects (controls) and then filtered single nucleotide variants by incorporating association results from meta-analyses of CD GWAS and in silico mutation effect prediction algorithms. We then genotyped 9348 subjects with CD, 2868 subjects with ulcerative colitis, and 14,567 control subjects and associated variants analyzed in functional studies using materials from subjects and controls and in vitro model systems.

RESULTS

We identified rare missense mutations in PR domain-containing 1 (PRDM1) and associated these with CD. These mutations increased proliferation of T cells and secretion of cytokines on activation and increased expression of the adhesion molecule L-selectin. A common CD risk allele, identified in GWAS, correlated with reduced expression of PRDM1 in ileal biopsy specimens and peripheral blood mononuclear cells (combined P = 1.6 × 10(-8)). We identified an association between CD and a common missense variant, Val248Ala, in nuclear domain 10 protein 52 (NDP52) (P = 4.83 × 10(-9)). We found that this variant impairs the regulatory functions of NDP52 to inhibit nuclear factor κB activation of genes that regulate inflammation and affect the stability of proteins in Toll-like receptor pathways.

CONCLUSIONS

We have extended the results of GWAS and provide evidence that variants in PRDM1 and NDP52 determine susceptibility to CD. PRDM1 maps adjacent to a CD interval identified in GWAS and encodes a transcription factor expressed by T and B cells. NDP52 is an adaptor protein that functions in selective autophagy of intracellular bacteria and signaling molecules, supporting the role of autophagy in the pathogenesis of CD.

Dense genotyping of immune-related disease regions identifies nine new risk loci for primary sclerosing cholangitis

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.

Relationship between proximal Crohn's disease location and disease behavior and surgery: a cross-sectional study of the IBD Genetics Consortium

OBJECTIVES

In classifying Crohn's disease (CD) location, proximal (L4) disease includes esophagogastroduodenal (EGD) and jejunal disease. Our aim was to determine the influence of proximal disease on outcomes of behavior and need for surgery and to determine if there was significant clinical heterogeneity between EGD and jejunal disease.

METHODS

We performed a cross-sectional query of the NIDDK (National Institute of Diabetes and Digestive and Kidney Disease) Inflammatory Bowel Disease Genetics Consortium (IBDGC) database of patients with a confirmed diagnosis of CD and phenotyped per the IBDGC manual. Presence of any L4, L4-EGD, L4-jejunal, and non-L4 disease (L1-ileal, L2-colonic, and L3-ileocolonic) was compared with demographic features including age, race, ethnicity, smoking and inflammatory bowel disease (IBD) family history, diagnosis age, disease duration, clinical outcomes of inflammatory, stricturing or penetrating behavior, and CD abdominal surgeries. Univariate and multivariable analyses were performed with R.

RESULTS

Among 2,105 patients with complete disease location data, 346 had L4 disease (175 L4-EGD, 115 L4-jejunal, and 56 EGD and jejunal) with 321 having concurrent L1-L3 disease. In all, 1,759 had only L1-L3 disease. L4 vs. non-L4 patients were more likely (P<0.001) to be younger at diagnosis, non-smokers, have coexisting ileal involvement, and have stricturing disease. L4-jejunal vs. L4-EGD patients were at least twice as likely (P<0.001) to have had ileal disease, stricturing behavior, and any or multiple abdominal surgeries. Remarkably, L4-jejunal patients had more (P<0.001) stricturing behavior and multiple abdominal surgeries than non-L4 ileal disease patients. Logistic regression showed stricturing risks were ileal (without proximal) site (odds ratio (OR) 3.18; 95% confidence interval 2.23-4.64), longer disease duration (OR 1.33/decade; 1.19-1.49), jejunal site (OR 2.90; 1.89-4.45), and older age at diagnosis (OR 1.21/decade; 1.10-1.34). Multiple surgery risks were disease duration (OR 3.74/decade; 3.05-4.64), penetrating disease (OR 2.60; 1.64-4.21), and jejunal site (OR 2.39; 1.36-4.20), with short duration from diagnosis to first surgery protective (OR 0.87/decade to first surgery; 0.84-0.90).

CONCLUSIONS

Jejunal disease is a significantly greater risk factor for stricturing disease and multiple abdominal surgeries than either EGD or ileal (without proximal) disease. The Montreal site classification should be revised to include separate designations for jejunal and EGD disease.

PTGER4 expression-modulating polymorphisms in the 5p13.1 region predispose to Crohn's disease and affect NF-κB and XBP1 binding sites

Background

Genome-wide association studies identified a PTGER4 expression-modulating region on chromosome 5p13.1 as Crohn's disease (CD) susceptibility region. The study aim was to test this association in a large cohort of patients with inflammatory bowel disease (IBD) and to elucidate genotypic and phenotypic interactions with other IBD genes.

Methodology/Principal Findings

A total of 7073 patients and controls were genotyped: 844 CD and 471 patients with ulcerative colitis and 1488 controls were analyzed for the single nucleotide polymorphisms (SNPs) rs4495224 and rs7720838 on chromosome 5p13.1. The study included two replication cohorts of North American (CD: n = 684; controls: n = 1440) and of German origin (CD: n = 1098; controls: n = 1048). Genotype-phenotype, epistasis and transcription factor binding analyses were performed. In the discovery cohort, an association of rs4495224 (p = 4.10×10⁻⁵; 0.76 [0.67–0.87]) and of rs7720838 (p = 6.91×10⁻⁴; 0.81 [0.71–0.91]) with susceptibility to CD was demonstrated. These associations were confirmed in both replication cohorts. In silico analysis predicted rs4495224 and rs7720838 as essential parts of binding sites for the transcription factors NF-κB and XBP1 with higher binding scores for carriers of the CD risk alleles, providing an explanation of how these SNPs might contribute to increased PTGER4 expression. There was no association of the PTGER4 SNPs with IBD phenotypes. Epistasis detected between 5p13.1 and ATG16L1 for CD susceptibility in the discovery cohort (p = 5.99×10⁻⁷ for rs7720838 and rs2241880) could not be replicated in both replication cohorts arguing against a major role of this gene-gene interaction in the susceptibility to CD.

Conclusions/Significance

We confirmed 5p13.1 as a major CD susceptibility locus and demonstrate by in silico analysis rs4495224 and rs7720838 as part of binding sites for NF-κB and XBP1. Further functional studies are necessary to confirm the results of our in silico analysis and to analyze if changes in PTGER4 expression modulate CD susceptibility.

Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease

Crohn’s disease (CD) and ulcerative colitis (UC), 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 (GWAS) and subsequent meta-analyses of CD and UC^2,3 as separate phenotypes implicated previously unsuspected mechanisms, such as autophagy⁴, in pathogenesis and showed that some IBD loci are shared with other inflammatory diseases⁵. Here we expand knowledge of relevant pathways by undertaking a meta-analysis of CD and UC genome-wide association scans, with validation of significant findings in 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 and balancing 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 striking 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.

Contribution of higher risk genes and European admixture to Crohn's disease in African Americans

BACKGROUND

African Americans (AAs) are an admixed population of West African (WA) and European ancestry (EA). Crohn's disease (CD) susceptibility genes have not been established. We therefore evaluated the contribution of European admixture and major established risk genes to AA CD.

METHODS

Ninety-seven admixture informative markers were genotyped for ancestry estimates using STRUCTURE. Overall, 354 AA CD cases and 354 ethnicity-matched controls were genotyped for total 21 single nucleotide polymorphisms (SNPs) in ATG16L1, NOD2, IBD5, IL23R and IRGM by TaqMan or direct sequencing. Association was evaluated by logistic regression, adjusted for ancestry.

RESULTS

Mean EA was similar among the CD cases and controls (20.9% and 20.4%, respectively, P = 0.58). No significant admixture differences were observed among 211 to 227 cases stratified by phenotypic subclassifications including onset, surgery, site, and behavior. CD was associated with NOD2 carrier (6.93% CD, 2.15% Controls, P = 0.007), ATG16L1 Thr300Ala (36.1% CD, 29.3% Controls, P = 0.003), SLC22A4 and SLC22A5 (IBD5 locus) functional SNPs (Leu503Phe [10.5% CD, 7.6% Controls, P = 0.05] and g-207c [41.3% CD, 35.7% Controls, P = 0.03], respectively), and IL23R rs2201841 (18.2% CD, 13.8% Controls, P = 0.03), but not IRGM variants, nor three African ancestral NOD2 nonsynonymous variants. IBD5 risk was recessive. An all-minor allele IBD5 haplotype from EA was associated (P = 0.05), whereas a more common haplotype isolating g-207c was not.

CONCLUSIONS

Specific functional gene variations contribute significantly to AA CD risk. Established NOD2, SLC22A4-A5, and ATG16L1 variants show increased CD risk, with IBD5 recessive. Although CD is more common in whites, European admixture is similar among AA cases and controls.

NADPH oxidase complex and IBD candidate gene studies: identification of a rare variant in NCF2 that results in reduced binding to RAC2

OBJECTIVE

The NOX2 NADPH oxidase complex produces reactive oxygen species and plays a critical role in the killing of microbes by phagocytes. Genetic mutations in genes encoding components of the complex result in both X-linked and autosomal recessive forms of chronic granulomatous disease (CGD). Patients with CGD often develop intestinal inflammation that is histologically similar to Crohn's colitis, suggesting a common aetiology for both diseases. The aim of this study is to determine if polymorphisms in NOX2 NADPH oxidase complex genes that do not cause CGD are associated with the development of inflammatory bowel disease (IBD).

METHODS

Direct sequencing and candidate gene approaches were used to identify susceptibility loci in NADPH oxidase complex genes. Functional studies were carried out on identified variants. Novel findings were replicated in independent cohorts.

RESULTS

Sequence analysis identified a novel missense variant in the neutrophil cytosolic factor 2 (NCF2) gene that is associated with very early onset IBD (VEO-IBD) and subsequently found in 4% of patients with VEO-IBD compared with 0.2% of controls (p=1.3×10(-5), OR 23.8 (95% CI 3.9 to 142.5); Fisher exact test). This variant reduced binding of the NCF2 gene product p67(phox) to RAC2. This study found a novel genetic association of RAC2 with Crohn's disease (CD) and replicated the previously reported association of NCF4 with ileal CD.

CONCLUSION

These studies suggest that the rare novel p67(phox) variant results in partial inhibition of oxidase function and are associated with CD in a subgroup of patients with VEO-IBD; and suggest that components of the NADPH oxidase complex are associated with CD.

Amino acid position 11 of HLA-DRβ1 is a major determinant of chromosome 6p association with ulcerative colitis

The major histocompatibility complex (MHC) on chromosome 6p is an established risk locus for ulcerative colitis (UC) and Crohn's disease (CD). We aimed to better define MHC association signals in UC and CD by combining data from dense single-nucleotide polymorphism (SNP) genotyping and from imputation of classical human leukocyte antigen (HLA) types, their constituent SNPs and corresponding amino acids in 562 UC, 611 CD and 1428 control subjects. Univariate and multivariate association analyses were performed, controlling for ancestry. In univariate analyses, absence of the rs9269955 C allele was strongly associated with risk for UC (P = 2.67 × 10(-13)). rs9269955 is a SNP in the codon for amino acid position 11 of HLA-DRβ1, located in the P6 pocket of the HLA-DR antigen binding cleft. This amino acid position was also the most significantly UC-associated amino acid in omnibus tests (P = 2.68 × 10(-13)). Multivariate modeling identified rs9269955-C and 13 other variants in best predicting UC vs control status. In contrast, there was only suggestive association evidence between the MHC and CD. Taken together, these data demonstrate that variation at HLA-DRβ1, amino acid 11 in the P6 pocket of the HLA-DR complex antigen binding cleft is a major determinant of chromosome 6p association with UC.

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn's disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD-susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2-4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.

Association of TNFSF15 polymorphism with irritable bowel syndrome

BACKGROUND

Irritable bowel syndrome (IBS) is the most common gastrointestinal disorder, affecting more than 10% of the general population worldwide. Although a genetic component is suspected, unambiguous susceptibility genes have so far not been identified. This study tested the hypothesis that genes contributing to epithelial barrier integrity, control of mucosal immune responses and interactions with bacteria in the gut are associated with IBS.

METHODS

Single nucleotide polymorphisms (SNPs) corresponding to top signals of association with Crohn's disease at 30 known susceptibility loci were tested for their effect on IBS risk in 1992 individuals from two independent case-control cohorts from Sweden and the USA. Association tests included a conservative Bonferroni correction for multiple comparisons, and were also performed on specific subgroups of patients characterised by constipation (IBS-C), diarrhoea (IBS-D) or alternating constipation and diarrhoea (IBS-A).

RESULTS

The Crohn's disease risk allele rs4263839 G in the TNFSF15 gene was significantly associated with an increased risk of both IBS (p=2.2×10(-5); OR 1.37) and more pronouncedly, IBS-C (p=8.7×10(-7); OR 1.79) in the entire sample. Similar associations and risk effects of the same magnitude were observed in the two cohorts analysed separately. A correlation between rs4263839 genotype and TNFSF15 mRNA expression was detected both in peripheral blood and in rectal mucosal biopsies from healthy individuals (combined p=0.0033).

CONCLUSIONS

TNFSF15 is a susceptibility gene for IBS and IBS constipation. As TL1A, the protein encoded by TNFSF15, contributes to the modulation of inflammatory responses, the results support a role of immune activation in IBS.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

Single nucleotide polymorphisms that increase expression of the guanosine triphosphatase RAC1 are associated with ulcerative colitis

BACKGROUND & AIMS

RAC1 is a guanosine triphosphatase that has an evolutionarily conserved role in coordinating immune defenses, from plants to mammals. Chronic inflammatory bowel diseases are associated with dysregulation of immune defenses. We studied the role of RAC1 in inflammatory bowel diseases using human genetic and functional studies and animal models of colitis.

METHODS

We used a candidate gene approach to HapMap-Tag single nucleotide polymorphisms in a discovery cohort; findings were confirmed in 2 additional cohorts. RAC1 messenger RNA expression was examined from peripheral blood cells of patients. Colitis was induced in mice with conditional disruption of Rac1 in phagocytes by administration of dextran sulfate sodium.

RESULTS

We observed a genetic association between RAC1 with ulcerative colitis in a discovery cohort, 2 independent replication cohorts, and in combined analysis for the single nucleotide polymorphisms rs10951982 (P(combined UC) = 3.3 × 10(-8), odds ratio = 1.43 [95% confidence interval: 1.26-1.63]) and rs4720672 (P(combined UC) = 4.7 × 10(-6), odds ratio = 1.36 [95% confidence interval: 1.19-1.58]). Patients with inflammatory bowel disease who had the rs10951982 risk allele had increased expression of RAC1 compared to those without this allele. Conditional disruption of Rac1 in macrophage and neutrophils of mice protected against dextran sulfate sodium-induced colitis.

CONCLUSIONS

Human studies and knockout mice demonstrated a role for the guanosine triphosphatase RAC1 in the development of ulcerative colitis; increased expression of RAC1 was associated with susceptibility to colitis.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies (GWAS) and candidate gene studies in ulcerative colitis (UC) have identified 18 susceptibility loci. We conducted a meta-analysis of 6 UC GWAS, comprising 6,687 cases and 19,718 controls, and followed-up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P<5×10^-8), increasing the number of UC associated loci to 47. After annotating associated regions using GRAIL, eQTL data and correlations with non-synonymous SNPs, we identified many candidate genes providing potentially important insights into disease pathogenesis, including IL1R2, IL8RA/B, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease (IBD) risk loci is now 99, including a minimum of 28 shared association signals between Crohn’s disease (CD) and UC.