Inflammation status modulates the effect of host genetic variation on intestinal gene expression in inflammatory bowel disease

More than 240 genetic risk loci have been associated with inflammatory bowel disease (IBD), but little is known about how they contribute to disease development in involved tissue. Here, we hypothesized that host genetic variation affects gene expression in an inflammation-dependent way, and investigated 299 snap-frozen intestinal biopsies from inflamed and non-inflamed mucosa from 171 IBD patients. RNA-sequencing was performed, and genotypes were determined using whole exome sequencing and genome wide genotyping. In total, 28,746 genes and 6,894,979 SNPs were included. Linear mixed models identified 8,881 independent intestinal cis-expression quantitative trait loci (cis-eQTLs) (FDR < 0.05) and interaction analysis revealed 190 inflammation-dependent intestinal cis-eQTLs (FDR < 0.05), including known IBD-risk genes and genes encoding immune-cell receptors and antibodies. The inflammation-dependent cis-eQTL SNPs (eSNPs) mainly interact with prevalence of immune cell types. Inflammation-dependent intestinal cis-eQTLs reveal genetic susceptibility under inflammatory conditions that can help identify the cell types involved in and the pathways underlying inflammation, knowledge that may guide future drug development and profile patients for precision medicine in IBD.

Inflammatory bowel diseases are heterogeneous, and little is known about how underlying genetic variation can affect their development. Here, the authors report that intestinal inflammation modulates the effect of host genetics on the gut mucosal expression of 190 genes in the context of inflammatory bowel diseases.

Whole exome sequencing analyses reveal gene-microbiota interactions in the context of IBD

OBJECTIVE

Both the gut microbiome and host genetics are known to play significant roles in the pathogenesis of IBD. However, the interaction between these two factors and its implications in the aetiology of IBD remain underexplored. Here, we report on the influence of host genetics on the gut microbiome in IBD.

DESIGN

To evaluate the impact of host genetics on the gut microbiota of patients with IBD, we combined whole exome sequencing of the host genome and whole genome shotgun sequencing of 1464 faecal samples from 525 patients with IBD and 939 population-based controls. We followed a four-step analysis: (1) exome-wide microbial quantitative trait loci (mbQTL) analyses, (2) a targeted approach focusing on IBD-associated genomic regions and protein truncating variants (PTVs, minor allele frequency (MAF) >5%), (3) gene-based burden tests on PTVs with MAF <5% and exome copy number variations (CNVs) with site frequency <1%, (4) joint analysis of both cohorts to identify the interactions between disease and host genetics.

RESULTS

We identified 12 mbQTLs, including variants in the IBD-associated genes IL17REL, MYRF, SEC16A and WDR78. For example, the decrease of the pathway acetyl-coenzyme A biosynthesis, which is involved in short chain fatty acids production, was associated with variants in the gene MYRF (false discovery rate <0.05). Changes in functional pathways involved in the metabolic potential were also observed in participants carrying rare PTVs or CNVs in CYP2D6, GPR151 and CD160 genes. These genes are known for their function in the immune system. Moreover, interaction analyses confirmed previously known IBD disease-specific mbQTLs in TNFSF15.

CONCLUSION

This study highlights that both common and rare genetic variants affecting the immune system are key factors in shaping the gut microbiota in the context of IBD and pinpoints towards potential mechanisms for disease treatment.

Performance of the Montreal classification for inflammatory bowel diseases

AIM: To validate the Montreal classification system for Crohn’s disease (CD) and ulcerative colitis (UC) within the Netherlands.

METHODS: A selection of 20 de-identified medical records with an appropriate representation of the inflammatory bowel disease (IBD) sub phenotypes were scored by 30 observers with different professions (gastroenterologist specialist in IBD, gastroenterologist in training and IBD-nurses) and experience level with IBD patient care. Patients were classified according to the Montreal classification. In addition, participants were asked to score extra-intestinal manifestations (EIM) and disease severity in CD based on their clinical judgment. The inter-observer agreement was calculated by percentages of correct answers (answers identical to the “expert evaluation”) and Fleiss-kappa (κ). Kappa cut-offs: < 0.4-poor; 0.41-0.6-moderate; 0.61-0.8-good; > 0.8 excellent.

RESULTS: The inter-observer agreement was excellent for diagnosis (κ = 0.96), perianal disease (κ = 0.92) and disease location in CD (κ = 0.82) and good for age of onset (κ = 0.67), upper gastrointestinal disease (κ = 0.62), disease behaviour in CD (κ = 0.79) and disease extent in UC (κ = 0.65). Disease severity in UC was scored poor (κ = 0.23). The additional items resulted in a good inter-observer agreement for EIM (κ = 0.68) and a moderate agreement for disease severity in CD (κ = 0.44). Percentages of correct answers over all Montreal items give a good reflection of the inter-observer agreement (> 80%), except for disease severity (48%-74%). IBD-nurses were significantly worse in scoring upper gastrointestinal disease in CD compared to gastroenterologists (P = 0.008) and gastroenterologists in training (P = 0.040). Observers with less than 10 years of experience were significantly better at scoring UC severity than observers with 10-20 years (P = 0.003) and more than 20 years (P = 0.003) of experience with IBD patient care. Observers with 10-20 years of experience with IBD patient care were significantly better at scoring upper gastrointestinal disease in CD than observers with less than 10 years (P = 0.007) and more than 20 years (P = 0.007) of experience with IBD patient care.

CONCLUSION: We found a good to excellent inter-observer agreement for all Montreal items except for disease severity in UC (poor).

Wnt signaling and Dupuytren's disease

BACKGROUND

Dupuytren's disease is a benign fibromatosis of the hands and fingers that leads to flexion contractures. We hypothesized that multiple genetic and environmental factors influence susceptibility to this disease and sought to identify susceptibility genes to better understand its pathogenesis.

METHODS

We conducted a genomewide association study of 960 Dutch persons with Dupuytren's disease and 3117 controls (the discovery set) to test for association between the disease and genetic markers. We tested the 35 single-nucleotide polymorphisms (SNPs) most strongly associated with Dupuytren's disease (P<1×10(-4)) in the discovery set in three additional, independent case series comprising a total of 1365 affected persons and 8445 controls from Germany, the United Kingdom, and The Netherlands.

RESULTS

Initially, we observed a significant genomewide association between Dupuytren's disease and 8 SNPs at three loci. Tests of replication and joint analysis of all data from 2325 patients with Dupuytren's disease and 11,562 controls yielded an association with 11 SNPs from nine different loci (P<5.0×10(-8)). Six of these loci contain genes known to be involved in the Wnt-signaling pathway: WNT4 (rs7524102) (P=2.8×10(-9); odds ratio, 1.28), SFRP4 (rs16879765) (P=5.6×10(-39); odds ratio, 1.98), WNT2 (rs4730775) (P=3.0×10(-8); odds ratio, 0.83), RSPO2 (rs611744) (P=7.9×10(-15); odds ratio, 0.75), SULF1 (rs2912522) (P=2.0×10(-13); odds ratio, 0.72), and WNT7B (rs6519955) (P=3.2×10(-33); odds ratio, 1.54).

CONCLUSIONS

This study implicates nine different loci involved in genetic susceptibility to Dupuytren's disease. The fact that six of these nine loci harbor genes encoding proteins in the Wnt-signaling pathway suggests that aberrations in this pathway are key to the process of fibromatosis in Dupuytren's disease.