Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the leading fatal malignancy with high morbidity and mortality worldwide. However, due to its complicated mechanism and lack of effective clinical therapeutics, early diagnosis and prognosis are still unsatisfactory. Most of the previous studies focused on cancer stem cells (CSCs), the relationship between cancer stemness (stem-like characteristics) and anti-tumor immunity has not been clearly revealed. Therefore, this study aimed to comprehensively analyze the role of cancer stemness and tumor microenvironment (TME) in LUAD using weighted gene co-expression network analysis (WGCNA). We constructed a gene co-expression network, identified key modules and hub genes, and further explored the relationship between hub gene expression and cancer immunological characteristics through a variety of algorithms, including Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and Gene Set Enrichment Analysis (GSEA). The hub genes were renamed stemness-immune related genes (SIRGs), whose functions were examined at the transcription and protein levels through survival analysis with additional samples, Oncomine database, immunohistochemistry, single cell RNA sequencing (scRNA-seq) and single-sample Gene Set Enrichment Analysis (ssGSEA). Subsequently, Tumor Immune Dysfunction and Exclusion (TIDE) and Connectivity Map (CMap) were implemented for treatment and prognosis analyses. As a result, 15 co-expressed SIRGs (CCNA2, CCNB1, CDC20, CDCA5, CDCA8, FEN1, KIF2C, KPNA2, MCM6, NUSAP1, RACGAP1. RRM2, SPAG5, TOP2A, TPX2) were identified. It was discovered that there was a general negative correlation between cancer stemness and immunity. The expression of SIRGs could greatly affect our prediction of tumor occurrence, progression, the efficacy of chemotherapy and immunotherapy, and clinical outcomes. In conclusion, the 15 SIRGs reported in our study may be used as potential candidate biomarkers for prognostic indicators and therapeutic targets after further validation.