The pan-cancer analysis of whole genomes The expansion of whole-genome sequencing studies from individual ICGC and TCGA working groups presented the opportunity to undertake a meta-analysis of genomic features across tumour types. To achieve this, the PCAWG Consortium was established. A Technical Working Group implemented the informatics analyses by aggregating the raw sequencing data from different working groups that studied individual tumour types, aligning the sequences to the human genome and delivering a set of high-quality somatic mutation calls for downstream analysis (Extended Data Fig. 1). Given the recent meta-analysis

Acquired mutations are pervasive across normal tissues. However, our understanding of the processes that drive transformation of certain clones to cancer is limited. Here we study this phenomenon in the context of clonal hematopoiesis (CH) and the development of therapy-related myeloid neoplasms (tMN). We find mutations are selected differentially based on exposures. Mutations in ASXL1 are enriched in current or former smokers, whereas cancer therapy with radiation, platinum and topoisomerase II inhibitors preferentially selects for mutations in DNA damage response (DDR) genes ( TP53, PPM1D, CHEK2 ). Sequential sampling provides definitive evidence that DDR clones outcompete other clones when exposed to certain therapies. Among cases where CH was previously detected, the CH mutation was present at tMN diagnosis. We identify the molecular characteristics of CH that increase risk of tMN. The increasing implementation of clinical sequencing at diagnosis provides an opportunity to identify patients at risk of tMN for prevention strategies.

The multiple myeloma (MM) genome is heterogeneous and evolves through preclinical and post-diagnosis phases. Here we report a catalog and hierarchy of driver lesions using sequences from 67 MM genomes serially collected from 30 patients together with public exome datasets. Bayesian clustering defines at least 7 genomic subgroups with distinct sets of co-operating events. Focusing on whole genome sequencing data, complex structural events emerge as major drivers, including chromothripsis and a novel replication-based mechanism of templated insertions, which typically occur early. Hyperdiploidy also occurs early, with individual trisomies often acquired in different chronological windows during evolution, and with a preferred order of acquisition. Conversely, positively selected point mutations, whole genome duplication and chromoplexy events occur in later disease phases. Thus, initiating driver events, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of evolution that are biologically relevant but heterogeneous across patients.

Summary Whole genome sequencing (WGS) brings comprehensive insights to cancer genome interpretation. To explore clinical value of WGS, we sequenced 254 triple negative breast cancers (TNBC) with associated treatment and outcome data collected between 2010-2015 via the population-based Sweden Cancerome Analysis Network-Breast (SCAN-B) project (ClinicalTrials.gov ID:NCT02306096). Applying the HRDetect mutational-signature-based algorithm to classify tumors, 59% were predicted to have Homologous-recombination-repair deficiency (HRDetect-high): 67% explained by germline/somatic mutations of BRCA1/BRCA2, BRCA1 promoter hypermethylation, RAD51C hypermethylation or biallelic loss of PALB2. A novel mechanism of BRCA1 abrogation was discovered via germline SINE-VNTR-Alu retrotransposition. HRDetect provided independent prognostic information, with HRDetect-high patients having better outcome on adjuvant chemotherapy for invasive disease-free survival (Hazard Ratio, HR=0.42, 95% confidence interval, CI=0.2-0.87), and distant relapse-free interval (HR=0.31, CI=0.13-0.76) compared to HRDetect-low, regardless of whether a genetic/epigenetic cause was identified. HRDetect-intermediate, some possessing potentially targetable biological abnormalities, had poorest outcomes. HRDetect-low cancers also had inadequate outcomes: ~4.7% were mismatch-repair-deficient - another targetable defect, not typically sought; and was enriched for (but not restricted to) PIK3CA/AKT1 pathway abnormalities. New treatment options need to be considered for now-discernible HRDetect-intermediate and HRDetect-low categories. This population-based study advocates for WGS of TNBC to better inform trial stratification and improve clinical decision-making.

Mutational signatures are patterns of mutations that arise during tumorigenesis. We present an enhanced, practical framework for mutational signature analyses. Applying these methods on 3,107 whole genome sequenced (WGS) primary cancers of 21 organs reveals known signatures and nine previously undescribed rearrangement signatures. We highlight inter-organ variability of

Multiple myeloma (MM) has a heterogeneous genome, evolving through both pre-clinical and post-diagnosis phases. Here, using sequences from 67 MM genomes serially collected from 30 patients together with public datasets, we establish a hierarchy of driver lesions. Point mutations, structural variants and copy number aberrations define at least 7 genomic subgroups of MM, each with distinct sets of co-operating driver mutations. Complex structural events are major drivers of MM, including chromothripsis, chromoplexy and a replication-based mechanism of templated insertions: these typically occur early. Hyperdiploidy also occurs early, with individual chromosomes often gained in more than one chronological epoch of MM evolution, showing a preferred order of acquisition. Positively selected point mutations frequently occur in later phases of disease development, as do structural variants involving MYC. Thus, initiating driver events of MM, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of subsequent evolution. MFAG (n.17658).

The landscape of structural variants (SV) in multiple myeloma remains poorly understood. Here, we performed comprehensive analysis of SVs in a large cohort of 752 patients with multiple myeloma by low-coverage long-insert whole-genome sequencing. We identified 68 SV hotspots involving 17 new candidate driver genes, including the therapeutic targets BCMA (TNFRSF17), SLAM7, and MCL1. Catastrophic complex rearrangements termed chromothripsis were present in 24% of patients and independently associated with poor clinical outcomes. Templated insertions were the second most frequent complex event (19%), mostly involved in super-enhancer hijacking and activation of oncogenes such as CCND1 and MYC. Importantly, in 31% of patients, two or more seemingly independent putative driver events were caused by a single structural event, demonstrating that the complex genomic landscape of multiple myeloma can be acquired through few key events during tumor evolutionary history. Overall, this study reveals the critical role of SVs in multiple myeloma pathogenesis. Significance: Previous genomic studies in multiple myeloma have largely focused on single-nucleotide variants, recurrent copy-number alterations, and recurrent translocations. Here, we demonstrate the crucial role of SVs and complex events in the development of multiple myeloma and highlight the importance of whole-genome sequencing to decipher its genomic complexity. See related commentary by Bergsagel and Kuehl, p. 221. This article is highlighted in the In This Issue feature, p. 215

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/DAXX trunc) 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/DAXX trunc 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.