Brain metastases frequently develop in patients with non-small cell lung cancer (NSCLC) and are a common cause of cancer-related deaths, yet our understanding of the underlying human biology is limited. Here we performed multimodal single-nucleus RNA and T cell receptor, single-cell spatial and whole-genome sequencing of brain metastases and primary tumors of patients with treatment-naive NSCLC. Chromosomal instability (CIN) is a distinguishing genomic feature of brain metastases compared with primary tumors, which we validated through integrated analysis of molecular profiling and clinical data in 4,869 independent patients, and a new cohort of 12,275 patients with NSCLC. Unbiased analyses revealed transcriptional neural-like programs that strongly enriched in cancer cells from brain metastases, including a recurring, CINhigh cell subpopulation that preexists in primary tumors but strongly enriched in brain metastases, which was also recovered in matched single-cell spatial transcriptomics. Using multiplexed immunofluorescence in an independent cohort of treatment-naive pairs of primary tumors and brain metastases from the same patients with NSCLC, we validated genomic and tumor-microenvironmental findings and identified a cancer cell population characterized by neural features strongly enriched in brain metastases. This comprehensive analysis provides insights into human NSCLC brain metastasis biology and serves as an important resource for additional discovery.