Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients.

Metastatic progression is the main cause of death in cancer patients, whereas the underlying genomic mechanisms driving metastasis remain largely unknown. Here, we assembled MSK-MET, a pan-cancer cohort of over 25,000 patients with metastatic diseases. By analyzing genomic and clinical data from this cohort, we identified associations between genomic alterations and patterns of metastatic dissemination across 50 tumor types. We found that chromosomal instability is strongly correlated with metastatic burden in some tumor types, including prostate adenocarcinoma, lung adenocarcinoma, and HR+/HER2+ breast ductal carcinoma, but not in others, including colorectal cancer and high-grade serous ovarian cancer, where copy-number alteration patterns may be established early in tumor development. We also identified somatic alterations associated with metastatic burden and specific target organs. Our data offer a valuable resource for the investigation of the biological basis for metastatic spread and highlight the complex role of chromosomal instability in cancer progression.

KRAS Allelic Imbalance Enhances Fitness and Modulates MAP Kinase Dependence in Cancer.

Investigating therapeutic "outliers" that show exceptional responses to anti-cancer treatment can uncover biomarkers of drug sensitivity. We performed preclinical trials investigating primary murine acute myeloid leukemias (AMLs) generated by retroviral insertional mutagenesis in KrasG12D "knockin" mice with the MEK inhibitor PD0325901 (PD901). One outlier AML responded and exhibited intrinsic drug resistance at relapse. Loss of wild-type (WT) Kras enhanced the fitness of the dominant clone and rendered it sensitive to MEK inhibition. Similarly, human colorectal cancer cell lines with increased KRAS mutant allele frequency were more sensitive to MAP kinase inhibition, and CRISPR-Cas9-mediated replacement of WT KRAS with a mutant allele sensitized heterozygous mutant HCT116 cells to treatment. In a prospectively characterized cohort of patients with advanced cancer, 642 of 1,168 (55%) with KRAS mutations exhibited allelic imbalance. These studies demonstrate that serial genetic changes at the Kras/KRAS locus are frequent in cancer and modulate competitive fitness and MEK dependency.

Phase and context shape the function of composite oncogenic mutations.

Cancers develop as a result of driver mutations1,2 that lead to clonal outgrowth and the evolution of disease3,4. The discovery and functional characterization of individual driver mutations are central aims of cancer research, and have elucidated myriad phenotypes5 and therapeutic vulnerabilities6. However, the serial genetic evolution of mutant cancer genes7,8 and the allelic context in which they arise is poorly understood in both common and rare cancer genes and tumour types. Here we find that nearly one in four human tumours contains a composite mutation of a cancer-associated gene, defined as two or more nonsynonymous somatic mutations in the same gene and tumour. Composite mutations are enriched in specific genes, have an elevated rate of use of less-common hotspot mutations acquired in a chronology driven in part by oncogenic fitness, and arise in an allelic configuration that reflects context-specific selective pressures. cis-acting composite mutations are hypermorphic in some genes in which dosage effects predominate (such as TERT), whereas they lead to selection of function in other genes (such as TP53). Collectively, composite mutations are driver alterations that arise from context- and allele-specific selective pressures that are dependent in part on gene and mutation function, and which lead to complex-often neomorphic-functions of biological and therapeutic importance.

Author Correction: Tumour lineage shapes BRCA-mediated phenotypes.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Tumour lineage shapes BRCA-mediated phenotypes.

Mutations in BRCA1 and BRCA2 predispose individuals to certain cancers1-3, and disease-specific screening and preventative strategies have reduced cancer mortality in affected patients4,5. These classical tumour-suppressor genes have tumorigenic effects associated with somatic biallelic inactivation, although haploinsufficiency may also promote the formation and progression of tumours6,7. Moreover, BRCA1/2-mutant tumours are often deficient in the repair of double-stranded DNA breaks by homologous recombination8-13, and consequently exhibit increased therapeutic sensitivity to platinum-containing therapy and inhibitors of poly-(ADP-ribose)-polymerase (PARP)14,15. However, the phenotypic and therapeutic relevance of mutations in BRCA1 or BRCA2 remains poorly defined in most cancer types. Here we show that in the 2.7% and 1.8% of patients with advanced-stage cancer and germline pathogenic or somatic loss-of-function alterations in BRCA1/2, respectively, selective pressure for biallelic inactivation, zygosity-dependent phenotype penetrance, and sensitivity to PARP inhibition were observed only in tumour types associated with increased heritable cancer risk in BRCA1/2 carriers (BRCA-associated cancer types). Conversely, among patients with non-BRCA-associated cancer types, most carriers of these BRCA1/2 mutation types had evidence for tumour pathogenesis that was independent of mutant BRCA1/2. Overall, mutant BRCA is an indispensable founding event for some tumours, but in a considerable proportion of other cancers, it appears to be biologically neutral-a difference predominantly conditioned by tumour lineage-with implications for disease pathogenesis, screening, design of clinical trials and therapeutic decision-making.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer.

PURPOSE: Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell. MATERIALS AND METHODS: We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival. RESULTS: Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high. CONCLUSION: Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.

The evolution of RET inhibitor resistance in RET-driven lung and thyroid cancers.

The efficacy of the highly selective RET inhibitor selpercatinib is now established in RET-driven cancers, and we sought to characterize the molecular determinants of response and resistance. We find that the pre-treatment genomic landscape does not shape the variability of treatment response except for rare instances of RAS-mediated primary resistance. By contrast, acquired selpercatinib resistance is driven by MAPK pathway reactivation by one of two distinct routes. In some patients, on- and off-target pathway reactivation via secondary RET solvent front mutations or MET amplifications are evident. In other patients, rare RET-wildtype tumor cell populations driven by an alternative mitogenic driver are selected for by treatment. Multiple distinct mechanisms are often observed in the same patient, suggesting polyclonal resistance may be common. Consequently, sequential RET-directed therapy may require combination treatment with inhibitors targeting alternative MAPK effectors, emphasizing the need for prospective characterization of selpercatinib-treated tumors at the time of monotherapy progression.

The Genetic Evolution of Treatment-Resistant Cutaneous, Acral, and Uveal Melanomas.

PURPOSE: Melanoma is a biologically heterogeneous disease composed of distinct clinicopathologic subtypes that frequently resist treatment. To explore the evolution of treatment resistance and metastasis, we used a combination of temporal and multilesional tumor sampling in conjunction with whole-exome sequencing of 110 tumors collected from 7 patients with cutaneous (n = 3), uveal (n = 2), and acral (n = 2) melanoma subtypes. EXPERIMENTAL DESIGN: Primary tumors, metastases collected longitudinally, and autopsy tissues were interrogated. All but 1 patient died because of melanoma progression. RESULTS: For each patient, we generated phylogenies and quantified the extent of genetic diversity among tumors, specifically among putative somatic alterations affecting therapeutic resistance. CONCLUSIONS: In 4 patients who received immunotherapy, we found 1-3 putative acquired and intrinsic resistance mechanisms coexisting in the same patient, including mechanisms that were shared by all tumors within each patient, suggesting that future therapies directed at overcoming intrinsic resistance mechanisms may be broadly effective.

The context-specific role of germline pathogenicity in tumorigenesis.

Human cancers arise from environmental, heritable and somatic factors, but how these mechanisms interact in tumorigenesis is poorly understood. Studying 17,152 prospectively sequenced patients with cancer, we identified pathogenic germline variants in cancer predisposition genes, and assessed their zygosity and co-occurring somatic alterations in the concomitant tumors. Two major routes to tumorigenesis were apparent. In carriers of pathogenic germline variants in high-penetrance genes (5.1% overall), lineage-dependent patterns of biallelic inactivation led to tumors exhibiting mechanism-specific somatic phenotypes and fewer additional somatic oncogenic drivers. Nevertheless, 27% of cancers in these patients, and most tumors in patients with pathogenic germline variants in lower-penetrance genes, lacked particular hallmarks of tumorigenesis associated with the germline allele. The dependence of tumors on pathogenic germline variants is variable and often dictated by both penetrance and lineage, a finding with implications for clinical management.

Germ Cell Tumor Molecular Heterogeneity Revealed Through Analysis of Primary and Metastasis Pairs.

PURPOSE: Although primary germ cell tumors (GCTs) have been extensively characterized, molecular analysis of metastatic sites has been limited. We performed whole-exome sequencing and targeted next-generation sequencing on paired primary and metastatic GCT samples in a patient cohort enriched for cisplatin-resistant disease. PATIENTS AND METHODS: Tissue sequencing was performed on 100 tumor specimens from 50 patients with metastatic GCT, and sequencing of plasma cell-free DNA was performed for a subset of patients. RESULTS: The mutational landscape of primary and metastatic pairs from GCT patients was highly discordant (68% of all somatic mutations were discordant). Whereas genome duplication was common and highly concordant between primary and metastatic samples, only 25% of primary-metastasis pairs had ≥ 50% concordance at the level of DNA copy number alterations (CNAs). Evolutionary-based analyses revealed that most mutations arose after CNAs at the respective loci in both primary and metastatic samples, with oncogenic mutations enriched in the set of early-occurring mutations versus variants of unknown significance (VUSs). TP53 pathway alterations were identified in nine cisplatin-resistant patients and had the highest degree of concordance in primary and metastatic specimens, consistent with their association with this treatment-resistant phenotype. CONCLUSION: Analysis of paired primary and metastatic GCT specimens revealed significant molecular heterogeneity for both CNAs and somatic mutations. Among loci demonstrating serial genetic evolution, most somatic mutations arose after CNAs, but oncogenic mutations were enriched in the set of early-occurring mutations as compared with VUSs. Alterations in TP53 were clonal when present and shared among primary-metastasis pairs.

The Evolutionary Origins of Recurrent Pancreatic Cancer.

Surgery is the only curative option for stage I/II pancreatic cancer; nonetheless, most patients will experience a recurrence after surgery and die of their disease. To identify novel opportunities for management of recurrent pancreatic cancer, we performed whole-exome or targeted sequencing of 10 resected primary cancers and matched intrapancreatic recurrences or distant metastases. We identified that recurrent disease after adjuvant or first-line platinum therapy corresponds to an increased mutational burden. Recurrent disease is enriched for genetic alterations predicted to activate MAPK/ERK and PI3K-AKT signaling and develops from a monophyletic or polyphyletic origin. Treatment-induced genetic bottlenecks lead to a modified genetic landscape and subclonal heterogeneity for driver gene alterations in part due to intermetastatic seeding. In 1 patient what was believed to be recurrent disease was an independent (second) primary tumor. These findings suggest routine post-treatment sampling may have value in the management of recurrent pancreatic cancer. SIGNIFICANCE: The biological features or clinical vulnerabilities of recurrent pancreatic cancer after pancreaticoduodenectomy are unknown. Using whole-exome sequencing we find that recurrent disease has a distinct genomic landscape, intermetastatic genetic heterogeneity, diverse clonal origins, and higher mutational burden than found for treatment-naïve disease.See related commentary by Bednar and Pasca di Magliano, p. 762.This article is highlighted in the In This Issue feature, p. 747.

Publisher Correction: The long tail of oncogenic drivers in prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Widespread Selection for Oncogenic Mutant Allele Imbalance in Cancer.

Driver mutations in oncogenes encode proteins with gain-of-function properties that enhance fitness. Heterozygous mutations are thus viewed as sufficient for tumorigenesis. We describe widespread oncogenic mutant allele imbalance in 13,448 prospectively characterized cancers. Imbalance was selected for through modest dosage increases of gain-of-fitness mutations. Negative selection targeted haplo-essential effectors of the spliceosome. Loss of the normal allele comprised a distinct class of imbalance driven by competitive fitness, which correlated with enhanced response to targeted therapies. In many cancers, an antecedent oncogenic mutation drove evolutionarily dependent allele-specific imbalance. In other instances, oncogenic mutations co-opted independent copy-number changes via the evolutionary process of exaptation. Oncogenic allele imbalance is a pervasive evolutionary innovation that enhances fitness and modulates sensitivity to targeted therapy.

Genome doubling shapes the evolution and prognosis of advanced cancers.

Ploidy abnormalities are a hallmark of cancer, but their impact on the evolution and outcomes of cancers is unknown. Here, we identified whole-genome doubling (WGD) in the tumors of nearly 30% of 9,692 prospectively sequenced advanced cancer patients. WGD varied by tumor lineage and molecular subtype, and arose early in carcinogenesis after an antecedent transforming driver mutation. While associated with TP53 mutations, 46% of all WGD arose in TP53-wild-type tumors and in such cases was associated with an E2F-mediated G1 arrest defect, although neither aberration was obligate in WGD tumors. The variability of WGD across cancer types can be explained in part by cancer cell proliferation rates. WGD predicted for increased morbidity across cancer types, including KRAS-mutant colorectal cancers and estrogen receptor-positive breast cancers, independently of established clinical prognostic factors. We conclude that WGD is highly common in cancer and is a macro-evolutionary event associated with poor prognosis across cancer types.

The long tail of oncogenic drivers in prostate cancer.

Comprehensive genomic characterization of prostate cancer has identified recurrent alterations in genes involved in androgen signaling, DNA repair, and PI3K signaling, among others. However, larger and uniform genomic analysis may identify additional recurrently mutated genes at lower frequencies. Here we aggregate and uniformly analyze exome sequencing data from 1,013 prostate cancers. We identify and validate a new class of E26 transformation-specific (ETS)-fusion-negative tumors defined by mutations in epigenetic regulators, as well as alterations in pathways not previously implicated in prostate cancer, such as the spliceosome pathway. We find that the incidence of significantly mutated genes (SMGs) follows a long-tail distribution, with many genes mutated in less than 3% of cases. We identify a total of 97 SMGs, including 70 not previously implicated in prostate cancer, such as the ubiquitin ligase CUL3 and the transcription factor SPEN. Finally, comparing primary and metastatic prostate cancer identifies a set of genomic markers that may inform risk stratification.

The Genomic Landscape of Endocrine-Resistant Advanced Breast Cancers.

We integrated the genomic sequencing of 1,918 breast cancers, including 1,501 hormone receptor-positive tumors, with detailed clinical information and treatment outcomes. In 692 tumors previously exposed to hormonal therapy, we identified an increased number of alterations in genes involved in the mitogen-activated protein kinase (MAPK) pathway and in the estrogen receptor transcriptional machinery. Activating ERBB2 mutations and NF1 loss-of-function mutations were more than twice as common in endocrine resistant tumors. Alterations in other MAPK pathway genes (EGFR, KRAS, among others) and estrogen receptor transcriptional regulators (MYC, CTCF, FOXA1, and TBX3) were also enriched. Altogether, these alterations were present in 22% of tumors, mutually exclusive with ESR1 mutations, and associated with a shorter duration of response to subsequent hormonal therapies.

Accelerating Discovery of Functional Mutant Alleles in Cancer.

Most mutations in cancer are rare, which complicates the identification of therapeutically significant mutations and thus limits the clinical impact of genomic profiling in patients with cancer. Here, we analyzed 24,592 cancers including 10,336 prospectively sequenced patients with advanced disease to identify mutant residues arising more frequently than expected in the absence of selection. We identified 1,165 statistically significant hotspot mutations of which 80% arose in 1 in 1,000 or fewer patients. Of 55 recurrent in-frame indels, we validated that novel AKT1 duplications induced pathway hyperactivation and conferred AKT inhibitor sensitivity. Cancer genes exhibit different rates of hotspot discovery with increasing sample size, with few approaching saturation. Consequently, 26% of all hotspots in therapeutically actionable oncogenes were novel. Upon matching a subset of affected patients directly to molecularly targeted therapy, we observed radiographic and clinical responses. Population-scale mutant allele discovery illustrates how the identification of driver mutations in cancer is far from complete.Significance: Our systematic computational, experimental, and clinical analysis of hotspot mutations in approximately 25,000 human cancers demonstrates that the long right tail of biologically and therapeutically significant mutant alleles is still incompletely characterized. Sharing prospective genomic data will accelerate hotspot identification, thereby expanding the reach of precision oncology in patients with cancer. Cancer Discov; 8(2); 174-83. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 127.

EGFR feedback-inhibition by Ran-binding protein 6 is disrupted in cancer.

Transport of macromolecules through the nuclear pore by importins and exportins plays a critical role in the spatial regulation of protein activity. How cancer cells co-opt this process to promote tumorigenesis remains unclear. The epidermal growth factor receptor (EGFR) plays a critical role in normal development and in human cancer. Here we describe a mechanism of EGFR regulation through the importin ß family member RAN-binding protein 6 (RanBP6), a protein of hitherto unknown functions. We show that RanBP6 silencing impairs nuclear translocation of signal transducer and activator of transcription 3 (STAT3), reduces STAT3 binding to the EGFR promoter, results in transcriptional derepression of EGFR, and increased EGFR pathway output. Focal deletions of the RanBP6 locus on chromosome 9p were found in a subset of glioblastoma (GBM) and silencing of RanBP6 promoted glioma growth in vivo. Our results provide an example of EGFR deregulation in cancer through silencing of components of the nuclear import pathway.

Opposing effects of cancer-type-specific SPOP mutants on BET protein degradation and sensitivity to BET inhibitors.

It is generally assumed that recurrent mutations within a given cancer driver gene elicit similar drug responses. Cancer genome studies have identified recurrent but divergent missense mutations affecting the substrate-recognition domain of the ubiquitin ligase adaptor SPOP in endometrial and prostate cancers. The therapeutic implications of these mutations remain incompletely understood. Here we analyzed changes in the ubiquitin landscape induced by endometrial cancer-associated SPOP mutations and identified BRD2, BRD3 and BRD4 proteins (BETs) as SPOP-CUL3 substrates that are preferentially degraded by endometrial cancer-associated SPOP mutants. The resulting reduction of BET protein levels sensitized cancer cells to BET inhibitors. Conversely, prostate cancer-specific SPOP mutations resulted in impaired degradation of BETs, promoting their resistance to pharmacologic inhibition. These results uncover an oncogenomics paradox, whereby mutations mapping to the same domain evoke opposing drug susceptibilities. Specifically, we provide a molecular rationale for the use of BET inhibitors to treat patients with endometrial but not prostate cancer who harbor SPOP mutations.