A correction has been applied to this article in:
Corrigendum: Improved diagnosis of inflammatory bowel disease and prediction and monitoring of response to anti-TNF alpha treatment based on measurement of signal transduction pathway activity
Wilbert Bouwman
Wim Verhaegh
Anja van de Stolpe- Philips Research, Eindhoven, Netherlands
Objective: Ulcerative colitis (UC) and Crohn’s disease (CD) are two subtypes of chronic inflammatory bowel disease (IBD). Differential diagnosis remains a challenge. Anti-TNFα treatment is an important treatment for IBD, yet resistance frequently occurs and cannot be predicted. Consequently, many patients receive ineffective therapy with potentially adverse effects. Novel assays are needed to improve diagnosis, and predict and monitor response to anti-TNF-α compounds.
Design: Signal transduction pathway (STP) technology was used to quantify activity of STPs (androgen and estrogen receptor, PI3K, MAPK, TGFβ, Notch, Hedgehog, Wnt, NFκB, JAK-STAT1/2, and JAK-STAT3 pathways) in colon mucosa samples of CD and UC patients, based on transcriptome analysis. Previously described STP assay technology is based on computational inference of STP activity from mRNA levels of target genes of the STP transcription factor.
Results: Results show that NFκB, JAK-STAT3, Wnt, MAPK, and androgen receptor pathways were abnormally active in CD and UC. Colon and ileum-localized CD differed with respect to STP activity, the JAK-STAT1/2 pathway being abnormally active in ileal CD. High activity of NFκB, JAK-STAT3, and TGFβ pathways was associated with resistance to anti-TNFα treatment in UC and colon-located CD, but not in ileal CD. Abnormal STP activity decreased with successful treatment.
Conclusion: We believe that measuring mucosal STP activity provides clinically relevant information to improve differential diagnosis of IBD and prediction of resistance to anti-TNFα treatment in patients with colon-localized IBD, and provides new targets for treatment and overcoming anti-TNFα resistance.
Introduction
Inflammatory bowel disease (IBD) is an inflammatory chronic intestinal disorder. The two major types include ulcerative colitis (UC) and Crohn’s disease (CD). Crohn’s disease affects the small intestine and large intestine, as well as the mouth, esophagus, stomach and anus, whereas ulcerative colitis primarily affects colon and rectum. Differential diagnosis with respect to CD and UC remains a challenge, and in 15% of cases, definitive discrimination is not possible. Correct diagnosis is important to decide on the most suitable therapy including systemic drug treatment. Systemic therapy is directed at inhibition of inflammation and the auto-immune process, healing of the intestinal wall, and correction of the microbiome (Hazel and O’Connor, 2020). Multiple types of drugs can be used, e.g., aminosalicylic acids, corticosteroids, and tumor necrosis factor alpha (TNFα) blockers, but around one-third of patients will be primary resistant to therapy and three-quarters of responders develop secondary resistance within 1 year of treatment (Feuerstein et al., 2021), (Colombel et al., 2010), (D’Haens et al., 2008), (D’Haens, 2010). As a consequence, many patients receive ineffective therapy with potentially severe adverse effects, including serious infections. Unfortunately, biomarker assays to predict response to anti-TNF-α drugs, to measure the extent of the response, and to identify emerging drug resistance are lacking, precluding a personalized treatment approach (D’Haens and Deventer, 2021). Novel assays to improve diagnosis, and to predict and monitor response to anti-TNF-α compounds are needed. New drug targets are also needed for those patients who are resistant to conventional treatment, including anti-TNFα therapy.
We show that our signal transduction pathway (STP) technology can be used to develop these clinically needed assays, and to identify new drug targets in patients with CD and UC.
Methods
Measurement of activity of signal transduction pathways on affymetrix microarray data from cell samples
Development and validation of computational models to quantitatively measure activity of the androgen receptor (AR), MAPK, NFκB, WNT, TGFβ, JAK-STAT1/2, and JAK-STAT3 pathways have been described before and is simplified shown in Supplementary Figure S1. The models can in principle be used to analyze pathway activity on any cell type, based on analysis of transcriptome data (Bouwman et al., 2020), (Verhaegh et al., 2014), (Stolpe et al., 2019), (Canté-Barrett et al., 2020), (Wesseling-Rozendaal et al., 2021), (Van Ooijen et al., 2018), (Bouwman et al., 2020), (Wesseling-Rozendaal et al., 2022), (Bouwman et al., 2021). For the current study, pathway activity scores (PAS) were measured on publicly available Affymetrix expression microarray data from preclinical and clinical IBD studies.
PAS were calculated from Affymetrix expression microarray data obtained from the GEO database (NCBI, 2021). For each signaling pathway, PAS represent the log2 of the odds of the pathway-associated transcription factor being active vs. inactive as calculated by a probabilistic computational model, as described before (Verhaegh et al., 2014), (Van Ooijen et al., 2018), (Stolpe et al., 2019).
Microarray data quality control
Quality control (QC) was performed on Affymetrix data of each individual sample, using standard methods as described before (Stolpe et al., 2019). In summary, QC parameters include: the average value of all probe intensities, presence of negative or extremely high (>16-bit) intensity values, poly-A RNA (sample preparation spike-ins) and labelled cRNA (hybridization spike ins) controls, GAPDH, and ACTB 3′/5′ ratio, center of intensity and values of positive and negative border controls determined by affyQCReport package in R, and an RNA degradation value determined by the AffyRNAdeg function from the Affymetrix package in R (Parman et al., 2018), (Gautier et al., 2004). Sample data that failed QC were removed prior to data analysis.
Brief description of the analyzed preclinical and clinical studies
GSE16879 (Arijs et al., 2009). Mucosal biopsies had been obtained at endoscopy from actively inflamed mucosa from IBD patients [24 ulcerative colitis (UC), 19 Crohn’s colitis (CDc) and 18 Crohn’s ileitis (CDi)], refractory to corticosteroids and/or immunosuppression, before and 4–6 weeks after intravenous infliximab (5 mg/kg) infusion, and from normal mucosa of 12 control patients. Biopsies had been taken at sites of active inflammation; in case of healing the repeat biopsy had been obtained in the area where lesions had been present prior to anti-TNFα therapy. Clinical response had been assessed 4–6 weeks after starting treatment. Of note, response assessment criteria used for CDc had been adapted to define (partial) response in CDi.
GSE23597 (Toedter et al., 2011). 48 UC patients with moderate-to-severe UC, refractory to corticosteroids, had been treated with 5 or 10 mg/kg infliximab or placebo at weeks 0, 2, 6 and every 8 weeks thereafter. Colon biopsy samples had been collected 15–20 cm from the anal verge prior to treatment and during treatment at weeks 8 and 30, that is, after induction and during maintenance therapy, respectively.
GSE52746 (Leal et al., 2015), Colon mucosa biopsy samples had been taken from CD patients and 17 non-inflammatory controls. CD patients included 10 patients with active CD without anti-TNF therapy, 5 with active CD under anti-TNF therapy (non-responders) and seven patients with inactive CD under anti-TNF therapy (responders). In patients with active CD, biopsies had been taken from the affected areas; in patients with inactive disease from the healed region. Patients had been either treated with adalimumab or infliximab. Adalimumab had been given subcutaneously (160/80 mg) at weeks 0 and 2 for induction, and (40 mg) every other week for maintenance treatment. Infliximab had been given intravenously (5 mg/kg) at weeks 0, 2, and 6 for induction, and every 8 weeks for maintenance treatment. For six patients paired biopsy data were available before and after 12 weeks of anti-TNFα treatment (5 responder patients and one non-responder).
Privacy
All analyzed data were obtained in agreement with GDPR regulations, as implemented within Philips.
Interpretation of signal transduction pathway activity score (PAS)
An important and unique advantage of the STP technology is that in principle analysis can be performed on each cell type. As described before (Stolpe et al., 2019), one should bear in mind a few considerations when interpreting log2 odds PAS:
1) On the same sample, log2 odds PAS cannot be compared between different signaling pathways, since each of the signaling pathways has its own range in log2 odds activity scores (Stolpe et al., 2019).
2) The log2 odds range for pathway activity (minimum to maximum activity) may vary depending on cell type. Once the range has been defined using samples with known pathway activity, on each new sample the absolute value can be directly interpreted against that reference. If the range has not been defined, only differences in log2 odds activity score between samples can be interpreted.
3) PAS are highly quantitative, and even small differences in log2 odds PAS can be reproducible and meaningful.
4) A negative log2 odds ratio does not necessarily mean that the pathway is inactive.
Statistics
Boxplots, individual sample plots and hierarchical clustering plots have been made using the Python data visualization library function seaborn; additional statistical annotations have been created using Python package statannot. Box-and-whisker plots show median values, interquartile range (IQR) between 25th and 75th percentile, extreme values, and outlier samples (diamond shape). Individual samples are presented as small circles (outlier samples are hence plotted twice, as diamond and as circle). Two-sided t-test (low sample numbers) or Mann-Whitney-Wilcoxon testing was used to compare PAS across groups. p-values are indicated in the figures, with p < 0.01 considered as significant. Mean log2 odds, mean delta between groups, p-values, u-stats and area under the curve (AUC) values have been summarized in Supplementary Table S1.
Results
STP activity scores from colon mucosa biopsy samples from clinical studies on CD and UC are shown in Figures 1–4 and Supplementary Figures S2, S3. In the description of results, unless otherwise specified, mentioned differences in PAS are significant, meaning p < 0.01.
FIGURE 1
FIGURE 1. Diagnosis of IBD. Datasets indicated on top of figure (A–D). STP activity analysis of intestinal mucosa samples for diagnosis of IBD. Mucosa samples from healthy control individuals and pre-treatment mucosa samples from IBD patients. For details, see Methods. Signal transduction pathway (STP) activity scores (PAS) are shown for the androgen receptor (AR), MAPK, NFκB, TGFβ (TGFbeta), JAK-STAT1/2, JAK-STAT3, and Wnt STPs. PAS on Y-axis on a log2 odds scale. Two sided Mann-Whitney-Wilcoxon statistical tests were performed; p-values are indicated in the figures as *p < 0.05, ***p < 0.001, ****p < 0.0001, ns: not significant.
FIGURE 2
FIGURE 2. Prediction and assessment of response to anti-TNFα remission-induction treatment in CD (A) and UC (B). Datasets GSE16879. STP activity analysis of intestinal mucosa samples for prediction of response to anti-TNFα remission induction treatment in patients with CD and UC. STP PAS before and after remission induction treatment. For details, see Methods. STP PAS are shown for the AR, MAPK, NFκB, TGFβ, JAK-STAT1/2, JAK-STAT3, and Wnt STPs. Two sided Mann–Whitney–Wilcoxon statistical tests were performed; p-values are indicated in the figures as *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001, ns: not significant. For supporting analysis results, see Supplementary Table S1.
FIGURE 3
FIGURE 3. Prediction and assessment of response to anti-TNFα remission-induction treatment. Dataset GSE16879. STP activity analysis of intestinal mucosa samples for prediction of response to anti-TNFα remission induction treatment in patients with CDc (A) and UC (B). STP PAS before and after remission induction treatment, for responder and non-responder patients. For details, see Methods. STP PAS are shown for the AR, MAPK, NFκB, TGFβ, JAK-STAT1/2, JAK-STAT3, and Wnt STPs. From left to right: CDc responders; CDc non-responders; UC responders; UC non-responders. Pre-treatment (left) and post-treatment (right) PAS are connected by a line for each individual patient. Paired t-tests were performed; p-values are indicated above figures. For associated barplots, see Figure 2. For supporting analysis results, see Supplementary Table S1.
FIGURE 4
FIGURE 4. Assessment of response to anti-TNFα remission-induction versus maintenance treatment. Dataset GSE23597. STP activity analysis of intestinal mucosa samples for assessment of response to low or high dose anti-TNFα remission maintenance treatment of patients with UC. Clinical response to treatment was assessed at 8 [(A), remission-induction treatment] and 30 [(B), maintenance treatment] weeks after start of treatment. Treatment with infliximab, 5 or 10 mg/kg, or placebo. Colon mucosa samples: 0, 8, and 30 weeks of treatment. For STP analysis, data of 5 and 10 mg/kg were combined. For details, see Methods. (A) STP activity scores (PAS) at 0, 8, and 30 weeks of treatment. Clinical response assessment at 8 weeks, following remission induction treatment. (B) PAS at 0, 8, and 30 weeks of treatment. Clinical response assessment at 30 weeks, following remission-induction and maintenance treatment. (C) Mean delta (delta PAS calculated per individual patient) PAS between 0 and 8 weeks, and between 0 and 30 weeks of treatment, for responders (R; red and turquoise) and non-responders (NR; green and purple). STP PAS shown for AR, MAPK, NFκB, TGFβ, JAK-STAT1/2, JAK-STAT3, and Wnt STPs. PAS on Y-axis on a log2 odds scale. Two sided Mann–Whitney–Wilcoxon statistical tests were performed; p-values are indicated in the figures as *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001, ns: not significant. For supporting analysis results, see Supplementary Table S1.
Diagnosis of CD and UC
To investigate which STPs are abnormally active in colon mucosa of IBD patients, and whether measuring STP activity can improve diagnostic accuracy, and can be used to facilitate differential diagnosis of CD and UC, STP activity was calculated for mucosa biopsy samples of IBD patients with different clinical manifestations. For CD, colon biopsy sample data were available from patients with specific involvement of either colon (CDc) or ileum (CDi) (dataset GSE16879), and of patients with disease present in colon, but not necessarily limited to colon (dataset GSE52746).
In untreated CD (dataset GSE16879, CDc and CDi combined), MAPK, NFκB, JAK-STAT1/2, JAK-STAT3, and Wnt PAS were increased compared to healthy controls (Figure 1A, left column, and Supplementary Figure S1). There were also clear differences in pathway activity between CDc and CDi (Figure 1B; Supplementary Table S1 and see Table 1 for summary): Compared to healthy mucosa, AR, NFκB, JAK-STAT3, and Wnt PAS were increased in CDc, while JAK-STAT1/2 PAS was specifically increased in CDi. This suggests different mechanisms of disease for colon- and ileum-localized CD. Analysis of the second clinical CD study (clinical disease present in, but not necessarily limited to, colon; dataset GSE52746) showed an abnormal pathway activity profile similar to the combined profile of CDc and CDi, complemented by increased TGFβ PAS (Figure 1D; Supplementary Table S1 and Table 1 for summary).
TABLE 1
TABLE 1. Abnormal STP activity in patients with various forms of IBD (diagnosis), and in responder versus non-responder patients, pre-treatment and during treatment. Diagnosis: Comparison between IBD and healthy control STP activity; indicated STPs have higher PAS in IBD group. Pre-treatment prediction of response to anti-TNFα treatment: comparison between pre-treatment STP activity of responder versus non-responder groups; indicated STPs have higher PAS in responder patient group. Individual patient response assessment during treatment: comparison between delta (for individual patients) in STP activity before - during treatment of responder versus delta in non-responder groups; indicated STPs have higher delta PAS in responder patients. Patient group response assessment during treatment: Comparison patient responder with non-responder groups during treatment; indicated STPs have higher PAS in responder group during treatment. Indicated STPs in the table: STPs with significantly higher PAS (p < 0.01); between brackets: STPs with marginal significance (p < 0.05). See Supplementary Table S1 for absolute p-values. For specifics of patient inclusion, dosing and sample type, see Methods for description of datasets. NA: Not applicable.
For untreated UC patients (dataset GSE16879) the same STP PAS were increased as found in patients with CDc, with addition of the TGFβ pathway (Figure 1C; Table 1).
Summarizing these results (Supplementary Table S1 and Table 1 for summary), activity of the MAPK pathway was reproducibly increased in colon mucosa samples of patients with all investigated forms of IBD (CD, CDc, CDi, UC), indicative for a role for this pathway in the core mechanism of disease of IBD. Furthermore, IBD restricted to the ileum (CDi) was specifically associated with increased JAK-STAT1/2 pathway activity, while colon location was associated with increased AR pathway activity. Differences in STP activity between CD, CDc, CDi, and UC provide evidence for the presence of multiple IBD subtypes with variant mechanisms of disease.
Prediction of response to anti-TNFα treatment using STP analysis of pre-treatment mucosa biopsy
To investigate whether it is possible to predict successful remission-induction by anti-TNFα treatment in IBD patients based on pre-treatment mucosa biopsy analysis, pre-treatment data from CD, CDc and CDi subtype (GSE16879; GSE527460) and UC (GSE16879, GSE23597) patients were analyzed.
In both CDc and UC patients (GSE16879), NFκB, TGFβ, and JAK-STAT3 PAS were higher in non-responder compared to responder patients (Figure 2; Supplementary Table S1 and summary Table 1). CDi patients are known to generally not respond to anti-TNFα therapy, and response criteria in this group were reported to have been adapted to allow partial responders to be included in the responder group. However, no STPs were identified of which PAS predicted partial response in CDi patients. In the CD patients with colonic/ileocolonic disease from GSE52746, patients had been treated with two types of TNFα inhibitors (subcutaneous adalimumab of intravenous infliximab) and matched samples before and after treatment were only available for six adalimumab-treated patients, of which 5 responder and one non-responder patient; thus statistical analysis was not possible (Supplementary Figure S3).
In the second UC study (GSE23597, lesion-independent biopsy) two different dosages of anti-TNFα had been administered, and data were combined for our analysis. Again, high pre-treatment NFκB, TGFβ and JAK-STAT3 PAS were associated with resistance to remission-induction treatment (but not related to failure of maintenance treatment) (Figure 4; Supplementary Table S1 and summary Table 1).
Summarizing these results, remission-induction response to the TNFα inhibitor Infliximab in patients with colon-localized CD and UC was reproducibly predicted by high pre-treatment NFκB, TGFβ, and JAK-STAT3 signa transduction pathway activity.
Response monitoring of remission-induction and maintenance anti-TNFα treatment
Comparison of STP activity profiles between post-treatment responder and non-responder groups (not per individual patient) revealed the same hyperactive STPs which had been identified as pathogenic by comparing healthy individuals with IBD patients (Figures 2–4; Supplementary Table S1 and summary Table 1, “Patient group response assessment during treatment”). Successful treatment resulted in return to (near) normal STP PAS values in the responder patient groups. Figure 4C illustrates the large differences in delta PAS induced by anti-TNFα treatment in responder patients compared to non-responder patients. These results support the pathogenic involvement of these STPs in IBD.
Clinically it is expected to be of clinical relevance to know whether a patient has a biochemical response, for example, to guide therapy decisions in case of an unclear clinical presentation. To check whether successful remission-induction can be assessed in an individual patient by comparing pre- and post-treatment STP activity profiles in sequentially taken mucosal biopsy samples, we calculated the delta STP PAS values for each individual patient (Figure 3; Supplementary Table S1; Table 1, “Individual patient response assessment during treatment”).
In individual UC patients who underwent successful remission-induction therapy (GSE16879), the decrease (delta) in AR and MAPK PAS was higher than in non-responder patients (Supplementary Table S1). In addition, successful maintenance therapy in individual UC patients (GSE23597) was also characterized by a larger decrease (delta) in MAPK PAS (Supplementary Table S1 and summary Table 1).
Comparison of STP analysis with data analysis results obtained by the investigators who performed the clinical studies and generated the datasets
We compared our STP analysis results with the data analysis results obtained by the investigators who had generated the Affymetrix data which we analyzed for the purpose of the current study. For their data analysis they used conventional public data analysis tools and compared groups of samples, aiming to identify differentially involved cellular biochemical processes. In contrast to the current analysis, their analysis was not performed on individual samples. Below we briefly summarize their results with the used data analysis tools, per dataset.
GSE16879. Data was analyzed using Bioconductor tools, with focus on differentially involved antimicrobial peptides (AMP) genes. In non-responders to infliximab, differences in AMP gene expression were not found, neither before nor during treatment.
GSE52746. Differential expression and clustering analysis was performed using Ingenuity Pathway Analysis (IPA). Top 10 biochemical processes identified were: Docosahexanoic Acid signaling; IL17A signaling in fibroblasts; tumoricidal function of hepatic natural killer cells; role of IL17F in allergic inflammatory airway disease; role of IL17A in psoriasis; granulocyte adhesion and diapedesis; differential regulation of cytokine production in intestinal epithelial cells by IL17A and IL17F; agranulocyte adhesion and diapedesis; role of Tissue Factor in cancer. Mechanistic networks: IL1B, IL17, NFkB, STAT3 are candidate upstream regulators. The investigators concluded that in non-responders an IL-17A-driven response was present which was not sensitive to anti-TNFα blockade. While involvement of NFκB and JAK-STAT3 pathways was identified, no information on activity status of these signaling pathways was obtained.
GSE23597. Data was analyzed using differential expression and clustering analysis with IPA. Top 6 biochemical processes identified were: cell-to-cell signaling and interaction; cellular movement; hematological system development and function; immune cell trafficking; inflammatory response; skeletal and muscular disorders. No information on involvement of specific signaling pathways and their activity was reported.
Summary of results
Despite many differences between the investigated studies, STP analysis results were generally consistent across patient groups, both for CD and UC. We identified the MAPK STP as a ‘core’ abnormally active STP in IBD, while NFκB, TGFβ, JAK-STAT1/2, JAK-STAT3, and possibly the Wnt pathway, seem to be differentially involved in the various clinical IBD presentations. High activities of NFκB, TGFβ and JAK-STAT3 pathways prior to treatment were predictive for resistance to anti-TNFα treatment in patients with colon-localized IBD, but not in ileal CD. Successful remission-induction treatment and maintenance was associated with a reset in pathway activity towards normal. For all STPs there was a large variation in STP PAS among patients within the various IBD subtypes.
Discussion
We previously reported a completely new method to quantitatively measure activity of clinically relevant signal transduction pathways (STP) on transcriptome data of individual cell-containing samples (Bouwman et al., 2020), (Verhaegh et al., 2014), (Stolpe et al., 2019), (Canté-Barrett et al., 2020), (Wesseling-Rozendaal et al., 2021), (Van Ooijen et al., 2018), (Bouwman et al., 2020), (Wesseling-Rozendaal et al., 2022), (Bouwman et al., 2021). Using this technology to analyze RNA expression data of colon mucosa biopsy samples, signaling pathway activity was investigated in CD and UC, and related to pathophysiology of IBD and the clinical response to anti-TNFα treatment.
Since abnormal activity of STPs is informative on the pathophysiology of a disease (the “mechanism of disease”) we investigated which STPs are abnormally active in different forms of IBD. Activity of the MAPK growth factor pathway was consistently increased in colon biopsies of patients with either UC or CD, irrespective of colonic or ileal disease location. Clinically effective anti-TNFα treatment (partly) corrected this abnormal signaling pathway activity, further supporting a core pathogenic role for the MAPK signaling pathway in IBD, independent of disease location. Indeed, the MAPK pathway has been described as an important signaling pathway in a number of related auto-immune diseases, such as IBD and rheumatoid arthritis (ten Hove et al., 2002), (Kumar et al., 2003), (Coskun et al., 2011), (Xu et al., 2019).
A few signaling pathways showed differential activity between colon- and ileum-located IBD. Activity of the NFκB, AR, JAK-STAT3, Wnt, and (less consistently) the TGFβ, signaling pathways was increased in patients with IBD disease localization in the colon, while the JAK-STAT1/2 pathway was specifically active in ileal localized CD, suggesting differences in pathophysiology depending on disease location. This is in line with current thinking regarding differences between ileal- and colon-localized IBD, including marked differences in clinical response to anti-TNFα treatment (Arijs et al., 2009). The inflammatory NFκB pathway, driving production of TNFα, is well known as a core mechanistic pathway in colon- but less in ileum-located IBD (Billmeier et al., 2016). Indeed, patients with ileum-located CD generally are known to respond less well to anti-TNFα treatment (Liu et al., 2017). Due to sequential pathway activation (e.g., an active NFκB pathway produces IL6, resulting in activation of the JAK-STAT3 pathway) and intracellular crosstalk between signaling pathways (e.g., cooperation between transcription factors of MAPK and TGFβ pathways), multiple pathways, i.e., the Wnt, MAPK, JAK-STAT and TGFβ pathways, are recruited in the pathogenic process underlying IBD (Moparthi and Koch, 2019), (Chen and Ten Dijke, 2016), (Massagué, 2012), (Salas et al., 2020), (Carey et al., 2008), (Na and Moon, 2019). While most of these pathways are thought to contribute to the inflammatory character of the disease, the TGFβ pathway may play an instrumental role in the fibrotic process associated with IBD (Yun et al., 2019).
Our results support the increasing evidence that the pathogenesis (or mechanism of disease) of ileal IBD is distinct from colon-located IBD (Pierre et al., 2020), (Naftali et al., 2016), (Kredel et al., 2019), (Howell et al., 2018), (Meinzer et al., 2005), (Wu et al., 2010), (van der Meer et al., 2020). Interferons are activators of the JAK-STAT1/2 pathway which was in our analysis identified as specifically active in ileal CD (Platanias, 2005). Indeed, interferon lambda concentrations have been reported as increased in mucosa of IBD-affected ileum, both in mouse models and in patients, supporting our findings (Günther et al., 2019). We believe that STP analysis contributes to our understanding of the mechanism of disease underlying the various manifestations of IBD and may improve diagnosis and differential diagnosis between colon- and ileum-located disease.
In patients with a diagnosis of IBD, either CD or UC, a biomarker assay to predict response to anti-TNFα treatment is needed because of potentially serious side effects of the treatment and significant treatment costs (D’Haens and Deventer, 2021). In our study, in IBD patients with colon involvement, high activity of NFκB, JAK-STAT3 and TGFβ pathways prior to treatment was associated with failure to induce remission by anti-TNF-α treatment. The NFκB pathway is the prime target for anti-TNFα treatment, since activity of this pathway is both induced by TNFα and responsible for production of TNFα (Liu et al., 2017), (Billmeier et al., 2016). The known secondary activation of JAK-STAT3 and TGFβ pathways by a (TNFα-) activated NFκB pathway readily explains their combined association with resistance to anti-TNFα treatment (Moparthi and Koch, 2019), (Chen and Ten Dijke, 2016), (Massagué, 2012).
NFκB pathway activity seemed to play a less dominant role in ileal-CD, clearly in line with anti-TNFα treatment being less effective in patients with ileal-restricted CD (Liu et al., 2017). Indeed, in the here analyzed clinical study, clinical response to anti-TNFα treatment had been separately defined on less strict criteria than used for clinical response in the other IBD groups (Arijs et al., 2009). Together our results provide initial evidence that STP analysis can be used to predict response to anti-TNFα treatment, except in patients with ileum-restricted CD.
While considered highly necessary from a health economic standpoint as well as from a clinical point of view, currently clinical response to anti-TNFα treatment cannot be predicted in IBD patients (van der Valk et al., 2014; Van der Valk et al., 2016), (Lopetuso et al., 2017), (D’Haens and Deventer, 2021), (Targownik et al., 2019). Literature data report that 20–40% of patients in clinical trials and 10–20% of patients in “real life” series are primary resistant to anti-TNFα induction therapy, while secondary resistance after a year treatment is even more common (Ben-Horin et al., 2014).
Our results suggest that to overcome resistance to anti-TNFα treatment, it may be needed to target the NFκB signal transduction pathway in a different manner than by inhibiting its ligand TNFα (van de Stolpe et al., 1994), (Miller et al., 2010). Indeed, the signaling pathways here identified as abnormally active in IBD patients provide new drug targets and opportunities for alternative treatment strategies. Development of MAPK pathway inhibitors for IBD has been explored without clinical success, probably due to the complex and not well understood role of this pathway in IBD (Kumar et al., 2003), (ten Hove et al., 2002), (Xu et al., 2019), (Coskun et al., 2011). JAK-STAT pathway inhibitors are being explored for IBD and the JAK-STAT3 pathway inhibitor tofacitinib has been approved for treatment of UC (Salas et al., 2020), (Ducreux et al., 2014), (Bouwman et al., 2020), (Akhurst and Hata, 2012), (Huynh et al., 2019), (Mitra et al., 2020). TGFβ inhibitors are predominantly being developed for cancer immunotherapy, but reducing intestinal fibrosis in IBD has been suggested as a potential clinical application (Huynh et al., 2019). For ileum-localized CD, based on our current results we hypothesize that JAK-STAT1/2 inhibitors may be more effective than TNFα inhibitors (Rogler, 2020), (Salas et al., 2020).
Potentially relevant for any treatment of IBD is the observation that activity of the signaling pathways showed large variations among patients with IBD. This suggests that clinical efficacy of a drug, including an anti-TNFα drug, in an individual patient may depend on the level of activity of the targeted signaling pathway, implicating that a personalized treatment approach is likely to be necessary to optimize treatment results. This requires an assay to measure STP activity in a sample of every individual patient to predict response to a drug. For example, a decision on treatment with an anti-TNFα drug such as infliximab ideally would be preceded by measuring NFκB pathway activity to guide this treatment decision.
For each Affymetrix dataset for which we performed STP analysis, data-analysis results obtained by the original investigators were compared to our results. While gene expression profiles had been obtained for CD and UC patient groups (instead of individual patients), no gene profile had been identified by the previous investigators to separate anti-TNFα responder and non-responder group, and no information on STP activity had been extracted from the Affymetrix data. This emphasizes the difference between our STP activity analysis, which can be performed on individual patient samples and provides clinically actionable STP information for an individual patient (e.g., responder or non-responder to infliximab), and other bioinformatics tools which are used to discover gene expression profiles as potential novel biomarkers for groups of patients.
In conclusion, STP analysis provided a first-of-a-kind insight into the mechanism of disease of IBD and its relation to treatment resistance. Our results provide new treatment targets, for example, the MAPK pathway in all types of IBD, in ileal Crohn’s disease the JAK-STAT1/2 pathway, and in patients refractory to anti-TNFα treatment the NFκB, TGFβ, and JAK-STAT3 pathways.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
WB: Data, results analysis, concept and writing. WV: Pathway model development, results analysis. AS: Pathway model development, concept, results analysis, and writing. All authors contributed to the article and approved the submitted version.
Acknowledgments
We wish to acknowledge all investigators who performed the clinical studies and generated the publicly available GEO datasets used in the current analysis.
Conflict of interest
Authors WB, WV, and AS were employed by the company Philips Research.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.1008976/full#supplementary-material
References
Akhurst, R. J., and Hata, A. (2012). Targeting the TGFβ signalling pathway in disease. Nat. Rev. Drug Discov. 11 (10), 790–811. doi:10.1038/nrd3810
Arijs, I., De Hertogh, G., Lemaire, K., Quintens, R., Van Lommel, L., Van Steen, K., et al. (2009). Mucosal gene expression of antimicrobial peptides in inflammatory bowel disease before and after first infliximab treatment. PLoS One 4 (11), e7984. doi:10.1371/journal.pone.0007984
Ben-Horin, S., Kopylov, U., and Chowers, Y. (2014). Optimizing anti-TNF treatments in inflammatory bowel disease. Autoimmun. Rev. 13 (1), 24–30. doi:10.1016/j.autrev.2013.06.002
Billmeier, U., Dieterich, W., Neurath, M. F., and Atreya, R. (2016). Molecular mechanism of action of anti-tumor necrosis factor antibodies in inflammatory bowel diseases. World J. Gastroenterol. 22 (42), 9300–9313. doi:10.3748/wjg.v22.i42.9300
Bouwman, W., Verhaegh, W., Holtzer, L, and van de Stolpe, A. (2020). Measurement of cellular immune response to viral infection and vaccination. Front. Immunol. 11, 575074. doi:10.3389/fimmu.2020.575074
Bouwman, W., Verhaegh, W., and Stolpe, A. Androgen receptor pathway activity assay for sepsis diagnosis and prediction of favorable prognosis. medRxiv. 2021.
Canté-Barrett, K., Holtzer, L., van Ooijen, H., Hagelaar, R., Cordo’, V., Verhaegh, W., et al. (2020). A molecular test for quantifying functional Notch signaling pathway activity in human cancer. Cancers (Basel) 12, E3142. doi:10.3390/cancers12113142
Carey, R., Jurickova, I., Ballard, E., Bonkowski, E., Han, X., Xu, H., et al. (2008). Activation of an IL-6:STAT3-dependent transcriptome in pediatric-onset inflammatory bowel disease. Inflamm. Bowel Dis. 14 (4), 446–457. doi:10.1002/ibd.20342
Chen, W., and Ten Dijke, P. (2016). Immunoregulation by members of the TGFβ superfamily. Nat. Rev. Immunol. 16 (12), 723–740. doi:10.1038/nri.2016.112
Colombel, J. F., Sandborn, W. J., Reinisch, W., Mantzaris, G. J., Kornbluth, A., Rachmilewitz, D., et al. (2010). Infliximab, azathioprine, or combination therapy for Crohn’s disease. N. Engl. J. Med. 362 (15), 1383–1395. doi:10.1056/NEJMoa0904492
Coskun, M., Olsen, J., Seidelin, J. B., and Nielsen, O. H. (2011). MAP kinases in inflammatory bowel disease. Clin. Chim. Acta. 412 (7–8), 513–520. doi:10.1016/j.cca.2010.12.020
D’Haens, G., Baert, F., Assche, G., Caenepeel, P., Vergauwe, P., Tuynman, H., et al. (2008). Early combined immunosuppression or conventional management in patients with newly diagnosed Crohn’s disease: An open randomised trial. Lancet 371 (9613), 660–667. doi:10.1016/S0140-6736(08)60304-9
D’Haens, G. R., and Deventer, S. (2021). 25 years of anti-TNF treatment for inflammatory bowel disease: Lessons from the past and a look to the future. Gut 70 (7), 1396–1405. doi:10.1136/gutjnl-2019-320022
D’Haens, G. R. (2010). Top-down therapy for IBD: Rationale and requisite evidence. Nat. Rev. Gastroenterol. Hepatol. 7 (2), 86–92. doi:10.1038/nrgastro.2009.222
Ducreux, J., Durez, P., Galant, C., Nzeusseu Toukap, A., Van den Eynde, B., Houssiau, F. A., et al. (2014). Global molecular effects of tocilizumab therapy in rheumatoid arthritis synovium. Arthritis Rheumatol. 66 (1), 15–23. doi:10.1002/art.38202
Feuerstein, J. D., Ho, E. Y., Shmidt, E., Singh, H., Falck-Ytter, Y., Sultan, S., et al. (2021). AGA clinical practice guidelines on the medical management of moderate to severe luminal and perianal fistulizing Crohn’s disease. Gastroenterology 160 (7), 2496–2508. doi:10.1053/j.gastro.2021.04.022
Gautier, L., Cope, L., Bolstad, B. M., and Irizarry, R. A. (2004). affyanalysis of Affymetrix GeneChip data at the probe level. Bioinformatics 20 (3), 307–315. doi:10.1093/bioinformatics/btg405
Günther, C., Ruder, B., Stolzer, I., Dorner, H., He, G. W., Chiriac, M. T., et al. (2019). Interferon lambda promotes paneth cell death via STAT1 signaling in mice and is increased in inflamed ileal tissues of patients with Crohn’s disease. Gastroenterology 157 (5), 1310–1322. doi:10.1053/j.gastro.2019.07.031
Hazel, K., and O’Connor, A. (2020). Emerging treatments for inflammatory bowel disease. Ther. Adv. Chronic Dis. 11, 2040622319899297. doi:10.1177/2040622319899297
Howell, K. J., Kraiczy, J., Nayak, K. M., Gasparetto, M., Ross, A., Lee, C., et al. (2018). DNA methylation and transcription patterns in intestinal epithelial cells from pediatric patients with inflammatory bowel diseases differentiate disease subtypes and associate with outcome. Gastroenterology 154 (3), 585–598. doi:10.1053/j.gastro.2017.10.007
Huynh, L. K., Hipolito, C. J., and Ten Dijke, P. (2019). A perspective on the development of TGF-β inhibitors for cancer treatment. Biomolecules 9 (11), E743. doi:10.3390/biom9110743
Kredel, L. I., Jödicke, L. J., Scheffold, A., Gröne, J., Glauben, R., Erben, U., et al. (2019). T-Cell composition in ileal and colonic creeping fat - separating ileal from colonic Crohn’s disease. J. Crohns Colitis 13 (1), 79–91. doi:10.1093/ecco-jcc/jjy146
Kumar, S., Boehm, J., and Lee, J. C. (2003). p38 MAP kinases: key signalling molecules as therapeutic targets for inflammatory diseases. Nat. Rev. Drug Discov. 2 (9), 717–726. doi:10.1038/nrd1177
Leal, R. F., Planell, N., Kajekar, R., Lozano, J. J., Ordás, I., Dotti, I., et al. (2015). Identification of inflammatory mediators in patients with Crohn’s disease unresponsive to anti-TNFα therapy. Gut 64 (2), 233–242. doi:10.1136/gutjnl-2013-306518
Liu, T., Zhang, L., Joo, D., and Sun, S. C. (2017). NF-κB signaling in inflammation. Signal Transduct. Target. Ther. 2 (1), 17023–17029. doi:10.1038/sigtrans.2017.23
Lopetuso, L. R., Gerardi, V., Papa, V., Scaldaferri, F., Rapaccini, G. L., Gasbarrini, A., et al. (2017). Can we predict the efficacy of anti-TNF-α agents? Int. J. Mol. Sci. 18 (9), 1973. doi:10.3390/ijms18091973
Massagué, J. (2012). TGFβ signalling in context. Nat. Rev. Mol. Cell Biol. 13 (10), 616–630. doi:10.1038/nrm3434
Meinzer, U., Ideström, M., Alberti, C., Peuchmaur, M., Belarbi, N., Bellaïche, M., et al. (2005). Ileal involvement is age dependent in pediatric Crohn’s disease. Inflamm. Bowel Dis. 11 (7), 639–644. doi:10.1097/01.mib.0000165114.10687.bf
Miller, S. C., Huang, R., Sakamuru, S., Shukla, S. J., Attene-Ramos, M. S., Shinn, P., et al. (2010). Identification of known drugs that act as inhibitors of NF-kappaB signaling and their mechanism of action. Biochem. Pharmacol. 79 (9), 1272–1280. doi:10.1016/j.bcp.2009.12.021
Mitra, M. S., Lancaster, K., Adedeji, A. O., Palanisamy, G. S., Dave, R. A., Zhong, F., et al. (2020). A potent pan-tgfβ neutralizing monoclonal antibody elicits cardiovascular toxicity in mice and cynomolgus monkeys. Toxicol. Sci. 175 (1), 24–34. doi:10.1093/toxsci/kfaa024
Moparthi, L., and Koch, S. (2019). Wnt signaling in intestinal inflammation. Differentiation. 108, 24–32. doi:10.1016/j.diff.2019.01.002
Na, S. Y., and Moon, W. (2019). Perspectives on current and novel treatments for inflammatory bowel disease. Gut Liver 13 (6), 604–616. doi:10.5009/gnl19019
Naftali, T., Reshef, L., Kovacs, A., Porat, R., Amir, I., Konikoff, F. M., et al. (2016). Distinct microbiotas are associated with ileum-restricted and colon-involving Crohn’s disease. Inflamm. Bowel Dis. 22 (2), 293–302. doi:10.1097/MIB.0000000000000662
NCBI. Home GEO DataSets NCBI. 2021, Available at: https://www.ncbi.nlm.nih.gov/gds/. [Accessed 2021 Feb 24].
Parman, C., Halling, C., Gentleman, R., Parman, C., Halling, C., and Gentleman, R. Data from: AffyQCReport: QC report generation for affyBatch objects. 2018. Version 1.68.0
Pierre, N., Salée, C., Massot, C., Blétard, N., Mazzucchelli, G., Smargiasso, N., et al. (2020). Proteomics highlights common and distinct pathophysiological processes associated with ileal and colonic ulcers in Crohn’s disease. J. Crohns Colitis 14 (2), 205–215. doi:10.1093/ecco-jcc/jjz130
Platanias, L. C. (2005). Mechanisms of type-I- and type-II-interferon-mediated signalling. Nat. Rev. Immunol. 5 (5), 375–386. doi:10.1038/nri1604
Rogler, G. (2020). Efficacy of JAK inhibitors in Crohn’s disease. J. Crohns Colitis 14, S746–S754. doi:10.1093/ecco-jcc/jjz186
Salas, A., Hernandez-Rocha, C., Duijvestein, M., Faubion, W., McGovern, D., Vermeire, S., et al. (2020). JAK-STAT pathway targeting for the treatment of inflammatory bowel disease. Nat. Rev. Gastroenterol. Hepatol. 17 (6), 323–337. doi:10.1038/s41575-020-0273-0
Stolpe, A., Holtzer, L., van Ooijen, H., de Inda, M. A., and Verhaegh, W. (2019). Enabling precision medicine by unravelling disease pathophysiology: Quantifying signal transduction pathway activity across cell and tissue types. Sci. Rep. 9 (1), 1603. doi:10.1038/s41598-018-38179-x
Targownik, L. E., Benchimol, E. I., Witt, J., Bernstein, C. N., Singh, H., Lix, L., et al. (2019). The effect of initiation of anti-TNF therapy on the subsequent direct health care costs of inflammatory bowel disease. Inflamm. Bowel Dis. 25 (10), 1718–1728. doi:10.1093/ibd/izz063
ten Hove, T., van den Blink, B., Pronk, I., Drillenburg, P., Peppelenbosch, M. P., and van Deventer, S. J. H. (2002). Dichotomal role of inhibition of p38 MAPK with SB 203580 in experimental colitis. Gut 50 (4), 507–512. doi:10.1136/gut.50.4.507
Toedter, G., Li, K., Marano, C., Ma, K., Sague, S., Huang, C. C., et al. (2011). Gene expression profiling and response signatures associated with differential responses to infliximab treatment in ulcerative colitis. Am. J. Gastroenterol. 106 (7), 1272–1280. doi:10.1038/ajg.2011.83
Van der Valk, M. E., Mangen, M. J. J., Severs, M., van der Have, M., Dijkstra, G., van Bodegraven, A. A., et al. (2016). Evolution of costs of inflammatory bowel disease over two years of follow-up. PLoS One 11 (4), e0142481. doi:10.1371/journal.pone.0142481
van de Stolpe, A., Caldenhoven, E., Stade, B. G., Koenderman, L., Raaijmakers, J. A., Johnson, J. P., et al. (1994). 12-O-tetradecanoylphorbol-13-acetate- and tumor necrosis factor alpha-mediated induction of intercellular adhesion molecule-1 is inhibited by dexamethasone. Functional analysis of the human intercellular adhesion molecular-1 promoter. J. Biol. Chem. 269 (8), 6185–6192. doi:10.1016/s0021-9258(17)37586-5
van der Meer, J. H. M., Wasmann, K. A. T. G. M., van der Bilt, J. D. W., Becker, M. A. J., Boermeester, M. A., Bemelman, W. A., et al. (2020). Anatomical variation in mesenteric macrophage phenotypes in Crohn’s disease. Clin. Transl. Gastroenterol. 11 (8), e00198. doi:10.14309/ctg.0000000000000198
van der Valk, M. E., Mangen, M., Leenders, M., Dijkstra, G., van Bodegraven, A. A., Fidder, H. H., et al. (2014). Healthcare costs of inflammatory bowel disease have shifted from hospitalisation and surgery towards anti-tnfα therapy: Results from the COIN study. Gut 63 (1), 72–79. doi:10.1136/gutjnl-2012-303376
Van Ooijen, H., Hornsveld, M., Dam-de Veen, C., Velter, R., Dou, M., Verhaegh, W., et al. (2018). Assessment of functional phosphatidylinositol 3-kinase pathway activity in cancer tissue using forkhead box-O target gene expression in a knowledge-based computational model. Am. J. Pathol. 188 (9), 1956–1972. doi:10.1016/j.ajpath.2018.05.020
Verhaegh, W., van Ooijen, H., Inda, M. A., Hatzis, P., Versteeg, R., Smid, M., et al. (2014). Selection of personalized patient therapy through the use of knowledge-based computational models that identify tumor-driving signal transduction pathways. Cancer Res. 74 (11), 2936–2945. doi:10.1158/0008-5472.CAN-13-2515
Wesseling-Rozendaal, Y., Holtzer, L., Verhaegh, W., and Stolpe, A. Measurement of activity of developmental signal transduction pathways to quantify stem cell pluripotency and phenotypically characterize differentiated cells. bioRxiv. 2021.
Wesseling-Rozendaal, Y., Van Doorn, A., Willard-Gallo, K., and Van de Stolpe, A. (2022). Characterization of immunoactive and immunotolerant CD4+ T cells in breast cancer by measuring activity of signaling pathways that determine immune cell function. Cancer Biol. 14 (3), 490. doi:10.3390/cancers14030490
Wu, F., Zhang, S., Dassopoulos, T., Harris, M. L., Bayless, T. M., Meltzer, S. J., et al. (2010). Identification of microRNAs associated with ileal and colonic Crohn’s disease. Inflamm. Bowel Dis. 16 (10), 1729–1738. doi:10.1002/ibd.21267
Xu, P., Elamin, E., Elizalde, M., Bours, P. P. H. A., Pierik, M. J., Masclee, A. A. M., et al. (2019). Modulation of intestinal epithelial permeability by plasma from patients with Crohn’s disease in a three-dimensional cell culture model. Sci. Rep. 9 (1), 2030. doi:10.1038/s41598-018-38322-8
Keywords: ulcerative colitis, crohn’s disease, IBD, signal transduction pathway, treatment response
Citation: Bouwman W, Verhaegh W and van de Stolpe A (2022) Improved diagnosis of inflammatory bowel disease and prediction and monitoring of response to anti-TNF alpha treatment based on measurement of signal transduction pathway activity. Front. Pharmacol. 13:1008976. doi: 10.3389/fphar.2022.1008976
Received: 01 August 2022; Accepted: 22 September 2022;
Published: 13 October 2022.
Edited by:
Kyle L. Flannigan, University of Calgary, CanadaReviewed by:
Amosy E. M’Koma, Meharry Medical College, United StatesSamira Salihovic, Örebro University, Sweden
Copyright © 2022 Bouwman, Verhaegh and van de Stolpe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wilbert Bouwman, V2lsYmVydC5Cb3V3bWFuQHBoaWxpcHMuY29t