- State Key Laboratory of Natural Medicines, Research Center of Biostatistics and Computational Pharmacy, China Pharmaceutical University, Nanjing, China
Background: Colorectal cancer is the fourth most deadly cancer worldwide. Although current treatment regimens have prolonged the survival of patients, the prognosis is still unsatisfactory. Inflammation and lncRNAs are closely related to tumor occurrence and development in CRC. Therefore, it is necessary to establish a new prognostic signature based on inflammation-related lncRNAs to improve the prognosis of patients with CRC.
Methods: LASSO-penalized Cox analysis was performed to construct a prognostic signature. Kaplan-Meier curves were used for survival analysis and ROC curves were used to measure the performance of the signature. Functional enrichment analysis was conducted to reveal the biological significance of the signature. The R package “maftool” and GISTIC2.0 algorithm were performed for analysis and visualization of genomic variations. The R package “pRRophetic”, CMap analysis and submap analysis were performed to predict response to chemotherapy and immunotherapy.
Results: An effective and independent prognostic signature, IRLncSig, was constructed based on sixteen inflammation-related lncRNAs. The IRLncSig was proved to be an independent prognostic indicator in CRC and was superior to clinical variables and the other four published signatures. The nomograms were constructed based on inflammation-related lncRNAs and detected by calibration curves. All samples were classified into two groups according to the median value, and we found frequent mutations of the TP53 gene in the high-risk group. We also found some significantly amplificated regions in the high-risk group, 8q24.3, 20q12, 8q22.3, and 20q13.2, which may regulate the inflammatory activity of cancer cells in CRC. Finally, we identified chemotherapeutic agents for high-risk patients and found that these patients were more likely to respond to immunotherapy, especially anti-CTLA4 therapy.
Conclusion: In short, we constructed a new signature based on sixteen inflammation-related lncRNAs to improve the outcomes of patients in CRC. Our findings have proved that the IRLncSig can be used as an effective and independent marker for predicting the survival of patients with CRC.
Introduction
Colorectal cancer (CRC) is the fourth deadliest cancer worldwide, its incidence is closely related to lifestyle, heredity and environmental factors, and the incidence and death rate in women is approximately 25% lower than those in men (Dekker et al., 2019). CRC accounts for approximately 10% of cancer-related deaths every year (Bray et al., 2018). And it is prone to metastasis and invasion, especially in advanced patients (Li et al., 2020). CRC can be divided into rectal adenocarcinoma and colon adenocarcinoma based on anatomic location (Zuo et al., 2019). The treatment of colon cancer and rectal cancer is mainly surgical resection, and rectal surgery is more difficult than colon surgery. For advanced patients with CRC, oxaliplatin in combination with fluorouracil is the standard treatment (Xie et al., 2020). Although the current treatment strategies have prolonged the overall survival of patients, the 5-years survival rate of patients remains low, at just over 50% (Frampton and Houlston, 2016).
Long noncoding RNAs (lncRNAs) are involved in a great diversity of biological processes, and various mechanisms of lncRNA transcriptional regulation of gene expression have been demonstrated (Navarro et al., 2006; Rinn et al., 2007; Gupta et al., 2010). Increasing evidence suggests that the abnormal expression of lncRNAs is closely related to tumor progression. For example, MIR17HG, SNHG14 and VIM-AS1 can promote the tumorigenesis and metastasis of tumor cells in CRC (Rezanejad et al., 2018; Xu et al., 2019a; Wu et al., 2019). Furthermore, inflammation is a marker of cancer and can be caused by diverse factors, including infection, environmental, cell senescence and obesity (Terzi et al., 2010). Inflammation is double-sided in cancer: on the one hand, cytotoxic T lymphocytes, which can target specific cancer cells, and regulatory T cells, which can inhibit nonspecific inflammation, may contribute to the antitumor response (Long et al., 2017). On the other hand, inflammatory cells and the chemokines and cytokines they produce can affect the biological processes of cells in the tumor microenvironment, thus promoting tumorigenesis, development, malignant transformation, invasion and metastasis of tumors (Balkwill and Mantovani, 2001; Werb and Coussens, 2006; Grivennikov et al., 2010; Greten and Grivennikov, 2019). Sustained use of low-dose anti-inflammatory drugs can slow the tumor progression in the long term, which also illustrates the tumor-promoting role of inflammation in CRC (Friis et al., 2015). In addition, lncRNAs are closely related to immunity and inflammation and can regulate inflammation and immune-related signaling pathways through a variety of mechanisms (Xu et al., 2019b).
Here, we constructed a new prognostic signature, IRLncSig and demonstrated its stability and reliability in validation cohorts. Functional enrichment analysis, somatic mutation and copy number alterations analysis were performed to reveal the biological significance of the IRLncSig. Then, we identified potential therapeutic agents for high-risk patients and found that high-risk group may respond to immunotherapy, especially anti-CTLA4 therapy. Overall, the IRLncSig was a robust and independent marker for predicting the outcomes of patients, which provides a new basis for improving outcomes in CRC.
Materials and methods
Study cohorts and data preprocessing
The data used in this study mainly include three publicly available datasets, TCGA-COADREAD (including TCGA-COAD and TCGA-READ) for training and testing, GSE38832, GSE39582 and GSE72970 for validation. For the TCGA-COADREAD cohort, the expression profiles and clinical information, including age, sex, pathological stage, status, and overall survival (OS), were downloaded from TCGA (https://portal.gdc.cancer.gov/) using the R package “TCGAbiolinks”, and the expression value was converted into TPM values with log2 (TPM+1) for subsequent analysis. Samples lacking clinical information were excluded. In addition, the cohort was divided into the training set and testing set at a ratio of 7:3 using the R package “caret”. For the GSE39582, GSE38832 and GSE72970 cohorts, the expression profiles and clinical information were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). Similarly, samples lacking clinical information were removed. Detailed sample information of these datasets is summarized (Table 1).
Construction and validation of the IRLncSig
Genes associated with inflammation were downloaded from the gene database in NCBI (www.ncbi.nlm.nih.gov/gene) (Brown et al., 2015; Choe et al., 2021). First, we conducted Pearson correlation analysis to identify lncRNAs associated with inflammation (Pearson correlation coefficient: r > 0.3). Subsequently, univariate Cox regression analysis was performed to screen out lncRNAs with prognostic values (p < 0.01). To minimize the risk of overfitting, we used a multivariate Cox regression model with the least absolute shrinkage and selection operator (LASSO) with tenfold cross-validation to tune the optimal value of penalty parameter 𝜆 using the R package “glmnet”. Finally, a new signature, IRLncSig, was constructed based on sixteen inflammation-related lncRNAs. The IRLncSig was defined as follows:
where
Construction and evaluation of nomograms and calibration curves
A nomogram model was constructed according to 16 inflammation-related lncRNAs to predict 1-,3-, and 5-years OS using the R package “rms”. Then, we use the calibration curves of 1-,3-, and 5-years to verify the performance of nomogram models by bootstrap method with 1,000 resamples.
Independence of prognostic effect of the IRLncSig
To evaluate the independence of the IRLncSig, we conducted multivariate regression analysis on age, sex, pathologic stage and the IRLncSig. We also performed stratification analysis by pathologic stage and age respectively. For the pathological stage, we divided samples into two groups, “stage I/II” and “stage III/IV” after removing samples with missing values. For age, we classified samples into two groups based on the median value. Kaplan–Meier curves and Log-rank test were used for evaluating the differences in OS using the R package “survival”.
Biological function and pathway analysis
To further understand the biological significance of the IRLncSig, we performed gene set enrichment analysis (GSEA) through the R package “clusterProfiler” (Yu et al., 2012). The functional annotation gene set (msigdb.v7.2. symbols.gmt) was downloaded from the MSigDB database (http://www.gsea-msigdb.org/gsea/msigdb). First, we obtained differentially expressed genes using the R package “limma” and sorted these genes by log2FoldChange (Ritchie et al., 2015). Then, we conducted GSEA based on the pre-ranked gene list (Subramanian et al., 2005). In addition, ten oncogenic pathways and their gene signatures were obtained from published literature (Sanchez-Vega et al., 2018). Gene set variation analysis (GSVA) was used to calculate the enrichment score of carcinogenic pathways.
Somatic mutation and copy number variation analysis
Somatic mutation data were obtained from cBioPortal (https://www.cbioportal.org/) and the R package “maftools” was used to analyze and visualize (Mayakonda et al., 2018). And the MutSigCV method was performed to identify driver genes in CRC (Lawrence et al., 2013).
Copy number variation data were obtained from Fire Browse (http://firebrowse.org/) and the Genomic Identification of Significant Targets in Cancer (GISTIC2.0, GenePattern) algorithm was performed to analyze and compare genomic alterations between two risk groups, including gains and losses (Mermel et al., 2011). The confidence level was set to 0.95 (Yang et al., 2020a; Luo et al., 2020).
Drug sensitivity analysis
We predicted the sensitivity of each patient to chemotherapeutic agents by the R package “pRRophetic” (Paul et al., 2014). All 138 chemotherapy drugs were included in this package. The estimated IC50 value for each patient treated with a specific chemotherapy drug was obtained through the function “pRRopheticPredict”. Then, we screened out drugs with higher sensitivity for high-risk group patients. In addition, we used CMap analysis (https://clue.io) to determine agents, with which gene expression value increased in high-risk patients but decreased by treatment (scores <0) (Chen et al., 2020). Finally, we set the intersection of agents obtained in the two steps to determine the final potential treatment agents for patients with high risk.
Characterization of immune microenvironment and prediction of immunotherapy response
To further reveal the differences in immune microenvironment between the two groups, we first estimated the abundance of 22 immune cells and 2 stromal cells using single sample gene set enrichment analysis (ssGSEA). We also obtained 10 immune checkpoints from published literature and compared the difference in expression between the two groups (Yang et al., 2020b). In addition, the submap algorithm was used to evaluate the response of patients to immunotherapy based on the dataset of 47 patients with cutaneous melanoma (Hoshida et al., 2007). The method was performed by the Submap module in GenePattern.
Performance comparison of the IRLncSig with other clinical variables and existing prognostic signatures
ROC curves were used to estimate the 1-year, 3-years and 5-years survival. The clinical variables age and pathological stage were considered. In detail, the IRLncSig and age were treated as continuous variables, while the pathologic stage was coded as binary, “Stage I/II” and “Stage III/IV”. We also compared the IRLncSig with four published prognostic signatures in CRC: a 15-gene signature derived from Yang’s study (Yang et al., 2021), a 6-lncRNA signature derived from Li’s study (Li et al., 2020), a 6-gene signature derived from Dai’s study (Dai et al., 2020), and a 7-gene signature derived from Lu’s study (Lu et al., 2021).
Statistical analysis
All statistical analyses were conducted by R version 4.0.2. Statistical test methods used in this study mainly include the Wilcoxon test for two groups, the Kruskal-Wallis test for more than two groups, and the Chi-Square test or Fisher’s exact test for contingency table variables (Hazra and Gogtay, 2016b). All statistical tests were two-sided, and it was considered statistically significant with p < 0.05.
Results
Construction and validation of the IRLncSig based on inflammation-related lncRNAs
A detailed workflow for prognostic signature construction and analysis was developed (Supplementary Figure S1). We downloaded 2,826 inflammation-related genes from NCBI and finally obtained 2,380 inflammation-related genes after the intersection with expression profile data. We screened out 934 lncRNAs strongly correlated with inflammatory genes (Pearson correlation coefficient: r > 0.3) and 34 lncRNAs with prognostic values (p < 0.01; Supplementary Figure S2). Next, sixteen lncRNAs with prognostic values were identified by univariate and LASSO-penalized regression analysis (Table 2; Figures 1A,B). The sixteen inflammation-related lncRNAs, which were independently correlated with survival, were preserved as prognostic factors. Of them, thirteen lncRNAs were related to worse prognosis and presented as “risk” factors, while the remaining three lncRNAs were related to better prognosis and presented as “protective” factors (Figure 1C).Figures 1D,E show the expression levels of the sixteen lncRNAs in the TCGA and GSE39582 datasets, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g001.jpg)
FIGURE 1. Construction of the prognostic signature, IRLncSig (A) Feature selection using the LASSO-penalized regression by 10-fold cross-validation with minimum criteria. (B) LASSO coefficient profiles of the 34 inflammation-related lncRNAs with nonzero coefficients (C) Forest plot for identified sixteen lncRNAs. (D,E) Boxviolin plot for the expression level of the sixteen lncRNAs.
In each cohort, we calculated the risk score according to the coefficients and expression values of the sixteen lncRNAs and divided samples into high-risk and low-risk groups based on the median value. Survival analysis indicated that high-risk group had a worse prognosis in the TCGA training cohort (Figure 2A). We found similar results in TCGA testing, TCGA entire, GSE38832, GSE39582 and GSE72970 datasets (Figures 2B–F). The AUC value was 0.79 in the first year and 0.78 both in the third and fifth year in the TCGA training dataset (Figure 2G), and the AUC values in other validation datasets are shown in Figures 2H–L. In addition, we found that the mortality was elevated with the increase of the risk score (Supplementary Figure S3). The “risk” lncRNAs were highly expressed in high-risk group and the protective lncRNAs were highly expressed in the low-risk group (Supplementary Figure S3). These results indicated that the IRLncSig could act as an effective prognostic marker in CRC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g002.jpg)
FIGURE 2. Validation of the performance of the IRLncSig for prognostic prediction (A–F) Survival analysis for two groups. The p values were calculated by the Log-rank test. (G–L) ROC curves for predicting 1-, 3- and 5-years OS.
Construction and evaluation of nomograms and calibration curves
The nomogram comprising the sixteen lncRNAs was fabricated to estimate the 1-year, 3-years and 5-years OS in the TCGA dataset (Figure 3A). The calibration curve for the observed versus predicted probability of the 1-year, 3-years and 5-years OS revealed consistency (Figure 3D). Similar results could be observed in GSE39582 and GSE38832 validation cohorts Figures 3B,C; Figures 3E,F). These results suggested that the nomogram model based on the sixteen lncRNAs could help predict the prognosis of patients in CRC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g003.jpg)
FIGURE 3. Construction and evaluation of prognostic nomograms and calibration curves (A–C) The nomogram predicts the probability of 1-, 3-, and 5-years OS. (D–F) The calibration plot of the nomogram predicts the probability of 1-, 3-, and 5-years OS.
Independence of prognostic effects of the IRLncSig
To evaluate the independence of the IRLncSig, we conducted multivariate regression analysis on age, sex, pathologic stage and the IRLncSig. The results showed that the IRLncSig could act as an independent prognostic signature in all cohorts (Figure 4A). In addition to the IRLncSig, the pathologic stage and age were observed to be highly significant (Figure 4A). Therefore, we performed stratification analysis by pathologic stage and age respectively. For the pathological stage, the results indicated that there were significant differences in OS between two groups “stage I/II” and “stage III/IV” in all datasets Figures 4B–E). Similar results were found in the stratification analysis of age (Supplementary Figure S4). These results all indicated that the IRLncSig could act as an independent and effective marker in CRC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g004.jpg)
FIGURE 4. Independence of prognostic effect of the IRLncSig (A) Forest plot for the IRLncSig and other clinical factors by multivariate regression analysis. (B–E) Stratification analysis by pathological stage. The p values were calculated by the Log-rank test.
Biological function and pathway analysis
To further understand the biological significance of the IRLncSig, differential expression analysis was first performed. Next, GSEA was performed based on Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Hallmark and Reactome gene sets. The results revealed that epithelial-mesenchymal transition (EMT), angiogenesis, inflammatory response, hypoxia, WNT and KRAS signaling pathways were enriched in the high-risk group (Figure 5A). Oxidative phosphorylation, MYC signaling, cell cycle, fatty acid metabolism, drug metabolism and TCA cycle were enriched in the low-risk group (Figure 5B). Next, GSVA was used to estimate the abundance of up-regulated, down-regulated, and carcinogenic pathways. We found that the up-regulated pathways exhibited a higher GSVA score in the high-risk group, while the down-regulated pathways exhibited the opposite (Figures 5C,D). In addition, HIPPO, NOTCH, RTK/RAS and WNT carcinogenic pathways had higher GSVA score in the high-risk group, while cell cycle and TP53 signaling had higher GSVA score in the low-risk group (Figures 5C,D). The results were similar in the GSE39582 validation set (Supplementary Figure S5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g005.jpg)
FIGURE 5. Function enrichment analysis in the TCGA cohort (A,B) Gene set enrichment analysis based on GO, KEGG, Hallmark and Reactome. (C) Heatmap for the GSVA score of upregulated, downregulated and oncogenic pathways (D) Boxplot for the GSVA score of upregulated, downregulated and oncogenic pathways between two subgroups. The p values were calculated by Wilcoxon test (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0,001).
Characterization of somatic mutation and copy number variations
We first identified driver genes in CRC by MutSigCV (q < 0.05; Supplementary Figure S6A). Then, we found that the mutation rate of gene TP53 was significantly higher in the high-risk group (Supplementary Figure S6B). To further verify the prognostic independence of the IRLncSig, we divided samples into four subgroups based on TP53 mutation status and risk group: Mutated in HRisk, Mutated in LRisk, Wild in HRisk and Wild in LRisk. Survival analysis showed that the OS of both mutated-type and wild-type in the high-risk group was significantly lower (Supplementary Figure S6C). These results again demonstrated the independence of the IRLncSig in predicting outcomes of patients in CRC.
GISTIC2.0 was used to analyze the copy number variation (Mermel et al., 2011). First, we compared genome alterations between two groups, including the fraction of genome altered (FGA), the fraction of genome lost (FGL) and the fraction of genome gained (FGG). Results showed that high-risk group exhibited more genome alterations, including FGA, FGL, and FGG (Figure 6A). We also found seven significantly amplificated regions (8q22.3, 8q24.3, 12p13.32, 17q11.1, 17q24.2, 20q12, 20q13.2) and five significantly deleted regions (4q31.3, 5q21.3, 8p21.3, 8p23.1, 15q12) in the high-risk group (Figure 5B). We further analyzed the correlation between the copy number of these amplified genes and their expression values and identified eleven inflammatory genes in the regions 8q24.3, 20q12, 8q22.3 and 20q13.2 with a strongly positive correlation (Table 3; Supplementary Figure S6D). These genes were highly expressed in the high-risk group except for the genes PLCG1, NFATC2, TOP1 and CYP24A1 (Supplementary Figure S6E). These results suggested that the four regions may regulate the inflammatory activity of cancer cells, further clarifying that alterations in specific chromosomal regions may influence the heterogeneity of CRC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g006.jpg)
FIGURE 6. Characterization of copy number variation in the entire TCGA cohort (A) The difference in genomic alterations between the two groups. The p values were calculated by Wilcoxon test (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0,001) (B) GISTIC cytoband of copy number alteration for all samples. (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0,001).
Identification of potential therapeutic agents for patients in CRC
We first identified five drugs with higher sensitivity for high-risk group patients from 138 chemotherapy drugs Figures 7A,B). Then, the correlation of IC50 values and target expression level with the IRLncSig for five drugs was analyzed (Supplementary Figure S7). We further tested the efficacy of these drugs by CMap analysis based on differentially expressed genes. The results showed that PD.173,074 (CMap score = −91.88) and Docetaxel (CMap score = −91.55) might be potential therapeutic drugs for patients with high risk (Figure 7C). PD-173074 is a fibroblast growth factor receptor (FGFR) inhibitor. FGFR1 is overexpressed and is considered a therapeutic target in CRC (Jang, 2005; Göke et al., 2013). Studies have also shown that its overexpression is associated with liver metastasis of CRC (Kwak et al., 2015). Docetaxel is a tubulin polymerization inhibitor, and its target is BCL-2, which is an oncogene and can inhibit apoptosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g007.jpg)
FIGURE 7. Identification of potential therapeutic agents for high-risk patients (A,B) Five agents with higher drug sensitivity in the high-risk group. The p values were calculated by the Wilcoxon test. (C) CMap analysis of five drugs with potential therapeutic agents for high risk patients.
Characterization of immune microenvironment between two risk groups
We first characterized the cellular interactions in the immune microenvironment, including 22 immune cells and 2 stromal cells, and the results showed that there was a strong positive correlation between these cells (Supplementary Figure S8A). CD8 T cells (HR: 0.84, 95% CI: 0.72–0.99, p = 0.038) and activated memory CD4 T cells (HR: 0.80, 95% CI: 0.67–0.95, p = 0.013) exhibited prognostic value in CRC (Supplementary Figure S8B). For activated memory CD4 T cells, the group with higher immune infiltration had a better prognosis (p = 0.004; Supplementary Figure S8C). However, there was no significant difference for CD8 T cells (p = 0.125; Supplementary Figure S8D). We also analyzed the correlation between immune cell infiltration and inflammation-related lncRNAs (Supplementary Figure S8E).
To further explore the response to immunotherapy of patients, we analyzed the infiltration levels of immune cells and stromal cells in the immune microenvironment. We found that the abundance of naive B cells, activated dendritic cells, macrophages M0, macrophages M2, monocytes, naive CD4 T cells, regulatory T cells and gamma delta T cells were higher in the high-risk group, while the abundance of memory B cells, eosinophils, NK cells, activated memory CD4 T cells and resting memory CD4 T cells were higher in the low-risk group (Figure 8A). We also found that high-risk group exhibited higher expression values of immune checkpoints, including CCL2, CD276, CD4, CXCR4, IL1A, and TGFB1 (Figure 8B). In addition, the submap analysis showed that high-risk group was more likely to respond to anti-CTLA-4 immunotherapy (nominal p = 0.02; Figure 8C). However, no difference was found in anti-PD1 immunotherapy, which may need further study. These results were verified in the GSE39582 cohort (Supplementary Figure S9). Collectively, high-risk group exhibited higher immune cell infiltration and expression of immune checkpoints, suggesting that patients in the high-risk group may be responsive to immunotherapy, especially anti-CTLA4 therapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/955240/fgene-13-955240-HTML/image_m/fgene-13-955240-g008.jpg)
FIGURE 8. Immune infiltration and the response to immunotherapy in TCGA cohort (A) Boxplot for the abundance of 22 immune cells and 2 stromal cells estimated by ssGSEA. The p values were calculated by Wilcoxon test (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0,001) (B) Boxplot for the expression level of 10 immune checkpoints. The p values were calculated by Wilcoxon test (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0,001). (C) Submap analysis for predicting the response to immunotherapy.
Performance comparison of the IRLncSig with other clinical variables and existing prognostic signatures
We compared the IRLncSig with four recently published prognostic signatures in CRC: a 15-gene signature derived from Yang’s study (Yang et al., 2021), a 6-lncRNA signature derived from Li’s study (Li et al., 2020), a 6-gene signature derived from Dai’s study (Dai et al., 2020), and a 7-gene signature derived from Lu’s study (Lu et al., 2021). The results showed that the AUC values at 1-year, 3-years and 5-years were 0.78, 0.78 and 0.79, respectively, which were all higher than those of the other four signatures (Figure 9). The results proved that the IRLncSig had better predictive performance than the other four published signatures in CRC. In addition, the performance of the signature was better than that of the clinical variables, age and pathologic stage (Supplementary Figure S10).
Finally, we associated the IRLncSig with the CMS subtypes established by previous studies (Guinney et al., 2015). The CMS4 subtypes had the worst prognosis and the highest risk score (Supplementary Figures S11A,B), and high-risk group exhibited a higher proportion of CMS4 subtype (p < 0.001; Supplementary Figure S11C). We also applied the signature to pan-cancer and found that it showed prognostic value in liver hepatocellular carcinoma (LIHC; HR: 1.47, 95% CI: 1.22–1.76, p < 0.001), uterine corpus endometrial carcinoma (UCEC; HR: 1.38, 95% CI: 1.14–1.67, p < 0.001), kidney renal clear cell carcinoma (KIRC; HR: 1.25, 95% CI: 1.11–1.41, p < 0.001), head and neck squamous cell carcinoma (HNSC; HR: 1.20, 95% CI: 1.05–1.36, p = 0.006), brain lower grade glioma (LGG; HR: 0.81, 95% CI: 0.68–0.97, p = 0.020), uveal melanoma (UVM; HR: 0.59, 95% CI: 0.40–0.89, p = 0.012) (Supplementary Figure S12).
Discussion
CRC is a common cancer of the digestive tract, with morbidity and mortality among the top five in the world (Kuipers et al., 2015; Qaderi et al., 2020). More than 1.4 million people have been diagnosed as CRC, and more than 0.5 million people have died of this tumor in a year (Jeffery et al., 2007; Siegel et al., 2020). Approximately 30% of patients have tumor recurrence within 3 years after surgery and adjuvant chemotherapy (Sargent et al., 2009). Many studies have shown that the construction of prognostic signatures could accelerate the prognostic assessment of patients in CRC (Huang et al., 2020). In this study, we constructed a robust prognostic signature, IRLncSig, based on sixteen inflammation-related lncRNAs to improve the prognosis of patients with CRC. The predictive performance of the IRLncSig was verified on TCGA, GSE38832, GSE39582, and GSE72970 cohorts. Then, samples were divided into two risk groups based on the median value in each dataset and Kaplan–Meier curves exhibited that patients in the low-risk group had a better prognosis. The performance of the IRLncSig was also proved to outperform traditional clinical predictors and four published signatures in CRC. Subgroup analyses showed that the signature could predict outcomes in different subgroups of patients. These results indicated that the IRLncSig was an independent and effective marker in predicting the prognosis of patients with CRC.
To better understand the biological significance of the IRLncSig, we conducted enrichment analysis and found that high-risk group mainly enriched in EMT, angiogenesis, inflammatory response, hypoxia, WNT signaling and KRAS signaling pathways. EMT can promote tumor progression, including CRC, especially when it involves invasion and metastasis (Arias, 2001; Bates and Mercurio, 2005). Angiogenesis is a hallmark of cancer and it is associated with tumor progression and metastasis in some cancers, such as CRC, skin melanoma, breast cancer, prostate cancer and lung cancer (Yuji et al., 1996). Moreover, carcinogenic pathways, including HIPPO, NOTCH and RTK/RAS were enriched in the high-risk group. Low-risk group enriched mainly in oxidative phosphorylation, MYC signaling, cell cycle, fatty acid metabolism, drug metabolism and TCA cycle. Restoration of oxidative phosphorylation levels is one of the characteristics of tumor development and progression (Maiuri and Kroemer, 2015). The tricarboxylic acid cycle, also known as the citric acid cycle, is a common metabolic pathway in the human body. It is the hub of the metabolic links between the three nutrients, including carbohydrates, lipids and amino acids. The dysfunction of TCA is one of the causes of human disease and tumor formation (Brière et al., 2007). Several studies have shown that inflammation has been implicated in EMT (Suarez-Carmona et al., 2017), angiogenesis (Ping, 2012; Aguilar-Cazares et al., 2019), hypoxia (Mamlouk and Wielockx, 2013), WNT signaling (Yang et al., 2008), KRAS signaling (Kitajima et al., 2016), MYC signaling (Sipos et al., 2016) and lipid metabolism (Zuo et al., 2020). Therefore, it is speculated that the IRLncSig constructed based on inflammation-related lncRNAs may impact the activity of these pathways through inflammatory responses.
Gene mutations are implicated in cancer. Among them, TP53 is an important tumor suppressor gene and the wild-type causes apoptosis of cancer cells preventing carcinogenesis, and TP53 mutation will increase the probability of carcinogenesis (Guimaraes and Hainaut, 2002). Pan-cancer studies have shown that TP53 mutations have been found in a variety of malignancies and the mutation rate is 43.2% in CRC (Olivier et al., 2010). A large number of studies have shown that TP53 mutation is related to poor survival in CRC (Iacopetta, 2003). Here, we found that the TP53 gene was frequently mutated in the high-risk group. We also found that the prognostic effect of the IRLncSig was independent of TP53 mutation status, which again demonstrated the independence of the IRLncSig in predicting outcomes of patients in CRC.
Copy number variation is an important component of genome alterations, which affects gene expression and the activity of signaling pathways. Here, we found that high-risk group had more genome alterations. We also found seven significantly amplificated regions in the high-risk group, 8q22.3, 8q24.3, 12p13.32, 17q11.1, 17q24.2, 20q12, and 20q13.2. The expression of gene SCRIB in region 8q24.3 was positively correlated with its copy number (Spearman’s r = 0.522). SCRIB gene is closely related to tumors and is considered to be a key factor of tumor development and metastasis (Zen et al., 2009; Royer and Lu, 2011). It is associated with proliferation, apoptosis, EMT and poor prognosis in CRC (Shen et al., 2021). SCRIB also interacts with the HIPPO signaling pathway in CRC (Enomoto and Igaki, 2011). These results clarified that alterations in specific chromosomal regions may affect the expression of genes, and further influence a series of signaling pathways and biological processes in CRC.
Next, we identified five agents, PD.173,074, Pazopanib, Docetaxel, KU-55933 and Embelin, which may have potential efficacy for high-risk patients in CRC. Of the five agents, PD.173,074 (CMap score = −91.88) and Docetaxel (CMap score = −91.55) were probably the most promising agents. PD-173074 is a fibroblast growth factor receptor (FGFR) inhibitor. The FGFR inhibitors have shown significant antitumor activity against multiple tumor cell lines, including stomach, lung, bladder, endometrium and breast (Rusnati MP and Presta, 2007). Many FGFR inhibitors have been developed to treat various cancers; for example, brivanib (Ayers et al., 2007), AZD4547 (Gavine et al., 2012), E-3810 (Bello et al., 2011), NP603 (Kammasud et al., 2007), etc. Docetaxel is a tubulin polymerization inhibitor. Tubulin polymerization inhibitors inhibit spindle formation by inducing depolymerization of microtubules, block cells in phase M, and eventually induce cell apoptosis (Morris and Fornier, 2008). Therefore, the two drugs are considered to have the most promising therapeutic potential for patients with high risk in CRC.
In addition to the tumorigenic effect, inflammation can also affect the immune system, thereby enhancing the response to chemotherapy (Apetoh et al., 2008). However, inflammation can blunt the beneficial effects of chemotherapy in some cases (Ammirante et al., 2010). Therefore, we explored the relationship between the IRLncSig and immunotherapy. We found high-risk group exhibited higher immune cell infiltration and expression value of immune checkpoints, suggesting that high-risk patients may be more likely to respond to immunotherapy. The submap analysis further confirmed this conclusion, which may provide some help for the development of clinical treatment strategies for patients with CRC.
There are some limitations in this study. Firstly, less clinical information was considered in this study considering the availability of data. It needs to verify the independence of the IRLncSig in a larger cohort containing more clinical information. Secondly, this study conducted a series of bioinformatics analyses based on public databases and may need further biological experimental verification. Finally, this study is a prospective study, and prospective studies may be needed to further examine the prognostic performance of the IRLncSig.
Conclusion
Collectively, we established a new signature, IRLncSig, based on sixteen inflammation-related lncRNAs. Our analysis suggested that the IRLncSig is an independent and effective prognostic marker and superior to the other four existing prognostic signatures in CRC. Functional enrichment analysis, somatic mutation and copy number alterations were conducted to reveal the biological significance of the IRLncSig. We performed drug sensitivity predictions and identified chemotherapeutic agents that might be effective for high-risk group patients. We also found that high-risk patients were more likely to respond to immunotherapy, especially anti-CTLA4 therapy. These findings could provide some help for the development of clinical treatment strategies in CRC.
Data availability statement
All data were downloaded from public databases. The TCGA-COADREAD cohort (including TCGA-COAD and TCGA-READ) was downloaded from the TCGA database (https://portal.gdc.cancer.gov/). The GSE39582 cohort, GSE38832 cohort and GSE72970 cohort were downloaded from the GEO database (GSE39582: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39582; GSE38832: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE38832; GSE72970: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72970). In addition, no new cohorts were generated in this study.
Author contributions
FY and XL guided this work. MH participated in data collection and processing, program implementation, and paper writing. YY, YC, JZ, LX, and WC contributed to the statistical analysis. All authors were involved in this study.
Funding
This study was supported by the National Natural Science Foundation of China [No.81973145, No. 82273735], and the Key R&D Program of Jiangsu Province (Social Development) (BE2020694).
Acknowledgments
The analysis work is all based on public databases. We thank the TCGA and GEO databases for providing the data analyzed for this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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/fgene.2022.955240/full#supplementary-material
Reference
Aguilar-Cazares, D., Chavez-Dominguez, R., Carlos-Reyes, A., Lopez-Camarillo, C., Hernadez de la Cruz, O.N., and Lopez-Gonzalez, J.S. (2019). Contribution of angiogenesis to inflammation and cancer. Front. Oncol. 9, 1399. doi:10.3389/fonc.2019.01399
Arias, A. M. (2001). Epithelial mesenchymal interactions in cancer and development. Cell 105 (4), 425–431. doi:10.1016/s0092-8674(01)00365-8
Ammirante, M., Luo, J. L., Grivennikov, S., Nedospasov, S., and Karin, M. (2010). B-cell-derived lymphotoxin promotes castration-resistant prostate cancer. Nature 464 (7300), 302–305. doi:10.1038/nature08782
Apetoh, L., Ghiringhelli, F., Kroemer, G., and Zitvogel, L. (2008). Immunological aspects of cancer chemotherapy. Nat. Rev. Immunol. 8 (1), 59–73. doi:10.1038/nri2216
Ayers, M., Fargnoli, J., Lewin, A., Wu, Q., and Platero, J. S. (2007). Discovery and validation of biomarkers that respond to treatment with brivanib alaninate, a small-molecule VEGFR-2/FGFR-1 antagonist. Cancer Res. 67 (14), 6899–6906. doi:10.1158/0008-5472.CAN-06-4555
Balkwill, F., and Mantovani, A. (2001). Inflammation and cancer: Back to virchow? Lancet 357, 539–545. doi:10.1016/S0140-6736(00)04046-0
Bates, R. C., and Mercurio, A. (2005). The epithelial-mesenchymal transition (EMT) and colorectal cancer progression. Cancer Biol. Ther. 4 (4), 365–370. doi:10.4161/cbt.4.4.1655
Bello, E., Colella, G., Scarlato, V., Oliva, P., Berndt, A., Valbusa, G., et al. (2011). E-3810 is a potent dual inhibitor of VEGFR and FGFR that exerts antitumor activity in multiple preclinical models. Cancer Res. 71 (4), 1396–1405. doi:10.1158/0008-5472.CAN-10-2700
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin. 68 (8), 394–424. doi:10.3322/caac.21492
Brière, J., Favier, J., Gimenez-Roqueplo, A. P., and Rustin, P. (2007). Tricarboxylic acid cycle dysfunction as a cause of human diseases and tumor formation. Am. J. Physiol. Cell Physiol. 291 (6), C1114–C1120. doi:10.1152/ajpcell.00216.2006
Brown, G. R., Vichet, H., Katz, K. S., Michael, O., Craig, W., Olga, E., et al. (2015). Gene: A gene-centered information resource at NCBI. Nucleic Acids Res. 43 (D1), 36–42. doi:10.1093/nar/gku1055
Chen, Y., Xiaowen, H., Yan, L., Junfei, C., Yuanyuan, L., and Shixue, D. (2020). Prognosis and personalized treatment prediction in TP53-mutant hepatocellular carcinoma: An in silico strategy towards precision oncology. Brief. Bioinform. 22, bbaa164. doi:10.1093/bib/bbaa164
Choe, E., Lee, S., Kim, S. Y., Shivakumar, M., Kim, D., Chai, Y. J., et al. (2021). Prognostic effect of inflammatory genes on stage I-III colorectal cancer—integrative analysis of TCGA data. Cancers 13 (4), 751. doi:10.3390/cancers13040751
Dai, S., Xu, S., Ye, Y., and Ding, K. (2020). Identification of an immune-related gene signature to improve prognosis prediction in colorectal cancer patients. Front. Genet. 11, 607009. doi:10.3389/fgene.2020.607009
Dekker, E., Tanis, P. J., Vleugels, J., Kasi, P. M., and Wallace, M. B. (2019). Colorectal cancer. Lancet 394 (10207), 1467–1480. doi:10.1016/S0140-6736(19)32319-0
Enomoto, M., and Igaki, T. (2011). Deciphering tumor-suppressor signaling in flies: Genetic link between scribble/dlg/lgl and the hippo pathways. J. Genet. Genomics 38 (10), 461–470. doi:10.1016/j.jgg.2011.09.005
Frampton, M., and Houlston, R. S. (2016). Modeling the prevention of colorectal cancer from the combined impact of host and behavioral risk factors. Genet. Med. 19 (3), 314–321. doi:10.1038/gim.2016.101
Friis, S., Riis, A. H., Erichsen, R., Baron, J. A., and Sørensen, H. (2015). Low-dose aspirin or nonsteroidal anti-inflammatory drug use and colorectal cancer risk: A population-based, case-control study. Ann. Intern. Med. 163 (5), 347–355. doi:10.7326/M15-0039
Göke, F., Göke, A., Mässenhausen, A. V., Franzen, A., Sharma, R., Kirsten, R., et al. (2013). Fibroblast growth factor receptor 1 as a putative therapy target in colorectal cancer. Digestion 88 (3), 172–181. doi:10.1159/000355018
Gavine, P. R., Mooney, L., Kilgour, E., Thomas, A. P., Klinowska, T., Beck, S., et al. (2012). AZD4547: An orally bioavailable, potent, and selective inhibitor of the fibroblast growth factor receptor tyrosine kinase family. Cancer Res. 72 (8), 2045–2056. doi:10.1158/0008-5472.CAN-11-3034
Greten, F. R., and Grivennikov, S.I. (2019). Inflammation and cancer: Triggers, mechanisms, and consequences. Immunity 51 (1), 27–41. doi:10.1016/j.immuni.2019.06.025
Grivennikov, S. I., Greten, F. R., and Karin, M. (2010). Immunity, inflammation, and cancer. Cell 140 (6), 883–899. doi:10.1016/j.cell.2010.01.025
Guimaraes, D. P., and Hainaut, P. (2002). TP53: A key gene in human cancer. Biochimie 84 (1), 83–93. doi:10.1016/s0300-9084(01)01356-6
Guinney, J., Dienstmann, R., Wang, X., Reyniès, A., Tejpar, S., Soneson, C., et al. (2015). The consensus molecular subtypes of colorectal cancer. Nat. Med. 21 (11), 1350–1356. doi:10.1038/nm.3967
Gupta, R. A., Shah, N., Wang, K.C., Kim, J., Horlings, H. M., Wong, D. J., et al. (2010). Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 464 (7291), 1071–1076. doi:10.1038/nature08975
Hazra, A., and Gogtay, N. (2016). Biostatistics series module 3: Comparing groups: Numerical variables. Indian J. dermatol. 61 (3), 251–260. doi:10.4103/0019-5154.182416
Hoshida, Y., Brunet, J. P., Tamayo, P., Golub, T. R., and Mesirov, J. P. (2007). Subclass mapping: Identifying common subtypes in independent disease data sets. PLoS ONE 2 (11), e1195. doi:10.1371/journal.pone.0001195
Huang, R., Zhou, L., Chi, Y., Wu, H., and Shi, L. (2020). LncRNA profile study reveals a seven-lncRNA signature predicts the prognosis of patients with colorectal cancer. Biomark. Res. 8 (1). doi:10.1186/s40364-020-00187-3
Iacopetta, B. (2003). TP53 mutation in colorectal cancer. Hum. Mutat. 21, 271–276. doi:10.1002/humu.10175
Jang, J.H. (2005). Reciprocal relationship in gene expression between FGFR1 and FGFR3: Implication for tumorigenesis. Oncogene 24 (5), 945–948. doi:10.1038/sj.onc.1208254
Jeffery, M., Hickey, B. E., and Hider, P. N. (2007). Follow-up strategies for patients treated for non-metastatic colorectal cancer. Cochrane Database Syst. Rev. 1 (1), CD002200. doi:10.1002/14651858.CD002200.pub2
Kammasud, N., Boonyarat, C., Tsunoda, S., Sakurai, H., Saiki, I., Grierson, D. S., et al. (2007). Novel inhibitor for fibroblast growth factor receptor tyrosine kinase. Bioorg. Med. Chem. Lett. 17 (17), 4812–4818. doi:10.1016/j.bmcl.2007.06.058
Kitajima, S., Thummalapalli, R., and Barbie, D. A. (2016). Inflammation as a driver and vulnerability of KRAS mediated oncogenesis. Semin. Cell Dev. Biol. 58, 127–135. doi:10.1016/j.semcdb.2016.06.009
Kuipers, E. J., Grady, W. M., Lieberman, D., Seufferlein, T., Watanabe, T., Boelens, P. G., et al. (2015). Colorectal cancer. Nat. Rev. Dis. Prim. 1, 15065. doi:10.1038/nrdp.2015.65
Kwak, Y., Nam, S. K., Seo, A. N., Kim, D. W., Kang, S. B., Kim, W. H., et al. (2015). Fibroblast growth factor receptor 1 gene copy number and mRNA expression in primary colorectal cancer and its clinicopathologic correlation. Pathobiology 82 (2), 76–83. doi:10.1159/000398807
Lawrence, M. S., Stojanov, P., Polak, P., Kryukov, G. V., Cibulskis, K., Sivachenko, A., et al. (2013). Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature 499 (7457), 214–218. doi:10.1038/nature12213
Li, S., Chen, S., Wang, B., Zhang, L., and Zhang, X. (2020). A robust 6-lncRNA prognostic signature for predicting the prognosis of patients with colorectal cancer metastasis. Front. Med. 7, 56. doi:10.3389/fmed.2020.00056
Long, A. G., Lundsmith, E. T., and Hamilton, K. E. (2017). Inflammation and colorectal cancer. Curr. Colorectal Cancer Rep. 13 (4), 341–351. doi:10.1007/s11888-017-0373-6
Lu, Y., Zhou, X., Liu, Z., Wang, W., and Fu, W. (2021). Prognosis prediction for colorectal cancer patients: A risk score based on the metabolic-related genes. Int. J. Med. Sci. 18 (3), 801–810. doi:10.7150/ijms.49576
Luo, H., Xu, X., Yang, J., Wang, K., Cai, H., Yang, P., et al. (2020). Genome-wide somatic copy number alteration analysis and database construction for cervical cancer. Mol. Genet. Genomics 295 (269), 765–773. doi:10.1007/s00438-019-01636-x
Maiuri, M. C., and Kroemer, G. (2015). Essential role for oxidative phosphorylation in cancer progression. Cell Metab. 21 (1), 11–12. doi:10.1016/j.cmet.2014.12.013
Mamlouk, S., and Wielockx, B. (2013). Hypoxia-inducible factors as key regulators of tumor inflammation. Int. J. Cancer 132 (12), 2721–2729. doi:10.1002/ijc.27901
Mayakonda, A., Lin, D. C., Assenov, Y., Plass, C., and Koeffler, H. P. (2018). Maftools: Efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 28 (11), 1747–1756. doi:10.1101/gr.239244.118
Mermel, C. H., Schumacher, S. E., Hill, B, Meyerson, M. L., Beroukhim, R., and Getz, G. (2011). GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 12 (4), R41–R. doi:10.1186/gb-2011-12-4-r41
Morris, P.G., and Fornier, M. N. (2008). Microtubule active agents: Beyond the taxane frontier. Clin. Cancer Res. 14 (22), 7167–7172. doi:10.1158/1078-0432.CCR-08-0169
Navarro, P., Page, D. R., Avner, P., and Rougeulle, C. (2006). Tsix-mediated epigenetic switch of a CTCF-flanked region of the Xist promoter determines the Xist transcription program. Genes Dev. 20 (20), 2787–2792. doi:10.1101/gad.389006
Olivier, M., Hollstein, M., and Hainaut, P. (2010). TP53 mutations in human cancers: Origins, consequences, and clinical use. Cold Spring Harb. Perspect. Biol. 2 (1), a001008. doi:10.1101/cshperspect.a001008
Paul, G., Nancy, C., Huang, R. S., and Barbour, J. D. (2014). pRRophetic: An R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. Plos One 9 (9), e107468. doi:10.1371/journal.pone.0107468
Qaderi, S. M., Swartjes, H., Custers, J., and Wilt, J. (2020). Health care provider and patient preparedness for alternative colorectal cancer follow-up; a review. Eur. J. Surg. Oncol. 46 (10), 1779–1788. doi:10.1016/j.ejso.2020.06.017
Rezanejad, H., Asadi, M. H., and Yaghoobi, M. M. (2018). Long noncoding RNA VIM-AS1 promotes colorectal cancer progression and metastasis by inducing EMT. Eur. J. Cell Biol. 97, 279–288. doi:10.1016/j.ejcb.2018.04.004
Rinn, J. L., Kertesz, M., Wang, J. K., Squazzo, S. L., Xiao, X., Brugmann, S. A., et al. (2007). Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell 129 (7), 1311–1323. doi:10.1016/j.cell.2007.05.022
Ritchie, M. E., Belinda, P., Wu, D., Hu, Y., Law, C. W., Wei, S., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Royer, C., and Lu, X. (2011). Epithelial cell polarity: A major gatekeeper against cancer?. Cell Death Differ. 18 (9), 1470–1477. doi:10.1038/cdd.2011.60
Rusnati Mp, M., and Presta, M. (2007). Fibroblast growth factors/fibroblast growth factor receptors as targets for the development of anti-angiogenesis strategies. Curr. Pharm. Des. 13 (20), 2025–2044. doi:10.2174/138161207781039689
Sanchez-Vega, F., Mina, M., Armenia, J., Chatila, W.K., Luna, A., La, K.C., et al. (2018). Oncogenic signaling pathways in the cancer genome atlas. Cell 173 (2), 321–337. doi:10.1016/j.cell.2018.03.035
Sargent, D., Sobrero, A., Grothey, A., O'Connell, M.J., Buyse, M., Andre, T., et al. (2009). Evidence for cure by adjuvant therapy in colon cancer: Observations based on individual patient data from 20, 898 patients on 18 randomized trials. J. Clin. Oncol. 27 (6), 872–877. doi:10.1200/JCO.2008.19.5362
Shen, H., Huang, C., Wu, J., Li, J., Hu, T., Wang, Z., et al. (2021). SCRIB promotes proliferation and metastasis by targeting hippo/YAP signalling in colorectal cancer. Front. Cell Dev. Biol. 9, 656359. doi:10.3389/fcell.2021.656359
Siegel, R.L., Miller, K.D., Goding Sauer, A., Fedewa, S.A., Butterly, L.F., Anderson, J.C., et al. (2020). Colorectal cancer statistics, 2020. Ca. Cancer J. Clin. 70 (3), 145–164. doi:10.3322/caac.21601
Sipos, F., Firneisz, G., and Muzes, G. (2016). Therapeutic aspects of c-MYC signaling in inflammatory and cancerous colonic diseases. World J. Gastroenterol. 22 (35), 7938–7950. doi:10.3748/wjg.v22.i35.7938
Suarez-Carmona, M., Lesage, J., Cataldo, D., and Gilles, C. (2017). EMT and inflammation: Inseparable actors of cancer progression. Mol. Oncol. 11 (7), 805–823. doi:10.1002/1878-0261.12095
Subramanian, A., Tamayo, P., Mootha, V.K., Mukherjee, S., Ebert, B.L., Gillette, M.A., et al. (2005). Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Terzi, J., Grivennikov, S., Karin, E., and Karin, M. (2010). Inflammation and colon cancer. Gastroenterology 138 (6), 2101–2114. e5. doi:10.1053/j.gastro.2010.01.058
Wu, D., Xue, W., Li, X., Yu, Z., Li, M., Jinxue, T., et al. (2019). Long noncoding RNA SNHG14 facilitates colorectal cancer metastasis through targeting EZH2-regulated EPHA7. Cell Death Dis. 10 (7), 514. doi:10.1038/s41419-019-1707-x
Xie, Y. H., Chen, Y. X., and Fang, J. Y. (2020). Comprehensive review of targeted therapy for colorectal cancer. Signal Transduct. Target. Ther. 5 (1), 22. doi:10.1038/s41392-020-0116-z
Xu, H., Jiang, Y., Xu, X., Su, X., Liu, Y., Ma, Y., et al. (2019). Inducible degradation of lncRNA Sros1 promotes IFN-gamma-mediated activation of innate immune responses by stabilizing Stat1 mRNA. Nat. Immunol. 20 (12), 1621–1630. doi:10.1038/s41590-019-0542-7
Xu, J., Meng, Q., Li, X., Yang, H., Chen, R., Gao, N., et al. (2019). Long noncoding RNA MIR17HG promotes colorectal cancer progression via miR-17-5p. Cancer Res. 79 (19), 4882–4895. doi:10.1158/0008-5472.CAN-18-3880
Yang, C., Huang, X., Liu, Z., Qin, W., and Wang, C. (2020). Metabolism-associated molecular classification of hepatocellular carcinoma. Mol. Oncol. 14 (4), 896–913. doi:10.1002/1878-0261.12639
Yang, J., Chen, Y., Luo, H., and Cai, H. (2020). The landscape of somatic copy number alterations in head and neck squamous cell carcinoma. Front. Oncol. 10, 321. doi:10.3389/fonc.2020.00321
Yang, K., Popova, N. V., Yang, W.C., Lozonschi, I., Tadesse, S., Kent, S., et al. (2008). Interaction of Muc2 and apc on wnt signaling and in intestinal tumorigenesis: Potential role of chronic inflammation. Cancer Res. 68 (18), 7313–7322. doi:10.1158/0008-5472.CAN-08-0598
Yang, Y., Feng, M., Bai, L. L., Zhang, M., Li, Q., Liao, W., et al. (2021). The effects of autophagy-related genes and lncRNAs in therapy and prognosis of colorectal cancer. Front. Oncol. 11, 582040. doi:10.3389/fonc.2021.582040
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics-a J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yuji, T., Shin-Ichi, A., Kazutaka, Y., Suminori, A., and Takashi, A. (1996). Angiogenesis as an unfavorable prognostic factor in human colorectal carcinoma. Cancer 78 (2), 226–231. doi:10.1002/(SICI)1097-0142(19960715)78:2<226::AID-CNCR6>3.0.CO;2-J
Zen, K., Yasui, K., Gen, Y., Dohi, O., Yoshikawa, T., MitSufuji, S., et al. (2009). Defective expression of polarity protein PAR-3 gene (PARD3) in esophageal squamous cell carcinoma. Oncogene 28 (32), 2910–2918. doi:10.1038/onc.2009.148
Zuo, M.S., Yang, J. Y., Wang, P. Y., Yang, H. O., and Ba, Y. (2020). Linking lipid metabolism with cell transformation and tumor progression. J. Nutr. Oncol. 5 (4), 153–161. doi:10.34175/jno202004001
Keywords: colorectal cancer, inflammation, lncRNAs, prognosis, immunotherapy
Citation: Huang M, Ye Y, Chen Y, Zhu J, Xu L, Cheng W, Lu X and Yan F (2022) Identification and validation of an inflammation-related lncRNAs signature for improving outcomes of patients in colorectal cancer. Front. Genet. 13:955240. doi: 10.3389/fgene.2022.955240
Received: 28 May 2022; Accepted: 29 August 2022;
Published: 30 September 2022.
Edited by:
Feng Wang, Affiliated Hospital of Nantong University, ChinaReviewed by:
Qiuxia Cui, Zhongnan Hospital, Wuhan University, ChinaMingqing Zhang, Tianjin Union Medical Center, China
Copyright © 2022 Huang, Ye, Chen, Zhu, Xu, Cheng, Lu and Yan. 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: Xiaofan Lu, xlu.cpu@foxmail.com; Fangrong Yan, f.r.yan@163.com