Lagging-strand replication shapes the mutational landscape of the genome

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Bulk synthesis and beyond: The roles of eukaryotic replicative DNA polymerases

An organism's genomic DNA must be accurately duplicated during each cell cycle. DNA synthesis is catalysed by DNA polymerase enzymes, which extend nucleotide polymers in a 5' to 3' direction. This inherent directionality necessitates that one strand is synthesised forwards (leading), while the other is synthesised backwards discontinuously (lagging) to couple synthesis to the unwinding of duplex DNA. Eukaryotic cells possess many diverse polymerases that coordinate to replicate DNA, with the three main replicative polymerases being Pol α, Pol δ and Pol ε. Studies conducted in yeasts and human cells utilising mutant polymerases that incorporate molecular signatures into nascent DNA implicate Pol ε in leading strand synthesis and Pol α and Pol δ in lagging strand replication. Recent structural insights have revealed how the spatial organization of these enzymes around the core helicase facilitates their strand-specific roles. However, various challenging situations during replication require flexibility in the usage of these enzymes, such as during replication initiation or encounters with replication-blocking adducts. This review summarises the roles of the replicative polymerases in bulk DNA replication and explores their flexible and dynamic deployment to complete genome replication. We also examine how polymerase usage patterns can inform our understanding of global replication dynamics by revealing replication fork directionality to identify regions of replication initiation and termination.

Unraveling the complexity of asymmetric DNA replication: Advancements in ribonucleotide mapping techniques and beyond

DNA replication is a fundamental process for cell proliferation, governed by intricate mechanisms involving leading and lagging strand synthesis. In eukaryotes, canonical DNA replication occurs during the S phase of the cell cycle, facilitated by various components of the replicative machinery at sites known as replication origins. Leading and lagging strands exhibit distinct replication dynamics, with leading strand replication being relatively straightforward compared to the complex synthesis of lagging strands involving Okazaki fragment maturation. Central to DNA synthesis are DNA polymerases, with Polα, Polε, and Polδ playing pivotal roles, each specializing in specific tasks during replication. Notably, leading and lagging strands are replicated by different polymerases, contributing to the division of labor in DNA replication. Understanding the enzymology of asymmetric DNA replication has been challenging, with methods relying on ribonucleotide incorporation and next-generation sequencing techniques offering comprehensive insights. These methodologies, such as HydEn-seq, PU-seq, ribose-seq, and emRiboSeq, offer insights into polymerase activity and strand synthesis, aiding in understanding DNA replication dynamics. Recent advancements include novel conditional mutants for ribonucleotide excision repair, enzymatic cleavage alternatives, and unified pipelines for data analysis. Further developments in adapting techniques to different organisms, studying non-canonical polymerases, and exploring new sequencing platforms hold promise for expanding our understanding of DNA replication dynamics. Integrating strand-specific information into single-cell studies could offer novel insights into enzymology, opening avenues for future research and applications in repair and replication biology.

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

Adaptive use of error-prone DNA polymerases provides flexibility in genome replication during tumorigenesis

Human cells possess many different polymerase enzymes, which collaborate in conducting DNA replication and genome maintenance to ensure faithful duplication of genetic material. Each polymerase performs a specialized role, together providing a balance of accuracy and flexibility to the replication process. Perturbed replication increases the requirement for flexibility to ensure duplication of the entire genome. Flexibility is provided via the use of error-prone polymerases, which maintain the progression of challenged DNA replication at the expense of mutagenesis, an enabling characteristic of cancer. This review describes our recent understanding of mechanisms that alter the usage of polymerases during tumorigenesis and examines the implications of this for cell survival and tumor progression. Although expression levels of polymerases are often misregulated in cancers, this does not necessarily alter polymerase usage since an additional regulatory step may govern the use of these enzymes. We therefore also examine how the regulatory mechanisms of DNA polymerases, such as Rad18-mediated PCNA ubiquitylation, may impact the functionalization of error-prone polymerases to tolerate oncogene-induced replication stress. Crucially, it is becoming increasingly evident that cancer cells utilize error-prone polymerases to sustain ongoing replication in response to oncogenic mutations which inactivate key DNA replication and repair pathways, such as BRCA deficiency. This accelerates mutagenesis and confers chemoresistance, but also presents a dependency that can potentially be exploited by therapeutics.

Monitoring and quantifying replication fork dynamics with high-throughput methods

Before each cell division, eukaryotic cells must replicate their chromosomes to ensure the accurate transmission of genetic information. Chromosome replication involves more than just DNA duplication; it also includes chromatin assembly, inheritance of epigenetic marks, and faithful resumption of all genomic functions after replication. Recent progress in quantitative technologies has revolutionized our understanding of the complexity and dynamics of DNA replication forks at both molecular and genomic scales. Here, we highlight the pivotal role of these novel methods in uncovering the principles and mechanisms of chromosome replication. These technologies have illuminated the regulation of genome replication programs, quantified the impact of DNA replication on genomic mutations and evolutionary processes, and elucidated the mechanisms of replication-coupled chromatin assembly and epigenome maintenance.

Strand-resolved mutagenicity of DNA damage and repair

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

S-phase checkpoint prevents leading strand degradation from strand-associated nicks at stalled replication forks

The S-phase checkpoint is involved in coupling DNA unwinding with nascent strand synthesis and is critical to maintain replication fork stability in conditions of replicative stress. However, its role in the specific regulation of leading and lagging strands at stalled forks is unclear. By conditionally depleting RNaseH2 and analyzing polymerase usage genome-wide, we examine the enzymology of DNA replication during a single S-phase in the presence of replicative stress and show that there is a differential regulation of lagging and leading strands. In checkpoint proficient cells, lagging strand replication is down-regulated through an Elg1-dependent mechanism. Nevertheless, when checkpoint function is impaired we observe a defect specifically at the leading strand, which was partially dependent on Exo1 activity. Further, our genome-wide mapping of DNA single-strand breaks reveals that strand discontinuities highly accumulate at the leading strand in HU-treated cells, whose dynamics are affected by checkpoint function and Exo1 activity. Our data reveal an unexpected role of Exo1 at the leading strand and support a model of fork stabilization through prevention of unrestrained Exo1-dependent resection of leading strand-associated nicks after fork stalling.

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.

Human primosome requires replication protein A when copying DNA with inverted repeats

The human primosome, a four-subunit complex of primase and DNA polymerase alpha (Polα), initiates DNA synthesis on both chromosome strands by generating chimeric RNA-DNA primers for loading DNA polymerases delta and epsilon (Polε). Replication protein A (RPA) tightly binds to single-stranded DNA strands, protecting them from nucleolytic digestion and unauthorized transactions. We report here that RPA plays a critical role for the human primosome during DNA synthesis across inverted repeats prone to hairpin formation. On other alternatively structured DNA forming a G-quadruplex, RPA provides no assistance for primosome. A stimulatory effect of RPA on DNA synthesis across hairpins was also observed for the catalytic domain of Polα but not of Polε. The important factors for an efficient hairpin bypass by primosome are the high affinity of RPA to DNA based on four DNA-binding domains and the interaction of the winged-helix-turn-helix domain of RPA with Polα. Binding studies indicate that this interaction stabilizes the RPA/Polα complex on the primed template. This work provides insight into a cooperative action of RPA and primosome on DNA, which is critical for DNA synthesis across inverted repeats.

Primase promotes the competition between transcription and replication on the same template strand resulting in DNA damage

Transcription-replication conflicts (TRCs), especially Head-On TRCs (HO-TRCs) can introduce R-loops and DNA damage, however, the underlying mechanisms are still largely unclear. We previously identified a chloroplast-localized RNase H1 protein AtRNH1C that can remove R-loops and relax HO-TRCs for genome integrity. Through the mutagenesis screen, we identify a mutation in chloroplast-localized primase ATH that weakens the binding affinity of DNA template and reduces the activities of RNA primer synthesis and delivery. This slows down DNA replication, and reduces competition of transcription-replication, thus rescuing the developmental defects of atrnh1c. Strand-specific DNA damage sequencing reveals that HO-TRCs cause DNA damage at the end of the transcription unit in the lagging strand and overexpression of ATH can boost HO-TRCs and exacerbates DNA damage. Furthermore, mutation of plastid DNA polymerase Pol1A can similarly rescue the defects in atrnh1c mutants. Taken together these results illustrate a potentially conserved mechanism among organisms, of which the primase activity can promote the occurrence of transcription-replication conflicts leading to HO-TRCs and genome instability.

A mutation rate model at the basepair resolution identifies the mutagenic effect of polymerase III transcription

De novo mutations occur at substantially different rates depending on genomic location, sequence context and DNA strand. The success of methods to estimate selection intensity, infer demographic history and map rare disease genes, depends strongly on assumptions about the local mutation rate. Here we present Roulette, a genome-wide mutation rate model at basepair resolution that incorporates known determinants of local mutation rate. Roulette is shown to be more accurate than existing models. We use Roulette to refine the estimates of population growth within Europe by incorporating the full range of human mutation rates. The analysis of significant deviations from the model predictions revealed a tenfold increase in mutation rate in nearly all genes transcribed by polymerase III (Pol III), suggesting a new mutagenic mechanism. We also detected an elevated mutation rate within transcription factor binding sites restricted to sites actively used in testis and residing in promoters.

A sensor complements the steric gate when DNA polymerase ϵ discriminates ribonucleotides

The cellular imbalance between high concentrations of ribonucleotides (NTPs) and low concentrations of deoxyribonucleotides (dNTPs), is challenging for DNA polymerases when building DNA from dNTPs. It is currently believed that DNA polymerases discriminate against NTPs through a steric gate model involving a clash between a tyrosine and the 2'-hydroxyl of the ribonucleotide in the polymerase active site in B-family DNA polymerases. With the help of crystal structures of a B-family polymerase with a UTP or CTP in the active site, molecular dynamics simulations, biochemical assays and yeast genetics, we have identified a mechanism by which the finger domain of the polymerase sense NTPs in the polymerase active site. In contrast to the previously proposed polar filter, our experiments suggest that the amino acid residue in the finger domain senses ribonucleotides by steric hindrance. Furthermore, our results demonstrate that the steric gate in the palm domain and the sensor in the finger domain are both important when discriminating NTPs. Structural comparisons reveal that the sensor residue is conserved among B-family polymerases and we hypothesize that a sensor in the finger domain should be considered in all types of DNA polymerases.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Impairment of the non-catalytic subunit Dpb2 of DNA Pol ɛ results in increased involvement of Pol δ on the leading strand

The generally accepted model assumes that leading strand synthesis is performed by Pol ε, while lagging-strand synthesis is catalyzed by Pol δ. Pol ε has been shown to target the leading strand by interacting with the CMG helicase [Cdc45 Mcm2-7 GINS(Psf1-3, Sld5)]. Proper functioning of the CMG-Pol ɛ, the helicase-polymerase complex is essential for its progression and the fidelity of DNA replication. Dpb2p, the essential non-catalytic subunit of Pol ε plays a key role in maintaining the correct architecture of the replisome by acting as a link between Pol ε and the CMG complex. Using a temperature-sensitive dpb2-100 mutant previously isolated in our laboratory, and a genetic system which takes advantage of a distinct mutational signature of the Pol δ-L612M variant which allows detection of the involvement of Pol δ in the replication of particular DNA strands we show that in yeast cells with an impaired Dpb2 subunit, the contribution of Pol δ to the replication of the leading strand is significantly increased.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

Generalised interrelations among mutation rates drive the genomic compliance of Chargaff's second parity rule

Chargaff's second parity rule (PR-2), where the complementary base and k-mer contents are matching within the same strand of a double stranded DNA (dsDNA), is a phenomenon that invited many explanations. The strict compliance of nearly all nuclear dsDNA to PR-2 implies that the explanation should also be similarly adamant. In this work, we revisited the possibility of mutation rates driving PR-2 compliance. Starting from the assumption-free approach, we constructed kinetic equations for unconstrained simulations. The results were analysed for their PR-2 compliance by employing symbolic regression and machine learning techniques. We arrived to a generalised set of mutation rate interrelations in place in most species that allow for their full PR-2 compliance. Importantly, our constraints explain PR-2 in genomes out of the scope of the prior explanations based on the equilibration under mutation rates with simpler no-strand-bias constraints. We thus reinstate the role of mutation rates in PR-2 through its molecular core, now shown, under our formulation, to be tolerant to previously noted strand biases and incomplete compositional equilibration. We further investigate the time for any genome to reach PR-2, showing that it is generally earlier than the compositional equilibrium, and well within the age of life on Earth.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

Structures of human primosome elongation complexes

The synthesis of RNA-DNA primer by primosome requires coordination between primase and DNA polymerase α subunits, which is accompanied by unknown architectural rearrangements of multiple domains. Using cryogenic electron microscopy, we solved a 3.6 Å human primosome structure caught at an early stage of RNA primer elongation with deoxynucleotides. The structure confirms a long-standing role of primase large subunit and reveals new insights into how primosome is limited to synthesizing short RNA-DNA primers.

DNA Polymerase Delta Exhibits Altered Catalytic Properties on Lysine Acetylation

DNA polymerase delta is the primary polymerase that is involved in undamaged nuclear lagging strand DNA replication. Our mass-spectroscopic analysis has revealed that the human DNA polymerase δ is acetylated on subunits p125, p68, and p12. Using substrates that simulate Okazaki fragment intermediates, we studied alterations in the catalytic properties of acetylated polymerase and compared it to the unmodified form. The current data show that the acetylated form of human pol δ displays a higher polymerization activity compared to the unmodified form of the enzyme. Additionally, acetylation enhances the ability of the polymerase to resolve complex structures such as G-quadruplexes and other secondary structures that might be present on the template strand. More importantly, the ability of pol δ to displace a downstream DNA fragment is enhanced upon acetylation. Our current results suggest that acetylation has a profound effect on the activity of pol δ and supports the hypothesis that acetylation may promote higher-fidelity DNA replication.

Okazaki fragment maturation: DNA flap dynamics for cell proliferation and survival

Unsuccessful processing of Okazaki fragments leads to the accumulation of DNA breaks which are associated with many human diseases including cancer and neurodegenerative disorders. Recently, Okazaki fragment maturation (OFM) has received renewed attention regarding how unprocessed Okazaki fragments are sensed and repaired, and how inappropriate OFM impacts on genome stability and cell viability, especially in cancer cells. We provide an overview of the highly efficient and faithful canonical OFM pathways and their regulation of genomic integrity and cell survival. We also discuss how cells induce alternative error-prone OFM processes to promote cell survival in response to environmental stresses. Such stress-induced OFM processes may be important mechanisms driving mutagenesis, cellular evolution, and resistance to radio/chemotherapy and targeted therapeutics in human cancers.

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.

Determination of the Ribonucleotide Content of mtDNA Using Alkaline Gels

Impaired mitochondrial DNA (mtDNA) maintenance, due to, e.g., defects in the replication machinery or an insufficient dNTP supply, underlies a number of mitochondrial disorders. The normal process of mtDNA replication leads to the incorporation of multiple single ribonucleotides (rNMPs) per mtDNA molecule. Given that embedded rNMPs alter the stability and properties of the DNA, they may have consequences for mtDNA maintenance and thereby for mitochondrial disease. They also serve as a readout of the intramitochondrial NTP/dNTP ratios. In this chapter, we describe a method for the determination of mtDNA rNMP content using alkaline gel electrophoresis and Southern blotting. This procedure is suited for the analysis of mtDNA in total genomic DNA preparations as well as in purified form. Moreover, it can be performed using equipment found in most biomedical laboratories, allows the simultaneous analysis of 10-20 samples depending on the gel system employed, and can be modified for the analysis of other mtDNA modifications.

Population Genetic Considerations Regarding Evidence for Biased Mutation Rates in Arabidopsis thaliana

It has recently been proposed that lower mutation rates in gene bodies compared with upstream and downstream sequences in Arabidopsis thaliana are the result of an "adaptive" modification of the rate of beneficial and deleterious mutations in these functional regions. This claim was based both on analyses of mutation accumulation lines and on population genomics data. Here, we show that several questionable assumptions were used in the population genomics analyses. In particular, we demonstrate that the difference between gene bodies and less selectively constrained sequences in the magnitude of Tajima's D can in principle be explained by the presence of sites subject to purifying selection and does not require lower mutation rates in regions experiencing selective constraints.

Physicochemical models of protein-DNA binding with standard and modified base pairs

DNA-binding proteins play important roles in various cellular processes, but the mechanisms by which proteins recognize genomic target sites remain incompletely understood. Functional groups at the edges of the base pairs (bp) exposed in the DNA grooves represent physicochemical signatures. As these signatures enable proteins to form specific contacts between protein residues and bp, their study can provide mechanistic insights into protein-DNA binding. Existing experimental methods, such as X-ray crystallography, can reveal such mechanisms based on physicochemical interactions between proteins and their DNA target sites. However, the low throughput of structural biology methods limits mechanistic insights for selection of many genomic sites. High-throughput binding assays enable prediction of potential target sites by determining relative binding affinities of a protein to massive numbers of DNA sequences. Many currently available computational methods are based on the sequence of standard Watson-Crick bp. They assume that the contribution of overall binding affinity is independent for each base pair, or alternatively include dinucleotides or short k-mers. These methods cannot directly expand to physicochemical contacts, and they are not suitable to apply to DNA modifications or non-Watson-Crick bp. These variations include DNA methylation, and synthetic or mismatched bp. The proposed method, DeepRec, can predict relative binding affinities as function of physicochemical signatures and the effect of DNA methylation or other chemical modifications on binding. Sequence-based modeling methods are in comparison a coarse-grain description and cannot achieve such insights. Our chemistry-based modeling framework provides a path towards understanding genome function at a mechanistic level.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Pol32, an accessory subunit of DNA polymerase delta, plays an essential role in genome stability and pathogenesis of Candida albicans

Candida albicans is a pathobiont that inflicts serious bloodstream fungal infections in individuals with compromised immunity and gut dysbiosis. Genomic diversity in the form of copy number alteration, ploidy variation, and loss of heterozygosity as an adaptive mechanism to adverse environments is frequently observed in C. albicans. Such genomic variations also confer a varied degree of fungal virulence and drug resistance, yet the factors propelling these are not completely understood. DNA polymerase delta (Polδ) is an essential replicative DNA polymerase in the eukaryotic cell and is yet to be characterized in C. albicans. Therefore, this study was designed to gain insights into the role of Polδ, especially its non-essential subunit Pol32, in the genome plasticity and life cycle of C. albicans. PCNA, the DNA clamp, recruits Polδ to the replication fork for processive DNA replication. Unlike in Saccharomyces cerevisiae, the PCNA interaction protein (PIP) motif of CaPol32 is critical for Polδ's activity during DNA replication. Our comparative genetic analyses and whole-genome sequencing of POL32 proficient and deficient C. albicans cells revealed a critical role of Pol32 in DNA replication, cell cycle progression, and genome stability as SNPs, indels, and repeat variations were largely accumulated in pol32 null strain. The loss of pol32 in C. albicans conferred cell wall deformity; Hsp90 mediated azoles resistance, biofilm development, and a complete attenuation of virulence in an animal model of systemic candidiasis. Thus, although Pol32 is dispensable for cell survival, its function is essential for C. albicans pathogenesis; and we discuss its translational implications in antifungal drugs and whole-cell vaccine development.

Human DNA polymerase α has a strong mutagenic potential at the initial steps of DNA synthesis

DNA polymerase α (Polα) is essential for DNA replication initiation and makes a notable contribution to genome mutagenesis. The activity and fidelity of Polα during the early steps of DNA replication have not been well studied. Here we show that at the beginning of DNA synthesis, when extending the RNA primer received from primase, Polα is more mutagenic than during the later DNA elongation steps. Kinetic and binding studies revealed substantially higher activity and affinity to the template:primer when Polα interacts with ribonucleotides of a chimeric RNA-DNA primer. Polα activity greatly varies during first six steps of DNA synthesis, and the bias in the rates of correct and incorrect dNTP incorporation leads to impaired fidelity, especially upon the second step of RNA primer extension. Furthermore, increased activity and stability of Polα/template:primer complexes containing RNA-DNA primers result in higher efficiency of mismatch extension.

Global landscape of replicative DNA polymerase usage in the human genome

The division of labour among DNA polymerase underlies the accuracy and efficiency of replication. However, the roles of replicative polymerases have not been directly established in human cells. We developed polymerase usage sequencing (Pu-seq) in HCT116 cells and mapped Polε and Polα usage genome wide. The polymerase usage profiles show Polε synthesises the leading strand and Polα contributes mainly to lagging strand synthesis. Combining the Polε and Polα profiles, we accurately predict the genome-wide pattern of fork directionality plus zones of replication initiation and termination. We confirm that transcriptional activity contributes to the pattern of initiation and termination and, by separately analysing the effect of transcription on co-directional and converging forks, demonstrate that coupled DNA synthesis of leading and lagging strands is compromised by transcription in both co-directional and convergent forks. Polymerase uncoupling is particularly evident in the vicinity of large genes, including the two most unstable common fragile sites, FRA3B and FRA3D, thus linking transcription-induced polymerase uncoupling to chromosomal instability. Together, our result demonstrated that Pu-seq in human cells provides a powerful and straightforward methodology to explore DNA polymerase usage and replication fork dynamics.

Mechanistic investigation of human maturation of Okazaki fragments reveals slow kinetics

The final steps of lagging strand synthesis induce maturation of Okazaki fragments via removal of the RNA primers and ligation. Iterative cycles between Polymerase δ (Polδ) and Flap endonuclease-1 (FEN1) remove the primer, with an intermediary nick structure generated for each cycle. Here, we show that human Polδ is inefficient in releasing the nick product from FEN1, resulting in non-processive and remarkably slow RNA removal. Ligase 1 (Lig1) can release the nick from FEN1 and actively drive the reaction toward ligation. These mechanisms are coordinated by PCNA, which encircles DNA, and dynamically recruits Polδ, FEN1, and Lig1 to compete for their substrates. Our findings call for investigating additional pathways that may accelerate RNA removal in human cells, such as RNA pre-removal by RNase Hs, which, as demonstrated herein, enhances the maturation rate ~10-fold. They also suggest that FEN1 may attenuate the various activities of Polδ during DNA repair and recombination.

DNA replication initiation shapes the mutational landscape and expression of the human genome

The interplay between active biological processes and DNA repair is central to mutagenesis. Here, we show that the ubiquitous process of replication initiation is mutagenic, leaving a specific mutational footprint at thousands of early and efficient replication origins. The observed mutational pattern is consistent with two distinct mechanisms, reflecting the two-step process of origin activation, triggering the formation of DNA breaks at the center of origins and local error-prone DNA synthesis in their immediate vicinity. We demonstrate that these replication initiation-dependent mutational processes exert an influence on phenotypic diversity in humans that is disproportionate to the origins' genomic size: By increasing mutational loads at gene promoters and splice junctions, the presence of an origin significantly influences both gene expression and mRNA isoform usage. Last, we show that mutagenesis at origins not only drives the evolution of origin sequences but also contributes to sculpting regulatory domains of the human genome.

Effects of replication domains on genome-wide UV-induced DNA damage and repair

Nucleotide excision repair is the primary repair mechanism that removes UV-induced DNA lesions in placentals. Unrepaired UV-induced lesions could result in mutations during DNA replication. Although the mutagenesis of pyrimidine dimers is reasonably well understood, the direct effects of replication fork progression on nucleotide excision repair are yet to be clarified. Here, we applied Damage-seq and XR-seq techniques and generated replication maps in synchronized UV-treated HeLa cells. The results suggest that ongoing replication stimulates local repair in both early and late replication domains. Additionally, it was revealed that lesions on lagging strand templates are repaired slower in late replication domains, which is probably due to the imbalanced sequence context. Asymmetric relative repair is in line with the strand bias of melanoma mutations, suggesting a role of exogenous damage, repair, and replication in mutational strand asymmetry.

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.

Insight into RNA-DNA primer length counting by human primosome

The human primosome, a four-subunit complex of primase and DNA polymerase alpha (Polα), synthesizes chimeric RNA–DNA primers of a limited length for DNA polymerases delta and epsilon to initiate DNA replication on both chromosome strands. Despite recent structural insights into the action of its two catalytic centers, the mechanism of DNA synthesis termination is still unclear. Here we report results of functional and structural studies revealing how the human primosome counts RNA–DNA primer length and timely terminates DNA elongation. Using a single-turnover primer extension assay, we defined two factors that determine a mature primer length (∼35-mer): (i) a tight interaction of the C-terminal domain of the DNA primase large subunit (p58_C) with the primer 5′-end, and (ii) flexible tethering of p58_C and the DNA polymerase alpha catalytic core domain (p180_core) to the primosome platform domain by extended linkers. The obtained data allow us to conclude that p58_C is a key regulator of all steps of RNA–DNA primer synthesis. The above-described findings provide a notable insight into the mechanism of DNA synthesis termination by a eukaryotic primosome, an important process for ensuring successful primer handover to replication DNA polymerases and for maintaining genome integrity.

Comprehensive Analysis of CRISPR-Cas9 Editing Outcomes in Yeast Xanthophyllomyces dendrorhous

DNA repair after Cas9 cutting can result in deletions/insertions, genomic rearrangements, and rare nucleotide substitutions. However, most work has only focused on deletions/insertions resulting from repair after CRISPR-Cas9 action. Here, we comprehensively analyzed the editing outcomes induced by CRISPR-Cas9 treatment in yeast Xanthophyllomyces dendrorhous by Sanger and Illumina sequencing and identified diverse DNA repair patterns, including DNA deletions, interchromosomal translocations, and on-target nucleotide substitutions (point mutations). Some deletions were observed repeatedly, and others, especially large deletions, varied in size. Genome sequencing and structural variation analysis showed that the interchromosomal translocations happened between Cas9 target sites and the endogenous ADH4 promoter. In contrast to previous studies, analysis revealed that the on-target point mutations were not random. Importantly, these point mutations showed strong sequence dependence that is not consistent with previous work in Hela cells, where CRISPR-mediated substitutions were found to lack sequence dependence and conversion preferences. Finally, we found that the non-homologous end joining components Ku70, Ku80, Mre11, or RAD50, and the overlapping roles of non-essential DNA polymerases were necessary for the production of both point mutations and deletions. This work expands our knowledge of CRISPR-Cas9 mediated DNA repair.

Structural and functional insight into mismatch extension by human DNA polymerase α

Human DNA polymerase α (Polα) does not possess proofreading ability and plays an important role in genome replication and mutagenesis. Polα extends the RNA primers generated by primase and provides a springboard for loading other replication factors. Here we provide the structural and functional analysis of the human Polα interaction with a mismatched template:primer. The structure of the human Polα catalytic domain in the complex with an incoming deoxycytidine triphosphate (dCTP) and the template:primer containing a T-C mismatch at the growing primer terminus was solved at a 2.9 Å resolution. It revealed the absence of significant distortions in the active site and in the conformation of the substrates, except the primer 3′-end. The T-C mismatch acquired a planar geometry where both nucleotides moved toward each other by 0.4 Å and 0.7 Å, respectively, and made one hydrogen bond. The binding studies conducted at a physiological salt concentration revealed that Polα has a low affinity to DNA and is not able to discriminate against a mispaired template:primer in the absence of deoxynucleotide triphosphate (dNTP). Strikingly, in the presence of cognate dNTP, Polα showed a more than 10-fold higher selectivity for a correct duplex versus a mismatched one. According to pre-steady-state kinetic studies, human Polα extends the T-C mismatch with a 249-fold lower efficiency due to reduction of the polymerization rate constant by 38-fold and reduced affinity to the incoming nucleotide by 6.6-fold. Thus, a mismatch at the postinsertion site affects all factors important for primer extension: affinity to both substrates and the rate of DNA polymerization.

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. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

PrimPol: A Breakthrough among DNA Replication Enzymes and a Potential New Target for Cancer Therapy

DNA replication can encounter blocking obstacles, leading to replication stress and genome instability. There are several mechanisms for evading this blockade. One mechanism consists of repriming ahead of the obstacles, creating a new starting point; in humans, PrimPol is responsible for carrying out this task. PrimPol is a primase that operates in both the nucleus and mitochondria. In contrast with conventional primases, PrimPol is a DNA primase able to initiate DNA synthesis de novo using deoxynucleotides, discriminating against ribonucleotides. In vitro, PrimPol can act as a DNA primase, elongating primers that PrimPol itself sythesizes, or as translesion synthesis (TLS) DNA polymerase, elongating pre-existing primers across lesions. However, the lack of evidence for PrimPol polymerase activity in vivo suggests that PrimPol only acts as a DNA primase. Here, we provide a comprehensive review of human PrimPol covering its biochemical properties and structure, in vivo function and regulation, and the processes that take place to fill the gap-containing lesion that PrimPol leaves behind. Finally, we explore the available data on human PrimPol expression in different tissues in physiological conditions and its role in cancer.

2-Dimensional Genetic Algorithm Exhibited an Essentiality of Gene Interaction for Evolution

Organisms consist of several genetic factors differing between species. However, the evolutionary effects of gene interactions on the evolutionary rate, adaptation, and divergence of organisms remain unknown. In a previous study, the 2-dimensional genetic algorithm (2DGA) program, including a gene interaction parameter, could simulate punctuated equilibrium under the disparity mode. Following this, we verified the effect of the number of gene interactions (gene cluster size) on evolution speed, adaptation, and divergence using the advanced 2DGA program. In this program, the population was replicated, mutated, and selected for 200,000 generations, and the fitness score, divergence, number of population, and genotype were output and plotted. The genotype data were used for evaluating the phylogenetic relations among the population. The gene cluster size 1) affected the disparity and parity mutagenesis modes differently, 2) determined the growth/exclusion rate and error threshold, and 3) accelerated or decelerated the population's speed of evolutionary advancement. In particular, when the gene cluster size expanded, the rate of increase in fitness scores decreased independently of the mutation rate and mode of mutation (disparity mode/parity mode). The mutation rate at the error threshold was also decreased by expanding the gene cluster size. Dendrograms traced the genotypes of the simulated population, indicating that the disparity mode caused the evolutionary process to enter 1) a stun mode, 2) an evolution mode, or 3) a divergence mode based on the mutation rate and gene cluster size, while the parity mode did not cause the population to enter a stun mode. Based on the above findings, we compared the predictions of the present study with evolution observed in the laboratory or the natural world and the processes of ongoing virus evolution, suggesting that our findings possibly explained the real evolution.

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.

Charles J, Dahl J, Farrell CE, Williams JS, Kunkel TA. 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.

Mutational Effects of Mobile Introns on the Mitochondrial Genomes of Metschnikowia Yeasts

It has been argued that DNA repair by homologous recombination in the context of endonuclease-mediated cleavage can cause mutations. To better understand this phenomenon, we examined homologous recombination following endonuclease cleavage in a native genomic context: the movement of self-splicing introns in the mitochondrial genomes of Metschnikowia yeasts. Self-splicing mitochondrial introns are mobile elements, which can copy and paste themselves at specific insertion sites in mitochondrial DNA using a homing endonuclease in conjunction with homologous recombination. Here, we explore the mutational effects of self-splicing introns by comparing sequence variation within the intron-rich cox1 and cob genes from 71 strains (belonging to 40 species) from the yeast genus Metschnikowia. We observed a higher density of single nucleotide polymorphisms around self-splicing-intron insertion sites. Given what is currently known about the movement of organelle introns, it is likely that their mutational effects result from the high binding affinity of endonucleases and their interference with repair machinery during homologous recombination (or, alternatively, via gene conversion occurring during the intron insertion process). These findings suggest that there are fitness costs to harbouring self-splicing, mobile introns and will help us better understand the risks associated with modern biotechnologies that use endonuclease-mediated homologous recombination, such as CRISPR-Cas9 gene editing.

RESCOT: Restriction Enzyme Set and Combination Optimization Tools for rNMP Capture Techniques

The incorporation of ribonucleoside monophosphates (rNMPs) in genomic DNA is a frequent phenomenon in many species, often associated with genome instability and disease. The ribose-seq technique is one of a few techniques designed to capture and map rNMPs embedded in genomic DNA. The first step of ribose-seq is restriction enzyme (RE) fragmentation, which cuts the genome into smaller fragments for subsequent rNMP capture. The RE selection chosen for genomic DNA fragmentation in the first step of the rNMP-capture techniques determines the genomic regions in which the rNMPs can be captured. Here, we designed a computational method, Restriction Enzyme Set and Combination Optimization Tools (RESCOT), to calculate the genomic coverage of rNMP-captured regions for a given RE set and to optimize the RE set to significantly increase the rNMP-captured-region coverage. Analyses of ribose-seq libraries for which the RESCOT tools were applied reveal that many rNMPs were captured in the expected genomic regions. Since different rNMP-mapping techniques utilize RE fragmentation and purification steps based on size-selection of the DNA fragments in the protocol, we discuss the possible usage of RESCOT for other rNMP-mapping techniques. In summary, RESCOT generates optimized RE sets for the fragmentation step of many rNMP capture techniques to maximize rNMP capture rate and thus enable researchers to better study characteristics of rNMP incorporation.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Frequency and patterns of ribonucleotide incorporation around autonomously replicating sequences in yeast reveal the division of labor of replicative DNA polymerases

Ribonucleoside triphosphate (rNTP) incorporation in DNA by DNA polymerases is a frequent phenomenon that results in DNA structural change and genome instability. However, it is unclear whether the rNTP incorporation into DNA follows any specific sequence patterns. We analyzed multiple datasets of ribonucleoside monophosphates (rNMPs) embedded in DNA, generated from three rNMP-sequencing techniques. These rNMP libraries were obtained from Saccharomyces cerevisiae cells expressing wild-type or mutant replicative DNA polymerase and ribonuclease H2 genes. We performed computational analyses of rNMP sites around early and late-firing autonomously replicating sequences (ARSs) of the yeast genome, where leading and lagging DNA synthesis starts bidirectionally. We found the preference of rNTP incorporation on the leading strand in wild-type DNA polymerase yeast cells. The leading/lagging-strand ratio of rNTP incorporation changes dramatically within the first 1,000 nucleotides from ARSs, highlighting the Pol δ - Pol ϵ handoff during early leading-strand synthesis. Furthermore, the pattern of rNTP incorporation is markedly distinct between the leading and lagging strands not only in mutant but also in wild-type polymerase cells. Such specific signatures of Pol δ and Pol ϵ provide a new approach to track the labor of these polymerases.

Genome-wide mapping of genomic DNA damage: methods and implications

Exposures from the external and internal environments lead to the modification of genomic DNA, which is implicated in the cause of numerous diseases, including cancer, cardiovascular, pulmonary and neurodegenerative diseases, together with ageing. However, the precise mechanism(s) linking the presence of damage, to impact upon cellular function and pathogenesis, is far from clear. Genomic location of specific forms of damage is likely to be highly informative in understanding this process, as the impact of downstream events (e.g. mutation, microsatellite instability, altered methylation and gene expression) on cellular function will be positional—events at key locations will have the greatest impact. However, until recently, methods for assessing DNA damage determined the totality of damage in the genomic location, with no positional information. The technique of “mapping DNA adductomics” describes the molecular approaches that map a variety of forms of DNA damage, to specific locations across the nuclear and mitochondrial genomes. We propose that integrated comparison of this information with other genome-wide data, such as mutational hotspots for specific genotoxins, tumour-specific mutation patterns and chromatin organisation and transcriptional activity in non-cancerous lesions (such as nevi), pre-cancerous conditions (such as polyps) and tumours, will improve our understanding of how environmental toxins lead to cancer. Adopting an analogous approach for non-cancer diseases, including the development of genome-wide assays for other cellular outcomes of DNA damage, will improve our understanding of the role of DNA damage in pathogenesis more generally.