Wei Li
Guiyun Yang2
Hao Dong
Tongjun Liu- 1Department of Colorectal and Anal Surgery, The Second Hospital of Jilin University, Changchun, China
- 2Department of Operating Room, The Second Hospital of Jilin University, Changchun, China
- 3Department of Gastrointestinal Nutrition and Hernia Surgery, The Second Hospital of Jilin University, Changchun, China
- 4Department of Radiology, China-Japan Union Hospital of Jilin University, Changchun, China
Cuprotosis, a newly proposed mechanism of cell death, can trigger acute oxidative stress that leads to cell death by mediating protein lipidation in the tricarboxylic acid cycle. However, cuprotosis-related long non-coding RNAs (CRLNCs) and their relationship with prognosis and the immunological landscape of colorectal cancer (CRC) are unclear. We have developed a lncRNA signature to predict survival time, immune infiltration, and sensitivity to chemotherapy. CRLNCs were screened using the Cor function of the R software and the differentially expressed lncRNAs were collected with the limma package. Differentially expressed long non-coding RNAs (lncRNAs) associated with prognosis were selected using univariate regression analysis. A prognostic signature was developed using the least absolute shrinkage and selection operator (LASSO) and multivariate regression analysis. Patients with CRC were divided into two groups based on the risk score. The low-risk group had a more favorable prognosis, higher expression of immune checkpoints, and a higher level of immune cell infiltration compared with the high-risk group. Furthermore, there was a close association between the risk score and the clinical stage, tumor mutational burden, cancer stem cell index, and microsatellite instability. We also assessed chemotherapy response in the two risk groups. Our study analyzed the role of CRLNCs in CRC and provided novel targets and strategies for CRC chemotherapy and immunotherapy.
Introduction
Colorectal cancer (CRC) is responsible for approximately 10% of cancer cases and related deaths worldwide (1). Only in developed countries does the incidence of CRC show a stable or declining trend, which is primarily due to the widespread use of large-scale screening and colonoscopy, as well as the continuous improvement of people’s living and eating habits. It is estimated that there will be 25 million new cases of CRC worldwide by 2035 (2). In addition to surgery, radiotherapy and chemotherapy are still widely applied to reduce recurrence and improve survival. Chemotherapy, which involves the application of chemical compounds to inhibit the growth of tumor cells, is an indispensable part of the treatment process. Currently, platinum-based chemotherapy in combination with 5-fluorouracil is the first-line treatment option in treating patients with CRC (3, 4). However, different patients respond differently to the same chemotherapy regimen, leading to large variations in patient prognoses (5, 6). Exploring new, specific, and effective targets related to chemotherapy sensitivity, as well as recognizing individualized and precise treatment, is therefore critical for CRC therapy.
Copper homeostasis is an ancient phenomenon in living organisms. Copper is an indispensable trace element for the homeostasis of the internal environment (7). Copper contributes to the progression of tumors, such as breast and lung cancer, where it is involved in tumor angiogenesis, epithelial-mesenchymal transition, and cell proliferation and metastasis (8, 9). Therefore, copper-chelating agents have been studied and reported to inhibit tumor growth in some clinical trials (10). Meanwhile, copper can also promote oxidative stress to mediate cell death (11, 12). Copper-specific ionophores can transport copper into cells at specific sites, increasing the copper level in tumor cells, and then mediating the toxicity of copper overload, which results in cell death (10). The role of copper in the treatment of tumors is complex and versatile. Mutations in lncRNAs are believed to mediate several forms of tumor development along with protein-coding genes (13).
LncRNAs can regulate immune and inflammatory responses at the transcriptional and posttranscriptional levels by interacting with proteins, RNA, and DNA (14). At the same time, lncRNAs have a close relationship with the tumor microenvironment (TME) (15). Several lncRNAs, including TUG1, MALAT1, H19, GAS5, LINC00152, UCA1, CUDR, and AA174084, have been identified as predictive biomarkers of CRC. Investigating such lncRNAs as potential targets for CRC therapy is of long-term value. GAS5 is involved in regulating chemotherapy resistance in CRC. The other lncRNAs require further investigation. To determine whether cuprotosis-related long non-coding RNAs (CRLNCs) play a role in CRC, a prognostic signature of the immune infiltration and survival of patients with CRC was developed. A different prognosis was revealed by the Kaplan-Meier analysis. Various methodologies, such as XCELL, TIMER, and ssGSEA, were also used to analyze the immune infiltration in patients with CRC. The analyses of immune checkpoints, clinicopathological data, tumor mutational burden (TMB), cancer stem cells (CSCs), microsatellite instability (MSI), and chemotherapy response were also performed.
Materials and methods
Datasets and samples
The transcriptome, mutation, and clinical data for COAD containing 32 healthy tissues and 375 tumors were downloaded from The Cancer Genome Atlas (TCGA) database. Six fresh frozen CRC and paracancerous paired tissues were obtained from the Second Hospital of Jilin University. The cuprotosis-related genes (CRGs) are shown in Supplementary Table 1.
Identification of differentially expressed cuprotosis-related long non-coding RNAs
To identify lncRNAs closely related to CRGs, we performed a screen using the Cor function of the R software, with the filter conditions set to require a correlation coefficient of >0.3 with a false discovery rate of <0.001. Subsequently, differentially expressed CRLNCs between the 32 normal and 375 tumor samples were selected using the limma package (| log Foldchange| > 1 and false discovery rate < 0.05).
Construction and validation of a prognostic long non-coding RNA signature
The Supplementary material provide details about the construction and validation of the prognostic lncRNA signature.
Gene set enrichment analysis and nomogram construction
The gene set enrichment analysis (GSEA) and nomogram are presented in the Supplementary material.
Immune landscape, immune checkpoints, and clinical data analyses
Analyses of the immune landscape, immune checkpoints, and clinical data are presented in the Supplementary material.
Analyses of tumor mutational burden, cancer stem cells, and microsatellite instability
Tumor mutational burden is an essential marker of immunotherapy response and prognosis. Therefore, we compared genetic mutations in samples from low-risk and high-risk groups. The mutational burdens from all samples were then calculated and compared. A correlation analysis was applied to determine the significant relationships between the risk scores, TMB, and immune infiltration. We also explored the link between CRGs and risk scores. MSI could reflect the effect of immunotherapy. Therefore, the association between MSI and risk score was analyzed. We compared patients’ survival times between MSI-H and MSS/MSI-L. We also integrated MSI into the signature for survival analysis.
Drug sensitivity analysis and identification of differential genes
The limma package was used to identify the differentially expressed genes (DEGs) (| log Foldchange| > 1 and false discovery rate < 0.05). To further search for the hub genes, we used the CytoNCA plugin in Cytoscape software. Based on the scores of Betweenness, Closeness, Degree, Eigenvector, LAC, and Network, we screened the DEGs twice to obtain core genes. Furthermore, Gene Ontology and Kyoto Encyclopedia Genes and Genomes pathway enrichment analyses were used to explore the functional pathways based on the DEGs. Finally, to investigate the differences in response to chemotherapy, we used the pRRophetic package to predict drug sensitivity.
Quantitative real-time PCR
Total RNA was extracted from CRC tissues using Trizol reagent (Invitrogen, Carlsbad, CA, United States). We used a reverse transcription kit (Takara, Tokyo, Japan) to synthesize cDNA. The SYBR Premix Ex Taq™ kit (Takara, Japan) was used to perform the quantitative real-time PCR (RT-qPCR). The expression level of LINC00412, AC016737.1, AC026782.2, AC090204.1, AC129507.1, and AC116914.2 was normalized using glyceraldehyde-3-phosphate dehydrogenase. The data were analyzed using the 2–ΔΔCt method. The primers of the seven genes are listed in Supplementary Table 2.
Statistical analyses
All statistical analyses were performed using R version 4.1.1. P < 0.05 was considered significant.
Results
Analysis of differentially expressed cuprotosis-related long non-coding RNAs
The study design is presented in Supplementary Figure 1. The Cor function was performed to select 880 CRLNCs. The association between CRGs and lncRNAs is shown in Figure 1A. We discovered 487 CRLNCs that were differentially expressed, with 445 being upregulated in CRC and 42 being downregulated (Figures 1B–C).
Figure 1. (A) Network of long noncoding RNAs and cuprotosis-related genes. (B) Differential expression of cuprotosis-related long non-coding RNAs (CRLNCs) between normal and tumor tissues. (C) The volcano plot of CRLNCs. Red represents upregulated CRLNCs; green represents downregulated CRLNCs. CRGs, cuprotosis-related genes; CRLNCs, cuprotosis-related long non-coding RNAs.
Construction and validation of the prognostic signature
As a result of the univariate regression analysis, six CRLNCs were discovered to be linked with the prognosis of patients with CRC (Figure 2A). We then selected the genes corresponding to the smallest lambda value for the multivariate Cox regression analysis (Figures 2B–C). Finally, LINC00412, AC016737.1, AC026782.2, AC090204.1, AC129507.1, and AC116914.2 were screened to construct the risk signature. The formula of the risk signature is as follows:
Figure 2. (A) Univariate Cox regression analysis. (B) The LASSO algorithm further selected the most crucial genes and LASSO coefficient profiles. (C) Expression of risk lncRNAs between normal and tumor tissues. (D) Sankey diagram showing the relationship between CRGs and lncRNAs. (E) Heatmap showing the expression of risk lncRNAs between low- and high-risk groups, and the ranked dot plot showing the risk score distribution in all samples. (*P < 0.05; **P < 0.01; ***P < 0.001; ns, not significant). LASSO, least absolute shrinkage and selection operator; CRGs, cuprotosis-related genes.
High-risk lncRNAs included AC016737.1, AC026782.2, AC090204.1, and AC129507.1. Low-risk lncRNAs included LINC00412 and AC116914.2 (Figure 2E). We found that LINC00412, AC016737.1, AC026782.2, and AC090204.1 were highly expressed in CRC. AC129507.1 was downregulated in CRC (Figure 2C). The relationship between CRGs and the lncRNAs is displayed in Figure 2D. The low-risk group showed better survival outcomes than the high-risk group (Figures 2E, 3A). The training group and test group confirmed this conclusion (Figures 3B–E). The area under the curve (AUC) values demonstrated that our prognostic signature had moderate performance (Figure 3F). The AUC values of age, gender, grade, and tumor stage were 0.575, 0.524, 0.558, and 0.586, respectively, indicating that the risk model had the best predictive ability (Figure 3G). The Kaplan-Meier survival curve further proved that this risk signature applied to patients of any age, gender, grade, and TNM stage (Figures 4A–G).
Figure 3. (A) Survival curve of low- and high-risk groups in all samples. The survival analysis of low- and high-risk groups in the training group (B) and the test group (C). Heatmap displaying the expression of risk lncRNAs between low- and high-risk groups, and the ranked dot plot showing the risk score distribution in the training group (D) and the test group (E). (F) The AUC values of 1-, 3-, and 5-year survival. (G) Comparison of AUC between risk signature and age, gender, grade, and stage. AUC, area under the curve.
Figure 4. Kaplan–Meier analyses of age (A), gender (B), grade (C), T stage (D), N stage (E), M stage (F), and stage (G). (H) GSEA of low- and high-risk groups. GSEA, gene set enrichment analysis.
Gene set enrichment analysis and nomogram construction
The results of the GSEA showed that the main functional pathways in the high-risk group were the calcium signaling pathway, GAP junction, extracellular matrix receptor interaction, and complement and coagulation cascades. The main functional pathways in the low-risk group were homologous recombination, neuroactive ligand-receptor interaction, retinoic acid-inducible gene-I-like receptor signaling pathway, and RNA degradation (Figure 4H). The univariate regression analysis showed that age (hazard ratio [HR]: 1.003–1.039; P < 0.05), stage (HR: 1.193–1.833; P < 0.001), and risk score (HR: 1.054–1.125; P < 0.001) were associated with prognosis (Figure 5A). It was demonstrated that age (HR: 1.011–1.049; P < 0.001), stage (HR: 1.233–2.075; P < 0.001), and risk score (HR: 1.059–1.129; P < 0.001) were found to be independent prognostic factors (Figure 5B). A prognostic nomogram was also developed for the prediction of the survival time (Figure 5C). The calibration curves for the 1-, 3-, and 5-year survival rates confirmed the accuracy of the nomogram (Figure 5D). The decision curve analysis indicated that the nomogram had a better predictive ability for survival time than the stage, age, and risk score (Figure 5E). The AUC values of the stage, age, risk score, and nomogram were 0.590, 0.571, 0.679, and 0.716, respectively (Figure 5F).
Figure 5. (A) Univariate Cox regression analysis. (B) Multivariate regression analysis. (C) Nomogram to predict the 1-, 3-, and 5-year survival probability for patients with CRC. (D) The calibration curves of 1-, 3-, and 5-year survival. (E) Decision curve analysis comparing stage, age, risk score, and the nomogram. (F) AUC values of the stage, age, risk score, and nomogram. CRC, colorectal cancer; AUC, area under the curve.
Analyses of immune infiltration, immune checkpoints, and clinical data
The low-risk group had more immune cell infiltration (Figure 6B). Additionally, we discovered that some immune cell types had positive correlations with risk scores while others had negative correlations (Figure 6A). The results of ssGSEA indicated that some activities, such as inflammation promotion, were upregulated in the high-risk group (Figures 6C–D). There was a differential expression of 17 immune checkpoints, of which 16 (94.12%) had higher expression levels in the low-risk group (Figure 7A). Meanwhile, higher risk scores were observed in late-stage CRC (Figure 7B).
Figure 6. (A) Correlation between risk score and immune cells. (B) XCELL, TIMER, QUANTISEQ, MCPCOUNTER, EPIC, and CIBERSORT to analyze the immune landscape of patients with CRC. (C,D) ssGSEA analysis. CRC, colorectal cancer; GSEA, gene set enrichment analysis.
Figure 7. (A) Expression of immune checkpoints in the two risk groups. (B) Clinical features of the model. Somatic mutations in the low-risk group (C) and the high-risk group (D). (E) Association between risk score and CRGs. (F) TMB of the two risk groups. (G) The relationship between TMB and risk score. (I,J) Correlation between CSC and risk score. (H) Circle picture displaying the relationship between TMB and immune cells. CRGs, cuprotosis-related genes; TMB, tumor mutational burden; CSC, cancer stem cell;
Characteristics of tumor mutational burden, cancer stem cells, and microsatellite instability
Both low- and high-risk groups had the same top five mutated genes. However, the low-risk group demonstrated a higher mutation probability (Figures 7C–D). TMB was high in the low-risk group (Figure 7F). TMB and risk score were inversely correlated (Figure 7G). Differences in TMB between the two groups may be related to endothelial cells and neutrophils (Figure 7H). Patients with high-risk scores had lower RNA and higher DNA markers in CSCs than patients with low-risk scores (Figures 7I–J). Figure 7E illustrates the relationship between risk scores and CRGs. MSS/MSI-L was strongly associated with higher risk scores (Figures 8A–B). There was no association between survival rates and MSS/MSI-L and MSI-H status. However, the MSI-H + low-risk score had the most favorable prognosis compared with the MSS/MSI-L + high-risk score, the MSS/MSI-L + low-risk score, and the MSI-H + high-risk score groups (Figure 8D).
Figure 8. (A,B) Relationship between MSI and risk score. (C) Correlation between survival time and MSS/MSI-L and MSI-H. (D) Survival analysis integrating MSI status and risk signature. (E) Differentially expressed genes (DEGs) between the two risk groups. (F) Volcano plot of DEGs. (G–I) Core DEGs in the two risk groups. (J) Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses of DEGs. MSI, microsatellite instability.
Analyses of chemotherapeutic drug sensitivity and differential genes
After running the limma package, we selected 140 DEGs (Figure 8E), including 131 overexpressed and 9 underexpressed genes (Figure 8F). Based on the scores of Betweenness, Closeness, Degree, Eigenvector, LAC, and Network, we performed two screenings and obtained eight core genes (Figures 8G–I). The Gene Ontology and Kyoto Encyclopedia Genes and Genomes analyses based on the 140 DEGs showed that the cGMP-PKG signaling pathway, heparin binding, and contractile fiber may be the main biological functions. Finally, we screened a total of 37 chemotherapeutic drugs to evaluate the differences in sensitivity between the two groups. Out of these 37 chemotherapeutic drugs, 32 (86.49%) had higher IC50 values in the low-risk group than in the high-risk group. This suggests that high-risk patients may be more sensitive to chemotherapy drugs, contributing to a more favorable prognosis (Figure 9).
Figure 9. Drug sensitivity analysis between high- and low-risk groups.
Discussion
Numerous studies have demonstrated that cell death is related to tumor occurrence and progression (16). The known mechanisms of cell death mainly include ferroptosis, pyroptosis, necroptosis, apoptosis, and autophagy (17–21). Tsvetkov et al. first described cuprotosis, a new mechanism of cell death (12). However, the relationship between cuprotosis and CRC is unclear, especially CRLNCs. In this study, we identified CRLNCs using bioinformatics studies. Six lncRNAs (i.e., LINC00412, AC016737.1, AC026782.2, AC090204.1, AC129507.1, and AC116914.2) were then used to construct a prognostic model for survival time prediction, immune infiltration, and chemotherapy drug sensitivity of CRC. LINC00412 was included in the 10 biomarkers for the construction of the cardia cancer prognostic model and contributed to the modification of the prognostic models by Xin et al. (22). Taniguchi-Ponciano et al. found that LINC00412 was upregulated in all kinds of pituitary tumors (23). Zhang et al. provided evidence that AC016737.1 was associated with the inflammatory response of CRC, and this association was later validated by constructing a prognostic model (24). Chen et al. used hypoxia-related lncRNAs, which included AC016737.1, to construct a model for predicting the outcome of CRC. AC016737.1 was also used in the m6A-modified lncRNA prognostic nomogram by Song et al. and the immune-related lncRNA pair model by Shenglei et al. for CRC prognosis. Zha et al. discovered that AC129507.1 played a role in predicting prognosis and multiple tumor-related pathways in CRC (25). AC129507.1 was also used in the exosome-related lncRNA CRC prognostic model by Li et al. as well as in the survival prediction model of gastroesophageal junction adenocarcinoma by Song et al. AC116914.2 was involved in autophagy, m6A RNA methylation, and hypoxia in head and neck squamous cell carcinoma. Meanwhile, AC116914.2 was associated with survival and immune cell infiltration (26–28). The results of q-RT PCR showed that AC016737.1, AC026782.2, AC090204.1, and AC129507.1 were highly expressed in tumor cells compared with normal cell, while LINC00412 and AC116914.2 were low expressed using Tukey’s HSD posttest as the method of multiple comparisons. This is consistent with the fact that AC016737.1, AC026782.2, AC090204.1, and AC129507.1 were high-risk genes, and LINC00412 and AC116914.2 were low-risk genes.
The tumor microenvironment is well known to be the site of tumor survival, with multiple components interacting to form a complex and polymorphic environment (29). Immune cell infiltration is one of the key components of TME. Comprehensive analysis of immunological signatures in the TME could facilitate the progress of native and effective immunotherapeutic strategies, as well as the discovery of highly effective biomarkers (30). TME also plays an important role in regulating tumor sensitivity to treatment (31). B cell is the most crucial humoral immune cell, mediating the antitumor response. It is associated with a favorable prognosis and immunotherapy response (32). Notably, the low-risk group had a higher infiltration level of B cell, B cell memory, B cell plasma, and naïve B cell. CD4 + T cells can kill tumors either directly by destroying the tumor cells or indirectly by mediating TME regulation. In addition, CD4 + T cells also can promote gene expression and differentiation of CD8 + T cells (33–35). As a result, we found that some T cell types, including T cell CD4 + memory and T cell gamma delta, were present in higher levels in the low-risk group than in the high-risk school. The first line of defense in identifying tumors is the ability of CD8 + T cells to recognize MHC class I molecules expressed by tumor cells. CD8 + T cells are the most efficient immune cells against cancer (36). In this study, the infiltration level of CD8 + T cells and CD8 + central memory T cells were higher in the low-risk group than in the high-risk group. Cancer-associated fibroblasts (CAFs) in the TME have been shown to promote the proliferation of multiple tumors by secreting a variety of biological factors to suppress the immune response (37). Various molecules, such as epidermal growth factor and interleukin-6, can be secreted by CAFs to enhance cell proliferation, tumor invasion and metastasis, and epithelial-mesenchymal transition. Notably, a higher infiltration level of CAFs was observed in the high-risk group than in the low-risk group, possibly resulting in the difference in prognosis between the two groups. Meanwhile, we also discovered that the patients with low-risk scores obtained higher immune scores and lower stromal scores than those with high-risk scores, which further gives a reasonable explanation for the difference in prognosis between the two groups.
Immune checkpoint inhibitors have become a promising treatment strategy in almost all kinds of malignant tumors. Several clinical trials involving nivolumab, pembrolizumab, ipilimumab, avelumab, and durvalumab have either been completed or are currently being conducted. In the low-risk group, we found the overexpression of 16 immune checkpoints, which could reveal potential immune therapy targets and help develop combination therapies and predictive biomarkers.
Cancer genomics studies have found that most cancers develop with the accumulation of somatic gene mutations (38). In this study, a higher mutation probability was detected in the low-risk group than in the high-risk group. It is widely known that TMB and MSI are predictive biomarkers of immunotherapy response. High TMB and MSI-H appear to be associated with favorable immunotherapy response and prognosis (39, 40). Our findings also confirmed this conclusion and may contribute to revealing potential therapeutic targets.
We also found that patients in the low-risk group were more sensitive to chemotherapy than patients in the high-risk group. Fluorouracil-based adjuvant chemotherapy is recommended for resected stage III and some stage II colon cancers to improve patient survival. Several studies have concentrated on the addition of oxaliplatin to fluorouracil as a novel standardized CRC chemotherapy (41–44). The standard course of adjuvant chemotherapy is 6 months. A major disadvantage of oxaliplatin chemotherapy is cumulative sensory neuropathy. In a clinical trial, 3-month adjuvant chemotherapy in low-risk stage III (not T4 or N2) colon cancer did not compromise treatment efficacy but reduced drug toxicity (such as neuropathy) (45). Chemotherapy sensitivity is vital for CRC treatment (46). Our prognostic signature can help make chemotherapy more effective or tailor treatment to each individual, which is critical for survival.
Our study also had several limitations. Comprehensive and detailed in vitro and in vivo experiments are still needed to further validate our conclusion. Also, more clinical samples need to be included.
Conclusion
Based on CRLNCs, a prognostic signature was constructed to predict the survival and chemotherapy sensitivity of patients with CRC. In summary, our study analyzed the role of CRLNCs in CRC and provided new targets and strategies for CRC therapy.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://portal.gdc.cancer.gov/repository.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Second Hospital of Jilin University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
TL: work concept or design. JZ: data collection. WL and HD: draft the manuscript. GY: make important revisions to the manuscript and approved the final manuscript. All authors read and approved the final manuscript.
Funding
This research was supported by the Jilin Provincial Department of Finance (No. 2020SCZT004).
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/fmed.2022.1055785/full#supplementary-material
References
1. Bray F, Ferlay J, Soerjomataram I, Siegel R, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424.
2. Ait Ouakrim D, Pizot C, Boniol M, Malvezzi M, Boniol M, Negri E, et al. Trends in colorectal cancer mortality in Europe: retrospective analysis of the WHO mortality database. BMJ. (2015) 351:h4970.
3. Vodenkova S, Buchler T, Cervena K, Veskrnova V, Vodicka P, Vymetalkova VJP, et al. 5-fluorouracil and other fluoropyrimidines in colorectal cancer: past, present and future. Pharmacol Ther. (2020) 206:107447. doi: 10.1016/j.pharmthera.2019.107447
6. Dekker E, Tanis P, Vleugels J, Kasi P, Wallace MJL. Colorectal cancer. Lancet. (2019) 394:1467–80.
7. Bhattacharjee A, Chakraborty K, Shukla A. Cellular copper homeostasis: current concepts on its interplay with glutathione homeostasis and its implication in physiology and human diseases. Metallomics. (2017) 9:1376–88. doi: 10.1039/c7mt00066a
8. da Silva DA, De Luca A, Squitti R, Rongioletti M, Rossi L, Machado CML, et al. Copper in tumors and the use of copper-based compounds in cancer treatment. J Inorg Biochem. (2022) 226:111634.
10. Denoyer D, Masaldan S, La Fontaine S, Cater MA. Targeting copper in cancer therapy: “copper that cancer”. Metallomics. (2015) 7: 1459–76.
11. Wang Y, Chen R, Gan Y, Ying S. The roles of PAD2- and PAD4-mediated protein citrullination catalysis in cancers. Int J Cancer. (2021) 148:267–76. doi: 10.1002/ijc.33205
12. Suthiphasilp V, Raksat A, Maneerat T, Hadsadee S, Jungsuttiwong S, Pyne SG, et al. Cytotoxicity and nitric oxide production inhibitory activities of compounds isolated from the plant pathogenic fungus Curvularia sp. J Fungi (Basel). (2021) 7:408. doi: 10.3390/jof7060408
13. Bhan A, Soleimani M, Mandal SS. Long noncoding RNA and cancer: a new paradigm. Cancer Res. (2017) 77:3965–81.
14. Robinson EK, Covarrubias S, Carpenter S. The how and why of lncRNA function: an innate immune perspective. Biochim Biophys Acta Gene Regul Mech. (2020) 1863:194419. doi: 10.1016/j.bbagrm.2019.194419
15. Liao CL, Hu N, Sun XY, Zhou Q, Tian M, Cao Y, et al. Identification and validation of tumor microenvironment-related lncRNA prognostic signature for uveal melanoma. Int J Ophthalmol. (2021) 14:1151–9. doi: 10.18240/ijo.2021.08.03
16. Strasser A, Vaux DL. Cell death in the origin and treatment of cancer. Mol Cell. (2020) 78:1045–54.
17. Mou Y, Wang J, Wu J, He D, Zhang C, Duan C, et al. Ferroptosis, a new form of cell death: opportunities and challenges in cancer. J Hematol Oncol. (2019) 12:34.
18. Tang R, Xu J, Zhang B, Liu J, Liang C, Hua J, et al. Ferroptosis, necroptosis, and pyroptosis in anticancer immunity. J Hematol Oncol. (2020) 13:110.
19. Hsu SK, Li CY, Lin IL, Syue WJ, Chen YF, Cheng KC, et al. Inflammation-related pyroptosis, a novel programmed cell death pathway, and its crosstalk with immune therapy in cancer treatment. Theranostics. (2021) 11:8813–35.
20. Kim C, Kim B. Anti-cancer natural products and their bioactive compounds inducing ER stress-mediated apoptosis: a review. Nutrients. (2018) 10:1021. doi: 10.3390/nu10081021
22. Tsvetkov PCS, Petrova B, Dreishpoon M, Verma A, Abdusamad M, Rossen J, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Science. (2022) 375:1254–61.
23. Taniguchi-Ponciano K, Portocarrero-Ortiz L, Guinto G, Moreno-Jimenez S, Gomez-Apo E, Chavez-Macias L, et al. The kinome, cyclins and cyclin-dependent kinases of pituitary adenomas, a look into the gene expression profile among tumors from different lineages. BMC Med Genomics. (2022) 15:52. doi: 10.1186/s12920-022-01206-y
24. Zhang S, Li X, Tang C, Kuang W. Inflammation-related long non-coding RNA signature predicts the prognosis of gastric carcinoma. Front Genet. (2021) 12:736766. doi: 10.3389/fgene.2021.736766
25. Zha Z, Zhang P, Li D, Liu G, Lu L. Identification and construction of a long noncoding RNA prognostic risk model for stomach adenocarcinoma patients. Dis Markers. (2021) 2021:8895723.
26. Shen L, Li N, Zhou Q, Li Z, Shen L. Development and validation of an autophagy-related LncRNA prognostic signature in head and neck squamous cell carcinoma. Front Oncol. (2021) 11:743611. doi: 10.3389/fonc.2021.743611
27. Feng ZY, Gao HY, Feng TD. Immune infiltrates of m6A RNA methylation-related lncRNAs and identification of PD-L1 in patients with primary head and neck squamous cell carcinoma. Front Cell Dev Biol. (2021) 9:672248. doi: 10.3389/fcell.2021.672248
28. Yang C, Zheng X. Identification of a hypoxia-related lncRNA biomarker signature for head and neck squamous cell carcinoma. J Oncol. (2022) 2022:6775496. doi: 10.1155/2022/6775496
29. Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res. (2019) 79:4557–66.
30. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. (2013) 14:1014–22.
31. Zhong C, Xie T, Chen L, Zhong X, Li X, Cai X, et al. Immune depletion of the methylated phenotype of colon cancer is closely related to resistance to immune checkpoint inhibitors. Front Immunol. (2022) 13:983636. doi: 10.3389/fimmu.2022.983636
32. Engelhard V, Conejo-Garcia JR, Ahmed R, Nelson BH, Willard-Gallo K, Bruno TC, et al. B cells and cancer. Cancer Cell. (2021) 39:1293–6.
33. Borst J, Ahrends T, Babala N, Melief CJM, Kastenmuller W. CD4(+) T cell help in cancer immunology and immunotherapy. Nat Rev Immunol. (2018) 18:635–47.
34. Oh DY, Fong L. Cytotoxic CD4(+) T cells in cancer: expanding the immune effector toolbox. Immunity. (2021) 54:2701–11. doi: 10.1016/j.immuni.2021.11.015
35. Zander R, Schauder D, Xin G, Nguyen C, Wu X, Zajac A, et al. CD4(+) T cell help is required for the formation of a cytolytic CD8(+) T cell subset that protects against chronic infection and cancer. Immunity. (2019) 51:1028–42.e4. doi: 10.1016/j.immuni.2019.10.009
36. Raskov H, Orhan A, Christensen JP, Gogenur I. Cytotoxic CD8(+) T cells in cancer and cancer immunotherapy. Br J Cancer. (2021) 124:359–67.
37. Li QQ, Hao JJ, Zhang Z, Krane LS, Hammerich KH, Sanford T, et al. Proteomic analysis of proteome and histone post-translational modifications in heat shock protein 90 inhibition-mediated bladder cancer therapeutics. Sci Rep. (2017) 7:201. doi: 10.1038/s41598-017-00143-6
38. Iranzo J, Martincorena I, Koonin EV. Cancer-mutation network and the number and specificity of driver mutations. Proc Natl Acad Sci USA. (2018) 115:E6010–9.
39. Luchini C, Bibeau F, Ligtenberg MJL, Singh N, Nottegar A, Bosse T, et al. ESMO recommendations on microsatellite instability testing for immunotherapy in cancer, and its relationship with PD-1/PD-L1 expression and tumour mutational burden: a systematic review-based approach. Ann Oncol. (2019) 30:1232–43. doi: 10.1093/annonc/mdz116
40. Liu L, Bai X, Wang J, Tang XR, Wu DH, Du SS, et al. Combination of TMB and CNA stratifies prognostic and predictive responses to immunotherapy across metastatic cancer. Clin Cancer Res. (2019) 25:7413–23. doi: 10.1158/1078-0432.CCR-19-0558
41. André T, Boni C, Navarro M, Tabernero J, Hickish T, Topham C, et al. Improved overall survival with oxaliplatin, fluorouracil, and leucovorin as adjuvant treatment in stage II or III colon cancer in the MOSAIC trial. J Clin Oncol. (2009) 27:3109–16.
42. Schmoll H, Tabernero J, Maroun J, de Braud F, Price T, Van Cutsem E, et al. Capecitabine Plus Oxaliplatin compared with fluorouracil/folinic acid as adjuvant therapy for stage III colon cancer: final results of the NO16968 randomized controlled phase III trial. J Clin Oncol. (2015) 33:3733–40.
43. Land S, Kopec J, Cecchini R, Ganz P, Wieand H, Colangelo L, et al. Neurotoxicity from oxaliplatin combined with weekly bolus fluorouracil and leucovorin as surgical adjuvant chemotherapy for stage II and III colon cancer: NSABP C-07. J Clin Oncol. (2007) 25:2205–11. doi: 10.1200/JCO.2006.08.6652
44. Yothers G, O’Connell M, Allegra C, Kuebler J, Colangelo L, Petrelli N, et al. Oxaliplatin as adjuvant therapy for colon cancer: updated results of NSABP C-07 trial, including survival and subset analyses. J Clin Oncol. (2011) 29:3768–74. doi: 10.1200/JCO.2011.36.4539
45. Grothey A, Sobrero AF, Shields AF, Yoshino T, Paul J, Taieb J, et al. Duration of adjuvant chemotherapy for stage III colon cancer. N Engl J Med. (2018) 378:1177–88.
Keywords: colorectal cancer, cuprotosis, prognostic signature, immune infiltration, immunotherapy, chemotherapy sensitivity
Citation: Li W, Yang G, Dong H, Zhu J and Liu T (2022) A prognostic signature based on cuprotosis-related long non-coding RNAs predicts the prognosis and sensitivity to chemotherapy in patients with colorectal cancer. Front. Med. 9:1055785. doi: 10.3389/fmed.2022.1055785
Received: 28 September 2022; Accepted: 17 October 2022;
Published: 16 November 2022.
Edited by:
Xingzhe Ma, Fresh Wind Biotechnologies, United StatesReviewed by:
Wangsheng Xue, Tongji University, ChinaGuiyu Wang, The Second Affiliated Hospital of Harbin Medical University, China
Copyright © 2022 Li, Yang, Dong, Zhu and Liu. 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: Jiajing Zhu, zhujj16@mails.jlu.edu.cn; Tongjun Liu, tongjunliu@163.com