Multiregional Tumor Trees Are Not Phylogenies

Global lineage evolution pattern of sars-cov-2 in Africa, America, Europe, and Asia: A comparative analysis of variant clusters and their relevance across continents

Objective

The objective of this study is to provide a comparative analysis of variant clusters and their relevance across Africa, America, Europe, and Asia, in order to understand the evolutionary patterns of the virus across different regions and to inform the development of targeted interventions and genomic surveillance eforts.

Methods

The study analyzed the global lineage evolution pattern of 74, 075 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from 32 countries across four continents, focusing on variant clusters and their relevance across regions. Variants were weighted according to their hierarchical level. The correlation between variants was visualized through Dimensionality reduction analysis and Pairwise Pearson's correlation. We presented a reconstructed phylogenetic tree based on correlation analysis and variant weights.

Results

The analysis revealed that each continent had distinct variant clusters and different evolutionary patterns. The Americas had two clustered variants before lineage divergence and a downstream confluence lineage, Europe had bifurcation into two global lineages with an early occurrence of certain cluster while Asia had a downstream confluence of two large lineages diverging by two distinct clusters. Based on the cluster patterns of shared variants of the SARS-CoV-2 virus, Africa demonstrated a relatively clear distinction among three distinct regions.

Conclusions

The study provides insights into the evolutionary patterns of SARS-CoV-2 and highlights the importance of international collaboration in tracking and responding to emerging variants. The study found that the global pandemic was driven by Omicron variants that evolved with significant differences between countries and regions, and with different patterns across continents.

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.

Late-Stage Metastatic Melanoma Emerges through a Diversity of Evolutionary Pathways

Understanding the evolutionary pathways to metastasis and resistance to immune-checkpoint inhibitors (ICI) in melanoma is critical for improving outcomes. Here, we present the most comprehensive intrapatient metastatic melanoma dataset assembled to date as part of the Posthumous Evaluation of Advanced Cancer Environment (PEACE) research autopsy program, including 222 exome sequencing, 493 panel-sequenced, 161 RNA sequencing, and 22 single-cell whole-genome sequencing samples from 14 ICI-treated patients. We observed frequent whole-genome doubling and widespread loss of heterozygosity, often involving antigen-presentation machinery. We found KIT extrachromosomal DNA may have contributed to the lack of response to KIT inhibitors of a KIT-driven melanoma. At the lesion-level, MYC amplifications were enriched in ICI nonresponders. Single-cell sequencing revealed polyclonal seeding of metastases originating from clones with different ploidy in one patient. Finally, we observed that brain metastases that diverged early in molecular evolution emerge late in disease. Overall, our study illustrates the diverse evolutionary landscape of advanced melanoma.

SIGNIFICANCE

Despite treatment advances, melanoma remains a deadly disease at stage IV. Through research autopsy and dense sampling of metastases combined with extensive multiomic profiling, our study elucidates the many mechanisms that melanomas use to evade treatment and the immune system, whether through mutations, widespread copy-number alterations, or extrachromosomal DNA. See related commentary by Shain, p. 1294. This article is highlighted in the In This Issue feature, p. 1275.

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.

State-dependent evolutionary models reveal modes of solid tumour growth

Spatial properties of tumour growth have profound implications for cancer progression, therapeutic resistance and metastasis. Yet, how spatial position governs tumour cell division remains difficult to evaluate in clinical tumours. Here, we demonstrate that faster division on the tumour periphery leaves characteristic genetic patterns, which become evident when a phylogenetic tree is reconstructed from spatially sampled cells. Namely, rapidly dividing peripheral lineages branch more extensively and acquire more mutations than slower-dividing centre lineages. We develop a Bayesian state-dependent evolutionary phylodynamic model (SDevo) that quantifies these patterns to infer the differential division rates between peripheral and central cells. We demonstrate that this approach accurately infers spatially varying birth rates of simulated tumours across a range of growth conditions and sampling strategies. We then show that SDevo outperforms state-of-the-art, non-cancer multi-state phylodynamic methods that ignore differential sequence evolution. Finally, we apply SDevo to single-time-point, multi-region sequencing data from clinical hepatocellular carcinomas and find evidence of a three- to six-times-higher division rate on the tumour edge. With the increasing availability of high-resolution, multi-region sequencing, we anticipate that SDevo will be useful in interrogating spatial growth restrictions and could be extended to model non-spatial factors that influence tumour progression.

Accurate tumor clonal structures require single-cell analysis

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

Estimation of cancer cell fractions and clone trees from multi-region sequencing of tumors

MOTIVATION

Multi-region sequencing of solid tumors can improve our understanding of intratumor subclonal diversity and the evolutionary history of mutational events. Due to uncertainty in clonal composition and the multitude of possible ancestral relationships between clones, elucidating the most probable relationships from bulk tumor sequencing poses statistical and computational challenges.

RESULTS

We developed a Bayesian hierarchical model called PICTograph to model uncertainty in assigning mutations to subclones, to enable posterior distributions of cancer cell fractions, and to visualize the most probable ancestral relationships between subclones. Compared to available methods, PICTograph provided more consistent and accurate estimates of cancer cell fractions and improved tree inference over a range of simulated clonal diversity. Application of PICTograph to multi-region whole exome sequencing of tumors from individuals with pancreatic cancer precursor lesions confirmed known early-occurring mutations and indicated substantial molecular diversity, including 6-12 distinct subclones and intra-sample mixing of subclones. Using ensemble-based visualizations, we highlight highly probable evolutionary relationships recovered in multiple models. PICTograph provides a useful approximation to evolutionary inference from cross-sectional multi-region sequencing, particularly for complex cases.

AVAILABILITY

https://github.com/KarchinLab/pictograph.

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.

Conditional prediction of consecutive tumor evolution using cancer progression models: What genotype comes next?

Accurate prediction of tumor progression is key for adaptive therapy and precision medicine. Cancer progression models (CPMs) can be used to infer dependencies in mutation accumulation from cross-sectional data and provide predictions of tumor progression paths. However, their performance when predicting complete evolutionary trajectories is limited by violations of assumptions and the size of available data sets. Instead of predicting full tumor progression paths, here we focus on short-term predictions, more relevant for diagnostic and therapeutic purposes. We examine whether five distinct CPMs can be used to answer the question "Given that a genotype with n mutations has been observed, what genotype with n + 1 mutations is next in the path of tumor progression?" or, shortly, "What genotype comes next?". Using simulated data we find that under specific combinations of genotype and fitness landscape characteristics CPMs can provide predictions of short-term evolution that closely match the true probabilities, and that some genotype characteristics can be much more relevant than global features. Application of these methods to 25 cancer data sets shows that their use is hampered by a lack of information needed to make principled decisions about method choice. Fruitful use of these methods for short-term predictions requires adapting method's use to local genotype characteristics and obtaining reliable indicators of performance; it will also be necessary to clarify the interpretation of the method's results when key assumptions do not hold.

Characterizing Intra-Tumor Heterogeneity From Somatic Mutations Without Copy-Neutral Assumption

Bulk samples of the same patient are heterogeneous in nature, comprising of different subpopulations (subclones) of cancer cells. Cells in a tumor subclone are characterized by unique mutational genotype profile. Resolving tumor heterogeneity by estimating the genotypes, cellular proportions and the number of subclones present in the tumor can help in understanding cancer progression and treatment. We present a novel method, ChaClone2, to efficiently deconvolve the observed variant allele fractions (VAFs), with consideration for possible effects from copy number aberrations at the mutation loci. Our method describes a state-space formulation of the feature allocation model, deconvolving the observed VAFs from samples of the same patient into three matrices: subclonal total and variant copy numbers for mutated genes, and proportions of subclones in each sample. We describe an efficient sequential Monte Carlo (SMC) algorithm to estimate these matrices. Extensive simulation shows that the ChaClone2 yields better accuracy when compared with other state-of-the-art methods for addressing similar problem and it offers scalability to large datasets. Also, ChaClone2 features that the model parameter estimates can be refined whenever new mutation data of freshly sequenced genomic locations are available. MATLAB code and datasets are available to download at: https://github.com/moyanre/method2.

DEVOLUTION-A method for phylogenetic reconstruction of aneuploid cancers based on multiregional genotyping data

Phylogenetic reconstruction of cancer cell populations remains challenging. There is a particular lack of tools that deconvolve clones based on copy number aberration analyses of multiple tumor biopsies separated in time and space from the same patient. This has hampered investigations of tumors rich in aneuploidy but few point mutations, as in many childhood cancers and high-risk adult cancer. Here, we present DEVOLUTION, an algorithm for subclonal deconvolution followed by phylogenetic reconstruction from bulk genotyping data. It integrates copy number and sequencing information across multiple tumor regions throughout the inference process, provided that the mutated clone fraction for each mutation is known. We validate DEVOLUTION on data from 56 pediatric tumors comprising 253 tumor biopsies and show a robust performance on simulations of bulk genotyping data. We also benchmark DEVOLUTION to similar bioinformatic tools using an external dataset. DEVOLUTION holds the potential to facilitate insights into the development, progression, and response to treatment, particularly in tumors with high burden of chromosomal copy number alterations.

Clonal evolution in liver cancer at single-cell and single-variant resolution

Genetic heterogeneity of tumor is closely related to its clonal evolution, phenotypic diversity and treatment resistance, and such heterogeneity has only been characterized at single-cell sub-chromosomal scale in liver cancer. Here we reconstructed the single-variant resolution clonal evolution in human liver cancer based on single-cell mutational profiles. The results indicated that key genetic events occurred early during tumorigenesis, and an early metastasis followed by independent evolution was observed in primary liver tumor and intrahepatic metastatic portal vein tumor thrombus. By parallel single-cell RNA-Seq, the transcriptomic phenotype of HCC was found to be related with genetic heterogeneity. For the first time we reconstructed the single-cell and single-variant clonal evolution in human liver cancer, and dissection of both genetic and phenotypic heterogeneity will facilitate better understanding of their relationship.

Quantifying the influence of mutation detection on tumour subclonal reconstruction

Whole-genome sequencing can be used to estimate subclonal populations in tumours and this intra-tumoural heterogeneity is linked to clinical outcomes. Many algorithms have been developed for subclonal reconstruction, but their variabilities and consistencies are largely unknown. We evaluate sixteen pipelines for reconstructing the evolutionary histories of 293 localized prostate cancers from single samples, and eighteen pipelines for the reconstruction of 10 tumours with multi-region sampling. We show that predictions of subclonal architecture and timing of somatic mutations vary extensively across pipelines. Pipelines show consistent types of biases, with those incorporating SomaticSniper and Battenberg preferentially predicting homogenous cancer cell populations and those using MuTect tending to predict multiple populations of cancer cells. Subclonal reconstructions using multi-region sampling confirm that single-sample reconstructions systematically underestimate intra-tumoural heterogeneity, predicting on average fewer than half of the cancer cell populations identified by multi-region sequencing. Overall, these biases suggest caution in interpreting specific architectures and subclonal variants.

Helicobacter pylori diversification during chronic infection within a single host generates sub-populations with distinct phenotypes

Helicobacter pylori chronically infects the stomach of approximately half of the world's population. Manifestation of clinical diseases associated with H. pylori infection, including cancer, is driven by strain properties and host responses; and as chronic infection persists, both are subject to change. Previous studies have documented frequent and extensive within-host bacterial genetic variation. To define how within-host diversity contributes to phenotypes related to H. pylori pathogenesis, this project leverages a collection of 39 clinical isolates acquired prospectively from a single subject at two time points and from multiple gastric sites. During the six years separating collection of these isolates, this individual, initially harboring a duodenal ulcer, progressed to gastric atrophy and concomitant loss of acid secretion. Whole genome sequence analysis identified 1,767 unique single nucleotide polymorphisms (SNPs) across isolates and a nucleotide substitution rate of 1.3x10-4 substitutions/site/year. Gene ontology analysis identified cell envelope genes among the genes with excess accumulation of nonsynonymous SNPs (nSNPs). A maximum likelihood tree based on genetic similarity clusters isolates from each time point separately. Within time points, there is segregation of subgroups with phenotypic differences in bacterial morphology, ability to induce inflammatory cytokines, and mouse colonization. Higher inflammatory cytokine induction in recent isolates maps to shared polymorphisms in the Cag PAI protein, CagY, while rod morphology in a subgroup of recent isolates mapped to eight mutations in three distinct helical cell shape determining (csd) genes. The presence of subgroups with unique genetic and phenotypic properties suggest complex selective forces and multiple niches within the stomach during chronic infection.

Malignant transformation and genetic alterations are uncoupled in early colorectal cancer progression

BACKGROUND

Colorectal cancer (CRC) development is generally accepted as a sequential process, with genetic mutations determining phenotypic tumor progression. However, matching genetic profiles with histological transition requires the analyses of temporal samples from the same patient at key stages of progression.

RESULTS

Here, we compared the genetic profiles of 34 early carcinomas with their respective adenomatous precursors to assess timing and heterogeneity of driver alterations accompanying the switch from benign adenoma to malignant carcinoma. In almost half of the cases, driver mutations specific to the carcinoma stage were not observed. In samples where carcinoma-specific alterations were present, TP53 mutations and chromosome 20 copy gains commonly accompanied the switch from adenomatous tissue to carcinoma. Remarkably, 40% and 50% of high-grade adenomas shared TP53 mutations and chromosome 20 gains, respectively, with their matched carcinomas. In addition, multi-regional analyses revealed greater heterogeneity of driver mutations in adenomas compared to their matched carcinomas.

CONCLUSION

Genetic alterations in TP53 and chromosome 20 occur at the earliest histological stage in colorectal carcinomas (pTis and pT1). However, high-grade adenomas can share these alterations despite their histological distinction. Based on the well-defined sequence of CRC development, we suggest that the timing of genetic changes during neoplastic progression is frequently uncoupled from histological progression.

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.

Genetic Alterations in Preinvasive Lung Synchronous Lesions

Purpose

Despite advances in treatment, lung cancer remains the leading cause of cancer mortality. This study aimed to characterise genome-wide tumorigenesis events and to understand the hypothesis of the multistep carcinogenesis of lung adenocarcinoma (LUAD).

Materials and Methods

We conducted multiregion whole-exome sequencing of LUAD with synchronous atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ, or minimally invasive adenocarcinoma of 19 samples from three patients to characterize genome-wide tumorigenesis events and validate the hypothesis of the multistep carcinogenesis of LUAD. We identified potential pathogenic mutations preserved in preinvasive lesions and supplemented the finding by allelic variant level from RNA sequencing.

Results

Overall, independent mutational profiles were observed per patient and between patients. Some shared mutations including epidermal growth factor receptor (EGFR, p.L858R) were present across synchronous lesions.

Conclusion

Here, we show that there are driver gene mutations in AAH, and they may exacerbate as a sequence in a histological continuum, supporting the Darwinian evolution model of cancer genome. The intertumoral and intratumoral heterogeneity of synchronous LUAD implies that multi-biomarker strategies might be necessary for appropriate treatment.

A new phylogenetic protocol: dealing with model misspecification and confirmation bias in molecular phylogenetics

Molecular phylogenetics plays a key role in comparative genomics and has increasingly significant impacts on science, industry, government, public health and society. In this paper, we posit that the current phylogenetic protocol is missing two critical steps, and that their absence allows model misspecification and confirmation bias to unduly influence phylogenetic estimates. Based on the potential offered by well-established but under-used procedures, such as assessment of phylogenetic assumptions and tests of goodness of fit, we introduce a new phylogenetic protocol that will reduce confirmation bias and increase the accuracy of phylogenetic estimates.

SCARLET: Single-cell tumor phylogeny inference with copy-number constrained mutation losses

A small number of somatic mutations drive the development of cancer, but all somatic mutations are markers of the evolutionary history of a tumor. Prominent methods to construct phylogenies from single-cell sequencing data use single-nucleotide variants (SNVs) as markers but fail to adequately account for copy-number aberrations (CNAs), which can overlap SNVs and result in SNV losses. Here, we introduce SCARLET, an algorithm that infers tumor phylogenies from single-cell DNA sequencing data while accounting for both CNA-driven loss of SNVs and sequencing errors. SCARLET outperforms existing methods on simulated data, with more accurate inference of the order in which mutations were acquired and the mutations present in individual cells. Using a single-cell dataset from a patient with colorectal cancer, SCARLET constructs a tumor phylogeny that is consistent with the observed CNAs and suggests an alternate origin for the patient's metastases. SCARLET is available at: github.com/raphael-group/scarlet.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Analysis of Telomere Lengths in p53 Signatures and Incidental Serous Tubal Intraepithelial Carcinomas Without Concurrent Ovarian Cancer

Telomere alterations represent one of the major molecular changes in the development of human cancer. We have previously reported that telomere lengths in most serous tubal intraepithelial carcinomas (STIC) are shorter than they are in ovarian high-grade serous carcinomas (HGSC) or in normal-appearing fallopian tube epithelium from the same patients. However, it remains critical to determine if similar telomere alterations occur in TP53-mutated but histologically unremarkable "p53 signature" lesions, as well as incidental STICs without concurrent HGSC. In this study, we quantitatively measured telomere lengths by performing telomere-specific fluorescence in situ hybridization in conjunction with p53 immunolabeling in 15 p53 signatures and 30 incidental STICs without concurrent HGSC. We compared these new results with our previous data in paired STICs and concurrent HGSCs. We found that most p53 signatures (80%) and incidental STICs without HGSC (77%) exhibited significant telomere shortening compared with adjacent normal-appearing fallopian tube epithelium (P<0.01). Interestingly, however, p53 signatures and incidental STICs without HGSC displayed longer telomeres and less cell-to-cell telomere length heterogeneity than STICs associated with HGSC (P<0.001). These findings indicate that telomere shortening occurs in p53 signatures, the earliest precancer lesion. Moreover, incidental STICs without concurrent HGSC are indeed similar to p53 signatures as they have less telomere shortening and less cell-to-cell telomere length heterogeneity than STICs associated with HGSC.

A Glioblastoma Genomics Primer for Clinicians

New discoveries in Glioblastoma (GBM) biology have been made using genomics data. Genomic markers are routinely integrated into clinical neurosurgical practice. In this manuscript, we review the fundamentals of genomics such as the differences between first, second, and third generation sequencing technology. We also review the impact of single cell genomics in understanding the complex heterogenous GBM microenvironment. Finally, we will discuss advances in epigenetics that have lent insights into treatment resistance. The integration of genomics into neuro-oncology clinical practice is routine and will continue to expand with the expansion of precision of medicine. We provide a primer for clinicians.

Extensive Clonal Branching Shapes the Evolutionary History of High-Risk Pediatric Cancers

Darwinian evolution of tumor cells remains underexplored in childhood cancer. We here reconstruct the evolutionary histories of 56 pediatric primary tumors, including 24 neuroblastomas, 24 Wilms tumors, and 8 rhabdomyosarcomas. Whole-genome copy-number and whole-exome mutational profiling of multiple regions per tumor were performed, followed by clonal deconvolution to reconstruct a phylogenetic tree for each tumor. Overall, 88% of the tumors exhibited genetic variation among primary tumor regions. This variability typically emerged through collateral phylogenetic branching, leading to spatial variability in the distribution of more than 50% (96/173) of detected diagnostically informative genetic aberrations. Single-cell sequencing of 547 individual cancer cells from eight solid pediatric tumors confirmed branching evolution to be a fundamental underlying principle of genetic variation in all cases. Strikingly, cell-to-cell genetic diversity was almost twice as high in aggressive compared with clinically favorable tumors (median Simpson index of diversity 0.45 vs. 0.88; P = 0.029). Similarly, a comparison of multiregional sampling data from a total of 274 tumor regions showed that new phylogenetic branches emerge at a higher frequency per sample and carry a higher mutational load in high-risk than in low-risk tumors. Timelines based on spatial genetic variation showed that the mutations most influencing relapse risk occur at initiation of clonal expansion in neuroblastoma and rhabdomyosarcoma, whereas in Wilms tumor, they are late events. Thus, from an evolutionary standpoint, some high-risk childhood cancers are born bad, whereas others grow worse over time. SIGNIFICANCE: Different pediatric cancers with a high risk of relapse share a common generic pattern of extensively branching evolution of somatic mutations. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/7/1512/F1.large.jpg.

Rapid evolution and biogeographic spread in a colorectal cancer

How and when tumoral clones start spreading to surrounding and distant tissues is currently unclear. Here we leveraged a model-based evolutionary framework to investigate the demographic and biogeographic history of a colorectal cancer. Our analyses strongly support an early monoclonal metastatic colonization, followed by a rapid population expansion at both primary and secondary sites. Moreover, we infer a hematogenous metastatic spread under positive selection, plus the return of some tumoral cells from the liver back to the colon lymph nodes. This study illustrates how sophisticated techniques typical of organismal evolution can provide a detailed, quantitative picture of the complex tumoral dynamics over time and space.

The clonal origins of metastases and the timing of dissemination remains an open question for most cancer types. Using primary and metastatic samples taken from one colorectal cancer patient, Alves et al. use Bayesian phylogenetics to reconstruct the history of metastasis.

Predicting clone genotypes from tumor bulk sequencing of multiple samples

Motivation

Analyses of data generated from bulk sequencing of tumors have revealed extensive genomic heterogeneity within patients. Many computational methods have been developed to enable the inference of genotypes of tumor cell populations (clones) from bulk sequencing data. However, the relative and absolute accuracy of available computational methods in estimating clone counts and clone genotypes is not yet known.

Results

We have assessed the performance of nine methods, including eight previously-published and one new method (CloneFinder), by analyzing computer simulated datasets. CloneFinder, LICHeE, CITUP and cloneHD inferred clone genotypes with low error (<5% per clone) for a majority of datasets in which the tumor samples contained evolutionarily-related clones. Computational methods did not perform well for datasets in which tumor samples contained mixtures of clones from different clonal lineages. Generally, the number of clones was underestimated by cloneHD and overestimated by PhyloWGS, and BayClone2, Canopy and Clomial required prior information regarding the number of clones. AncesTree and Canopy did not produce results for a large number of datasets. Overall, the deconvolution of clone genotypes from single nucleotide variant (SNV) frequency differences among tumor samples remains challenging, so there is a need to develop more accurate computational methods and robust software for clone genotype inference.

Availability and implementation

CloneFinder is implemented in Python and is available from https://github.com/gstecher/CloneFinderAPI.

Supplementary information

Supplementary data are available at Bioinformatics online.

Spatially constrained tumour growth affects the patterns of clonal selection and neutral drift in cancer genomic data

Quantification of the effect of spatial tumour sampling on the patterns of mutations detected in next-generation sequencing data is largely lacking. Here we use a spatial stochastic cellular automaton model of tumour growth that accounts for somatic mutations, selection, drift and spatial constraints, to simulate multi-region sequencing data derived from spatial sampling of a neoplasm. We show that the spatial structure of a solid cancer has a major impact on the detection of clonal selection and genetic drift from both bulk and single-cell sequencing data. Our results indicate that spatial constrains can introduce significant sampling biases when performing multi-region bulk sampling and that such bias becomes a major confounding factor for the measurement of the evolutionary dynamics of human tumours. We also propose a statistical inference framework that incorporates spatial effects within a growing tumour and so represents a further step forwards in the inference of evolutionary dynamics from genomic data. Our analysis shows that measuring cancer evolution using next-generation sequencing while accounting for the numerous confounding factors remains challenging. However, mechanistic model-based approaches have the potential to capture the sources of noise and better interpret the data.

Measuring Clonal Evolution in Cancer with Genomics

Cancers originate from somatic cells in the human body that have accumulated genetic alterations. These mutations modify the phenotype of the cells, allowing them to escape the homeostatic regulation that maintains normal cell number. Viewed through the lens of evolutionary biology, the transformation of normal cells into malignant cells is evolution in action. Evolution continues throughout cancer growth, progression, treatment resistance, and disease relapse, driven by adaptation to changes in the cancer's environment, and intratumor heterogeneity is an inevitable consequence of this evolutionary process. Genomics provides a powerful means to characterize tumor evolution, enabling quantitative measurement of evolving clones across space and time. In this review, we discuss concepts and approaches to quantify and measure this evolutionary process in cancer using genomics.

Genomic landscape and evolutionary trajectories of ovarian cancer precursor lesions

The clonal relationship between ovarian high-grade serous carcinoma (HGSC) and its presumed precursor lesion, serous tubal intraepithelial carcinoma (STIC), has been reported. However, when analyzing patients with concurrent ovarian carcinoma and tubal lesion, the extensive carcinoma tissues present at diagnosis may have effaced the natural habitat of precursor clone(s), obscuring tumor clonal evolutionary history, or may have disseminated to anatomically adjacent fimbriae ends, masquerading as precursor lesions. To circumvent these limitations, we analyzed the genomic landscape of incidental tubal precursor lesions including p53 signature, dormant STIC or serous tubal intraepithelial lesion (STIL) and proliferative STIC in women without ovarian carcinoma or any cancer diagnosis using whole-exome sequencing and amplicon sequencing. In three of the four cancer-free women with multiple discrete tubal lesions we observed non-identical TP53 mutations between precursor lesions from the same individual. In one of the four women with co-existing ovarian HGSC and tubal precursor lesion we found non-identical TP53 mutations and a lack of common mutations shared between her precursor lesion and carcinoma. Analyzing the evolutionary history of multiple tubal lesions in the same four patients with concurrent ovarian carcinoma indicated distinct evolution trajectories. Collectively, the results support diverse clonal origins of tubal precursor lesions at the very early stages of tumorigenesis. Mathematical modeling based on lesion-specific proliferation rates indicated that p53 signature and dormant STIC may take a prolonged time (two decades or more) to develop into STIC, whereas STIC may progress to carcinoma in a much shorter time (6 years). The above findings may have implications for future research aimed at prevention and early detection of ovarian cancer. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Clonal evolution through genetic bottlenecks and telomere attrition: Potential threats to in vitro data reproducibility

Tissue cultures of immortalized human cells, also known as established cell lines, are broadly accessible and cost-efficient tools for biomedical research. We here review potential genetic sources of systematic error in cell line experiments due to clonal evolution in vitro. In particular, the authors highlight alterations in telomere function over prolonged culture and population bottlenecks, respectively, as two commonly overlooked phenomena that can result in significant alterations in cell line genotypes over just one or a few passages in vitro. These alterations may include changes in mutation status of oncogenes and large scale chromosomal imbalances. We introduce a simple list of factors to be avoided in order to reduce the risk of data misinterpretation due to clonal evolution, including unacknowledged in vitro selection pressures, prolonged culture per se, harsh population size reductions, experiments at early phases after establishment, and the employment of cell lines not sufficiently analyzed by high resolution genetic techniques.

The Spatiotemporal Evolution of Lymph Node Spread in Early Breast Cancer

Purpose: The most significant prognostic factor in early breast cancer is lymph node involvement. This stage between localized and systemic disease is key to understanding breast cancer progression; however, our knowledge of the evolution of lymph node malignant invasion remains limited, as most currently available data are derived from primary tumors.Experimental Design: In 11 patients with treatment-naïve node-positive early breast cancer without clinical evidence of distant metastasis, we investigated lymph node evolution using spatial multiregion sequencing (n = 78 samples) of primary and lymph node deposits and genomic profiling of matched longitudinal circulating tumor DNA (ctDNA).Results: Linear evolution from primary to lymph node was rare (1/11), whereas the majority of cases displayed either early divergence between primary and nodes (4/11) or no detectable divergence (6/11), where both primary and nodal cells belonged to a single recent expansion of a metastatic clone. Divergence of metastatic subclones was driven in part by APOBEC. Longitudinal ctDNA samples from 2 of 7 subjects with evaluable plasma taken perioperatively reflected the two major evolutionary patterns and demonstrate that private mutations can be detected even from early metastatic nodal deposits. Moreover, node removal resulted in disappearance of private lymph node mutations in ctDNA.Conclusions: This study sheds new light on a crucial evolutionary step in the natural history of breast cancer, demonstrating early establishment of axillary lymph node metastasis in a substantial proportion of patients. Clin Cancer Res; 24(19); 4763-70. ©2018 AACR.

Inferring parsimonious migration histories for metastatic cancers

Metastasis is the migration of cancerous cells from a primary tumor to other anatomical sites. Although metastasis was long thought to result from monoclonal seeding, or single cellular migrations, recent phylogenetic analyses of metastatic cancers have reported complex patterns of cellular migrations between sites, including polyclonal migrations and reseeding. However, accurate determination of migration patterns from somatic mutation data is complicated by intratumor heterogeneity and discordance between clonal lineage and cellular migration. We introduce MACHINA, a multi-objective optimization algorithm that jointly infers clonal lineages and parsimonious migration histories of metastatic cancers from DNA sequencing data. MACHINA analysis of data from multiple cancers shows that migration patterns are often not uniquely determined from sequencing data alone and that complicated migration patterns among primary tumors and metastases may be less prevalent than previously reported. MACHINA's rigorous analysis of migration histories will aid in studies of the drivers of metastasis.

Sensitivity to sequencing depth in single-cell cancer genomics

BACKGROUND

Querying cancer genomes at single-cell resolution is expected to provide a powerful framework to understand in detail the dynamics of cancer evolution. However, given the high costs currently associated with single-cell sequencing, together with the inevitable technical noise arising from single-cell genome amplification, cost-effective strategies that maximize the quality of single-cell data are critically needed. Taking advantage of previously published single-cell whole-genome and whole-exome cancer datasets, we studied the impact of sequencing depth and sampling effort towards single-cell variant detection.

METHODS

Five single-cell whole-genome and whole-exome cancer datasets were independently downscaled to 25, 10, 5, and 1× sequencing depth. For each depth level, ten technical replicates were generated, resulting in a total of 6280 single-cell BAM files. The sensitivity of variant detection, including structural and driver mutations, genotyping, clonal inference, and phylogenetic reconstruction to sequencing depth was evaluated using recent tools specifically designed for single-cell data.

RESULTS

Altogether, our results suggest that for relatively large sample sizes (25 or more cells) sequencing single tumor cells at depths > 5× does not drastically improve somatic variant discovery, characterization of clonal genotypes, or estimation of single-cell phylogenies.

CONCLUSIONS

We suggest that sequencing multiple individual tumor cells at a modest depth represents an effective alternative to explore the mutational landscape and clonal evolutionary patterns of cancer genomes.

Evolution of Barrett's esophagus through space and time at single-crypt and whole-biopsy levels

The low risk of progression of Barrett's esophagus (BE) to esophageal adenocarcinoma can lead to over-diagnosis and over-treatment of BE patients. This may be addressed through a better understanding of the dynamics surrounding BE malignant progression. Although genetic diversity has been characterized as a marker of malignant development, it is still unclear how BE arises and develops. Here we uncover the evolutionary dynamics of BE at crypt and biopsy levels in eight individuals, including four patients that experienced malignant progression. We assay eight individual crypts and the remaining epithelium by SNP array for each of 6-11 biopsies over 2 time points per patient (358 samples in total). Our results indicate that most Barrett's segments are clonal, with similar number and inferred rates of alterations observed for crypts and biopsies. Divergence correlates with geographical location, being higher near the gastro-esophageal junction. Relaxed clock analyses show that genomic instability precedes and is enhanced by genome doubling. These results shed light on the clinically relevant evolutionary dynamics of BE.

Pathological Bases and Clinical Impact of Intratumor Heterogeneity in Clear Cell Renal Cell Carcinoma

PURPOSE OF REVIEW

Intratumor heterogeneity is an inherent event in tumor development that is receiving much attention in the last years since it is responsible for most failures of current targeted therapies. The purpose of this review is to offer clinicians an updated insight of the multiple manifestations of a complex event that impacts significantly patient's life.

RECENT FINDINGS

Clear cell renal cell carcinoma is the most common renal tumor and a paradigmatic example of a heterogeneous neoplasm. Next-generation sequencing has demonstrated that intratumor heterogeneity encompasses genetic, epigenetic, and microenvironmental variability. Currently accepted protocols of tumor sampling seem insufficient in unveiling intratumor heterogeneity with reliability and need to be updated. This variability challenges the precise morphological diagnosis, its molecular characterization, and the selection of optimal personalized therapies in clear cell renal cell carcinoma, a neoplasm traditionally considered chemo- and radio-resistant. We review the state of the art of the different approaches to intratumor heterogeneity in clear cell renal cell carcinomas, from the simple morphology to the most sophisticated massive sequencing tools.