OKseqHMM: a genome-wide replication fork directionality analysis toolkit

Motivation: During each cell division, tens of thousands of DNA replication origins are coordinately activated to ensure the complete duplication of the entire human genome. However, the progression of replication forks can be challenged by numerous factors. One such factor is transcription-replication conflicts (TRC), which can either be co-directional or head-on with the latter being revealed as more dangerous for genome integrity. Results: In order to study the direction of replication fork movement and TRC, we developed a bioinformatics tool, called OKseqHMM, to directly measure the genome-wide replication fork directionality (RFD) as well as replication initiation and termination from data obtained by Okazaki fragment sequencing (OK-Seq) and related techniques. Availability and Implementation: We have gathered and analyzed OK-seq data from a large number of organisms including yeast, mouse and human, to generate high-quality RFD profiles and determine initiation zones and termination zones by using Hidden Markov Model (HMM) algorithm (all tools and data are available at https://github.com/CL-CHEN-Lab/OK-Seq). In addition, we have extended our analysis to data obtained by related techniques, such as eSPAN and TrAEL-seq, which also contain RFD information. Our works, therefore, provide an important tool and resource for the community to further study TRC and genome instability, in a wide range of cell line models and growth conditions, which is of prime importance for human health.

Single-molecule mapping of replisome progression

Fundamental aspects of DNA replication, such as the anatomy of replication stall sites, how replisomes are influenced by gene transcription and whether the progression of sister replisomes is coordinated are poorly understood. Available techniques do not allow the precise mapping of the positions of individual replisomes on chromatin. We have developed a new method called Replicon-seq that entails the excision of full-length replicons by controlled nuclease cleavage at replication forks. Replicons are sequenced using Nanopore, which provides a single molecule readout of long DNA molecules. Using Replicon-seq, we have investigated replisome movement along chromatin. We found that sister replisomes progress with remarkable consistency from the origin of replication but function autonomously. Replication forks that encounter obstacles pause for a short duration but rapidly resume synthesis. The helicase Rrm3 plays a critical role both in mitigating the effect of protein barriers and facilitating efficient termination. Replicon-seq provides an unprecedented means of defining replisome movement across the genome.

Chl1 helicase controls replication fork progression by regulating dNTP pools

Eukaryotic cells have evolved a replication stress response that helps to overcome stalled/collapsed replication forks and ensure proper DNA replication. The replication checkpoint protein Mrc1 plays important roles in these processes, although its functional interactions are not fully understood. Here, we show that MRC1 negatively interacts with CHL1, which encodes the helicase protein Chl1, suggesting distinct roles for these factors during the replication stress response. Indeed, whereas Mrc1 is known to facilitate the restart of stalled replication forks, we uncovered that Chl1 controls replication fork rate under replication stress conditions. Chl1 loss leads to increased RNR1 gene expression and dNTP levels at the onset of S phase likely without activating the DNA damage response. This in turn impairs the formation of RPA-coated ssDNA and subsequent checkpoint activation. Thus, the Chl1 helicase affects RPA-dependent checkpoint activation in response to replication fork arrest by ensuring proper intracellular dNTP levels, thereby controlling replication fork progression under replication stress conditions.

The Bacillus subtilis PriA winged helix domain is critical for surviving DNA damage

DNA replication forks regularly encounter lesions or other impediments that result in a blockage to fork progression. PriA is one of the key proteins used by virtually all eubacteria to survive conditions that result in a blockage to replication fork movement. PriA directly binds stalled replication forks and initiates fork restart allowing for chromosomes to be fully duplicated under stressful conditions. We used a CRISPR-Cas gene editing approach to map PriA residues critical for surviving DNA damage induced by several antibiotics in B. subtilis . We find that the winged helix (WH) domain in B. subtilis PriA is critical for surviving DNA damage and participates in DNA binding. The critical in vivo function of the WH domain mapped to distinct surfaces that were also conserved among several Gram-positive human pathogens. In addition, we identified an amino acid linker neighboring the WH domain that is greatly extended in B. subtilis due to an insertion. Shortening this linker induced a hypersensitive phenotype to DNA damage, suggesting that its extended length is critical for efficient replication fork restart in vivo . Because the WH domain is dispensable in E. coli PriA, our findings demonstrate an important difference in the contribution of the WH domain during fork restart in B. subtilis . Further, with our results we suggest that this highly variable region in PriA could provide different functions across diverse bacterial organisms. IMPORTANCE PriA is an important protein found in virtually all bacteria that recognizes stalled replication forks orchestrating fork restart. PriA homologs contain a winged helix (WH) domain which is dispensable in E. coli and functions in a fork restart pathway that is not conserved outside of E. coli and closely related proteobacteria. We analyzed the importance of the WH domain and an associated linker in B. subtilis and found that both are critical for surviving DNA damage. This function mapped to a small motif at the C-terminal end of the WH domain, which is also conserved in pathogenic bacteria. The motif was not required for DNA binding and therefore may perform a novel function in the replication fork restart pathway.

Genetic Study of Four Candidate Holliday Junction Processing Proteins in the Thermophilic Crenarchaeon Sulfolobus acidocaldarius

Homologous recombination (HR) is thought to be important for the repair of stalled replication forks in hyperthermophilic archaea. Previous biochemical studies identified two branch migration helicases (Hjm and PINA) and two Holliday junction (HJ) resolvases (Hjc and Hje) as HJ-processing proteins; however, due to the lack of genetic evidence, it is still unclear whether these proteins are actually involved in HR in vivo and how their functional relation is associated with the process. To address the above questions, we constructed hjc-, hje-, hjm-, and pina single-knockout strains and double-knockout strains of the thermophilic crenarchaeon Sulfolobus acidocaldarius and characterized the mutant phenotypes. Notably, we succeeded in isolating the hjm- and/or pina-deleted strains, suggesting that the functions of Hjm and PINA are not essential for cellular growth in this archaeon, as they were previously thought to be essential. Growth retardation in Δpina was observed at low temperatures (cold sensitivity). When deletion of the HJ resolvase genes was combined, Δpina Δhjc and Δpina Δhje exhibited severe cold sensitivity. Δhjm exhibited severe sensitivity to interstrand crosslinkers, suggesting that Hjm is involved in repairing stalled replication forks, as previously demonstrated in euryarchaea. Our findings suggest that the function of PINA and HJ resolvases is functionally related at lower temperatures to support robust cellular growth, and Hjm is important for the repair of stalled replication forks in vivo.

Mandatory coupling of zygotic transcription to DNA replication in early Drosophila embryos

Collisions between transcribing RNA polymerases and DNA replication forks are disruptive. The threat of collisions is particularly acute during the rapid early embryonic cell cycles of Drosophila when S phase occupies the entirety of interphase. We hypothesized that collision-avoidance mechanisms safeguard the onset of zygotic transcription in these cycles. To explore this hypothesis, we used real-time imaging of transcriptional events at the onset of each interphase. Endogenously tagged RNA polymerase II (RNAPII) abruptly formed clusters before nascent transcripts accumulated, indicating recruitment prior to transcriptional engagement. Injection of inhibitors of DNA replication prevented RNAPII clustering, blocked formation of foci of the pioneer factor Zelda, and largely prevented expression of transcription reporters. Knockdown of Zelda or the histone acetyltransferase CBP prevented RNAPII cluster formation except at the replication-dependent (RD) histone gene locus. We suggest a model in which the passage of replication forks allows Zelda and a distinct pathway at the RD histone locus to reconfigure chromatin to nucleate RNAPII clustering and promote transcriptional initiation. The replication dependency of these events defers initiation of transcription and ensures that RNA polymerases transcribe behind advancing replication forks. The resulting coordination of transcription and replication explains how early embryos circumvent collisions and promote genome stability.

SUMO-Based Regulation of Nuclear Positioning to Spatially Regulate Homologous Recombination Activities at Replication Stress Sites

DNA lesions have properties that allow them to escape their nuclear compartment to achieve DNA repair in another one. Recent studies uncovered that the replication fork, when its progression is impaired, exhibits increased mobility when changing nuclear positioning and anchors to nuclear pore complexes, where specific types of homologous recombination pathways take place. In yeast models, increasing evidence points out that nuclear positioning is regulated by small ubiquitin-like modifier (SUMO) metabolism, which is pivotal to maintaining genome integrity at sites of replication stress. Here, we review how SUMO-based pathways are instrumental to spatially segregate the subsequent steps of homologous recombination during replication fork restart. In particular, we discussed how routing towards nuclear pore complex anchorage allows distinct homologous recombination pathways to take place at halted replication forks.

Delineation of the Ancestral Tus-Dependent Replication Fork Trap

In Escherichia coli, DNA replication termination is orchestrated by two clusters of Ter sites forming a DNA replication fork trap when bound by Tus proteins. The formation of a ‘locked’ Tus–Ter complex is essential for halting incoming DNA replication forks. However, the absence of replication fork arrest at some Ter sites raised questions about their significance. In this study, we examined the genome-wide distribution of Tus and found that only the six innermost Ter sites (TerA–E and G) were significantly bound by Tus. We also found that a single ectopic insertion of TerB in its non-permissive orientation could not be achieved, advocating against a need for ‘back-up’ Ter sites. Finally, examination of the genomes of a variety of Enterobacterales revealed a new replication fork trap architecture mostly found outside the Enterobacteriaceae family. Taken together, our data enabled the delineation of a narrow ancestral Tus-dependent DNA replication fork trap consisting of only two Ter sites.

Mechanism of replication fork reversal and protection by human RAD51 and RAD51 paralogs

SMARCAL1, ZRANB3 and HLTF are all required for the remodeling of replication forks upon stress. Using reconstituted reactions, we show that the motor proteins have unequal biochemical capacities, explaining why they have non-redundant functions. Whereas SMARCAL1 uniquely anneals RPA-coated ssDNA, suggesting an initial function in fork reversal, it becomes comparatively inefficient in subsequent branch migration. We also show that low concentrations of RAD51 and the RAD51 paralog complex, RAD51B-RAD51C-RAD51D-XRCC2 (BCDX2), directly stimulate SMARCAL1 and ZRANB3 but not HLTF, providing a mechanism underlying previous cellular data implicating these factors in fork reversal. Upon reversal, RAD51 protects replication forks from degradation by MRE11, DNA2 and EXO1 nucleases. We show that the protective function of RAD51 unexpectedly depends on its binding to double-stranded DNA, and higher RAD51 concentrations are required for DNA protection compared to reversal. Together, we define the non-canonical functions of RAD51 and its paralogs in replication fork reversal and protection.