Dissecting the clonal origins of childhood acute lymphoblastic leukemia by single-cell genomics. Proc Natl Acad Sci U S A. 2014;111:17947-17952. [PMID: 25425670 DOI: 10.1073/pnas.1420822111]

Single-Cell Sequencing: Biological Insight and Potential Clinical Implications in Pediatric Leukemia

Single-cell sequencing (SCS) provides high-resolution insight into the genomic, epigenomic, and transcriptomic landscape of oncohematological malignancies including pediatric leukemia, the most common type of childhood cancer. Besides broadening our biological understanding of cellular heterogeneity, sub-clonal architecture, and regulatory network of tumor cell populations, SCS can offer clinically relevant, detailed characterization of distinct compartments affected by leukemia and identify therapeutically exploitable vulnerabilities. In this review, we provide an overview of SCS studies focused on the high-resolution genomic and transcriptomic scrutiny of pediatric leukemia. Our aim is to investigate and summarize how different layers of single-cell omics approaches can expectedly support clinical decision making in the future. Although the clinical management of pediatric leukemia underwent a spectacular improvement during the past decades, resistant disease is a major cause of therapy failure. Currently, only a small proportion of childhood leukemia patients benefit from genomics-driven therapy, as 15-20% of them meet the indication criteria of on-label targeted agents, and their overall response rate falls in a relatively wide range (40-85%). The in-depth scrutiny of various cell populations influencing the development, progression, and treatment resistance of different disease subtypes can potentially uncover a wider range of driver mechanisms for innovative therapeutic interventions.

Deciphering intratumoral heterogeneity using integrated clonal tracking and single-cell transcriptome analyses

Cellular heterogeneity is a major cause of treatment resistance in cancer. Despite recent advances in single-cell genomic and transcriptomic sequencing, it remains difficult to relate measured molecular profiles to the cellular activities underlying cancer. Here, we present an integrated experimental system that connects single cell gene expression to heterogeneous cancer cell growth, metastasis, and treatment response. Our system integrates single cell transcriptome profiling with DNA barcode based clonal tracking in patient-derived xenograft models. We show that leukemia cells exhibiting unique gene expression respond to different chemotherapies in distinct but consistent manners across multiple mice. In addition, we uncover a form of leukemia expansion that is spatially confined to the bone marrow of single anatomical sites and driven by cells with distinct gene expression. Our integrated experimental system can interrogate the molecular and cellular basis of the intratumoral heterogeneity underlying disease progression and treatment resistance.

DNA barcoding is a promising technology for the simultaneous analysis of genetic and phenotypic heterogeneity. Here, the authors combine DNA barcoding and single-cell RNA-seq to study heterogeneity, progression and response to therapy in B-cell acute lymphoblastic leukaemia patient-derived xenografts.

Somatic variant analysis suite: copy number variation clonal visualization online platform for large-scale single-cell genomics

The recent advance of single-cell copy number variation (CNV) analysis plays an essential role in addressing intratumor heterogeneity, identifying tumor subgroups and restoring tumor-evolving trajectories at single-cell scale. Informative visualization of copy number analysis results boosts productive scientific exploration, validation and sharing. Several single-cell analysis figures have the effectiveness of visualizations for understanding single-cell genomics in published articles and software packages. However, they almost lack real-time interaction, and it is hard to reproduce them. Moreover, existing tools are time-consuming and memory-intensive when they reach large-scale single-cell throughputs. We present an online visualization platform, single-cell Somatic Variant Analysis Suite (scSVAS), for real-time interactive single-cell genomics data visualization. scSVAS is specifically designed for large-scale single-cell genomic analysis that provides an arsenal of unique functionalities. After uploading the specified input files, scSVAS deploys the online interactive visualization automatically. Users may conduct scientific discoveries, share interactive visualizations and download high-quality publication-ready figures. scSVAS provides versatile utilities for managing, investigating, sharing and publishing single-cell CNV profiles. We envision this online platform will expedite the biological understanding of cancer clonal evolution in single-cell resolution. All visualizations are publicly hosted at https://sc.deepomics.org.

Single-cell DNA and RNA sequencing reveals the dynamics of intra-tumor heterogeneity in a colorectal cancer model

BACKGROUND

Intra-tumor heterogeneity (ITH) encompasses cellular differences in tumors and is related to clinical outcomes such as drug resistance. However, little is known about the dynamics of ITH, owing to the lack of time-series analysis at the single-cell level. Mouse models that recapitulate cancer development are useful for controlled serial time sampling.

RESULTS

We performed single-cell exome and transcriptome sequencing of 200 cells to investigate how ITH is generated in a mouse colorectal cancer model. In the model, a single normal intestinal cell is grown into organoids that mimic the intestinal crypt structure. Upon RNAi-mediated downregulation of a tumor suppressor gene APC, the transduced organoids were serially transplanted into mice to allow exposure to in vivo microenvironments, which play relevant roles in cancer development. The ITH of the transcriptome increased after the transplantation, while that of the exome decreased. Mutations generated during organoid culture did not greatly change at the bulk-cell level upon the transplantation. The RNA ITH increase was due to the emergence of new transcriptional subpopulations. In contrast to the initial cells expressing mesenchymal-marker genes, new subpopulations repressed these genes after the transplantation. Analyses of colorectal cancer data from The Cancer Genome Atlas revealed a high proportion of metastatic cases in human subjects with expression patterns similar to the new cell subpopulations in mouse. These results suggest that the birth of transcriptional subpopulations may be a key for adaptation to drastic micro-environmental changes when cancer cells have sufficient genetic alterations at later tumor stages.

CONCLUSIONS

This study revealed an evolutionary dynamics of single-cell RNA and DNA heterogeneity in tumor progression, giving insights into the mesenchymal-epithelial transformation of tumor cells at metastasis in colorectal cancer.

Applications of Single-Cell DNA Sequencing

Over the past decade, genomic analyses of single cells-the fundamental units of life-have become possible. Single-cell DNA sequencing has shed light on biological questions that were previously inaccessible across diverse fields of research, including somatic mutagenesis, organismal development, genome function, and microbiology. Single-cell DNA sequencing also promises significant future biomedical and clinical impact, spanning oncology, fertility, and beyond. While single-cell approaches that profile RNA and protein have greatly expanded our understanding of cellular diversity, many fundamental questions in biology and important biomedical applications require analysis of the DNA of single cells. Here, we review the applications and biological questions for which single-cell DNA sequencing is uniquely suited or required. We include a discussion of the fields that will be impacted by single-cell DNA sequencing as the technology continues to advance.

Studying the History of Tumor Evolution from Single-Cell Sequencing Data by Exploring the Space of Binary Matrices

Single-cell sequencing (SCS) data have great potential in reconstructing the evolutionary history of tumors. Rapid advances in SCS technology in the past decade were followed by the design of various computational methods for inferring trees of tumor evolution. Some of the earliest methods were based on the direct search in the space of trees with the goal of finding the maximum likelihood tree. However, it can be shown that instead of searching directly in the tree space, we can perform a search in the space of binary matrices and obtain maximum likelihood tree directly from the maximum likelihood matrix. The potential of the latter tree search strategy has recently been recognized by different research groups and several related methods were published in the past 2 years. Here we provide a review of the theoretical background of these methods and a detailed discussion, which are largely missing in the available publications, of the correlation between the two tree search strategies. We also discuss each of the existing methods based on the search in the space of binary matrices and summarize the best-known single-cell DNA sequencing data sets, which can be used in the future for assessing performance on real data of newly developed methods.

doubletD: detecting doublets in single-cell DNA sequencing data

Motivation

While single-cell DNA sequencing (scDNA-seq) has enabled the study of intratumor heterogeneity at an unprecedented resolution, current technologies are error-prone and often result in doublets where two or more cells are mistaken for a single cell. Not only do doublets confound downstream analyses, but the increase in doublet rate is also a major bottleneck preventing higher throughput with current single-cell technologies. Although doublet detection and removal are standard practice in scRNA-seq data analysis, options for scDNA-seq data are limited. Current methods attempt to detect doublets while also performing complex downstream analyses tasks, leading to decreased efficiency and/or performance.

Results

We present doubletD, the first standalone method for detecting doublets in scDNA-seq data. Underlying our method is a simple maximum likelihood approach with a closed-form solution. We demonstrate the performance of doubletD on simulated data as well as real datasets, outperforming current methods for downstream analysis of scDNA-seq data that jointly infer doublets as well as standalone approaches for doublet detection in scRNA-seq data. Incorporating doubletD in scDNA-seq analysis pipelines will reduce complexity and lead to more accurate results.

Availability and implementation

https://github.com/elkebir-group/doubletD.

Supplementary information

Supplementary data are available at Bioinformatics online.

A BAYESIAN NONPARAMETRIC MODEL FOR INFERRING SUBCLONAL POPULATIONS FROM STRUCTURED DNA SEQUENCING DATA

There are distinguishing features or "hallmarks" of cancer that are found across tumors, individuals, and types of cancer, and these hallmarks can be driven by specific genetic mutations. Yet, within a single tumor there is often extensive genetic heterogeneity as evidenced by single-cell and bulk DNA sequencing data. The goal of this work is to jointly infer the underlying genotypes of tumor subpopulations and the distribution of those subpopulations in individual tumors by integrating single-cell and bulk sequencing data. Understanding the genetic composition of the tumor at the time of treatment is important in the personalized design of targeted therapeutic combinations and monitoring for possible recurrence after treatment. We propose a hierarchical Dirichlet process mixture model that incorporates the correlation structure induced by a structured sampling arrangement and we show that this model improves the quality of inference. We develop a representation of the hierarchical Dirichlet process prior as a Gamma-Poisson hierarchy and we use this representation to derive a fast Gibbs sampling inference algorithm using the augment-and-marginalize method. Experiments with simulation data show that our model outperforms standard numerical and statistical methods for decomposing admixed count data. Analyses of real acute lymphoblastic leukemia cancer sequencing dataset shows that our model improves upon state-of-the-art bioinformatic methods. An interpretation of the results of our model on this real dataset reveals co-mutated loci across samples.

Exploiting Single-Cell Tools in Gene and Cell Therapy

Single-cell molecular tools have been developed at an incredible pace over the last five years as sequencing costs continue to drop and numerous molecular assays have been coupled to sequencing readouts. This rapid period of technological development has facilitated the delineation of individual molecular characteristics including the genome, transcriptome, epigenome, and proteome of individual cells, leading to an unprecedented resolution of the molecular networks governing complex biological systems. The immense power of single-cell molecular screens has been particularly highlighted through work in systems where cellular heterogeneity is a key feature, such as stem cell biology, immunology, and tumor cell biology. Single-cell-omics technologies have already contributed to the identification of novel disease biomarkers, cellular subsets, therapeutic targets and diagnostics, many of which would have been undetectable by bulk sequencing approaches. More recently, efforts to integrate single-cell multi-omics with single cell functional output and/or physical location have been challenging but have led to substantial advances. Perhaps most excitingly, there are emerging opportunities to reach beyond the description of static cellular states with recent advances in modulation of cells through CRISPR technology, in particular with the development of base editors which greatly raises the prospect of cell and gene therapies. In this review, we provide a brief overview of emerging single-cell technologies and discuss current developments in integrating single-cell molecular screens and performing single-cell multi-omics for clinical applications. We also discuss how single-cell molecular assays can be usefully combined with functional data to unpick the mechanism of cellular decision-making. Finally, we reflect upon the introduction of spatial transcriptomics and proteomics, its complementary role with single-cell RNA sequencing (scRNA-seq) and potential application in cellular and gene therapy.

Distinguishing linear and branched evolution given single-cell DNA sequencing data of tumors

Background

Cancer arises from an evolutionary process where somatic mutations give rise to clonal expansions. Reconstructing this evolutionary process is useful for treatment decision-making as well as understanding evolutionary patterns across patients and cancer types. In particular, classifying a tumor’s evolutionary process as either linear or branched and understanding what cancer types and which patients have each of these trajectories could provide useful insights for both clinicians and researchers. While comprehensive cancer phylogeny inference from single-cell DNA sequencing data is challenging due to limitations with current sequencing technology and the complexity of the resulting problem, current data might provide sufficient signal to accurately classify a tumor’s evolutionary history as either linear or branched.

Results

We introduce the Linear Perfect Phylogeny Flipping (LPPF) problem as a means of testing two alternative hypotheses for the pattern of evolution, which we prove to be NP-hard. We develop Phyolin, which uses constraint programming to solve the LPPF problem. Through both in silico experiments and real data application, we demonstrate the performance of our method, outperforming a competing machine learning approach.

Conclusion

Phyolin is an accurate, easy to use and fast method for classifying an evolutionary trajectory as linear or branched given a tumor’s single-cell DNA sequencing data.

Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation, or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108 188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (>35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had >2 NOTCH1 mutations, some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells and clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Therefore, single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

The Bourque distances for mutation trees of cancers

Background

Mutation trees are rooted trees in which nodes are of arbitrary degree and labeled with a mutation set. These trees, also referred to as clonal trees, are used in computational oncology to represent the mutational history of tumours. Classical tree metrics such as the popular Robinson–Foulds distance are of limited use for the comparison of mutation trees. One reason is that mutation trees inferred with different methods or for different patients often contain different sets of mutation labels.

Results

We generalize the Robinson–Foulds distance into a set of distance metrics called Bourque distances for comparing mutation trees. We show the basic version of the Bourque distance for mutation trees can be computed in linear time. We also make a connection between the Robinson–Foulds distance and the nearest neighbor interchange distance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13015-021-00188-3.

Accurate genomic variant detection in single cells with primary template-directed amplification

Improvements in whole genome amplification (WGA) would enable new types of basic and applied biomedical research, including studies of intratissue genetic diversity that require more accurate single-cell genotyping. Here, we present primary template-directed amplification (PTA), an isothermal WGA method that reproducibly captures >95% of the genomes of single cells in a more uniform and accurate manner than existing approaches, resulting in significantly improved variant calling sensitivity and precision. To illustrate the types of studies that are enabled by PTA, we developed direct measurement of environmental mutagenicity (DMEM), a tool for mapping genome-wide interactions of mutagens with single living human cells at base-pair resolution. In addition, we utilized PTA for genome-wide off-target indel and structural variant detection in cells that had undergone CRISPR-mediated genome editing, establishing the feasibility for performing single-cell evaluations of biopsies from edited tissues. The improved precision and accuracy of variant detection with PTA overcomes the current limitations of accurate WGA, which is the major obstacle to studying genetic diversity and evolution at cellular resolution.

GRMT: Generative Reconstruction of Mutation Tree From Scratch Using Single-Cell Sequencing Data

Single-cell sequencing (SCS) now promises the landscape of genetic diversity at single cell level, and is particularly useful to reconstruct the evolutionary history of tumor. There are multiple types of noise that make the SCS data notoriously error-prone, and significantly complicate tumor tree reconstruction. Existing methods for tumor phylogeny estimation suffer from either high computational intensity or low-resolution indication of clonal architecture, giving a necessity of developing new methods for efficient and accurate reconstruction of tumor trees. We introduce GRMT (Generative Reconstruction of Mutation Tree from scratch), a method for inferring tumor mutation tree from SCS data. GRMT exploits the k-Dollo parsimony model to allow each mutation to be gained once and lost at most k times. Under this constraint on mutation evolution, GRMT searches for mutation tree structures from a perspective of tree generation from scratch, and implements it to an iterative process that gradually increases the tree size by introducing a new mutation per time until a complete tree structure that contains all mutations is obtained. This enables GRMT to efficiently recover the chronological order of mutations and scale well to large datasets. Extensive evaluations on simulated and real datasets suggest GRMT outperforms the state-of-the-arts in multiple performance metrics. The GRMT software is freely available at https://github.com/qasimyu/grmt.

Direct Assessment of Single-Cell DNA Using Crudely Purified Live Cells: A Proof of Concept for Noninvasive Prenatal Definitive Diagnosis

Noninvasive testing techniques are often used for fetal diagnosis of genetic abnormalities but are limited by certain characteristics, including noninformative results. Thus, novel methods of noninvasive definitive diagnosis of fetal genetic abnormalities are needed. The aim of this study was to develop a single-cell DNA analysis method with high sensitivity and specificity that enables direct extraction of genetic information from live fetal cells in a crude mixture for simultaneous evaluation. Genomic DNA from circulating fetal CD45^-CD14^- cells, an extremely rare cell type, extracted from 10-mL samples of maternal peripheral blood, was extracted using a single-cell-based droplet digital (sc-dd) PCR system with a modified amount of polymerase. A hexachloro-6-carboxyfluorescein-labeled RPP30 probe was used as an internal control and a 6-carboxyfluorescein-labeled SRY probe as a target. The results indicated that no droplets generated with samples from pregnant women carrying female fetuses were positive for both probe signals, whereas droplets prepared with samples from pregnant women carrying male fetuses were positive for both probe signals. The latter was considered a direct assessment of genetic information from single circulating male fetal cells. Thus, the modified sc-ddPCR system allows the detection of genetic information from rare target cells in a crudely purified cell population. This research also serves as a proof of concept for noninvasive prenatal definitive diagnosis.

Single-cell technologies and analyses in hematopoiesis and hematological malignancies

In recent years, single-cell technologies have emerged as breakthrough techniques that enable the characterization of hematopoietic cell populations of normal and malignant tissue samples and will be combined in the near future with bulk technologies, currently used in clinical practice, to improve diagnosis, prognosis, and the search for novel molecular targets. These single-cell methods have the advantage of not masking cell-to-cell variation features and involve the study of genetic, epigenetic, transcriptional, and proteomic landscapes from a single-cell perspective. Latest advances in this field have enabled the development of novel strategies that significantly increase both sensitivity and high throughput. In this review, we emphasize emerging techniques aimed at assessing individual or multiomic parameters at single-cell resolution and analyze how these technologies have helped us understand hematopoietic variability and identify unknown and/or rare subpopulations. We also summarize the impact of these single-cell profiling strategies on the characterization of cell diversity within the tumor and the clonal evolution of multiple hematological malignancies in samples from untreated and treated patients, which provide valuable information for diagnosis, prognosis, and future treatments and explain why current therapies may fail. However, despite these improvements, new challenges lie ahead.

The evolution of leukaemia from pre-leukaemic and leukaemic stem cells

Haematopoietic stem and progenitor cells (HSPCs) are defined as unspecialized cells that give rise to more differentiated cells. In a similar way, leukaemic stem and progenitor cells (LSPCs) are defined as unspecialized leukaemic cells, which can give rise to more differentiated cells. Leukaemic cells carry leukaemic mutations/variants and have clear differentiation abnormalities. Pre-leukaemic HSPCs (PreL-HSPCs) carry pre-leukaemic mutations/variants (pLMs) and are capable of producing mature functional cells, which will carry the same variants. Under the roof of LSPCs, one can find a broad range of cell types genetic and disease phenotypes. Present-day knowledge suggests that this phenotypic heterogeneity is the result of interactions between the cell of origin, the genetic background and the microenvironment background. The combination of these attributes will define the LSPC phenotype, frequency, differentiation capacity and evolutionary trajectory. Importantly, as LSPCs are leukaemia-initiating cells that sustain clinical remission and are the source of relapse, an improved understanding of LSPCs phenotype would offer better clinical opportunities for the treatment and hopefully prevention of human leukaemia. The current review will focus on LSPCs attributes in the context of human haematologic malignancies.

Inferring cancer progression from Single-Cell Sequencing while allowing mutation losses

Motivation

In recent years, the well-known Infinite Sites Assumption has been a fundamental feature of computational methods devised for reconstructing tumor phylogenies and inferring cancer progressions. However, recent studies leveraging single-cell sequencing (SCS) techniques have shown evidence of the widespread recurrence and, especially, loss of mutations in several tumor samples. While there exist established computational methods that infer phylogenies with mutation losses, there remain some advancements to be made.

Results

We present Simulated Annealing Single-Cell inference (SASC): a new and robust approach based on simulated annealing for the inference of cancer progression from SCS datasets. In particular, we introduce an extension of the model of evolution where mutations are only accumulated, by allowing also a limited amount of mutation loss in the evolutionary history of the tumor: the Dollo-k model. We demonstrate that SASC achieves high levels of accuracy when tested on both simulated and real datasets and in comparison with some other available methods.

Availability and implementation

The SASC tool is open source and available at https://github.com/sciccolella/sasc.

Supplementary information

Supplementary data are available at Bioinformatics online.

Precision Oncology Beyond Genomics: The Future Is Here-It Is Just Not Evenly Distributed

Cancer is a multifactorial disease with increasing incidence. There are more than 100 different cancer types, defined by location, cell of origin, and genomic alterations that influence oncogenesis and therapeutic response. This heterogeneity between tumors of different patients and also the heterogeneity within the same patient's tumor pose an enormous challenge to cancer treatment. In this review, we explore tumor heterogeneity on the longitudinal and the latitudinal axis, reviewing current and future approaches to study this heterogeneity and their potential to support oncologists in tailoring a patient's treatment regimen. We highlight how the ideal of precision oncology is reaching far beyond the knowledge of genetic variants to inform clinical practice and discuss the technologies and strategies already available to improve our understanding and management of heterogeneity in cancer treatment. We will focus on integrating multi-omics technologies with suitable in vitro models and their proficiency in mimicking endogenous tumor heterogeneity.

Moss enables high sensitivity single-nucleotide variant calling from multiple bulk DNA tumor samples

Intra-tumor heterogeneity renders the identification of somatic single-nucleotide variants (SNVs) a challenging problem. In particular, low-frequency SNVs are hard to distinguish from sequencing artifacts. While the increasing availability of multi-sample tumor DNA sequencing data holds the potential for more accurate variant calling, there is a lack of high-sensitivity multi-sample SNV callers that utilize these data. Here we report Moss, a method to identify low-frequency SNVs that recur in multiple sequencing samples from the same tumor. Moss provides any existing single-sample SNV caller the ability to support multiple samples with little additional time overhead. We demonstrate that Moss improves recall while maintaining high precision in a simulated dataset. On multi-sample hepatocellular carcinoma, acute myeloid leukemia and colorectal cancer datasets, Moss identifies new low-frequency variants that meet manual review criteria and are consistent with the tumor's mutational signature profile. In addition, Moss detects the presence of variants in more samples of the same tumor than reported by the single-sample caller. Moss' improved sensitivity in SNV calling will enable more detailed downstream analyses in cancer genomics.

Time to Move to the Single-Cell Level: Applications of Single-Cell Multi-Omics to Hematological Malignancies and Waldenström's Macroglobulinemia-A Particularly Heterogeneous Lymphoma

Single-cell sequencing techniques have become a powerful tool for characterizing intra-tumor heterogeneity, which has been reflected in the increasing number of studies carried out and reported. We have rigorously reviewed and compiled the information about these techniques inasmuch as they are relative to the area of hematology to provide a practical view of their potential applications. Studies show how single-cell multi-omics can overcome the limitations of bulk sequencing and be applied at all stages of tumor development, giving insights into the origin and pathogenesis of the tumors, the clonal architecture and evolution, or the mechanisms of therapy resistance. Information at the single-cell level may help resolve questions related to intra-tumor heterogeneity that have not been previously explained by other techniques. With that in mind, we review the existing knowledge about a heterogeneous lymphoma called Waldenström's macroglobulinemia and discuss how single-cell studies may help elucidate the underlying causes of this heterogeneity.

CACTUS: integrating clonal architecture with genomic clustering and transcriptome profiling of single tumor cells

Background

Drawing genotype-to-phenotype maps in tumors is of paramount importance for understanding tumor heterogeneity. Assignment of single cells to their tumor clones of origin can be approached by matching the genotypes of the clones to the mutations found in RNA sequencing of the cells. The confidence of the cell-to-clone mapping can be increased by accounting for additional measurements. Follicular lymphoma, a malignancy of mature B cells that continuously acquire mutations in parallel in the exome and in B cell receptor loci, presents a unique opportunity to join exome-derived mutations with B cell receptor sequences as independent sources of evidence for clonal evolution.

Methods

Here, we propose CACTUS, a probabilistic model that leverages the information from an independent genomic clustering of cells and exploits the scarce single cell RNA sequencing data to map single cells to given imperfect genotypes of tumor clones.

Results

We apply CACTUS to two follicular lymphoma patient samples, integrating three measurements: whole exome, single-cell RNA, and B cell receptor sequencing. CACTUS outperforms a predecessor model by confidently assigning cells and B cell receptor-based clusters to the tumor clones.

Conclusions

The integration of independent measurements increases model certainty and is the key to improving model performance in the challenging task of charting the genotype-to-phenotype maps in tumors. CACTUS opens the avenue to study the functional implications of tumor heterogeneity, and origins of resistance to targeted therapies. CACTUS is written in R and source code, along with all supporting files, are available on GitHub (https://github.com/LUMC/CACTUS).

Supplementary Information

The online version contains supplementary material available at (10.1186/s13073-021-00842-w).

Genetic and Non-Genetic Mechanisms Underlying Cancer Evolution

Simple Summary

Our manuscript summarizes the up-to-date data on the complex and dynamic nature of adaptation mechanisms and evolutionary processes taking place during cancer initiation, development and progression. Although for decades cancer has been viewed as a process governed by genetic mechanisms, it is becoming more and more clear that non-genetic mechanisms may play an equally important role in cancer evolution. In this review, we bring together these fundamental concepts and discuss how those tightly interconnected mechanisms lead to the establishment of highly adaptive quickly evolving cancers. Furthermore, we argue that in depth understanding of cancer progression from the evolutionary perspective may allow the prediction and direction of the evolutionary path of cancer populations towards drug sensitive phenotypes and thus facilitate the development of more effective anti-cancer approaches.

Abstract

Cancer development can be defined as a process of cellular and tissular microevolution ultimately leading to malignancy. Strikingly, though this concept has prevailed in the field for more than a century, the precise mechanisms underlying evolutionary processes occurring within tumours remain largely uncharacterized and rather cryptic. Nevertheless, although our current knowledge is fragmentary, data collected to date suggest that most tumours display features compatible with a diverse array of evolutionary paths, suggesting that most of the existing macro-evolutionary models find their avatar in cancer biology. Herein, we discuss an up-to-date view of the fundamental genetic and non-genetic mechanisms underlying tumour evolution with the aim of concurring into an integrated view of the evolutionary forces at play throughout the emergence and progression of the disease and into the acquisition of resistance to diverse therapeutic paradigms. Our ultimate goal is to delve into the intricacies of genetic and non-genetic networks underlying tumour evolution to build a framework where both core concepts are considered non-negligible and equally fundamental.

BnpC: Bayesian non-parametric clustering of single-cell mutation profiles

Motivation

The high resolution of single-cell DNA sequencing (scDNA-seq) offers great potential to resolve intratumor heterogeneity (ITH) by distinguishing clonal populations based on their mutation profiles. However, the increasing size of scDNA-seq datasets and technical limitations, such as high error rates and a large proportion of missing values, complicate this task and limit the applicability of existing methods.

Results

Here, we introduce BnpC, a novel non-parametric method to cluster individual cells into clones and infer their genotypes based on their noisy mutation profiles. We benchmarked our method comprehensively against state-of-the-art methods on simulated data using various data sizes, and applied it to three cancer scDNA-seq datasets. On simulated data, BnpC compared favorably against current methods in terms of accuracy, runtime and scalability. Its inferred genotypes were the most accurate, especially on highly heterogeneous data, and it was the only method able to run and produce results on datasets with 5000 cells. On tumor scDNA-seq data, BnpC was able to identify clonal populations missed by the original cluster analysis but supported by Supplementary Experimental Data. With ever growing scDNA-seq datasets, scalable and accurate methods such as BnpC will become increasingly relevant, not only to resolve ITH but also as a preprocessing step to reduce data size.

Availability and implementation

BnpC is freely available under MIT license at https://github.com/cbg-ethz/BnpC.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-cell sequencing in translational cancer research and challenges to meet clinical diagnostic needs

The ability to capture alterations in the genome or transcriptome by next-generation sequencing has provided critical insight into molecular changes and programs underlying cancer biology. With the rapid technological development in single-cell sequencing, it has become possible to study individual cells at the transcriptional, genetic, epigenetic, and protein level. Using single-cell analysis, an increased resolution of fundamental processes underlying cancer development is obtained, providing comprehensive insights otherwise lost by sequencing of entire (bulk) samples, in which molecular signatures of individual cells are averaged across the entire cell population. Here, we provide a concise overview on the application of single-cell analysis of different modalities within cancer research by highlighting key articles of their respective fields. We furthermore examine the potential of existing technologies to meet clinical diagnostic needs and discuss current challenges associated with this translation.

Melanoma Single-Cell Biology in Experimental and Clinical Settings

Cellular heterogeneity is regarded as a major factor for treatment response and resistance in a variety of malignant tumors, including malignant melanoma. More recent developments of single-cell sequencing technology provided deeper insights into this phenomenon. Single-cell data were used to identify prognostic subtypes of melanoma tumors, with a special emphasis on immune cells and fibroblasts in the tumor microenvironment. Moreover, treatment resistance to checkpoint inhibitor therapy has been shown to be associated with a set of differentially expressed immune cell signatures unraveling new targetable intracellular signaling pathways. Characterization of T cell states under checkpoint inhibitor treatment showed that exhausted CD8+ T cell types in melanoma lesions still have a high proliferative index. Other studies identified treatment resistance mechanisms to targeted treatment against the mutated BRAF serine/threonine protein kinase including repression of the melanoma differentiation gene microphthalmia-associated transcription factor (MITF) and induction of AXL receptor tyrosine kinase. Interestingly, treatment resistance mechanisms not only included selection processes of pre-existing subclones but also transition between different states of gene expression. Taken together, single-cell technology has provided deeper insights into melanoma biology and has put forward our understanding of the role of tumor heterogeneity and transcriptional plasticity, which may impact on innovative clinical trial designs and experimental approaches.

Cytologic and Molecular Diagnostics for Vitreoretinal Lymphoma: Current Approaches and Emerging Single-Cell Analyses

Vitreoretinal lymphoma (VRL) is a rare ocular malignancy that manifests as diffuse large B-cell lymphoma. Early and accurate diagnosis is essential to prevent mistreatment and to reduce the high morbidity and mortality associated with VRL. The disease can be diagnosed using various methods, including cytology, immunohistochemistry, cytokine analysis, flow cytometry, and molecular analysis of bulk vitreous aspirates. Despite these options, VRL diagnosis remains challenging, as samples are often confounded by low cellularity, the presence of debris and non-target immunoreactive cells, and poor cytological preservation. As such, VRL diagnostic accuracy is limited by both false-positive and false-negative outcomes. Missed or inappropriate diagnosis may cause delays in treatment, which can have life-threatening consequences for patients with VRL. In this review, we summarize current knowledge and the diagnostic modalities used for VRL diagnosis. We also highlight several emerging molecular techniques, including high-resolution single cell-based analyses, which may enable more comprehensive and precise VRL diagnoses.

Applying Single-Cell Technology in Uveal Melanomas: Current Trends and Perspectives for Improving Uveal Melanoma Metastasis Surveillance and Tumor Profiling

Uveal melanoma (UM) is the most common primary adult intraocular malignancy. This rare but devastating cancer causes vision loss and confers a poor survival rate due to distant metastases. Identifying clinical and molecular features that portend a metastatic risk is an important part of UM workup and prognostication. Current UM prognostication tools are based on determining the tumor size, gene expression profile, and chromosomal rearrangements. Although we can predict the risk of metastasis fairly accurately, we cannot obtain preclinical evidence of metastasis or identify biomarkers that might form the basis of targeted therapy. These gaps in UM research might be addressed by single-cell research. Indeed, single-cell technologies are being increasingly used to identify circulating tumor cells and profile transcriptomic signatures in single, drug-resistant tumor cells. Such advances have led to the identification of suitable biomarkers for targeted treatment. Here, we review the approaches used in cutaneous melanomas and other cancers to isolate single cells and profile them at the transcriptomic and/or genomic level. We discuss how these approaches might enhance our current approach to UM management and review the emerging data from single-cell analyses in UM.

Integrated decoding hematopoiesis and leukemogenesis using single-cell sequencing and its medical implication

Single-cell RNA sequencing provides exciting opportunities to unbiasedly study hematopoiesis. However, our understanding of leukemogenesis was limited due to the high individual differences. Integrated analyses of hematopoiesis and leukemogenesis potentially provides new insights. Here we analyzed ~200,000 single-cell transcriptomes of bone marrow mononuclear cells (BMMCs) and its subsets from 23 clinical samples. We constructed a comprehensive cell atlas as hematopoietic reference. We developed counterpart composite index (CCI; available at GitHub: https://github.com/pengfeeei/cci) to search for the healthy counterpart of each leukemia cell subpopulation, by integrating multiple statistics to map leukemia cells onto reference hematopoietic cells. Interestingly, we found leukemia cell subpopulations from each patient had different healthy counterparts. Analysis showed the trajectories of leukemia cell subpopulations were similar to that of their healthy counterparts, indicating that developmental termination of leukemia initiating cells at different phases leads to different leukemia cell subpopulations thus explained the origin of leukemia heterogeneity. CCI further predicts leukemia subtypes, cellular heterogeneity, and cellular stemness of each leukemia patient. Analyses of leukemia patient at diagnosis, refractory, remission and relapse vividly presented dynamics of cell population during leukemia treatment. CCI analyses showed the healthy counterparts of relapsed leukemia cells were closer to the root of hematopoietic tree than that of other leukemia cells, although single-cell transcriptomic genetic variants and haplotype tracing analyses showed the relapsed leukemia cell were derived from an early minor leukemia cell population. In summary, this study developed a unified framework for understanding leukemogenesis with hematopoiesis reference, which provided novel biological and medical implication.

Single-cell omics analyses with single molecular detection: challenges and perspectives

The ultimate goal of single-cell analyses is to obtain the biomolecular content for each cell in unicellular and multicellular organisms at different points of their life cycle under variable environmental conditions. These require an assessment of: a) the total number of cells, b) the total number of cell types, and c) the complete and quantitative single molecular detection and identification for all classes of biopolymers, and organic and inorganic compounds, in each individual cell. For proteins, glycans, lipids, and metabolites, whose sequences cannot be amplified by copying as in the case of nucleic acids, the detection limit by mass spectrometry is about 10⁵ molecules. Therefore, proteomic, glycomic, lipidomic, and metabolomic analyses do not yet permit the assembly of the complete single-cell omes. The construction of novel nanoelectrophoretic arrays and nano in microarrays on a single 1-cm-diameter chip has shown proof of concept for a high throughput platform for parallel processing of thousands of individual cells. Combined with dynamic secondary ion mass spectrometry, with 3D scanning capability and lateral resolution of 50 nm, the sensitivity of single molecular quantification and identification for all classes of biomolecules could be reached. Further development and routine application of such technological and instrumentation solution would allow assembly of complete omes with a quantitative assessment of structural and functional cellular diversity at the molecular level.

Red panda: a novel method for detecting variants in single-cell RNA sequencing

BACKGROUND

Single-cell sequencing enables us to better understand genetic diseases, such as cancer or autoimmune disorders, which are often affected by changes in rare cells. Currently, no existing software is aimed at identifying single nucleotide variations or micro (1-50 bp) insertions and deletions in single-cell RNA sequencing (scRNA-seq) data. Generating high-quality variant data is vital to the study of the aforementioned diseases, among others.

RESULTS

In this study, we report the design and implementation of Red Panda, a novel method to accurately identify variants in scRNA-seq data. Variants were called on scRNA-seq data from human articular chondrocytes, mouse embryonic fibroblasts (MEFs), and simulated data stemming from the MEF alignments. Red Panda had the highest Positive Predictive Value at 45.0%, while other tools-FreeBayes, GATK HaplotypeCaller, GATK UnifiedGenotyper, Monovar, and Platypus-ranged from 5.8-41.53%. From the simulated data, Red Panda had the highest sensitivity at 72.44%.

CONCLUSIONS

We show that our method provides a novel and improved mechanism to identify variants in scRNA-seq as compared to currently existing software. However, methods for identification of genomic variants using scRNA-seq data can be still improved.

Innate Lymphoid Cells and Celiac Disease: Current Perspective

Celiac disease (CD) is a common autoimmune disorder triggered by the ingestion of gluten in genetically susceptible individuals. Although the mechanisms underlying gliadin-mediated activation of adaptive immunity in CD have been well-characterized, regulation of innate immune responses and the functions of certain immune cell populations within the epithelium and lamina propria are not well-understood at present. Innate lymphoid cells (ILCs) are types of innate immune cells that have lymphoid morphology, lack antigen-specific receptors, and play important roles in tissue homeostasis, inflammation, and protective immune responses against pathogens. Information regarding the diversity and functions of ILCs in lymphoid organs and at mucosal sites has grown over the past decade, and roles of different ILC subsets in the pathogenesis of some inflammatory intestinal diseases have been proposed. However, our understanding of the contribution of ILCs toward the initiation and progression of CD is still limited. In this review, we discuss current pathophysiological aspects of ILCs within the gastrointestinal tract, findings of recent investigations characterizing ILC alterations in CD and refractory CD, and suggest avenues for future research.

gpps: an ILP-based approach for inferring cancer progression with mutation losses from single cell data

BACKGROUND

Cancer progression reconstruction is an important development stemming from the phylogenetics field. In this context, the reconstruction of the phylogeny representing the evolutionary history presents some peculiar aspects that depend on the technology used to obtain the data to analyze: Single Cell DNA Sequencing data have great specificity, but are affected by moderate false negative and missing value rates. Moreover, there has been some recent evidence of back mutations in cancer: this phenomenon is currently widely ignored.

RESULTS

We present a new tool, gpps, that reconstructs a tumor phylogeny from Single Cell Sequencing data, allowing each mutation to be lost at most a fixed number of times. The General Parsimony Phylogeny from Single cell (gpps) tool is open source and available at https://github.com/AlgoLab/gpps .

CONCLUSIONS

gpps provides new insights to the analysis of intra-tumor heterogeneity by proposing a new progression model to the field of cancer phylogeny reconstruction on Single Cell data.

Digital-WGS: Automated, highly efficient whole-genome sequencing of single cells by digital microfluidics

Single-cell whole-genome sequencing (WGS) is critical for characterizing dynamic intercellular changes in DNA. Current sample preparation technologies for single-cell WGS are complex, expensive, and suffer from high amplification bias and errors. Here, we describe Digital-WGS, a sample preparation platform that streamlines high-performance single-cell WGS with automatic processing based on digital microfluidics. Using the method, we provide high single-cell capture efficiency for any amount and types of cells by a wetted hydrodynamic structure. The digital control of droplets in a closed hydrophobic interface enables the complete removal of exogenous DNA, sufficient cell lysis, and lossless amplicon recovery, achieving the low coefficient of variation and high coverage at multiple scales. The single-cell genomic variations profiling performs the excellent detection of copy number variants with the smallest bin of 150 kb and single-nucleotide variants with allele dropout rate of 5.2%, holding great promise for broader applications of single-cell genomics.

The Trends of Single-Cell Analysis: A Global Study

Objective

The field of single-cell analysis has rapidly grown worldwide, and a bibliometric analysis and visualization of data and publications pertaining to such single-cell research has the potential to offer insights into the development of this field over the past two decades while also highlighting future avenues of research.

Methods

Single-cell analysis-related studies published from 2000-2019 were identified through searches of the Web of Science, Scopus, and PubMed databases, and corresponding bibliometric data were systematically compiled. Extracted data from each study included author names, country of origin, and affiliations. GraphPad Prism was used to analyze these data, while VOSviewer was used to perform global analyses of bibliographic coupling, coauthorship, cocitation, and co-occurrence.

Results

In total, 4,071 relevant studies were included in this analysis. The number of publications increased substantially with time, suggesting that single-cell analyses are becoming increasingly more prevalent in recent years. Studies from the USA had the greatest impact in this field, with higher H-index values and numbers of citations relative to other countries, whereas Israel exhibited the highest average number of citations per publication. Bibliographic coupling, coauthorship, cocitation, and co-occurrence analyses revealed that Analytical Chemistry was associated with the highest number of publications in this field, and the University of Stanford contributed the most to this field. The most cited study included in this analysis was published by Macosko et al. in 2015 in Cell. Co-occurrence analyses revealed that the most common single-cell research topics included “mechanistic studies,” “in vitro studies,” “in vivo studies,” and “fabrication studies.”

Conclusions

Single-cell analyses are a rapidly growing area of scientific interest, and higher volumes of publications in this field are expected in the coming years, particularly for studies conducting fabrication and in vivo single-cell analyses.

Tumor Phylogeny Topology Inference via Deep Learning

Principled computational approaches for tumor phylogeny reconstruction via single-cell sequencing typically aim to build the most likely perfect phylogeny tree from the noisy genotype matrix - which represents genotype calls of single cells. This problem is NP-hard, and as a result, existing approaches aim to solve relatively small instances of it through combinatorial optimization techniques or Bayesian inference. As expected, even when the goal is to infer basic topological features of the tumor phylogeny, rather than reconstructing the topology entirely, these approaches could be prohibitively slow. In this paper, we introduce fast deep learning solutions to the problems of inferring whether the most likely tree has a linear (chain) or branching topology and whether a perfect phylogeny is feasible from a given genotype matrix. We also present a reinforcement learning approach for reconstructing the most likely tumor phylogeny. This preliminary work demonstrates that data-driven approaches can reconstruct key features of tumor evolution.

Inference of mutability landscapes of tumors from single cell sequencing data

One of the hallmarks of cancer is the extremely high mutability and genetic instability of tumor cells. Inherent heterogeneity of intra-tumor populations manifests itself in high variability of clone instability rates. Analogously to fitness landscapes, the instability rates of clonal populations form their mutability landscapes. Here, we present MULAN (MUtability LANdscape inference), a maximum-likelihood computational framework for inference of mutation rates of individual cancer subclones using single-cell sequencing data. It utilizes the partial information about the orders of mutation events provided by cancer mutation trees and extends it by inferring full evolutionary history and mutability landscape of a tumor. Evaluation of mutation rates on the level of subclones rather than individual genes allows to capture the effects of genomic interactions and epistasis. We estimate the accuracy of our approach and demonstrate that it can be used to study the evolution of genetic instability and infer tumor evolutionary history from experimental data. MULAN is available at https://github.com/compbel/MULAN.

RobustClone: a robust PCA method for tumor clone and evolution inference from single-cell sequencing data

MOTIVATION

Single-cell sequencing (SCS) data provide unprecedented insights into intratumoral heterogeneity. With SCS, we can better characterize clonal genotypes and reconstruct phylogenetic relationships of tumor cells/clones. However, SCS data are often error-prone, making their computational analysis challenging.

RESULTS

To infer the clonal evolution in tumor from the error-prone SCS data, we developed an efficient computational framework, termed RobustClone. It recovers the true genotypes of subclones based on the extended robust principal component analysis, a low-rank matrix decomposition method, and reconstructs the subclonal evolutionary tree. RobustClone is a model-free method, which can be applied to both single-cell single nucleotide variation (scSNV) and single-cell copy-number variation (scCNV) data. It is efficient and scalable to large-scale datasets. We conducted a set of systematic evaluations on simulated datasets and demonstrated that RobustClone outperforms state-of-the-art methods in large-scale data both in accuracy and efficiency. We further validated RobustClone on two scSNV and two scCNV datasets and demonstrated that RobustClone could recover genotype matrix and infer the subclonal evolution tree accurately under various scenarios. In particular, RobustClone revealed the spatial progression patterns of subclonal evolution on the large-scale 10X Genomics scCNV breast cancer dataset.

AVAILABILITY AND IMPLEMENTATION

RobustClone software is available at https://github.com/ucasdp/RobustClone.

CONTACT

lwan@amss.ac.cn or maliang@ioz.ac.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Highly Scalable Method for Joint Whole-Genome Sequencing and Gene-Expression Profiling of Single Cells

To address how genetic variation alters gene expression in complex cell mixtures, we developed direct nuclear tagmentation and RNA sequencing (DNTR-seq), which enables whole-genome and mRNA sequencing jointly in single cells. DNTR-seq readily identified minor subclones within leukemia patients. In a large-scale DNA damage screen, DNTR-seq was used to detect regions under purifying selection and identified genes where mRNA abundance was resistant to copy-number alteration, suggesting strong genetic compensation. mRNA sequencing (mRNA-seq) quality equals RNA-only methods, and the low positional bias of genomic libraries allowed detection of sub-megabase aberrations at ultra-low coverage. Each cell library is individually addressable and can be re-sequenced at increased depth, allowing multi-tiered study designs. Additionally, the direct tagmentation protocol enables coverage-independent estimation of ploidy, which can be used to identify cell singlets. Thus, DNTR-seq directly links each cell's state to its corresponding genome at scale, enabling routine analysis of heterogeneous tumors and other complex tissues.

PhyDOSE: Design of follow-up single-cell sequencing experiments of tumors

The combination of bulk and single-cell DNA sequencing data of the same tumor enables the inference of high-fidelity phylogenies that form the input to many important downstream analyses in cancer genomics. While many studies simultaneously perform bulk and single-cell sequencing, some studies have analyzed initial bulk data to identify which mutations to target in a follow-up single-cell sequencing experiment, thereby decreasing cost. Bulk data provide an additional untapped source of valuable information, composed of candidate phylogenies and associated clonal prevalence. Here, we introduce PhyDOSE, a method that uses this information to strategically optimize the design of follow-up single cell experiments. Underpinning our method is the observation that only a small number of clones uniquely distinguish one candidate tree from all other trees. We incorporate distinguishing features into a probabilistic model that infers the number of cells to sequence so as to confidently reconstruct the phylogeny of the tumor. We validate PhyDOSE using simulations and a retrospective analysis of a leukemia patient, concluding that PhyDOSE's computed number of cells resolves tree ambiguity even in the presence of typical single-cell sequencing errors. We also conduct a retrospective analysis on an acute myeloid leukemia cohort, demonstrating the potential to achieve similar results with a significant reduction in the number of cells sequenced. In a prospective analysis, we demonstrate the advantage of selecting cells to sequence across multiple biopsies and that only a small number of cells suffice to disambiguate the solution space of trees in a recent lung cancer cohort. In summary, PhyDOSE proposes cost-efficient single-cell sequencing experiments that yield high-fidelity phylogenies, which will improve downstream analyses aimed at deepening our understanding of cancer biology.

Accurate and efficient cell lineage tree inference from noisy single cell data: the maximum likelihood perfect phylogeny approach

MOTIVATION

Cells in an organism share a common evolutionary history, called cell lineage tree. Cell lineage tree can be inferred from single cell genotypes at genomic variation sites. Cell lineage tree inference from noisy single cell data is a challenging computational problem. Most existing methods for cell lineage tree inference assume uniform uncertainty in genotypes. A key missing aspect is that real single cell data usually has non-uniform uncertainty in individual genotypes. Moreover, existing methods are often sampling based and can be very slow for large data.

RESULTS

In this article, we propose a new method called ScisTree, which infers cell lineage tree and calls genotypes from noisy single cell genotype data. Different from most existing approaches, ScisTree works with genotype probabilities of individual genotypes (which can be computed by existing single cell genotype callers). ScisTree assumes the infinite sites model. Given uncertain genotypes with individualized probabilities, ScisTree implements a fast heuristic for inferring cell lineage tree and calling the genotypes that allow the so-called perfect phylogeny and maximize the likelihood of the genotypes. Through simulation, we show that ScisTree performs well on the accuracy of inferred trees, and is much more efficient than existing methods. The efficiency of ScisTree enables new applications including imputation of the so-called doublets.

AVAILABILITY AND IMPLEMENTATION

The program ScisTree is available for download at: https://github.com/yufengwudcs/ScisTree.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pancreatic circulating tumor cell detection by targeted single-cell next-generation sequencing

BACKGROUND AND AIMS

Single-cell next-generation sequencing (scNGS) technology has been widely used in genomic profiling, which relies on whole-genome amplification (WGA). However, WGA introduces errors and is especially less accurate when applied to single nucleotide variant (SNV) analysis. Targeted scNGS for SNV without WGA has not been described. We aimed to develop a method to detect circulating tumor cells (CTCs) with DNA SNVs.

METHODS

We tested this targeted scNGS method with three driver mutant genes (KRAS/TP53/SMAD4) on one pancreatic cancer cell line AsPC-1 and then applied it to patients with metastatic PDAC for the validation.

RESULTS

All single-cell of AsPC-1 and spiked-in AsPC-1 cells in healthy donor blood, which were isolated by the filtration with size or by flow cytometry, were detected by targeted scNGS method. All blood samples from six patients with metastatic PDAC, for the validation of target scNGS method, showed CTCs with SNVs of KRAS/TP53/SMAD4 and the positive confirmation of immunofluorescent stainings with Pan-CK/Vimentin/CD45. Four patients with early stage disease, one patient with benign pancreatic cyst and a healthy control sample all showed concordant results between targeted scNGS and CTC enumeration.

CONCLUSIONS

The novel technique of targeted scNGS for SNV analysis, without pre-amplification, is a promising method for identifying and characterizing circulating tumor cells.

Stem cells in cancer initiation and progression

Bajaj et al. review how cancers originate, how their heterogeneity is linked to cancer stem cells, and the signals fundamental to driving these processes.

While standard therapies can lead to an initial remission of aggressive cancers, they are often only a transient solution. The resistance and relapse that follows is driven by tumor heterogeneity and therapy-resistant populations that can reinitiate growth and promote disease progression. There is thus a significant need to understand the cell types and signaling pathways that not only contribute to cancer initiation, but also those that confer resistance and drive recurrence. Here, we discuss work showing that stem cells and progenitors may preferentially serve as a cell of origin for cancers, and that cancer stem cells can be key in driving the continued growth and  functional heterogeneity of established cancers. We also describe emerging evidence for the role of developmental signals in cancer initiation, propagation, and therapy resistance and discuss how targeting these pathways may be of therapeutic value.

Unraveling the Complexity of the Cancer Microenvironment With Multidimensional Genomic and Cytometric Technologies

Cancers are characterized by extensive heterogeneity that occurs intratumorally, between lesions, and across patients. To study cancer as a complex biological system, multidimensional analyses of the tumor microenvironment are paramount. Single-cell technologies such as flow cytometry, mass cytometry, or single-cell RNA-sequencing have revolutionized our ability to characterize individual cells in great detail and, with that, shed light on the complexity of cancer microenvironments. However, a key limitation of these single-cell technologies is the lack of information on spatial context and multicellular interactions. Investigating spatial contexts of cells requires the incorporation of tissue-based techniques such as multiparameter immunofluorescence, imaging mass cytometry, or in situ detection of transcripts. In this Review, we describe the rise of multidimensional single-cell technologies and provide an overview of their strengths and weaknesses. In addition, we discuss the integration of transcriptomic, genomic, epigenomic, proteomic, and spatially-resolved data in the context of human cancers. Lastly, we will deliberate on how the integration of multi-omics data will help to shed light on the complex role of cell types present within the human tumor microenvironment, and how such system-wide approaches may pave the way toward more effective therapies for the treatment of cancer.

Parallel bimodal single-cell sequencing of transcriptome and chromatin accessibility

Joint profiling of transcriptome and chromatin accessibility within single cells allows for the deconstruction of the complex relationship between transcriptional states and upstream regulatory programs determining different cell fates. Here, we developed an automated method with high sensitivity, assay for single-cell transcriptome and accessibility regions (ASTAR-seq), for simultaneous measurement of whole-cell transcriptome and chromatin accessibility within the same single cell. To show the utility of ASTAR-seq, we profiled 384 mESCs under naive and primed pluripotent states as well as a two-cell like state, 424 human cells of various lineage origins (BJ, K562, JK1, and Jurkat), and 480 primary cord blood cells undergoing erythroblast differentiation. With the joint profiles, we configured the transcriptional and chromatin accessibility landscapes of discrete cell states, uncovered linked sets of cis-regulatory elements and target genes unique to each state, and constructed interactome and transcription factor (TF)–centered upstream regulatory networks for various cell states.

Biological and Medical Importance of Cellular Heterogeneity Deciphered by Single-Cell RNA Sequencing

The present review discusses recent progress in single-cell RNA sequencing (scRNA-seq), which can describe cellular heterogeneity in various organs, bodily fluids, and pathologies (e.g., cancer and Alzheimer's disease). We outline scRNA-seq techniques that are suitable for investigating cellular heterogeneity that is present in cell populations with very high resolution of the transcriptomic landscape. We summarize scRNA-seq findings and applications of this technology to identify cell types, activity, and other features that are important for the function of different bodily organs. We discuss future directions for scRNA-seq techniques that can link gene expression, protein expression, cellular function, and their roles in pathology. We speculate on how the field could develop beyond its present limitations (e.g., performing scRNA-seq in situ and in vivo). Finally, we discuss the integration of machine learning and artificial intelligence with cutting-edge scRNA-seq technology, which could provide a strong basis for designing precision medicine and targeted therapy in the future.

Unraveling Heterogeneity in Transcriptome and Its Regulation Through Single-Cell Multi-Omics Technologies

Cellular heterogeneity plays a pivotal role in tissue homeostasis and the disease development of multicellular organisms. To deconstruct the heterogeneity, a multitude of single-cell toolkits measuring various cellular contents, including genome, transcriptome, epigenome, and proteome, have been developed. More recently, multi-omics single-cell techniques enable the capture of molecular footprints with a higher resolution by simultaneously profiling various cellular contents within an individual cell. Integrative analysis of multi-omics datasets unravels the relationships between cellular modalities, builds sophisticated regulatory networks, and provides a holistic view of the cell state. In this review, we summarize the major developments in the single-cell field and review the current state-of-the-art single-cell multi-omic techniques and the bioinformatic tools for integrative analysis.

Advances in histone deacetylase inhibitors in targeting glioblastoma stem cells

Glioblastoma multiforme (GBM) is a lethal grade IV glioma (WHO classification) and widely prevalent primary brain tumor in adults. GBM tumors harbor cellular heterogeneity with the presence of a small subpopulation of tumor cells, described as GBM cancer stem cells (CSCs) that pose resistance to standard anticancer regimens and eventually mediate aggressive relapse or intractable progressive GBM. Existing conventional anticancer therapies for GBM do not target GBM stem cells and are mostly palliative; therefore, exploration of new strategies to target stem cells of GBM has to be prioritized for the development of effective GBM therapy. Recent developments in the understanding of GBM pathophysiology demonstrated dysregulation of epigenetic mechanisms along with the genetic changes in GBM CSCs. Altered expression/activity of key epigenetic regulators, especially histone deacetylases (HDACs) in GBM stem cells has been associated with poor prognosis; inhibiting the activity of HDACs using histone deacetylase inhibitors (HDACi) has been promising as mono-therapeutic in targeting GBM and in sensitizing GBM stem cells to an existing anticancer regimen. Here, we review the development of pan/selective HDACi as potential anticancer agents in targeting the stem cells of glioblastoma as a mono or combination therapy.

Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.