- 1Department of General Surgery, Xiangya Hospital, Central South University, Changsha, China
- 2Department of Pathology, Xiangya Hospital, Central South University, Changsha, China
- 3National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
- 4Department of Pharmacy, Xiangya Hospital, Central South University, Changsha, China
Background: Pyroptosis is a newly recognized form of cell death. Emerging evidence has suggested the crucial role of long non-coding RNAs (lncRNAs) in the tumorigenesis and progression of ovarian cancer (OC). However, there is still poor understanding of pyroptosis-related lncRNAs in OC.
Methods: The TCGA database was accessed for gene expression and clinical data of 377 patients with OC. Two cohorts for training and validation were established by random allocation. Correlation analysis and Cox regression analysis were performed to identify pyroptosis-related lncRNAs and construct a risk model.
Results: Six pyroptosis-related lncRNAs were included in the final signature with unfavorable survival data. Subsequent ROC curves showed promising predictive value of patient prognosis. Further multivariate regression analyses confirmed the signature as an independent risk factor in the training (HR: 2.242, 95% CI: 1.598–3.145) and validation (HR: 1.884, 95% CI: 1.204–2.95) cohorts. A signature-based nomogram was also established with a C-index of.684 (95% CI: 0.662–0.705). Involvement of the identified signature in multiple immune-related pathways was revealed by functional analysis. Moreover, the signature was also associated with higher expression of three immune checkpoints (PD-1, B7-H3, and VSIR), suggesting the potential of the signature as an indicator for OC immunotherapies.
Conclusion: This study suggests that the identified pyroptosis-related lncRNA signature and signature-based nomogram may serve as methods for risk stratification of OC. The signature is also associated with the tumor immune microenvironment, potentially providing an indicator for patient selection of immunotherapy in OC.
Background
Ovarian cancer (OC) is a lethal gynecologic malignancy among all gynecological tumors (1). Surgery is recognized as the first-line treatment for OC; however, 4 patients out of 5 are unable to receive surgery because of advanced disease. Although new therapies other than chemotherapy are developed or under-developing for advanced OC, the 5-year survival rate improves with a relatively slow speed (2). The recognition of high- and low-risk patient groups may help individualize treatments, thus improving patient prognosis. Moreover, although immunotherapies improve patient prognosis of several cancers (3), there is currently no significant breakthrough in the development of immunotherapy in treating OC during the decade. A promising indicator is also needed to screen patients with OC who may benefit from certain immunotherapies
Pyroptosis, an inflammatory form of cell death triggered by certain inflammasomes (4), has been found to be related to multiple human diseases such as OC (5). Many genes are identified to play an important role in the processes of pyroptosis, such as NLRP3, which consists of the NLRP3/caspase-1 signaling pathway (6, 7). In the meantime, long non-coding RNAs (lncRNAs), members of the non-coding RNA family, also participate in the development of a variety of cancers (8). It is confirmed by many researchers that lncRNAs may have a complex impact on the development of OC (9). However, the importance of lncRNAs, which are associated with pyroptosis-related genes, is poorly investigated in OC.
The aim of this study is to identify prognostic pyroptosis-related lncRNAs, therefore establishing a promising signature and a signature-based nomogram for risk stratification in OC. The association between the signature and immune microenvironment is also investigated, as well as its role as an immunotherapy indicator.
Materials and Methods
Data Acquisition
Transcriptome and clinical data of 377 patients with OC were retrieved from the OV project of the TCGA database (http://cancergenome.nih.gov/), while the data of 88 normal tissues were retrieved from GTEx. Patients without adequate clinical data were excluded from the analyses.
Identification of Pyroptosis-Related lncRNAs
A total of 33 pyroptosis-related genes were obtained from reports of Ye (Supplementary Table 1) (10). Pearson correlation test was performed to calculate correlations between lncRNAs and pyroptosis-related genes. A pyroptosis-related lncRNA was identified with Pearson correlation coefficient >0.3 and p < 0.001.
Construction and Validation of the Signature
The cohort was randomly divided into a training cohort and a validation cohort at a 2:1 ratio. The data from the training cohort were used for constructing the prognostic pyroptosis-related lncRNAs signature, while the other was for validation. Univariate Cox regression analysis was performed to identify prognostic pyroptosis-related lncRNAs. Subsequently, the least absolute shrinkage and selection operator (LASSO) Cox regression was used to construct the signature, presented as follows: risk score = expression of lncRNA1 × β1lncRNA1 + expression of lncRNA2 × β2lncRNA2 + …expression of lncRNAn × βnlncRNAn. The two cohorts were further divided into the low-risk and high-risk groups, respectively. Survival analyses and time-dependent ROC curves were performed to investigate the prognostic value. Moreover, multivariate Cox regression of available patient characteristics was performed to reconfirm the prognostic value of the signature. In addition, a nomogram was constructed to predict patient prognosis more precisely.
The mRNA-lncRNA Co-Expression Network
To better demonstrate the associations between pyroptosis-related genes and pyroptosis-related lncRNAs, a co-expression network was constructed. A Sankey diagram was used to illustrate mRNA-lncRNA relationships.
Gene Set Enrichment Analysis and Subsequent Functional Enrichment Analyses
Tumor tissues were divided into the low-risk and high-risk groups based on risk scores. Differentially expressed genes between the groups were identified by the “DEseq2” package with cut-off criteria of false discovery rate <0.05 and |log2foldchange| > 1. After that, the differentially expressed genes were uploaded for GSEA analysis (http://www.broadinstitute.org/gsea) (11). The protein-protein interaction (PPI) chart was realized with the STRING database (12).
Immunological Analysis
The abundance of tumor-infiltrating immune cells in OC tissues was investigated using the MCPcounter (13) and ssGSEA algorithms (14).
Cell Culture
The human ovarian epithelial cell IOSE80 and ovarian cancer cells A2780 and TOV112D were obtained from Center for Molecular Medicine, Xiangya Hospital, Central South University. All cell lines were cultured in Dulbecco's Modified Eagle's Medium (DMEM; Gibco, Grand Island, NY, USA) with 10% fetal bovine serum (FBS; Gibco, Grand Island, NY, USA) and 1% penicillin/streptomycin (Gibco, Grand Island, NY, USA). The cultures were placed in a sterile incubator maintained at 37 C with 5% CO2.
Quantitative PCR (qPCR)
Total RNA was extracted from disparate cell lines using a TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the protocol of the manufacturer and then converted to cDNA using a PrimeScriptTMRT reagent kit (6210; Takara, Dalian, China). A qPCR assay was performed with iTaqTM Universal SYBR Green Supermix (1725121; Bio-Rad, Hercules, CA, USA) to determine relative RNA levels. β-actin was used as an internal control for quantification of each gene. The sequences of gene primers are displayed in Supplementary Table 2. Relative expression levels of RNAs were determined using the 2−ΔΔCT method.
Statistical Analyses
Statistical analyses were performed with the R 3.3.0 platform. One-way analysis of variance (ANOVA) and Welch's ANOVA were performed as appropriate. Generally, p < 0.05 was considered statistically significant.
Results
The Training and Validation Cohorts
After excluding 2 patients without adequate clinical data (Figure 1A), the training cohort with 250 patients and the validation cohort with 125 patients were established by random allocation. Table 1 shows the patient and tumor characteristics of the cohorts. Except for neoplasm histologic grade, no difference was shown to have statistical significance.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g001.jpg)
Figure 1. Identification of prognostic pyroptosis-associated long noncoding ribonucleic acids (lncRNAs) in patients with ovarian cancer (OC). (A) Flow chart of constructing the pyroptosis-related lncRNAs signature. (B) Heatmap of 6 prognostic pyroptosis-related lncRNAs in OC tissues and normal ovarian tissues. (C) Univariate Cox regression of 6 prognostic pyroptosis-related lncRNAs.
Identification of Prognostic Pyroptosis-Related lncRNAs
A total of 4,222 pyroptosis-related lncRNAs were identified by Pearson correlation test, and subsequent univariate Cox regression and Kaplan-Meier curves were performed to determine the prognostic value of each lncRNA in the training cohort, suggesting 71 prognostic pyroptosis-related lncRNAs. After excluding 8 lncRNAs without significance in expression analyses, 63 candidate lncRNAs were included in the LASSO Cox regression analysis. Six pyroptosis-related lncRNAs (AC006001.2, LINC02585, AL136162.1, AC005041.3, AL023583.1, and LINC02881) were finally screened in the prognostic signature: risk score = (0.0624*AC006001.2 expression) + (−0.1014*LINC02585 expression) + (−0.1389*AL136162.1 expression) + (0.2398*AC005041.3 expression) + (−0.1734*AL023583.1 expression) + (2.1483*LINC02881 expression). The heatmap showed the expression level of the 6 lncRNAs in OC (Figure 1B), and the results of univariate regression analyses showed that the unfavorable prognostic value of AC006001.2, AC005041.3, LINC02881, and favorable prognostic value of LINC02585, AL136162.1, and AL023583.1 (Figure 1C).
Associations between the pyroptosis-related genes and lncRNAs were also investigated by co-expression analyses (Supplementary Table 3). Figure 2A shows the constructed network including 6 pyroptosis-related genes (CASP9, GSDMA, NLRP1, NOD1, PJVK, and PLCG1). It was noteworthy that AC005041.3 was associated with three pyroptosis-related genes, namely, CASP9, NLRP1, and PLCG1. Additionally, the Sankey diagram was also adopted to illustrate associations between genes and lncRNAs with their risk types (Figure 2B). Biological effects of 33 pyroptosis-related genes were also investigated, including GO (Figure 2C), KEGG (Figure 2D), and PPI chart (Figure 2E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g002.jpg)
Figure 2. The messenger RNA (mRNA)-lncRNA co-expression network. (A) mRNA-lncRNA co-expression network of the pyroptosis-related genes and selected pyroptosis-related lncRNAs. (B) Sankey diagram showing the connection degree between the pyroptosis-related lncRNAs and the pyroptosis-related genes. (C) Gene Ontology (GO) analysis of 33 pyroptosis-related genes. (D) Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of 33 pyroptosis-related genes. (E) PPI chart of 33 pyroptosis-related genes.
Validation of Prognostic Pyroptosis-Related lncRNA Signature
Next, the prognostic value of the signature was further examined in both the training and validation cohorts by Kaplan-Meier curves and time-dependent receiver operating characteristic (ROC) curves. After dividing the patients into the low-risk and high-risk groups (Table 2), worse survival data were shown among patients in the high-risk group in both cohorts (Figures 3A–D). Survival analyses further confirmed the unfavorable prognostic value of the signature in both cohorts (Figures 3E,F). In terms of ROC curves (Figures 3G,H), the area under curve (AUC) reached 0.681 at 1-year, 0.712 at 3-year, and 0.753 at 5-year in the training cohort, and 0.7 at 1-year, 0.573 at 3-year, and 0.621 at 5-year in the validation cohort. Survival analyses on the whole cohort also confirmed the unfavorable prognostic value of the signature (Figure 3I). Subsequently, the results of multivariate analyses confirmed the signature as an independent prognostic factor in both cohorts (Tables 3, 4). These results suggest a promising value of the signature in predicting prognosis of patients with OC. Additionally, the results of multivariate analysis of the risk score and the clinical characteristics in the whole cohort are shown in Table 5, while the clinical correlations between the clinical characteristics and risk score are shown in Table 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-t002.jpg)
Table 2. Associations between risk score and characteristics of patients with ovarian cancer in the training cohort.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g003.jpg)
Figure 3. Prognostic analysis of pyroptosis-related lncRNA signature in the training and validation cohorts. (A) Distribution of risk scores in the training cohort. (B) Distribution of risk scores in the validation cohort. (C) Distributions of overall survival status, overall survival, and risk score in the training cohort. (D) Distributions of overall survival status, overall survival, and risk score in the validation cohort. (E) Kaplan-Meier curves for the overall survival of patients in the high- and low-risk groups in the training cohort. (F) Kaplan-Meier curves for the overall survival of patients in the high- and low-risk groups in the validation cohort. (G) AUC of time-dependent ROC curves verified the prognostic accuracy of the risk score in the training cohort. (H) AUC of time-dependent ROC curves verified the prognostic accuracy of the risk score in the validation cohort. (I) Kaplan-Meier curves for the overall survival of patients in the high- and low-risk groups in the whole cohort.
Signature-Based Nomogram
Considering the contributions of other characteristics in patient prognosis, an additional endeavor was made to construct a signature-based nomogram. After screening prognostic factors by univariate Cox regression analyses (Figure 4A), a nomogram, including pyroptosis-related lncRNA signature, age, race, clinical stage, and chemotherapy, was constructed to predict the 1-, 3-, and 5-year survival probability with a C-index of 0.684 (95% CI: 0.662–0.705) in OC (Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g004.jpg)
Figure 4. Prognostic values of the pyroptosis-related lncRNA signature. (A) Multivariate Cox regression of patient characteristics and the signature in the whole cohort. (B) Nomogram constructed using patient characteristics and the signature.
Associations Between the Signature and Immune-Related Pathways
Gene set enrichment analysis (GSEA) was adopted to investigate the potential mechanism of the signature (Supplementary Table 4). Figure 5 shows the top 9 immune-associated signaling pathways, including immunoregulatory interactions between a lymphoid and a non-lymphoid cell, inflammatory response pathway, antigen processing and presentation, CD8 TCR downstream pathway, complement activation, complement and coagulation cascades, NK cell pathway, MHC class II antigen presentation, and CD8 TCR pathway. These findings suggest the potential associations between the identified signature and immune regulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g005.jpg)
Figure 5. Gene set enrichment analysis (GSEA) of the pyroptosis-related lncRNA prognostic signature. (A) Immunoregulatory interactions between a lymphoid cell and a non-lymphoid cell. (B) Inflammatory response pathway. (C) Antigen processing and presentation. (D) CD8TCR downstream pathway. (E) Complement activation. (F) Complement and coagulation cascades. (G) NK cell pathway. (H) MHC class 1II antigen presentation. (I) CD8TCR pathway.
Associations Between the Signature and Immune Infiltration
To further investigate associations between the signature and the immune microenvironment, the MCPcounter algorithm was used to calculate immune infiltration cells in OC. The difference in proportions of immune cells is illustrated in Figures 6A–I, and the relationships among different cells is illustrated in Figure 6J. T cells, CD8+ T cells, NK cell, B cell, macrophage/monocyte, myeloid dendritic cell, neutrophil, and cancer-associated fibroblast were shown to be upregulated in the high-risk group. Furthermore, we also applied the ssGSEA algorithm to recalculate immune infiltration cells in OC, which confirmed that proportions of NK cell, macrophage, and neutrophil were indeed upregulated in the high-risk group (Figures 6K–M). Moreover, comparisons of several common immune checkpoints expression levels showed that PD-1, B7-H3, and VSIR were expressed higher in the high-risk group than in the low-risk group (Figure 6N), suggesting the potential value of the signature as an immunotherapy indicator in OC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g006.jpg)
Figure 6. Interactions between the pyroptosis-related lncRNA signature and immune regulation in patients with OC. (A) Heatmap of the tumor-infiltrating cell in low-risk and high-risk patients. (B) Comparisons of T cell between the low-risk group and the high-risk group. (C) Comparisons of CD8 T cell between the low-risk group and the high-risk group. (D) Comparisons of NK cell between the low-risk group and the high-risk group. (E) Comparisons of B cell between the low-risk group and the high-risk group. (F) Comparisons of macrophage/monocyte between the low-risk group and the high-risk group. (G) Comparisons of myeloid dendritic cell between the low-risk group and the high-risk group. (H) Comparisons of neutrophil between the low-risk group and the high-risk group. (I) Comparisons of cancer-associated fibroblast between the low-risk group and the high-risk group. (J) Correlation matrix of immune cells in OC. (K) Comparisons of NK cell between the low-risk group and the high-risk group using the ssGSEA algorithm. (L) Comparisons of macrophage between the low-risk group and the high-risk group using the ssGSEA algorithm. (M) Comparisons of neutrophil between the low-risk group and the high-risk group using the ssGSEA algorithm. (N) Comparisons of multiple immune checkpoints between the low-risk group and the high-risk group, including PD-1, PD-L1, CTLA4, B7-H3, VSIR, LAG3, and TIGIT.
qRT-PCR
Quantitative real-time PRC (qRT-PCR) was performed to confirm the expression pattern of the six identified pyroptosis-related lncRNAs in two OC cell lines (Figure 7). The results confirmed that the expression levels of AC006001.2, LINC02585, AL136162.1, AC005041.3, and AL023583.1 were upregulated in OC, and that LINC02881 was downregulated, which was consistent with the results of the RNA-seq.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/793515/fmed-08-793515-HTML/image_m/fmed-08-793515-g007.jpg)
Figure 7. qRT-PCR of the six identified pyroptosis-related lncRNAs. (A) AC006001.2. (B) LINC02585. (C) AL136162.1. (D) AC005041.2. (E) AL023583.1. (F) LINC02881.
Discussion
In this study, comprehensive analyses were performed to identify and investigate the prognostic pyroptosis-related lncRNAs in OC. A novel prognostic signature consisting of 6 pyroptosis-related lncRNAs (AC006001.2, LINC02585, AL136162.1, AC005041.3, AL023583.1, and LINC02881) and a signature-based nomogram were constructed and showed promising and validated results. Subsequent GSEA and immune microenvironment analyses all indicated the potential involvement of the signature in the immunology of OC. This study may provide new insights into risk stratification and immunotherapy development of patients with OC.
Long non-coding RNAs (lncRNAs) are recognized as key players in epigenetic regulation among patients with cancer (15), including patients with OC (16). However, few studies have focused on the role of lncRNAs in the pyroptosis of OC. Li et al. have previously reported that lncRNA GAS5 could suppress ovarian cancer by inducing inflammasome formation and pyroptosis (17). Moreover, the study by Tan et al. also indicated lncRNA HOTTIP targeted the miR-148a-3p/AKT2 axis, thus inhibiting ASK1/JNK signaling and NLRP1 inflammasome-mediated pyroptosis in ovarian cancer (18). In this study, 6 pyroptosis-related lncRNAs (AC006001.2, LINC02585, AL136162.1, AC005041.3, AL023583.1, and LINC02881) were identified to be associated with progression of OC, and none of them has been investigated previously. Future studies should be encouraged to reveal the underlying mechanisms of these candidate lncRNAs in OC biology, especially cell pyroptosis. Additionally, the 6 identified lncRNAs were found to be not associated with the other 27 pyroptosis-related genes. However, it should be pointed out that the false negative rate could also potentially cause the situation. Thus, future studies should also put efforts to validate the associations between the 6 identified lncRNAs and other pyroptosis-related genes.
Several immune checkpoint inhibitors (ICIs) have been applied to various cancers with promising results; however, unfortunately, there is currently no approved immunotherapy for OC. There are strong theoretical rationales to apply immunotherapy in OC, which presents an elevated number of tumor-infiltrating lymphocytes and neoantigen load (19). However, clinical results of two anti-PD-1 drugs (pembrolizumab and nivolumab) and three anti-PD-L1 drugs (atezolizumab, avelumab, and durvalumab) were all disappointing. Therefore, some strategies were developed aiming at sensitizing OC to immunotherapy by combining it with chemotherapy, anti-angiogenetics, poly ADP-ribose polymerase inhibitor, and radiotherapy (20). Patient selection can be another significant issue when applying immunotherapy in OC. Proper patient selection for each type of immunotherapy or proper immunotherapy selection for each patient is crucial for improving patient prognosis in cancers, including OC. In this study, three immune checkpoints (PD-1, B7-H3, and VSIR) were shown to be expressed higher among patients with high-risk score. B7-H3, also known as CD276, is an immune checkpoint molecule and immunoregulatory protein, which participates in tumor microenvironment shaping and development (21). VSIR, also known as VISTA, is a well-established immune regulatory receptor functioning like a homeostatic regulator that actively normalizes immune responses (22). On the one hand, since the immunotherapy of B7-H3 and VSIR has not been developed in OC, further studies developing immunotherapies in OC can focus on these two immune checkpoints. On the other hand, the identified signature in this study can serve as an indicator to select patients with OC who may benefit from immunotherapies of PD-1, B7-H3, and VSIR.
Pyroptosis is proven to be highly related to immune response, and recent evidence also reported the immunostimulatory function of pyroptosis, as well as a significant role in promoting the efficacy of cancer immunotherapy (23, 24). The roles of some lncRNAs in cancer immunotherapy are also well-elucidated by various studies (25). However, there has been no previous robust research reporting the role of the 6 identified lncRNAs (AC006001.2, LINC02585, AL136162.1, AC005041.3, AL023583.1, and LINC02881) in immune therapies. For the very first time, this study preliminarily identified the potential contributions of these lncRNAs in immune therapies. Future studies should validate the function of these lncRNAs in immunotherapies, and focus on the specific mechanisms of lncRNAs in influencing immunotherapies, especially among patients with OC.
However, this study has certain limitations. First, the cohorts in this study were mainly based on the TCGA database. Thus, practical data were needed for further validations of the pyroptosis-related lncRNA signature in patients with OC. Second, this study failed to investigate the underlying mechanisms for regulating the pyroptosis-related lncRNA signature in tumor immune response. Moreover, further studies should focus on the potential of the pyroptosis-related lncRNA signature as an indicator of immunotherapies.
Conclusion
This study suggests that the identified pyroptosis-related lncRNA signature and the signature-based nomogram may serve as methods for risk stratification of OC. The signature is also associated with the tumor immune microenvironment, potentially providing an indicator for patient selection of immunotherapy in OC.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
YY and ZX conceived the design of the study and edited the manuscript. ZZ and YY performed the study, collected the data, and contributed to the design of the study. ZZ, YY, and ZX prepared the manuscript. All the authors made substantive intellectual contributions to this study and read and approved the final version of the manuscript.
Funding
This study was supported by grants from the National Natural Science Foundation of China (81803035), China Postdoctoral Science Foundation (2021T140754, 2020M672521), Natural Science Foundation of Hunan Province (2019JJ50932, 2020JJ5934), and Postdoctoral Science Foundation of Central South University (248485).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.793515/full#supplementary-material
Abbreviations
OC, ovarian cancer; lncRNAs, long non-coding RNAs; TPM, transcripts per kilobase million; LASSO, least absolute shrinkage and selection operator; GSEA, gene set enrichment analysis; ANOVA, one-way analysis of variance; AUC, area under curve.
References
1. El-Sherif A, El-Sherif S, Taylor AH, Ayakannu T. Ovarian cancer: lifestyle, diet and nutrition. Nutr Cancer. (2021) 73:1092–107. doi: 10.1080/01635581.2020.1792948
2. Kuroki L, Guntupalli SR. Treatment of epithelial ovarian cancer. BMJ. (2020) 371:m3773. doi: 10.1136/bmj.m3773
3. Topalian SL, Hodi FS, Brahmer JR, Gettinger SN, Smith DC, McDermott DF, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med. (2012) 366:2443–54. doi: 10.1056/NEJMoa1200690
4. Shi C, Wang Q, Rao Z, Shi Y, Wei S, Wang H, et al. Diabetes induces hepatocyte pyroptosis by promoting oxidative stress-mediated NLRP3 inflammasome activation during liver ischaemia and reperfusion injury. Ann Transl Med. (2020) 8:739. doi: 10.21037/atm-20-1839
5. Qiao L, Wu X, Zhang J, Liu L, Sui X, Zhang R, et al. α-NETA induces pyroptosis of epithelial ovarian cancer cells through the GSDMD/caspase-4 pathway. FASEB J. (2019) 33:12760–7. doi: 10.1096/fj.201900483RR
6. Wang D, Weng Y, Zhang Y, Wang R, Wang T, Zhou J, et al. Exposure to hyperandrogen drives ovarian dysfunction and fibrosis by activating the NLRP3 inflammasome in mice. Sci Total Environ. (2020) 745:141049. doi: 10.1016/j.scitotenv.2020.141049
7. Hou J, Lei Z, Cui L, Hou Y, Yang L, An R, et al. Polystyrene microplastics lead to pyroptosis and apoptosis of ovarian granulosa cells via NLRP3/Caspase-1 signaling pathway in rats. Ecotox Environ Safe. (2021) 212:112012. doi: 10.1016/j.ecoenv.2021.112012
8. Chen W, Yang J, Fang H, Li L, Sun J. Relevance function of Linc-ROR in the pathogenesis of cancer. Front Cell Dev Biol. (2020) 8:696. doi: 10.3389/fcell.2020.00696
9. Oncul S, Amero P, Rodriguez-Aguayo C, Calin GA, Sood AK, Lopez-Berestein G. Long non-coding RNAs in ovarian cancer: expression profile and functional spectrum. Rna Biol. (2020) 17:1523–34. doi: 10.1080/15476286.2019.1702283
10. Ye Y, Dai Q, Qi H. A novel defined pyroptosis-related gene signature for predicting the prognosis of ovarian cancer. Cell Death Discov. (2021) 7:71. doi: 10.1038/s41420-021-00451-x
11. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
12. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et 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
13. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. (2016) 17:218. doi: 10.1186/s13059-016-1070-5
14. Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. (2009) 462:108–12. doi: 10.1038/nature08460
15. Vishnubalaji R, Shaath H, Elango R, Alajez NM. Noncoding RNAs as potential mediators of resistance to cancer immunotherapy. Semin Cancer Biol. (2020) 65:65–79. doi: 10.1016/j.semcancer.2019.11.006
16. Calanca N, Abildgaard C, Rainho CA, Rogatto SR. The interplay between long noncoding RNAs and proteins of the epigenetic machinery in ovarian cancer. Cancers. (2020) 12:2701. doi: 10.3390/cancers12092701
17. Li J, Yang C, Li Y, Chen A, Li L, You Z. LncRNA GAS5 suppresses ovarian cancer by inducing inflammasome formation. Biosci Rep. (2018) 38:BSR20171150. doi: 10.1042/BSR20171150
18. Tan C, Liu W, Zheng Z, Wan X. LncRNA HOTTIP inhibits cell pyroptosis by targeting miR-148a-3p/AKT2 axis in ovarian cancer. Cell Biol Int. (2021) 45:1487–97. doi: 10.1002/cbin.11588
19. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Programmed cell death 1 ligand 1 and tumor-infiltrating CD8+ T lymphocytes are prognostic factors of human ovarian cancer. Proc Natl Acad Sci USA. (2007) 104:3360–5. doi: 10.1073/pnas.0611533104
20. Macherla S, Laks S, Naqash AR, Bulumulle A, Zervos E, Muzaffar M. Emerging role of immune checkpoint inhibitors in the treatment of ovarian cancer. Expert Opin Emerg Dr. (2020) 25:445–53. doi: 10.1080/14728214.2020.1836155
21. Zhou W, Jin W. B7-H3/CD276: an emerging cancer immunotherapy. Front Immunol. (2021) 12:701006. doi: 10.3389/fimmu.2021.701006
22. Huang X, Zhang X, Li E, Zhang G, Wang X, Tang T, et al. VISTA: an immune regulatory protein checking tumor and immune cells in cancer immunotherapy. J Hematol Oncol. (2020) 13:83. doi: 10.1186/s13045-020-00917-y
23. Hsu SK, Li CY, Lin IL, Syue WJ, Chen YF, Cheng KC, et al. Inflammation-related pyroptosis, a novel programmed cell death pathway, and its crosstalk with immune therapy in cancer treatment. Theranostics. (2021) 11:8813–35. doi: 10.7150/thno.62521
24. Chen X, Chen H, Yao H, Zhao K, Zhang Y, He D, et al. Turning up the heat on non-immunoreactive tumors: pyroptosis influences the tumor immune microenvironment in bladder cancer. Oncogene. (2021) 40:6381–93. doi: 10.1038/s41388-021-02024-9
Keywords: ovarian cancer, pyroptosis, lncRNA, prognostic signature, tumor immune microenvironment
Citation: Zhang Z, Xu Z and Yan Y (2022) Role of a Pyroptosis-Related lncRNA Signature in Risk Stratification and Immunotherapy of Ovarian Cancer. Front. Med. 8:793515. doi: 10.3389/fmed.2021.793515
Received: 12 October 2021; Accepted: 13 December 2021;
Published: 12 January 2022.
Edited by:
Alice Chen, National Cancer Institute (NCI), United StatesReviewed by:
Yanyong Yang, Second Military Medical University, ChinaHaicheng Tang, Fudan University, China
Copyright © 2022 Zhang, Xu and Yan. 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: Yuanliang Yan, yanyuanliang@csu.edu.cn