A new method to accurately identify single nucleotide variants using small FFPE breast samples.

Most tissue collections of neoplasms are composed of formalin-fixed and paraffin-embedded (FFPE) excised tumor samples used for routine diagnostics. DNA sequencing is becoming increasingly important in cancer research and clinical management; however it is difficult to accurately sequence DNA from FFPE samples. We developed and validated a new bioinformatic pipeline to use existing variant-calling strategies to robustly identify somatic single nucleotide variants (SNVs) from whole exome sequencing using small amounts of DNA extracted from archival FFPE samples of breast cancers. We optimized this strategy using 28 pairs of technical replicates. After optimization, the mean similarity between replicates increased 5-fold, reaching 88% (range 0-100%), with a mean of 21.4 SNVs (range 1-68) per sample, representing a markedly superior performance to existing tools. We found that the SNV-identification accuracy declined when there was less than 40 ng of DNA available and that insertion-deletion variant calls are less reliable than single base substitutions. As the first application of the new algorithm, we compared samples of ductal carcinoma in situ of the breast to their adjacent invasive ductal carcinoma samples. We observed an increased number of mutations (paired-samples sign test, P < 0.05), and a higher genetic divergence in the invasive samples (paired-samples sign test, P < 0.01). Our method provides a significant improvement in detecting SNVs in FFPE samples over previous approaches.

RecPhyloXML: a format for reconciled gene trees.

Motivation: A reconciliation is an annotation of the nodes of a gene tree with evolutionary events-for example, speciation, gene duplication, transfer, loss, etc.-along with a mapping onto a species tree. Many algorithms and software produce or use reconciliations but often using different reconciliation formats, regarding the type of events considered or whether the species tree is dated or not. This complicates the comparison and communication between different programs. Results: Here, we gather a consortium of software developers in gene tree species tree reconciliation to propose and endorse a format that aims to promote an integrative-albeit flexible-specification of phylogenetic reconciliations. This format, named recPhyloXML, is accompanied by several tools such as a reconciled tree visualizer and conversion utilities. Availability and implementation: http://phylariane.univ-lyon1.fr/recphyloxml/.

SimPhy: Phylogenomic Simulation of Gene, Locus, and Species Trees.

We present a fast and flexible software package--SimPhy--for the simulation of multiple gene families evolving under incomplete lineage sorting, gene duplication and loss, horizontal gene transfer--all three potentially leading to species tree/gene tree discordance--and gene conversion. SimPhy implements a hierarchical phylogenetic model in which the evolution of species, locus, and gene trees is governed by global and local parameters (e.g., genome-wide, species-specific, locus-specific), that can be fixed or be sampled from a priori statistical distributions. SimPhy also incorporates comprehensive models of substitution rate variation among lineages (uncorrelated relaxed clocks) and the capability of simulating partitioned nucleotide, codon, and protein multilocus sequence alignments under a plethora of substitution models using the program INDELible. We validate SimPhy's output using theoretical expectations and other programs, and show that it scales extremely well with complex models and/or large trees, being an order of magnitude faster than the most similar program (DLCoal-Sim). In addition, we demonstrate how SimPhy can be useful to understand interactions among different evolutionary processes, conducting a simulation study to characterize the systematic overestimation of the duplication time when using standard reconciliation methods. SimPhy is available at https://github.com/adamallo/SimPhy, where users can find the source code, precompiled executables, a detailed manual and example cases.

A Bayesian Supertree Model for Genome-Wide Species Tree Reconstruction.

Current phylogenomic data sets highlight the need for species tree methods able to deal with several sources of gene tree/species tree incongruence. At the same time, we need to make most use of all available data. Most species tree methods deal with single processes of phylogenetic discordance, namely, gene duplication and loss, incomplete lineage sorting (ILS) or horizontal gene transfer. In this manuscript, we address the problem of species tree inference from multilocus, genome-wide data sets regardless of the presence of gene duplication and loss and ILS therefore without the need to identify orthologs or to use a single individual per species. We do this by extending the idea of Maximum Likelihood (ML) supertrees to a hierarchical Bayesian model where several sources of gene tree/species tree disagreement can be accounted for in a modular manner. We implemented this model in a computer program called guenomu whose inputs are posterior distributions of unrooted gene tree topologies for multiple gene families, and whose output is the posterior distribution of rooted species tree topologies. We conducted extensive simulations to evaluate the performance of our approach in comparison with other species tree approaches able to deal with more than one leaf from the same species. Our method ranked best under simulated data sets, in spite of ignoring branch lengths, and performed well on empirical data, as well as being fast enough to analyze relatively large data sets. Our Bayesian supertree method was also very successful in obtaining better estimates of gene trees, by reducing the uncertainty in their distributions. In addition, our results show that under complex simulation scenarios, gene tree parsimony is also a competitive approach once we consider its speed, in contrast to more sophisticated models.

Diversity and distribution of unicellular opisthokonts along the European coast analysed using high-throughput sequencing.

The opisthokonts are one of the major super groups of eukaryotes. It comprises two major clades: (i) the Metazoa and their unicellular relatives and (ii) the Fungi and their unicellular relatives. There is, however, little knowledge of the role of opisthokont microbes in many natural environments, especially among non-metazoan and non-fungal opisthokonts. Here, we begin to address this gap by analysing high-throughput 18S rDNA and 18S rRNA sequencing data from different European coastal sites, sampled at different size fractions and depths. In particular, we analyse the diversity and abundance of choanoflagellates, filastereans, ichthyosporeans, nucleariids, corallochytreans and their related lineages. Our results show the great diversity of choanoflagellates in coastal waters as well as a relevant representation of the ichthyosporeans and the uncultured marine opisthokonts (MAOP). Furthermore, we describe a new lineage of marine fonticulids (MAFO) that appears to be abundant in sediments. Taken together, our work points to a greater potential ecological role for unicellular opisthokonts than previously appreciated in marine environments, both in water column and sediments, and also provides evidence of novel opisthokont phylogenetic lineages. This study highlights the importance of high-throughput sequencing approaches to unravel the diversity and distribution of both known and novel eukaryotic lineages.

Stronger Together: Cancer Clones Cooperate to Alleviate Growth Barriers in Critical Cancer Progression Transitions.

Hershey and colleagues recently showed how clones in a triple-negative breast cancer cell line cooperate for their mutual fitness benefit. In this system, clones exchange soluble metabolites to increase their in vitro growth rate at low population densities, therefore mitigating the documented growth barrier that reduces individual fitness in small tumor cell populations (Allee effect). Such cooperation could aid important transitions in cancer progression in which cancer cell populations are small, like invasion or metastasis. Using orthotopic transplantation, the authors demonstrate that this cooperation is functional in one such transition in vivo, increasing the metastatic load and number of metastases, which are usually polyclonal. Together, these findings highlight the need to consider ecologic interactions to properly understand tumor growth dynamics, and how they complement the standing evolutionary model of cancer progression in our quest to understand and treat cancer.

Growth Dynamics of Ductal Carcinoma in Situ Recapitulate Normal Breast Development.

Ductal carcinoma in situ (DCIS) and invasive breast cancer share many morphologic, proteomic, and genomic alterations. Yet in contrast to invasive cancer, many DCIS tumors do not progress and may remain indolent over decades. To better understand the heterogenous nature of this disease, we reconstructed the growth dynamics of 18 DCIS tumors based on the geo-spatial distribution of their somatic mutations. The somatic mutation topographies revealed that DCIS is multiclonal and consists of spatially discontinuous subclonal lesions. Here we show that this pattern of spread is consistent with a new 'Comet' model of DCIS tumorigenesis, whereby multiple subclones arise early and nucleate the buds of the growing tumor. The discontinuous, multiclonal growth of the Comet model is analogous to the branching morphogenesis of normal breast development that governs the rapid expansion of the mammary epithelium during puberty. The branching morphogenesis-like dynamics of the proposed Comet model diverges from the canonical model of clonal evolution, and better explains observed genomic spatial data. Importantly, the Comet model allows for the clinically relevant scenario of extensive DCIS spread, without being subjected to the selective pressures of subclone competition that promote the emergence of increasingly invasive phenotypes. As such, the normal cell movement inferred during DCIS growth provides a new explanation for the limited risk of progression in DCIS and adds biologic rationale for ongoing clinical efforts to reduce DCIS overtreatment.

Cancer Prevalence Across Vertebrates.

Cancer is pervasive across multicellular species. Are there any patterns that can explain differences in cancer prevalence across species? Using 16,049 necropsy records for 292 species spanning three clades (amphibians, sauropsids and mammals) we found that neoplasia and malignancy prevalence increases with adult weight and decreases with gestation time, contrary to Peto’s Paradox. Evolution of cancer susceptibility appears to have undergone sudden shifts followed by stabilizing selection. Outliers for neoplasia prevalence include the common porpoise (<1.3%), the Rodrigues fruit bat (<1.6%) the black-footed penguin (<0.4%), ferrets (63%) and opossums (35%). Discovering why some species have particularly high or low levels of cancer may lead to a better understanding of cancer syndromes and novel strategies for the management and prevention of cancer.

Cancer Prevalence Across Vertebrates.

Cancer is pervasive across multicellular species, but what explains differences in cancer prevalence across species? Using 16,049 necropsy records for 292 species spanning three clades (amphibians, sauropsids and mammals) we found that neoplasia and malignancy prevalence increases with adult weight (contrary to Peto's Paradox) and somatic mutation rate, but decreases with gestation time. Evolution of cancer susceptibility appears to have undergone sudden shifts followed by stabilizing selection. Outliers for neoplasia prevalence include the common porpoise (<1.3%), the Rodrigues fruit bat (<1.6%) the black-footed penguin (<0.4%), ferrets (63%) and opossums (35%). Discovering why some species have particularly high or low levels of cancer may lead to a better understanding of cancer syndromes and novel strategies for the management and prevention of cancer.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer.

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

Minimal barriers to invasion during human colorectal tumor growth.

Intra-tumoral heterogeneity (ITH) could represent clonal evolution where subclones with greater fitness confer more malignant phenotypes and invasion constitutes an evolutionary bottleneck. Alternatively, ITH could represent branching evolution with invasion of multiple subclones. The two models respectively predict a hierarchy of subclones arranged by phenotype, or multiple subclones with shared phenotypes. We delineate these modes of invasion by merging ancestral, topographic, and phenotypic information from 12 human colorectal tumors (11 carcinomas, 1 adenoma) obtained through saturation microdissection of 325 small tumor regions. The majority of subclones (29/46, 60%) share superficial and invasive phenotypes. Of 11 carcinomas, 9 show evidence of multiclonal invasion, and invasive and metastatic subclones arise early along the ancestral trees. Early multiclonal invasion in the majority of these tumors indicates the expansion of co-evolving subclones with similar malignant potential in absence of late bottlenecks and suggests that barriers to invasion are minimal during colorectal cancer growth.

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.

Natural Selection in Cancer Biology: From Molecular Snowflakes to Trait Hallmarks.

Evolution by natural selection is the conceptual foundation for nearly every branch of biology and increasingly also for biomedicine and medical research. In cancer biology, evolution explains how populations of cells in tumors change over time. It is a fundamental question whether this evolutionary process is driven primarily by natural selection and adaptation or by other evolutionary processes such as founder effects and drift. In cancer biology, as in organismal evolutionary biology, there is controversy about this question and also about the use of adaptation through natural selection as a guiding framework for research. In this review, we discuss the differences and similarities between evolution among somatic cells versus evolution among organisms. We review what is known about the parameters and rate of evolution in neoplasms, as well as evidence for adaptation. We conclude that adaptation is a useful framework that accurately explains the defining characteristics of cancer. Further, convergent evolution through natural selection provides the only satisfying explanation both for how a group of diverse pathologies have enough in common to usefully share the descriptive label of "cancer" and for why this convergent condition becomes life-threatening.

Multilocus inference of species trees and DNA barcoding.

The unprecedented amount of data resulting from next-generation sequencing has opened a new era in phylogenetic estimation. Although large datasets should, in theory, increase phylogenetic resolution, massive, multilocus datasets have uncovered a great deal of phylogenetic incongruence among different genomic regions, due both to stochastic error and to the action of different evolutionary process such as incomplete lineage sorting, gene duplication and loss and horizontal gene transfer. This incongruence violates one of the fundamental assumptions of the DNA barcoding approach, which assumes that gene history and species history are identical. In this review, we explain some of the most important challenges we will have to face to reconstruct the history of species, and the advantages and disadvantages of different strategies for the phylogenetic analysis of multilocus data. In particular, we describe the evolutionary events that can generate species tree-gene tree discordance, compare the most popular methods for species tree reconstruction, highlight the challenges we need to face when using them and discuss their potential utility in barcoding. Current barcoding methods sacrifice a great amount of statistical power by only considering one locus, and a transition to multilocus barcodes would not only improve current barcoding methods, but also facilitate an eventual transition to species-tree-based barcoding strategies, which could better accommodate scenarios where the barcode gap is too small or inexistent.This article is part of the themed issue 'From DNA barcodes to biomes'.

When (distant) relatives stay too long: implications for cancer medicine.

Whole-genome analyses of human medulloblastomas show that the dominant clone at relapse is present as a rare subclone at primary diagnosis.