Characterizing Mutational Signatures in Human Cancer Cell Lines Reveals Episodic APOBEC Mutagenesis.

Multiple signatures of somatic mutations have been identified in cancer genomes. Exome sequences of 1,001 human cancer cell lines and 577 xenografts revealed most common mutational signatures, indicating past activity of the underlying processes, usually in appropriate cancer types. To investigate ongoing patterns of mutational-signature generation, cell lines were cultured for extended periods and subsequently DNA sequenced. Signatures of discontinued exposures, including tobacco smoke and ultraviolet light, were not generated in vitro. Signatures of normal and defective DNA repair and replication continued to be generated at roughly stable mutation rates. Signatures of APOBEC cytidine deaminase DNA-editing exhibited substantial fluctuations in mutation rate over time with episodic bursts of mutations. The initiating factors for the bursts are unclear, although retrotransposon mobilization may contribute. The examined cell lines constitute a resource of live experimental models of mutational processes, which potentially retain patterns of activity and regulation operative in primary human cancers.

Mitotic tethering enables inheritance of shattered micronuclear chromosomes.

Chromothripsis, the shattering and imperfect reassembly of one (or a few) chromosome(s)1, is an ubiquitous2 mutational process generating localized and complex chromosomal rearrangements that drive genome evolution in cancer. Chromothripsis can be initiated by mis-segregation errors in mitosis3,4 or DNA metabolism5-7 that lead to entrapment of chromosomes within micronuclei and their subsequent fragmentation in the next interphase or following mitotic entry6,8-10. Here we use inducible degrons to demonstrate that chromothriptically produced pieces of a micronucleated chromosome are tethered together in mitosis by a protein complex consisting of mediator of DNA damage checkpoint 1 (MDC1), DNA topoisomerase II-binding protein 1 (TOPBP1) and cellular inhibitor of PP2A (CIP2A), thereby enabling en masse segregation to the same daughter cell. Such tethering is shown to be crucial for the viability of cells undergoing chromosome mis-segregation and shattering after transient inactivation of the spindle assembly checkpoint. Transient, degron-induced reduction in CIP2A following chromosome micronucleation-dependent chromosome shattering is shown to drive acquisition of segmental deletions and inversions. Analyses of pancancer tumour genomes showed that expression of CIP2A and TOPBP1 was increased overall in cancers with genomic rearrangements, including copy number-neutral chromothripsis with minimal deletions, but comparatively reduced in cancers with canonical chromothripsis in which deletions were frequent. Thus, chromatin-bound tethers maintain the proximity of fragments of a shattered chromosome enabling their re-encapsulation into, and religation within, a daughter cell nucleus to form heritable, chromothriptically rearranged chromosomes found in the majority of human cancers.

Chromosomal fragile site breakage by EBV-encoded EBNA1 at clustered repeats.

Epstein-Barr virus (EBV) is an oncogenic herpesvirus associated with several cancers of lymphocytic and epithelial origin1-3. EBV encodes EBNA1, which binds to a cluster of 20 copies of an 18-base-pair palindromic sequence in the EBV genome4-6. EBNA1 also associates with host chromosomes at non-sequence-specific sites7, thereby enabling viral persistence. Here we show that the sequence-specific DNA-binding domain of EBNA1 binds to a cluster of tandemly repeated copies of an EBV-like, 18-base-pair imperfect palindromic sequence encompassing a region of about 21 kilobases at human chromosome 11q23. In situ visualization of the repetitive EBNA1-binding site reveals aberrant structures on mitotic chromosomes characteristic of inherently fragile DNA. We demonstrate that increasing levels of EBNA1 binding trigger dose-dependent breakage at 11q23, producing a fusogenic centromere-containing fragment and an acentric distal fragment, with both mis-segregated into micronuclei in the next cell cycles. In cells latently infected with EBV, elevating EBNA1 abundance by as little as twofold was sufficient to trigger breakage at 11q23. Examination of whole-genome sequencing of EBV-associated nasopharyngeal carcinomas revealed that structural variants are highly enriched on chromosome 11. Presence of EBV is also shown to be associated with an enrichment of chromosome 11 rearrangements across 2,439 tumours from 38 cancer types. Our results identify a previously unappreciated link between EBV and genomic instability, wherein EBNA1-induced breakage at 11q23 triggers acquisition of structural variations in chromosome 11.

Extrachromosomal DNA in the cancerous transformation of Barrett's oesophagus.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

Mechanisms of APOBEC3 mutagenesis in human cancer cells.

The APOBEC3 family of cytosine deaminases has been implicated in some of the most prevalent mutational signatures in cancer1-3. However, a causal link between endogenous APOBEC3 enzymes and mutational signatures in human cancer genomes has not been established, leaving the mechanisms of APOBEC3 mutagenesis poorly understood. Here, to investigate the mechanisms of APOBEC3 mutagenesis, we deleted implicated genes from human cancer cell lines that naturally generate APOBEC3-associated mutational signatures over time4. Analysis of non-clustered and clustered signatures across whole-genome sequences from 251 breast, bladder and lymphoma cancer cell line clones revealed that APOBEC3A deletion diminished APOBEC3-associated mutational signatures. Deletion of both APOBEC3A and APOBEC3B further decreased APOBEC3 mutation burdens, without eliminating them. Deletion of APOBEC3B increased APOBEC3A protein levels, activity and APOBEC3A-mediated mutagenesis in some cell lines. The uracil glycosylase UNG was required for APOBEC3-mediated transversions, whereas the loss of the translesion polymerase REV1 decreased overall mutation burdens. Together, these data represent direct evidence that endogenous APOBEC3 deaminases generate prevalent mutational signatures in human cancer cells. Our results identify APOBEC3A as the main driver of these mutations, indicate that APOBEC3B can restrain APOBEC3A-dependent mutagenesis while contributing its own smaller mutation burdens and dissect mechanisms that translate APOBEC3 activities into distinct mutational signatures.

Mapping clustered mutations in cancer reveals APOBEC3 mutagenesis of ecDNA.

Clustered somatic mutations are common in cancer genomes and previous analyses reveal several types of clustered single-base substitutions, which include doublet- and multi-base substitutions1-5, diffuse hypermutation termed omikli6, and longer strand-coordinated events termed kataegis3,7-9. Here we provide a comprehensive characterization of clustered substitutions and clustered small insertions and deletions (indels) across 2,583 whole-genome-sequenced cancers from 30 types of cancer10. Clustered mutations were highly enriched in driver genes and associated with differential gene expression and changes in overall survival. Several distinct mutational processes gave rise to clustered indels, including signatures that were enriched in tobacco smokers and homologous-recombination-deficient cancers. Doublet-base substitutions were caused by at least 12 mutational processes, whereas most multi-base substitutions were generated by either tobacco smoking or exposure to ultraviolet light. Omikli events, which have previously been attributed to APOBEC3 activity6, accounted for a large proportion of clustered substitutions; however, only 16.2% of omikli matched APOBEC3 patterns. Kataegis was generated by multiple mutational processes, and 76.1% of all kataegic events exhibited mutational patterns that are associated with the activation-induced deaminase (AID) and APOBEC3 family of deaminases. Co-occurrence of APOBEC3 kataegis and extrachromosomal DNA (ecDNA), termed kyklonas (Greek for cyclone), was found in 31% of samples with ecDNA. Multiple distinct kyklonic events were observed on most mutated ecDNA. ecDNA containing known cancer genes exhibited both positive selection and kyklonic hypermutation. Our results reveal the diversity of clustered mutational processes in human cancer and the role of APOBEC3 in recurrently mutating and fuelling the evolution of ecDNA.

Signatures of copy number alterations in human cancer.

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Genome sequencing and analysis of the Tasmanian devil and its transmissible cancer.

The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

Mutational processes molding the genomes of 21 breast cancers.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

Integrated mutational landscape analysis of poorly differentiated high-grade neuroendocrine carcinoma of the uterine cervix.

High-grade neuroendocrine cervical cancers (NETc) are exceedingly rare, highly aggressive tumors. We analyzed 64 NETc tumor samples by whole-exome sequencing (WES). Human papillomavirus DNA was detected in 65.6% (42/64) of the tumors. Recurrent mutations were identified in PIK3CA, KMT2D/MLL2, K-RAS, ARID1A, NOTCH2, and RPL10. The top mutated genes included RB1, ARID1A, PTEN, KMT2D/MLL2, and WDFY3, a gene not yet implicated in NETc. Somatic CNV analysis identified two copy number gains (3q27.1 and 19q13.12) and five copy number losses (1p36.21/5q31.3/6p22.2/9q21.11/11p15.5). Also, gene fusions affecting the ACLY-CRHR1 and PVT1-MYC genes were identified in one of the eight samples subjected to RNA sequencing. To resolve evolutionary history, multiregion WES in NETc admixed with adenocarcinoma cells was performed (i.e., mixed-NETc). Phylogenetic analysis of mixed-NETc demonstrated that adenocarcinoma and neuroendocrine elements derive from a common precursor with mutations typical of adenocarcinomas. Over one-third (22/64) of NETc demonstrated a mutator phenotype of C > T at CpG consistent with deficiencies in MBD4, a member of the base excision repair (BER) pathway. Mutations in the PI3K/AMPK pathways were identified in 49/64 samples. We used two patient-derived-xenografts (PDX) (i.e., NET19 and NET21) to evaluate the activity of pan-HER (afatinib), PIK3CA (copanlisib), and ATR (elimusertib) inhibitors, alone and in combination. PDXs harboring alterations in the ERBB2/PI3K/AKT/mTOR/ATR pathway were sensitive to afatinib, copanlisib, and elimusertib (P < 0.001 vs. controls). However, combinations of copanlisib/afatinib and copanlisib/elimusertib were significantly more effective in controlling NETc tumor growth. These findings define the genetic landscape of NETc and suggest that a large subset of these highly lethal malignancies might benefit from existing targeted therapies.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Author Correction: The evolutionary history of lethal metastatic prostate cancer.

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

Author Correction: Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

Somatic 9p24.1 alterations in HPV- head and neck squamous cancer dictate immune microenvironment and anti-PD-1 checkpoint inhibitor activity.

Somatic copy number alterations (SCNAs), generally (1) losses containing interferons and interferon-pathway genes, many on chromosome 9p, predict immune-cold, immune checkpoint therapy (ICT)-resistant tumors (2); however, genomic regions mediating these effects are unclear and probably tissue specific. Previously, 9p21.3 loss was found to be an early genetic driver of human papillomavirus-negative (HPV-) head and neck squamous cancer (HNSC), associated with an immune-cold tumor microenvironment (TME) signal, and recent evidence suggested that this TME-cold phenotype was greatly enhanced with 9p21 deletion size, notably encompassing band 9p24.1 (3). Here, we report multi-omic, -threshold and continuous-variable dissection of 9p21 and 9p24 loci (including depth and degree of somatic alteration of each band at each locus, and each gene at each band) and TME of four HPV- HNSC cohorts. Preferential 9p24 deletion, CD8 T-cell immune-cold associations were observed, driven by 9p24.1 loss, and in turn by an essential telomeric regulatory gene element, JAK2-CD274. Surprisingly, same genetic region gains were immune hot. Related 9p21-TME analyses were less evident. Inherent 9p-band-level influences on anti-PD1 ICT survival rates, coincident with TME patterns, were also observed. At a 9p24.1 whole-transcriptome expression threshold of 60th percentile, ICT survival rate exceeded that of lower expression percentiles and of chemotherapy; below this transcript threshold, ICT survival was inferior to chemotherapy, the latter unaffected by 9p24.1 expression level (P-values < 0.01, including in a PD-L1 immunohistochemistry-positive patient subgroup). Whole-exome analyses of 10 solid-tumor types suggest that these 9p-related ICT findings could be relevant to squamous cancers, in which 9p24.1 gain/immune-hot associations exist.

Assigning mutational signatures to individual samples and individual somatic mutations with SigProfilerAssignment.

MOTIVATION: Analysis of mutational signatures is a powerful approach for understanding the mutagenic processes that have shaped the evolution of a cancer genome. To evaluate the mutational signatures operative in a cancer genome, one first needs to quantify their activities by estimating the number of mutations imprinted by each signature. RESULTS: Here we present SigProfilerAssignment, a desktop and an online computational framework for assigning all types of mutational signatures to individual samples. SigProfilerAssignment is the first tool that allows both analysis of copy-number signatures and probabilistic assignment of signatures to individual somatic mutations. As its computational engine, the tool uses a custom implementation of the forward stagewise algorithm for sparse regression and nonnegative least squares for numerical optimization. Analysis of 2700 synthetic cancer genomes with and without noise demonstrates that SigProfilerAssignment outperforms four commonly used approaches for assigning mutational signatures. AVAILABILITY AND IMPLEMENTATION: SigProfilerAssignment is available under the BSD 2-clause license at https://github.com/AlexandrovLab/SigProfilerAssignment with a web implementation at https://cancer.sanger.ac.uk/signatures/assignment/.

Arsenic and UVR co-exposure results in unique gene expression profile identifying key co-carcinogenic mechanisms.

Changes in gene expression underlie many pathogenic endpoints including carcinogenesis. Metals, like arsenic, alter gene expression; however, the consequences of co-exposures of metals with other stressors are less understood. Although arsenic acts as a co-carcinogen by enhancing the development of UVR skin cancers, changes in gene expression in arsenic UVR co-carcinogenesis have not been investigated. We performed RNA-sequencing analysis to profile changes in gene expression distinct from arsenic or UVR exposures alone. A large number of differentially expressed genes (DEGs) were identified after arsenic exposure alone, while after UVR exposure alone fewer genes were changed. A distinct increase in the number of DEGs was identified after exposure to combined arsenic and UVR exposure that was synergistic rather than additive. In addition, a majority of these DEGs were unique from arsenic or UVR alone suggesting a distinct response to combined arsenic-UVR exposure. Globally, arsenic alone and arsenic plus UVR exposure caused a global downregulation of genes while fewer genes were upregulated. Gene Ontology analysis using the DEGs revealed cellular processes related to chromosome instability, cell cycle, cellular transformation, and signaling were targeted by combined arsenic and UVR exposure, distinct from UVR alone and arsenic alone, while others were related to epigenetic mechanisms such as the modification of histones. This result suggests the cellular functions we identified in this study may be key in understanding how arsenic enhances UVR carcinogenesis and that arsenic-enhanced gene expression changes may drive co-carcinogenesis of UVR exposure.

Intra-tumour diversification in colorectal cancer at the single-cell level.

Every cancer originates from a single cell. During expansion of the neoplastic cell population, individual cells acquire genetic and phenotypic differences from each other. Here, to investigate the nature and extent of intra-tumour diversification, we characterized organoids derived from multiple single cells from three colorectal cancers as well as from adjacent normal intestinal crypts. Colorectal cancer cells showed extensive mutational diversification and carried several times more somatic mutations than normal colorectal cells. Most mutations were acquired during the final dominant clonal expansion of the cancer and resulted from mutational processes that are absent from normal colorectal cells. Intra-tumour diversification of DNA methylation and transcriptome states also occurred; these alterations were cell-autonomous, stable, and followed the phylogenetic tree of each cancer. There were marked differences in responses to anticancer drugs between even closely related cells of the same tumour. The results indicate that colorectal cancer cells experience substantial increases in somatic mutation rate compared to normal colorectal cells, and that genetic diversification of each cancer is accompanied by pervasive, stable and inherited differences in the biological states of individual cancer cells.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.