- 1Department of Clinical Pharmacy, 920th Hospital of Joint Logistics Support Force of People’s Liberation Army, Kunming, China
- 2College of Pharmacy, Dali University, Dali, China
- 3Department of Pharmacy, The Second People’s Hospital of Quzhou Zhejiang, Quzhou, China
- 4Department of Cardiothoracic Surgery, 920th Hospital of Joint Logistics Support Force of People’s Liberation Army, Kunming, China
- 5Research Center of Clinical Pharmacology, Yunnan Provincial Hospital of Traditional Chinese Medicine, Kunming, China
Background: Pyroptosis is a novel inflammatory form of programmed cell death and a prospective target for cancer therapy. Nevertheless, little is known about the association between pyroptosis-related genes (PRGs) and acute myeloid leukemia (AML) prognosis. Herein, we systematically investigated the specific functions and clinical prognostic value of multiple PRGs in AML.
Methods: Univariate and LASSO Cox regression analyses based on TCGA and GTEx databases were used to generate the PRG signature, whose predictive efficacy of survival was evaluated using survival analysis, ROC, univariate and multivariate Cox analyses as well as subgroup analysis. The BeatAML cohort was used for data validation. The association between risk score and immune cell infiltration, HLA, immune checkpoints, cancer stem cell (CSC), tumor mutation burden (TMB), and therapeutic drug sensitivity were also analyzed.
Results: Six -PRG signatures, namely, CASP3, ELANE, GSDMA, NOD1, PYCARD, and VDR were generated. The high-risk score represented a poorer prognosis and the PRG risk score was also validated as an independent predictor of prognosis. A nomogram including the cytogenetic risk, age, and risk score was constructed for accurate prediction of 1-, 3-, and 5-year survival probabilities. Meanwhile, this risk score was significantly associated with the tumor immune microenvironment (TIME). A high-risk score is characterized by high immune cell infiltration, HLA, and immune checkpoints, as well as low CSC and TMB. In addition, patients with low-risk scores presented significantly lower IC50 values for ATRA, cytarabine, midostaurin, doxorubicin, and etoposide.
Conclusion: Our findings might contribute to further understanding of PRGs in the prognosis and development of AML and provide novel and reliable biomarkers for its precise prevention and treatment.
Introduction
Acute myeloid leukemia (AML) is the second most common type of leukemia diagnosed in adults and accounted for approximately 30% of all adult leukemia cases (Perna et al., 2017), which is primarily characterized by heterogeneity of molecular abnormalities and aberrant accumulation of immature myeloid progenitors in bone marrow and peripheral blood owing to impaired differentiation of hematopoietic progenitors (Sperlazza et al., 2015). The invasive infiltration of AML is mainly represented by a malignant extramedullary infiltration which involves the skin, lymph nodes, liver, spleen, and even central nervous system, and shows poorer prognosis in clinic (Stölzel et al., 2011). In recent years, although the survival and prognosis of AML patients were relatively prolonged with the development and application of molecular targeted therapy and combination therapy in clinical practice, the 5-year survival rate remained at 27% for AML patients over 20 years old and even 69% for patients younger than 20 partially due to the absence of reliable prognostic biomarkers (Ball and Stein, 2019; Schlenk et al., 2019; Zheng et al., 2020). Thus, identifying novel and effective prognostic biomarkers is crucial for improving the AML prognosis and better understanding the pathogenesis of AML.
Recently, much attention was paid to the influences of the tumor microenvironment (TME) on tumorigenesis and development, and alterations in TME components were found at every stage of malignant development in almost all carcinomas (Wu and Dai, 2017; Arneth, 2019). As an important component of TME, the tumor immune microenvironment (TIME) was also found to play prominent roles in tumor cell proliferation, invasion, and metastasis (Fu et al., 2019). Numerous studies showed that TIME was a key determinant of diagnosis and therapeutic response in tumor patients (Xu B. et al., 2021; Yang et al., 2021; Zeng et al., 2021; Wang et al., 2022). However, TIME is complex and variable mainly because of the multiple interaction networks among tumor, immune, stromal, and mesenchymal cells along with various soluble factors and changes in extracellular matrix (ECM) components as well (Wu and Dai, 2017; Xu B. et al., 2021). Therefore, identification of certain potential biomarkers related to TIME would eventually contribute to better understanding of tumor development and further identification of candidate therapeutic targets.
As a highly specific inflammatory programmed cell death, pyroptosis was reported to be significantly cross-correlated with TIME according to previous studies (Orning et al., 2019; Erkes et al., 2020). When persistent inflammation is present, initial activation and assembly of inflammasome started within host cells (Balahura et al., 2020). Subsequently, the caspase was further activated and produced inflammatory cytokines, eventually resulting in pyroptotic cell death (Broz et al., 2010; Moujalled et al., 2021). In recent years, a growing number of studies demonstrated that pyroptosis played crucial roles in pathogenesis and progression of various types of cancers including AML. It was reported that activation of NLR family pyrin domain containing 1 b (NLRP1b) by small molecule inhibitors of serine dipeptidase 8/9 (DPP8/9) induced caspase-1 dependent pyroptosis, which, in turn, suppressed the development of AML (Johnson et al., 2018). Meanwhile, pyridoxine was also found to induce death of primary AML cell in AML patients and prevent disease progression by activating caspase-3/8 and promoting the release of inflammatory factors (Yang et al., 2020). These findings suggested that pyroptosis provided a tumor-suppressive microenvironment and played an immunomodulatory role in AML. However, the prognostic influence of pyroptosis on AML patients was still largely unknown. Moreover, due to technical limitations, most previous studies were limited to a small number of pyroptosis-related genes (PRGs), whereas the involved PRGs might be far more in numbers and their antitumor effects were very likely to interact with each other in a highly coordinated manner. Therefore, a comprehensive analysis regarding the features of TIME cell infiltration mediated by multiple PRGs may provide a relatively whole profile of the function of PRGs and also further insights into the underlying mechanisms of AML occurrence and progression, which, however, has not been investigated so far as we know.
Based on these backgrounds, we, herein, systematically analyzed the differential expression of PRGs and the prognostic value of these genes in clinical practice between AML and normal samples and established an independent prognostic PRG signature. Subsequently, we also explored the relationship between pyroptosis and TIME as well as evaluated the sensitivity of therapeutic drugs for AML patients according to PRG prognostic signature. This study identified reliable prognostic biomarkers for AML patients and provided a novel scientific basis for future immunotherapy in AML.
Materials and methods
Acquisition of data
The specific analysis process of the present study is illustrated in Supplementary Figure S1. The RNA sequencing (RNA-seq) data of 151 AML patients’ bone marrow (BM) samples and 755 normal peripheral whole blood samples were acquired from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) and the Genotype-Tissue Expression (GTEx) database (https://www.gtexportal.org/). Meanwhile, we obtained the corresponding clinical features of 151 AML patients from TCGA and eventually included 132 samples after excluding 19 samples without survival time. Furthermore, the gene expression profile and relevant clinical characteristic data of 91 AML patients were downloaded from the BeatAML database (Tyner et al., 2018) (http://www.vizome.org/aml/) as a validation cohort. All patients’ clinical features are detailed in Supplementary Table S1. The RNA-seq data of raw count normalized from the aforementioned databases were utilized for differential expression analysis. Then, the count values were converted to transcripts per kilobase million (TPM) values, which were further transformed to log2 (TPM + 1) for subsequent analysis. In addition, the somatic mutation and copy number variation (CNV) data were also retrieved from TCGA database.
Pyroptosis-related genes
A total of 44 PRGs were retrieved from the GeneCards (https://www.genecards.org/) and previously published literature works (Man and Kanneganti, 2015; Wang and Yin, 2017; Karki and Kanneganti, 2019; Xia et al., 2019), all of which were protein-coding genes. The full details of these genes are listed in Supplementary Table S2. However, as the expression profile of GSDME and PTVK could not be acquired from the GTEx database, we finally selected 42 PRGs for further analysis. The location of PRGs on the chromosome was plotted by the “RCircos” R package. The “limma” R package was utilized to identify the differentially expressed pyroptosis–related genes (DEPRGs) with a p value <0.05 and |log2FC| > 0 between tumor and normal tissues. The PRG somatic mutation landscape was presented via the “maftools” R packages. The frequencies of genetic amplification and deletion were also summarized. The univariate Cox regression was performed to identify PRGs significantly associated with overall survival (OS) in TCGA cohort. Simultaneously, the Spearman correlation test was applied to evaluate the associations across all PRGs and the comprehensive results were visualized using the “igraph” R package.
Consensus clustering
Consensus clustering was adopted to identify the distinct pyroptosis-related patterns pertaining to the expression of pyroptosis regulators using k-means algorithms (Hartigan and Wong, 1979). The numbers and stability of clusters were determined by the consensus clustering algorithms of the “ConsensuClusterPlus” R package (Wilkerson and Hayes, 2010). We conducted 1,000 times repetitions to guarantee the classification stability.
Construction and validation of a prognostic gene signature by prognostic DEPRGs
The significant prognostic-related DEPRGs were presented via the “VennDiagram” R package. Subsequently, the least absolute shrinkage and selection operator (LASSO) Cox regression analysis using the “glmnet” R package was applied to screen out the optimal candidate gene combination to construct the prognostic gene signature (Tibshirani, 1997). The optimal value of the penalty parameter λ was determined by 10-fold cross-validation based on the minimum criteria. According to the coefficient calculated by LASSO regression and the standardized and normalized TCGA AML expression level, the individual risk score of each AML patient could be calculated using the following formula:
Simultaneously, the TCGA AML patients were separated into high- and low-risk categories according to the median risk score. The principal component analysis (PCA) and t-distributed stochastic neighbor–embedding (t-SNE) analysis were conducted using the “Rtsne” and “ggplot2” R packages to investigate the distribution of various groups in terms of gene expression levels in the constructed model. Thereafter, we performed Kaplan–Meier analysis via the “survminer” R package to assess the survival difference between the two categories. The “survival” and “timeROC” R packages were utilized to perform the time-dependent receiver operating characteristic (ROC) curve analysis, which was applied to evaluate the prognostic gene signature’s predictive value. The prognostic significances of gene signature and other clinical characteristics were further investigated using univariate and multivariate Cox regression analyses. Moreover, the same formula and statistical methods were used to further validate the prognostic capacity of the gene signature in the BeatAML cohort.
The chi-squared test was adopted to explore the association of gene signature and clinicopathological characteristics, which was visualized with a heatmap using the “pheatmap” R package. The Wilcoxon signed-rank test and Kruskal–Wallis H-test were utilized to compare the risk score among various categories of these clinicopathological characteristics, and the visualization of results was presented via the boxplots.
Establishment of the predictive nomogram
The independent clinical features (cytogenetic risk, age, and risk score) validated using univariate and multivariate Cox regression analyses were enrolled to construct a predictive prognosis nomogram using the “rms” and “survival” R packages. Time-dependent ROC curves for 1, 3, and 5 years were used to assess the performance of the nomogram. The calibration curves for 1-, 3-, and 5-year prediction were utilized to depict the consistency between predicted and actual survival. Furthermore, an alluvial diagram was drawn to show the changes in pyroptosis-related clusters, risk score, age, and cytogenetic risk using the “ggalluvial” R package.
Functional enrichment analysis of DEGs based on the high- and low-risk groups
AML patients in TCGA cohort were batched into two groups according to the median risk score, respectively. The same approach was also performed in the BeatAML cohort. Afterward, the differentially expressed genes (DEGs) were extracted by utilizing the “limma” R package with the criteria of FDR <0.05 and |log2FC| > 1. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were performed using the “clusterProfiler” R package (Wu et al., 2021).
Assessment of the tumor immune environment
We applied the ESTIMATE algorithm to calculate the stromal, immune scores, and ESTIMATE scores of each patient in TCGA and BeatAML cohorts (Yoshihara et al., 2013). Furthermore, considering the significant roles of the immune cell infiltration in the TIME, the infiltrating scores of 22 kinds of human immune cells in each AML sample were computed by the CIBERSORT algorithm with 1,000 permutation. Furthermore, the single-sample gene set enrichment analysis (ssGSEA) algorithm using the “gsva” R package was utilized to calculate infiltration abundance of 29 immune signatures in the AML TIME. Subsequently, we explored the association between the enrichment scores of 22 kinds of immune cells and the risk score or expression levels of the aforementioned identified optimal candidate genes. In addition, the Wilcoxon signed-rank test was applied to estimate the differences in expression of human lymphocyte antigen (HLA) signature and immune checkpoint genes between high- and low-risk groups.
Analysis of CSC, TMB, and drug susceptibility between high- and low-risk groups
The gene expression–based stemness index was acquired from the previous study (Malta et al., 2018). We explored the correlation between risk score and cancer stemness cell (CSC). Simultaneously, we also computed the tumor mutation burden (TMB) score for each AML sample based on the “maftools” R package and analyzed the differences in the TMB score between high- and low-risk groups, as well as the associations between the TMB score and risk score. To investigate the differences in efficacy of therapeutic drugs in patients between high- and low-risk categories, we calculated the semi-inhibitory concentration (IC50) values of drugs commonly used for the treatment of AML via the “pRRophetic” R package.
Statistical analysis
All statistical analyses were performed by R software (version 4.0.4). Association coefficients were calculated by the Spearman correlation test. Log-rank tests were used for identifying the significance of differences in Kaplan–Meier analysis curves. p values of less than 0.05 were considered statistically significant (*p < 0.05) in all analyses.
Results
Landscape of expression and genetic alterations of PRGs in AML
This study first summarized the incidence of CNVs and somatic mutations of 42 PRGs in AML. The exploration of CNVs demonstrated prevalent CNV alterations in all 42 PRGs, among which the CNV of TIRAP was significantly increased while the CNVs of CEBPB, PLCG1, and VDR were significantly decreased (Figure 1A). Furthermore, the location of CNV variation in the PRGs on their respective chromosomes is displayed in Figure 1B. In the following assessment of the genetic mutation of PRGs in depth, only CASP3 (1%), NLRC4 (1%), NLRP1 (1%), NLRP2 (1%), NLRP3 (1%), TFAM (1%), and TXNIP (1%) showed the genetic mutation in AML patients (Figure 1C). To determine whether these gene variants influence the expression of PRGs in patients with AML, we further calculated the mRNA expression levels of 42 PRG between normal and tumor specimens and then identified the statistically significant DEGs (p < 0.05) that correspond to them. As a result, 19 downregulated and 16 upregulated genes were observed in the tumor group (Figure 1D). Subsequently, we discovered that the CNV alteration may be a major factor resulting in perturbations on PRG expression levels. In comparison to normal samples, the expressions of genes with deleted CNVs were found to be significantly decreased in AML samples (such as CEBPB, PLCG1, and VDR), and vice versa (such as TIRAP). Based on the aforementioned results, the composite landscape of interaction and correlation of the 42 PRGs and their prognostic value for AML patients was further comprehensively displayed in a pyroptosis network, which showed that 12 pyroptosis genes exhibited a significant prognostic value (Figure 1E; Supplementary Tables S3, S4).
FIGURE 1. Genetic and transcriptional alteration of PRGs in AML. (A) Frequencies of CNV gain, loss, and non-CNV for 42 PRGs in TCGA cohort. (B) Locations of CNV alterations in PRGs on 23 chromosomes. (C) Mutation frequency of PRGs in TCGA cohort. (D) Heatmap of expression patterns of PRGs in normal and tumor samples. (E) Network of correlations including PRGs in TCGA cohort. PRGs, pyroptosis-related genes; AML, acute myeloid leukemia; CNV, copy number variant; TCGA, The Cancer Genome Atlas.
Tumor classification based on the DEPRGs
To further investigate the expression features of the PRGs in AML, we applied the consensus clustering algorithm to classify the patients with AML according to the expression spectrums of 35 DEPRGs (Supplementary Figures S2A–H). We found that k = 3 seemed to be the optimum alternative for categorizing the entire cohort into three subtypes (Figure 2A and Supplementary Figures S2I–J), based on which satisfied separation across the three clusters was achieved according to the PCA and t-SNE plots (Figures 2B,C). The OS time was also compared between the three groups, but no significant difference was observed (p = 0.11, Figure 2D). Furthermore, the heatmap displayed the PRG expression profile and clinicopathological characteristic such as race, gender, age, FAB classification, cytogenetic risk, FLT3-ITD mutation, and NPM1 mutation, with NPM1 mutation (chi-squared test: p < 0.01) and FAB classification (p < 0.001) demonstrating significant differences among the three groups (Figure 2E).
FIGURE 2. Tumor classification based on the DEPRGs. (A) 132 AML patients were classified into three clusters according to the consensus clustering matrix (k = 3). (B) PCA and (C) t-SNE analysis exhibiting a remarkable difference in transcriptomes among three clusters. (D) Kaplan–Meier OS curves for three clusters. (E) Heatmap and clinicopathological features of three clusters. DEPRGs, differentially expressed pyroptosis-related genes; AML, acute myeloid leukemia; PCA, principal component analysis; t-SNE, t-distributed stochastic neighbor embedding; OS, overall survival.
Gene signature construction and validation from prognostic DEPRGs
To further explore the prognostic value of the DEPRGs in AML patients, the gene signature was constructed. As shown in Figure 3A, 12 of 42 PRGs showed a significant prognostic value according to univariate Cox regression analysis, among which 10 genes were overlapped between DEPRGs and prognostic genes and were identified as the prognostic DEPRGs. The specific prognostic value of 10 genes is further depicted in Figure 3B and Supplementary Table S5 with five risk and five favorable factors, respectively. Subsequently, based on 1,000 times 10-fold cross-validation in LASSO Cox regression analysis, the minimum of the λ value was selected as the optimum λ value (0.0278). Then, six of 10 genes with not-zero coefficients were screened to construct the gene signature according to the optimum λ (Figures 3C,D). The risk score of each patient is computed as follows: risk score = (−0.0993153386641694 × expression of CASP3) + (−0.0669072839409446 × expression of ELANE) + (1.29694083393916 × expression of GSDMA) + (−0.164405100033389 × expression of NOD1) + (0.440274519545682 × expression of PYCARD) + (0.00707428724075096 × expression of VDR). Then, the patients were separated into the high-risk group (n = 66) and low-risk (n = 66) group according to the median risk score (Figure 4A). The distribution of the risk score indicated that the OS status in the high-risk category was significantly worse than that in the low-risk category (Figure 4B). Further PCA and t-SNE analyses showed identifiable dimensions between the high-risk and low-risk categories (Figures 4C,D). Meanwhile, the Kaplan–Meier survival curves demonstrated significantly superior OS in patients with low scores than those with high scores (p < 0.0001; Figure 4E). Furthermore, time-dependent ROC analysis also revealed that this gene signature exhibited a favorable prognostic performance with AUCs of 0.736, 0.769, and 0.815 at 1-, 3-, and 5-year, respectively (Figure 4F).
FIGURE 3. Construction of a PRG signature for AML in TCGA cohort. (A) Ten prognostic DEPRGs were identified via the Venn diagram. (B) Forest plots show the results of univariate Cox analysis of OS for 10 prognostic DEPRGs. (C) LASSO regression of the 10 prognostic DEPRGs. (D) Cross-validation for turning the parameter selection in the LASSO regression. PRG, pyroptosis-related gene; AML, acute myeloid leukemia; TCGA, The Cancer Genome Atlas; DEPRGs, differentially expressed pyroptosis-related genes; OS, overall survival; LASSO, the least absolute shrinkage and selection operator.
FIGURE 4. Prognostic study of TCGA cohort via the PRG signature model. (A) Risk scores and (B) survival status of AML patients. (C) PCA and (D) t-SNE analysis showing the different gene expression of samples. (E) Kaplan–Meier OS curves for high-risk and low-risk groups. (F) The 1-, 3-, and 5-year ROC curve to predict the survival status. TCGA, The Cancer Genome Atlas; PRG, pyroptosis-related gene; PCA, principal component analysis; t-SNE, t-distributed stochastic neighbor embedding; OS, overall survival; ROC, receiver operating characteristic.
To verify the capability of the 6-gene signature, the BeatAML dataset was downloaded as an external validation cohort. We calculated the risk scores of patients with the same formula and stratified patients into high-risk (n = 45) and low-risk (n = 46) groups (Supplementary Figure S3A). Similarly, patients in the two groups of the BeatAML cohort were distributed in various directions based on the PCA and t-SNE analyses (Supplementary Figures S2C–D). Similar to the results of TCGA cohort, patients with low scores were found to have significantly longer OS and a favorable prognostic value (p = 0.0034) in the validation cohort (Supplementary Figures S3B,E). In addition, the results of 1-, 3-, and 5-year ROC curves also possessed relatively higher AUC values (0.654, 0.800, and 0.659), indicating that the signature had excellent predictive capability for survival of AML patients (Supplementary Figure S3F).
Independent prognostic value of the 6-gene signature and evaluation of clinical characteristics
We performed univariate and multivariate Cox analyses to assess the possibility of the risk score functioning as an independent prognostic factor. The results of univariate Cox regression analysis demonstrated that the risk score was a significant independent predictor of poor survival in both TCGA (p < 0.001, HR = 3.176, 95% CI = 2.122–4.754; Figure 5A) and BeatAML cohorts (p = 0.005, HR = 3.108, 95% CI = 1.407–6.867; Figure 5C). Multivariate analysis also revealed that the risk score was a critical prognostic factor for AML patients in both cohorts after accounting for other confounders (p < 0.001, HR = 3.356, 95% CI = 1.958–5.753 for TCGA cohort and p = 0.031, HR = 2.602, 95% CI = 1.092–6.201 for the BeatAML cohort; Figures 5B,D).
FIGURE 5. Univariate and multivariate Cox regression analyses for the risk score. (A) Univariate and (B) multivariate Cox regression of the risk score and other clinical characteristics associated with OS in the TCGA cohort. (C) Univariate and (D) multivariate Cox regression of the risk score and other clinical characteristics associated with OS in the BeatAML cohort. (E) Heatmap for association between clinicopathologic characteristics and risk groups. OS, overall survival.
The chi-squared test of correlations between the risk subgroups and clinicopathological features indicated that NPM1 mutation (p < 0.05), cytogenetic risk (p < 0.001), FAB classification (p < 0.001), and PRG clusters (p < 0.001) were significantly associated with risk subgroups (Figure 5E). Subsequently, the Wilcoxon signed-rank test was applied to compare the differences in risk scores across respective ;groups of the aforementioned clinicopathological characteristics, also demonstrating a remarkable relationshipof these clinicopathological characteristics with risk scores (Supplementary Figures S4A–D).
In addition, we further assessed whether the risk score could still maintain a good prediction of survival across different subgroups. The results showed that low-risk patients also showed a better prognosis than high-risk patients in all subgroups, with statistically significant results in the following subgroups: race (white: p < 0.0001), gender (male: p = 0.018; female: p = 0.0001), age (≤55 years: p = 0.024; >55 years: p = 0.00045), FAB classification (M0: p = 0.017; M2: p = 0.011), cytogenetic risk (intermediate: p = 0.036), NPM1 mutation (positive: p = 0.0051; negative: p = 0.0029), and FLT3-ITD mutation (positive: p = 0.034; negative: p = 0.0019) (Supplementary Figures S5A–R).
Development of a prognostic nomogram for AML
In consideration of the fact that the risk score alone was not sufficient to predict OS in AML patients, a nomogram incorporating the risk score and clinicopathological features was constructed to forecast 1-, 3-, and 5-year OS of AML patients according to the significant results of multivariate Cox regression analysis (Figure 6A). Furthermore, ROC curves indicated that this nomogram exhibited a good prognostic performance with AUCs of 0.738, 0.768, and 0.815 at 1-, 3-, and 5-year, respectively (Figure 6B). The subsequent calibration plot showed the proposed nomogram operated in a manner consistent with an ideal model (Figure 6C). Moreover, the alluvial diagram was applied to visualize variations in the aforementioned characteristics of AML patients (Figure 6D).
FIGURE 6. Construction and evaluation of a PRG signature-based nomogram. (A) Nomogram incorporating cytogenetic risk, age, and risk score was constructed to predict 1-, 3-, and 5-year survival probabilities. (B) ROC curves for 1-, 3-, and 5-year were used to assess the performance of nomogram. (C) Calibration curves for 1-, 3-, and 5-year prediction were utilized to depict the consistency between predicted and actual survival. (D) Alluvial diagram of subgroup distributions with different risk scores and survival outcomes. PRG, pyroptosis-related gene; ROC, receiver operating characteristic.
Functional annotation of the 6-gene signature
To further investigate the potential biological functions and pathways of the 6-gene signature, the DEGs across the high-risk and low-risk categories were applied to perform GO and KEGG analyses. The specific details of results were presented in Supplementary Tables S6, S7. Furthermore, these DEGs of TCGA cohort were significantly enriched in biological processes and molecular functions with immunity, such as neutrophil activation involved in the immune response, leukocyte cell–cell adhesion, regulation of leukocyte proliferation, immune receptor activity, immunoglobulin binding, and IgG binding (Figure 7A). Several of these biological processes and molecular functions were validated in the BeatAML cohort, including neutrophil activation involved in the immune response, leukocyte cell–cell adhesion, and immune receptor activity (Figure 7C). The results of KEGG analysis also showed enrichment of immune-related pathways, which included cytokine–cytokine receptor interaction, phagosome, viral protein interaction with cytokine and cytokine receptor, intestinal immune network for IgA production, and leukocyte transendothelial migration in both TCGA and BeatAML cohorts (Figures 7B,D). Moreover, pathways regarding B-cell receptor signaling and Th1 and Th2 cell differentiation were also found in TCGA (Figure 7B) and BeatAML cohorts (Figure 7D), respectively. In addition, several cancer-related pathways were simultaneously identified in both cohorts, such as transcriptional misregulation in cancer, hematopoietic cell lineage, ECM-receptor interaction, and proteoglycans in cancer. These results revealed that the pyroptosis-related 6-gene signature was significantly associated with cancer progression and particularly possessed an important influence on the immunoregulation of TME.
FIGURE 7. GO and KEGG functional enrichment between high- and low-risk groups. Top 10 results of (A) GO and (B) KEGG pathway enrichment of DEGs among high- and low-risk groups in TCGA cohort. Top 10 results of (C) GO and (D) KEGG pathway enrichment of DEGs among high- and low-risk groups in the BeatAML cohort. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes.
Tumor immune microenvironment analysis of 6-gene signature
Next, we further explored the association between the pyroptosis-related gene signature and tumor immune microenvironment and observed that the low-risk patients were significantly correlated with inferior immune and ESTIMATE scores in the TCGA cohort (Figure 8A). Subsequently, based on the CIBERSORT algorithm, we compared the distribution of 22 kinds of immune cells in diverse risk subgroups. A significant difference in the distribution of immune cells was observed in high-risk patients with superior infiltration of monocytes and M2 macrophages, but inferior infiltration of plasma cells, resting memory CD4+ T cells, follicular helper T cells, activated mast cells, and resting mast cells (Figure 8C and Supplementary Table S8). Concurrently, we analyzed the connection between the risk score and the infiltration score of immune cells, which further suggested that the risk score was significantly associated with the six kinds of immune cells (Supplementary Figure S6). We also evaluated the relationship between the six pyroptosis-related genes in the proposed signature and abundance of immune cells. The results showed that partial immune cells were significantly associated with the six genes (Figure 8B). Furthermore, comparisons of 29 immune signatures provided by the ssGSEA algorithm revealed that high-risk patients exhibited higher infiltration scores of APC co-inhibition, APC co-stimulation, B cells, CCR, checkpointa, DCs, HLA, iDCs, increased inflammation , neutrophils, parainflammation, pDCs, T-cell co-inhibition, T helper cells, Tfh, TIL, and type I IFN response, whereas lower infiltration scores of mast cells (Figure 8D and Supplementary Table S9) indicated higher immune infiltration among high-risk AML patients.
FIGURE 8. Immune characteristics analysis in TCGA cohort. (A) Associations between risk score and both immune and stromal scores. (B) Correlations between the abundance of immune cells and six genes in the proposed signature. Comparison of the enrichment scores of (C) 22 kinds of immune cells and (D) 29 types of immune signatures between high- and low-risk groups. (E) Expression of HLA in the high- and low-risk groups. (F) Expression of immune checkpoints in the high- and low-risk groups. p values are shown as: ns, not significant. *p = < 0.05; **p = < 0.01; ***p = < 0.001. HLA, human lymphocyte antigen.
Considering that HLA-related genes play a critical role in regulating the immune response, we then compared the expression of HLA-related genes between different subgroups and observed that most of the HLA-related genes were upregulated in the high-risk group (Figure 8E). Furthermore, we investigated the correlation between 33 immune checkpoints and the 6-PRG signature. Figure 8F demonstrated that high-risk patients exhibited significantly higher expression of PDCD1, CD200R1, CAG3, and LGALS9 as well as lower expression of CD160, NRP1, and TMIGD2 compared with low-risk patients. In addition, similar results were also observed in the BeatAML cohort (Supplementary Figures S7, S8; Supplementary Tables S10, S11).
Analysis of CSC index, TMB, and drug susceptibility
Considering that the CSC index and TMB play a critical role in the pathogenesis and immunotherapy of AML (Eppert et al., 2011; Snyder et al., 2014), we further explored the potential correlation between them and the PRG signature. As shown in Figure 9A, the risk score was negatively associated with the CSC index (R = −0.29, p = 0.0011), suggesting that AML cells with a lower risk score exhibited more significant stem cell characteristics and a lower degree of cell differentiation. Meanwhile, we also found that TMB of the low-risk group was significantly higher than that of the high-risk group (p = 0.039; Figure 9B) and was negatively associated with the PRG risk score (R = −0.25, p = 0.019; Figure 9C). Furthermore, the Kaplan–Meier survival curve for combining the PRG risk score and TMB revealed significant differences in survival outcomes and patients with high TMB and low PRG risk scores exhibited a more pronounced survival advantage (Figure 9D).
FIGURE 9. Comprehensive analysis of the risk score in AML. (A) Associations between the risk score and CSC index. (B) Comparison of TMB between high- and low-risk groups. (C) Spearman correlation analysis of risk score and TMB. (D) Kaplan–Meier OS curves for different TMB and risk score subgroups. (E–I) Relationships between risk score and therapeutic drugs of AML sensitivity. AML, acute myeloid leukemia; CSC, cancer stemness cell; TMB, tumor mutation burden; OS, overall survival.
Moreover, given the impact of drug susceptibility on patients with AML, we further selected the drugs currently used in the treatment of AML to assess the sensitivity of these drugs to patients in both high-risk and low-risk subgroups. Interestingly, we observed that patients in the low-risk group presented significantly lower IC50 values for ATRA, cytarabine, midostaurin, doxorubicin, and etoposide than those in the high-risk group (Figures 9E–I). Therefore, patients with a low PRG risk score might exhibit better treatment benefits when administrating these drugs. Nevertheless, the effect of these drugs in the treatment of AML patients remains to be further proven in future clinical studies.
Discussion
AML is a rapidly progressive hematologic tumor with poor prognosis, and its development can be influenced by multiple factors (Coombs et al., 2016), among which genes associated with programmed cell death (PCD), such as autophagy-related genes and ferroptosis-related genes, were demonstrated to serve as reliable prognostic biomarkers for AML (Fu et al., 2021; Shao et al., 2021). However, studies regarding the role of PRGs in AML are restricted to individual PRG and have never been investigated systematically and comprehensively. In this study, we detected the global alterations in PRGs at transcriptional and genetic levels in AML and identified a total of 42 currently available PRGs, the majority of which were differentially expressed between AML and normal samples and were associated with prognostic of this disease. Subsequently, we first constructed a reliable and valid PRG signature for AML based on six PRGs. In both training and validation cohorts, the PRG signature exhibited robust capabilities in predicting survival outcomes in AML patients. Furthermore, patients with low and high PRG scores showed significantly different clinicopathological features. In addition, univariate and multivariate Cox regression analyses also indicated that the PRG signature was an independent prognostic factor. These findings confirm that the present PRG signature can be used as a potentially reliable prognostic biomarker in AML patients with various clinicopathological features.
Pyroptosis, a novel form of PCD, is activated via the classical caspase-1 inflammasome or non-classical caspase-4, caspase-5, and caspase-11–mediated pathways (Tang et al., 2020). Recently, the dual roles of pyroptosis in tumor progression gained substantial attention. On one hand, it was reported that inflammatory molecules released from cancer cells underwent pyroptosis and gradually converted the surrounding normal cells into cancer cells by altering the microenvironment, which in turn promoted tumor development. On the other hand, induction of tumor cell pyroptosis was also found to inhibit tumor progression and was demonstrated to be a potential therapeutic target for drug development (Tang et al., 2020). This situation raised great uncertainties regarding the exact functions of pyroptosis in different kinds of tumors. Furthermore, as a matter of fact, previous studies paid much more attention to the roles of PRGs in solid carcinomas [e.g., hepatocellular carcinoma (Deng et al., 2022), breast cancer (Xu D. et al., 2021), lung adenocarcinoma (Lin W. et al., 2021), and colon adenocarcinoma (Luo et al., 2021)] with few studies, however, focusing on non-solid tumor such as AML. Therefore, the present study provided a novel signature featuring six PRGs (CASP3, ELANE, GSDMA, PYCARD, VDR, and NOD1) that exactly and independently predict OS in AML patients, which will further contribute to driving the progress of individualized prevention and treatment of AML.
Among these PRGs, caspase-3 (CASP3) is an important member of the caspase family, whose activation degrades structural and functional proteins within cells, thereby inducing cell death (Yuan et al., 1993; Jiang et al., 2020). In the present study, we further revealed that this gene was a favorable predictor of survival outcome and was associated with increased sensitivity to chemotherapeutic drug–induced pyroptosis in AML, which was in accordance with the previous study indicating that CASP3 activated by chemotherapeutic drugs initiated pyroptosis (Wang et al., 2017). Moreover, as one of the primary serine proteases secreted by neutrophils, ELANE is another known promoter of pyroptosis, which activates inflammatory factors such as TNF-α, IL-1β, and IL-18, and induces neutrophils to develop pyroptosis (Kambara et al., 2018; Mirea et al., 2020). Consistently, our study demonstrated that the expression of ELANE was significantly higher while the neutrophil infiltration score was remarkably lower in the low-risk group than that in the high-risk group, which was very likely to be ascribed to its pyroptosis activating effect in neutrophils. In addition, the following three genes (GSMDA, PYCARD, and VDR) that were previously identified as possible executors of pyroptosis and usually exhibited tumor-suppressive effects (Ding et al., 2016; Šutić et al., 2019; Ling et al., 2022) were also identified and included in the present PRG signature, further confirming its reliability.
Interestingly, as a cytoplasmic pattern recognition receptor, NOD1 was initially identified as a cancer-promoting factor and might cause tumor recurrence and metastasis, resulting in a poorer prognosis through pyroptosis (Fernández-García et al., 2022; Nomoto et al., 2022). However, in our research, NOD1 showed a significant cancer suppressive effect and acted as a protective factor against AML. This discrepancy may be attributed to the specialized tumor microenvironment in non-solid tumors and also the antitumor immune activity generated by the combined action of multiple PRGs although the exact mechanism still needs to be further explored.
Previous studies showed that the pro-inflammatory effects of pyroptosis are strongly associated with the regulation of the TIME (Tsuchiya, 2021). This study hence further evaluated the association between the risk score of PRG signature and TIME and found that patients in the high-risk group showed significantly higher immune scores than those in the low-risk group. The abundance of infiltration of tumor immune cells also differed between the high- and low-risk groups. Compared to the high-risk group, resting CD4+ memory T cells and follicular helper T cells, both of which were well acknowledged to exert an important antitumor immune response (Watanabe, 2021), infiltrated at higher levels in the low-risk group; whereas tumorigenesis-, angiogenesis- and immune suppressing-related cell types, such as M2 macrophages, monocytes, antigen-presenting cells, and dendritic cells (Nahas et al., 2019; Fu and Song, 2021), showed higher infiltration levels in the high-risk group. Furthermore, this study also found that patients with high-risk scores showed a worse prognosis than those with low-risk scores. These findings indicated that the immunosuppressive microenvironment played important roles during the genesis and development of AML. In fact, it was reported that the formation of an immunosuppressive microenvironment usually prevented the clearance of tumor cells by tumor killer cells, resulting in an increased risk of malignant progression and death (Fridman et al., 2022). Therefore, treatment targeting the immunosuppressive microenvironment may be a more effective and feasible strategy for patients with a poor prognosis of AML.
In addition, we found that most HLA-related genes and PDCD1 expressed at higher levels in the high-risk group, which was in accordance with the current increasing evidence regarding solid tumors suggesting that more HLA presentation increased the recognition of tumor-associated antigens in HLA and in turn increased the success of immune checkpoint inhibitor therapy (Rizvi et al., 2015; Lin W. Y. et al., 2021). Therefore, patients with high-risk scores for AML might benefit more from immunotherapy, especially with immune checkpoint inhibitors PD-1.
Although we analyzed the effects of pyroptosis on the prognosis of AML and the immune microenvironment as comprehensively as possible, the following points remained inadequate. First, all analyses in this study were performed based on retrospective data from public databases, and large prospective studies and additional in vivo and in vitro experimental studies are still needed to confirm our findings. Second, the proposed gene signature model was validated only through public databases. Therefore, further clinical trials are still necessary to confirm its clinical utility. In addition, the potential mechanisms of the present six key genes used for model construction need to be further explored in order to better understand the tumorigenesis and development of AML. Finally, subgroup analyses for ethnicity were not performed in our study mainly due to the limited availability of current original data, which should also be validated in future studies based on additional risk models.
In summary, we comprehensively analyzed the expression and genetic changes of PGRs in AML, their prognostic value in the clinic, their important role in TIME, and constructed a signature consisting of six PRGs, which was confirmed to be an independent predictor for OS in AML patients. The results of this study will contribute to further understanding of the important role of pyroptosis in the prognosis and development of AML and provide novel and reliable biomarkers for its precise prevention and treatment.
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.
Author contributions
G-HH conceived and designed the study, obtained the funding. TZ performed the majority of data analysis and wrote the manuscript. W-KC, S-JY, and Y-YL participated in checking the statistic method and helped to draft the manuscript. KQ and PW prepared the figures. All authors revised and approved the final manuscript.
Funding
This work was supported by the Grants from the National Science Foundation of China (No. 81960664), the Applied Basic Research Program Yunnan Province of China (Joint Special Project of Kunming Medical University) (No. 202101AY070001-300), and the Science and Technology Planning Project of 920th Hospital of Joint Logistics Support Force (No. 2020YGB11).
Acknowledgments
All authors sincerely acknowledge the contributions from the TCGA, GTEx, and BeatAML projects.
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/fphar.2022.951480/full#supplementary-material
References
Arneth, B. (2019). Tumor microenvironment. Med. Kaunas. Lith. 56 (1), 15. doi:10.3390/medicina56010015
Balahura, L. R., Selaru, A., Dinescu, S., and Costache, M. (2020). Inflammation and inflammasomes: pros and cons in tumorigenesis. J. Immunol. Res. 2020, 2549763. doi:10.1155/2020/2549763
Ball, B., and Stein, E. M. (2019). Which are the most promising targets for minimal residual disease-directed therapy in acute myeloid leukemia prior to allogeneic stem cell transplant? Haematologica 104 (8), 1521–1531. doi:10.3324/haematol.2018.208587
Broz, P., von Moltke, J., Jones, J. W., Vance, R. E., and Monack, D. M. (2010). Differential requirement for Caspase-1 autoproteolysis in pathogen-induced cell death and cytokine processing. Cell Host Microbe 8 (6), 471–483. doi:10.1016/j.chom.2010.11.007
Coombs, C. C., Tallman, M. S., and Levine, R. L. (2016). Molecular therapy for acute myeloid leukaemia. Nat. Rev. Clin. Oncol. 13 (5), 305–318. doi:10.1038/nrclinonc.2015.210
Deng, M., Sun, S., Zhao, R., Guan, R., Zhang, Z., Li, S., et al. (2022). The pyroptosis-related gene signature predicts prognosis and indicates immune activity in hepatocellular carcinoma. Mol. Med. 28 (1), 16. doi:10.1186/s10020-022-00445-0
Ding, J., Wang, K., Liu, W., She, Y., Sun, Q., Shi, J., et al. (2016). Pore-forming activity and structural autoinhibition of the gasdermin family. Nature 535 (7610), 111–116. doi:10.1038/nature18590
Eppert, K., Takenaka, K., Lechman, E. R., Waldron, L., Nilsson, B., van Galen, P., et al. (2011). Stem cell gene expression programs influence clinical outcome in human leukemia. Nat. Med. 17 (9), 1086–1093. doi:10.1038/nm.2415
Erkes, D. A., Cai, W., Sanchez, I. M., Purwin, T. J., Rogers, C., Field, C. O., et al. (2020). Mutant BRAF and MEK inhibitors regulate the tumor immune microenvironment via pyroptosis. Cancer Discov. 10 (2), 254–269. doi:10.1158/2159-8290.CD-19-0672
Fernández-García, V., González-Ramos, S., Avendaño-Ortiz, J., Martín-Sanz, P., Delgado, C., Castrillo, A., et al. (2022). NOD1 splenic activation confers ferroptosis protection and reduces macrophage recruitment under pro-atherogenic conditions. Biomed. Pharmacother. 148, 112769. doi:10.1016/j.biopha.2022.112769
Fridman, W. H., Meylan, M., Petitprez, F., Sun, C. M., Italiano, A., and Sautès-Fridman, C. (2022). B cells and tertiary lymphoid structures as determinants of tumour immune contexture and clinical outcome. Nat. Rev. Clin. Oncol. 15, 1–17. Advance online publication. doi:10.1038/s41571-022-00619-z
Fu, D., Zhang, B., Wu, S., Zhang, Y., Xie, J., Ning, W., et al. (2021). Prognosis and characterization of immune microenvironment in acute myeloid leukemia through identification of an autophagy-related signature. Front. Immunol. 12, 695865. doi:10.3389/fimmu.2021.695865
Fu, X. W., and Song, C. Q. (2021). Identification and validation of pyroptosis-related gene signature to predict prognosis and reveal immune infiltration in hepatocellular carcinoma. Front. Cell Dev. Biol. 9, 748039. doi:10.3389/fcell.2021.748039
Fu, Y., Liu, S., Zeng, S., and Shen, H. (2019). From bench to bed: the tumor immune microenvironment and current immunotherapeutic strategies for hepatocellular carcinoma. J. Exp. Clin. Cancer Res. 38 (1), 396. doi:10.1186/s13046-019-1396-4
Hartigan, J. A., and Wong, M. A. (1979). Algorithm as 136: a k-means clustering algorithm. Appl. Stat. 28 (1), 100. doi:10.2307/2346830
Jiang, M., Qi, L., Li, L., and Li, Y. (2020). The caspase-3/GSDME signal pathway as a switch between apoptosis and pyroptosis in cancer. Cell Death Discov. 6, 112. doi:10.1038/s41420-020-00349-0
Johnson, D. C., Taabazuing, C. Y., Okondo, M. C., Chui, A. J., Rao, S. D., Brown, F. C., et al. (2018). DPP8/DPP9 inhibitor-induced pyroptosis for treatment of acute myeloid leukemia. Nat. Med. 24 (8), 1151–1156. doi:10.1038/s41591-018-0082-y
Kambara, H., Liu, F., Zhang, X., Liu, P., Bajrami, B., Teng, Y., et al. (2018). Gasdermin D exerts anti-inflammatory effects by promoting neutrophil death. Cell Rep. 22 (11), 2924–2936. doi:10.1016/j.celrep.2018.02.067
Karki, R., and Kanneganti, T. D. (2019). Diverging inflammasome signals in tumorigenesis and potential targeting. Nat. Rev. Cancer 19 (4), 197–214. doi:10.1038/s41568-019-0123-y
Lin, W., Chen, Y., Wu, B., Chen, Y., and Li, Z. (2021). Identification of the pyroptosis-related prognostic gene signature and the associated regulation axis in lung adenocarcinoma. Cell Death Discov. 7 (1), 161. doi:10.1038/s41420-021-00557-2
Lin, W. Y., Fordham, S. E., Hungate, E., Sunter, N. J., Elstob, C., Xu, Y., et al. (2021). Genome-wide association study identifies susceptibility loci for acute myeloid leukemia. Nat. Commun. 12 (1), 6233. doi:10.1038/s41467-021-26551-x
Ling, Y., Xu, F., Xia, X., Dai, D., Sun, R., Xie, Z., et al. (2022). Vitamin D receptor regulates proliferation and differentiation of thyroid carcinoma via the E-cadherin-β-catenin complex. J. Mol. Endocrinol. 68 (3), 137–151. doi:10.1530/JME-21-0167
Luo, B., Lin, J., Cai, W., and Wang, M. (2021). Identification of the pyroptosis-related gene signature and risk score model for colon adenocarcinoma. Front. Genet. 12, 771847. doi:10.3389/fgene.2021.771847
Malta, T. M., Sokolov, A., Gentles, A. J., Burzykowski, T., Poisson, L., Weinstein, J. N., et al. (2018). Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell 173 (2), 338–354. e15. doi:10.1016/j.cell.2018.03.034
Man, S. M., and Kanneganti, T. D. (2015). Regulation of inflammasome activation. Immunol. Rev. 265 (1), 6–21. doi:10.1111/imr.12296
Mirea, A. M., Stienstra, R., Kanneganti, T. D., Tack, C. J., Chavakis, T., Toonen, E., et al. (2020). Mice deficient in the IL-1β activation genes Prtn3, elane, and Casp1 are protected against the development of obesity-induced NAFLD. Inflammation 43 (3), 1054–1064. doi:10.1007/s10753-020-01190-4
Moujalled, D., Strasser, A., and Liddell, J. R. (2021). Molecular mechanisms of cell death in neurological diseases. Cell Death Differ. 28 (7), 2029–2044. doi:10.1038/s41418-021-00814-y
Nahas, M. R., Stroopinsky, D., Rosenblatt, J., Cole, L., Pyzer, A. R., Anastasiadou, E., et al. (2019). Hypomethylating agent alters the immune microenvironment in acute myeloid leukaemia (AML) and enhances the immunogenicity of a dendritic cell/AML vaccine. Br. J. Haematol. 185 (4), 679–690. doi:10.1111/bjh.15818
Nomoto, D., Baba, Y., Liu, Y., Tsutsuki, H., Okadome, K., Harada, K., et al. (2022). Fusobacterium nucleatum promotes esophageal squamous cell carcinoma progression via the NOD1/RIPK2/NF-κB pathway. Cancer Lett. 530, 59–67. doi:10.1016/j.canlet.2022.01.014
Orning, P., Lien, E., and Fitzgerald, K. A. (2019). Gasdermins and their role in immunity and inflammation. J. Exp. Med. 216 (11), 2453–2465. doi:10.1084/jem.20190545
Perna, F., Berman, S. H., Soni, R. K., Mansilla-Soto, J., Eyquem, J., Hamieh, M., et al. (2017). Integrating proteomics and transcriptomics for systematic combinatorial chimeric antigen receptor therapy of AML. Cancer Cell 32 (4), 506–519. doi:10.1016/j.ccell.2017.09.004
Rizvi, N. A., Hellmann, M. D., Snyder, A., Kvistborg, P., Makarov, V., Havel, J. J., et al. (2015). Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science 348 (6230), 124–128. doi:10.1126/science.aaa1348
Schlenk, R. F., Jaramillo, S., and Müller-Tidow, C. (2019). What's new in consolidation therapy in AML? Semin. Hematol. 56 (2), 96–101. doi:10.1053/j.seminhematol.2018.08.005
Shao, R., Wang, H., Liu, W., Wang, J., Lu, S., Tang, H., et al. (2021). Establishment of a prognostic ferroptosis-related gene profile in acute myeloid leukaemia. J. Cell. Mol. Med. 25 (23), 10950–10960. doi:10.1111/jcmm.17013
Snyder, A., Makarov, V., Merghoub, T., Yuan, J., Zaretsky, J. M., Desrichard, A., et al. (2014). Genetic basis for clinical response to CTLA-4 blockade in melanoma. N. Engl. J. Med. 371 (23), 2189–2199. doi:10.1056/NEJMoa1406498
Sperlazza, J., Rahmani, M., Beckta, J., Aust, M., Hawkins, E., Wang, S. Z., et al. (2015). Depletion of the chromatin remodeler CHD4 sensitizes AML blasts to genotoxic agents and reduces tumor formation. Blood 126 (12), 1462–1472. doi:10.1182/blood-2015-03-631606
Stölzel, F., Röllig, C., Radke, J., Mohr, B., Platzbecker, U., Bornhäuser, M., et al. (2011). 18F-FDG-PET/CT for detection of extramedullary acute myeloid leukemia. Haematologica 96 (10), 1552–1556. doi:10.3324/haematol.2011.045047
Šutić, M., Motzek, A., Bubanović, G., Linke, M., Sabol, I., Vugrek, O., et al. (2019). Promoter methylation status of ASC/TMS1/PYCARD is associated with decreased overall survival and TNM status in patients with early-stage non-small cell lung cancer (NSCLC). Transl. Lung Cancer Res. 8 (6), 1000–1015. doi:10.21037/tlcr.2019.12.08
Tang, R., Xu, J., Zhang, B., Liu, J., Liang, C., Hua, J., et al. (2020). Ferroptosis, necroptosis, and pyroptosis in anticancer immunity. J. Hematol. Oncol. 13 (1), 110. doi:10.1186/s13045-020-00946-7
Tibshirani, R. (1997). The lasso method for variable selection in the Cox model. Stat. Med. 16 (4), 385–395. doi:10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3
Tsuchiya, K. (2021). Switching from apoptosis to pyroptosis: gasdermin-elicited inflammation and antitumor immunity. Int. J. Mol. Sci. 22 (1), 426. doi:10.3390/ijms22010426
Tyner, J. W., Tognon, C. E., Bottomly, D., Wilmot, B., Kurtz, S. E., Savage, S. L., et al. (2018). Functional genomic landscape of acute myeloid leukaemia. Nature 562 (7728), 526–531. doi:10.1038/s41586-018-0623-z
Wang, B., and Yin, Q. (2017). AIM2 inflammasome activation and regulation: a structural perspective. J. Struct. Biol. 200 (3), 279–282. doi:10.1016/j.jsb.2017.08.001
Wang, S., Yuan, P., Mao, B., Li, N., Ying, J., Tao, X., et al. (2022). Genomic features and tumor immune microenvironment alteration in NSCLC treated with neoadjuvant PD-1 blockade. NPJ Precis. Oncol. 6 (1), 2. doi:10.1038/s41698-021-00244-6
Wang, Y., Gao, W., Shi, X., Ding, J., Liu, W., He, H., et al. (2017). Chemotherapy drugs induce pyroptosis through caspase-3 cleavage of a gasdermin. Nature 547 (7661), 99–103. doi:10.1038/nature22393
Watanabe, T. (2021). The tumor microenvironment in follicular lymphoma: its pro-malignancy role with therapeutic potential. Int. J. Mol. Sci. 22 (10), 5352. doi:10.3390/ijms22105352
Wilkerson, M. D., and Hayes, D. N. (2010). ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics 26 (12), 1572–1573. doi:10.1093/bioinformatics/btq170
Wu, T., and Dai, Y. (2017). Tumor microenvironment and therapeutic response. Cancer Lett. 387, 61–68. doi:10.1016/j.canlet.2016.01.043
Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation 2 (3), 100141. doi:10.1016/j.xinn.2021.100141
Xia, X., Wang, X., Cheng, Z., Qin, W., Lei, L., Jiang, J., et al. (2019). The role of pyroptosis in cancer: pro-cancer or pro-"host" Cell Death Dis. 10 (9), 650. doi:10.1038/s41419-019-1883-8
Xu, B., Lu, M., Yan, L., Ge, M., Ren, Y., Wang, R., et al. (2021). A pan-cancer analysis of predictive methylation signatures of response to cancer immunotherapy. Front. Immunol. 12, 796647. doi:10.3389/fimmu.2021.796647
Xu, D., Ji, Z., and Qiang, L. (2021). Molecular characteristics, clinical implication, and cancer immunity interactions of pyroptosis-related genes in breast cancer. Front. Med. 8, 702638. doi:10.3389/fmed.2021.702638
Yang, J., Hong, S., Zhang, X., Liu, J., Wang, Y., Wang, Z., et al. (2021). Tumor immune microenvironment related gene-based model to predict prognosis and response to compounds in ovarian cancer. Front. Oncol. 11, 807410. doi:10.3389/fonc.2021.807410
Yang, W., Liu, S., Li, Y., Wang, Y., Deng, Y., Sun, W., et al. (2020). Pyridoxine induces monocyte-macrophages death as specific treatment of acute myeloid leukemia. Cancer Lett. 492, 96–105. doi:10.1016/j.canlet.2020.08.018
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612
Yuan, J., Shaham, S., Ledoux, S., Ellis, H. M., and Horvitz, H. R. (1993). The C. elegans cell death gene ced-3 encodes a protein similar to mammalian interleukin-1 beta-converting enzyme. Cell 75 (4), 641–652. doi:10.1016/0092-8674(93)90485-9
Zeng, T., Cui, L., Huang, W., Liu, Y., Si, C., Qian, T., et al. (2021). The establishment of a prognostic scoring model based on the new tumor immune microenvironment classification in acute myeloid leukemia. BMC Med. 19 (1), 176. doi:10.1186/s12916-021-02047-9
Keywords: pyroptosis, acute myeloid leukemia, prognosis, immune microenvironment, TCGA
Citation: Zhou T, Qian K, Li Y-Y, Cai W-K, Yin S-J, Wang P and He G-H (2022) The pyroptosis-related gene signature predicts prognosis and reveals characterization of the tumor immune microenvironment in acute myeloid leukemia. Front. Pharmacol. 13:951480. doi: 10.3389/fphar.2022.951480
Received: 24 May 2022; Accepted: 13 July 2022;
Published: 10 August 2022.
Edited by:
Victor Mangas Sanjuan, University of Valencia, SpainReviewed by:
Jorge Duconge, University of Puerto Rico, Puerto RicoBo Zheng, Naval Medical Center, China
Copyright © 2022 Zhou, Qian, Li, Cai, Yin, Wang and He. 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: Gong-Hao He, Z29uZ2hvd0Bob3RtYWlsLmNvbQ==
†These authors have contributed equally to this work and share first authorship.