Stuart Horswell

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.).

Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing.

Clear cell renal carcinomas (ccRCCs) can display intratumor heterogeneity (ITH). We applied multiregion exome sequencing (M-seq) to resolve the genetic architecture and evolutionary histories of ten ccRCCs. Ultra-deep sequencing identified ITH in all cases. We found that 73–75% of identified ccRCC driver aberrations were subclonal, confounding estimates of driver mutation prevalence. ITH increased with the number of biopsies analyzed, without evidence of saturation in most tumors. Chromosome 3p loss and VHL aberrations were the only ubiquitous events. The proportion of C>T transitions at CpG sites increased during tumor progression. M-seq permits the temporal resolution of ccRCC evolution and refines mutational signatures occurring during tumor development.

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

Tracking the Evolution of Non–Small-Cell Lung Cancer

BACKGROUND Among patients with non‐small‐cell lung cancer (NSCLC), data on intratumor heterogeneity and cancer genome evolution have been limited to small retrospective cohorts. We wanted to prospectively investigate intratumor heterogeneity in relation to clinical outcome and to determine the clonal nature of driver events and evolutionary processes in early‐stage NSCLC. METHODS In this prospective cohort study, we performed multiregion whole‐exome sequencing on 100 early‐stage NSCLC tumors that had been resected before systemic therapy. We sequenced and analyzed 327 tumor regions to define evolutionary histories, obtain a census of clonal and subclonal events, and assess the relationship between intratumor heterogeneity and recurrence‐free survival. RESULTS We observed widespread intratumor heterogeneity for both somatic copy‐number alterations and mutations. Driver mutations in EGFR, MET, BRAF, and TP53 were almost always clonal. However, heterogeneous driver alterations that occurred later in evolution were found in more than 75% of the tumors and were common in PIK3CA and NF1 and in genes that are involved in chromatin modification and DNA damage response and repair. Genome doubling and ongoing dynamic chromosomal instability were associated with intratumor heterogeneity and resulted in parallel evolution of driver somatic copy‐number alterations, including amplifications in CDK4, FOXA1, and BCL11A. Elevated copy‐number heterogeneity was associated with an increased risk of recurrence or death (hazard ratio, 4.9; P=4.4×10‐4), which remained significant in multivariate analysis. CONCLUSIONS Intratumor heterogeneity mediated through chromosome instability was associated with an increased risk of recurrence or death, a finding that supports the potential value of chromosome instability as a prognostic predictor. (Funded by Cancer Research UK and others; TRACERx ClinicalTrials.gov number, NCT01888601.)

Rho-actin signaling to the MRTF coactivators dominates the immediate transcriptional response to serum in fibroblasts

The transcription factor SRF (serum response factor) recruits two families of coactivators, the MRTFs (myocardin-related transcription factors) and the TCFs (ternary complex factors), to couple gene transcription to growth factor signaling. Here we investigated the role of the SRF network in the immediate transcriptional response of fibroblasts to serum stimulation. SRF recruited its cofactors in a gene-specific manner, and virtually all MRTF binding was directed by SRF. Much of SRF DNA binding was serum-inducible, reflecting a requirement for MRTF-SRF complex formation in nucleosome displacement. We identified 960 serum-responsive SRF target genes, which were mostly MRTF-controlled, as assessed by MRTF chromatin immunoprecipitation (ChIP) combined with deep sequencing (ChIP-seq) and/or sensitivity to MRTF-linked signals. MRTF activation facilitates RNA polymerase II (Pol II) recruitment or promoter escape according to gene context. MRTF targets encode regulators of the cytoskeleton, transcription, and cell growth, underpinning the role of SRF in cytoskeletal dynamics and mechanosensing. Finally, we show that specific activation of either MRTFs or TCFs can reset the circadian clock.

Tracking genomic cancer evolution for precision medicine: the lung TRACERx study.

TRACERx, a prospective study of patients with primary non-small cell lung cancer, aims to map the genomic landscape of lung cancer by tracking clonal heterogeneity and tumour evolution from diagnosis to relapse.

Tracking the genomic evolution of esophageal adenocarcinoma through neoadjuvant chemotherapy

UNLABELLED Esophageal adenocarcinomas are associated with a dismal prognosis. Deciphering the evolutionary history of this disease may shed light on therapeutically tractable targets and reveal dynamic mutational processes during the disease course and following neoadjuvant chemotherapy (NAC). We exome sequenced 40 tumor regions from 8 patients with operable esophageal adenocarcinomas, before and after platinum-containing NAC. This revealed the evolutionary genomic landscape of esophageal adenocarcinomas with the presence of heterogeneous driver mutations, parallel evolution, early genome-doubling events, and an association between high intratumor heterogeneity and poor response to NAC. Multiregion sequencing demonstrated a significant reduction in thymine to guanine mutations within a CpTpT context when comparing early and late mutational processes and the presence of a platinum signature with enrichment of cytosine to adenine mutations within a CpC context following NAC. Esophageal adenocarcinomas are characterized by early chromosomal instability leading to amplifications containing targetable oncogenes persisting through chemotherapy, providing a rationale for future therapeutic approaches. SIGNIFICANCE This work illustrates dynamic mutational processes occurring during esophageal adenocarcinoma evolution and following selective pressures of platinum exposure, emphasizing the iatrogenic impact of therapy on cancer evolution. Identification of amplifications encoding targetable oncogenes maintained through NAC suggests the presence of stable vulnerabilities, unimpeded by cytotoxics, suitable for therapeutic intervention.

Recurrent chromosomal gains and heterogeneous driver mutations characterise papillary renal cancer evolution

Papillary renal cell carcinoma (pRCC) is an important subtype of kidney cancer with a problematic pathological classification and highly variable clinical behaviour. Here we sequence the genomes or exomes of 31 pRCCs, and in four tumours, multi-region sequencing is undertaken. We identify BAP1, SETD2, ARID2 and Nrf2 pathway genes (KEAP1, NHE2L2 and CUL3) as probable drivers, together with at least eight other possible drivers. However, only ~10% of tumours harbour detectable pathogenic changes in any one driver gene, and where present, the mutations are often predicted to be present within cancer sub-clones. We specifically detect parallel evolution of multiple SETD2 mutations within different sub-regions of the same tumour. By contrast, large copy number gains of chromosomes 7, 12, 16 and 17 are usually early, monoclonal changes in pRCC evolution. The predominance of large copy number variants as the major drivers for pRCC highlights an unusual mode of tumorigenesis that may challenge precision medicine approaches.

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.

Mutation of cancer driver MLL2 results in transcription stress and genome instability

Genome instability is a recurring feature of tumorigenesis. Mutation in MLL2, encoding a histone methyltransferase, is a driver in numerous different cancer types, but the mechanism is unclear. Here, we present evidence that MLL2 mutation results in genome instability. Mouse cells in which MLL2 gene deletion can be induced display elevated levels of sister chromatid exchange, gross chromosomal aberrations, 53BP1 foci, and micronuclei. Human MLL2 knockout cells are characterized by genome instability as well. Interestingly, MLL2 interacts with RNA polymerase II (RNAPII) and RECQL5, and, although MLL2 mutated cells have normal overall H3K4me levels in genes, nucleosomes in the immediate vicinity of RNAPII are hypomethylated. Importantly, MLL2 mutated cells display signs of substantial transcription stress, and the most affected genes overlap with early replicating fragile sites, show elevated levels of γH2AX, and suffer frequent mutation. The requirement for MLL2 in the maintenance of genome stability in genes helps explain its widespread role in cancer and points to transcription stress as a strong driver in tumorigenesis.