Identification of a Pyroptosis-Related Gene Signature for Prediction of Overall Survival in Lung Adenocarcinoma

Dong, Zheng; Bian, Lv; Wang, Minglang; Wang, Luoqing; Wang, Yilian

doi:https://doi.org/10.1155/2021/6365459

Journal of Oncology

On this page

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

Special Issue

Bioinformatics in Cancer and Immune Microenvironment

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 6365459 | https://doi.org/10.1155/2021/6365459

Identification of a Pyroptosis-Related Gene Signature for Prediction of Overall Survival in Lung Adenocarcinoma

Zheng Dong,¹Lv Bian,²Minglang Wang,¹Luoqing Wang,¹and Yilian Wang¹

Academic Editor: Jimei Wang

Received08 Aug 2021

Accepted13 Sept 2021

Published30 Sept 2021

Abstract

Pyroptosis is a kind of programmed cell death that is characterized by inflammation. However, the expression of pyroptosis-related genes and their connection with prognosis in lung adenocarcinoma (LUAD) remain unknown. The aim of this study is to create and validate a LUAD prediction signature based on genes associated with pyroptosis. The TCGA and GEO were used to collect gene sequencing data and clinical information for LUAD samples. To identify patients with LUAD from the TCGA cohort, consensus clustering by pyroptosis-related genes was employed. Our prognostic model was constructed using LASSO-Cox analysis after Cox regression using differentially expressed genes. To predict patient survival, we created a seven-mRNA signature. Additionally, reliability and validity were established in the GEO cohort. To assess its diagnostic and prognostic usefulness, an integrated bioinformatics method was used. Using a risk score with varying overall survival (OS) in two cohorts (all ), a seven-gene signature was developed to categorize patients into two risk categories. The signature was shown to be an independent predictor of LUAD using multivariate regression analysis. The signature was linked to a variety of immune cell subtypes according to a study of immune cell infiltration. We constructed a signature consisting of seven genes as a robust biomarker with potential for clinical use in risk stratification and OS prediction in LUAD patients, as well as a potential indicator of immunotherapy in LUAD.

1. Introduction

Worldwide, lung cancer causes a significant percentage of cancer mortality and has a bad prognosis [1, 2]. Lung cancer accounted for more than a quarter (27%) of all cancer fatalities in 2015 [3, 4]. Clinically, LUAD is most often diagnosed in patients with nonsmall cell lung cancer (NSCLC), while the most prevalent histologic form of NSCLC accounts for 40% of all lung cancer cases [5]. Surgical excision, chemotherapy, and radiation are the primary clinical therapies for LUAD at the moment. Despite advances in medical therapies in recent years, only 15% of LUAD patients survive for more than 5 years [6, 7]. Numerous studies have shown that LUAD is a very diverse illness with unique genetic and transcriptome features among individual individuals, and predicting LUAD prognosis remains difficult. As a result, it is critical to discover novel prognostic gene signatures that may be utilized to generate prognostic predictions and act as new therapeutic targets for LUAD patients.

Pyroptosis is the cleavage of gasdermins through conventional and nonclassical mechanisms, which may cause cells to expand indefinitely until the cell membrane ruptures, allowing the contents of the cell to escape, resulting in a severe inflammatory reaction [8, 9]. Pyroptosis is critical for combating infection and endogenous danger signals. Pyroptosis, by releasing inflammatory chemicals, provides a tumor-suppressive environment; yet, it may also impair the body’s immunological response to tumor cells and promote tumor development in certain malignancies [10, 11]. However, it is unknown that how pyroptosis affects the prognosis of LUAD.

The use of next-generation sequencing to classify LUAD patients is a new technique for rapidly identifying cancer features and determining the best therapy options. Although EGFR tyrosine kinase inhibitors (TKIs) and PD-L1/PD1 immune checkpoint targeting treatments have been shown effective, the underlying mechanisms that contribute to the lack of efficacy in patients with LUAD or even complete unresponsiveness to these therapies have still to be identified [12, 13]. According to the current understanding, molecular subtypes cannot be used as guides in clinical therapy due to the absence of subgroup classifications. As a result, the creation of a reliable gene signature to predict prognosis and guide therapeutic treatment, particularly in the areas of targeted therapy and immunotherapy, is critical.

The goal of this research was to develop a model for predicting prognosis and guiding therapeutic therapy by categorizing LUAD patients based on genes associated with pyroptosis. After grouping 414 people with LUAD based on genes associated with pyroptosis, we identified two subgroups that were connected to prognosis. The risk score may then be calculated by utilizing the LASSO-Cox technique to build a model associated with pyroptosis. This risk score may forecast prognosis and immunological infiltration. Our results suggest a possible link between pyroptosis, prognosis, and the immunological milieu of individuals with LUAD.

2. Materials and Methods

2.1. Data Collection

We conducted a comprehensive search of publicly accessible transcriptome cohorts for LUAD with matching clinical data. The TCGA database was used to obtain normalized gene expression data from RNA sequencing (FPKM), somatic mutation data, and clinical data from TCGA-LUAD cohort. Using the total number of mutations identified in the exome content, the tumor mutation burden (TMB) was then calculated per megabase for each sample. The GEO database was also used to gather clinical information and normalized gene expression data; the accession number is GSE13507. R was used to examine the data (version 4.0.2).

2.2. Identification of Pyroptosis-Related Genes with Variable Expression

From previous studies [14–17], we retrieved 52 pyroptosis-related genes, which are listed in Table S1. In order to discover the genes that were differentially expressed between normal and tumor tissues, we collected 19 normal lung tissues as well as 414 LUAD tissues. DEGs with a value of 0.05 were identified using the “limma” package. A heatmap was then utilized to compare the gene expression associated with pyroptosis in LUAD versus normal lung tissues.

2.3. Consensus Clustering

Consensus clustering using k-means algorithms was used to discover unique pyroptosis-related patterns pertaining to gene expression [18]. The number and stability of clusters were established using the consensus clustering method, which was implemented in the “ConsensuClusterPlus” package. We repeated our categorization 1,000 times to ensure its stability [19].

2.4. The Development and Validation of a Predictive Pyroptosis-Related Pattern Gene Signature

With an adjusted value of 0.05 and an absolute value of |log₂FC| ≥ 1, in TCGA cohort, DEGs between pyroptosis-related patterns were identified. Using the univariate Cox method, the DEGs with prognostic significance were identified, and the value was corrected using the Benjamini & Hochberg (BH) correction approach. A prognostic model was built using LASSO-Cox regression analysis to minimize the risk of overfitting [20, 21]. The LASSO method was used to select and minimize variables. A regression analysis was performed using the normalized expression matrix of potential prognostic DEGs as an independent variable, with the OS and status of patients in the TCGA cohort serving as the dependent variables. Patients’ risk scores were computed using a method that included finding a patient’s level of expression for each pattern gene associated with pyroptosis and the regression coefficient associated with that level. The formula was devised as follows:

Patients were classified into high-and low-risk categories based on their median risk score. In order to investigate the distribution of various groups in terms of gene expression levels in the constructed model, PCA and t-SNE analyses were carried out using the “Rtsne” and “ggplot2” R packages, respectively [22]. Using the “survminer” R package, a survival analysis was performed to compare the OS of high-and low-risk groups. Survival and TimeROC packages were used to perform time-dependent ROC curve analysis, which was used to assess the prognostic signature’s predictive value. The signature’s prognostic significance was further investigated using univariate and multivariate Cox analysis.

2.5. Analysis of Functional Enrichment of DEGs between Patterns

Consensus clustering classified LUAD patients in the TCGA cohort into two categories. The software “clusterProfiler” was used to perform GO and KEGG analysis on the DEGs. The ssGSEA was performed using the “gsva” package to compute infiltrating immune cell scores and assess immune-related pathway activity [23, 24].

2.6. Predictive Nomogram Establishment

A nomogram was developed by taking into consideration the patient’s risk score and other clinicopathological features and then using that information to provide a reliable clinical prediction tool for LUAD patients, especially with regard to their 1-, 3-, and 5-year survival. Following that we conducted calibration curve analyses to determine the suitability of our nomogram for clinical use.

2.7. Statistical Analysis

For statistical analyses, R software version 4.0.2 and various R packages were utilized, with a 2-tailed value of 0.05 indicating statistical significance. The “survival” package was used to run univariate and multivariate Cox regression analysis. The “glmnet” package was used to perform the LASSO-Cox regression analysis, and ten times cross-validations were employed to find the optimum penalty parameter lambda. The “survival” package was used to create Kaplan–Meier analyses and survival curves. The nomogram and calibration curve were created using the “rms” package. The “timeROC” software was used to analyze time-dependent ROC curves. Based on FDR, the Benjamini–Hochberg technique was employed to identify differentially expressed genes. Using “GSVA,” the ssGSEA-normalized DEGs were compared with a genome.

3. Results

3.1. An Overview of Genetic Alterations and Expression Variations in Pyroptosis-Related Genes in LUAD

In order to compare the pyroptosis-related gene expression levels found in the TCGA data from 19 normal and 414 tumor tissues, we first calculated the 52 pyroptosis-related gene expression levels and then looked for the DEGs that corresponded to them. 6 genes were found to be downregulated, while 23 genes were found to be enriched in the tumor group. Figure 1(a) shows heatmaps of these genes’ RNA levels. Figure 1(b) depicts the correlation network including all pyroptosis-related genes. At the genetic level, 337 of the 561 samples (approximately 60.07%) had mutations in pyroptosis-related genes (Figure 1(c)). We also discovered CNVs in pyroptosis-related genes, which were frequent alterations with a focus on copy number amplification (Figure 1(e)). We discovered changes in pyroptosis-related genes with CNVs on the chromosome (Figure 1(d)). We predicted that alterations in CNV might play a significant role in causing aberrant gene expression. Finding that pyroptosis-related genes had a relationship with the LUAD suggests that they may reflect unique features in patients in our research.

(a)

(b)

(c)

(d)

(e)

Figure 1

Genes associated with pyroptosis in LUAD have distinct characteristics and differences. (a) Heatmap depicting the variations in pyroptosis-related gene expression between the normal and tumor samples. (b) A network of correlations including pyroptosis-related genes in the TCGA cohort. (c) The mutation landscape in 561 LUAD patients from the TCGA-LUAD cohort. Each pyroptosis-related gene’s mutation information was displayed by a waterfall plot. Annotations at the bottom of the corresponding hues denoted various mutation types. The barplot above depicted the mutation load. Individual mutation frequencies were indicated by the right values. (d) The TCGA cohort’s location of CNV modification of genes associated with pyroptosis on chromosomes. (e) Pyroptosis-related gene CNV frequency in the TCGA cohort. The height of the columns revealed various proportions. Copy number variants, or CNVs, are a kind of copy number variation.

3.2. Tumor Classification Based on the DEGs

We did a consensus clustering analysis of all the LUAD cases in the TCGA cohort, with a focus on all the 29 pyroptosis-related DEGs, in order to investigate the links between those expression patterns and LUAD subtypes. In order to see whether two groups can be successfully distinguished using the DEGs, we increased the clustering variable (k) from 2 to 10. As we predicted, at k = 2, the intragroup correlations were low, suggesting that the LUAD patients could be well split into two clusters (Figure 2(a)). A heatmap depicts the gene expression profile as well as clinical characteristics such as stage, grade, age, and gender, with grade indicating a difference between the two groups (Figure 2(c)). The OS time of the two clusters was also examined, and a substantial difference was discovered (; Figure 2(b)).

(a)

(b)

(c)

3.3. Gene Signature Constructed from Pyroptosis-Related Clustering

On the basis of these findings, we created a gene signature that may be used in the treatment of LUAD and can be used to calculate an individual patient’s LUAD score. The gene signature was also assessed in order to diagnose and treat each patient on an individual basis. First, 1458 DEGs with |log₂FC| ≥ 1 and were shown to be linked with the two clusters (Figure S1; Table S2). 13 genes were identified as independent prognostic signatures by univariate regression analysis (Table S3). Applying the LASSO-Cox regression model with a minimum of λ, 7 of the 13 DEGs were maintained in order to construct a model capable of quantifying each patient. We utilized these to create a signature that was associated with pyroptosis (Figures 3(a) and 3(b)). The next step was to try to evaluate the signature’s worth by forecasting patient prognosis. In accordance with the median cut-off value, the patients were split into two groups (Figure 4(c)). Patients in various risk categories were dispersed in two directions according to PCA and t-SNE analyses (Figures 4(e) and 4(g)). Additionally, the scatter chart showed that individuals at high risk died sooner than those at low risk (Figure 4(a)). The Kaplan–Meier curve consistently indicated that high-risk patients had a substantially shorter OS than low-risk patients (Figure 3(d); ). For examination of the prognostic model’s survival prediction, time-dependent ROC curves were constructed, and the AUC reached 0.688 at one year, 0.665 at three years, and 0.688 at five years (Figure 3(c)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3

The development of a gene signature to predict patient OS based on clusters associated with pyroptosis. (a-b) Ten-fold cross-validation was used to determine the value of the super parameter in the TCGA LASSO-Cox model. (c) Time-dependent ROC analysis of the risk score in the TCGA cohort. (d) OS curves for the TCGA cohort’s various risk score subgroups. (e) Time-dependent ROC analysis of the risk score in the GEO cohort. (f) OS curves for the GEO cohort’s various risk score subgroups.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.4. Validation of the 7-Gene Signature in the GEO Cohort

Patients in the GEO cohort were also divided into high-risk and low-risk groups based on the median value from the TCGA cohort in order to assess the stability of the model built from the TCGA dataset. PCA and t-SNE analysis revealed a distinct distribution of patients in the two categories, similar to the TCGA cohort findings (Figures 4(f) and 4(h)). Additionally, individuals in the high-risk group died sooner (Figure 4(b)) and had a lower OS time when compared with those in the low-risk group (Figure 3(f)). Additionally, the AUC for the seven-gene signature was 0.743 at one year, 0.736 at three years, and 0.710 at five years (Figure 3(e)).

3.5. The 7-Gene Signature’s Independent Prognostic Value

The Cox analyses of factors, both univariate and multivariate, were used to evaluate whether or not the risk score was an independent prognostic predictor for OS. Using univariate Cox analysis, it was discovered that risk ratings in both the TCGA and GEO cohorts were substantially associated with OS (TCGA cohort: HR = 3.663, 95%CI = 2.303 − 5.826, ; GEO cohort: HR = 4.422, 95%CI = 2.364 − 8.271, ) (Figures 5(a) and 5(b)). Multivariate Cox analysis showed that the risk score remained a significant predictor when additional confounding factors were taken into consideration. (TCGA cohort: HR = 2.717, 95%CI = 2.175 − 5.677, ; GEO cohort: HR = 1.185 − 6.227, ) (Figures 5(c) and 5(d)). Furthermore, we created a heatmap of clinical characteristics for the TCGA cohort (Figure 5(e)) and discovered that gender, grade, stage, and so on were distributed differently across the low- and high-risk categories .

(a)

(b)

(c)

(d)

(e)

3.6. Developing and Validating a Nomogram That Incorporates Clinical Features

Because the risk score was strongly associated with high malignancy, we integrated clinical variables and created a nomogram in the TCGA cohort to expand the clinical applicability and use of the signature (Figure 6(a)). By combining the points for each prognostic criterion, each patient was given a total point value. Patients with a higher overall point had a poorer clinical outcome. Additionally, the calibration plot demonstrated that the nomogram operated in a manner consistent with an ideal model (Figure 6(b)).

(a)

(b)

3.7. Functional Analyses Based on the Risk Model

To investigate further the variations in gene functions and pathways across the risk model subgroups, we used the “limma” R package to extract DEGs using the FDR < 0.05 and |log₂FC| ≥ 1 criterion. According to the TCGA cohort, a total of 317 DEGs were found between those who were at low and high risk of LUAD. As a result, 161 genes were upregulated, whereas 156 genes were downregulated in the high-risk group (Table S4). Based on these DEGs, GO and KEGG analyses were carried out. The DEGs were mostly associated with immunological and cell differentiation signaling pathways according to the findings (Figures 7(a) and 7(b)).

(a)

(b)

3.8. Comparative Analysis of Immunological Activation among Subgroups

On the basis of the functional analyses, we used the single-sample gene set enrichment analysis (ssGSEA) to compare the enrichment scores for 16 different types of immune cells, as well as the activity of 13 different immune-related pathways, between the low-and high-risk groups in both the TCGA and GEO cohorts. The high-risk subgroups in the TCGA cohort (Figure 8(a)) exhibited higher levels of immune cell infiltration, particularly CD8+ T cells, macrophages, and T helper cells (Tfh, Th1, and Th2 cells) than the low-risk category. Except for the Type II IFN response pathway, the other 12 immune pathways showed more activity in the high-risk group than in the low-risk group during immune function analysis (Figure 8(c)). Similar results were found when evaluating the immunological state of the GEO cohort (Figures 8(b) and 8(d)).

(a)

(b)

(c)

(d)

4. Discussion

Lung cancer is widely recognized to be the leading cause of mortality around the globe. Lung cancer patients with NSCLC, of whom roughly 50% have LUAD, make up almost 80 percent of all lung cancer patients [25]. Despite advances in treatment regimens, the survival rate for LUAD patients remains low. The high-level variability of LUAD, as well as the complicated etiologic variables, makes prognostic prediction difficult. As a result, the development of new prognostic models is critical.

Pyroptosis is a kind of programmed cell death that happens in cells attacked by pathogens, triggering the body’s inflammatory response [15, 26]. Apoptosis may therefore be changed to pyroptosis in response to pathogen activation. Pyroptosis is involved in a variety of malignancies [27]. It inhibits tumor development in colorectal, liver, and skin cancers but has a bidirectional impact on breast cancer [28–31]. As a result, we cannot assess the predictive usefulness of LUAD only on the basis of the expression of various gasdermins. As a result, we investigated all pathways associated with pyroptosis and developed a predictive signature. Pyroptosis is now being explored for antitumor treatment, and our research suggests that pyroptosis in conjunction with immunotherapy may be a feasible therapeutic approach for patients with poor prognoses.

Classification of samples according to specified gene expression criteria is a well-established technique. On the basis of the expression of pyroptosis-related genes, we developed a subtyping approach for LUAD patients. We demonstrated a substantial association between the expression of these pyroptosis-related genes and various survival risks. Then, using the two clusters, we built a risk score model to quantify the prognostic risk. Our research established a solid case for LUAD clinical treatment. To begin, the risk score takes the patient’s heterogeneity into consideration. Second, this score may be used to establish a connection between pyroptosis and prognosis. A high-risk score was associated with poor clinical characteristics and a shorter projected survival time. The results on TME cell infiltration showed that the risk score is critical for immunotherapy. Most active immune cell infiltration resulted in a superior response to immunotherapy in individuals with a high score. A predictive model for LUAD under hypoxia or modified m6A conditions was investigated as was the immunoscore. A major emphasis of our study has been on the factors that directly contribute to tumor cell death and change the tumor microenvironment. As a consequence, our method is more advantageous in terms of facilitating treatment.

Seven genes were included in our model (CPT1C, PAQR8, CD109, GOLT1A, ADCY7, SERPINB2, and KRT1). CPT1C is a member of the CPT1 family of enzymes that catalyzes the acylation of fatty acids and their entrance into the mitochondria for oxidation. CPT1C has a role in the regulation of energy balance, ceramide metabolism, and the hypothalamic control of food intake [32]. CPT1C levels were recently implicated in the poor prognosis and metastatic development of human malignancies, which are intimately linked to fatty acid absorption and metabolism. CD109 is a glycoprotein with a glycosylphosphatidylinositol (GPI) anchor that is overexpressed in squamous cell carcinomas of the lung, esophagus, uterus, and oral cavity. CD109 inhibits TGF-β signaling in keratinocytes via direct modulation of receptor activation [33]. According to Kazuhiro, miR-378a-3p regulates tamoxifen sensitivity in breast cancer MCF-7 cells via targeting GOLT1A [34]. According to Li et al., ADCY7 promotes the development of acute myeloid leukemia and that inhibiting ADCY7 may be a new approach for treating leukemia [35]. SERPINB2 is a stress protein because it is highly expressed in activated monocytes and macrophages, as well as differentiating keratinocytes [36]. According to Han et al., the transcript level of KRT1 may serve as a possible prognostic biomarker in patients with melanoma [37]. PAQR8, one of the sex steroid hormones, is supposed to be potential prognostic biomarker of endometrial cancer [38].

Our research sought to categorize patients with LUAD, discover DEGs, develop a predictive model, and establish a connection between pyroptosis and patient prognosis. There are limits to this study despite our attempts to verify it from many perspectives and using multiple databases. Further, in vitro and in vivo testing of our model is needed to get a deeper understanding of how the risk score and cell pyroptosis are related. These have not only exacerbated the difficulties but also provided us with optimism, which motivates us to continue digging.

5. Conclusion

In addition to providing some hints for future study on the mechanism of pyroptosis-related genes, our model may also offer additional resources for a better understanding of immune cell-specific genes implicated in cancer control.

Abbreviations

DEGs:	Differentially expressed genes
TCGA:	The cancer genome atlas
GEO:	Gene expression omnibus
GO:	Gene ontology
BP:	Biological processes
MF:	Molecular function
CC:	Cellular components
KEGG:	Kyoto encyclopedia of genes and genomes
FDR:	False discovery rate
ssGSEA:	Single-sample gene set enrichment analysis (ssGSEA).

Data Availability

Publicly available datasets were used in this study. These data can be found in the gene expression omnibus (GEO) database and in the cancer genome atlas (TCGA) database.

Ethical Approval

TCGA and GEO belong to public datasets. The patients involved in the databases have obtained ethical approval. Users can download relevant data for free for research and publish relevant articles. Our study is based on open-source data, so there are no ethical issues and other conflicts of interest.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Zheng Dong and Lv Bian designed and analyzed the research study. Minglang Wang and Yilian Wang wrote and revised the manuscript. Luoqing Wang collected and analyzed the data, and all authors have read and approved the manuscript. Zheng Dong and Lv Bian contributed equally to this work.

Acknowledgments

The authors acknowledge TCGA and GEO database for providing the platform for uploading the meaningful datasets.

Supplementary Materials

Supplementary File Table S1: 52 genes associated with pyroptosis from prior reviews. Supplementary File Table S2: 1458 DEGs associated with cluster 1 and cluster 2. Supplementary File Table S3: 13 genes identified by univariate regression. Supplementary File Table S4: 317 DEGs between low- and high-risk groups in TCGA cohort. Supplementary File Figure S1: An overview of the differential gene expression between the two pyroptosis-related clusters in TCGA cohort. (Supplementary Materials)

References

F. Nasim, B. F. Sabath, and G. A. Eapen, “Lung cancer,” Medical Clinics of North America, vol. 103, no. 3, pp. 463–473, 2019.
View at: Publisher Site | Google Scholar
C. E. DeSantis, K. D. Miller, A. Goding Sauer, A. Jemal, and R. L. Siegel, “Cancer statistics for african Americans, 2019,” CA: A Cancer Journal for Clinicians, vol. 69, no. 3, pp. 211–233, 2019.
View at: Publisher Site | Google Scholar
G. Absenger, J. Terzic, and A. Bezan, “ASCO update: lung cancer,” Memo-Magazine of European Medical Oncology, vol. 10, no. 4, pp. 224–227, 2017.
View at: Publisher Site | Google Scholar
Y. Song, D. Chen, X. Zhang, Y. Luo, and S. Li, “Integrating genetic mutations and expression profiles for survival prediction of lung adenocarcinoma,” Thoracic Cancer, vol. 10, no. 5, pp. 1220–1228, 2019.
View at: Publisher Site | Google Scholar
F. Guerrera, L. Errico, A. Evangelista et al., “Exploring stage I non-small-cell lung cancer: development of a prognostic model predicting 5-year survival after surgical resection†,” European Journal of Cardio-Thoracic Surgery, vol. 47, no. 6, pp. 1037–1043, 2015.
View at: Publisher Site | Google Scholar
J.-P. Pignon, H. Tribodet, G. V. Scagliotti et al., “Lung adjuvant cisplatin evaluation: a pooled analysis by the LACE Collaborative Group,” Journal of Clinical Oncology, vol. 26, no. 21, pp. 3552–3559, 2008.
View at: Publisher Site | Google Scholar
G. M. Strauss, J. E. Herndon, M. A. Maddaus et al., “Adjuvant paclitaxel plus carboplatin compared with observation in stage IB non-small-cell lung cancer: CALGB 9633 with the cancer and leukemia group B, radiation therapy oncology group, and north central cancer treatment group study groups,” Journal of Clinical Oncology, vol. 26, no. 31, pp. 5043–5051, 2008.
View at: Publisher Site | Google Scholar
P. Broz, P. Pelegrín, and F. Shao, “The gasdermins, a protein family executing cell death and inflammation,” Nature Reviews Immunology, vol. 20, no. 3, pp. 143–157, 2020.
View at: Publisher Site | Google Scholar
J. Ding, K. Wang, W. Liu et al., “Pore-forming activity and structural autoinhibition of the gasdermin family,” Nature, vol. 535, no. 7610, pp. 111–116, 2016.
View at: Publisher Site | Google Scholar
C. Rogers, T. Fernandes-Alnemri, L. Mayes, D. Alnemri, G. Cingolani, and E. S. Alnemri, “Cleavage of DFNA5 by caspase-3 during apoptosis mediates progression to secondary necrotic/pyroptotic cell death,” Nature Communications, vol. 8, Article ID 14128, 2017.
View at: Publisher Site | Google Scholar
J. Shi, Y. Zhao, K. Wang et al., “Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death,” Nature, vol. 526, no. 7575, pp. 660–665, 2015.
View at: Publisher Site | Google Scholar
K. Yoneda, N. Imanishi, Y. Ichiki, and F. Tanaka, “Treatment of non-small cell lung cancer with EGFR-mutations,” Journal of UOEH, vol. 41, no. 2, pp. 153–163, 2019.
View at: Publisher Site | Google Scholar
T. Abu Hejleh, M. Furqan, Z. Ballas, and G. Clamon, “The clinical significance of soluble PD-1 and PD-L1 in lung cancer,” Critical Reviews in Oncology, vol. 143, pp. 148–152, 2019.
View at: Publisher Site | Google Scholar
R. Karki and T.-D. Kanneganti, “Diverging inflammasome signals in tumorigenesis and potential targeting,” Nature Reviews Cancer, vol. 19, no. 4, pp. 197–214, 2019.
View at: Publisher Site | Google Scholar
X. Xia, X. Wang, Z. Cheng et al., “The role of pyroptosis in cancer: pro-cancer or pro-“host”?” Cell Death & Disease, vol. 10, no. 9, Article ID 650, 2019.
View at: Publisher Site | Google Scholar
B. Wang and Q. Yin, “AIM2 inflammasome activation and regulation: a structural perspective,” Journal of Structural Biology, vol. 200, no. 3, pp. 279–282, 2017.
View at: Publisher Site | Google Scholar
S. M. Man and T.-D. Kanneganti, “Regulation of inflammasome activation,” Immunological Reviews, vol. 265, no. 1, pp. 6–21, 2015.
View at: Publisher Site | Google Scholar
J. A. H. A. Wong, “Algorithm AS 136: a K-means clustering algorithm,” Journal of the Royal Statistical Society, vol. 28, no. 1, pp. 100–108, 1979.
View at: Google Scholar
M. D. Wilkerson and D. N. Hayes, “ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking,” Bioinformatics, vol. 26, no. 12, pp. 1572-1573, 2010.
View at: Publisher Site | Google Scholar
N. Simon, J. Friedman, T. Hastie, and R. Tibshirani, “Regularization paths for cox’s proportional hazards model via coordinate descent,” Journal of Statistical Software, vol. 39, no. 5, pp. 1–13, 2011.
View at: Publisher Site | Google Scholar
R. Tibshirani, “The LASSO method for variable selection in the Cox Model,” Statistics in Medicine, vol. 16, no. 4, pp. 385–395, 1997.
View at: Publisher Site | Google Scholar
G. E. Hinton, “Visualizing high-dimensional data using t-SNE,” Vigiliae Christianae, vol. 9, pp. 2579–2605, 2008.
View at: Google Scholar
V. K. Mootha, C. M. Lindgren, K.-F. Eriksson et al., “PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes,” Nature Genetics, vol. 34, no. 3, pp. 267–273, 2003.
View at: Publisher Site | Google Scholar
A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences, vol. 102, no. 43, pp. 15545–15550, 2005.
View at: Publisher Site | Google Scholar
H. X. Dong, R. Wang, X. Y. Jin, J. Zeng, and J. Pan, “LncRNA DGCR5 promotes lung adenocarcinoma (LUAD) progression via inhibiting hsa‐mir‐22‐3p,” Journal of Cellular Physiology, vol. 233, no. 5, pp. 4126–4136, 2018.
View at: Publisher Site | Google Scholar
Y. Fang, S. Tian, Y. Pan et al., “Pyroptosis: a new frontier in cancer,” Biomedicine & pharmacotherapy = Biomedecine & pharmacotherapie, vol. 121, Article ID 109595, 2020.
View at: Publisher Site | Google Scholar
Y. Y. Wang, X. L. Liu, and R. Zhao, “Induction of pyroptosis and its implications in cancer management,” Frontiers in oncology, vol. 9, Article ID 971, 2019.
View at: Publisher Site | Google Scholar
Z. Tang, L. Ji, M. Han et al., “Pyroptosis is involved in the inhibitory effect of FL118 on growth and metastasis in colorectal cancer,” Life Sciences, vol. 257, Article ID 118065, 2020.
View at: Publisher Site | Google Scholar
J. Wu, S. Lin, B. Wan, B. Velani, and Y. Zhu, “Pyroptosis in liver disease: new insights into disease mechanisms,” Aging and Disease, vol. 10, no. 5, pp. 1094–1108, 2019.
View at: Publisher Site | Google Scholar
B. Zhou, J.-Y. Zhang, X.-S. Liu et al., “Tom20 senses iron-activated ROS signaling to promote melanoma cell pyroptosis,” Cell Research, vol. 28, no. 12, pp. 1171–1185, 2018.
View at: Publisher Site | Google Scholar
X. Wu, X. Mao, Y. Huang, Q. Zhu, J. Guan, and L. Wu, “Detection of proteins associated with the pyroptosis signaling pathway in breast cancer tissues and their significance,” International Journal of Clinical and Experimental Pathology, vol. 13, no. 6, pp. 1408–1414, 2020.
View at: Google Scholar
I. R. Schlaepfer and M. Joshi, “CPT1A-mediated fat oxidation, mechanisms, and therapeutic potential,” Endocrinology, vol. 161, no. 2, 2020.
View at: Publisher Site | Google Scholar
C. Li, M. A. Hancock, P. Sehgal, S. Zhou, D. P. Reinhardt, and A. Philip, “Soluble CD109 binds TGF-β and antagonizes TGF-β signalling and responses,” Biochemical Journal, vol. 473, no. 5, Article ID 537, 2016.
View at: Publisher Site | Google Scholar
K. Ikeda, K. Horie-Inoue, T. Ueno et al., “miR-378a-3p modulates tamoxifen sensitivity in breast cancer MCF-7 cells through targeting GOLT1A,” Scientific Reports, vol. 5, Article ID 13170, 2015.
View at: Publisher Site | Google Scholar
C. Li, J. Xie, Z. Lu et al., “ADCY7 supports development of acute myeloid leukemia,” Biochemical and Biophysical Research Communications, vol. 465, no. 1, pp. 47–52, 2015.
View at: Publisher Site | Google Scholar
D. R. Croucher, D. N. Saunders, S. Lobov, and M. Ranson, “Revisiting the biological roles of PAI2 (SERPINB2) in cancer,” Nature Reviews Cancer, vol. 8, no. 7, pp. 535–545, 2008.
View at: Publisher Site | Google Scholar
W. Han, C. Hu, Z. J. Fan, and G. L. Shen, “Transcript levels of keratin 1/5/6/14/15/16/17 as potential prognostic indicators in melanoma patients,” Scientific Reports, vol. 11, no. 1, Article ID 1023, 2021.
View at: Publisher Site | Google Scholar
M. Sinreih, T. Knific, P. Thomas, S. Frković Grazio, and T. L. Rižner, “Membrane progesterone receptors β and γ have potential as prognostic biomarkers of endometrial cancer,” The Journal of Steroid Biochemistry and Molecular Biology, vol. 178, pp. 303–311, 2018.
View at: Publisher Site | Google Scholar

Copyright

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

PDF Download Citation

Download other formats

Order printed copies

Views

898

Downloads

842

Citations