Next Article in Journal
A Review of Biomimetic and Biodegradable Magnetic Scaffolds for Bone Tissue Engineering and Oncology
Next Article in Special Issue
Contribution of Blood Vessel Activation, Remodeling and Barrier Function to Inflammatory Bowel Diseases
Previous Article in Journal
circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain
Previous Article in Special Issue
Natural NADH and FAD Autofluorescence as Label-Free Biomarkers for Discriminating Subtypes and Functional States of Immune Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models

by
Andrey V. Markov
*,
Innokenty A. Savin
,
Marina A. Zenkova
and
Aleksandra V. Sen’kova
Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences, Lavrent’ev Ave., 8, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4311; https://doi.org/10.3390/ijms24054311
Submission received: 3 February 2023 / Revised: 18 February 2023 / Accepted: 20 February 2023 / Published: 21 February 2023
(This article belongs to the Special Issue Solving the Puzzle: Molecular Research in Inflammatory Bowel Diseases)

Abstract

:
Inflammatory bowel disease (IBD) is a complex and multifactorial systemic disorder of the gastrointestinal tract and is strongly associated with the development of colorectal cancer. Despite extensive studies of IBD pathogenesis, the molecular mechanism of colitis-driven tumorigenesis is not yet fully understood. In the current animal-based study, we report a comprehensive bioinformatics analysis of multiple transcriptomics datasets from the colon tissue of mice with acute colitis and colitis-associated cancer (CAC). We performed intersection of differentially expressed genes (DEGs), their functional annotation, reconstruction, and topology analysis of gene association networks, which, when combined with the text mining approach, revealed that a set of key overexpressed genes involved in the regulation of colitis (C3, Tyrobp, Mmp3, Mmp9, Timp1) and CAC (Timp1, Adam8, Mmp7, Mmp13) occupied hub positions within explored colitis- and CAC-related regulomes. Further validation of obtained data in murine models of dextran sulfate sodium (DSS)-induced colitis and azoxymethane/DSS-stimulated CAC fully confirmed the association of revealed hub genes with inflammatory and malignant lesions of colon tissue and demonstrated that genes encoding matrix metalloproteinases (acute colitis: Mmp3, Mmp9; CAC: Mmp7, Mmp13) can be used as a novel prognostic signature for colorectal neoplasia in IBD. Finally, using publicly available transcriptomics data, translational bridge interconnecting of listed colitis/CAC-associated core genes with the pathogenesis of ulcerative colitis, Crohn’s disease, and colorectal cancer in humans was identified. Taken together, a set of key genes playing a core function in colon inflammation and CAC was revealed, which can serve both as promising molecular markers and therapeutic targets to control IBD and IBD-associated colorectal neoplasia.

1. Introduction

Colorectal cancer (CRC) is the third most common malignancy and the second leading cause of cancer-related deaths worldwide [1,2]. Colon inflammation, along with the particular host and environmental factors, plays a crucial role in the initiation and progression of CRC [3]. Colitis-associated cancer (CAC) is a type of CRC, which is preceded by clinically detectable inflammatory bowel disease (IBD), including Crohn’s disease (CD) and ulcerative colitis (UC), two highly heterogeneous, incurable, persistent, relapsing/worsening, and immune-arbitrated inflammatory pathologies of the digestive system [4,5]. Epidemiologic studies have showed that patients with IBD have a predisposition to CRC, and cancer risk is highly correlated with the duration and severity of colon inflammation [6,7].
In IBD, chronic long-term colon inflammation accompanied by oxidative stress can alter the expression patterns of key carcinogenesis-associated genes [8]. Moreover, persistent stimulation of epithelial proliferation in the colon by the pro-inflammatory stimuli and excessive cell damage with increased epithelial cell turnover result in detrimental genetic and immunological alterations, making patients with IBD prone to developing CRC [9]. Despite the proven involvement of “inflammation-dysplasia-carcinoma” axis in the malignant transformation of cells in IBD-related CRC [10], the molecular mechanism underlying this process is not yet fully understood. In particular, it remains rather unclear which core genes are involved in the regulation of acute colitis and how markedly their profiles change during colitis-associated malignant transformation of the colon tissue. In addition, the proven complexity of the colitis/CAC-related regulome underlies the low efficacy of conventional IBD/CRC therapy, making it inevitable that surgery is recommended for treating these pathologies [5,11]. Given the known adverse impact of the surgical management of colonic diseases on the quality of life, mental health, and work productivity of patients [5,11], the search for novel key genes involved in the inflammation-related tumor transformation, which can be used as potential molecular targets for IBD therapy, is urgently needed. Moreover, such regulatory genes can be considered as biomarkers of inflammation-driven tumorigenesis and serve as predictors for surveillance strategies and chemoprevention of colitis-related dysplasia and CRC in IBD patients.
To date, extensive exploration of colitis- and CAC-associated regulomes has been performed using transcriptomics-based approaches [12,13,14,15,16,17,18,19,20,21]. Reported bioinformatics studies have revealed some candidate biomarker genes and key signaling pathways susceptible to the development of the mentioned disorders [14,15,16,17,18,19,20,21,22,23], colitis-induced changes in the landscape of immune infiltration of colon tissue [14,16], and a range of hub genes probably involved in the development of CAC [15,20,23]. Despite a plethora of published studies, obtained results are still uncertain and are not well correlated with each other, probably due to insufficient usage of a multiple microarray analysis algorithm (the exploration of three or more independent microarray datasets in the same study, which gives more valid results [18,20,21]), ineffective manual searching of the published literature on the topic of study [14,15,16,17,18,19,20,21], and, in some cases, the absence of proper experimental validation [18]. Since the obtained data still remain insufficient for a thorough understanding of colitis/CAC-associated gene signature, further comprehensive bioinformatics analysis of colitis/CAC-related core genes is required.
In this study, deep re-analysis of multiple microarray datasets related to murine acute colitis (GSE42768, GSE35609, GSE64658, GSE71920, GSE35609) and CAC (GSE31106, GSE5605, GSE64658, GSE42768) was performed. Firstly, the differentially expressed genes (DEGs) were computed between injured and healthy colon tissues, followed by their functional annotation and Venn diagram analysis to identify acute colitis- and CAC-associated core genes. Next, the changes in the sets of core genes associated with the transition from colon inflammation to CRC were identified. Further reconstruction and analysis of gene association networks revealed a range of hub regulators among core genes, subsequent exploration of which by the text mining approach identified a list of candidate genes, which can be used as novel promising biomarkers and therapeutic targets for colitis and CAC. The obtained results were finally validated using an in vivo model of dextran sulfate sodium (DSS)-induced acute colitis and azoxymethane (AOM)/DSS-induced CAC. Furthermore, the role of identified core genes in the colonic carcinogenesis in the backstage of chronic long-term inflammation was analyzed with respect to IBD and CRC in humans.

2. Results

2.1. Identification of Core Genes Related to Colitis and Colitis-Associated Cancer

To reveal key genes involved in the regulation of acute colitis and its transformation to CAC in mice, a range of independent expression profiles of murine colon tissue were retrieved from the GEO database, including samples of mice of both sexes and different strains with acute colitis stimulated by DSS (GSE42768, GSE35609, GSE64658, GSE71920) or dinitrobenzene sulfonic acid (DNBS) (GSE35609), or chronic colitis driven by azoxymethane (AOM)/DSS accompanied by the development of colorectal cancer (GSE31106, GSE5605, GSE64658, GSE42768). The analysis of selected transcriptomic datasets using the GEO2R tool revealed the sets of differentially expressed genes (DEGs) (colitis vs. control and CAC vs. control) susceptible to the mentioned pathologies, further overlapping of which identified 54 and 109 common DEGs specific to colitis and CAC, respectively (hereafter referred to as core genes) (Figure 1A).

2.1.1. Hierarchical Clustering and Functional Analysis of DEGs

Hierarchical clustering of the expression profiles of identified colitis-associated core genes revealed two main clades separating up- and down-regulated DEGs from each other (Figure 1B). The sub-clade of the most overexpressed DEGs included genes related to immune response (Ccl3, S100a9, S100a8, Cxcl2) and heme metabolism (Hp), whereas the most suppressed core genes in the colitis group were Hao2 and Slc26a3, associated with fatty acid metabolism and chloride ion transport, respectively (Figure 1B). Further functional analysis of colitis-specific core genes revealed high enrichment of inflammatory-related terms, including the production of pro-inflammatory cytokines IL-1 and TNF-α, IL-17, IGF1-Akt and Tyrobp signaling pathways, antiviral response, matrix metalloproteinases (MMPs), lung fibrosis, and rheumatoid arthritis (Figure 1C, upper panel).
Hierarchical clustering of CAC-specific core genes (Figure 1D) revealed two main clades, grouping activated and suppressed DEGs separately, and one outgroup consisted of the most overexpressed CAC-associated DEGs, notably, regulators of host-microbiota interplay (Reg3b, Reg3g), immune response (S100a9), and extracellular matrix (ECM) remodeling (Mmp7). In turn, the most suppressed core genes in the CAC group were involved in the regulation of cell adhesion (Zan), pH homeostasis (Car4), and ion transport (Slc26a3, Slc37a2, Aqp8) (Figure 1D). Performed gene set enrichment analysis revealed that CAC-specific core genes are tightly associated with cell invasiveness (wound healing involved in inflammatory response and MMPs), immune response (acute inflammatory response, antimicrobial peptides, etc.), redox imbalance, ion transport, bile secretion, and numerous metabolic processes (Figure 1C, lower panel). Interestingly, the retrieved functional annotation map specific for CAC was significantly less interconnected compared with the acute colitis-associated GO term/pathways network (Figure 1C, upper panel), which can be explained by the more discrete disposition of identified core genes in the CAC-related regulome.

2.1.2. Analysis of Interconnection between Acute Colitis- and CAC-Specific Core Genes

To explore how strongly identified core genes are interconnected in acute and chronic (CAC) phases of colitis, their Venn diagram analysis and the reconstruction of the gene association network were performed. Overlapping of acute colitis- and CAC-related genes demonstrated that 22 of the core genes, playing a regulatory role in acute inflammation, were involved in CAC pathogenesis (Figure 1E), including immune genes (Ifitm1, Ifitm3, Il1a, Lcn2, S100a9, Saa3, Tnf), genes encoding protease inhibitors (Serpina3n, Slpi, Wfdc18), ion transporters (Slc26a2, Slc26a3, Trpm6), ECM remodeling proteins (Mmp10, Timp1, Mep1a), signal transduction components (Igfbp4, Lrg1) and regulators of cell motility (Capg), fatty acid homeostasis (Hao2), host-microbiota interplay (Sult1a1), and heme metabolism (Hp).
Analysis of the gene association network generated from acute colitis- and CAC-associated core genes using the STRING database [24] demonstrated their relatively high interconnection: 72 of 141 uploaded core genes (51%) formed interactions with each other within the network (Figure 1F). Interestingly, only 34 of 87 CAC-specific genes (39%) were involved in the network, whereas the shares of acute colitis-specific and common genes in the reconstructed interactome were 66% (21 of 32 genes) and 73% (16 of 22 genes), respectively. Considering that highly interconnected genes can be involved in the same or similar biological processes [25], revealed low enrichment of the analyzed network by CAC-specific genes (Figure 1F) was in line with the discrete structure of the CAC-related functional annotation map shown above (Figure 1C).
Further computing of degree centrality scores of explored core genes revealed a range of genes occupying hub positions in the analyzed network (Figure 1F). It was found that the most interconnected nodes were acute colitis-specific or common genes involved in immune response (C3, Cxcl2, Il1b, Tnf) and ECM remodeling (Mmp9, Timp1). Among CAC-specific genes, the highest degree was identified for stabilizer of endoplasmic reticulum structure Ckap4, regulator of cell–cell interaction Cd44, and gene Lyz1 encoding lysozyme (Figure 1F). Given the hub position of Mmp9 and Timp1 and the formation of a highly connected cluster of MMPs in the core gene-retrieved network (Figure 1F), the changes in the MMPs profile can be involved in the regulation of malignant transformation of colon tissue during chronic colitis. This pattern needs further clarification.

2.2. Identification of Novel Acute Colitis- and CAC-Specific Hub Genes

2.2.1. Acute Colitis-Associated Hub Genes

To identify novel candidate genes for acute colitis and CAC, which can be used as both diagnostic markers and promising therapeutic targets, next we questioned how strongly evaluated core genes can be involved in the regulation of the mentioned pathologies and how well these genes have been studied in the field of inflammatory and neoplastic disorders of the colon.
To address the first issue, the degree centrality scores of the core genes in gene association networks created for each analyzed transcriptomic dataset were computed. Given that hub genes can exert key regulatory functions in reconstructed gene networks [26], the top 20 acute colitis-specific hub genes were identified and are shown in Figure 2A. The obtained results demonstrated that the most interconnected genes associated with acute colitis included genes encoding cytokines (Tnf, Il1a, Il1b), chemokines and its receptors (Ccl2, Cxcl2, Ccl3, Ccr5), growth factors and signal transduction components (Igf1, Tyrobp, Arrb2), ECM remodeling regulators (Mmp3, Mmp9, Timp1), and immune (C3, Clec7a, H2-Aa, Sell, Selp) and protective (Hp, Ugt2b35) proteins.
Next, to select genes poorly characterized for their role in colitis and colitis-associated disorders, a text mining approach was performed. Analysis of the mention of acute colitis-related core genes (Figure 1B) alongside the keywords “Colitis”, “Crohn’s”, “Dysplasia”, and “Colon cancer” in scientific texts deposited in the MEDLINE database revealed the most studied genes in the field of colitis (Tnf, Il1b, Mmp9, Igf1, Ccl2, Slc26a2, Timp1, Lcn2, Il1a, and Sell); the majority of them occupied hub positions in retrieved colitis-associated gene networks (key nodes) (Figure 2A). The rest of the genes were found to be less explored as colitis-related ones (Figure 2B), and, therefore, could be used as a source of novel promising markers/regulators of colitis. To experimentally verify the obtained data, Mmp3, C3, and Tyrobp, displaying, on the one hand, little connection with colitis in the published reports (Figure 2B), and, on the other hand, high degree centrality scores in colitis-associated gene networks (Figure 2A), were selected for further qRT-PCR analysis. Since the profile of MMPs was identified as hypothetically susceptible to transforming acute colitis into CAC (Figure 1F), expressions of Mmp9 and Timp1 (known inhibitor of MMPs) were also further validated.

2.2.2. CAC-Associated Hub Genes

The ranking of CAC-specific core genes according to their degree centrality scores in CAC-related gene networks identified the top 20 genes occupying hub positions, including genes encoding cyto- and chemokines (Tnf, Il1a, Cxcl16), regulators of ECM remodeling (Timp1, Mmp7, Mmp13, Gusb), immune (Ctla4, Cyba) and protective (Gstt1, Hp, Clu, Cyp2s1) response, lipid homeostasis (Acss2, Chpt1), ROS production (Maoa), cell–cell interaction (Cd44), membrane fusion (Snap25), and signal transduction (Plce1, Lgr5) (Figure 2C). Further text mining study, combined with the computing of the association of CAC-related core genes with the overall survival of patients with colon (COAD) and rectal (READ) adenocarcinomas, clearly confirmed the credibility of our bioinformatics analysis: the most reported CAC-related genes (Tnf, Cd44, Timp1, Mmp7, Ctla4, Clu, Il1a, Hp) were not only associated with poor prognosis in COAD and READ patients but also occupied the hub positions in the networks retrieved from CAC-associated DEGs (Figure 2D). These results indicate a probable important regulatory function of the listed core genes in colitis-associated neoplastic transformation of colon tissue.
To identify novel candidate genes for CAC, our attention was centered on the core genes that are, on the one hand, poorly characterized in the field of CAC, and, on the other hand, associated with ECM remodeling susceptible to “inflammation-dysplasia-carcinoma” axis (Figure 1C,F), notably, Mmp13 (key node) and Adam8 (extracellular metalloprotease-disintegrin involved in ECM digestion and markedly associated with pathogenesis of gastrointestinal malignancies [27]) (Figure 2D). In addition, the key nodes Timp1 and Mmp7 previously reported as probable regulators of CRC were also selected for qRT-PCR analysis.

2.3. Validation of Novel Candidate Genes for Colitis and CAC

2.3.1. Murine Model of DSS-Induced Colitis and CAC

Acute colitis was induced in mice by administration of 2.5% DSS solution in drinking water for 7 days, followed by a 3 day recovery (Figure 3A). CAC was induced in mice by single intraperitoneal (i.p.) injection of AOM 1 week before DSS administration. Furthermore, mice were exposed to 3 consecutive cycles of 1.5% DSS instillations for 7 days, followed by 2 weeks of recovery (Figure 3A). After the experiment termination, the colons were separated from the proximal rectum, mechanically cleaned with saline buffer, and collected for subsequent histological analysis and qRT-PCR.
Gross morphological analysis of healthy colons revealed the normal thickness of the colonic wall and mucosa structure (Figure 3B). Administration of 2.5% DSS for one week led to acute inflammatory changes in the colonic tissues, clearly demonstrating the development of acute colitis and represented by thickening of the colonic wall, hyperemia, hemorrhages, and scattered ulcers (Figure 3B). Long-term cyclic administration of 1.5% DSS with prior injections of carcinogen AOM caused the development of multiple adenomas in the distal part of mice colons with a significant decrease in the intensity of acute inflammatory changes in the colonic tissues (Figure 3B).
Histologically, the colon tissue of healthy mice demonstrated intact colon architecture, non-disrupted crypts, and goblet cells with active mucus vacuoles (Figure 3C). Acute administration of DSS caused severe colon tissue damage, represented by massive epithelium disruption with erosions and ulcerations, diffuse destruction of crypts, and loss of mucosal architecture (Figure 3C). Pronounced inflammatory infiltration through the whole colonic wall, due to neutrophils and lymphocytes as well as mucosa edema, was revealed (Figure 3C). In the case of CAC, chronic administration of DSS after AOM injection caused adenomatous transformation of the colon mucosa, represented by multiple adenomas in the colonic tissue with epithelial hyperproliferation and hyperplastic crypts (Figure 3C). Residual inflammatory infiltration located in the mucosa and submucosa of colon tissue with adenomas and represented by lymphocytes and macrophages was detected (Figure 3C). In the colon tissue adjacent to adenomas (colitis in CAC), signs of chronic colonic inflammation with moderate destruction of the mucosal architecture and crypt damage were found (Figure 3C).
Thus, we reproduced the process of colon carcinogenesis, starting with acute inflammation in the colon tissue, transitioning to chronic inflammation, and eventually ending up with the colonic tumor formation.

2.3.2. Core Genes Expression in the Colonic Tissue of Mice with Acute Colitis and CAC

Finally, the expression of the revealed hub genes related to acute colitis (C3, Tyrobp, Mmp3, Mmp9, Timp1) and CAC (Timp1, Adam8, Mmp7, Mmp13) was validated by qRT-PCR in the colon tissue of mice with acute colitis and colitis-driven adenomas (Figure 3D). As expected, the expression of colitis-related genes C3, Tyrobp, Mmp3, Mmp9, and Timp1 was significantly up-regulated in inflamed colon tissue compared with healthy controls; among them, Mmp3 and Timp1 were found to be the most susceptible to acute colitis induction, demonstrating 306.3- and 110.6-fold increases in the expression, respectively, in DSS-treated mice compared with healthy controls (Figure 3D). The chronification of colonic inflammation led to significant reduction in the expression of C3, Tyrobp, Mmp3, and Timp1 in the adjacent to adenomas colonic tissue by 17.5, 6.6, 2.8 and 46.1 times compared with the samples from acute colitis group, and, moreover, the expression of C3 and Tyrobp in this compartment decreased to the healthy level (Figure 3D). Interestingly, chronification of colitis had no obvious effect on the expression of Mmp9: comparable induction of this gene in both DSS- and AOM/DSS-inflamed colon tissues was observed (Figure 3D), which could indicate the important role of Mmp9 in both acute and chronic colon inflammation, agreeing with [28]. The analysis of colonic adenomatous nodes revealed low expression of all the explored acute colitis-associated key genes: the expression levels of C3, Tyrobp, Mmp3, Mmp9, and Timp1 in adenoma tissue were 26.3, 3.4, 20.8, 11.9, and 27.7 times lower than those in the samples with acute colitis (Figure 3D). Note that adenomatous and adjacent tissues in mice with CAC mainly differed in the expression of the following genes: Tyrobp and Timp1 were found to be 1.9 and 1.7 times overexpressed in adenomas compared with the adjacent counterparts, respectively, whereas Mmp3 and Mmp9 were 7.6 and 13.9 times suppressed in tumor tissue, respectively (Figure 3D). Taken together, the obtained results clearly demonstrated that selected key genes associated with acute colitis indeed reached the maximum expression in the acute phase of colon inflammation, whereas chronification of the latter led to a marked decline in this parameter.
As expected, all CAC-associated hub genes (Adam8, Mmp7, Mmp13, and Timp1 (mentioned above)) were characterized by significant overexpression in tumor nodes compared with healthy tissue, which confirms the expediency of their further exploration as CAC-related marker genes (Figure 3D). Interestingly, only Mmp7 and Mmp13 displayed a significantly higher level of activation in colon adenomas compared with both the adjacent tissue (13.8- and 13.4-fold increase, respectively) and colon tissue with acute colitis (186.4- and 19.6-fold increase, respectively). Adam8 and Timp1 mentioned above were also up-regulated in adenomas by 8.3 and 1.7 times compared with the adjacent tissue; however, the maximum of their expression was revealed in the acute colitis samples (86.6- and 110.6-fold increase compared with the healthy group, respectively) (Figure 3D).
Thus, the performed qRT-PCR analysis successfully confirmed the expression of the acute colitis- and CAC-related hub genes identified by the in silico analysis in corresponding murine tissues and clearly demonstrated that colitis-driven colonic adenomatous transformation is accompanied by significant changes in the expression profiles of matrix metalloproteinases, which can be used as a novel prognostic signature for colorectal neoplasia in IBD.

3. Discussion

Despite the large collection of transcriptomics data from IBD and CAC studies, molecular regulators of the transition of colonic inflammatory lesions to cancer have not yet been clearly defined. The current study aimed to reveal core genes involved in the regulation of acute colitis and CAC development in mice and to explore how far their expression profiles changed during the chronification of colon inflammation.
Performed bioinformatics analysis of multiple cDNA microarray datasets of acute colitis and CAC identified a range of core genes associated with the explored pathologies, further functional annotation of which clearly confirmed the reliability of the obtained data. Indeed, high enrichment of acute colitis-related functional terms with pro-inflammatory cytokines and IGF1-Akt signaling pathway (Figure 1C, upper network) agrees well with the proven regulatory role of the latter in colon inflammation and inflammation-induced mucosal injury [29,30]. Along with this, CAC-related core genes were associated with the processes which markedly changed during colitis-driven tumorigenesis (Figure 1C, lower network): it is known that dysplastic and malignant lesions of colon tissue markedly dysregulate sodium transport [31], bile acid secretion [32], and metabolic [33] and oxidative [34] homeostasis.
Interestingly, the analysis of core genes common for both acute colitis and CAC (Figure 1E) also demonstrated the credibility of the performed in silico study. According to the published reports, the acute-phase genes Hp, Lcn2, Lrg1, and Serpina3n included in this list are not only activated in response to inflammatory stimuli, but their aberrant expression is also strongly implicated in tumorigenesis: high levels of Hp and Lcn2 resulted in glucose metabolic dysfunction, angiogenesis, and metastasis in different tumor types [35,36], and Lrg1 and Serpina3n were associated with epithelial–mesenchymal transition in colorectal cancer [37,38]. In addition, the interferon-responsive gene Ifitm3 is critical to early colon cancer development [12,39], along with S100a9 and Slpi, which, when highly expressed in inflamed colon tissues in mice and patients with colitis and IBD, respectively, can be considered as potent amplifiers of tumor invasion [40,41].
Analysis of gene association networks with subsequent processing of obtained results using the text mining approach revealed a range of core genes occupied hub positions in the acute colitis- and CAC-associated regulomes, which had not yet been extensively studied in relation to the explored diseases (acute colitis: C3, Tyrobp, Mmp3; CAC: Adam8, Mmp13) (Figure 2). Further qRT-PCR analysis clearly confirmed the overexpression of the mentioned hub genes in the colon tissue of mice with acute colitis and CAC (Figure 3D) that indicated the expediency of further exploration of these genes as promising novel biomarkers of colon inflammation and colon tumorigenesis.
To independently examine how tightly revealed hub genes were associated with inflammation and colorectal cancer, their sub-networks with first gene neighbors from rodent inflammatome [42] and the gene network related to malignant tumors of the colon (DisGeNET ID: C00071202) were reconstructed and analyzed. As depicted in Figure 4, all explored hub genes, except for Adam8, indeed form tight modules with gene partners within the evaluated regulomes, and are related to diverse processes and signaling pathways important for the pathogenesis of colitis and CAC. For instance, the detection of the functional group “Interleukin-4 and 13 signaling” is in accordance with [43]: a marked IL-13 response from CD4+ natural killer T cells was previously detected in mice with oxazolone-induced colitis and its blockage was found to ameliorate intestinal inflammation and injury. The members of the integrin family (Figure 4A, Timp1-, C3- and Mmp3-centered sub-networks) play a crucial role in the intestinal homing of immune cells and in supporting the inflammatory mechanisms in the gut [44]. uPA-mediated signaling (Figure 4, Timp1-, Mmp3-, Mmp9-centered sub-networks) controls macrophage phagocytosis in intestinal inflammation, and uPA receptor deficiency leads to marked aggravation of experimental colitis in mice [45]. Moreover, uPA-/- mice demonstrated more severe colorectal neoplasia compared with their wild-type littermates [46]. In addition, remodeling of the extracellular matrix is a hallmark of both colitis/IBD [47] and CAC [48], and prostaglandin signaling is involved in the malignant transformation of inflamed intestinal tissue [49].
The detailed comparison of obtained results revealed a group of MMPs as key participants of acute colon inflammation and its transition to malignancy: functional term “Matrix Metalloproteinases” was identified as statistically significant in both acute colitis- and CAC-associated functional annotation maps (Figure 1C), the highly interconnected cluster of MMPs related to different phases of colitis was revealed in the gene network retrieved from computed core genes (Figure 1F), and MMPs occupied hub positions in all analyzed regulomes related to both acute colitis (Figure 2A: Mmp3, Mmp9) and CAC (Figure 2B: Mmp7, Mmp13). Interestingly, the tissue inhibitor of matrix metalloproteinase-1 (Timp1) was also detected as a hub gene specific to both acute colitis and CAC (Figure 1E and Figure 2A,B) and tightly interconnected with MMPs module (Figure 1F), which clearly indicated the importance of Timp1/MMPs balance in colitis-induced tumorigenesis. Indeed, Timp1 is a known regulator of colitis, knockout of which markedly attenuated fibrosis in DSS-inflamed colon tissue [50], and, according to the recent report of Niu et al. [51], a hub gene in colorectal cancer regulome. High expression of MMP3 and MMP9 in mucosa-resident macrophages/neutrophils and IgG plasma cells was detected in patients with IBD [52,53]. According to Pedersen et al. [54], MMP3 and MMP9 are two key enzymes involved in the degradation of intestinal tissue during IBD. Interestingly, the silencing of Mmp3 by siRNA markedly ameliorated DSS-induced colitis in mice [55], whereas knockout of Mmp9 or its pharmacological inhibition surprisingly had no obvious effect on the progression of DSS- and TNBS-stimulated colitis in the murine model [56]. Thus, the master regulatory functions of MMPs in colitis pathogenesis require further clarification: in some cases, their overexpression can be considered as a consequence rather than a cause of intestinal inflammation [56]. In the case of CAC-associated MMPs (Mmp7, Mmp13) revealed in this study, focal high expression of Mmp7 was previously observed in CAC-related dysplastic lesions [48] and its overexpression was associated with tumor growth, metastasis, and worse overall survival in patients with colon cancer [57]. According to Wernicke et al. [58], the up-regulation of MMP-13 was considered as an early predictive cancer biomarker in patients with colon adenoma, which agrees well with the results of our qRT-PCR analysis (Figure 3D). Despite the extensive studies of MMPs as candidate marker genes of colitis and CAC, to the best of our knowledge, the complex evaluation of the expression of Mmp3, Mmp7, Mmp9, and Mmp13 in acutely inflamed, adenomatous, and adjacent colon tissues has not yet been reported. Revealed marked changes in their expression profiles during chronification of colitis (Figure 3D) can be considered as a novel gene signature for predicting CAC.
Besides MMPs, another ECM remodeling player, Adam8, a member of a disintegrin and metalloproteinase family (ADAMs), was identified as a core gene associated with CAC development (Figure 1F, Figure 2D and Figure 3D). Surprisingly, high expression of Adam8 was detected not only in CAC but also in DSS-inflamed colon tissue (Figure 3D). Along with the reorganization of ECM, ADAMs are engaged in the processing of various substrates, including cytokines, growth factors, cell adhesion molecules, and receptors, that determines their important role in a range of pathological processes [59]. The most studied ADAMs in IBD was Adam17, associated with EGFR and STAT3 signaling pathways crucial for the pathogenesis of colitis [60], high epithelial expression of which positively correlated with cell proliferation and goblet cell number in UC patients [61]. To the best of our knowledge, the involvement of Adam8 in the regulation of acute colitis and colitis-induced adenomatous transformation of colon tissue had not yet been reported. Only Christophi et al. and Guo et al. have discussed the overexpression of Adam8 in IBD patients [62] and AOM/DSS-induced colitis in mice [63]. Given the recently demonstrated ability of Adam8 to control neutrophil transmigration [64] and NLRP3 inflammasome activation [65], the processes tightly associated with colon inflammation [66,67], Adam8 can be considered as a novel promising master regulator of colitis and CAC; this requires further clarification. Interestingly, despite the revealed low interconnection of Adam8 with the colon cancer-associated gene network retrieved from DisGeNET (Figure 4), this gene seems to play an important role in the pathogenesis of CAC: Adam8 is involved in the activation of integrin, FAK, ERK1/2, and Akt/PKB signaling pathways related to cancer progression [68], its overexpression was identified in colorectal cancer compared with adjacent normal tissues [69], and the suppression of the expression of Adam8 by knockout or siRNA approaches resulted in reduced proliferation and invasiveness of colon cancer cells [69,70].
Finally, C3 and Tyrobp were also revealed as colitis-specific hub genes (Figure 2A and Figure 3D), which is in line with published reports. Previously, a high level of C3 in the serum and jejunal secretion of IBD patients was identified [71,72]. Moreover, C3 was found to be up-regulated in intestinal epithelial cells in the DSS-induced colitis model [73], and its ablation promoted inflammatory responses in the mid colon [74] and significantly reinforced DSS-induced colitis in C3 knockout mice compared with wild-type littermates [72]. Tyrobp is a known regulator of the production of pro-inflammatory mediators in macrophages and neutrophils [75], and, thus, is implicated in pathogenesis of various inflammation-associated diseases [75,76,77]. According to recent studies, Tyrobp was identified as a probable upstream regulator of UC [78], and its knockout robustly attenuated the severity of DSS-induced colitis in mice, whereas its overexpression resulted in a striking exacerbation of colon damage caused by DSS [79].
The published works discussed above demonstrated the involvement of the revealed core genes in the regulation of inflammation and malignant lesion of the colon, not only in murine models but also in patients. To independently confirm the translational bridge between our findings and the pathogenesis of colitis/CAC in humans, expression of core genes (acute colitis: C3, Tyrobp, Mmp3, Mmp9, Timp1; CAC: Timp1, Mmp7, Mmp13, Adam8) was further evaluated in the transcriptomics profiles of colon tissue from patients with UC and CD collected from GEO (Figure 5A) and colorectal cancer retrieved from The Cancer Genome Atlas (TCGA) (Figure 5B). As depicted in Figure 5A, the majority of the explored key genes were overexpressed in IBD and demonstrated more pronounced susceptibility to the induction of UC compared with CD, except for TYROBP, expression of which was more up-regulated in CD patients. Interestingly, despite the proven association with CAC (Figure 2B,D), TIMP1, MMP7, and ADAM8 were activated in IBD-affected colon tissues (Figure 5A), which is fully in line with our data: the high expression of these genes was demonstrated in DSS-inflamed and adjacent to adenomas colon tissues in mice (Figure 3D). In addition, similar to our results (Figure 3D), CAC-specific MMP13 was found to be slightly associated with IBD: its low activation in two of the four analyzed UC transcriptomics datasets and unchanged levels in CD samples were observed (Figure 5A). Presumably, Mmp13 plays a minor role in ECM remodeling in colitis, whereas CAC was associated with significant up-regulation of its expression, which makes Mmp13 a promising gene candidate for the predicting of colitis-associated tumorigenesis; this requires further detailed study. TCGA analysis of the identified CAC-related core genes revealed a significant association between high expression of TIMP1 and ADAM8 with low overall survival of patients with both colon (COAD) and rectal (READ) adenocarcinomas (Figure 5B). Despite the finding that Timp1 and Adam8 can play important regulatory functions in CAC, this supposition requires further detailed confirmation, since TCGA analysis was performed without consideration of the ratio of UC- and CD-associated CAC patients in COAD and READ cohorts. In addition, given recently reported sex disparities in the association of Timp1 expression with cancer progression [80], further exploration of its regulatory role in CAC in mice of both sexes is needed.
The obtained results were finally summarized in the scheme depicted in Figure 5C. According to our findings, (a) revealed core genes not only occupy hub positions within explored acute colitis- and CAC-specific regulomes, but also are interconnected with each other, (b) Timp1 is identified as a hub node in gene association networks retrieved for both acute colitis and CAC, which can indicate its crucial role in colitis-associated tumorigenesis, (c) chronification of colonic inflammation is accompanied by a switch in MMPs profile (acute colitis: Mmp3, Mmp9; CAC: Mmp7, Mmp13), which can serve as a gene signature panel for prognosis of malignant transformation of inflamed colon tissue; and (d) identified core genes are overexpressed in the colon tissue of patients with IBD (all explored genes) and highly aggressive colorectal cancer (TIMP1, ADAM8), confirming the interest in studying these genes within the framework of intestinal pathologies in humans (Figure 5C).

Limitations of the Study

The limitations of the study are as follows: First, given the relatively low number of mice used for experimental validation of the obtained data (n = 6), and their belonging to only one sex (female) and one strain (C57Bl6), further study is required to validate the results using a larger sample size obtained from mice of both sexes and different strains. Second, considering that our findings are predominantly animal-based, to more clearly elucidate how closely (if at all) the identified core genes are involved in the regulation of intestinal pathologies in humans, revealed translational bridge needs further large-scale verification study, using clinical samples of patients with UC, CD, and UC/CD-associated colorectal cancer. Third, despite the identification of high degree centrality scores of the explored key genes and their tight association with crucial colitis/CAC-related signaling pathways, the master regulatory functions of these genes in colitis and CAC should be further verified experimentally (for instance, using knockout models).

4. Materials and Methods

4.1. Microarray Data Collection and Differential Expression Analysis

The gene expression profiles associated with murine acute colitis and CAC, as well as ulcerative colitis and Crohn’s disease, in patients were acquired from the Gene Expression Omnibus database [81] (Table 1). The fold changes between the mean expression values of the genes in the experimental (pathology) versus control groups were computed using the GEO2R tool [82]. The Benjamini–Hochberg false discovery rate method was selected for adjusting p-values. The genes with a p-value < 0.05 and |fold change| > 1.5 were identified as differentially expressed genes (DEGs) and were collected for further analysis. Overlapping of the DEGs from different datasets was performed using the InteractiVenn tool [83]. Hierarchical clustering of DEGs according to their expression profiles was carried out using the Euclidean distance metric, using the Morpheus tool (https://software.broadinstitute.org/morpheus, accessed on 12 December 2022).

4.2. Functional Analysis of DEGs

Functional annotation of acute colitis- and CAC-associated DEGs was performed using the ClueGO 2.5.7 plugin in Cytoscape 3.7.2, using the latest updates of Gene Ontology (Biological Processes), Kyoto Encyclopedia of Genes and Genomes (KEGG), WikiPathways, and REACTOME databases. The GO Tree interval was ranged from 3 to 8 and the minimum number of genes per cluster was set to 3. Enrichment of functional terms was tested using the two-sided hypergeometric test corrected using the Bonferroni method, followed by selecting significantly enriched terms with a p-value < 0.05. To cluster similar functional groups retrieved from different databases in the common pathway-specific modules, the GO Term Fusion was used. Functional grouping of finally selected functional terms was performed using kappa statistics (kappa score ≥ 0.4). Functional annotation of gene modules, consisting of core genes and their first gene partners extracted from murine inflammatome and colon cancer-related regulome, was performed using the ToppFun tool (databases: KEGG, REACTOME, MSigDB C2 BIOCARTA, BioSystems: Pathway Interaction Database, Pathway Ontology; Bonferroni adjustment) [84].

4.3. Reconstruction of Gene Association Networks

Gene association networks were reconstructed from the genes of interest using the Search Tool for the Retrieval of Interaction Genes (STRING) database, using the stringApp 1.5.1 tool [85], and were visualized using Cytoscape 3.7.2. The cutoff criterion of the confidence score was set as >0.7 to eliminate inconsistent “gene–gene” pairs from the dataset. The number of neighbors of a gene of interest within reconstructed networks was calculated using the NetworkAnalyzer plugin [86] and visualized using the Morpheus platform [87].

4.4. Data Mining Analysis

The search for the co-occurrence of the names of core genes with various colitis- and CAC-related terms in the same sentences in abstracts of published reports deposited in the MEDLINE database was performed using the GenCLiP3 tool [88], with the following settings: impact factor of 0–50 and year of publication of 1992–2022. The results were visualized using Circos [89].

4.5. Murine Models of Acute Colitis and Colitis-Associated Cancer (CAC)

Eight-week-old female C57Bl6 mice with an average weight of 22–24 g were obtained from the Vivarium of the Institute of Chemical Biology and Fundamental Medicine SB RAS (Novosibirsk, Russia). Mice were housed in plastic cages (7 animals per cage) under normal daylight conditions. Water and food were provided ad libitum. Experiments were carried out in accordance with the European Communities Council Directive 86/609/CEE. The experimental protocols were approved by the Committee on the Ethics of Animal Experiments at the Institute of Cytology and Genetics SB RAS (Novosibirsk, Russia) (protocol No. 56 from 10 August 2019).
Acute colitis was induced in mice (n = 10) by administration of 2.5% DSS solution in drinking water for 7 days, followed by 3 days of recovery. Mice were sacrificed on day 10 after colitis initiation. CAC was induced in mice (n = 10) by a single intraperitoneal (i.p.) injection of carcinogen AOM (10 mg/kg) 1 week before DSS administration, as described in [90]. Furthermore, mice were exposed to 3 consecutive cycles of 1.5% DSS instillations with drinking water for 7 days, followed by 2 weeks of recovery. The mice were sacrificed 10 weeks after the start of the experiment. At the end of the study, the colons were separated from the proximal rectum, mechanically cleaned with saline buffer, and were then collected. Only 8 of 10 samples had well-formed adenomas in the colon, which were selected for the subsequent gross examination, histological analysis, and qRT–PCR.

4.6. Histology

For the histological study, colon specimens were fixed in 10% neutral-buffered formalin (BioVitrum, Moscow, Russia), dehydrated in ascending ethanol and xylols, and embedded in HISTOMIX paraffin (BioVitrum, Moscow, Russia). The paraffin sections (5 μm) were sliced on a Microm HM 355S microtome (Thermo Fisher Scientific, Waltham, MA, USA) and stained with haematoxylin and eosin. The images were examined and scanned using an Axiostar Plus microscope equipped with an Axiocam MRc5 digital camera (Zeiss, Oberkochen, Germany) at magnifications of ×100.

4.7. Quantitative Real-Time PCR (qRT-PCR)

Total RNA was isolated from the colons of experimental animals using TRIzol reagent (Ambion, Austin, TX, USA) according to the manufacturer’s instructions. Briefly, colon tissue was collected in 1.5 mL capped tubes, filled with 1 g of lysing matrix D (MP Biomedicals, Irvine, CA, USA) and 1 mL of TRIzol reagent, then homogenized using a FastPrep-24 TM 5G homogenizer (MP Biomedicals, Irvine, CA, USA) with QuickPrep 24 adapter. The homogenization was performed at 6.0 m/s for 40 s. After homogenization, the content of the tubes was transferred to the new 1.5 mL tubes without lysing matrix. Total RNA extraction was performed according to the TRIzol reagent protocol.
Due to the known ability of DSS to linger in the RNA extracted from the colon tissue, and, thus, interfere with both reverse transcription and PCR reactions, the extracted total RNA was diluted to a volume of 250 μL and purified using Microcon Centrifugal Filter Devices (MilliPore, Burlington, MA, USA) by centrifuging for 1 h at 14,000× g.
The first strand of cDNA was synthesized from total RNA (n = 6 per group, the samples with the highest RNA purity and integrity) in 100 μL of reaction mixture containing 2.5 μg of total RNA, 20 μL of 5× RT buffer (Biolabmix, Novosibirsk, Russia), 250 U of M-MuLV-RH revertase (Biolabmix, Novosibirsk, Russia), and 100 μM of dT(15) diluted to a volume of 100 μL. Reverse transcription was performed at 25 °C for 10 min followed by the incubation at 42 °C for 60 min with subsequent termination at 70 °C for 10 min.
Amplification of cDNA was performed in a 25 μL PCR reaction mixture containing 5 μL of cDNA, 12.5 μL of HS-qPCR (2×) master mix (Biolabmix, Novosibirsk, Russia), 0.25 μM each of the forward and reverse primers to Hprt and Hprt specific ROX-labeled probe, 0.25 μM each of the forward and reverse gene-specific primers, and FAM-labeled probe (Table 2). Amplification was performed as follows: (1) 94 °C, 2 min; (2) 94 °C, 10 s; (3) 60 °C, 30 s (steps 2–3: 50 cycles). The relative level of gene expression was normalized to the level of Hprt expression according to the ΔΔCt method.
Amplification was performed using a C1000 Touch with CFX96 module Real-Time system (BioRad, Hercules, CA, USA), and the relative level of gene expression was calculated using BioRad CFX manager software (BioRad, Hercules, CA, USA). Three to five samples from each experimental group were analyzed in triplicate. The sequences of the primers used in the study are listed in Table 2.

4.8. The Association of DEGs Expression with Survival Rates of Patients with Colorectal Cancer

To explore the association of revealed core genes with the progression of colon (COAD) and rectal (READ) adenocarcinomas, analysis of the survival rates and their correlation with the expression of studied genes was performed using The Cancer Genome Atlas (TCGA) clinical data for patients with COAD and READ. Kaplan–Meier survival curves for COAD and READ patients depending on the mRNA expression level of core genes were constructed using the OncoLnc tool [91].

4.9. Statistical Analysis

The statistical analysis was performed using Benjamini–Hochberg false discovery rate method (identification of DEGs; GEO2R tool), two-sided hypergeometric test with Bonferroni correction (functional analysis of DEGs; ClueGO plugin and ToppFunn tool), and two-tailed unpaired Student’s t-test (qRT-PCR analysis; Microsoft Excel). p-values of less than 0.05 were considered statistically significant.

5. Conclusions

In summary, this animal-based research revealed a range of core genes associated with acute colitis (C3, Tyrobp, Mmp3, Mmp9, Timp1) and CAC (Timp1, Mmp7, Mmp13) in mice. The observed high rate of interconnection of these genes with gene networks retrieved for intestinal inflammation and malignancy, their significant association with key colitis/CAC-related signaling pathways, and probable involvement in the pathogenesis of IBD and colorectal cancer in patients demonstrated the expediency of further detailed studies of identified core genes as novel master regulators and promising therapeutic targets for colitis and CAC.

Author Contributions

Conceptualization, A.V.M.; methodology, A.V.M. and A.V.S.; software, A.V.M.; validation, I.A.S. and A.V.S.; formal analysis, A.V.M., I.A.S., and A.V.S.; investigation, A.V.M., I.A.S., and A.V.S.; resources, A.V.M. and M.A.Z.; data curation, A.V.M. and A.V.S.; writing—original draft preparation, A.V.M. and A.V.S.; writing—review and editing, M.A.Z.; visualization, A.V.M. and A.V.S.; supervision, M.A.Z.; project administration, A.V.M.; funding acquisition, M.A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation (Grant No. 19–74–30011) and the Russian state-funded budget project of ICBFM SB RAS No. 121031300044–5.

Institutional Review Board Statement

Animal studies were carried out in accordance with the European Communities Council Directive 86/609/CEE. The experimental protocols were approved by the Committee on the Ethics of Animal Experiments at the Institute of Cytology and Genetics SB RAS (Novosibirsk, Russia) (protocol No. 56 from 10 August 2019).

Informed Consent Statement

Not applicable since only public cDNA microarray datasets deposited in the Gene Expression Omnibus database were used.

Data Availability Statement

The data presented in this study are openly available on the Gene Expression Omnibus database. Reference numbers: GSE31106, GSE35609, GSE42768, GSE64658, GSE71920, GSE31106, GSE42768, GSE56050, GSE64658, GSE73661, GSE75214, GSE87466, GSE92415, GSE9686, GSE10616, GSE16879, and GSE75214.

Acknowledgments

The authors gratefully thank Anton V. Filatov (Institute of Chemical Biology and Fundamental Medicine SB RAS, Novosibirsk, Russia) for help with the GEO datasets analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bardou, M.; Rouland, A.; Martel, M.; Loffroy, R.; Barkun, A.N.; Chapelle, N. Review: Obesity and colorectal cancer. Aliment. Pharmacol. Ther. 2022, 56, 407–418. [Google Scholar] [CrossRef]
  2. Dai, Z.H.; Tang, M.; Chen, Y.L.; Zhang, T.L.; Li, J.; Lv, G.H.; Yan, Y.G.; Ouyang, Z.H.; Huang, W.; Zou, M.X. Incidence and Risk Factors for Cerebrovascular-Specific Mortality in Patients with Colorectal Cancer: A Registry-Based Cohort Study Involving 563,298 Patients. Cancers 2022, 14, 2053. [Google Scholar] [CrossRef]
  3. Majumder, S.; Shivaji, U.N.; Kasturi, R.; Sigamani, A.; Ghosh, S.; Iacucci, M. Inflammatory bowel disease-related colorectal cancer: Past, present and future perspectives. World J. Gastrointest. Oncol. 2022, 14, 547. [Google Scholar] [CrossRef] [PubMed]
  4. Grivennikov, S.I. Inflammation and colorectal cancer: Colitis-associated neoplasia. Semin. Immunopathol. 2013, 35, 229–244. [Google Scholar] [CrossRef] [PubMed]
  5. M’koma, A.E. Inflammatory Bowel Disease: Clinical Diagnosis and Surgical Treatment-Overview. Medicina 2022, 58, 567. [Google Scholar] [CrossRef] [PubMed]
  6. Hsiao, S.-W.; Yen, H.-H.; Chen, Y.-Y. Chemoprevention of Colitis-Associated Dysplasia or Cancer in Inflammatory Bowel Disease. Gut Liver 2022, 16, 840–848. [Google Scholar] [CrossRef] [PubMed]
  7. Iwanaga, K.; Nakamura, T.; Maeda, S.; Aritake, K.; Hori, M.; Urade, Y.; Ozaki, H.; Murata, T. Mast cell-derived prostaglandin D2 inhibits colitis and colitis-associated colon cancer in mice. Cancer Res. 2014, 74, 3011–3019. [Google Scholar] [CrossRef] [Green Version]
  8. Shah, S.C.; Itzkowitz, S.H. Colorectal Cancer in Inflammatory Bowel Disease: Mechanisms and Management. Gastroenterology 2022, 162, 715–730.e3. [Google Scholar] [CrossRef]
  9. Frigerio, S.; Lartey, D.A.; D’haens, G.R.; Grootjans, J. The Role of the Immune System in IBD-Associated Colorectal Cancer: From Pro to Anti-Tumorigenic Mechanisms. Int. J. Mol. Sci. 2021, 22, 12739. [Google Scholar] [CrossRef]
  10. Rogler, G. Chronic ulcerative colitis and colorectal cancer. Cancer Lett. 2014, 345, 235–241. [Google Scholar] [CrossRef] [Green Version]
  11. Matson, J.; Ramamoorthy, S.; Lopez, N.E. The Role of Biomarkers in Surgery for Ulcerative Colitis: A Review. J. Clin. Med. 2021, 10, 3362. [Google Scholar] [CrossRef] [PubMed]
  12. Suzuki, R.; Miyamoto, S.; Yasui, Y.; Sugie, S.; Tanaka, T. Global gene expression analysis of the mouse colonic mucosa treated with azoxymethane and dextran sodium sulfate. BMC Cancer 2007, 7, 84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shi, L.; Han, X.; Li, J.-X.; Liao, Y.-T.; Kou, F.-S.; Wang, Z.-B.; Shi, R.; Zhao, X.-J.; Sun, Z.-M.; Hao, Y. Identification of differentially expressed genes in ulcerative colitis and verification in a colitis mouse model by bioinformatics analyses. World J. Gastroenterol. 2020, 26, 5983–5996. [Google Scholar] [CrossRef] [PubMed]
  14. Xiu, M.; Liu, Y.; Chen, G.; Hu, C.; Kuang, B. Identifying Hub Genes, Key Pathways and Immune Cell Infiltration Characteristics in Pediatric and Adult Ulcerative Colitis by Integrated Bioinformatic Analysis. Dig. Dis. Sci. 2021, 66, 3002–3014. [Google Scholar] [CrossRef] [PubMed]
  15. Lu, C.; Zhang, X.; Luo, Y.; Huang, J.; Yu, M. Identification of CXCL10 and CXCL11 as the candidate genes involving the development of colitis-associated colorectal cancer. Front. Genet. 2022, 13, 945414. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, M.; Kong, Y.; Chen, N.; Peng, W.; Zi, R.; Jiang, M.; Zhu, J.; Wang, Y.; Yue, J.; Lv, J.; et al. Identification of Immune-Related Gene Signature and Prediction of CeRNA Network in Active Ulcerative Colitis. Front. Immunol. 2022, 13, 855645. [Google Scholar] [CrossRef]
  17. Cheng, C.; Hua, J.; Tan, J.; Qian, W.; Zhang, L.; Hou, X. Identification of differentially expressed genes, associated functional terms pathways, and candidate diagnostic biomarkers in inflammatory bowel diseases by bioinformatics analysis. Exp. Ther. Med. 2019, 18, 278–288. [Google Scholar] [CrossRef] [Green Version]
  18. Hu, W.; Fang, T.; Chen, X. Identification of Differentially Expressed Genes and miRNAs for Ulcerative Colitis Using Bioinformatics Analysis. Front. Genet. 2022, 13, 914384. [Google Scholar] [CrossRef]
  19. Zhang, J.; Wang, X.; Xu, L.; Zhang, Z.; Wang, F.; Tang, X. Investigation of Potential Genetic Biomarkers and Molecular Mechanism of Ulcerative Colitis Utilizing Bioinformatics Analysis. Biomed. Res. Int. 2020, 2020, 4921387. [Google Scholar] [CrossRef] [Green Version]
  20. Shi, W.; Zou, R.; Yang, M.; Mai, L.; Ren, J.; Wen, J.; Liu, Z.; Lai, R. Analysis of Genes Involved in Ulcerative Colitis Activity and Tumorigenesis Through Systematic Mining of Gene Co-expression Networks. Front. Physiol. 2019, 10, 662. [Google Scholar] [CrossRef]
  21. Chen, H.-M.; MacDonald, J.A. Molecular Network Analyses Implicate Death-Associated Protein Kinase 3 (DAPK3) as a Key Factor in Colitis-Associated Dysplasia Progression. Inflamm. Bowel Dis. 2022, 28, 1485–1496. [Google Scholar] [CrossRef]
  22. Ding, H.; Liu, X.; Jian-ming, X.; Qiao, M. Identification of Crucial Genes and Related Transcription Factors in Ulcerative Colitis. Ann. Clin. Lab. Sci. 2021, 51, 245–254. [Google Scholar] [PubMed]
  23. Huang, Y.; Zhang, X.; Wang, P.; Li, Y.; Yao, J. Identification of hub genes and pathways in colitis-associated colon cancer by integrated bioinformatic analysis. BMC Genom. Data 2022, 23, 48. [Google Scholar] [CrossRef] [PubMed]
  24. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  25. Markov, A.V.; Kel, A.E.; Salomatina, O.V.; Salakhutdinov, N.F.; Zenkova, M.A.; Logashenko, E.B. Deep insights into the response of human cervical carcinoma cells to a new cyano enone-bearing triterpenoid soloxolone methyl: A transcriptome analysis. Oncotarget 2019, 10, 5267–5297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yu, D.; Lim, J.; Wang, X.; Liang, F.; Xiao, G. Enhanced construction of gene regulatory networks using hub gene information. BMC Bioinform. 2017, 18, 186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Łukaszewicz-Zając, M.; Pączek, S.; Mroczko, B. A Disintegrin and Metalloproteinase (ADAM) Family—Novel Biomarkers of Selected Gastrointestinal (GI) Malignancies? Cancers 2022, 14, 2307. [Google Scholar] [CrossRef]
  28. Marônek, M.; Marafini, I.; Gardlík, R.; Link, R.; Troncone, E.; Monteleone, G. Metalloproteinases in Inflammatory Bowel Diseases. J. Inflamm. Res. 2021, 14, 1029–1041. [Google Scholar] [CrossRef] [PubMed]
  29. Liu, D.; Saikam, V.; Skrada, K.A.; Merlin, D.; Iyer, S.S. Inflammatory bowel disease biomarkers. Med. Res. Rev. 2022, 42, 1856–1887. [Google Scholar] [CrossRef]
  30. Guijarro, L.G.; Cano-Martínez, D.; Toledo-Lobo, M.V.; Salinas, P.S.; Chaparro, M.; Gómez-Lahoz, A.M.; Zoullas, S.; Rodríguez-Torres, R.; Román, I.D.; Monasor, L.S.; et al. Relationship between IGF-1 and body weight in inflammatory bowel diseases: Cellular and molecular mechanisms involved. Biomed. Pharmacother. 2021, 144, 112239. [Google Scholar] [CrossRef]
  31. Prasad, H.; Visweswariah, S.S. Impaired Intestinal Sodium Transport in Inflammatory Bowel Disease: From the Passenger to the Driver’s Seat. Cell. Mol. Gastroenterol. Hepatol. 2021, 12, 277–292. [Google Scholar] [CrossRef]
  32. Camilleri, M. Bile acid detergency: Permeability, inflammation, and effects of sulfation. Am. J. Physiol. Liver Physiol. 2022, 322, G480–G488. [Google Scholar] [CrossRef] [PubMed]
  33. Shen, Y.; Sun, M.; Zhu, J.; Wei, M.; Li, H.; Zhao, P.; Wang, J.; Li, R.; Tian, L.; Tao, Y.; et al. Tissue metabolic profiling reveals major metabolic alteration in colorectal cancer. Mol. Omics 2021, 17, 464–471. [Google Scholar] [CrossRef] [PubMed]
  34. Cao, Y.; Deng, S.; Yan, L.; Gu, J.; Mao, F.; Xue, Y.; Zheng, C.; Yang, M.; Liu, H.; Liu, L.; et al. An Oxidative Stress Index-Based Score for Prognostic Prediction in Colorectal Cancer Patients Undergoing Surgery. Oxid. Med. Cell. Longev. 2021, 2021, 6693707. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, J.; Cheuk, I.W.-Y.; Siu, M.-T.; Yang, W.; Cheng, A.S.; Shin, V.Y.; Kwong, A. Human haptoglobin contributes to breast cancer oncogenesis through glycolytic activity modulation. Am. J. Cancer Res. 2020, 10, 2865–2877. [Google Scholar]
  36. Santiago-Sánchez, G.S.; Pita-Grisanti, V.; Quiñones-Díaz, B.; Gumpper, K.; Cruz-Monserrate, Z.; Vivas-Mejía, P.E. Biological Functions and Therapeutic Potential of Lipocalin 2 in Cancer. Int. J. Mol. Sci. 2020, 21, 4365. [Google Scholar] [CrossRef]
  37. Camilli, C.; Hoeh, A.E.; De Rossi, G.; Moss, S.E.; Greenwood, J. LRG1: An emerging player in disease pathogenesis. J. Biomed. Sci. 2022, 29, 6. [Google Scholar] [CrossRef]
  38. Soman, A.; Asha Nair, S. Unfolding the cascade of SERPINA3: Inflammation to cancer. Biochim. Biophys. Acta Rev. Cancer 2022, 1877, 188760. [Google Scholar] [CrossRef]
  39. Fan, J.; Peng, Z.; Zhou, C.; Qiu, G.; Tang, H.; Sun, Y.; Wang, X.; Li, Q.; Le, X.; Xie, K. Gene-expression profiling in Chinese patients with colon cancer by coupling experimental and bioinformatic genomewide gene-expression analyses. Cancer 2008, 113, 266–275. [Google Scholar] [CrossRef]
  40. Shabani, F.; Farasat, A.; Mahdavi, M.; Gheibi, N. Calprotectin (S100A8/S100A9): A key protein between inflammation and cancer. Inflamm. Res. 2018, 67, 801–812. [Google Scholar] [CrossRef]
  41. Nugteren, S.; Samsom, J.N. Secretory Leukocyte Protease Inhibitor (SLPI) in mucosal tissues: Protects against inflammation, but promotes cancer. Cytokine Growth Factor Rev. 2021, 59, 22–35. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, I.M.; Zhang, B.; Yang, X.; Zhu, J.; Stepaniants, S.; Zhang, C.; Meng, Q.; Peters, M.; He, Y.; Ni, C.; et al. Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers. Mol. Syst. Biol. 2012, 8, 594. [Google Scholar] [CrossRef] [PubMed]
  43. May, R.D.; Fung, M. Strategies targeting the IL-4/IL-13 axes in disease. Cytokine 2015, 75, 89–116. [Google Scholar] [CrossRef] [PubMed]
  44. Garlatti, V.; Lovisa, S.; Danese, S.; Vetrano, S. The Multiple Faces of Integrin–ECM Interactions in Inflammatory Bowel Disease. Int. J. Mol. Sci. 2021, 22, 10439. [Google Scholar] [CrossRef]
  45. Genua, M.; D’Alessio, S.; Cibella, J.; Gandelli, A.; Sala, E.; Correale, C.; Spinelli, A.; Arena, V.; Malesci, A.; Rutella, S.; et al. The urokinase plasminogen activator receptor (uPAR) controls macrophage phagocytosis in intestinal inflammation. Gut 2015, 64, 589–600. [Google Scholar] [CrossRef]
  46. Karamanavi, E.; Angelopoulou, K.; Lavrentiadou, S.; Tsingotjidou, A.; Abas, Z.; Taitzoglou, I.; Vlemmas, I.; Erdman, S.E.; Poutahidis, T. Urokinase-Type Plasminogen Activator Deficiency Promotes Neoplasmatogenesis in the Colon of Mice. Transl. Oncol. 2014, 7, 174–187.e5. [Google Scholar] [CrossRef] [Green Version]
  47. Derkacz, A.; Olczyk, P.; Olczyk, K.; Komosinska-Vassev, K. The Role of Extracellular Matrix Components in Inflammatory Bowel Diseases. J. Clin. Med. 2021, 10, 1122. [Google Scholar] [CrossRef]
  48. Walter, L.; Harper, C.; Garg, P. Role of Matrix Metalloproteinases in Inflammation/Colitis-Associated Colon Cancer. ImmunoGastroenterology 2013, 2, 22. [Google Scholar] [CrossRef]
  49. Wang, Q.; Morris, R.J.; Bode, A.M.; Zhang, T. Prostaglandin Pathways: Opportunities for Cancer Prevention and Therapy. Cancer Res. 2022, 82, 949–965. [Google Scholar] [CrossRef]
  50. Breynaert, C.; de Bruyn, M.; Arijs, I.; Cremer, J.; Martens, E.; Van Lommel, L.; Geboes, K.; De Hertogh, G.; Schuit, F.; Ferrante, M.; et al. Genetic Deletion of Tissue Inhibitor of Metalloproteinase-1/TIMP-1 Alters Inflammation and Attenuates Fibrosis in Dextran Sodium Sulphate-induced Murine Models of Colitis. J. Crohn’s Colitis 2016, 10, 1336–1350. [Google Scholar] [CrossRef] [Green Version]
  51. Niu, L.; Gao, C.; Li, Y. Identification of potential core genes in colorectal carcinoma and key genes in colorectal cancer liver metastasis using bioinformatics analysis. Sci. Rep. 2021, 11, 23938. [Google Scholar] [CrossRef] [PubMed]
  52. O’Sullivan, S.; Gilmer, J.F.; Medina, C. Matrix Metalloproteinases in Inflammatory Bowel Disease: An Update. Mediat. Inflamm. 2015, 2015, 964131. [Google Scholar] [CrossRef] [PubMed]
  53. Nakov, R. New markers in ulcerative colitis. Clin. Chim. Acta 2019, 497, 141–146. [Google Scholar] [CrossRef] [PubMed]
  54. Pedersen, G.; Saermark, T.; Kirkegaard, T.; Brynskov, J. Spontaneous and cytokine induced expression and activity of matrix metalloproteinases in human colonic epithelium. Clin. Exp. Immunol. 2009, 155, 257–265. [Google Scholar] [CrossRef]
  55. Kobayashi, K.; Arimura, Y.; Goto, A.; Okahara, S.; Endo, T.; Shinomura, Y.; Imai, K. Therapeutic implications of the specific inhibition of causative matrix metalloproteinases in experimental colitis induced by dextran sulphate sodium. J. Pathol. 2006, 209, 376–383. [Google Scholar] [CrossRef]
  56. de Bruyn, M.; Breynaert, C.; Arijs, I.; De Hertogh, G.; Geboes, K.; Thijs, G.; Matteoli, G.; Hu, J.; Van Damme, J.; Arnold, B.; et al. Inhibition of gelatinase B/MMP-9 does not attenuate colitis in murine models of inflammatory bowel disease. Nat. Commun. 2017, 8, 15384. [Google Scholar] [CrossRef] [Green Version]
  57. Liao, H.-Y.; Da, C.-M.; Liao, B.; Zhang, H.-H. Roles of matrix metalloproteinase-7 (MMP-7) in cancer. Clin. Biochem. 2021, 92, 9–18. [Google Scholar] [CrossRef]
  58. Wernicke, A.-K.; Churin, Y.; Sheridan, D.; Windhorst, A.; Tschuschner, A.; Gattenlöhner, S.; Roderfeld, M.; Roeb, E. Matrix metalloproteinase-13 refines pathological staging of precancerous colorectal lesions. Oncotarget 2016, 7, 73552–73557. [Google Scholar] [CrossRef] [Green Version]
  59. Lambrecht, B.N.; Vanderkerken, M.; Hammad, H. The emerging role of ADAM metalloproteinases in immunity. Nat. Rev. Immunol. 2018, 18, 745–758. [Google Scholar] [CrossRef]
  60. Calligaris, M.; Cuffaro, D.; Bonelli, S.; Spanò, D.P.; Rossello, A.; Nuti, E.; Scilabra, S.D. Strategies to Target ADAM17 in Disease: From Its Discovery to the iRhom Revolution. Molecules 2021, 26, 944. [Google Scholar] [CrossRef]
  61. Blaydon, D.C.; Biancheri, P.; Di, W.-L.; Plagnol, V.; Cabral, R.M.; Brooke, M.A.; van Heel, D.A.; Ruschendorf, F.; Toynbee, M.; Walne, A.; et al. Inflammatory Skin and Bowel Disease Linked to ADAM17 Deletion. N. Engl. J. Med. 2011, 365, 1502–1508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Christophi, G.P.; Rong, R.; Holtzapple, P.G.; Massa, P.T.; Landas, S.K. Immune Markers and Differential Signaling Networks in Ulcerative Colitis and Crohn’s Disease. Inflamm. Bowel Dis. 2012, 18, 2342–2356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Guo, Y.; Su, Z.-Y.; Zhang, C.; Gaspar, J.M.; Wang, R.; Hart, R.P.; Verzi, M.P.; Kong, A.-N.T. Mechanisms of colitis-accelerated colon carcinogenesis and its prevention with the combination of aspirin and curcumin: Transcriptomic analysis using RNA-seq. Biochem. Pharmacol. 2017, 135, 22–34. [Google Scholar] [CrossRef] [PubMed]
  64. Conrad, C.; Yildiz, D.; Cleary, S.J.; Margraf, A.; Cook, L.; Schlomann, U.; Panaretou, B.; Bowser, J.L.; Karmouty-Quintana, H.; Li, J.; et al. ADAM8 signaling drives neutrophil migration and ARDS severity. JCI Insight 2022, 7, e87489. [Google Scholar] [CrossRef] [PubMed]
  65. Lu, H.; Meng, Y.; Han, X.; Zhang, W. ADAM8 Activates NLRP3 Inflammasome to Promote Cerebral Ischemia-Reperfusion Injury. J. Healthc. Eng. 2021, 2021, 3097432. [Google Scholar] [CrossRef] [PubMed]
  66. Wang, Y.; Wang, K.; Han, G.-C.; Wang, R.-X.; Xiao, H.; Hou, C.-M.; Guo, R.-F.; Dou, Y.; Shen, B.-F.; Li, Y.; et al. Neutrophil infiltration favors colitis-associated tumorigenesis by activating the interleukin-1 (IL-1)/IL-6 axis. Mucosal Immunol. 2014, 7, 1106–1115. [Google Scholar] [CrossRef] [Green Version]
  67. Perera, A.P.; Sajnani, K.; Dickinson, J.; Eri, R.; Körner, H. NLRP3 inflammasome in colitis and colitis-associated colorectal cancer. Mamm. Genome 2018, 29, 817–830. [Google Scholar] [CrossRef]
  68. Conrad, C.; Benzel, J.; Dorzweiler, K.; Cook, L.; Schlomann, U.; Zarbock, A.; Slater, E.P.; Nimsky, C.; Bartsch, J.W. ADAM8 in invasive cancers: Links to tumor progression, metastasis, and chemoresistance. Clin. Sci. 2019, 133, 83–99. [Google Scholar] [CrossRef]
  69. Yang, Z.; Bai, Y.; Huo, L.; Chen, H.; Huang, J.; Li, J.; Fan, X.; Yang, Z.; Wang, L.; Wang, J. Expression of A disintegrin and metalloprotease 8 is associated with cell growth and poor survival in colorectal cancer. BMC Cancer 2014, 14, 568. [Google Scholar] [CrossRef] [Green Version]
  70. Park Bin, G.; Choi, S.; Yoon Sang, Y.; Kim, D. TrkB/C-induced HOXC6 activation enhances the ADAM8-mediated metastasis of chemoresistant colon cancer cells. Mol. Med. Rep. 2021, 23, 423. [Google Scholar] [CrossRef]
  71. Okada, K.; Itoh, H.; Ikemoto, M. Serum complement C3 and α2-macroglobulin are potentially useful biomarkers for inflammatory bowel disease patients. Heliyon 2021, 7, e06554. [Google Scholar] [CrossRef] [PubMed]
  72. Sina, C.; Kemper, C.; Derer, S. The intestinal complement system in inflammatory bowel disease: Shaping intestinal barrier function. Semin. Immunol. 2018, 37, 66–73. [Google Scholar] [CrossRef] [PubMed]
  73. Sünderhauf, A.; Skibbe, K.; Preisker, S.; Ebbert, K.; Verschoor, A.; Karsten, C.M.; Kemper, C.; Huber-Lang, M.; Basic, M.; Bleich, A.; et al. Regulation of epithelial cell expressed C3 in the intestine—Relevance for the pathophysiology of inflammatory bowel disease? Mol. Immunol. 2017, 90, 227–238. [Google Scholar] [CrossRef] [PubMed]
  74. Choi, Y.J.; Kim, J.E.; Lee, S.J.; Gong, J.E.; Jin, Y.J.; Lee, H.; Hwang, D.Y. Promotion of the inflammatory response in mid colon of complement component 3 knockout mice. Sci. Rep. 2022, 12, 1700. [Google Scholar] [CrossRef]
  75. Subramanian, S.; Pallati, P.K.; Sharma, P.; Agrawal, D.K.; Nandipati, K.C. Significant association of TREM-1 with HMGB1, TLRs and RAGE in the pathogenesis of insulin resistance in obese diabetic populations. Am. J. Transl. Res. 2017, 9, 3224–3244. [Google Scholar]
  76. Fan, D.; He, X.; Bian, Y.; Guo, Q.; Zheng, K.; Zhao, Y.; Lu, C.; Liu, B.; Xu, X.; Zhang, G.; et al. Triptolide Modulates TREM-1 Signal Pathway to Inhibit the Inflammatory Response in Rheumatoid Arthritis. Int. J. Mol. Sci. 2016, 17, 498. [Google Scholar] [CrossRef] [Green Version]
  77. Tammaro, A.; Stroo, I.; Rampanelli, E.; Blank, F.; Butter, L.M.; Claessen, N.; Takai, T.; Colonna, M.; Leemans, J.C.; Florquin, S.; et al. Role of TREM1-DAP12 in Renal Inflammation during Obstructive Nephropathy. PLoS ONE 2013, 8, e82498. [Google Scholar] [CrossRef]
  78. Li, X.; Lee, E.J.; Gawel, D.R.; Lilja, S.; Schäfer, S.; Zhang, H.; Benson, M. Meta-Analysis of Expression Profiling Data Indicates Need for Combinatorial Biomarkers in Pediatric Ulcerative Colitis. J. Immunol. Res. 2020, 2020, 8279619. [Google Scholar] [CrossRef] [Green Version]
  79. Biagioli, M.; Mencarelli, A.; Carino, A.; Cipriani, S.; Marchianò, S.; Fiorucci, C.; Donini, A.; Graziosi, L.; Baldelli, F.; Distrutti, E.; et al. Genetic and Pharmacological Dissection of the Role of Spleen Tyrosine Kinase (Syk) in Intestinal Inflammation and Immune Dysfunction in Inflammatory Bowel Diseases. Inflamm. Bowel Dis. 2018, 24, 123–135. [Google Scholar] [CrossRef] [Green Version]
  80. Hermann, C.D.; Schoeps, B.; Eckfeld, C.; Munkhbaatar, E.; Kniep, L.; Prokopchuk, O.; Wirges, N.; Steiger, K.; Häußler, D.; Knolle, P.; et al. TIMP1 expression underlies sex disparity in liver metastasis and survival in pancreatic cancer. J. Exp. Med. 2021, 218, e20210911. [Google Scholar] [CrossRef]
  81. Clough, E.; Barrett, T. The Gene Expression Omnibus Database BT—Statistical Genomics: Methods and Protocols; Mathé, E., Davis, S., Eds.; Springer: New York, NY, USA, 2016; pp. 93–110. ISBN 978-1-4939-3578-9. [Google Scholar]
  82. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets—Update. Nucleic Acids Res. 2013, 41, D991–D995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Heberle, H.; Meirelles, V.G.; da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Chen, J.; Bardes, E.E.; Aronow, B.J.; Jegga, A.G. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009, 37, W.305–W.311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef]
  86. Assenov, Y.; Ramírez, F.; Schelhorn, S.E.S.E.; Lengauer, T.; Albrecht, M. Computing topological parameters of biological networks. Bioinformatics 2008, 24, 282–284. [Google Scholar] [CrossRef] [Green Version]
  87. MORPHEUS. Versatile Matrix Visualization and Analysis Software. Available online: https://software.broadinstitute.org/morpheus/ (accessed on 12 November 2022).
  88. Wang, J.-H.; Zhao, L.-F.; Wang, H.-F.; Wen, Y.-T.; Jiang, K.-K.; Mao, X.-M.; Zhou, Z.-Y.; Yao, K.-T.; Geng, Q.-S.; Guo, D.; et al. GenCLiP 3: Mining human genes’ functions and regulatory networks from PubMed based on co-occurrences and natural language processing. Bioinformatics 2019, 36, 1973–1975. [Google Scholar] [CrossRef]
  89. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [Green Version]
  90. Tanaka, T.; Kohno, H.; Suzuki, R.; Yamada, Y.; Sugie, S.; Mori, H. A novel inflammation-related mouse colon carcinogenesis model induced by azoxymethane and dextran sodium sulfate. Cancer Sci. 2003, 94, 965–973. [Google Scholar] [CrossRef]
  91. Anaya, J. OncoLnc: Linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Comput. Sci. 2016, 2016, e67. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Core genes involved in the development of acute colitis and CAC in mice revealed by bioinformatics analysis. (A) Venn diagrams overlapping differentially expressed genes (DEGs) identified by re-analysis of cDNA microarray datasets of colon tissue of mice with colitis and CAC. (B) Heatmap demonstrating expression levels of DEGs in acute colitis-related GSE datasets (acute colitis vs. healthy control). (C) Functional analysis of overlapped DEGs associated with acute colitis and CAC. Enrichment for Gene Ontology (Biological Processes), KEGG, REACTOME, and WikiPathways terms performed using ClueGO plugin in Cytoscape. Only pathways with p < 0.05 after Bonferroni correction for multiple testing were included in the functional annotation map. (D) Heatmap demonstrating expression levels of DEGs in CAC-related GSE datasets (CAC vs. healthy control). (E) Venn diagrams overlapping revealed acute colitis- and CAC-associated core genes. (F) Gene association network retrieved from identified core genes associated with acute colitis and CAC reconstructed using STRING database (confidence score ≥ 0.7). Degree—the number of first neighbors (gene partners) within the gene network.
Figure 1. Core genes involved in the development of acute colitis and CAC in mice revealed by bioinformatics analysis. (A) Venn diagrams overlapping differentially expressed genes (DEGs) identified by re-analysis of cDNA microarray datasets of colon tissue of mice with colitis and CAC. (B) Heatmap demonstrating expression levels of DEGs in acute colitis-related GSE datasets (acute colitis vs. healthy control). (C) Functional analysis of overlapped DEGs associated with acute colitis and CAC. Enrichment for Gene Ontology (Biological Processes), KEGG, REACTOME, and WikiPathways terms performed using ClueGO plugin in Cytoscape. Only pathways with p < 0.05 after Bonferroni correction for multiple testing were included in the functional annotation map. (D) Heatmap demonstrating expression levels of DEGs in CAC-related GSE datasets (CAC vs. healthy control). (E) Venn diagrams overlapping revealed acute colitis- and CAC-associated core genes. (F) Gene association network retrieved from identified core genes associated with acute colitis and CAC reconstructed using STRING database (confidence score ≥ 0.7). Degree—the number of first neighbors (gene partners) within the gene network.
Ijms 24 04311 g001
Figure 2. Core genes characterized by hub position within acute colitis- and CAC-related regulomes. (A,C) Heatmaps demonstrating interconnection of core genes in gene association networks reconstructed for each acute colitis- (A) and CAC-related (C) dataset using the STRING database (confidence score ≥ 0.7). Only the top 20 of the most interconnected hub genes are represented. (B,D) Co-occurrence of revealed hub genes with keywords associated with colitis (B) and CAC (D) in the scientific literature deposited in the MEDLINE database. Analysis was performed using the GenClip3 tool. Data were visualized using Circos. COAD and READ-association of the aberrant expression of hub genes with poor prognosis in patients with colon and rectal adenocarcinomas, respectively, according to the TCGA database.
Figure 2. Core genes characterized by hub position within acute colitis- and CAC-related regulomes. (A,C) Heatmaps demonstrating interconnection of core genes in gene association networks reconstructed for each acute colitis- (A) and CAC-related (C) dataset using the STRING database (confidence score ≥ 0.7). Only the top 20 of the most interconnected hub genes are represented. (B,D) Co-occurrence of revealed hub genes with keywords associated with colitis (B) and CAC (D) in the scientific literature deposited in the MEDLINE database. Analysis was performed using the GenClip3 tool. Data were visualized using Circos. COAD and READ-association of the aberrant expression of hub genes with poor prognosis in patients with colon and rectal adenocarcinomas, respectively, according to the TCGA database.
Ijms 24 04311 g002
Figure 3. Morphological changes and expression levels of key genes identified by the bioinformatics analysis in the colon tissues of healthy mice and mice with acute colitis and colitis-associated cancer (CAC). (A) Experimental setup. Acute colitis was induced by the administration of 2.5% dextran sulfate sodium (DSS) solution in drinking water for 7 days followed by a 3 day recovery (upper panel) (n = 10). CAC was induced in mice by single intraperitoneal (i.p.) injection of azoxymethane (AOM) 1 week before DSS administration. Mice were then exposed to 3 consecutive cycles of 1.5% DSS instillations for 7 days followed by 2 weeks recovery (lower panel) (n = 10). At the end of the experiment, the colons were collected for subsequent histological analysis and qRT-PCR. (B,C) Gross morphology (B) and histology (C) of the colon tissue of healthy mice and mice with acute colitis and CAC. Hematoxylin and eosin staining. Original magnification ×100. Black arrows indicate inflammatory infiltration. Green arrows indicate adenomatous transformation of colon epithelium. (D) Expression levels of identified genes in the colon tissues measured by qRT-PCR. HPRT was used as a reference gene. Healthy—healthy colon tissue, colitis—colon tissue with acute colitis, colitis in CAC—colon tissue with chronic colitis adjacent to adenomas. The data are expressed as the mean ± SD (n = 6). * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 3. Morphological changes and expression levels of key genes identified by the bioinformatics analysis in the colon tissues of healthy mice and mice with acute colitis and colitis-associated cancer (CAC). (A) Experimental setup. Acute colitis was induced by the administration of 2.5% dextran sulfate sodium (DSS) solution in drinking water for 7 days followed by a 3 day recovery (upper panel) (n = 10). CAC was induced in mice by single intraperitoneal (i.p.) injection of azoxymethane (AOM) 1 week before DSS administration. Mice were then exposed to 3 consecutive cycles of 1.5% DSS instillations for 7 days followed by 2 weeks recovery (lower panel) (n = 10). At the end of the experiment, the colons were collected for subsequent histological analysis and qRT-PCR. (B,C) Gross morphology (B) and histology (C) of the colon tissue of healthy mice and mice with acute colitis and CAC. Hematoxylin and eosin staining. Original magnification ×100. Black arrows indicate inflammatory infiltration. Green arrows indicate adenomatous transformation of colon epithelium. (D) Expression levels of identified genes in the colon tissues measured by qRT-PCR. HPRT was used as a reference gene. Healthy—healthy colon tissue, colitis—colon tissue with acute colitis, colitis in CAC—colon tissue with chronic colitis adjacent to adenomas. The data are expressed as the mean ± SD (n = 6). * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 24 04311 g003
Figure 4. The involvement of the identified acute colitis- and CAC-related hub genes in murine inflammatome (A) and the gene network associated with malignant tumor of the colon, retrieved from DisGeNET (B). The modules containing the explored genes with their gene partners were extracted from gene networks created using the STRING database and visualized using Cytoscape, according to Section 4.3. Functional analysis of gene modules was performed using the ToppFun tool.
Figure 4. The involvement of the identified acute colitis- and CAC-related hub genes in murine inflammatome (A) and the gene network associated with malignant tumor of the colon, retrieved from DisGeNET (B). The modules containing the explored genes with their gene partners were extracted from gene networks created using the STRING database and visualized using Cytoscape, according to Section 4.3. Functional analysis of gene modules was performed using the ToppFun tool.
Ijms 24 04311 g004
Figure 5. Expression levels of identified DEGs in human patients with inflammatory bowel diseases (ulcerative colitis and Crohn’s disease) and their association with the survival of patients with colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ). (A) Heatmap demonstrating the expression levels of identified acute colitis- and CAC-related hub genes in datasets of human patients with ulcerative colitis and Crohn’s disease. The differential analysis of cDNA microarray data was performed according to Section 4.1. The heat map was constructed using Morpheus. Log2FC = Log2(fold change). In all datasets, biopsy samples were taken from the colon at the sites of active inflammation. (B) Survival of patients with COAD and READ depending on the expression levels of identified DEGs (TIMP1 and ADAM8) in the colon tissue. Kaplan–Meier survival curves were constructed using TCGA data according to Section 4.8. (C) The final scheme of core genes associated with acute colitis and CAC revealed in the current study, their expression profiles in murine models, and association with IBD and colorectal cancer in humans.
Figure 5. Expression levels of identified DEGs in human patients with inflammatory bowel diseases (ulcerative colitis and Crohn’s disease) and their association with the survival of patients with colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ). (A) Heatmap demonstrating the expression levels of identified acute colitis- and CAC-related hub genes in datasets of human patients with ulcerative colitis and Crohn’s disease. The differential analysis of cDNA microarray data was performed according to Section 4.1. The heat map was constructed using Morpheus. Log2FC = Log2(fold change). In all datasets, biopsy samples were taken from the colon at the sites of active inflammation. (B) Survival of patients with COAD and READ depending on the expression levels of identified DEGs (TIMP1 and ADAM8) in the colon tissue. Kaplan–Meier survival curves were constructed using TCGA data according to Section 4.8. (C) The final scheme of core genes associated with acute colitis and CAC revealed in the current study, their expression profiles in murine models, and association with IBD and colorectal cancer in humans.
Ijms 24 04311 g005
Table 1. The GEO microarray datasets used in the study.
Table 1. The GEO microarray datasets used in the study.
ObjectPathologyGEO IDMurine StrainSexAnalyzed GroupsNumber of DSS Cycles
MiceAcute colitisGSE31106ICRm3 untreated mice, 3 AOM/DSS-treated mice1
GSE35609ICRf4 untreated mice, 4 TNBS-treated mice
GSE42768C57Bl6f3 untreated mice, 3 DSS-treated mice1
GSE64658C57Bl6f6 untreated mice, 3 AOM/DSS-treated mice1
GSE71920C57Bl6m3 untreated mice, 3 DSS-treated mice1
CACGSE31106ICRm3 untreated mice, 3 AOM/DSS-treated mice3
GSE42768C57Bl6f3 untreated mice, 3 DSS-treated mice3
GSE56050Lrg5-lacZm2 untreated mice, 2 AOM/DSS-treated mice3
GSE64658C57Bl6f6 untreated mice, 6 AOM/DSS-treated mice3
HumansUlcerative colitis (UC)GSE73661f, m12 healthy samples, 67 UC samples
GSE75214f, m11 healthy samples, 74 UC samples
GSE87466f, m21 healthy samples, 27 UC samples
GSE92415f, m21 healthy samples, 87 UC samples
Crohn’s disease (CD)GSE9686f, m8 healthy samples, 11 CD samples
GSE10616f, m26 healthy samples, 18 CD samples
GSE16879f, m6 healthy samples, 18 CD samples
GSE75214f, m11 healthy samples, 59 CD samples
Table 2. Primers used in the study.
Table 2. Primers used in the study.
GeneTypeSequence
Adam8Forward5′-TATGCAACCACAAGAGGGAG-3′
Probe5′-((5,6)-FAM)-TCATCTGATACATCTGCCAGCCGC-3′–BHQ1
Reverse5′-ACCAAGACCACAACCACAC-3′
C3Forward5′-GTTTATTCCTTCATTTCGCCTGG-3′
Probe5′-((5,6)-FAM)-ACACCCTGATTGGAGCTAGTGGC-3′–BHQ1
Reverse5′-GATGGTTATCTCTTGGGTCACC-3′
Timp1Forward5′-CTCAAAGACCTATAGTGCTGGC-3′
Probe5′-((5,6)-FAM)-ACTCACTGTTTGTGGACGGATCAGG-3′–BHQ1
Reverse5′-CAAAGTGACGGCTCTGGTAG-3′
TyrobpForward5′-GGTGACTTGGTGTTGACTCTG-3′
Probe5′-((5,6)-FAM)-CCTTCCGCTGTCCCTTGACCTC-3′–BHQ1
Reverse5′-GACCCTGAAGCTCCTGATAAG-3′
Mmp3Forward5′-TGCATATGAGGTTACTAACAGAGAC-3′
Probe((5,6)-FAM)-5′-AATCAGTTCTGGGCTATACGAGGGC-3′-BHQ1
Reverse5′-CAGGGTGTGAATGCTTTTAGG-3′
Mmp7Forward5′-CATAATTGGCTTCGCAAGGAG-3′
Probe((5,6)-FAM)-5′-TACTGGACTGATGGTGAGGACGCA-3′-BHQ1
Reverse5′-CAAATTCATGGGTGGCAGC-3′
Mmp9Forward5′-ACCTGAAAACCTCCAACCTC-3′
Probe((5,6)-FAM)-5′-TAGCGGTACAAGTATGCCTCTGCC-3′-BHQ1
Reverse5′-TCGAATGGCCTTTAGTGTCTG-3′
Mmp13Forward5′-GATTATCCCCGCCTCATAGAAG-3′
Probe((5,6)-FAM)-5′-CAGCATCTACTTTGTTGCCAATTCCAGG-3′-BHQ1
Reverse5′-CCCACCCCATACATCTGAAAG-3′
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Markov, A.V.; Savin, I.A.; Zenkova, M.A.; Sen’kova, A.V. Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models. Int. J. Mol. Sci. 2023, 24, 4311. https://doi.org/10.3390/ijms24054311

AMA Style

Markov AV, Savin IA, Zenkova MA, Sen’kova AV. Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models. International Journal of Molecular Sciences. 2023; 24(5):4311. https://doi.org/10.3390/ijms24054311

Chicago/Turabian Style

Markov, Andrey V., Innokenty A. Savin, Marina A. Zenkova, and Aleksandra V. Sen’kova. 2023. "Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models" International Journal of Molecular Sciences 24, no. 5: 4311. https://doi.org/10.3390/ijms24054311

APA Style

Markov, A. V., Savin, I. A., Zenkova, M. A., & Sen’kova, A. V. (2023). Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models. International Journal of Molecular Sciences, 24(5), 4311. https://doi.org/10.3390/ijms24054311

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop