Comprehensive evaluation of smoking exposures and their interactions on DNA methylation.

BACKGROUND: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS). METHODS: We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C). FINDINGS: Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs. INTERPRETATION: Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases. FUNDING: Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

Synergistic roles of tristetraprolin family members in myeloid cells in the control of inflammation.

Members of the tristetraprolin (TTP) family of RNA-binding proteins can bind to and promote the decay of specific transcripts containing AU-rich motifs. ZFP36 (TTP) is best known for regulating pro-inflammatory cytokine expression in myeloid cells; however, its mammalian paralogues ZFP36L1 and ZFP36L2 have not been viewed as important in controlling inflammation. We knocked out these genes in myeloid cells in mice, singly and together. Single-gene myeloid-specific knockouts resulted in almost no spontaneous phenotypes. In contrast, mice with myeloid cell deficiency of all three genes developed severe inflammation, with a median survival of 8 wk. Macrophages from these mice expressed many more stabilized transcripts than cells from myeloid-specific TTP knockout mice; many of these encoded pro-inflammatory cytokines and chemokines. The failure of weight gain, arthritis, and early death could be prevented completely by two normal alleles of any of the three paralogues, and even one normal allele of Zfp36 or Zfp36l2 was enough to prevent the inflammatory phenotype. Our findings emphasize the importance of all three family members, acting in concert, in myeloid cell function.

Vanadium Pentoxide Exposure Causes Strain-Dependent Changes in Mitochondrial DNA Heteroplasmy, Copy Number, and Lesions, but Not Nuclear DNA Lesions.

Interstitial lung diseases (ILDs) are lethal lung diseases characterized by pulmonary inflammation and progressive lung interstitial scarring. We previously developed a mouse model of ILD using vanadium pentoxide (V2O5) and identified several gene candidates on chromosome 4 associated with pulmonary fibrosis. While these data indicated a significant genetic contribution to ILD susceptibility, they did not include any potential associations and interactions with the mitochondrial genome that might influence disease risk. To conduct this pilot work, we selected the two divergent strains we previously categorized as V2O5-resistant C57BL6J (B6) and -responsive DBA/2J (D2) and compared their mitochondrial genome characteristics, including DNA variants, heteroplasmy, lesions, and copy numbers at 14- and 112-days post-exposure. While we did not find changes in the mitochondrial genome at 14 days post-exposure, at 112 days, we found that the responsive D2 strain exhibited significantly fewer mtDNA copies and more lesions than control animals. Alongside these findings, mtDNA heteroplasmy frequency decreased. These data suggest that mice previously shown to exhibit increased susceptibility to pulmonary fibrosis and inflammation sustain damage to the mitochondrial genome that is evident at 112 days post-V2O5 exposure.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans.

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Race/ethnicity-stratified fine-mapping of the MHC locus reveals genetic variants associated with late-onset asthma.

Introduction: Asthma is a chronic disease of the airways that impairs normal breathing. The etiology of asthma is complex and involves multiple factors, including the environment and genetics, especially the distinct genetic architecture associated with ancestry. Compared to early-onset asthma, little is known about genetic predisposition to late-onset asthma. We investigated the race/ethnicity-specific relationship among genetic variants within the major histocompatibility complex (MHC) region and late-onset asthma in a North Carolina-based multiracial cohort of adults. Methods: We stratified all analyses by self-reported race (i.e., White and Black) and adjusted all regression models for age, sex, and ancestry. We conducted association tests within the MHC region and performed fine-mapping analyses conditioned on the race/ethnicity-specific lead variant using whole-genome sequencing (WGS) data. We applied computational methods to infer human leukocyte antigen (HLA) alleles and residues at amino acid positions. We replicated findings in the UK Biobank. Results: The lead signals, rs9265901 on the 5' end of HLA-B, rs55888430 on HLA-DOB, and rs117953947 on HCG17, were significantly associated with late-onset asthma in all, White, and Black participants, respectively (OR = 1.73, 95%CI: 1.31 to 2.14, p = 3.62 × 10-5; OR = 3.05, 95%CI: 1.86 to 4.98, p = 8.85 × 10-6; OR = 19.5, 95%CI: 4.37 to 87.2, p = 9.97 × 10-5, respectively). For the HLA analysis, HLA-B*40:02 and HLA-DRB1*04:05, HLA-B*40:02, HLA-C*04:01, and HLA-DRB1*04:05, and HLA-DRB1*03:01 and HLA-DQB1 were significantly associated with late-onset asthma in all, White, and Black participants. Conclusion: Multiple genetic variants within the MHC region were significantly associated with late-onset asthma, and the associations were significantly different by race/ethnicity group.

DUX4 double whammy: The transcription factor that causes a rare muscular dystrophy also kills the precursors of the human nose.

SMCHD1 mutations cause congenital arhinia (absent nose) and a muscular dystrophy called FSHD2. In FSHD2, loss of SMCHD1 repressive activity causes expression of double homeobox 4 (DUX4) in muscle tissue, where it is toxic. Studies of arhinia patients suggest a primary defect in nasal placode cells (human nose progenitors). Here, we show that upon SMCHD1 ablation, DUX4 becomes derepressed in H9 human embryonic stem cells (hESCs) as they differentiate toward a placode cell fate, triggering cell death. Arhinia and FSHD2 patient-derived induced pluripotent stem cells (iPSCs) express DUX4 when converted to placode cells and demonstrate variable degrees of cell death, suggesting an environmental disease modifier. HSV-1 may be one such modifier as herpesvirus infection amplifies DUX4 expression in SMCHD1 KO hESC and patient iPSC. These studies suggest that arhinia, like FSHD2, is due to compromised SMCHD1 repressive activity in a cell-specific context and provide evidence for an environmental modifier.

Pif1 family helicases promote mutation avoidance during DNA replication.

Pif1 family 5' → 3' DNA helicases are important for replication fork progression and genome stability. The budding yeast Saccharomyces cerevisiae encodes two Pif1 family helicases, Rrm3 and Pif1, both of which are multi-functional. Here we describe novel functions for Rrm3 in promoting mutation avoidance during DNA replication. We show that loss of RRM3 results in elevated spontaneous mutations made by DNA polymerases Pols ϵ and δ, which are subject to DNA mismatch repair. The absence of RRM3 also causes higher mutagenesis by the fourth B-family DNA polymerase Pol ζ. By genome-wide analysis, we show that the mutational consequences due to loss of RRM3 vary depending on the genomic locus. Rrm3 promotes the accuracy of DNA replication by Pols ϵ and δ across the genome, and it is particularly important for preventing Pol ζ-dependent mutagenesis at tRNA genes. In addition, mutation avoidance by Rrm3 depends on its helicase activity, and Pif1 serves as a backup for Rrm3 in suppressing mutagenesis. We present evidence that the sole human Pif1 family helicase in human cells likely also promotes replication fidelity, suggesting that a role for Pif1 family helicases in mutation avoidance may be evolutionarily conserved, a possible underlying mechanism for its potential tumor-suppressor function.

How asymmetric DNA replication achieves symmetrical fidelity.

Accurate DNA replication of an undamaged template depends on polymerase selectivity for matched nucleotides, exonucleolytic proofreading of mismatches, and removal of remaining mismatches via DNA mismatch repair (MMR). DNA polymerases (Pols) δ and ε have 3'-5' exonucleases into which mismatches are partitioned for excision in cis (intrinsic proofreading). Here we provide strong evidence that Pol δ can extrinsically proofread mismatches made by itself and those made by Pol ε, independently of both Pol δ's polymerization activity and MMR. Extrinsic proofreading across the genome is remarkably efficient. We report, with unprecedented accuracy, in vivo contributions of nucleotide selectivity, proofreading, and MMR to the fidelity of DNA replication in Saccharomyces cerevisiae. We show that extrinsic proofreading by Pol δ improves and balances the fidelity of the two DNA strands. Together, we depict a comprehensive picture of how nucleotide selectivity, proofreading, and MMR cooperate to achieve high and symmetrical fidelity on the two strands.

A post-transcriptional regulon controlled by TtpA, the single tristetraprolin family member expressed in Dictyostelium discoideum.

Post-transcriptional processes mediated by mRNA binding proteins represent important control points in gene expression. In eukaryotes, mRNAs containing specific AU-rich motifs are regulated by binding of tristetraprolin (TTP) family tandem zinc finger proteins, which promote mRNA deadenylation and decay, partly through interaction of a conserved C-terminal CNOT1 binding (CNB) domain with CCR4-NOT protein complexes. The social ameba Dictyostelium discoideum shared a common ancestor with humans more than a billion years ago, and expresses only one TTP family protein, TtpA, in contrast to three members expressed in humans. Evaluation of ttpA null-mutants identified six transcripts that were consistently upregulated compared to WT during growth and early development. The 3'-untranslated regions (3'-UTRs) of all six 'TtpA-target' mRNAs contained multiple TTP binding motifs (UUAUUUAUU), and one 3'-UTR conferred TtpA post-transcriptional stability regulation to a heterologous mRNA that was abrogated by mutations in the core TTP-binding motifs. All six target transcripts were upregulated to similar extents in a C-terminal truncation mutant, in contrast to less severe effects of analogous mutants in mice. All six target transcripts encoded probable membrane proteins. In Dictyostelium, TtpA may control an 'RNA regulon', where a single RNA binding protein, TtpA, post-transcriptionally co-regulates expression of several functionally related proteins.

UV-exposure, endogenous DNA damage, and DNA replication errors shape the spectra of genome changes in human skin.

Human skin is continuously exposed to environmental DNA damage leading to the accumulation of somatic mutations over the lifetime of an individual. Mutagenesis in human skin cells can be also caused by endogenous DNA damage and by DNA replication errors. The contributions of these processes to the somatic mutation load in the skin of healthy humans has so far not been accurately assessed because the low numbers of mutations from current sequencing methodologies preclude the distinction between sequencing errors and true somatic genome changes. In this work, we sequenced genomes of single cell-derived clonal lineages obtained from primary skin cells of a large cohort of healthy individuals across a wide range of ages. We report here the range of mutation load and a comprehensive view of the various somatic genome changes that accumulate in skin cells. We demonstrate that UV-induced base substitutions, insertions and deletions are prominent even in sun-shielded skin. In addition, we detect accumulation of mutations due to spontaneous deamination of methylated cytosines as well as insertions and deletions characteristic of DNA replication errors in these cells. The endogenously induced somatic mutations and indels also demonstrate a linear increase with age, while UV-induced mutation load is age-independent. Finally, we show that DNA replication stalling at common fragile sites are potent sources of gross chromosomal rearrangements in human cells. Thus, somatic mutations in skin of healthy individuals reflect the interplay of environmental and endogenous factors in facilitating genome instability and carcinogenesis.

Increased Burden of Rare Sequence Variants in GnRH-Associated Genes in Women With Hypothalamic Amenorrhea.

CONTEXT: Functional hypothalamic amenorrhea (HA) is a common, acquired form of hypogonadotropic hypogonadism that occurs in the setting of energy deficits and/or stress. Variability in individual susceptibility to these stressors, HA heritability, and previous identification of several rare sequence variants (RSVs) in genes associated with the rare disorder, isolated hypogonadotropic hypogonadism (IHH), in individuals with HA suggest a possible genetic contribution to HA susceptibility. OBJECTIVE: We sought to determine whether the burden of RSVs in IHH-related genes is greater in women with HA than controls. DESIGN: We compared patients with HA to control women. SETTING: The study was conducted at secondary referral centers. PATIENTS AND OTHER PARTICIPANTS: Women with HA (n = 106) and control women (ClinSeq study; n = 468). INTERVENTIONS: We performed exome sequencing in all patients and controls. MAIN OUTCOME MEASURE(S): The frequency of RSVs in 53 IHH-associated genes was determined using rare variant burden and association tests. RESULTS: RSVs were overrepresented in women with HA compared with controls (P = .007). Seventy-eight heterozygous RSVs in 33 genes were identified in 58 women with HA (36.8% of alleles) compared to 255 RSVs in 41 genes among 200 control women (27.2%). CONCLUSIONS: Women with HA are enriched for RSVs in genes that cause IHH, suggesting that variation in genes associated with gonadotropin-releasing hormone neuronal ontogeny and function may be a major determinant of individual susceptibility to developing HA in the face of diet, exercise, and/or stress.

Association between Mitochondrial DNA Sequence Variants and V˙O2 max Trainability.

PURPOSE: We designed the study to determine whether mitochondrial DNA (mtDNA) haplogroup, sequence, and heteroplasmy differed between individuals previously characterized as low (LR) or high responders (HR) as defined by their maximal oxygen uptake response to a standardized aerobic exercise training program. METHODS: DNA was isolated from whole blood in subjects from the HERITAGE Family Study that were determined to be either HR (n = 15) or LR (n = 15). mtDNA was amplified by long-range polymerase chain reaction, then tagged with Nextera libraries and sequenced on a MiSeq instrument. RESULTS: Different mtDNA haplogroup subtypes were found in HR and LR individuals. Compared with HR subjects, significantly more LR subjects had variants in 13 sites, including 7 in hypervariable (HV) regions: HV2 (G185A: 0 vs 6, P = 0.02; G228A: 0 vs 5, P = 0.04; C295T: 0 vs 6; P = 0.04), HV3 (C462T: 0 vs 5, P = 0.04; T489C: 0 vs 5; P = 0.04), and HV1 (C16068T: 0 vs 6, P = 0.02; T16125C: 0 vs 6, P = 0.02). Remaining variants were in protein coding genes, mtND1 (1 vs 8, P = 0.02), mtND3 (A10397G: 0 vs 5, P = 0.04), mtND4 (A11250G: 1 vs 8, P = 0.02), mtND5 (G13707A: 0 vs 5, P = 0.04), and mtCYTB (T14797C: 0 vs 5, P = 0.04; C15451A: 1 vs 8, P = 0.02). Average total numbers of heteroplasmies (P = 0.83) and frequency of heteroplasmies (P = 0.05) were similar between the groups. CONCLUSIONS: Our findings provide specific sites across the mitochondrial genome that may be related to maximal oxygen uptake trainability.

H/ACA snoRNA levels are regulated during stem cell differentiation.

H/ACA small nucleolar RNAs (snoRNAs) guide pseudouridylation as part of a small nucleolar ribonucleoprotein complex (snoRNP). Disruption of H/ACA snoRNA levels in stem cells impairs pluripotency, yet it remains unclear how H/ACA snoRNAs contribute to differentiation. To determine if H/ACA snoRNA levels are dynamic during differentiation, we comprehensively profiled H/ACA snoRNA abundance in multiple murine cell types and during differentiation in three cellular models, including mouse embryonic stem cells and mouse myoblasts. We determined that the profiles of H/ACA snoRNA abundance are cell-type specific, and we identified a subset of snoRNAs that are specifically regulated during differentiation. Additionally, we demonstrated that a decrease in Snora27 abundance upon differentiation corresponds to a decrease in pseudouridylation of its target site within the E-site transfer RNA (tRNA) binding region of the 28S ribosomal RNA (rRNA) in the large ribosomal subunit. Together, these data point toward a potential model in which H/ACA snoRNAs are specifically regulated during differentiation to alter pseudouridylation and fine tune ribosome function.

Mutation signatures specific to DNA alkylating agents in yeast and cancers.

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.

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

[This corrects the article DOI: 10.1371/journal.pbio.3000464.].

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.

Roles for DNA polymerase δ in initiating and terminating leading strand DNA replication.

Most current evidence indicates that DNA polymerases ε and δ, respectively, perform the bulk of leading and lagging strand replication of the eukaryotic nuclear genome. Given that ribonucleotide and mismatch incorporation rates by these replicases influence somatic and germline patterns of variation, it is important to understand the details and exceptions to this overall division of labor. Using an improved method to map where these replicases incorporate ribonucleotides during replication, here we present evidence that DNA polymerase δ universally participates in initiating leading strand synthesis and that nascent leading strand synthesis switches from Pol ε to Pol δ during replication termination. Ribonucleotide maps from both the budding and fission yeast reveal conservation of these processes. These observations of replisome dynamics provide important insight into the mechanisms of eukaryotic replication and genome maintenance.

Multi-walled carbon nanotubes upregulate mitochondrial gene expression and trigger mitochondrial dysfunction in primary human bronchial epithelial cells.

Nanomaterials are a relatively new class of materials that acquire novel properties based on their reduced size. While these materials have widespread use in consumer products and industrial applications, the potential health risks associated with exposure to them remain to be fully characterized. Carbon nanotubes are among the most widely used nanomaterials and have high potential for human exposure by inhalation. These nanomaterials are known to penetrate the cell membrane and interact with intracellular molecules, resulting in a multitude of documented effects, including oxidative stress, genotoxicity, impaired metabolism, and apoptosis. While the capacity for carbon nanotubes to damage nuclear DNA has been established, the effect of exposure on mitochondrial DNA (mtDNA) is relatively unexplored. In this study, we investigated the potential of multi-walled carbon nanotubes (MWCNTs) to impair mitochondrial gene expression and function in human bronchial epithelial cells (BECs). Primary BECs were exposed to sub-cytotoxic doses (up to 3 µg/ml) of MWCNTs for 5 d and assessed for changes in expression of all mitochondrial protein-coding genes, heteroplasmies, and insertion/deletion mutations (indels). Exposed cells were also measured for cytotoxicity, metabolic function, mitochondrial abundance, and mitophagy. We found that MWCNTs upregulated mitochondrial gene expression, while significantly decreasing oxygen consumption rate and mitochondrial abundance. Confocal microscopy revealed induction of mitophagy by 2 hours of exposure. Mitochondrial DNA heteroplasmy and insertion/deletion mutations were not significantly affected by any treatment. We conclude that carbon nanotubes cause mitochondrial dysfunction that leads to mitophagy in exposed BECs via a mechanism unrelated to its reported genotoxicity.

DNA Polymerase Delta Synthesizes Both Strands during Break-Induced Replication.

Break-induced replication (BIR) is a pathway of homology-directed repair that repairs one-ended DNA breaks, such as those formed at broken replication forks or uncapped telomeres. In contrast to conventional S phase DNA synthesis, BIR proceeds by a migrating D-loop and results in conservative synthesis of the nascent strands. DNA polymerase delta (Pol δ) initiates BIR; however, it is not known whether synthesis of the invading strand switches to a different polymerase or how the complementary strand is synthesized. By using alleles of the replicative DNA polymerases that are permissive for ribonucleotide incorporation, thus generating a signature of their action in the genome that can be identified by hydrolytic end sequencing, we show that Pol δ replicates both the invading and the complementary strand during BIR. In support of this conclusion, we show that depletion of Pol δ from cells reduces BIR, whereas depletion of Pol ε has no effect.

Genome-wide mutagenesis resulting from topoisomerase 1-processing of unrepaired ribonucleotides in DNA.

Ribonucleotides are the most common non-canonical nucleotides incorporated into DNA during replication, and their processing leads to mutations and genome instability. Yeast mutation reporter systems demonstrate that 2-5 base pair deletions (Δ2-5bp) in repetitive DNA are a signature of unrepaired ribonucleotides, and that these events are initiated by topoisomerase 1 (Top1) cleavage. However, a detailed understanding of the frequency and locations of ribonucleotide-dependent mutational events across the genome has been lacking. Here we present the results of genome-wide mutational analysis of yeast strains deficient in Ribonucleotide Excision Repair (RER). We identified mutations that accumulated over thousands of generations in strains expressing either wild-type or variant replicase alleles (M644G Pol ε, L612M Pol δ, L868M Pol α) that confer increased ribonucleotide incorporation into DNA. Using a custom-designed mutation-calling pipeline called muver (for mutationes verificatae), we observe a number of surprising mutagenic features. This includes a 24-fold preferential elevation of AG and AC relative to AT dinucleotide deletions in the absence of RER, suggesting specificity for Top1-initiated deletion mutagenesis. Moreover, deletion rates in di- and trinucleotide repeat tracts increase exponentially with tract length. Consistent with biochemical and reporter gene mutational analysis, these deletions are no longer observed upon deletion of TOP1. Taken together, results from these analyses demonstrate the global impact of genomic ribonucleotide processing by Top1 on genome integrity.