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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Nanogrid single-nucleus RNA sequencing reveals phenotypic diversity in breast cancer.

Single cell RNA sequencing has emerged as a powerful tool for resolving transcriptional diversity in tumors, but is limited by throughput, cost and the ability to process archival frozen tissue samples. Here we develop a high-throughput 3' single-nucleus RNA sequencing approach that combines nanogrid technology, automated imaging, and cell selection to sequence up to ~1800 single nuclei in parallel. We compare the transcriptomes of 485 single nuclei to 424 single cells in a breast cancer cell line, which shows a high concordance (93.34%) in gene levels and abundance. We also analyze 416 nuclei from a frozen breast tumor sample and 380 nuclei from normal breast tissue. These data reveal heterogeneity in cancer cell phenotypes, including angiogenesis, proliferation, and stemness, and a minor subpopulation (19%) with many overexpressed cancer genes. Our studies demonstrate the utility of nanogrid single-nucleus RNA sequencing for studying the transcriptional programs of tumor nuclei in frozen archival tissue samples.Single cell RNA sequencing is a powerful tool for understanding cellular diversity but is limited by cost, throughput and sample preparation. Here the authors use nanogrid technology with integrated imaging to sequence thousands of cancer nuclei in parallel from fresh or frozen tissue.

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.

Punctuated copy number evolution and clonal stasis in triple-negative breast cancer.

Aneuploidy is a hallmark of breast cancer; however, knowledge of how these complex genomic rearrangements evolve during tumorigenesis is limited. In this study, we developed a highly multiplexed single-nucleus sequencing method to investigate copy number evolution in patients with triple-negative breast cancer. We sequenced 1,000 single cells from tumors in 12 patients and identified 1-3 major clonal subpopulations in each tumor that shared a common evolutionary lineage. For each tumor, we also identified a minor subpopulation of non-clonal cells that were classified as metastable, pseudodiploid or chromazemic. Phylogenetic analysis and mathematical modeling suggest that these data are unlikely to be explained by the gradual accumulation of copy number events over time. In contrast, our data challenge the paradigm of gradual evolution, showing that the majority of copy number aberrations are acquired at the earliest stages of tumor evolution, in short punctuated bursts, followed by stable clonal expansions that form the tumor mass.