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]

Advances in single-cell DNA sequencing enable insights into human somatic mosaicism

DNA sequencing from bulk or clonal human tissues has shown that genetic mosaicism is common and contributes to both cancer and non-cancerous disorders. However, single-cell resolution is required to understand the full genetic heterogeneity that exists within a tissue and the mechanisms that lead to somatic mosaicism. Single-cell DNA-sequencing technologies have traditionally trailed behind those of single-cell transcriptomics and epigenomics, largely because most applications require whole-genome amplification before costly whole-genome sequencing. Now, recent technological and computational advances are enabling the use of single-cell DNA sequencing to tackle previously intractable problems, such as delineating the genetic landscape of tissues with complex clonal patterns, of samples where cellular material is scarce and of non-cycling, postmitotic cells. Single-cell genomes are also revealing the mutational patterns that arise from biological processes or disease states, and have made it possible to track cell lineage in human tissues. These advances in our understanding of tissue biology and our ability to identify disease mechanisms will ultimately transform how disease is diagnosed and monitored.

Measuring Longitudinal Genome-wide Clonal Evolution of Pediatric Acute Lymphoblastic Leukemia at Single-Cell Resolution

Over 80% of children with acute lymphoblastic leukemia (pALL) can be cured by treating them with multiple chemotherapeutic agents administered over several years, whereas pALL is incurable with 1-3 medications, suggesting significant variation in drug susceptibility across clonal populations. While bulk sequencing studies indicate that pALL cells contain relatively few genetic variants compared to other cancers, the true extent of genetic diversity at the single-cell level remains unknown. Here, we used three complementary approaches to investigate pALL genetic heterogeneity: error-corrected bulk sequencing, single-cell exome sequencing, and primary template-directed amplification (PTA)-enabled single-cell genome sequencing. We discovered that some ETV6-RUNX1 samples harbor multiple independent ras clones and that individual pALL cells harbor substantially more mutations (mean 3,553 per cell) than detected in bulk samples (mean 965 mutations), with variant signatures suggesting both early and late APOBEC-driven mutagenesis in ETV6-RUNX1 patients. Using PTA-based phylogenetic analysis of over 150 single-cell genomes from four pALL patients, we identified heritable phenotypes associated with specific genetic alterations, including some low-frequency clones that are preferentially selected for during chemotherapy treatment. Our findings reveal previously undetected genetic diversity in pALL and suggest that pre-existing mutations influence treatment response, with implications for future therapeutic strategies. This study provides a high-resolution framework for understanding cancer clonal evolution during treatment, yielding important new insights for developing more effective therapeutic approaches for pALL.

Recent Innovations and Technical Advances in High-Throughput Parallel Single-Cell Whole-Genome Sequencing Methods

Single-cell whole-genome sequencing (scWGS) detects cell heterogeneity at the aspect of genomic variations, which are inheritable and play an important role in life processes such as aging and cancer progression. The recent explosive development of high-throughput single-cell sequencing methods has enabled high-performance heterogeneity detection through a vast number of novel strategies. Despite the limitation on total cost, technical advances in high-throughput single-cell whole-genome sequencing methods are made for higher genome coverage, parallel throughput, and level of integration. This review highlights the technical advancements in high-throughput scWGS in the aspects of strategies design, data efficiency, parallel handling platforms, and their applications on human genome. The experimental innovations, remaining challenges, and perspectives are summarized and discussed.

Research progress and the prospect of using single-cell sequencing technology to explore the characteristics of the tumor microenvironment

In precision cancer therapy, addressing intra-tumor heterogeneity poses a significant obstacle. Due to the heterogeneity of each cell subtype and between cells within the tumor, the sensitivity and resistance of different patients to targeted drugs, chemotherapy, etc., are inconsistent. Concerning a specific tumor type, many feasible treatments or combinations can be used by specifically targeting the tumor microenvironment. To solve this problem, it is necessary to further study the tumor microenvironment. Single-cell sequencing techniques can dissect distinct tumor cell populations by isolating cells and using statistical computational methods. This technology may assist in the selection of targeted combination therapy, and the obtained cell subset information is crucial for the rational application of targeted therapy. In this review, we summarized the research and application advances of single-cell sequencing technology in the tumor microenvironment, including the most commonly used single-cell genomic and transcriptomic sequencing, and their future development direction was proposed. The application of single-cell sequencing technology has been expanded to include epigenomics, proteomics, metabolomics, and microbiome analysis. The integration of these different omics approaches has significantly advanced the development of single-cell multiomics sequencing technology. This innovative approach holds immense potential for various fields, such as biological research and medical investigations. Finally, we discussed the advantages and disadvantages of using single-cell sequencing to explore the tumor microenvironment.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Gene doping detection in the era of genomics

Recent progress in gene editing has enabled development of gene therapies for many genetic diseases, but also made gene doping an emerging risk in sports and competitions. By delivery of exogenous transgenes into human body, gene doping not only challenges competition fairness but also places health risk on athletes. World Anti-Doping Agency (WADA) has clearly inhibited the use of gene and cell doping in sports, and many techniques have been developed for gene doping detection. In this review, we will summarize the main tools for gene doping detection at present, highlight the main challenges for current tools, and elaborate future utilizations of high-throughput sequencing for unbiased, sensitive, economic and large-scale gene doping detections. Quantitative real-time PCR assays are the widely used detection methods at present, which are useful for detection of known targets but are vulnerable to codon optimization at exon-exon junction sites of the transgenes. High-throughput sequencing has become a powerful tool for various applications in life and health research, and the era of genomics has made it possible for sensitive and large-scale gene doping detections. Non-biased genomic profiling could efficiently detect new doping targets, and low-input genomics amplification and long-read third-generation sequencing also have application potentials for more efficient and straightforward gene doping detection. By closely monitoring scientific advancements in gene editing and sport genetics, high-throughput sequencing could play a more and more important role in gene detection and hopefully contribute to doping-free sports in the future.

Unlocking the Complexity: Exploration of Acute Lymphoblastic Leukemia at the Single Cell Level

Acute lymphoblastic leukemia (ALL) is the most common cancer in children. ALL originates from precursor lymphocytes that acquire multiple genomic changes over time, including chromosomal rearrangements and point mutations. While a large variety of genomic defects was identified and characterized in ALL over the past 30 years, it was only in recent years that the clonal heterogeneity was recognized. Thanks to the latest advancements in single-cell sequencing techniques, which have evolved from the analysis of a few hundred cells to the analysis of thousands of cells simultaneously, the study of tumor heterogeneity now becomes possible. Different modalities can be explored at the single-cell level: DNA, RNA, epigenetic modifications, and intracellular and cell surface proteins. In this review, we describe these techniques and highlight their advantages and limitations in the study of ALL biology. Moreover, multiomics technologies and the incorporation of the spatial dimension can provide insight into intercellular communication. We describe how the different single-cell sequencing technologies help to unravel the molecular complexity of ALL, shedding light on its development, its heterogeneity, its interaction with the leukemia microenvironment and possible relapse mechanisms.

Unveiling aging dynamics in the hematopoietic system insights from single-cell technologies

As the demographic structure shifts towards an aging society, strategies aimed at slowing down or reversing the aging process become increasingly essential. Aging is a major predisposing factor for many chronic diseases in humans. The hematopoietic system, comprising blood cells and their associated bone marrow microenvironment, intricately participates in hematopoiesis, coagulation, immune regulation and other physiological phenomena. The aging process triggers various alterations within the hematopoietic system, serving as a spectrum of risk factors for hematopoietic disorders, including clonal hematopoiesis, immune senescence, myeloproliferative neoplasms and leukemia. The emerging single-cell technologies provide novel insights into age-related changes in the hematopoietic system. In this review, we summarize recent studies dissecting hematopoietic system aging using single-cell technologies. We discuss cellular changes occurring during aging in the hematopoietic system at the levels of the genomics, transcriptomics, epigenomics, proteomics, metabolomics and spatial multi-omics. Finally, we contemplate the future prospects of single-cell technologies, emphasizing the impact they may bring to the field of hematopoietic system aging research.

Transient Differentiation-State Plasticity Occurs during Acute Lymphoblastic Leukemia Initiation

Leukemia is characterized by oncogenic lesions that result in a block of differentiation, whereas phenotypic plasticity is retained. A better understanding of how these two phenomena arise during leukemogenesis in humans could help inform diagnosis and treatment strategies. Here, we leveraged the well-defined differentiation states during T-cell development to pinpoint the initiation of T-cell acute lymphoblastic leukemia (T-ALL), an aggressive form of childhood leukemia, and study the emergence of phenotypic plasticity. Single-cell whole genome sequencing of leukemic blasts was combined with multiparameter flow cytometry to couple cell identity and clonal lineages. Irrespective of genetic events, leukemia-initiating cells altered their phenotypes by differentiation and dedifferentiation. The construction of the phylogenies of individual leukemias using somatic mutations revealed that phenotypic diversity is reflected by the clonal structure of cancer. The analysis also indicated that the acquired phenotypes are heritable and stable. Together, these results demonstrate a transient period of plasticity during leukemia initiation, where phenotypic switches seem unidirectional. Significance: A method merging multicolor flow cytometry with single-cell whole genome sequencing to couple cell identity with clonal lineages uncovers differentiation-state plasticity in leukemia, reconciling blocked differentiation with phenotypic plasticity in cancer.

Single-cell lineage tracing approaches to track kidney cell development and maintenance

The kidney is a complex organ consisting of various cell types. Previous studies have aimed to elucidate the cellular relationships among these cell types in developing and mature kidneys using Cre-loxP-based lineage tracing. However, this methodology falls short of fully capturing the heterogeneous nature of the kidney, making it less than ideal for comprehensively tracing cellular progression during kidney development and maintenance. Recent technological advancements in single-cell genomics have revolutionized lineage tracing methods. Single-cell lineage tracing enables the simultaneous tracing of multiple cell types within complex tissues and their transcriptomic profiles, thereby allowing the reconstruction of their lineage tree with cell state information. Although single-cell lineage tracing has been successfully applied to investigate cellular hierarchies in various organs and tissues, its application in kidney research is currently lacking. This review comprehensively consolidates the single-cell lineage tracing methods, divided into 4 categories (clustered regularly interspaced short palindromic repeat [CRISPR]/CRISPR-associated protein 9 [Cas9]-based, transposon-based, Polylox-based, and native barcoding methods), and outlines their technical advantages and disadvantages. Furthermore, we propose potential future research topics in kidney research that could benefit from single-cell lineage tracing and suggest suitable technical strategies to apply to these topics.

Heterogeneity and molecular landscape of melanoma: implications for targeted therapy

Uveal cancer (UM) offers a complex molecular landscape characterized by substantial heterogeneity, both on the genetic and epigenetic levels. This heterogeneity plays a critical position in shaping the behavior and response to therapy for this uncommon ocular malignancy. Targeted treatments with gene-specific therapeutic molecules may prove useful in overcoming radiation resistance, however, the diverse molecular makeups of UM call for a patient-specific approach in therapy procedures. We need to understand the intricate molecular landscape of UM to develop targeted treatments customized to each patient's specific genetic mutations. One of the promising approaches is using liquid biopsies, such as circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA), for detecting and monitoring the disease at the early stages. These non-invasive methods can help us identify the most effective treatment strategies for each patient. Single-cellular is a brand-new analysis platform that gives treasured insights into diagnosis, prognosis, and remedy. The incorporation of this data with known clinical and genomics information will give a better understanding of the complicated molecular mechanisms that UM diseases exploit. In this review, we focused on the heterogeneity and molecular panorama of UM, and to achieve this goal, the authors conducted an exhaustive literature evaluation spanning 1998 to 2023, using keywords like "uveal melanoma, "heterogeneity". "Targeted therapies"," "CTCs," and "single-cellular analysis".

Comparison of high-throughput single-cell RNA-seq methods for ex vivo drug screening

Functional precision medicine (FPM) aims to optimize patient-specific drug selection based on the unique characteristics of their cancer cells. Recent advancements in high throughput ex vivo drug profiling have accelerated interest in FPM. Here, we present a proof-of-concept study for an integrated experimental system that incorporates ex vivo treatment response with a single-cell gene expression output enabling barcoding of several drug conditions in one single-cell sequencing experiment. We demonstrate this through a proof-of-concept investigation focusing on the glucocorticoid-resistant acute lymphoblastic leukemia (ALL) E/R+ Reh cell line. Three different single-cell transcriptome sequencing (scRNA-seq) approaches were evaluated, each exhibiting high cell recovery and accurate tagging of distinct drug conditions. Notably, our comprehensive analysis revealed variations in library complexity, sensitivity (gene detection), and differential gene expression detection across the methods. Despite these differences, we identified a substantial transcriptional response to fludarabine, a highly relevant drug for treating high-risk ALL, which was consistently recapitulated by all three methods. These findings highlight the potential of our integrated approach for studying drug responses at the single-cell level and emphasize the importance of method selection in scRNA-seq studies. Finally, our data encompassing 27 327 cells are freely available to extend to future scRNA-seq methodological comparisons.

Single Nucleotide Polymorphism (SNP) and Antibody-based Cell Sorting (SNACS): A tool for demultiplexing single-cell DNA sequencing data

Motivation

Recently, single-cell DNA sequencing (scDNA-seq) and multi-modal profiling with the addition of cell-surface antibodies (scDAb-seq) have provided key insights into cancer heterogeneity. Scaling these technologies across large patient cohorts, however, is cost and time prohibitive. Multiplexing, in which cells from unique patients are pooled into a single experiment, offers a possible solution. While multiplexing methods exist for scRNAseq, accurate demultiplexing in scDNAseq remains an unmet need.

Results

Here, we introduce SNACS: Single-Nucleotide Polymorphism (SNP) and Antibody-based Cell Sorting. SNACS relies on a combination of patient-level cell-surface identifiers and natural variation in genetic polymorphisms to demultiplex scDNAseq data. We demonstrated the performance of SNACS on a dataset consisting of multi-sample experiments from patients with leukemia where we knew truth from single-sample experiments from the same patients. Using SNACS, accuracy ranged from 0.948 - 0.991 vs 0.552 - 0.934 using demultiplexing methods from the single-cell literature.

Decoding leukemia at the single-cell level: clonal architecture, classification, microenvironment, and drug resistance

Leukemias are refractory hematological malignancies, characterized by marked intrinsic heterogeneity which poses significant obstacles to effective treatment. However, traditional bulk sequencing techniques have not been able to effectively unravel the heterogeneity among individual tumor cells. With the emergence of single-cell sequencing technology, it has bestowed upon us an unprecedented resolution to comprehend the mechanisms underlying leukemogenesis and drug resistance across various levels, including the genome, epigenome, transcriptome and proteome. Here, we provide an overview of the currently prevalent single-cell sequencing technologies and a detailed summary of single-cell studies conducted on leukemia, with a specific focus on four key aspects: (1) leukemia's clonal architecture, (2) frameworks to determine leukemia subtypes, (3) tumor microenvironment (TME) and (4) the drug-resistant mechanisms of leukemia. This review provides a comprehensive summary of current single-cell studies on leukemia and highlights the markers and mechanisms that show promising clinical implications for the diagnosis and treatment of leukemia.

Computational immunogenomic approaches to predict response to cancer immunotherapies

Cancer immunogenomics is an emerging field that bridges genomics and immunology. The establishment of large-scale genomic collaborative efforts along with the development of new single-cell transcriptomic techniques and multi-omics approaches have enabled characterization of the mutational and transcriptional profiles of many cancer types and helped to identify clinically actionable alterations as well as predictive and prognostic biomarkers. Researchers have developed computational approaches and machine learning algorithms to accurately obtain clinically useful information from genomic and transcriptomic sequencing data from bulk tissue or single cells and explore tumours and their microenvironment. The rapid growth in sequencing and computational approaches has resulted in the unmet need to understand their true potential and limitations in enabling improvements in the management of patients with cancer who are receiving immunotherapies. In this Review, we describe the computational approaches currently available to analyse bulk tissue and single-cell sequencing data from cancer, stromal and immune cells, as well as how best to select the most appropriate tool to address various clinical questions and, ultimately, improve patient outcomes.

Inferring Markov Chains to Describe Convergent Tumor Evolution With CIMICE

The field of tumor phylogenetics focuses on studying the differences within cancer cell populations. Many efforts are done within the scientific community to build cancer progression models trying to understand the heterogeneity of such diseases. These models are highly dependent on the kind of data used for their construction, therefore, as the experimental technologies evolve, it is of major importance to exploit their peculiarities. In this work we describe a cancer progression model based on Single Cell DNA Sequencing data. When constructing the model, we focus on tailoring the formalism on the specificity of the data. We operate by defining a minimal set of assumptions needed to reconstruct a flexible DAG structured model, capable of identifying progression beyond the limitation of the infinite site assumption. Our proposal is conservative in the sense that we aim to neither discard nor infer knowledge which is not represented in the data. We provide simulations and analytical results to show the features of our model, test it on real data, show how it can be integrated with other approaches to cope with input noise. Moreover, our framework can be exploited to produce simulated data that follows our theoretical assumptions. Finally, we provide an open source R implementation of our approach, called CIMICE, that is publicly available on BioConductor.

Intratumor heterogeneity: models of malignancy emergence and evolution

Cancer is a complex and heterogeneous disease characterized by the accumulation of genetic alterations that drive uncontrolled cell growth and proliferation. Evolutionary dynamics plays a crucial role in the emergence and development of tumors, shaping the heterogeneity and adaptability of cancer cells. From the perspective of evolutionary theory, tumors are complex ecosystems that evolve through a process of microevolution influenced by genetic mutations, epigenetic changes, tumor microenvironment factors, and therapy-induced changes. This dynamic nature of tumors poses significant challenges for effective cancer treatment, and understanding it is essential for developing effective and personalized therapies. By uncovering the mechanisms that determine tumor heterogeneity, researchers can identify key genetic and epigenetic changes that contribute to tumor progression and resistance to treatment. This knowledge enables the development of innovative strategies for targeting specific tumor clones, minimizing the risk of recurrence and improving patient outcomes. To investigate the evolutionary dynamics of cancer, researchers employ a wide range of experimental and computational approaches. Traditional experimental methods involve genomic profiling techniques such as next-generation sequencing and fluorescence in situ hybridization. These techniques enable the identification of somatic mutations, copy number alterations, and structural rearrangements within cancer genomes. Furthermore, single-cell sequencing methods have emerged as powerful tools for dissecting intratumoral heterogeneity and tracing clonal evolution. In parallel, computational models and algorithms have been developed to simulate and analyze cancer evolution. These models integrate data from multiple sources to predict tumor growth patterns, identify driver mutations, and infer evolutionary trajectories. In this paper, we set out to describe the current approaches to address this evolutionary complexity and theories of its occurrence.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.

Understanding blood development and leukemia using sequencing-based technologies and human cell systems

Our current understanding of human hematopoiesis has undergone significant transformation throughout the years, challenging conventional views. The evolution of high-throughput technologies has enabled the accumulation of diverse data types, offering new avenues for investigating key regulatory processes in blood cell production and disease. In this review, we will explore the opportunities presented by these advancements for unraveling the molecular mechanisms underlying normal and abnormal hematopoiesis. Specifically, we will focus on the importance of enhancer-associated regulatory networks and highlight the crucial role of enhancer-derived transcription regulation. Additionally, we will discuss the unprecedented power of single-cell methods and the progression in using in vitro human blood differentiation system, in particular induced pluripotent stem cell models, in dissecting hematopoietic processes. Furthermore, we will explore the potential of ever more nuanced patient profiling to allow precision medicine approaches. Ultimately, we advocate for a multiparameter, regulatory network-based approach for providing a more holistic understanding of normal hematopoiesis and blood disorders.

Assessing the performance of methods for cell clustering from single-cell DNA sequencing data

BACKGROUND

Many cancer genomes have been known to contain more than one subclone inside one tumor, the phenomenon of which is called intra-tumor heterogeneity (ITH). Characterizing ITH is essential in designing treatment plans, prognosis as well as the study of cancer progression. Single-cell DNA sequencing (scDNAseq) has been proven effective in deciphering ITH. Cells corresponding to each subclone are supposed to carry a unique set of mutations such as single nucleotide variations (SNV). While there have been many studies on the cancer evolutionary tree reconstruction, not many have been proposed that simply characterize the subclonality without tree reconstruction. While tree reconstruction is important in the study of cancer evolutionary history, typically they are computationally expensive in terms of running time and memory consumption due to the huge search space of the tree structure. On the other hand, subclonality characterization of single cells can be converted into a cell clustering problem, the dimension of which is much smaller, and the turnaround time is much shorter. Despite the existence of a few state-of-the-art cell clustering computational tools for scDNAseq, there lacks a comprehensive and objective comparison under different settings.

RESULTS

In this paper, we evaluated six state-of-the-art cell clustering tools-SCG, BnpC, SCClone, RobustClone, SCITE and SBMClone-on simulated data sets given a variety of parameter settings and a real data set. We designed a simulator specifically for cell clustering, and compared these methods' performances in terms of their clustering accuracy, specificity and sensitivity and running time. For SBMClone, we specifically designed an ultra-low coverage large data set to evaluate its performance in the face of an extremely high missing rate.

CONCLUSION

From the benchmark study, we conclude that BnpC and SCG's clustering accuracy are the highest and comparable to each other. However, BnpC is more advantageous in terms of running time when cell number is high (> 1500). It also has a higher clustering accuracy than SCG when cluster number is high (> 16). SCClone's accuracy in estimating the number of clusters is the highest. RobustClone and SCITE's clustering accuracy are the lowest for all experiments. SCITE tends to over-estimate the cluster number and has a low specificity, whereas RobustClone tends to under-estimate the cluster number and has a much lower sensitivity than other methods. SBMClone produced reasonably good clustering (V-measure > 0.9) when coverage is > = 0.03 and thus is highly recommended for ultra-low coverage large scDNAseq data sets.

Novel precision medicine approaches and treatment strategies in hematological malignancies

Genetic testing has been applied for decades in clinical routine diagnostics of hematological malignancies to improve disease (sub)classification, prognostication, patient management, and survival. In recent classifications of hematological malignancies, disease subtypes are defined by key recurrent genetic alterations detected by conventional methods (i.e., cytogenetics, fluorescence in situ hybridization, and targeted sequencing). Hematological malignancies were also one of the first disease areas in which targeted therapies were introduced, the prime example being BCR::ABL1 inhibitors, followed by an increasing number of targeted inhibitors hitting the Achilles' heel of each disease, resulting in a clear patient benefit. Owing to the technical advances in high-throughput sequencing, we can now apply broad genomic tests, including comprehensive gene panels or whole-genome and whole-transcriptome sequencing, to identify clinically important diagnostic, prognostic, and predictive markers. In this review, we give examples of how precision diagnostics has been implemented to guide treatment selection and improve survival in myeloid (myelodysplastic syndromes and acute myeloid leukemia) and lymphoid malignancies (acute lymphoblastic leukemia, diffuse large B-cell lymphoma, and chronic lymphocytic leukemia). We discuss the relevance and potential of monitoring measurable residual disease using ultra-sensitive techniques to assess therapy response and detect early relapses. Finally, we bring up the promising avenue of functional precision medicine, combining ex vivo drug screening with various omics technologies, to provide novel treatment options for patients with advanced disease. Although we are only in the beginning of the field of precision hematology, we foresee rapid development with new types of diagnostics and treatment strategies becoming available to the benefit of our patients.

Barcoded multiple displacement amplification for high coverage sequencing in spatial genomics

Determining mutational landscapes in a spatial context is essential for understanding genetically heterogeneous cell microniches. Current approaches, such as Multiple Displacement Amplification (MDA), offer high genome coverage but limited multiplexing, which hinders large-scale spatial genomic studies. Here, we introduce barcoded MDA (bMDA), a technique that achieves high-coverage genomic analysis of low-input DNA while enhancing the multiplexing capabilities. By incorporating cell barcodes during MDA, bMDA streamlines library preparation in one pot, thereby overcoming a key bottleneck in spatial genomics. We apply bMDA to the integrative spatial analysis of triple-negative breast cancer tissues by examining copy number alterations, single nucleotide variations, structural variations, and kataegis signatures for each spatial microniche. This enables the assessment of subclonal evolutionary relationships within a spatial context. Therefore, bMDA has emerged as a scalable technology with the potential to advance the field of spatial genomics significantly.

A weighted distance-based approach for deriving consensus tumor evolutionary trees

MOTIVATION

The acquisition of somatic mutations by a tumor can be modeled by a type of evolutionary tree. However, it is impossible to observe this tree directly. Instead, numerous algorithms have been developed to infer such a tree from different types of sequencing data. But such methods can produce conflicting trees for the same patient, making it desirable to have approaches that can combine several such tumor trees into a consensus or summary tree. We introduce The Weighted m-Tumor Tree Consensus Problem (W-m-TTCP) to find a consensus tree among multiple plausible tumor evolutionary histories, each assigned a confidence weight, given a specific distance measure between tumor trees. We present an algorithm called TuELiP that is based on integer linear programming which solves the W-m-TTCP, and unlike other existing consensus methods, allows the input trees to be weighted differently.

RESULTS

On simulated data we show that TuELiP outperforms two existing methods at correctly identifying the true underlying tree used to create the simulations. We also show that the incorporation of weights can lead to more accurate tree inference. On a Triple-Negative Breast Cancer dataset, we show that including confidence weights can have important impacts on the consensus tree identified.

AVAILABILITY

An implementation of TuELiP and simulated datasets are available at https://bitbucket.org/oesperlab/consensus-ilp/src/main/.

Single-Cell Next-Generation Sequencing to Monitor Hematopoietic Stem-Cell Transplantation: Current Applications and Future Perspectives

Acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS) are genetically complex and diverse diseases. Such complexity makes challenging the monitoring of response to treatment. Measurable residual disease (MRD) assessment is a powerful tool for monitoring response and guiding therapeutic interventions. This is accomplished through targeted next-generation sequencing (NGS), as well as polymerase chain reaction and multiparameter flow cytometry, to detect genomic aberrations at a previously challenging leukemic cell concentration. A major shortcoming of NGS techniques is the inability to discriminate nonleukemic clonal hematopoiesis. In addition, risk assessment and prognostication become more complicated after hematopoietic stem-cell transplantation (HSCT) due to genotypic drift. To address this, newer sequencing techniques have been developed, leading to more prospective and randomized clinical trials aiming to demonstrate the prognostic utility of single-cell next-generation sequencing in predicting patient outcomes following HSCT. This review discusses the use of single-cell DNA genomics in MRD assessment for AML/MDS, with an emphasis on the HSCT time period, including the challenges with current technologies. We also touch on the potential benefits of single-cell RNA sequencing and analysis of accessible chromatin, which generate high-dimensional data at the cellular resolution for investigational purposes, but not currently used in the clinical setting.

RobustTree: An adaptive, robust PCA algorithm for embedded tree structure recovery from single-cell sequencing data

Robust Principal Component Analysis (RPCA) offers a powerful tool for recovering a low-rank matrix from highly corrupted data, with growing applications in computational biology. Biological processes commonly form intrinsic hierarchical structures, such as tree structures of cell development trajectories and tumor evolutionary history. The rapid development of single-cell sequencing (SCS) technology calls for the recovery of embedded tree structures from noisy and heterogeneous SCS data. In this study, we propose RobustTree, a unified framework to reconstruct the inherent topological structure underlying high-dimensional data with noise. By extending RPCA to handle tree structure optimization, RobustTree leverages data denoising, clustering, and tree structure reconstruction. It solves the tree optimization problem with an adaptive parameter selection scheme that we proposed. In addition to recovering real datasets, RobustTree can reconstruct continuous topological structure and discrete-state topological structure of underlying SCS data. We apply RobustTree on multiple synthetic and real datasets and demonstrate its high accuracy and robustness when analyzing high-noise SCS data with embedded complex structures. The code is available at https://github.com/ucasdp/RobustTree.

Application of high-throughput single-nucleus DNA sequencing in pancreatic cancer

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high-throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimize a highly automated nuclei extraction workflow that achieves fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic ductal adenocarcinoma patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals genomic evolution patterns of pancreatic ductal adenocarcinoma as well.

Application of single-cell transcriptome sequencing in gastric cancer

Gastric cancer (GC) is the fifth most common cancer worldwide and the third leading cause of cancer death. With the development of single-cell RNA-sequencing (scRNA-seq) technology, the research on GC has gradually developed from the histopathological level to the transcriptional level. In this paper, we discuss the principle of scRNA-seq technology and its application in GC research, including the transcriptional characteristics and origin of GC precancerous lesions, intratumor heterogeneity of primary tumors, tumor microenvironment, and metastatic dissemination.

Single-cell analysis of acute lymphoblastic and lineage-ambiguous leukemia: approaches and molecular insights

Despite recent progress in identifying the genetic drivers of acute lymphoblastic leukemia (ALL), prognosis remains poor for those individuals who experience disease recurrence. Moreover, acute leukemias of ambiguous lineage lack a biologically informed framework to guide classification and therapy. These needs have driven the adoption of multiple complementary single-cell sequencing approaches to explore key issues in the biology of these leukemias, including cell of origin, developmental hierarchy and ontogeny, and the molecular heterogeneity driving pathogenesis, progression, and therapeutic responsiveness. There are multiple single-cell techniques for profiling a specific modality, including RNA, DNA, chromatin accessibility and methylation; and an expanding range of approaches for simultaneous analysis of multiple modalities. Single-cell sequencing approaches have also enabled characterization of cell-intrinsic and -extrinsic features of ALL biology. In this review we describe these approaches and highlight the extensive heterogeneity that underpins ALL gene expression, cellular differentiation, and clonal architecture throughout disease pathogenesis and treatment resistance. In addition, we discuss the importance of the dynamic interactions that occur between leukemia cells and the nonleukemia microenvironment. We discuss potential opportunities and limitations of single-cell sequencing for the study of ALL biology and treatment responsiveness.

Reconstructing mutational lineages in breast cancer by multi-patient-targeted single-cell DNA sequencing

Single-cell DNA sequencing (scDNA-seq) methods are powerful tools for profiling mutations in cancer cells; however, most genomic regions sequenced in single cells are non-informative. To overcome this issue, we developed a multi-patient-targeted (MPT) scDNA-seq method. MPT involves first performing bulk exome sequencing across a cohort of cancer patients to identify somatic mutations, which are then pooled together to develop a single custom targeted panel for high-throughput scDNA-seq using a microfluidics platform. We applied MPT to profile 330 mutations across 23,500 cells from 5 patients with triple negative-breast cancer (TNBC), which showed that 3 tumors were monoclonal and 2 tumors were polyclonal. From these data, we reconstructed mutational lineages and identified early mutational and copy-number events, including early TP53 mutations that occurred in all five patients. Collectively, our data suggest that MPT can overcome a major technical obstacle for studying tumor evolution using scDNA-seq by profiling information-rich mutation sites.

Single-cell technologies uncover intra-tumor heterogeneity in childhood cancers

Childhood cancer is the second leading cause of death in children aged 1 to 14. Although survival rates have vastly improved over the past 40 years, cancer resistance and relapse remain a significant challenge. Advances in single-cell technologies enable dissection of tumors to unprecedented resolution. This facilitates unraveling the heterogeneity of childhood cancers to identify cell subtypes that are prone to treatment resistance. The rapid accumulation of single-cell data from different modalities necessitates the development of novel computational approaches for processing, visualizing, and analyzing single-cell data. Here, we review single-cell approaches utilized or under development in the context of childhood cancers. We review computational methods for analyzing single-cell data and discuss best practices for their application. Finally, we review the impact of several studies of childhood tumors analyzed with these approaches and future directions to implement single-cell studies into translational cancer research in pediatric oncology.

bmVAE: a variational autoencoder method for clustering single-cell mutation data

MOTIVATION

Genetic intra-tumor heterogeneity (ITH) characterizes the differences in genomic variations between tumor clones, and accurately unmasking ITH is important for personalized cancer therapy. Single-cell DNA sequencing now emerges as a powerful means for deciphering underlying ITH based on point mutations of single cells. However, detecting tumor clones from single-cell mutation data remains challenging due to the error-prone and discrete nature of the data.

RESULTS

We introduce bmVAE, a bioinformatics tool for learning low-dimensional latent representation of single cell based on a variational autoencoder and then clustering cells into subpopulations in the latent space. bmVAE takes single-cell binary mutation data as inputs, and outputs inferred cell subpopulations as well as their genotypes. To achieve this, the bmVAE framework is designed to consist of three modules including dimensionality reduction, cell clustering and genotype estimation. We assess the method on various synthetic datasets where different factors including false negative rate, data size and data heterogeneity are considered in simulation, and further demonstrate its effectiveness on two real datasets. The results suggest bmVAE is highly effective in reasoning ITH, and performs competitive to existing methods.

AVAILABILITY AND IMPLEMENTATION

bmVAE is freely available at https://github.com/zhyu-lab/bmvae.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Accurate tumor clonal structures require single-cell analysis

Tumor clonal structure is closely related to future progression, which has been mainly investigated as mutation abundance clustering in bulk samples. With relatively limited studies at single-cell resolution, a systematic comparison of the two approaches is still lacking. Here, using bulk and single-cell mutational data from the liver and colorectal cancers, we checked whether co-mutations determined by single-cell analysis had corresponding bulk variant allele frequency (VAF) peaks. While bulk analysis suggested the absence of subclonal peaks and, possibly, neutral evolution in some cases, the single-cell analysis identified coexisting subclones. The overlaps of bulk VAF ranges for co-mutations from different subclones made it difficult to separate them. Complex subclonal structures and dynamic evolution could be hidden under the seemingly clonal neutral pattern at the bulk level, suggesting single-cell analysis is necessary to avoid underestimation of tumor heterogeneity.

Integrating Micro and Nano Technologies for Cell Engineering and Analysis: Toward the Next Generation of Cell Therapy Workflows

The emerging field of cell therapy offers the potential to treat and even cure a diverse array of diseases for which existing interventions are inadequate. Recent advances in micro and nanotechnology have added a multitude of single cell analysis methods to our research repertoire. At the same time, techniques have been developed for the precise engineering and manipulation of cells. Together, these methods have aided the understanding of disease pathophysiology, helped formulate corrective interventions at the cellular level, and expanded the spectrum of available cell therapeutic options. This review discusses how micro and nanotechnology have catalyzed the development of cell sorting, cellular engineering, and single cell analysis technologies, which have become essential workflow components in developing cell-based therapeutics. The review focuses on the technologies adopted in research studies and explores the opportunities and challenges in combining the various elements of cell engineering and single cell analysis into the next generation of integrated and automated platforms that can accelerate preclinical studies and translational research.

Duplex Sequencing Uncovers Recurrent Low-frequency Cancer-associated Mutations in Infant and Childhood KMT2A-rearranged Acute Leukemia

Infant acute lymphoblastic leukemia (ALL) with KMT2A-gene rearrangements (KMT2A-r) have few mutations and a poor prognosis. To uncover mutations that are below the detection of standard next-generation sequencing (NGS), a combination of targeted duplex sequencing and NGS was applied on 20 infants and 7 children with KMT2A-r ALL, 5 longitudinal and 6 paired relapse samples. Of identified nonsynonymous mutations, 87 had been previously implicated in cancer and targeted genes recurrently altered in KMT2A-r leukemia and included mutations in KRAS, NRAS, FLT3, TP53, PIK3CA, PAX5, PIK3R1, and PTPN11, with infants having fewer such mutations. Of identified cancer-associated mutations, 62% were below the resolution of standard NGS. Only 33 of 87 mutations exceeded 2% of cellular prevalence and most-targeted PI3K/RAS genes (31/33) and typically KRAS/NRAS. Five patients only had low-frequency PI3K/RAS mutations without a higher-frequency signaling mutation. Further, drug-resistant clones with FLT3 D835H or NRAS G13D/G12S mutations that comprised only 0.06% to 0.34% of diagnostic cells, expanded at relapse. Finally, in longitudinal samples, the relapse clone persisted as a minor subclone from diagnosis and through treatment before expanding during the last month of disease. Together, we demonstrate that infant and childhood KMT2A-r ALL harbor low-frequency cancer-associated mutations, implying a vast subclonal genetic landscape.

ipaQTL-atlas: an atlas of intronic polyadenylation quantitative trait loci across human tissues

Functional interpretation of disease-associated non-coding variants remains a significant challenge in the post-GWAS era. Our recent study has identified 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs) and connects APA events with QTLs as a major driver of human traits and diseases. Besides 3'UTR, APA events can also occur in intron regions, and increasing evidence has connected intronic polyadenylation with disease risk. However, systematic investigation of the roles of intronic polyadenylation in human diseases remained challenging due to the lack of a comprehensive database across a variety of human tissues. Here, we developed ipaQTL-atlas (http://bioinfo.szbl.ac.cn/ipaQTL) as the first comprehensive portal for intronic polyadenylation. The ipaQTL-atlas is based on the analysis of 15 170 RNA-seq data from 838 individuals across 49 Genotype-Tissue Expression (GTEx v8) tissues and contains ∼0.98 million SNPs associated with intronic APA events. It provides an interface for ipaQTLs search, genome browser, boxplots, and data download, as well as the visualization of GWAS and ipaQTL colocalization results. ipaQTL-atlas provides a one-stop portal to access intronic polyadenylation information and could significantly advance the discovery of APA-associated disease susceptibility genes.

Single-cell mutation calling and phylogenetic tree reconstruction with loss and recurrence

Motivation

Tumours evolve as heterogeneous populations of cells, which may be distinguished by different genomic aberrations. The resulting intra-tumour heterogeneity plays an important role in cancer patient relapse and treatment failure, so that obtaining a clear understanding of each patient’s tumour composition and evolutionary history is key for personalized therapies. Single-cell sequencing (SCS) now provides the possibility to resolve tumour heterogeneity at the highest resolution of individual tumour cells, but brings with it challenges related to the particular noise profiles of the sequencing protocols as well as the complexity of the underlying evolutionary process.

Results

By modelling the noise processes and allowing mutations to be lost or to reoccur during tumour evolution, we present a method to jointly call mutations in each cell, reconstruct the phylogenetic relationship between cells, and determine the locations of mutational losses and recurrences. Our Bayesian approach allows us to accurately call mutations as well as to quantify our certainty in such predictions. We show the advantages of allowing mutational loss or recurrence with simulated data and present its application to tumour SCS data.

Availability and implementation

SCIΦN is available at https://github.com/cbg-ethz/SCIPhIN.

Supplementary information

Supplementary data are available at Bioinformatics online.

Frontiers in single cell analysis: multimodal technologies and their clinical perspectives

Single cell multimodal analysis is at the frontier of single cell research: it defines the roles and functions of distinct cell types through simultaneous analysis to provide unprecedented insight into cellular processes. Current single cell approaches are rapidly moving toward multimodal characterizations. It replaces one-dimensional single cell analysis, for example by allowing for simultaneous measurement of transcription and post-transcriptional regulation, epigenetic modifications and/or surface protein expression. By providing deeper insights into single cell processes, multimodal single cell analyses paves the way to new understandings in various cellular processes such as cell fate decisions, physiological heterogeneity or genotype-phenotype linkages. At the forefront of this, microfluidics is key for high-throughput single cell analysis. Here, we present an overview of the recent multimodal microfluidic platforms having a potential in biomedical research, with a specific focus on their potential clinical applications.

Concordant and Heterogeneity of Single-Cell Transcriptome in Cardiac Development of Human and Mouse

Normal heart development is vital for maintaining its function, and the development process is involved in complex interactions between different cell lineages. How mammalian hearts develop differently is still not fully understood. In this study, we identified several major types of cardiac cells, including cardiomyocytes (CMs), fibroblasts (FBs), endothelial cells (ECs), ECs/FBs, epicardial cells (EPs), and immune cells (macrophage/monocyte cluster, MACs/MONOs), based on single-cell transcriptome data from embryonic hearts of both human and mouse. Then, species-shared and species-specific marker genes were determined in the same cell type between the two species, and the genes with consistent and different expression patterns were also selected by constructing the developmental trajectories. Through a comparison of the development stage similarity of CMs, FBs, and ECs/FBs between humans and mice, it is revealed that CMs at e9.5 and e10.5 of mice are most similar to those of humans at 7 W and 9 W, respectively. Mouse FBs at e10.5, e13.5, and e14.5 are correspondingly more like the same human cells at 6, 7, and 9 W. Moreover, the e9.5-ECs/FBs of mice are most similar to that of humans at 10W. These results provide a resource for understudying cardiac cell types and the crucial markers able to trace developmental trajectories among the species, which is beneficial for finding suitable mouse models to detect human cardiac physiology and related diseases.

The Promise of Single-cell Technology in Providing New Insights Into the Molecular Heterogeneity and Management of Acute Lymphoblastic Leukemia

Drug resistance and treatment failure in pediatric acute lymphoblastic leukemia (ALL) are in part driven by tumor heterogeneity and clonal evolution. Although bulk tumor genomic analyses have provided some insight into these processes, single-cell sequencing has emerged as a powerful technique to profile individual cells in unprecedented detail. Since the introduction of single-cell RNA sequencing, we now have the capability to capture not only transcriptomic, but also genomic, epigenetic, and proteomic variation between single cells separately and in combination. This rapidly evolving field has the potential to transform our understanding of the fundamental biology of pediatric ALL and guide the management of ALL patients to improve their clinical outcome. Here, we discuss the impact single-cell sequencing has had on our understanding of tumor heterogeneity and clonal evolution in ALL and provide examples of how single-cell technology can be integrated into the clinic to inform treatment decisions for children with high-risk disease.

Monitoring of Leukemia Clones in B-cell Acute Lymphoblastic Leukemia at Diagnosis and During Treatment by Single-cell DNA Amplicon Sequencing

Acute lymphoblastic leukemia (ALL) is characterized by the presence of chromosomal changes, including numerical changes, translocations, and deletions, which are often associated with additional single-nucleotide mutations. In this study, we used single cell–targeted DNA sequencing to evaluate the clonal heterogeneity of B-ALL at diagnosis and during chemotherapy treatment. We designed a custom DNA amplicon library targeting mutational hotspot regions (in 110 genes) present in ALL, and we measured the presence of mutations and small insertions/deletions (indels) in bone marrow or blood samples from 12 B-ALL patients, with a median of 7973 cells per sample. Nine of the 12 cases showed at least 1 subclonal mutation, of which cases with PAX5 alterations or high hyperdiploidy (with intermediate to good prognosis) showed a high number of subclones (1 to 7) at diagnosis, defined by a variety of mutations in the JAK/STAT, RAS, or FLT3 signaling pathways. Cases with RAS pathway mutations had multiple mutations in FLT3, NRAS, KRAS, or BRAF in various clones. For those cases where we detected multiple mutational clones at diagnosis, we also studied blood samples during the first weeks of chemotherapy treatment. The leukemia clones disappeared during treatment with various kinetics, and few cells with mutations were easily detectable, even at low frequency (<0.1%). Our data illustrate that about half of the B-ALL cases show >2 subclones at diagnosis and that even very rare mutant cells can be detected at diagnosis or during treatment by single cell–targeted DNA sequencing.

Distinct genomic landscape of Chinese pediatric acute myeloid leukemia impacts clinical risk classification

Studies have revealed key genomic aberrations in pediatric acute myeloid leukemia (AML) based on Western populations. It is unknown to what extent the current genomic findings represent populations with different ethnic backgrounds. Here we present the genomic landscape of driver alterations of Chinese pediatric AML and discover previously undescribed genomic aberrations, including the XPO1-TNRC18 fusion. Comprehensively comparing between the Chinese and Western AML cohorts reveal a substantially distinct genomic alteration profile. For example, Chinese AML patients more commonly exhibit mutations in KIT and CSF3R, and less frequently mutated of genes in the RAS signaling pathway. These differences in mutation frequencies lead to the detection of previously uncharacterized co-occurring mutation pairs. Importantly, the distinct driver profile is clinical relevant. We propose a refined prognosis risk classification model which better reflected the adverse event risk for Chinese AML patients. These results emphasize the importance of genetic background in precision medicine.

ARCANE-ROG: Algorithm for Reconstruction of Cancer Evolution from single-cell data using Robust Graph Learning

Tumor heterogeneity, marked by the presence of divergent clonal subpopulations of tumor cells, impedes the treatment response in cancer patients. Single-cell sequencing technology provides substantial prospects to gain an in-depth understanding of the cellular phenotypic variability driving tumor progression. A comprehensive insight into the intra-tumor heterogeneity may further assist in dealing with the treatment-resistant clones in cancer patients, thereby improving their overall survival. However, this task is hampered due to the challenges associated with the single-cell data, such as false positives, false negatives and missing bases, and the increase in their size. As a result, the computational cost of the existing methods increases, thereby limiting their usage. In this work, we propose a robust graph learning-based method, ARCANE-ROG (Algorithm for Reconstruction of CANcer Evolution via RObust Graph learning), for inferring clonal evolution from single-cell datasets. The first step of the proposed method is a joint framework of denoising with data imputation for the noisy and incomplete matrix while simultaneously learning an adjacency graph. Both the operations in the joint framework boost each other such that the overall performance of the denoising algorithm is improved. In the second step, an optimal number of clusters are identified via the Leiden method. In the last step, clonal evolution trees are inferred via a minimum spanning tree algorithm. The method has been benchmarked against a state-of-the-art method, RobustClone, using simulated datasets of varying sizes and five real datasets. The performance of our proposed method is found to be significantly superior (p-value < 0.05) in terms of reconstruction error, False Positive to False Negative (FPFN) ratio, tree distance error and V-measure compared to the other method. Overall, the proposed method is an improvement over the existing methods as it enhances cluster assignment and inference on clonal hierarchies.

SCClone: Accurate Clustering of Tumor Single-Cell DNA Sequencing Data

Single-cell DNA sequencing (scDNA-seq) enables high-resolution profiling of genetic diversity among single cells and is especially useful for deciphering the intra-tumor heterogeneity and evolutionary history of tumor. Specific technical issues such as allele dropout, false-positive errors, and doublets make scDNA-seq data incomplete and error-prone, giving rise to a severe challenge of accurately inferring clonal architecture of tumor. To effectively address these issues, we introduce a new computational method called SCClone for reasoning subclones from single nucleotide variation (SNV) data of single cells. Specifically, SCClone leverages a probability mixture model for binary data to cluster single cells into distinct subclones. To accurately decipher underlying clonal composition, a novel model selection scheme based on inter-cluster variance is employed to find the optimal number of subclones. Extensive evaluations on various simulated datasets suggest SCClone has strong robustness against different technical noises in scDNA-seq data and achieves better performance than the state-of-the-art methods in reasoning clonal composition. Further evaluations of SCClone on three real scDNA-seq datasets show that it can effectively find the underlying subclones from severely disturbed data. The SCClone software is freely available at https://github.com/qasimyu/scclone.

Integration of multiple lineage measurements from the same cell reconstructs parallel tumor evolution

Organoid evolution models complemented with integrated single-cell sequencing technology provide a powerful platform to characterize intra-tumor heterogeneity (ITH) and tumor evolution. Here, we conduct a parallel evolution experiment to mimic the tumor evolution process by evolving a colon cancer organoid model over 100 generations, spanning 6 months in time. We use single-cell whole-genome sequencing (WGS) in combination with viral lineage tracing at 12 time points to simultaneously monitor clone size, CNV states, SNV states, and viral lineage barcodes for 1,641 single cells. We integrate these measurements to construct clonal evolution trees with high resolution. We characterize the order of events in which chromosomal aberrations occur and identify aberrations that recur multiple times within the same tumor sub-population. We observe recurrent sequential loss of chromosome 4 after loss of chromosome 18 in four unique tumor clones. SNVs and CNVs identified in our organoid experiments are also frequently reported in colorectal carcinoma samples, and out of 334 patients with chromosome 18 loss in a Memorial Sloan Kettering colorectal cancer cohort, 99 (29.6%) also harbor chromosome 4 loss. Our study reconstructs tumor evolution in a colon cancer organoid model at high resolution, demonstrating an approach to identify potentially clinically relevant genomic aberrations in tumor evolution.

Bringing precision oncology to cellular resolution with single-cell genomics

Single-cell sequencing technologies have undergone rapid development and adoption by the scientific community in the past 5 years, fueling discoveries about the etiology, pathogenesis, and treatment responsiveness of individual tumor cells within cancer ecosystems. Most of the advancements in our understanding of cancer with these new technologies have focused on basic tumor biology. However, the knowledge produced by these and other studies are beginning to provide biomarkers and drug targets for clinically-relevant subpopulations within a tumor, creating opportunities for the development of biologically-informed, clone-specific combination treatment strategies. Here we provide an overview of the development of the field of single-cell cancer sequencing and provide a roadmap for shepherding these technologies from research tools to diagnostic instruments that provide high-resolution, treatment-directing details of tumors to clinical oncologists.

Natural Barcodes for Longitudinal Single Cell Tracking of Leukemic and Immune Cell Dynamics

Blood malignancies provide unique opportunities for longitudinal tracking of disease evolution following therapeutic bottlenecks and for the monitoring of changes in anti-tumor immunity. The expanding development of multi-modal single-cell sequencing technologies affords newer platforms to elucidate the mechanisms underlying these processes at unprecedented resolution. Furthermore, the identification of molecular events that can serve as in-vivo barcodes now facilitate the tracking of the trajectories of malignant and of immune cell populations over time within primary human samples, as these permit unambiguous identification of the clonal lineage of cell populations within heterogeneous phenotypes. Here, we provide an overview of the potential for chromosomal copy number changes, somatic nuclear and mitochondrial DNA mutations, single nucleotide polymorphisms, and T and B cell receptor sequences to serve as personal natural barcodes and review technical implementations in single-cell analysis workflows. Applications of these methodologies include the study of acquired therapeutic resistance and the dissection of donor- and host cellular interactions in the context of allogeneic hematopoietic stem cell transplantation.

GraPhyC: Using Consensus to Infer Tumor Evolution

We consider the problem of finding a consensus tumor evolution tree from a set of conflicting input trees. In contrast to traditional phylogenetic trees, the tumor trees we consider do not have the same set of labels applied to the leaves of each tree. We describe several distance measures between these tumor trees. Our GraPhyC algorithm solves the consensus problem using a weighted directed graph where vertices are sets of mutations and edges are weighted based on the number of times a parental relationship is observed between their constituent mutations in the input trees. We find a minimum weight spanning arborescence in this graph and prove that it minimizes the total distance to all input trees for one of our distance measures. We also describe several extensions of our GraPhyC approach. On simulated data we show that GraPhyC outperforms a baseline method and demonstrate that GraPhyC can be an effective means of computing centroids in k-medians clustering. We analyze two real sequencing datasets and find that GraPhyC is able to identify a tree not included in the set of input trees, but that contains characteristics supported by other reported evolutionary reconstructions of this tumor.

Immunoglobulin Heavy Chain High-Throughput Sequencing in Pediatric B-Precursor Acute Lymphoblastic Leukemia: Is the Clonality of the Disease at Diagnosis Related to Its Prognosis?

High-throughput sequencing (HTS) of the immunoglobulin heavy chain (IgH) locus is a recent very efficient technique to monitor minimal residual disease of B-cell precursor acute lymphoblastic leukemia (BCP-ALL). It also reveals the sequences of clonal rearrangements, therefore, the multiclonal structure, of BCP-ALL. In this study, we performed IgH HTS on the diagnostic bone marrow of 105 children treated between 2004 and 2008 in Belgium for BCP-ALL in the European Organization for Research and Treatment of Cancer (EORTC)-58951 clinical trial. Patients were included irrespectively of their outcome. We described the patterns of clonal complexity at diagnosis and investigated its association with patients' characteristics. Two indicators of clonal complexity were used, namely, the number of foster clones, described as clones with similar D-N₂-J rearrangements but other V-rearrangement and N₁-joining, and the maximum across all foster clones of the number of evolved clones from one foster clone. The maximum number of evolved clones was significantly higher in patients with t(12;21)/ETV6:RUNX1. A lower number of foster clones was associated with a higher risk group after prephase and t(12;21)/ETV6:RUNX1 genetic type. This study observes that clonal complexity as accessed by IgH HTS is linked to prognostic factors in childhood BCP-ALL, suggesting that it may be a useful diagnostic tool for BCP-ALL status and prognosis.

RDAClone: Deciphering Tumor Heterozygosity through Single-Cell Genomics Data Analysis with Robust Deep Autoencoder

Rapid advances in single-cell genomics sequencing (SCGS) have allowed researchers to characterize tumor heterozygosity with unprecedented resolution and reveal the phylogenetic relationships between tumor cells or clones. However, high sequencing error rates of current SCGS data, i.e., false positives, false negatives, and missing bases, severely limit its application. Here, we present a deep learning framework, RDAClone, to recover genotype matrices from noisy data with an extended robust deep autoencoder, cluster cells into subclones by the Louvain-Jaccard method, and further infer evolutionary relationships between subclones by the minimum spanning tree. Studies on both simulated and real datasets demonstrate its robustness and superiority in data denoising, cell clustering, and evolutionary tree reconstruction, particularly for large datasets.