C6 and KLRG2 are pyroptosis subtype-related prognostic biomarkers and correlated with tumor-infiltrating lymphocytes in lung adenocarcinoma

Abstract

Abstract Pyroptosis plays an important role in lung adenocarcinoma (LUAD). In this study, we aimed to explore the pyroptosis-related gene (PRG) expression pattern and to identify promising pyroptosis-related biomarkers to improve the prognosis of LUAD. The gene expression profiles and clinical information of LUAD patients were downloaded from the Cancer Genome Atlas (TCGA), and validation cohort information was extracted from the Gene Expression Omnibus database. Gene expression data were analyzed using the limma package and visualized using the ggplot2 package as well as the pheatmap package in R software. Functional enrichment analysis was also performed for the 44 differentially expressed PRGs (DEPRGs). Then, consensus clustering revealed pyroptosis-related tumor subtypes, and differentially expressed genes (DEGs) were screened according to the subtypes. Next, univariate Cox and multivariate Cox regression analyses were used to identify independent prognostic PRGs. After overlapping DEGs and the Lasso regression analysis-based prognostic genes, the predictive risk model was established and validated. Correlation analysis between PRGs and clinicopathological variables was also explored. Finally, the TIMER and TISIDB databases were used to further explore the correlation analysis between immune cell infiltration levels, the risk score, and clinicopathological variables in the predictive risk model. A total of 52 genes from the PubMed were identified as PRGs, and 44 of the 52 genes were pooled as DEPRGs. The most significant GO term was “collagen trimer” (P = 2.46E-13), and KEGG analysis results indicated that 44 DEPRGs were significantly enriched in Salmonella infection (P < 0.001). Then, consensus clustering analysis divided LUAD patients into two clusters, and a total of 79 DEGs were identified according to these cluster subtypes. Subsequently, univariate and multivariate Cox regression analyses were used to identify 12 genes that could serve as independent prognostic indicators and we also performed Lasso regression analysis and screened 23 DEGs. After overlapping 23 DEGs and 12 genes, only 4 (KLRG2, MAPK4, C6 and SFRP5) of 12 genes were selected for the further exploration of the prognostic pattern. Survival analysis results indicated that this risk model effectively predicted the prognosis (P < 0.001). Combined with the correlation analysis results between the 4 genes and clinicopathological variables, C6 and KLRG2 were screened as prognostic genes. In this study, we constructed a predictive risk model and identified two pyroptosis subtype-related gene expression patterns to improve the prognosis of LUAD. Understanding the subtypes of LUAD is helpful for accurately characterizing the LUAD and developing personalized treatment.

Keywords

• Lung adenocarcinoma
• Pyroptosis
• C6
• KLRG2
• Biomarkers
• Risk model