Constructing a Prognostic Gene Signature for Lung Adenocarcinoma Based on Weighted Gene Co-Expression Network Analysis and Single-Cell Analysis

Biqian Fu,* Lin Lu,* Haifu Huang Internal Medicine-Oncology, Shenzhen Hospital of Guangzhou University of Traditional Chinese Medicine, Shenzhen, People’s Republic of China*These authors contributed equally to this workCorrespondence: Haifu Huang, Internal Medicine-Oncology, Shenzhen Hospital of Guangzhou University of traditional Chinese Medicine, 6001 North Ring Road, Futian District, Shenzhen, 440300, People’s Republic of China, Tel +0755-83548611, Email huanghaifu2004@163.comPurpose: Lung adenocarcinoma (LUAD) has a high degree of intratumor heterogeneity. Advanced single-cell RNA sequencing (scRNA-seq) technologies have offered tools to analyze intratumor heterogeneity, which improves the accuracy of identifying biomarkers based on single-cell expression data, and thus helps in predicting prognosis of cancer patients and assisting decision-makings for cancer treatment.Patients and Methods: ScRNA-seq data containing two LUAD and two para-cancerous tissue samples were included to identify different cell clusters in tumor tissues. To identify the most relevant modules and important cell subpopulations (clusters) in LUAD tissues, weighted gene co-expression network analysis (WGCNA) was performed. Subsequently, LUAD molecular subtypes were constructed by unsupervised consensus clustering based on genes in key modules. Using differential analysis, univariate Cox regression analysis, and least absolute shrinkage and selection operator (LASSO) regression analysis, a prognostic model of LUAD was established.Results: A total of 14 cell clusters belonging to 10 cell types in LUAD were identified. The turquoise module was the most relevant to LUAD among all the modules; cluster 10 (C10, lung epithelial cells) was found to be the most strongly associated with the turquoise module. LUAD samples were divided into two groups of distinct molecular subtypes. Based on the 165 shared genes between the turquoise module and C10, 511 DEGs between the two molecular subtypes were obtained, and five of them were selected to construct the gene signature, which was validated to be an independent prognostic marker of LUAD.Conclusion: Fourteen cell clusters co-existed in LUAD, which contributed to its intratumor heterogeneity. Two molecular subtypes of LUAD were identified and a five-gene signature was developed and validated to be significantly associated with prognostic and clinical characteristics of LUAD patients.Keywords: single-cell RNA sequencing, lung adenocarcinoma, intratumor heterogeneity, molecular subtypes, prognosis, five-gene signature