THetA: inferring intra-tumor heterogeneity from high-throughput DNA sequencing data

Tumor samples are typically heterogeneous, containing admixture by normal, non-cancerous cells and one or more subpopulations of cancerous cells. Whole-genome sequencing of a tumor sample yields reads from this mixture, but does not directly reveal the cell of origin for each read. We introduce THetA (Tumor Heterogeneity Analysis), an algorithm that infers the most likely collection of genomes and their proportions in a sample, for the case where copy number aberrations distinguish subpopulations. THetA successfully estimates normal admixture and recovers clonal and subclonal copy number aberrations in real and simulated sequencing data. THetA is available at http://compbio.cs.brown.edu/software/.     

Whole-Genome Sequencing Reveals Diverse Models of Structural Variations in Esophageal Squamous Cell Carcinoma

Comprehensive identification of somatic structural variations (SVs) and understanding their mutational mechanisms in cancer might contribute to understanding biological differences and help to identify new therapeutic targets. Unfortunately, characterization of complex SVs across the whole genome and the mutational mechanisms underlying esophageal squamous cell carcinoma (ESCC) is largely unclear. To define a comprehensive catalog of somatic SVs, affected target genes, and their underlying mechanisms in ESCC, we re-analyzed whole-genome sequencing (WGS) data from 31 ESCCs using Meerkat algorithm to predict somatic SVs and Patchwork to determine copy-number changes. We found deletions and translocations with NHEJ and alt-EJ signature as the dominant SV types, and 16% of deletions were complex deletions. SVs frequently led to disruption of cancer-associated genes (e.g., CDKN2A and NOTCH1) with different mutational mechanisms. Moreover, chromothripsis, kataegis, and breakage-fusion-bridge (BFB) were identified as contributing to locally mis-arranged chromosomes that occurred in 55% of ESCCs. These genomic catastrophes led to amplification of oncogene through chromothripsis-derived double-minute chromosome formation (e.g., FGFR1 and LETM2) or BFB-affected chromosomes (e.g., CCND1, EGFR, ERBB2, MMPs, and MYC), with approximately 30% of ESCCs harboring BFB-derived CCND1 amplification. Furthermore, analyses of copy-number alterations reveal high frequency of whole-genome duplication (WGD) and recurrent focal amplification of CDCA7 that might act as a potential oncogene in ESCC. Our findings reveal molecular defects such as chromothripsis and BFB in malignant transformation of ESCCs and demonstrate diverse models of SVs-derived target genes in ESCCs. These genome-wide SV profiles and their underlying mechanisms provide preventive, diagnostic, and therapeutic implications for ESCCs.     

Recurrent patterns of DNA copy number alterations in tumors reflect metabolic selection pressures

Copy number alteration (CNA) profiling of human tumors has revealed recurrent patterns of DNA amplifications and deletions across diverse cancer types. These patterns are suggestive of conserved selection pressures during tumor evolution but cannot be fully explained by known oncogenes and tumor suppressor genes. Using a pan‐cancer analysis of CNA data from patient tumors and experimental systems, here we show that principal component analysis‐defined CNA signatures are predictive of glycolytic phenotypes, including ¹⁸F‐fluorodeoxy‐glucose (FDG) avidity of patient tumors, and increased proliferation. The primary CNA signature is enriched for p53 mutations and is associated with glycolysis through coordinate amplification of glycolytic genes and other cancer‐linked metabolic enzymes. A pan‐cancer and cross‐species comparison of CNAs highlighted 26 consistently altered DNA regions, containing 11 enzymes in the glycolysis pathway in addition to known cancer‐driving genes. Furthermore, exogenous expression of hexokinase and enolase enzymes in an experimental immortalization system altered the subsequent copy number status of the corresponding endogenous loci, supporting the hypothesis that these metabolic genes act as drivers within the conserved CNA amplification regions. Taken together, these results demonstrate that metabolic stress acts as a selective pressure underlying the recurrent CNAs observed in human tumors, and further cast genomic instability as an enabling event in tumorigenesis and metabolic evolution.     

Discovering tumor suppressor genes through genome-wide copy number analysis

Classical tumor suppressor gene discovery has largely involved linkage analysis and loss-of-heterozygosity (LOH) screens, followed by detailed mapping of relatively large chromosomal regions. Subsequent efforts made use of genome-wide PCR-based methods to detect rare homozygous deletions. More recently, high-resolution genomic arrays have been applied to cancer gene discovery. However, accurate characterization of regions of genomic loss is particularly challenging due to sample heterogeneity, the small size of deleted regions and the high frequency of germline copy number polymorphisms. Here, we review the application of genome-wide copy number analysis to the specific problem of identifying tumor suppressor genes.     

Interrogating cell division errors using random and chromosome-specific missegregation approaches

Accurate segregation of the duplicated genome in mitosis is essential for maintaining genetic stability. Errors in this process can cause numerical and/or structural chromosome abnormalities – hallmark genomic features commonly associated with both tumorigenesis and developmental disorders. A cell-based approach was recently developed permitting inducible missegregation of the human Y chromosome by selectively disrupting kinetochore assembly onto the Y centromere. Although this strategy initially requires several steps of genetic manipulation, it is easy to use, highly efficient and specific for the Y without affecting the autosomes or the X, and does not require cell cycle synchronization or mitotic perturbation. Here we describe currently available tools for studying chromosome segregation errors, aneuploidy, and micronuclei, as well as discuss how the Y-specific missegregation system has been used to elucidate how chromosomal micronucleation can trigger a class of extensive rearrangements termed chromothripsis. The combinatorial use of these different tools will allow unresolved aspects of cell division defects and chromosomal instability to be experimentally explored.     

Genome instability: Does genetic diversity amplification drive tumorigenesis?

Recent data show that catastrophic events during one cell cycle can cause massive genome damage producing viable clones with unstable genomes. This is in contrast with the traditional view that tumorigenesis requires a long‐term process in which mutations gradually accumulate over decades. These sudden events are likely to result in a large increase in genomic diversity within a relatively short time, providing the opportunity for selective advantages to be gained by a subset of cells within a population. This genetic diversity amplification, arising from a single aberrant cell cycle, may drive a population conversion from benign to malignant. However, there is likely a period of relative genome stability during the clonal expansion of tumors – this may provide an opportunity for therapeutic intervention, especially if mechanisms that limit tolerance of aneuploidy are exploited.     

Polyploid mitosis and depolyploidization promote chromosomal instability and tumor progression in a Notch-induced tumor model

1. Download : Download high-res image (210KB)
2. Download : Download full-size image

     

Integrated Analysis of Tumor Samples Sheds Light on Tumor Heterogeneity

The heterogeneity of tumor samples is a major challenge in the analysis of high-throughput profiling of tumor biopsies and cell lines. The measured aggregate signals of multigenerational progenies often represent an average of several tumor subclones with varying genomic aberrations and different gene expression levels. The goal of the present study was to integrate copy number analyses from SNP-arrays and karyotyping, gene expression profiling, and pathway analyses to detect heterogeneity, identify driver mutations, and explore possible mechanisms of tumor evolution. We showed the heterogeneity of the studied samples, characterized the global copy number alteration profiles, and identified genes whose copy number status and expression levels were aberrant. In particular, we identified a recurrent association between two BRAF^V600E and BRAF^V600K mutations and changes in DKK1 gene expression levels, which might indicate an association between the BRAF and WNT pathways. These findings show that the integrated approaches used in the present study can robustly address the challenging issue of tumor heterogeneity in high-throughput profiling.