Intestinal Inflammation Modulates the Expression of ACE2 and TMPRSS2 and Potentially Overlaps With the Pathogenesis of SARS-CoV-2-related Disease.

BACKGROUND AND AIMS: The presence of gastrointestinal symptoms and high levels of viral RNA in the stool suggest active severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication within enterocytes.
METHODS: Here, in multiple, large cohorts of patients with inflammatory bowel disease (IBD), we have studied the intersections between Coronavirus Disease 2019 (COVID-19), intestinal inflammation, and IBD treatment.
RESULTS: A striking expression of ACE2 on the small bowel enterocyte brush border supports intestinal infectivity by SARS-CoV-2. Commonly used IBD medications, both biologic and nonbiologic, do not significantly impact ACE2 and TMPRSS2 receptor expression in the uninflamed intestines. In addition, we have defined molecular responses to COVID-19 infection that are also enriched in IBD, pointing to shared molecular networks between COVID-19 and IBD.
CONCLUSIONS: These data generate a novel appreciation of the confluence of COVID-19- and IBD-associated inflammation and provide mechanistic insights supporting further investigation of specific IBD drugs in the treatment of COVID-19. Preprint doi: https://doi.org/10.1101/2020.05.21.109124.

Background and context

ACE2 and TMPRSS2, receptors that mediate SARS-CoV-2 entry are abundantly expressed in the intestines. The effect of IBD and its treatment on ACE2 and TMPRSS2 is unclear. Additionally, intersections between IBD and COVID-19 are not well defined.

New findings

Common IBD medications have region-specific effect on ACE2 and TMPRSS2 expression. Patients with COVID-19 share overlapping gene signatures with sites of inflammation in IBD.

Impact

These data support the continued clinically indicated use of IBD medications during the pandemic and suggest a potential role for these medications in the treatment of COVID-19.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing coronavirus disease 2019 (COVID-19) have evolved into a global pandemic of unprecedented proportions.

Angiotensin-converting enzyme-2 (ACE2) is a carboxypeptidase that catalyzes the conversion of angiotensin I into angiotensin 1–9, and angiotensin II into angiotensin 1–7.3, 4, 5, 6 ACE2 can also be cleaved by serine proteases such as transmembrane serine protease (TMPRSS) 2, TMPRSS11D, and TMPRSS1. Early events in the pathogenesis of SARS-CoV-2 infection include attachment of the receptor binding domain of the viral spike (S) protein to epithelial ACE2.7, 8, 9, 10 The S protein is then cleaved by TMPRSS2, which facilitates viral entry into the cytoplasm of the host cell. Following infection with SARS-CoV, ACE2 is downregulated in the lungs, resulting in unopposed renin-angiotensin-aldosterone system and contributing to disease severity. ,

Inflammatory bowel diseases (IBD) encompassing Crohn’s disease (CD) and ulcerative colitis (UC) are chronic, inflammatory disorders of the gastrointestinal (GI) tract that are treated with conventional immunosuppressive drugs such as corticosteroids, biologic therapies and immunomodulatory drugs. , Given that SARS-CoV-2 co-opts receptors expressed by intestinal epithelial cells, COVID-19 has the potential to intersect with the pathogenesis of IBD and by extension, its treatment, at a number of points. , For example, ACE2 expression may be potentially altered during gut inflammation or by IBD medications. Further, immunomodulatory drugs used in IBD therapeutics , could potentially be used in COVID-19 patients to manage the “cytokine storm” associated with severe disease.

Therefore, in this study we systematically examined potential areas of intersection between the uninflamed and inflamed GI tract and COVID-19 disease. The results of this study may improve our molecular understanding of how COVID-19 intersects IBD and may provide a rationale for further investigation of drugs used in IBD therapeutics for use in patients with COVID-19.

Methods

Immunofluorescence Microscopy

Specimens were obtained via clinical endoscopy during routine care (Supplementary Tables 1 and 2). Tissue was formalin fixed and paraffin embedded by the clinical pathology core at our institution. Primary antibodies used included ACE2 (abcam-ab15348, 1:1000), EPCAM (abcam-ab228023, prediluted), and mouse anti-TMPRSS2 (Millipore-MABF2158, 1:500) and staining was performed as detailed in supplementary methods.

Study Cohorts

Cross-sectional cohorts

The Mount Sinai Crohn’s and Colitis Registry (MSCCR): Peripheral blood and biopsy whole transcriptome sequencing data were obtained from a cross-sectional cohort (∼1200 patients) that was enrolled in the MSCCR (December 2013–September 2016); protocol approved by the Icahn School of Medicine at Mount Sinai Institutional Review Board. Detailed in supplementary methods.

The RISK cohort: ACE2 and TMPRSS2 in patients with treatment-free pediatric CD (<17 years of age) was studied using RNA-sequencing (RNA-seq) expression profiles from GSE57945, which includes ileal biopsies from endoscopically defined inflamed samples (n = 160), non-inflamed (n = 53), and non-IBD controls (n = 42).

Longitudinal cohorts

The GSE100833 series, which includes expression profiles from the gut of 80 patients with anti-tumor necrosis factor α refractory CD and the blood from 226 patients enrolled in a phase 2b crossover trial (CERTIFI trial) with ustekinumab (detailed in supplementary methods).

The GSE73661 series includes gene expression profiles (Affymetrix Human Gene 1.0 ST arrays) from colonic biopsies from patients with moderate-to-severe UC enrolled in 2 Vedolizumab efficacy trials (GEMINI-I and GEMINI LTS) (detailed in supplementary methods). GSE73661 series also included 12 non-IBD colonic biopsies and colonic biopsies from 23 patients with UC before and 4 to 6 weeks after first infliximab treatment. Response was defined as endoscopic mucosal healing.

MSCCR Bayesian Gene Regulatory Network Generation

Bayesian Gene Regulatory Networks (BGRNs) can capture fundamental properties of complex systems in states that give rise to complex (diseased) phenotypes. BGRNs were generated from intestinal biopsy RNA sequence data (MSCCR, using intestinal expression QTL information as priors). The BGRNs were region- (ileum or colon/rectum) and disease- (CD, UC, and control) specific and included both inflamed and uninflamed biopsies and were constructed using RIMBAnet software and visualized using Cytoscape 3.7. We also used 2 publicly available BGRNs from the RISK and the CERTIFI cohort (Supplementary Methods).

Bayesian Gene Regulatory Subnetwork Generation

ACE2 and TMPRSS2 subnetworks: Gene-centric subnetworks were generated by selecting either ACE2 or TMPRSS2 from various BGRNs and expanding out 3 to 5 layers (undirected) to obtain the nearest ACE2 or TMPRSS2 neighbors. The connected subnetworks obtained were generally between 200 and 500 genes in total.

IBD Inflammation, COVID-19, and IBD Drug Response-Subnetwork Generation

We curated RNA-seq–based IBD and COVID-19 response signatures by identifying differentially expressed genes (supplementary methods). Genes differentially expressed in blood, lung NHBE/A549, or human small intestinal organoids (hSIO) following SARS-CoV-2 infection; IBD inflammation; or response to medications were separately projected onto various BGRNs allowing for 1 or 2 nearest neighbors depending on the signature sizes. The most connected subnetworks were then extracted to generate model-specific SARS-CoV-2 infection–, IBD inflammation–, or drug-response–associated subnetworks (supplementary methods).

Pathway and Geneset Enrichment Analysis of Subnetworks

Gene subnetworks were tested for functional enrichment using a Fisher’s exact test with Benjamini-Hochberg (BH) multiple test correction on a collection of genesets. The collection of genesets included (1) Reactome pathways sourced from Enrichr, (2) gene sets from Smillie et al, (3) Huang et al, (4) various macrophage perturbations (eg, cytokines), (5) ACE2 coexpressed genes, and (6) reported IBD GWAS genes (see supplementary methods). Pathway and geneset enrichment, as well as intersection between networks were tested using a Fisher’s exact test and P values were adjusted for multiple hypothesis with BH correction.

Key Driver Gene Analysis

Key driver analysis identifies key or “master” driver genes for a given gene set in a given BGRN (supplementary methods). Genesets for key driver analysis included those associated with NHBE-COVID-19 infection or IBD inflammation. Key driver genes (KDGs) were summarized by frequency across the networks.

Geneset Variation Analysis of SARS-CoV-2 Infection Gene Expression Signatures

COVID-19 response gene expression (from whole blood or epithelial models) was evaluated in the context of IBD-related inflammation using gene set variation analysis (GSVA). For each COVID-19 response signature, a sample-wise enrichment score was quantified from each transcriptomic profile in the MSCCR and CERTIFI cohorts using GSVA. COVID-19 response GSVA scores were then modeled to test the association with patient-derived phenotypic information.

Results

Healthy Gut Segments Express ACE2 and TMPRSS2 Proteins

To define the localization and distribution of ACE2, immunofluorescence microscopy was performed on histologically normal GI tissue in 20 adults (9 men, 11 women) and 11 children (7 boys, 4 girls) (Figure 1 , Supplementary Table 1). ACE2 expression was observed on the small intestinal surface epithelium in all subjects in a continuous distribution with the exception of occasional breaks representing the mucin from goblet cells (Figure 1 A and B). In all examined small intestinal segments, ACE2 could be detected on the crypt epithelium, although to a lesser extent than on the surface epithelium. TMPRSS2 expression was less abundant in the small bowel, and when detectable, was exclusively found on crypt epithelium (Supplementary Figure 1 B and B). There was no observable age or sex dependence of ACE2 or TMPRSS2 protein expression in the small bowel.

Figure 1

Robust ACE2 expression is found on small bowel surface epithelium in both children and adults. Representative immunofluorescence images of ACE2 (green) and EPCAM (red) counterstained with DAPI (blue) in intestinal biopsies of healthy patients, including magnified images of surface epithelium (se) and crypt epithelium (ce). (A) Duodenal biopsies, (B) terminal ileum biopsies (Goblet cell indicated with arrow and (C) colonic biopsies from indicated sites. Patient age (years) and sex (M, male; F, female) as indicated. Isotype controls and no primary controls for each segment are included on the far right of each panel. Scale bar, 100 μm.

Supplementary Figure 1

TMPRSS2 and ACE2 distribution and localization in the intestine in children and adults. Representative immunofluorescence images of ACE2 (green) and TMPRSS2 (red) counterstained with DAPI (blue) in intestinal biopsies of non-IBD patients. Magnified images of surface epithelium (se) and crypt epithelium (ce) showing only TMPRSS2 and DAPI. (A) Duodenal biopsies. (B) Terminal ileum biopsies. (C) Biopsies from indicated colonic segments. Patient age (years) and sex (M, male; F, female) indicated in the top left corner of each image. Isotype controls and no primary controls for each segment are included on the far right of each panel. Scale bar, 100 μm.

In the colon, ACE2 expression was patchy and could not be identified in every subject, in contrast to the small bowel (Figure 1 C). This inconsistency across different donors could not be readily associated with age, sex, ACE inhibitor treatment, or the colonic segment being examined. In contrast, TMPRSS2 expression was more robust in the colon and was readily detectable on both surface and crypt epithelia (Supplementary Figure 1 C). Thus, ACE2 in the healthy gut is higher in the small bowel than the large bowel and inversely, expresses more TMPRSS2 protein in the colon compared to the small bowel.

ACE2 and TMPRSS2 messenger RNA Expression Varies in Healthy or Inflamed Gut Segments by Region

Next, we examined ACE2 and TMPRSS2 messenger RNA (mRNA) levels in the intestine of non-IBD controls and patients with IBD with active and inactive disease enrolled in MSCCR (Supplementary Table 3). As the number of samples for the colon non-rectum locations was low (Supplementary Figure 3 A and B) and no discernable differences were observed, colon nonrectum biopsies were grouped together to increase statistical power. Consistent with the protein data, ACE2 gene expression was higher in the uninflamed ileum compared with the uninflamed colon or rectum. With inflammation, ileal ACE2 mRNA expression was significantly decreased compared with either uninflamed biopsies from patients with IBD or non-IBD controls. In contrast, in the rectum ACE2 mRNA expression was increased with inflammation when compared to the uninflamed patients with IBD or non-IBD controls (Figure 2 A, top). There were no significant differences by disease location or between patients with UC versus CD (Supplementary Figure 4). We further validated these results using the pediatric IBD RISK cohort where ACE2 mRNA was significantly decreased in the ileum of patients with active IBD as well (Figure 2 B).

Supplementary Figure 3

ACE2 and TMPRSS2 gut gene expression in MSCCR CD and UC patients versus healthy controls. Plots summarize the expression level of ACE2 (A) and TMPRSS2 (B) across various gut regions from CD (top) or UC (bottom), which were either endoscopically noninflamed (Non-Inf) or inflamed (Inf). Each IBD region is compared with healthy controls. Numbers at the bottom represent the number of samples (#Non-Inf vs #Control / #Inf vs #Control). Level of significance is indicated by ∗, ∗∗, or ∗∗∗ for P < .05, < .01 or < .005, respectively.

Figure 2

Expression of ACE2 and TMPRSS2 in intestinal biopsies of adult and pediatric patients with IBD and controls. (A) Box plot summarizing the normalized expression level of ACE2 (top) and TMPRSS2 (bottom) as measured in ileum or colonic and rectal biopsies from control and IBD MSCCR patients, which were either endoscopically inflamed (IBD.I) or noninflamed (Non.I). (B) Box plot summarizing the normalized expression level of ACE2 (left) and TMPRSS2 (right), as measured in ileum CD samples from the RISK pediatric cohort. Clinical characteristics of MSCCR patients with IBD and biopsies are summarized in Supplementary Tables 3 and 4.

Supplementary Figure 4

ACE2 and TMPRSS2 gene expression according to location of disease in CD patients and effects of age and gender. (A) Normalized gene expression of ACE2 (top) and TMPRSS2 (bottom) summarized according to patients with CD with L1 (purple) or L2/L3 (orange) disease separated by region. Biopsies from endoscopically defined noninflamed (left) or inflamed (right) areas are examined separately. (B) The effect of gender and age effect on ACE2 TMPRSS2 gene expression was estimated using a multivariable model with smoking, age, gender, and 2 interaction of age and gender and IBD subtype and Tissue. The coefficient for Age effect is presented for CD and UC at different regions (ileum, colon, and rectum). P values denote significance of the age effect for each disease and tissue group, +P < .1, ∗P < .05, ∗∗P < .01. (C) We estimated the marginal mean expression of ACE2 and TMPRSS2 for male and female at each region (ileum, colon, and rectum) and disease (CD and UC) group. P values are presented where significant differences between males and females were found.

The expression of TMPRSS2 in the MSCCR cohort was moderately higher in the colon compared to ileum. In both ileum and colon biopsies TMPRSS2 expression was found significantly increased in the inflamed relative to noninflamed samples, although the effect sizes were small (Figure 2 A, bottom).

With immunofluorescence microscopy (Supplementary Table 2), we could not appreciate differences in ACE2 expression in the ileum, which possibly stemmed from the elevated physiological expression of ACE2. In the colon, patchy epithelial ACE2 expression from control non-IBD controls increased in patients with IBD with inflammation and this increase was mostly evident on the crypt epithelium (Supplementary Figure 2 A).

Supplementary Figure 2

Expression of ACE2 and TMPRSS2 protein in the intestine of patients with IBD. (A) Representative immunofluorescence images of ACE2 (green) and DAPI (blue) on the left and TMPRSS2 (red) and DAPI (blue) on the right in paired inflamed and uninflamed IBD intestinal specimens. Terminal ileum from a CD patient pre-IFX (inflamed) and after IFX treatment (uninflamed). Inflamed left colon and uninflamed sigmoid colon from a patient with UC pre-biologic. Inflamed rectum from a CD patient pre-IFX and uninflamed rectum post-IFX therapy. Clinical characteristics of IBD patients and biopsies are summarized in Supplementary Table 3.

In the ileum, low intensity TMPRSS2 expression on the crypt epithelium was also comparable in inflamed or uninflamed mucosal segments from patients with IBD and in noninflamed mucosal segments from non-IBD controls. In the colon, TMPRSS2 expression appeared comparable in patients with IBD and non-IBD controls. In the rectum, TMPRSS2 expression was enhanced by inflammation (Supplementary Figure 2 A).

Age and Gender But Not Smoking Increases ACE2 mRNA in the Colon

ACE2 mRNA was higher with age in the uninflamed rectum samples. However, these effects were essentially nullified in the presence of inflammation. A positive association with age in inflamed CD ileum and a negative association in inflamed UC rectum biopsies was observed. ACE2 mRNA in the uninflamed rectum was significantly lower in male versus female individuals but no gender associations in TMPRSS2 mRNA levels were found (Supplementary Figure 4 B and C). The expression of ACE2 and TMPRSS2 was similar when comparing active smokers with nonsmokers, either between healthy controls or patients with IBD (data not shown) and no significant interactions with inflammation status, region, or other covariates were found. Thus, age and gender, but not smoking, modulates ACE2 but not TMPRSS2 mRNA expression in the IBD colon.

Nonbiological Medications: Corticosteroids, Thiopurines, and 5-aminosalicylates Reduce ACE2 and TMPRSS2 Gene Expression in the Inflamed Colon and Rectum But Not in the Ileum

We further evaluated the impact of nonbiologic and biologic medication use (self-reported) on the expression of ACE2 and TMPRSS2 mRNA (Supplementary Table 3) by propensity matching the MSCCR cohort (see supplementary methods).

In the ileum of patients with IBD (CD), corticosteroid, thiopurine, or 5-aminosalicylate use had no impact on ACE2 mRNA expression in either inflamed or uninflamed biopsies (Figure 3 A–C). In the rectum; however, a significant decrease in ACE2 mRNA expression was observed with corticosteroid use in inflamed biopsies. A similar decrease of ACE2 mRNA was noticed in thiopurine-treated noninflamed samples from the rectum. The use of corticosteroids, thiopurine, or 5-aminosalicylate did not significantly affect TMPRSS2 mRNA expression in ileum samples. However, each of these 3 medications significantly decreased TMPRSS2 mRNA expression in inflamed rectum or colon samples. Thus, corticosteroid, thiopurine, or 5-aminosalicylate attenuated ACE2 and TMPRSS2 mRNA expression in inflamed colon and rectum.

Figure 3

The effect of IBD medication use on expression of ACE2 and TMPRSS2 in intestinal biopsies of IBD and control MSCCR patients. Propensity-matched cohorts of UC and CD MSCCR patients receiving or not receiving corticosteroids (A), thiopurines (B), 5-aminosalicylates (C), or anti-TNF (D) were used to estimate the mean (±SEM) differences in ACE2 (top) and TMPRSS2 (bottom) expression between the medicated and nonmedicated groups. P values <0.1 are reported. Under the model, we estimated the change in ACE2 or TMPRSS2 gene expression between the medicated and nonmedicated group according to disease subtype (CD, UC), region (ileum, colon, rectum), and tissue type (inflamed, noninflamed). Sample sizes are in Supplementary Table 5.

Among Biologic Medications, Infliximab Reduced ACE2 Expression in the Inflamed Colon While Ustekinumab Increased ACE2 and TMPRSS2 in the Inflamed Colon in Treatment Responders

We also defined the effect of current anti–tumor necrosis factor (anti-TNF) therapy (either adalimumab or infliximab) use on ACE2 and TMPRSS2 (Figure 3 D) expression. In the ileum, patients taking anti-TNF biologics did not show significantly different ACE2 or TMPRSS2 mRNA expression compared with those not on anti-TNF medication. In the large intestine, anti-TNF users showed increased ACE2 and TMPRSS expression, particularly in the inflamed rectum. Because the use of a cross-sectional cohort to study treatment effect has its limitations, we used published datasets from clinical trial cohorts, where longitudinal patient sampling was performed and information on endoscopic responses to treatment was available.

The results from a patient cohort treated with infliximab or vedolizumab on the GEMINI LTS trial are summarized (Figure 4 A, Supplementary Figure 5). Compared with baseline ACE2 expression, post-infliximab (week 6) colonic ACE2 gene expression was significantly lower and this decrease was observed predominantly in endoscopic responders (Supplementary Figure 5). In contrast, post-vedolizumab (week 6) ACE2 gene expression did not significantly change, although a nominally significant decrease was observed in endoscopic responders (Figure 4 A, Supplementary Figure 5). Finally, neither vedolizumab nor infliximab modified TMPRSS2 expression.

Figure 4

A summary of the effects of demographics, age, gender, inflammation, and medication use on ACE2 and TMPRSS2 gene expression in the intestine. (A) Changes on colonic gene expression profiles for ACE2 (top) and TMPRSS2 (bottom) on patients treated with vedolizumab (VDZ) or infliximab (IFX). (A) Estimated marginal means for the expression (mean ± SEM) at baseline and week 4–6. (B) Changes in gut expression of ACE2 (top) and TMPRSS2 (bottom) in patients with CD treated with ustekinumab (CERTIFI cohort). Estimated marginal means for the expression (mean ± SEM) at baseline and week 6 for inflamed and noninflamed biopsies. P values denote significance of each time point compared with screening visit +(P < .1), ∗(P < .05), ∗∗(P < .01). Sample sizes are in Supplementary Table 6. (C) A summary of various contrasts of interest on ACE2 and TMPRSS2 expression in the ileum, colon, or rectum. Arrows are used to indicate the effect of various conditions on ACE2 or TMPRSS2 gene expression, with Non.I indicating noninflamed biopsy vs control and Inf indicating inflamed vs control. N/a indicates data not available. Lighter colored arrows represent nominal significance (P between 0.1 and 0.05) and darker arrows represent significance of P < .05 in the various contrasts.

Supplementary Figure 5

The effect of infliximab and vedolizumab on ACE2 and TMPRSS2 expression. (A) Changes on colonic gene expression profiles for ACE2 (top) and TMPRSS2 (bottom) on patients treated with vedolizumab (VDZ) or infliximab (IFX). Differences in endoscopic responders and nonresponders at week 4 to 6 vs baseline samples. P values denote significance of each time point compared with screening visit +P < .1, ∗P < . 05, ∗∗P < . 01. Samples sizes are in Supplementary Table 6.

To study the impact of interleukin (IL)-12/IL-23-targeting ustekinumab on ACE2 and TMPRSS2 mRNA levels in IBD, we used the CERTIFI trial cohort. , We first confirmed that, as in our MSCCR cohort, ACE2 gene expression was higher in uninflamed ileum compared with the colon regions (Supplementary Figure 6 A). Next, we observed that ileal ACE2 gene expression was significantly decreased with inflammation as was observed in the ileal samples from both MSCCR and RISK cohorts (Supplementary Figure 6 A). In colonic and rectal samples, a trend to increase ACE2 expression in inflamed biopsies as compared with non-inflamed was observed (Supplementary Figure 6 A).

Supplementary Figure 6

The effect of ustekinumab (CERTIFI cohort) on ACE2 and TMPRSS2 gene expression in the intestine. (A) Baseline differences in expression of ACE2 and TMPRSS2 between inflamed and noninflamed tissue was estimated across different regions (ileum, colon, and rectum) using a mixed-effect model with tissue and region and its interactions as fixed effects and random intercepts for each patient. P values indicated the significance of the inflamed vs noninflamed comparison. (B) Changes in gut expression of ACE2 (top) and TMPRSS2 (bottom) in patients with CD treated with ustekinumab (CERTIFI cohort). Treatment changes in expression of ACE2 and TMPRSS2 were modeled using a mixed-effect model with visit, region, tissue, and treatment and its interactions as fixed effects. Marginal estimated means are presented for patients treated with ustekinumab and placebo group at baseline and week 6 across different gut regions. P values denote significance of change at week 6 from screening, +P < . 1, ∗P < . 05. (C) As no change was observed in noninflamed biopsies, treatment effect on inflamed biopsies was compared between week 22 clinical responders and nonresponders vs baseline. +P < . 1, ∗P < . 05. The samples sizes are summarized in Supplementary Table 6.

Following ustekinumab treatment, ACE2 gene expression was increased (nominal significance) in the inflamed tissue (both small and large intestinal) after week 6 post ustekinumab as compared with the screening biopsy (Figure 4 B, Supplementary Figure 6 B) contrasting with a decrease in expression in placebo-treated patients. On examining expression by ustekinumab response at week 22, the increase of ACE2 and TMPRSS2 expression was mainly observed in the colon and was stronger in responders (fold-change = 1.79, P = .055 for ACE2 and fold-change= 2.25, P = .024 for TMPRSS2, Supplementary Figure 6 C). Thus, TNF-targeting biologics attenuated ACE2 mRNA expression in the inflamed colon, whereas IL-12/IL-23-targeting biologics increased both ACE2 and TMPRSS2 mRNA expression in the inflamed colon, particularly in therapy responders.

Overall, the effects of IBD medications on ACE2 and TMRSS2 mRNA is complex, region-specific, and drug-specific. A summary of the medication effects is provided in Figure 4 C.

Intestinal ACE2 Gene Regulatory Subnetworks Are Enriched in Metabolic Functions and Interferon Signaling in the Inflamed Colon

To identify potential functions associated with ACE2 and TMPRSS2 in the gut, we studied these genes in the context of BGRNs. These probabilistic graphical models consider all available trait data (gene expression and genotype) simultaneously to derive gene:gene casual relationships among thousands of intermediate molecular traits.

The nearest neighbors of ACE2 were extracted including genes within either 3 to 5 path lengths in each BGRN network (5 networks total) keeping the subnetwork sizes relatively similar (∼200–300 genes, Figure 5 A and D, Supplementary Table 7). Immediate neighbors of ACE2 in the ileum CD network included SLC6A19 (a known interacting partner of ACE2), and other SLC transporters but also other viral-associated receptor proteins, like DPP4 and EZR. A summary of the recurring genes across subnetworks is shown (Figure 5 B and E).

Figure 5

BGRN analysis of ACE2 reveals metabolic function and interferon pathway associations. The ACE2-associated subnetworks extracted from the (A) ileum CD BGRN (235 of 8458 nodes total) or from the (D) Colon CD BGRN (235 of 8549 total nodes). Genes (nodes) found within the first, second, or third/fourth layer are colored blue, yellow, or gray, respectively. Five ACE-associated subnetworks were generated (Supplementary Table 7). (B and E) A summary of the genes found in both ileum and colon ACE2-associated networks. Subnetwork sizes are: Control_ileum = 221; CD_ileum= 235; Control_colon = 229; CD_colon = 235 and UC_colon = 346. (C and F) Reactome pathway enrichment analysis of ACE2-associated subnetworks (with BH adj P < .1). The heatmap depicts fold enrichment. Level of significance is indicated by ∗, ∗∗, or ∗∗∗ for P value of < .1, < .05, or < .01, respectively (see Supplementary Table 8 for full results).

Functional enrichment of ACE2-associated subnetwork was interrogated in the ileum (Figure 5 C) or colon (Figure 5 F) using the Reactome database. Identified metabolic pathways included SLC-mediated transport; xenobiotic metabolism; vitamins and cofactors; and hexose transport. Additional ileum-associated pathways included amino acid and oligopeptide SLC transporters, whereas colon-associated networks included interferon and immune signaling.

To support our network approaches, we verified a significant overlap was observed between the colonic ACE2-associated subnetworks and genes reported to be correlated with ACE2 expression in colonocytes (Supplementary Figure 7 A). In addition, ACE2-subnetworks were significantly enriched in expression profiles associated with epithelial cell types, including enterocytes and absorptive cells , (Supplementary Figure 7 B and C). Colonic ACE2-associated subnetworks also co-enriched in immune cell types as well as macrophage gene signatures following various cytokine perturbations including interferon (IFN)γ/β (Supplementary Figure 7 D). In summary, our ACE2-subnetworks are a novel source of insight into the regulation and function of ACE2.

Supplementary Figure 7

ACE2-associated subnetworks reveal gut epithelial cell type enrichment. Each ACE2-associated subnetwork was interrogated for enrichment with: (A) Genes co-correlated with ACE2 in colonocytes was curated from Wang et al; (B) gut cell type associated signatures from (Huang et al, 2019) and (C) (Smillie et al, 2019). (D) Gene sets associated with various perturbations in macrophages. The number of genes within each ACE2-associated subnetwork are: control ileum = 221; CD ileum = 235; control colon = 229; CD colon = 235; and UC colon = 346. The heatmap coloring depicts the fold enrichment and the level of significance based on BH adj P value of ∗P < .1) ∗∗P < . 05, or ∗∗∗P < . 01. The full enrichment results are available in Supplementary Tables 9 to 12.

Intestinal TMPRSS2 Gene Regulatory Subnetworks Are Predominantly Enriched for Metabolic Functions

To address the function of TMRPSS2, the subnetworks associated with TMPRSS2 were extracted from the BGRNs (Supplementary Figure 8 A) allowing 3 or 4 layers to obtain similar subnetworks sizes (∼200–300 genes). TMPRSS2 was not found in the BGRN from the ileum of controls, likely due to a low expression variance, which is a filter used before BGRN network construction (Supplementary Table 13). Genes recurring in 4 of 4 subnetworks were identified (Supplementary Figure 8 B). TMPRSS2-associated subnetwork genes are enriched for cell-cell communication, tight junction interaction, O-linked glycosylation of mucins and membrane trafficking-associated pathways (Supplementary Figure 8 C). Consistent with these functions, we observed enrichment of the TMPRSS2 subnetworks in genesets associated with enterocytes, and goblet and secretory cells (Supplementary Figure 9 A and B).

Supplementary Figure 8

Bayesian gene regulatory network (BGRN) analysis of TMPRSS2 reveals gut barrier function pathway associations. (A) The TMPRSS2 subnetwork extracted from the UC Colon BGRN (490 of 8557 nodes total) was generated by including 4 additional layers of genes (undirected network expansion). Genes or nodes found within the first, second, or third/fourth layer are colored blue, yellow, or gray, respectively. In total, 4 TMPRSS2 subnetworks were generated, 1 from ileum gene expression data from MSCCR CD patients and 3 from colon (and rectum) gene expression data from MSCCR control, UC, or CD patients (Supplementary Table 13). Both inflamed and noninflamed biopsies used for the IBD networks are shown. (B) A summary of the genes found in 4 of the 4 TMPRSS2-associated networks. The number of genes within each TMPRSS2 associated subnetworks are: CD ileum = 397; Control colon = 268; CD colon = 366; and UC colon= 490. (C) Pathway enrichment analysis (Fisher’s exact test) of each TMPRSS2-associated subnetwork according to Reactome pathways was performed. Only pathways that were found significantly enriched in at least 1 TMPRSS2-associated subnetwork (at BH adj P < .1) are presented in the heatmaps in (C). The heatmap coloring depicts the fold enrichment and the level of significance based on BH adj P value of ∗ P < . 1, ∗∗ P < . 05, or ∗∗∗ P < . 01. The full enrichment results are available in Supplementary Table 14.

Supplementary Figure 9

TMPRSS2-associated subnetworks reveal gut cell type enrichment. Each of the 4 TMPRSS2-associated subnetworks were examined for enrichment with various gut cell type related signatures (Smillie et al, 2019) (A) (Huang et al, 2019). (B) The number of genes within each TMPRSS2-associated subnetworks are as follows: CD ileum = 397; control colon = 268; CD colon = 366; and UC colon = 490. The heatmap coloring depicts the fold enrichment and the level of significance based on BH adj P value of ∗P < .1, ∗∗P < . 05, or ∗∗∗P < . 01). The full enrichment results are available in Supplementary Tables 15 and 16.

A Subset of Pathways Associated With SARS-CoV2 Response and IBD Inflammation Overlap

COVID-19 is a multisystem disorder in which innate and adaptive immune cells as well as nonimmune cells likely play a role in disease pathogenesis. Therefore, apart from alterations in the receptor expression, we investigated additional areas of overlap between COVID-19 responsive pathways and pathways associated with IBD. As the first step, we examined a host molecular response signature generated following SARS-CoV-2 infection of a primary human lung epithelium (NHBE) cell and a transformed lung alveolar cell line (A549). Using GSVA, we generated a per-sample score summarizing expression of either up or down-regulated COVID-19-responsive genes, and evaluated differences in these scores according to intestinal region, disease and inflammation status (Figure 6 A, Supplementary Figure 10 A and B). Genes up-regulated by SARS-CoV-2 infection show significantly higher expression in inflamed regions as compared with uninflamed regions or non-IBD control subjects, an observation confirmed in the CERTIFI CD cohort (Supplementary Figure 10 C and D). We directly compared the genes associated with response to lung cell SARS-CoV-2 infection and various IBD-centric genesets generated in our MSCCR cohort or IBD GWAS genes. We observed that genes up-regulated with inflammation, or positively associated with macroscopic or microscopic measures of disease, or associated with the risk of IBD, were significantly enriched with genes up-regulated by SARS-CoV2 infection of lung epithelial cells (Supplementary Figure 10 E).

Figure 6

Molecular genes associated with IBD inflammation overlap with genes responsive to COVID-19 infection in lung cell model and whole blood. (A) Boxplot of the overall activity of COVID-19–responsive genes (up-regulated) as determined in A549 (top) or NHBE (bottom) lung epithelial cell models in the MSCCR biopsy samples using GSVA. Box plots summarize the mean expression (± SEM) of the 2 signatures (SF11 and sample sizes in Supplementary Table 4, ∗∗∗P < .001). (B) Box plot summarizing the expression of the COVID-19 blood signature (up-regulated genes) in the blood transcriptome data of the MSCCR showing association between IBD disease (left) and clinical severity (right) (SF11). (C) Three NHBE-COVID-19– and 3 IBD-Inf–associated subnetworks were generated one each for ileum CD, colon UC, and colon CD BGRNs (Supplementary Tables 17 and 18). The 3 ileum-associated CD networks are shown for example. Venn diagrams of the overlap in subnetwork genes including the fold enrichment (FE) and a P value of the enrichment test are shown. Using the BGRN drug response, subnetworks were generated (see supplementary methods) for infliximab and ustekinumab response genes. The bar graph summarizes the FE for the overlaps between the drug response and the (1) COVID-19 subnetworks; (2) IBD inflammation subnetworks, or (3) the intersecting nodes between COVID-19 and IBD-Inf subnetworks. A full table of network nodes and enrichment results can be found in Supplementary Tables 23 and 24. (D) The 3 sets of genes in the intersecting subnetworks from (C) were interrogated for enrichment in the Reactome database. Heatmap depicting the FE in pathways (BH adj P < .1 and minimum 5-fold enrichment, Supplementary Table 19). (E) Expression of a molecular signature consisting of genes responsive to COVID-19 infection as determined in blood transcriptome data from adult patients with IBD with CD (CERTIFI cohort). (Top) Estimated marginal means (mean ± SEM) for the activity (GSVA scores) of the COVID-19 blood signature on the blood transcriptome of CERTIFI patients with CD at baseline and after weeks 4 to 22 of treatment with ustekinumab or placebo. (Bottom) Ustekinumab-induced changes in COVID-19 blood signature between week 6 responders and nonresponders (defined as change in CDAI). P values denote significance of each time point compared with screening visit ∗P < .1, ∗∗P < .05, ∗∗∗P < .01.

Supplementary Figure 10

Geneset variation analysis of lung COVID-19 responsive genes as determined in MSCCR and CERTIFI cohort gut biopsies. (A) We evaluated expression of COVID-19–responsive genes as determined in NHBE (A) or A549 (B) lung epithelial cell models in the MSCCR biopsy samples using GSVA. Mixed-effect linear models with region, tissue type, and its interaction were used to compare COVID-19 scores among control, noninflamed, and inflamed samples (sample sizes in Supplementary Table 4, ∗∗∗P < .001). (C, D) Two molecular expression signatures reflecting a host’s transcriptional response to SARS-CoV-2 infection were curated from Blanco-Melo et al. (C) Transformed lung alveolar cell line (A549) and (D) primary human lung epithelium (NHBE) were profiled following SARS-CoV-2 exposure. Boxplots for the GSVA scores of COVID-19–responsive genes in inflamed and noninflamed biopsies across gut regions for patients in CERTIFI cohort. Mixed-effect models with region and tissue as fixed-effects were used to compare expression between noninflamed and inflamed samples. Sample sizes are presented in Supplementary Table 6. (E) A heatmap summarizing the significance (-log adj P value) for the enrichment of genes up- or down-regulated following lung cell SARS-CoV2 infection in various IBD disease–associated genesets derived from the MSCCR cohort analysis and IBD GWAS genes.

Next, we examined for the congruence of COVID-19–related peripheral blood gene responses and active IBD inflammation. We observed that genes up-regulated in the blood of COVID-19–infected patients have significantly higher expression in the blood of patients with IBD as compared with healthy control blood (Figure 6 B, Supplementary Figure 11 A) as well as in patients with active IBD versus quiescent IBD (Figure 6 B, Supplementary Figure 11 B).

Supplementary Figure 11

Expression of a blood signature of genes responsive to COVID-19 infection in blood transcriptome of adult patients with IBD and controls (MSCCR cohort). Box plot summarizing the expression of the blood signature identified in COVID-19–infected patients vs healthy controls as determined in the blood transcriptome data of (A) the MSCCR CD (n = 432) and UC (n = 389) patients and controls (n = 209) and (B) between clinically defined inactive (n = 288 CD, n = 340 UC) vs active (n = 72 CD, n = 49 UC) MSCCR IBD patients (right). GSVA scores associated with up-regulated genes are shown on the left and down-regulated genes are shown on the right. P values denote significance with ∗P < .05, ∗∗P < .01, and ∗∗∗P < .001.

To further interrogate the molecular congruence between SARS-CoV-2 infection and IBD inflammatory responses, we generated NHBE-COVID-19–associated subnetworks (NHBE-COVID_subnet) and IBD inflammation–associated subnetworks (IBD_Inf_subnet) from 3 IBD BGRNs (Supplementary Table 17). We then determined the overlap between them, with the rationale that common gene membership implies similar molecular pathobiologies. We observed a significant overlap between IBD-Inf and NHBE-COVID_subnets across all 3 independent networks tested (Figure 6 C). Interestingly, the common genes (148 in ileum CD, 213 in colon UC, and 170 in colon CD networks) were significantly enriched in candidate IBD GWAS genes (Supplementary Table 18).

To determine the underlying shared pathobiological mechanisms, we evaluated enrichment of the 3 intersecting subnetworks against Reactome pathways. Interestingly, although subnetworks were generated by projecting a COVID-19 response signature generated in an epithelial-based model, most of the pathways and cell-type enrichments were immune oriented with the most striking related to innate immune signaling via IFN and IL-6 (Figure 6 D, Supplementary Figure 12, Supplementary Table 19). This is consistent with the fact that although SARS-CoV-2 infection initiates within the epithelium, COVID-19 is a multisystem disorder in which innate and adaptive immune cells as well as nonimmune cells likely play a role in disease pathogenesis. Our COVID-19–associated gut subnetwork findings are consistent with a large body of recent data that describe a dramatic up-regulation of proinflammatory cytokines in patients with COVID-19, including induction of IFN-stimulated genes.

Supplementary Figure 12

Lung model COVID-19–associated gene subnetworks and the shared subnetwork genes with IBD-inflamed genes are enriched with immune cell types. NHBE COVID-19 subnetworks and the intersecting genes found between them and the IBD-inflamed subnetwork genes were interrogated for enrichment with gut cell type signatures. The heatmap coloring depicts the fold enrichment and the level of significance based on BH adj P value of ∗P < .1, ∗∗P < .05, or ∗∗∗P < .01.

We confirmed these observations using the recent dataset in which SARS-CoV-2 was shown to infect hSIOs (Supplementary Figure 13 A and B). Direct overlap of genesets from the SARS-CoV-2–infected hSIOs and various IBD-centric genesets showed significant gene enrichments, similar to those observed between the lung-COVID-19 model and IBD (Supplementary Figure 13 C). Finally, subnetworks generated in the ileum CD using hSIO-COVID-19 responsive genesets (Supplementary Table 20) overlapped with IBD-Inf associated subnetworks. The intersecting genes were similar to those found between IBD-Inf and NHBE-COVID19-associated ileum CD subnetworks, containing many IFN-stimulated genes (Supplementary Figure 13 D and E).

Supplementary Figure 13

Subnetworks associated with COVID-19 response in hSIOs overlaps with NHBE-COVID-19 responsive subnetworks. Two molecular expression signatures reflecting transcriptional responses to SARS-CoV-2 infection in human small intestinal organoids (hSIOs) were curated. (A) hSIOs cultured with differentiation media (DIF) contain predominantly enterocytes, goblet cells, and low numbers of enteroendocrine cells. (B) hSIOs grown in Wnt high expansion medium (EXP) consist mainly of stem cells and enterocyte progenitors. We determined expression of the hSIO COVID-19–responsive genes (at false discovery rate <0.1) in the ileum, colon, and rectum MSCCR cohort samples using geneset variation analysis (GSVA) and mixed-effect model with region, tissue type, and its interaction as fixed effects was performed to compare expression between control, noninflamed, and inflamed samples. P values are as indicated. (C) A heatmap summarizing the significance (-log adj P value) for the enrichment of genes up- or down-regulated following hSIO SARS-CoV2 infection in various IBD disease-associated genesets derived from the MSCCR cohort analysis and IBD GWAS genes. (D) A Venn diagram summarizing the overlaps of 3 subnetworks generated using the ileum CD BGRN, from projecting signatures associated with (1) IBD-Inf; (2) NHBE_COVID-19 response; or (3) hSIO (DIF media)-COVID-19 response, allowing 1 or 2 nearest neighboring genes. (E) A Venn diagram summarizing the overlaps of 3 subnetworks generated using the ileum CD BGRN, from projecting signatures associated with (1) IBD-Inf; (2) NHBE_COVID-19 response; or (3) hSIO (EXP media)-COVID-19 response, allowing 1 nearest neighboring gene. The hSIO(DIF)- and hSIO(EXP)- COVID-19 associated subnetworks are found in Supplementary Table 20. The intersecting genes are shown. Note many genes are also identified as KDGs in Figure 6E. The significance of the overlaps of various genesets is presented in each panel.

Next, as a data reduction approach, we evaluated each of the NHBE-COVID-19 and IBD-Inflammation molecular response subnetworks for KDGs. We then summarized the KDGs across the 3 BGRNs and determined those shared by IBD inflammation and NHBE-COVID-19 molecular responses (Supplementary Figure 14, Supplementary Table 21). The shared KDGs identified were interferon-stimulated genes such as IRF1, GBP1/2, and PARP14/9, and several of these were replicated in BGRNs from the RISK and CERTIFI cohort, including CXCL1, GBP4, and PARP9, among others (data not shown).

Supplementary Figure 14

Summary of KDGs for the COVID-19 and IBD-inflammation–associated colonic and ileum IBD subnetworks. KDGs were determined, separately, for NHBE-COVID-19 or IBD-Inf signature genes in each Bayesian gene regulatory network (BGRN) from Figure 6C and a summary of the top recurring KDGs are presented (see Supplementary Table 21 for full results). A gene was considered a KDG if its subnetwork (within 3 layers) was significantly enriched (adj P < .05) in signature genes.

We also compared the intersections between COVID-19–associated genes and those associated with 3 different murine models of intestinal injury or inflammation. We observed a significant overlap between genes up-regulated with COVID-19 response in the NHBE lung model with genes up-regulated in a (1) Dextran Sodium Sulfate–induced intestinal injury model; (2) Trinitrobenzene sulfonic acid-intestinal injury model, especially 2 days post TNBS administration; or (3) adoptive T-cell transfer colitis model following a 6-week time period (ie, W0IBD-mouse models and NHBE-COVID-19 response included C3, IFITM3, IL1B, S100A9, TGM2 (transglutaminase 2), and PLAUR (plasminogen activator, urokinase receptor) (Supplementary Table 22).

Finally, we generated colonic or ileal gene subnetworks associated with response to infliximab or ustekinumab (Supplementary Table 23) and tested their enrichment in the tissue and IBD-type specific subnetworks and observed significant overlap between many of the genes associated with IBD therapy and either IBD-inflammation (used as positive control) or COVID infection (Figure 6 C), indicative of commonality between COVID-19 and response to IBD treatment. We also evaluated changes in the activity of the blood COVID-19 response signature in the blood of CERTIFI CD patients treated with ustekinumab and observed that COVID-19 up-regulated genes significantly decreased after 4 weeks of ustekinumab treatment (Figure 6 F).

Altogether, through analyses of multiple genesets, we observed a significant overlap between COVID-19 response genes and genes associated with IBD response.

Discussion

The objective of this study was to systematically determine the molecular intersections among COVID-19–associated inflammation, IBD, and immunomodulatory drugs. Our data provide mRNA- and protein-level evidence of the regional distribution of ACE2 and TMPRSS2 across different parts of the GI tract and the impact of the commonly used IBD medications on ACE2 and TMPRSS2 expression in the inflamed and uninflamed intestines. In addition, our data highlight an overlap between COVID-responsive pathways and pathways associated with IBD inflammation. These findings generate the possibility that some of the current and emerging therapies in IBD may be of benefit in patients with COVID-19.

In exploring the intersections between COVID-19 and IBD, we initially considered (1) impact of inflammation on ACE2 expression and TMPRSS2 expression in the intestines; and (2) impact of IBD medications on ACE2 and TMPRSS2 expression. ACE2 has a less appreciated renin-angiotensin-aldosterone system–independent role in the intestine by promoting amino acid absorption. Consistent with this function, we observed high levels of ACE2 expression in the small bowel brush border38, 39, 40 supportive of the role of ACE2 in mucosal homeostasis. Remarkably, the analyses of gene:regulatory networks empirically derived from the terminal ileum, showed that ACE2 was coregulated with SLC6A19, an amino acid transporter that physically interacts with ACE2. The colon-derived ACE2 subnetwork would also suggest a potential metabolic role for ACE2 in the colon including solute carrier dependent processes.

Although the physiological function of TMPRSS2 remains largely elusive, it has been linked to epithelial sodium channel (ENaC) regulation, and possibly in regulating sperm function. In agreement with literature, our BGRN analysis showed SCNN1A, which corresponds to a subunit of the ENaC, associated with TMPRSS2. Remarkably, gene products colocalized with TMPRSS2 were consistently enriched in fundamental epithelial functions. For example, F11R, which encodes a JAM-A protein, mediates tight junction formation. Interestingly, we noted connections to viral responses. JAM-A is a reovirus receptor and TRIM31 and BAIAIP2L1 associate with mitochondrial antiviral-associated proteins.

We have observed a reduction in ACE2 expression in the inflamed ileum. Because ACE2 appears to be a brush-border enzyme, its reduction with inflammation in the ileum is consistent with loss of expression of other brush border enzymes during enteritis. That said, even during inflammation, the expression of ACE2 remains significant in the small intestines. In contrast, inflammation associated with IBD enhances the expression of ACE2 and TMPRSS2 in the rectum. Overall, modulation of ACE2 and TMPRSS2 by IBD-associated inflammation is complex and appears to be region-specific in the intestines.

Akin to the impact of inflammation on intestinal ACE2 and TMPRSS2 expression, the effect of IBD medications was complex. Corticosteroids, thiopurines, or 5-aminosalicylates did not significantly affect TMPRSS2 mRNA expression in the ileum. However, each of these 3 medications significantly decreased TMPRSS2 mRNA expression in inflamed rectum or colon samples. In addition, the impact of IBD medications like TNF inhibitors may vary by the stage of treatment. One may speculate that in early stages of therapy, anti-TNF drugs, similar to anti-IL-12/23 agents, may increase the expression of ACE2 in the intestines and potentially have a detrimental effect; however, these medications could have a beneficial effect in the long run because of their ability to reduce inflammation.

Importantly, although SARS-CoV-2 infection initiates with the viral attachment to ACE2 and its cleavage by TMPRSS2 within the epithelial surfaces,7, 8, 9, 10 COVID-19 is a multisystem disorder in which innate and adaptive immune cells as well as nonimmune cells likely play a role in disease pathogenesis. Therefore, apart from alterations in the receptor expression, we investigated additional areas of overlap between COVID-19 responsive pathways and pathways associated with IBD.

A number of such intersections appeared. IL6, CXCL1/2/5, PDPN, and S100A8/A9, which were up-regulated following SARS-CoV-2 infection of primary human lung epithelium (NHBE) were also significantly up-regulated in the inflamed intestines of patients in our IBD cohort (MSCCR). Similarly, a peripheral blood gene signature from patients with COVID-19, which included the up-regulated genes, CLEC4D, S100A8/A9, and FCAR, were also significantly up-regulated in the blood of MSCCR patients with IBD (as compared with controls) and in patients with active IBD (as compared with patients with quiescent disease). In addition, a number of SARS-CoV-2–associated genes were up-regulated in murine models of intestinal injury with DSS (DAPP1, PDPN, IL1RN, DUOX2, IL1B, S100A9, CXCL2, CXCL3) or TNBS (MARCKSL1, IFITIM3, IFITIM1, C3, AGR2, REG4) or adoptive T-cell transfer colitis model (TAP2, MARCKSL1, SLP1, PARP9, IFITIM3, MMP13, IL1B, S100A8). These genes relate to a number of IBD-relevant pathways including those associated with inflammatory cytokine signaling (including IL-6, IL-1, IFN-G), chemokine signaling, but also with interferon-associated pathways, regulation of complement cascade, G protein coupled receptor signaling, as well as collagen degradation.

Having observed significant molecular intersections between COVID-19– and IBD-associated pathways, we next examined the impact of biologic medications used in IBD therapeutics where pre- and post-treatment transcriptomic data was available. Interestingly, our network analyses had identified a number of shared COVID-19 and IBD-associated “key driver genes,” including CXCL1, GBP4, SOCS3, PARP14, and PARP9. Importantly, we observed that these KDGs, which were up-regulated with COVID-19 and IBD inflammation, were all down-regulated following infliximab treatment. Thus, the impact of IBD medications on attenuating some of the key inflammatory genes and pathways would be independent of their complex effects on the expression of ACE2 and TMPRSS2 on enterocytes.

Given the unregulated inflammatory responses in patients with severe COVID-19, it has been argued that targeted use of anti-inflammatory medications be considered as a therapeutic option. This approach has been met with variable success. Although the use of dexamethasone has provided a striking mortality benefit, the results of trials using anti-IL-6 have been equivocal (NCT04320615) and NCT04372186. An anti-TNF trial is currently under way in the United Kingdom (NCT04425538). Our data suggest that anti-IL-12/23 therapy could also be considered in the therapeutic armamentarium in patients with COVID-19. Reassuringly, real-time data from a registry of patients with IBD-COVID did not report significant adverse outcomes associated with the use of biologic medications, including anti-TNF and anti-IL-12/23 inhibitors. To the contrary, the risk of severe COVID-19 was found to be reduced in patients with IBD on anti-TNF inhibitor medication. This is consistent with the observation that individuals treated with cytokine inhibitors had lower rates of SARS-CoV-2 seroprevalence compared with controls.

In summary, through detailed analyses of intestinal tissues in health and IBD, we conclude that high expression of ACE2 and TMPRSS2 potentially supports local, GI-associated replication of SARS-COV-2. Further, a number of overlapping inflammatory pathways between COVID-19 and IBD are noted. These data support the use of specific anti-inflammatory agents in the treatment of patients with COVID-19.