SUMMARY Sequencing studies of breast tumor cohorts have identified many prevalent mutations, but provide limited insight into the genomic diversity within tumors. Here, we developed a whole-genome and exome single cell sequencing approach called Nuc-Seq that utilizes G2/M nuclei to achieve 91% mean coverage breadth. We applied this method to sequence single normal and tumor nuclei from an estrogen-receptor positive breast cancer and a triple-negative ductal carcinoma. In parallel, we performed single nuclei copy number profiling. Our data show that aneuploid rearrangements occurred early in tumor evolution and remained highly stable as the tumor masses clonally expanded. In contrast, point mutations evolved gradually, generating extensive clonal diversity. Many of the diverse mutations were shown to occur at low frequencies (<10%) in the tumor mass by targeted single-molecule sequencing. Using mathematical modeling we found that the triple-negative tumor cells had an increased mutation rate (13.3X) while the ER+ tumor cells did not. These findings have important implications for the diagnosis, therapeutic treatment and evolution of chemoresistance in breast cancer.

Cancer progression in humans is difficult to infer because we do not routinely sample patients at multiple stages of their disease. However, heterogeneous breast tumors provide a unique opportunity to study human tumor progression because they still contain evidence of early and intermediate subpopulations in the form of the phylogenetic relationships. We have developed a method we call Sector-Ploidy-Profiling (SPP) to study the clonal composition of breast tumors. SPP involves macro-dissecting tumors, flow-sorting genomic subpopulations by DNA content, and profiling genomes using comparative genomic hybridization (CGH). Breast carcinomas display two classes of genomic structural variation: (1) monogenomic and (2) polygenomic. Monogenomic tumors appear to contain a single major clonal subpopulation with a highly stable chromosome structure. Polygenomic tumors contain multiple clonal tumor subpopulations, which may occupy the same sectors, or separate anatomic locations. In polygenomic tumors, we show that heterogeneity can be ascribed to a few clonal subpopulations, rather than a series of gradual intermediates. By comparing multiple subpopulations from different anatomic locations, we have inferred pathways of cancer progression and the organization of tumor growth.

The deadliest anaplastic thyroid cancer (ATC) often transforms from indolent differentiated thyroid cancer (DTC); however, the complex intra-tumor transformation process is poorly understood.We investigated an anaplastic transformation model by dissecting both cell lineage and cell fate transitions using single cell transcriptomes and genetic alteration data from patients with different subtypes of thyroid cancer. The resulting spectrum of ATC transformation included stressresponsive DTC cells, inflammatory ATC cells (iATCs), mitotic-defective ATC cells and extended all the way to mesenchymal ATC cells (mATCs). Further, our analysis identified two important milestones: 1) a diploid stage, where iATC cells were diploids with inflammatory phenotypes, and 2) an aneuploid stage, where mATCs gained aneuploid genomes and mesenchymal phenotypes producing excessive collagens and collagen-interacting receptors. In parallel, cancer-associatedfibroblasts showed strong interactions among mesenchymal cell-types, macrophages shifted from M1 to M2 states, and T cells reprogrammed from cytotoxic to exhausted states, highlighting new therapeutic opportunities for ATC.

http://prober.cshl.edu

Pure ductal carcinoma in situ (DCIS) is being diagnosed more frequently through breast screening programmes and is associated with an increased risk of developing invasive breast cancer. We assessed the clonal relatedness of 143 cases of pure DCIS and their subsequent events using a combination of whole exome, targeted and copy number sequencing, supplemented by single cell analysis. Unexpectedly, 18% of all invasive events after DCIS were clonally unrelated to the primary DCIS. Single cell sequencing of selected pairs confirmed our findings. In contrast, synchronous DCIS and invasive disease (n=44) were almost always (93%) clonally related. This challenges the dogma that most invasive events after DCIS represent invasive transformation of the initial DCIS and suggests that DCIS could be an independent risk factor for developing invasive disease as well as a precursor lesion.