Advances in understanding tumour evolution through single-cell sequencing

Isotype-aware inference of B cell clonal lineage trees from single-cell sequencing data

Single-cell RNA sequencing (scRNA-seq) enables comprehensive characterization of the micro-evolutionary processes of B cells during an adaptive immune response, capturing features of somatic hypermutation (SHM) and class switch recombination (CSR). Existing phylogenetic approaches for reconstructing B cell evolution have primarily focused on the SHM process alone. Here, we present tree inference of B cell clonal lineages (TRIBAL), an algorithm designed to optimally reconstruct the evolutionary history of B cell clonal lineages undergoing both SHM and CSR from scRNA-seq data. Through simulations, we demonstrate that TRIBAL produces more comprehensive and accurate B cell lineage trees compared to existing methods. Using real-world datasets, TRIBAL successfully recapitulates expected biological trends in a model affinity maturation system while reconstructing evolutionary histories with more parsimonious class switching than state-of-the-art methods. Thus, TRIBAL significantly improves B cell lineage tracing, useful for modeling vaccine responses, disease progression, and the identification of therapeutic antibodies.

Weighted centroid trees: a general approach to summarize phylogenies in single-labeled tumor mutation tree inference

MOTIVATION

Tumor trees, which depict the evolutionary process of cancer, provide a backbone for discovering recurring evolutionary processes in cancer. While they are not the primary information extracted from genomic data, they are valuable for this purpose. One such extraction method involves summarizing multiple trees into a single representative tree, such as consensus trees or supertrees.

RESULTS

We define the "weighted centroid tree problem" to find the centroid tree of a set of single-labeled rooted trees through the following steps: (i) mapping the given trees into the Euclidean space, (ii) computing the weighted centroid matrix of the mapped trees, and (iii) finding the nearest mapped tree (NMTP) to the centroid matrix. We show that this setup encompasses previously studied parent-child and ancestor-descendent metrics as well as the GraPhyC and TuELiP consensus tree algorithms. Moreover, we show that, while the NMTP problem is polynomial-time solvable for the adjacency embedding, it is NP-hard for ancestry and distance mappings. We introduce integer linear programs for NMTP in different setups where we also provide a new algorithm for the case of ancestry embedding called 2-AncL2, that uses a novel weighting scheme for ancestry signals. Our experimental results show that 2-AncL2 has a superior performance compared to available consensus tree algorithms. We also illustrate our setup's application on providing representative trees for a large real breast cancer dataset, deducing that the cluster centroid trees summarize reliable evolutionary information about the original dataset.

AVAILABILITY AND IMPLEMENTATION

https://github.com/vasei/WAncILP.

Oncotree2vec - a method for embedding and clustering of tumor mutation trees

MOTIVATION

Understanding the genomic heterogeneity of tumors is an important task in computational oncology, especially in the context of finding personalized treatments based on the genetic profile of each patient's tumor. Tumor clustering that takes into account the temporal order of genetic events, as represented by tumor mutation trees, is a powerful approach for grouping together patients with genetically and evolutionarily similar tumors and can provide insights into discovering tumor subtypes, for more accurate clinical diagnosis and prognosis.

RESULTS

Here, we propose oncotree2vec, a method for clustering tumor mutation trees by learning vector representations of mutation trees that capture the different relationships between subclones in an unsupervised manner. Learning low-dimensional tree embeddings facilitates the visualization of relations between trees in large cohorts and can be used for downstream analyses, such as deep learning approaches for single-cell multi-omics data integration. We assessed the performance and the usefulness of our method in three simulation studies and on two real datasets: a cohort of 43 trees from six cancer types with different branching patterns corresponding to different modes of spatial tumor evolution and a cohort of 123 AML mutation trees.

AVAILABILITY AND IMPLEMENTATION

https://github.com/cbg-ethz/oncotree2vec.

Enhanced Complement Expression in the Tumor Microenvironment Following Neoadjuvant Therapy: Implications for Immunomodulation and Survival in Pancreatic Ductal Adenocarcinoma

Background

Neoadjuvant therapy (NAT) is increasingly being used for pancreatic ductal adenocarcinoma (PDAC) treatment. However, its specific effects on carcinoma cells and the tumor microenvironment (TME) are not fully understood. This study aims to investigate how NAT differentially impacts PDAC's carcinoma cells and TME.

Methods

Spatial transcriptomics was used to compare gene expression profiles in carcinoma cells and the TME between 23 NAT-treated and 13 NAT-naïve PDAC patients, correlating with their clinicopathologic features. Analysis of an online single-nucleus RNA sequencing (snRNA-seq) dataset was performed for validation of the specific cell types responsible for NAT-induced gene expression alterations.

Results

NAT not only induces apoptosis and inhibits proliferation in carcinoma cells but also significantly remodels the TME. Notably, NAT induces a coordinated upregulation of multiple key complement genes (C3, C1S, C1R, C4B and C7) in the TME, making the complement pathway one of the most significantly affected pathways by NAT. Patients with higher TME complement expression following NAT exhibit improved overall survival. These patients also exhibit increased immunomodulatory and neurotrophic cancer-associated fibroblasts (CAFs); more CD4+ T cells, monocytes, and mast cells; and reduced immune exhaustion gene expression. snRNA-seq analysis demonstrates C3 complement was specifically upregulated in CAFs but not in other stroma cell types.

Conclusions

NAT can enhance complement production and signaling within the TME, which is associated with reduced immunosuppression in PDAC. These findings suggest that local complement dynamics could serve as a novel biomarker for prognosis, evaluating treatment response and resistance, and guiding therapeutic strategies in NAT-treated PDAC patients.

scGAL: unmask tumor clonal substructure by jointly analyzing independent single-cell copy number and scRNA-seq data

BACKGROUND

Accurately deciphering clonal copy number substructure can provide insights into the evolutionary mechanism of cancer, and clustering single-cell copy number profiles has become an effective means to unmask intra-tumor heterogeneity (ITH). However, copy numbers inferred from single-cell DNA sequencing (scDNA-seq) data are error-prone due to technically confounding factors such as amplification bias and allele-dropout, and this makes it difficult to precisely identify the ITH.

RESULTS

We introduce a hybrid model called scGAL to infer clonal copy number substructure. It combines an autoencoder with a generative adversarial network to jointly analyze independent single-cell copy number profiles and gene expression data from same cell line. Under an adversarial learning framework, scGAL exploits complementary information from gene expression data to relieve the effects of noise in copy number data, and learns latent representations of scDNA-seq cells for accurate inference of the ITH. Evaluation results on three real cancer datasets suggest scGAL is able to accurately infer clonal architecture and surpasses other similar methods. In addition, assessment of scGAL on various simulated datasets demonstrates its high robustness against the changes of data size and distribution. scGAL can be accessed at: https://github.com/zhyu-lab/scgal .

CONCLUSIONS

Joint analysis of independent single-cell copy number and gene expression data from a same cell line can effectively exploit complementary information from individual omics, and thus gives more refined indication of clonal copy number substructure.

CoT: a transformer-based method for inferring tumor clonal copy number substructure from scDNA-seq data

Single-cell DNA sequencing (scDNA-seq) has been an effective means to unscramble intra-tumor heterogeneity, while joint inference of tumor clones and their respective copy number profiles remains a challenging task due to the noisy nature of scDNA-seq data. We introduce a new bioinformatics method called CoT for deciphering clonal copy number substructure. The backbone of CoT is a Copy number Transformer autoencoder that leverages multi-head attention mechanism to explore correlations between different genomic regions, and thus capture global features to create latent embeddings for the cells. CoT makes it convenient to first infer cell subpopulations based on the learned embeddings, and then estimate single-cell copy numbers through joint analysis of read counts data for the cells belonging to the same cluster. This exploitation of clonal substructure information in copy number analysis helps to alleviate the effect of read counts non-uniformity, and yield robust estimations of the tumor copy numbers. Performance evaluation on synthetic and real datasets showcases that CoT outperforms the state of the arts, and is highly useful for deciphering clonal copy number substructure.

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.

Advancements in MDM2 inhibition: Clinical and pre-clinical investigations of combination therapeutic regimens

Cancer cells often depend on multiple pathways for their growth and survival, resulting in therapeutic resistance and the limited effectiveness of treatments. Combination therapy has emerged as a favorable approach to enhance treatment efficacy and minimize acquired resistance and harmful side effects. The murine double minute 2 (MDM2) protein regulates cellular proliferation and promotes cancer-related activities by negatively regulating the tumor suppressor protein p53. MDM2 aberrations have been reported in a variety of human cancers, making it an appealing target for cancer therapy. As a result, several small-molecule MDM2 inhibitors have been developed and are currently being investigated in clinical studies. Nevertheless, it has been shown that the inhibition of MDM2 alone is inadequate to achieve long-term suppression of tumor growth, thus prompting the need for further investigation into combination therapeutic strategies. In this review, possible clinical and preclinical MDM2 combination inhibitor regimens are thoroughly analyzed and discussed. It provides a rationale for combining MDM2 inhibitors with other therapeutic approaches in the management of cancer, taking into consideration ongoing clinical trials that evaluate the combination of MDM2 inhibitors. The review explores the current status of MDM2 inhibitors in combination with chemotherapy or targeted therapy, as well as promising approach of combining MDM2 inhibitors with immunotherapy. In addition, it investigates the function of PROTACs as MDM2 degraders in cancer treatment. A comprehensive examination of these combination regimens highlights the potential for advancing MDM2-inhibitor therapy and improving clinical outcomes for cancer patients and establishes the foundation for future research and development in this promising area of study.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Pharmacogenomics: Driving Personalized Medicine

Personalized medicine tailors therapies, disease prevention, and health maintenance to the individual, with pharmacogenomics serving as a key tool to improve outcomes and prevent adverse effects. Advances in genomics have transformed pharmacogenetics, traditionally focused on single gene-drug pairs, into pharmacogenomics, encompassing all "-omics" fields (e.g., proteomics, transcriptomics, metabolomics, and metagenomics). This review summarizes basic genomics principles relevant to translation into therapies, assessing pharmacogenomics' central role in converging diverse elements of personalized medicine. We discuss genetic variations in pharmacogenes (drug-metabolizing enzymes, drug transporters, and receptors), their clinical relevance as biomarkers, and the legacy of decades of research in pharmacogenetics. All types of therapies, including proteins, nucleic acids, viruses, cells, genes, and irradiation, can benefit from genomics, expanding the role of pharmacogenomics across medicine. Food and Drug Administration approvals of personalized therapeutics involving biomarkers increase rapidly, demonstrating the growing impact of pharmacogenomics. A beacon for all therapeutic approaches, molecularly targeted cancer therapies highlight trends in drug discovery and clinical applications. To account for human complexity, multicomponent biomarker panels encompassing genetic, personal, and environmental factors can guide diagnosis and therapies, increasingly involving artificial intelligence to cope with extreme data complexities. However, clinical application encounters substantial hurdles, such as unknown validity across ethnic groups, underlying bias in health care, and real-world validation. This review address the underlying science and technologies germane to pharmacogenomics and personalized medicine, integrated with economic, ethical, and regulatory issues, providing insights into the current status and future direction of health care. SIGNIFICANCE STATEMENT: Personalized medicine aims to optimize health care for the individual patients with use of predictive biomarkers to improve outcomes and prevent adverse effects. Pharmacogenomics drives biomarker discovery and guides the development of targeted therapeutics. This review addresses basic principles and current trends in pharmacogenomics, with large-scale data repositories accelerating medical advances. The impact of pharmacogenomics is discussed, along with hurdles impeding broad clinical implementation, in the context of clinical care, ethics, economics, and regulatory affairs.

Joint inference of exclusivity patterns and recurrent trajectories from tumor mutation trees

Cancer progression is an evolutionary process shaped by both deterministic and stochastic forces. Multi-region and single-cell sequencing of tumors enable high-resolution reconstruction of the mutational history of each tumor and highlight the extensive diversity across tumors and patients. Resolving the interactions among mutations and recovering recurrent evolutionary processes may offer greater opportunities for successful therapeutic strategies. To this end, we present a novel probabilistic framework, called TreeMHN, for the joint inference of exclusivity patterns and recurrent trajectories from a cohort of intra-tumor phylogenetic trees. Through simulations, we show that TreeMHN outperforms existing alternatives that can only focus on one aspect of the task. By analyzing datasets of blood, lung, and breast cancers, we find the most likely evolutionary trajectories and mutational patterns, consistent with and enriching our current understanding of tumorigenesis. Moreover, TreeMHN facilitates the prediction of tumor evolution and provides probabilistic measures on the next mutational events given a tumor tree, a prerequisite for evolution-guided treatment strategies.

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Reconstructing clonal tree for phylo-phenotypic characterization of cancer using single-cell transcriptomics

Functional characterization of the cancer clones can shed light on the evolutionary mechanisms driving cancer's proliferation and relapse mechanisms. Single-cell RNA sequencing data provide grounds for understanding the functional state of cancer as a whole; however, much research remains to identify and reconstruct clonal relationships toward characterizing the changes in functions of individual clones. We present PhylEx that integrates bulk genomics data with co-occurrences of mutations from single-cell RNA sequencing data to reconstruct high-fidelity clonal trees. We evaluate PhylEx on synthetic and well-characterized high-grade serous ovarian cancer cell line datasets. PhylEx outperforms the state-of-the-art methods both when comparing capacity for clonal tree reconstruction and for identifying clones. We analyze high-grade serous ovarian cancer and breast cancer data to show that PhylEx exploits clonal expression profiles beyond what is possible with expression-based clustering methods and clear the way for accurate inference of clonal trees and robust phylo-phenotypic analysis of cancer.

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.

Identification of immune cell function in breast cancer by integrating multiple single-cell data

Breast cancer has now become the most commonly diagnosed cancer worldwide. It is a highly complex and heterogeneous disease that comprises distinct histological features and treatment response. With the development of molecular biology and immunology, immunotherapy has become a new field of breast cancer treatment. Identifying cell-type-specific genes critical to the immune microenvironment contributes to breast cancer treatment. Single-cell RNA sequencing (scRNA-seq) technology could serve as a powerful tool to analyze cellular genetic information at single-cell resolution and to uncover the gene expression status of each cell, thus allowing comprehensive assessment of intercellular heterogeneity. Because of the influence of sample size and sequencing depth, the specificity of genes in different cell types for breast cancer cannot be fully revealed. Therefore, the present study integrated two public breast cancer scRNA-seq datasets aiming to investigate the functions of different type of immune cells in tumor microenvironment. We identified total five significant differential expressed genes of B cells, T cells and macrophage and explored their functions and immune mechanisms in breast cancer. Finally, we performed functional annotation analyses using the top fifteen differentially expressed genes in each immune cell type to discover the immune-related pathways and gene ontology (GO) terms.

Single cell RNA-sequencing: A powerful yet still challenging technology to study cellular heterogeneity

Almost all biomedical research to date has relied upon mean measurements from cell populations, however it is well established that what it is observed at this macroscopic level can be the result of many interactions of several different single cells. Thus, the observable macroscopic 'average' cannot outright be used as representative of the 'average cell'. Rather, it is the resulting emerging behaviour of the actions and interactions of many different cells. Single-cell RNA sequencing (scRNA-Seq) enables the comparison of the transcriptomes of individual cells. This provides high-resolution maps of the dynamic cellular programmes allowing us to answer fundamental biological questions on their function and evolution. It also allows to address medical questions such as the role of rare cell populations contributing to disease progression and therapeutic resistance. Furthermore, it provides an understanding of context-specific dependencies, namely the behaviour and function that a cell has in a specific context, which can be crucial to understand some complex diseases, such as diabetes, cardiovascular disease and cancer. Here, we provide an overview of scRNA-Seq, including a comparative review of emerging technologies and computational pipelines. We discuss the current and emerging applications and focus on tumour heterogeneity a clear example of how scRNA-Seq can provide new understanding of a complex disease. Additionally, we review the limitations and highlight the need of powerful computational pipelines and reproducible protocols for the broader acceptance of this technique in basic and clinical research.

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.

Scelestial: Fast and accurate single-cell lineage tree inference based on a Steiner tree approximation algorithm

Single-cell genome sequencing provides a highly granular view of biological systems but is affected by high error rates, allelic amplification bias, and uneven genome coverage. This creates a need for data-specific computational methods, for purposes such as for cell lineage tree inference. The objective of cell lineage tree reconstruction is to infer the evolutionary process that generated a set of observed cell genomes. Lineage trees may enable a better understanding of tumor formation and growth, as well as of organ development for healthy body cells. We describe a method, Scelestial, for lineage tree reconstruction from single-cell data, which is based on an approximation algorithm for the Steiner tree problem and is a generalization of the neighbor-joining method. We adapt the algorithm to efficiently select a limited subset of potential sequences as internal nodes, in the presence of missing values, and to minimize cost by lineage tree-based missing value imputation. In a comparison against seven state-of-the-art single-cell lineage tree reconstruction algorithms—BitPhylogeny, OncoNEM, SCITE, SiFit, SASC, SCIPhI, and SiCloneFit—on simulated and real single-cell tumor samples, Scelestial performed best at reconstructing trees in terms of accuracy and run time. Scelestial has been implemented in C++. It is also available as an R package named RScelestial.

SECEDO: SNV-based subclone detection using ultra-low coverage single-cell DNA sequencing

MOTIVATION

Several recently developed single-cell DNA sequencing technologies enable whole-genome sequencing of thousands of cells. However, the ultra-low coverage of the sequenced data (<0.05× per cell) mostly limits their usage to the identification of copy number alterations in multi-megabase segments. Many tumors are not copy number-driven, and thus single-nucleotide variant (SNV)-based subclone detection may contribute to a more comprehensive view on intra-tumor heterogeneity. Due to the low coverage of the data, the identification of SNVs is only possible when superimposing the sequenced genomes of hundreds of genetically similar cells. Thus, we have developed a new approach to efficiently cluster tumor cells based on a Bayesian filtering approach of relevant loci and exploiting read overlap and phasing.

RESULTS

We developed Single Cell Data Tumor Clusterer (SECEDO, lat. 'to separate'), a new method to cluster tumor cells based solely on SNVs, inferred on ultra-low coverage single-cell DNA sequencing data. We applied SECEDO to a synthetic dataset simulating 7250 cells and eight tumor subclones from a single patient and were able to accurately reconstruct the clonal composition, detecting 92.11% of the somatic SNVs, with the smallest clusters representing only 6.9% of the total population. When applied to five real single-cell sequencing datasets from a breast cancer patient, each consisting of ≈2000 cells, SECEDO was able to recover the major clonal composition in each dataset at the original coverage of 0.03×, achieving an Adjusted Rand Index (ARI) score of ≈0.6. The current state-of-the-art SNV-based clustering method achieved an ARI score of ≈0, even after merging cells to create higher coverage data (factor 10 increase), and was only able to match SECEDOs performance when pooling data from all five datasets, in addition to artificially increasing the sequencing coverage by a factor of 7. Variant calling on the resulting clusters recovered more than twice as many SNVs as would have been detected if calling on all cells together. Further, the allelic ratio of the called SNVs on each subcluster was more than double relative to the allelic ratio of the SNVs called without clustering, thus demonstrating that calling variants on subclones, in addition to both increasing sensitivity of SNV detection and attaching SNVs to subclones, significantly increases the confidence of the called variants.

AVAILABILITY AND IMPLEMENTATION

SECEDO is implemented in C++ and is publicly available at https://github.com/ratschlab/secedo. Instructions to download the data and the evaluation code to reproduce the findings in this paper are available at: https://github.com/ratschlab/secedo-evaluation. The code and data of the submitted version are archived at: https://doi.org/10.5281/zenodo.6516955.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Accurate Identification of Subclones in Tumor Genomes

Understanding intratumor heterogeneity is critical for studying tumorigenesis and designing personalized treatments. To decompose the mixed cell population in a tumor, subclones are inferred computationally based on variant allele frequency (VAF) from bulk sequencing data. In this study, we showed that sequencing depth, mean VAF, and variance of VAF of a subclone are confounded. Without considering this effect, current methods require deep-sequencing data (>300× depth) to reliably infer subclones. Here, we present a novel algorithm that incorporates depth-variance and mean-variance dependencies in a clustering error model and successfully identifies subclones in tumors sequenced at depths of as low as 30×. We implemented the algorithm as a model-based adaptive grouping of subclones (MAGOS) method. Analyses of computer simulated data and empirical sequencing data showed that MAGOS outperformed existing methods on minimum sequencing depth, decomposition accuracy, and computation efficiency. The most prominent improvements were observed in analyzing tumors sequenced at depths between 30× and 200×, whereas the performance was comparable between MAGOS and existing methods on deeply sequenced tumors. MAGOS supports analysis of single-nucleotide variants and copy number variants from a single sample or multiple samples of a tumor. We applied MAGOS to whole-exome data of late-stage liver cancers and discovered that high subclone count in a tumor was a significant risk factor of poor prognosis. Lastly, our analysis suggested that sequencing multiple samples of the same tumor at standard depth is more cost-effective and robust for subclone characterization than deep sequencing a single sample. MAGOS is available at github (https://github.com/liliulab/magos).

Developments of Conventional and Microfluidic Flow Cytometry Enabling High-Throughput Characterization of Single Cells

This article first reviews scientific meanings of single-cell analysis by highlighting two key scientific problems: landscape reconstruction of cellular identities during dynamic immune processes and mechanisms of tumor origin and evolution. Secondly, the article reviews clinical demands of single-cell analysis, which are complete blood counting enabled by optoelectronic flow cytometry and diagnosis of hematologic malignancies enabled by multicolor fluorescent flow cytometry. Then, this article focuses on the developments of optoelectronic flow cytometry for the complete blood counting by comparing conventional counterparts of hematology analyzers (e.g., DxH 900 of Beckman Coulter, XN-1000 of Sysmex, ADVIA 2120i of Siemens, and CELL-DYN Ruby of Abbott) and microfluidic counterparts (e.g., microfluidic impedance and imaging flow cytometry). Future directions of optoelectronic flow cytometry are indicated where intrinsic rather than dependent biophysical parameters of blood cells must be measured, and they can replace blood smears as the gold standard of blood analysis in the near future.

Accounting for errors in data improves divergence time estimates in single-cell cancer evolution

Single-cell sequencing provides a new way to explore the evolutionary history of cells. Compared to traditional bulk sequencing, where a population of heterogeneous cells is pooled to form a single observation, single-cell sequencing isolates and amplifies genetic material from individual cells, thereby preserving the information about the origin of the sequences. However, single-cell data is more error-prone than bulk sequencing data due to the limited genomic material available per cell. Here, we present error and mutation models for evolutionary inference of single-cell data within a mature and extensible Bayesian framework, BEAST2. Our framework enables integration with biologically informative models such as relaxed molecular clocks and population dynamic models. Our simulations show that modeling errors increase the accuracy of relative divergence times and substitution parameters. We reconstruct the phylogenetic history of a colorectal cancer patient and a healthy patient from single-cell DNA sequencing data. We find that the estimated times of terminal splitting events are shifted forward in time compared to models which ignore errors. We observed that not accounting for errors can overestimate the phylogenetic diversity in single-cell DNA sequencing data. We estimate that 30-50% of the apparent diversity can be attributed to error. Our work enables a full Bayesian approach capable of accounting for errors in the data within the integrative Bayesian software framework BEAST2.

CONET: copy number event tree model of evolutionary tumor history for single-cell data

Copy number alterations constitute important phenomena in tumor evolution. Whole genome single-cell sequencing gives insight into copy number profiles of individual cells, but is highly noisy. Here, we propose CONET, a probabilistic model for joint inference of the evolutionary tree on copy number events and copy number calling. CONET employs an efficient, regularized MCMC procedure to search the space of possible model structures and parameters. We introduce a range of model priors and penalties for efficient regularization. CONET reveals copy number evolution in two breast cancer samples, and outperforms other methods in tree reconstruction, breakpoint identification and copy number calling.

Single-cell sequencing: a promising approach for uncovering the mechanisms of tumor metastasis

Single-cell sequencing (SCS) is an emerging high-throughput technology that can be used to study the genomics, transcriptomics, and epigenetics at a single cell level. SCS is widely used in the diagnosis and treatment of various diseases, including cancer. Over the years, SCS has gradually become an effective clinical tool for the exploration of tumor metastasis mechanisms and the development of treatment strategies. Currently, SCS can be used not only to analyze metastasis-related malignant biological characteristics, such as tumor heterogeneity, drug resistance, and microenvironment, but also to construct metastasis-related cell maps for predicting and monitoring the dynamics of metastasis. SCS is also used to identify therapeutic targets related to metastasis as it provides insights into the distribution of tumor cell subsets and gene expression differences between primary and metastatic tumors. Additionally, SCS techniques in combination with artificial intelligence (AI) are used in liquid biopsy to identify circulating tumor cells (CTCs), thereby providing a novel strategy for treating tumor metastasis. In this review, we summarize the potential applications of SCS in the field of tumor metastasis and discuss the prospects and limitations of SCS to provide a theoretical basis for finding therapeutic targets and mechanisms of metastasis.

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.

CellPhy: accurate and fast probabilistic inference of single-cell phylogenies from scDNA-seq data

We introduce CellPhy, a maximum likelihood framework for inferring phylogenetic trees from somatic single-cell single-nucleotide variants. CellPhy leverages a finite-site Markov genotype model with 16 diploid states and considers amplification error and allelic dropout. We implement CellPhy into RAxML-NG, a widely used phylogenetic inference package that provides statistical confidence measurements and scales well on large datasets with hundreds or thousands of cells. Comprehensive simulations suggest that CellPhy is more robust to single-cell genomics errors and outperforms state-of-the-art methods under realistic scenarios, both in accuracy and speed. CellPhy is freely available at https://github.com/amkozlov/cellphy .

Independent Drug Action in Combination Therapy: Implications for Precision Oncology

Combination therapies are superior to monotherapy for many cancers. This advantage was historically ascribed to the ability of combinations to address tumor heterogeneity, but synergistic interaction is now a common explanation as well as a design criterion for new combinations. We review evidence that independent drug action, described in 1961, explains the efficacy of many practice-changing combination therapies: it provides populations of patients with heterogeneous drug sensitivities multiple chances of benefit from at least one drug. Understanding response heterogeneity could reveal predictive or pharmacodynamic biomarkers for more precise use of existing drugs and realize the benefits of additivity or synergy.

OUP accepted manuscript

MOTIVATION

Cancer develops through a process of clonal evolution in which an initially healthy cell gives rise to progeny gradually differentiating through the accumulation of genetic and epigenetic mutations. These mutations can take various forms, including single-nucleotide variants (SNVs), copy number alterations (CNAs) or structural variations (SVs), with each variant type providing complementary insights into tumor evolution as well as offering distinct challenges to phylogenetic inference.

RESULTS

In this work, we develop a tumor phylogeny method, TUSV-ext, which incorporates SNVs, CNAs and SVs into a single inference framework. We demonstrate on simulated data that the method produces accurate tree inferences in the presence of all three variant types. We further demonstrate the method through application to real prostate tumor data, showing how our approach to coordinated phylogeny inference and clonal construction with all three variant types can reveal a more complicated clonal structure than is suggested by prior work, consistent with extensive polyclonal seeding or migration.

AVAILABILITY AND IMPLEMENTATION

https://github.com/CMUSchwartzLab/TUSV-ext.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

OUP accepted manuscript

Motivation

Identifying cell types and their abundances and how these evolve during tumor progression is critical to understanding the mechanisms of metastasis and identifying predictors of metastatic potential that can guide the development of new diagnostics or therapeutics. Single-cell RNA sequencing (scRNA-seq) has been especially promising in resolving heterogeneity of expression programs at the single-cell level, but is not always feasible, e.g. for large cohort studies or longitudinal analysis of archived samples. In such cases, clonal subpopulations may still be inferred via genomic deconvolution, but deconvolution methods have limited ability to resolve fine clonal structure and may require reference cell type profiles that are missing or imprecise. Prior methods can eliminate the need for reference profiles but show unstable performance when few bulk samples are available.

Results

In this work, we develop a new method using reference scRNA-seq to interpret sample collections for which only bulk RNA-seq is available for some samples, e.g. clonally resolving archived primary tissues using scRNA-seq from metastases. By integrating such information in a Quadratic Programming framework, our method can recover more accurate cell types and corresponding cell type abundances in bulk samples. Application to a breast tumor bone metastases dataset confirms the power of scRNA-seq data to improve cell type inference and quantification in same-patient bulk samples.

Availability and implementation

Source code is available on Github at https://github.com/CMUSchwartzLab/RADs.

Statistical tests for intra-tumour clonal co-occurrence and exclusivity

Tumour progression is an evolutionary process in which different clones evolve over time, leading to intra-tumour heterogeneity. Interactions between clones can affect tumour evolution and hence disease progression and treatment outcome. Intra-tumoural pairs of mutations that are overrepresented in a co-occurring or clonally exclusive fashion over a cohort of patient samples may be suggestive of a synergistic effect between the different clones carrying these mutations. We therefore developed a novel statistical testing framework, called GeneAccord, to identify such gene pairs that are altered in distinct subclones of the same tumour. We analysed our framework for calibration and power. By comparing its performance to baseline methods, we demonstrate that to control type I errors, it is essential to account for the evolutionary dependencies among clones. In applying GeneAccord to the single-cell sequencing of a cohort of 123 acute myeloid leukaemia patients, we find 1 clonally co-occurring and 8 clonally exclusive gene pairs. The clonally exclusive pairs mostly involve genes of the key signalling pathways.

Tumours typically display high levels of heterogeneity, not only between different tumours but also within a single one. Intra-tumour heterogeneity results from an evolutionary process, giving rise to different populations of cancer cells known as clones. How clones interact may affect tumour evolution, which in turn determines disease progression and treatment outcome. In practice, we may observe pairs of mutations that co-occur in clones or exclude each other more often than we would expect for a given cohort of patient samples. Exclusive pairs are suggestive that clones carrying one or the other mutation may cooperate in the evolutionary process. Targeting only one of them may then suffice to alter the tumour evolution. Therefore it is critical to have statistical methods which allow us to identify such pairs. GeneAccord is a novel statistical testing framework we developed especially to identify pairs of genes altered in distinct clones of the same tumour. Accounting for the evolutionary dependencies among clones emerged as critical to adequately control testing errors. In a cohort of 123 acute myeloid leukaemia patients, GeneAccord identified one clonally co-occurring and eight clonally exclusive gene pairs. The latter predominantly involved genes of key signalling pathways.

Accurate and scalable variant calling from single cell DNA sequencing data with ProSolo

Accurate single cell mutational profiles can reveal genomic cell-to-cell heterogeneity. However, sequencing libraries suitable for genotyping require whole genome amplification, which introduces allelic bias and copy errors. The resulting data violates assumptions of variant callers developed for bulk sequencing. Thus, only dedicated models accounting for amplification bias and errors can provide accurate calls. We present ProSolo for calling single nucleotide variants from multiple displacement amplified (MDA) single cell DNA sequencing data. ProSolo probabilistically models a single cell jointly with a bulk sequencing sample and integrates all relevant MDA biases in a site-specific and scalable-because computationally efficient-manner. This achieves a higher accuracy in calling and genotyping single nucleotide variants in single cells in comparison to state-of-the-art tools and supports imputation of insufficiently covered genotypes, when downstream tools cannot handle missing data. Moreover, ProSolo implements the first approach to control the false discovery rate reliably and flexibly. ProSolo is implemented in an extendable framework, with code and usage at: https://github.com/prosolo/prosolo.

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.

VDAC Modulation of Cancer Metabolism: Advances and Therapeutic Challenges

Most anionic metabolites including respiratory substrates, glycolytic adenosine triphosphate (ATP), and small cations that enter mitochondria, and mitochondrial ATP moving to the cytosol, cross the outer mitochondrial membrane (OMM) through voltage dependent anion channels (VDAC). The closed states of VDAC block the passage of anionic metabolites, and increase the flux of small cations, including calcium. Consequently, physiological or pharmacological regulation of VDAC opening, by conditioning the magnitude of both anion and cation fluxes, is a major contributor to mitochondrial metabolism. Tumor cells display a pro-proliferative Warburg phenotype characterized by enhanced aerobic glycolysis in the presence of partial suppression of mitochondrial metabolism. The heterogeneous and flexible metabolic traits of most human tumors render cells able to adapt to the constantly changing energetic and biosynthetic demands by switching between predominantly glycolytic or oxidative phenotypes. Here, we describe the biological consequences of changes in the conformational state of VDAC for cancer metabolism, the mechanisms by which VDAC-openers promote cancer cell death, and the advantages of VDAC opening as a valuable pharmacological target. Particular emphasis is given to the endogenous regulation of VDAC by free tubulin and the effects of VDAC-tubulin antagonists in cancer cells. Because of its function and location, VDAC operates as a switch to turn-off mitochondrial metabolism (closed state) and increase aerobic glycolysis (pro-Warburg), or to turn-on mitochondrial metabolism (open state) and decrease glycolysis (anti-Warburg). A better understanding of the role of VDAC regulation in tumor progression is relevant both for cancer biology and for developing novel cancer chemotherapies.

Conifer: clonal tree inference for tumor heterogeneity with single-cell and bulk sequencing data

BACKGROUND

Genetic heterogeneity of a cancer tumor that develops during clonal evolution is one of the reasons for cancer treatment failure, by increasing the chance of drug resistance. Clones are cell populations with different genotypes, resulting from differences in somatic mutations that occur and accumulate during cancer development. An appropriate approach for identifying clones is determining the variant allele frequency of mutations that occurred in the tumor. Although bulk sequencing data can be used to provide that information, the frequencies are not informative enough for identifying different clones with the same prevalence and their evolutionary relationships. On the other hand, single-cell sequencing data provides valuable information about branching events in the evolution of a cancerous tumor. However, the temporal order of mutations may be determined with ambiguities using only single-cell data, while variant allele frequencies from bulk sequencing data can provide beneficial information for inferring the temporal order of mutations with fewer ambiguities.

RESULT

In this study, a new method called Conifer (ClONal tree Inference For hEterogeneity of tumoR) is proposed which combines aggregated variant allele frequency from bulk sequencing data with branching event information from single-cell sequencing data to more accurately identify clones and their evolutionary relationships. It is proven that the accuracy of clone identification and clonal tree inference is increased by using Conifer compared to other existing methods on various sets of simulated data. In addition, it is discussed that the evolutionary tree provided by Conifer on real cancer data sets is highly consistent with information in both bulk and single-cell data.

CONCLUSIONS

In this study, we have provided an accurate and robust method to identify clones of tumor heterogeneity and their evolutionary history by combining single-cell and bulk sequencing data.

Computational Reconstruction of Clonal Hierarchies From Bulk Sequencing Data of Acute Myeloid Leukemia Samples

Acute myeloid leukemia is an aggressive cancer of the blood forming system. The malignant cell population is composed of multiple clones that evolve over time. Clonal data reflect the mechanisms governing treatment response and relapse. Single cell sequencing provides most direct insights into the clonal composition of the leukemic cells, however it is still not routinely available in clinical practice. In this work we develop a computational algorithm that allows identifying all clonal hierarchies that are compatible with bulk variant allele frequencies measured in a patient sample. The clonal hierarchies represent descendance relations between the different clones and reveal the order in which mutations have been acquired. The proposed computational approach is tested using single cell sequencing data that allow comparing the outcome of the algorithm with the true structure of the clonal hierarchy. We investigate which problems occur during reconstruction of clonal hierarchies from bulk sequencing data. Our results suggest that in many cases only a small number of possible hierarchies fits the bulk data. This implies that bulk sequencing data can be used to obtain insights in clonal evolution.

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.

Leveraging Single-Cell Approaches in Cancer Precision Medicine

Cancer precision medicine aims to improve patient outcomes by tailoring treatment to the unique genomic background of a tumor. However, efforts to develop prognostic and drug response biomarkers largely rely on bulk 'omic' data, which fails to capture intratumor heterogeneity (ITH) and deconvolve signals from normal versus tumor cells. These shortcomings in measuring clinically relevant features are being addressed with single-cell technologies, which provide a fine-resolution map of the genetic and phenotypic heterogeneity in tumors and their microenvironment, as well as an improved understanding of the patterns of subclonal tumor populations. Here we present recent advances in the application of single-cell technologies, towards gaining a deeper understanding of ITH and evolution, and potential applications in developing personalized therapeutic strategies.

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).

PhISCS-BnB: a fast branch and bound algorithm for the perfect tumor phylogeny reconstruction problem

MOTIVATION

Recent advances in single-cell sequencing (SCS) offer an unprecedented insight into tumor emergence and evolution. Principled approaches to tumor phylogeny reconstruction via SCS data are typically based on general computational methods for solving an integer linear program, or a constraint satisfaction program, which, although guaranteeing convergence to the most likely solution, are very slow. Others based on Monte Carlo Markov Chain or alternative heuristics not only offer no such guarantee, but also are not faster in practice. As a result, novel methods that can scale up to handle the size and noise characteristics of emerging SCS data are highly desirable to fully utilize this technology.

RESULTS

We introduce PhISCS-BnB (phylogeny inference using SCS via branch and bound), a branch and bound algorithm to compute the most likely perfect phylogeny on an input genotype matrix extracted from an SCS dataset. PhISCS-BnB not only offers an optimality guarantee, but is also 10-100 times faster than the best available methods on simulated tumor SCS data. We also applied PhISCS-BnB on a recently published large melanoma dataset derived from the sublineages of a cell line involving 20 clones with 2367 mutations, which returned the optimal tumor phylogeny in <4 h. The resulting phylogeny agrees with and extends the published results by providing a more detailed picture on the clonal evolution of the tumor.

AVAILABILITY AND IMPLEMENTATION

https://github.com/algo-cancer/PhISCS-BnB.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-cell RNA sequencing in cancer research

Single-cell RNA sequencing (scRNA-seq), a technology that analyzes transcriptomes of complex tissues at single-cell levels, can identify differential gene expression and epigenetic factors caused by mutations in unicellular genomes, as well as new cell-specific markers and cell types. scRNA-seq plays an important role in various aspects of tumor research. It reveals the heterogeneity of tumor cells and monitors the progress of tumor development, thereby preventing further cellular deterioration. Furthermore, the transcriptome analysis of immune cells in tumor tissue can be used to classify immune cells, their immune escape mechanisms and drug resistance mechanisms, and to develop effective clinical targeted therapies combined with immunotherapy. Moreover, this method enables the study of intercellular communication and the interaction of tumor cells and non-malignant cells to reveal their role in carcinogenesis. scRNA-seq provides new technical means for further development of tumor research and is expected to make significant breakthroughs in this field. This review focuses on the principles of scRNA-seq, with an emphasis on the application of scRNA-seq in tumor heterogeneity, pathogenesis, and treatment.

Analysis of Intratumoral Heterogeneity in Myelodysplastic Syndromes with Isolated del(5q) Using a Single Cell Approach

Simple Summary

Myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic stem cell malignancies characterized by ineffective differentiation of one or more bone marrow cell lineages. Only 50% of patients with de novo MDS will be found to have cytogenetic abnormalities, of which del(5q) is the most common. In 10% of MDS cases, del(5q) is found as a sole abnormality. In this work, a single cell approach was used to analyze intratumoral heterogeneity in four patients with MDS with isolated del(5q). We were able to observe that an ancestral event in one patient can appear as a secondary hit in another one, thus reflecting the high intratumoral heterogeneity in MDS with isolated del(5q) and the importance of patient-specific molecular characterization.

Abstract

Myelodysplastic syndromes (MDS) are a heterogeneous group of hematological diseases. Among them, the most well characterized subtype is MDS with isolated chromosome 5q deletion (MDS del(5q)), which is the only one defined by a cytogenetic abnormality that makes these patients candidates to be treated with lenalidomide. During the last decade, single cell (SC) analysis has emerged as a powerful tool to decipher clonal architecture and to further understand cancer and other diseases at higher resolution level compared to bulk sequencing techniques. In this study, a SC approach was used to analyze intratumoral heterogeneity in four patients with MDS del(5q). Single CD34+CD117+CD45+CD19- bone marrow hematopoietic stem progenitor cells were isolated using the C1 system (Fluidigm) from diagnosis or before receiving any treatment and from available follow-up samples. Selected somatic alterations were further analyzed in SC by high-throughput qPCR (Biomark HD, Fluidigm) using specific TaqMan assays. A median of 175 cells per sample were analyzed. Inferred clonal architectures were relatively simple and either linear or branching. Similar to previous studies based on bulk sequencing to infer clonal architecture, we were able to observe that an ancestral event in one patient can appear as a secondary hit in another one, thus reflecting the high intratumoral heterogeneity in MDS del(5q) and the importance of patient-specific molecular characterization.

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.

Cancer progression and the invisible phase of metastatic colonization

Metastatic dissemination occurs very early in the malignant progression of a cancer but the clinical manifestation of metastases often takes years. In recent decades, 5-year survival of patients with many solid cancers has increased due to earlier detection, local disease control and adjuvant therapies. As a consequence, we are confronted with an increase in late relapses as more antiproliferative cancer therapies prolong disease courses, raising questions about how cancer cells survive, evolve or stop growing and finally expand during periods of clinical latency. I argue here that the understanding of early metastasis formation, particularly of the currently invisible phase of metastatic colonization, will be essential for the next stage in adjuvant therapy development that reliably prevents metachronous metastasis.

Genomic Analysis of Hematopoietic Stem Cell at the Single-Cell Level: Optimization of Cell Fixation and Whole Genome Amplification (WGA) Protocol

Single-cell genomics has become the method of choice for the study of heterogeneous cell populations and represents an elective application in defining the architecture and clonal evolution in hematological neoplasms. Reconstructing the clonal evolution of a neoplastic population therefore represents the main way to understand more deeply the pathogenesis of the neoplasm, but it is also a potential tool to understand the evolution of the tumor population with respect to its response to therapy. Pre-analytical phase for single-cell genomics analysis is crucial to obtain a cell population suitable for single-cell sorting, and whole genome amplification is required to obtain the necessary amount of DNA from a single cell in order to proceed with sequencing. Here, we evaluated the impact of different methods of cellular immunostaining, fixation and whole genome amplification on the efficiency and yield of single-cell sequencing.

Towards multi-omics characterization of tumor heterogeneity: a comprehensive review of statistical and machine learning approaches

The multi-omics molecular characterization of cancer opened a new horizon for our understanding of cancer biology and therapeutic strategies. However, a tumor biopsy comprises diverse types of cells limited not only to cancerous cells but also to tumor microenvironmental cells and adjacent normal cells. This heterogeneity is a major confounding factor that hampers a robust and reproducible bioinformatic analysis for biomarker identification using multi-omics profiles. Besides, the heterogeneity itself has been recognized over the years for its significant prognostic values in some cancer types, thus offering another promising avenue for therapeutic intervention. A number of computational approaches to unravel such heterogeneity from high-throughput molecular profiles of a tumor sample have been proposed, but most of them rely on the data from an individual omics layer. Since the heterogeneity of cells is widely distributed across multi-omics layers, methods based on an individual layer can only partially characterize the heterogeneous admixture of cells. To help facilitate further development of the methodologies that synchronously account for several multi-omics profiles, we wrote a comprehensive review of diverse approaches to characterize tumor heterogeneity based on three different omics layers: genome, epigenome and transcriptome. As a result, this review can be useful for the analysis of multi-omics profiles produced by many large-scale consortia. Contact:sunkim.bioinfo@snu.ac.kr.

Low levels of AMPK promote epithelial-mesenchymal transition in lung cancer primarily through HDAC4- and HDAC5-mediated metabolic reprogramming

AMP-activated protein kinase (AMPK) serves as a "supermetabolic regulator" that helps maintain cellular energy homeostasis. However, the role of AMPK in glucose metabolism reprogramming in lung cancer remains unclear. Here, our study shows that low AMPK expression correlates with metastasis and clinicopathologic parameters of non-small-cell lung cancer. Low AMPK significantly enhances the Warburg effect in HBE and A549 cells, which in turn induces the expression of mesenchymal markers and enhances their invasion and migration. At the mechanistic level, low AMPK up-regulates HK2 expression and glycolysis levels through HDAC4 and HDAC5. Collectively, our findings demonstrate that low AMPK-induced metabolism can promote epithelial-mesenchymal transition progression in normal bronchial epithelial cells and lung cancer cells, and increase the risk for tumour metastasis.

Inference of clonal selection in cancer populations using single-cell sequencing data

Summary

Intra-tumor heterogeneity is one of the major factors influencing cancer progression and treatment outcome. However, evolutionary dynamics of cancer clone populations remain poorly understood. Quantification of clonal selection and inference of fitness landscapes of tumors is a key step to understanding evolutionary mechanisms driving cancer. These problems could be addressed using single-cell sequencing (scSeq), which provides an unprecedented insight into intra-tumor heterogeneity allowing to study and quantify selective advantages of individual clones. Here, we present Single Cell Inference of FItness Landscape (SCIFIL), a computational tool for inference of fitness landscapes of heterogeneous cancer clone populations from scSeq data. SCIFIL allows to estimate maximum likelihood fitnesses of clone variants, measure their selective advantages and order of appearance by fitting an evolutionary model into the tumor phylogeny. We demonstrate the accuracy our approach, and show how it could be applied to experimental tumor data to study clonal selection and infer evolutionary history. SCIFIL can be used to provide new insight into the evolutionary dynamics of cancer.

Availability and implementation

Its source code is available at https://github.com/compbel/SCIFIL.

Tumor Copy Number Deconvolution Integrating Bulk and Single-Cell Sequencing Data

Characterizing intratumor heterogeneity (ITH) is crucial to understanding cancer development, but it is hampered by limits of available data sources. Bulk DNA sequencing is the most common technology to assess ITH, but involves the analysis of a mixture of many genetically distinct cells in each sample, which must then be computationally deconvolved. Single-cell sequencing is a promising alternative, but its limitations-for example, high noise, difficulty scaling to large populations, technical artifacts, and large data sets-have so far made it impractical for studying cohorts of sufficient size to identify statistically robust features of tumor evolution. We have developed strategies for deconvolution and tumor phylogenetics combining limited amounts of bulk and single-cell data to gain some advantages of single-cell resolution with much lower cost, with specific focus on deconvolving genomic copy number data. We developed a mixed membership model for clonal deconvolution via non-negative matrix factorization balancing deconvolution quality with similarity to single-cell samples via an associated efficient coordinate descent algorithm. We then improve on that algorithm by integrating deconvolution with clonal phylogeny inference, using a mixed integer linear programming model to incorporate a minimum evolution phylogenetic tree cost in the problem objective. We demonstrate the effectiveness of these methods on semisimulated data of known ground truth, showing improved deconvolution accuracy relative to bulk data alone.