Genomic and Expression Analyses Identify a Disease-Modifying Variant for Fibrostenotic Crohn's Disease.

Background and Aims: Crohn's disease [CD] is a chronic inflammatory disease with unpredictable behaviour. More than half of CD patients eventually develop complications such as stenosis, for which they then require endoscopic dilatation or surgery, as no anti-fibrotic drugs are currently available. We aim to identify disease-modifying genes associated with fibrostenotic CD. Methods: We performed a within-case analysis comparing 'extreme phenotypes' using the Immunochip and replication of the top single nucleotide polymorphisms [SNPs] with Agena Bioscience in two independent case-control cohorts totalling 322 cases with fibrostenotis [recurrent after surgery] and 619 cases with purely inflammatory CD. Results: Combined meta-analysis resulted in a genome-wide significant signal for SNP rs11861007 [p = 6.0910-11], located on chromosome 16, in lncRNA RP11-679B19.1, an lncRNA of unknown function, and close to exon 9 of the WWOX gene, which codes for WW domain-containing oxidoreductase. We analysed mRNA expression of TGF-ß and downstream genes in ileocecal resection material from ten patients with and without the WWOX risk allele. Patients carrying the risk allele [A] showed enhanced colonic expression of TGF-ß compared to patients homozygous for the wild-type [G] allele [p = 0.0079]. Conclusion: We have identified a variant in WWOX and in lncRNA RP11-679B19.1 as a disease-modifying genetic variant associated with recurrent fibrostenotic CD and replicated this association in an independent cohort. WWOX can potentially play a crucial role in fibrostenosis in CD, being positioned at the crossroads of inflammation and fibrosis.

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.

Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations.

Ulcerative colitis and Crohn's disease are the two main forms of inflammatory bowel disease (IBD). Here we report the first trans-ancestry association study of IBD, with genome-wide or Immunochip genotype data from an extended cohort of 86,640 European individuals and Immunochip data from 9,846 individuals of East Asian, Indian or Iranian descent. We implicate 38 loci in IBD risk for the first time. For the majority of the IBD risk loci, the direction and magnitude of effect are consistent in European and non-European cohorts. Nevertheless, we observe genetic heterogeneity between divergent populations at several established risk loci driven by differences in allele frequency (NOD2) or effect size (TNFSF15 and ATG16L1) or a combination of these factors (IL23R and IRGM). Our results provide biological insights into the pathogenesis of IBD and demonstrate the usefulness of trans-ancestry association studies for mapping loci associated with complex diseases and understanding genetic architecture across diverse populations.

Genome-wide association scan in north Indians reveals three novel HLA-independent risk loci for ulcerative colitis.

OBJECTIVE: Over 100 ulcerative colitis (UC) loci have been identified by genome-wide association studies (GWASs) primarily in Caucasians (CEUs). Many of them have weak effects on disease susceptibility, and the bulk of the heritability cannot be ascribed to these loci. Very little is known about the genetic background of UC in non-CEU groups. Here we report the first GWAS on UC in a genetically distinct north Indian (NI) population. DESIGN: A genome-wide scan was performed on 700 cases and 761 controls. 18 single-nucleotide polymorphisms (SNPs) (p<5×10(-5)) were genotyped in an independent cohort of 733 cases and 1148 controls. A linear mixed model was used for case-control association tests. RESULTS: Seven novel human leucocyte antigen (HLA)-independent SNPs from chromosome 6, located in 3.8-1, BAT2, MSH5, HSPA1L, SLC44A4, CFB and NOTCH4, exceeded p<5×10(-8) in the combined analysis. To assess the independent biological contribution of such genes from the extended HLA region, we determined the percentage alternative pathway activity of complement factor B (CFB), the top novel hit. The activity was significantly different (p=0.01) between the different genotypes at rs12614 in UC cases. Transethnic comparisons revealed a shared contribution of a fraction of UC risk genes between NI and CEU populations, in addition to genetic heterogeneity. CONCLUSIONS: This study shows varying contribution of the HLA region to UC in different populations. Different environmental exposures and the characteristic genetic structure of the HLA locus across ethnic groups collectively make it amenable to the discovery of causative alleles by transethnic resequencing. This may lead to an improved understanding of the molecular mechanisms underlying UC.

Copper metabolism domain-containing 1 represses genes that promote inflammation and protects mice from colitis and colitis-associated cancer.

BACKGROUND & AIMS: Activation of the transcription factor nuclear factor-κB (NF-κB) has been associated with the development of inflammatory bowel disease (IBD). Copper metabolism MURR1 domain containing 1 (COMMD1), a regulator of various transport pathways, has been shown to limit NF-κB activation. We investigated the roles of COMMD1 in the pathogenesis of colitis in mice and IBD in human beings. METHODS: We created mice with a specific disruption of Commd1 in myeloid cells (Mye-knockout [K/O] mice); we analyzed immune cell populations and functions and expression of genes regulated by NF-κB. Sepsis was induced in Mye-K/O and wild-type mice by cecal ligation and puncture or intraperitoneal injection of lipopolysaccharide (LPS), colitis was induced by administration of dextran sodium sulfate, and colitis-associated cancer was induced by administration of dextran sodium sulfate and azoxymethane. We measured levels of COMMD1 messenger RNA in colon biopsy specimens from 29 patients with IBD and 16 patients without (controls), and validated findings in an independent cohort (17 patients with IBD and 22 controls). We searched for polymorphisms in or near COMMD1 that were associated with IBD using data from the International IBD Genetics Consortium and performed quantitative trait locus analysis. RESULTS: In comparing gene expression patterns between myeloid cells from Mye-K/O and wild-type mice, we found that COMMD1 represses expression of genes induced by LPS. Mye-K/O mice had more intense inflammatory responses to LPS and developed more severe sepsis and colitis, with greater mortality. More Mye-K/O mice with colitis developed colon dysplasia and tumors than wild-type mice. We observed a reduced expression of COMMD1 in colon biopsy specimens and circulating leukocytes from patients with IBD. We associated single-nucleotide variants near COMMD1 with reduced expression of the gene and linked them with increased risk for ulcerative colitis. CONCLUSIONS: Expression of COMMD1 by myeloid cells has anti-inflammatory effects. Reduced expression or function of COMMD1 could be involved in the pathogenesis of IBD.

Extraintestinal manifestations and complications in inflammatory bowel disease: from shared genetics to shared biological pathways.

BACKGROUND: The clinical presentation of the inflammatory bowel diseases (IBD) is extremely heterogenous and is characterized by various extraintestinal manifestations and complications (EIM). Increasing genetic insight for IBD and EIM shows multiple shared susceptibility loci. We hypothesize that, next to these overlapping genetic risk loci, distinct disease pathways are shared between IBD and EIM. METHODS: The overlapping genetic risk loci for IBD and its EIM were searched in literature. We assessed shared disease pathways by performing an extensive pathway analysis by protein-protein interaction and cotranscriptional analysis, using both publicly available and newly developed databases. RESULTS: Reliable genetic data were available for primary sclerosing cholangitis, ankylosing spondylitis, decreased bone mineral density, colorectal carcinoma, gallstones, kidney stones, and deep venous thrombosis. We found an extensive overlap in genetic risk loci, especially for IBD and primary sclerosing cholangitis and ankylosing spondylitis. We identified 370 protein-protein interactions, of which 108 are statistically specific. We identified 446 statistically specific cotranscribed gene pairs. The interactions are shown to cluster in specific biological pathways. CONCLUSIONS: We show that the pathogenetic overlap between IBD and its EIM extends beyond shared risk genes to distinctive shared biological pathways. We define genetic background as a risk factor for IBD-EIM alongside known mechanisms such as malabsorption and medication. Clustering patients based on distinctive pathways may enable stratification of patients to predict development of EIM.

Correlation of genetic risk and messenger RNA expression in a Th17/IL23 pathway analysis in inflammatory bowel disease.

BACKGROUND: The Th17/IL23 pathway has both genetically and biologically been implicated in the pathogenesis of the inflammatory bowel diseases (IBD), Crohn's disease, and ulcerative colitis. So far, it is unknown whether and how associated risk variants affect expression of the genes encoding for Th17/IL23 pathway proteins. METHODS: Ten IBD-associated SNPs residing near Th17/IL23 genes were used to construct a genetic risk model in 753 Dutch IBD cases and 1045 controls. In an independent cohort of 40 Crohn's disease, 40 ulcerative colitis, and 40 controls, the genetic risk load and presence of IBD were correlated to quantitative PCR-generated messenger RNA (mRNA) expression of 9 representative Th17/IL23 genes in both unstimulated and PMA/CaLo stimulated peripheral blood mononuclear cells. In 1240 individuals with various immunological diseases with whole genome genotype and mRNA-expression data, we also assessed correlation between genetic risk load and differential mRNA expression and sought for SNPs affecting expression of all currently known Th17/IL23 pathway genes (cis-expression quantitative trait locus). RESULTS: The presence of IBD, but not the genetic risk load, was correlated to differential mRNA expression for IL6 in unstimulated peripheral blood mononuclear cells and to IL23A and RORC in response to stimulation. The cis-expression quantitative trait locus analysis showed little evidence for correlation between genetic risk load and mRNA expression of Th17/IL23 genes, because we identified for only 2 of 22 Th17/IL23 genes a cis-expression quantitative trait locus single nucleotide polymorphism that is also associated to IBD (STAT3 and CCR6). CONCLUSIONS: Our results suggest that only the presence of IBD and not the genetic risk load alters mRNA expression levels of IBD-associated Th17/IL23 genes.

Complex host genetics influence the microbiome in inflammatory bowel disease.

BACKGROUND: Human genetics and host-associated microbial communities have been associated independently with a wide range of chronic diseases. One of the strongest associations in each case is inflammatory bowel disease (IBD), but disease risk cannot be explained fully by either factor individually. Recent findings point to interactions between host genetics and microbial exposures as important contributors to disease risk in IBD. These include evidence of the partial heritability of the gut microbiota and the conferral of gut mucosal inflammation by microbiome transplant even when the dysbiosis was initially genetically derived. Although there have been several tests for association of individual genetic loci with bacterial taxa, there has been no direct comparison of complex genome-microbiome associations in large cohorts of patients with an immunity-related disease. METHODS: We obtained 16S ribosomal RNA (rRNA) gene sequences from intestinal biopsies as well as host genotype via Immunochip in three independent cohorts totaling 474 individuals. We tested for correlation between relative abundance of bacterial taxa and number of minor alleles at known IBD risk loci, including fine mapping of multiple risk alleles in the Nucleotide-binding oligomerization domain-containing protein 2 (NOD2) gene exon. We identified host polymorphisms whose associations with bacterial taxa were conserved across two or more cohorts, and we tested related genes for enrichment of host functional pathways. RESULTS: We identified and confirmed in two cohorts a significant association between NOD2 risk allele count and increased relative abundance of Enterobacteriaceae, with directionality of the effect conserved in the third cohort. Forty-eight additional IBD-related SNPs have directionality of their associations with bacterial taxa significantly conserved across two or three cohorts, implicating genes enriched for regulation of innate immune response, the JAK-STAT cascade, and other immunity-related pathways. CONCLUSIONS: These results suggest complex interactions between genetically altered host functional pathways and the structure of the microbiome. Our findings demonstrate the ability to uncover novel associations from paired genome-microbiome data, and they suggest a complex link between host genetics and microbial dysbiosis in subjects with IBD across independent cohorts.

Functional genomics identifies type I interferon pathway as central for host defense against Candida albicans.

Candida albicans is the most common human fungal pathogen causing mucosal and systemic infections. However, human antifungal immunity remains poorly defined. Here by integrating transcriptional analysis and functional genomics, we identified Candida-specific host defence mechanisms in humans. Candida induced significant expression of genes from the type I interferon pathway in human peripheral blood mononuclear cells. This unexpectedly prominent role of type I interferon pathway in anti-Candida host defence was supported by additional evidence. Polymorphisms in type I interferon genes modulated Candida-induced cytokine production and were correlated with susceptibility to systemic candidiasis. In in vitro experiments, type I interferons skewed Candida-induced inflammation from a Th17 response towards a Th1 response. Patients with chronic mucocutaneous candidiasis displayed defective expression of genes in the type I interferon pathway. These findings indicate that the type I interferon pathway is a main signature of Candida-induced inflammation and has a crucial role in anti-Candida host defence in humans.

HNF4α and CDH1 are associated with ulcerative colitis in a Dutch cohort.

BACKGROUND: Inflammatory bowel diseases (IBDs), consisting of ulcerative colitis (UC) and Crohn's disease (CD), are complex disorders with multiple genes contributing to disease pathogenesis. A recent genome-wide association scan identified three novel susceptibility loci for UC: HNF4α, CDH1, and LAMB1. We performed an analysis of these three loci in an independent cohort. METHODS: In all, 821 UC patients and 1260 healthy controls of central European Caucasian descent were genotyped for single nucleotide polymorphisms (SNPs): rs6017342 (HNF4α), rs1728785 (CDH1), and rs6949033 (LAMB1). Differences in allele and genotype distribution in cases and controls were tested for significance with the χ² test. RESULTS: Allelic association analysis showed that SNP rs6017342 in the HNF4α locus was strongly associated with UC (P = 1,04 × 10(-11) , odds ratio [OR] = 1.56, 95% confidence interval [CI] = 1.37-1.77) and SNP rs1728785 (CDH1) was associated with P = 0.01 (OR = 1.23, 95% CI = 1.05-1.44). SNP rs6949033 in LAMB1 was not associated in our cohort (P = 0.12, OR = 1.11, 95% CI = 0.97-1.26). We found an association for SNP rs6949033 (LAMB1) for disease limited to the rectum (P = 0.02). However, this association was lost after correcting for multiple testing. No further specific subphenotype associations were identified. CONCLUSIONS: This is the first independent study to replicate the HNF4α and CDH1 loci as susceptibility loci for UC. The main candidate genes in these risk loci play important roles in the maintenance of the integrity of the epithelial barrier, highlighting the importance of the mucosal barrier function for UC pathogenesis.

Differential association of two PTPN22 coding variants with Crohn's disease and ulcerative colitis.

BACKGROUND: The PTPN22 gene is an important risk factor for human autoimmunity. The aim of this study was to evaluate for the first time the role of the R263Q PTPN22 polymorphism in ulcerative colitis (UC) and Crohn's disease (CD), and to reevaluate the association of the R620W PTPN22 polymorphism with both diseases. METHODS: A total of 1677 UC patients, 1903 CD patients, and 3111 healthy controls from an initial case-control set of Spanish Caucasian ancestry and two independent sample sets of European ancestry (Dutch and New Zealand) were included in the study. Genotyping was performed using TaqMan SNP assays for the R263Q (rs33996649) and R620W (rs2476601) PTPN22 polymorphisms. Meta-analysis was performed on 6977 CD patients, 5695 UC patients, and 9254 controls to test the overall effect of the minor allele of R620W and R263Q polymorphisms. RESULTS: The PTPN22 263Q loss-of-function variant showed initial evidence of association with UC in the Spanish cohort (P = 0.026, odds ratio [OR] = 0.61, 95% confidence interval [CI]: 0.39-0.95), which was confirmed in the meta-analysis (P = 0.013 pooled, OR = 0.69, 95% CI: 0.51-0.93). In contrast, the 263Q allele showed no association with CD (P = 0.22 pooled, OR = 1.16, 95% CI: 0.91-1.47). We found in the pooled analysis that the PTPN22 620W gain-of-function variant was associated with reduced risk of CD (P = 7.4E-06 pooled OR = 0.81, 95% CI: 0.75-0.89) but not of UC (P = 0.88 pooled, OR = 0.98, 95% CI: 0.85-1.15). CONCLUSIONS: Our data suggest that two autoimmunity-associated polymorphisms of the PTPN22 gene are differentially associated with CD and UC. The R263Q polymorphism only associated with UC, whereas the R620W was significantly associated with only CD.

Analysis of SNPs with an effect on gene expression identifies UBE2L3 and BCL3 as potential new risk genes for Crohn's disease.

Genome-wide association studies (GWAS) for Crohn's disease (CD) have identified loci explaining approximately 20% of the total genetic risk of CD. Part of the other genetic risk loci is probably partly hidden among signals discarded by the multiple testing correction needed in the analysis of GWAS data. Strategies for finding these hidden loci require large replication cohorts and are costly to perform. We adopted a strategy of selecting SNPs for follow-up that showed a correlation to gene expression [cis-expression quantitative trait loci (eQTLs)] since these have been shown more likely to be trait-associated. First we show that there is an overrepresentation of cis-eQTLs in the known CD-associated loci. Then SNPs were selected for follow-up by screening the top 500 SNP hits from a CD GWAS data set. We identified 10 cis-eQTL SNPs. These 10 SNPs were tested for association with CD in two independent cohorts of Dutch CD patients (1539) and healthy controls (2648). In a combined analysis, we identified two cis-eQTL SNPs that were associated with CD rs2298428 in UBE2L3 (P=5.22x10(-5)) and rs2927488 in BCL3 (P=2.94x10(-4)). After adding additional publicly available data from a previously reported meta-analysis, the association with rs2298428 almost reached genome-wide significance (P=2.40x10(-7)) and the association with rs2927488 was corroborated (P=6.46x10(-4)). We have identified UBE2L3 and BCL3 as likely novel risk genes for CD. UBE2L3 is also associated with other immune-mediated diseases. These results show that eQTL-based pre-selection for follow-up is a useful approach for identifying risk loci from a moderately sized GWAS.