BackgroundLung adenocarcinoma (LUAD) with lymph node (LN) metastasis is linked to poor prognosis, yet the underlying mechanisms remain largely undefined. This study aimed to elucidate the immunogenomic landscape associated with LN metastasis in LUAD.MethodsWe employed broad-panel next-generation sequencing (NGS) on a cohort of 257 surgically treated LUAD patients to delineate the molecular landscape of primary tumors and identify actionable driver-gene alterations. Additionally, we used multiplex immunohistochemistry (mIHC) on a propensity score-matched cohort, which enabled us to profile the immune microenvironment of primary tumors in detail while preserving cellular metaclusters, interactions, and neighborhood functional units. By integrating data from NGS and mIHC, we successfully identified spatial immunogenomic patterns and developed a predictive model for LN metastasis, which was subsequently validated independently.ResultsOur analysis revealed distinct immunogenomic alteration patterns associated with LN metastasis stages. Specifically, we observed increased mutation frequencies in genes such as PIK3CG and ATM in LN metastatic primary tumors. Moreover, LN positive primary tumors exhibited a higher presence of macrophage and regulatory T cell metaclusters, along with their enriched neighborhood units (p < 0.05), compared to LN negative tumors. Furthermore, we developed a novel predictive model for LN metastasis likelihood, designed to inform non-surgical treatment strategies, optimize personalized therapy plans, and potentially improve outcomes for patients who are ineligible for surgery.ConclusionsThis study offers a comprehensive analysis of the genetic and immune profiles in LUAD primary tumors with LN metastasis, identifying key immunogenomic patterns linked to metastatic progression. The predictive model derived from these insights marks a substantial advancement in personalized treatment, underscoring its potential to improve patient management.