key: cord-0910709-w315prob
authors: Shen, Lin; Li, Yanyan; Li, Na; Zhao, Yajie; Zhou, Qin; Shen, Liangfang; Li, Zhanzhan
title: Integrative analysis reveals the functional implications and clinical relevance of pyroptosis in low-grade glioma
date: 2022-03-16
journal: Sci Rep
DOI: 10.1038/s41598-022-08619-w
sha: 087cc10432582b7a34a17b8d1c368e226fc194f5
doc_id: 910709
cord_uid: w315prob

Using the Chinese Glioma Genome Atlas (training dataset) and The Cancer Genome Atlas (validation dataset), we found that low-grade gliomas can be divided into two molecular subclasses based on 30 pyroptosis genes. Cluster 1 presented higher immune cell and immune function scores and poorer prognosis than Cluster 2. We established a prognostic model based on 10 pyroptosis genes; the model could predict overall survival in glioma and was well validated in an independent dataset. The high-risk group had relatively higher immune cell and immune function scores and lower DNA methylation levels in pyroptosis genes than the low-risk group. There were no marked differences in pyroptosis gene alterations between the high- and low-risk groups. The competing endogenous RNA (ceRNA) regulatory network uncovered the lncRNA–miRNA–mRNA regulation patterns of the different risk groups in low-grade glioma. Five pairs of target genes and drugs were identified. In vitro, CASP8 silencing inhibited the migration and invasion of glioma cells. The expression of pyroptosis genes can reflect the molecular biological and clinical features of low-grade glioma subclasses. The developed prognostic model can predict overall survival and distinguish molecular alterations in patients. Our integrated analyses could provide valuable guidelines for improving risk management and therapy for low-grade glioma patients.

We identified the molecular subclasses using the consensus clustering method, which can determine the optimal cluster number via a cumulative distribution function. We also used a t-distributed stochastic neighbour embedding-based approach to validate the clustering consensus. We used the GSVA R package to determine the potential pathway enrichment by calculating the enrichment scores of each sample.

Development and validation of the prognostic model. We established a prognostic model of overall survival in the CGGA training dataset using pyroptosis genes. LASSO regression was performed to identify the potential genes in the prognostic signature, and a Cox regression model was developed. Then, we calculated the risk score of each sample using the following equation: risk score = coef1*gene 1 expression + coef 2 * gene 2 expression + …coef n *gene n expression. The patients in the training and validation cohorts were classified into the high-risk (risk score > median risk score) and low-risk (risk score ≤ median risk score) groups. We calculated the 1-year, 2-year, and 3-year areas under the curves (AUCs) to evaluate the predictive ability of the prognostic model. Finally, we developed a quantized assignment operator to calculate the 1-year, 3-year, and 5-year survival probabilities of individual patients based on the prognostic model, which were evaluated by consensus calibration analysis.

Gene Ontology and KEGG pathway analyses were carried out using the TCGA and CGGA datasets separately 14, 15 . To calculate the enrichment scores of stem cell-like properties (RNAss and DNAss) and the TME (stromal score, immune score, and ESTIMATE score), we carried out single-sample gene set enrichment analysis in the TCGA dataset. In addition, the immune-related cell and function scores were calculated for each sample (downloaded from https:// www. gsea-msigdb. org/). Somatic copy number alteration, mutation, and DNA methylation analysis. We further compared the copy number alteration frequencies and DNA methylation levels between the high-and low-risk groups using the "limma" R package.

Construction of the ceRNA network and analysis of drug sensitivity. To investigate the potential lncRNA-miRNA-mRNA regulatory network, we identified the differentially expressed lncRNAs, miRNAs, and mRNAs between the high-and low-risk groups. Then, we constructed the lncRNA-miRNA-mRNA regulatory network using Cytoscape (version 3.8.2). Based on the correlation coefficient (|R| > 0.25, P < 0.05), we identified the potential small molecule compounds related to the pyroptosis genes included in the prognostic model.

In vitro experimental verification. We further performed quantitative polymerase chain reaction, Western blot analysis, migration assays, invasion assays, and clonogenic assays to explore the effects of silencing CASP3 on cell phenotypes. The descriptions of the in vitro experimental processes are provided in Supplementary material 1.

Statistical analysis. We adopted the log-rank test to compare the overall survival curves for the two groups.

Univariate and multivariate Cox regression analysis was conducted to assess the correlation of the risk score with prognosis in low-grade glioma, and hazard ratios (HRs) and confidence intervals (CIs) were calculated. We performed all statistical analyses using R software (version 4.0). A P value < 0.05 was considered significant unless specifically defined otherwise.

Identification of molecular subclasses of low-grade glioma. The flow diagram of the integrated analysis in the study is shown in Fig. 1A . We obtained 630 samples with mRNA expression data from the CGGA dataset, and 30 pyroptosis genes were extracted (MAD > 0.5). Supplementary material 2 shows the gene symbols and Entrez IDs of the pyroptosis genes used in this study. Using the STRING database, we constructed a proteinprotein interaction (PPI) network of all pyroptosis genes, and the PPI network showed that several pyroptosis genes, including CASP1, CASP3, CASP4, CASP8, NLRP1, NLRP3, and NLRC4, were connected to more nodes than were other genes (Fig. 1B www.nature.com/scientificreports/ The expression correlations among these pyroptosis genes are presented in Fig. 1C (red: positive correlations; dark colour: negative correlations). The consensus clustering analysis indicated that all samples of low-grade glioma could be divided into two subclasses. The consensus matrix exhibited a relatively sharp and clear boundary, which indicated stable and robust clustering (Fig. 1D) . We further found a two-dimensional t-sensitivity distribution that supported the subtype clustering (Fig. 1E ). The specific subclass information of each sample is provided in Supplementary material 2: Table S2 . Finally, we compared the overall survival curves for the two subclasses using Kaplan-Meier analysis. Our results indicated that the median overall survival time was significantly longer in Cluster 2 than in Cluster 1. Kaplan-Meier analysis indicated that the median survival time was significantly shorter in Cluster 2 than in Cluster 1 (MST: 9.62 vs. 4.35 years, P < 0.001; Fig. 1F ). This result showed that the two subclasses had markedly different prognostic patterns.

Correlations of the molecular subclasses with pyroptosis genes. To explore the signalling pathway enrichment in the two subclasses, we performed gene set variation analysis (GSVA) by transforming the expression data from a gene-by-sample matrix to a gene set by subclass matrix for both subclasses. The results indicated that the two subclasses had different pathway enrichment profiles. Unlike Cluster 2, Cluster 1 had 78 kinds of significantly different signalling pathways (Supplementary material 2: Table S3 ). The upregulated pathways mainly included type 1 diabetes mellitus, intestinal immune network for IGA production, asthma, autoimmune thyroid disease, graft versus host disease, allograft rejection, and ribosome. In addition, the DNA replication, mismatch repair, ECM receptor interaction, cell cycle, and p53 signalling pathways, and some metabolic pathways were significantly enriched ( Fig. 2A) .

Clinical characteristics and pyroptosis gene expression patterns of the molecular subclasses. We investigated the correlations of the two subclasses with clinical parameters (Fig. 2B ). Compared with patients in Cluster 1, who had poor prognoses, patients from Cluster 2 tended to have a lower WHO stage (P < 0.001), lower rate of chemotherapy (P < 0.05), 1p19q non-codeletion status (P < 0.001), and IDH mutation status (P < 0.001). No significant differences were observed in PRS type, sex, age, or radiotherapy status (P > 0.05). Regarding the expression levels of pyroptosis genes, significant differential expression was observed between the two subclasses for all genes except AIM2 and TIRAP. All of these differentially expressed genes were upregulated in Cluster 2 (Fig. 2B) . We also compared the differences in pyroptosis gene expression by WHO stage, IDH mutation status, 1p19q status, and chemotherapy history. Compared with the WHO II group, the WHO III group had sixteen upregulated genes (CASP1, CASP3, CASP4, CASP5, CASP6, CASP8, GSDMB, GSDMD, IL18, IL6, NLRC4, NOD1, NOD2, PLCG1, PRKACA, PYCARD) (Supplementary material 3: Fig. S1A ). For IDH mutation status, 21 differentially expressed genes (DEGs) were found (Supplementary material 3: Fig. S1B ). A total of 27 pyroptosis DEGs were found for 1p191 status (Supplementary material 3: Fig. S1C ). Thirteen DEGs were found for chemotherapy history (Supplementary material 3: Fig. S1D) .

We further performed differential expression analysis between Cluster 1 and Cluster 2. A total of 1053 DEGs were found; 67 genes were upregulated and 986 were downregulated in Cluster 2 (Supplementary material 2: Table S4 ). We performed GO and KEGG enrichment analyses using these DEGs (Supplementary material 2: Table S5 and Table S6 ). A total of 1033 differential functional terms were enriched, namely, 856 biological process, 116 cellular component and 61 molecular function terms. The top 30 enriched terms are presented in Supplementary material 3: Fig. S2 . Most of these functions were associated with immunity. In addition, we identified 56 significant pathways by KEGG analysis (Supplementary material 3: Fig. S3 ), and the top five were phagosome, Staphylococcus aureus infection, coronavirus disease, antigen processing and presentation, and tuberculosis.

To explore tumour heterogeneity between the two subclasses, we compared immune cell and immune function differences. Compared with Cluster 2, Cluster 1 had higher aDC, B cell, CD8+ T cell, DC, iDC, macrophage, mast cell, neutrophil, NK cell, pDC, T helper cell, Tfh cell, Th2 cell, Tfh cell, TIL, and Treg cell levels (all P < 0.001, Supplementary material 3: Fig. S4A ). Similarly, Cluster 1 had higher immune function scores, including APC coinhibition, APC costimulation, CCR, checkpoint, cytolytic activity, HLA, inflammation promotion, MHC class I, parainflammation, T cell coinhibition, type I IFN response and type II IFN response scores, than Cluster 2 (all P < 0.001, Supplementary material 3: Fig. S4B ).

Establishment of the pyroptosis-related prognostic model. Initially, we performed univariate Cox regression analysis to identify prognosis-related pyroptosis genes in the CGGA cohort (Supplementary material 3: Fig. S5A ). In total, 21 pyroptosis genes were found to be associated with the overall survival of low-grade glioma patients. Kaplan-Meier analysis indicated that high expression of all 21 genes (CASP1, CASP3, CASP4, CASP5, CASP6, CASP8, GSDMA, GSDMB, GSDMC, GSDMD, IL18, IL1B, IL6, NLRC4, NLRP3, NOD1, NOD2, PLCG1, PRKACA, PYCARD, SCAF11) was associated with poor OS in low-grade glioma. Furthermore, using LASSO regression in the CGGA training cohort, we identified 10 pyroptosis genes (CASP5, CASP6, CASP8, GSDMC, IL1B, IL6, NLRP3, NOD2, PLCG1, SCAF11) to establish the prognostic model (Supplementary material 3: Fig. S5B ,C). We calculated the risk score for each sample using the regression coefficients of the 10 genes (Supplementary material 2: Table S7 ). Patients with risk scores higher than the median value were classified into the high-risk group, and the other patients were classified into the low-risk group. Compared with the patients in the low-risk group, the patients in the high-risk group tended to have advanced WHO stage (P < 0.001), recurrence (P < 0.001), a history of chemotherapy (P < 0.001), IDH mutation status (P < 0.001), and 1p19q codeletion status (P < 0.001) (Supplementary material 3: Fig. S5D ). Finally, we found that glioma patients in Cluster 1, with Fig. S6 ). www.nature.com/scientificreports/ Kaplan-Meier analysis indicated that the high-risk group had a significantly worse OS than the low-risk group (Fig. 3A,B) . Univariate Cox regression analysis indicated that the risk score was positively associated with poor OS (HR = 3. 878, 95% CI: 3.012-4.992, P < 0.001; Fig. 3C ). Multivariate Cox regression analysis showed that an elevated risk score was an independent predictor of unfavourable prognosis in low-grade glioma (HR = 2.419, 95% CI: 1.823-3.211, P < 0.001; Fig. 3D ). In addition, PRS type (recurrence: HR = 1.901, 95% CI: 1.419-2.546, P < 0.001) and WHO III grade (HR = 2.188, 95% CI: 1.581-3.028, P < 0.001) were positively associated with poor OS, and 1p19q non-co-deletion status (HR = 0.450, 95% CI: 0.295-0.686, P < 0.001) was negatively associated with OS. PCA also indicated that the high-and low-risk groups showed two markedly different distributions (Fig. 3E) . Time-dependent receiver operating characteristic analysis was performed to evaluate the predictive ability of the prognostic model. Our results showed that the AUCs at 1 year, 2 years, and 3 years were 0.670, 0.734 and 0.723 (Fig. 3F) , respectively. We further compared the OS status among the WHO stage, sex, age, IDH status, 1p19q codeletion status, radiotherapy history, and chemotherapy history subgroups. All results indicated that the OS in the high-risk group was still poorer than that in the low-risk group (Supplementary material 3: Fig. S7 , all P < 0.001).

Validation of the prognostic model. To further validate the pyroptosis gene model, we also calculated the risk scores of low-grade glioma patients in the TCGA cohort. Kaplan-Meier analysis indicated a significant correlation with poor OS in the high-risk group (Fig. 4A,B) . Univariate Cox regression analysis showed that an increased risk score was significantly associated with poor OS in the TCGA cohort (HR = 2.011, 95% CI: 1.714-2.360, P < 0.001; Fig. 4C ). An elevated risk score was also an independent prognostic indicator in multivariate Cox regression analysis (HR = 1.837, 95% CI: 1.538-2.193, P < 0.001; Fig. 4D ). PCA also validated the high-and low-risk distribution patterns of all patients (Fig. 4E) . Furthermore, the AUCs of the risk score were 0.832 at 1 year, 0.793 at 2 years, and 0.804 at 3 years (Fig. 4F ).

Clinical application of the prognostic model. To further evaluate the clinical predictive value of the prognostic model, we developed a nomogram based on multivariate Cox regression analysis that included significant clinical parameters in the CGGA dataset (Fig. 5A) . The calibration curves indicated that the clinical nomogram could precisely predict the 1-year, 3-year and 5-year OS of glioma patients (C-index = 0.799, Fig. 5B-D) . The predictive accuracy of this nomogram was well validated in the TCGA dataset (C-index = 0.841, Fig. 5E-G) .

To explore the potential biological functions that affect the overall survival of low-grade glioma patients, we performed GO term and KEGG pathway enrichment analyses. We first identified DEGs between the high-and low-risk groups and then annotated the functions of the DEGs in terms of biological processes, cellular components, and molecular functions. We found 1571 DEGs in the CGGA cohort (Supplementary material 2: Table S8 ) and 609 DEGs in the TCGA cohort (Supplementary material 2: Table S9 ). The GO term and KEGG pathway enrichment analyses indicated that the CGGA and TCGA cohorts shared some enriched terms and pathways, such as extracellular matrix organization, extracellular structure organization, immune response, ECM-receptor interaction, cell adhesion molecules, PI3K-Akt signalling pathway, and Epstein-Barr virus infection (Fig. 6A-D) .

We further compared the differences in immune cells and immune functions between the high-and low-risk groups in the CGGA cohort (Fig. 6E ,G) and TCGA cohort (Fig. 6F,H) . As shown in the box plots, the immune cell score showed a similar trend in the CGGA and TCGA cohorts. Most immune cell scores showed a tendency to be increased in the high-risk group. The differences in immune function between the different risk groups were similar in the two datasets (all P < 0.001). All immune function scores were significantly increased in the high-risk group. The 10 pyroptosis genes included in the prognostic model showed significant differences among the four immune subtypes (C3, C4, C5, C6) (Supplementary material 3: Fig. S8 ). We finally investigated the correlations of the signature genes with cancer stem cell-like properties (RNAss and DNAss) and the TME (stromal score, immune score, and ESTIMATE score). We found that all ten genes were negatively associated with the RNAss and the expression of CASP5, CASP6, CASP8, GSDMC, IL18, IL6m, NLRP3, and NOD2 and that all genes except PLCG1 and SCAF11 showed positive associations with the DNAss, stromal score, immune score, and ESTIMATE score. PLCGa and SCAF11 were negatively associated with the stromal score, immune score, and ESTIMATE score. SCAF11 was the only gene that was negatively associated with the DNAss (Supplementary material 3: Fig. S9 ).

Profiling of pyroptosis gene alterations. Molecular alterations in pyroptosis-related genes were also evaluated based on the high-and low-risk groups in the TCGA dataset. PLCG1 was the only gene with alteration in the low-risk group, and NLRP2, GSDMC and NLRP3 were the genes with alterations in the high-risk group. All altered genes had an alteration frequency of less than or equal to 2% (Fig. 7A,B) . Somatic copy number alteration analysis indicated significant differences among the pyroptosis genes. Among these genes, the copy number variation in IL18, NLRP6, IL6, CASP5, CASP4, CASP1, and NOD1 was significantly increased and that of TNF, NLRP2, NLRP7, CASP3, and CASP6 was significantly decreased in the high-risk group (Fig. 7C ). In addition, differences in the DNA methylation levels of the pyroptosis genes between the high-and low-risk groups were assessed. The results showed that the overall DNA methylation levels of all pyroptosis genes were significantly higher in the high-risk group than in the low-risk group (Fig. 7D ).

Construction of a ceRNA network based on the prognostic signature. To explore differences in the lncRNA-miRNA-mRNA regulatory network between the high-and low-risk groups, we constructed a ceRNA network based on the differentially expressed mRNAs, lncRNAs and miRNAs between the high-and low-risk Table S13 ). There were 5 pairs with a correlation coefficient > 0.5 or < − 0.5. The NOD2-isotretinoin, IL1B-rebimastat, NLRP3-rebimastat, and NOD2-elesclomol pairs showed drug sensitivity. The GSDMC-ixazomib citrate pair showed drug resistance (Fig. 9 ).

We used the following criteria to select the validated genes: (1) genes differentially expressed between tumour and normal tissue; (2) genes included in the prognostic model; and (3) genes that were rarely reported previously (Supplementary material 3: Fig. S13 ). We first detected the expression of CASP8 in glioma cell lines using Western blot analysis (Fig. 10A,B) and found that CASP8 expression was highest in LN299 cells. We generated CASP8-si LN229, H4 and U87 glioma cells. qPCR analysis indicated that the mRNA expression of CASP8 was significantly downregulated in CASP8-si H4 and LN229 cells (Fig. 10C) . Furthermore, silencing CASP8 inhibited the migration and invasion of glioma cells (Fig. 10D-G) . The clonogenic assay showed that CASP8 silencing also inhibited cell growth (Fig. 10F-G) . These results suggested that CASP8 silencing suppressed the malignant progression of glioma cells.

In the present study, we found that low-grade glioma can be categorized into two molecular subtypes based on the expression levels of 30 pyroptosis genes, and markedly different survival outcomes were observed between the two subclasses. The traditional classification system based on clinical parameters have been updated several times. However, the traditional classification based on histology also has some limitations. The primary limitation is mainly due to interobserver heterogeneity 2 . It was previously reported that there was only approximately 50% agreement between different neuropathologists in reviewing cases, specifically astrocytoma and www.nature.com/scientificreports/ oligodendroglioma cases 16 . Genetic studies and the associated findings enable us to better understand how tumours differ in clinical outcomes and molecular patterns and to facilitate the effective treatment of tumour subtypes based on gene expression features. We also found that some immune-related pathways were enriched in the poor prognosis group. Furthermore, the poor prognosis group showed elevated immune cell levels and immune-related functions, such as immune checkpoint activity, inflammation promotion, and parainflammation. It has been reported that pyroptosis represents an antitumour immune function in cancers, which means that pyroptosis can induce inflammation, triggering robust antitumour immunity and synergizing with immune checkpoint blockade 17 . Moreover, some key pathways, such as oxidative phosphorylation, PPAR signalling pathway, primary immunodeficiency, and citrate cycle (TCA cycle), which have been reported to be involved in glioma progression, were highly enriched in the poor prognosis group [18] [19] [20] [21] . These results suggested that pyroptosis genes can accurately differentiate low-grade glioma patients into two-dimensional distributions.

We developed a prognostic model based on 10 pyroptosis genes (CASP5, CASP6, CASP8, GSDMC, IL1B, IL6, NLRP3, NOD2, PLCG1, SCAF11). This prognostic model was well validated in an independent cohort, and www.nature.com/scientificreports/ the 1-year, 2-year, and 3-year AUCs were 0.670, 0.734 and 0.723, respectively. This predictive ability showed the moderate discernibility of the model. Furthermore, using the clinical parameters and risk scores obtained from the 10 genes included in the model, we developed a nomogram for estimating an individual's overall survival probability. The results from the training dataset and validation dataset showed high consistency. Our results suggested that the established prognostic model has clinical value. The cell death via pyroptosis involves two biological mechanisms. The classical pyroptosis mechanism involves the assembly of inflammasomes. Inflammasomes are macromolecular protein complexes that are necessary for inflammation in the cytoplasm and recognize danger signal molecules such as those released by bacteria and viruses 22 . Inflammasomes are mainly composed of pattern recognition receptors (PRRs), apoptosis-associated speck-like protein (ASC) and pro-caspase-1 precursors 23 . PRRs are receptor proteins responsible for recognizing different signal stimuli in cells. They are mainly composed of nucleotide-binding oligomerization domainlike receptor protein (NLRP) 1, NLRP3, nucleotide-binding oligomerization domain-like receptor protein C4 (NLRC4), absent in melanoma 2 (AIM2) and other components 24 . ASC is an adaptor protein that is mainly composed of an N-terminal pyrin domain (PYD) and a C-terminal caspase activation and recruitment domain (CARD) 25 . Procaspase-1 is an effector molecule that can specifically cleave GSDMD after activation. After the danger signal sensor NLR1, NLRP3 or AIM2 recognizes a danger signal molecule, the N-terminal PYD interacts with the N-terminal PYD of the adaptor protein. ASC then recruits Caspase-1 through interaction with the CARD domain to complete the assembly of the inflammasome 26 . This method of cell death mediated by Caspase-1 is called the classical pathway of pyroptosis 27 . The nonclassical pyroptosis pathway is mainly mediated by Caspases-4, -5, and -11. After cells are stimulated by bacterial LPS, Caspases-4, -5, and -11 directly bind to bacterial LPS and are activated 28 . Activated Caspases-4, -5, and -11 specifically cleave GSDMD and alleviate the intramolecular inhibition 29 . The interaction of the GSDMD-N-terminus with cell membrane phospholipids causes pore formation in the cell membrane, cell swelling and cell rupture and induces pyroptosis; the GSDMD-N-terminus can also activate Caspase-1 by activating the NLRP3 inflammasome 30 . Activated Caspase-1 stimulates the maturation of the IL-18 and IL-1β precursors, and IL-18 and IL-1β are secreted into the extracellular space and amplify the inflammatory response. Yang et al. found that in the nonclassical Caspase-11-dependent pathway, gap junction protein-1 (Pannexin-1) can be cleaved and that cleavage of Pannexin-1 can activate its own channel and release ATP, which induces pyroptosis 31 . Lamkanfi et al. found that in the nonclassical Caspase-11-dependent pathway, Pannexin-1 cleavage can also activate the NLRP3 inflammasome, which in turn activates Caspase-1 and induces pyroptosis 32 . In our study, we found that the prognostic model included CASP5, CASP6, and CASP8 along with IL18, NLRP3, and GSDMC.

These results indicated that two biological mechanisms are involved in the development and progression of low-grade glioma. Further research is required to illustrate the specific molecular mechanisms.

We divided low-grade glioma patients into high-and low-risk groups using the estimated risk score and investigated the differences in genomic patterns and clinical features between the two risk groups. Our results indicated that the enriched functions and pathways were markedly different between the high-and low-risk groups. However, the results in two independents datasets showed that these risk groups shared some similar functional enrichment and signalling pathways, such as extracellular matrix organization, focal adhesion, ECMreceptor interaction, and GABAergic synapse. Immune infiltration is very closely associated with the progression of tumours 33 . Our results showed that the immune cell and immune function scores were significantly elevated in the high-risk group and that the poor prognosis group also had a higher risk score than the favourable prognosis group. These results showed that the molecular subtypes and risk classifications were stable and robust. Subsequently, we explored the differences in alterations, CNVs, and DNA methylation levels of pyroptosis genes between the high-and low-risk groups. We found that these genomic patterns were significantly different between the two risk groups, providing some molecular basis for the changes corresponding to poor clinical outcomes. Finally, we constructed a ceRNA regulatory network that identified several key lncRNA-miRNA-mRNA regulatory axes. The ceRNA network identified several key lncRNA-miRNA-mRNA regulatory axes: FAM181A-AS1-miR-21-(MAP2K3, JAG1, TGFBI, and FAM24A), CRNDE-miR-155-(DPYA1L1, NAMPT, TRIP13, IKBIP, SPL1, EGFR, WEE1), NEAT1-miR-200a-(EPHA2, DPY19L1, PTPRD, LATS2, and ELAVL2). Survival analysis further suggested regulatory correlations: elevated expression of FAM18A-AS1 and miR-21 was associated with poor prognosis in low-grade glioma, and high expression of EGFR, WEE1, MAP2K3, JA1, and EPHA2 was associated with poor prognosis. Previous experiments have reported the promoting role of miR-21 in glioma 34, 35 , and upregulation of EGFR and EPHA2 is associated with the development and progression of glioma [36] [37] [38] [39] . Drug sensitivity analysis indicated that NOD2, ELANE, NLRP3, CASP3, and PRKACA showed sensitivity to small molecule drugs, and that PRKACA and IL6 showed resistance to some compounds. A previous study reported that the migration ability of glioma cells was reduced after inhibition of the NLRP3 inflammasome by betahydroxybutyrate 40 . These results may provide some guidelines for clinical therapy.

In conclusion, our results indicate that pyroptosis genes can be used to classify low-grade glioma patients into two subgroups. The established prognostic model based on 10 pyroptosis genes can not only predict the prognosis of low-grade glioma patients but also reflect the molecular alterations, immune status, and stem cell-like properties of the high-and low-risk groups. Classification based on the risk score of the prognostic signature revealed a lncRNA-miRNA-mRNA regulatory network. The correlations of the signature genes with drug sensitivity may provide a rationale for clinical treatment. Finally, our study uncovers the biological function and clinical relevance of pyroptosis in the occurrence and progression of low-grade glioma and may offer some risk management and treatment strategies. Future research should focus on molecular mechanisms and therapeutic targets.

Glioma subclassifications and their clinical significance

Low-grade glioma radiotherapy treatment and trials

Management of pediatric low-grade glioma

Molecular markers and targeted therapy in pediatric low-grade glioma

Programmed death, cells on the last train to glory

Coordinated host responses during pyroptosis: Caspase-1-dependent lysosome exocytosis and inflammatory cytokine maturation

Caspase-1-dependent pore formation during pyroptosis leads to osmotic lysis of infected host macrophages

Pillars article: Caspase-1-dependent pore formation during pyroptosis leads to osmotic lysis of infected host macrophages

Diverging inflammasome signals in tumorigenesis and potential targeting

The role of pyroptosis in cancer: Pro-cancer or pro-"host

AIM2 Inflammasome activation and regulation: A structural perspective

Regulation of inflammasome activation

KEGG as a reference resource for gene and protein annotation

Kyoto Encyclopedia of Genes and Genomes

Improving diagnostic accuracy and interobserver concordance in the classification and grading of primary gliomas

A bioorthogonal system reveals antitumour immune function of pyroptosis

Lactate reprograms energy and lipid metabolism in glucose-deprived oxidative glioma stem cells

Correlation between TSP-1, TGF-beta and PPAR-gamma expression levels and glioma microvascular density

Clinical characteristics and prognostic analysis of glioma in human immunodeficiency virus-infected patients

Emerging roles and therapeutic interventions of aerobic glycolysis in glioma

Inflammasome activation triggers blood clotting and host death through pyroptosis

Inflammasome-associated cell death: Pyroptosis, apoptosis, and physiological implications

Pattern recognition receptors and inflammation

The inflammasome adapter ASC assembles into filaments with integral participation of its two death domains PYD and CARD

Unified polymerization mechanism for the assembly of ASC-dependent inflammasomes

Analyzing caspase-1 activation during legionella pneumophila infection in macrophages

Human caspase-4 detects tetra-acylated LPS and cytosolic Francisella and functions differently from murine caspase-11

Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death

Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling

Caspase-11 requires the pannexin-1 channel and the purinergic P2X7 pore to mediate pyroptosis and endotoxic shock

Mechanisms and functions of inflammasomes

Tumor immune microenvironment and immune checkpoint inhibitors in esophageal squamous cell carcinoma

MiR-21 is linked to glioma angiogenesis: A co-localization study

MiR-21 enhanced glioma cells resistance to carmustine via decreasing Spry2 expression

vascular endothelial growth factor, and vascular endothelial growth factor correlate with adverse outcomes and poor survival in patients with glioma

Updated insights on EGFR signaling pathways in glioma

Microglial-stimulation of glioma invasion involves the EGFR ligand amphiregulin

Involvement of EGFR, ERK-1,2 and AKT-1,2 activity on human glioma cell growth. Asian Pac

The beta-hydroxybutyrate suppresses the migration of glioma cells by inhibition of NLRP3 inflammasome

The authors declare no competing interests.

Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-08619-w.Correspondence and requests for materials should be addressed to Z.L.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.