Abstract

The molecular distinction between benign and malignant pulmonary nodules remains a significant diagnostic challenge. While genomic drivers are well studied, multiomic integration of the epigenetic–transcriptional landscape and its translation into noninvasive tools are lacking. We performed a multiomic characterization (genomic, epigenomic, and transcriptomic) of 158 pulmonary nodules. Unsupervised factor analysis integrated these layers to identify core regulatory axes. A 9-gene cell-free DNA (cfDNA) methylation classifier was developed and validated in blood and tissue cohorts. Genomic profiling revealed EGFR mutations (exclusive to malignant nodules) and MYC amplification as fundamental initiators of malignancy. Multiomic factor analysis (Factor 1) revealed profound genetic‒epigenetic synergy, in which these alterations dictate a permissive methylome, leading to aberrant epigenetic programming of chromatin accessibility, as well as epigenetic-transcriptional effects: hypomethylation at the promoters of cell cycle genes that augments their expression, and hypermethylation at immune related pathways gene loci that silences their transcription. This effect orchestrates formation of proproliferative (E2F target/G2M checkpoint) and “immune-cold” malignant phenotype, characterized by elevated Treg/CD8+ ratios and fibroblast recruitment. Notably, we observed a gradual accumulation of methylation aberrations along the premalignant-to-invasive continuum (adenocarcinoma in situ [AIS]→minimally invasive adenocarcinoma [MIA]→adenocarcinoma [ADC]), identifying progressive epigenetic dysregulation as a hallmark of tumor aggressiveness. Global methylome remodeling drives ADC progression through hypermethylation-mediated silencing of tumor suppressors (RASA3 and PPARG) and hypomethylation-activated oncogenic axes, specifically the GDF15 axis, which independently predict poor survival in patients with lung ADC in the TCGA cohort.