African Ancestry Proportion Influences Ileal Gene Expression in Inflammatory Bowel Disease

Oxidative stress gene expression, DNA methylation, and gut microbiota interaction trigger Crohn's disease: a multi-omics Mendelian randomization study

BACKGROUND

Oxidative stress (OS) is a key pathophysiological mechanism in Crohn's disease (CD). OS-related genes can be affected by environmental factors, intestinal inflammation, gut microbiota, and epigenetic changes. However, the role of OS as a potential CD etiological factor or triggering factor is unknown, as differentially expressed OS genes in CD can be either a cause or a subsequent change of intestinal inflammation. Herein, we used a multi-omics summary data-based Mendelian randomization (SMR) approach to identify putative causal effects and underlying mechanisms of OS genes in CD.

METHODS

OS-related genes were extracted from the GeneCards database. Intestinal transcriptome datasets were collected from the Gene Expression Omnibus (GEO) database and meta-analyzed to identify differentially expressed genes (DEGs) related to OS in CD. Integration analyses of the largest CD genome-wide association study (GWAS) summaries with expression quantitative trait loci (eQTLs) and DNA methylation QTLs (mQTLs) from the blood were performed using SMR methods to prioritize putative blood OS genes and their regulatory elements associated with CD risk. Up-to-date intestinal eQTLs and fecal microbial QTLs (mbQTLs) were integrated to uncover potential interactions between host OS gene expression and gut microbiota through SMR and colocalization analysis. Two additional Mendelian randomization (MR) methods were used as sensitivity analyses. Putative results were validated in an independent multi-omics cohort from the First Affiliated Hospital of Sun Yat-sen University (FAH-SYS).

RESULTS

A meta-analysis from six datasets identified 438 OS-related DEGs enriched in intestinal enterocytes in CD from 817 OS-related genes. Five genes from blood tissue were prioritized as candidate CD-causal genes using three-step SMR methods: BAD, SHC1, STAT3, MUC1, and GPX3. Furthermore, SMR analysis also identified five putative intestinal genes, three of which were involved in gene-microbiota interactions through colocalization analysis: MUC1, CD40, and PRKAB1. Validation results showed that 88.79% of DEGs were replicated in the FAH-SYS cohort. Associations between pairs of MUC1-Bacillus aciditolerans and PRKAB1-Escherichia coli in the FAH-SYS cohort were consistent with eQTL-mbQTL colocalization.

CONCLUSIONS

This multi-omics integration study highlighted that OS genes causal to CD are regulated by DNA methylation and host-microbiota interactions. This provides evidence for future targeted functional research aimed at developing suitable therapeutic interventions and disease prevention.

Wnt signaling is boosted during intestinal regeneration by a CD44-positive feedback loop

Enhancement of Wnt signaling is fundamental for stem cell function during intestinal regeneration. Molecular modules control Wnt activity by regulating signal transduction. CD44 is such a positive regulator and a Wnt target gene. While highly expressed in intestinal crypts and used as a stem cell marker, its role during intestinal homeostasis and regeneration remains unknown. Here we propose a CD44 positive-feedback loop that boosts Wnt signal transduction, thus impacting intestinal regeneration. Excision of Cd44 in Cd44^fl/fl;VillinCreER^T2 mice reduced Wnt target gene expression in intestinal crypts and affected stem cell functionality in organoids. Although the integrity of the intestinal epithelium was conserved in mice lacking CD44, they were hypersensitive to dextran sulfate sodium, and showed more severe inflammation and delayed regeneration. We localized the molecular function of CD44 at the Wnt signalosome, and identified novel DVL/CD44 and AXIN/CD44 complexes. CD44 thus promotes optimal Wnt signaling during intestinal regeneration.

Stratification of risk of progression to colectomy in ulcerative colitis via measured and predicted gene expression

An important goal of clinical genomics is to be able to estimate the risk of adverse disease outcomes. Between 5% and 10% of individuals with ulcerative colitis (UC) require colectomy within 5 years of diagnosis, but polygenic risk scores (PRSs) utilizing findings from genome-wide association studies (GWASs) are unable to provide meaningful prediction of this adverse status. By contrast, in Crohn disease, gene expression profiling of GWAS-significant genes does provide some stratification of risk of progression to complicated disease in the form of a transcriptional risk score (TRS). Here, we demonstrate that a measured TRS based on bulk rectal gene expression in the PROTECT inception cohort study has a positive predictive value approaching 50% for colectomy. Single-cell profiling demonstrates that the genes are active in multiple diverse cell types from both the epithelial and immune compartments. Expression quantitative trait locus (QTL) analysis identifies genes with differential effects at baseline and week 52 follow-up, but for the most part, differential expression associated with colectomy risk is independent of local genetic regulation. Nevertheless, a predicted polygenic transcriptional risk score (PPTRS) derived by summation of transcriptome-wide association study (TWAS) effects identifies UC-affected individuals at 5-fold elevated risk of colectomy with data from the UK Biobank population cohort studies, independently replicated in an NIDDK-IBDGC dataset. Prediction of gene expression from relatively small transcriptome datasets can thus be used in conjunction with TWASs for stratification of risk of disease complications.

Altered splicing associated with the pathology of inflammatory bowel disease

Background

Aberrant splicing of individual genes is a well-known mechanism promoting pathology for a wide range of conditions, but disease is less commonly attributed to global disruption of exon usage. To explore the possible association of aberrant splicing with inflammatory bowel disease, we developed a pipeline for quantifying transcript abundance and exon inclusion transcriptome-wide and applied it to a dataset of ileal and rectal biopsies, both obtained in duplicate from 34 pediatric or young adult cases of ulcerative colitis and Crohn’s disease.

Results

Expression and splicing covary to some extent, and eight individuals exhibited aberrant profiles that can be explained by altered ratios of epithelial to stromal and immune cells. Ancestry-related biases in alternative splicing accounting for 5% of the variance were also observed, in part also related to cell-type proportions. In addition, two individuals were identified who had 284 exons with significantly divergent percent spliced in exons, including in the established IBD risk gene CEACAM1, which caused their ileal samples to resemble the rectum.

Conclusions

These results imply that quantitative differences in splice usage contribute to the pathology of inflammatory bowel disease in a previously unrecognized manner.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00347-y.

Increased colonic expression of ACE2 associates with poor prognosis in Crohn's disease

The host receptor for SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2), is highly expressed in small intestine. Our aim was to study colonic ACE2 expression in Crohn's disease (CD) and non-inflammatory bowel disease (non-IBD) controls. We hypothesized that the colonic expression levels of ACE2 impacts CD course. We examined the expression of colonic ACE2 in 67 adult CD and 14 NIBD control patients using RNA-seq and quantitative (q) RT-PCR. We validated ACE2 protein expression and localization in formalin-fixed, paraffin-embedded matched colon and ileal tissues using immunohistochemistry. The impact of increased ACE2 expression in CD for the risk of surgery was evaluated by a multivariate regression analysis and a Kaplan–Meier estimator. To provide critical support for the generality of our findings, we analyzed previously published RNA-seq data from two large independent cohorts of CD patients. Colonic ACE2 expression was significantly higher in a subset of adult CD patients which was defined as the ACE2-high CD subset. IHC in a sampling of ACE2-high CD patients confirmed high ACE2 protein expression in the colon and ileum compared to ACE2-low CD and NIBD patients. Notably, we found that ACE2-high CD patients are significantly more likely to undergo surgery within 5 years of CD diagnosis, and a Cox regression analysis found that high ACE2 levels is an independent risk factor for surgery (OR 2.17; 95% CI, 1.10–4.26; p = 0.025). Increased intestinal expression of ACE2 is associated with deteriorated clinical outcomes in CD patients. These data point to the need for molecular stratification that can impact CD disease-related outcomes.

In vivo negative regulation of SARS-CoV-2 receptor, ACE2, by interferons and its genetic control

Background: Angiotensin I converting enzyme 2 (ACE2) is a receptor for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and differences in its expression may affect susceptibility to infection. Methods: We performed a genome-wide expression quantitative trait loci (eQTL) analysis using hepatitis C virus-infected liver tissue from 190 individuals. Results: We discovered that polymorphism in a type III interferon gene (IFNL4), which eliminates IFN-λ4 production, is associated with a two-fold increase in ACE2 RNA expression. Conversely, among genes negatively correlated with ACE2 expression, IFN-signalling pathways were highly enriched and ACE2 was downregulated after IFN-α treatment. Negative correlation was also found in the gastrointestinal tract where inflammation driven IFN-stimulated genes were negatively correlated with ACE2 expression and in lung tissue from a murine model of SARS-CoV-1 infection suggesting conserved regulation of ACE2 across tissue and species. Conclusions: We conclude that ACE2 is likely a negatively-regulated interferon-stimulated gene (ISG) and carriage of IFNL4 gene alleles which modulates ISGs expression in viral infection may play a role in SARS-CoV-2 pathogenesis with implications for therapeutic interventions.

Increased Colonic Expression of ACE2 Associates with Poor Prognosis in Crohn's disease

Background and Aims

The host receptor for SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2), is highly expressed in small intestine. Our aim was to study colonic ACE2 expression in Crohn's disease (CD) and non-inflammatory bowel disease (non-IBD) controls. We hypothesized that the colonic expression levels of ACE2 impacts CD course.

Methods

We examined the expression of colon ACE2 using RNA-seq and quantitative (q) RT-PCR from 69 adult CD and 14 NIBD control patients. In a subset of this cohort we validated ACE2 protein expression and localization in formalin-fixed, paraffin-embedded matched colon and ileal tissues using immunohistochemistry. The impact of increased ACE2 expression in CD for the risk of surgery was evaluated by a multivariate regression analysis and a Kaplan-Meier estimator. To provide critical support for the generality of our findings, we analyzed previously published RNA-seq data from two large independent cohorts of CD patients.

Results

Colonic ACE2 expression was significantly higher in a subset of adult CD patients (ACE2-high CD). IHC in a sampling of ACE2-high CD patients confirmed high ACE2 protein expression in the colon and ileum compared to ACE2-low CD and NIBD patients. Notably, we found that ACE2-high CD patients are significantly more likely to undergo surgery within 5 years of diagnosis, with a Cox regression analysis finding that high ACE2 levels is an independent risk factor (OR 2.18; 95%CI, 1.05-4.55; p=0.037).

Conclusion

Increased intestinal expression of ACE2 is associated with deteriorated clinical outcomes in CD patients. These data point to the need for molecular stratification that may impact CD disease-related outcomes.