Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Framework for quality assessment of whole genome cancer sequences.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .

Genomic footprints of activated telomere maintenance mechanisms in cancer.

Cancers require telomere maintenance mechanisms for unlimited replicative potential. They achieve this through TERT activation or alternative telomere lengthening associated with ATRX or DAXX loss. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, we dissect whole-genome sequencing data of over 2500 matched tumor-control samples from 36 different tumor types aggregated within the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium to characterize the genomic footprints of these mechanisms. While the telomere content of tumors with ATRX or DAXX mutations (ATRX/DAXXtrunc) is increased, tumors with TERT modifications show a moderate decrease of telomere content. One quarter of all tumor samples contain somatic integrations of telomeric sequences into non-telomeric DNA. This fraction is increased to 80% prevalence in ATRX/DAXXtrunc tumors, which carry an aberrant telomere variant repeat (TVR) distribution as another genomic marker. The latter feature includes enrichment or depletion of the previously undescribed singleton TVRs TTCGGG and TTTGGG, respectively. Our systematic analysis provides new insight into the recurrent genomic alterations associated with telomere maintenance mechanisms in cancer.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Defective homologous recombination DNA repair as therapeutic target in advanced chordoma.

Chordomas are rare bone tumors with few therapeutic options. Here we show, using whole-exome and genome sequencing within a precision oncology program, that advanced chordomas (n = 11) may be characterized by genomic patterns indicative of defective homologous recombination (HR) DNA repair and alterations affecting HR-related genes, including, for example, deletions and pathogenic germline variants of BRCA2, NBN, and CHEK2. A mutational signature associated with HR deficiency was significantly enriched in 72.7% of samples and co-occurred with genomic instability. The poly(ADP-ribose) polymerase (PARP) inhibitor olaparib, which is preferentially toxic to HR-incompetent cells, led to prolonged clinical benefit in a patient with refractory chordoma, and whole-genome analysis at progression revealed a PARP1 p.T910A mutation predicted to disrupt the autoinhibitory PARP1 helical domain. These findings uncover a therapeutic opportunity in chordoma that warrants further exploration, and provide insight into the mechanisms underlying PARP inhibitor resistance.

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma.

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

Enhancer hijacking activates oncogenic transcription factor NR4A3 in acinic cell carcinomas of the salivary glands.

The molecular pathogenesis of salivary gland acinic cell carcinoma (AciCC) is poorly understood. The secretory Ca-binding phosphoprotein (SCPP) gene cluster at 4q13 encodes structurally related phosphoproteins of which some are specifically expressed at high levels in the salivary glands and constitute major components of saliva. Here we report on recurrent rearrangements [t(4;9)(q13;q31)] in AciCC that translocate active enhancer regions from the SCPP gene cluster to the region upstream of Nuclear Receptor Subfamily 4 Group A Member 3 (NR4A3) at 9q31. We show that NR4A3 is specifically upregulated in AciCCs, and that active chromatin regions and gene expression signatures in AciCCs are highly correlated with the NR4A3 transcription factor binding motif. Overexpression of NR4A3 in mouse salivary gland cells increases expression of known NR4A3 target genes and has a stimulatory functional effect on cell proliferation. We conclude that NR4A3 is upregulated through enhancer hijacking and has important oncogenic functions in AciCC.

The mutational landscape of Burkitt-like lymphoma with 11q aberration is distinct from that of Burkitt lymphoma.

The new recently described provisional lymphoma category Burkitt-like lymphoma with 11q aberration comprises cases similar to Burkitt lymphoma (BL) on morphological, immunophenotypic and gene-expression levels but lacking the IG-MYC translocation. They are characterized by a peculiar imbalance pattern on chromosome 11, but the landscape of mutations is not yet described. Thus, we investigated 15 MYC-negative Burkitt-like lymphoma with 11q aberration (mnBLL,11q,) cases by copy-number analysis and whole-exome sequencing. We refined the regions of 11q imbalance and identified the INO80 complex-associated gene NFRKB as a positional candidate in 11q24.3. Next to recurrent gains in 12q13.11-q24.32 and 7q34-qter as well as losses in 13q32.3-q34, we identified 47 genes recurrently affected by protein-changing mutations (each ≥3 of 15 cases). Strikingly, we did not detect recurrent mutations in genes of the ID3-TCF3 axis or the SWI/SNF complex that are frequently altered in BL, or in genes frequently mutated in germinal center-derived B-cell lymphomas like KMT2D or CREBBP An exception is GNA13, which was mutated in 7 of 15 cases. We conclude that the genomic landscape of mnBLL,11q, differs from that of BL both at the chromosomal and mutational levels. Our findings implicate that mnBLL,11q, is a lymphoma category distinct from BL at the molecular level.

Whole genome sequencing puts forward hypotheses on metastasis evolution and therapy in colorectal cancer.

Incomplete understanding of the metastatic process hinders personalized therapy. Here we report the most comprehensive whole-genome study of colorectal metastases vs. matched primary tumors. 65% of somatic mutations originate from a common progenitor, with 15% being tumor- and 19% metastasis-specific, implicating a higher mutation rate in metastases. Tumor- and metastasis-specific mutations harbor elevated levels of BRCAness. We confirm multistage progression with new components ARHGEF7/ARHGEF33. Recurrently mutated non-coding elements include ncRNAs RP11-594N15.3, AC010091, SNHG14, 3' UTRs of FOXP2, DACH2, TRPM3, XKR4, ANO5, CBL, CBLB, the latter four potentially dual protagonists in metastasis and efferocytosis-/PD-L1 mediated immunosuppression. Actionable metastasis-specific lesions include FAT1, FGF1, BRCA2, KDR, and AKT2-, AKT3-, and PDGFRA-3' UTRs. Metastasis specific mutations are enriched in PI3K-Akt signaling, cell adhesion, ECM and hepatic stellate activation genes, suggesting genetic programs for site-specific colonization. Our results put forward hypotheses on tumor and metastasis evolution, and evidence for metastasis-specific events relevant for personalized therapy.

IG-MYC + neoplasms with precursor B-cell phenotype are molecularly distinct from Burkitt lymphomas.

The WHO Classification of Tumours of Haematopoietic and Lymphoid Tissue notes instances of Burkitt lymphoma/leukemia (BL) with IG-MYC rearrangement displaying a B-cell precursor immunophenotype (termed herein "preBLL"). To characterize the molecular pathogenesis of preBLL, we investigated 13 preBLL cases (including 1 cell line), of which 12 were analyzable using genome, exome, and targeted sequencing, imbalance mapping, and DNA methylation profiling. In 5 patients with reads across the IG-MYC breakpoint junctions, we found evidence that the translocation derived from an aberrant VDJ recombination, as is typical for IG translocations arising in B-cell precursors. Genomic changes like biallelic IGH translocations or VDJ rearrangements combined with translocation into the VDJ region on the second allele, potentially preventing expression of a productive immunoglobulin, were detected in 6 of 13 cases. We did not detect mutations in genes frequently altered in BL, but instead found activating NRAS and/or KRAS mutations in 7 of 12 preBLLs. Gains on 1q, recurrent in BL and preB lymphoblastic leukemia/lymphoma (pB-ALL/LBL), were detected in 7 of 12 preBLLs. DNA methylation profiling showed preBLL to cluster with precursor B cells and pB-ALL/LBL, but apart from BL. We conclude that preBLL genetically and epigenetically resembles pB-ALL/LBL rather than BL. Therefore, we propose that preBLL be considered as a pB-ALL/LBL with recurrent genetic abnormalities.

Author Correction: The landscape of genomic alterations across childhood cancers.

In this Article, author Benedikt Brors was erroneously associated with affiliation number '8' (Department of Developmental Neurobiology, St Jude Children's Research Hospital, Memphis, Tennessee, USA); the author's two other affiliations (affiliations '3' and '7', both at the German Cancer Research Center (DKFZ)) were correct. This has been corrected online.

NRG1 Fusions in KRAS Wild-Type Pancreatic Cancer.

We used whole-genome and transcriptome sequencing to identify clinically actionable genomic alterations in young adults with pancreatic ductal adenocarcinoma (PDAC). Molecular characterization of 17 patients with PDAC enrolled in a precision oncology program revealed gene fusions amenable to pharmacologic inhibition by small-molecule tyrosine kinase inhibitors in all patients with KRAS wild-type (KRASWT) tumors (4 of 17). These alterations included recurrent NRG1 rearrangements predicted to drive PDAC development through aberrant ERBB receptor-mediated signaling, and pharmacologic ERBB inhibition resulted in clinical improvement and remission of liver metastases in 2 patients with NRG1-rearranged tumors that had proved resistant to standard treatment. Our findings demonstrate that systematic screening of KRASWT tumors for oncogenic fusion genes will substantially improve the therapeutic prospects for a sizeable fraction of patients with PDAC.Significance: Advanced PDAC is a malignancy with few treatment options that lacks molecular mechanism-based therapies. Our study uncovers recurrent gene rearrangements such as NRG1 fusions as disease-driving events in KRASwt tumors, thereby providing novel insights into oncogenic signaling and new therapeutic options in this entity. Cancer Discov; 8(9); 1087-95. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Genomic features of renal cell carcinoma with venous tumor thrombus.

A venous tumor thrombus (VTT) is a potentially lethal complication of renal cell carcinoma (RCC) but virtually nothing is known about the underlying natural history. Based on our observation that venous thrombi contain significant numbers of viable tumor cells, we applied multiregion whole exome sequencing to a total of 37 primary tumor and VTT samples including normal tissue specimens from five consecutive patients. Our findings demonstrate mutational heterogeneity between primary tumor and VTT with 106 of 483 genes (22%) harboring functional SNVs and/or indels altered in either primary tumor or thrombus. Reconstruction of the clonal phylogeny showed clustering of tumor samples and VTT samples, respectively, in the majority of tumors. However, no new subclones were detected suggesting that pre-existing subclones of the primary tumor drive VTT formation. Importantly, we found several lines of evidence for "BRCAness" in a subset of tumors. These included mutations in genes that confer "BRCAness", a mutational signature and an increase of small indels. Re-analysis of SNV calls from the TCGA KIRC-US cohort confirmed a high frequency of the "BRCAness" mutational signature AC3 in clear cell RCC. Our findings warrant further pre-clinical experiments and may lead to novel personalized therapies for RCC patients.

The landscape of genomic alterations across childhood cancers.

Pan-cancer analyses that examine commonalities and differences among various cancer types have emerged as a powerful way to obtain novel insights into cancer biology. Here we present a comprehensive analysis of genetic alterations in a pan-cancer cohort including 961 tumours from children, adolescents, and young adults, comprising 24 distinct molecular types of cancer. Using a standardized workflow, we identified marked differences in terms of mutation frequency and significantly mutated genes in comparison to previously analysed adult cancers. Genetic alterations in 149 putative cancer driver genes separate the tumours into two classes: small mutation and structural/copy-number variant (correlating with germline variants). Structural variants, hyperdiploidy, and chromothripsis are linked to TP53 mutation status and mutational signatures. Our data suggest that 7-8% of the children in this cohort carry an unambiguous predisposing germline variant and that nearly 50% of paediatric neoplasms harbour a potentially druggable event, which is highly relevant for the design of future clinical trials.

Integrative genomic and transcriptomic analysis of leiomyosarcoma.

Leiomyosarcoma (LMS) is an aggressive mesenchymal malignancy with few therapeutic options. The mechanisms underlying LMS development, including clinically actionable genetic vulnerabilities, are largely unknown. Here we show, using whole-exome and transcriptome sequencing, that LMS tumors are characterized by substantial mutational heterogeneity, near-universal inactivation of TP53 and RB1, widespread DNA copy number alterations including chromothripsis, and frequent whole-genome duplication. Furthermore, we detect alternative telomere lengthening in 78% of cases and identify recurrent alterations in telomere maintenance genes such as ATRX, RBL2, and SP100, providing insight into the genetic basis of this mechanism. Finally, most tumors display hallmarks of "BRCAness", including alterations in homologous recombination DNA repair genes, multiple structural rearrangements, and enrichment of specific mutational signatures, and cultured LMS cells are sensitive towards olaparib and cisplatin. This comprehensive study of LMS genomics has uncovered key biological features that may inform future experimental research and enable the design of novel therapies.

Screening drug effects in patient-derived cancer cells links organoid responses to genome alterations.

Cancer drug screening in patient-derived cells holds great promise for personalized oncology and drug discovery but lacks standardization. Whether cells are cultured as conventional monolayer or advanced, matrix-dependent organoid cultures influences drug effects and thereby drug selection and clinical success. To precisely compare drug profiles in differently cultured primary cells, we developed DeathPro, an automated microscopy-based assay to resolve drug-induced cell death and proliferation inhibition. Using DeathPro, we screened cells from ovarian cancer patients in monolayer or organoid culture with clinically relevant drugs. Drug-induced growth arrest and efficacy of cytostatic drugs differed between the two culture systems. Interestingly, drug effects in organoids were more diverse and had lower therapeutic potential. Genomic analysis revealed novel links between drug sensitivity and DNA repair deficiency in organoids that were undetectable in monolayers. Thus, our results highlight the dependency of cytostatic drugs and pharmacogenomic associations on culture systems, and guide culture selection for drug tests.