OKseqHMM: a genome-wide replication fork directionality analysis toolkit

Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression

BACKGROUND

DNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression.

RESULTS

To characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density.

CONCLUSIONS

Altogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

Linear interaction between replication and transcription shapes DNA break dynamics at recurrent DNA break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unravel the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA reveals co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity uncovers that head-to-head interactions between replication and transcription machinery result in 60% DSB contribution to the head-on compared to 40% for co-directional. Taken together we reveal TTR as a fragile class and show how the linear interaction between transcription and replication impacts genome stability.

A Decade of Discovery-Eukaryotic Replisome Disassembly at Replication Termination

The eukaryotic replicative helicase (CMG complex) is assembled during DNA replication initiation in a highly regulated manner, which is described in depth by other manuscripts in this Issue. During DNA replication, the replicative helicase moves through the chromatin, unwinding DNA and facilitating nascent DNA synthesis by polymerases. Once the duplication of a replicon is complete, the CMG helicase and the remaining components of the replisome need to be removed from the chromatin. Research carried out over the last ten years has produced a breakthrough in our understanding, revealing that replication termination, and more specifically replisome disassembly, is indeed a highly regulated process. This review brings together our current understanding of these processes and highlights elements of the mechanism that are conserved or have undergone divergence throughout evolution. Finally, we discuss events beyond the classic termination of DNA replication in S-phase and go over the known mechanisms of replicative helicase removal from chromatin in these particular situations.

Fork coupling directs DNA replication elongation and termination

DNA replication is initiated at multiple loci to ensure timely duplication of eukaryotic genomes. Sister replication forks progress bidirectionally, and replication terminates when two convergent forks encounter one another. To investigate the coordination of replication forks, we developed a replication-associated in situ HiC method to capture chromatin interactions involving nascent DNA. We identify more than 2000 fountain-like structures of chromatin contacts in human and mouse genomes, indicative of coupling of DNA replication forks. Replication fork interaction not only occurs between sister forks but also involves forks from two distinct origins to predetermine replication termination. Termination-associated chromatin fountains are sensitive to replication stress and lead to coupled forks-associated genomic deletions in cancers. These findings reveal the spatial organization of DNA replication forks within the chromatin context.

Linear Interaction Between Replication and Transcription Shapes DNA Break Dynamics at Recurrent DNA Break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unraveled the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA revealed co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity revealed that head-to-head interactions between replication and transcription machinery resulted in 60% DSB contribution to the head-on compared to 40% for co-directional. Our findings revealed TTR as a novel fragile class and highlighted how the linear interaction between transcription and replication impacts genome stability.

Asymmetric distribution of parental H3K9me3 in S phase silences L1 elements

In eukaryotes, repetitive DNA sequences are transcriptionally silenced through histone H3 lysine 9 trimethylation (H3K9me3). Loss of silencing of the repeat elements leads to genome instability and human diseases, including cancer and ageing1-3. Although the role of H3K9me3 in the establishment and maintenance of heterochromatin silencing has been extensively studied4-6, the pattern and mechanism that underlie the partitioning of parental H3K9me3 at replicating DNA strands are unknown. Here we report that H3K9me3 is preferentially transferred onto the leading strands of replication forks, which occurs predominantly at long interspersed nuclear element (LINE) retrotransposons (also known as LINE-1s or L1s) that are theoretically transcribed in the head-on direction with replication fork movement. Mechanistically, the human silencing hub (HUSH) complex interacts with the leading-strand DNA polymerase Pol ε and contributes to the asymmetric segregation of H3K9me3. Cells deficient in Pol ε subunits (POLE3 and POLE4) or the HUSH complex (MPP8 and TASOR) show compromised H3K9me3 asymmetry and increased LINE expression. Similar results were obtained in cells expressing a MPP8 mutant defective in H3K9me3 binding and in TASOR mutants with reduced interactions with Pol ε. These results reveal an unexpected mechanism whereby the HUSH complex functions with Pol ε to promote asymmetric H3K9me3 distribution at head-on LINEs to suppress their expression in S phase.

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.