Abstract
Multiple pulmonary lung cancers (MPLCs) are frequently encountered on computed tomography (CT) scanning of chest, yet their intrinsic characteristics associated with genomic features and radiological or pathological textures that may lead to distinct clinical outcomes remain largely unexplored. A total of 27 pulmonary nodules covering different radiological or pathological textures as well as matched adjacent normal tissues and blood samples were collected from patients diagnosed with MPLCs. Whole-exome sequencing (WES) and whole-transcriptome sequencing were performed. The molecular and immune features of MPLCs associated with distinct radiological or pathological textures were comprehensively investigated. Genomics analysis unveiled the distinct branches of pulmonary nodules originating independently within the same individual. EGFR and KRAS mutations were found to be prevalent in MPLCs, exhibiting mutual exclusivity. The group with KRAS mutations exhibited stronger immune signatures compared to the group with EGFR mutations. Additionally, MPLCs exhibited a pronounced immunosuppressive microenvironment, which was particularly distinct when compared with normal tissues. The expression of the FDSCP gene was specifically observed in MPLCs. When categorizing MPLCs based on radiological or pathological characteristics, a progressive increase in mutation accumulation was observed, accompanied by heightened chromatin-level instability as ground-glass opacity component declined or invasive progression occurred. A close association with the immunosuppressive microenvironment was also observed during the progression of pulmonary nodules. Notably, the upregulation of B cell and regulatory T cell marker genes occurred progressively. Immune cell abundance analysis further demonstrated a marked increase in exhausted cells and regulatory T cells during the progression of pulmonary nodules. These results were further validated by independent datasets including nCounter RNA profiling, single-cell RNA sequencing, and spatial transcriptomic datasets. Our study provided a comprehensive representation of the diverse landscape of MPLCs originating within the same individual and emphasized the significant influence of the immunosuppressive microenvironment in the occurrence and development of pulmonary nodules. These findings hold great potential for enhancing the clinical diagnosis and treatment strategies for MPLCs.