Yun Chen- Department of TCM Gynecology, Hangzhou TCM Hospital Affiliated to Zhejiang Chinese Medical University, Hangzhou, China
Background: Uterine corpus endometrial carcinoma (UCEC) is a common gynecological malignancy globally with high recurrence and mortality rates. Cuproptosis is a new type of programmed cell death involved in tumor cell proliferation and growth, angiogenesis, and metastasis.
Methods: The difference in cuproptosis-related genes (CRGs) between UCEC tissues and normal tissues deposited in The Cancer Genome Atlas database was calculated using the “limma” R package. LASSO Cox regression analysis was conducted to construct a prognostic cuproptosis–related signature. Kaplan–Meier analysis was conducted to compare the survival of UCEC patients. A ceRNA network was constructed to identify the lncRNA–miRNA–mRNA regulatory axis. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to verify CRG expression in UCEC.
Results: The expression of FDX1, LIAS, DLAT, and CDKN2A were upregulated, whereas the expression of LIPT1, DLD, PDHB, MTF1, and GLS were downregulated in UCEC versus normal tissues. The genetic mutation landscape of CRGs in UCEC was also summarized. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that these CRGs were enriched in the tricarboxylic acid (TCA) cycle, glycolysis, and HIF-1 signaling pathway. LASSO Cox regression analysis was performed and identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1). UCEC patients with high risk scores had a poor prognosis with an area under the curve of 0.782 and 0.764 on 3- and 5-year receiver operating characteristic curves. Further analysis demonstrated a significant correlation between CDKN2A and pTNM stage, tumor grade, immune cell infiltration, drug sensitivity, tumor mutational burden (TMB) score, and microsatellite instable (MSI) score. The data validation of qRT-PCR further demonstrated the upregulation of CDKN2A and the downregulation of LIPT1 and GLS in UCEC versus normal tissues. The ceRNA network also identified lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis for UCEC.
Conclusion: The current study identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1) for UCEC. The ceRNA network also identified that lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis may be involved in the progression of UCEC. Further in vivo and in vitro studies should be conducted to verify these results.
Introduction
Uterine corpus endometrial carcinoma (UCEC), originating from the endometrium, is a common gynecological malignancy globally (Yang et al., 2020a). In 2020, there is an estimated 417,367 new cases of UCEC and an estimated 97,370 disease-related deaths worldwide (Sung et al., 2021). Although the cure rate of early UCEC can reach 95%, the prognosis of UCEC was still poor, because of its high recurrence rate and high mortality rate (Wang et al., 2020; Cui et al., 2021). Furthermore, there are few satisfactory biomarkers or approaches that could predict the prognosis of UCEC patients accurately. Despite a lot of effort having been made to explore the mechanism of the tumorigenesis and progression of UCEC, the specific mechanism was not fully elucidated. Thus, identifying the prognostic signature and associated potential mechanism of UCEC is greatly significant.
Cuproptosis is a new type of programmed cell death identified by Tsvetkov et al. (2022). Excess copper could result in mitochondrial protein aggregation and cause a distinct form of cell death (Kahlson and Dixon, 2022). Copper is fundamental material in various materials in biological processes including mitochondrial respiration, iron uptake, and antioxidant/detoxification processes (Ruiz et al., 2021). Increasing pieces of evidence demonstrated the involvement of copper in tumor cell proliferation and growth, angiogenesis, and metastasis (Oliveri, 2022). Different from apoptosis, ferroptosis, pyroptosis, and necroptosis, limited studies have been conducted on the role of cuproptosis-related genes (CRGs) in the progression and prognosis of cancers, including UCEC. Hence, the present study focus on the role of CRGs in UCEC.
Competing endogenous RNAs (ceRNAs) could regulate each other via the competition for shared miRNAs at the posttranscriptional level (Qi et al., 2015). The lncRNA–miRNA–mRNA ceRNA network has been suggested to play a vital role in many types of cancers, including LUAD (Wu et al., 2020). More and more studies suggested the importance of lncRNA–miRNA–mRNA ceRNA network in the study of the molecular mechanism of malignancies (Wang et al., 2019; Chen et al., 2021a).
The Cancer Genome Atlas (TCGA, https://cancergenome.nih.gov/) is a landmark cancer genomics program that contains a lot of clinical and bioinformatics data on cancers (Yang et al., 2013). Increasing pieces of evidence have highlighted database mining as one of the promising approaches to studying the molecular mechanism and prognostic signatures of cancer. Thus, bioinformatics analyses were applied in the current study to identify cuproptosis-related prognostic signature and associated regulatory axis in UCEC.
Materials and Methods
Dataset and Preprocessing
The flowchart of the current study is shown in Supplementary Figure S1. Based on a previous study, a total of 10 CRGs were obtained, including FDX1, LIAS, LIPT1, DLD, DLAT, PDHA1, PDHB, MTF1, GLS, and CDKN2A (Yang et al., 2020a). Using TCGA (https://portal.gdc.cancer.gov/) database (Tomczak et al., 2015), the gene expression profile of CRGs of UCEC was isolated and R (version 4.0.5) (Kang et al., 2021) was used to perform data preprocessing with associated packages. The RNA-seq data were then normalized to transcripts per kilobase million value.
Mutation Landscape, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes Analysis
After obtaining the single nucleotide variation (SNV) and copy number variation (CNV) data of UCEC from the TCGA database, we then explore their mutation landscape in UCEC using R software with a “maftools” package (Mayakonda et al., 2018). Seven types of mutation were included in this analysis: Missense_Mutation, Nonsense_Mutation, Frame_Shift_Ins, Splice_Site, Frame_Shift_Del, In_Frame_Del, and In_Frame_Ins. In Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, the “ggplot2” package (Maag, 2018) was utilized to perform this analysis, and p < 0.05 was set.
Cuproptosis-Related Prognostic Signature Analysis
Prognosis analysis, including overall survival (OS), progression-free survival (PFS), and disease-free survival (DFS), was conducted to identify potential prognostic biomarkers. Log-rank test was performed to calculate the p-values, hazard ratio (HR), and 95% confidence interval. Using these potential prognostic biomarkers, LASSO Cox regression analysis was performed to construct a cuproptosis-related prognostic signature. After calculating the risk score of each patient with UCEC [Risk score =
Single Gene Analysis
After obtaining the abundance of immune cells using TIMER (https://cistrome.shinyapps.io/timer/), the correlation between cuproptosis-related prognostic genes and immune cells was analyzed using the Pearson correlation test (Li et al., 2017). The IC50 of 265 small molecules and its corresponding gene expression were isolated from Genomics of Drug Sensitivity in Cancer (GDSC) (Yang et al., 2013). The correlation between cuproptosis-related prognostic gene expression and drug IC50 was analyzed using the Pearson correlation test. Moreover, correlation analysis between cuproptosis-related prognostic gene expression and TMB/MSI was performed using the Spearman method.
Cell Culture and qRT-PCR
The cells and culture conditions were as follows: Human endometriosis cell line hEM15A (Dulbecco’s modified eagle medium (DMEM) containing 10% FBS) and UCEC cell lines KLE (MEM containing 10% FBS) and ANC3A (DMEM containing 10% FBS) in a 5% CO2 incubator at 37°C. The total RNA of these cells was isolated from TRIzol (Invitrogen) and RNA extraction kits. The PrimeScript RT and miRNA PrimescripTM RT kits (E047-01B; Novoprotein) were used to reverse transcribe RNA into complementary DNA (cDNA). This was followed by qRT-PCR for gene expression analysis. The mRNA expression levels were normalized against glyceraldehyde 3-phosphate dehydrogenase expression.
lncRNA–miRNA–mRNA ceRNA Network Analysis
The miRNA target of CDKN2A was explored with StarBase (http://starbase.sysu.edu.cn/) (Li et al., 2014), TargetScan (https://www.targetscan.org/) (Lin et al., 2021), and miRWalk (http://mirwalk.umm.uni-heidelberg.de/) (Sticht et al., 2018). Furthermore, the lncRNA target interacting with the miRNA target was explored with lncBase (http://carolina.imis.athena-innovation.gr/) (Karagkouni et al., 2020) and StarBase (http://starbase.sysu.edu.cn/) (Li et al., 2014). Furthermore, Student t-test and KM test were performed to evaluate the expression and prognosis of miRNA and lncRNA in UCEC.
Results
Expression of Cuproptosis-Related Genes in Uterine Corpus Endometrial Carcinoma
Analysis expression revealed that 9 of 10 CRGs were differentially expressed in UCEC (Figure 1A). In a more specific sense, the expression of FDX1, LIAS, DLAT, and CDKN2A was upregulated, whereas the expression of LIPT1, DLD, PDHB, MTF1, and GLS were downregulated in UCEC versus normal tissues (Figure 1A, all p < 0.05). Further correlation analysis revealed that most of these CRGs were positively correlated with each other (Figure 1B).
FIGURE 1. Expression of cuproptosis-related genes (CRGs) in uterine corpus endometrial carcinoma (UCEC). (A) The mRNA level of CRG in UCEC versus normal tissues. (B) The correlation between each member of CRG in UCEC. *p < 0.05; **p < 0.01; and ***p < 0.001;—p > 0.05.
Mutation Landscape and Functional Enrichment Analysis
Figures 2A,B show the mutation landscape of CRGs in UCEC, revealing that 80 UCEC samples had a genetic mutation. Among these 10 CRGs, MTF1 was the gene with the highest frequency of mutation followed by PDHA1 and GLS (Figures 2A,B). In variant classification and SNV class, missense mutation and C > T ranked at the top, respectively (Figure 2B). CNV analysis suggested that more than half of 10 PRGs had homozygous amplification, whereas FDX1 and DLAT had a widespread homozygous deletion (Figure 2C). In functional enrichment analysis, the result of GO analysis was shown in Figure 3A, which reveals that these CRGs were enriched in sulfur compound biosynthetic process, acetyl-CoA biosynthetic and metabolic process, mitochondrial matrix, and metal cluster binding. Moreover, these CRGs were enriched in the TCA cycle, pyruvate metabolism, glycolysis, central carbon metabolism in cancer, HIF-1 signaling pathway, and miRNAs in cancer in the KEGG pathway analysis (Figure 3B).
FIGURE 2. Mutation landscape of cuproptosis-related genes (CRGs) in uterine corpus endometrial carcinoma (UCEC). (A,B) Single nucleotide variation analysis of CRG in UCEC. (C) Copy number variation analysis of CRG in UCEC.
FIGURE 3. Enrichment analysis of cuproptosis-related genes in uterine corpus endometrial carcinoma. (A) Enriched items in gene ontology analysis. (B) Enriched items in Kyoto Encyclopedia of Genes and Genomes. BP, biological process; CC, cellular component; MF molecular function.
Development of a Cuproptosis-Related Prognostic Signature in Uterine Corpus Endometrial Carcinoma
OS, PFS, and DFS analyses were performed to identify a potential cuproptosis-related prognostic marker for UCEC. The result of OS analysis is shown in Figure 4A, which reveals that UCEC patients with high level of CDKN2A (p < 0.0001, HR = 2.32) and GLS (p = 0.001, HR = 2.03) had a poor OS rate. As shown in Figure 4B, PFS analysis suggested a poor OS rate in UCEC patients with a high level of CDKN2A (p < 0.0001, HR = 2.12) and GLS (p = 0.019, HR = 1.53). Furthermore, the result of PFS analysis demonstrated that high expression of CDKN2A (p < 0.0001, HR = 3.58), GLS (p < 0.0001, HR = 3.23), and LIPT1 (p = 0.013, HR = 1.95) had a poor prognosis in UCEC (Figure 4C). All in all, CDKN2A, GLS, and LIPT1 were potential prognostic biomarkers for UCEC. Based on these potential prognosis biomarkers, LASSO Cox regression analysis was performed and a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1) was constructed. Figures 5A,B reveal the coefficient and partial likelihood deviance of cuproptosis-related prognostic signature. After calculating the risk score of UCEC patients using the formula ((0.2725)*CDKN2A+(0.416)*GLS+(0.0143)*LIPT1), UCEC patients were divided into high- and low-risk groups and the risk score, survival status, and gene expression of prognostic signature was shown in Figure 5C. Moreover, the data revealed that UCEC patients with high risk scores had a poor prognosis (Figure 5D, p = 2.05e−6, HR = 4.45), with an area under the curve of 0.782 and 0.764 in a 3- and 5-year ROC curve (Figure 5E), demonstrating a good performance of this signature in predicting the prognosis of UCEC patients. Considering clinicopathologic characters and prognostic signature, univariate and multivariate analyses were performed to further identify the prognostic factors, which demonstrated CDKN2A, LIPT1, and pTNM stage as independent prognosis factors for UCEC patients (Figures 6A,B). A predictive nomogram was constructed based on these results, suggesting the good performance of this predictive nomogram in the 3‐ and 5‐year survival rates compared with an ideal model (Figures 6C,D).
FIGURE 4. Prognostic significance of cuproptosis-related genes (CRGs) in uterine corpus endometrial carcinoma (UCEC). (A) Result of overall survival analysis of CRGs in UCEC. (B) Result of progression-free survival analysis of CRGs in UCEC. (C) Result of disease-free survival analysis of CRGs in UCEC. OS, overall survival; PFS, progression-free survival; DFS, disease-free survival.
FIGURE 5. A cuproptosis-related prognostic signature in uterine corpus endometrial carcinoma (UCEC). (A,B) Coefficient and partial likelihood deviance of prognostic signature. (C) The risk score distribution, patient survival status, and cuproptosis-related gene expression profile of prognostic signature. (D,E) Survival curve of UCEC patients with high/low-risk score and ROC curve of prognostic signature.
FIGURE 6. Development of a predictive nomogram. (A,B) Univariate and multivariate Cox regression considering clinical parameters and cuproptosis-related prognostic biomarkers. (C,D) Predictive nomogram to predict the 3- and 5-year survival of uterine corpus endometrial carcinoma patients. Calibration curve for the survival nomogram model in the discovery group. A dashed diagonal line represents the ideal nomogram.
Risk Model Genes Analysis
The correlation between cuproptosis-related prognostic genes and tumor grade as well as pTNM stage was also analyzed. The expression of CDKN2A (Figure 7A, p = 1.5e−5), GLS (Figure 7B, p = 3.1e−6) increased as the pTNM stage increased in UCEC. Moreover, the expression of CDKN2A (Figure 7D, p = 3.5e−10) and GLS (Figure 7E, p = 5.2e−17) in UCEC patients with low tumor grade were higher than in those with high tumor grade. However, there is no significant correlation between LIPT1 expression and pTNM stage (Figure 7C, p = 0.62) and tumor grade (Figure 7F, p = 0.25). In immune infiltration analysis, CDKN2A expression increased as the abundance of CD8+ T cells decreased and the abundance of CD4+ T cells and neutrophils increased (Figure 8A). The result suggested a positive correlation between GLS expression and the abundance of CD8+ T cells, neutrophils, and dendritic cells (Figure 8B). Moreover, the result found that LIPT1 expression was positively correlated with the abundance of CD8+ T cells, macrophages, and neutrophils (Figure 8C). Somatic cell copy number alteration of CDKN2A, GLS, and LIPT1 could inhibit the infiltration level of some immune cells (Figures 8D–F). Drug sensitivity analysis was performed to explore the potential of CDKN2A, GLS, and LIPT1 as drug scanning targets for UCEC, showing that high CDKN2A expression and low GLS expression were positively correlated with drug resistance of GDSC (Figure 8G). Increasing pieces of evidence suggested TMB and MSI as predictive markers for tumor immunotherapy efficacy in cancer (Liu et al., 2019; Samstein et al., 2019). In the current study, the data indicated a negative correlation between TMB score and CDNK2A expression (Figure 9A, p = 7.14e−17) but not the expression of GLS (Figure 9B, p = 0.954) and LIPT1 (Figure 9C, p = 0.536). Moreover, the MSI score of UCEC patients increased as CDNK2A expression (Figure 9D, p = 2.74e−7) decreased and the expression of GLS (Figure 9E, p = 4.23e−6) and LIPT1 (Figure 9F, p = 0.003) increased. Above all, among these three genes, CDNK2A was the most possible gene involved in the progression of UCEC.
FIGURE 7. Expression of cuproptosis-related prognostic biomarkers in different groups of uterine corpus endometrial carcinoma (UCEC) patients. (A–C) Expression of CDKN2A, GLS, and LIPT1 in UCEC patients in different pTNM stages. (D–F) Expression of CDKN2A, GLS, and LIPT1 in UCEC patients in different tumor grades.
FIGURE 8. Immune cell infiltration and drug sensitivity analysis of cuproptosis-related prognostic biomarkers in uterine corpus endometrial carcinoma (UCEC). (A–C) Correlation between the expression of CDKN2A, GLS, and LIPT1 and immune cell infiltration in UCEC. (D–F) Correlation between copy number variation of CDKN2A, GLS, and LIPT1 and immune cell infiltration in UCEC. (G) Correlation between the expression of CDKN2A, GLS, and LIPT1 and drug sensitivity in UCEC. GDSC, genomics of drug sensitivity in cancer.
FIGURE 9. TMB and MSI of cuproptosis-related prognostic biomarkers in uterine corpus endometrial carcinoma (UCEC). (A-C) Correlation between the expression of CDKN2A, GLS, and LIPT1 and TMB score in UCEC. (D–F) Correlation between the expression of CDKN2A, GLS, and LIPT1 and MSI score in UCEC. TMB, tumor mutation burden; MSI, microsatellite instability.
Validation of the Expression of Risk Model Genes
As expected, the expression of CDKN2A (Supplementary Figure S2A) was upregulated, whereas the expression of GLS (Supplementary Figure S2B) and LIPT1 (Supplementary Figure S2C) was downregulated in UCEC cells versus normal human endometriosis cells. These were consistent with the above results.
lncRNA–miRNA–mRNA ceRNA Network Analysis
To further clarify the hub gene-associated mechanism in UCEC, a lncRNA–miRNA–mRNA ceRNA network analysis was performed. Combined with the predicted result of StarBase, miRWalk, and TargetScan, a total of four miRNAs (miR-125a-5p, miR-125b-5p, miR-4319, and miR-485-5p) were identified as the most potential miRNA targets of CDNK2A (Figure 10A). However, only two of four miRNAs were differently expressed in UCEC (Figures 10B,C). To be more specific, the expression of miR-125a-5p (Figure 10B, p = 0.014) and miR-125b-5p (Figure 10C, p = 3e−9) were downregulated in UCEC versus normal tissues. Further OS revealed that UCEC patients with high miR-125a-5p levels had better survival (Figure 10D, p = 0.018). Thus, miR-125a-5p was the most possible miRNA target of CDKN2A. The lncRNA target of miR-125a-5p was further explored. Combined with the predicted result of lncBase V3 and StarBase, lncRNA XIST was suggested as the most potential target of miR-125a-5p (Figure 10E). Further analysis revealed that XIST expression was downregulated in UCEC (Figure 10F, p = 0.004) and high XIST levels had a better survival (Figure 10G, p = 0.048). Thus, lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis may play a vital role in the progression of UCEC and further in vivo and in vitro studies should be conducted to verify this hypothesis.
FIGURE 10. lncRNA–miRNA–mRNA regulatory axis in uterine corpus endometrial carcinoma (UCEC). (A) The miRNA targets predicted by TargetScan, miRWalk, and StarBase. (B–D) The expression and prognosis significance of miRNA target in UCEC. (F) The lncRNA targets predicted by lncBase and StarBase. (F,G) Expression and prognosis significance of lncRNA XIST in UCEC.
Discussion
Characterized by biological invasiveness and potent metastasis, UCEC is still a major problem affecting public women’s health (Li et al., 2021). Moreover, the disadvantage of it was that the incidence and mortality of UCEC are increasing (Amant et al., 2005). There are few satisfactory biomarkers or approaches that could predict the prognosis of UCEC patients accurately. Since the significant role exerted by cuproptosis in the oncogenesis, severity, and development of malignancy, it is important to explore the significance of CRGs in the development and prognosis of UCEC.
Gene expression profile suggested that the expression of FDX1, LIAS, DLAT, and CDKN2A were upregulated, whereas the expression of LIPT1, DLD, PDHB, MTF1, and GLS were downregulated in UCEC versus normal tissues. These data were consistent with the results of previous studies. A previous study revealed that CDKN2A was upregulated in UCEC and it might be implicated in the pathogenesis of UCEC (Su et al., 2015). However, no study had been conducted to elucidate the expression and function of the other CRGs in UCEC. Correlation analysis demonstrated that the majority of CRGs were positively correlated with each other. GLS and PDHA1 played a synergistic role in promoting greater glutamine dependence in prostate cancer patients (Li et al., 2016). Moreover, PDHA1, PDHB, DLAT, and DLD played a synergistic role in pyruvate dehydrogenase complex deficiency (Inui et al., 2022). GO and KEGG analyses revealed that these CRGs were enriched in the TCA cycle, glycolysis, and HIF-1 signaling pathway. These pathways have been found to be involved in the progression and therapy of UCEC. By regulating the HIF-1 signaling pathway, AGR2-induced glucose metabolism could promote the progression of UCEC (Gong et al., 2020). HIF-1 signaling pathway-related hypoxia and hyperglycemia were associated with metformin resistance in endometrial cancer (Sivalingam et al., 2020). TCA cycle and its enzyme components play a vital role in basal cell metabolism, thus affecting tumor cell proliferation and invasion in UCEC (Iplik et al., 2018).
Further prognosis revealed that UCEC patients with high levels of CDKN2A and GLS had a poor OS rate. PFS analysis suggested a poor OS rate in UCEC patients with a high level of CDKN2A and GLS. Moreover, the result of the PFS analysis demonstrated that high expression of CDKN2A, GLS, and LIPT1 had a poor prognosis in UCEC. All in all, CDKN2A, GLS, and LIPT1 were potential prognostic biomarkers for UCEC. Another study also suggested CDKN2A as a prognostic biomarker in human endometrial cancer (Zhang et al., 2020). Chen et al. (2021b) found that LIPT1 was associated with survival in urothelial cancer. Moreover, GLS differentially could modulate the prognosis of human cancer (Saha et al., 2019). Furthermore, LASSO Cox regression analysis was performed and identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1), which had a good performance in predicting the prognosis of UCEC patients. As far as we know, the TNM system was the most reliable predictive method for the prognosis of UCEC patients. Although some prognostic signatures had been identified for the prognosis of UCEC patients (Hu et al., 2021; Weijiao et al., 2021), these prognostic signatures had not been used for the prognosis of UCEC patients clinically. The current study identified the first cuproptosis-related prognostic signature for UCEC patients, which may be used for the prognosis prediction of UCEC patients clinically.
In immune infiltration analysis, the result revealed that the expression of CDKN2A, GLS, and LIPT1 was positively correlated with the abundance of certain immune cells, including CD8+ T cells and neutrophils. Similar results had been reported in other types of cancers. CDKN2A expression was positively correlated with infiltrating levels into CD8+ T cells, CD4+ T cells, and neutrophils in hepatocellular carcinoma (Luo et al., 2021). Another integrative bioinformatics analysis revealed that CDKN2A shows a significant correlation with immune signatures in multiple myeloma, including CD4+ regulatory T cells, T cell exhaustion, and neutrophils (Tuerxun et al., 2022). A previous study revealed that STAT5 and CDKN2A/CDKN2B could promote the proliferation and terminal differentiation of CD8+ T cells (Grange et al., 2015).
Another important finding of the current study was that we identified that lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis may be involved in the progression of UCEC. Increasing pieces of evidence had revealed that lncRNA XIST played a vital role in many biological processes in human cancer (Liu et al., 2018; Zhu et al., 2018; Yang et al., 2020b). Overexpression of miR-125a-5p could inhibit tumor cell apoptosis in UCEC (Shi et al., 2019). CDKN2A was involved in tumor cell growth and drug resistance in human breast cancer (Lubecka et al., 2018). A previous study revealed that miR-125a-5p/CDKN2A regulatory axis exerts a vital role in prognostic value in cervical cancer (Wang et al., 2021). The development and progression of UCEC was a complex process with continuous progression and multiple steps. The specific molecular mechanism of UCEC had not been fully elucidated. In the current study, a lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis was identified for UCEC and it has not been studied before, which provided another theoretical basis for the further study of the mechanism of oncogenesis and progression of UCEC.
The current study also had some limitations. Cuproptosis-related prognostic signature should be verified by another dataset. Further in vivo and in vitro studies should be conducted to verify the lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis.
Conclusion
In conclusion, the current study identified a cuproptosis-related prognostic signature including these three prognostic biomarkers (CDKN2A, GLS, and LIPT1) for UCEC. The ceRNA network also identified that lncRNA XIST/miR-125a-5p/CDKN2A regulatory axis may be involved in the progression of UCEC. Further in vivo and in vitro studies should be conducted to verify these results.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
YC designed the study, performed data analysis work, and aided in writing and editing the manuscript.
Funding
This work was financially supported through grants from the Program for Medical Key Discipline of Hangzhou-Gynecology of Traditional Chinese Medicine (2020SJZDXK02).
Conflict of Interest
The author declares 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.912037/full#supplementary-material
References
Amant, F., Moerman, P., Neven, P., Timmerman, D., Van Limbergen, E., and Vergote, I. (2005). Endometrial Cancer. Lancet 366 (9484), 491–505. doi:10.1016/s0140-6736(05)67063-8
Bai, H., Wang, Y., Liu, H., and Lu, J. (2021). Development of a Four-mRNA Expression-Based Prognostic Signature for Cutaneous Melanoma. Front. Genet. 12, 680617. doi:10.3389/fgene.2021.680617
Chen, L., Zou, W., Zhang, L., Shi, H., Li, Z., and Ni, C. (2021). ceRNA Network Development and Tumor-Infiltrating Immune Cell Analysis in Hepatocellular Carcinoma. Med. Oncol. 38 (7), 85. doi:10.1007/s12032-021-01534-6
Chen, Y., Xu, T., Xie, F., Wang, L., Liang, Z., Li, D., et al. (2021). Evaluating the Biological Functions of the Prognostic Genes Identified by the Pathology Atlas in Bladder Cancer. Oncol. Rep. 45 (1), 191–201. doi:10.3892/or.2020.7853
Cui, Y.-Q., Xiang, Y., Meng, F., Ji, C.-H., Xiao, R., Li, J.-P., et al. (2021). ALDH2 Promotes Uterine Corpus Endometrial Carcinoma Proliferation and Construction of Clinical Survival Prognostic Model. Aging 13 (20), 23588–23602. doi:10.18632/aging.203605
Gong, W., Ekmu, B., Wang, X., Lu, Y., and Wan, L. (2020). AGR2-induced Glucose Metabolism Facilitated the Progression of Endometrial Carcinoma via Enhancing the MUC1/HIF-1α Pathway. Hum. Cell. 33 (3), 790–800. doi:10.1007/s13577-020-00356-4
Grange, M., Giordano, M., Mas, A., Roncagalli, R., Firaguay, G., Nunes, J. A., et al. (2015). Control of CD8 T Cell Proliferation and Terminal Differentiation by Active STAT5 and CDKN2A/CDKN2B. Immunology 145 (4), 543–557. doi:10.1111/imm.12471
Hu, Y., Zheng, M., Zhang, D., Gou, R., Liu, O., Wang, S., et al. (2021). Identification of the Prognostic Value of a 2-gene Signature of the WNT Gene Family in UCEC Using Bioinformatics and Real-World Data. Cancer Cell. Int. 21 (1), 516. doi:10.1186/s12935-021-02215-0
Inui, T., Wada, Y., Shibuya, M., Arai-Ichinoi, N., Okubo, Y., Endo, W., et al. (2022). Intravenous Ketogenic Diet Therapy for Neonatal-Onset Pyruvate Dehydrogenase Complex Deficiency. Brain Dev. 44 (3), 244–248. doi:10.1016/j.braindev.2021.11.005
Iplik, E. S., Catmakas, T., and Cakmakoglu, B. (2018). A New Target for the Treatment of Endometrium Cancer by Succinic Acid. Cell. Mol. Biol. (Noisy-le-grand) 64 (1), 60–63. doi:10.14715/cmb/2018.64.1.11
Kahlson, M. A., and Dixon, S. J. (2022). Copper-induced Cell Death. Science 375 (6586), 1231–1232. doi:10.1126/science.abo3959
Kang, Z., Li, W., Yu, Y.-H., Che, M., Yang, M.-L., Len, J.-J., et al. (2021). Identification of Immune-Related Genes Associated with Bladder Cancer Based on Immunological Characteristics and Their Correlation with the Prognosis. Front. Genet. 12, 763590. doi:10.3389/fgene.2021.763590
Karagkouni, D., Paraskevopoulou, M. D., Tastsoglou, S., Skoufos, G., Karavangeli, A., Pierros, V., et al. (2020). DIANA-LncBase V3: Indexing Experimentally Supported miRNA Targets on Non-coding Transcripts. Nucleic Acids Res. 48 (D1), D101–d110. doi:10.1093/nar/gkz1036
Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucleic Acids Res. 42 (Database issue), D92–D97. doi:10.1093/nar/gkt1248
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307
Li, Y., Li, X., Li, X., Zhong, Y., Ji, Y., Yu, D., et al. (2016). PDHA1 Gene Knockout in Prostate Cancer Cells Results in Metabolic Reprogramming towards Greater Glutamine Dependence. Oncotarget 7 (33), 53837–53852. doi:10.18632/oncotarget.10782
Li, Y., Yu, S., Li, Y., Liang, X., Su, M., and Li, R. (2021). Medical Significance of Uterine Corpus Endometrial Carcinoma Patients Infected with SARS-CoV-2 and Pharmacological Characteristics of Plumbagin. Front. Endocrinol. 12, 714909. doi:10.3389/fendo.2021.714909
Lin, W., Chen, Y., Wu, B., Chen, Y., and Li, Z. (2021). Identification of the Pyroptosis-related P-rognostic G-ene S-ignature and the A-ssociated R-egulation axis in L-ung A-denocarcinoma. Cell. Death Discov. 7 (1), 161. doi:10.1038/s41420-021-00557-2
Liu, H., Deng, H., Zhao, Y., Li, C., and Liang, Y. (2018). LncRNA XIST/miR-34a axis Modulates the Cell Proliferation and Tumor Growth of Thyroid Cancer through MET-Pi3k-AKT Signaling. J. Exp. Clin. Cancer Res. 37 (1), 279. doi:10.1186/s13046-018-0950-9
Liu, L., Bai, X., Wang, J., Tang, X.-R., Wu, D.-H., Du, S.-S., et al. (2019). Combination of TMB and CNA Stratifies Prognostic and Predictive Responses to Immunotherapy across Metastatic Cancer. Clin. Cancer Res. 25 (24), 7413–7423. doi:10.1158/1078-0432.ccr-19-0558
Lubecka, K., Kaufman-Szymczyk, A., and Fabianowska-Majewska, K. (2018). Inhibition of Breast Cancer Cell Growth by the Combination of Clofarabine and Sulforaphane Involves Epigenetically Mediated CDKN2A Upregulation. Nucleosides, Nucleotides Nucleic Acids 37 (5), 280–289. doi:10.1080/15257770.2018.1453075
Luo, J. P., Wang, J., and Huang, J. H. (2021). CDKN2A Is a Prognostic Biomarker and Correlated with Immune Infiltrates in Hepatocellular Carcinoma. Biosci. Rep. 41 (10), BSR20211103. doi:10.1042/BSR20211103
Maag, J. L. V. (2018). Gganatogram: An R Package for Modular Visualisation of Anatograms and Tissues Based on Ggplot2. F1000Res 7, 1576. doi:10.12688/f1000research.16409.2
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
Oliveri, V. (2022). Selective Targeting of Cancer Cells by Copper Ionophores: An Overview. Front. Mol. Biosci. 9, 841814. doi:10.3389/fmolb.2022.841814
Qi, X., Zhang, D.-H., Wu, N., Xiao, J.-H., Wang, X., and Ma, W. (2015). ceRNA in Cancer: Possible Functions and Clinical Implications. J. Med. Genet. 52 (10), 710–718. doi:10.1136/jmedgenet-2015-103334
Ruiz, L. M., Libedinsky, A., and Elorza, A. A. (2021). Role of Copper on Mitochondrial Function and Metabolism. Front. Mol. Biosci. 8, 711227. doi:10.3389/fmolb.2021.711227
Saha, S. K., Islam, S. M. R., Abdullah-Al-Wadud, M., Islam, S., Ali, F., and Park, K. S. (2019). Multiomics Analysis Reveals that GLS and GLS2 Differentially Modulate the Clinical Outcomes of Cancer. J. Clin. Med. 8 (3), 355. doi:10.3390/jcm8030355
Samstein, R. M., Lee, C.-H., Shoushtari, A. N., Hellmann, M. D., Shen, R., Janjigian, Y. Y., et al. (2019). Tumor Mutational Load Predicts Survival after Immunotherapy across Multiple Cancer Types. Nat. Genet. 51 (2), 202–206. doi:10.1038/s41588-018-0312-8
Shi, F., Wang, T., Liu, Z., Zhang, Y., Wang, J., Zhang, K., et al. (2019). LncRNA miR143HG Up-Regulates P53 in Endometrial Carcinoma by Sponging miR-125a. Cmar Vol. 11, 10117–10123. doi:10.2147/cmar.s222215
Sivalingam, V. N., Latif, A., Kitson, S., McVey, R., Finegan, K. G., Marshall, K., et al. (2020). Hypoxia and Hyperglycaemia Determine Why Some Endometrial Tumours Fail to Respond to Metformin. Br. J. Cancer 122 (1), 62–71. doi:10.1038/s41416-019-0627-y
Sticht, C., De La Torre, C., Parveen, A., and Gretz, N. (2018). miRWalk: An Online Resource for Prediction of microRNA Binding Sites. PLoS One 13 (10), e0206239. doi:10.1371/journal.pone.0206239
Su, L., Wang, H., Miao, J., and Liang, Y. (2015). Clinicopathological Significance and Potential Drug Target of CDKN2A/p16 in Endometrial Carcinoma. Sci. Rep. 5, 13238. doi:10.1038/srep13238
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tomczak, K., Czerwińska, P., and Wiznerowicz, M. (2015). The Cancer Genome Atlas (TCGA): an Immeasurable Source of Knowledge. Contemp. Oncol. Pozn. 19 (1a), A68–A77. doi:10.5114/wo.2014.47136
Tsvetkov, P., Coy, S., Petrova, B., Dreishpoon, M., Verma, A., Abdusamad, M., et al. (2022). Copper Induces Cell Death by Targeting Lipoylated TCA Cycle Proteins. Science 375 (6586), 1254–1261. doi:10.1126/science.abf0529
Tuerxun, N., Wang, J., Qin, Y.-t., Zhao, F., Wang, H., Qu, J.-h., et al. (2022). Identification of Key Genes and miRNA-mRNA Regulatory Networks Associated with Bone Marrow Immune Microenvironment Regulations in Multiple Myeloma by Integrative Bioinformatics Analysis. Hematology 27 (1), 506–517. doi:10.1080/16078454.2022.2068873
Wang, L.-X., Wan, C., Dong, Z.-B., Wang, B.-H., Liu, H.-Y., and Li, Y. (2019). Integrative Analysis of Long Noncoding RNA (lncRNA), microRNA (miRNA) and mRNA Expression and Construction of a Competing Endogenous RNA (ceRNA) Network in Metastatic Melanoma. Med. Sci. Monit. 25, 2896–2907. doi:10.12659/msm.913881
Wang, P., Zeng, Z., Shen, X., Tian, X., and Ye, Q. (2020). Identification of a Multi-RNA-Type-Based Signature for Recurrence-free Survival Prediction in Patients with Uterine Corpus Endometrial Carcinoma. DNA Cell. Biol. 39 (4), 615–630. doi:10.1089/dna.2019.5148
Wang, T., Zhang, X.-D., and Hua, K.-Q. (2021). A ceRNA Network of BBOX1-AS1-hsa-miR-125b-5p/hsa-miR-125a-5p-CDKN2A Shows Prognostic Value in Cervical Cancer. Taiwan. J. Obstetrics Gynecol. 60 (2), 253–261. doi:10.1016/j.tjog.2020.12.006
Weijiao, Y., Fuchun, L., Mengjie, C., Xiaoqing, Q., Hao, L., Yuan, L., et al. (2021). Immune Infiltration and a Ferroptosis-Associated Gene Signature for Predicting the Prognosis of Patients with Endometrial Cancer. Aging 13 (12), 16713–16732. doi:10.18632/aging.203190
Wu, X., Sui, Z., Zhang, H., Wang, Y., and Yu, Z. (2020). Integrated Analysis of lncRNA-Mediated ceRNA Network in Lung Adenocarcinoma. Front. Oncol. 10, 554759. doi:10.3389/fonc.2020.554759
Yang, J., Li, H., Hu, S., and Zhou, Y. (2020). ACE2 Correlated with Immune Infiltration Serves as a Prognostic Biomarker in Endometrial Carcinoma and Renal Papillary Cell Carcinoma: Implication for COVID-19. Aging 12 (8), 6518–6535. doi:10.18632/aging.103100
Yang, W., Soares, J., Greninger, P., Edelman, E. J., Lightfoot, H., Forbes, S., et al. (2013). Genomics of Drug Sensitivity in Cancer (GDSC): a Resource for Therapeutic Biomarker Discovery in Cancer Cells. Nucleic Acids Res. 41 (Database issue), D955–D961. doi:10.1093/nar/gks1111
Yang, X., Zhang, S., He, C., Xue, P., Zhang, L., He, Z., et al. (2020). METTL14 Suppresses Proliferation and Metastasis of Colorectal Cancer by Down-Regulating Oncogenic Long Non-coding RNA XIST. Mol. Cancer 19 (1), 46. doi:10.1186/s12943-020-1146-4
Zhang, J., Wang, Z., Zhao, R., An, L., Zhou, X., Zhao, Y., et al. (2020). An Integrated Autophagy-Related Gene Signature Predicts Prognosis in Human Endometrial Cancer. BMC Cancer 20 (1), 1030. doi:10.1186/s12885-020-07535-4
Keywords: uterine corpus endometrial carcinoma, cuproptosis, prognostic signature, CDKN2A, immunotherapy
Citation: Chen Y (2022) Identification and Validation of Cuproptosis-Related Prognostic Signature and Associated Regulatory Axis in Uterine Corpus Endometrial Carcinoma. Front. Genet. 13:912037. doi: 10.3389/fgene.2022.912037
Received: 03 April 2022; Accepted: 13 June 2022;
Published: 22 July 2022.
Edited by:
Jinhui Liu, Nanjing Medical University, ChinaReviewed by:
Simin Li, Southern Medical University, ChinaKaushik Kumar Dey, St. Jude Children’s Research Hospital, United States
Copyright © 2022 Chen. 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: Yun Chen, ZmlvbmNlZUAxNjMuY29t