Coexisting genomic aberrations associated with lymph node metastasis in breast cancer

Single cancer cell-sequencing studies currently use randomly selected cells, limiting correlations among genomic aberrations, morphology, and spatial localization. We laser-captured microdissected single cells from morphologically distinct areas of primary breast cancer and corresponding lymph node metastasis and performed whole-exome or deep-target sequencing of more than 100 such cells. Two major subclones coexisted in different areas of the primary tumor, and the lymph node metastasis originated from a minor subclone in the invasive front of the primary tumor, with additional copy number changes, including chr8q gain, but no additional point mutations in driver genes. Lack of metastasis-specific driver events led us to assess whether other clonal and subclonal genomic aberrations preexisting in primary tumors contribute to lymph node metastasis. Gene mutations and copy number variations analyzed in 5 breast cancer tissue sample sets revealed that copy number variations in several genomic regions, including areas within chr1p, chr8q, chr9p, chr12q, and chr20q, harboring several metastasis-associated genes, were consistently associated with lymph node metastasis. Moreover, clonal expansion was observed in an area of morphologically normal breast epithelia, likely driven by a driver mutation and a subsequent amplification in chr1q. Our study illuminates the molecular evolution of breast cancer and genomic aberrations contributing to metastases.