High-resolution transcriptomics analysis of CXCL13⁺ EPSTI1⁺ CDK1⁺ cells with a specific focus on lung adenocarcinoma

Programmed cell death ligand 1 (PD-L1) blocking therapy has transformed the treatment of lung adenocarcinoma (LUAD), which has significantly changed the landscape of immunotherapy. We aimed to explore specific cell subpopulations to understand tumor progression and identify markers of response to PD-L1 blocking therapy.

Bulk, fluorescence-activated cell sorting (FACS), and single-cell RNA (scRNA) sequencing were used to profile CXCL13, EPSTI1, and CDK1. The gene set variation analysis (GSVA) R package was utilized for score calculation, and prognostic analyses included receiver operating characteristic (ROC) curves, Cox proportional hazard models, and meta-analysis. Additionally, we analyzed tumor microenvironment (TME), genomics, compound perturbations, and clinical indicators. The high-dimensional analysis captured the intrinsic characteristics of the subpopulation. Furthermore, subpopulation differential genes were used for enrichment analysis of transcription factors and compounds.

Literature and website analyses supported the essential role of CXCL13, CDK1, and EPSTI1 in immunotherapy. This led us to focus specifically on LUAD by representing a pan-cancer profile of immune-sensitive genes. Logically, the high-characteristic population may consist of samples positive for CXCL13, EPSTI1, and CDK1. The three-gene signature was a favorable indicator of immunotherapy response in the Stand Up to Cancer-Mark Foundation (SU2C-MARK) LUAD cohort but showed a poor prognosis before treatment in the Lung Cancer Explorer (LCE) database. Further mechanistic exploration revealed specific mutations associated with the three-gene signature in SU2C-MARK LUAD, such as STK11. In The Cancer Genome Atlas (TCGA)-LUAD cohort, the high-scoring group exhibited a higher tumor mutational burden (TMB) and global methylation but a lower fraction genome altered (FGA) and estimated tumor purity. Moreover, dasatinib demonstrated sensitivity in the high-scoring group. The co-localization of the CXCL13, EPSTI1, and CDK1 subpopulation was validated through spatial transcriptome and immunohistochemical databases. Assessment of the subpopulation depicted high-resolution intercellular communication. Maintenance of specific pathways, such as TNF, CD74, and CD44, contributed to immunotherapy sensitivity. Finally, the subpopulation-enriched targets and drugs were confirmed through ConnectivityMap (CMAP) analysis and multi-omics, respectively.

In this study, positive samples for CXCL13, EPSTI1, and CDK1 exhibited poor prognostic significance in treatment-naïve LUAD cases but demonstrated benefits from PD-L1 blockade and dasatinib therapies.