Circulating tumor cells' (CTCs) heterogeneity contributes to counteract their introduction in clinical practice. Through single-cell sequencing we aim at exploring CTC heterogeneity in metastatic breast cancer (MBC) patients. Single CTCs were isolated using DEPArray NxT. After whole genome amplification, libraries were prepared for copy number aberration (CNA) and single nucleotide variant (SNV) analysis and sequenced using Ion GeneStudio S5 and Illumina MiSeq, respectively. CTCs demonstrate distinctive mutational signatures but retain molecular traces of their common origin. CNA profiling identifies frequent aberrations involving critical genes in pathogenesis: gains of 1q (CCND1) and 11q (WNT3A), loss of 22q (CHEK2). The longitudinal single-CTC analysis allows tracking of clonal selection and the emergence of resistance-associated aberrations, such as gain of a region in 12q (CDK4). A group composed of CTCs from different patients sharing common traits emerges. Further analyses identify losses of 15q and enrichment of terms associated with pseudopodium formation as frequent and exclusive events. CTCs from MBC patients are heterogeneous, especially concerning their mutational status. The single-cell analysis allows the identification of aberrations associated with resistance, and is a candidate tool to better address treatment strategy. The translational significance of the group populated by similar CTCs should be elucidated.
Keywords: circulating tumor cells; copy number aberrations; liquid biopsy; metastatic breast cancer; next generation sequencing; single cell sequencing; single nucleotide variants.