Strand specificity of ribonucleotide excision repair in Escherichia coli.

In Escherichia coli, replication of both strands of genomic DNA is carried out by a single replicase-DNA polymerase III holoenzyme (pol III HE). However, in certain genetic backgrounds, the low-fidelity TLS polymerase, DNA polymerase V (pol V) gains access to undamaged genomic DNA where it promotes elevated levels of spontaneous mutagenesis preferentially on the lagging strand. We employed active site mutants of pol III (pol IIIα_S759N) and pol V (pol V_Y11A) to analyze ribonucleotide incorporation and removal from the E. coli chromosome on a genome-wide scale under conditions of normal replication, as well as SOS induction. Using a variety of methods tuned to the specific properties of these polymerases (analysis of lacI mutational spectra, lacZ reversion assay, HydEn-seq, alkaline gel electrophoresis), we present evidence that repair of ribonucleotides from both DNA strands in E. coli is unequal. While RNase HII plays a primary role in leading-strand Ribonucleotide Excision Repair (RER), the lagging strand is subject to other repair systems (RNase HI and under conditions of SOS activation also Nucleotide Excision Repair). Importantly, we suggest that RNase HI activity can also influence the repair of single ribonucleotides incorporated by the replicase pol III HE into the lagging strand.

5'-End Mapping in Human Mitochondrial DNA.

Here, we describe an assay that enables mapping of 5'-ends across the genome using next-generation sequencing on an Illumina platform, 5'-End-sequencing (5'-End-seq). We use this method to map free 5'-ends in mtDNA isolated from fibroblasts. This method can be used to answer key questions regarding DNA integrity, DNA replication mechanisms and to identify priming events, primer processing, nick processing, and double strand break processing on the entire genome.

Mammalian RNase H1 directs RNA primer formation for mtDNA replication initiation and is also necessary for mtDNA replication completion.

The in vivo role for RNase H1 in mammalian mitochondria has been much debated. Loss of RNase H1 is embryonic lethal and to further study its role in mtDNA expression we characterized a conditional knockout of Rnaseh1 in mouse heart. We report that RNase H1 is essential for processing of RNA primers to allow site-specific initiation of mtDNA replication. Without RNase H1, the RNA:DNA hybrids at the replication origins are not processed and mtDNA replication is initiated at non-canonical sites and becomes impaired. Importantly, RNase H1 is also needed for replication completion and in its absence linear deleted mtDNA molecules extending between the two origins of mtDNA replication are formed accompanied by mtDNA depletion. The steady-state levels of mitochondrial transcripts follow the levels of mtDNA, and RNA processing is not altered in the absence of RNase H1. Finally, we report the first patient with a homozygous pathogenic mutation in the hybrid-binding domain of RNase H1 causing impaired mtDNA replication. In contrast to catalytically inactive variants of RNase H1, this mutant version has enhanced enzyme activity but shows impaired primer formation. This finding shows that the RNase H1 activity must be strictly controlled to allow proper regulation of mtDNA replication.

Tracking Escherichia coli DNA polymerase V to the entire genome during the SOS response.

Ribonucleotides are frequently incorporated into DNA and can be used as a marker of DNA replication enzymology. To investigate on a genome-wide scale, how E. coli pol V accesses undamaged chromosomal DNA during the SOS response, we mapped the location of ribonucleotides incorporated by steric gate variants of pol V across the entire E. coli genome. To do so, we used strains that are deficient in ribonucleotide excision repair (ΔrnhB), deficient in pol IV DNA polymerase, constitutively express all SOS-regulated genes [lexA(Def)] and constitutively "activated" RecA* (recA730). The strains also harbor two steric gate variants of E. coli pol V (Y11A, or F10L), or a homolog of pol V, (pol VR391-Y13A). Ribonucleotides are frequently incorporated by the pol V-Y11A and pol VR391-Y13A variants, with a preference to the lagging strand. In contrast, the pol V-F10L variant incorporates less ribonucleotides and no strand preference is observed. Sharp transitions in strand specificity are observed at the replication origin (oriC), while a gradient is observed at the termination region. To activate RecA* in a recA+ strain, we treated the strains with ciprofloxacin and genome-wide mapped the location of the incorporated ribonucleotides. Again, the pol V-Y11A steric gate variant exhibited a lagging strand preference. Our data are consistent with a specific role for pol V in lagging strand DNA synthesis across the entire E. coli genome during the SOS response.

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.

RNase H1 directs origin-specific initiation of DNA replication in human mitochondria.

Human mitochondrial DNA (mtDNA) replication is first initiated at the origin of H-strand replication. The initiation depends on RNA primers generated by transcription from an upstream promoter (LSP). Here we reconstitute this process in vitro using purified transcription and replication factors. The majority of all transcription events from LSP are prematurely terminated after ~120 nucleotides, forming stable R-loops. These nascent R-loops cannot directly prime mtDNA synthesis, but must first be processed by RNase H1 to generate 3'-ends that can be used by DNA polymerase γ to initiate DNA synthesis. Our findings are consistent with recent studies of a knockout mouse model, which demonstrated that RNase H1 is required for R-loop processing and mtDNA maintenance in vivo. Both R-loop formation and DNA replication initiation are stimulated by the mitochondrial single-stranded DNA binding protein. In an RNase H1 deficient patient cell line, the precise initiation of mtDNA replication is lost and DNA synthesis is initiated from multiple sites throughout the mitochondrial control region. In combination with previously published in vivo data, the findings presented here suggest a model, in which R-loop processing by RNase H1 directs origin-specific initiation of DNA replication in human mitochondria.

DNA polymerase η contributes to genome-wide lagging strand synthesis.

DNA polymerase η (pol η) is best known for its ability to bypass UV-induced thymine-thymine (T-T) dimers and other bulky DNA lesions, but pol η also has other cellular roles. Here, we present evidence that pol η competes with DNA polymerases α and δ for the synthesis of the lagging strand genome-wide, where it also shows a preference for T-T in the DNA template. Moreover, we found that the C-terminus of pol η, which contains a PCNA-Interacting Protein motif is required for pol η to function in lagging strand synthesis. Finally, we provide evidence that a pol η dependent signature is also found to be lagging strand specific in patients with skin cancer. Taken together, these findings provide insight into the physiological role of DNA synthesis by pol η and have implications for our understanding of how our genome is replicated to avoid mutagenesis, genome instability and cancer.

Elevated pyrimidine dimer formation at distinct genomic bases underlies promoter mutation hotspots in UV-exposed cancers.

Sequencing of whole cancer genomes has revealed an abundance of recurrent mutations in gene-regulatory promoter regions, in particular in melanoma where strong mutation hotspots are observed adjacent to ETS-family transcription factor (TF) binding sites. While sometimes interpreted as functional driver events, these mutations are commonly believed to be due to locally inhibited DNA repair. Here, we first show that low-dose UV light induces mutations preferably at a known ETS promoter hotspot in cultured cells even in the absence of global or transcription-coupled nucleotide excision repair (NER). Further, by genome-wide mapping of cyclobutane pyrimidine dimers (CPDs) shortly after UV exposure and thus before DNA repair, we find that ETS-related mutation hotspots exhibit strong increases in CPD formation efficacy in a manner consistent with tumor mutation data at the single-base level. Analysis of a large whole genome cohort illustrates the widespread contribution of this effect to recurrent mutations in melanoma. While inhibited NER underlies a general increase in somatic mutation burden in regulatory elements including ETS sites, our data supports that elevated DNA damage formation at specific genomic bases is at the core of the prominent promoter mutation hotspots seen in skin cancers, thus explaining a key phenomenon in whole-genome cancer analyses.

Mapping Ribonucleotides Incorporated into DNA by Hydrolytic End-Sequencing.

Ribonucleotides embedded within DNA render the DNA sensitive to the formation of single-stranded breaks under alkali conditions. Here, we describe a next-generation sequencing method called hydrolytic end sequencing (HydEn-seq) to map ribonucleotides inserted into the genome of Saccharomyce cerevisiae strains deficient in ribonucleotide excision repair. We use this method to map several genomic features in wild-type and replicase variant yeast strains.

Simultaneous Mapping and Quantitation of Ribonucleotides in Human Mitochondrial DNA.

Established approaches to estimate the number of ribonucleotides present in a genome are limited to the quantitation of incorporated ribonucleotides using short synthetic DNA fragments or plasmids as templates and then extrapolating the results to the whole genome. Alternatively, the number of ribonucleotides present in a genome may be estimated using alkaline gels or Southern blots. More recent in vivo approaches employ Next-generation sequencing allowing genome-wide mapping of ribonucleotides, providing the position and identity of embedded ribonucleotides. However, they do not allow quantitation of the number of ribonucleotides which are incorporated into a genome. Here we describe how to simultaneously map and quantitate the number of ribonucleotides which are incorporated into human mitochondrial DNA in vivo by Next-generation sequencing. We use highly intact DNA and introduce sequence specific double strand breaks by digesting it with an endonuclease, subsequently hydrolyzing incorporated ribonucleotides with alkali. The generated ends are ligated with adapters and these ends are sequenced on a Next-generation sequencing machine. The absolute number of ribonucleotides can be calculated as the number of reads outside the recognition site per average number of reads at the recognition site for the sequence specific endonuclease. This protocol may also be utilized to map and quantitate free nicks in DNA and allows adaption to map other DNA lesions that can be processed to 5´-OH ends or 5´-phosphate ends. Furthermore, this method can be applied to any organism, given that a suitable reference genome is available. This protocol therefore provides an important tool to study DNA replication, 5´-end processing, DNA damage, and DNA repair.

Ribonucleotides incorporated by the yeast mitochondrial DNA polymerase are not repaired.

Incorporation of ribonucleotides into DNA during genome replication is a significant source of genomic instability. The frequency of ribonucleotides in DNA is determined by deoxyribonucleoside triphosphate/ribonucleoside triphosphate (dNTP/rNTP) ratios, by the ability of DNA polymerases to discriminate against ribonucleotides, and by the capacity of repair mechanisms to remove incorporated ribonucleotides. To simultaneously compare how the nuclear and mitochondrial genomes incorporate and remove ribonucleotides, we challenged these processes by changing the balance of cellular dNTPs. Using a collection of yeast strains with altered dNTP pools, we discovered an inverse relationship between the concentration of individual dNTPs and the amount of the corresponding ribonucleotides incorporated in mitochondrial DNA, while in nuclear DNA the ribonucleotide pattern was only altered in the absence of ribonucleotide excision repair. Our analysis uncovers major differences in ribonucleotide repair between the two genomes and provides concrete evidence that yeast mitochondria lack mechanisms for removal of ribonucleotides incorporated by the mtDNA polymerase. Furthermore, as cytosolic dNTP pool imbalances were transmitted equally well into the nucleus and the mitochondria, our results support a view of the cytosolic and mitochondrial dNTP pools in frequent exchange.

Nucleotide pools dictate the identity and frequency of ribonucleotide incorporation in mitochondrial DNA.

Previous work has demonstrated the presence of ribonucleotides in human mitochondrial DNA (mtDNA) and in the present study we use a genome-wide approach to precisely map the location of these. We find that ribonucleotides are distributed evenly between the heavy- and light-strand of mtDNA. The relative levels of incorporated ribonucleotides reflect that DNA polymerase γ discriminates the four ribonucleotides differentially during DNA synthesis. The observed pattern is also dependent on the mitochondrial deoxyribonucleotide (dNTP) pools and disease-causing mutations that change these pools alter both the absolute and relative levels of incorporated ribonucleotides. Our analyses strongly suggest that DNA polymerase γ-dependent incorporation is the main source of ribonucleotides in mtDNA and argues against the existence of a mitochondrial ribonucleotide excision repair pathway in human cells. Furthermore, we clearly demonstrate that when dNTP pools are limiting, ribonucleotides serve as a source of building blocks to maintain DNA replication. Increased levels of embedded ribonucleotides in patient cells with disturbed nucleotide pools may contribute to a pathogenic mechanism that affects mtDNA stability and impair new rounds of mtDNA replication.

Gene duplications and losses among vertebrate deoxyribonucleoside kinases of the non-TK1 Family.

Deoxyribonucleoside kinases (dNKs) salvage deoxyribonucleosides (dNs) and catalyze the rate limiting step of this salvage pathway by converting dNs into corresponding monophosphate forms. These enzymes serve as an excellent model to study duplicated genes and their evolutionary history. So far, among vertebrates only four mammalian dNKs have been studied for their substrate specificity and kinetic properties. However, some vertebrates, such as fish, frogs, and birds, apparently possess a duplicated homolog of deoxycytidine kinase (dCK). In this study, we characterized a family of dCK/deoxyguanosine kinase (dGK)-like enzymes from a frog Xenopus laevis and a bird Gallus gallus. We showed that X. laevis has a duplicated dCK gene and a dGK gene, whereas G. gallus has a duplicated dCK gene but has lost the dGK gene. We cloned, expressed, purified, and subsequently determined the kinetic parameters of the dCK/dGK enzymes encoded by these genes. The two dCK enzymes in G. gallus have broader substrate specificity than their human or X. laevis counterparts. Additionally, the duplicated dCK enzyme in G. gallus might have become mitochondria. Based on our study we postulate that changing and adapting substrate specificities and subcellular localization are likely the drivers behind the evolution of vertebrate dNKs.

Measuring ribonucleotide incorporation into DNA in vitro and in vivo.

Ribonucleotides are incorporated into genomes by DNA polymerases, they can be removed, and if not removed, they can have deleterious and beneficial consequences. Here, we describe an assay to quantify stable ribonucleotide incorporation by DNA polymerases in vitro, and an assay to probe for ribonucleotides in each of the two DNA strands of the yeast nuclear genome.

Evidence that processing of ribonucleotides in DNA by topoisomerase 1 is leading-strand specific.

Ribonucleotides incorporated during DNA replication are removed by RNase H2-dependent ribonucleotide excision repair (RER). In RER-defective yeast, topoisomerase 1 (Top1) incises DNA at unrepaired ribonucleotides, initiating their removal, but this is accompanied by RNA-DNA-damage phenotypes. Here we show that these phenotypes are incurred by a high level of ribonucleotides incorporated by a leading strand-replicase variant, DNA polymerase (Pol) ɛ, but not by orthologous variants of the lagging-strand replicases, Pols α or δ. Moreover, loss of both RNases H1 and H2 is lethal in combination with increased ribonucleotide incorporation by Pol ɛ but not by Pols α or δ. Several explanations for this asymmetry are considered, including the idea that Top1 incision at ribonucleotides relieves torsional stress in the nascent leading strand but not in the nascent lagging strand, in which preexisting nicks prevent the accumulation of superhelical tension.

Tracking replication enzymology in vivo by genome-wide mapping of ribonucleotide incorporation.

Ribonucleotides are frequently incorporated into DNA during replication in eukaryotes. Here we map genome-wide distribution of these ribonucleotides as markers of replication enzymology in budding yeast, using a new 5' DNA end-mapping method, hydrolytic end sequencing (HydEn-seq). HydEn-seq of DNA from ribonucleotide excision repair-deficient strains reveals replicase- and strand-specific patterns of ribonucleotides in the nuclear genome. These patterns support the roles of DNA polymerases α and δ in lagging-strand replication and of DNA polymerase ɛ in leading-strand replication. They identify replication origins, termination zones and variations in ribonucleotide incorporation frequency across the genome that exceed three orders of magnitude. HydEn-seq also reveals strand-specific 5' DNA ends at mitochondrial replication origins, thus suggesting unidirectional replication of a circular genome. Given the conservation of enzymes that incorporate and process ribonucleotides in DNA, HydEn-seq can be used to track replication enzymology in other organisms.

Heterogeneous polymerase fidelity and mismatch repair bias genome variation and composition.

Mutational heterogeneity must be taken into account when reconstructing evolutionary histories, calibrating molecular clocks, and predicting links between genes and disease. Selective pressures and various DNA transactions have been invoked to explain the heterogeneous distribution of genetic variation between species, within populations, and in tissue-specific tumors. To examine relationships between such heterogeneity and variations in leading- and lagging-strand replication fidelity and mismatch repair, we accumulated 40,000 spontaneous mutations in eight diploid yeast strains in the absence of selective pressure. We found that replicase error rates vary by fork direction, coding state, nucleosome proximity, and sequence context. Further, error rates and DNA mismatch repair efficiency both vary by mismatch type, responsible polymerase, replication time, and replication origin proximity. Mutation patterns implicate replication infidelity as one driver of variation in somatic and germline evolution, suggest mechanisms of mutual modulation of genome stability and composition, and predict future observations in specific cancers.

Plants salvage deoxyribonucleosides in mitochondria.

Deoxyribonucleoside kinases phosphorylate deoxyribonucleosides into the corresponding 5'-monophosphate deoxyribonucleosides to supply the cell with nucleic acid precursors. In mitochondrial fractions of the model plant Arabidopsis thaliana, we detected deoxyadenosine and thymidine kinase activities, while the cytosol fraction contained six-fold lower activity and chloroplasts contained no measurable activities. In addition, a mitochondrial fraction isolated from the potato Solanum tuberosum contained thymidine kinase and deoxyadenosine kinase activities. We conclude that an active salvage of deoxyribonucleosides in plants takes place in their mitochondria. In general, the observed localization of the plant dNK activities in the mitochondrion suggests that plants have a different organization of the deoxyribonucleoside salvage compared to mammals.

Structure-function analysis of ribonucleotide bypass by B family DNA replicases.

Ribonucleotides are frequently incorporated into DNA during replication, they are normally removed, and failure to remove them results in replication stress. This stress correlates with DNA polymerase (Pol) stalling during bypass of ribonucleotides in DNA templates. Here we demonstrate that stalling by yeast replicative Pols δ and ε increases as the number of consecutive template ribonucleotides increases from one to four. The homologous bacteriophage RB69 Pol also stalls during ribonucleotide bypass, with a pattern most similar to that of Pol ε. Crystal structures of an exonuclease-deficient variant of RB69 Pol corresponding to multiple steps in single ribonucleotide bypass reveal that increased stalling is associated with displacement of Tyr391 and an unpreferred C2'-endo conformation for the ribose. Even less efficient bypass of two consecutive ribonucleotides in DNA correlates with similar movements of Tyr391 and displacement of one of the ribonucleotides along with the primer-strand DNA backbone. These structure-function studies have implications for cellular signaling by ribonucleotides, and they may be relevant to replication stress in cells defective in ribonucleotide excision repair, including humans suffering from autoimmune disease associated with RNase H2 defects.

Ribonucleotides are signals for mismatch repair of leading-strand replication errors.

To maintain genome stability, mismatch repair of nuclear DNA replication errors must be directed to the nascent strand, likely by DNA ends and PCNA. Here we show that the efficiency of mismatch repair in Saccharomyces cerevisiae is reduced by inactivating RNase H2, which nicks DNA containing ribonucleotides incorporated during replication. In strains encoding mutator polymerases, this reduction is preferential for repair of mismatches made by leading-strand DNA polymerase ε as compared to lagging-strand DNA polymerase δ. The results suggest that RNase-H2-dependent processing of ribonucleotides transiently present in DNA after replication may direct mismatch repair to the continuously replicated nascent leading strand.