Archival single-cell genomics reveals persistent subclones during DCIS progression.

Ductal carcinoma in situ (DCIS) is a common precursor of invasive breast cancer. Our understanding of its genomic progression to recurrent disease remains poor, partly due to challenges associated with the genomic profiling of formalin-fixed paraffin-embedded (FFPE) materials. Here, we developed Arc-well, a high-throughput single-cell DNA-sequencing method that is compatible with FFPE materials. We validated our method by profiling 40,330 single cells from cell lines, a frozen tissue, and 27 FFPE samples from breast, lung, and prostate tumors stored for 3-31 years. Analysis of 10 patients with matched DCIS and cancers that recurred 2-16 years later show that many primary DCIS had already undergone whole-genome doubling and clonal diversification and that they shared genomic lineages with persistent subclones in the recurrences. Evolutionary analysis suggests that most DCIS cases in our cohort underwent an evolutionary bottleneck, and further identified chromosome aberrations in the persistent subclones that were associated with recurrence.

Chemoresistance Evolution in Triple-Negative Breast Cancer Delineated by Single-Cell Sequencing.

Triple-negative breast cancer (TNBC) is an aggressive subtype that frequently develops resistance to chemotherapy. An unresolved question is whether resistance is caused by the selection of rare pre-existing clones or alternatively through the acquisition of new genomic aberrations. To investigate this question, we applied single-cell DNA and RNA sequencing in addition to bulk exome sequencing to profile longitudinal samples from 20 TNBC patients during neoadjuvant chemotherapy (NAC). Deep-exome sequencing identified 10 patients in which NAC led to clonal extinction and 10 patients in which clones persisted after treatment. In 8 patients, we performed a more detailed study using single-cell DNA sequencing to analyze 900 cells and single-cell RNA sequencing to analyze 6,862 cells. Our data showed that resistant genotypes were pre-existing and adaptively selected by NAC, while transcriptional profiles were acquired by reprogramming in response to chemotherapy in TNBC patients.

Multiclonal Invasion in Breast Tumors Identified by Topographic Single Cell Sequencing.

Ductal carcinoma in situ (DCIS) is an early-stage breast cancer that infrequently progresses to invasive ductal carcinoma (IDC). Genomic evolution has been difficult to delineate during invasion due to intratumor heterogeneity and the low number of tumor cells in the ducts. To overcome these challenges, we developed Topographic Single Cell Sequencing (TSCS) to measure genomic copy number profiles of single tumor cells while preserving their spatial context in tissue sections. We applied TSCS to 1,293 single cells from 10 synchronous patients with both DCIS and IDC regions in addition to exome sequencing. Our data reveal a direct genomic lineage between in situ and invasive tumor subpopulations and further show that most mutations and copy number aberrations evolved within the ducts prior to invasion. These results support a multiclonal invasion model, in which one or more clones escape the ducts and migrate into the adjacent tissues to establish the invasive carcinomas.

A spatially resolved single-cell genomic atlas of the adult human breast.

The adult human breast is comprised of an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue1-3. Although most previous studies have focused on the breast epithelial system4-6, many of the non-epithelial cell types remain understudied. Here we constructed the comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics study profiled 714,331 cells from 126 women, and 117,346 nuclei from 20 women, identifying 12 major cell types and 58 biological cell states. These data reveal abundant perivascular, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Spatial mapping using four different technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide a reference of the adult normal breast tissue for studying mammary biology and diseases such as breast cancer.

Simple oligonucleotide-based multiplexing of single-cell chromatin accessibility.

Microdroplet single-cell ATAC-seq is widely used to measure chromatin accessibility, however, highly scalable and simple sample multiplexing procedures are not available. Here, we present a transposome-assisted single nucleus barcoding approach for ATAC-seq (SNuBar-ATAC) that utilizes a single oligonucleotide adaptor for multiplexing samples during the existing tagmentation step and does not require a pre-labeling procedure. The accuracy and scalability of SNuBar-ATAC was evaluated using cell line mixture experiments. We applied SNuBar-ATAC to investigate treatment-induced chromatin accessibility dynamics by multiplexing 28 mice with lung tumors that received different combinations of chemo, radiation, and targeted immunotherapy. We also applied SNuBar-ATAC to study spatial epigenetic heterogeneity by multiplexing 32 regions from a human breast tissue. Additionally, we show that SNuBar can multiplex single cell ATAC and RNA multiomic assays in cell lines and human breast tissue samples. Our data show that SNuBar is a highly accurate, easy-to-use, and scalable system for multiplexing scATAC-seq and scATAC and RNA co-assay experiments.

Breast tumours maintain a reservoir of subclonal diversity during expansion.

Our knowledge of copy number evolution during the expansion of primary breast tumours is limited1,2. Here, to investigate this process, we developed a single-cell, single-molecule DNA-sequencing method and performed copy number analysis of 16,178 single cells from 8 human triple-negative breast cancers and 4 cell lines. The results show that breast tumours and cell lines comprise a large milieu of subclones (7-22) that are organized into a few (3-5) major superclones. Evolutionary analysis suggests that after clonal TP53 mutations, multiple loss-of-heterozygosity events and genome doubling, there was a period of transient genomic instability followed by ongoing copy number evolution during the primary tumour expansion. By subcloning single daughter cells in culture, we show that tumour cells rediversify their genomes and do not retain isogenic properties. These data show that triple-negative breast cancers continue to evolve chromosome aberrations and maintain a reservoir of subclonal diversity during primary tumour growth.

Single-cell DNA sequencing reveals a late-dissemination model in metastatic colorectal cancer.

Metastasis is a complex biological process that has been difficult to delineate in human colorectal cancer (CRC) patients. A major obstacle in understanding metastatic lineages is the extensive intra-tumor heterogeneity at the primary and metastatic tumor sites. To address this problem, we developed a highly multiplexed single-cell DNA sequencing approach to trace the metastatic lineages of two CRC patients with matched liver metastases. Single-cell copy number or mutational profiling was performed, in addition to bulk exome and targeted deep-sequencing. In the first patient, we observed monoclonal seeding, in which a single clone evolved a large number of mutations prior to migrating to the liver to establish the metastatic tumor. In the second patient, we observed polyclonal seeding, in which two independent clones seeded the metastatic liver tumor after having diverged at different time points from the primary tumor lineage. The single-cell data also revealed an unexpected independent tumor lineage that did not metastasize, and early progenitor clones with the "first hit" mutation in APC that subsequently gave rise to both the primary and metastatic tumors. Collectively, these data reveal a late-dissemination model of metastasis in two CRC patients and provide an unprecedented view of metastasis at single-cell genomic resolution.

Cabazitaxel plus carboplatin for the treatment of men with metastatic castration-resistant prostate cancers: a randomised, open-label, phase 1-2 trial.

BACKGROUND: Taxane-platinum combinations have shown promising activity in metastatic castration-resistant prostate cancers in single-group clinical studies but not in randomised trials. Distinct biological subsets of the disease might derive the greatest benefit from the addition of platinum. We aimed to determine whether adding carboplatin to cabazitaxel would improve the outcomes of men with metastatic castration-resistant prostate cancer. METHODS: We did a phase 1-2, open label, randomised study at two centres in men with progressive metastatic castration-resistant prostate cancer. In phase 1, patients received intravenous cabazitaxel 20-25 mg/m2 and intravenous carboplatin area under the curve (AUC) 3-4 mg/mL per min every 21 days. The maximum tolerated dose was defined as the highest dose cohort studied in which one of six or fewer patients experienced a dose-limiting toxicity. In phase 2, patients were randomly assigned (1:1) centrally by a computerised algorithm to intravenous cabazitaxel 25 mg/m2 with or without intravenous carboplatin AUC 4 mg/mL per min. All patients received growth factor support and oral prednisone 10 mg daily. The primary endpoints were the maximum tolerated dose of the combination in phase 1 and investigator-assessed progression-free survival in phase 2. This trial is registered at ClinicalTrials.gov, number NCT01505868. FINDINGS: Between Aug 17, 2012, and May 11, 2015, nine patients completed phase 1 as planned, and 160 were randomly assigned to cabazitaxel (n=79) or cabazitaxel plus carboplatin (n=81) in phase 2. During phase I, grade 3 adverse events were anaemia (n=2), fatigue (n=1), thrombocytopenia (n=1), hypomagnesaemia (n=1), diarrhoea (n=1), hypokalaemia (n=1), anorexia (n=1), and dehydration (n=1), and no grade 4 adverse events occurred. No dose-limiting toxicities were observed, therefore, a maximum tolerated dose of cabazitaxel of 25 mg/m2 and carboplatin of AUC 4 mg/mL per min was selected for phase 2. At a median follow-up of 31·0 months (IQR 20·5-37·1), the combination improved the median progression-free survival from 4·5 months (95% CI 3·5-5·7) to 7·3 months (95% CI 5·5-8·2; hazard ratio 0·69, 95% CI 0·50-0·95, p=0·018). In the phase 2 study, the most common grade 3-5 adverse events were fatigue (7 [9%] of 79 in the cabazitaxel group vs 16 [20%] of 81 in the combination group), anaemia (3 [4%] vs 19 [23%]), neutropenia (3 [4%] vs 13 [16%]), and thrombocytopenia (1 [1%] vs 11 [14%]). There were no treatment-related deaths. INTERPRETATION: Carboplatin added to cabazitaxel showed improved clinical efficacy compared with cabazitaxel alone for men with metastatic castration-resistant prostate cancer. Although adverse events were more common with the combination, the treatment was safe and generally well tolerated. Our data suggest that taxane-platinum combinations have a clinically beneficial role in advanced prostate cancer and a randomised phase 3 study is planned. FUNDING: Sanofi Genzyme, University of Texas MD Anderson Cancer Center Prostate Cancer Moon Shot Program, and Solon Scott III Prostate Cancer Research Fund.

HITS-CLIP analysis uncovers a link between the Kaposi's sarcoma-associated herpesvirus ORF57 protein and host pre-mRNA metabolism.

The Kaposi's sarcoma associated herpesvirus (KSHV) is an oncogenic virus that causes Kaposi's sarcoma, primary effusion lymphoma (PEL), and some forms of multicentric Castleman's disease. The KSHV ORF57 protein is a conserved posttranscriptional regulator of gene expression that is essential for virus replication. ORF57 is multifunctional, but most of its activities are directly linked to its ability to bind RNA. We globally identified virus and host RNAs bound by ORF57 during lytic reactivation in PEL cells using high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP). As expected, ORF57-bound RNA fragments mapped throughout the KSHV genome, including the known ORF57 ligand PAN RNA. In agreement with previously published ChIP results, we observed that ORF57 bound RNAs near the oriLyt regions of the genome. Examination of the host RNA fragments revealed that a subset of the ORF57-bound RNAs was derived from transcript 5' ends. The position of these 5'-bound fragments correlated closely with the 5'-most exon-intron junction of the pre-mRNA. We selected four candidates (BTG1, EGR1, ZFP36, and TNFSF9) and analyzed their pre-mRNA and mRNA levels during lytic phase. Analysis of both steady-state and newly made RNAs revealed that these candidate ORF57-bound pre-mRNAs persisted for longer periods of time throughout infection than control RNAs, consistent with a role for ORF57 in pre-mRNA metabolism. In addition, exogenous expression of ORF57 was sufficient to increase the pre-mRNA levels and, in one case, the mRNA levels of the putative ORF57 targets. These results demonstrate that ORF57 interacts with specific host pre-mRNAs during lytic reactivation and alters their processing, likely by stabilizing pre-mRNAs. These data suggest that ORF57 is involved in modulating host gene expression in addition to KSHV gene expression during lytic reactivation.

Chromatin immunoprecipitation and microarray analysis suggest functional cooperation between Kaposi's Sarcoma-associated herpesvirus ORF57 and K-bZIP.

The Kaposi's sarcoma-associated herpesvirus (KSHV) open reading frame 57 (ORF57)-encoded protein (Mta) is a multifunctional regulator of viral gene expression. ORF57 is essential for viral replication, so elucidation of its molecular mechanisms is important for understanding KSHV infection. ORF57 has been implicated in nearly every aspect of viral gene expression, including transcription, RNA stability, splicing, export, and translation. Here we demonstrate that ORF57 interacts with the KSHV K-bZIP protein in vitro and in cell extracts from lytically reactivated infected cells. To further test the biological relevance of the interaction, we performed a chromatin immunoprecipitation and microarray (ChIP-chip) analysis using anti-ORF57 antibodies and a KSHV tiling array. The results revealed four specific areas of enrichment, including the ORF4 and K8 (K-bZIP) promoters, as well as oriLyt, all of which interact with K-bZIP. In addition, ORF57 associated with DNA corresponding to the PAN RNA transcribed region, a known posttranscriptional target of ORF57. All of the peaks were RNase insensitive, demonstrating that ORF57 association with the viral genome is unlikely to be mediated exclusively by an RNA tether. Our data demonstrate that ORF57 associates with the viral genome by using at least two modes of recruitment, and they suggest that ORF57 and K-bZIP coregulate viral gene expression during lytic infection.

Reconstructing mutational lineages in breast cancer by multi-patient-targeted single-cell DNA sequencing.

Single-cell DNA sequencing (scDNA-seq) methods are powerful tools for profiling mutations in cancer cells; however, most genomic regions sequenced in single cells are non-informative. To overcome this issue, we developed a multi-patient-targeted (MPT) scDNA-seq method. MPT involves first performing bulk exome sequencing across a cohort of cancer patients to identify somatic mutations, which are then pooled together to develop a single custom targeted panel for high-throughput scDNA-seq using a microfluidics platform. We applied MPT to profile 330 mutations across 23,500 cells from 5 patients with triple negative-breast cancer (TNBC), which showed that 3 tumors were monoclonal and 2 tumors were polyclonal. From these data, we reconstructed mutational lineages and identified early mutational and copy-number events, including early TP53 mutations that occurred in all five patients. Collectively, our data suggest that MPT can overcome a major technical obstacle for studying tumor evolution using scDNA-seq by profiling information-rich mutation sites.

Targeting T cell checkpoints 41BB and LAG3 and myeloid cell CXCR1/CXCR2 results in antitumor immunity and durable response in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is considered non-immunogenic, with trials showing its recalcitrance to PD1 and CTLA4 immune checkpoint therapies (ICTs). Here, we sought to systematically characterize the mechanisms underlying de novo ICT resistance and to identify effective therapeutic options for PDAC. We report that agonist 41BB and antagonist LAG3 ICT alone and in combination, increased survival and antitumor immunity, characterized by modulating T cell subsets with antitumor activity, increased T cell clonality and diversification, decreased immunosuppressive myeloid cells and increased antigen presentation/decreased immunosuppressive capability of myeloid cells. Translational analyses confirmed the expression of 41BB and LAG3 in human PDAC. Since single and dual ICTs were not curative, T cell-activating ICTs were combined with a CXCR1/2 inhibitor targeting immunosuppressive myeloid cells. Triple therapy resulted in durable complete responses. Given similar profiles in human PDAC and the availability of these agents for clinical testing, our findings provide a testable hypothesis for this lethal disease.

A spatially resolved single cell genomic atlas of the adult human breast.

The adult human breast comprises an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue. While previous studies have mainly focused on the breast epithelial system, many of the non-epithelial cell types remain understudied. Here, we constructed a comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics data profiled 535,941 cells from 62 women, and 120,024 nuclei from 20 women, identifying 11 major cell types and 53 cell states. These data revealed abundant pericyte, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Our spatial mapping using three technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells in the ducts and lobules, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide an unprecedented reference of adult normal breast tissue for studying mammary biology and disease states such as breast cancer.

Spatial charting of single-cell transcriptomes in tissues.

Single-cell RNA sequencing methods can profile the transcriptomes of single cells but cannot preserve spatial information. Conversely, spatial transcriptomics assays can profile spatial regions in tissue sections, but do not have single-cell resolution. Here, we developed a computational method called CellTrek that combines these two datasets to achieve single-cell spatial mapping through coembedding and metric learning approaches. We benchmarked CellTrek using simulation and in situ hybridization datasets, which demonstrated its accuracy and robustness. We then applied CellTrek to existing mouse brain and kidney datasets and showed that CellTrek can detect topological patterns of different cell types and cell states. We performed single-cell RNA sequencing and spatial transcriptomics experiments on two ductal carcinoma in situ tissues and applied CellTrek to identify tumor subclones that were restricted to different ducts, and specific T cell states adjacent to the tumor areas. Our data show that CellTrek can accurately map single cells in diverse tissue types to resolve their spatial organization.

Spatial intra-tumor heterogeneity is associated with survival of lung adenocarcinoma patients.

Intra-tumor heterogeneity (ITH) of human tumors is important for tumor progression, treatment response, and drug resistance. However, the spatial distribution of ITH remains incompletely understood. Here, we present spatial analysis of ITH in lung adenocarcinomas from 147 patients using multi-region mass spectrometry of >5,000 regions, single-cell copy number sequencing of ~2,000 single cells, and cyclic immunofluorescence of >10 million cells. We identified two distinct spatial patterns among tumors, termed clustered and random geographic diversification (GD). These patterns were observed in the same samples using both proteomic and genomic data. The random proteomic GD pattern, which is characterized by decreased cell adhesion and lower levels of tumor-interacting endothelial cells, was significantly associated with increased risk of recurrence or death in two independent patient cohorts. Our study presents comprehensive spatial mapping of ITH in lung adenocarcinoma and provides insights into the mechanisms and clinical consequences of GD.

Single-Cell Sequencing Reveals Trajectory of Tumor-Infiltrating Lymphocyte States in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) has few effective treatments. Immunotherapy, an attractive alternative strategy, remains challenging with the lack of knowledge on the tumor-infiltrating lymphocyte (TIL) landscape in PDAC. To generate a reference of T-cell subpopulations, we profiled 80,000 T cells from 57 PDAC samples, 22 uninvolved/normal samples, and cultured TIL using single-cell transcriptomic and T-cell receptor analysis. These data revealed 20 cell states and heterogeneous distributions of TIL populations. The CD8+ TIL contained a putative transitional GZMK+ population based on T-cell receptor clonotype sharing, and cell-state trajectory analysis showed similarity to a GZMB+PRF1+ cytotoxic and a CXCL13+ dysfunctional population. Statistical analysis suggested that certain TIL states, such as dysfunctional and inhibitory populations, often occurred together. Finally, analysis of cultured TIL revealed that high-frequency clones from effector populations were preferentially expanded. These data provide a framework for understanding the PDAC TIL landscape for future TIL use in immunotherapy for PDAC. SIGNIFICANCE: To improve the efficacy of immunotherapy in PDAC, there is a great need to understand the PDAC TIL landscape. This study represents a reference of PDAC TIL subpopulations and their relationships and provides a foundation upon which to base future immunotherapeutic efforts. This article is highlighted in the In This Issue feature, p. 2221.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer.

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

Delineating copy number and clonal substructure in human tumors from single-cell transcriptomes.

Single-cell transcriptomic analysis is widely used to study human tumors. However, it remains challenging to distinguish normal cell types in the tumor microenvironment from malignant cells and to resolve clonal substructure within the tumor. To address these challenges, we developed an integrative Bayesian segmentation approach called copy number karyotyping of aneuploid tumors (CopyKAT) to estimate genomic copy number profiles at an average genomic resolution of 5 Mb from read depth in high-throughput single-cell RNA sequencing (scRNA-seq) data. We applied CopyKAT to analyze 46,501 single cells from 21 tumors, including triple-negative breast cancer, pancreatic ductal adenocarcinoma, anaplastic thyroid cancer, invasive ductal carcinoma and glioblastoma, to accurately (98%) distinguish cancer cells from normal cell types. In three breast tumors, CopyKAT resolved clonal subpopulations that differed in the expression of cancer genes, such as KRAS, and signatures, including epithelial-to-mesenchymal transition, DNA repair, apoptosis and hypoxia. These data show that CopyKAT can aid in the analysis of scRNA-seq data in a variety of solid human tumors.

Decoding the evolutionary response to prostate cancer therapy by plasma genome sequencing.

BACKGROUND: Investigating genome evolution in response to therapy is difficult in human tissue samples. To address this challenge, we develop an unbiased whole-genome plasma DNA sequencing approach that concurrently measures genomic copy number and exome mutations from archival cryostored plasma samples. This approach is applied to study longitudinal blood plasma samples from prostate cancer patients, where longitudinal tissue biopsies from the bone and other metastatic sites have been challenging to collect. RESULTS: A molecular characterization of archival plasma DNA from 233 patients and genomic profiling of 101 patients identifies clinical correlations of aneuploid plasma DNA profiles with poor survival, increased plasma DNA concentrations, and lower plasma DNA size distributions. Deep-exome sequencing and genomic copy number profiling are performed on 23 patients, including 9 patients with matched metastatic tissues and 12 patients with serial plasma samples. These data show a high concordance in genomic alterations between the plasma DNA and metastatic tissue samples, suggesting the plasma DNA is highly representative of the tissue alterations. Longitudinal sequencing of 12 patients with 2-5 serial plasma samples reveals clonal dynamics and genome evolution in response to hormonal and chemotherapy. By performing an integrated evolutionary analysis, minor subclones are identified in 9 patients that expanded in response to therapy and harbored mutations associated with resistance. CONCLUSIONS: This study provides an unbiased evolutionary approach to non-invasively delineate clonal dynamics and identify clones with mutations associated with resistance in prostate cancer.

Single-Cell Circulating Tumor Cell Analysis Reveals Genomic Instability as a Distinctive Feature of Aggressive Prostate Cancer.

PURPOSE: Aggressive variant prostate cancer (AVPC) represents a clinical subset distinguished by therapy resistance and poor prognosis, linked to combined losses of the tumor suppressor genes (TSG) PTEN, RB1, and TP53. Circulating tumor cells (CTC) provide a minimally invasive opportunity for identification and molecular characterization of AVPC. We aimed to evaluate the incidence and clinical significance of compound (2+)TSG losses and genomic instability in prostate cancer CTC, and to expand the set genomic biomarkers relevant to AVPC. EXPERIMENTAL DESIGN: Genomic analysis of chromosomal copy-number alterations (CNA) at single-cell resolution was performed in CTC from patients with and without AVPC before initiating chemotherapy with cabazitaxel or cabazitaxel and carboplatin. We evaluated associations between single-CTC genomics and clinical features, progression-free survival, and overall survival. RESULTS: A total of 257 individual CTC were sequenced from 47 patients (1-22 CTC/patient). Twenty patients (42.6%) had concurrent 2+TSG losses in at least one CTC in association with poor survival and increased genomic instability, inferred by high large-scale transitions scores. Higher LST in CTC were independent of CTC enumerated, clinically more indicative of aggressive behavior than co-occurring TSG losses, and molecularly associated with gains in chromosomal regions including PTK2, Myc, and NCOA2; increased androgen receptor expression; and BRCA2 loss. In 57 patients with matched cell-free tumor DNA data, CTC were more frequently detectable and evaluable for CNA analysis (in 73.7% vs. 42.1%, respectively). CONCLUSIONS: Our findings suggest that genomic instability in CTC is a hallmark of advanced prostate cancer aggressiveness, and support single-CTC sequencing as a compelling tool to noninvasively characterize cancer heterogeneity.