- 1Department of Gynecology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
- 2Department of Biostatistics, School of Public Heath, Nanjing Medical University, Nanjing, China
Background: As the fourth most common malignant tumors in women, uterine corpus endometrial carcinoma (UCEC) requires novel and reliable biomarkers for prognosis prediction to improve the overall survival. Oxidative phosphorylation (OXPHOS) is found to be strongly correlated with the progression of tumor. Here, we aimed to construct an OXPHOS-related and immune microenvironment prognostic signature to stratify UCEC patients for optimization of treatment strategies.
Method: Prognosis-associated OXPHOS-related differentially expressed genes were identified by multivariable Cox regression from TCGA–UCEC cohort. Based on the candidate genes, an OXPHOS-related prognostic signature was constructed by the train set data and verified by the entire set. When integrated with relevant clinical characteristics, a nomogram was also created for clinical application. Through comparison of tumor microenvironment between different risk groups, the underlying mechanism of the model and the inner correlation between immune microenvironment and energy metabolism were further investigated.
Results: An OXPHOS-related signature containing ATP5IF1, COX6B1, FOXP3, and NDUFB11 was constructed and had better predictive ability compared with other recently published signatures in UCEC. Patients with lower risk score showed higher immune cell infiltration, higher ESTIMATE score (p = 2.808E−18), lower tumor purity (p = 2.808E−18), higher immunophenoscores (IPSs) (p < 0.05), lower expression of mismatch repair (MMR) proteins (p < 0.05), higher microsatellite instability (MSI), lower expression of markers of N6-methyladenosine (m6A) mRNA methylation regulators, higher tumor mutation burden (TMB) (p = 1.278E−9), and more sensitivity to immune checkpoint blockade (ICB) (p < 0.001) and chemotherapy drugs, thus, possessing improved prognosis.
Conclusion: An OXPHOS-related and immune microenvironment prognostic signature classifying EC patients into different risk subsets was constructed in our study, which could be used to predict the prognosis of patients and help to select a specific subset of patients who might benefit from immunotherapy and chemotherapy, thus, improving the overall survival rate of UCEC. These findings may contribute to the discovery of novel and robust biomarkers or target therapy in UCEC and give new insights into the molecular mechanism of tumorigenesis and progression of UCEC.
1 Introduction
Uterine corpus endometrial carcinoma (UCEC) is the fourth most common malignancy among women in the United States, with a trend of increasing morbidity and mortality worldwide (Cancer Genome Atlas Research et al., 2013). According to the GLOBOCAN, there would be an estimated 382,069 (2.1% of the total cancer cases) new cases, with 89,929 (0.9% of total cancer deaths) deaths in 2018 (Bray et al., 2018). Due to the advance in techniques for early diagnosis, the patients diagnosed at an early stage account for approximately 75% of UCEC, and most of the patients in stage I can be nearly cured by surgery, with a 5-years overall survival rate of 47–69%. In comparison, advanced-stage patients (stage III or IV) have a poor prognosis due to high risk for recurrence and limited therapeutic strategies, with a 5-years overall survival rate of 15–17% (Lee et al., 2017). Therefore, attention should be paid to discovering novel and reliable biomarkers for prognosis prediction, which can also work as a sensitive classifier for a specific subset of UCEC patients who will benefit from immunotherapy and have improved survival.
Recently, increasing attention has been paid into the exploration of the relationship between cancer cell metabolic plasticity and migration and metastasis (Mosier et al., 2021). Cancer cells with a more invasive and distal metastasis phenotype are found to be strongly correlated with the upregulated expression of peroxisome proliferator-associated receptor gamma and coactivator 1-alpha, which are hallmarks of active mitochondrial biogenesis and oxidative phosphorylation (OXPHOS) (LeBleu et al., 2014). OXPHOS is supposed to be a potential biomarker of tumor progression (Commander et al., 2020). In addition, emerging evidence shows that OXPHOS inhibition can be a valuable target in cancer treatment (Ashton et al., 2018; Molina et al., 2018; Xu et al., 2020a; Cardenas et al., 2020). OXPHOS provides the required energy for cancer cells to strive, and cancer stem cells with primary or acquired resistance against chemotherapy or tyrosine kinase inhibitors are characterized with markedly enhanced OXPHOS dependency (Sica et al., 2020). Metformin and thiazolidinediones can inhibit the mitochondrial electron transport chain, which provides energy for ovarian cancer growth. Thus, interfering with the process of OXPHOS may be considered as a new target for cancer therapy (Nayak et al., 2018). With the development of molecular medicine and next-generation sequencing technology, OXPHOS has been employed in the cancer risk prediction model in lung adenocarcinoma (Xu et al., 2020b). OXPHOS has also attached great importance to UCEC. As a regulator participating in oxidative phosphorylation and glycolysis, PKM2 can function as a biomarker for malignant and premalignant endometrial lesions, and the presence of PKM2high tumor cells in UCEC tissue also indicates a poor prognosis (Lai et al., 2019). To the best of our knowledge, OXPHOS has not been included in the UCEC prognostic model in previous studies, which deserves to be explored for subgrouping of the patients with distinct prognosis for personal treatment and identifying potential novel therapeutic targets.
Here, we identified prognosis-associated OXPHOS-related differentially expressed genes (DEGs) from TCGA–UCEC cohort (Tomczak et al., 2015). Based on the candidate genes, we constructed an OXPHOS-related prognostic signature model using the train dataset and verified it using the entire dataset. To further investigate the underlying mechanism of the model and its correlation with immune microenvironment and energy metabolism, then between high- and low-risk groups, we compared tumor-infiltrating immune cells (TIICs), ESTIMATE score, tumor purity, immunophenoscores (IPSs), the expression of mismatch repair (MMR) proteins, mRNA expression-based stemness index (mRNAsi), microsatellite instability expression (MSI), the expression of markers of N6-methyladenosine (m6A) mRNA methylation regulators, tumor mutation burden (TMB) and the response to immune checkpoint blockade (ICB) and chemotherapy agents, which may provide evidence for clinical application of the signature for immunotherapy and chemotherapy option (Figure 1). The findings in this study not only provided a useful predictive model for UCEC but also offered evidence for mutual impacts of OXPHOS, immune microenvironment, and tumor malignant biology.
2 Materials and Methods
2.1 Data Collection
The gene expression profiles were collected from TCGA–UCEC database (https://tcga-data.nci.nih.gov/tcga/), which contained 552 tumor and 23 normal tissues. Meanwhile, the corresponding clinical and pathological information was also obtained from TCGA database, including patient age, tumor grade, tumor stage, and histological type of each sample. After removing the samples with unknown survival time, <30 days, and no survival status, we finally integrated the transcriptome and complete clinical data, enrolling a total of 511 patients in this study (Table 1).
2.2 Oxidative Phosphorylation-Related Survival-Related Differentially Expressed Genes and Functional Annotation
We used the R package limma to define the differentially expressed genes among 289 OXPHOS-related genes acquired from the National Center for Biotechnology Information—GENE and the Metabolic Atlas database (Xu et al., 2020b). The genes with p < 0.05 and |log2 (fold change)| >1 were regarded as DEGs. We used the survival package to show the association between OXPHOS-related DEGs and survival time. We defined the DEGs with p less than 0.05 as survival-related DEGs (sDEGs). To explore the potential biological functions of the OXPHOS-related sDEGs, we implemented the ClusterProfiler package of R software for Gene Ontology (GO) term enrichment analysis (Ashburner et al., 2000), which included three categories: biological process (BP), molecular functions (MF), cellular components (CC), as well as the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (Kanehisa and Goto, 2000).
2.3 Construction of an Oxidative Phosphorylation-Related Prognostic Signature
We randomly assigned 256 patients as the train set and took the entire set for validation. First, we utilized the univariate Cox regression in the train set to identify candidate sDEGs with a significant level of p < 0.05. Then, we performed the least absolute shrinkage and selection operator (LASSO) Cox regression analysis by the glmnet package (Qu et al., 2018). Finally, we employed multivariate Cox proportional hazard regression analysis to construct an OXPHOS-related prognostic signature. The formula for calculating risk score was as follows:
We treated the median of risk score in the train set as the cutoff value to divide the validation data into the high- and low-risk groups and plotted the Kaplan–Meier curves, as well as the ROC curves by survminer, survivalROC package (Heagerty et al., 2000). To confirm the independence of the signature, a conjoint univariate and multivariate Cox analyses were conducted. In addition, we compared the subgroup survival analysis in tumor stage (stages I and II, stages III and IV), tumor grade (grades 1 and 2, grades 3 and 4), histological type (endometrial, mixed, and serous), and age (≤60, >60) to further evaluate the value of the signature combined with clinical factors in the survival prediction of UCEC. We used the rms package to construct a baseline nomogram for clinical application in UCEC patients (Iasonos et al., 2008).
2.4 RNA Isolation and Quantitative Real-Time PCR
UCEC and normal endometrial tissues were obtained from patients of the First Affiliated Hospital of Nanjing Medical University with the approval of the ethics committee and with informed consent. The clinicopathological parameters are shown in Supplementary Table S1. Total RNA was isolated from tissue samples with TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, United States), and the integrity of the extracted RNA was estimated by the Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, United States). RNA was reversely transcribed into cDNA using high-capacity reverse transcription kits (TaKaRa, Shiga, Japan), and quantitative real-time PCR (qRT-PCR) was conducted on Light Cycler 480 (Roche, Switzerland) using SYBR Green PCR Kit (Thermo Fisher Scientific) with the 2−ΔΔCt method, in which GAPDH was the endogenous control. All program procedures of qRT-PCR were performed according to the protocol of the manufacturer. Primer sequences for GAPDH and four OXPHOS-related genes are presented in Supplementary Table S2.
2.5 Gene Set Enrichment Analysis and Immunity Analyses
Gene set enrichment analysis (GSEA) was utilized to elucidate the molecular mechanisms of the OXPHOS-related sDEGs. We separated the samples in the entire set into the high- and low-risk groups based on scores calculated by the signature and then compared the enriched BP between the two groups. The value of p < 0.05 was regarded as the cutoff criterion (Zhu et al., 2016). Furthermore, we employed the ESTIMATE algorithm by the ESTIMATE package to calculate the stromal, immune, and estimate scores of each sample and evaluate the relationship between the OXPHOS-related prognostic signature and immunization (https://bioinformatics.mdanderson.org/public-software/estimate/) (Yoshihara et al., 2013).
2.6 Analysis of Tumor-Infiltrating Immune Cells
We employed CIBERSORT (https://cibersort.stanford.edu/) to explore the immune cell infiltration of each sample based on the RNA-seq data (Newman et al., 2015). To further analyze the association between risk score and tumor-infiltrating immune cells (TIICs), single-sample gene-set enrichment analysis (ssGSEA) was also performed to quantify the immune activity in different risk groups by exploring 29 immune-related genes (He et al., 2018).
2.7 Analysis of Immune Status Between Different Risk Groups
The immunogenicity of tumor is determined by the four dominant components containing immunomodulators, immunosuppressive cells, effector cells as well as MHC molecules. Immunophenoscores (IPSs) (including IPS, IPS-CTLA4, IPS-PD1-PD-L1-PD-L2, IPS-PD1-PD-L1-PD-L2-CTLA4 scores) are calculated based on the representative cell-type gene expression z-scores to evaluate and compare the potential response to immune checkpoint inhibitor between different risk groups (Charoentong et al., 2017). IPS data of each UCEC patient in the entire set was downloaded from The Cancer Immunome Atlas (TCIA) (https://tcia.at/home). Besides, the immunotherapy response could be predicted with high accuracy by an Immune Cell Abundance Identifier (ImmuCellAI) result-based model, which precisely estimates the abundance of 24 immune cell types including 18 T-cell subsets from gene expression data (http://bioinfo.life.hust.edu.cn/ImmuCellAI) (Wang et al., 2020a; Toh et al., 2021). Microsatellite instability (MSI) (Toh et al., 2021) and the expression of markers of N6-methyladenosine (m6A) mRNA methylation regulators (Wang et al., 2020a) of each sample were also collected from TCGA, which can predict potential response to immune checkpoint inhibitors (ICIs) and help explain the association between the signature and immunogenicity. Tumor cells sharing more similarities with stem cells represents cancer progression. Malta et al. developed an innovative mRNAsi algorithm to calculate the degree of similarities for UCEC patients in TCGA database (Malta et al., 2018). We also downloaded the mRNAsi information from the research for analyses.
2.8 Analysis of Tumor Mutation Burden
The mutation data collected from TCGA were analyzed by R package maftools (Robinson et al., 2017). The tumor mutation burden (TMB) was calculated by the formula: TMB = (total mutation/total covered bases) *
2.9 Potential Chemotherapeutic Response
The response to temsirolimus, roscovitine, AZD6244, PD.0325901, RDEA119, PF.02341066, AKT. inhibitor.VIII, BMS.509744, vinblastine, bryostatin.1, metformin, AZ628, nutlin.3a, bortezomib, bicalutamide, AZD6482, cytarabine, rapamycin, camptothecin, mitomycin. C, S. Trityl.L.cysteine, tipifarnib, parthenolide, sorafenib, and methotrexate, the 25 common chemo drugs, were predicted by the Genomics of Drug Sensitivity in Cancer (GDSC) (https://www.cancerrxgene.org/) to analyze the relationship between the signature and chemotherapeutic response (Yang et al., 2013). We used R package pRRophetic to estimate and compare the half-maximal inhibitory concentration (IC50) between different risk groups (Geeleher et al., 2014).
2.10 Statistical Analysis
All statistical analyses were applied by R version 3.6.1 (Package: limma, survival, ClusterProfiler, glmnet, survminer, survivalROC, rms, ESTIMATE, pRRophetic). Categorical variables were presented by counts and percentages. Continuous variables in normal distribution were analyzed using Student’s t-test and presented as mean ± standard deviation, while in abnormal distribution, they were presented as median (range). Multiple groups of continuous variables were analyzed by one-way ANOVA. The hazard ratio and 95% confidence interval were calculated to identify genes associated with overall survival. Unless with special explanation, p < 0.05 was considered statistically significant.
3 Results and Discussion
3.1 Results
3.1.1 Identification of Oxidative Phosphorylation-Related Differentially Expressed Genes
We conducted the differential expression analysis of 289 OXPHOS-related genes between tumor and normal endometrium (Xu et al., 2020b). Sixty-six OXPHOS-related DEGs were identified, including 54 upregulated genes and 12 downregulated genes (Supplementary Figures S1A,B and Supplementary Table S3). The GO term enrichment analysis presented that they were related to ATP metabolic process (p < 0.0005) and purine ribonucleoside triphosphate metabolic process (p < 0.0005) in the BP group, mitochondrial inner membrane (p < 0.0005), and organelle inner membrane (p < 0.0005) in the CC group, oxidoreductase activity, acting on NAD(P)H and electron transfer activity (p < 0.0005) in the MF group (Supplementary Figures S1C); the KEGG results manifested that they mainly took part in pathways of oxidative phosphorylation (p < 0.05), endocannabinoid signaling (p < 0.05), and HIF-1 signaling (p < 0.05) in cancer (Supplementary Figures S1D).
3.1.2 Establishment of Oxidative Phosphorylation-Related Prognostic Signature
Using univariate Cox regression from the train data, we defined seven OXPHOS-related sDEGs, including ATP5IF1, MRPL12, FOXP3, NDUFA13, ATP5F1E, and NDUFB11 (Supplementary Table S4). ATP5IF1 and FOXP3 were protective elements of UCEC, while the others were dangerous ones. Then we conducted LASSO Cox regression analysis and multivariate Cox regression analysis (Supplementary Figures S2A,B) to construct a reliable OXPHOS-related prognostic signature. Consequently, ATP5IF1, COX6B1, FOXP3, and NDUFB11 were included in the signature as key prognosis-associated ones (Supplementary Figures S2C). Furthermore, the results of the expression analyses of the four genes showed that all of them were highly expressed in tumor compared with normal tissue, but only ATP5IF1 and FOXP3 were highly associated with overall survival (OS) of UCEC patients (Supplementary Figures S3A–F). The qRT-PCR results of our own samples also demonstrated that the expression level of ATP5IF1 (p = 0.0311), COX6B1 (p = 0.0212), FOXP3 (p = 0.1082), and NDUFB11 (p = 0.0450) were significantly higher in UCEC tissues than in normal endometrial tissues (Supplementary Figures S3G–J). The formula to calculate the risk score of patients is shown as follows: risk score = (−0.011678 * expression value of ATP5IF1) + (0.0014649 * expression value of COX6B1) + (−0.309569 * expression value of FOXP3) + (0.0029791 * expression value of NDUFB11).
3.1.3 Validation of Prognostic Signature in Uterine Corpus Endometrial Carcinoma
We used the validation data to evaluate the performance of the risk model in survival prediction. Using the preceding equation, we calculated the risk score and divided the samples of the train and entire sets separately into high- and low-risk groups by taking the median risk score in the train set as the cutoff line. The distribution of risk score (Figure 2A), survival status of each sample (Figure 2B), and the key gene expression profiles in each set (Figure 2C) are exhibited in Figure 2. By comparison with the low-risk group, Kaplan–Meier survival curves revealed that the OS of the high-risk group was significantly worse (log rank test, p = 1.35E−4) (Figure 2D). Besides, a survival ROC curve analysis was conducted, and the area under the curve (AUC) at 1 year in the train set was 0.748 (Figure 2E), which verified the predictive value of the signature, and the AUC value in the entire set also suggested a similar potential for prediction of survival with the above results (Figure 2E). The PCA plots indicated that patients in different risk groups tended to distribute differently (Figure 2F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g002.jpg)
FIGURE 2. Risk score analysis, time-dependent ROC analysis, Kaplan–Meier analysis, and PCA plots for the validation of the oxidative phosphorylation (OXPHOS)-related prognostic signature in the train and entire sets based on overall survival (OS). (A) Rank of risk score and distribution of groups based on the signature. Patients in each group were divided into the high- or low-risk groups. (B) Survival status of patients in each group. (C) Heatmap of expression profiles of selected four OXPHOS-related prognostic differentially expressed genes (DEGs). (D) Kaplan–Meier survival curves analyses of risk score in TCGA cohort. In the train set, in the low-risk group, N = 128; in the high-risk group, N = 128. In the entire set, in the low-risk group, N = 244; in the high-risk group, N = 267. (E) Survival ROC analysis of risk score in predicting prognoses in the train and entire sets. (F) PCA plots of the train and entire sets. In the train set, N = 256; in the entire set, N = 511.
3.1.4 Clinical Utility of the Oxidative Phosphorylation-Related Prognostic Signature
The conjoint univariate and multivariate analyses in both sets confirmed the independence of the OXPHOS-related prognostic signature (Figure 3). Integrated with clinical factors, we performed subgroup survival analysis in stages I and II, stages III and IV, grades 1 and 2, grades 3 and 4, endometrial, mixed, and serous, age ≤6 0, age >60, and we found that a higher risk score was related to poorer prognosis in all subgroups (Supplementary Figures S4A–H). The correlation between risk core and clinical characteristics is shown in Supplementary Figures S4I–L and a heatmap (Supplementary Figure S5), indicating that patients with lower stage, lower grade, endometrial histological type, and lower age tended to have a lower risk score. Aimed at creating a feasible method for clinical application, we constructed a nomogram cooperating risk core, patient age, histological type, tumor grade, and tumor stage (Figure 4A). By comparing the AUC values at 1, 3, and 5 years plotted by the signature and clinical factors, we discovered that the gene-based signature was superior to others in the OS prediction. What is more, when the gene-based signature was combined with the clinical characteristics, it had better predictive ability than the signature alone (Figure 4B). The calibration curves also demonstrated that the results predicted by the nomogram were in high consistency with actual survival, which suggested that the nomogram had a high clinical value (Figure 4C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g003.jpg)
FIGURE 3. The regression analyses of the independent predictor of uterine corpus endometrial carcinoma (UCEC). In the train set, N = 256; in the entire set, N = 511.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g004.jpg)
FIGURE 4. Construction of nomogram for predicting survival of patients with UCEC (N = 511). (A) The baseline nomogram including five factors: patient age, histological type, tumor grade, tumor stage, and risk score. (B) ROC curves of 1-, 3-, and 5-years OS indicated that the gene-based signature had better predictive ability compared with other clinical characteristics. In addition, the nomogram composed of the signature and clinical characteristics showed superior prediction value for the prognosis. (C) Calibration plots of 1-, 3-, and 5-years OS suggested that the prognosis predicted by the nomogram was consistent with the actual outcome.
3.1.5 Gene Set Enrichment Analyses
To ascertain involved biological processes, we employed GSEA to analyze the transcript message of patients in different risk groups. In the high-risk group, representative KEGG pathways were axon guidance, glycosaminoglycan biosynthesis ghondroitin sulfate, glycosaminoglycan biosynthesis heparan sulfate, RNA polymerase, and tight junction (Figure 5A), while in the low-risk group, representative KEGG pathways were cytokine–cytokine receptor interaction, chemokine signaling pathway, and so on (Figure 5B). Through comparison, we discovered that in the low-risk group, the representative biological processes were mainly related to immunization. Furthermore, we implemented the ESTIMATE algorithm and found that the higher the risk score, the lower the immune score, the stromal score, and the ESTIMATE score, and the higher the tumor purity (Figures 5C–F), which further proved that the signature was significantly related to immunization.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g005.jpg)
FIGURE 5. The gene set enrichment analysis, the immune, stromal, and ESTIMATE scores, and the tumor purity of different risk groups. (A,B) The gene set enrichment analysis showed the significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in TCGA. (C–F) The AVONA results of the immune, stromal, and ESTIMATE scores and tumor purity in the high- and low-risk groups. Patients in the high-risk group had lower immune score, stromal score, ESTIMATE score, and higher tumor purity. In the high-risk group, N = 267; in the low-risk group, N = 244.
3.1.6 Difference of Tumor-Infiltrating Immune Cells in Different Risk Groups
To further explore the relationship between the OXPHOS-related prognostic signature and immunization, we found the difference of TIICs between different risk groups. M0 macrophages, M2 macrophages, and activated mast cells were significantly elevated in the high-risk group (p < 0.05). In contrast, plasma cells, CD8+ T cells, activated memory CD4+ T cells, and regulatory T cells (Tregs) were evidently elevated in the low-risk group (p < 0.05) (Figure 6A). We analyzed the correlation of risk score and TIICs; as we can see in Figures 6B,C, M0 macrophages, M2 macrophages, plasma cells, activated memory CD4+ T cells, CD8+ T cells, gamma delta T cells, and regulatory T cells (Tregs) were the seven most relevant types of immune cells with risk core. There was a significant difference in the enrichment scores of diverse immune cell subpopulations and related functions or pathways with ssGSEA between different risk groups, such as TIL and HLA, which showed coherence with the GO and KEGG analyses results (Figures 6D,E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g006.jpg)
FIGURE 6. The relationship between risk score and the tumor-infiltrating immune cells (TIICs). (A) Comparison of the TIICs between different risk groups. (B) The correlation of risk score and TIICs. (C) The associations between M0 macrophages, M2 macrophages, plasma cells, activated memory CD4+ T cells, CD8+ T cells, gamma delta T cells, regulatory T cells (Tregs), and risk score. Comparison of the single-sample gene-set enrichment analysis (ssGSEA) scores between different risk groups in TCGA cohort, including the scores of 16 immune cells (D) and 13 immune-related functions (E) displayed in boxplots. CCR, cytokine–cytokine receptor. Adjusted p-values were shown as ns, not significant; *p < 0.05; **p < 0.01; ***p < 0.001. In the high-risk group, N = 267; in the low-risk group, N = 244.
3.1.7 Comparison of Immune Status Between Different Risk Groups
IPS and the expression of immune modulators are able to predict the potential response of the patients to the immune checkpoint inhibitor. The scores of IPS (p = 0.004), IPS-PD1-PD-L1-PD-L2 (p = 1.297E−8), IPS-CTLA4 (7.334E−7), and IPS-PD1-PD-L1-PD-L2-CTLA4 (p = 3.776E−11) in the high-risk group were clearly higher than that in the low-risk group (Figures 7A–D). The expression of immune modulators including CD27, CTLA4, ICOS, PD-L2, TIGIT, PD-1, LAG3, TIM-3, CD86, PD-L1, CD70, and CD270 (p < 0.001) in the low-risk group was also remarkably higher than in the high-risk group (Figure 7E). We further discovered that the genes related to immune checkpoint inhibitor (ICI), including TIGIT, ICOS, CTLA4, CD27, CD58, PD-1, CD86, TIM-3, IDO1, and PD-L1 were negatively related to risk score (Figures 7F–Q), which suggested that the patients in the low-risk group promised to respond to immunotherapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g007.jpg)
FIGURE 7. Immunophenoscore (IPS) and gene expression analyses of immune checkpoint. (A–D) Patients in the low-risk group had higher scores of IPS, IPS-PD1-PD-L1-PD-L2, IPS-CTLA4, and IPS-PD1-PD-L1-PD-L2-CTLA4 than patients in the high-risk group. (E) Comparison of the expression of immune checkpoint between the different risk groups. (F–Q) The associations among TIGIT, ICOS, CTLA4, CD27, CD270, CD58, PD-1, CD86, TIM-3, IDO1, PD-L1, and risk score. In the high-risk group, N = 267; in the low-risk group, N = 244.
By comparing the reaction to immune checkpoint blockade (ICB), we found that the percentage of patients with response to ICB in the high-risk group (53.18%, 142/267) was lower than that in the low-risk group (75.41%, 184/244) (Figure 8A). The risk scores of the responders were also lower than those of the nonresponders (Figure 8B). These results could demonstrate that the OXPHOS-related prognostic signature had the potential to predict the possible effect of ICB in UCEC patients and work as a classifier for therapy selection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g008.jpg)
FIGURE 8. The response to immune checkpoint blockade (ICB) therapy and the microsatellite instability (MSI) status of different risk groups. (A) Comparison of the response to immune checkpoint blockade (ICB) therapy between the low- and the high-risk groups. Patients in the low-risk group showed more sensitive response to ICB therapy than those in the high-risk group. (B–G) Patients in the low-risk group had lower expression of the dominant mismatch repair (MMR) proteins including MLH1, MSH2, MSH6, PMS2, and EPCAM than patients in the high-risk group. (H) The mRNAsi of the low-risk group is lower than that of the high-risk group. (I) The correlation between the MSI expression and the risk score. (J) The expression of the markers of m6A mRNA methylation regulators in the high- and low-risk groups. m6A, N6-methyladenosine. In the high-risk group, N = 267; in the low-risk group, N = 244.
Moreover, we compared the MSI status of different risk groups and found, as shown in Figures 8C–G, the expression level of dominant mismatch repair (MMR) proteins, such as MLH1, MSH2, MSH6, PMS2, and EPCAM, in the low-risk group were lower, and risk score was negatively correlated with the MSI expression (Figure 8I). The patients in the low-risk group had a lower score of mRNAsi (Figure 8H), which can partly explain the poor prognosis in the high-risk group. In addition, significantly decreased gene expression of KIAA1429, FMR1, HINRNPC, ZC3H13, YTHDF1, and YTHDF3, which are key factors of m6A mRNA methylation could be found in Figure 8J. These findings all provided evidence for the efficiency of immunotherapy for patients in the low-risk group and helped to elucidate the mechanisms of different outcomes.
3.1.8 Relationship Between Tumor Mutational Burden and Risk Score
To further clarify the mechanism of tumorigenesis and development of UCEC patients, we compared the tumor mutational burden (TMB) between different risk groups based on somatic mutation data. We discovered that the low-risk group has remarkably higher TMB than the high-risk group (p = 1.278E−8) (Figure 9A). The correlation between risk score and TMB was negative (Cor = −0.269, p = 6.178E−11) (Figure 9B). We also found that only the mutation of TP53 was found in the high-risk group, whereas the others were found in the low-risk group (Figure 9C). Next, we divided UCEC patients into H-TMB and L-TMB groups based on the median TMB; as shown in Figure 9D, the OS of the H-TMB group was better than that of the L-TMB group (log-rank test, p < 0.001). Stratified survival analysis was conducted to assess the synergistic effect of these factors, which corroborated that TMB would not influence the effect of the signature in prognosis prediction. Significant survival differences were observed in both TMB-stratified subgroups [log-rank test, H-TMB and high risk score (HH) versus H-TMB and low risk score (HL), p < 0.001; L-TMB and high risk score (LH) versus L-TMB and low risk score (LL), p < 0.001; Figure 9E], which indicated that the signature could function as an independent predictive indicator and effectively evaluate the potential response to immunotherapy. The most frequent somatic mutations in the high-risk group were as follows: PTEN > TP53 > PIK3CA > ARID1A > TTN > PIK3R1 > CTNNB1 > KMT2D > CHD4 > CSMD3 (Figure 9F). While the most frequent somatic mutations in the low-risk group followed the order PTEN > ARID1A > PIK3CA > TTN > PIK3R1 > CTCF > MUC16 > ZFHX3 > KMT2D > MUC5B (Figure 9G). Other mutation details, including variant classification, mutation type, SNV class, mutation load, and top 10 mutated genes in UCEC between different risk groups, are shown in Supplementary Figure S6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g009.jpg)
FIGURE 9. The relationship between risk score and tumor mutational burden (TMB). (A) Comparison of TMB between the low- and high-risk groups. (B) The correlation between risk score and TMB was negative. (C) The regression analyses of the mutated genes. (D) Kaplan–Meier curves for the high and low TMB groups of the TCGA–UCEC cohort. Log-rank test, p < 0.001. (E) Kaplan–Meier curves for patients in the TCGA–UCEC cohort stratified by both TMB and risk score. Log-rank test, p < 0.001. (F) Oncoplot displaying the somatic landscape of UCEC in the high-risk group. (G) Oncoplot displaying the somatic landscape of UCEC in the low-risk group. In the high-risk group, N = 267; in the low-risk group, N = 244.
3.1.9 The Chemotherapeutic Response Between Different Risk Groups
The high-risk group had a higher estimated IC50 of temsirolimus, roscovitine, AZD6244, PD.0325901, RDEA119, PF.02341066, AKT. inhibitor.VIII, BMS.509744, vinblastine, bryostatin.1, metformin, AZ628, nutlin.3a, bortezomib, bicalutamide, AZD6482, cytarabine, rapamycin, camptothecin, mitomycin. C, S. Trityl.L.cysteine, tipifarnib, parthenolide, sorafenib, and methotrexate than the low-risk group (p < 0.0001), which demonstrated that patients with a low risk score were more sensitive to these drugs (Figure 10). Therefore, the OXPHOS-related prognostic signature could also be utilized as a robust and valuable classifier for chemotherapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/753004/fcell-09-753004-HTML/image_m/fcell-09-753004-g010.jpg)
FIGURE 10. Prediction of the response to common chemotherapy drugs between the low- and high-risk groups. Patients in the high-risk group possessed higher estimated IC50 than patients in the low-risk group. In the high-risk group, N = 267; in the low-risk group, N = 244.
4 Discussion
Though the overall survival state of UCEC patients has improved a lot due to the development of early diagnostic techniques (Braun et al., 2016), there are still 15% of patients with UCEC diagnosed at an advanced stage, thus, suffering from tumor recurrence owing to a limited response to unsuitable therapies (Li and Wan, 2020). Therefore, developing valid and reliable biomarkers for survival prediction and treatment selection to avoid overtreatment in patients who will not relapse and propose individual adjuvant treatment to patients who will recur counts (Lee et al., 2017), and efforts have been put into developing a robust prognostic signature based on DEGs for identification of biologic subsets to guide treatment strategies for UCEC patients, such as a cell cycle-related prognostic signature (Liu et al., 2020a), an autophagy-related prognostic signature (Wang et al., 2020b), and so on. However, there is no consensus on predictive biomarkers for scientifically validated therapy.
In recent years, emerging evidence has proved that OXPHOS could function as a new target in cancer therapy (Ashton et al., 2018). However, OXPHOS has not been involved in developing a prognostic signature for UCEC, thus far, and the potential underlying mechanism and possible clinical application of OXPHOS-related signature in cancer are worth exploring.
Based on the 289 OXPHOS-related gene matrix obtained from the research of Zihao Xu, we used differential expression analysis and identified 66 OXPHOS-related DEGs, including 54 upregulated genes and 12 downregulated genes; most of them have been reported to be associated with cancer. For instance, EPAS1 could enhance the effect of paclitaxel on breast cancer cells by inhibiting growth and promoting apoptosis of MCF-7/TAX cells (Song et al., 2020). Furthermore, functional annotation of the OXPHOS-related DEGs also indicated that the candidate genes were significantly related to ATP metabolic process, purine ribonucleoside triphosphate metabolic process, and the essential pathways reported in many cancers (Matsuda et al., 2017; Fallah and Rini, 2019), including retrograde endocannabinoid signaling and HIF-1 signaling pathways.
Then, we randomly selected samples from the entire set in half into the train set and constructed the model based on it. After univariate Cox analysis and LASSO Cox regression analysis, finally, a prognostic model composed of four OXPHOS-related genes (ATP5IF1, COX6B1, FOXP3, and NDUFB11) was built. Except for ATP5IF1, other candidate genes included in the prognostic signature have been proven to play a vital role in cancer development. COX6B1, as one of the core proteins involved in the molecular basis of the below-background radiation stress response, might inhibit the proliferation of laryngeal squamous cell carcinoma cells (Liu et al., 2020b). FOXP3 can activate the Wnt/beta-catenin signaling pathway, thus, inducing epithelial–mesenchymal transition and promoting the progression of non-small cell lung cancer (Yang et al., 2017). In addition, the regulation of the expression of the NDUFB11 gene can regulate the programmed cell death process, so it promises to be a new target for cancer therapy (Panelli et al., 2013).
Furthermore, Kaplan–Meier survival curves in the train and entire set indicated that the patients with high risk score tend to have a worse prognosis. Besides, a survival ROC curve analysis demonstrated the predictive value of the OXPHOS-related prognostic signature. In addition, by conjoint univariate and multivariate analyses, the independence of the OXPHOS-related signature for prognosis prediction was verified. By combining the clinical features, we performed subgroup survival analysis in different risk groups and discovered that high risk was related to poor prognosis in all subgroups. Primarily, we observed that the risk score was highly correlated to tumor stage, tumor grade, histological type, and patient age. What is more, multi-ROC curve analyses and the evaluation of the nomogram provided evidence for the potential clinical application of the signature in prognosis prediction.
The results of GSEA indicated that the representative pathways in the low-risk group were mainly related to immunization. What is more, the immune score, stromal score, and ESTIMATE score of the low-risk group were higher, and the tumor purity was lower, which implied that patients with a low risk score appeared to be more immunogenic and might preferentially benefit from immunotherapy. To penetrate the immune mechanism of the signature, we compared TIICs between the two different risk groups and analyzed the relationship between the infiltrating immune cells and risk score. We found that plasma cells, CD8+ T cells, activated memory CD4+ T cells, and regulatory T cells (Tregs) were remarkably elevated in the low-risk group, while M0 macrophages, M2 macrophages, and activated mast cells were significantly elevated in the high-risk group. Immune cell profiles in the tumor microenvironment (TME) of UCEC tissues could be utilized to predict the survival of patients with UCEC (Li and Wan, 2020). Under the influence of environmental stimuli, tumor-associated macrophages (TAMs) can differentiate into two types of macrophage with opposite functions, in which M1 macrophages are protective immune cells that can remove pathogens and malignant cells. In contrast, M2 macrophages are dangerous factors, which can promote angiogenesis and the growth of tumor. High M2 macrophage infiltration in most tumor types is correlated with an adverse prognosis (Zhao et al., 2017; Wang et al., 2018), which is consistent with our results. A meta-analysis proved that UCEC patients with a high CD8+ T-cell density had a favorable prognosis (Guo et al., 2020). Increased expression of activated CD4+ memory T cells had been confirmed to be beneficial for prognosis in bladder cancer (Li et al., 2020), and the previous study of our team also provided evidence for the opinion (Liu et al., 2020c). The regulatory T cell (Treg) could suppress the antitumor immune response, thus, becoming a target of immunotherapy (Tanaka and Sakaguchi, 2017; Antomarchi et al., 2019), which is different from our results. However, another study showed that FOXP3 lymphocytic infiltration had no significant impact on the survival state of UCEC patients (Giatromanolaki et al., 2008). Our results suggested that CD8+ T cells, activated memory CD4+ T cells, regulatory T cells (Tregs), M0 macrophages, and M2 macrophages were crucial to the development of UCEC, which had the potential to be targeted for immunotherapies. The immune cell subpopulation enrichment analyses, as well as the ssGSEA-related functions or pathway analyses, showed that the scores of all the items except for aDCs and type I IFN response were higher in the low-risk group, which further proved that UCEC patients with low risk scores appeared to be more immunogenic.
Furthermore, patients with low risk scores had significantly higher scores of IPS, IPS-PD1-PD-L1-PD-L2, IPS-CTLA4, and IPS-PD1-PD-L1-PD-L2-CTLA4; thereby, they tend to respond more sensitively to immune checkpoint inhibitors. The Food and Drug Administration (FDA) has approved a series of monoclonal antibodies (mAbs) targeting PD-1, PD-L1, and CTLA-4 in several malignancy therapies (Arora et al., 2018). Besides, several clinical trials demonstrated the clinical value of anti-PD-1 antibody, pembrolizumab, which can inhibit tumor immune escaping in endometrial cancer. It has been approved as an alternative therapy for patients in the advanced stage (Musacchio et al., 2020). Also, the expression of CD27, CTLA4, ICOS, PD-L2, TIGIT, PD-1, LAG3, TIM-3, CD86, PD-L1, CD70, and CD270 were also remarkably higher in the low-risk group, which deserved to be explored as a new target for immunotherapy. TIGIT could enhance T-cell function by inhibiting PVRIG, so TIGIT-PVRIG pathways play a vital role in human cancers (Whelan et al., 2019). LAG-3 has been identified on various tumor cells, including UCEC, as a critical inducer for the malignancy progression (Friedman et al., 2020). There is also increasing evidence for TIM-3 as an emerging target for malignancy immunotherapy (Moore et al., 2019). The difference found in the comparison of the response to immune checkpoint blockade (ICB) therapy between different risk groups showed the same trend with the above findings. We also compared the risk score of the patients showing a different response to ICB, and found that the responders tended to possess a low risk score. These results demonstrated that the OXPHOS-related prognostic signature had the potential to be the classifier for a specific subset of UCEC patients, who would have a favorable effect on ICB therapy. Compared with UCEC patients in the advanced stage, patients with primary tumor tend to have more somatic loss of mismatch repair (MMR) protein expression; thus, the loss of MMR should be considered as guidance for immunotherapy (Ta et al., 2018). We compared the MSI status of different risk groups and found that the expression of dominant MMR proteins (MLH1, MSH2, MSH6, PMS2, and EPCAM) and mRNAsi in the low-risk group was downregulated, which illustrated that the patients with low risk score were burdened with high-MSI and showed fewer stemness features associated with oncogenic dedifferentiation. Haoya Xu et al. found that cell stemness can predict patient prognosis and constructed an mRNAsi-related prognostic signature in endometrial cancer (Xu et al., 2021). Accumulating evidence revealed that m6A mRNA methylation participates in the pathogenesis and progression of cancers through molecular mechanisms, such as inhibiting the antitumor response of CD8+ T cells (Chen et al., 2019; Sun et al., 2019). Patients with a low risk score had decreased gene expression of the m6A mRNA methylation, thus, showing a favorable response to immunotherapy.
To explore the cause of tumorigenesis, we compared the TMB of each group, and we discovered that the low-risk group had remarkably higher TMB, and the risk score was negatively correlated with TMB. It is reported that patients with high TMB had a more sensitive response to immune checkpoint blockade, such as anti-PD-1 agents (Fancello et al., 2019). A multicohort study demonstrated the value of tissue TMB (tTMB) as a predictive biomarker for response to pembrolizumab monotherapy (Marabelle et al., 2020). In our present study, the most frequent somatic mutations in the low-risk group were PTEN, ARID1A, PIK3CA, TTN, PIK3R1, CTCF, MUC16, ZFHX3, KMT2D, and MUC5B. PTEN functions as a tumor inhibitor by decreasing the activity of phosphoinositide 3-kinase (PI3K) (Bian et al., 2018), and its expression diminished with the progression of UCEC (Xi et al., 2019). It has been proven that ARID1A alterations indicated improved prognosis of tumor patients for compromising MMR and enhancing the infiltration of lymphocytes (Jiang et al., 2020). UCEC patients with high PIK3CA had shorter survival time because PIK3CA could change the tumor immune microenvironment by altering the fraction of tumor-associated neutrophils (Pan et al., 2019). Eighty-five percent of lung squamous cell carcinoma patients with TTN mutant type had missense variations and favorable overall survival (Cheng et al., 2019). In addition, patients in the high TMB group showed better overall survival than in the low TMB group, and the combined Kaplan–Meier analyses proved that risk score was a reliable predictive indicator independent of TMB. Our results were consistent with previous studies, indicating that UCEC patients with low risk score tend to have improved response to immunotherapy and favorable prognosis.
We compared the predicted effectiveness of chemotherapy drugs between different risk groups by using the GDSC dataset. Patients with UCEC in the low-risk group were more sensitive to temsirolimus, roscovitine, AZD6244, PD.0325901, RDEA119, PF.02341066, AKT. inhibitor.VIII, BMS.509744, vinblastine, bryostatin.1, metformin, AZ628, nutlin.3a, bortezomib, bicalutamide, AZD6482, cytarabine, rapamycin, camptothecin, mitomycin. C, S. Trityl.L.cysteine, tipifarnib, parthenolide, sorafenib, and methotrexate, which might be a potential helpful treatment option for UCEC.
These findings indicated that patients with low risk scores possessed higher IPS, higher TMB, and responded more sensitively to ICI and chemotherapy, which partly explained the reason for the improved prognosis of the low-risk patients.
However, there are still some deficiencies with our study. First of all, our data were based entirely on TCGA databases and lacked an external database for validation. Second, our molecular results lack the evidence of the experiments in vitro or in vivo. Third, these OXPHOS-related prognostic genes are supposed to be further examined by cell function assay for biological function and pathway analyses.
5 Conclusion
In conclusion, an OXPHOS-related and immune microenvironment prognostic signature classifying EC patients into different risk subsets was constructed in our study. This model could predict the prognosis of patients and help to select a specific subset of patients who might benefit from immunotherapy and chemotherapy, thus, improving the overall survival rate. These results may also contribute to the identification of novel immune biomarkers or target therapy in UCEC and give a new insight into the tumorigenesis and progression of UCEC.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the Nanjing Medical University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
JB, HY, JL, TC, and ZZ conceived the study and participated in the study design, performance, and manuscript writing. MY, SN, and SY conducted the bioinformatics analysis. LC and FS revised the manuscript. All authors read and approved the final manuscript.
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.
Acknowledgments
We would like to thank the researchers and study participants for their contributions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.753004/full#supplementary-material
Supplementary Figure S1 | Selected OXPHOS-related DEGs from uterine corpus endometrial carcinoma samples. Heatmap (A) and volcano plot (B) of OXPHOS-related DEGs between UCEC and normal endometrium. Red dots represented up-regulated OXPHOS-related DEGs, green dots represented down-regulated OXPHOS-related DEGs and black dots represented no OXPHOS-related DEGs. Functional annotation of the OXPHOS-related DEGs included Gene ontology analysis (C) and Kyoto Encyclopedia of Genes and Genomes pathway analysis (D). (E) The PPI network indicated the interactions between OXPHOS-related DEGs.
Supplementary Figure S2 | Construction of a prognostic model based on OXPHOS-related DEGs. (A,B) LASSO Cox regression analysis and multivariate Cox regression analysis. (C) Forest plots present the multivariate Cox proportional hazards regression analysis of the prognosis-associated OXPHOS-related DEGs.
Supplementary Figure S3 | (A–D) The expression of the four prognostic OXPHOS-related DEGs in tumor (N = 511) and normal tissue (N = 23). (E,F) The relationship between ATP5IF1, FOXP3 and the overall survival of UCEC patients. (G–J) Expression levels of the four OXPHOS-related DEGs mRNA in clinical tissue samples.
Supplementary Figure S4 | Stratified analysis. (A–H) The Kaplan-Meier survival analyses based on the prognostic signature stratified by tumor stage, tumor grade, histological type, patient age in the TCGA-UCEC cohort. (I–L) Relationship between the prognostic signature and tumor stage, tumor grade, histological type and age in the TCGA-UCEC cohort.
Supplementary Figure S5 | The heatmap shows the OXPHOS-related prognostic genes expression and clinical characteristics of the TCGA-UCEC cohort.
Supplementary Figure S6 | Variant classification, mutation type, SNV class, mutation load and top 10 mutated genes in UCEC between low (A) and high risk groups (B). In the high risk group, N = 267; in the low risk group, N = 244.
References
Antomarchi, J., Ambrosetti, D., Cohen, C., Delotte, J., Chevallier, A., Karimdjee-Soilihi, B., et al. (2019). Immunosuppressive Tumor Microenvironment Status and Histological Grading of Endometrial Carcinoma. Cancer Microenvironment 12, 169–179. doi:10.1007/s12307-019-00225-1
Arora, E., Masab, M., Mittar, P., Jindal, V., Gupta, S., and Dourado, C. (2018). Role of Immune Checkpoint Inhibitors in Advanced or Recurrent Endometrial Cancer. Cureus 10, e2521. doi:10.7759/cureus.2521
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: Tool for the Unification of Biology. Nat. Genet. 25, 25–29. doi:10.1038/75556
Ashton, T. M., McKenna, W. G., Kunz-Schughart, L. A., and Higgins, G. S. (2018). Oxidative Phosphorylation as an Emerging Target in Cancer Therapy. Clin. Cancer Res. 24, 2482–2490. doi:10.1158/1078-0432.ccr-17-3070
Bian, X., Gao, J., Luo, F., Rui, C., Zheng, T., Wang, D., et al. (2018). PTEN Deficiency Sensitizes Endometrioid Endometrial Cancer to Compound PARP-Pi3k Inhibition but Not PARP Inhibition as Monotherapy. Oncogene 37, 341–351. doi:10.1038/onc.2017.326
Braun, M. M., Overbeek-Wager, E. A., and Grumbo, R. J. (2016). Diagnosis and Management of Endometrial Cancer. Am. Fam. Physician 93, 468–474.
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: a Cancer J. clinicians 68, 394–424. doi:10.3322/caac.21492
Cancer Genome Atlas Research, N., Kandoth, C., Schultz, N., Cherniack, A. D., Akbani, R., Liu, Y., et al. (2013). Integrated Genomic Characterization of Endometrial Carcinoma. Nature 497, 67–73. doi:10.1038/nature12113
Cardenas, C., Lovy, A., Silva-Pavez, E., Urra, F., Mizzoni, C., Ahumada-Castro, U., et al. (2020). Cancer Cells with Defective Oxidative Phosphorylation Require Endoplasmic Reticulum-To-Mitochondria Ca2+ Transfer for Survival. Sci. Signal. 13, eaay1212. doi:10.1126/scisignal.aay1212
Charoentong, P., Finotello, F., Angelova, M., Mayer, C., Efremova, M., Rieder, D., et al. (2017). Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cel Rep. 18, 248–262. doi:10.1016/j.celrep.2016.12.019
Chen, X.-Y., Zhang, J., and Zhu, J.-S. (2019). The Role of m6A RNA Methylation in Human Cancer. Mol. Cancer 18, 103. doi:10.1186/s12943-019-1033-z
Cheng, X., Yin, H., Fu, J., Chen, C., An, J., Guan, J., et al. (2019). Aggregate Analysis Based on TCGA: TTN Missense Mutation Correlates with Favorable Prognosis in Lung Squamous Cell Carcinoma. J. Cancer Res. Clin. Oncol. 145, 1027–1035. doi:10.1007/s00432-019-02861-y
Commander, R., Wei, C., Sharma, A., Mouw, J. K., Burton, L. J., Summerbell, E., et al. (2020). Subpopulation Targeting of Pyruvate Dehydrogenase and GLUT1 Decouples Metabolic Heterogeneity during Collective Cancer Cell Invasion. Nat. Commun. 11, 1533. doi:10.1038/s41467-020-15219-7
Fallah, J., and Rini, B. I. (2019). HIF Inhibitors: Status of Current Clinical Development. Curr. Oncol. Rep. 21, 6. doi:10.1007/s11912-019-0752-z
Fancello, L., Gandini, S., Pelicci, P. G., and Mazzarella, L. (2019). Tumor Mutational burden Quantification from Targeted Gene Panels: Major Advancements and Challenges. J. Immunotherapy Cancer 7, 183. doi:10.1186/s40425-019-0647-4
Friedman, L. A., Ring, K. L., and Mills, A. M. (2020). LAG-3 and GAL-3 in Endometrial Carcinoma. Int. J. Gynecol. Pathol. 39, 203–212. doi:10.1097/pgp.0000000000000608
Geeleher, P., Cox, N. J., and Huang, R. (2014). Clinical Drug Response Can Be Predicted Using Baseline Gene Expression Levels and In Vitro Drug Sensitivity in Cell Lines. Genome Biol. 15, R47. doi:10.1186/gb-2014-15-3-r47
Giatromanolaki, A., Bates, G. J., Koukourakis, M. I., Sivridis, E., Gatter, K. C., Harris, A. L., et al. (2008). The Presence of Tumor-Infiltrating FOXP3+ Lymphocytes Correlates with Intratumoral Angiogenesis in Endometrial Cancer. Gynecol. Oncol. 110, 216–221. doi:10.1016/j.ygyno.2008.04.021
Guo, F., Dong, Y., Tan, Q., Kong, J., and Yu, B. (2020). Tissue Infiltrating Immune Cells as Prognostic Biomarkers in Endometrial Cancer: A Meta-Analysis. Dis. Markers 2020, 1805764. doi:10.1155/2020/1805764
He, Y., Jiang, Z., Chen, C., and Wang, X. (2018). Classification of Triple-Negative Breast Cancers Based on Immunogenomic Profiling. J. Exp. Clin. Cancer Res. 37, 327. doi:10.1186/s13046-018-1002-1
Heagerty, P. J., Lumley, T., and Pepe, M. S. (2000). Time-dependent ROC Curves for Censored Survival Data and a Diagnostic Marker. Biometrics 56, 337–344. doi:10.1111/j.0006-341x.2000.00337.x
Iasonos, A., Schrag, D., Raj, G. V., and Panageas, K. S. (2008). How to Build and Interpret a Nomogram for Cancer Prognosis. J. Clin. Oncol. 26, 1364–1370. doi:10.1200/jco.2007.12.9791
Jiang, T., Chen, X., Su, C., Ren, S., and Zhou, C. (2020). Pan-cancer Analysis of ARID1A Alterations as Biomarkers for Immunotherapy Outcomes. J. Cancer 11, 776–780. doi:10.7150/jca.41296
Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Lai, Y. J., Chou, Y. C., Lin, Y. J., Yu, M. H., Ou, Y. C., Chu, P. W., et al. (2019). Pyruvate Kinase M2 Expression: A Potential Metabolic Biomarker to Differentiate Endometrial Precancer and Cancer that Is Associated with Poor Outcomes in Endometrial Carcinoma. Int. J. Environ. Res. Public Health 16, 4589. doi:10.3390/ijerph16234589
LeBleu, V. S., O’Connell, J. T., Gonzalez Herrera, K. N., Wikman, H., Pantel, K., Haigis, M. C., et al. (2014). PGC-1α Mediates Mitochondrial Biogenesis and Oxidative Phosphorylation in Cancer Cells to Promote Metastasis. Nat. Cel Biol 16, 9921–100315. doi:10.1038/ncb3039
Lee, Y. C., Lheureux, S., and Oza, A. M. (2017). Treatment Strategies for Endometrial Cancer: Current Practice and Perspective. Curr. Opin. Obstet. Gynecol. 29, 47–58. doi:10.1097/gco.0000000000000338
Li, B. L., and Wan, X. P. (2020). Prognostic Significance of Immune Landscape in Tumour Microenvironment of Endometrial Cancer. J. Cel Mol Med 24, 7767–7777. doi:10.1111/jcmm.15408
Li, W., Zeng, J., Luo, B., Mao, Y., Liang, Y., Zhao, W., et al. (2020). High Expression of Activated CD4+ Memory T Cells and CD8+ T Cells and Low Expression of M0 Macrophage Are Associated with Better Clinical Prognosis in Bladder Cancer Patients. Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi 36, 97–103.
Liu, J., Mei, J., Li, S., Wu, Z., and Zhang, Y. (2020a). Establishment of a Novel Cell Cycle-Related Prognostic Signature Predicting Prognosis in Patients with Endometrial Cancer. Cancer Cel Int 20, 329. doi:10.1186/s12935-020-01428-z
Liu, J., Ma, T., Gao, M., Liu, Y., Liu, J., Wang, S., et al. (2020b). Proteomic Characterization of Proliferation Inhibition of Well-Differentiated Laryngeal Squamous Cell Carcinoma Cells under Below-Background Radiation in a Deep Underground Environment. Front. Public Health 8, 584964. doi:10.3389/fpubh.2020.584964
Liu, J., Chen, X., Jiang, Y., and Cheng, W. (2020c). Development of an Immune Gene Prognostic Classifier for Survival Prediction and Respond to Immunocheckpoint Inhibitor Therapy/chemotherapy in Endometrial Cancer. Int. immunopharmacology 86, 106735. doi:10.1016/j.intimp.2020.106735
Malta, T. M., Sokolov, A., Gentles, A. J., Burzykowski, T., Poisson, L., Weinstein, J. N., et al. (2018). Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation. Cell 173, 338–354. doi:10.1016/j.cell.2018.03.034
Marabelle, A., Fakih, M., Lopez, J., Shah, M., Shapira-Frommer, R., Nakagawa, K., et al. (2020). Association of Tumour Mutational burden with Outcomes in Patients with Advanced Solid Tumours Treated with Pembrolizumab: Prospective Biomarker Analysis of the Multicohort, Open-Label, Phase 2 KEYNOTE-158 Study. Lancet Oncol. 21, 1353–1365. doi:10.1016/s1470-2045(20)30445-9
Matsuda, F., Toya, Y., and Shimizu, H. (2017). Learning from Quantitative Data to Understand central Carbon Metabolism. Biotechnol. Adv. 35, 971–980. doi:10.1016/j.biotechadv.2017.09.006
Molina, J. R., Sun, Y., Protopopova, M., Gera, S., Bandi, M., Bristow, C., et al. (2018). An Inhibitor of Oxidative Phosphorylation Exploits Cancer Vulnerability. Nat. Med. 24, 1036–1046. doi:10.1038/s41591-018-0052-4
Moore, M., Ring, K. L., and Mills, A. M. (2019). TIM-3 in Endometrial Carcinomas: an Immunotherapeutic Target Expressed by Mismatch Repair-Deficient and Intact Cancers. Mod. Pathol. 32, 1168–1179. doi:10.1038/s41379-019-0251-7
Mosier, J. A., Schwager, S. C., Boyajian, D. A., and Reinhart-King, C. A. (2021). Cancer Cell Metabolic Plasticity in Migration and Metastasis. Clin. Exp. metastasis 38 (4), 343–359. doi:10.1007/s10585-021-10102-1
Musacchio, L., Boccia, S. M., Caruso, G., Santangelo, G., Fischetti, M., Tomao, F., et al. (2020). Immune Checkpoint Inhibitors: A Promising Choice for Endometrial Cancer Patients? J. Clin. Med. 9, 1721. doi:10.3390/jcm9061721
Nayak, A. P., Kapur, A., Barroilhet, L., and Patankar, M. S. (2018). Oxidative Phosphorylation: A Target for Novel Therapeutic Strategies against Ovarian Cancer. Cancers (Basel) 10, 337. doi:10.3390/cancers10090337
Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust Enumeration of Cell Subsets from Tissue Expression Profiles. Nat. Methods 12, 453–457. doi:10.1038/nmeth.3337
Pan, Y., Jia, L. P., Liu, Y., Han, Y., and Deng, Q. (2019). Alteration of Tumor Associated Neutrophils by PIK3CA Expression in Endometrial Carcinoma from TCGA Data. J. Ovarian Res. 12, 81. doi:10.1186/s13048-019-0557-6
Panelli, D., Lorusso, F. P., Papa, F., Panelli, P., Stella, A., Caputi, M., et al. (2013). The Mechanism of Alternative Splicing of the X-Linked NDUFB11 Gene of the Respiratory Chain Complex I, Impact of Rotenone Treatment in Neuroblastoma Cells. Biochim. Biophys. Acta (Bba) - Gene Regul. Mech. 1829, 211–218. doi:10.1016/j.bbagrm.2012.12.001
Qu, L., Wang, Z.-l., Chen, Q., Li, Y.-m., He, H.-w., Hsieh, J. J., et al. (2018). Prognostic Value of a Long Non-coding RNA Signature in Localized Clear Cell Renal Cell Carcinoma. Eur. Urol. 74, 756–763. doi:10.1016/j.eururo.2018.07.032
Robinson, D. R., Wu, Y.-M., Lonigro, R. J., Vats, P., Cobain, E., Everett, J., et al. (2017). Integrative Clinical Genomics of Metastatic Cancer. Nature 548, 297–303. doi:10.1038/nature23306
Sica, V., Bravo‐San Pedro, J. M., Stoll, G., and Kroemer, G. (2020). Oxidative Phosphorylation as a Potential Therapeutic Target for Cancer Therapy. Int. J. Cancer 146, 10–17. doi:10.1002/ijc.32616
Song, Y., Zhang, M., Lu, M. M., Qu, L. Y., Xu, S. G., Li, Y. Z., et al. (2020). EPAS1 Targeting by miR-152-3p in Paclitaxel-Resistant Breast Cancer. J. Cancer 11, 5822–5830. doi:10.7150/jca.46898
Sun, T., Wu, R., and Ming, L. (2019). The Role of m6A RNA Methylation in Cancer. Biomed. Pharmacother. 112, 108613. doi:10.1016/j.biopha.2019.108613
Ta, R. M., Hecht, J. L., and Lin, D. I. (2018). Discordant Loss of Mismatch Repair Proteins in Advanced Endometrial Endometrioid Carcinoma Compared to Paired Primary Uterine Tumors. Gynecol. Oncol. 151, 401–406. doi:10.1016/j.ygyno.2018.10.012
Tanaka, A., and Sakaguchi, S. (2017). Regulatory T Cells in Cancer Immunotherapy. Cell Res 27, 109–118. doi:10.1038/cr.2016.151
Toh, J. W. T., Phan, K., Reza, F., Chapuis, P., and Spring, K. J. (2021). Rate of Dissemination and Prognosis in Early and Advanced Stage Colorectal Cancer Based on Microsatellite Instability Status: Systematic Review and Meta-Analysis. Int. J. colorectal Dis. 36 (8), 1573–1596. doi:10.1007/s00384-021-03874-1
Tomczak, K., Czerwińska, P., and Wiznerowicz, M. (2015). The Cancer Genome Atlas (TCGA): an Immeasurable Source of Knowledge. Contemp. Oncol. (Pozn) 19, A68–A77. doi:10.5114/wo.2014.47136
Wang, N., Liu, W., Zheng, Y., Wang, S., Yang, B., Li, M., et al. (2018). CXCL1 Derived from Tumor-Associated Macrophages Promotes Breast Cancer Metastasis via Activating NF-Κb/sox4 Signaling. Cell Death Dis 9, 880. doi:10.1038/s41419-018-0876-3
Wang, T., Kong, S., Tao, M., and Ju, S. (2020a). The Potential Role of RNA N6-Methyladenosine in Cancer Progression. Mol. Cancer 19, 88. doi:10.1186/s12943-020-01204-7
Wang, H., Ma, X., Liu, J., Wan, Y., Jiang, Y., Xia, Y., et al. (2020b). Prognostic Value of an Autophagy-Related Gene Expression Signature for Endometrial Cancer Patients. Cancer Cel Int 20, 306. doi:10.1186/s12935-020-01413-6
Whelan, S., Ophir, E., Kotturi, M. F., Levy, O., Ganguly, S., Leung, L., et al. (2019). PVRIG and PVRL2 Are Induced in Cancer and Inhibit CD8+ T-Cell Function. Cancer Immunol. Res. 7, 257–268. doi:10.1158/2326-6066.cir-18-0442
Xi, Z., Jing, L., Le-Ni, K., Zhu, L., Ze-Wen, D., Hui, Y., et al. (2019). Evaluation of PTEN and CD4+FOXP3+ T Cell Expressions as Diagnostic and Predictive Factors in Endometrial Cancer. Medicine 98, e16345. doi:10.1097/md.0000000000016345
Xu, H., Zou, R., Liu, J., and Zhu, L. (2021). A Risk Signature with Nine Stemness Index-Associated Genes for Predicting Survival of Patients with Uterine Corpus Endometrial Carcinoma. J. Oncol. 2021, 6653247. doi:10.1155/2021/6653247
Xu, Y., Xue, D., Bankhead, A., and Neamati, N. (2020a). Why All the Fuss about Oxidative Phosphorylation (OXPHOS)? J. Med. Chem. 63, 14276–14307. doi:10.1021/acs.jmedchem.0c01013
Xu, Z., Wu, Z., Zhang, J., Zhou, R., Ye, L., Yang, P., et al. (2020b). Development and Validation of an Oxidative Phosphorylation-Related Gene Signature in Lung Adenocarcinoma. Epigenomics 12, 1333–1348. doi:10.2217/epi-2020-0217
Yang, S., Liu, Y., Li, M.-Y., Ng, C. S. H., Yang, S.-l., Wang, S., et al. (2017). FOXP3 Promotes Tumor Growth and Metastasis by Activating Wnt/β-Catenin Signaling Pathway and EMT in Non-small Cell Lung Cancer. Mol. Cancer 16, 124. doi:10.1186/s12943-017-0700-1
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, D955–D961. doi:10.1093/nar/gks1111
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612
Zhao, X., Qu, J., Sun, Y., Wang, J., Liu, X., Wang, F., et al. (2017). Prognostic Significance of Tumor-Associated Macrophages in Breast Cancer: a Meta-Analysis of the Literature. Oncotarget 8, 30576–30586. doi:10.18632/oncotarget.15736
Keywords: uterine corpus endometrial carcinoma, oxidative phosphorylation, prognosis, tumor microenvironment, immunotherapy
Citation: Liu J, Chen T, Yang M, Zhong Z, Ni S, Yang S, Shao F, Cai L, Bai J and Yu H (2021) Development of an Oxidative Phosphorylation-Related and Immune Microenvironment Prognostic Signature in Uterine Corpus Endometrial Carcinoma. Front. Cell Dev. Biol. 9:753004. doi: 10.3389/fcell.2021.753004
Received: 04 August 2021; Accepted: 25 October 2021;
Published: 25 November 2021.
Edited by:
Na Luo, Nankai University, ChinaReviewed by:
Chun Wang, Thomas Jefferson University, United StatesShizhi Wang, Southeast University, China
Fangrong Yan, China Pharmaceutical University, China
Copyright © 2021 Liu, Chen, Yang, Zhong, Ni, Yang, Shao, Cai, Bai and Yu. 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: Hao Yu, haoyu@njmu.edu.cn; Jianling Bai, baijianling@njmu.edu.cn
†These authors have contributed equally to this work