Integrated bulk and single-cell RNA sequencing identifies an aneuploidy-based gene signature to predict sensitivity of lung adenocarcinoma to traditional chemotherapy drugs and patients' prognosis.

Background: Patients with lung adenocarcinoma (LUAD) often develop a poor prognosis. Currently, researches on prognostic and immunotherapeutic capacity of aneuploidy-related genes in LUAD are limited. Methods: Genes related to aneuploidy were screened based on bulk RNA sequencing data from public databases using Spearman method. Next, univariate Cox and Lasso regression analyses were performed to establish an aneuploidy-related riskscore (ARS) model. Results derived from bioinformatics analysis were further validated using cellular experiments. In addition, typical LUAD cells were identified by subtype clustering, followed by SCENIC and intercellular communication analyses. Finally, ESTIMATE, ssGSEA and CIBERSORT algorithms were employed to analyze the potential relationship between ARS and tumor immune environment. Results: A five-gene ARS signature was developed. These genes were abnormally high-expressed in LUAD cell lines, and in particular the high expression of CKS1B promoted the proliferative, migratory and invasive phenotypes of LUAD cell lines. Low ARS group had longer overall survival time, higher degrees of inflammatory infiltration, and could benefit more from receiving immunotherapy. Patients in low ASR group responded more actively to traditional chemotherapy drugs (Erlotinib and Roscovitine). The scRNA-seq analysis annotated 17 cell subpopulations into seven cell clusters. Core transcription factors (TFs) such as CREB3L1 and CEBPD were enriched in high ARS cell group, while TFs such as BCLAF1 and UQCRB were enriched in low ARS cell group. CellChat analysis revealed that high ARS cell groups communicated with immune cells Conclusion: In this research, integrative analysis based on the ARS model provided a potential direction for improving the diagnosis and treatment of LUAD.