Global landscape of replicative DNA polymerase usage in the human 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.

Strategic targeting of Cas9 nickase induces large segmental duplications

Gene/segmental duplications play crucial roles in genome evolution and variation. Here, we introduce paired nicking-induced amplification (PNAmp) for their experimental induction. PNAmp strategically places two Cas9 nickases upstream and downstream of a replication origin on opposite strands. This configuration directs the sister replication forks initiated from the origin to break at the nicks, generating a pair of one-ended double-strand breaks. If homologous sequences flank the two break sites, then end resection converts them to single-stranded DNAs that readily anneal to drive duplication of the region bounded by the homologous sequences. PNAmp induces duplication of segments as large as ∼1 Mb with efficiencies exceeding 10% in the budding yeast Saccharomyces cerevisiae. Furthermore, appropriate splint DNAs allow PNAmp to duplicate/multiplicate even segments not bounded by homologous sequences. We also provide evidence for PNAmp in mammalian cells. Therefore, PNAmp provides a prototype method to induce structural variations by manipulating replication fork progression.

Replicating chromatin in the nucleus: A histone variant perspective

In eukaryotes, chromatin and DNA replication are intimately linked, whereby chromatin impacts DNA replication control while genome duplication involves recovery of chromatin organisation. Here, we review recent advances in this area using a histone variant lens. We highlight how nucleosomal features interplay with origin definition and how the order of origin firing links with chromatin states in early mammalian development. We next discuss histone recycling and de novo deposition at the fork to finally open on the post-replicative recovery of the chromatin landscape to promote maintenance of cell identity.

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.

Distinct features of ribonucleotides within genomic DNA in Aicardi-Goutières syndrome ortholog mutants of Saccharomyces cerevisiae

Ribonucleoside monophosphates (rNMPs) are abundantly found within genomic DNA of cells. The embedded rNMPs alter DNA properties and impact genome stability. Mutations in ribonuclease (RNase) H2, a key enzyme for rNMP removal, are associated with the Aicardi-Goutières syndrome (AGS), a severe neurological disorder. Here, we engineered orthologs of the human RNASEH2A-G37S and RNASEH2C-R69W AGS mutations in yeast Saccharomyces cerevisiae: rnh201-G42S and rnh203-K46W. Using the ribose-seq technique and the Ribose-Map bioinformatics toolkit, we unveiled rNMP abundance, composition, hotspots, and sequence context in these AGS-ortholog mutants. We found a high rNMP presence in the nuclear genome of rnh201-G42S-mutant cells, and an elevated rCMP content in both mutants, reflecting preferential cleavage of RNase H2 at rGMP. We discovered unique rNMP patterns in each mutant, showing differential activity of the AGS mutants on the leading or lagging replication strands. This study guides future research on rNMP characteristics in human genomes with AGS mutations.

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.

Dormant origin firing promotes head-on transcription-replication conflicts at transcription termination sites in response to BRCA2 deficiency

BRCA2 is a tumor suppressor protein responsible for safeguarding the cellular genome from replication stress and genotoxicity, but the specific mechanism(s) by which this is achieved to prevent early oncogenesis remains unclear. Here, we provide evidence that BRCA2 acts as a critical suppressor of head-on transcription-replication conflicts (HO-TRCs). Using Okazaki-fragment sequencing (Ok-seq) and computational analysis, we identified origins (dormant origins) that are activated near the transcription termination sites (TTS) of highly expressed, long genes in response to replication stress. Dormant origins are a source for HO-TRCs, and drug treatments that inhibit dormant origin firing led to a reduction in HO-TRCs, R-loop formation, and DNA damage. Using super-resolution microscopy, we showed that HO-TRC events track with elongating RNA polymerase II, but not with transcription initiation. Importantly, RNase H2 is recruited to sites of HO-TRCs in a BRCA2-dependent manner to help alleviate toxic R-loops associated with HO-TRCs. Collectively, our results provide a mechanistic basis for how BRCA2 shields against genomic instability by preventing HO-TRCs through both direct and indirect means occurring at predetermined genomic sites based on the pre-cancer transcriptome.

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.

rNMPID: a database for riboNucleoside MonoPhosphates in DNA

Motivation

Ribonucleoside monophosphates (rNMPs) are the most abundant non-standard nucleotides embedded in genomic DNA. If the presence of rNMP in DNA cannot be controlled, it can lead to genome instability. The actual regulatory functions of rNMPs in DNA remain mainly unknown. Considering the association between rNMP embedment and various diseases and cancer, the phenomenon of rNMP embedment in DNA has become a prominent area of research in recent years.

Results

We introduce the rNMPID database, which is the first database revealing rNMP-embedment characteristics, strand bias, and preferred incorporation patterns in the genomic DNA of samples from bacterial to human cells of different genetic backgrounds. The rNMPID database uses datasets generated by different rNMP-mapping techniques. It provides the researchers with a solid foundation to explore the features of rNMP embedded in the genomic DNA of multiple sources, and their association with cellular functions, and, in future, disease. It also significantly benefits researchers in the fields of genetics and genomics who aim to integrate their studies with the rNMP-embedment data.

Availability and implementation

rNMPID is freely accessible on the web at https://www.rnmpid.org.

Insights into common fragile site instability: DNA replication challenges at DNA repeat sequences

Common fragile sites (CFS) are specific genomic regions prone to chromosomal instability under conditions of DNA replication stress. CFSs manifest as breaks, gaps, and constrictions on metaphase chromosomes under mild replication stress. These replication-sensitive CFS regions are preferentially unstable during cancer development, as reflected by their association with copy number variants (CNVs) frequently arise in most tumor types. Over the years, it became clear that a combination of different characteristics underlies the enhanced sensitivity of CFSs to replication stress. As of today, there is a strong evidence that the core fragility regions along CFSs overlap with actively transcribed large genes with delayed replication timing upon replication stress. Recently, the mechanistic basis for CFS instability was further extended to regions which span topologically associated domain (TAD) boundaries, generating a fragility signature composed of replication, transcription and genome organization. The presence of difficult-to-replicate AT-rich repeats was one of the early features suggested to characterize a subgroup of CFSs. These long stretches of AT-dinucleotide have the potential to fold into stable secondary structures which may impede replication fork progression, leaving the region under-replicated. Here, we focus on the molecular mechanisms underlying repeat instability at CFSs and on the proteins involved in the resolution of secondary structure impediments arising along repetitive sequence elements which are essential for the maintenance of genome stability.

Transcription-Replication Conflicts as a Source of Genome Instability

Transcription and replication both require large macromolecular complexes to act on a DNA template, yet these machineries cannot simultaneously act on the same DNA sequence. Conflicts between the replication and transcription machineries (transcription-replication conflicts, or TRCs) are widespread in both prokaryotes and eukaryotes and have the capacity to both cause DNA damage and compromise complete, faithful replication of the genome. This review will highlight recent studies investigating the genomic locations of TRCs and the mechanisms by which they may be prevented, mitigated, or resolved. We address work from both model organisms and mammalian systems but predominantly focus on multicellular eukaryotes owing to the additional complexities inherent in the coordination of replication and transcription in the context of cell type-specific gene expression and higher-order chromatin organization.

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.

Nontraditional Roles of DNA Polymerase Eta Support Genome Duplication and Stability

DNA polymerase eta (Pol η) is a Y-family polymerase and the product of the POLH gene. Autosomal recessive inheritance of POLH mutations is the cause of the xeroderma pigmentosum variant, a cancer predisposition syndrome. This review summarizes mounting evidence for expanded Pol η cellular functions in addition to DNA lesion bypass that are critical for maintaining genome stability. In vitro, Pol η displays efficient DNA synthesis through difficult-to-replicate sequences, catalyzes D-loop extensions, and utilizes RNA-DNA hybrid templates. Human Pol η is constitutively present at the replication fork. In response to replication stress, Pol η is upregulated at the transcriptional and protein levels, and post-translational modifications regulate its localization to chromatin. Numerous studies show that Pol η is required for efficient common fragile site replication and stability. Additionally, Pol η can be recruited to stalled replication forks through protein-protein interactions, suggesting a broader role in replication fork recovery. During somatic hypermutations, Pol η is recruited by mismatch repair proteins and is essential for V_H gene A:T basepair mutagenesis. Within the global context of repeat-dense genomes, the recruitment of Pol η to perform specialized functions during replication could promote genome stability by interrupting pure repeat arrays with base substitutions. Alternatively, not engaging Pol η in genome duplication is costly, as the absence of Pol η leads to incomplete replication and increased chromosomal instability.

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.