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

Genomic, transcriptomic, and viral integration profiles associated with recurrent/metastatic progression in high-risk human papillomavirus cervical carcinomas

Acquisition of recurrent/metastatic potential by a tumor cell defines a critical step in malignant progression. However, understanding of metastatic progression at the molecular level is scarce for cervical carcinomas (CES). In this study, we performed genomic, transcriptomic, and viral profiling of five pairs of primary (CES‐P) and matched recurrent/metastatic tumors (CES‐R/M) with high risk human papillomavirus. Whole exome sequencing revealed mutation features of CES‐R/M including elevated mutation burdens and prevalent copy number alterations compared to their matched CES‐P. A relative deficit of APOBEC‐related mutation signatures accompanying the transcriptional downregulation of APOBEC3A was observed for CES‐R/M. Mutations in genes encoding epigenetic regulators were commonly observed as CES‐R/M‐specific alterations. Immunoprofiling and gene set analysis revealed CES‐Ps were enriched with transcripts representing activated anticancer immunity such as interferon‐gamma pathway, while CES‐R/M exhibited upregulation of genes involved in epithelial‐mesenchymal transition and angiogenesis. Viral capture sequencing revealed that integration sites remained enriched in viral E1 protein domain during malignant progression. Moreover, we found transcriptional upregulation of POSTN and downregulation of APOBEC3A were associated with unfavorable clinical outcomes in CES. Comprehensive genomic and transcriptomic profiling of a rare cohort including CES‐R/M identified metastases‐specific features to advance the molecular understanding into CES metastatic progression with potential clinical implications.

Molecular origins of APOBEC-associated mutations in cancer

The APOBEC family of cytidine deaminases has been proposed to represent a major enzymatic source of mutations in cancer. Here, we summarize available evidence that links APOBEC deaminases to cancer mutagenesis. We also highlight newly identified human cell models of APOBEC mutagenesis, including cancer cell lines with suspected endogenous APOBEC activity and a cell system of telomere crisis-associated mutations. Finally, we draw on recent data to propose potential causes of APOBEC misregulation in cancer, including the instigating factors, the relevant mutator(s), and the mechanisms underlying generation of the genome-dispersed and clustered APOBEC-induced mutations.

Single-nucleotide polymorphism of the DNA cytosine deaminase APOBEC3H haplotype I leads to enzyme destabilization and correlates with lung cancer

A number of APOBEC family DNA cytosine deaminases can induce mutations in tumor cells. APOBEC3H haplotype I is one of the deaminases that has been proposed to cause mutations in lung cancer. Here, we confirmed that APOBEC3H haplotype I can cause uracil-induced DNA damage in lung cancer cells that results in γH2AX foci. Interestingly, the database of cancer biomarkers in DNA repair genes (DNArCdb) identified a single-nucleotide polymorphism (rs139298) of APOBEC3H haplotype I that is involved in lung cancer. While we thought this may increase the activity of APOBEC3H haplotype I, instead we found through computational modeling and cell-based experiments that this single-nucleotide polymorphism causes the destabilization of APOBEC3H Haplotype I. Computational analysis suggests that the resulting K121E change affects the structure of APOBEC3H leading to active site disruption and destabilization of the RNA-mediated dimer interface. A K117E mutation in a K121E background stabilized the APOBEC3H haplotype I, thus enabling biochemical study. Subsequent analysis showed that K121E affected catalytic activity, single-stranded DNA binding and oligomerization on single-stranded DNA. The destabilization of a DNA mutator associated with lung cancer supports the model that too much APOBEC3-induced mutation could result in immune recognition or death of tumor cells.

DNA mismatch repair promotes APOBEC3-mediated diffuse hypermutation in human cancers

Certain mutagens, including the APOBEC3 (A3) cytosine deaminase enzymes, can create multiple genetic changes in a single event. Activity of A3s results in striking 'mutation showers' occurring near DNA breakpoints; however, less is known about the mechanisms underlying the majority of A3 mutations. We classified the diverse patterns of clustered mutagenesis in tumor genomes, which identified a new A3 pattern: nonrecurrent, diffuse hypermutation (omikli). This mechanism occurs independently of the known focal hypermutation (kataegis), and is associated with activity of the DNA mismatch-repair pathway, which can provide the single-stranded DNA substrate needed by A3, and contributes to a substantial proportion of A3 mutations genome wide. Because mismatch repair is directed towards early-replicating, gene-rich chromosomal domains, A3 mutagenesis has a high propensity to generate impactful mutations, which exceeds that of other common carcinogens such as tobacco smoke and ultraviolet exposure. Cells direct their DNA repair capacity towards more important genomic regions; thus, carcinogens that subvert DNA repair can be remarkably potent.

The LL-100 Cell Lines Panel: Tool for Molecular Leukemia-Lymphoma Research

Certified cell line models provide ideal experimental platforms to answer countless scientific questions. The LL-100 panel is a cohort of cell lines that are broadly representative of all leukemia–lymphoma entities (including multiple myeloma and related diseases), rigorously authenticated and validated, and comprehensively annotated. The process of the assembly of the LL-100 panel was based on evidence and experience. To expand the genetic characterization across all LL-100 cell lines, we performed whole-exome sequencing and RNA sequencing. Here, we describe the conception of the panel and showcase some exemplary applications with a focus on cancer genomics. Due diligence was paid to exclude cross-contaminated and non-representative cell lines. As the LL-100 cell lines are so well characterized and readily available, the panel will be a valuable resource for identifying cell lines with mutations in cancer genes, providing superior model systems. The data also add to the current knowledge of the molecular pathogenesis of leukemia–lymphoma. Additional efforts to expand the breadth of available high-quality cell lines are clearly warranted.

Human genes differ by their UV sensitivity estimated through analysis of UV-induced silent mutations in melanoma

We hypothesized that human genes differ by their sensitivity to ultraviolet (UV) exposure. We used somatic mutations detected by genome-wide screens in melanoma and reported in the Catalog Of Somatic Mutations In Cancer. As a measure of UV sensitivity, we used the number of silent mutations generated by C>T transitions in pyrimidine dimers of a given transcript divided by the number of potential sites for this type of mutations in the transcript. We found that human genes varied by UV sensitivity by two orders of magnitude. We noted that the melanoma-associated tumor suppressor gene CDKN2A was among the top five most UV-sensitive genes in the human genome. Melanoma driver genes have a higher UV-sensitivity compared with other genes in the human genome. The difference was more prominent for tumor suppressors compared with oncogene. The results of this study suggest that differential sensitivity of human transcripts to UV light may explain melanoma specificity of some driver genes. Practical significance of the study relates to the fact that differences in UV sensitivity among human genes need to be taken into consideration whereas predicting melanoma-associated genes by the number of somatic mutations detected in a given gene.

Pervasive lesion segregation shapes cancer genome evolution

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

The duality of human oncoproteins: drivers of cancer and congenital disorders

Human oncoproteins promote transformation of cells into tumours by dysregulating the signalling pathways that are involved in cell growth, proliferation and death. Although oncoproteins were discovered many years ago and have been widely studied in the context of cancer, the recent use of high-throughput sequencing techniques has led to the identification of cancer-associated mutations in other conditions, including many congenital disorders. These syndromes offer an opportunity to study oncoprotein signalling and its biology in the absence of additional driver or passenger mutations, as a result of their monogenic nature. Moreover, their expression in multiple tissue lineages provides insight into the biology of the proto-oncoprotein at the physiological level, in both transformed and unaffected tissues. Given the recent paradigm shift in regard to how oncoproteins promote transformation, we review the fundamentals of genetics, signalling and pathogenesis underlying oncoprotein duality.

Matching cell lines with cancer type and subtype of origin via mutational, epigenomic, and transcriptomic patterns

Cell lines are commonly used as cancer models. The tissue of origin provides context for understanding biological mechanisms and predicting therapy response. We therefore systematically examined whether cancer cell lines exhibit features matching the presumed cancer type of origin. Gene expression and DNA methylation classifiers trained on ~9000 tumors identified 35 (of 614 examined) cell lines that better matched a different tissue or cell type than the one originally assigned. Mutational patterns further supported most reassignments. For instance, cell lines identified as originating from the skin often exhibited a UV mutational signature. We cataloged 366 "golden set" cell lines in which transcriptomic and epigenomic profiles strongly resemble the cancer type of origin, further proposing their assignments to subtypes. Accounting for the uncertain tissue of origin in cell line panels can change the interpretation of drug screening and genetic screening data, revealing previously unknown genomic determinants of sensitivity or resistance.

Hypermutation in single-stranded DNA

Regions of genomic DNA can become single-stranded in the course of normal replication and transcription as well as during DNA repair. Abnormal repair and replication intermediates can contain large stretches of persistent single-stranded DNA, which is extremely vulnerable to DNA damaging agents and hypermutation. Since such single-stranded DNA spans only a fraction of the genome at a given instance, hypermutation in these regions leads to tightly-spaced mutation clusters. This phenomenon of hypermutation in single-stranded DNA has been documented in several experimental models as well as in cancer genomes. Recently, hypermutated single-stranded RNA viral genomes also have been documented. Moreover, indications of hypermutation in single-stranded DNA may also be found in the human germline. This review will summarize key current knowledge and the recent developments in understanding the diverse mechanisms and sources of ssDNA hypermutation.

The mutREAD method detects mutational signatures from low quantities of cancer DNA

Mutational processes acting on cancer genomes can be traced by investigating mutational signatures. Because high sequencing costs limit current studies to small numbers of good-quality samples, we propose a robust, cost- and time-effective method, called mutREAD, to detect mutational signatures from small quantities of DNA, including degraded samples. We show that mutREAD recapitulates mutational signatures identified by whole genome sequencing, and will ultimately allow the study of mutational signatures in larger cohorts and, by compatibility with formalin-fixed paraffin-embedded samples, in clinical settings.

MTOR signaling orchestrates stress-induced mutagenesis, facilitating adaptive evolution in cancer

In microorganisms, evolutionarily conserved mechanisms facilitate adaptation to harsh conditions through stress-induced mutagenesis (SIM). Analogous processes may underpin progression and therapeutic failure in human cancer. We describe SIM in multiple in vitro and in vivo models of human cancers under nongenotoxic drug selection, paradoxically enhancing adaptation at a competing intrinsic fitness cost. A genome-wide approach identified the mechanistic target of rapamycin (MTOR) as a stress-sensing rheostat mediating SIM across multiple cancer types and conditions. These observations are consistent with a two-phase model for drug resistance, in which an initially rapid expansion of genetic diversity is counterbalanced by an intrinsic fitness penalty, subsequently normalizing to complete adaptation under the new conditions. This model suggests synthetic lethal strategies to minimize resistance to anticancer therapy.

POLE Mutation Spectra Are Shaped by the Mutant Allele Identity, Its Abundance, and Mismatch Repair Status

Human tumors with exonuclease domain mutations in the gene encoding DNA polymerase ε (POLE) have incredibly high mutation burdens. These errors arise in four unique mutation signatures occurring in different relative amounts, the etiologies of which remain poorly understood. We used CRISPR-Cas9 to engineer human cell lines expressing POLE tumor variants, with and without mismatch repair (MMR). Whole-exome sequencing of these cells after defined numbers of population doublings permitted analysis of nascent mutation accumulation. Unlike an exonuclease active site mutant that we previously characterized, POLE cancer mutants readily drive signature mutagenesis in the presence of functional MMR. Comparison of cell line and human patient data suggests that the relative abundance of mutation signatures partitions POLE tumors into distinct subgroups dependent on the nature of the POLE allele, its expression level, and MMR status. These results suggest that different POLE mutants have previously unappreciated differences in replication fidelity and mutagenesis.

Distinct immune evasion in APOBEC-enriched, HPV-negative HNSCC

Immune checkpoint inhibition leads to response in some patients with head and neck squamous cell carcinoma (HNSCC). Robust biomarkers are lacking to date. We analyzed viral status, gene expression signatures, mutational load and mutational signatures in whole exome and RNA-sequencing data of the HNSCC TCGA dataset (n = 496) and a validation set (DKTK MASTER cohort, n = 10). Public single-cell gene expression data from 17 HPV-negative HNSCC were separately reanalyzed. APOBEC3-associated TCW motif mutations but not total single nucleotide variant burden were significantly associated with inflammation. This association was restricted to HPV-negative HNSCC samples. An APOBEC-enriched, HPV-negative subgroup was identified, that showed higher T-cell inflammation and immune checkpoint expression, as well as expression of APOBEC3 genes. Mutations in immune-evasion pathways were also enriched in these tumors. Analysis of single-cell sequencing data identified expression of APOBEC3B and 3C genes in malignant cells. We identified an APOBEC-enriched subgroup of HPV-negative HNSCC with a distinct immunogenic phenotype, potentially mediating response to immunotherapy.

Vulvar Squamous Cell Carcinoma: Comprehensive Genomic Profiling of HPV+ Versus HPV- Forms Reveals Distinct Sets of Potentially Actionable Molecular Targets

PURPOSE

Vulvar squamous cell carcinoma (vSCC) encompasses two predominant variants: one associated with detectable high-risk strains of human papillomavirus (hrHPV) and a second form often occurring in the context of chronic dermatitis in postmenopausal women. Genomic assessment of a large-scale cohort of patients with aggressive vSCC may identify distinct mutational signatures.

MATERIALS AND METHODS

Tumor samples from a total of 280 patients with vSCC underwent hybridization capture with analysis of up to 406 cancer-related genes. Human papillomavirus (HPV) sequences were detected by de novo assembly of nonhuman sequencing reads and aligned to the RefSeq database. Immunohistochemistry for programmed death-ligand 1 (PD-L1) was assessed.

RESULTS

One hundred two of 280 vSCCs (36%) contained hrHPV sequences, predominantly HPV 16 (88%). The HPV-positive (HPV+) group was significantly younger (median age, 59 v 64 years; P = .001). Compared with HPV-negative (HPV–) vSCCs, HPV+ tumors showed more frequent pathogenic alterations in PIK3CA (31% v 16%; P = .004), PTEN (14% v 2%; P < .0001), EP300 (14% v 1%; P < .0001), STK11 (14% v 1%; P < .0001), AR (5% v 0%; P = .006), and FBXW7 (10% v 3%; P = .03). In contrast, HPV– vSCCs showed more alterations in TP53 (83% v 6%; P < .0001), TERTp (71% v 9%; P < .0001), CDKN2A (55% v 2%; P < .0001), CCND1 amplification (22% v 2%; P < .0001), FAT1 (25% v 4%; P < .0001), NOTCH1 (19% v 6%; P = .002), and EGFR amplification (11% v 0%; P < .0001), as well as a higher rate of 9p24.1 (PDL1/PDL2) amplification (5% v 1%) and PD-L1 immunohistochemistry high-positive tumor staining (33% v 9%; P = .04).

CONCLUSION

Comprehensive molecular profiles of vSCC vary considerably with hrHPV status and may inform patient selection into clinical trials. Sixty-one percent of HPV+ vSCCs had a pathogenic alteration in the PI3K/mTOR pathway, whereas HPV– vSCCs showed alterations in TP53, TERTp, CDKN2A, CCND1, and EGFR, and biomarkers associated with responsiveness to immunotherapy.

Somatic genetic aberrations in benign breast disease and the risk of subsequent breast cancer

It is largely unknown how the development of breast cancer (BC) is transduced by somatic genetic alterations in the benign breast. Since benign breast disease is an established risk factor for BC, we established a case-control study of women with a history of benign breast biopsy (BBB). Cases developed BC at least one year after BBB and controls did not develop BC over an average of 17 years following BBB. 135 cases were matched to 69 controls by age and type of benign change: non-proliferative or proliferation without atypia (PDWA). Whole-exome sequencing (WES) was performed for the BBB. Germline DNA (available from n = 26 participants) was utilized to develop a mutation-calling pipeline, to allow differentiation of somatic from germline variants. Among the 204 subjects, two known mutational signatures were identified, along with a currently uncatalogued signature that was significantly associated with triple negative BC (TNBC) (p = 0.007). The uncatalogued mutational signature was validated in 109 TNBCs from TCGA (p = 0.001). Compared to non-proliferative samples, PDWA harbors more abundant mutations at PIK3CA pH1047R (p < 0.001). Among the 26 BBB whose somatic copy number variation could be assessed, deletion of MLH3 is significantly associated with the mismatch repair mutational signature (p < 0.001). Matched BBB-cancer pairs were available for ten cases; several mutations were shared between BBB and cancers. This initial study of WES of BBB shows its potential for the identification of genetic alterations that portend breast oncogenesis. In future larger studies, robust personalized breast cancer risk indicators leading to novel interception paradigms can be assessed.

Quantification of ongoing APOBEC3A activity in tumor cells by monitoring RNA editing at hotspots

APOBEC3A is a cytidine deaminase driving mutagenesis, DNA replication stress and DNA damage in cancer cells. While the APOBEC3A-induced vulnerability of cancers offers an opportunity for therapy, APOBEC3A protein and mRNA are difficult to quantify in tumors due to their low abundance. Here, we describe a quantitative and sensitive assay to measure the ongoing activity of APOBEC3A in tumors. Using hotspot RNA mutations identified from APOBEC3A-positive tumors and droplet digital PCR, we develop an assay to quantify the RNA-editing activity of APOBEC3A. This assay is superior to APOBEC3A protein- and mRNA-based assays in predicting the activity of APOBEC3A on DNA. Importantly, we demonstrate that the RNA mutation-based APOBEC3A assay is applicable to clinical samples from cancer patients. Our study presents a strategy to follow the dysregulation of APOBEC3A in tumors, providing opportunities to investigate the role of APOBEC3A in tumor evolution and to target the APOBEC3A-induced vulnerability in therapy.

The DNA cytosine deaminases APOBEC3A and APOBEC3B have emerged from cancer genomics studies as drivers of mutation in cancers and tumor heterogeneity. Here the authors present a computational approach to identify the RNA mutations specifically driven by APOBEC3A, and developed an RNA mutation-based assay to quantify ongoing APOBEC3A activity in tumor cells.

Genome aging: somatic mutation in the brain links age-related decline with disease and nominates pathogenic mechanisms

Aging is a mysterious process, not only controlled genetically but also subject to random damage that can accumulate over time. While DNA damage and subsequent mutation in somatic cells were first proposed as drivers of aging more than 60 years ago, whether and to what degree these processes shape the neuronal genome in the human brain could not be tested until recent technological breakthroughs related to single-cell whole-genome sequencing. Indeed, somatic single-nucleotide variants (SNVs) increase with age in the human brain, in a somewhat stochastic process that may nonetheless be controlled by underlying genetic programs. Evidence from the literature suggests that in addition to demonstrated increases in somatic SNVs during aging in normal brains, somatic mutation may also play a role in late-onset, sporadic neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease. In this review, we will discuss somatic mutation in the human brain, mechanisms by which somatic mutations occur and can be controlled, and how this process can impact human health.

Whole-genome sequencing of 508 patients identifies key molecular features associated with poor prognosis in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a poor-prognosis cancer type with limited understanding of its molecular etiology. Using 508 ESCC genomes, we identified five novel significantly mutated genes and uncovered mutational signature clusters associated with metastasis and patients' outcomes. Several functional assays implicated that NFE2L2 may act as a tumor suppressor in ESCC and that mutations in NFE2L2 probably impaired its tumor-suppressive function, or even conferred oncogenic activities. Additionally, we found that the NFE2L2 mutations were significantly associated with worse prognosis of ESCC. We also identified potential noncoding driver mutations including hotspot mutations in the promoter region of SLC35E2 that were correlated with worse survival. Approximately 5.9% and 15.2% of patients had high tumor mutation burden or actionable mutations, respectively, and may benefit from immunotherapy or targeted therapies. We found clinically relevant coding and noncoding genomic alterations and revealed three major subtypes that robustly predicted patients' outcomes. Collectively, we report the largest dataset of genomic profiling of ESCC useful for developing ESCC-specific biomarkers for diagnosis and treatment.

Mutational signatures are jointly shaped by DNA damage and repair

Cells possess an armamentarium of DNA repair pathways to counter DNA damage and prevent mutation. Here we use C. elegans whole genome sequencing to systematically quantify the contributions of these factors to mutational signatures. We analyse 2,717 genomes from wild-type and 53 DNA repair defective backgrounds, exposed to 11 genotoxins, including UV-B and ionizing radiation, alkylating compounds, aristolochic acid, aflatoxin B1, and cisplatin. Combined genotoxic exposure and DNA repair deficiency alters mutation rates or signatures in 41% of experiments, revealing how different DNA alterations induced by the same genotoxin are mended by separate repair pathways. Error-prone translesion synthesis causes the majority of genotoxin-induced base substitutions, but averts larger deletions. Nucleotide excision repair prevents up to 99% of point mutations, almost uniformly across the mutation spectrum. Our data show that mutational signatures are joint products of DNA damage and repair and suggest that multiple factors underlie signatures observed in cancer genomes.

Experimental investigations of carcinogen-induced mutation spectra: Innovation, challenges and future directions

Recent years have witnessed an expansion of mutagenesis research focusing on experimentally modeled genome-scale mutational signatures of carcinogens and of endogenous processes. Experimental mutational signatures can explain etiologic links to patterns found in human tumors that may be linked to same exposures, and can serve as biomarkers of exposure history and may even provide insights on causality. A number of innovative exposure models have been employed and reported, based on cells cultured in monolayers or in 3-D, on organoids, induced pluripotent stem cells, non-mammalian organisms, microorganisms and rodent bioassays. Here we discuss some of the latest developments and pros and cons of these experimental systems used in mutational signature analysis. Integrative designs that bring together multiple exposure systems (in vitro, in vivo and in silico pan-cancer data mining) started emerging as powerful tools to identify robust mutational signatures of the tested cancer risk agents. We further propose that devising a new generation of cell-based models is warranted to streamline systematic testing of carcinogen effects on the cell genomes, while seeking to increasingly supplant animal with non-animal systems to address relevant ethical issues and accentuate the 3R principles. We conclude that the knowledge accumulating from the growing body of signature modelling investigations has considerable power to advance cancer etiology studies and to support cancer prevention efforts through streamlined characterization of cancer-causing agents and the recognition of their specific effects.

Low rates of mutation in clinical grade human pluripotent stem cells under different culture conditions

The occurrence of repetitive genomic changes that provide a selective growth advantage in pluripotent stem cells is of concern for their clinical application. However, the effect of different culture conditions on the underlying mutation rate is unknown. Here we show that the mutation rate in two human embryonic stem cell lines derived and banked for clinical application is low and not substantially affected by culture with Rho Kinase inhibitor, commonly used in their routine maintenance. However, the mutation rate is reduced by >50% in cells cultured under 5% oxygen, when we also found alterations in imprint methylation and reversible DNA hypomethylation. Mutations are evenly distributed across the chromosomes, except for a slight increase on the X-chromosome, and an elevation in intergenic regions suggesting that chromatin structure may affect mutation rate. Overall the results suggest that pluripotent stem cells are not subject to unusually high rates of genetic or epigenetic alterations.

Exome sequencing identifies novel mutation signatures of UV radiation and trichostatin A in primary human keratinocytes

Canonical ultraviolet (UV) mutation type and spectra are traditionally defined by direct sequencing-based approaches to map mutations in a limited number of representative DNA elements. To obtain an unbiased view of genome wide UV mutation features, we performed whole exome-sequencing (WES) to profile single nucleotide substitutions in UVB-irradiated primary human keratinocytes. Cross comparison of UV mutation profiles under different UVB radiation conditions revealed that T > C transition was highly prevalent in addition to C > T transition. We also identified 5'-ACG-3' as a common sequence motif of C > T transition. Furthermore, our analyses uncovered several recurring UV mutations following acute UVB radiation affecting multiple genes including HRNR, TRIOBP, KCNJ12, and KMT2C, which are frequently mutated in skin cancers, indicating their potential role as founding mutations in UV-induced skin tumorigenesis. Pretreatment with trichostatin A, a pan-histone deacetylase inhibitor that renders chromatin decondensation, significantly decreased the number of mutations in UVB-irradiated keratinocytes. Unexpectedly, we found trichostatin A to be a mutagen that caused DNA damage and mutagenesis at least partly through increased reactive oxidation. In summary, our study reveals new UV mutation features following acute UVB radiation and identifies novel UV mutation hotspots that may potentially represent founding driver mutations in skin cancer development.

Two main mutational processes operate in the absence of DNA mismatch repair

The analysis of tumour genome sequences has demonstrated high rates of base substitution mutagenesis upon the inactivation of DNA mismatch repair (MMR), and the resulting somatic mutations in MMR deficient tumours appear to significantly enhance the response to immune therapy. A handful of different algorithmically derived base substitution mutation signatures have been attributed to MMR deficiency in tumour somatic mutation datasets. In contrast, mutation data obtained from whole genome sequences of isogenic wild type and MMR deficient cell lines in this study, as well as from published sources, show a more uniform experimental mutation spectrum of MMR deficiency. In order to resolve this discrepancy, we reanalysed mutation data from MMR deficient tumour whole exome and whole genome sequences. We derived two base substitution signatures using non-negative matrix factorisation, which together adequately describe mutagenesis in all tumour and cell line samples. The two new signatures broadly resemble COSMIC signatures 6 and 20, but perform better than existing COSMIC signatures at identifying MMR deficient tumours in mutation signature deconstruction. We show that the contribution of the two identified signatures, one of which is dominated by C to T mutations at CpG sites, is biased by the different sequence composition of the exome and the whole genome. We further show that the identity of the inactivated MMR gene, the tissue type, the mutational burden or the patient's age does not influence the mutation spectrum, but that a tendency for a greater contribution by the CpG mutational process is observed in tumours as compared to cultured cells. Our analysis suggest that two separable mutational processes operate in the genomes of MMR deficient cells.

Mutational signatures: experimental design and analytical framework

Mutational signatures provide a powerful alternative for understanding the pathophysiology of cancer. Currently, experimental efforts aimed at validating and understanding the etiologies of cancer-derived mutational signatures are underway. In this review, we highlight key aspects of mutational signature experimental design and describe the analytical framework. We suggest guidelines and quality control measures for handling whole-genome sequencing data for mutational signature analyses and discuss pitfalls in interpretation. We envision that improved next-generation sequencing technologies and molecular cell biology approaches will usher in the next generation of studies into the etiologies and mechanisms of mutational patterns uncovered in cancers.

Tobacco smoking and somatic mutations in human bronchial epithelium

Tobacco smoking causes lung cancer1-3, a process that is driven by more than 60 carcinogens in cigarette smoke that directly damage and mutate DNA4,5. The profound effects of tobacco on the genome of lung cancer cells are well-documented6-10, but equivalent data for normal bronchial cells are lacking. Here we sequenced whole genomes of 632 colonies derived from single bronchial epithelial cells across 16 subjects. Tobacco smoking was the major influence on mutational burden, typically adding from 1,000 to 10,000 mutations per cell; massively increasing the variance both within and between subjects; and generating several distinct mutational signatures of substitutions and of insertions and deletions. A population of cells in individuals with a history of smoking had mutational burdens that were equivalent to those expected for people who had never smoked: these cells had less damage from tobacco-specific mutational processes, were fourfold more frequent in ex-smokers than current smokers and had considerably longer telomeres than their more-mutated counterparts. Driver mutations increased in frequency with age, affecting 4-14% of cells in middle-aged subjects who had never smoked. In current smokers, at least 25% of cells carried driver mutations and 0-6% of cells had two or even three drivers. Thus, tobacco smoking increases mutational burden, cell-to-cell heterogeneity and driver mutations, but quitting promotes replenishment of the bronchial epithelium from mitotically quiescent cells that have avoided tobacco mutagenesis.

Low-abundance mutations in colorectal cancer patients and healthy adults

Detecting low-abundance mutations is very important for cancer diagnosis and treatment. Here we describe an improved targeted sequencing analysis that dramatically increases sequencing depth. Seven colorectal cancer (CRC) patients and seven healthy adults were enrolled in this study. We examined genetic mutations in tissue samples from the central and peripheral regions of tumors from the CRC patients and in blood cells from the healthy adults. We observed that each CRC carried larger numbers of mutations more than previously estimated. These included numerous deletion mutations in the tumor tissue. While the cellular morphology in the surrounding normal colonic tissues was healthy, these cells also carried many mutations. Similarly, the blood cells from the healthy donors carried numerous mutations. These findings shed new light on the processes of tumorigenesis and aging, and also present a potentially effective method for detecting low-abundance mutations for cancer diagnosis and targeted treatments.

Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

The genomics of oxidative DNA damage, repair, and resulting mutagenesis

Reactive oxygen species are a constant threat to DNA as they modify bases with the risk of disrupting genome function, inducing genome instability and mutation. Such risks are due to primary oxidative DNA damage and also mediated by the repair process. This leads to a delicate decision process for the cell as to whether to repair a damaged base at a specific genomic location or better leave it unrepaired. Persistent DNA damage can disrupt genome function, but on the other hand it can also contribute to gene regulation by serving as an epigenetic mark. When such processes are out of balance, pathophysiological conditions could get accelerated, because oxidative DNA damage and resulting mutagenic processes are tightly linked to ageing, inflammation, and the development of multiple age-related diseases, such as cancer and neurodegenerative disorders. Recent technological advancements and novel data analysis strategies have revealed that oxidative DNA damage, its repair, and related mutations distribute heterogeneously over the genome at multiple levels of resolution. The involved mechanisms act in the context of genome sequence, in interaction with genome function and chromatin. This review addresses what we currently know about the genome distribution of oxidative DNA damage, repair intermediates, and mutations. It will specifically focus on the various methodologies to measure oxidative DNA damage distribution and discuss the mechanistic conclusions derived from the different approaches. It will also address the consequences of oxidative DNA damage, specifically how it gives rise to mutations, genome instability, and how it can act as an epigenetic mark.

Perturbation of base excision repair sensitizes breast cancer cells to APOBEC3 deaminase-mediated mutations

Abundant APOBEC3 (A3) deaminase-mediated mutations can dominate the mutational landscape ('mutator phenotype') of some cancers, however, the basis of this sporadic vulnerability is unknown. We show here that elevated expression of the bifunctional DNA glycosylase, NEIL2, sensitizes breast cancer cells to A3B-mediated mutations and double-strand breaks (DSBs) by perturbing canonical base excision repair (BER). NEIL2 usurps the canonical lyase, APE1, at abasic sites in a purified BER system, rendering them poor substrates for polymerase β. However, the nicked NEIL2 product can serve as an entry site for Exo1 in vitro to generate single-stranded DNA, which would be susceptible to both A3B and DSBs. As NEIL2 or Exo1 depletion mitigates the DNA damage caused by A3B expression, we suggest that aberrant NEIL2 expression can explain certain instances of A3B-mediated mutations.

Targeted drugs ramp up cancer mutability

Drug-induced mutagenesis may accelerate cancer evolution and resistance

Mutation signature analysis identifies increased mutation caused by tobacco smoke associated DNA adducts in larynx squamous cell carcinoma compared with oral cavity and oropharynx

Squamous cell carcinomas of the head and neck (HNSCC) arise from mucosal keratinocytes of the upper aero-digestive tract. Despite a common cell of origin and similar driver-gene mutations which divert cell fate from differentiation to proliferation, HNSCC are considered a heterogeneous group of tumors categorized by site of origin within the aero-digestive mucosa, and the presence or absence of HPV infection. Tobacco use is a major driver of carcinogenesis in HNSCC and is a poor prognosticator that has previously been associated with poor immune cell infiltration and higher mutation numbers. Here, we study patterns of mutations in HNSCC that are derived from the specific nucleotide changes and their surrounding nucleotide context (also known as mutation signatures). We identify that mutations linked to DNA adducts associated with tobacco smoke exposure are predominantly found in the larynx. Presence of this class of mutation, termed COSMIC signature 4, is responsible for the increased burden of mutation in this anatomical sub-site. In addition, we show that another mutation pattern, COSMIC signature 5, is positively associated with age in HNSCC from non-smokers and that larynx SCC from non-smokers have a greater number of signature 5 mutations compared with other HNSCC sub-sites. Immunohistochemistry demonstrates a significantly lower Ki-67 proliferation index in size matched larynx SCC compared with oral cavity SCC and oropharynx SCC. Collectively, these observations support a model where larynx SCC are characterized by slower growth and increased susceptibility to mutations from tobacco carcinogen DNA adducts.

FIREVAT: finding reliable variants without artifacts in human cancer samples using etiologically relevant mutational signatures

Background

Accurate identification of real somatic variants is a primary part of cancer genome studies and precision oncology. However, artifacts introduced in various steps of sequencing obfuscate confidence in variant calling. Current computational approaches to variant filtering involve intensive interrogation of Binary Alignment Map (BAM) files and require massive computing power, data storage, and manual labor. Recently, mutational signatures associated with sequencing artifacts have been extracted by the Pan-cancer Analysis of Whole Genomes (PCAWG) study. These spectrums can be used to evaluate refinement quality of a given set of somatic mutations.

Results

Here we introduce a novel variant refinement software, FIREVAT (FInding REliable Variants without ArTifacts), which uses known spectrums of sequencing artifacts extracted from one of the largest publicly available catalogs of human tumor samples. FIREVAT performs a quick and efficient variant refinement that accurately removes artifacts and greatly improves the precision and specificity of somatic calls. We validated FIREVAT refinement performance using orthogonal sequencing datasets totaling 384 tumor samples with respect to ground truth. Our novel method achieved the highest level of performance compared to existing filtering approaches. Application of FIREVAT on additional 308 The Cancer Genome Atlas (TCGA) samples demonstrated that FIREVAT refinement leads to identification of more biologically and clinically relevant mutational signatures as well as enrichment of sequence contexts associated with experimental errors. FIREVAT only requires a Variant Call Format file (VCF) and generates a comprehensive report of the variant refinement processes and outcomes for the user.

Conclusions

In summary, FIREVAT facilitates a novel refinement strategy using mutational signatures to distinguish artifactual point mutations called in human cancer samples. We anticipate that FIREVAT results will further contribute to precision oncology efforts that rely on accurate identification of variants, especially in the context of analyzing mutational signatures that bear prognostic and therapeutic significance. FIREVAT is freely available at https://github.com/cgab-ncc/FIREVAT

Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans

Reducing or eliminating persistent disparities in lung cancer incidence and survival has been challenging because our current understanding of lung cancer biology is derived primarily from populations of European descent. Here we show results from a targeted sequencing panel using NCI-MD Case Control Study patient samples and reveal a significantly higher prevalence of PTPRT and JAK2 mutations in lung adenocarcinomas among African Americans compared with European Americans. This increase in mutation frequency was validated with independent WES data from the NCI-MD Case Control Study and TCGA. We find that patients carrying these mutations have a concomitant increase in IL-6/STAT3 signaling and miR-21 expression. Together, these findings suggest the identification of these potentially actionable mutations could have clinical significance for targeted therapy and the enrollment of minority populations in clinical trials.

A Subset of Colorectal Cancers with Cross-Sensitivity to Olaparib and Oxaliplatin

PURPOSE

Defects in the homologous recombination (HR) repair pathway are of clinical interest due to sensitivity of HR-deficient cells to PARP inhibitors. We were interested in defining PARP vulnerability in patients with metastatic colorectal cancer (mCRC) carrying KRAS and BRAF mutations who display poor prognosis, have limited therapeutic options, and represent an unmet clinical need.

EXPERIMENTAL DESIGN

We tested colorectal cancer cell lines, patient-derived organoids (PDO), and patient-derived xenografts (PDX) enriched for KRAS and BRAF mutations for sensitivity to the PARP inhibitor olaparib, and the chemotherapeutic agents oxaliplatin and 5-fluorouracil (5-FU). Genomic profiles and DNA repair proficiency of colorectal cancer models were compared with pharmacologic response.

RESULTS

Thirteen of 99 (around 13%) colorectal cancer cell lines were highly sensitive to clinically active concentrations of olaparib and displayed functional deficiency in HR. Response to PARP blockade was positively correlated with sensitivity to oxaliplatin in colorectal cancer cell lines as well as patient-derived organoids. Treatment of PDXs with olaparib impaired tumor growth and maintenance therapy with PARP blockade after initial oxaliplatin response delayed disease progression in mice.

CONCLUSIONS

These results indicate that a colorectal cancer subset characterized by poor prognosis and limited therapeutic options is vulnerable to PARP inhibition and suggest that PDO-based drug-screening assays can be used to identify patients with colorectal cancer likely to benefit from olaparib. As patients with mCRC almost invariably receive therapies based on oxaliplatin, "maintenance" treatment with PARP inhibitors warrants further clinical investigation.

APOBEC3A is a prominent cytidine deaminase in breast cancer

APOBEC cytidine deaminases are the second-most prominent source of mutagenesis in sequenced tumors. Previous studies have proposed that APOBEC3B (A3B) is the major source of mutagenesis in breast cancer (BRCA). We show that APOBEC3A (A3A) is the only APOBEC whose expression correlates with APOBEC-induced mutation load and that A3A expression is responsible for cytidine deamination in multiple BRCA cell lines. Comparative analysis of A3A and A3B expression by qRT-PCR, RSEM-normalized RNA-seq, and unambiguous RNA-seq validated the use of RNA-seq to measure APOBEC expression, which indicates that A3A is the primary correlate with APOBEC-mutation load in primary BRCA tumors. We also demonstrate that A3A has >100-fold more cytidine deamination activity than A3B in the presence of cellular RNA, likely explaining why higher levels of A3B expression contributes less to mutagenesis in BRCA. Our findings identify A3A as a major source of cytidine deaminase activity in breast cancer cells and possibly a prominent contributor to the APOBEC mutation signature.

The case for BK polyomavirus as a cause of bladder cancer

In 2014, the International Agency for Research on Cancer judged Merkel cell polyomavirus (MCPyV) to be a probable human carcinogen. BK polyomavirus (BKPyV, a distant cousin of MCPyV) was ruled a possible carcinogen. In this review, we argue that it has recently become reasonable to view both of these viruses as known human carcinogens. In particular, several complementary lines of evidence support a causal role for BKPyV in the development of bladder carcinomas affecting organ transplant patients. The expansion of inexpensive deep sequencing has opened new approaches to investigating the important question of whether BKPyV causes urinary tract cancers in the general population.

Epidermolysis Bullosa-Associated Squamous Cell Carcinoma: From Pathogenesis to Therapeutic Perspectives

Epidermolysis bullosa (EB) is a heterogeneous group of inherited skin disorders determined by mutations in genes encoding for structural components of the cutaneous basement membrane zone. Disease hallmarks are skin fragility and unremitting blistering. The most disabling EB (sub)types show defective wound healing, fibrosis and inflammation at lesional skin. These features expose patients to serious disease complications, including the development of cutaneous squamous cell carcinomas (SCCs). Almost all subjects affected with the severe recessive dystrophic EB (RDEB) subtype suffer from early and extremely aggressive SCCs (RDEB-SCC), which represent the first cause of death in these patients. The genetic determinants of RDEB-SCC do not exhaustively explain its unique behavior as compared to low-risk, ultraviolet-induced SCCs in the general population. On the other hand, a growing body of evidence points to the key role of tumor microenvironment in initiation, progression and spreading of RDEB-SCC, as well as of other, less-investigated, EB-related SCCs (EB-SCCs). Here, we discuss the recent advances in understanding the complex series of molecular events (i.e., fibrotic, inflammatory, and immune processes) contributing to SCC development in EB patients, cross-compare tumor features in the different EB subtypes and report the most promising therapeutic approaches to counteract or delay EB-SCCs.

Mouse APOBEC1 cytidine deaminase can induce somatic mutations in chromosomal DNA

Background

APOBEC1 (A1) enzymes are cytidine deaminases involved in RNA editing. In addition to this activity, a few A1 enzymes have been shown to be active on single stranded DNA. As two human ssDNA cytidine deaminases APOBEC3A (A3A), APOBEC3B (A3B) and related enzymes across the spectrum of placental mammals have been shown to introduce somatic mutations into nuclear DNA of cancer genomes, we explored the mutagenic threat of A1 cytidine deaminases to chromosomal DNA.

Results

Molecular cloning and expression of various A1 enzymes reveal that the cow, pig, dog, rabbit and mouse A1 have an intracellular ssDNA substrate specificity. However, among all the enzymes studied, mouse A1 appears to be singular, being able to introduce somatic mutations into nuclear DNA with a clear 5’TpC editing context, and to deaminate 5-methylcytidine substituted DNA which are characteristic features of the cancer related mammalian A3A and A3B enzymes. However, mouse A1 activity fails to elicit formation of double stranded DNA breaks, suggesting that mouse A1 possess an attenuated nuclear DNA mutator phenotype reminiscent of human A3B.

Conclusions

At an experimental level mouse APOBEC1 is remarkable among 12 mammalian A1 enzymes in that it represents a source of somatic mutations in mouse genome, potentially fueling oncogenesis. While the order Rodentia is bereft of A3A and A3B like enzymes it seems that APOBEC1 may well substitute for it, albeit remaining much less active. This modifies the paradigm that APOBEC3 and AID enzymes are the sole endogenous mutator enzymes giving rise to off-target editing of mammalian genomes.

Adaptive mutability of colorectal cancers in response to targeted therapies

The emergence of drug resistance limits the efficacy of targeted therapies in human tumors. The prevalent view is that resistance is a fait accompli: when treatment is initiated, cancers already contain drug-resistant mutant cells. Bacteria exposed to antibiotics transiently increase their mutation rates (adaptive mutability), thus improving the likelihood of survival. We investigated whether human colorectal cancer (CRC) cells likewise exploit adaptive mutability to evade therapeutic pressure. We found that epidermal growth factor receptor (EGFR)/BRAF inhibition down-regulates mismatch repair (MMR) and homologous recombination DNA-repair genes and concomitantly up-regulates error-prone polymerases in drug-tolerant (persister) cells. MMR proteins were also down-regulated in patient-derived xenografts and tumor specimens during therapy. EGFR/BRAF inhibition induced DNA damage, increased mutability, and triggered microsatellite instability. Thus, like unicellular organisms, tumor cells evade therapeutic pressures by enhancing mutability.

The spectrum of APOBEC3 activity: From anti-viral agents to anti-cancer opportunities

The APOBEC3 family of cytosine deaminases are part of the innate immune response to viral infection, but also have the capacity to damage cellular DNA. Detection of mutational signatures consistent with APOBEC3 activity, together with elevated APOBEC3 expression in cancer cells, has raised the possibility that these enzymes contribute to oncogenesis. Genome deamination by APOBEC3 enzymes also elicits DNA damage response signaling and presents therapeutic vulnerabilities for cancer cells. Here, we discuss implications of APOBEC3 activity in cancer and the potential to exploit their mutagenic activity for targeted cancer therapies.

Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation

A single cancer genome can harbor thousands of clustered mutations. Mutation signature analyses have revealed that the origin of clusters are lesions in long tracts of single-stranded (ss) DNA damaged by apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases, raising questions about molecular mechanisms that generate long ssDNA vulnerable to hypermutation. Here, we show that ssDNA intermediates formed during the repair of gamma-induced bursts of double-strand breaks (DSBs) in the presence of APOBEC3A in yeast lead to multiple APOBEC-induced clusters similar to cancer. We identified three independent pathways enabling cluster formation associated with repairing bursts of DSBs: 5′ to 3′ bidirectional resection, unidirectional resection, and break-induced replication (BIR). Analysis of millions of mutations in APOBEC-hypermutated cancer genomes revealed that cancer tolerance to formation of hypermutable ssDNA is similar to yeast and that the predominant pattern of clustered mutagenesis is the same as in resection-defective yeast, suggesting that cluster formation in cancers is driven by a BIR-like mechanism. The phenomenon of genome-wide burst of clustered mutagenesis revealed by our study can play an important role in generating somatic hypermutation in cancers as well as in noncancerous cells.

This study uses yeast expressing a human cytidine deaminase to reveal simultaneous stretches of long single-strand DNA and multiple vast mutation clusters in a single eukaryotic cell repairing multiple double-strand breaks. This is reminiscent of the phenomenon of “kataegis” or hypermutation observed in cancer genomes, suggesting that a similar mechanism is involved.

Drug-adapted cancer cell lines as preclinical models of acquired resistance

Acquired resistance formation limits the efficacy of anti-cancer therapies. Acquired and intrinsic resistance differ conceptually. Acquired resistance is the consequence of directed evolution, whereas intrinsic resistance depends on the (stochastic) presence of pre-existing resistance mechanisms. Preclinical model systems are needed to study acquired drug resistance because they enable: (1) in depth functional studies; (2) the investigation of non-standard treatments for a certain disease condition (which is necessary to identify small groups of responders); and (3) the comparison of multiple therapies in the same system. Hence, they complement data derived from clinical trials and clinical specimens, including liquid biopsies. Many groups have successfully used drug-adapted cancer cell lines to identify and elucidate clinically relevant resistance mechanisms to targeted and cytotoxic anti-cancer drugs. Hence, we argue that drug-adapted cancer cell lines represent a preclinical model system in their own right that is complementary to other preclinical model systems and clinical data.

PCR-free whole exome sequencing: Cost-effective and efficient in detecting rare mutations

In this study, we describe the development of a PCR-free whole exome sequencing method. Using this method, 2 μg DNA was sufficient for library preparation for whole exome sequencing. Furthermore, the method is simple and makes use of a commercial kit, with additional step of concentrating the captured library by ethanol precipitation. The accuracy of the PCR-free method was found to be equivalent to that of unique molecular identifier-corrected analysis method, which is the commonly used method to detect rare mutations. Thus, the PCR-free whole exome sequencing method is cost-effective as well as efficient in detecting rare mutations.