Inferring parsimonious migration histories for metastatic cancers

Clonal Lineage Tracing with Somatic Delivery of Recordable Barcodes Reveals Migration Histories of Metastatic Prostate Cancer

The patterns by which primary tumors spread to metastatic sites remain poorly understood. Here, we define patterns of metastatic seeding in prostate cancer using a novel injection-based mouse model-EvoCaP (Evolution in Cancer of the Prostate), featuring aggressive metastatic cancer to bone, liver, lungs, and lymph nodes. To define migration histories between primary and metastatic sites, we used our EvoTraceR pipeline to track distinct tumor clones containing recordable barcodes. We detected widespread intratumoral heterogeneity from the primary tumor in metastatic seeding, with few clonal populations instigating most migration. Metastasis-to-metastasis seeding was uncommon, as most cells remained confined within the tissue. Migration patterns in our model were congruent with human prostate cancer seeding topologies. Our findings support the view of metastatic prostate cancer as a systemic disease driven by waves of aggressive clones expanding their niche, infrequently overcoming constraints that otherwise keep them confined in the primary or metastatic site. Significance: Defining the kinetics of prostate cancer metastasis is critical for developing novel therapeutic strategies. This study uses CRISPR/Cas9-based barcoding technology to accurately define tumor clonal patterns and routes of migration in a novel somatically engineered mouse model (EvoCaP) that recapitulates human prostate cancer using an in-house developed analytical pipeline (EvoTraceR).

Defining heritability, plasticity, and transition dynamics of cellular phenotypes in somatic evolution

Single-cell sequencing has characterized cell state heterogeneity across diverse healthy and malignant tissues. However, the plasticity or heritability of these cell states remains largely unknown. To address this, we introduce PATH (phylogenetic analysis of trait heritability), a framework to quantify cell state heritability versus plasticity and infer cell state transition and proliferation dynamics from single-cell lineage tracing data. Applying PATH to a mouse model of pancreatic cancer, we observed heritability at the ends of the epithelial-to-mesenchymal transition spectrum, with higher plasticity at more intermediate states. In primary glioblastoma, we identified bidirectional transitions between stem- and mesenchymal-like cells, which use the astrocyte-like state as an intermediary. Finally, we reconstructed a phylogeny from single-cell whole-genome sequencing in B cell acute lymphoblastic leukemia and delineated the heritability of B cell differentiation states linked with genetic drivers. Altogether, PATH replaces qualitative conceptions of plasticity with quantitative measures, offering a framework to study somatic evolution.

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.

Inferring cancer type-specific patterns of metastatic spread

The metastatic spread of a cancer can be reconstructed from DNA sequencing of primary and metastatic tumours, but doing so requires solving a challenging combinatorial optimization problem. This problem often has multiple solutions that cannot be distinguished based on current maximum parsimony principles alone. Current algorithms use ad hoc criteria to select among these solutions, and decide, a priori, what patterns of metastatic spread are more likely, which is itself a key question posed by studies of metastasis seeking to use these tools. Here we introduce Metient, a freely available open-source tool which proposes multiple possible hypotheses of metastatic spread in a cohort of patients and rescores these hypotheses using independent data on genetic distance of metastasizing clones and organotropism. Metient is more accurate and is up to 50x faster than current state-of-the-art. Given a cohort of patients, Metient can calibrate its parsimony criteria, thereby identifying shared patterns of metastatic dissemination in the cohort. Reanalyzing metastasis in 169 patients based on 490 tumors, Metient automatically identifies cancer type-specific trends of metastatic dissemination in melanoma, high-risk neuroblastoma and non-small cell lung cancer. Metient's reconstructions usually agree with semi-manual expert analysis, however, in many patients, Metient identifies more plausible migration histories than experts, and further finds that polyclonal seeding of metastases is more common than previously reported. By removing the need for hard constraints on what patterns of metastatic spread are most likely, Metient introduces a way to further our understanding of cancer type-specific metastatic spread.

Enforcing Temporal Consistency in Migration History Inference

In addition to undergoing evolution, members of biological populations may also migrate between locations. Examples include the spread of tumor cells from the primary tumor to distant metastases or the spread of pathogens from one host to another. One may represent migration histories by assigning a location label to each vertex of a given phylogenetic tree such that an edge connecting vertices with distinct locations represents a migration. Some biological populations undergo comigration, a phenomenon where multiple taxa from distinct lineages simultaneously comigrate from one location to another. In this work, we show that a previous problem statement for inferring migration histories that are parsimonious in terms of migrations and comigrations may lead to temporally inconsistent solutions. To remedy this deficiency, we introduce precise definitions of temporal consistency of comigrations in a phylogenetic tree, leading to three successive problems. First, we formulate the temporally consistent comigration problem to check if a set of comigrations is temporally consistent and provide a linear time algorithm for solving this problem. Second, we formulate the parsimonious consistent comigrations (PCC) problem, which aims to find comigrations given a location labeling of a phylogenetic tree. We show that PCC is NP-hard. Third, we formulate the parsimonious consistent comigration history (PCCH) problem, which infers the migration history given a phylogenetic tree and locations of its extant vertices only. We show that PCCH is NP-hard as well. On the positive side, we propose integer linear programming models to solve the PCC and PCCH problems. We demonstrate our algorithms on simulated and real data.

Decoding the interplay between genetic and non-genetic drivers of metastasis

Metastasis is a multistep process by which cancer cells break away from their original location and spread to distant organs, and is responsible for the vast majority of cancer-related deaths. Preventing early metastatic dissemination would revolutionize the ability to fight cancer. Unfortunately, the relatively poor understanding of the molecular underpinnings of metastasis has hampered the development of effective anti-metastatic drugs. Although it is now accepted that disseminating tumour cells need to acquire multiple competencies to face the many obstacles they encounter before reaching their metastatic site(s), whether these competencies are acquired through an accumulation of metastasis-specific genetic alterations and/or non-genetic events is often debated. Here we review a growing body of literature highlighting the importance of both genetic and non-genetic reprogramming events during the metastatic cascade, and discuss how genetic and non-genetic processes act in concert to confer metastatic competencies. We also describe how recent technological advances, and in particular the advent of single-cell multi-omics and barcoding approaches, will help to better elucidate the cross-talk between genetic and non-genetic mechanisms of metastasis and ultimately inform innovative paths for the early detection and interception of this lethal process.

Maximum Likelihood Inference of Time-scaled Cell Lineage Trees with Mixed-type Missing Data

Recent dynamic lineage tracing technologies combine CRISPR-based genome editing with single-cell sequencing to track cell divisions during development. A key computational problem in dynamic lineage tracing is to infer a cell lineage tree from the measured CRISPR-induced mutations. Three features of dynamic lineage tracing data distinguish this problem from standard phylogenetic tree inference. First, the CRISPR-editing process modifies a genomic location exactly once. This non-modifiable property is not well described by the time-reversible models commonly used in phylogenetics. Second, as a consequence of non-modifiability, the number of mutations per time unit decreases over time. Third, CRISPR-based genome-editing and single-cell sequencing results in high rates of both heritable and non-heritable (dropout) missing data. To model these features, we introduce the Probabilistic Mixed-type Missing (PMM) model. We describe an algorithm, LAML (Lineage Analysis via Maximum Likelihood), to search for the maximum likelihood (ML) tree under the PMM model. LAML combines an Expectation Maximization (EM) algorithm with a heuristic tree search to jointly estimate tree topology, branch lengths and missing data parameters. We derive a closed-form solution for the M-step in the case of no heritable missing data, and a block coordinate ascent approach in the general case which is more efficient than the standard General Time Reversible (GTR) phylogenetic model. On simulated data, LAML infers more accurate tree topologies and branch lengths than existing methods, with greater advantages on datasets with higher ratios of heritable to non-heritable missing data. We show that LAML provides unbiased time-scaled estimates of branch lengths. In contrast, we demonstrate that maximum parsimony methods for lineage tracing data not only underestimate branch lengths, but also yield branch lengths which are not proportional to time, due to the nonlinear decay in the number of mutations on branches further from the root. On lineage tracing data from a mouse model of lung adenocarcinoma, we show that LAML infers phylogenetic distances that are more concordant with gene expression data compared to distances derived from maximum parsimony. The LAML tree topology is more plausible than existing published trees, with fewer total cell migrations between distant metastases and fewer reseeding events where cells migrate back to the primary tumor. Crucially, we identify three distinct time epochs of metastasis progression, which includes a burst of metastasis events to various anatomical sites during a single month.

Loss of Pip4k2c confers liver-metastatic organotropism through insulin-dependent PI3K-AKT pathway activation

Liver metastasis (LM) confers poor survival and therapy resistance across cancer types, but the mechanisms of liver-metastatic organotropism remain unknown. Here, through in vivo CRISPR-Cas9 screens, we found that Pip4k2c loss conferred LM but had no impact on lung metastasis or primary tumor growth. Pip4k2c-deficient cells were hypersensitized to insulin-mediated PI3K/AKT signaling and exploited the insulin-rich liver milieu for organ-specific metastasis. We observed concordant changes in PIP4K2C expression and distinct metabolic changes in 3,511 patient melanomas, including primary tumors, LMs and lung metastases. We found that systemic PI3K inhibition exacerbated LM burden in mice injected with Pip4k2c-deficient cancer cells through host-mediated increase in hepatic insulin levels; however, this circuit could be broken by concurrent administration of an SGLT2 inhibitor or feeding of a ketogenic diet. Thus, this work demonstrates a rare example of metastatic organotropism through co-optation of physiological metabolic cues and proposes therapeutic avenues to counteract these mechanisms.

Startle: A star homoplasy approach for CRISPR-Cas9 lineage tracing

CRISPR-Cas9-based genome editing combined with single-cell sequencing enables the tracing of the history of cell divisions, or cellular lineage, in tissues and whole organisms. Although standard phylogenetic approaches may be applied to reconstruct cellular lineage trees from this data, the unique features of the CRISPR-Cas9 editing process motivate the development of specialized models that describe the evolution of CRISPR-Cas9-induced mutations. Here, we introduce the "star homoplasy" evolutionary model that constrains a phylogenetic character to mutate at most once along a lineage, capturing the "non-modifiability" property of CRISPR-Cas9 mutations. We derive a combinatorial characterization of star homoplasy phylogenies and use this characterization to develop an algorithm, "Startle", that computes a maximum parsimony star homoplasy phylogeny. We demonstrate that Startle infers more accurate phylogenies on simulated lineage tracing data compared with existing methods and finds parsimonious phylogenies with fewer metastatic migrations on lineage tracing data from mouse metastatic lung adenocarcinoma.

TRIBAL: Tree Inference of B cell Clonal Lineages

B cells are a critical component of the adaptive immune system, responsible for producing antibodies that help protect the body from infections and foreign substances. Single cell RNA-sequencing (scRNA-seq) has allowed for both profiling of B cell receptor (BCR) sequences and gene expression. However, understanding the adaptive and evolutionary mechanisms of B cells in response to specific stimuli remains a significant challenge in the field of immunology. We introduce a new method, TRIBAL, which aims to infer the evolutionary history of clonally related B cells from scRNA-seq data. The key insight of TRIBAL is that inclusion of isotype data into the B cell lineage inference problem is valuable for reducing phylogenetic uncertainty that arises when only considering the receptor sequences. Consequently, the TRIBAL inferred B cell lineage trees jointly capture the somatic mutations introduced to the B cell receptor during affinity maturation and isotype transitions during class switch recombination. In addition, TRIBAL infers isotype transition probabilities that are valuable for gaining insight into the dynamics of class switching. Via in silico experiments, we demonstrate that TRIBAL infers isotype transition probabilities with the ability to distinguish between direct versus sequential switching in a B cell population. This results in more accurate B cell lineage trees and corresponding ancestral sequence and class switch reconstruction compared to competing methods. Using real-world scRNA-seq datasets, we show that TRIBAL recapitulates expected biological trends in a model affinity maturation system. Furthermore, the B cell lineage trees inferred by TRIBAL were equally plausible for the BCR sequences as those inferred by competing methods but yielded lower entropic partitions for the isotypes of the sequenced B cell. Thus, our method holds the potential to further advance our understanding of vaccine responses, disease progression, and the identification of therapeutic antibodies.

Integrated Molecular Characterization of Intraductal Papillary Mucinous Neoplasms: An NCI Cancer Moonshot Precancer Atlas Pilot Project

Intraductal papillary mucinous neoplasms (IPMN) are cystic precursor lesions to pancreatic ductal adenocarcinoma (PDAC). IPMNs undergo multistep progression from low-grade (LG) to high-grade (HG) dysplasia, culminating in invasive neoplasia. While patterns of IPMN progression have been analyzed using multiregion sequencing for somatic mutations, there is no integrated assessment of molecular events, including copy-number alterations (CNA) and transcriptional changes that accompany IPMN progression. We performed laser capture microdissection on surgically resected IPMNs of varying grades of histologic dysplasia obtained from 23 patients, followed by whole-exome and whole-transcriptome sequencing. Overall, HG IPMNs displayed a significantly greater aneuploidy score than LG lesions, with chromosome 1q amplification being associated with HG progression and with cases that harbored co-occurring PDAC. Furthermore, the combined assessment of single-nucleotide variants (SNV) and CNAs identified both linear and branched evolutionary trajectories, underscoring the heterogeneity in the progression of LG lesions to HG and PDAC. At the transcriptome level, upregulation of MYC-regulated targets and downregulation of transcripts associated with the MHC class I antigen presentation machinery as well as pathways related to glycosylation were a common feature of progression to HG. In addition, the established PDAC transcriptional subtypes (basal-like and classical) were readily apparent within IPMNs. Taken together, this work emphasizes the role of 1q copy-number amplification as a putative biomarker of high-risk IPMNs, underscores the importance of immune evasion even in noninvasive precursor lesions, and reinforces that evolutionary pathways in IPMNs are heterogenous, comprised of both SNV and CNA-driven events.

SIGNIFICANCE

Integrated molecular analysis of genomic and transcriptomic alterations in the multistep progression of IPMNs, which are bona fide precursors of pancreatic cancer, identifies features associated with progression of low-risk lesions to high-risk lesions and cancer, which might enable patient stratification and cancer interception strategies.

Genomic and evolutionary characteristics of metastatic gastric cancer by routes

BACKGROUND

In gastric cancer (GC) patients, metastatic progression through the lymphatic, hematogenous, peritoneal, and ovarian routes, is the ultimate cause of death. However, the genomic and evolutionary characteristics of metastatic GC have not been widely evaluated.

METHODS

Whole-exome sequencing data were analyzed for 99 primary and paired metastatic gastric cancers from 15 patients who underwent gastrectomy and metastasectomy.

RESULTS

Hematogenous metastatic tumors were associated with increased chromosomal instability and de novo gain/amplification in cancer driver genes, whereas peritoneal/ovarian metastasis was linked to sustained chromosomal stability and de novo somatic mutations in driver genes. The genomic distance of the hematogenous and peritoneal metastatic tumors was found to be closer to the primary tumors than lymph node (LN) metastasis, while ovarian metastasis was closer to LN and peritoneal metastasis than the primary tumor. Two migration patterns for metastatic GCs were identified; branched and diaspora. Both molecular subtypes of the metastatic tumors, rather than the primary tumor, and their migration patterns were related to patient survival.

CONCLUSIONS

Genomic characteristics of metastatic gastric cancer is distinctive by routes and associated with patients' prognosis along with genomic evolution pattenrs, indicating that both primary and metastatic gastric cancers require genomic evaluation.

Single-cell methylation sequencing data reveal succinct metastatic migration histories and tumor progression models

Recent studies exploring the impact of methylation in tumor evolution suggest that although the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Because changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single-cell methylation lineage tree and for jointly identifying lineage-informative CpG sites that harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole-genome sequencing data of multiregionally sampled tumor cells from nine metastatic colorectal cancer patients, as well as multiregionally sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We show that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and with more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Insights into the metastatic cascade through research autopsies

Metastasis is a complex process and the leading cause of cancer-related death globally. Recent studies have demonstrated that genomic sequencing data from paired primary and metastatic tumours can be used to trace the evolutionary origins of cells responsible for metastasis. This approach has yielded new insights into the genomic alterations that engender metastatic potential, and the mechanisms by which cancer spreads. Given that the reliability of these approaches is contingent upon how representative the samples are of primary and metastatic tumour heterogeneity, we review insights from studies that have reconstructed the evolution of metastasis within the context of their cohorts and designs. We discuss the role of research autopsies in achieving the comprehensive sampling necessary to advance the current understanding of metastasis.

Bootstrap confidence for molecular evolutionary estimates from tumor bulk sequencing data

Bulk sequencing is commonly used to characterize the genetic diversity of cancer cell populations in tumors and the evolutionary relationships of cancer clones. However, bulk sequencing produces aggregate information on nucleotide variants and their sample frequencies, necessitating computational methods to predict distinct clone sequences and their frequencies within a sample. Interestingly, no methods are available to measure the statistical confidence in the variants assigned to inferred clones. We introduce a bootstrap resampling approach that combines clone prediction and statistical confidence calculation for every variant assignment. Analysis of computer-simulated datasets showed the bootstrap approach to work well in assessing the reliability of predicted clones as well downstream inferences using the predicted clones (e.g., mapping metastatic migration paths). We found that only a fraction of inferences have good bootstrap support, which means that many inferences are tentative for real data. Using the bootstrap approach, we analyzed empirical datasets from metastatic cancers and placed bootstrap confidence on the estimated number of mutations involved in cell migration events. We found that the numbers of driver mutations involved in metastatic cell migration events sourced from primary tumors are similar to those where metastatic tumors are the source of new metastases. So, mutations with driver potential seem to keep arising during metastasis. The bootstrap approach developed in this study is implemented in software available at https://github.com/SayakaMiura/CloneFinderPlus.

The evolution of non-small cell lung cancer metastases in TRACERx

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.

Testing for Phylogenetic Signal in Single-Cell RNA-Seq Data

Phylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing or single-cell DNA sequencing and are based on calling copy number alterations and single nucleotide variants (SNVs). Single-cell RNA sequencing (scRNA-seq) is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this article, we demonstrate for the first time that scRNA-seq data contain sufficient evolutionary signal and can also be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer.

Orchestration of Collective Migration and Metastasis by Tumor Cell Clusters

Metastatic dissemination has lethal consequences for cancer patients. Accruing evidence supports the hypothesis that tumor cells can migrate and metastasize as clusters of cells while maintaining contacts with one another. Collective metastasis enables tumor cells to colonize secondary sites more efficiently, resist cell death, and evade the immune system. On the other hand, tumor cell clusters face unique challenges for dissemination particularly during systemic dissemination. Here, we review recent progress toward understanding how tumor cell clusters overcome these disadvantages as well as mechanisms they utilize to gain advantages throughout the metastatic process. We consider useful models for studying collective metastasis and reflect on how the study of collective metastasis suggests new opportunities for eradicating and preventing metastatic disease.

Single-cell dissection of remodeled inflammatory ecosystem in primary and metastatic gallbladder carcinoma

Gallbladder carcinoma (GBC) is the most common biliary tract malignancy with the lowest survival rate, primarily arising from chronic inflammation. To better characterize the progression from inflammation to cancer to metastasis, we performed single-cell RNA sequencing across samples of 6 chronic cholecystitis, 12 treatment-naive GBCs, and 6 matched metastases. Benign epithelial cells from inflamed gallbladders displayed resting, immune-regulating, and gastrointestinal metaplastic phenotypes. A small amount of PLA2G2A⁺ epithelial cells with copy number variation were identified from a histologically benign sample. We validated significant overexpression of PLA2G2A across in situ GBCs, together with increased proliferation and cancer stemness in PLA2G2A-overexpressing GBC cells, indicating an important role for PLA2G2A during early carcinogenesis. Malignant epithelial cells displayed pervasive cancer hallmarks and cellular plasticity, differentiating into metaplastic, inflammatory, and mesenchymal subtypes with distinct transcriptomic, genomic, and prognostic patterns. Chronic cholecystitis led to an adapted microenvironment characterized by MDSC-like macrophages, CD8⁺ T_RM cells, and CCL2⁺ immunity-regulating fibroblasts. By contrast, GBC instigated an aggressive and immunosuppressive microenvironment, featured by tumor-associated macrophages, Treg cells, CD8⁺ T_EX cells, and STMN1⁺ tumor-promoting fibroblasts. Single-cell and bulk RNA-seq profiles consistently showed a more suppressive immune milieu for GBCs with inflammatory epithelial signatures, coupled with strengthened epithelial-immune crosstalk. We further pinpointed a subset of senescence-like fibroblasts (FN1⁺TGM2⁺) preferentially enriched in metastatic lesions, which promoted GBC migration and invasion via their secretory phenotype. Collectively, this study provides comprehensive insights into epithelial and microenvironmental reprogramming throughout cholecystitis-propelled carcinogenesis and metastasis, laying a new foundation for the precision therapy of GBC.

Reconstructing cancer phylogenies using Pairtree, a clone tree reconstruction algorithm

Pairtree is a clone tree reconstruction algorithm that uses somatic point mutations to build clone trees describing the evolutionary history of individual cancers. Using the Pairtree software package, we describe steps to preprocess somatic mutation data, cluster mutations into subclones, search for clone trees, and visualize clone trees. Pairtree builds clone trees using up to 100 samples from a single cancer with at least 30 subclonal populations. For complete details on the use and execution of this protocol, please refer to Wintersinger et al. (2022).

Clone Phylogenetics Reveals Metastatic Tumor Migrations, Maps, and Models

Simple Summary

Metastasis is the spread of cancer cells across organs and is a major cause of cancer mortality. Analysis of tumor sequencing data provides a means toward the reconstruction of routes of metastatic cell migrations. Our reconstructions demonstrated that many metastases were likely seeded from pre-existing metastasis of primary tumors. Additionally, multiple clone exchanges between tumor sites were common. In conclusion, the pattern of cancer cell migrations is often complex and is highly variable among patients.

Abstract

Dispersal routes of metastatic cells are not medically detected or even visible. A molecular evolutionary analysis of tumor variation provides a way to retrospectively infer metastatic migration histories and answer questions such as whether the majority of metastases are seeded from clones within primary tumors or seeded from clones within pre-existing metastases, as well as whether the evolution of metastases is generally consistent with any proposed models. We seek answers to these fundamental questions through a systematic patient-centric retrospective analysis that maps the dynamic evolutionary history of tumor cell migrations in many cancers. We analyzed tumor genetic heterogeneity in 51 cancer patients and found that most metastatic migration histories were best described by a hybrid of models of metastatic tumor evolution. Synthesizing across metastatic migration histories, we found new tumor seedings arising from clones of pre-existing metastases as often as they arose from clones from primary tumors. There were also many clone exchanges between the source and recipient tumors. Therefore, a molecular phylogenetic analysis of tumor variation provides a retrospective glimpse into general patterns of metastatic migration histories in cancer patients.

A phylogenetic approach to study the evolution of somatic mutational processes in cancer

Cancer cell genomes change continuously due to mutations, and mutational processes change over time in patients, leaving dynamic signatures in the accumulated genomic variation in tumors. Many computational methods detect the relative activities of known mutation signatures. However, these methods may produce erroneous signatures when applied to individual branches in cancer cell phylogenies. Here, we show that the inference of branch-specific mutational signatures can be improved through a joint analysis of the collections of mutations mapped on proximal branches of the cancer cell phylogeny. This approach reduces the false-positive discovery rate of branch-specific signatures and can sometimes detect faint signatures. An analysis of empirical data from 61 lung cancer patients supports trends based on computer-simulated datasets for which the correct signatures are known. In lung cancer somatic variation, we detect a decreasing trend of smoking-related mutational processes over time and an increasing influence of APOBEC mutational processes as the tumor evolution progresses. These analyses also reveal patterns of conservation and divergence of mutational processes in cell lineages within patients.

Model-based inference of metastatic seeding rates in de novo metastatic breast cancer reveals the impact of secondary seeding and molecular subtype

We present a stochastic network model of metastasis spread for de novo metastatic breast cancer, composed of tumor to metastasis (primary seeding) and metastasis to metastasis spread (secondary seeding), parameterized using the SEER (Surveillance, Epidemiology, and End Results) database. The model provides a quantification of tumor cell dissemination rates between the tumor and metastasis sites. These rates were used to estimate the probability of developing a metastasis for untreated patients. The model was validated using tenfold cross-validation. We also investigated the effect of HER2 (Human Epidermal Growth Factor Receptor 2) status, estrogen receptor (ER) status and progesterone receptor (PR) status on the probability of metastatic spread. We found that dissemination rate through secondary seeding is up to 300 times higher than through primary seeding. Hormone receptor positivity promotes seeding to the bone and reduces seeding to the lungs and primary seeding to the liver, while HER2 expression increases dissemination to the bone, lungs and primary seeding to the liver. Secondary seeding from the lungs to the liver seems to be hormone receptor-independent, while that from the lungs to the brain appears HER2-independent.

Deciphering clonal dynamics and metastatic routines in a rare patient of synchronous triple-primary tumors and multiple metastases with MPTevol

Multiple primary tumor (MPT) is a special and rare cancer type, defined as more than two primary tumors presenting at the diagnosis in a single patient. The molecular characteristics and tumorigenesis of MPT remain unclear due to insufficient approaches. Here, we present MPTevol, a practical computational framework for comprehensively exploring the MPT from multiregion sequencing (MRS) experiments. To verify the utility of MPTevol, we performed whole-exome MRS for 33 samples of a rare patient with triple-primary tumors and three metastatic sites and systematically investigated clonal dynamics and metastatic routines. MPTevol assists in comparing genomic profiles across samples, detecting clonal evolutionary history and metastatic routines and quantifying the metastatic history. All triple-primary tumors were independent origins and their genomic characteristics were consistent with corresponding sporadic tumors, strongly supporting their independent tumorigenesis. We further showed two independent early monoclonal seeding events for the metastases in the ovary and uterus. We revealed that two ovarian metastases were disseminated from the same subclone of the primary tumor through undergoing whole-genome doubling processes, suggesting metastases-to-metastases seeding occurred when tumors had similar microenvironments. Surprisingly, according to the metastasis timing model of MPTevol, we found that primary tumors of about 0.058–0.124 cm diameter have been disseminating to distant organs, which is much earlier than conventional clinical views. We developed MPT-specialized analysis framework MPTevol and demonstrated its utility in explicitly resolving clonal evolutionary history and metastatic seeding routines with a rare MPT case. MPTevol is implemented in R and is available at https://github.com/qingjian1991/MPTevol under the GPL v3 license.

SimSCSnTree: a simulator of single-cell DNA sequencing data

SUMMARY

We report on a new single-cell DNA sequence simulator, SimSCSnTree, which generates an evolutionary tree of cells and evolves single nucleotide variants (SNVs) and copy number aberrations (CNAs) along its branches. Data generated by the simulator can be used to benchmark tools for single-cell genomic analyses, particularly in cancer where SNVs and CNAs are ubiquitous.

AVAILABILITY AND IMPLEMENTATION

SimSCSnTree is now on BioConda and also is freely available for download at https://github.com/compbiofan/SimSCSnTree.git with detailed documentation.

The translational challenges of precision oncology

The translational challenges in the field of precision oncology are in part related to the biological complexity and diversity of this disease. Technological advances in genomics have facilitated large sequencing efforts and discoveries that have further supported this notion. In this review, we reflect on the impact of these discoveries on our understanding of several concepts: cancer initiation, cancer prevention, early detection, adjuvant therapy and minimal residual disease monitoring, cancer drug resistance, and cancer evolution in metastasis. We discuss key areas of focus for improving cancer outcomes, from biological insights to clinical application, and suggest where the development of these technologies will lead us. Finally, we discuss practical challenges to the wider adoption of molecular profiling in the clinic and the need for robust translational infrastructure.

The Genetic Evolution of Metastasis

Cancer is an evolutionary process that is characterized by the emergence of multiple genetically distinct populations or clones within the primary tumor. Intratumor heterogeneity provides a substrate for the selection of adaptive clones, such as those that lead to metastasis. Comparative molecular studies of primary tumors and metastases have identified distinct genomic features associated with the development of metastases. In this review, we discuss how these insights could inform clinical decision-making and uncover rational antimetastasis treatment strategies.

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.

Genomic evolution and diverse models of systemic metastases in colorectal cancer

OBJECTIVE

The systemic spread of colorectal cancer (CRC) is dominated by the portal system and exhibits diverse patterns of metastasis without systematical genomic investigation. Here, we evaluated the genomic evolution of CRC with multiorgan metastases using multiregion sequencing.

DESIGN

Whole-exome sequencing was performed on multiple regions (n=74) of matched primary tumour, adjacent non-cancerous mucosa, liver metastasis and lung metastasis from six patients with CRC. Phylogenetic reconstruction and evolutionary analyses were used to investigate the metastatic seeding pattern and clonal origin. Recurrent driver gene mutations were analysed across patients and validated in two independent cohorts. Metastatic assays were performed to examine the effect of the novel driver gene on the malignant behaviour of CRC cells.

RESULTS

Based on the migration patterns and clonal origins, three models were revealed (sequential, branch-off and diaspora), which not only supported the anatomic assumption that CRC cells spread to lung after clonally expanding in the liver, but also illustrated the direct seeding of extrahepatic metastases from primary tumours independently. Unlike other cancer types, polyphyletic seeding occurs in CRC, which may result in late metastases with intermetastatic driver gene heterogeneity. In cases with rapid dissemination, we found recurrent trunk loss-of-function mutations in ZFP36L2, which is enriched in metastatic CRC and associated with poor overall survival. CRISPR/Cas9-mediated knockout of ZFP36L2 enhances the metastatic potential of CRC cells.

CONCLUSION

Our results provide genomic evidence for metastatic evolution and indicate that biopsy/sequencing of metastases may be considered for patients with CRC with multiorgan or late postoperative metastasis.

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.

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.

Evolution of metastasis: new tools and insights

Metastasis is an evolutionary process occurring across multiple organs and timescales. Due to its continuous and dynamic nature, this multifaceted process has been challenging to investigate and remains incompletely understood, in part due to the lack of tools capable of probing genomic evolution at high enough resolution. However, technological advances in genetic sequencing and editing have provided new and powerful methods to refine our understanding of the complex series of events that lead to metastatic dissemination. In this review, we summarize the latest genetic and lineage-tracing approaches developed to unravel the genetic evolution of metastasis. The findings that have emerged have enhanced our comprehension of the mechanistic trajectories and timescales of metastasis and could provide new strategies for therapy.

Tumors Are Evolutionary Island-Like Ecosystems

Integration of ecological and evolutionary features has begun to understand the interplay of tumor heterogeneity, microenvironment, and metastatic potential. Developing a theoretical framework is intrinsic to deciphering tumors' tremendous spatial and longitudinal genetic variation patterns in patients. Here, we propose that tumors can be considered evolutionary island-like ecosystems, that is, isolated systems that undergo evolutionary and spatiotemporal dynamic processes that shape tumor microenvironments and drive the migration of cancer cells. We examine attributes of insular systems and causes of insularity, such as physical distance and connectivity. These properties modulate migration rates of cancer cells through processes causing spatial and temporal isolation of the organs and tissues functioning as a supply of cancer cells for new colonizations. We discuss hypotheses, predictions, and limitations of tumors as islands analogy. We present emerging evidence of tumor insularity in different cancer types and discuss their relevance to the islands model. We suggest that the engagement of tumor insularity into conceptual and mathematical models holds promise to illuminate cancer evolution, tumor heterogeneity, and metastatic potential of cells.

Migrations of cancer cells through the lens of phylogenetic biogeography

Malignant cells leave their initial tumor of growth and disperse to other tissues to form metastases. Dispersals also occur in nature when individuals in a population migrate from their area of origin to colonize other habitats. In cancer, phylogenetic biogeography is concerned with the source and trajectory of cell movements. We examine the suitability of primary features of organismal biogeography, including genetic diversification, dispersal, extinction, vicariance, and founder effects, to describe and reconstruct clone migration events among tumors. We used computer-simulated data to compare fits of seven biogeographic models and evaluate models’ performance in clone migration reconstruction. Models considering founder effects and dispersals were often better fit for the clone phylogenetic patterns, especially for polyclonal seeding and reseeding of metastases. However, simpler biogeographic models produced more accurate estimates of cell migration histories. Analyses of empirical datasets of basal-like breast cancer had model fits consistent with the patterns seen in the analysis of computer-simulated datasets. Our analyses reveal the powers and pitfalls of biogeographic models for modeling and inferring clone migration histories using tumor genome variation data. We conclude that the principles of molecular evolution and organismal biogeography are useful in these endeavors but that the available models and methods need to be applied judiciously.

The somatic molecular evolution of cancer: Mutation, selection, and epistasis

Cancer progression has been attributed to somatic changes in single-nucleotide variants, copy-number aberrations, loss of heterozygosity, chromosomal instability, epistatic interactions, and the tumor microenvironment. It is not entirely clear which of these changes are essential and which are ancillary to cancer. The dynamic nature of cancer evolution in a patient can be illuminated using several concepts and tools from classical evolutionary biology. Neutral mutation rates in cancer cells are calculable from genomic data such as synonymous mutations, and selective pressures are calculable from rates of fixation occurring beyond the expectation by neutral mutation and drift. However, these cancer effect sizes of mutations are complicated by epistatic interactions that can determine the likely sequence of gene mutations. In turn, longitudinal phylogenetic analyses of somatic cancer progression offer an opportunity to identify key moments in cancer evolution, relating the timing of driver mutations to corresponding landmarks in the clinical timeline. These analyses reveal temporal aspects of genetic and phenotypic change during tumorigenesis and across clinical timescales. Using a related framework, clonal deconvolution, physical locations of clones, and their phylogenetic relations can be used to infer tumor migration histories. Additionally, genetic interactions with the tumor microenvironment can be analyzed with longstanding approaches applied to organismal genotype-by-environment interactions. Fitness landscapes for cancer evolution relating to genotype, phenotype, and environment could enable more accurate, personalized therapeutic strategies. An understanding of the trajectories underlying the evolution of neoplasms, primary, and metastatic tumors promises fundamental advances toward accurate and personalized predictions of therapeutic response.

Genetic and non-genetic clonal diversity in cancer evolution

The observation and analysis of intra-tumour heterogeneity (ITH), particularly in genomic studies, has advanced our understanding of the evolutionary forces that shape cancer growth and development. However, only a subset of the variation observed in a single tumour will have an impact on cancer evolution, highlighting the need to distinguish between functional and non-functional ITH. Emerging studies highlight a role for the cancer epigenome, transcriptome and immune microenvironment in functional ITH. Here, we consider the importance of both genetic and non-genetic ITH and their role in tumour evolution, and present the rationale for a broad research focus beyond the cancer genome. Systems-biology analytical approaches will be necessary to outline the scale and importance of functional ITH. By allowing a deeper understanding of tumour evolution this will, in time, encourage development of novel therapies and improve outcomes for patients.

Drivers of dynamic intratumor heterogeneity and phenotypic plasticity

Cancer is a clonal disease, i.e., all tumor cells within a malignant lesion trace their lineage back to a precursor somatic cell that acquired oncogenic mutations during development and aging. And yet, those tumor cells tend to have genetic and nongenetic variations among themselves-which is denoted as intratumor heterogeneity. Although some of these variations are inconsequential, others tend to contribute to cell state transition and phenotypic heterogeneity, providing a substrate for somatic evolution. Tumor cell phenotypes can dynamically change under the influence of genetic mutations, epigenetic modifications, and microenvironmental contexts. Although epigenetic and microenvironmental changes are adaptive, genetic mutations are usually considered permanent. Emerging reports suggest that certain classes of genetic alterations show extensive reversibility in tumors in clinically relevant timescales, contributing as major drivers of dynamic intratumor heterogeneity and phenotypic plasticity. Dynamic heterogeneity and phenotypic plasticity can confer resistance to treatment, promote metastasis, and enhance evolvability in cancer. Here, we first highlight recent efforts to characterize intratumor heterogeneity at genetic, epigenetic, and microenvironmental levels. We then discuss phenotypic plasticity and cell state transition by tumor cells, under the influence of genetic and nongenetic determinants and their clinical significance in classification of tumors and therapeutic decision-making.

PathFinder: Bayesian inference of clone migration histories in cancer

SUMMARY

Metastases cause a vast majority of cancer morbidity and mortality. Metastatic clones are formed by dispersal of cancer cells to secondary tissues, and are not medically detected or visible until later stages of cancer development. Clone phylogenies within patients provide a means of tracing the otherwise inaccessible dynamic history of migrations of cancer cells. Here, we present a new Bayesian approach, PathFinder, for reconstructing the routes of cancer cell migrations. PathFinder uses the clone phylogeny, the number of mutational differences among clones, and the information on the presence and absence of observed clones in primary and metastatic tumors. By analyzing simulated datasets, we found that PathFinder performes well in reconstructing clone migrations from the primary tumor to new metastases as well as between metastases. It was more challenging to trace migrations from metastases back to primary tumors. We found that a vast majority of errors can be corrected by sampling more clones per tumor, and by increasing the number of genetic variants assayed per clone. We also identified situations in which phylogenetic approaches alone are not sufficient to reconstruct migration routes.In conclusion, we anticipate that the use of PathFinder will enable a more reliable inference of migration histories and their posterior probabilities, which is required to assess the relative preponderance of seeding of new metastasis by clones from primary tumors and/or existing metastases.

AVAILABILITY AND IMPLEMENTATION

PathFinder is available on the web at https://github.com/SayakaMiura/PathFinder.

Elements and evolutionary determinants of genomic divergence between paired primary and metastatic tumors

Can metastatic-primary (M-P) genomic divergence measured from next generation sequencing reveal the natural history of metastatic dissemination? This remains an open question of utmost importance in facilitating a deeper understanding of metastatic progression, and thereby, improving its prevention. Here, we utilize mathematical and computational modeling to tackle this question as well as to provide a framework that illuminates the fundamental elements and evolutionary determinants of M-P divergence. Our framework facilitates the integration of sequencing detectability of somatic variants, and hence, paves the way towards bridging the measurable between-tumor heterogeneity with analytical modeling and interpretability. We show that the number of somatic variants of the metastatic seeding cell that are experimentally undetectable in the primary tumor, can be characterized as the path of the phylogenetic tree from the last appearing variant of the seeding cell back to the most recent detectable variant. We find that the expected length of this path is principally determined by the decay in detectability of the variants along the seeding cell's lineage; and thus, exhibits a significant dependence on the underlying tumor growth dynamics. A striking implication of this fact, is that dissemination from an advanced detectable subclone of the primary tumor can lead to an abrupt drop in the expected measurable M-P divergence, thereby breaking the previously assumed monotonic relation between seeding time and M-P divergence. This is emphatically verified by our single cell-based spatial tumor growth simulation, where we find that M-P divergence exhibits a non-monotonic relationship with seeding time when the primary tumor grows under branched and linear evolution. On the other hand, a monotonic relationship holds when we condition on the dynamics of progressive diversification, or by restricting the seeding cells to always originate from undetectable subclones. Our results highlight the fact that a precise understanding of tumor growth dynamics is the sine qua non for exploiting M-P divergence to reconstruct the chronology of metastatic dissemination. The quantitative models presented here enable further careful evaluation of M-P divergence in association with crucial evolutionary and sequencing parameters.

Single-cell lineages reveal the rates, routes, and drivers of metastasis in cancer xenografts

Detailed phylogenies of tumor populations can recount the history and chronology of critical events during cancer progression, such as metastatic dissemination. We applied a Cas9-based, single-cell lineage tracer to study the rates, routes, and drivers of metastasis in a lung cancer xenograft mouse model. We report deeply resolved phylogenies for tens of thousands of cancer cells traced over months of growth and dissemination. This revealed stark heterogeneity in metastatic capacity, arising from preexisting and heritable differences in gene expression. We demonstrate that these identified genes can drive invasiveness and uncovered an unanticipated suppressive role for KRT17 We also show that metastases disseminated via multidirectional tissue routes and complex seeding topologies. Overall, we demonstrate the power of tracing cancer progression at subclonal resolution and vast scale.

A practical guide to cancer subclonal reconstruction from DNA sequencing

Subclonal reconstruction from bulk tumor DNA sequencing has become a pillar of cancer evolution studies, providing insight into the clonality and relative ordering of mutations and mutational processes. We provide an outline of the complex computational approaches used for subclonal reconstruction from single and multiple tumor samples. We identify the underlying assumptions and uncertainties in each step and suggest best practices for analysis and quality assessment. This guide provides a pragmatic resource for the growing user community of subclonal reconstruction methods.

Inference of mutability landscapes of tumors from single cell sequencing data

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

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

CONTACT

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

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Prostate cancer evolution from multilineage primary to single lineage metastases with implications for liquid biopsy

The evolutionary progression from primary to metastatic prostate cancer is largely uncharted, and the implications for liquid biopsy are unexplored. We infer detailed reconstructions of tumor phylogenies in ten prostate cancer patients with fatal disease, and investigate them in conjunction with histopathology and tumor DNA extracted from blood and cerebrospinal fluid. Substantial evolution occurs within the prostate, resulting in branching into multiple spatially intermixed lineages. One dominant lineage emerges that initiates and drives systemic metastasis, where polyclonal seeding between sites is common. Routes to metastasis differ between patients, and likely genetic drivers of metastasis distinguish the metastatic lineage from the lineage that remains confined to the prostate within each patient. Body fluids capture features of the dominant lineage, and subclonal expansions that occur in the metastatic phase are non-uniformly represented. Cerebrospinal fluid analysis reveals lineages not detected in blood-borne DNA, suggesting possible clinical utility.

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

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

Reconstruction of clone- and haplotype-specific cancer genome karyotypes from bulk tumor samples

Many cancer genomes are extensively rearranged with aberrant chromosomal karyotypes. Deriving these karyotypes from high-throughput DNA sequencing of bulk tumor samples is complicated because most tumors are a heterogeneous mixture of normal cells and subpopulations of cancer cells, or clones, that harbor distinct somatic mutations. We introduce a new algorithm, Reconstructing Cancer Karyotypes (RCK), to reconstruct haplotype-specific karyotypes of one or more rearranged cancer genomes from DNA sequencing data from a bulk tumor sample. RCK leverages evolutionary constraints on the somatic mutational process in cancer to reduce ambiguity in the deconvolution of admixed sequencing data into multiple haplotype-specific cancer karyotypes. RCK models mixtures containing an arbitrary number of derived genomes and allows the incorporation of information both from short-read and long-read DNA sequencing technologies. We compare RCK to existing approaches on 17 primary and metastatic prostate cancer samples. We find that RCK infers cancer karyotypes that better explain the DNA sequencing data and conform to a reasonable evolutionary model. RCK's reconstructions of clone- and haplotype-specific karyotypes will aid further studies of the role of intra-tumor heterogeneity in cancer development and response to treatment. RCK is freely available as open source software.

Accurate quantification of copy-number aberrations and whole-genome duplications in multi-sample tumor sequencing data

Copy-number aberrations (CNAs) and whole-genome duplications (WGDs) are frequent somatic mutations in cancer but their quantification from DNA sequencing of bulk tumor samples is challenging. Standard methods for CNA inference analyze tumor samples individually; however, DNA sequencing of multiple samples from a cancer patient has recently become more common. We introduce HATCHet (Holistic Allele-specific Tumor Copy-number Heterogeneity), an algorithm that infers allele- and clone-specific CNAs and WGDs jointly across multiple tumor samples from the same patient. We show that HATCHet outperforms current state-of-the-art methods on multi-sample DNA sequencing data that we simulate using MASCoTE (Multiple Allele-specific Simulation of Copy-number Tumor Evolution). Applying HATCHet to 84 tumor samples from 14 prostate and pancreas cancer patients, we identify subclonal CNAs and WGDs that are more plausible than previously published analyses and more consistent with somatic single-nucleotide variants (SNVs) and small indels in the same samples.

Multi-site clonality analysis uncovers pervasive heterogeneity across melanoma metastases

Metastatic melanoma carries a poor prognosis despite modern systemic therapies. Understanding the evolution of the disease could help inform patient management. Through whole-genome sequencing of 13 melanoma metastases sampled at autopsy from a treatment naïve patient and by leveraging the analytical power of multi-sample analyses, we reveal evidence of diversification among metastatic lineages. UV-induced mutations dominate the trunk, whereas APOBEC-associated mutations are found in the branches of the evolutionary tree. Multi-sample analyses from a further seven patients confirmed that lineage diversification was pervasive, representing an important mode of melanoma dissemination. Our analyses demonstrate that joint analysis of cancer cell fraction estimates across multiple metastases can uncover previously unrecognised levels of tumour heterogeneity and highlight the limitations of inferring heterogeneity from a single biopsy.