- 1Morphology Laboratory, Medical College of Nantong University, Nantong, Jiangsu, China
- 2Department of Clinical Medicine, Ningbo College of Health Science, Ningbo, Zhejiang, China
Purpose: Low-grade gliomas (LGG), which are malignant primary brain tumors, are more prevalent in young adults. Pyroptosis, an inflammatory form of programmed cell death, has been shown in recent years to be directly associated with tumor growth and tumor microenvironment (TME). However, the correlation between LGG and pyroptosis remained to be explored. In this research, we explored pyroptosis-related gene expression patterns and their prognostic significance based on transcriptome profiles and clinical data in LGG.
Methods: We identified 31 pyroptosis-related genes differentially expressed at the mRNA level between the data of LGG patients from TCGA and the data of normal brain tissues from GTEx. Univariate Cox regression analysis was used to screen 16 differentially expressed genes (DEGs) based on survival data. Next, the prognostic model was established using LASSO Cox regression, which divided LGG patients into high- and low- risk subgroups and showed an independent prognostic value for overall survival (OS) combined with clinical factors in the CGGA test cohort. Pyroptosis and immune cells were correlated through the CIBERSORT R package and the TIMER database.
Results: Based on the analyses of 523 LGG and 1152 normal tissues, nine significant differential genes were identified. The AUC remained at about 0.74 when combined with the risk score and clinical factors. Enrichment analyses revealed that DEGs were mainly enriched in cytokine-cytokine receptor interactions, immune response and chemokine signaling pathways. Immune cell enrichment analysis demonstrated that scores for most immune cell types differed significantly between the high-and low-risk groups, and further infiltrating analysis showed obvious differences between these two risk subgroups.
Conclusion: Pyroptosis-related genes play a pivotal role in LGG and are associated with tumor immunity, which may be beneficial to the prognosis and immunotherapy of LGG.
Introduction
According to the World Health Organization (WHO), low-grade gliomas (LGG) account for about 20% of gliomas and are often diagnosed in young patients (1, 2). While most LGGs are inactive in patients, some cases rapidly deteriorate to neurological disorders and death in the short term (3). Based on prior clinical experience, an accurate and sensitive prognostic model may serve as the foundation for more effective diagnostic and therapeutic approaches.
According to the WHO classification guide, many molecular biomarkers are used to detect tumor occurrence, proliferation, infiltration and migration. These biomarkers provide clues by which to assess the comprehensive status of gliomas and then formulate personalized precision treatment (4). Previous studies on the biomarkers of glioma mainly focused on epidermal growth factor receptor (EFGR) amplification, chromosome 1p/19q deletion, O6-methylguanine DNA methyltransferase (MGMT) promoter methylation, and isocitrate dehydrogenase (IDH) mutations (5, 6). The literature confirms many studies with an appreciable quantity of reporting on prognostic models related to LGG, such as those involving metabolism (7), autophagy (8), ferroptosis (9), cuproptosis (10), and so on. While such markers and models are present in a wide variety of gliomas, they are insufficient for predicting the prognosis of LGG. Even though there have been many studies on prognostic models for LGG, no unified prognostic models have been applied in clinical practice. The complexity of glioma development and biological pathways call for the discovery of novel and specific LGG biomarkers.
Studies on malignancies of the nervous system have consistently highlighted the phenomena of programmed cell death (PCD) (11, 12). First, the formation of nerve tissue depends on the regulation of PCD. Abnormal death of nerve cells can promote tumor growth (11, 12). Second, tumor cells can resist or escape from cell-killing effects mediated by immune cells through a variety of adaptive mechanisms (13, 14). Thus, the detection and regulation of PCD performed an essential function in the diagnosis and treatment of nervous system tumors. Newly discovered close relationships among different kinds of tumors have made pyroptosis a hotspot as a new kind of PCD. Inflammasomes have been found in a large number of tumor cells, which activate the caspase pathways to initiate pyroptosis. Furthermore, pyroptosis can regulate the TME through proinflammatory effects to play a dual role in carcinomatosis and cancer promotion. This all indicates the potential of pyroptosis as a new biomarker of nervous system tumors (15, 16). Li et al. have constructed a prognostic model of glioblastoma (GBM, grade IV glioma) based on pyroptosis-related genes to provide accurate one, three, and five years overall survival rates (OS) of malignant gliomas (17). However, the indicators of LGG patients are still unclear, making it worthwhile to investigate the role of pyroptosis-related genes. Therefore, we developed a risk model for pyrolysis-related genes, which provides a new perspective for improving the prognosis and treatment of LGG.
Material and methods
Datasets
The transcriptome profile and clinical data of LGG patients were obtained from The Cancer Genome Atlas (18) (TCGA, https://www.tcga.org/) database and the Chinese Glioma Genome Atlas Project (19) (CGGA, http://www.cgga.org.cn/, Supplementary Table 1). Transcriptome profiles of 1152 normal brain tissues were collected from the Genotype-Tissue Expression Project (20) (GTEx, https://xenabrowser.net/) database. We acquired 523 low-grade gliomas in the TCGA database and 422 patients in the CGGA database to construct a validation set.
DEGs identification
The gene expression data were normalized in all sets to fragments per thousand base million (FPKM) (21). The differentially expressed genes (DEGs) with P< 0.05 were identified by using the ‘limma’ R package (22). There are 33 pyroptosis-related genes were retrieved in Pubmed with ‘pyroptosis’ as the key word (23–25). In addition, a protein interaction network (PPI) was construct by using an interactive gene/protein search tool (String, http://www.string-db.org/) to analyze the correlation among DEGs. A Pearson correlation analysis diagram was drawn among DEGs.
Development and validation of prognostic models for pyroptosis-related gene
In the training set, univariate COX regression analysis was adopted to determine the correlation between DEGs and overall survival (OS) to evaluate the prognostic value of DEGs in LGGs with the truncation condition as P< 0.05 and HR unequal to 1 (26). Next, Lasso were used to adjust COX proportional hazard regression to avoid overfitting (27). The above operations were performed through ‘glmnet’ R package, and penalty parameter λ was determined through minimum criterion (28). The prognostic model was constructed based on multivariate COX regression analysis (29). Based on the median risk score, LGG patients were assigned into low-risk group and high-risk group. The risk score formula following: risk score = Xi ×Yi, where n is the number of surviving genes after Lasso regression, Y is gene expression level, and X is coefficients. The OS time between subgroups was compared by using Kaplan-Meier analysis. The 3-year and 5-year ROC analyses were performed using the “survivalroc” R package and the area under the curve (AUC) was calculated. We used CGGA data as validation set to calculate the risk score. The Kaplan-Meier curves and ROC (Receiver Operating Characteristic) were plotted.
Analysis of functional enrichment
The GO (Gene Ontology) (30) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis (31) were performed by using ‘ClusterProfiler’ R package based on the difference between risk subgroups (|log2FC|≥mean(|log2FC|) + 4*sd(|logFC|);P<0.05).
Estimation of tumor−infiltrating immune cells
According to the transcription profile of LGGs in TCGA, the relative proportion of immune infiltrating cells in all tumor samples was calculated by using cell type identification analysis in the ‘CIBERSORT’ R package (32, 33). We used the Wilcoxon rank-sum test to evaluate differences in immune cell infiltration levels between risk subgroups.
TIMER database and GDSC database
The TIMER database provided a reliable estimate of the level of immune infiltration for tumor-immune interactions (34). The level of tumor immune infiltration and the correlation between gene expression were calculated by using ‘GENE’ module of the TIMER database (35). We screened a wide range of drugs from the GDSC database (36) and calculated the IC50 value of drugs using the pRRophetic algorithm (37) for patients with LGG.
Quantitative real-time polymerase chain reaction
The HEB (human normal glial cell line) was obtained from Mingzhoubio (Ningbo, China) and was cultured H-DMEM medium with 10% FBS (fetal bovine serum). Human glioma cell lines, including U251 and U87, were obtained from the ATCC (American Type Culture Collection; Manassas, VA, USA) and were cultured in Roswell Park Memorial Institute (RPMI)-1640 with 10% fetal bovine serum (FBS). Total RNA from the cell lines was extracted by applying RNA simple Total RNA Kit (Tiangen, China). Subsequently, Total RNA was reverse transcribed with PrimeScript RT reagent kit (Takara, Otsu, Japan) to obtain cDNA. Using 2 μL cDNA and SYBR Premix Ex Taq (Takara, Otsu, Japan) and primers, the expression of the target genes was determined using the Biosystems StepOne Plus real-time PCR system (Life Technologies, Grand Island, NY, USA). The Primers of the target genes were obtained from Sangon Biotechnology (Shanghai, China), and the sequences are shown in Supplementary Table 2.
Statistical analysis
Normally distributed variables were compared using Student’s t-test. Wilcoxon test was used to compare non-normally distributed data. Utilizing the Survminer package from R, we estimated the OS status between the two subgroups by Cox regression analysis and Kaplan–Meier curves. Ggplot2 package from R were used to plot the figures, thus visualizing our data. Test level was set at both sides α =0. 05, P < 0.05 was considered statistically significant unless otherwise specified.
Results
The flow chart of this study is shown in Figure 1.
Identification of DEGs between tumor and normal tissues
First, we compared 33 pyroptosis-related genes in the pooled GTEx (Genotype-Tissue Expression) and TCGA (The Cancer Genome Atlas) data from 523 LGG and 1152 normal tissues and then identified 31 differentially expressed genes (DEGs) under the condition of P<0.05 (Figure 2A and Supplementary Table 3) besides ELANE and GSDMD. Among these DEGs, upregulated of 27 genes, while downregulated of 4 genes, in LGG tissues. The expression levels of DEGs were presented as heatmaps in Figure S1. We then constructed protein interaction networks (PPI) to investigate the interactions of DEGs (Figure 2B). As a result of the analysis, most DEGs showed obvious positive correlations, among which CASP4 and SCAF11 were significantly positively correlated (Cor = 0.87), while PRKACA was not connected to other DEGs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g002.jpg)
Figure 2 Expressions and screening of pyroptosis-related genes in low-grade glioma. (A) Box plot demonstrating 31 differentially expressed pyroptosis-related genes between normal tissue (n=1152) and lower-grade glioma (n=523). **P< 0.01; ∗∗∗∗P< 0.0001, ns: not significant. The green box shows normal tissue, and the red box shows tumor tissue. (B) PPI network showing the interactions of DEGs (interaction score=0.9). The darker the color of a node in a PPI network, the closer its connection to other nodes. (C) Univariate COX analysis of 16 DEGs (CASP1, CASP3, CASP4, CASP5, CASP6, CASP8, CASP9, GSDMC, IL18, IL6, NLRC4, NOD1, PLCG1, PRKACA, PYCARD, and SCAF11). P<0.05; HR, Hazard Ratio. (D) Cross-validation for tuning parameter selection in LASSO regression. (E) LASSO analysis of 9 prognostic pyroptosis-related genes.
Development of a prognostic gene model in the training and testing sets
Next, A total of 388 LGG patients were matched based on their survival data. We used univariate Cox regression for preliminary analysis of DEGs, and 16 survival-related genes were screened at P<0.05 (Figure 2C). According to the optimum Λ value (Figures 2D, E), a prognostic model of 9 genes (CASP3, CASP4, CASP8, CASP9, GSDMC, IL18, IL6, PLCG1 and PRKACA) was construct by using Least absolute shrinkage and selection operator (LASSO) Cox regression analysis, thus calculating the risk score for each LGG patient.
A total of 388 LGG patients were divided into high- (n=194) and low-risk (n=194) subgroups according to the median score calculated by the risk score formula(Figures 3A, B). Time-dependent Kaplan-Meier curves and receiver operating characteristic (ROC) curves were constructed to evaluate the sensitivity of the prognostic model (Figures 3C, E). The results showed differences in the survival curve between the low-risk group and the high-risk group (P<0.001) as the one year, three years, and five years AUC values were 0.856, 0.832 and 0.742, respectively. In addition, we matched 420 LGG patients in CGGA (Chinese Glioma Genome Atlas) as the validation set to further verify the model, resulting in one year, three years, and five years AUC values of 0.596, 0.636 and 0.659, respectively (Figures 3B, D, F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g003.jpg)
Figure 3 Subgroups of LGGs based on the construction of risk signature. (A, B) Distribution of risk score, differences in survival between the high- and low-risk groups and heatmap of the expression patterns of 9 pyroptosis-related genes in the training set (A) and testing set (B). Red dots indicate dead, whereas blue dots indicate live. (C, D) Kaplan-Meier curve for the OS of LGG patients in the high- (n=194) and low-risk (n=194) groups in the TCGA training cohort (C) and the CGGA (n=420) testing cohort (D). The red curve represents the high-risk group, and the blue curve represents the low-risk group. (E, F) Time-dependent ROC curve based on the prognostic model regarding OS and survival status in the TCGA training cohort (E) and the CGGA testing cohort (F). Green, blue and red curves represent 1-year, 3-year and 5-year, respectively.
Clinical evaluation of the prognostic risk model
In order to verify practical value of prognostic model, univariate and multivariate COX analysis were used to evaluate whether the risk score of the model could also be an independent prognostic point compared with such clinical factors as grade, gender, age and radiation therapy. Univariate Cox regression analysis showed that risk scores were not associated with poor survival in both TCGA and CGGA groups of patients (Figure S2A). Interestingly, multivariate COX analysis showed similar trends to univariate COX regression results, suggesting that the prognosis of LGG is closely related with pyroptosis. According to ROC curve analysis, the AUC value remained at about 0.8 when combined with the risk score and clinical factors (Figures 4A–C). In addition, we used the TCGA and CGGA cohorts to analyze clinical factors and found that age, grade, and radiotherapy differed in the distribution of low- and high-risk categories (Figures S2B, C). We then integrated clinical variables and created a nomogram in the TCGA cohort (Figure 4D). The total score for each LGG patient was acquired by combining the scores of each prognostic standard in the nomogram. A higher total score would indicate patients with a worse prognosis. As a comparison, the predicted survival rate showed more consistency with the observed survival rate (Figures 4E–G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g004.jpg)
Figure 4 Assessment of the risk model combined with clinical factors. (A-C) Time-dependent ROC to evaluate prognostic power based on risk score and clinical factors at (A) 1-year, (B) 3-year, and (C) 5-year. Different colored curves represent different clinical factors. (D) A nomogram consisting of risk score and other clinical indicators for predicting 1-, 3-, and 5-year OS of primary LGG based on the TCGA cohort. (E-G) Nomogram calibration plots for predicting OS in the TCGA cohort at (E) 1-year, (F) 3-year, and (G) 5-year. The Red line indicates actual survival.
Functional and mutational analyses based on the risk model
Subsequently, the DEGs were extracted and further explore the signal pathways related to different subgroups in the risk model. There are 181 differential genes were identified in the TCGA cohort, among which up-regulated of 109 genes and down-regulated of 72 genes in the high-risk group. Based on the differences between the risk subgroups, GO enrichment (Figure 5A) and KEGG (Figure 5B) pathway analysis was conducted. The results revealed the differential genes were mainly enriched in chemokine signaling pathways, cytokine-cytokine receptor interactions and immune response (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g005.jpg)
Figure 5 Functional and mutational analyses based on the risk model. (A) The enriched item in gene ontology analysis. Blue represents BP (biological process), orange represents CC (cellular component), and the green represents MF (molecular function). (B) The enriched item in Kyoto Encyclopedia of Genes and Genomes analysis. The color indicates the size of the P-value, and the size of the circle indicates the number of genes. (C) Mutation profile in the low-risk group. (D) Mutation profile in the low-risk group. The small figure above shows the TMB, whereas the number on the left shows the mutation frequency of each gene, and the figure on the right shows the proportion of each variant.
Combined with the results of gene enrichment in immune-related pathways, as suggested by functional analysis, we first compared the effects of cell mutations that tumors may induce. We used the ‘maftools’ R package on the high-risk group (Figure 5C) and the low-risk group (Figure 5D) to achieve a visual analysis of differences in the distribution of somatic mutations. Results showed the mutation rates among high-risk groups (46.06%) and low- (48.82%) to be relatively close.
Immune activity and PCR analysis based on subgroups of the risk model
In order to further investigate the correlation among risk score of prognostic model and immune cell infiltration, ssGSEA module in the ‘GSVAR’ R package was used. This module quantifies the function of immune cells and the signal pathways in tumor samples. The scores for most immune cell types differed significantly between the low-risk and high-risk groups, which were confirmed in the validation set. The relative proportion of 22 types of immune cells associated with each LGG patient was calculated by using CIBERSOR and Estimate algorithms. The analysis of the correlation between the risk score and the degree of immune cell infiltration showed that many immune cells differed in the degree of infiltration between subgroups (Figure 6A); especially, plasma cells, TFH cells, and M2 macrophages were significantly upregulated (P<0.05). As a comparison, the Estimate results suggested the scores of patients were significantly higher in the high-risk group (Figures 6B–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g006.jpg)
Figure 6 Immunoassay in the risk subgroups. (A) Boxplot showing the differential abundance of 22 infiltrating immune cells calculated by CIBERSORT between the high- and low-risk groups in LGG. ∗P< 0.05; ∗∗∗P< 0.001; ∗∗∗∗P< 0.0001; ns: not significant. (B-D) The higher expression level of Estimate Scores (B), Immune Score (C), and Stroma Score (D) correlated with the high-risk group. The green box shows the low-risk group, and the red box shows the high-risk group.
We used TIMER database to analyses the relationship among 9 prognostic genes and immune cell abundance, and the results revealed that 9 genes of prediction model were all closely related to immune cell abundance (P<0.01, Figures 7A-C and Figure S4). To be more specific, the expression of PLCG1 was positively correlated with neutrophils, macrophages, myeloid dendritic cells, B cells and CD4+ T cells (Figure 7A), but negatively correlated with CD8+ T cells (Figure 7A). The expressions of CASP3 (Figure S4A), CASP4 (Figure S4B), CASP8 (Figure S4C), IL18 (Figure 7B), and IL6 (Figure 7C) were positively correlated with neutrophils. At the level of cellular immune infiltration, it is worth noting that CASP3 (Figure S4A) and CASP4 (Figure S4B) were positively correlated with macrophages and dendritic cells, but negatively correlated with CD4+ T cells and CD8+ T cells.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1015850/fonc-12-1015850-HTML-r1/image_m/fonc-12-1015850-g007.jpg)
Figure 7 Expression of pyroptosis-related genes in immune cells and cell lines. (A-C) Association between the abundance of immune cells and the expression of three pyroptosis-related genes for PLCG1 (A), IL18 (B), and IL6 (D). (D-F). Quantitative Real-time PCR in cell lines for PLCG1 (D), IL18 (E), and IL6 (F). ∗P< 0.05; ∗∗∗P< 0.001; ∗∗∗∗P< 0.0001.
Furthermore, we randomly selected three pyroptosis-related genes (PLCG1, IL18, and IL6) in the prognostic model and confirmed the expression of three genes in cell lines through qPCR assays. Compared with the human normal glial cell line (HEB), PLCG1 (Figure 7D), IL18 (Figure 7E), and IL6 (Figure 7F) were highly expressed in cancer cell lines (U251 and U87). These three genes were therefore tightly associated with immunological activity in addition to being substantially expressed in cancer cell lines.
Discussion
In this study, a prognostic model of 9 genes (CASP3, CASP4, CASP8, CASP9, GSDMC, IL18, IL6, PLCG1, and PRKACA) was constructed based on bioinformatics analysis with the survival data of LGG patients in the TCGA and CGGA cohorts and validated it by both clinical prognosis and correlation with tumor immunity. Although there are studies that have reported prognostic models related to LGG (38–42), these studies only used data from a single database - TCGA, thus the accuracy of the constructed models was limited. On the contrary, the present study used an external database (CGGA), which makes it more reliable for modeling. Furthermore, the constructed prognostic model was based on a newly discovered PCD, pyroptosis, making this current research fundamentally different from the others.
Initially, pyroptosis was thought to be PCD-activated by Caspase-1 only. However, Caspase-3/4/5/6/8/9/11 were successively identified to cause pyroptosis in different cell types; notably, Caspase-3 was found to be an essential molecule for membrane blistering and activation of apoptosis (43). Caspase-3 can cleave GSDME to make an alternative pathway to induce pyroptosis (44), which may be rapidly triggered by the Caspase-4/9 complex. Caspase-4 can also be turned on under infectious conditions to make a non-classical pyroptotic pathway by cleaving GASMD, along with caspase-5/11 (45, 46). Meanwhile, GASMD can stimulate macrophages to induce pyroptosis in a PLCG1 (phospholipase)-dependent pathway (47). Besides cleaving GSDMD directly, caspase-8 cleaves GSDMC when TNF-α is released to promote pyroptosis (48). By comparison, Xiao et al. (49) found that Caspase-3 influences dying glioma cells after chemotherapy by promoting tumor angiogenesis, leading to the recurrence of gliomas. Caspase-8 can also stimulate the NF-kB pathway and upregulate the secretion of IL-8, IL-1β, IL-6, VEGF and MCP-1, improving temozolomide resistance in gliomas (50). Caspase-9 had been reported to repair mitochondrial damage by maintaining autophagy, which then avoided influencing NLRP3 inflammasome to produce pyroptosis (51). These investigations all suggested that via inducing pyroptosis, the caspase family plays a dual role in gliomas.
Previous studies have revealed that proinflammatory effects of pyroptosis are closely related to the regulation of the tumor immune microenvironment (TIME) (52). In the TIME, tumor cells are protected by chemokines, cytokines, stromal cells, and metabolites that allow tumor cells to survive (53). Through damage-associated molecular patterns (DAMPs) released after cellular osmotic lysis, pyroptosis can reprogram TIME into an immune stimulatory state, so as to inhibit tumor cell growth and metastasis. Under the action of inflammatory factors, it can also promote the growth of tumor cells (54). Pyroptosis appeared to have a specific effect on TIME and facilitate immune surveillance (54). There has research has revealed NLRP3 inflammasome activation could affect pyroptosis or hyperactivity, resulting in a cascade of immune or inflammatory responses through the release of interleukins IL-β and IL-18, affecting antitumor immunity (55). In addition, IL-18 is upregulated in tumor-infiltrating lymphocytes (TILs) secreted by inflammasomes (56, 57). It has a protective pro-inflammatory effect, stimulating the generation of MDSCs to accelerate tumor progression (58). IL-6 performed an essential function in the regulation of macrophages and lymphocytes in the TME. Past studies reported that IL-6 directly influences the invasion of gliomas by activating the STAT3 pathway (59). PLCG1 and PRKACA are traditionally considered to be the executive molecules of apoptosis by signal transduction of phosphokinases. These molecules play regulatory roles in macrophage differentiation and inflammatory response to regulate the tumor microenvironment (60). Our prognostic model quantified the function of immune cells and signal pathways in tumor samples such that the scores for most immune cell types differed significantly between the low-risk and high-risk groups. In summary, our studies revealed different pyroptosis pathways (classical/non-classical/alternative/kinase-related), objectively reflecting the various factors influencing the development of LGG with a considerable reference value. Interestingly, Shao et al. (61) reported non-expression of GSDMD in the vast majority of tumor cells, suggesting that tumor cells can induce pyroptosis through more than one pathway, which is consistent with differential gene expression that forms the basis of our model. Though there is a slight lack of accuracy in judging the 5-year survival rate of LGGs, we still have a considerable advantage in the prognosis of LGG by considering the possibility of high-risk patients deteriorating into high-grade glioma within a short timeframe.
In order to analyze the immune infiltration of LGG tissue, our model compared different risk subgroups and found that naive B cells, plasma cells, TFH cells, and M2 macrophages were significantly upregulated in the high-risk group. Studies have shown accumulation of M2 macrophages can induce chronic inflammation, which is beneficial to the growth of invasive tumors. The intratumoral density of macrophages is highest in malignant gliomas, which is indicative of a negative correlation with the survival rate of patients (62). On the other hand, when combined with the analysis of the TIMER database, we further observed that the risk subgroups showed differences in the level of CD8+T cell infiltration. This is consistent with J. Robert Kane et al. whose study revealed that the accumulation of CD8+T cells in the TME indicates a better prognosis, while their absence is conducive to the growth of gliomas (63). In summary, we can speculate that pyroptosis regulates the TIME, which can be quantitated through our pyroptosis-related risk model, giving, in turn, a better prognostic picture than previous markers and models. Our further studies will elucidate the specific mechanisms.
Current tumor treatment strategies include, for example, immune checkpoint inhibitors, CAR-T, cytokines, virus lysis, and tumor vaccines, all widely used as representatives of immunotherapies in recent years (64). These methods have shown obvious carcinostatic effects in several tumors, but challenges remain in applying them to the treatment of glioma. One such challenge involves the physical blocking effect of the blood-brain barrier on drugs and cells. Another challenge arises from most gliomas that appear as a ‘cold tumor’ environment showing a higher level of immunosuppression, but fewer immune targets (65). Hung et al. addressed these challenges and found that large numbers of PD-L1 molecules enter the nucleus under tumor hypoxia. The complex of PD-L1 and phosphorylated STAT3 can mediate the expression of GSDMC and then induce tumor cell pyroptosis. In addition, pyroptosis combined with PD-L1 inhibitor can promote the conversion of ‘cold tumors’ into ‘hot tumors’. Liu et al. (66) reported that CAR-T cells activate the pyroptosis pathway of Caspase-3-GSDME by releasing GzmB, while co-culture experiments in vitro found that CAR-T can stimulate macrophages to release IL-6 and IL-1β by activating the Caspase-1-GSDMD pathway (67). Based on the studies above, our LGG prognostic model suggests a new approach to the treatment of gliomas by using it as a guide for the management of current mainstream immunotherapy.
There are still several limitations in our study. The analytical data were from public databases (TCGA and CGGA). To further confirmed the accuracy of pyroptosis-related gene prognostic model in predicting prognosis, the LGG tissues need to be collected and closely followed. Due to the lack of time and money for follow-up, the specific mechanism of pyroptosis-related genes has not been explored yet.
Conclusion
Activation of inflammasome-related sensing proteins causes pyroptosis, a form of immune cell death. There are several that detect different substances. We identified differently expressed pyroptosis-related genes that may be involved in LGG. Pyroptosis may be an alternative therapeutic target since differentially expressed genes associated with it have significant predictive values for patient survival.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.tcga.org/, http://www.cgga.org.cn/ and https://xenabrowser.net/.
Author contributions
TT conceived, designed, and supervised the study. HL performed data analysis, arranged the figures and drafted the manuscript. All authors reviewed and approved the final manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Zhejiang Provincial Natural Science Foundation of China (Grant Number LY19H100001), Zhejiang Province Traditional Chinese Medicine Science and Technology Project (Grant Number 2022RC255) and Zhejiang Province Medical and Health Science and Technology Project (Grant Number 2022ZB331).
Acknowledgments
We appreciate the information provided by The Cancer Genome Atlas.
Conflict of interests
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/fonc.2022.1015850/full#supplementary-material
Supplementary Figure 1 | Expressions of pyroptosis-related genes in LGGs. (A) Heatmap of RNA expression levels of pyroptosis-related genes. Green represents a low expression level, and red represents a high expression level. (B) Pearson correlation analysis of the 33 pyroptosis-related genes in LGGs. Red indicates the degree of positive correlation, and blue indicates the degree of negative correlation.
Supplementary Figure 2 | Analysis of clinical factors in subgroups. (A) Univariate analysis and multivariate analysis for hazard ratio values of risk score and clinical characters in the training set (CI, confidence interval). (B-C) Relationship between prognostic gene expression and clinical factors. Green represents a low expression level, and red represents a high expression level.
Supplementary Figure 3 | Comparison of ssGSEA scores in high- and low-risk groups. (A) Training set. (B) Testing set. The green box shows the low-risk group, and the red box shows the high-risk group. ∗P< 0.05; ∗∗P< 0.01; ∗∗∗P< 0.001; ns: not significant.
Supplementary Figure 4 | Association between the abundance of immune cells and the expression in LGGs. (A) CASP3. (B) CASP4. (C) CASP8. (D) CASP9. (E) GSDMC. (F) PRKACA.
References
1. Ostrom QT, Gittleman H, Fulop J, Liu M, Blanda R, Kromer C, et al. CBTRUS statistical report: Primary brain and central nervous system tumors diagnosed in the united states in 2008-2012. Neuro Oncol (2015) 17 Suppl 4:iv1–iv62. doi: 10.1093/neuonc/nov189
2. Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, et al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol (2007) 114:97–109. doi: 10.1007/s00401-007-0243-4
3. Schiff D, Brown PD, Giannini C. Outcome in adult low-grade glioma: The impact of prognostic factors and treatment. Neurology (2007) 69:1366–73. doi: 10.1212/01.wnl.0000277271.47601.a1
4. Louis DN, Perry A, Wesseling P, Brat DJ, Cree IA, Figarella-Branger D, et al. The 2021 WHO classification of tumors of the central nervous system: A summary. Neuro Oncol (2021) 23:1231–51. doi: 10.1093/neuonc/noab106
5. Appin CL, Brat DJ. Biomarker-driven diagnosis of diffuse gliomas. Mol Aspects Med (2015) 45:87–96. doi: 10.1016/j.mam.2015.05.002
6. Jones DTW, Bandopadhayay P, Jabado N. The power of human cancer genetics as revealed by low-grade gliomas. Annu Rev Genet (2019) 53:483–503. doi: 10.1146/annurev-genet-120417-031642
7. Li G, Wu Z, Gu J, Zhu Y, Zhang T, Wang F, et al. Metabolic signature-based subtypes may pave novel ways for low-grade glioma prognosis and therapy. Front Cell Dev Biol (2021) 9:755776. doi: 10.3389/fcell.2021.755776
8. Feng S, Liu H, Dong X, Du P, Guo H, Pang Q. Identification and validation of an autophagy-related signature for predicting survival in lower-grade glioma. Bioengineered (2021) 12:9692–708. doi: 10.1080/21655979.2021.1985818
9. Zhao J, Liu Z, Zheng X, Gao H, Li L. Prognostic model and nomogram construction based on a novel ferroptosis-related gene signature in lower-grade glioma. Front Genet (2021) 12:753680. doi: 10.3389/fgene.2021.753680
10. Yan X, Wang N, Dong J, Wang F, Zhang J, Hu X, et al. A cuproptosis-related lncRNAs signature for prognosis, chemotherapy, and immune checkpoint blockade therapy of low-grade glioma. Front Mol Biosci (2022) 9:966843. doi: 10.3389/fmolb.2022.966843
11. Bögler O, Mikkelsen T. Angiogenesis and apoptosis in glioma: Two arenas for promising new therapies. J Cell Biochem (2005) 96:16–24. doi: 10.1002/jcb.20475
12. Amirlak B, Couldwell WT. Apoptosis in glioma cells: review and analysis of techniques used for study with focus on the laser scanning cytometer. J Neurooncol (2003) 63:129–45. doi: 10.1023/a:1023906316524
13. Roesch S, Rapp C, Dettling S, Herold-Mende C. When immune cells turn bad-Tumor-Associated Microglia/Macrophages in glioma. Int J Mol Sci (2018) 19(2):436. doi: 10.3390/ijms19020436
14. Grabowski MM, Sankey EW, Ryan KJ, Chongsathidkiet P, Lorrey SJ, Wilkinson DS, et al. Immune suppression in gliomas. J Neurooncol (2021) 151:3–12. doi: 10.1007/s11060-020-03483-y
15. Minton K. Pyroptosis heats tumour immunity. Nat Rev Drug Discovery (2020) 19:309. doi: 10.1038/d41573-020-00062-8
16. Wang YY, Liu XL, Zhao R. Induction of pyroptosis and its implications in cancer management. Front Oncol (2019) 9:971. doi: 10.3389/fonc.2019.00971
17. Li XY, Zhang LY, Li XY, Yang XT, Su LX. A pyroptosis-related gene signature for predicting survival in glioblastoma. Front Oncol (2021) 11:697198. doi: 10.3389/fonc.2021.697198
18. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD, et al. An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell (2018) 173:400–416.e411. doi: 10.1016/j.cell.2018.02.052
19. Zhao Z, Zhang KN, Wang Q, Li G, Zeng F, Zhang Y, et al. Chinese Glioma genome atlas (CGGA): A comprehensive resource with functional genomic data from Chinese glioma patients. Genomics Proteomics Bioinf (2021) 19:1–12. doi: 10.1016/j.gpb.2020.10.005
20. Burgess DJ. Reaching completion for GTEx. Nat Rev Genet (2020) 21:717. doi: 10.1038/s41576-020-00296-7
21. Zhao Y. TPM, FPKM, or normalized counts? A comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository. J Transl Med (2021) 19:269. doi: 10.1186/s12967-021-02936-w
22. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43:e47. doi: 10.1093/nar/gkv007
23. Shi J, Gao W, Shao F. Pyroptosis: Gasdermin-mediated programmed necrotic cell death. Trends Biochem Sci (2017) 42:245–54. doi: 10.1016/j.tibs.2016.10.004
24. Bergsbaken T, Fink SL, Cookson BT. Pyroptosis: Host cell death and inflammation. Nat Rev Microbiol (2009) 7:99–109. doi: 10.1038/nrmicro2070
25. Xia X, Wang H, Lengerich BJ, Aragam B, Xing EP. The role of pyroptosis in cancer: pro-cancer or pro-"host"? Cell Death Dis (2019) 10:650. doi: 10.1038/s41419-019-1883-8
26. van Dijk PC, Jager KJ, Zwinderman AH, Zoccali C, Dekker FW. The analysis of survival data in nephrology: Basic concepts and methods of cox regression. Kidney Int (2008) 74:705–9. doi: 10.1038/ki.2008.294
27. Tibshirani R. The lasso method for variable selection in the cox model. Stat Med (1997) 16:385–95. doi: 10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3
28. Wang H, Lengerich BJ, Aragam B, Xing EP. Precision lasso: accounting for correlations and linear dependencies in high-dimensional genomic data. Bioinformatics (2019) 35:1181–7. doi: 10.1093/bioinformatics/bty750
29. Scosyrev E, Glimm E. Power analysis for multivariable cox regression models. Stat Med (2019) 38:88–99. doi: 10.1002/sim.7964
30. Gene ontology consortium: Going forward. Nucleic Acids Res (2015) 43:D1049–1056. doi: 10.1093/nar/gku1179
31. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res (2017) 45:D353–d361. doi: 10.1093/nar/gkw1092
32. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods (2015) 12:453–7. doi: 10.1038/nmeth.3337
33. Meng Z, Ren D, Zhang K, Zhao J, Jin X, Wu H. Using ESTIMATE algorithm to establish an 8-mRNA signature prognosis prediction system and identify immunocyte infiltration-related genes in pancreatic adenocarcinoma. Aging (Albany NY) (2020) 12:5048–70. doi: 10.18632/aging.102931
34. Li T, Fu J, Wang B, Traugh N, Chen Q, Liu J, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res (2017) 77:e108–10. doi: 10.1158/0008-5472.Can-17-0307
35. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res (2020) 48:W509–w514. doi: 10.1093/nar/gkaa407
36. Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of drug sensitivity in cancer (GDSC): A resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res (2013) 41:D955–961. doi: 10.1093/nar/gks1111
37. Geeleher P, Cox N, Huang RS. pRRophetic: An r package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PloS One (2014) 9:e107468. doi: 10.1371/journal.pone.0107468
38. Han T, et al. Comprehensive analysis of inflammatory response-related genes, and prognosis and immune infiltration in patients with low-grade glioma. Front Pharmacol (2021) 12:748993. doi: 10.3389/fphar.2021.748993
39. Lin T, Cheng H, Liu D, Wen L, Kang J, Xu L, et al. A novel six autophagy-related genes signature associated with outcomes and immune microenvironment in lower-grade glioma. Front Genet (2021) 12:698284. doi: 10.3389/fgene.2021.698284
40. Lin JZ, Lin N, Zhao WJ. Identification and validation of a six-lncRNA prognostic signature with its ceRNA networks and candidate drugs in lower-grade gliomas. Genomics (2020) 112:2990–3002. doi: 10.1016/j.ygeno.2020.05.016
41. Li ZZ, Liu PF, An TT, Yang HC, Zhang W, Wang JX. Construction of a prognostic immune signature for lower grade glioma that can be recognized by MRI radiomics features to predict survival in LGG patients. Transl Oncol (2021) 14:101065. doi: 10.1016/j.tranon.2021.101065
42. Wu D, Wang X, Xue Y, Sun C, Zhang M. Identification of a novel immune-related lncRNA signature to predict prognostic outcome and therapeutic efficacy of LGG. J Integr Neurosci (2022) 21:55. doi: 10.31083/j.jin2102055
43. McIlwain DR, Berger T, Mak TW. Caspase functions in cell death and disease. Cold Spring Harb Perspect Biol (2013) 5:a008656. doi: 10.1101/cshperspect.a008656
44. Jiang M, Qi L, Li L, Li Y. The caspase-3/GSDME signal pathway as a switch between apoptosis and pyroptosis in cancer. Cell Death Discovery (2020) 6:112. doi: 10.1038/s41420-020-00349-0
45. Ding J, Wang K, Liu W, She Y, Sun Q, Shi J, et al. Pore-forming activity and structural autoinhibition of the gasdermin family. Nature (2016) 535:111–6. doi: 10.1038/nature18590
46. Tang R, Xu J, Zhang B, Liu J, Liang C, Hua J, et al. Ferroptosis, necroptosis, and pyroptosis in anticancer immunity. J Hematol Oncol (2020) 13:110. doi: 10.1186/s13045-020-00946-7
47. Kang R, Zeng L, Zhu S, Xie Y, Liu J, Wen Q, et al. Lipid peroxidation drives gasdermin d-mediated pyroptosis in lethal polymicrobial sepsis. Cell Host Microbe (2018) 24:97–108.e104. doi: 10.1016/j.chom.2018.05.009
48. Hou J, Zhao R, Xia W, Chang CW, You Y, Hsu JM, et al. PD-L1-mediated gasdermin c expression switches apoptosis to pyroptosis in cancer cells and facilitates tumour necrosis. Nat Cell Biol (2020) 22:1264–75. doi: 10.1038/s41556-020-0575-z
49. Feng X, Yu Y, He S, Cheng J, Gong Y, Zhang Z, et al. Dying glioma cells establish a proangiogenic microenvironment through a caspase 3 dependent mechanism. Cancer Lett (2017) 385:12–20. doi: 10.1016/j.canlet.2016.10.042
50. Fianco G, Mongiardi MP, Levi A, De Luca T, Desideri M, Trisciuoglio D, et al. Caspase-8 contributes to angiogenesis and chemotherapy resistance in glioblastoma. Elife (2017) 6. doi: 10.7554/eLife.22593
51. An HK, Chung KM, Park H, Hong J, Gim JE, Choi H, et al. CASP9 (caspase 9) is essential for autophagosome maturation through regulation of mitochondrial homeostasis. Autophagy (2020) 16:1598–617. doi: 10.1080/15548627.2019.1695398
52. Fernández-García V, González-Ramos S, Avendaño-Ortiz J, Martín-Sanz P, Delgado C, Castrillo A, et al. NOD1 splenic activation confers ferroptosis protection and reduces macrophage recruitment under pro-atherogenic conditions. BioMed Pharmacother (2022) 148:112769. doi: 10.1016/j.biopha.2022.112769
53. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol (2013) 14:1014–22. doi: 10.1038/ni.2703
54. Du T, Gao J, Li P, Wang Y, Qi Q, Liu X, et al. Pyroptosis, metabolism, and tumor immune microenvironment. Clin Transl Med (2021) 11:e492. doi: 10.1002/ctm2.492
55. Brydges SD, Broderick L, McGeough MD, Pena CA, Mueller JL, Hoffman HM. Divergence of IL-1, IL-18, and cell death in NLRP3 inflammasomopathies. J Clin Invest (2013) 123:4695–705. doi: 10.1172/jci71543
56. Zhou T, Damsky W, Weizman OE, McGeary MK, Hartmann KP, Rosen CE, et al. IL-18BP is a secreted immune checkpoint and barrier to IL-18 immunotherapy. Nature (2020) 583:609–14. doi: 10.1038/s41586-020-2422-6
57. Kaplanski G. Interleukin-18: Biological properties and role in disease pathogenesis. Immunol Rev (2018) 281:138–53. doi: 10.1111/imr.12616
58. Fabbi M, Carbotti G, Ferrini S. Context-dependent role of IL-18 in cancer biology and counter-regulation by IL-18BP. J Leukoc Biol (2015) 97:665–75. doi: 10.1189/jlb.5RU0714-360RR
59. Jiang Y, Han S, Cheng W, Wang Z, Wu A. NFAT1-regulated IL6 signalling contributes to aggressive phenotypes of glioma. Cell Commun Signal (2017) 15:54. doi: 10.1186/s12964-017-0210-1
60. Liu B, Tan J, Wang X, Liu X. Identification of recurrent risk-related genes and establishment of support vector machine prediction model for gastric cancer. Neoplasma (2018) 65:360–6. doi: 10.4149/neo_2018_170507N326
61. Wang Y, Gao W, Shi X, Ding J, Liu W, He H, et al. Chemotherapy drugs induce pyroptosis through caspase-3 cleavage of a gasdermin. Nature (2017) 547:99–103. doi: 10.1038/nature22393
62. Gieryng A, Pszczolkowska D, Walentynowicz KA, Rajan WD, Kaminska B. Immune microenvironment of gliomas. Lab Invest (2017) 97:498–518. doi: 10.1038/labinvest.2017.19
63. Kane JR, Zhao J, Tsujiuchi T, Laffleur B, Arrieta VA, Mahajan A, et al. CD8(+) T-cell-Mediated immunoediting influences genomic evolution and immune evasion in murine gliomas. Clin Cancer Res (2020) 26:4390–401. doi: 10.1158/1078-0432.Ccr-19-3104
64. Zhang Y, Zhang Z. The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications. Cell Mol Immunol (2020) 17:807–21. doi: 10.1038/s41423-020-0488-6
65. Xu S, Tang L, Li X, Fan F, Liu Z. Immunotherapy for glioma: Current management and future application. Cancer Lett (2020) 476:1–12. doi: 10.1016/j.canlet.2020.02.002
66. Liu Y, Fang Y, Chen X, Wang Z, Liang X, Zhang T, et al. Gasdermin e-mediated target cell pyroptosis by CAR T cells triggers cytokine release syndrome. Sci Immunol (2020) 5(43):eaax7969. doi: 10.1126/sciimmunol.aax7969
Keywords: LGG, pyroptosis, prognostic model, bioinformatics, immunity
Citation: Liu H and Tao T (2022) Prognosis and immune features of pyroptosis-related RNA patterns in low-grade glioma. Front. Oncol. 12:1015850. doi: 10.3389/fonc.2022.1015850
Received: 10 August 2022; Accepted: 08 November 2022;
Published: 20 December 2022.
Edited by:
Tian Li, Independent Researcher, Xi’an, ChinaReviewed by:
Siyuan Song, Affiliated Hospital of Nanjing University of Chinese Medicine, ChinaNan Zhang, Harbin Medical University, China
Zewei Tu, The Second Affiliated Hospital of Nanchang University, China
Copyright © 2022 Liu and Tao. 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: Tao Tao, Taotao2017@vip.163.com; Hanzhang Liu, ntdxlhz@ntu.edu.cn