Detection of ribonucleotides embedded in DNA by Nanopore sequencing.

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Light-strand bias and enriched zones of embedded ribonucleotides are associated with DNA replication and transcription in the human-mitochondrial genome.

Abundant ribonucleoside-triphosphate (rNTP) incorporation into DNA by DNA polymerases in the form of ribonucleoside monophosphates (rNMPs) is a widespread phenomenon in nature, resulting in DNA-structural change and genome instability. The rNMP distribution, characteristics, hotspots and association with DNA metabolic processes in human mitochondrial DNA (hmtDNA) remain mostly unknown. Here, we utilize the ribose-seq technique to capture embedded rNMPs in hmtDNA of six different cell types. In most cell types, the rNMPs are preferentially embedded on the light strand of hmtDNA with a strong bias towards rCMPs; while in the liver-tissue cells, the rNMPs are predominately found on the heavy strand. We uncover common rNMP hotspots and conserved rNMP-enriched zones across the entire hmtDNA, including in the control region, which links the rNMP presence to the frequent hmtDNA replication-failure events. We show a strong correlation between coding-sequence size and rNMP-embedment frequency per nucleotide on the non-template, light strand in all cell types, supporting the presence of transient RNA-DNA hybrids preceding light-strand replication. Moreover, we detect rNMP-embedment patterns that are only partly conserved across the different cell types and are distinct from those found in yeast mtDNA. The study opens new research directions to understand the biology of hmtDNA and genomic rNMPs.

Synthetic lethal mutants in Escherichia coli define pathways necessary for survival with RNase H deficiency.

Ribonucleotides frequently contaminate DNA and, if not removed, cause genomic instability. Consequently, all organisms are equipped with RNase H enzymes to remove RNA-DNA hybrids (RDHs). Escherichia coli lacking RNase HI (rnhA) and RNase HII (rnhB) enzymes, the ∆rnhA ∆rnhB double mutant, accumulates RDHs in its DNA. These RDHs can convert into RNA-containing DNA lesions (R-lesions) of unclear nature that compromise genomic stability. The ∆rnhAB double mutant has severe phenotypes, like growth inhibition, replication stress, sensitivity to ultraviolet radiation, SOS induction, increased chromosomal fragmentation, and defects in nucleoid organization. In this study, we found that RNase HI deficiency also alters wild-type levels of DNA supercoiling. Despite these severe chromosomal complications, ∆rnhAB double mutant survives, suggesting that dedicated pathways operate to avoid or repair R-lesions. To identify these pathways, we systematically searched for mutants synthetic lethal (colethal) with the rnhAB defect using an unbiased color screen and a candidate gene approach. We identified both novel and previously reported rnhAB-colethal and -coinhibited mutants, characterized them, and sorted them into avoidance or repair pathways. These mutants operate in various parts of nucleic acid metabolism, including replication fork progression, R-loop prevention and removal, nucleoid organization, tRNA modification, recombinational repair, and chromosome-dimer resolution, demonstrating the pleiotropic nature of RNase H deficiency. IMPORTANCE Ribonucleotides (rNs) are structurally very similar to deoxyribonucleotides. Consequently, rN contamination of DNA is common and pervasive across all domains of life. Failure to remove rNs from DNA has severe consequences, and all organisms are equipped with RNase H enzymes to remove RNA-DNA hybrids. RNase H deficiency leads to complications in bacteria, yeast, and mouse, and diseases like progressive external ophthalmoplegia (mitochondrial defects in RNASEH1) and Aicardi-Goutières syndrome (defects in RNASEH2) in humans. Escherichia coli ∆rnhAB mutant, deficient in RNases H, has severe chromosomal complications. Despite substantial problems, nearly half of the mutant population survives. We have identified novel and previously confirmed pathways in various parts of nucleic acid metabolism that ensure survival with RNase H deficiency.

NME6: ribonucleotide salvage sustains mitochondrial transcription.

The building blocks for RNA and DNA are made in the cytosol, meaning mitochondria depend on the import and salvage of ribonucleoside triphosphates (rNTPs) and deoxyribonucleoside triphosphates (dNTPs) for the synthesis of their own genetic material. While extensive research has focused on mitochondrial dNTP homeostasis due to its defects being associated with various mitochondrial DNA (mtDNA) depletion and deletion syndromes, the investigation of mitochondrial rNTP homeostasis has received relatively little attention. In this issue of the EMBO Journal, Grotehans et al provide compelling evidence of a major role for NME6, a mitochondrial nucleoside diphosphate kinase, in the conversion of pyrimidine ribonucleoside diphosphates into the corresponding triphosphates. These data also suggest a significant physiological role for NME6, as its absence results in the depletion of mitochondrial transcripts and destabilization of the electron transport chain (Grotehans et al, 2023).

Ribonucleotide synthesis by NME6 fuels mitochondrial gene expression.

Replication of the mitochondrial genome and expression of the genes it encodes both depend on a sufficient supply of nucleotides to mitochondria. Accordingly, dysregulated nucleotide metabolism not only destabilises the mitochondrial genome, but also affects its transcription. Here, we report that a mitochondrial nucleoside diphosphate kinase, NME6, supplies mitochondria with pyrimidine ribonucleotides that are necessary for the transcription of mitochondrial genes. Loss of NME6 function leads to the depletion of mitochondrial transcripts, as well as destabilisation of the electron transport chain and impaired oxidative phosphorylation. These deficiencies are rescued by an exogenous supply of pyrimidine ribonucleosides. Moreover, NME6 is required for the maintenance of mitochondrial DNA when the access to cytosolic pyrimidine deoxyribonucleotides is limited. Our results therefore reveal an important role for ribonucleotide salvage in mitochondrial gene expression.

Promoter-independent synthesis of chemically modified RNA by human DNA polymerase θ variants.

Synthetic RNA oligonucleotides composed of canonical and modified ribonucleotides are highly effective for RNA antisense therapeutics and RNA-based genome engineering applications utilizing CRISPR-Cas9. Yet, synthesis of synthetic RNA using phosphoramidite chemistry is highly inefficient and expensive relative to DNA oligonucleotides, especially for relatively long RNA oligonucleotides. Thus, new biotechnologies are needed to significantly reduce costs, while increasing synthesis rates and yields of synthetic RNA. Here, we engineer human DNA polymerase theta (Polθ) variants and demonstrate their ability to synthesize long (95-200 nt) RNA oligonucleotides with canonical ribonucleotides and ribonucleotide analogs commonly used for stabilizing RNA for therapeutic and genome engineering applications. In contrast to natural promoter-dependent RNA polymerases, Polθ variants synthesize RNA by initiating from DNA or RNA primers, which enables the production of RNA without short abortive byproducts. Remarkably, Polθ variants show the lower capacity to misincorporate ribonucleotides compared to T7 RNA polymerase. Automation of this enzymatic RNA synthesis technology can potentially increase yields while reducing costs of synthetic RNA oligonucleotide production.

Ribodysgenesis: sudden genome instability in the yeast Saccharomyces cerevisiae arising from RNase H2 cleavage at genomic-embedded ribonucleotides.

Ribonucleotides can be incorporated into DNA during replication by the replicative DNA polymerases. These aberrant DNA subunits are efficiently recognized and removed by Ribonucleotide Excision Repair, which is initiated by the heterotrimeric enzyme RNase H2. While RNase H2 is essential in higher eukaryotes, the yeast Saccharomyces cerevisiae can survive without RNase H2 enzyme, although the genome undergoes mutation, recombination and other genome instability events at an increased rate. Although RNase H2 can be considered as a protector of the genome from the deleterious events that can ensue from recognition and removal of embedded ribonucleotides, under conditions of high ribonucleotide incorporation and retention in the genome in a RNase H2-negative strain, sudden introduction of active RNase H2 causes massive DNA breaks and genome instability in a condition which we term 'ribodysgenesis'. The DNA breaks and genome instability arise solely from RNase H2 cleavage directed to the ribonucleotide-containing genome. Survivors of ribodysgenesis have massive loss of heterozygosity events stemming from recombinogenic lesions on the ribonucleotide-containing DNA, with increases of over 1000X from wild-type. DNA breaks are produced over one to two divisions and subsequently cells adapt to RNase H2 and ribonucleotides in the genome and grow with normal levels of genome instability.

Ribonucleotide Incorporation by Eukaryotic B-Family Replicases and Its Implications for Genome Stability.

Our current view of how DNA-based genomes are efficiently and accurately replicated continues to evolve as new details emerge on the presence of ribonucleotides in DNA. Ribonucleotides are incorporated during eukaryotic DNA replication at rates that make them the most common noncanonical nucleotide placed into the nuclear genome, they are efficiently repaired, and their removal impacts genome integrity. This review focuses on three aspects of this subject: the incorporation of ribonucleotides into the eukaryotic nuclear genome during replication by B-family DNA replicases, how these ribonucleotides are removed, and the consequences of their presence or removal for genome stability and disease.

Signatures of TOP1 transcription-associated mutagenesis in cancer and germline.

The mutational landscape is shaped by many processes. Genic regions are vulnerable to mutation but are preferentially protected by transcription-coupled repair1. In microorganisms, transcription has been demonstrated to be mutagenic2,3; however, the impact of transcription-associated mutagenesis remains to be established in higher eukaryotes4. Here we show that ID4-a cancer insertion-deletion (indel) mutation signature of unknown aetiology5 characterized by short (2 to 5 base pair) deletions -is due to a transcription-associated mutagenesis process. We demonstrate that defective ribonucleotide excision repair in mammals is associated with the ID4 signature, with mutations occurring at a TNT sequence motif, implicating topoisomerase 1 (TOP1) activity at sites of genome-embedded ribonucleotides as a mechanistic basis. Such TOP1-mediated deletions occur somatically in cancer, and the ID-TOP1 signature is also found in physiological settings, contributing to genic de novo indel mutations in the germline. Thus, although topoisomerases protect against genome instability by relieving topological stress6, their activity may also be an important source of mutations in the human genome.

Mapping ribonucleotides embedded in genomic DNA to single-nucleotide resolution using Ribose-Map.

The most common nonstandard nucleotides found in genomic DNA are ribonucleotides. Although ribonucleotides are frequently incorporated into DNA by replicative DNA polymerases, very little is known about the distribution and signatures of ribonucleotides incorporated into DNA. Recent advances in high-throughput ribonucleotide sequencing can capture the exact locations of ribonucleotides in genomic DNA. Ribose-Map is a user-friendly, standardized bioinformatics toolkit for the comprehensive analysis of ribonucleotide sequencing experiments. It allows researchers to map the locations of ribonucleotides in DNA to single-nucleotide resolution and identify biological signatures of ribonucleotide incorporation. In addition, it can be applied to data generated using any currently available high-throughput ribonucleotide sequencing technique, thus standardizing the analysis of ribonucleotide sequencing experiments and allowing direct comparisons of results. This protocol describes in detail how to use Ribose-Map to analyze ribonucleotide sequencing data, including preparing the reads for analysis, locating the genomic coordinates of ribonucleotides, exploring the genome-wide distribution of ribonucleotides, determining the nucleotide sequence context of ribonucleotides and identifying hotspots of ribonucleotide incorporation. Ribose-Map does not require background knowledge of ribonucleotide sequencing analysis and assumes only basic command-line skills. The protocol requires less than 3 h of computing time for most datasets and ~30 min of hands-on time. Ribose-Map is available at https://github.com/agombolay/ribose-map .

Protein Kinase D-Dependent Downregulation of Immediate Early Genes through Class IIA Histone Deacetylases in Acute Lymphoblastic Leukemia.

Acute lymphoblastic leukemia (ALL) is a leading cause of cancer-related death in children and adolescents, and cure rates for relapsed/refractory ALL remain dismal, highlighting the need for novel targeted therapies. To identify genome-wide metabolic-stress regulated genes, we used RNA-sequencing in ALL cells treated with AICAR, an AMPK activator. RNA-sequencing identified the immediate early genes (IEGs) as a subset of genes downregulated by AICAR. We show that AICAR-induced IEGs downregulation was blocked by an adenosine uptake inhibitor indicating AICAR was responsible for IEGs reprogramming. Using pharmacologic and genetic models we established this mechanism was AMPK-independent. Further investigations using kinase assays, PKD/PKC inhibitors and rescue experiments, demonstrated that AICAR directly inhibited PKD kinase activity and identified PKD as responsible for IEGs downregulation. Mechanistically, PKD inhibition suppressed phosphorylation and nuclear export of class IIa HDACs, which lowered histone H3 acetylation and decreased NFκB(p65) recruitment to IEGs promoters. Finally, PKD inhibition induced apoptosis via DUSP1/DUSP6 downregulation eliciting a DNA damage response. More importantly, ALL patient cells exhibited the same PKD-HDACs-IEGs-mediated mechanism. As proof of principle of the therapeutic potential of targeting PKD, we established the in vivo relevance of our findings using an NSG ALL mouse model. In conclusion, we identified a previously unreported PKD-dependent survival mechanism in response to AICAR-induced cellular stress in ALL through regulation of DUSPs and IEGs' expression. IMPLICATIONS: PKD mediates early transcriptional responses in ALL cells as an adaptive survival mechanism to overcome cellular stress.

Evaluation of repair activity by quantification of ribonucleotides in the genome.

Ribonucleotides incorporated in the genome are a source of endogenous DNA damage and also serve as signals for repair. Although recent advances of ribonucleotide detection by sequencing, the balance between incorporation and repair of ribonucleotides has not been elucidated. Here, we describe a competitive sequencing method, Ribonucleotide Scanning Quantification sequencing (RiSQ-seq), which enables absolute quantification of misincorporated ribonucleotides throughout the genome by background normalization and standard adjustment within a single sample. RiSQ-seq analysis of cells harboring wild-type DNA polymerases revealed that ribonucleotides were incorporated nonuniformly in the genome with a 3'-shifted distribution and preference for GC sequences. Although ribonucleotide profiles in wild-type and repair-deficient mutant strains showed a similar pattern, direct comparison of distinct ribonucleotide levels in the strains by RiSQ-seq enabled evaluation of ribonucleotide excision repair activity at base resolution and revealed the strand bias of repair. The distinct preferences of ribonucleotide incorporation and repair create vulnerable regions associated with indel hotspots, suggesting that repair at sites of ribonucleotide misincorporation serves to maintain genome integrity and that RiSQ-seq can provide an estimate of indel risk.

Ribonucleotide incorporation into DNA during DNA replication and its consequences.

Ribonucleotides are the most abundant non-canonical nucleotides in the genome. Their vast presence and influence over genome biology is becoming increasingly appreciated. Here we review the recent progress made in understanding their genomic presence, incorporation characteristics and usefulness as biomarkers for polymerase enzymology. We also discuss ribonucleotide processing, the genetic consequences of unrepaired ribonucleotides in DNA and evidence supporting the significance of their transient presence in the nuclear genome.

Tyrosyl-DNA phosphodiesterases are involved in mutagenic events at a ribonucleotide embedded into DNA in human cells.

Ribonucleoside triphosphates are often incorporated into genomic DNA during DNA replication. The accumulation of unrepaired ribonucleotides is associated with genomic instability, which is mediated by DNA topoisomerase 1 (Top1) processing of embedded ribonucleotides. The cleavage initiated by Top1 at the site of a ribonucleotide leads to the formation of a Top1-DNA cleavage complex (Top1cc), occasionally resulting in a DNA double-strand break (DSB). In humans, tyrosyl-DNA phosphodiesterases (TDPs) are essential repair enzymes that resolve the trapped Top1cc followed by downstream repair factors. However, there is limited cellular evidence of the involvement of TDPs in the processing of incorporated ribonucleotides in mammals. We assessed the role of TDPs in mutagenesis induced by a single ribonucleotide embedded into DNA. A supF shuttle vector site-specifically containing a single riboguanosine (rG) was introduced into the human lymphoblastoid TK6 cell line and its TDP1-, TDP2-, and TDP1/TDP2-deficient derivatives. TDP1 and TDP2 insufficiency remarkably decreased the mutant frequency caused by an embedded rG. The ratio of large deletion mutations induced by rG was also substantially lower in TDP1/TDP2-deficient cells than wild-type cells. Furthermore, the disruption of TDPs reduced the length of rG-mediated large deletion mutations. The recovery ratio of the propagated plasmid was also increased in TDP1/TDP2-deficient cells after the transfection of the shuttle vector containing rG. The results suggest that TDPs-mediated ribonucleotide processing cascade leads to unfavorable consequences, whereas in the absence of these repair factors, a more error-free processing pathway might function to suppress the ribonucleotide-induced mutagenesis. Furthermore, base substitution mutations at sites outside the position of rG were detected in the supF gene via a TDPs-independent mechanism. Overall, we provide new insights into the mechanism of mutagenesis induced by an embedded ribonucleotide in mammalian cells, which may lead to the fatal phenotype in the ribonucleotide excision repair deficiency.

Novel ribonucleotide discrimination in the RNA polymerase-like two-barrel catalytic core of Family D DNA polymerases.

Family D DNA polymerase (PolD) is the essential replicative DNA polymerase for duplication of most archaeal genomes. PolD contains a unique two-barrel catalytic core absent from all other DNA polymerase families but found in RNA polymerases (RNAPs). While PolD has an ancestral RNA polymerase catalytic core, its active site has evolved the ability to discriminate against ribonucleotides. Until now, the mechanism evolved by PolD to prevent ribonucleotide incorporation was unknown. In all other DNA polymerase families, an active site steric gate residue prevents ribonucleotide incorporation. In this work, we identify two consensus active site acidic (a) and basic (b) motifs shared across the entire two-barrel nucleotide polymerase superfamily, and a nucleotide selectivity (s) motif specific to PolD versus RNAPs. A novel steric gate histidine residue (H931 in Thermococcus sp. 9°N PolD) in the PolD s-motif both prevents ribonucleotide incorporation and promotes efficient dNTP incorporation. Further, a PolD H931A steric gate mutant abolishes ribonucleotide discrimination and readily incorporates a variety of 2' modified nucleotides. Taken together, we construct the first putative nucleotide bound PolD active site model and provide structural and functional evidence for the emergence of DNA replication through the evolution of an ancestral RNAP two-barrel catalytic core.

RNases H1 and H2: guardians of the stability of the nuclear genome when supply of dNTPs is limiting for DNA synthesis.

RNA/DNA hybrids are processed by RNases H1 and H2, while single ribonucleoside-monophosphates (rNMPs) embedded in genomic DNA are removed by the error-free, RNase H2-dependent ribonucleotide excision repair (RER) pathway. In the absence of RER, however, topoisomerase 1 (Top1) can cleave single genomic rNMPs in a mutagenic manner. In RNase H2-deficient mice, the accumulation of genomic rNMPs above a threshold of tolerance leads to catastrophic genomic instability that causes embryonic lethality. In humans, deficiencies in RNase H2 induce the autoimmune disorders Aicardi-Goutières syndrome and systemic lupus erythematosus, and cause skin and intestinal cancers. Recently, we reported that in Saccharomyces cerevisiae, the depletion of Rnr1, the major catalytic subunit of ribonucleotide reductase (RNR), which converts ribonucleotides to deoxyribonucleotides, leads to cell lethality in absence of RNases H1 and H2. We hypothesized that under replicative stress and compromised DNA repair that are elicited by an insufficient supply of deoxyribonucleoside-triphosphates (dNTPs), cells cannot survive the accumulation of persistent RNA/DNA hybrids. Remarkably, we found that cells lacking RNase H2 accumulate ~ 5-fold more genomic rNMPs in absence than in presence of Rnr1. When the load of genomic rNMPs is further increased in the presence of a replicative DNA polymerase variant that over-incorporates rNMPs in leading or lagging strand, cells missing both Rnr1 and RNase H2 suffer from severe growth defects. These are reversed in absence of Top1. Thus, in cells lacking RNase H2 and containing a limiting supply of dNTPs, there is a threshold of tolerance for the accumulation of genomic ribonucleotides that is tightly associated with Top1-mediated DNA damage. In this mini-review, we describe the implications of the loss of RNase H2, or RNases H1 and H2, on the integrity of the nuclear genome and viability of budding yeast cells that are challenged with a critically low supply of dNTPs. We further propose that our findings in budding yeast could pave the way for the study of the potential role of mammalian RNR in RNase H2-related diseases.

Validation strategies for antibodies targeting modified ribonucleotides.

Chemical modifications are found on almost all RNAs and affect their coding and noncoding functions. The identification of m6A on mRNA and its important role in gene regulation stimulated the field to investigate whether additional modifications are present on mRNAs. Indeed, modifications including m1A, m5C, m7G, 2'-OMe, and Ψ were detected. However, since their abundances are low and tools used for their corroboration are often not well characterized, their physiological relevance remains largely elusive. Antibodies targeting modified nucleotides are often used but have limitations such as low affinity or specificity. Moreover, they are not always well characterized and due to the low abundance of the modification, particularly on mRNAs, generated data sets might resemble noise rather than specific modification patterns. Therefore, it is critical that the affinity and specificity is rigorously tested using complementary approaches. Here, we provide an experimental toolbox that allows for testing antibody performance prior to their use.

Specificity of end resection pathways for double-strand break regions containing ribonucleotides and base lesions.

DNA double-strand break repair by homologous recombination begins with nucleolytic resection of the 5' DNA strand at the break ends. Long-range resection is catalyzed by EXO1 and BLM-DNA2, which likely have to navigate through ribonucleotides and damaged bases. Here, we show that a short stretch of ribonucleotides at the 5' terminus stimulates resection by EXO1. Ribonucleotides within a 5' flap are resistant to cleavage by DNA2, and extended RNA:DNA hybrids inhibit both strand separation by BLM and resection by EXO1. Moreover, 8-oxo-guanine impedes EXO1 but enhances resection by BLM-DNA2, and an apurinic/apyrimidinic site stimulates resection by BLM-DNA2 and DNA strand unwinding by BLM. Accordingly, depletion of OGG1 or APE1 leads to greater dependence of DNA resection on DNA2. Importantly, RNase H2A deficiency impairs resection overall, which we attribute to the accumulation of long RNA:DNA hybrids at DNA ends. Our results help explain why eukaryotic cells possess multiple resection nucleases.

The conformational landscape of transcription intermediates involved in the regulation of the ZMP-sensing riboswitch from Thermosinus carboxydivorans.

Recently, prokaryotic riboswitches have been identified that regulate transcription in response to change of the concentration of secondary messengers. The ZMP (5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR))-sensing riboswitch from Thermosinus carboxydivorans is a transcriptional ON-switch that is involved in purine and carbon-1 metabolic cycles. Its aptamer domain includes the pfl motif, which features a pseudoknot, impeding rho-independent terminator formation upon stabilization by ZMP interaction. We herein investigate the conformational landscape of transcriptional intermediates including the expression platform of this riboswitch and characterize the formation and unfolding of the important pseudoknot structure in the context of increasing length of RNA transcripts. NMR spectroscopic data show that even surprisingly short pre-terminator stems are able to disrupt ligand binding and thus metabolite sensing. We further show that the pseudoknot structure, a prerequisite for ligand binding, is preformed in transcription intermediates up to a certain length. Our results describe the conformational changes of 13 transcription intermediates of increasing length to delineate the change in structure as mRNA is elongated during transcription. We thus determine the length of the key transcription intermediate to which addition of a single nucleotide leads to a drastic drop in ZMP affinity.

Elimination of rNMPs from mitochondrial DNA has no effect on its stability.

Ribonucleotides (rNMPs) incorporated in the nuclear genome are a well-established threat to genome stability and can result in DNA strand breaks when not removed in a timely manner. However, the presence of a certain level of rNMPs is tolerated in mitochondrial DNA (mtDNA) although aberrant mtDNA rNMP content has been identified in disease models. We investigated the effect of incorporated rNMPs on mtDNA stability over the mouse life span and found that the mtDNA rNMP content increased during early life. The rNMP content of mtDNA varied greatly across different tissues and was defined by the rNTP/dNTP ratio of the tissue. Accordingly, mtDNA rNMPs were nearly absent in SAMHD1-/- mice that have increased dNTP pools. The near absence of rNMPs did not, however, appreciably affect mtDNA copy number or the levels of mtDNA molecules with deletions or strand breaks in aged animals near the end of their life span. The physiological rNMP load therefore does not contribute to the progressive loss of mtDNA quality that occurs as mice age.