A general framework for analyzing tumor subclonality using SNP array and DNA sequencing data

Intra-tumor heterogeneity reflects cancer genome evolution and provides key information for diagnosis and treatment. When bulk tumor tissues are profiled for somatic copy number alterations (sCNA) and point mutations, it may be difficult to estimate their cellular fractions when a mutation falls within a sCNA. We present the Clonal Heterogeneity Analysis Tool, which estimates cellular fractions for both sCNAs and mutations, and uses their distributions to inform macroscopic clonal architecture. In a set of approximately 700 breast tumors, more than half appear to contain multiple recognizable aneuploid tumor clones, and many show subtype-specific differences in clonality for known cancer genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0473-4) contains supplementary material, which is available to authorized users.     

肿瘤遗传异质性

尽管大多数自发性肿瘤源于单细胞, 但人们逐渐认识到个体肿瘤的内部存在异质性, 这种差异涉及分化程度、细胞增殖率、侵袭和转移能力及治疗反应等众多方面。分子生物学研究也证实肿瘤演进过程中不断产生新的突变, 为肿瘤内部存在异质性增加了更有力的佐证。文章综述了关于肿瘤遗传异质性研究的主要进展。鉴于遗传异质性分析可为揭示肿瘤细胞的产生、扩散、转移的时间段提供重要信息, 文章以瘤内异质性为主线, 列举了遗传异质性存在的实验依据; 阐述了遗传多样性在人类个体肿瘤发展进化史中的研究价值; 介绍了遗传异质性产生的两种模式, 包括肿瘤干细胞模型和克隆进化模型; 总结了遗传异质性在肿瘤转移、治疗中的重要临床意义。文章最后展示了常用的遗传异质性的研究方法, 包括特定基因的分析方法和基于基因组水平的方法, 并对各自的优缺点进行了分析。     

Tumor heterogeneity: Causes and consequences

With rare exceptions, spontaneous tumors originate from a single cell. Yet, at the time of clinical diagnosis, the majority of human tumors display startling heterogeneity in many morphological and physiological features, such as expression of cell surface receptors, proliferative and angiogenic potential. To a substantial extent, this heterogeneity might be attributed to morphological and epigenetic plasticity, but there is also strong evidence for the co-existence of genetically divergent tumor cell clones within tumors. In this perspective, we summarize the sources of intra-tumor phenotypic heterogeneity with emphasis on genetic heterogeneity. We review experimental evidence for the existence of both intra-tumor clonal heterogeneity as well as frequent evolutionary divergence between primary tumors and metastatic outgrowths. Furthermore, we discuss potential biological and clinical implications of intra-tumor clonal heterogeneity.     

Intra-tumor heterogeneity of MLH1 promoter methylation revealed by deep single molecule bisulfite sequencing

A single tumor may contain cells with different somatic mutations. By characterizing this genetic heterogeneity within tumors, advances have been made in the prognosis, treatment and understanding of tumorigenesis. In contrast, the extent of epigenetic intra-tumor heterogeneity and how it influences tumor biology is under-explored. We have characterized epigenetic heterogeneity within individual tumors using next-generation sequencing. We used deep single molecule bisulfite sequencing and sample-specific DNA barcodes to determine the spectrum of MLH1 promoter methylation across an average of 1000 molecules in each of 33 individual samples in parallel, including endometrial cancer, matched blood and normal endometrium. This first glimpse, deep into each tumor, revealed unexpectedly heterogeneous patterns of methylation at the MLH1 promoter within a subset of endometrial tumors. This high-resolution analysis allowed us to measure the clonality of methylation in individual tumors and gain insight into the accumulation of aberrant promoter methylation on both alleles during tumorigenesis.     

Development of synchronous VHL syndrome tumors reveals contingencies and constraints to tumor evolution

Background

Genomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution.

Results

We perform whole exome sequencing on four clear cell renal cell carcinomas removed from both kidneys of a patient with a germline VHL mutation. We report that tumors arising in this context are clonally independent and harbour distinct secondary events exemplified by loss of chromosome 3p, despite an identical genetic background and tissue microenvironment. We propose that divergent mutational and copy number anomalies are contingent upon the nature of 3p loss of heterozygosity occurring early in tumorigenesis. However, despite distinct 3p events, genomic, proteomic and immunohistochemical analyses reveal evidence for convergence upon the PI3K-AKT-mTOR signaling pathway. Four germline tumors in this young patient, and in a second, older patient with VHL syndrome demonstrate minimal intra-tumor heterogeneity and mutational burden, and evaluable tumors appear to follow a linear evolutionary route, compared to tumors from patients with sporadic clear cell renal cell carcinoma.

Conclusions

In tumors developing from a germline VHL mutation, the evolutionary principles of contingency and convergence in tumor development are complementary. In this small set of patients with early stage VHL-associated tumors, there is reduced mutation burden and limited evidence of intra-tumor heterogeneity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0433-z) contains supplementary material, which is available to authorized users.     

Intra-tumor Genetic Heterogeneity and Mortality in Head and Neck Cancer: Analysis of Data from The Cancer Genome Atlas

Background

Although the involvement of intra-tumor genetic heterogeneity in tumor progression, treatment resistance, and metastasis is established, genetic heterogeneity is seldom examined in clinical trials or practice. Many studies of heterogeneity have had prespecified markers for tumor subpopulations, limiting their generalizability, or have involved massive efforts such as separate analysis of hundreds of individual cells, limiting their clinical use. We recently developed a general measure of intra-tumor genetic heterogeneity based on whole-exome sequencing (WES) of bulk tumor DNA, called mutant-allele tumor heterogeneity (MATH). Here, we examine data collected as part of a large, multi-institutional study to validate this measure and determine whether intra-tumor heterogeneity is itself related to mortality.

Methods and Findings

Clinical and WES data were obtained from The Cancer Genome Atlas in October 2013 for 305 patients with head and neck squamous cell carcinoma (HNSCC), from 14 institutions. Initial pathologic diagnoses were between 1992 and 2011 (median, 2008). Median time to death for 131 deceased patients was 14 mo; median follow-up of living patients was 22 mo. Tumor MATH values were calculated from WES results. Despite the multiple head and neck tumor subsites and the variety of treatments, we found in this retrospective analysis a substantial relation of high MATH values to decreased overall survival (Cox proportional hazards analysis: hazard ratio for high/low heterogeneity, 2.2; 95% CI 1.4 to 3.3). This relation of intra-tumor heterogeneity to survival was not due to intra-tumor heterogeneity’s associations with other clinical or molecular characteristics, including age, human papillomavirus status, tumor grade and TP53 mutation, and N classification. MATH improved prognostication over that provided by traditional clinical and molecular characteristics, maintained a significant relation to survival in multivariate analyses, and distinguished outcomes among patients having oral-cavity or laryngeal cancers even when standard disease staging was taken into account. Prospective studies, however, will be required before MATH can be used prognostically in clinical trials or practice. Such studies will need to examine homogeneously treated HNSCC at specific head and neck subsites, and determine the influence of cancer therapy on MATH values. Analysis of MATH and outcome in human-papillomavirus-positive oropharyngeal squamous cell carcinoma is particularly needed.

Conclusions

To our knowledge this study is the first to combine data from hundreds of patients, treated at multiple institutions, to document a relation between intra-tumor heterogeneity and overall survival in any type of cancer. We suggest applying the simply calculated MATH metric of heterogeneity to prospective studies of HNSCC and other tumor types.     

Deep sequencing of multiple regions of glial tumors reveals spatial heterogeneity for mutations in clinically relevant genes

Background

The extent of intratumoral mutational heterogeneity remains unclear in gliomas, the most common primary brain tumors, especially with respect to point mutation. To address this, we applied single molecule molecular inversion probes targeting 33 cancer genes to assay both point mutations and gene amplifications within spatially distinct regions of 14 glial tumors.

Results

We find evidence of regional mutational heterogeneity in multiple tumors, including mutations in TP53 and RB1 in an anaplastic oligodendroglioma and amplifications in PDGFRA and KIT in two glioblastomas (GBMs). Immunohistochemistry confirms heterogeneity of TP53 mutation and PDGFRA amplification. In all, 3 out of 14 glial tumors surveyed have evidence for heterogeneity for clinically relevant mutations.

Conclusions

Our results underscore the need to sample multiple regions in GBM and other glial tumors when devising personalized treatments based on genomic information, and furthermore demonstrate the importance of measuring both point mutation and copy number alteration while investigating genetic heterogeneity within cancer samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0530-z) contains supplementary material, which is available to authorized users.     

Microenvironmental Heterogeneity Parallels Breast Cancer Progression: A Histology–Genomic Integration Analysis

Background

The intra-tumor diversity of cancer cells is under intense investigation; however, little is known about the heterogeneity of the tumor microenvironment that is key to cancer progression and evolution. We aimed to assess the degree of microenvironmental heterogeneity in breast cancer and correlate this with genomic and clinical parameters.

Methods and Findings

We developed a quantitative measure of microenvironmental heterogeneity along three spatial dimensions (3-D) in solid tumors, termed the tumor ecosystem diversity index (EDI), using fully automated histology image analysis coupled with statistical measures commonly used in ecology. This measure was compared with disease-specific survival, key mutations, genome-wide copy number, and expression profiling data in a retrospective study of 510 breast cancer patients as a test set and 516 breast cancer patients as an independent validation set. In high-grade (grade 3) breast cancers, we uncovered a striking link between high microenvironmental heterogeneity measured by EDI and a poor prognosis that cannot be explained by tumor size, genomics, or any other data types. However, this association was not observed in low-grade (grade 1 and 2) breast cancers. The prognostic value of EDI was superior to known prognostic factors and was enhanced with the addition of TP53 mutation status (multivariate analysis test set, p = 9 × 10⁻⁴, hazard ratio = 1.47, 95% CI 1.17–1.84; validation set, p = 0.0011, hazard ratio = 1.78, 95% CI 1.26–2.52). Integration with genome-wide profiling data identified losses of specific genes on 4p14 and 5q13 that were enriched in grade 3 tumors with high microenvironmental diversity that also substratified patients into poor prognostic groups. Limitations of this study include the number of cell types included in the model, that EDI has prognostic value only in grade 3 tumors, and that our spatial heterogeneity measure was dependent on spatial scale and tumor size.

Conclusions

To our knowledge, this is the first study to couple unbiased measures of microenvironmental heterogeneity with genomic alterations to predict breast cancer clinical outcome. We propose a clinically relevant role of microenvironmental heterogeneity for advanced breast tumors, and highlight that ecological statistics can be translated into medical advances for identifying a new type of biomarker and, furthermore, for understanding the synergistic interplay of microenvironmental heterogeneity with genomic alterations in cancer cells.     

Expression profiling of human renal carcinomas with functional taxonomic analysis

Background  

Molecular characterization has contributed to the understanding of the inception, progression, treatment and prognosis of cancer. Nucleic acid array-based technologies extend molecular characterization of tumors to thousands of gene products. To effectively discriminate between tumor sub-types, reliable laboratory techniques and analytic methods are required.     

A minimally invasive multiple marker approach allows highly efficient detection of meningioma tumors

Background  

The development of effective frameworks that permit an accurate diagnosis of tumors, especially in their early stages, remains a grand challenge in the field of bioinformatics. Our approach uses statistical learning techniques applied to multiple antigen tumor antigen markers utilizing the immune system as a very sensitive marker of molecular pathological processes. For validation purposes we choose the intracranial meningioma tumors as model system since they occur very frequently, are mostly benign, and are genetically stable.     

Deficiency of Pten accelerates mammary oncogenesis in MMTV-Wnt-1 transgenic mice

Background  

Germline mutations in the tumor suppressor PTEN predispose human beings to breast cancer, and genetic and epigenetic alterations of PTEN are also detected in sporadic human breast cancer. Germline Pten mutations in mice lead to the development of a variety of tumors, but mammary carcinomas are infrequently found, especially in mice under the age of six months.     

Clonal Analysis in Recurrent Astrocytic, Oligoastrocytic and Oligodendroglial Tumors Implicates IDH1- Mutation as Common Tumor Initiating Event

Background

To investigate the dynamics of inter- and intratumoral molecular alterations during tumor progression in recurrent gliomas.

Methodology/Principal Findings

To address intertumoral heterogeneity we investigated non- microdissected tumor tissue of 106 gliomas representing 51 recurrent tumors. To address intratumoral heterogeneity a set of 16 gliomas representing 7 tumor pairs with at least one recurrence, and 4 single mixed gliomas were investigated by microdissection of distinct oligodendroglial and astrocytic tumor components. All tumors and tumor components were analyzed for allelic loss of 1p/19q (LOH 1p/19q), for TP53- mutations and for R132 mutations in the IDH1 gene. The investigation of non- microdissected tumor tissue revealed clonality in 75% (38/51). Aberrant molecular alterations upon recurrence were detected in 25% (13/51). 64% (9/14) of these were novel and associated with tumor progression. Loss of previously detected alterations was observed in 36% (5/14). One tumor pair (1/14; 7%) was significant for both. Intratumoral clonality was detected in 57% (4/7) of the microdissected tumor pairs and in 75% (3/4) of single microdissected tumors. 43% (3/7) of tumor pairs and one single tumor (25%) revealed intratumoral heterogeneity. While intratumoral heterogeneity affected both the TP53- mutational status and the LOH1p/19q status, all tumors with intratumoral heterogeneity shared the R132 IDH1- mutation as a common feature in both their microdissected components.

Conclusions/Significance

The majority of recurrent gliomas are of monoclonal origin. However, the detection of divertive tumor cell clones in morphological distinct tumor components sharing IDH1- mutations as early event may provide insight into the tumorigenesis of true mixed gliomas.     

Quantifying Metabolic Heterogeneity in Head and Neck Tumors in Real Time: 2-DG Uptake Is Highest in Hypoxic Tumor Regions

Purpose

Intratumoral metabolic heterogeneity may increase the likelihood of treatment failure due to the presence of a subset of resistant tumor cells. Using a head and neck squamous cell carcinoma (HNSCC) xenograft model and a real-time fluorescence imaging approach, we tested the hypothesis that tumors are metabolically heterogeneous, and that tumor hypoxia alters patterns of glucose uptake within the tumor.

Experimental Design

Cal33 cells were grown as xenograft tumors (n = 16) in nude mice after identification of this cell line''s metabolic response to hypoxia. Tumor uptake of fluorescent markers identifying hypoxia, glucose import, or vascularity was imaged simultaneously using fluorescent molecular tomography. The variability of intratumoral 2-deoxyglucose (IR800-2-DG) concentration was used to assess tumor metabolic heterogeneity, which was further investigated using immunohistochemistry for expression of key metabolic enzymes. HNSCC tumors in patients were assessed for intratumoral variability of ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) uptake in clinical PET scans.

Results

IR800-2-DG uptake in hypoxic regions of Cal33 tumors was 2.04 times higher compared to the whole tumor (p = 0.0001). IR800-2-DG uptake in tumors containing hypoxic regions was more heterogeneous as compared to tumors lacking a hypoxic signal. Immunohistochemistry staining for HIF-1α, carbonic anhydrase 9, and ATP synthase subunit 5β confirmed xenograft metabolic heterogeneity. We detected heterogeneous ¹⁸F-FDG uptake within patient HNSCC tumors, and the degree of heterogeneity varied amongst tumors.

Conclusion

Hypoxia is associated with increased intratumoral metabolic heterogeneity. ¹⁸F-FDG PET scans may be used to stratify patients according to the metabolic heterogeneity within their tumors, which could be an indicator of prognosis.     

Tumor Phosphatidylinositol-3-Kinase Signaling and Development of Metastatic Disease in Locally Advanced Rectal Cancer

Background

Recognizing EGFR as key orchestrator of the metastatic process in colorectal cancer, but also the substantial heterogeneity of responses to anti-EGFR therapy, we examined the pattern of composite tumor kinase activities governed by EGFR-mediated signaling that might be implicated in development of metastatic disease.

Patients and Methods

Point mutations in KRAS, BRAF, and PIK3CA and ERBB2 amplification were determined in primary tumors from 63 patients with locally advanced rectal cancer scheduled for radical treatment. Using peptide arrays with tyrosine kinase substrates, ex vivo phosphopeptide profiles were generated from the same baseline tumor samples and correlated to metastasis-free survival.

Results

Unsupervised clustering analysis of the resulting phosphorylation of 102 array substrates defined two tumor classes, both consisting of cases with and without KRAS/BRAF mutations. The smaller cluster group of patients, with tumors generating high ex vivo phosphorylation of phosphatidylinositol-3-kinase-related substrates, had a particularly aggressive disease course, with almost a half of patients developing metastatic disease within one year of follow-up.

Conclusion

High phosphatidylinositol-3-kinase-mediated signaling activity of the primary tumor, rather than KRAS/BRAF mutation status, was identified as a hallmark of poor metastasis-free survival in patients with locally advanced rectal cancer undergoing radical treatment of the pelvic cavity.     

Reporting Tumor Molecular Heterogeneity in Histopathological Diagnosis

Background

Detection of molecular tumor heterogeneity has become of paramount importance with the advent of targeted therapies. Analysis for detection should be comprehensive, timely and based on routinely available tumor samples.

Aim

To evaluate the diagnostic potential of targeted multigene next-generation sequencing (TM-NGS) in characterizing gastrointestinal cancer molecular heterogeneity.

Methods

35 gastrointestinal tract tumors, five of each intestinal type gastric carcinomas, pancreatic ductal adenocarcinomas, pancreatic intraductal papillary mucinous neoplasms, ampulla of Vater carcinomas, hepatocellular carcinomas, cholangiocarcinomas, pancreatic solid pseudopapillary tumors were assessed for mutations in 46 cancer-associated genes, using Ion Torrent semiconductor-based TM-NGS. One ampulla of Vater carcinoma cell line and one hepatic carcinosarcoma served to assess assay sensitivity. TP53, PIK3CA, KRAS, and BRAF mutations were validated by conventional Sanger sequencing.

Results

TM-NGS yielded overlapping results on matched fresh-frozen and formalin-fixed paraffin-embedded (FFPE) tissues, with a mutation detection limit of 1% for fresh-frozen high molecular weight DNA and 2% for FFPE partially degraded DNA. At least one somatic mutation was observed in all tumors tested; multiple alterations were detected in 20/35 (57%) tumors. Seven cancers displayed significant differences in allelic frequencies for distinct mutations, indicating the presence of intratumor molecular heterogeneity; this was confirmed on selected samples by immunohistochemistry of p53 and Smad4, showing concordance with mutational analysis.

Conclusions

TM-NGS is able to detect and quantitate multiple gene alterations from limited amounts of DNA, moving one step closer to a next-generation histopathologic diagnosis that integrates morphologic, immunophenotypic, and multigene mutational analysis on routinely processed tissues, essential for personalized cancer therapy.     

Deep sequencing of the X chromosome reveals the proliferation history of colorectal adenomas

Background

Mismatch repair deficient colorectal adenomas are composed of transformed cells that descend from a common founder and progressively accumulate genomic alterations. The proliferation history of these tumors is still largely unknown. Here we present a novel approach to rebuild the proliferation trees that recapitulate the history of individual colorectal adenomas by mapping the progressive acquisition of somatic point mutations during tumor growth.

Results

Using our approach, we called high and low frequency mutations acquired in the X chromosome of four mismatch repair deficient colorectal adenomas deriving from male individuals. We clustered these mutations according to their frequencies and rebuilt the proliferation trees directly from the mutation clusters using a recursive algorithm. The trees of all four lesions were formed of a dominant subclone that co-existed with other genetically heterogeneous subpopulations of cells. However, despite this similar hierarchical organization, the growth dynamics varied among and within tumors, likely depending on a combination of tumor-specific genetic and environmental factors.

Conclusions

Our study provides insights into the biological properties of individual mismatch repair deficient colorectal adenomas that may influence their growth and also the response to therapy. Extended to other solid tumors, our novel approach could inform on the mechanisms of cancer progression and on the best treatment choice.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0437-8) contains supplementary material, which is available to authorized users.     

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors

Background

Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations).

Methods

Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas.

Results

We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10–15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair.

Conclusion

We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.