Unraveling the clonal hierarchy of somatic genomic aberrations

Defining the chronology of molecular alterations may identify milestones in carcinogenesis. To unravel the temporal evolution of aberrations from clinical tumors, we developed CLONET, which upon estimation of tumor admixture and ploidy infers the clonal hierarchy of genomic aberrations. Comparative analysis across 100 sequenced genomes from prostate, melanoma, and lung cancers established diverse evolutionary hierarchies, demonstrating the early disruption of tumor-specific pathways. The analyses highlight the diversity of clonal evolution within and across tumor types that might be informative for risk stratification and patient selection for targeted therapies. CLONET addresses heterogeneous clinical samples seen in the setting of precision medicine. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0439-6) contains supplementary material, which is available to authorized users

Clonal architecture of secondary acute myeloid leukemia defined by single-cell sequencing

Next-generation sequencing has been used to infer the clonality of heterogeneous tumor samples. These analyses yield specific predictions-the population frequency of individual clones, their genetic composition, and their evolutionary relationships-which we set out to test by sequencing individual cells from three subjects diagnosed with secondary acute myeloid leukemia, each of whom had been previously characterized by whole genome sequencing of unfractionated tumor samples. Single-cell mutation profiling strongly supported the clonal architecture implied by the analysis of bulk material. In addition, it resolved the clonal assignment of single nucleotide variants that had been initially ambiguous and identified areas of previously unappreciated complexity. Accordingly, we find that many of the key assumptions underlying the analysis of tumor clonality by deep sequencing of unfractionated material are valid. Furthermore, we illustrate a single-cell sequencing strategy for interrogating the clonal relationships among known variants that is cost-effective, scalable, and adaptable to the analysis of both hematopoietic and solid tumors, or any heterogeneous population of cells

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1 %) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones, We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.Peer reviewe

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data

RECONSTRUCTING MUTATIONAL LINEAGES IN BREAST CANCER BY MULTI-PATIENT-TARGETED SINGLE CELL DNA SEQUENCING

Triple negative breast cancer (TNBC) is an aggressive subtype of breast cancer with high rates of metastasis and recurrence, where TNBC patients have a poor 5-year survival and ~50% are non-responsive to chemotherapy. Aneuploidy is a cancer hallmark that is pervasive in over 90% of breast cancer patients and is indicative of complex genomic rearrangements that are acquired during tumor initiation. Although copy number aberrations have been extensively studied in relation to aneuploidy and TNBC initiation, little is currently known regarding the timing and impact of single nucleotide variants (SNVs) contributing to these early transformative genomic events. Paramount to novel treatment options is understanding the underlying biology of initiation in the early stages of TNBC development, where inferring clonal lineages and mutational chronologies can help characterize the order, timing, and potential impact of single point mutations, such as TP53, in association with early mechanisms of genomic instability, genome doubling, and aneuploid transformation. Single cell DNA sequencing (scDNA-seq) methods are powerful tools for profiling these mutations in cancer cells, however most genomic regions sequenced in single cells are non-informative. To overcome this issue, we developed a Multi-Patient-Targeted (MPT) scDNA-seq method. MPT involves first performing bulk exome sequencing across a cohort of cancer patients to identify somatic mutations, which are then pooled together to develop a single custom targeted panel for high-throughput scDNA-seq using a microfluidics platform. We applied MPT to profile 330 mutations across 23,500 cells from 5 TNBC patients, which showed that 3 tumors were monoclonal and 2 tumors were polyclonal. From this data, we reconstructed mutational lineages and identified early mutational and copy number events, including early TP53 mutations that occurred in all five patients. Collectively, our data suggests that MPT can overcome a major technical obstacle for studying tumor evolution using scDNA-seq by profiling information-rich mutation sites. By applying MPT to resolve the order, timing, and evolutionary relationships of early point mutations, in association with TP53, we discovered new insights into the compelling biology behind transformative genomic events in TNBC initiation and progression, thereby inspiring dynamic basic and translational genomics research to further develop early detection approaches, preventative treatments, and precision medicine therapies in breast cancer

Evolutionary Trajectories of IDH Glioblastomas Reveal a Common Path of Early Tumorigenesis Instigated Years ahead of Initial Diagnosis

We studied how intratumoral genetic heterogeneity shapes tumor growth and therapy response for isocitrate dehydrogenase (IDH)-wild-type glioblastoma, a rapidly regrowing tumor. We inferred the evolutionary trajectories of matched pairs of primary and relapsed tumors based on deep whole-genome-sequencing data. This analysis suggests both a distant origin of de novo glioblastoma, up to 7 years before diagnosis, and a common path of early tumorigenesis, with one or more of chromosome 7 gain, 9p loss, or 10 loss, at tumor initiation. TERT promoter mutations often occurred later as a prerequisite for rapid growth. In contrast to this common early path, relapsed tumors acquired no stereotypical pattern of mutations and typically regrew from oligoclonal origins, suggesting sparse selective pressure by therapeutic measures

Clonality and evolutionary history of rhabdomyosarcoma.

To infer the subclonality of rhabdomyosarcoma (RMS) and predict the temporal order of genetic events for the tumorigenic process, and to identify novel drivers, we applied a systematic method that takes into account germline and somatic alterations in 44 tumor-normal RMS pairs using deep whole-genome sequencing. Intriguingly, we find that loss of heterozygosity of 11p15.5 and mutations in RAS pathway genes occur early in the evolutionary history of the PAX-fusion-negative-RMS (PFN-RMS) subtype. We discover several early mutations in non-RAS mutated samples and predict them to be drivers in PFN-RMS including recurrent mutation of PKN1. In contrast, we find that PAX-fusion-positive (PFP) subtype tumors have undergone whole-genome duplication in the late stage of cancer evolutionary history and have acquired fewer mutations and subclones than PFN-RMS. Moreover we predict that the PAX3-FOXO1 fusion event occurs earlier than the whole genome duplication. Our findings provide information critical to the understanding of tumorigenesis of RMS

Longitudinal Molecular Trajectories of Diffuse Glioma in Adults

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear ¹² . Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of difuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specifc gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at diferent rates across the glioma subtypes, and hypermutation was not associated with diferences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner

Aromatase inhibition remodels the clonal architecture of estrogen-receptor-positive breast cancers

Resistance to oestrogen-deprivation therapy is common in oestrogen-receptor-positive (ER+) breast cancer. To better understand the contributions of tumour heterogeneity and evolution to resistance, here we perform comprehensive genomic characterization of 22 primary tumours sampled before and after 4 months of neoadjuvant aromatase inhibitor (NAI) treatment. Comparing whole-genome sequencing of tumour/normal pairs from the two time points, with coincident tumour RNA sequencing, reveals widespread spatial and temporal heterogeneity, with marked remodelling of the clonal landscape in response to NAI. Two cases have genomic evidence of two independent tumours, most obviously an ER− ‘collision tumour', which was only detected after NAI treatment of baseline ER+ disease. Many mutations are newly detected or enriched post treatment, including two ligand-binding domain mutations in ESR1. The observed clonal complexity of the ER+ breast cancer genome suggests that precision medicine approaches based on genomic analysis of a single specimen are likely insufficient to capture all clinically significant information