A Novel Defined Pyroptosis-Related Gene Signature for Predicting the Prognosis of Endometrial Cancer

Liu, Shuguang; Zeng, Chao; Lv, Huaisheng; Zhang, Yan; Xiong, Hong; Tang, Hongping

doi:https://doi.org/10.1155/2022/7570494

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 7570494 | https://doi.org/10.1155/2022/7570494

A Novel Defined Pyroptosis-Related Gene Signature for Predicting the Prognosis of Endometrial Cancer

Shuguang Liu,¹Chao Zeng,¹Huaisheng Lv,²Yan Zhang,^3,4Hong Xiong,⁵and Hongping Tang²

Academic Editor: Mallikarjuna Korivi

Received28 Mar 2022

Revised15 Oct 2022

Accepted24 Nov 2022

Published16 Dec 2022

Abstract

Endometrial carcinoma (EC) is the second major female genital malignancy. Genetic signatures may be an improved choice to predict the prognosis of EC patients. The relationship between pyroptosis and tumours has attracted much attention in recent years. Here, we constructed a new pyroptosis-related gene (PRG) signature for predicting the prognosis of EC. In this study, gene data and clinical information of EC patients were obtained from The Cancer Genome Atlas (TCGA). Following the identification of PRGs correlated with EC prognosis, we further investigate the bioinformatics functions of these PRGs by univariate Cox regression analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Then, we used the least absolute contraction and selection operator (LASSO) regression and multiple Cox regression analysis to construct a new PRG signature that contains seven PRGs (NFKB1, EEF2K, CTSV, MDM2, GZMB, PANX1, and PTEN) and performed the Kaplan-Meier (K-M) analysis, receiver operating characteristic curve (ROC) analysis, and principal component analysis (PCA) to evaluate the prognostic value of our novel PRG signature. Finally, we assessed the correlations between pyroptosis and immune cells/checkpoints through the CIBERSORT tool and single-sample gene set enrichment analysis (ssGSEA). The result suggested that our signature was powerful in predicting EC prognosis and may play a part in assessing response to immunotherapy in EC patients. In conclusion, our study established a novel PRG signature for EC, which can be used as an effective prognostic marker in clinical practice in the future.

1. Introduction

Endometrial cancer (EC) severely threatens women’s health, ranking second in female genital malignancies. The incidence rate of endometrial cancer has been increasing significantly in recent years. According to the statistics, there were 417,000 new cases and 97,000 deaths in 2020 worldwide [1]. Although the therapeutic effect of EC is relatively good, a large number of patients still suffer from tumour recurrence, and approximately 18% of patients eventually die [2, 3]. Molecular classification and emerging targeted therapies can provide a basis for selecting clinical treatment strategies in EC [4, 5]. Therefore, exploring biomarkers to predict the prognosis is very important for the clinical diagnosis and treatment of EC.

Pyroptosis is a novel characterized form of programmed cell death (PCD) discovered after apoptosis and necrosis. Pyroptosis was first found in macrophages after pathogen infection and then confirmed to play a critical role in both inflammatory and immune defenses [6–9]. Recently, the role of pyroptosis in tumourigenesis, development, and treatment of tumours has attracted more and more attention. It was indicated that pyroptosis could promote inflammatory cell death of cancer and inhibit the proliferation and migration of tumour cells, thus affecting the progress and prognosis of tumour [10, 11]. For example, the expression of pyroptotic inflammasome GSDMD is decreased in gastric cancer (GC) cells, which promotes the proliferation of GC cells [10]. Recent studies have shown that PRG signatures can be used as effective prognostic markers in ovarian cancer [12], lung cancer [13], bladder cancer [14], and squamous cell cancer in the head and neck [15]. However, the prognostic value of PRGs in EC has not yet been elucidated.

In this study, we performed a systematic study to construct a prognostic multigene signature with PRGs and validated it in the TCGA cohort. This novel seven-PRG signature provides ideas for predicting the prognosis of EC and shares intriguing insights for more in-depth explorative research upon the mechanism of tumour immunity. Also, it may provide a basis for developing targeted anticancer therapy of EC.

2. Materials and Methods

2.1. Datasets and Preprocessing

The profiling gene data and corresponding clinical data of patients from 35 normal and 515 tumour tissues were obtained from TCGA (https://tcga-data.nci.nih.gov/tcga/) on November 12, 2021. Only samples from the primary tumour and RNA-seq data from fresh frozen samples were chosen in this study, and criteria were excluded as follows: data from FFPE tissue samples, duplicate RNA-seq data in the same patient sample, and cases without survival information. Furthermore, to develop a prognostic model, EC samples from TCGA database (after removing samples with patient’s survival days) were randomly divided into a training set () and a testing set (), and the ratio was set at 3 : 1.

2.2. Identification of PRGs in EC

A total of 133 PRGs were extracted from the GeneCards database (https://www.genecards.org/) to identify differentially expressed genes (DEGs, listed in Supplementary Table S1). The raw read counts were normalized, and DEGs related to pyroptosis were identified using the R package “DEseq2” [16]. The protein-protein interaction (PPI) network for the DEGs mentioned above was constructed by search tool for the retrieval of interacting genes (STRING, https://string-db.org/) [17].

2.3. Functional Enrichment Analysis

The Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/) pathway enrichment analyses and Gene Ontology (GO, http://www.geneontology.org) annotation were used to explore the biological functions of DEGs associated with pyroptosis through “org.Hs.eg.db,” “clusterProfiler,” “ggplot2,” and “enrichplot” packages in R. Adjusted value < 0.05 was established as the statistical significance threshold.

2.4. Construction of Prognostic PRG Signature

Univariate Cox regression analysis was performed first to assess the correlation between PRGs and the survival status of EC patients in the training set. Then, to reduce the dimensionality of the candidate genes according to the best penalty factor (λ), the LASSO regression analysis was applied. After multivariate Cox regression analysis was performed for these selected DEGs and their coefficients, a risk signature was constructed using the Akaike information criterion for the stepwise backward/forward model. The risk score (RS) for each patient was calculated using the risk formula: (expression of signature gene coefficient ).

2.5. Evaluation and Validation

According to the median risk score of our PRG signature, EC patients were divided into low-risk and high-risk groups. Survival curves were generated through the Kaplan-Meier method for the survival difference. Receiver operation characteristic (ROC) curves and area under the curve (AUC) values corresponding to 1, 2, and 3 years were used to evaluate the predictive efficiency of the signature. Principal component analysis (PCA) was adopted via the “prcomp” function in the “stats” R package. Correlation analysis was performed to verify the connections between the clinicopathological characteristics of the patients and the RS. The independent prediction ability was demonstrated by univariate and multivariate Cox regression analyses. In addition, stratified analysis was used to examine the precision of prognostic prediction based on different clinicopathological groups.

2.6. Construction and Evaluation by Nomogram

Based on the clinicopathological features and the RS, a nomogram model was then established to predict 1-, 2-, and 3-year mortality with the rms package of R. The calibration curve with the foreign package is used to determine the predictive performance of the nomogram.

2.7. Tumour Immunity Analysis

The CIBERSORT algorithm and the Spearman correlation were used to assess different distributions of 22 types of tumour-infiltrating immune cells (TIICs) with the variation of the risk score. The correlation between the signature-based risk score and the expression of the immune checkpoint genes (programmed cell death protein 1 (PD1), PD ligand 1 (PDL1), and cytotoxic T-lymphocyte-associated protein 4 (CTLA4)) was assessed by Pearson’s test.

2.8. Statistical Analysis

All statistical analyses were performed with the R (version 4.1.2, http://www.r-project.org) and R Bioconductor packages. The Kaplan-Meier curve with a log-rank test from the survival package and Cox proportional risk regression models were used for survival analysis. ROC analysis was used to detect the sensitivity and specificity of the gene signature risk score to predict survival. The AUC values can be used as an index of prognostic precision. In all analyses, value < 0.05 stands for a statistically significant difference.

3. Results

3.1. Defining the Expression of PRGs in EC

We first explored the expression of 133 PRG expressions from 35 normal and 515 tumour tissues in TCGA, using the DESeq2 package in R. A total of 97 PRGs were identified (all ), including 64 upregulated genes and 33 downregulated genes. Details of these PRGs were listed in Table S2. A heatmap of these genes is visualized in Figure 1(a). A PPI analysis with the minimum required interaction score of 0.9 (the highest confidence) was constructed to explore the interactions of the PRGs, which revealed that TP53, NLRP3, CASP8, NFKB1, AKT1, BIRC2, TNF, CASP3, JUN, and PYCARD (all node ) were the core genes (Figure 1(b)). The correlation network containing PRGs with is shown in Figure 1(c).

(a)

(b)

(c)

3.2. Functional Enrichment of the PRGs in EC

Function annotation analyses of the 97 PRGs were performed. GO enrichment (listed in Supplementary Table S3) suggested that these genes were mainly involved in “pyroptosis,” “interleukin-1 beta production,” “regulation of interleukin-1 beta production,” “response to virus,” and “cellular response to biotic stimulus” (Figure 2(a)). KEGG pathway enrichment analysis (listed in Supplementary Table S4) suggested that these PRGs were mainly involved in “NOD-like receptor signaling pathway,” “Salmonella infection,” “shigellosis,” “lipid and atherosclerosis,” and “apoptosis” (Figure 2(b)).

(a)

(b)

3.3. Identification and Validation of the Prognostic PRG Signature

Univariate Cox regression analysis was first performed to identify the prognostic-associated PRGs. Nineteen genes (MDM2, ALK, CRTAC1, IRF2, UBE2D2, GZMB, NFKB1, PTEN, EEF2K, CTSV, GBP5, ANXA2, IL32, SERPINB1, PRDM1, CASP3, IL13RA2, TP53, and PANX1) that meet the criteria are obtained for further analysis (Figure 3(a)). Then, by performing a least absolute shrinkage and selection operator (LASSO) (Figure 3(b)), fourteen PRGs were sought out according to the optimal value (Figure 3(c)). Finally, multivariate Cox regression analysis revealed that seven PRGs (NFKB1, CTSV, PANX1, PTEN, MDM2, GZMB, and EEF2K) were independent prognostic factors for patients with EC (Figure 3(d)).

(a)

(b)

(c)

(d)

As calculated in the following, the .

Based on the median risk score, 387 EC patients in the training set were divided into low-risk () and high-risk groups (). The distribution patterns of risk scores and survival status are illustrated in Figures 4(a) and 4(b). A heatmap exhibited seven-gene expression between the high- and low-risk groups (Figure 4(c)). PCA showed that patients with different risks were satisfactorily divided into two clusters (Figure 4(d)). The Kaplan-Meier curve showed that patients in the high-risk group had more deaths and shorter survival times than those in the low-risk group (Figure 4(e)). The areas under the ROC curve for predicting 1-, 3-, and 5-year survival rates were 0.732, 0.763, and 0.793, respectively (Figure 4(f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4

Construction of the risk signature in the training cohort. (a) Risk score map of patients in the high- and low-risk groups. (b) Survival status map of patients in the death (red spheres) and alive (blue spheres) groups. (c) Heatmap of gene expression in the high- and low-risk groups. (d) PCA plot showing different clustering between EC (red spheres) and normal endometrial tissues (blue spheres). (e) The Kaplan-Meier survival curves of patients in the high- and low-risk groups. (f) ROC curve evaluating the predictive performance of the risk signature.

To validate the efficiency of the signature, based on the same median risk score of the training set, 128 patients in the testing set were divided into low-risk () and high-risk groups (). The distribution patterns of risk scores and survival status are illustrated in Figures 5(a) and 5(b). A heatmap exhibited seven-gene expression between the high- and low-risk groups (Figure 5(c)). PCA showed that patients with different risks were satisfactorily divided into two clusters (Figure 5(d)). The Kaplan-Meier curve showed that patients in the high-risk group had more deaths and shorter survival times (Figure 5(e)). The AUCs for predicting 1-, 3-, and 5-year survival rates were 0.92, 0.789, and 0.739, respectively (Figure 5(f)). Meanwhile, we performed the same analyses in the entire set. The distribution patterns of risk scores and survival status are illustrated in Figures 6(a) and 6(b), and a heatmap exhibiting seven-gene expression between the high- and low-risk groups is shown in Figure 6(c). PCA showed that patients with different risks were satisfactorily divided into two clusters (Figure 6(d)). The Kaplan-Meier curve showed that patients in the high-risk group had more deaths and shorter survival times (Figure 6(e)). The AUCs for predicting 1-, 3-, and 5-year survival rates were 0.711, 0.729, and 0.754, respectively. (Figure 6(f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5

Construction of the risk signature in the test cohort. (a) Risk score map of patients in the high- and low-risk groups. (b) Survival status map of patients in the death (red spheres) and alive (blue spheres) groups. (c) Heatmap of gene expression in the high- and low-risk groups. (d) PCA plot showing different clustering between EC (red spheres) and normal endometrial tissues (blue spheres). (e) The Kaplan-Meier survival curves of patients in the high- and low-risk groups. (f) ROC curve evaluating the predictive performance of the risk signature.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6

Construction of the risk signature in the entire cohort. (a) Risk score map of patients in the high- and low-risk groups. (b) Survival status map of patients in the death (red spheres) and alive (blue spheres) groups. (c) Heatmap of gene expression in the high- and low-risk groups. (d) PCA plot showing different clustering between EC (red spheres) and normal endometrial tissues (blue spheres). (e) The Kaplan-Meier survival curves of patients in the high- and low-risk groups. (f) ROC curve evaluating the predictive performance of the risk signature.

3.4. Correlation of the PRG Prognostic Signature with Clinical Features

The correlation analysis showed significant relationship between the risk scores and clinicopathological features such as age (Figure 7(a), ), histology (Figure 7(b), ), and stage (Figure 7(c), ), while there was no significant relationship between the risk scores and grade (Figure 7(d), all ). These results suggested that our PRGs may be associated with EC progress.

(a)

(b)

(c)

(d)

The stratified analysis showed that high risk was an independent index of poorer survival outcome in comparison to low risk based on different subgroups of clinicopathological features, such as (, Figure 8(a)) and >60 (, Figure 8(b)), endometrioid (, Figure 8(c)) and serous (, Figure 8(d)), stages I-II (, Figure 8(e)) and stages III-IV (, Figure 8(f)), grade 1 (, Figure 8(g)), grade 2 (, Figure 8(h)), and grade 3 (, Figure 8(i)). The result suggested that our PRG signature was powerful to predict the prognosis in different subgroups of EC patients.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.5. Independent Predictive Value of PRG Risk Signature

Univariate and multivariate analyses were performed to verify the PRG risk signature as an independent prognostic factor in EC patients. Univariate analysis (listed in Supplementary Table S5) revealed that PRG risk score, as well as age, histology, and stage, was significantly correlated with the patient’s survival time (all , Figure 9(a)). Multivariate analysis (listed in Supplementary Table S6) demonstrated that PRG risk score, age, and stage predict the prognosis independently in patients with EC (all , Figure 9(b)). A heatmap for the connection between the predictive gene signature and clinicopathological features is shown in Figure 9(c). ROC curve analysis showed that the AUC value of the risk score was 0.716, which was significantly higher than that of patients’ age () and patients’ weight (), histology (), grade (), and stage () (Figure 9(d)). These results suggest that this PRG signature was an effectively independent predictor of the prognosis of EC patients.

(a)

(b)

(c)

(d)

3.6. Nomogram Building and Validation

To further estimate EC patients’ survival probability, a comprehensive prognostic nomogram was built based on six parameters, including age, weight, histology, grade, stage, and risk score (Figure 10(a)). The calibration plots showed excellent consistency between the nomogram predictions and actual observations in terms of the 1-year (Figure 10(b)), 3-year (Figure 10(c)), and 5-year (Figure 10(d)) survival rates in the entire TCGA cohort (listed in Supplementary Table S7). It was also confirmed the efficiency of prognosis prediction of our PRG risk signature.

(a)

(b)

(c)

(d)

3.7. Immunological Activity between Different Risk Subgroups

CIBERSORT was implemented to evaluate the association between the risk score and the 22 types of tumour-infiltrating immune cells (TIICs) in EC. We observed significantly higher proportions of B cell naïve, T cell CD4 memory resting, NK cells activated, and dendritic cells activated and lower proportions of plasma cells, CD8+ T cells, T cell regulatory (Tregs), NK cell resting, and dendritic cell resting in the high-risk group (Figure 11(a)). Then, the enrichment scores of 16 types of immune cells and the activity enrichment analysis of 13 immune-related pathways were compared between TCGA dataset using ssGSEA. Figure 11(b) shows that there are significantly higher proportions of immune cells, such as DCs and neutrophils in the high-risk group. Moreover, four immune pathways, such as APC_co_inhibition, checkpoint, HLA, and T_cell_coinhibition, were also significantly higher in the high-risk group (Figure 11(c)). The results indicated that the infiltration of these immune cells and pathways might significantly influence the prognosis of EC patients.

(a)

(b)

(c)

The expression levels of immune checkpoint genes and correlation analysis between the high- and low-risk groups showed significant decrease of CTLA4, PD1, and PDL1 in the high-risk group: the expression levels of CTLA4 (, Figure 12(a)) and the correlation between risk score and CTLA4 (, , Figure 12(b)), the expression levels of PD1 (, Figure 12(c)) and the correlation between risk score and PD1 (, , Figure 12(d)), and the expression levels of PDL1 (, Figure 12(e)) and the correlation between risk score and PDL1 (, , Figure 12(f)). These data indicated that the PRG signature might play a part in assessing response to immunotherapy in EC patients.

(a)

(b)

(c)

(d)

(e)

(f)

4. Discussion

EC is one of the most common female genital malignancies and threatens women’s lives worldwide [18]. Because of the lack of effective screening methods, there is still a large part of EC that is difficult to diagnose at the early stage [19]. The current clinicopathological features, such as age, tumour size, and tumour histological type, are not sufficient to accurately predict the outcome of patients with EC. Finding new molecular mechanisms is very necessary for the clinical diagnosis and treatment of EC.

As a novel immunogenic form of inflammatory programmed cell death pathway [20], pyroptosis can participate in local inflammation and attract immune cell infiltration, which provides a good opportunity to alleviate the immunosuppression of the tumour microenvironment [21]. It has been proven that pyroptosis may play various roles in tumours and PRGs have the effect of inhibiting tumour growth and tumour immunity [10]. Ye et al. [12] showed a novel signature of seven PRGs (AIM2, PJVK, PLCG1, GSDMA, ELANE, CASP3, and CASP6) for predicting the prognosis of ovarian cancer (OC). Zhang et al. [22] identified a new signature featuring seven PRGs (BAK1, CHMP4B, NOD2, NLRP6, GSDMC, PLCG1, and SCAF11) in predicting the prognosis of hepatocellular carcinoma (HCC) patients. Chen et al. [14] constructed eight-PRG signature (AIM2, BAK1, GZMA, GZMB, IRF1, NOD2, TNF, and TP63) to predict the survival in bladder cancer (BC) patients. Qian et al. [15] performed a novel PRG signature (GSDME, IL1B, NLRP1, and NLRP6) for prognostic prediction of head and neck squamous cell carcinoma. In EC, Chen et al. [23] recently analyzed 33 PRGs between tumour samples and normal samples and obtained six pyroptosis-related prognostic DEGs (GPX4, GSDMD, GSDME, IL6, NOD2, and PYCARD).

In this study, using bioinformatics, we compared the expression of 133 currently known PRGs between EC and normal endometrial tissues in TCGA dataset and identified 19 DEGs associated with pyroptosis in EC. In order to further assess the prognostic value of these PRGs, we constructed 7-gene signature (NFKB1, CTSV, PANX1, PTEN, MDM2, GZMB, and EEF2K) via Cox univariate analysis and LASSO Cox regression analysis and then subsequently performed GO and KEGG to validate our PRG risk signature.

The nuclear factor-κB (NF-κB) signaling pathway plays an important role in inflammation and immune response [24]. As a subunit of NF-κB, NFKB1 has been shown to be a pathway-specific suppressor of inflammation, aging, and tumours such as hematological malignancies [25]. Upregulation of NFKB1 can promote the invasiveness of breast cancer cells in vitro [26]. Cathepsin V (CTSV) is a member of the cathepsin family, which is highly expressed in activated macrophages and is involved in inflammatory diseases like myasthenia gravis [27, 28]. Studies showed that high expression of CTSV was associated with poor prognosis of breast cancer [29], and ectopic expression of CTSV increased the number of migrated and invaded colorectal cancer cells in vivo [30].

Pannexin 1 (PANX1) is a critical ATP-releasing channel and widely participates in the regulation of the tumour immune microenvironment by infiltrating multiple immune cells, such as cancer-associated fibroblast cells, macrophage, lymphocytes, and neutrophil cells [31]. High expression of PANX1 was significantly related to the poor outcome of multiple cancers, especially in pancreatic adenocarcinoma (PAAD) [32]. In addition, PANX1 overexpression is associated with the EMT transformation of breast cancer cells and poorer clinical outcomes in breast cancer patients [33]. As a tumour suppressor gene, phosphatase and tensin homolog (PTEN) protein plays an important role in PI3K/AKT pathway [34]. The inactivation of PTEN may be considered a key factor for early endometrial tumourigenesis [35]. Similar to our study, Kanamori et al. found that PTEN was a significant prognostic indicator of EC [36]. As a critical negative regulator of the tumour suppressor p53, mouse double minute 2 (MDM2) can bind mutant and wild-type p53 protein and plays a key role in controlling its function [37]. MDM2 is abnormally upregulated in several types of tumours, especially those of mesenchymal origin [38]. Ambros et al. first reported that there was a selective correlation of expression between MDM2 and p53 expression in EC [39]. Buchynska et al. further proved that high expression of p53 with low expression of MDM2 might be the characteristic features of low-differentiated EC [40]. In our study, as a prognostic PRG, the expression level of MDM2 in EC is also lower than in normal endometrial tissues. Granzyme B (GZMB) is a component of cytolytic granules within cytotoxic T lymphocytes (CTLs) and natural killer (NK) cells, which are involved in several pathologies, including the formation of the tumour microenvironment [41]. GZMB expression was associated with worse clinical outcomes in non-small-cell lung cancer [42] and with early signs of metastasis in colorectal cancer [43]. Eukaryotic elongation factor 2 kinase (EEF2K) is a new member of atypical protein kinase and acts as a negative regulator of cell growth [44]. It was shown that EEF2K is activated and overexpressed in many cancers. EEF2K can positively help to improve breast cancer cells’ survival ability under nutrient deprivation or insufficient growth factor [45]. EEF2K is also a clinical indicator of metastasis and prognosis of stomach adenocarcinoma and could serve as a potential therapeutic target [46].

In our study, this PRG prognostic signature’s association with EC progression can be preferably illustrated through the results of correlation analysis and stratified analysis, and this finding exerts satisfactory predicting power in terms of the survival period of EC patients. The development of univariate and multivariate analyses facilitates our verification, which suggests that PRG’s risk can serve as an independent prognostic factor in EC patients.

Furthermore, in-depth exploring the mechanisms of this genes signature, CIBERSORT functional analysis of immune cell subsets revealed significantly higher proportions of B cells naïve, T cell CD4 memory resting, NK cells activated, and dendritic cells activated and lower proportions of plasma cells, CD8+ T cells, Tregs, NK cell resting, and dendritic cell resting in the high-risk group of EC. Then, the enrichment score analyses showed that there were significantly higher proportions of immune cells, such as DCs and neutrophils, and four immune pathways, such as APC_co_inhibition, checkpoint, HLA, and T_cell_coinhibition, were also observed to be significantly higher in the high-risk group (Figure 10(c)). These results indicated that reduced levels of antitumour immunity might lead to a poor prognosis. Additionally, our study showed a significant decrease in immune checkpoint genes of CTLA4, PD1, and PDL1 in the high-risk group. Current studies have shown that immunotherapy, especially immune checkpoint inhibitors (ICIS), may produce lasting therapeutic effects. The combined application of induced blepharoptosis and ICIS can enhance anticancer activity [47]. The PD-1 pathway is a critical pathway of immunosuppression in the tumour microenvironment. The generation of endogenous antitumour immunity to inhibit cancer development can be ideally achieved through the inhibition of PD-1 and PD-L1 [48]. Further research on the pyroptosis-related immune mechanism will be beneficial to guide clinical diagnosis and treatment of EC.

In summary, more and more studies had shown that PRGs are closely related to tumour prognosis, but as a complex gene network, there may be differences in key genes related to prognosis in different types of tumours and even in different research groups of the same type of tumours. Here, we developed a novel PRG signature using different platforms and validated this signature model. Our study suggested that using this seven-PRG signature is able to predict prognosis in patients with EC. It may provide new insight into the molecular mechanism and a classification tool in the clinical practice of EC.

Data Availability

The datasets and code generated and analyzed in this study are available from the corresponding author on request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Shuguang Liu, Chao Zeng, and Hongping Tang designed this study. Shuguang Liu, Chao Zeng, Huaisheng Lv, and Yan Zhang collected and analyzed the data. Shuguang Liu, Hong Xiong, and Hongping Tang wrote and revised the manuscript. All authors have read and approved the manuscript. Shuguang Liu and Chao Zeng contributed equally to this work.

Acknowledgments

This study was supported by grants from the Futian District Health Public Welfare Scientific Research Project of Shenzhen (Grant No. FTWS2019014). We acknowledge TCGA database for providing the platform for uploading meaningful datasets.

Supplementary Materials

Supplementary 1. Table S1: 133 pyroptosis-related genes retrieved from the GeneCards.

Supplementary 2. Table S2: 97 pyroptosis-related DEGs in EC.

Supplementary 3. Table S3: GO enrichment analysis.

Supplementary 4. Table S4: KEGG pathway analysis.

Supplementary 5. Table S5: univariate Cox regression analysis.

Supplementary 6. Table S6: multivariate Cox regression analysis.

Supplementary 7. Table S7: clinicopathological features as well as risk score.

References

H. Sung, J. Ferlay, R. L. Siegel et al., “Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 71, no. 3, pp. 209–249, 2021.
View at: Publisher Site | Google Scholar
G. Trojano, C. Olivieri, R. Tinelli, G. R. Damiani, A. Pellegrino, and E. Cicinelli, “Conservative treatment in early stage endometrial cancer: a review,” Acta Bio-Medica, vol. 90, no. 4, pp. 405–410, 2019.
View at: Publisher Site | Google Scholar
M. G. Del Carmen, D. M. Boruta, and J. O. Schorge, “Recurrent endometrial cancer,” Clinical Obstetrics and Gynecology, vol. 54, no. 2, pp. 266–277, 2011.
View at: Publisher Site | Google Scholar
M. E. Urick and D. W. Bell, “Clinical actionability of molecular targets in endometrial cancer,” Nature Reviews Cancer, vol. 19, no. 9, pp. 510–521, 2019.
View at: Publisher Site | Google Scholar
T. T. Yen, T. L. Wang, A. N. Fader, I. M. Shih, and S. Gaillard, “Molecular classification and emerging targeted therapy in endometrial cancer,” International Journal of Gynecological Pathology, vol. 39, no. 1, pp. 26–35, 2020.
View at: Publisher Site | Google Scholar
M. Wellington, K. Koselny, F. S. Sutterwala, and D. J. Krysan, “Candida albicans triggers NLRP3-mediated pyroptosis in macrophages,” Eukaryotic Cell, vol. 13, no. 2, pp. 329–340, 2014.
View at: Publisher Site | Google Scholar
I. Jorgensen and E. A. Miao, “Pyroptotic cell death defends against intracellular pathogens,” Immunological Reviews, vol. 265, no. 1, pp. 130–142, 2015.
View at: Publisher Site | Google Scholar
C. Y. Taabazuing, M. C. Okondo, and D. A. Bachovchin, “Pyroptosis and apoptosis pathways engage in bidirectional crosstalk in monocytes and macrophages,” Cell Chemical Biology, vol. 24, no. 4, pp. 507–514.e4, 2017.
View at: Publisher Site | Google Scholar
D. A. Erkes, W. Cai, I. M. Sanchez et al., “Mutant BRAF and MEK inhibitors regulate the tumor immune microenvironment via pyroptosis,” Cancer Discovery, vol. 10, no. 2, pp. 254–269, 2020.
View at: Publisher Site | Google Scholar
Y. Fang, S. Tian, Y. Pan et al., “Pyroptosis: a new frontier in cancer,” Biomedicine & Pharmacotherapy, vol. 121, article 109595, 2020.
View at: Publisher Site | Google Scholar
Y. F. Tan, M. Wang, Z. Y. Chen, L. Wang, and X. H. Liu, “Inhibition of BRD4 prevents proliferation and epithelial-mesenchymal transition in renal cell carcinoma via NLRP3 inflammasome-induced pyroptosis,” Cell Death & Disease, vol. 11, no. 4, p. 239, 2020.
View at: Publisher Site | Google Scholar
Y. Ye, Q. Dai, and H. Qi, “A novel defined pyroptosis-related gene signature for predicting the prognosis of ovarian cancer,” Cell Death Discovery, vol. 7, no. 1, p. 71, 2021.
View at: Publisher Site | Google Scholar
J. Gao, X. Qiu, G. Xi et al., “Downregulation of GSDMD attenuates tumor proliferation via the intrinsic mitochondrial apoptotic pathway and inhibition of EGFR/Akt signaling and predicts a good prognosis in non-small cell lung cancer,” Oncology Reports, vol. 40, no. 4, pp. 1971–1984, 2018.
View at: Publisher Site | Google Scholar
W. Chen, W. Zhang, T. Zhou, J. Cai, Z. Yu, and Z. Wu, “A newly defined pyroptosis-related gene signature for the prognosis of bladder cancer,” International Journal of General Medicine, vol. 12, no. 14, pp. 8109–8120, 2021.
View at: Publisher Site | Google Scholar
X. Qian, J. Tang, Y. Chu et al., “A novel pyroptosis-related gene signature for prognostic prediction of head and neck squamous cell carcinoma,” International Journal of General Medicine, vol. 14, no. 14, pp. 7669–7679, 2021.
View at: Publisher Site | Google Scholar
M. I. Love, W. Huber, and S. Anders, “Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2,” Genome Biology, vol. 15, no. 12, p. 550, 2014.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
J. N. McAlpine, S. M. Temkin, and H. J. Mackay, “Endometrial cancer: not your grandmother’s cancer,” Cancer, vol. 122, no. 18, pp. 2787–2798, 2016.
View at: Publisher Site | Google Scholar
S. Bendifallah, M. Ballester, and E. Darai, “Endometrial cancer: predictive models and clinical impact,” Bulletin du Cancer, vol. 104, no. 12, pp. 1022–1031, 2017.
View at: Publisher Site | Google Scholar
S. M. Man, R. Karki, and T. D. Kanneganti, “Molecular mechanisms and functions of pyroptosis, inflammatory caspases and inflammasomes in infectious diseases,” Immunological Reviews, vol. 277, no. 1, pp. 61–75, 2017.
View at: Publisher Site | Google Scholar
M. Wang, S. Jiang, Y. Zhang, P. Li, and K. Wang, “The multifaceted roles of pyroptotic cell death pathways in cancer,” Cancers, vol. 11, no. 9, p. 1313, 2019.
View at: Publisher Site | Google Scholar
S. Zhang, X. Li, X. Zhang, S. Zhang, C. Tang, and W. Kuang, “The pyroptosis-related gene signature predicts the prognosis of hepatocellular carcinoma,” Frontiers in Molecular Biosciences, vol. 8, 2021.
View at: Publisher Site | Google Scholar
Y. Chen, Y. Liao, Q. Du et al., “Roles of pyroptosis-related gene signature in prediction of endometrial cancer outcomes,” Frontiers in Medicine, vol. 9, 2022.
View at: Publisher Site | Google Scholar
Q. Li and I. M. Verma, “NF-κB regulation in the immune system,” Nature Reviews Immunology, vol. 2, no. 10, pp. 725–734, 2002.
View at: Publisher Site | Google Scholar
T. Cartwright, N. D. Perkins, and C. L. Wilson, “NFKB1: a suppressor of inflammation, ageing and cancer,” The FEBS Journal, vol. 283, no. 10, pp. 1812–1822, 2016.
View at: Publisher Site | Google Scholar
G. C. Kim, H. K. Kwon, C. G. Lee et al., “Upregulation of Ets1 expression by NFATc2 and NFKB1/RELA promotes breast cancer cell invasiveness,” Oncogene, vol. 7, no. 11, p. 91, 2018.
View at: Publisher Site | Google Scholar
Y. Yasuda, Z. Li, D. Greenbaum, M. Bogyo, E. Weber, and D. Brömme, “Cathepsin V, a novel and potent elastolytic activity expressed in activated macrophages,” The Journal of Biological Chemistry, vol. 279, no. 35, pp. 36761–36770, 2004.
View at: Publisher Site | Google Scholar
E. Tolosa, W. Li, Y. Yasuda et al., “Cathepsin V is involved in the degradation of invariant chain in human thymus and is overexpressed in myasthenia gravis,” The Journal of Clinical Investigation, vol. 112, no. 4, pp. 517–526, 2003.
View at: Publisher Site | Google Scholar
N. Sereesongsaeng, S. H. McDowell, J. F. Burrows, C. J. Scott, and R. E. Burden, “Cathepsin V suppresses GATA3 protein expression in luminal A breast cancer,” Breast Cancer Research, vol. 22, no. 1, p. 139, 2020.
View at: Publisher Site | Google Scholar
C. H. Wang, L. K. Wang, C. C. Wu et al., “Cathepsin V mediates the tazarotene-induced gene 1-induced reduction in invasion in colorectal cancer cells,” Cell Biochemistry and Biophysics, vol. 78, no. 4, pp. 483–494, 2020.
View at: Publisher Site | Google Scholar
S. Penuela, L. Harland, J. Simek, and D. W. Laird, “Pannexin channels and their links to human disease,” The Biochemical Journal, vol. 461, no. 3, pp. 371–381, 2014.
View at: Publisher Site | Google Scholar
L. Bao, K. Sun, and X. Zhang, “PANX1 is a potential prognostic biomarker associated with immune infiltration in pancreatic adenocarcinoma- a pan-cancer analysis,” Channels, vol. 15, no. 1, pp. 680–696, 2021.
View at: Publisher Site | Google Scholar
N. Jalaleddine, L. El-Hajjar, H. Dakik et al., “Pannexin 1 is associated with enhanced epithelial-to-mesenchymal transition in human patient breast cancer tissues and in breast cancer cell lines,” Cancers, vol. 11, p. 1967, 2019.
View at: Publisher Site | Google Scholar
T. L. Yuan and L. C. Cantley, “PI3K pathway alterations in cancer: variations on a theme,” Oncogene, vol. 27, no. 41, pp. 5497–5510, 2008.
View at: Publisher Site | Google Scholar
B. Konopka, Z. Paszko, A. Janiec-Jankowska, and M. Goluda, “Assessment of the quality and frequency of mutations occurrence in PTEN gene in endometrial carcinomas and hyperplasias,” Cancer Letters, vol. 178, no. 1, pp. 43–51, 2002.
View at: Publisher Site | Google Scholar
Y. Kanamori, J. Kigawa, H. Itamochi et al., “PTEN expression is associated with prognosis for patients with advanced endometrial carcinoma undergoing postoperative chemotherapy,” International Journal of Cancer, vol. 100, no. 6, pp. 686–689, 2002.
View at: Publisher Site | Google Scholar
J. D. Oliner, A. Y. Saiki, and S. Caenepeel, “The role of MDM2 amplification and overexpression in tumorigenesis,” Cold Spring Harbor Perspectives in Medicine, vol. 6, no. 6, 2016.
View at: Publisher Site | Google Scholar
M. Mendoza, G. Mandani, and J. Momand, “The MDM2 gene family,” Biomolecular Concepts, vol. 5, no. 1, pp. 9–19, 2014.
View at: Publisher Site | Google Scholar
R. A. Ambros, C. E. Sheehan, B. V. Kallakury et al., “MDM2 and p 53 protein expression in the histologic subtypes of endometrial carcinoma,” Modern Pathology, vol. 9, no. 12, pp. 1165–1169, 1996.
View at: Google Scholar
L. G. Buchynska, I. P. Nesina, and E. V. Kashuba, “Different trends of p 53, MDM2 and p 14 ARF expression patterns in endometrial adenocarcinomas versus hyperplasia,” Experimental Oncology, vol. 29, no. 4, pp. 287–294, 2007.
View at: Google Scholar
M. Konishi, S. S. Erdem, R. Weissleder, A. H. Lichtman, J. R. McCarthy, and P. Libby, “Imaging granzyme B activity assesses immune-mediated myocarditis,” Circulation Research, vol. 117, no. 6, pp. 502–512, 2015.
View at: Publisher Site | Google Scholar
D. P. Hurkmans, E. A. Basak, N. Schepers et al., “Granzyme B is correlated with clinical outcome after PD-1 blockade in patients with stage IV non-small-cell lung cancer,” Journal for Immunotherapy of Cancer, vol. 8, no. 1, article e000586, 2020.
View at: Publisher Site | Google Scholar
P. Salama, M. Phillips, C. Platell, and B. Iacopetta, “Low expression of granzyme B in colorectal cancer is associated with signs of early metastastic invasion,” Histopathology, vol. 59, no. 2, pp. 207–215, 2011.
View at: Publisher Site | Google Scholar
J. W. Kenney, C. E. Moore, X. Wang, and C. G. Proud, “Eukaryotic elongation factor 2 kinase, an unusual enzyme with multiple roles,” Advances in Biological Regulation, vol. 55, pp. 15–27, 2014.
View at: Publisher Site | Google Scholar
G. Leprivier, B. Rotblat, D. Khan, E. Jan, and P. H. Sorensen, “Stress-mediated translational control in cancer cells,” Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms, vol. 1849, no. 7, pp. 845–860, 2015.
View at: Publisher Site | Google Scholar
M. Jiang, L. Qi, K. Jin et al., “eEF2K as a novel metastatic and prognostic biomarker in gastric cancer patients,” Pathology-Research and Practice, vol. 225, article 153568, 2021.
View at: Publisher Site | Google Scholar
Q. Wang, Y. Wang, J. Ding et al., “A bioorthogonal system reveals antitumour immune function of pyroptosis,” Nature, vol. 579, no. 7799, pp. 421–426, 2020.
View at: Publisher Site | Google Scholar
A. D. Garg and P. Agostinis, “Cell death and immunity in cancer: from danger signals to mimicry of pathogen defense responses,” Immunological Reviews, vol. 280, no. 1, pp. 126–148, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Shuguang Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

366

Downloads

491

Citations