Mechanisms of base substitution mutagenesis in cancer genomes

Histone mutations in cancer

Genes encoding histone proteins are recurrently mutated in tumor samples, and these mutations may impact nucleosome stability, histone post-translational modification, or chromatin dynamics. The prevalence of histone mutations across diverse cancer types suggest that normal chromatin structure is a barrier to tumorigenesis. Oncohistone mutations disrupt chromatin structure and gene regulatory mechanisms, resulting in aberrant gene expression and the development of cancer phenotypes. Examples of oncohistones include the histone H3 K27M mutation found in pediatric brain cancers that blocks post-translational modification of the H3 N-terminal tail and the histone H2B E76K mutation found in some solid tumors that disrupts nucleosome stability. Oncohistones may comprise a limited fraction of the total histone pool yet cause global effects on chromatin structure and drive cancer phenotypes. Here, we survey histone mutations in cancer and review their function and role in tumorigenesis.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Estimating the Vertical Ionization Potential of Single-Stranded DNA Molecules

The electronic properties of DNA molecules, defined by the sequence-dependent ionization potentials of nucleobases, enable long-range charge transport along the DNA stacks. This has been linked to a range of key physiological processes in the cells and to the triggering of nucleobase substitutions, some of which may cause diseases. To gain molecular-level understanding of the sequence dependence of these phenomena, we estimated the vertical ionization potential (vIP) of all possible nucleobase stacks in B-conformation, containing one to four Gua, Ade, Thy, Cyt, or methylated Cyt. To do this, we used quantum chemistry calculations and more precisely the second-order Møller-Plesset perturbation theory (MP2) and three double-hybrid density functional theory methods, combined with several basis sets for describing atomic orbitals. The calculated vIP of single nucleobases were compared to experimental data and those of nucleobase pairs, triplets, and quadruplets, to observed mutability frequencies in the human genome, reported to be correlated with vIP values. This comparison selected MP2 with the 6-31G* basis set as the best of the tested calculation levels. These results were exploited to set up a recursive model, called vIPer, which estimates the vIP of all possible single-stranded DNA sequences of any length based on the calculated vIPs of overlapping quadruplets. vIPer's vIP values correlate well with oxidation potentials measured by cyclic voltammetry and activities obtained through photoinduced DNA cleavage experiments, further validating our approach. vIPer is freely available on the github.com/3BioCompBio/vIPer repository.

Genotype-Phenotype Correlations in Autosomal Dominant and Recessive APC Mutation-Negative Colorectal Adenomatous Polyposis

The most prevalent type of intestinal polyposis, colorectal adenomatous polyposis (CAP), is regarded as a precancerous lesion of colorectal cancer with obvious genetic characteristics. Early screening and intervention can significantly improve patients' survival and prognosis. The adenomatous polyposis coli (APC) mutation is believed to be the primary cause of CAP. There is, however, a subset of CAP with undetectable pathogenic mutations in APC, known as APC (-)/CAP. The genetic predisposition to APC (-)/CAP has largely been associated with germline mutations in some susceptible genes, including the human mutY homologue (MUTYH) gene and the Nth-like DNA glycosylase 1 (NTHL1) gene, and DNA mismatch repair (MMR) can cause autosomal recessive APC (-)/CAP. Furthermore, autosomal dominant APC (-)/CAP could occur as a result of DNA polymerase epsilon (POLE)/DNA polymerase delta 1 (POLD1), axis inhibition protein 2 (AXIN2), and dual oxidase 2 (DUOX2) mutations. The clinical phenotypes of these pathogenic mutations vary greatly depending on their genetic characteristics. Therefore, in this study, we present a comprehensive review of the association between autosomal recessive and dominant APC (-)/CAP genotypes and clinical phenotypes and conclude that APC (-)/CAP is a disease caused by multiple genes with different phenotypes and interaction exists in the pathogenic genes.

Multiple roles of apolipoprotein B mRNA editing enzyme catalytic subunit 3B (APOBEC3B) in human tumors: a pan-cancer analysis

Although there have been some recent cell and animal experiments indicating that expression of the gene encoding apolipoprotein B mRNA editing enzyme catalytic subunit 3B (APOBEC3B) is closely related to cancer, it still lacks pan-cancer analysis. Here we analyzed the potential carcinogenic role of APOBEC3B in 33 tumors based on The Cancer Genome Atlas (TCGA). APOBEC3B was highly expressed in most tumors and weakly expressed in a few. Differences in expression level were significantly correlated with the pathological tumor stage and prognosis of affected patients. The high-frequency APOBEC3B changes were principally mutations and amplifications in some tumors, such as uterine corpus endometrial carcinomas or cutaneous melanomas. In testicular germ cell tumors and invasive breast carcinomas, APOBEC3B expression and CD8⁺ T lymphocyte counts were correlated. In other cancers, such as human papilloma virus (HPV)-related head and neck squamous cell carcinomas or esophageal adenocarcinomas, there was also cancer-associated fibroblast infiltration. The APOBEC3B enzyme acts in the mitochondrial respiratory electron transport chain and in oxidative phosphorylation. This first pan-cancer study provides a comprehensive understanding of the multiple roles of APOBEC3B in different tumor types.

Structures of LIG1 that engage with mutagenic mismatches inserted by polβ in base excision repair

DNA ligase I (LIG1) catalyzes the ligation of the nick repair intermediate after gap filling by DNA polymerase (pol) β during downstream steps of the base excision repair (BER) pathway. However, how LIG1 discriminates against the mutagenic 3'-mismatches incorporated by polβ at atomic resolution remains undefined. Here, we determine the X-ray structures of LIG1/nick DNA complexes with G:T and A:C mismatches and uncover the ligase strategies that favor or deter the ligation of base substitution errors. Our structures reveal that the LIG1 active site can accommodate a G:T mismatch in the wobble conformation, where an adenylate (AMP) is transferred to the 5'-phosphate of a nick (DNA-AMP), while it stays in the LIG1-AMP intermediate during the initial step of the ligation reaction in the presence of an A:C mismatch at the 3'-strand. Moreover, we show mutagenic ligation and aberrant nick sealing of dG:T and dA:C mismatches, respectively. Finally, we demonstrate that AP-endonuclease 1 (APE1), as a compensatory proofreading enzyme, removes the mismatched bases and interacts with LIG1 at the final BER steps. Our overall findings provide the features of accurate versus mutagenic outcomes coordinated by a multiprotein complex including polβ, LIG1, and APE1 to maintain efficient repair.

SomaticSiMu: A mutational signature simulator

SUMMARY

SomaticSiMu is an in silico simulator of single and double base substitutions, and single base insertions and deletions in an input genomic sequence to mimic mutational signatures. SomaticSiMu outputs simulated DNA sequences and mutational catalogues with imposed mutational signatures. The tool is the first mutational signature simulator featuring a graphical user interface, control of mutation rates, and built-in visualization tools of the simulated mutations. Simulated datasets are useful as a ground truth to test the accuracy and sensitivity of DNA sequence classification tools and mutational signature extraction tools under different experimental scenarios. The reliability of SomaticSiMu was affirmed by 1) supervised machine learning classification of simulated sequences with different mutation types and burdens, and 2) mutational signature extraction from simulated mutational catalogues.

AVAILABILITY AND IMPLEMENTATION

SomaticSiMu is written in Python 3.8.3. The open-source code, documentation, and tutorials are available at https://github.com/HillLab/SomaticSiMu under the terms of the Creative Commons Attribution 4.0 International License.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Rare Tautomers of Artificially Expanded Genetic Letters and their Effects on the Base pair Stabilities

To expand the existing genetic letters, it is necessary to design robust nucleotides that can function naturally in living cells. Therefore, it is desirable to examine the roles of recently proposed second-generation artificially expanded genetic letters in producing stable duplex DNA. Here, a reliable dispersion-corrected density functional theory method is used to understand the electronic structures and properties of different rare tautomers of proposed expanded genetic letters and their effects on the base pair stabilities in the duplex DNA. It is found that the rare tautomers are not only stable in the aqueous medium but can also base pair with natural bases to produce stable mispairs. Except for J and V, all the artificial genetic letters are found to produce mispairs that are about 1-7 kcal/mol more stable than their complementary counterparts. They are also appreciably more stable than the naturally occurring G:C, A:T, and G:T pairs. The higher base pair stabilities are found to be mainly because of the polarity of monomers and attractive electrostatic interactions.

DNA Repair Syndromes and Cancer: Insights Into Genetics and Phenotype Patterns

DNA damage response is essential to human physiology. A broad spectrum of pathologies are displayed by individuals carrying monoallelic or biallelic loss-of-function mutations in DNA damage repair genes. DNA repair syndromes with biallelic disturbance of essential DNA damage response pathways manifest early in life with multi-systemic involvement and a high propensity for hematologic and solid cancers, as well as bone marrow failure. In this review, we describe classic biallelic DNA repair cancer syndromes arising from faulty single- and double-strand DNA break repair, as well as dysfunctional DNA helicases. These clinical entities include xeroderma pigmentosum, constitutional mismatch repair deficiency, ataxia telangiectasia, Nijmegen breakage syndrome, deficiencies of DNA ligase IV, NHEJ/Cernunnos, and ERCC6L2, as well as Bloom, Werner, and Rothmund-Thompson syndromes. To give an in-depth understanding of these disorders, we provide historical overview and discuss the interplay between complex biology and heterogeneous clinical manifestations.

Relation between DNA ionization potentials, single base substitutions and pathogenic variants

Background

It is nowadays clear that single base substitutions that occur in the human genome, of which some lead to pathogenic conditions, are non-random and influenced by their flanking nucleobase sequences. However, despite recent progress, the understanding of these "non-local" effects is still far from being achieved.

Results

To advance this problem, we analyzed the relationship between the base mutability in specific gene regions and the electron hole transport along the DNA base stacks, as it is one of the mechanisms that have been suggested to contribute to these effects. More precisely, we studied the connection between the normalized frequency of single base substitutions and the vertical ionization potential of the base and its flanking sequence, estimated using MP2/6-31G* ab initio quantum chemistry calculations. We found a statistically significant overall anticorrelation between these two quantities: the lower the vIP value, the more probable the substitution. Moreover, the slope of the regression lines varies. It is larger for introns than for exons and untranslated regions, and for synonymous than for missense substitutions. Interestingly, the correlation appears to be more pronounced when considering the flanking sequence of the substituted base in the 3’ rather than in the 5’ direction, which corresponds to the preferred direction of charge migration. A weaker but still statistically significant correlation is found between the ionization potentials and the pathogenicity of the base substitutions. Moreover, pathogenicity is also preferentially associated with larger changes in ionization potentials upon base substitution.

Conclusions

With this analysis we gained new insights into the complex biophysical mechanisms that are at the basis of mutagenesis and pathogenicity, and supported the role of electron-hole transport in these matters.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5867-y) contains supplementary material, which is available to authorized users.

Mutational landscape in genetically engineered, carcinogen-induced, and radiation-induced mouse sarcoma

Cancer development is influenced by hereditary mutations, somatic mutations due to random errors in DNA replication, or external factors. It remains unclear how distinct cell-intrinsic and -extrinsic factors impact oncogenesis within the same tissue type. We investigated murine soft tissue sarcomas generated by oncogenic alterations (KrasG12D activation and p53 deletion), carcinogens (3-methylcholanthrene [MCA] or ionizing radiation), and in a novel model combining both factors (MCA plus p53 deletion). Whole-exome sequencing demonstrated distinct mutational signatures in individual sarcoma cohorts. MCA-induced sarcomas exhibited high mutational burden and predominantly G-to-T transversions, while radiation-induced sarcomas exhibited low mutational burden and a distinct genetic signature characterized by C-to-T transitions. The indel to substitution ratio and amount of gene copy number variations were high for radiation-induced sarcomas. MCA-induced tumors generated on a p53-deficient background showed the highest genomic instability. MCA-induced sarcomas harbored mutations in putative cancer-driver genes that regulate MAPK signaling (Kras and Nf1) and the Hippo pathway (Fat1 and Fat4). In contrast, radiation-induced sarcomas and KrasG12Dp53-/- sarcomas did not harbor recurrent oncogenic mutations, rather they exhibited amplifications of specific oncogenes: Kras and Myc in KrasG12Dp53-/- sarcomas, and Met and Yap1 for radiation-induced sarcomas. These results reveal that different initiating events drive oncogenesis through distinct mechanisms.

Early divergence of mutational processes in human fetal tissues

A developing human fetus needs to balance rapid cellular expansion with maintaining genomic stability. Here, we accurately quantified and characterized somatic mutation accumulation in fetal tissues by analyzing individual stem cells from human fetal liver and intestine. Fetal mutation rates were about fivefold higher than in tissue-matched adult stem cells. The mutational landscape of fetal intestinal stem cells resembled that of adult intestinal stem cells, while the mutation spectrum of fetal liver stem cells is distinct from stem cells of the fetal intestine and the adult liver. Our analyses indicate that variation in mutational mechanisms, including oxidative stress and spontaneous deamination of methylated cytosines, contributes to the observed divergence in mutation accumulation patterns and drives genetic mosaicism in humans.

Influence of Electron-Holes on DNA Sequence-Specific Mutation Rates

Biases in mutation rate can influence molecular evolution, yielding rates of evolution that vary widely in different parts of the genome and even among neighboring nucleotides. Here, we explore one possible mechanism of influence on sequence-specific mutation rates, the electron–hole, which can localize and potentially trigger a replication mismatch. A hole is a mobile site of positive charge created during one-electron oxidation by, for example, radiation, contact with a mutagenic agent, or oxidative stress. Its quantum wavelike properties cause it to localize at various sites with probabilities that vary widely, by orders of magnitude, and depend strongly on the local sequence. We find significant correlations between hole probabilities and mutation rates within base triplets, observed in published mutation accumulation experiments on four species of bacteria. We have also computed hole probability spectra for hypervariable segment I of the human mtDNA control region, which contains several mutational hotspots, and for heptanucleotides in noncoding regions of the human genome, whose polymorphism levels have recently been reported. We observe significant correlations between hole probabilities, and context-specific mutation and substitution rates. The correlation with hole probability cannot be explained entirely by CpG methylation in the heptanucleotide data. Peaks in hole probability tend to coincide with mutational hotspots, even in mtDNA where CpG methylation is rare. Our results suggest that hole-enhanced mutational mechanisms, such as oxidation-stabilized tautomerization and base deamination, contribute to molecular evolution.

NF-kappaB Regulates Redox Status in Breast Cancer Subtypes

Oxidative stress (OS) is an indispensable condition to ensure genomic instability in cancer cells. In breast cancer (BC), redox alterations have been widely characterized, but since this process results from a chain of inflammatory events, the causal molecular triggers remain to be identified. In this context, we used a microarray approach to investigate the role of the main pro-oxidant transcription factor, nuclear factor-kappa B (NF-κB), in gene profiles of BC subtypes. Our results showed that NF-κB knockdown in distinct BC subtypes led to differential expression of relevant factors involved in glutathione metabolism, prostaglandins, cytochrome P450 and cyclooxygenase, suggesting a relationship between the redox balance and NF-κB in such cells. In addition, we performed biochemical analyses to validate the microarray dataset focusing on OS and correlated these parameters with normal expression or NF-κB inhibition. Our data showed a distinct oxidative status pattern for each of the three studied BC subtype models, consistent with the intrinsic characteristics of each BC subtype. Thus, our findings suggest that NF-κB may represent an additional mechanism related to OS maintenance in BC, operating in various forms to mediate other important predominant signaling components of each BC subtype.

A site specific model and analysis of the neutral somatic mutation rate in whole-genome cancer data

Background

Detailed modelling of the neutral mutational process in cancer cells is crucial for identifying driver mutations and understanding the mutational mechanisms that act during cancer development. The neutral mutational process is very complex: whole-genome analyses have revealed that the mutation rate differs between cancer types, between patients and along the genome depending on the genetic and epigenetic context. Therefore, methods that predict the number of different types of mutations in regions or specific genomic elements must consider local genomic explanatory variables. A major drawback of most methods is the need to average the explanatory variables across the entire region or genomic element. This procedure is particularly problematic if the explanatory variable varies dramatically in the element under consideration.

Results

To take into account the fine scale of the explanatory variables, we model the probabilities of different types of mutations for each position in the genome by multinomial logistic regression. We analyse 505 cancer genomes from 14 different cancer types and compare the performance in predicting mutation rate for both regional based models and site-specific models. We show that for 1000 randomly selected genomic positions, the site-specific model predicts the mutation rate much better than regional based models.

We use a forward selection procedure to identify the most important explanatory variables. The procedure identifies site-specific conservation (phyloP), replication timing, and expression level as the best predictors for the mutation rate. Finally, our model confirms and quantifies certain well-known mutational signatures.

Conclusion

We find that our site-specific multinomial regression model outperforms the regional based models. The possibility of including genomic variables on different scales and patient specific variables makes it a versatile framework for studying different mutational mechanisms. Our model can serve as the neutral null model for the mutational process; regions that deviate from the null model are candidates for elements that drive cancer development.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2141-2) contains supplementary material, which is available to authorized users.

Different mutational rates and mechanisms in human cells at pregastrulation and neurogenesis

Somatic mosaicism in the human brain may alter function of individual neurons. We analyzed genomes of single cells from the forebrains of three human fetuses (15 to 21 weeks postconception) using clonal cell populations. We detected 200 to 400 single-nucleotide variations (SNVs) per cell. SNV patterns resembled those found in cancer cell genomes, indicating a role of background mutagenesis in cancer. SNVs with a frequency of >2% in brain were also present in the spleen, revealing a pregastrulation origin. We reconstructed cell lineages for the first five postzygotic cleavages and calculated a mutation rate of ~1.3 mutations per division per cell. Later in development, during neurogenesis, the mutation spectrum shifted toward oxidative damage, and the mutation rate increased. Both neurogenesis and early embryogenesis exhibit substantially more mutagenesis than adulthood.

*K-means and cluster models for cancer signatures

We present *K-means clustering algorithm and source code by expanding statistical clustering methods applied in https://ssrn.com/abstract=2802753 to quantitative finance. *K-means is statistically deterministic without specifying initial centers, etc. We apply *K-means to extracting cancer signatures from genome data without using nonnegative matrix factorization (NMF). *K-means' computational cost is a fraction of NMF's. Using 1389 published samples for 14 cancer types, we find that 3 cancers (liver cancer, lung cancer and renal cell carcinoma) stand out and do not have cluster-like structures. Two clusters have especially high within-cluster correlations with 11 other cancers indicating common underlying structures. Our approach opens a novel avenue for studying such structures. *K-means is universal and can be applied in other fields. We discuss some potential applications in quantitative finance.

Base Excision Repair Variants in Cancer

Base excision repair (BER) is a key genome maintenance pathway that removes endogenously damaged DNA bases that arise in cells at very high levels on a daily basis. Failure to remove these damaged DNA bases leads to increased levels of mutagenesis and chromosomal instability, which have the potential to drive carcinogenesis. Next-generation sequencing of the germline and tumor genomes of thousands of individuals has uncovered many rare mutations in BER genes. Given that BER is critical for genome maintenance, it is important to determine whether BER genomic variants have functional phenotypes. In this chapter, we present our in silico methods for the identification and prioritization of BER variants for further study. We also provide detailed instructions and commentary on the initial cellular assays we employ to dissect potentially important phenotypes of human BER variants and highlight the strengths and weaknesses of our approaches. BER variants possessing interesting functional phenotypes can then be studied in more detail to provide important mechanistic insights regarding the role of aberrant BER in carcinogenesis.

Somatic mutations in cancer: Stochastic versus predictable

The origins of human cancers remain unclear except for a limited number of potent environmental mutagens, such as tobacco and UV light, and in rare cases, familial germ line mutations that affect tumor suppressor genes or oncogenes. A significant component of cancer etiology has been deemed stochastic and correlated with the number of stem cells in a tissue, the number of times the stem cells divide and a low incidence of random DNA polymerase errors that occur during each cell division. While somatic mutations occur during each round of DNA replication, mutations in cancer driver genes are not stochastic. Out of a total of 2843 codons, 1031 can be changed to stop codons by a single base substitution in the tumor suppressor APC gene, which is mutated in 76% of colorectal cancers (CRC). However, the nonsense mutations, which comprise 65% of all the APC driver mutations in CRC, are not random: 43% occur at Arg CGA codons, although they represent <3% of the codons. In TP53, CGA codons comprise <3% of the total 393 codons but they account for 72% and 39% of the mutations in CRC and ovarian cancer OVC, respectively. This mutation pattern is consistent with the kinetically slow, but not stochastic, hydrolytic deamination of 5-methylcytosine residues at specific methylated CpG sites to afford T·G mismatches that lead to C→T transitions and stop codons at CGA. Analysis of nonsense mutations in CRC, OVC and a number of other cancers indicates the need to expand the predictable risk factors for cancer to include, in addition to random polymerase errors, the methylation status of gene body CGA codons in tumor suppressor genes.

Piecing together the puzzle: nanopore technology in detection and quantification of cancer biomarkers

This mini-review paper is a comprehensive outline of nanopore technology applications in the detection and study of various cancer causal factors.

Mutation tendency of mutator Plasmodium berghei with proofreading-deficient DNA polymerase δ

In this study, we investigated the mutation tendency of a mutator rodent malaria parasite, Plasmodium berghei, with proofreading-deficient DNA polymerase δ. Wild-type and mutator parasites were maintained in mice for over 24 weeks, and the genome-wide accumulated mutations were determined by high-throughput sequencing. The mutator P. berghei had a significant preference for C/G to A/T substitutions; thus, its genome had a trend towards a higher AT content. The mutation rate was influenced by the sequence context, and mutations were markedly elevated at TCT. Some genes mutated repeatedly in replicate passage lines. In particular, knockout mutations of the AP2-G gene were frequent, which conferred strong growth advantages on parasites during the blood stage but at the cost of losing the ability to form gametocytes. This is the first report to demonstrate a biased mutation tendency in malaria parasites, and its results help to promote our basic understanding of Plasmodium genetics.

Clusters of Multiple Mutations: Incidence and Molecular Mechanisms

It has been long understood that mutation distribution is not completely random across genomic space and in time. Indeed, recent surprising discoveries identified multiple simultaneous mutations occurring in tiny regions within chromosomes while the rest of the genome remains relatively mutation-free. Mechanistic elucidation of these phenomena, called mutation showers, mutation clusters, or kataegis, in parallel with findings of abundant clustered mutagenesis in cancer genomes, is ongoing. So far, the combination of factors most important for clustered mutagenesis is the induction of DNA lesions within unusually long and persistent single-strand DNA intermediates. In addition to being a fascinating phenomenon, clustered mutagenesis also became an indispensable tool for identifying a previously unrecognized major source of mutation in cancer, APOBEC cytidine deaminases. Future research on clustered mutagenesis may shed light onto important mechanistic details of genome maintenance, with potentially profound implications for human health.

Mutation Processes in 293-Based Clones Overexpressing the DNA Cytosine Deaminase APOBEC3B

Molecular, cellular, and clinical studies have combined to demonstrate a contribution from the DNA cytosine deaminase APOBEC3B (A3B) to the overall mutation load in breast, head/neck, lung, bladder, cervical, ovarian, and other cancer types. However, the complete landscape of mutations attributable to this enzyme has yet to be determined in a controlled human cell system. We report a conditional and isogenic system for A3B induction, genomic DNA deamination, and mutagenesis. Human 293-derived cells were engineered to express doxycycline-inducible A3B-eGFP or eGFP constructs. Cells were subjected to 10 rounds of A3B-eGFP exposure that each caused 80–90% cell death. Control pools were subjected to parallel rounds of non-toxic eGFP exposure, and dilutions were done each round to mimic A3B-eGFP induced population fluctuations. Targeted sequencing of portions of TP53 and MYC demonstrated greater mutation accumulation in the A3B-eGFP exposed pools. Clones were generated and microarray analyses were used to identify those with the greatest number of SNP alterations for whole genome sequencing. A3B-eGFP exposed clones showed global increases in C-to-T transition mutations, enrichments for cytosine mutations within A3B-preferred trinucleotide motifs, and more copy number aberrations. Surprisingly, both control and A3B-eGFP clones also elicited strong mutator phenotypes characteristic of defective mismatch repair. Despite this additional mutational process, the 293-based system characterized here still yielded a genome-wide view of A3B-catalyzed mutagenesis in human cells and a system for additional studies on the compounded effects of simultaneous mutation mechanisms in cancer cells.

Computational DNA hole spectroscopy: A new tool to predict mutation hotspots, critical base pairs, and disease 'driver' mutations

We report on a new technique, computational DNA hole spectroscopy, which creates spectra of electron hole probabilities vs. nucleotide position. A hole is a site of positive charge created when an electron is removed. Peaks in the hole spectrum depict sites where holes tend to localize and potentially trigger a base pair mismatch during replication. Our studies of mitochondrial DNA reveal a correlation between L-strand hole spectrum peaks and spikes in the human mutation spectrum. Importantly, we also find that hole peak positions that do not coincide with large variant frequencies often coincide with disease-implicated mutations and/or (for coding DNA) encoded conserved amino acids. This enables combining hole spectra with variant data to identify critical base pairs and potential disease 'driver' mutations. Such integration of DNA hole and variance spectra could ultimately prove invaluable for pinpointing critical regions of the vast non-protein-coding genome. An observed asymmetry in correlations, between the spectrum of human mtDNA variations and the L- and H-strand hole spectra, is attributed to asymmetric DNA replication processes that occur for the leading and lagging strands.

The somatic autosomal mutation matrix in cancer genomes

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Disruption of Transcriptional Coactivator Sub1 Leads to Genome-Wide Re-distribution of Clustered Mutations Induced by APOBEC in Active Yeast Genes

Mutations in genomes of species are frequently distributed non-randomly, resulting in mutation clusters, including recently discovered kataegis in tumors. DNA editing deaminases play the prominent role in the etiology of these mutations. To gain insight into the enigmatic mechanisms of localized hypermutagenesis that lead to cluster formation, we analyzed the mutational single nucleotide variations (SNV) data obtained by whole-genome sequencing of drug-resistant mutants induced in yeast diploids by AID/APOBEC deaminase and base analog 6-HAP. Deaminase from sea lamprey, PmCDA1, induced robust clusters, while 6-HAP induced a few weak ones. We found that PmCDA1, AID, and APOBEC1 deaminases preferentially mutate the beginning of the actively transcribed genes. Inactivation of transcription initiation factor Sub1 strongly reduced deaminase-induced can1 mutation frequency, but, surprisingly, did not decrease the total SNV load in genomes. However, the SNVs in the genomes of the sub1 clones were re-distributed, and the effect of mutation clustering in the regions of transcription initiation was even more pronounced. At the same time, the mutation density in the protein-coding regions was reduced, resulting in the decrease of phenotypically detected mutants. We propose that the induction of clustered mutations by deaminases involves: a) the exposure of ssDNA strands during transcription and loss of protection of ssDNA due to the depletion of ssDNA-binding proteins, such as Sub1, and b) attainment of conditions favorable for APOBEC action in subpopulation of cells, leading to enzymatic deamination within the currently expressed genes. This model is applicable to both the initial and the later stages of oncogenic transformation and explains variations in the distribution of mutations and kataegis events in different tumor cells.

Genomes of tumors are heavily enriched with mutations. Some of these mutations are distributed non-randomly, forming mutational clusters. Editing cytosine deaminases from APOBEC superfamily are responsible for the formation of many of these clusters. We have expressed APOBEC enzyme in diploid yeast cells and found that most of the mutations occur in the beginning of the active genes, where transcription starts. Clusters of mutations overlapped with promoters/transcription start sites. This is likely due to the weaker protection of ssDNA, an ultimate APOBEC deaminase enzyme target, in the beginning of the genes. This hypothesis was reinforced by the finding that inactivation of Sub1 transcription initiation factor, which is found predominantly in the regions of transcription initiation, leads to further increase in mutagenesis in the beginning of the genes. Interestingly, the total number of mutations in the genomes of Sub1-deficient clones did not change, despite the 100-fold decrease in frequency of mutants in a reporter gene. Thus, the drastic change in genome-wide distribution of mutations can be caused by inactivation of a single gene. We propose that the loss of ssDNA protection factors causes formation of mutation clusters in human cancer.

APOBEC family mutational signatures are associated with poor prognosis translocations in multiple myeloma

We have sequenced 463 presenting cases of myeloma entered into the UK Myeloma XI study using whole exome sequencing. Here we identify mutations induced as a consequence of misdirected AID in the partner oncogenes of IGH translocations, which are activating and associated with impaired clinical outcome. An APOBEC mutational signature is seen in 3.8% of cases and is linked to the translocation-mediated deregulation of MAF and MAFB, a known poor prognostic factor. Patients with this signature have an increased mutational load and a poor prognosis. Loss of MAF or MAFB expression results in decreased APOBEC3B and APOBEC4 expression, indicating a transcriptional control mechanism. Kataegis, a further mutational pattern associated with APOBEC deregulation, is seen at the sites of the MYC translocation. The APOBEC mutational signature seen in myeloma is, therefore, associated with poor prognosis primary and secondary translocations and the molecular mechanisms involved in generating them.

Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes

Single base substitutions (SBSs) and insertions/deletions are critical for generating population diversity and can lead both to inherited disease and cancer. Whereas on a genome-wide scale SBSs are influenced by cellular factors, on a fine scale SBSs are influenced by the local DNA sequence-context, although the role of flanking sequence is often unclear. Herein, we used bioinformatics, molecular dynamics and hybrid quantum mechanics/molecular mechanics to analyze sequence context-dependent mutagenesis at mononucleotide repeats (A-tracts and G-tracts) in human population variation and in cancer genomes. SBSs and insertions/deletions occur predominantly at the first and last base-pairs of A-tracts, whereas they are concentrated at the second and third base-pairs in G-tracts. These positions correspond to the most flexible sites along A-tracts, and to sites where a ‘hole’, generated by the loss of an electron through oxidation, is most likely to be localized in G-tracts. For A-tracts, most SBSs occur in the direction of the base-pair flanking the tracts. We conclude that intrinsic features of local DNA structure, i.e. base-pair flexibility and charge transfer, render specific nucleotides along mononucleotide runs susceptible to base modification, which then yields mutations. Thus, local DNA dynamics contributes to phenotypic variation and disease in the human population.

Molecular mechanism and clinical impact of APOBEC3B-catalyzed mutagenesis in breast cancer

Cancer genomic DNA sequences enable identification of all mutations and suggest targets for precision medicine. The identities and patterns of the mutations themselves also provide critical information for deducing the originating DNA damaging agents, causal molecular mechanisms, and thus additional therapeutic targets. A classic example is ultraviolet light, which crosslinks adjacent pyrimidines and leads to C-to-T transitions. A new example is the DNA cytosine deaminase APOBEC3B, which was identified recently as a source of DNA damage and mutagenesis in breast, head/neck, cervix, bladder, lung, ovary, and to lesser extents additional cancer types. This enzyme is normally an effector protein in the innate immune response to virus infection but upregulation in these cancer types causes elevated levels of genomic C-to-U deamination events, which manifest as C-to-T transitions and C-to-G transversions within distinct DNA trinucleotide contexts (preferentially 5’-TCA and 5’-TCG). Genomic C-to-U deamination events within the same trinucleotide contexts also lead to cytosine mutation clusters (kataegis), and may precipitate visible chromosomal aberrations such as translocations. Clinical studies indicate that APOBEC3B upregulation correlates with poorer outcomes for estrogen receptor-positive breast cancer patients, including shorter durations of disease-free survival and overall survival after surgery. APOBEC3B may therefore have both diagnostic and prognostic potential. APOBEC3B may also be a candidate for therapeutic targeting because inhibition of this non-essential enzyme is predicted to decrease tumor mutation rates and diminish the likelihood of undesirable mutation-dependent outcomes such as recurrence, metastasis, and the development of therapy resistant tumors.

Base excision repair capacity in informing healthspan

Base excision repair (BER) is a frontline defense mechanism for dealing with many common forms of endogenous DNA damage, several of which can drive mutagenic or cell death outcomes. The pathway engages proteins such as glycosylases, abasic endonucleases, polymerases and ligases to remove substrate modifications from DNA and restore the genome back to its original state. Inherited mutations in genes related to BER can give rise to disorders involving cancer, immunodeficiency and neurodegeneration. Studies employing genetically defined heterozygous (haploinsufficient) mouse models indicate that partial reduction in BER capacity can increase vulnerability to both spontaneous and exposure-dependent pathologies. In humans, measurement of BER variation has been imperfect to this point, yet tools to assess BER in epidemiological surveys are steadily evolving. We provide herein an overview of the BER pathway and discuss the current efforts toward defining the relationship of BER defects with disease susceptibility.

Hypermutation in human cancer genomes: footprints and mechanisms

A role for somatic mutations in carcinogenesis is well accepted, but the degree to which mutation rates influence cancer initiation and development is under continuous debate. Recently accumulated genomic data have revealed that thousands of tumour samples are riddled by hypermutation, broadening support for the idea that many cancers acquire a mutator phenotype. This major expansion of cancer mutation data sets has provided unprecedented statistical power for the analysis of mutation spectra, which has confirmed several classical sources of mutation in cancer, highlighted new prominent mutation sources (such as apolipoprotein B mRNA editing enzyme catalytic polypeptide-like (APOBEC) enzymes) and empowered the search for cancer drivers. The confluence of cancer mutation genomics and mechanistic insight provides great promise for understanding the basic development of cancer through mutations.

Nanopores discriminate among five C5-cytosine variants in DNA

Individual DNA molecules can be read at single nucleotide precision using nanopores coupled to processive enzymes. Discrimination among the four canonical bases has been achieved, as has discrimination among cytosine, 5-methylcytosine (mC), and 5-hydroxymethylcytosine (hmC). Two additional modified cytosine bases, 5-carboxylcytosine (caC) and 5-formylcytosine (fC), are produced during enzymatic conversion of hmC to cytosine in mammalian cells. Thus, an accurate picture of the cytosine epigenetic status in target cells should also include these C5-cytosine variants. In the present study, we used a patch clamp amplifier to acquire ionic current traces caused by phi29 DNA polymerase-controlled translocation of DNA templates through the M2MspA pore. Decision boundaries based on three consecutive ionic current states were implemented to call mC, hmC, caC, fC, or cytosine at CG dinucleotides in ∼4400 individual DNA molecules. We found that the percentage of correct base calls for single pass reads ranged from 91.6% to 98.3%. This accuracy depended upon the identity of nearest neighbor bases surrounding the CG dinucleotide.

The spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β

To provide molecular-level insights into the spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β (polβ), we report four crystal structures of polβ complexed with dG•dTTP and dA•dCTP mismatches in the presence of Mg2+ or Mn2+. The Mg(2+)-bound ground-state structures show that the dA•dCTP-Mg2+ complex adopts an 'intermediate' protein conformation while the dG•dTTP-Mg2+ complex adopts an open protein conformation. The Mn(2+)-bound 'pre-chemistry-state' structures show that the dA•dCTP-Mn2+ complex is structurally very similar to the dA•dCTP-Mg2+ complex, whereas the dG•dTTP-Mn2+ complex undergoes a large-scale conformational change to adopt a Watson-Crick-like dG•dTTP base pair and a closed protein conformation. These structural differences, together with our molecular dynamics simulation studies, suggest that polβ increases replication fidelity via a two-stage mismatch discrimination mechanism, where one is in the ground state and the other in the closed conformation state. In the closed conformation state, polβ appears to allow only a Watson-Crick-like conformation for purine•pyrimidine base pairs, thereby discriminating the mismatched base pairs based on their ability to form the Watson-Crick-like conformation. Overall, the present studies provide new insights into the spontaneous replication error and the replication fidelity mechanisms of polβ.

Melanoma's high C>T mutation rate: is deamination playing a role?

The majority of melanoma mutations are C>T transitions, and most bear UV signatures. However, other process may contribute to the high C>T mutation rate. Okura et al., have demonstrated immunohistochemical evidence of deaminating enzymes, activation-induced cytidine deaminase (AID) and apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3B (APOBEC3B) in melanoma. Both have been implicated in cancer. While further validation is necessary, these findings warrant consideration of a role for deamination in melanomagenesis. Deamination primarily drives C>T transitions. Compared with trunk/extremity melanomas, acral melanomas display a significantly higher percentage of 'spontaneous' and 'AID' mutation signature events suggesting deamination may be particularly important in this subgroup.

Replicative Stress and the FHIT Gene: Roles in Tumor Suppression, Genome Stability and Prevention of Carcinogenesis

The fragile FHIT gene, encompassing the chromosomal fragile site FRA3B, is an early target of DNA damage in precancerous cells. While vulnerable to DNA damage itself, FHIT protein expression is essential to protect from DNA damage-induced cancer initiation and progression by modulating genome stability, oxidative stress and levels of accumulating DNA damage. Thus, FHIT, whose expression is lost or reduced in many human cancers, is a tumor suppressor and genome caretaker whose loss initiates genome instability in preneoplastic lesions. Ongoing studies are seeking more detailed understanding of the role of FHIT in the cellular response to oxidative damage. This review discusses the relationship between FHIT, reactive oxygen species production, and DNA damage in the context of cancer initiation and progression.