PCA-Plus: Enhanced principal component analysis with illustrative applications to batch effects and their quantitation.

Background: Principal component analysis (PCA), a standard approach to analysis and visualization of large datasets, is commonly used in biomedical research for detecting similarities and differences among groups of samples. We initially used conventional PCA as a tool for critical quality control of batch and trend effects in multi-omic profiling data produced by The Cancer Genome Atlas (TCGA) project of the NCI. We found, however, that conventional PCA visualizations were often hard to interpret when inter-batch differences were moderate in comparison with intra-batch differences; it was also difficult to quantify batch effects objectively. We, therefore, sought enhancements to make the method more informative in those and analogous settings. Results: We have developed algorithms and a toolbox of enhancements to conventional PCA that improve the detection, diagnosis, and quantitation of differences between or among groups, e.g., groups of molecularly profiled biological samples. The enhancements include (i) computed group centroids; (ii) sample-dispersion rays; (iii) differential coloring of centroids, rays, and sample data points; (iii) trend trajectories; and (iv) a novel separation index (DSC) for quantitation of differences among groups. Conclusions: PCA-Plus has been our most useful single tool for analyzing, visualizing, and quantitating batch effects, trend effects, and class differences in molecular profiling data of many types: mRNA expression, microRNA expression, DNA methylation, and DNA copy number. An early version of PCA-Plus has been used as the central graphical visualization in our MBatch package for near-real-time surveillance of data for analysis working groups in more than 70 TCGA, PanCancer Atlas, PanCancer Analysis of Whole Genomes, and Genome Data Analysis Network projects of the NCI. The algorithms and software are generic, hence applicable more generally to other types of multivariate data as well. PCA-Plus is freely available in a down-loadable R package at our MBatch website.

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.

Learning to distinguish progressive and non-progressive ductal carcinoma in situ.

Ductal carcinoma in situ (DCIS) is a non-invasive breast neoplasia that accounts for 25% of all screen-detected breast cancers diagnosed annually. Neoplastic cells in DCIS are confined to the ductal system of the breast, although they can escape and progress to invasive breast cancer in a subset of patients. A key concern of DCIS is overtreatment, as most patients screened for DCIS and in whom DCIS is diagnosed will not go on to exhibit symptoms or die of breast cancer, even if left untreated. However, differentiating low-risk, indolent DCIS from potentially progressive DCIS remains challenging. In this Review, we summarize our current knowledge of DCIS and explore open questions about the basic biology of DCIS, including those regarding how genomic events in neoplastic cells and the surrounding microenvironment contribute to the progression of DCIS to invasive breast cancer. Further, we discuss what information will be needed to prevent overtreatment of indolent DCIS lesions without compromising adequate treatment for high-risk patients.

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.

Retinal patterns and the cellular repertoire of neuropsin (Opn5) retinal ganglion cells.

Obtaining a parts list of the sensory components of the retina is vital to understand the effects of light in behavior, health, and disease. Rods, cones, and intrinsically photosensitive retinal ganglion cells (ipRGCs) are the best described photoreceptors in the mammalian retina, but recent functional roles have been proposed for retinal neuropsin (Opn5)-an atypical opsin. However, little is known about the pattern of Opn5 expression in the retina. Using cre (Opn5cre ) and cre-dependent reporters, we uncover patterns of Opn5 expression and find that Opn5 is restricted to retinal ganglion cells (RGCs). Opn5-RGCs are nonhomogenously distributed through the retina, with greater densities of cells located in the dorsotemporal quadrant. In addition to the local topology of these cells, using cre-dependent AAV viral tracing, we surveyed their central targets and found that they are biased towards image-forming and image-stabilizing regions. Finally, molecular and electrophysiological profiling reveal that Opn5-RGCs comprise previously defined RGC types that respond optimally to edges and object-motion (F-mini-ONs, HD2, HD1, LEDs, ooDSRGCs, etc.). Together, these data describe the second collection of RGCs that express atypical opsins in the mouse, and expand the roles of image-forming cells in retinal physiology and function.

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.

LKB1 coordinates neurite remodeling to drive synapse layer emergence in the outer retina.

Structural changes in pre and postsynaptic neurons that accompany synapse formation often temporally and spatially overlap. Thus, it has been difficult to resolve which processes drive patterned connectivity. To overcome this, we use the laminated outer murine retina. We identify the serine/threonine kinase LKB1 as a key driver of synapse layer emergence. The absence of LKB1 in the retina caused a marked mislocalization and delay in synapse layer formation. In parallel, LKB1 modulated postsynaptic horizontal cell refinement and presynaptic photoreceptor axon growth. Mislocalized horizontal cell processes contacted aberrant cone axons in LKB1 mutants. These defects coincided with altered synapse protein organization, and horizontal cell neurites were misdirected to ectopic synapse protein regions. Together, these data suggest that LKB1 instructs the timing and location of connectivity in the outer retina via coordinate regulation of pre and postsynaptic neuron structure and the localization of synapse-associated proteins.

Spatiotemporal gene expression patterns reveal molecular relatedness between retinal laminae.

In several areas of the central nervous system, neurons are regionally organized into groups or layers that carry out specific activities. In this form of patterning, neurons of distinct types localize their cell bodies to just one or a few of the layers within a structure. However, little is known about whether diverse neuron types within a lamina share molecular features that coordinate their organization. To begin to identify such candidates, we used the laminated murine retina to screen 92 lacZ reporter lines available through the Knockout Mouse Project. Thirty-two of these displayed reporter expression in restricted subsets of inner retina neurons. We then identified the spatiotemporal expression patterns of these genes at key developmental stages. This uncovered several that were heavily enriched in development but reduced in adulthood, including the transcriptional regulator Hmga1. An additional set of genes displayed maturation associated laminar enrichment. Among these, we identified Bbox1 as a novel gene that specifically labels all neurons in the ganglion cell layer but is largely excluded from otherwise molecularly similar neurons in the inner retina. Finally, we established Dbn1 as a new marker enriched in amacrines and Fmnl3 as a marker for subsets of αRGCs. Together, these data provide a spatiotemporal map for laminae-specific molecules and suggest that diverse neuron types within a lamina share coordinating molecular features that may inform their fate or function.

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.

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.

Genome evolution in ductal carcinoma in situ: invasion of the clones.

Ductal carcinoma in situ (DCIS) is the most frequently diagnosed early-stage breast cancer. Only a subset of patients progress to invasive ductal carcinoma (IDC), and this presents a formidable clinical challenge for determining which patients to treat aggressively and which patients to monitor without therapeutic intervention. Understanding the molecular and genomic basis of invasion has been difficult to study in DCIS cancers due to several technical obstacles, including low tumour cellularity, lack of fresh-frozen tissues, and intratumour heterogeneity. In this review, we discuss the role of intratumour heterogeneity in the progression of DCIS to IDC in the context of three evolutionary models: independent lineages, evolutionary bottlenecks, and multiclonal invasion. We examine the evidence in support of these models and their relevance to the diagnosis and treatment of patients with DCIS. We also discuss how emerging technologies, such as single-cell sequencing, STAR-FISH, and imaging mass spectrometry, are likely to provide new insights into the evolution of this enigmatic disease. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

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.