Eva Grönroos

Intratumor heterogeneity and branched evolution revealed by multiregion sequencing.

BACKGROUND Intratumor heterogeneity may foster tumor evolution and adaptation and hinder personalized-medicine strategies that depend on results from single tumor-biopsy samples. METHODS To examine intratumor heterogeneity, we performed exome sequencing, chromosome aberration analysis, and ploidy profiling on multiple spatially separated samples obtained from primary renal carcinomas and associated metastatic sites. We characterized the consequences of intratumor heterogeneity using immunohistochemical analysis, mutation functional analysis, and profiling of messenger RNA expression. RESULTS Phylogenetic reconstruction revealed branched evolutionary tumor growth, with 63 to 69% of all somatic mutations not detectable across every tumor region. Intratumor heterogeneity was observed for a mutation within an autoinhibitory domain of the mammalian target of rapamycin (mTOR) kinase, correlating with S6 and 4EBP phosphorylation in vivo and constitutive activation of mTOR kinase activity in vitro. Mutational intratumor heterogeneity was seen for multiple tumor-suppressor genes converging on loss of function; SETD2, PTEN, and KDM5C underwent multiple distinct and spatially separated inactivating mutations within a single tumor, suggesting convergent phenotypic evolution. Gene-expression signatures of good and poor prognosis were detected in different regions of the same tumor. Allelic composition and ploidy profiling analysis revealed extensive intratumor heterogeneity, with 26 of 30 tumor samples from four tumors harboring divergent allelic-imbalance profiles and with ploidy heterogeneity in two of four tumors. CONCLUSIONS Intratumor heterogeneity can lead to underestimation of the tumor genomics landscape portrayed from single tumor-biopsy samples and may present major challenges to personalized-medicine and biomarker development. Intratumor heterogeneity, associated with heterogeneous protein function, may foster tumor adaptation and therapeutic failure through Darwinian selection. (Funded by the Medical Research Council and others.).

Spatial and temporal diversity in genomic instability processes defines lung cancer evolution

Spatial and temporal dissection of the genomic changes occurring during the evolution of human non–small cell lung cancer (NSCLC) may help elucidate the basis for its dismal prognosis. We sequenced 25 spatially distinct regions from seven operable NSCLCs and found evidence of branched evolution, with driver mutations arising before and after subclonal diversification. There was pronounced intratumor heterogeneity in copy number alterations, translocations, and mutations associated with APOBEC cytidine deaminase activity. Despite maintained carcinogen exposure, tumors from smokers showed a relative decrease in smoking-related mutations over time, accompanied by an increase in APOBEC-associated mutations. In tumors from former smokers, genome-doubling occurred within a smoking-signature context before subclonal diversification, which suggested that a long period of tumor latency had preceded clinical detection. The regionally separated driver mutations, coupled with the relentless and heterogeneous nature of the genome instability processes, are likely to confound treatment success in NSCLC. Different regions of a human lung tumor harbor different mutations, possibly explaining why the disease is so tough to treat. [Also see Perspective by Govindan] Space, time, and the lung cancer genome Lung cancer poses a formidable challenge to clinical oncologists. It is often detected at a late stage, and most therapies work for only a short time before the tumors resume their relentless growth. Two independent analyses of the human lung cancer genome may help explain why this disease is so resilient (see the Perspective by Govindan). Rather than take a single “snapshot” of the cancer genome, de Bruin et al. and Zhang et al. identified genomic alterations in spatially distinct regions of single lung tumors and used this information to infer the tumors evolutionary history. Each tumor showed tremendous spatial and temporal diversity in its mutational profiles. Thus, the efficacy of drugs may be short-lived because they destroy only a portion of the tumor. Science, this issue p. 251, p. 256; see also p. 169

Replication stress links structural and numerical cancer chromosomal instability.

Cancer chromosomal instability (CIN) results in an increased rate of change of chromosome number and structure and generates intratumour heterogeneity. CIN is observed in most solid tumours and is associated with both poor prognosis and drug resistance. Understanding a mechanistic basis for CIN is therefore paramount. Here we find evidence for impaired replication fork progression and increased DNA replication stress in CIN+ colorectal cancer (CRC) cells relative to CIN− CRC cells, with structural chromosome abnormalities precipitating chromosome missegregation in mitosis. We identify three new CIN-suppressor genes (PIGN (also known as MCD4), MEX3C (RKHD2) and ZNF516 (KIAA0222)) encoded on chromosome 18q that are subject to frequent copy number loss in CIN+ CRC. Chromosome 18q loss was temporally associated with aneuploidy onset at the adenoma–carcinoma transition. CIN-suppressor gene silencing leads to DNA replication stress, structural chromosome abnormalities and chromosome missegregation. Supplementing cells with nucleosides, to alleviate replication-associated damage, reduces the frequency of chromosome segregation errors after CIN-suppressor gene silencing, and attenuates segregation errors and DNA damage in CIN+ cells. These data implicate a central role for replication stress in the generation of structural and numerical CIN, which may inform new therapeutic approaches to limit intratumour heterogeneity.

Tolerance of Whole-Genome Doubling Propagates Chromosomal Instability and Accelerates Cancer Genome Evolution

UNLABELLED The contribution of whole-genome doubling to chromosomal instability (CIN) and tumor evolution is unclear. We use long-term culture of isogenic tetraploid cells from a stable diploid colon cancer progenitor to investigate how a genome-doubling event affects genome stability over time. Rare cells that survive genome doubling demonstrate increased tolerance to chromosome aberrations. Tetraploid cells do not exhibit increased frequencies of structural or numerical CIN per chromosome. However, the tolerant phenotype in tetraploid cells, coupled with a doubling of chromosome aberrations per cell, allows chromosome abnormalities to evolve specifically in tetraploids, recapitulating chromosomal changes in genomically complex colorectal tumors. Finally, a genome-doubling event is independently predictive of poor relapse-free survival in early-stage disease in two independent cohorts in multivariate analyses [discovery data: hazard ratio (HR), 4.70, 95% confidence interval (CI), 1.04-21.37; validation data: HR, 1.59, 95% CI, 1.05-2.42]. These data highlight an important role for the tolerance of genome doubling in driving cancer genome evolution. SIGNIFICANCE Our work sheds light on the importance of whole-genome–doubling events in colorectal cancer evolution. We show that tetraploid cells undergo rapid genomic changes and recapitulate the genetic alterations seen in chromosomally unstable tumors. Furthermore, we demonstrate that a genome-doubling event is prognostic of poor relapse-free survival in this disease type.

Requirement for Interaction of PI3-Kinase p110α with RAS in Lung Tumor Maintenance

Summary RAS proteins directly activate PI3-kinases. Mice bearing a germline mutation in the RAS binding domain of the p110α subunit of PI3-kinse are resistant to the development of RAS-driven tumors. However, it is unknown whether interaction of RAS with PI3-kinase is required in established tumors. The need for RAS interaction with p110α in the maintenance of mutant Kras-driven lung tumors was explored using an inducible mouse model. In established tumors, removal of the ability of p110α to interact with RAS causes long-term tumor stasis and partial regression. This is a tumor cell-autonomous effect, which is improved significantly by combination with MEK inhibition. Total removal of p110α expression or activity has comparable effects, albeit with greater toxicities.

SETD2 loss-of-function promotes renal cancer branched evolution through replication stress and impaired DNA repair.

Defining mechanisms that generate intratumour heterogeneity and branched evolution may inspire novel therapeutic approaches to limit tumour diversity and adaptation. SETD2 (Su(var), Enhancer of zeste, Trithorax-domain containing 2) trimethylates histone-3 lysine-36 (H3K36me3) at sites of active transcription and is mutated in diverse tumour types, including clear cell renal carcinomas (ccRCCs). Distinct SETD2 mutations have been identified in spatially separated regions in ccRCC, indicative of intratumour heterogeneity. In this study, we have addressed the consequences of SETD2 loss-of-function through an integrated bioinformatics and functional genomics approach. We find that bi-allelic SETD2 aberrations are not associated with microsatellite instability in ccRCC. SETD2 depletion in ccRCC cells revealed aberrant and reduced nucleosome compaction and chromatin association of the key replication proteins minichromosome maintenance complex component (MCM7) and DNA polymerase δ hindering replication fork progression, and failure to load lens epithelium-derived growth factor and the Rad51 homologous recombination repair factor at DNA breaks. Consistent with these data, we observe chromosomal breakpoint locations are biased away from H3K36me3 sites in SETD2 wild-type ccRCCs relative to tumours with bi-allelic SETD2 aberrations and that H3K36me3-negative ccRCCs display elevated DNA damage in vivo. These data suggest a role for SETD2 in maintaining genome integrity through nucleosome stabilization, suppression of replication stress and the coordination of DNA repair.

Transforming Growth Factor β Inhibits Bone Morphogenetic Protein-Induced Transcription through Novel Phosphorylated Smad1/5-Smad3 Complexes

ABSTRACT In vivo cells receive simultaneous signals from multiple extracellular ligands and must integrate and interpret them to respond appropriately. Here we investigate the interplay between pathways downstream of two transforming growth factor β (TGF-β) superfamily members, bone morphogenetic protein (BMP) and TGF-β. We show that in multiple cell lines, TGF-β potently inhibits BMP-induced transcription at the level of both BMP-responsive reporter genes and endogenous BMP target genes. This inhibitory effect requires the TGF-β type I receptor ALK5 and is independent of new protein synthesis. Strikingly, we show that Smad3 is required for TGF-βs inhibitory effects, whereas Smad2 is not. We go on to demonstrate that TGF-β induces the formation of complexes comprising phosphorylated Smad1/5 and Smad3, which bind to BMP-responsive elements in vitro and in vivo and mediate TGF-β-induced transcriptional repression. Furthermore, loss of Smad3 confers on TGF-β the ability to induce transcription via BMP-responsive elements. Our results therefore suggest that not only is Smad3 important for mediating TGF-βs inhibitory effects on BMP signaling but it also plays a critical role in restricting the transcriptional output in response to TGF-β.

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

BackgroundGenomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution.ResultsWe 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.ConclusionsIn 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.

BCL9L Dysfunction Impairs Caspase-2 Expression Permitting Aneuploidy Tolerance in Colorectal Cancer

Summary Chromosomal instability (CIN) contributes to cancer evolution, intratumor heterogeneity, and drug resistance. CIN is driven by chromosome segregation errors and a tolerance phenotype that permits the propagation of aneuploid genomes. Through genomic analysis of colorectal cancers and cell lines, we find frequent loss of heterozygosity and mutations in BCL9L in aneuploid tumors. BCL9L deficiency promoted tolerance of chromosome missegregation events, propagation of aneuploidy, and genetic heterogeneity in xenograft models likely through modulation of Wnt signaling. We find that BCL9L dysfunction contributes to aneuploidy tolerance in both TP53-WT and mutant cells by reducing basal caspase-2 levels and preventing cleavage of MDM2 and BID. Efforts to exploit aneuploidy tolerance mechanisms and the BCL9L/caspase-2/BID axis may limit cancer diversity and evolution.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal

Summary The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.