Abstract
Background:
Uterine corpus endometrial carcinoma (UCEC) is a gynecological malignant tumor with low survival rate and poor prognosis. The traditional clinicopathological staging is insufficient to estimate the prognosis of UCEC. It is necessary to select a more effective prognostic signature of UCEC to predict the prognosis and immunotherapy effect of UCEC.
Methods:
CIBERSORT and weighted correlation network analysis (WGCNA) algorithms were combined to screen modules related to regulatory T (Treg) cells. Subsequently, univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were used to identify the genes in key modules. The difference in overall survival (OS) between high- and low-risk patients was analyzed by Kaplan–Meier analysis. The Tregs-related risk signature (TRRS) was screened by uni- and multivariate Cox analyses. Afterward, we analyzed the expression difference of TRRS and verified its ability to predict the prognosis of UCEC and the effect of immunotherapy.
Results:
Red module has the highest correlation with Tregs among all clustered modules. Pathways enrichment indicated that the related processes of UCEC were primarily associated to the immune system. Eight genes (ZSWIM1, NPRL3, GOLGA7, ST6GALNAC4, CDC16, ITPK1, PCSK4, and CORO1B) were selected to construct TRRS. We found that this TRRS is a significantly independent prognostic factor of UCEC. Low-risk patients have higher overall survival than high-risk patients. The immune status of different groups was different, and tumor-related pathways were enriched in patients with higher risk score. Low-risk patients are more likely take higher tumor mutation burden (TMB). Meanwhile, they are more sensitive to chemotherapy than patients with high-risk score, which indicated a superior prognosis. Immune checkpoints such as PD-1, CTLA4, PD-L1, and PD-L2 all had a higher expression level in low-risk group. TRRS expression really has a relevance with the sensitivity of UCEC patients to chemotherapeutic drugs.
Conclusion:
We developed and validated a TRRS to estimate the prognosis and reflect the immune status of UCEC, which could accurately assess the prognosis of patients with UCEC and supply personalized treatments for them.
Introduction
Uterine corpus endometrial carcinoma (UCEC) is the third most common gynecological malignant cancer globally (). In 2019, the incidence and mortality of UCEC have been estimated at 61,880 and 12,160, respectively, merely in the United States (). Obesity can increase the risk of uterine corpus endometrial carcinoma (). At present, multiple treatment such as surgery, chemotherapy, radiation therapy, and hormone therapy are always applied to UCEC treatment, but the incidence and disease-related mortality are still increasing annually (–). Effective treatment is based on accurate assessment of prognosis. Nevertheless, patients in the same clinical stage may present different clinical characteristics, indicating that the prognosis of UCEC according to the traditional clinicopathological staging is not fully accurate (). Consequently, using effective biomarkers to accurately define different UCEC stage is helpful to precise treatment. In recent years, immunotherapy has become an effective therapy for cancer (), especially in melanoma, lung cancer, and liver cancer (). Tregs, a category of CD4+ T cells, can maintain immune homeostasis through regulating antimicrobial resistance, allergy, and transplantation rejection, and suppressing protective immune responses (). It has been confirmed that Tregs could be employed to predict the outcomes of solid tumors like breast and ovarian (, ). Checkpoint activity of T cells can be suppressed by immune checkpoint inhibitor (ICI). In particularly, CTLA‐4 and PD‐1 are important ICI (, ) and have shown good efficacy in cancer treatment (). More recently, olaparib has shown a clinical effect on UCEC to certain degree (); the relationships between olaparib exposure and UCEC biomarkers are still unknown (). The expression of RNF183, which has been considered as a good prognostic marker of UCEC, has been confirmed to be related to the markers of different subsets of Tregs (). Therefore, the identification of TRRS will be helpful to explore the immunotherapy of UCEC.
With the continuous development of bioinformation technology, a great deal of methods has been employed to define biomarkers. WGCNA algorithm can identify highly correlated modules and genes for cancer based on the network construct of genes expression (). This study applied WGCNA to identify relevant modules and genes of UCEC, and the red module was selected for further analysis. In addition, we construct a prognostic model related to Tregs and analyzed its relationship with immune microenvironment, chemotherapy, and immunotherapy. Furthermore, the TRRS might be used as a novel tool to diagnose UCEC patients and provide more effective personalized treatment.
Materials and Methods
Data Acquisition
A total of 552 UCEC cases and 23 normal samples were obtained from The Cancer Genome Atlas (TCGA) data portal (https://portal.gdc.cancer.gov/) (). Clinical factors data such as age, grade, and histological type were downloaded from TCGA data portal.
Estimation Tregs Proportion and Differential Analysis in UCEC
CIBERSORT was employed to estimate the proportion of 22 immune cells for normal and tumor tissues (). Therefore, we identified 309 samples with the standards of p < 0.05 from 575 cases of samples.
Construction of Co-expression Network and Select Hub Module
Intersecting Gene and Low-Tregs Groups. Second, we used Limma to define the differentially expressed genes with a p-value lower than 0.05 and log foldchange larger than 1 as the threshold (). We obtained 4,703 genes. Then, these genes were used to construct a WGCNA (). First, according to the Pearson’s correlation value, a weighted matrix was constructed. Next, amn = |cmn| β (cmn means Pearson’s correlation value of paired genes; amn is adjacency between paired genes) was used to construct a weighted adjacency matrix. Parameter β is a soft threshold, has the function of strengthening correlations, and reduce weakening correlations between genes. The value of β was defined as 6. For the purpose of dividing the genes with resembling expression levels into different modules. We cluster these genes with minimum size genes dendrogram of 50. The key module was identified through the analysis of the correlation between these genes and clinical factors. Eigengene dendrogram and adjacency heatmap also confirmed that the red module has a higher correlation coefficient. A protein–protein interaction (PPI) network was constructed base on STRING database, a tool that integrates all the connections between proteins of interacting genes (), intuitively showing the interacting nodes of each genes in the red module.
Gene Set Enrichment Analysis
Enrichr (https://maayanlab.cloud/Enrichr/) is a web-based program that has plenty information about protein and gene, in which investigators always use to do pathway and process enrichment analysis (). For further analysis, Enrichr was applied to conduct Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to analyze the functions of genes in the red module.
Construction of TRRS
We obtained 537 samples with complete gene expression profiles and OS time. We randomly assigned 269 patients as the train set based on a computer-generated allocation sequence and took the entire set as the validation data. The train set was employed to construct TRRS, while the entire set was chosen for validating the predicting value. LASSO analysis based on R package “glmnet” combined with multivariate Cox regression analyses had been employed to select genes that has significant connections with OS of UCEC. Afterwards, these genes were applied to construct TRRS (). The coefficients calculated by LASSO regression were applied to obtain the formula as follows: risk score = sum of coefficients × 8 TRRS expression level (). All the samples were split into two groups in the light of the median risk score. Through Kaplan–Meier curves, receiver operating characteristic curves (ROC), and principal component analysis, we assess the accuracy of TRRS. Risk score distribution, survival status, and genes expression were also taken into consideration. In addition, univariate and multivariate Cox regression analyses were applied to calculate the predictive ability of TRRS.
Quantitative Real-Time RT-PCR
Total RNA from 16 UCEC samples and 16 normal tissues was extracted on the basis of Trizol reagent (Invitrogen). Before reverse transcription to cDNA, 4× GDNA wiper mix (vazymer323-01) was employed to remove residual genomic DNA from total RNA. Complementary RNA was synthesized by using PrimeScript RT reagent kit. Real-time quantification was performed using the SYBR Premix Ex Taq Kit (TaKaRa DRR041). The relative expression level of the target gene was standardized by GAPDH and 2−△△Ct method. Quantitative real-time RT-PCR (qRT-PCR) primers are listed in Supplementary Table S1.
Construction of a Nomogram
Five clinical characteristics such as age, stage, grade, histological type, and risk score were combined to establish a nomogram aims to calculate the OS of 1, 3, and 5 years of UCEC in the entire set. Calibration curves were employed to evaluate whether the established nomogram is reliable ().
Immune and Clinical Characteristic Identification
A total of 552 UCEC samples and 23 normal samples were obtained from TCGA database, which was employed to verify the expression distinction between tumor and normal tissues. Furthermore, boxplots were drawn to analyze the relationship between gene expression and clinical factors. Finally, according to the level of expression, patients were split into two groups. Then, we verified the connection between TRRS expression and survival probability based on K–M analysis.
Genome-Wide Analysis of Genes
Gene Set Cancer Analysis (GSCALite), a web server, aims to flexibly calculate the expression, mutation, and interaction of genes in different cancers (). GSCALite has the function to analyze different gene expression levels, survival time, mutations, methylation, chemotherapeutics sensitivity, and so on. We analyze mutation distribution and global activity of eight genes. After that, we identified the connection between genes expression levels and copy number variations (CNVs) and methylation by using GSCALite ().
Survival Analyses of Genes in TRRS
The areas under the curve (AUCs) of 1, 3, and 5 years on account of risk score and different clinical factors were calculated, respectively. Furthermore, the curve was employed to identify the total influence of risk score combined with clinical characteristics on survival probability. Overall survival rates of different clinical characteristics were calculated as well.
IPS Analysis
Immunomodulators, major histocompatibility complex (MHC) molecules, effector cells, and immunosuppressive cells are the four main components to evaluate tumor immunogenicity. Through calculating the expression values of four kinds of immune genes, we derived a patient’s immunophenoscore (IPS). The IPS was achieved based on a scale range from 0 to 10. The higher the score, the stronger the immunogenicity. The IPSs (including IPS, IPS-CTLA4, IPS-CTLA4/PD-L1/PD1/PD-L2, and IPS-PD1/PD-L1/PD-L2 scores), acquired from The Cancer Immunome Atlas [TCIA (https://tcia.at/home)] (), were applied to assess the response of UCEC patients for ICI.
Estimation of ICIs Response
The connection between ICI expression and risk score of UCEC patients was identified based on Pearson correlation coefficient method. Boxplots display the results.
GSEA and ESTIMATE and ssGSEA
Function and pathways for different risk groups were investigated by GSEA, respectively (). Pathways with nominal p < 0.05 were considered significantly enriched. Immune scores, stromal scores, and estimate scores of each sample were received based on the “Estimation of Stromal and Immune cells in Malignant Tumors using Expression data” (ESTIMATE) algorithm (). Additionally, the difference in immune cell and immune function between two risk groups were achieved by CIBERSORT method. Single-sample gene set enrichment (ssGSEA) is another method to identify the distinction of immunological status between patients in different risks. ssGSEA utilized the enrichment scores to indicate the degree of absolute enrichment in each sample. Standardized enrichment scores for each immune category could be calculated ().
Mutation Analysis
The mutation materials of UCEC were obtained from TCGA. Somatic mutation data are stored in mutation annotation format (MAF) (). The TMB scores of every patient were calculated based on the following formula: TMB = (total mutation/total covered bases) × 106.
Chemotherapeutic Response Prediction
A total of 537 samples were split into high- and low-risk groups on the basis of median risk score, and each group was given six kinds of chemotherapeutics including cisplatin, docetaxel, doxorubicin, gemcitabine, methotrexate, and paclitaxel. The sensitivity of each sample to chemotherapy was predicted by using Genomics of Drug Sensitivity in Cancer (GDSC) (https://www.cancerrxgene.org/) (). The calculation was conducted through R package “pRRophetic”, where the samples’ half-maximal inhibitory concentration (IC50) was achieved by ridge regression. Pearson correlation coefficient was applied to explore the relevance between gene expression in TRRS and their sensitivity to pharmacotherapy.
Immunoassay for ZSWIM1
For the purpose of further studying the correlations between ZSWIM1 and immunity, tumor–immune system interactions and drugbank (TISIDB, http://cis.hku.hk/TISIDB/index.php) database were applied and inferred the correlation between ZSWIM1 expression and the TRRS of UCEC ().
Statistical Analysis
R (version 4.0.5) was applied to perform statistical analysis in our study. All statistical tests were two-sided, and p < 0.05 was considered as statistical difference. Student’s t-test was used to make a comparison between the normally distributed variables in the two groups, and Wilcoxon test was used to calculate the continuous variables.
Results
Tregs in Uterine Corpus Endometrial Carcinoma
The study indicated that a lower expression of Tregs was related to higher clinical grade and more serious pathological morphology (p = 4.756e−09, p = 6.05e−04, respectively) (Figures 1A, B). Meanwhile, patients in advanced stage had a lower expression of Tregs in general (p = 0.001) (Figure 1C). Additionally, the prognosis of UCEC patients with lower Tregs expression is poorer than those with higher Tregs expression (Figure 1D).
Figure 1
Construct Weighted Co-Expression Network
A total of 4,703 genes were applied to set up a weighted co-expression network based on WGCNA. First, combining with clinical factors, we clustered 298 samples and constructed samples dendrogram and trait heatmap for them (Figure 2A). Then, we identified the value of β = 6 (R2 = 0.91) as the threshold (Supplementary Figures S1A, B). Supplementary Figures S1C, D indicate positive result of rationality test. In addition, a hierarchical clustering tree had been established, and 12 modules were generated (Figure 2B). Genes having relative expression were clustered together to form a branch that constructed a module. Among the 12 modules formed by clustering, the correlation between the red module and Tregs is higher than that between other modules, which have 251 genes (Figure 2C). Figures 2D, E show that the modules are not independent of each other. Finally, PPI network was constructed by using STRING (Figure 2F), which indicates that most genes are connected with the others in the red module. Figure 2G shows the top 30 genes that have more connections with other genes.
Figure 2
Pathway and Process Enrichment Analysis
We employed GO and KEGG analysis to conduct enrichment analysis and find out the functions and pathways that related to the red module. The outcome indicated that the functions and pathways of UCEC were primarily bound with immune-related physiological processes (Supplementary Figure S2).
Key Genes Identification and TRRS Construction
In order to estimate the prognostic performance of DE genes, univariate Cox regression analysis was employed to the 251 genes (Supplementary Table S2). In the train set, 12 genes were closely related to the OS of UCEC patients (p < 0.05). Twelve genes were analyzed by LASSO, and eight of them were selected (Supplementary Figures S3A, B). After this, a multivariate Cox regression analysis was implemented, and five genes had significant statistical correlation with the hazard ratio of UCEC patients in the train set (Supplementary Figure S3C). We utilized these genes to establish the TRRS. Eight genes were weighted by relative coefficient, and the formula is as follows: risk score = (−0.070*CDC16) + (0.010*ZSWIM1) + (−0.045*ITPK1) + (0.099*NPRL3) + (0.274*GOLGA7) + (0.022*ST6GALNAC4) + (−0.201*PCSK4) + (−0.026*CORO1B). These eight genes were associated with high risk, including ZSWIM1 [hazards ratio (HR) =1.104(1.009–1.208), p = 0.031], NPRL3 [HR = 1.104(1.023–1.191), p = 0.010), GOLGA7 [HR = 1.027(0.995–1.061), p = 0.096], ST6GALNAC4 [HR = 1.022(0.998–1.046), p = 0.072], CDC16 [HR = 0.932(0.881–0.996), p = 0.015], ITPK1 [HR = 0.956(0.922–0.991), p = 0.015], PCSK4 (HR = 0.818[0.656–1.020], p = 0.072], and CORO1B [HR = 0.974(0.951–0.998), p = 0.033) (Table 1). The risk score of samples in the train set were computed according to the above formula. People in the train set were split into high-risk group (n = 134) and low-risk group (n = 135) on the basis of their median-risk score. The comparison showed that there exists a significant difference in the OS of different groups (p = 2.74e−06, log-rank test) (Figure 3A). The AUC for the signature of OS in 5 years is 0.753, and it is 0.836 in 3 years and 0.781 in 1 year (Figure 3B). We sort the patient’s risk scores, and their distribution is shown in Figure 3C. The living condition of UCEC patients is shown in the dot plot (Figure 3D). Gene expression pattern between the two groups of patients with different prognosis is presented in the heatmap expression (Figure 3E). Low- and high-risk patients showed significant difference using principal component analysis (PCA) (Figure 3F).
Table 1
| Gene | HR (95% CI) | p value |
|---|---|---|
| CDC16 | 0.93 (0.88, 0.99) | 0.01 |
| ZSWIM1 | 1.10 (1.01, 1.21) | 0.03 |
| ITPK1 | 0.96 (0.92, 1.00) | 0.02 |
| NPRL3 | 1.10 (1.02, 1.19) | 0.01 |
| GOLGA7 | 1.03 (1.00, 1.06) | 0.10 |
| ST6GALNAC4 | 1.02 (1.00, 1.05) | 0.07 |
| PCSK4 | 0.81 (0.66, 1.02) | 0.07 |
| CORO1B | 0.97 (0.95, 1.00) | 0.03 |
Multivariate Cox regression eight genes weighted by their relative coefficient.
Figure 3
The Relationships Between TRRS Expression and Clinical Factors
We achieved the expression profiles of TRRS from the TCGA. GOLGA7, ITPK1, and ST6GALNAC4 genes had lower expression in tumor tissues than in healthy persons (p < 0.05) by means of Wilcoxon signed-rank test (Supplementary Figure S4A). After that, qRT-PCR has been used to compare the difference in TRRS expression between normal and tumor tissue (Supplementary Figure S4B). The results showed that the expression levels of CORO1B, GOLGA7, PCSK4, ST6GALNAC4, and ZSWIM1 in normal tissues were significantly higher than those in tumor tissues, which were similar to the trends detected in the TCGA dataset. Then, the stratified studies were utilized to see whether the genes have different expression level in different clinical characteristics (Supplementary Figure S5). It is not hard to find that most of the eight TRRS have lower expression levels in advanced stage. The result showed that higher expression of GOLGA7 [HR =2.28(1.46–3.56), p < 0.001] and ZSWIM1 [HR = 1.7(1.09–2.67), p = 0.011] were combined with poorer OS of UCEC patients; however, CORO1B [HR = 0.51(0.33–0.79), p = 0.001], CDC16 [HR = 0.52(0.33–0.81), p = 0.001], PCSK4 [HR = 0.51 (0.32–0.81), p = 0.001], and ITPK1 [HR = 0.52(0.34–0.8), p = 0.002] expression level had a positively correlated with overall survival (Supplementary Figure S6).
Genome-Wide Analysis of TRRS
Genome-wide analysis of TRRS was carried out by utilizing the GSCALite. The results indicated that PCSK4 was the gene with the highest mutation frequency, followed by CDC16 and ITPK1, while GOLGA7 had the lowest (Figure 4A). The expression of the eight genes in TRRS had a positive correlation with copy number variations (CNVs), which is performed in the bubble diagram (Figure 4B). CNV frequency was positively correlated with gene expression. Spearman correlation coefficient between methylation and gene expression was performed in the methylation difference bubble chart. We found that the GOLGA7, PCSK4, ZSWIM1, ST6GALNAC4, and ITPK1’s methylation was downregulated in UCEC (Figure 4C). This means that methylation had a significant effect on gene expression. The exploration of pathway activity showed that ZSWIM1, PCSK4, CDC16, and CORO1B had positive correlations to DNA damage response pathway activation; ITPK1, NPRL3, PCSK4, ST6GALNAC4, and ZSWIM1 were related to the RTK pathway inhibition; and ITPK1 activated the apoptosis and inhibit PI3K/AKT pathway significantly (Figure 4D).
Figure 4
Verifying the Predictive Capability of the Eight TRRS
For the purpose of exploring the predictive power of TRRS, we constructed an entire set. In the entire set, the risk score was computed by using the median risk score. Each patient in the entire set was split into two groups according to the risk scores. High-risk groups have 267 cases and low-risk groups have 270. There was significant statistical distinction between the two groups’ Kaplan–Meier survival curves (Figure 3G, P =7.44e-07). In the entire set, the AUC was 0.797 in 1 year, 0.702 in 3 years, and 0.7 in 5 years (Figure 3H). The situation of risk score, survival status, and expression of eight TRRS in the entire set are performed in Figures 3I–K. Principal component analysis also showed a discrepancy between the groups (Figure 3L).
Using Nomogram to Predict the Survival Rates
By using the known risk score and some clinical characteristics, multivariate logistic regression was applied to structure a nomogram that may predict the survival rates of UCEC patients accurately. Age, stage, grade, histological type, and risk score were considered predictors of survival rates, which were combined into the nomogram (Figure 5A). It showed that risk score is the most influential factors of the nomogram total score. The calibration curve of the constructed nomogram (Figures 5B–D) performed that the survival rates predicted by TRRS was almost consistent with the actually survival rates we observed. This means the TRRS has good clinical practicability.
Figure 5
Clinical Factors and Risk Score Will Affect the Prognosis
AUC was employed to calculate the prognostic capability of clinical factor and risk score. The higher the AUC, the more precise the TRRS. Risk score and four clinical factors are shown on Figures 6A–C. The AUCs of the clinical factors combined with the risk score of 1-, 3- and 5-year survival rates were 0.751, 0.758, and 0.758, respectively (Figures 6D–F). The result showed that using the risk score combined with clinical characteristics to evaluate the prognosis resulted in high sensitivity and specificity. After that, Cox regression analyses including univariate and multivariate regression were applied to calculate the prognosis capacity of risk score and clinical factors. As performed in Figures 6G, H, the outcome of univariate cox regression analysis in the entire set showed that stage (HR, 3.881; 95%CI, 2.561–5.883, p < 0.001), histological type (HR, 2.836; 95%CI, 1.874–4.291; p < 0.001), and risk score (HR, 1.050; 95%CI, 1.028–1.072; p < 0.001) were related to UCEC prognosis. However, the outcome of multivariate Cox regression indicated that stage (HR, 2.836; 95%CI, 1.874–4.291; p < 0.001) is an independent prognostic factor of UCEC. Besides, the stage (HR, 2.785; 95% CI, 1.559–4.973; p < 0.001), grade (HR, 3.729; 95% CI, 1.739–7.997; p < 0.001), and risk score (HR, 1.152; 95%CI, 1.106–1.200; p < 0.001) were related to the prognosis of UCEC. However, in the train set, only the risk score (HR, 1.171; 95%CI, 1.066–1.170; p < 0.001) was independently correlated with the prognosis of UCEC patients (Figures 6I, J). Then we used prognostic stratification to analyze the function of risk score to judge the prognosis of UCEC. The result indicated that with different clinical characteristics (age, stage, grade, and endometrial histological type), UCEC patients with different risks have different overall survival rates (Supplementary Figure S7). This result also proved that risk score can affect the prognosis of UCEC patients independently.
Figure 6
GSEA and ESTIMATE and ssGSEA Analysis
Through GSEA, we concluded that the high risk is related to tumor-related pathways (Figure 7A), which can explain why high-risk groups have poor prognosis. Low risk in patients is associated with the pathways that related to immunity (Figure 7B). Therefore, we use ESTIMATE and ssGSEA analysis to verify the difference in immunological status between the groups. The immune score, stromal score, and estimate score were achieved by ESTIMATE algorithm through R “estimate” package (Figures 7C–H). The study on the association between risk score and immunity based on ssGSEA also confirmed that the risk score was negatively correlated with its immune ability (Figures 7I, J).
Figure 7
The Relationship Between the TRRS and Immune Response in UCEC
According to our previous studies, we comprehensively assess the values of 22 immune cells by using CIBERSORT. Figure 8A shows the results we got from 251 UCEC samples. The result showed that dendritic cells activation, macrophage M0, resting NK cells, T-follicular helper cells, and gamma delta T cells were obviously stronger in the high-risk group. Resting dendritic cells, monocytes, neutrophils, CD8 T cells, activating NK cells, and Tregs were obviously higher in low-risk patients (Figure 8B). The results of the relationship analysis were consistent with this. Activated dendritic cells (R = 0.23, p < 0.001), macrophage M1 (R = 0.14, p = 0.029), and macrophage M2 (R = 0.15, p = 0.021), T follicular helper cells (R = 0.15, p = 0.020), and gamma delta T cells (R = 0.18, p = 0.005) were positively associated with the risk score. While resting dendritic cells (R = −0.17, p = 0.009), monocytes (R = −0.16, p = 0.014), activating NK cells activated (R = −0.19, p = 0.003), CD8 T cells (R = −0.14, p = 0.023) were positively associated with the risk score. Specially, the risk scores were strongly negatively related to the expression of Tregs (R = −0.33, p = 8.5e−08). This is consistent with the previous results (Figure 8C).
Figure 8
IPS and Immune Checkpoints
Immune checkpoints (CTLA4 and PD1) are able to assess the response of patients to immunotherapy. The result showed that the expression of CTLA4 (entire set: R = −0.16, p = 0.00025; train set, R = −0.16, p = 0.0098) and PD1 (entire set: R = −0.15, p = 0.00042; train set: R = −0.19, p = 0.0023) were negatively correlated with TRRS in both train set and entire set (Figures 9A–D). In addition, patients in lower risk presented higher gene expression of CTLA4 and PD1 in both train set and entire set (p < 0.05) (Figures 9E–H). Consequently, we supposed that immune checkpoints related to Tregs might be used in immunotherapy of UCEC. In this study, we further analyzed the correlation between IPS and TRRS in UCEC. The IPS, IPS-CTLA4, IPS-CTLA4/PD-L1/PD1/PD-L2, and IPS-PD1/PD-L1/PD-L2 scores were employed to assess the probability of ICI. Low-risk patients had significantly higher scores (Figures 9I–L: IPS, p < 0.001; IPS-CTLA4, p = 0.003; IPS-CTLA4/PD-L1/PD1/PD-L2, p = 0.008; and IPS-PD1/PD-L1/PD-L2, p = 0.027). Consequently, we inferred that patients with low risk are more likely to trigger an immune response.
Figure 9
The TRRS and Mutation Profile
Tumor burden has always been an important factor affecting immunotherapy (). In our study, TMBs are negatively correlated with TRRS (Figure 10A). Genes that had the most frequent mutation in each groups are shown in Figures 10B, C. Low-risk patients had heavier TMB (p = 0.032) (Figure 10D). Furthermore, survival probability of lower TMB patients was significantly lower than those of patients with higher TMB (Figure 10E). Meanwhile, patients with lower tumor mutational burden and high risk had the lowest survival probability compared with other groups (p < 0.001) (Figure 10F).
Figure 10
Relationship Between TRRS and Chemotherapy Sensitivity
Recently, chemotherapy is an ordinary treatment for UCEC; we analyzed the response of two groups to six chemotherapeutic drugs including cisplatin, docetaxel, doxorubicin, gemcitabine, methotrexate, and paclitaxel. We calculated IC50 for each samples using the TRRS. The results showed that half of those drugs have significant distinction between the groups. Doxorubicin (p = 0.001) and gemcitabine (p = 0.004) have higher sensitivity in the low-risk group, which indicated that chemotherapeutic drugs may have better curative effect in the high-risk group (Figure 11A). After that, we further analyzed the relationship between expression of genes in TRRS and chemotherapy sensitivity. The outcome indicated that seven genes were strongly correlated to the sensitivity of some chemotherapeutic drug (p < 0.01) (Figure 11B). For example, GOLGA7, ITPK1, and NPRL3 were correlated with increased drug resistance of cancer cells to Dasatinib, Zoledronate, PF-06463922, Brigatinib, LDK-378, Vinorelbine, Carfilzomib, and Bortezomib, respectively. Meanwhile, increased expression of ST6GALNAC4 and PCSK4 was related to more sensitivity of tumor cells to a number of chemotherapy drugs like Temsirolimus, Bleomycin, Nelarabine, and Cladribine. In addition, CDC16 was positively correlated with ARRY-162 and Selumetinib but was negatively correlated with Everolimus. The mechanism needs further study.
Figure 11
Immunoassay for ZSWIM1
Among the eight genes, only the expression of ZSWIM1 was consistent with the prognosis, so we chose ZSWIM1 for further analysis. Tumor-infiltrating lymphocytes (TILs) can undertake the role of an independent predictor in several cancers. Consequently, TISIDB was applied to analyze the relationship between TILs and ZSWIM1expression. The correlation between ZSWIM1 expression and TILs in different descriptions of cancer is performed in Figure 12A. The expression of ZSWIM1 was negatively correlated with Tregs in UCEC (R = −0.161, p = 0.000168). Immunopotentiators, immunosuppressants, MHC molecules (Figure 12B), chemokines, and receptors (Figure 12C) had a negative relationship with the expression of ZSWIM1. There exist significant distinction in the expression of ZSWIM1 among different immune and molecular subtypes of UCEC. The relationship between ZSWIM1 expression and human cancer immune subtypes is shown in Figure 12D. Specifically, the expression of ZSWIM1 in wound healing, interferon gamma (IFN-γ) dominance, inflammation, lymphocyte depletion, and transforming growth factor alpha (TGF-α) was increased in turn. However, it was negatively correlated with molecular subtypes.
Figure 12
Discussion
UCEC is a common gynecological malignant tumor () with low survival rate and poor prognosis (37–39). However, if the disease is correctly diagnosed in the early stage, the 5-year survival rate can be as high as 90% (40). At present, the treatment of EC is based on surgical operation (41), postoperative radiotherapy, and chemotherapy, which can be performed according to the type and stage of tumor (42). However, patients in same clinical stage may manifest different clinical characteristics, indicating that it is insufficient to estimate the prognosis of UCEC according to the traditional clinicopathological staging (). In recent years, CA125 and HE4 have been used as serum markers of UCEC. However, the accuracy of serum prediction is relatively low (43). Clinical trials of UCEC patient showed that the effective rate of mammalian target of rapamycin (mTOR) inhibitors treatment was <10% (44). Immunotherapy has replaced traditional radiotherapy and chemotherapy as a remedy method for cancer. Especially, CTLA-4, PD-1, and PD-L1 antibodies have curative effect on tumor treatment (45). In this study, The Cancer Genome Atlas (TCGA) database was employed to construct TRRS for patients with UCEC, which is expected to be used to monitor prognosis and immunotherapy response in UCEC.
In the immunotherapy of different kinds of cancer, blocking immune checkpoint has shown a wide range of effects (46). Antibodies against checkpoint molecule, namely, CTLA4, PD1, and PD-L1 have shown clinical efficacy and persistence in more than 15 kinds of human malignant caners (47, 48). This has been approved by the Food and Drug Administration (FDA) as a treatment strategy for a variety of cancers. In particular, anti-PD-1 combination therapy can be applied as a remedy for advanced melanoma, which has good curative effect compared with any single drug therapy (49).
At the beginning, CIBERSORT was employed to assess the proportion of T cells, which indicated that the more advanced the clinical stage and grade, the lower the expression of Tregs. This means that TRRS is a tumor suppressor gene. By using the WGCNA algorithm and correlation analysis, the red module was selected as the most relevant module. GO and KEGG enrichment analysis showed that the red module has the closest relationship with Tregs. Univariate, LASSO, and multivariate Cox regression were applied to calculate prognostic significances in risk score and some clinical factors.
Genome-wide analysis of these eight genes indicated that these high frequency mutations and methylation had close relationships with gene expression, and they were involved in regulating the activity of injury response pathway. The high expression of these genes in tumor samples and early clinical stage indicates their potential use as biomarkers. The outcomes of qRT-PCR were similar to the trends detected in the TCGA dataset, which confirmed the prediction ability of the model to a certain degree. The AUC of the risk score showed the high sensitivity and specificity of TRRS. The AUCs of the risk score merged with clinical characteristics were higher, which suggests that using the TRRS together with clinical factors, we can give patients a better treatment guidance. In line with previous studies, the results of this study showed that age, stage, and grade were also connected with the prognosis of patients with UCEC. The nomogram shows that the TRRS can accurately assess the prognosis of UCEC (). Studies on several common chemotherapeutic drug reactions have found that low-risk patients had higher sensitivity to adriamycin and gemcitabine; furthermore, the expression of CTLA-4, PD1, and PDL1 was higher in the low-risk group, which means that chemotherapeutic drugs could have better efficacy in low-risk patients. TMB was used to assess the ability of TRRS to identify patients who have higher response to ICI (). TMB was higher in low-risk groups with UCEC. This means that low-risk patients have the possibility to be identified by immune cells and can benefit more from immunotherapy. In addition, the survival probability of low TMB cases was significantly lighter than that of cases with high TMB.
Eight genes were selected for further analysis and had been used as prognostic markers in other diseases. CDC16 has a connection with multiple neurodevelopmental disorders (50, 51). It has been found that CDC16 has potential therapeutic function in melanoma (52). ITPK1 may serve as biomarkers for GC pathogenesis (53, 54). The expression of PCSK4 mRNA is decreased in non-islet-cell tumor hypoglycemia (NICTH) patients, which has been demonstrated to be related to serum big IGF2 increase (55). Coronin 1B (Coro1B) is one of the actin binding proteins that can regulate platelet-derived growth factor (PDGF)-induced vascular smooth muscle cell (VSMC) migration, suggesting a new therapeutic target for vasculopathies (56). ZSWIM1 can be used as a biomarker of T helper cell differentiation (57). Nitrogen permease regulator-like 3 took part in the construct of GATOR1 complex (Supplementary Figure S1), which can regulate the mTOR pathway (58). Chang Soo Ryu found that NPRL3 is a common biomarker for ischemic stroke (59). Linghui Zhou found that rs11337 (G > T) in GOLGA7 is related to survival of glioma patients (60). ST6GALNAC4 expression is related to glycosphingolipids synthesis, which has a connection with breast cancer (61, 62). By using GSEA, tumor-related pathways like ERBB, TGF-BETA, and WNT were significantly enriched; these pathways are deemed to be associated with tumor, which could be used as novel therapeutic targets (63–65). Immune-related responses were significantly enriched in low-risk patients, which further validated the difference in immune status between the two risk groups. Besides, high-risk cases had higher fractions of monocytes, NK cells, CD8 T cells, neutrophils, and Tregs. Researchers have confirmed that NK cells, CD8 T cells, neutrophils, and Tregs were significantly associated with survival of UCEC patients (66).
Genome-wide analysis of the genes showed that the mutation frequency of PCSK4 was the highest. The expression of the TRRS was positively correlated with CNV. Besides, the methylation of GOLGA7, PCSK4, ZSWIM1, ST6GALNAC4, and ITPK1 was downregulated. DNA methylation is an epigenetic mechanism to control the expression of oncogenic or tumor-suppressive genes. Scientists have studied the potential effectiveness of some methylated biomarkers in predicting cancer prognosis (67). However, there has been seldom studies on these genes’ methylation in UCEC, so it makes sense to do further research. Further study on ZSWIM1 based on TISIDB showed that the expression level of ZSWIM1 in different immune subtypes and molecular subtypes of UCEC was significantly different. The immunosuppressant was effective on ZSWIM1 gene. Therefore, ZSWIM1 can be considered as a target for UCEC treatment.
Previous studies have also found gene prognostic markers closely related to UCEC (68–70), but it is the first time to construct a model in Treg cells to predict the prognosis of UCEC. AUC is a crucial standard to judge whether a prediction model has good discrimination. In another article that also identified a signature for predicting the prognosis of patients with UCEC, the AUC of the prognostic model using 10 immune genes was 0.756 (). In this study, the AUC for the signature that we constructed in 1 year is 0.781, and it is 0.836 in 3 years and 0.753 in 5 year. It is higher than that of a previous study that showed high sensitivity and specificity of our model, suggesting that the model has a better ability to predict the probability of disease occurrence. Nonetheless, this study still has some deficiencies. The conclusion of our study is mainly on account of bioinformatics analysis, and further clinical research is needed. Furthermore, risk factors of UCEC, such as obesity and smoking, were not discussed in this study.
Conclusion
All in all, through a series of bioinformatics analysis, we constructed a TRRS as potential biomarkers and targets for immunotherapy of UCEC. Low-risk patients had better prognosis and higher response rate to ICI. In the future, TRRS is expected to help predict prognosis and assess the efficacy of immunotherapy for UCEC patients, which can provide individualized treatment.
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.
Statements
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.
Author contributions
JB, JL, and RG conceived the study and participated in the study design, performance, and manuscript writing. JL, SY, and FS conducted the bioinformatics analysis. ZZ, MY, SN, and LC revised the manuscript. All authors contributed to the article and approved the submitted version.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.788431/full#supplementary-material
Supplementary Figure 1Soft-thresholding parameters of WGCNA. (A) Analyze the scale-free fit index for various soft-thresholding parameters. (B) Analyze the mean connectivity for various soft-thresholding parameters. (C) Histogram of connectivity distribution when β=6. (D) Check of scale-free topology when β=6.
Supplementary Figure 2GO and KEGG analysis (A) Top ten biological process related to key gene using GO analysis. (B) Top ten cellular component related to key genes using GO analysis. (C) Top ten molecular function related to key genes using GO analysis. (D) Top ten human pathway related to key genes using KEGG analysis.
Supplementary Figure 3LASSO COX regression and multivariate Cox analysis. (A, B) LASSO COX regression for OS of TRRS in the TCGA cohort. (C) Forest plot illustrating the multivariate Cox model results of eight gene related to Tregs.
Supplementary Figure 4Comparison of the expression of the eight genes. (A) Expression levels of the eight genes between EC samples and normal tissues evaluated by means of Wilcoxon signed-rank test. (B) Expression levels of the eight genes between normal tissues and tumor tissue evaluated by using qRT-PCR.
Supplementary Figure 5The relationships between gene expression and clinical characteristic.
Supplementary Figure 6The association of expression of eight gene and OS in TCGA UCEC dataset. (A) GOLGA7; (B) CORO1B; (C) CDC16; (D) ZSWIM1; (E)PCSK4; (F)ITPK1.
Supplementary Figure 7Overall survival result between patients with high and low risk. (A) age>60; (B) age<60; (C) grade G1&G2; (D) grade G3&G4; (E)histological-type Mixed &Serous; (F) Stage I& Stage II; (G) Stage III& Stage IV.
References
1
RomeroIRubioMJMedinaMMatias-GuiuXSantacanaMSchoenenbergerJAet al. An Olaparib Window-of-Opportunity Trial in Patients With Early-Stage Endometrial Carcinoma: POLEN Study. Gynecol Oncol (2020) 159:721–31. doi: 10.1016/j.ygyno.2020.09.013
2
SiegelRLMillerKDJemalA. Cancer Statistics, 2019. CA Cancer J Clin (2019) 69:7–34. doi: 10.3322/caac.21551
3
LehrerSRheinsteinPH. Endometrial Cancer, BRCA1, and BRCA2 in the UK Biobank Cohort. Am J Clin Oncol (2021) 44:299–300. doi: 10.1097/COC.0000000000000833
4
LiBLWanXP. Prognostic Significance of Immune Landscape in Tumour Microenvironment of Endometrial Cancer. J Cell Mol Med (2020) 24:7767–77. doi: 10.1111/jcmm.15408
5
MakkerVTaylorMHAghajanianCOakninAMierJCohnALet al. Lenvatinib Plus Pembrolizumab in Patients With Advanced Endometrial Cancer. J Clin Oncol (2020) 38:2981–92. doi: 10.1200/JCO.19.02627
6
SoslowRATornosCParkKJMalpicaAMatias-GuiuXOlivaEet al. Endometrial Carcinoma Diagnosis: Use of FIGO Grading and Genomic Subcategories in Clinical Practice: Recommendations of the International Society of Gynecological Pathologists. Int J Gynecol Pathol (2019) 38(Suppl 1):S64–74. doi: 10.1097/PGP.0000000000000518
7
ZhouCLiCYanFZhengY. Identification of an Immune Gene Signature for Predicting the Prognosis of Patients With Uterine Corpus Endometrial Carcinoma. Cancer Cell Int (2020) 20:541. doi: 10.1186/s12935-020-01560-w
8
MillerKDNogueiraLMariottoABRowlandJHYabroffKRAlfanoCMet al. Cancer Treatment and Survivorship Statistics, 2019. CA Cancer J Clin (2019) 69:363–85. doi: 10.3322/caac.21565
9
HavelJJChowellDChanTA. The Evolving Landscape of Biomarkers for Checkpoint Inhibitor Immunotherapy. Nat Rev Cancer (2019) 19:133–50. doi: 10.1038/s41568-019-0116-x
10
ChenYHChenSHHouJKeZBWuYPLinTTet al. Identifying Hub Genes of Clear Cell Renal Cell Carcinoma Associated With the Proportion of Regulatory T Cells by Weighted Gene Co-Expression Network Analysis. Aging (Albany NY) (2019) 11:9478–91. doi: 10.18632/aging.102397
11
Kucukgoz GulecUKilic BagirEPaydasSGuzelABGumurduluDVardarMA. Programmed Death-1 (PD-1) and Programmed Death-Ligand 1 (PD-L1) Expressions in Type 2 Endometrial Cancer. Arch Gynecol Obs (2019) 300:377–82. doi: 10.1007/s00404-019-05180-2
12
ForschnerABattkeFHadaschikDSchulzeMWeissgraeberSHanCTet al. Tumor Mutation Burden and Circulating Tumor DNA in Combined CTLA-4 and PD-1 Antibody Therapy in Metastatic Melanoma - Results of a Prospective Biomarker Study. J Immunother Cancer (2019) 7:180. doi: 10.1186/s40425-019-0659-0
13
de MiguelMCalvoE. Clinical Challenges of Immune Checkpoint Inhibitors. Cancer Cell (2020) 38:326–33. doi: 10.1016/j.ccell.2020.07.004
14
GockleyAAKolinDLAwtreyCSLindemanNIMatulonisUAKonstantinopoulosPA. Durable Response in a Woman With Recurrent Low-Grade Endometrioid Endometrial Cancer and a Germline BRCA2 Mutation Treated With a PARP Inhibitor. Gynecol Oncol (2018) 150:219–26. doi: 10.1016/j.ygyno.2018.05.028
15
GengRZhengYZhaoLHuangXQiangRZhangRet al. RNF183 Is a Prognostic Biomarker and Correlates With Tumor Purity, Immune Infiltrates in Uterine Corpus Endometrial Carcinoma. Front Genet (2020) 11:595733. doi: 10.3389/fgene.2020.595733
16
LangfelderPHorvathS. WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
17
BlumAWangPZenklusenJC. SnapShot: TCGA-Analyzed Tumors. Cell (2018) 173:530. doi: 10.1016/j.cell.2018.03.059
18
WaksAGStoverDGGuerrieroJLDillonDBarryWTGjiniEet al. The Immune Microenvironment in Hormone Receptor-Positive Breast Cancer Before and After Preoperative Chemotherapy. Clin Cancer Res (2019) 25:4644–55. doi: 10.1158/1078-0432.CCR-19-0173
19
LiuJZhouSLiSJiangYWanYMaXet al. Eleven Genes Associated With Progression and Prognosis of Endometrial Cancer (EC) Identified by Comprehensive Bioinformatics Analysis. Cancer Cell Int (2019) 19:136. doi: 10.1186/s12935-019-0859-1
20
LiuJHWuZSunRNieSPMengHYZhongYet al. Using Mrnasi to Identify Prognostic-Related Genes in Endometrial Carcinoma Based on WGCNA. Life Sci (2020) 258:118231. doi: 10.1016/j.lfs.2020.118231
21
SzklarczykDGableALNastouKCLyonDKirschRPyysaloSet al. The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/Measurement Sets. Nucleic Acids Res (2021) 49:D605–12. doi: 10.1093/nar/gkaa1074
22
ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10:1523. doi: 10.1038/s41467-019-09234-6
23
McEligotAJPoynorVSharmaRPanangadanA. Logistic Lasso Regression for Dietary Intakes and Breast Cancer. Nutrients (2020) 12. doi: 10.3390/nu12092652
24
LeeSGornitzNXingEPHeckermanDLippertC. Ensembles of Lasso Screening Rules. IEEE Trans Pattern Anal Mach Intell (2018) 40:2841–52. doi: 10.1109/TPAMI.2017.2765321
25
HoshinoNHidaKSakaiYOsadaSIdaniHSatoTet al. Nomogram for Predicting Anastomotic Leakage After Low Anterior Resection for Rectal Cancer. Int J Color Dis (2018) 33:411–8. doi: 10.1007/s00384-018-2970-5
26
LiuCJHuFFXiaMXHanLZhangQGuoAY. GSCALite: A Web Server for Gene Set Cancer Analysis. Bioinformatics (2018) 34:3771–2. doi: 10.1093/bioinformatics/bty411
27
ZhangWZhaoHChenJZhongXZengWLiZet al. Mining Database for the Expression and Gene Regulation Network of JAK2 in Skin Cutaneous Melanoma. Life Sci (2020) 253:117600. doi: 10.1016/j.lfs.2020.117600
28
KirbyJPriorFPetrickNHadjiskiLFarahaniKDrukkerKet al. Introduction to Special Issue on Datasets Hosted in The Cancer Imaging Archive (TCIA). Med Phys (2020) 47:6026–8. doi: 10.1002/mp.14595
29
CanzlerSHackermüllerJ. multiGSEA: A GSEA-Based Pathway Enrichment Analysis for Multi-Omics Data. BMC Bioinf (2020) 21:561. doi: 10.1186/s12859-020-03910-x
30
YoshiharaKShahmoradgoliMMartinezEVegesnaRKimHTorres-GarciaWet al. Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
31
ShenSWangGZhangRZhaoYYuHWeiYet al. Development and Validation of an Immune Gene-Set Based Prognostic Signature in Ovarian Cancer. EBioMedicine (2019) 40:318–26. doi: 10.1016/j.ebiom.2018.12.054
32
WangCLiangC. MSIpred: A Python Package for Tumor Microsatellite Instability Classification From Tumor Mutation Annotation Data Using a Support Vector Machine. Sci Rep (2018) 8:17546. doi: 10.1038/s41598-018-35682-z
33
CokelaerTChenEIorioFMendenMLightfootHSaez-RodriguezJet al. GDSC Tools for Mining Pharmacogenomic Interactions in Cancer. Bioinformatics (2087) 34:1226–8. doi: 10.1093/bioinformatics/btx744
34
RuBWongCNTongYZhongJYZhongSSWWuWCet al. TISIDB: An Integrated Repository Portal for Tumor-Immune System Interactions. Bioinformatics (2019) 35:4200–2. doi: 10.1093/bioinformatics/btz210
35
MerinoDMMcShaneLMFabrizioDFunariVChenSJWhiteJRet al. Establishing Guidelines to Harmonize Tumor Mutational Burden (TMB): In Silico Assessment of Variation in TMB Quantification Across Diagnostic Platforms: Phase I of the Friends of Cancer Research TMB Harmonization Project. J Immunother Cancer (2020) 8(1):e000147. doi: 10.1136/jitc-2019-000147
36
ZhouMZhangZZhaoHBaoSSunJ. A Novel lncRNA-Focus Expression Signature for Survival Prediction in Endometrial Carcinoma. BMC Cancer (2018) 18:39. doi: 10.1186/s12885-017-3983-0
37
BrooksRAFlemingGFLastraRRLeeNKMoroneyJWSonCHet al. Current Recommendations and Recent Progress in Endometrial Cancer. CA Cancer J Clin (2019) 69:258–79. doi: 10.3322/caac.21561
38
TravaglinoARaffoneASacconeGDe LucaCMolloAMascoloMet al. Immunohistochemical Nuclear Expression of Beta-Catenin as a Surrogate of CTNNB1 Exon 3 Mutation in Endometrial Cancer. Am J Clin Pathol (2019) 151:529–38. doi: 10.1093/ajcp/aqy178
39
RaffoneATravaglinoAMascoloMCarboneLGuidaMInsabatoLet al. TCGA Molecular Groups of Endometrial Cancer: Pooled Data About Prognosis. Gynecol Oncol (2019) 155:374–83. doi: 10.1016/j.ygyno.2019.08.019
40
ColomboNCreutzbergCAmantFBosseTGonzalez-MartinALedermannJet al. ESMO-ESGO-ESTRO Consensus Conference on Endometrial Cancer: Diagnosis, Treatment and Follow-Up. Radiother Oncol (2015) 117:559–81. doi: 10.1016/j.radonc.2015.11.013
41
MoricePLearyACreutzbergCAbu-RustumNDaraiE. Endometrial Cancer. Lancet (2016) 387:1094–108. doi: 10.1016/S0140-6736(15)00130-0
42
LeeYCLheureuxSOzaAM. Treatment Strategies for Endometrial Cancer: Current Practice and Perspective. Curr Opin Obs Gynecol (2017) 29:47–58. doi: 10.1097/GCO.0000000000000338
43
TewariKSBurgerRAEnserroDNorquistBMSwisherEMBradyMFet al. Final Overall Survival of a Randomized Trial of Bevacizumab for Primary Treatment of Ovarian Cancer. J Clin Oncol (2019) 37:2317–28. doi: 10.1200/JCO.19.01009
44
OzaAMPignataSPovedaAMcCormackMClampASchwartzBet al. Randomized Phase II Trial of Ridaforolimus in Advanced Endometrial Carcinoma. J Clin Oncol (2015) 33:3576–82. doi: 10.1200/JCO.2014.58.8871
45
HugoWZaretskyJMSunLSongCMorenoBHHu-LieskovanSet al. Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma. Cell (2017) 168:542. doi: 10.1016/j.cell.2017.01.010
46
KalbasiARibasA. Tumour-Intrinsic Resistance to Immune Checkpoint Blockade. Nat Rev Immunol (2020) 20:25–39. doi: 10.1038/s41577-019-0218-4
47
KrugerSIlmerMKoboldSCadilhaBLEndresSOrmannsSet al. Advances in Cancer Immunotherapy 2019 - Latest Trends. J Exp Clin Cancer Res (2019) 38:268. doi: 10.1186/s13046-019-1266-0
48
MartinsFSofiyaLSykiotisGPLamineFMaillardMFragaMet al. Adverse Effects of Immune-Checkpoint Inhibitors: Epidemiology, Management and Surveillance. Nat Rev Clin Oncol (2019) 16:563–80. doi: 10.1038/s41571-019-0218-0
49
WeiSCDuffyCRAllisonJP. Fundamental Mechanisms of Immune Checkpoint Blockade Therapy. Cancer Discov (2018) 8:1069–86. doi: 10.1158/2159-8290.CD-18-0367
50
CastroflorioEden HoedJSvistunovaDFinelliMJCebrian-SerranoACorrochanoSet al. The Ncoa7 Locus Regulates V-ATPase Formation and Function, Neurodevelopment and Behaviour. Cell Mol Life Sci (2021) 78:3503–24. doi: 10.1007/s00018-020-03721-6
51
HarmsFLParthasarathyPZorndtDAlawiMFuchsSHallidayBJet al. Biallelic Loss-of-Function Variants in TBC1D2B Cause a Neurodevelopmental Disorder With Seizures and Gingival Overgrowth. Hum Mutat (2020) 41:1645–61. doi: 10.1002/humu.24071
52
TangLPengCZhuSSZhouZLiuHChengQet al. Tre2-Bub2-Cdc16 Family Proteins Based Nomogram Serve as a Promising Prognosis Predicting Model for Melanoma. Front Oncol (2020) 10:579625. doi: 10.3389/fonc.2020.579625
53
HuZYangDTangYZhangXWeiZFuHet al. Five-Long non-Coding RNA Risk Score System for the Effective Prediction of Gastric Cancer Patient Survival. Oncol Lett (2019) 17:4474–86. doi: 10.3892/ol.2019.10124
54
WhitfieldHWhiteGSpriggCRileyAMPotterBVLHemmingsAMet al. An ATP-Responsive Metabolic Cassette Comprised of Inositol Tris/Tetrakisphosphate Kinase 1 (ITPK1) and Inositol Pentakisphosphate 2-Kinase (IPK1) Buffers Diphosphosphoinositol Phosphate Levels. Biochem J (2020) 477:2621–38. doi: 10.1042/BCJ20200423
55
KawaiSAriyasuHUrakiSTakeshimaKMoritaSInabaHet al. Imbalanced Expression of IGF2 and PCSK4 Is Associated With Overproduction of Big IGF2 in SFT With NICTH: A Pilot Study. J Clin Endocrinol Metab (2018) 103:2728–34. doi: 10.1210/jc.2018-00593
56
WilliamsHCSan MartinAAdamoCMSeidel-RogolBPounkovaLDatlaSRet al. Role of Coronin 1B in PDGF-Induced Migration of Vascular Smooth Muscle Cells. Circ Res (2012) 111:56–65. doi: 10.1161/CIRCRESAHA.111.255745
57
KoKKPowellMSHogarthPM. ZSWIM1: A Novel Biomarker in T Helper Cell Differentiation. Immunol Lett (2014) 160:133–8. doi: 10.1016/j.imlet.2014.01.016
58
NandagopalNRouxPP. Regulation of Global and Specific mRNA Translation by the mTOR Signaling Pathway. Transl (2015) 3:e983402. doi: 10.4161/21690731.2014.983402
59
RyuCSBaeJKimIJKimJOhSHKimOJet al. MPG and NPRL3 Polymorphisms are Associated With Ischemic Stroke Susceptibility and Post-Stroke Mortality. Diagn (Basel) (2020) 10. doi: 10.3390/diagnostics10110947
60
ZhouLDongSDengYYangPZhengYYaoLet al. GOLGA7 Rs11337, a Polymorphism at the MicroRNA Binding Site, Is Associated With Glioma Prognosis. Mol Ther Nucleic Acids (2019) 18:56–65. doi: 10.1016/j.omtn.2019.08.006
61
DongYZhangTLiXYuFGuoY. Comprehensive Analysis of Coexpressed Long Noncoding RNAs and Genes in Breast Cancer. J Obs Gynaecol Res (2019) 45:428–37. doi: 10.1111/jog.13840
62
MiaoXJiaLZhouHSongXZhouMXuJet al. miR-4299 Mediates the Invasive Properties and Tumorigenicity of Human Follicular Thyroid Carcinoma by Targeting ST6GALNAC4. IUBMB Life (2016) 68:136–44. doi: 10.1002/iub.1467
63
KoslaKOrzechowskaMJedroszkaDBarylaIBednarekAKPluciennikE. A Novel Set of WNT Pathway Effectors as a Predictive Marker of Uterine Corpus Endometrial Carcinoma-Study Based on Weighted Co-Expression Matrices. Front Oncol (2019) 9:360. doi: 10.3389/fonc.2019.00360
64
ZhangZLiLLiMWangX. The SARS-CoV-2 Host Cell Receptor ACE2 Correlates Positively With Immunotherapy Response and Is a Potential Protective Factor for Cancer Progression. Comput Struct Biotechnol J (2020) 18:2438–44. doi: 10.1016/j.csbj.2020.08.024
65
KolendaTGuglasKKopczynskaMSobocinskaJTeresiakABlizniakRet al. Good or Not Good: Role of miR-18a in Cancer Biology. Rep Pr Oncol Radiother (2020) 25:808–19. doi: 10.1016/j.rpor.2020.07.006
66
YangJLiHHuSZhouY. ACE2 Correlated With Immune Infiltration Serves as a Prognostic Biomarker in Endometrial Carcinoma and Renal Papillary Cell Carcinoma: Implication for COVID-19. Aging (Albany NY) (2020) 12:6518–35. doi: 10.18632/aging.103100
67
PadillaJJenkinsNTThornePKMartinJSRectorRSDavisJWet al. Identification of Genes Whose Expression Is Altered by Obesity Throughout the Arterial Tree. Physiol Genomics (2014) 46:821–32. doi: 10.1152/physiolgenomics.00091.2014
68
BianJXuYWuFPanQLiuY. Identification of a Five-Gene Signature for Predicting the Progression and Prognosis of Stage I Endometrial Carcinoma. Oncol Lett (2020) 20:2396–410. doi: 10.3892/ol.2020.11798
69
ZouRZhengMTanMXuHLuanNZhuL. Decreased PTGDS Expression Predicting Poor Survival of Endometrial Cancer by Integrating Weighted Gene Co-Expression Network Analysis and Immunohistochemical Validation. Cancer Manag Res (2020) 12:5057–75. doi: 10.2147/CMAR.S255753
70
ZhouHZouXLiHLiTChenLChengX. Decreased Secretoglobin Family 2A Member 1expression Is Associated With Poor Outcomes in Endometrial Cancer. Oncol Lett (2020) 20:24. doi: 10.3892/ol.2020.11885
Summary
Keywords
uterine corpus endometrial carcinoma, regulatory T cells, immunotherapy, prognostic markers, TCGA
Citation
Liu J, Geng R, Yang S, Shao F, Zhong Z, Yang M, Ni S, Cai L and Bai J (2021) Development and Clinical Validation of Novel 8-Gene Prognostic Signature Associated With the Proportion of Regulatory T Cells by Weighted Gene Co-Expression Network Analysis in Uterine Corpus Endometrial Carcinoma. Front. Immunol. 12:788431. doi: 10.3389/fimmu.2021.788431
Received
02 October 2021
Accepted
22 November 2021
Published
14 December 2021
Volume
12 - 2021
Edited by
Dominik Wolf, Innsbruck Medical University, Austria
Reviewed by
Yifei Liu, Affiliated Hospital of Nantong University, China; Nupur Mukherjee, National Institute for Research in Reproductive Health (ICMR), India
Updates
Copyright
© 2021 Liu, Geng, Yang, Shao, Zhong, Yang, Ni, Cai and Bai.
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: Jianling Bai, baijianling@njmu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Cancer Immunity and Immunotherapy, a section of the journal Frontiers in Immunology
Disclaimer
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.