Cancer Stemness-Based Prognostic Immune-Related Gene Signatures in Lung Adenocarcinoma and Lung Squamous Cell Carcinoma

A novel prognostic signature based on cancer stemness and metabolism-related genes for cervical squamous cell carcinoma and endocervical adenocarcinoma

BACKGROUND

CESC is the second most commonly diagnosed gynecological malignancy. Given the pivotal involvement of metabolism-related genes (MRGs) in the etiology of multiple tumors, our investigation aims to devise a prognostic risk signature rooted in cancer stemness and metabolism.

METHODS

The stemness index based on mRNA expression (mRNAsi) of samples from the TCGA dataset was computed using the One-class logistic regression (OCLR) algorithm. Furthermore, potential metabolism-related genes related to mRNAsi were identified through weighted gene co-expression network analysis (WGCNA). We construct a stemness-related metabolic gene signature through shrinkage estimation and univariate analysis, thereby calculating the corresponding risk scores. Moreover, we selected corresponding DEGs between groups with high- and low-risk score and conducted routine bioinformatic analyses. Furthermore, we validated the expression of four hub genes at the protein level through immunohistochemistry (IHC) in samples obtained from our patient cohort.

RESULTS

According to the findings, it was found that six genes-AKR1B10, GNA15, ALDH1B1, PLOD2, LPCAT1, and GPX8- were differentially expressed in both TCGA-CSEC and GEO datasets among 23 differentially expressed metabolism-related genes (DEMRGs). mRNAsi exhibited a notable association with the extent of key oncogene mutation. The results showed that the AUC values for forecasting survival at 1, 3, and 5 years are 0.715, 0.689, and 0.748, individually. We observed a notable association between the risk score and different immune cell populations, along with enrichment in crucial signaling pathways in CESC. Four genes differentially expressed between different risk score groups were validated by IHC to be highly expressed in the CESC samples at the protein level.

CONCLUSION

The current investigation indicated that a 3-gene signature based on stemness-related metabolic and 4 hub genes with differential expression between high and low-risk score subgroups may serve as valuable prognostic markers and potential therapeutic targets in CESC.

Exercise potentially prevents colorectal cancer liver metastases by suppressing tumor epithelial cell stemness via RPS4X downregulation

Background

Colorectal cancer (CRC) is the third most prevalent tumor globally. The liver is the most common site for CRC metastasis, and the involvement of the liver is a common cause of death in patients with late-stage CRC. Consequently, mitigating CRC liver metastasis (CRLM) is key to improving CRC prognosis and increasing survival. Exercise has been shown to be an effective method of improving the prognosis of many tumor types. However, the ability of exercise to inhibit CRLM is yet to be thoroughly investigated.

Methods

The GSE157600 and GSE97084 datasets were used for analysis. A pan-cancer dataset which was uniformly normalized was downloaded and analyzed from the UCSC database: TCGA, TARGET, GTEx (PANCAN, n = 19,131, G = 60,499). Several advanced bioinformatics analyses were conducted, including single-cell sequencing analysis, correlation algorithm, and prognostic screen. CRC tumor microarray (TMA) as well as cell/animal experiments are used to further validate the results of the analysis.

Results

The greatest variability was found in epithelial cells from the tumor group. RPS4X was generally upregulated in all types of CRC, while exercise downregulated RPS4X expression. A lowered expression of RPS4X may prolong tumor survival and reduce CRC metastasis. RPS4X and tumor stemness marker-CD44 were highly positively correlated and knockdown of RPS4X expression reduced tumor stemness both in vitro and in vivo.

Conclusion

RPS4X upregulation may enhance CRC stemness and increase the odds of metastasis. Exercise may reduce CRC metastasis through the regulation of RPS4X.

The endoplasmic reticulum stress-related genes and molecular typing predicts prognosis and reveals characterization of tumor immune microenvironment in lung squamous cell carcinoma

BACKGROUND

Endoplasmic reticulum stress (ERS) acts critical roles on cell growth, proliferation, and metastasis in various cancers. However, the relationship between ERs and lung squamous cell carcinoma (LUSC) prognoses still remains unclear.

METHODS

The consensus clustering analysis of ERS-related genes and the differential expression analysis between clusters were investigated in LUSC based on TCGA database. Furthermore, ERS-related prognostic risk models were constructed by LASSO regression and Cox regression analyses. Then, the predictive effect of the risk model was evaluated by Kaplan-Meier, Cox regression, and ROC Curve analyses, as well as validated in the GEO cohort. According to the optimal threshold, patients with LUSC were divided into high- and low- risk groups, and somatic mutations, immune cell infiltration, chemotherapy response and immunotherapy effect were systematically analyzed.

RESULTS

Two ERS-related clusters were identified in patients with LUSC that had distinct patterns of immune cell infiltration. A 5-genes ERS-related prognostic risk model and nomogram were constructed and validated. Kaplan-Meier curves and Cox regression analysis showed that ERS risk score was an independent prognostic factor (p < 0.001, HR = 1.317, 95% CI = 1.159-1.496). Patients with low-risk scores presented significantly lower TIDE scores and significantly lower IC50 values for common chemotherapy drugs such as cisplatin and gemcitabine.

CONCLUSION

ERS-related risk signature has certain prognostic value and may be a potential therapeutic target and prognostic biomarker for LUSC patients.

Investigating the clinical role and prognostic value of genes related to insulin-like growth factor signaling pathway in thyroid cancer

BACKGROUND

Thyroid cancer (THCA) is the most common endocrine malignancy having a female predominance. The insulin-like growth factor (IGF) pathway contributed to the unregulated cell proliferation in multiple malignancies. We aimed to explore the IGF-related signature for THCA prognosis.

METHOD

The TCGA-THCA dataset was collected from the Cancer Genome Atlas (TCGA) for screening of key prognostic genes. The limma R package was applied for differentially expressed genes (DEGs) and the clusterProfiler R package was used for the Gene Ontology (GO) and KEGG analysis of DEGs. Then, the un/multivariate and least absolute shrinkage and selection operator (Lasso) Cox regression analysis was used for the establishment of RiskScore model. Receiver Operating Characteristic (ROC) analysis was used to verify the model's predictive performance. CIBERSORT and MCP-counter algorithms were applied for immune infiltration analysis. Finally, we analyzed the mutation features and the correlation between the RiskScore and cancer hallmark pathway by using the GSEA.

RESULT

We obtained 5 key RiskScore model genes for patient's risk stratification from the 721 DEGs. ROC analysis indicated that our model is an ideal classifier, the high-risk patients are associated with the poor prognosis, immune infiltration, high tumor mutation burden (TMB), stronger cancer stemness and stronger correlation with the typical cancer-activation pathways. A nomogram combined with multiple clinical features was developed and exhibited excellent performance upon long-term survival quantitative prediction.

CONCLUSIONS

We constructed an excellent prognostic model RiskScore based on IGF-related signature and concluded that the IGF signal pathway may become a reliable prognostic phenotype in THCA intervention.

Development and validation of a novel stemness-related prognostic model for neuroblastoma using integrated machine learning and bioinformatics analyses

Background

Neuroblastoma (NB) is a common solid tumor in children, with a dismal prognosis in high-risk cases. Despite advancements in NB treatment, the clinical need for precise prognostic models remains critical, particularly to address the heterogeneity of cancer stemness which plays a pivotal role in tumor aggressiveness and patient outcomes. By utilizing machine learning (ML) techniques, we aimed to explore the cancer stemness features in NB and identify stemness-related hub genes for future investigation and potential targeted therapy.

Methods

The public dataset GSE49710 was employed as the training set for acquire gene expression data and NB sample information, including age, stage, and MYCN amplification status and survival. The messenger RNA (mRNA) expression-based stemness index (mRNAsi) was calculated and patients were grouped according to their mRNAsi value. Stemness-related hub genes were identified from the differentially expressed genes (DEGs) to construct a gene signature. This was followed by evaluating the relationship between cancer stemness and the NB immune microenvironment, and the development of a predictive nomogram. We assessed the prognostic outcomes including overall survival (OS) and event-free survival, employing machine learning methods to measure predictive accuracy through concordance indices and validation in an independent cohort E-MTAB-8248.

Results

Based on mRNAsi, we categorized NB patients into two groups to explore the association between varying levels of stemness and their clinical outcomes. High mRNAsi was linked to the advanced International Neuroblastoma Staging System (INSS) stage, amplified MYCN, and elder age. High mRNAsi patients had a significantly poorer prognosis than low mRNAsi cases. According to the multivariate Cox analysis, the mRNAsi was an independent risk factor of prognosis in NB patients. After least absolute shrinkage and selection operator (LASSO) regression analysis, four key genes (ERCC6L, DUXAP10, NCAN, DIRAS3) most related to mRNAsi scores were discovered and a risk model was built. Our model demonstrated a significant prognostic capacity with hazard ratios (HR) ranging from 18.96 to 41.20, P values below 0.0001, and area under the receiver operating characteristic curve (AUC) values of 0.918 in the training set, suggesting high predictive accuracy which was further confirmed by external verification. Individuals with a low four-gene signature score had a favorable outcome and better immune responses. Finally, a nomogram for clinical practice was constructed by integrating the four-gene signature and INSS stage.

Conclusions

Our findings confirm the influence of CSC features in NB prognosis. The newly developed NB stemness-related four-gene signature prognostic signature could facilitate the prognostic prediction, and the identified hub genes may serve as promising targets for individualized treatments.

Identification of molecular subtypes and prognostic model to reveal immune infiltration and predict prognosis based on immunogenic cell death-related genes in lung adenocarcinoma

Immunogenic cell death (ICD) has been increasingly indicated to be related to caners. However, ICD's role in Lung adenocarcinoma (LUAD) is still not well investigated. Clinical data along with associated mRNA expression profiles from LUAD cases were collected in TCGA and GEO databases. 13 ICD-related genes were identified. Relations of ICD-related genes expression with prognosis of patients, tumor immune microenvironment (TIME) was analyzed. Then, candidate genes were identified and the prognostic signature were constructed. Afterwards, one nomogram incorporating those chosen clinical data together with risk scores were built. Finally, the effect of HSP90AA1, one gene of the prognostic signature, on LUAD cell were analyzed. Two clusters were identified, which were designated as the ICD-high or -low subtype according to ICD-related genes levels. ICD-high subgroup showed good prognosis, high immune cell infiltration degrees, and enhanced immune response signaling activity compared with ICD-low subtype. Moreover, we established and verified the risk signature based on ICD-related genes. High risk group predicted poor prognosis of LUAD independently and presented negative association with immune score and immune status. Furthermore, nomogram contributed to the accurate prediction of LUAD prognostic outcome. Finally, HSP90AA1 levels were remarkably elevated within tumor cells in comparison with healthy pulmonary epithelial cells. HSP90α, HSP90AA1 protein product, promoted growth, migration, and invasion of LUAD cells. Molecular subtypes and prognostic model were identified by incorporating ICD-related genes, and it was related to TIME and might be adopted for the accurate prediction of LUAD prognosis.

Identification of specific prognostic markers for lung squamous cell carcinoma based on tumor progression, immune infiltration, and stem index

Introduction

Lung squamous cell carcinoma (LUSC) is a unique subform of nonsmall cell lung cancer (NSCLC). The lack of specific driver genes as therapeutic targets leads to worse prognoses in patients with LUSC, even with chemotherapy, radiotherapy, or immune checkpoint inhibitors. Furthermore, research on the LUSC-specific prognosis genes is lacking. This study aimed to develop a comprehensive LUSC-specific differentially expressed genes (DEGs) signature for prognosis correlated with tumor progression, immune infiltration,and stem index.

Methods

RNA sequencing data for LUSC and lung adenocarcinoma (LUAD) were extracted from The Cancer Genome Atlas (TCGA) data portal, and DEGs analyses were conducted in TCGA-LUSC and TCGA-LUAD cohorts to identify specific DEGs associated with LUSC. Functional analysis and protein-protein interaction network were performed to annotate the roles of LUSC-specific DEGs and select the top 100 LUSC-specific DEGs. Univariate Cox regression and least absolute shrinkage and selection operator regression analyses were performed to select prognosis-related DEGs.

Results

Overall, 1,604 LUSC-specific DEGs were obtained, and a validated seven-gene signature was constructed comprising FGG, C3, FGA, JUN, CST3, CPSF4, and HIST1H2BH. FGG, C3, FGA, JUN, and CST3 were correlated with poor LUSC prognosis, whereas CPSF4 and HIST1H2BH were potential positive prognosis markers in patients with LUSC. Receiver operating characteristic analysis further confirmed that the genetic profile could accurately estimate the overall survival of LUSC patients. Analysis of immune infiltration demonstrated that the high risk (HR) LUSC patients exhibited accelerated tumor infiltration, relative to low risk (LR) LUSC patients. Molecular expressions of immune checkpoint genes differed significantly between the HR and LR cohorts. A ceRNA network containing 19 lncRNAs, 50 miRNAs, and 7 prognostic DEGs was constructed to demonstrate the prognostic value of novel biomarkers of LUSC-specific DEGs based on tumor progression, stemindex, and immune infiltration. In vitro experimental models confirmed that LUSC-specific DEG FGG expression was significantly higher in tumor cells and correlated with immune tumor progression, immune infiltration, and stem index. In vitro experimental models confirmed that LUSC-specific DEG FGG expression was significantly higher in tumor cells and correlated with immune tumor progression, immune infiltration, and stem index.

Conclusion

Our study demonstrated the potential clinical implication of the 7- DEGs signature for prognosis prediction of LUSC patients based on tumor progression, immune infiltration, and stem index. And the FGG could be an independent prognostic biomarker of LUSC promoting cell proliferation, migration, invasion, THP-1 cell infiltration, and stem cell maintenance.

Immunologic constant of rejection as a predictive biomarker of immune checkpoint inhibitors efficacy in non-small cell lung cancer

BACKGROUND

Anti-PD1/PDL1 immune checkpoint inhibitors (ICI) transformed the prognosis of patients with advanced non-small cell lung cancer (NSCLC). However, the response rate remains disappointing and toxicity may be life-threatening, making urgent identification of biomarkers predictive for efficacy. Immunologic Constant of Rejection signature (ICR) is a 20-gene expression signature of cytotoxic immune response with prognostic value in some solid cancers. Our objective was to assess its predictive value for benefit from anti-PD1/PDL1 in patients with advanced NSCLC.

METHODS

We retrospectively profiled 44 primary tumors derived from NSCLC patients treated with ICI as single-agent in at least the second-line metastatic setting. Transcriptomic analysis was performed using the nCounter® analysis system and the PanCancer Immune Profiling Panel. We then pooled our data with clinico-biological data from four public gene expression data sets, leading to a total of 162 NSCLC patients treated with single-agent anti-PD1/PDL1. ICR was applied to all samples and correlation was searched between ICR classes and the Durable Clinical Benefit (DCB), defined as stable disease or objective response according to RECIST 1.1 for a minimum of 6 months after the start of ICI.

RESULTS

The DCB rate was 29%; 22% of samples were classified as ICR1, 30% ICR2, 22% ICR3, and 26% ICR4. These classes were not associated with the clinico-pathological variables, but showed enrichment from ICR1 to ICR4 in quantitative/qualitative markers of immune response. ICR2-4 class was associated with a 5.65-fold DCB rate when compared with ICR1 class. In multivariate analysis, ICR classification remained associated with DCB, independently from PDL1 expression and other predictive immune signatures. By contrast, it was not associated with disease-free survival in 556 NSCLC TCGA patients untreated with ICI.

CONCLUSION

The 20-gene ICR signature was independently associated with benefit from anti-PD1/PDL1 ICI in patients with advanced NSCLC. Validation in larger retrospective and prospective series is warranted.

A Comprehensive Pan-Cancer Analysis of the Regulation and Prognostic Effect of Coat Complex Subunit Zeta 1

The Coatomer protein complex Zeta 1 (COPZ1) has been reported to play an essential role in maintaining the survival of some types of tumors. In this study, we sought to explore the molecular characteristics of COPZ1 and its clinical prognostic value through a pan-cancers bioinformatic analysis. We found that COPZ1 was extremely prevalent in a variety of cancer types, and high expression of COPZ1 was linked to poor overall survival in many cancers, while low expression in LAML and PADC was correlated with tumorigenesis. Besides, the CRISPR Achilles' knockout analysis revealed that COPZ1 was vital for many tumor cells' survival. We further demonstrated that the high expression level of COPZ1 in tumors was regulated in multi-aspects, including abnormal CNV, DNA-methylation, transcription factor and microRNAs. As for the functional exploration of COPZ1, we found a positive relationship between COPZ1's expression and stemness and hypoxia signature, especially the contribution of COPZ1 on EMT ability in SARC. GSEA analysis revealed that COPZ1 was associated with many immune response pathways. Further investigation demonstrated that COPZ expression was negatively correlated with immune score and stromal score, and low expression of COPZ1 has been associated to more antitumor immune cell infiltration and pro-inflammatory cytokines. The further analysis of COPZ1 expression and anti-inflammatory M2 cells showed a consistent result. Finally, we verified the expression of COPZ1 in HCC cells, and proved its ability of sustaining tumor growth and invasion with biological experiments. Our study provides a multi-dimensional pan-cancer analysis of COPZ and demonstrates that COPZ1 can serve as both a prospective target for the treatment of cancer and a prognostic marker for a variety of cancer types.

Construction of a novel mRNAsi-related risk model for predicting prognosis and immunotherapy response in osteosarcoma

Background

Targeting cancer stem cells (CSC) may represent a future therapeutic direction for osteosarcoma (OS), which mainly relies on the identification of CSC markers. This study aimed to classify OS based on messenger ribonucleic acid (mRNA) stemness indices (mRNAsi) and construct a mRNAsi-related risk model to predict the prognosis of OS.

Methods

The one-class logistic regression (OCLR) algorithm was applied to the RNA- sequencing (seq) data of human embryonic stem cells (hESC) and induced pluripotent stem cell (iPSC) lines to calculate mRNAsi. Weighted gene co-expression network analysis (WGCNA) was performed on data obtained from the TARGET database to screen the mRNAsi-related genes. Univariate Cox regression analysis was implemented to screen mRNAsi-related genes with prognostic significance for consensus clustering of OS. The least absolute shrinkage and selection operator (LASSO) and COX regression analysis were conducted to construct a risk model based on mRNAsi-related genes.

Results

Six gene modules were identified in the TARGET database. The yellow module showed the strongest negative correlation with mRNAsi and the strongest significant positive correlation with the immune score and stromal score. OS was divided into three molecular subtypes with significant survival differences based on 73 mRNAsi-related genes with prognostic value for OS. The survival rate was ranked as C3 < C1 < C2 from low to high. The levels of immune components in C2 was significantly higher than those in C1 and C3. HSD11B2, GBP1, RNF130, APBB1IP, and NPC2 in the yellow module were used as variables for building the mRNAsi-related risk model. The survival rate of the high-risk group (as defined by this model) was significantly higher than that of the low-risk group, and it had significant survival prediction ability in 28 types of cancer. In addition, the mRNAsi-related risk model was superior to the Tumor Immune Dysfunction and Exclusion (TIDE) model in predicting the prognosis and immunotherapy response in all three immunotherapy cohorts.

Conclusions

This study classified OS and constructed a mRNAsi-related risk model based on mRNAsi-related genes, which provides a potential tool for more accurate risk stratification of OS and prediction of immunotherapy response.

Identification and validation of molecular subtype and prognostic signature for lung adenocarcinoma based on neutrophil extracellular traps

Background: Neutrophil Extracellular Traps (NETs) are fibrous networks made of DNA-histone complexes and proteins protruded from activated neutrophils. Accumulating evidences have highlighted the vital role of NETs in tumor progression and diffusion. However, limited systematic studies regarding the role of NETs in LUAD have been performed. Methods: Differentially expressed NETs-related genes and their mutation landscape were identified with TCGA database. Consensus clustering analysis was performed to determine the NETs-related subtypes of LUAD. LASSO algorithm was employed to construct a prognostic signature. Moreover, GSE30219 and GSE31210 were used as independent validation. We also constructed a lncRNA-miRNA-mRNA regulatory axis with several miRNA and lncRNA databases. Results: Consensus clustering identified two NETs-related clusters in LUAD. High NETs score was correlated with a favorable overall survival, abundant immune cell infiltration, and high activity of immune response signal pathways. Six NET-related genes (G0S2, KCNJ15, S100A12, AKT2, CTSG, and HMGB1) with significant prognostic value were screened to develop a prognostic signature. LUAD patients with low-risk had a significantly favorable overall survival both in the training set and validation set. Moreover, NETs-related risk score and clinical stage could act as an independent prognostic factor for LUAD patients. Significant correlation was obtained between risk score and tumor immune microenvironment. We also identified lncRNA BCYRN1/miR-3664-5p/CTSG regulatory axis that may be involved in the progression of LUAD. Conclusion: We developed two molecular subtypes and a prognostic signature for LUAD based on NETs-related genes. This stratification could provide more evidences for estimating the prognosis and immunotherapy of LAUD patients.

Immune checkpoints expression patterns in early-stage triple-negative breast cancer predict prognosis and remodel the tumor immune microenvironment

Background

Currently, targeting immune checkpoint molecules holds great promise for triple-negative breast cancer (TNBC). However, the expression landscape of immune checkpoint genes (ICGs) in TNBC remains largely unknown.

Method

Herein, we systematically investigated the ICGs expression patterns in 422 TNBC samples. We evaluated the ICGs molecular typing based on the ICGs expression profile and explored the associations between ICGs molecular subtypes and tumor immune characteristics, clinical significance, and response to immune checkpoint inhibitors (ICIs).

Results

Two ICGs clusters and two ICGs-related gene clusters were determined, which were involved in different survival outcomes, biological roles and infiltration levels of immune cells. We established a quantification system ICGs riskscore (named IRS) to assess the ICGs expression patterns for individuals. TNBC patients with lower IRS were characterized by increased immune cell infiltration, favorable clinical outcomes and high sensitivity to ICIs therapy. We also developed a nomogram model combining clinicopathological variables to predict overall survival in TNBC. Genomic feature analysis revealed that high IRS group presented an increased tumor mutation burden compared with the low IRS group.

Conclusion

Collectively, dissecting the ICGs expression patterns not only provides a new insight into TNBC subtypes but also deepens the understanding of ICGs in the tumor immune microenvironment.

Immunotherapy response and microenvironment provide biomarkers of immunotherapy options for patients with lung adenocarcinoma

Background: Immunotherapy has been a promising approach option for lung cancer.

Method: All the open-accessed data was obtained from the Cancer Genome Atlas (TCGA) database. All the analysis was conducted using the R software analysis.

Results: Firstly, the genes differentially expressed in lung cancer immunotherapy responders and non-responders were identified. Then, the lung adenocarcinoma immunotherapy-related genes were determined by LASSO logistic regression and SVM-RFE, respectively. A total of 18 immunotherapy response-related genes were included in our investigation. Subsequently, we constructed the logistics score model. Patients with high logistics score had a better clinical effect on immunotherapy, with 63.2% of patients responding to immunotherapy, while only 12.1% of patients in the low logistics score group responded to immunotherapy. Moreover, we found that pathways related to immunotherapy were mainly enriched in metabolic pathways such as fatty acid metabolism, bile acid metabolism, oxidative phosphorylation, and carcinogenic pathways such as KRAS signaling. Logistics score was positively correlated with NK cells activated, Mast cells resting, Monocytes, Macrophages M2, dendritic cells resting, dendritic cells activated and eosinophils, while was negatively related to Tregs, macrophages M0, macrophages M1, and mast cells activated. In addition, ERVH48-1 was screened for single-cell exploration. The expression of ERVH48-1 increased in patients with distant metastasis, and ERVH48-1 was associated with pathways such as pancreas beta cells, spermatogenesis, G2M checkpoints and KRAS signaling. The result of quantitative real-time PCR showed that ERVH48-1 was upregulated in lung cancer cells.

Conclusion: Our study developed an effective signature to predict the immunotherapy response of lung cancer patients.

Integrated Bioinformatics Analysis Identifies a New Stemness Index-Related Survival Model for Prognostic Prediction in Lung Adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is one of the most lethal malignancies and is currently lacking in effective biomarkers to assist in diagnosis and therapy. The aim of this study is to investigate hub genes and develop a risk signature for predicting prognosis of LUAD patients.

METHODS

RNA-sequencing data and relevant clinical data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes associated with mRNA expression-based stemness indices (mRNAsi) in TCGA. We utilized LASSO Cox regression to assemble our predictive model. To validate our predictive model, me applied it to an external cohort.

RESULTS

mRNAsi index was significantly associated with the tissue type of LUAD, and high mRNAsi scores may have a protective influence on survival outcomes seen in LUAD patients. WGCNA indicated that the turquoise module was significantly correlated with the mRNAsi. We identified a 9-gene signature (CENPW, MCM2, STIL, RACGAP1, ASPM, KIF14, ANLN, CDCA8, and PLK1) from the turquoise module that could effectively identify a high-risk subset of these patients. Using the Kaplan-Meier survival curve, as well as the time-dependent receiver operating characteristic (tdROC) analysis, we determined that this gene signature had a strong predictive ability (AUC = 0.716). By combining the 9-gene signature with clinicopathological features, we were able to design a predictive nomogram. Finally, we additionally validated the 9-gene signature using two external cohorts from GEO and the model proved to be of high value.

CONCLUSION

Our study shows that the 9-gene mRNAsi-related signature can predict the prognosis of LUAD patient and contribute to decisions in the treatment and prevention of LUAD patients.