Magnetically functionalized hydrogels for high-throughput genomic applications.

Single-cell genomics has revolutionized tissue analysis by revealing the genetic program of individual cells. The key aspect of the technology is the use of barcoded beads to unambiguously tag sequences originating from a single cell. The generation of unique barcodes on beads is mainly achieved by split-pooling methods, which are labor-intensive due to repeated washing steps. Towards the automation of the split-pooling method, we developed a simple method to magnetize hydrogel beads. We show that these hydrogel beads provide increased yields and washing efficiencies for purification procedures. They are also fully compatible with single-cell sequencing using the BAG-Seq workflow. Our work opens the automation of the split-pooling technique, which will improve single-cell genomic workflows.

FABP5 Inhibition against PTEN-Mutant Therapy Resistant Prostate Cancer.

Resistance to standard of care taxane and androgen deprivation therapy (ADT) causes the vast majority of prostate cancer (PC) deaths worldwide. We have developed RapidCaP, an autochthonous genetically engineered mouse model of PC. It is driven by the loss of PTEN and p53, the most common driver events in PC patients with life-threatening diseases. As in human ADT, surgical castration of RapidCaP animals invariably results in disease relapse and death from the metastatic disease burden. Fatty Acid Binding Proteins (FABPs) are a large family of signaling lipid carriers. They have been suggested as drivers of multiple cancer types. Here we combine analysis of primary cancer cells from RapidCaP (RCaP cells) with large-scale patient datasets to show that among the 10 FABP paralogs, FABP5 is the PC-relevant target. Next, we show that RCaP cells are uniquely insensitive to both ADT and taxane treatment compared to a panel of human PC cell lines. Yet, they share an exquisite sensitivity to the small-molecule FABP5 inhibitor SBFI-103. We show that SBFI-103 is well tolerated and can strongly eliminate RCaP tumor cells in vivo. This provides a pre-clinical platform to fight incurable PC and suggests an important role for FABP5 in PTEN-deficient PC.

Ordered and deterministic cancer genome evolution after p53 loss.

Although p53 inactivation promotes genomic instability1 and presents a route to malignancy for more than half of all human cancers2,3, the patterns through which heterogenous TP53 (encoding human p53) mutant genomes emerge and influence tumorigenesis remain poorly understood. Here, in a mouse model of pancreatic ductal adenocarcinoma that reports sporadic p53 loss of heterozygosity before cancer onset, we find that malignant properties enabled by p53 inactivation are acquired through a predictable pattern of genome evolution. Single-cell sequencing and in situ genotyping of cells from the point of p53 inactivation through progression to frank cancer reveal that this deterministic behaviour involves four sequential phases-Trp53 (encoding mouse p53) loss of heterozygosity, accumulation of deletions, genome doubling, and the emergence of gains and amplifications-each associated with specific histological stages across the premalignant and malignant spectrum. Despite rampant heterogeneity, the deletion events that follow p53 inactivation target functionally relevant pathways that can shape genomic evolution and remain fixed as homogenous events in diverse malignant populations. Thus, loss of p53-the 'guardian of the genome'-is not merely a gateway to genetic chaos but, rather, can enable deterministic patterns of genome evolution that may point to new strategies for the treatment of TP53-mutant tumours.

Parity-induced changes to mammary epithelial cells control NKT cell expansion and mammary oncogenesis.

Pregnancy reprograms mammary epithelial cells (MECs) to control their responses to pregnancy hormone re-exposure and carcinoma progression. However, the influence of pregnancy on the mammary microenvironment is less clear. Here, we used single-cell RNA sequencing to profile the composition of epithelial and non-epithelial cells in mammary tissue from nulliparous and parous female mice. Our analysis indicates an expansion of γδ natural killer T-like immune cells (NKTs) following pregnancy and upregulation of immune signaling molecules in post-pregnancy MECs. We show that expansion of NKTs following pregnancy is due to elevated expression of the antigen-presenting molecule CD1d on MECs. Loss of CD1d expression on post-pregnancy MECs, or overall lack of activated NKTs, results in mammary oncogenesis. Collectively, our findings illustrate how pregnancy-induced changes modulate the communication between MECs and the immune microenvironment and establish a causal link between pregnancy, the immune microenvironment, and mammary oncogenesis.

Chromosomal instability accelerates the evolution of resistance to anti-cancer therapies.

Aneuploidy is a ubiquitous feature of human tumors, but the acquisition of aneuploidy typically antagonizes cellular fitness. To investigate how aneuploidy could contribute to tumor growth, we triggered periods of chromosomal instability (CIN) in human cells and then exposed them to different culture environments. We discovered that transient CIN reproducibly accelerates the acquisition of resistance to anti-cancer therapies. Single-cell sequencing revealed that these resistant populations develop recurrent aneuploidies, and independently deriving one chromosome-loss event that was frequently observed in paclitaxel-resistant cells was sufficient to decrease paclitaxel sensitivity. Finally, we demonstrated that intrinsic levels of CIN correlate with poor responses to numerous therapies in human tumors. Our results show that, although CIN generally decreases cancer cell fitness, it also provides phenotypic plasticity to cancer cells that can allow them to adapt to diverse stressful environments. Moreover, our findings suggest that aneuploidy may function as an under-explored cause of therapy failure.

Intraductal Transplantation Models of Human Pancreatic Ductal Adenocarcinoma Reveal Progressive Transition of Molecular Subtypes.

Pancreatic ductal adenocarcinoma (PDAC) is the most lethal common malignancy, with little improvement in patient outcomes over the past decades. Recently, subtypes of pancreatic cancer with different prognoses have been elaborated; however, the inability to model these subtypes has precluded mechanistic investigation of their origins. Here, we present a xenotransplantation model of PDAC in which neoplasms originate from patient-derived organoids injected directly into murine pancreatic ducts. Our model enables distinction of the two main PDAC subtypes: intraepithelial neoplasms from this model progress in an indolent or invasive manner representing the classical or basal-like subtypes of PDAC, respectively. Parameters that influence PDAC subtype specification in this intraductal model include cell plasticity and hyperactivation of the RAS pathway. Finally, through intratumoral dissection and the direct manipulation of RAS gene dosage, we identify a suite of RAS-regulated secreted and membrane-bound proteins that may represent potential candidates for therapeutic intervention in patients with PDAC. SIGNIFICANCE: Accurate modeling of the molecular subtypes of pancreatic cancer is crucial to facilitate the generation of effective therapies. We report the development of an intraductal organoid transplantation model of pancreatic cancer that models the progressive switching of subtypes, and identify stochastic and RAS-driven mechanisms that determine subtype specification.See related commentary by Pickering and Morton, p. 1448.This article is highlighted in the In This Issue feature, p. 1426.

Novel insights into breast cancer copy number genetic heterogeneity revealed by single-cell genome sequencing.

Copy number alterations (CNAs) play an important role in molding the genomes of breast cancers and have been shown to be clinically useful for prognostic and therapeutic purposes. However, our knowledge of intra-tumoral genetic heterogeneity of this important class of somatic alterations is limited. Here, using single-cell sequencing, we comprehensively map out the facets of copy number alteration heterogeneity in a cohort of breast cancer tumors. Ou/var/www/html/elife/12-05-2020/backup/r analyses reveal: genetic heterogeneity of non-tumor cells (i.e. stroma) within the tumor mass; the extent to which copy number heterogeneity impacts breast cancer genomes and the importance of both the genomic location and dosage of sub-clonal events; the pervasive nature of genetic heterogeneity of chromosomal amplifications; and the association of copy number heterogeneity with clinical and biological parameters such as polyploidy and estrogen receptor negative status. Our data highlight the power of single-cell genomics in dissecting, in its many forms, intra-tumoral genetic heterogeneity of CNAs, the magnitude with which CNA heterogeneity affects the genomes of breast cancers, and the potential importance of CNA heterogeneity in phenomena such as therapeutic resistance and disease relapse.

Cells in the body remain healthy by tightly preventing and repairing random changes, or mutations, in their genetic material. In cancer cells, however, these mechanisms can break down. When these cells grow and multiply, they can then go on to accumulate many mutations. As a result, cancer cells in the same tumor can each contain a unique combination of genetic changes. This genetic heterogeneity has the potential to affect how cancer responds to treatment, and is increasingly becoming appreciated clinically. For example, if a drug only works against cancer cells carrying a specific mutation, any cells lacking this genetic change will keep growing and cause a relapse. However, it is still difficult to quantify and understand genetic heterogeneity in cancer. Copy number alterations (or CNAs) are a class of mutation where large and small sections of genetic material are gained or lost. This can result in cells that have an abnormal number of copies of the genes in these sections. Here, Baslan et al. set out to explore how CNAs might vary between individual cancer cells within the same tumor. To do so, thousands of individual cancer cells were isolated from human breast tumors, and a technique called single-cell genome sequencing used to screen the genetic information of each of them. These experiments confirmed that CNAs did differ ­ sometimes dramatically ­ between patients and among cells taken from the same tumor. For example, many of the cells carried extra copies of well-known cancer genes important for treatment, but the exact number of copies varied between cells. This heterogeneity existed for individual genes as well as larger stretches of DNA: this was the case, for instance, for an entire section of chromosome 8, a region often affected in breast and other tumors. The work by Baslan et al. captures the sheer extent of genetic heterogeneity in cancer and in doing so, highlights the power of single-cell genome sequencing. In the future, a finer understanding of the genetic changes present at the level of an individual cancer cell may help clinicians to manage the disease more effectively.

Integrated Computational Pipeline for Single-Cell Genomic Profiling.

PURPOSE: Copy-number profiling of multiple individual cells from sparse sequencing may be used to reveal a detailed picture of genomic heterogeneity and clonal organization in a tissue biopsy specimen. We sought to provide a comprehensive computational pipeline for single-cell genomics, to facilitate adoption of this molecular technology for basic and translational research. MATERIALS AND METHODS: The pipeline comprises software tools programmed in Python and in R and depends on Bowtie, HISAT2, Matplotlib, and Qt. It is installed and used with Anaconda. RESULTS: Here we describe a complete pipeline for sparse single-cell genomic data, encompassing all steps of single-nucleus DNA copy-number profiling, from raw sequence processing to clonal structure analysis and visualization. For the latter, a specialized graphical user interface termed the single-cell genome viewer (SCGV) is provided. With applications to cancer diagnostics in mind, the SCGV allows for zooming and linkage to the University of California at Santa Cruz Genome Browser from each of the multiple integrated views of single-cell copy-number profiles. The latter can be organized by clonal substructure or by any of the associated metadata such as anatomic location and histologic characterization. CONCLUSION: The pipeline is available as open-source software for Linux and OS X. Its modular structure, extensive documentation, and ease of deployment using Anaconda facilitate its adoption by researchers and practitioners of single-cell genomics. With open-source availability and Massachusetts Institute of Technology licensing, it provides a basis for additional development by the cancer bioinformatics community.

Multiplex accurate sensitive quantitation (MASQ) with application to minimal residual disease in acute myeloid leukemia.

Measuring minimal residual disease in cancer has applications for prognosis, monitoring treatment and detection of recurrence. Simple sequence-based methods to detect nucleotide substitution variants have error rates (about 10-3) that limit sensitive detection. We developed and characterized the performance of MASQ (multiplex accurate sensitive quantitation), a method with an error rate below 10-6. MASQ counts variant templates accurately in the presence of millions of host genomes by using tags to identify each template and demanding consensus over multiple reads. Since the MASQ protocol multiplexes 50 target loci, we can both integrate signal from multiple variants and capture subclonal response to treatment. Compared to existing methods for variant detection, MASQ achieves an excellent combination of sensitivity, specificity and yield. We tested MASQ in a pilot study in acute myeloid leukemia (AML) patients who entered complete remission. We detect leukemic variants in the blood and bone marrow samples of all five patients, after induction therapy, at levels ranging from 10-2 to nearly 10-6. We observe evidence of sub-clonal structure and find higher target variant frequencies in patients who go on to relapse, demonstrating the potential for MASQ to quantify residual disease in AML.

Single-Chromosomal Gains Can Function as Metastasis Suppressors and Promoters in Colon Cancer.

High levels of cancer aneuploidy are frequently associated with poor prognosis. To examine the relationship between aneuploidy and cancer progression, we analyzed a series of congenic cell lines that harbor single extra chromosomes. We found that across 13 different trisomic cell lines, 12 trisomies suppressed invasiveness or were largely neutral, while a single trisomy increased metastatic behavior by triggering a partial epithelial-mesenchymal transition. In contrast, we discovered that chromosomal instability activates cGAS/STING signaling but strongly suppresses invasiveness. By analyzing patient copy-number data, we demonstrate that specific aneuploidies are associated with distinct outcomes, and the acquisition of certain aneuploidies is in fact linked with a favorable prognosis. Thus, aneuploidy is not a uniform driver of malignancy, and different aneuploidies can uniquely influence tumor progression. At the same time, the gain of a single chromosome is capable of inducing a profound cell state transition, thereby linking genomic plasticity, phenotypic plasticity, and metastasis.

Copolymerization of single-cell nucleic acids into balls of acrylamide gel.

We show the use of 5'-Acrydite oligonucleotides to copolymerize single-cell DNA or RNA into balls of acrylamide gel (BAGs). Combining this step with split-and-pool techniques for creating barcodes yields a method with advantages in cost and scalability, depth of coverage, ease of operation, minimal cross-contamination, and efficient use of samples. We perform DNA copy number profiling on mixtures of cell lines, nuclei from frozen prostate tumors, and biopsy washes. As applied to RNA, the method has high capture efficiency of transcripts and sufficient consistency to clearly distinguish the expression patterns of cell lines and individual nuclei from neurons dissected from the mouse brain. By using varietal tags (UMIs) to achieve sequence error correction, we show extremely low levels of cross-contamination by tracking source-specific SNVs. The method is readily modifiable, and we will discuss its adaptability and diverse applications.

Utility of Single-Cell Genomics in Diagnostic Evaluation of Prostate Cancer.

A distinction between indolent and aggressive disease is a major challenge in diagnostics of prostate cancer. As genetic heterogeneity and complexity may influence clinical outcome, we have initiated studies on single tumor cell genomics. In this study, we demonstrate that sparse DNA sequencing of single-cell nuclei from prostate core biopsies is a rich source of quantitative parameters for evaluating neoplastic growth and aggressiveness. These include the presence of clonal populations, the phylogenetic structure of those populations, the degree of the complexity of copy-number changes in those populations, and measures of the proportion of cells with clonal copy-number signatures. The parameters all showed good correlation to the measure of prostatic malignancy, the Gleason score, derived from individual prostate biopsy tissue cores. Remarkably, a more accurate histopathologic measure of malignancy, the surgical Gleason score, agrees better with these genomic parameters of diagnostic biopsy than it does with the diagnostic Gleason score and related measures of diagnostic histopathology. This is highly relevant because primary treatment decisions are dependent upon the biopsy and not the surgical specimen. Thus, single-cell analysis has the potential to augment traditional core histopathology, improving both the objectivity and accuracy of risk assessment and inform treatment decisions.Significance: Genomic analysis of multiple individual cells harvested from prostate biopsies provides an indepth view of cell populations comprising a prostate neoplasm, yielding novel genomic measures with the potential to improve the accuracy of diagnosis and prognosis in prostate cancer. Cancer Res; 78(2); 348-58. ©2017 AACR.

Effect of Blood Collection Tube Type and Time to Processing on the Enumeration and High-Content Characterization of Circulating Tumor Cells Using the High-Definition Single-Cell Assay.

CONTEXT: - As circulating tumor cell (CTC) assays gain clinical relevance, it is essential to address preanalytic variability and to develop standard operating procedures for sample handling in order to successfully implement genomically informed, precision health care. OBJECTIVE: - To evaluate the effects of blood collection tube (BCT) type and time-to-assay (TTA) on the enumeration and high-content characterization of CTCs by using the high-definition single-cell assay (HD-SCA). DESIGN: - Blood samples of patients with early- and advanced-stage breast cancer were collected into cell-free DNA (CfDNA), EDTA, acid-citrate-dextrose solution, and heparin BCTs. Time-to-assay was evaluated at 24 and 72 hours, representing the fastest possible and more routine domestic shipping intervals, respectively. RESULTS: - We detected the highest CTC levels and the lowest levels of negative events in CfDNA BCT at 24 hours. At 72 hours in this BCT, all CTC subpopulations were decreased with the larger effect observed in high-definition CTCs and cytokeratin-positive cells smaller than white blood cells. Overall cell retention was also optimal in CfDNA BCT at 24 hours. Whole-genome copy number variation profiles were generated from single cells isolated from all BCT types and TTAs. Cells from CfDNA BCT at 24-hour TTA exhibited the least noise. CONCLUSIONS: - Circulating tumor cells can be identified and characterized under a variety of collection, handling, and processing conditions, but the highest quality can be achieved with optimized conditions. We quantified performance differences of the HD-SCA for specific preanalytic variables that may be used as a guide to develop best practices for implementation into patient care and/or research biorepository processes.

Early Detection of Cancer in Blood Using Single-Cell Analysis: A Proposal.

Here, we explore the potential of single-cell genomic analysis in blood for early detection of cancer; we consider a method that screens the presence of recurrent patterns of copy number (CN) alterations using sparse single-cell sequencing. We argue for feasibility, based on in silico analysis of existing single-cell data and cancer CN profiles. Sampling procedures from existing diploid single cells can render data for a cell with any given profile. Sampling from multiple published tumor profiles can interrogate cancer clonality via an algorithm that tests the multiplicity of close pairwise similarities among single-cell cancer genomes. The majority of common solid cancers would be detectable in this manner. As any early detection method must be verifiable and actionable, we describe how further analysis of suspect cells can aid in determining risk and anatomic origin. Future affordability rests on currently available procedures for tumor cell enrichment and inexpensive methods for single-cell analysis.

Whole-genome single-cell copy number profiling from formalin-fixed paraffin-embedded samples.

A substantial proportion of tumors consist of genotypically distinct subpopulations of cancer cells. This intratumor genetic heterogeneity poses a substantial challenge for the implementation of precision medicine. Single-cell genomics constitutes a powerful approach to resolve complex mixtures of cancer cells by tracing cell lineages and discovering cryptic genetic variations that would otherwise be obscured in tumor bulk analyses. Because of the chemical alterations that result from formalin fixation, single-cell genomic approaches have largely remained limited to fresh or rapidly frozen specimens. Here we describe the development and validation of a robust and accurate methodology to perform whole-genome copy-number profiling of single nuclei obtained from formalin-fixed paraffin-embedded clinical tumor samples. We applied the single-cell sequencing approach described here to study the progression from in situ to invasive breast cancer, which revealed that ductal carcinomas in situ show intratumor genetic heterogeneity at diagnosis and that these lesions may progress to invasive breast cancer through a variety of evolutionary processes.

TGF-ß reduces DNA ds-break repair mechanisms to heighten genetic diversity and adaptability of CD44+/CD24- cancer cells.

Many lines of evidence have indicated that both genetic and non-genetic determinants can contribute to intra-tumor heterogeneity and influence cancer outcomes. Among the best described sub-population of cancer cells generated by non-genetic mechanisms are cells characterized by a CD44+/CD24- cell surface marker profile. Here, we report that human CD44+/CD24- cancer cells are genetically highly unstable because of intrinsic defects in their DNA-repair capabilities. In fact, in CD44+/CD24- cells, constitutive activation of the TGF-beta axis was both necessary and sufficient to reduce the expression of genes that are crucial in coordinating DNA damage repair mechanisms. Consequently, we observed that cancer cells that reside in a CD44+/CD24- state are characterized by increased accumulation of DNA copy number alterations, greater genetic diversity and improved adaptability to drug treatment. Together, these data suggest that the transition into a CD44+/CD24- cell state can promote intra-tumor genetic heterogeneity, spur tumor evolution and increase tumor fitness.

SMASH, a fragmentation and sequencing method for genomic copy number analysis.

Copy number variants (CNVs) underlie a significant amount of genetic diversity and disease. CNVs can be detected by a number of means, including chromosomal microarray analysis (CMA) and whole-genome sequencing (WGS), but these approaches suffer from either limited resolution (CMA) or are highly expensive for routine screening (both CMA and WGS). As an alternative, we have developed a next-generation sequencing-based method for CNV analysis termed SMASH, for short multiply aggregated sequence homologies. SMASH utilizes random fragmentation of input genomic DNA to create chimeric sequence reads, from which multiple mappable tags can be parsed using maximal almost-unique matches (MAMs). The SMASH tags are then binned and segmented, generating a profile of genomic copy number at the desired resolution. Because fewer reads are necessary relative to WGS to give accurate CNV data, SMASH libraries can be highly multiplexed, allowing large numbers of individuals to be analyzed at low cost. Increased genomic resolution can be achieved by sequencing to higher depth.

Interactive analysis and assessment of single-cell copy-number variations.

We present Ginkgo (http://qb.cshl.edu/ginkgo), a user-friendly, open-source web platform for the analysis of single-cell copy-number variations (CNVs). Ginkgo automatically constructs copy-number profiles of cells from mapped reads and constructs phylogenetic trees of related cells. We validated Ginkgo by reproducing the results of five major studies. After comparing three commonly used single-cell amplification techniques, we concluded that degenerate oligonucleotide-primed PCR is the most consistent for CNV analysis.

Optimizing sparse sequencing of single cells for highly multiplex copy number profiling.

Genome-wide analysis at the level of single cells has recently emerged as a powerful tool to dissect genome heterogeneity in cancer, neurobiology, and development. To be truly transformative, single-cell approaches must affordably accommodate large numbers of single cells. This is feasible in the case of copy number variation (CNV), because CNV determination requires only sparse sequence coverage. We have used a combination of bioinformatic and molecular approaches to optimize single-cell DNA amplification and library preparation for highly multiplexed sequencing, yielding a method that can produce genome-wide CNV profiles of up to a hundred individual cells on a single lane of an Illumina HiSeq instrument. We apply the method to human cancer cell lines and biopsied cancer tissue, thereby illustrating its efficiency, reproducibility, and power to reveal underlying genetic heterogeneity and clonal phylogeny. The capacity of the method to facilitate the rapid profiling of hundreds to thousands of single-cell genomes represents a key step in making single-cell profiling an easily accessible tool for studying cell lineage.