Transcription shapes DNA replication initiation to preserve genome integrity
Genome Biology volume 22, Article number: 176 (2021)
Early DNA replication occurs within actively transcribed chromatin compartments in mammalian cells, raising the immediate question of how early DNA replication coordinates with transcription to avoid collisions and DNA damage.
We develop a high-throughput nucleoside analog incorporation sequencing assay and identify thousands of early replication initiation zones in both mouse and human cells. The identified early replication initiation zones fall in open chromatin compartments and are mutually exclusive with transcription elongation. Of note, early replication initiation zones are mainly located in non-transcribed regions adjacent to transcribed regions. Mechanistically, we find that RNA polymerase II actively redistributes the chromatin-bound mini-chromosome maintenance complex (MCM), but not the origin recognition complex (ORC), to actively restrict early DNA replication initiation outside of transcribed regions. In support of this finding, we detect apparent MCM accumulation and DNA replication initiation in transcribed regions due to anchoring of nuclease-dead Cas9 at transcribed genes, which stalls RNA polymerase II. Finally, we find that the orchestration of early DNA replication initiation by transcription efficiently prevents gross DNA damage.
RNA polymerase II redistributes MCM complexes, but not the ORC, to prevent early DNA replication from initiating within transcribed regions. This RNA polymerase II-driven MCM redistribution spatially separates transcription and early DNA replication events and avoids the transcription-replication initiation collision, thereby providing a critical regulatory mechanism to preserve genome stability.
The mammalian genome is categorized into active and inactive compartments linked to key cellular processes according to chromatin activity and histone modifications, and transcription mainly occurs in active chromatin compartments [1,2,3,4]. Early DNA replication is initiated within active chromatin compartments, followed by elongation into inactive compartments . Therefore, early DNA replication and gene transcription both occur within active chromatin compartments, raising the pivotal question of how these processes are spatially and temporally coordinated to avoid replication-transcription collisions and subsequent DNA damage [6,7,8,9]. Previous studies suggest that transcription might affect DNA replication initiation, but the mechanism underlying this interaction in mammalian cells remains elusive [10,11,12,13,14,15].
DNA replication occurs throughout the entire genome, in contrast to transcription, and tracing early DNA replication requires precise identification of genomic DNA replication origins and initiation zones. Replication origins have been well characterized in bacteria and yeast, but their locations in mammalian cells remain unclear [16, 17]. Many attempts have been made to identify replication origins and initiation zones in mammalian cells. Chromatin immunoprecipitation-sequencing (ChIP-seq) of pre-replication complex (pre-RC) components is difficult because high-quality antibodies against the origin recognition complex (ORC) and mini-chromosome maintenance complex (MCM) are not available . Repli-seq has been applied in studies exploring the determination of replication timing, which have suggested that replication domains range from 400 to 800 kb in size [5, 19]. Recently, high-resolution Repli-seq was employed to profile initiation zones by dissecting the S phase into 16 fractions . The dNTP synthesis inhibitor hydroxyurea (HU) can induce replication stress to slow replication elongation and thus greatly improve the resolution of EdU-marked DNA replication initiation (EdU-seq-HU) experiments [13, 15]. Small nascent DNA strands (SNS) have been used to map DNA replication origins [21,22,23]. Moreover, Okazaki fragments (OK) generated during DNA replication have also been employed to identify initiation zones that harbor one or more origin(s) within each zone [10, 24]. These methods reveal that early DNA replication initiation occurs in active chromatin compartments in which transcription also occurs .
The process of replication initiation requires origin recognition by the ORC, MCM loading, and MCM activation [6, 25]. MCM complexes are loaded as inactive double hexamers in a process dependent on the ORC during the late M and G1 phases. When a cell enters the S phase, a proportion of the MCM complexes are activated to begin unwinding double-stranded (ds)DNA, initiating two replication forks that move bidirectionally [26,27,28]. The ORC stably binds chromatin, while ring-shaped MCM double hexamers encircling dsDNA may slide along the chromatin, uncoupling the ORC and MCM [29,30,31,32]. MCM redistribution has been reported in Drosophila melanogaster, Xenopus egg extract, and asynchronized human HeLa cells and 2fTGH cells, but the mechanism underlying this process remains unexplored [33,34,35]. RNA polymerases are among the most likely candidates responsible for MCM redistribution. T7 RNA polymerase can force yeast MCM complexes to fall off linear dsDNA in vitro . In addition, RNA polymerase II is also able to cause MCM relocation at ribosomal DNA loci in budding yeast . However, MCM redistribution by the transcription machinery in mammalian cells has not been confirmed experimentally [32, 35].
Here, we developed a new assay, nucleoside analog incorporation loci sequencing (NAIL-seq), to precisely map genome-wide early replication initiation zones (ERIZs). ERIZs are primarily located in non-transcribed regions in open chromatin compartments, which are mutually exclusive with transcription elongation regions. Furthermore, inhibition of transcription leads to MCM redistribution and early replication, but not ORC re-localization, in transcribed regions. Failure by cells to prevent DNA replication initiation in transcribed regions leads to gross DNA damage. Therefore, we propose that transcription regulates early DNA replication initiation by redistributing MCM complexes to non-transcribed regions to avoid replication-transcription collisions and thereby preserve genome stability.
NAIL-seq identifies early DNA replication initiation zones
To identify replication initiation zones with the aim of tracing early DNA replication, we synchronized human GM12878 (primary-like) and K562 (leukemia) cells separately in the G1 phase using the CDK4/6 inhibitor palbociclib . More than 94% of the treated cells were arrested in the G1 phase after treatment for 36 h (Fig. 1a; Additional file 1: Figure S1a). After a release of 2–3 h from the G1 phase, the synchronized K562 and GM12878 cells reached the G1/S transition (Additional file 1: Figure S1a and b). The K562 and GM12878 cells were labeled after a release of 2.5 or 3 h, respectively, with EdU and then BrdU, each for a 15-min pulse (Fig. 1a; Additional file 1: Figure S1a and b). Theoretically, the two thymidine analog signals can mark replication initiation sites and adjacent fork elongation regions for tens of kilobases. The EdU- and BrdU-incorporated fragments can be distinguished by using streptavidin C1 beads after the Click reaction or an anti-BrdU antibody, respectively (Additional file 1: Figure S1c). Of note, EdU cross-reacted with the anti-BrdU antibody; nevertheless, this cross-reaction was undetectable after the Click reaction (Additional file 1: Figure S1d). Incorporation of EdU or BrdU from the sequential dual-labeling revealed strong early DNA replication signals in the early replication domains of synchronized K562 cells (Fig. 1b; Additional file 1: Figure S2a). Moreover, the “E-B” signal obtained by subtracting the BrdU signal from the EdU signal was further narrowed to the middle of the early replication domains (Fig. 1b; Additional file 1: Figure S2a). The resolution of the E-B peaks was more than 2-fold higher than that of the individual EdU and BrdU peaks, with a median width of 90 kb versus 245 and 290 kb, respectively, in K562 cells (Fig. 1c; Additional file 1: Figure S2b).
Next, we employed a second assay, EdU-seq-HU, by slowly incorporating EdU under HU treatment to map early replication initiation [13, 15]. G1-arrested cells were released into fresh medium containing HU and EdU for an additional 12 h prior to harvesting (Fig. 1a; Additional file 1: Figure S1c). Approximately 71.8% of EdU/HU hotspots in K562 cells overlapped with 78.9% of the E-B peaks (Fig. 1d). The EdU/HU hotspots showed a median width of 55 kb, exhibiting a slightly higher resolution than that of the E-B peaks (Fig. 1c; Additional file 1: Figure S2b). HU treatment yields high-resolution replication-associated peaks [13, 15], but it may induce DNA double-stranded breaks (DSBs) and early utilization of dormant DNA replication origins [39,40,41,42]. E-B signals are obtained from cells undergoing the G1/S transition without replication stress; therefore, E-B signals can be used to exclude early firing of dormant origins and potential DNA damage in the EdU/HU libraries. The EdU/HU peaks overlapping the E-B peaks showed earlier replication timing in comparison with that of the non-overlapping EdU/HU peaks (Fig. 1e). Therefore, we defined the EdU/HU peaks overlapping the E-B peaks as ERIZs and subjected them to further analysis (Fig. 1a, b, and d; Additional file 1: Figure S2a). Since EdU-seq-HU also employed a nucleoside analog for sequencing, this method was termed nucleoside analog incorporation loci sequencing (NAIL-seq). NAIL-seq identified 2265 and 2874 ERIZs with a median size of 71 kb and 76 kb in K562 and GM12878 cells, respectively (Fig. 1f; Additional file 2: Table S1; Additional file 3: Table S2). As expected, more than 95% of the identified ERIZs were found in early replication regions in K562 and GM12878 cells (Additional file 1: Figure S2c). Of note, only 0.04% (1 of 2265) of ERIZs in K562 cells and 1% (30 of 2874) of those in GM12878 cells occurred in typical late replication regions.
NAIL-seq has higher resolution than Repli-seq
To determine whether NAIL-seq is superior to typical single-signal incorporation assays, we used published Repli-seq data from asynchronized K562 cells to perform a direct comparison of the methods . The individual EdU and BrdU peaks from NAIL-seq had a markedly higher resolution in comparison with that of the early replication peaks identified by Repli-seq, which had a median peak width of 340 kb (Fig. 1b, c; Additional file 1: Figure S2a), suggesting that cell synchronization improved the resolution of the replication initiation zones identified by the experiments. With regard to ERIZs, the narrowest ERIZs identified by NAIL-seq were approximately 10 kb in size, while the narrowest early replication domains from the Repli-seq results were more than 100 kb in size (Fig. 1c, f). OK-seq and SNS-seq have also been widely used to identify initiation zones (IZs) and replication origins, respectively, associated with replication events throughout the entire S phase in unsynchronized cells. In K562 cells, nearly 50% of ERIZs identified by NAIL-seq overlapped with 43.4% of the IZs identified by OK-seq , and nearly 70% of these ERIZs overlapped with 16.3% of the replication origins identified by SNS-seq  (Additional file 1: Figure S2d). Of note, we also found that two or more replication origins identified by SNS-seq fell into a single ERIZ (highlighted in cyan in the Additional file 1: Figure S2a). Moreover, we performed EdU-seq-HU in stimulated mouse splenic B cells and the identified early replication zones were highly correlated with previous results  (Additional file 1: Figure S2e).
Early DNA replication is present in non-transcribed regions
Next, we sought to investigate the relationship between transcription and early replication initiation. Roughly 94−97% of the clustered ERIZs were present in transcription-occupied A compartments  in K562 and GM12878 cells (Fig. 2a), consistent with previous reports . To further characterize the chromatin localization of ERIZs, we built a logistic regression model to demonstrate the correlation between ERIZs and histone modification markers, as well as chromatin structural proteins, based on ChIP-seq data from ENCODE. The most significant predictive marker for ERIZs was H2A.Z (Fig. 2b; Additional file 1: Figure S3a), in line with a recent report that H2A.Z recruits ORC1 to initiate DNA replication . Cohesin, CTCF, and other histone markers for transcription enhancers or promoters, including H3K4me1, H3K4me3, and H3K27ac, also had positive prediction scores for ERIZs (Fig. 2b; Additional file 1: Figure S3a). However, transcription elongation markers H3K79me2 and H3K36me3, as well as transcription silencer marker H3K9me3, had negative prediction scores for ERIZs (Fig. 2b).
To further explore the relationship between transcription and early DNA replication, we examined the zoomed-in profile of global run-on (GRO)-seq transcription signals and ERIZ-associated early replication signals. Remarkably, we found that ERIZs are mainly located in non-transcribed regions, mutually exclusive with transcribed gene bodies, despite being in the same A compartments (Fig. 2c). The genome-wide analysis confirmed that ERIZ signals accumulated in non-transcribed regions flanked by transcription signals of GRO-seq  in GM12878 and K562 cells (Fig. 2d). We also identified ERIZs from two primary cell types, mouse embryonic stem cells (mESCs) and mouse splenic B cells (mSBCs), and obtained similar findings (Additional file 1: Figure S3b). Of note, the mESCs were arrested by thymidine and then nocodazole , after which they were released into medium containing HU for 3 h to allow the cells to enter the early S phase before harvesting (Additional file 1: Figure S3c and d). The G0-phase mSBCs were activated by cytokines in medium containing HU for 28 h to allow the cells to enter the early S phase  before harvesting (Additional file 1: Figure S3c). These results explain the negative correlation between ERIZs and transcription elongation markers and suggest that early DNA replication initiation occurs preferentially at non-transcribed regions within active chromatin compartments.
Transcription shapes early DNA replication initiation
To investigate whether transcription interacts with DNA replication initiation, we used α-amanitin to inhibit transcription  in synchronized K562 cells, after which we measured the impact of transcription perturbation on the distribution of ERIZs. We released G1-arrested K562 cells into fresh medium containing 2 or 10 μg/mL α-amanitin, in addition to HU and EdU, for 12 h before harvesting (Additional file 1: Figure S4a). Transcription blockade by α-amanitin slightly reduced EdU incorporation, but robust EdU signals were detected without significant elevation in DNA damage marked by γ-H2AX (Additional file 1: Figure S4b and c). Enrichment of early replication signals in ERIZ-associated non-transcribed regions (as defined in Fig. 2d) decreased as the α-amanitin concentration was increased (Fig. 3a; Additional file 1: Figure S4d). Early replication signals also showed increased enrichment at transcription start sites (TSSs) following α-amanitin treatment (Figs. 3a and 4f; Additional file 1: Figure S4d), possibly because TSS regions where ORC is located support DNA replication initiation, as previously reported [10, 12, 23, 48]. The ratio of early replication signals in non-transcribed regions to that of early replication signals in adjacent transcribed regions significantly decreased as the concentration of α-amanitin was increased, suggesting that early replication penetrated transcribed regions (Fig. 3b; Additional file 1: Figure S4d). Moreover, stronger inhibition of transcription resulted in more dramatic early replication redistribution on the genome (Fig. 3b and exemplified in Additional file 1: Figure S4d). Similarly, early replication signals also relocated from non-transcribed regions to neighboring transcribed regions in low-dose α-amanitin-treated mSBCs (Additional file 1: Figure S4e and S4f). We also applied another transcription inhibitor, 5,6-dichloro-1-β-d-ribofuranosyl benzimidazole (DRB), at concentrations ranging from 5 to 30 μM to interrupt transcription elongation in G1-arrested K562 cells and obtained similar, albeit milder, findings (Fig. 3c, d; Additional file 1: Figure S4g).
Since transcription inhibitors may have non-specific effects on DNA metabolism, we employed the auxin-inducible degron system to directly eliminate RNA polymerase II from mESCs in an acute manner as previously reported . Briefly, the catalytic subunit of RNA polymerase II was tagged with a mAID tag, and indole-3-acetic acid (IAA) was added to induce degradation of the fusion protein. mESCs were arrested at the M phase and then released into the early S phase (Additional file 1: Figure S4h). In order to eliminate RNA polymerase II before the G1/S transition, IAA was added 1 h before release from the M phase. RNA polymerase II was degraded almost completely after 2 h of induction by IAA, while the number of MCM complexes remained constant (Additional file 1: Figure S4h and i). In addition, rapid degradation of RNA polymerase II did not influence the entry of mESCs into the S phase (Additional file 1: Figuer S4j). IAA treatment also did not significantly change the distribution of ERIZ-associated early replication signals in wild-type cells (Fig. 3f; Additional file 1: Figure S4k; Additional file 4: Table S3). However, the absence of RNA polymerase II resulted in a remarkable decrease in the number of early replication signals in non-transcribed regions (Fig. 3e), in line with the findings obtained using transcription inhibitors in K562 cells. In addition, the number of early replication signals was declined in non-transcribed regions and increased in gene bodies, as exemplified in the Dnajc7 gene (Fig. 3f, g). Collectively, our data indicate that transcription actively orchestrates early DNA replication initiation in both human and mouse cells.
Transcription redistributes the MCM, but not the ORC, to non-transcribed regions
It has been reported that the ORC is distributed widely in both transcribed and non-transcribed regions, with enrichment at TSSs, and this finding was confirmed by re-analysis of published data from an ORC2 ChIP-seq study  in asynchronized K562 cells (Additional file 1: Figure S5a). In contrast, MCM complexes accumulate in non-transcribed regions, as shown by the ChIP-seq data for MCM3, MCM5, and MCM7 from ENCODE in asynchronized K562 cells, indicating that MCM complexes are uncoupled from the ORC (Additional file 1: Figure S5a).
To explore the mechanism by which transcription influences early DNA replication initiation, we examined the distribution of the ORC and MCM via ChIP-seq following α-amanitin-induced transcription perturbation in the G1 phase prior to DNA replication initiation. The protein levels of the ORC and MCM remained constant in the cytosol and on chromatin despite the degradation of RNA polymerase II in G1-arrested K562 cells following treatment with 10 μg/mL α-amanitin (Fig. 4a; Additional file 1: S5b). Of note, ORC2 ChIP-seq was achieved by expressing Avi-tagged ORC2 at an expression level much lower than that of endogenous ORC2, which had no impact on chromatin-bound MCM5 in the G1 phase (Additional file 1: Figure S5c). ORC2 and MCM5 showed distinct distribution patterns in the G1 phase without α-amanitin treatment, which were similar to those of asynchronized K562 cells (Fig. 4b left and 4d left; Additional file 1: Figure S5a, d, and e). Following α-amanitin treatment, ORC2 primarily remained at its original binding sites, showing similar distribution patterns in both non-transcribed and transcribed regions regardless of α-amanitin treatment (Fig. 4b, c). However, the MCM5 signal was diminished in non-transcribed regions, and the ratio of the MCM5 signal in non-transcribed regions to that of the signal in transcribed regions decreased significantly (Fig. 4d, e).
To better present changes in early replication, ORC, and MCM signals in transcribed regions, we aligned ERIZ-adjacent transcribed genes at both TSS and transcription termination sites (TTS), which revealed that early replication signals were rarely detected within transcribed genes before treatment (Fig. 4f). Following α-amanitin treatment, early replication signals were detected in transcribed regions, with obvious enrichment at TSSs (Fig. 4f), while the abundance of ORC2 remained unchanged within transcribed genes and was enriched at TSSs (Fig. 4b, c, and g). However, the MCM5 signal within transcribed genes was intensified, and mild signal pileups were observed within gene bodies (Fig. 4e, h). In particular, α-amanitin treatment resulted in similar distribution patterns for ERIZs and ORC2 (Fig. 4f right versus 4g right). These findings imply that RNA polymerase II redistributes MCM complexes outside of transcribed regions, driving uncoupling of the ORC and MCM complexes and thus preventing early DNA replication initiation at transcribed regions.
MCM accumulates at RNA polymerase II blockade sites to initiate DNA replication
The driving of MCM complexes by RNA polymerase II requires the coupling of these two complexes in transcribed regions in the G1 phase. To detect the coupling of RNA polymerase II and MCM double hexamers, we utilized a transcription barrier to stall RNA polymerase II and thereby trap MCM. In this assay, trapped MCM accumulates and may induce DNA replication initiation in the arrival direction of transcription upstream of the barrier. We designed four guide (g)RNA-target sites spanning a 162-bp exonic region of the CMIP gene without ORC binding capability to recruit nuclease-dead Cas9 (dCas9) as a transcription roadblock (Fig. 5a). The efficiency of transcription blockade by dCas9 reached 60−70%, and the upstream to downstream transcription levels were significantly increased (Additional file 1: Figure S6a-c). Both the upstream and downstream ERIZs neighboring CMIP still supported substantial early DNA replication following the dCas9-mediated transcription blockade as a result of the presence of the ORC in neighboring non-transcribed regions and inefficient transcription inhibition (Additional file 1: Figure S6b-d). Interestingly, we found enrichment of EdU signals within the 2-kb region upstream of the dCas9 binding sites in the arrival direction of transcription, but not in the other direction (Fig. 5a–c; Additional file 1: Figure S6e). To evaluate whether these EdU signals were produced by DNA replication initiation or DNA damage due to transcription stalling, we performed EdU-seq-HU in the G1-arrested K562 cells. We found no EdU enrichment within the same 2-kb region following transcription blockade (Fig. 5d, e; Additional file 1: Figure S6b and e). Given that both transcription and DNA replication occurred in the S phase, but only transcription occurred in the G1 phase, the enriched EdU signals in the early S phase must have originated from new DNA replication initiation. In addition, no significant early replication signals were observed within the gene body downstream of the dCas9 block site, possibly because ORC was not present (Fig. 5b, c; Additional file 1: Figure S6d and e). We also performed dCas9-induced transcription blockade at GALNT10. Although the blockade efficiency was lower than that observed at CMIP (Additional file 1: Figure S6b and c), we observed abrupt DNA replication initiation following dCas9-induced transcription blockade (Additional file 1: Figure S6f-h).
To detect MCM accumulation, we performed anti-MCM5 ChIP-qPCR on CMIP in the G1 phase following dCas9-induced transcription blockade, because ChIP-seq is not suitable for detecting minor changes at individual sites. dCas9-mediated stalling of RNA polymerase II did not alter the amount of MCM5 in the upstream or downstream non-transcribed regions of CMIP (Fig. 5f). Remarkably, MCM5 significantly accumulated in the 2-kb region upstream of the blockade sites in three biological replicates, consistent with the distribution pattern of EdU signals in the early S phase following dCas9-induced transcription blockade (Fig. 5b, c, and f; Additional file 1: Figure S6i). These findings suggest that RNA polymerase II may drive MCM along active genes to transcription-poor regions or transcription obstacles to initiate early DNA replication (Fig. 5g).
Collisions between transcription and DNA replication initiation induce genome instability
It is conceivable that the exclusion of early DNA replication initiation from transcribed regions may represent a mechanism for preserving genome integrity, especially in transcribed regions. Therefore, we investigated the impact of collisions between dysregulated early DNA replication and transcription on genome stability. We used 30 μM DRB to interrupt transcription in G1-arrested K562 cells, which allowed DNA replication initiation to occur within actively transcribed genes (Fig. 3c, d). Subsequently, the DRB treatment was ceased for 3.5 h to restart transcription and induce transcription-replication collisions, and γ-H2AX was used to mark DNA damage (Fig. 6a). Treatment with 30 μM DRB caused similar levels of damage in cells in the G1 and early S phases in comparison with untreated cells (Fig. 6b, c). After DRB treatment was ceased for 3.5 h, the intensity of the γ-H2AX signal in G1-arrested K562 cells increased slightly in comparison with that of the untreated cells (Fig. 6b, c). However, in the DRB-released early S phase cells, in which DNA replication initiation occurred in transcribed regions (Fig. 3c, d), the γ-H2AX signal intensity was significantly higher than that of the G1-phase or untreated cells (Fig. 6b, c).
R-loops are important deleterious structures produced during the collision of replication and transcription, which are sensitive to RNase H1 [8, 9]. In this context, we overexpressed RNase H1 in the DRB-treated-and-released K562 cells by plasmid transfection (Fig. 6a, d; Additional file 1: Figure S7a). Overexpression of RNase H1 caused DNA damage in asynchronized cells  and slightly, but significantly, increased the γ-H2AX signal in the G1-arrested cells before, but not after, DRB release (Fig. 6e; Additional file 1: Figure S7a). The γ-H2AX signal was elevated after DRB release in the GFP-overexpressed K562 cells in the early S phase, consistent with the findings described above. Overexpression of RNase H1 resulted in a significant decrease in the number of γ-H2AX foci in comparison with that of the GFP-overexpressed DRB-released S phase cells (Fig. 6d, e), indicating that RNase H1 alleviated DNA damage in DRB-released cells undergoing replication-transcription collisions. Collectively, these results suggest that genome instability is increased when transcription is restarted within transcribed regions with dysregulated DNA replication initiation.
It is well established that thymidine analogs can be incorporated into the genome to map DNA replication loci. CldU and IdU have been used in DNA combing assays to analyze DNA replication events at a single-molecule level using two distinct antibodies . Here, we used two different enriching strategies to detect sequentially incorporated EdU and BrdU signals and thus map replication loci via high-throughput sequencing. To ensure that EdU and BrdU-labeled replication initiation zones and adjacent elongation regions, cells were synchronized in the G1 phase using a CDK4/6 inhibitor, released at the G1/S transition stage, and sequentially labeled by EdU and BrdU for short pulses. Since no extra replication stress is present, NAIL-seq with dual-labeling of EdU and BrdU can more accurately identify early replication initiation in comparison with EdU-seq-HU (Fig. 1b, e). Moreover, NAIL-seq has a much higher resolution than Repli-seq, even when S phase cells are divided into 16 fractions for Repli-seq library generation . However, NAIL-seq with dual-labeling has a slightly lower resolution than EdU-seq-HU (Fig. 1c; Additional file 1: Figure S2b), and library preparation for the dual-labeling assay is more time-consuming. Ultimately, we combined dual-labeling of EdU and BrdU with EdU-seq-HU to precisely locate bona fide ERIZs in this study.
Both transcription and early DNA replication occur in open chromatin compartments, and each of these two pivotal cellular processes is indispensable to the other [1, 5, 7, 8]. Transcription plays a dual role in regulating early DNA replication. On the one hand, high-level transcription is coincident with robust DNA replication initiation in the same active compartments, and enhanced transcription can induce late replicating regions to replicate earlier, indicating that transcription and DNA replication initiation are collaborative [12, 14, 21, 23, 48, 52]. On the other hand, transcription can inhibit replication initiation at origins, as reported in budding yeast [53,54,55]. In this study, we propose a “transcription bulldozing” model to describe the RNA Pol II-driven MCM redistribution mechanism that spatially separates transcription and early DNA replication (Fig. 5g). We concluded that RNA polymerase II, analogous to a bulldozer in this model, pushes away MCM to suppress DNA replication initiation in transcribed regions (Fig. 4). Our results explain the observation that early DNA replication initiation is rarely detected at transcribed genes [10, 15]. According to the “transcription bulldozing” model, DNA replication initiation is predisposed to start at transcription barriers or other chromatin structural obstacles such as TSSs, G-quadruplex stretches, Z-DNA, and hairpins [10, 21,22,23, 48, 56,57,58]. Consistently, ERIZs are also enriched in native transcription barriers including G-quadruplex stretches (Additional file 1: Figure S7b). Here, we introduced a dCas9-induced transcription blockade at the transcribed regions of CMIP and GALNT10, providing another piece of evidence suggesting that RNA polymerase II may drive MCM complexes along transcribed genes (Fig. 5). When encountering the dCas9 blockade, RNA polymerase II-coupled MCM complexes were trapped upstream of dCas9-binding sites. As expected, we detected an early replication signal correlated with the trapped MCM complexes at the block (Fig. 5; Additional file 1: Figure S6). Regarding regions downstream from the block, loaded MCM double hexamers would not be expected to be cleared and might also initiate DNA replication. However, the lack of the ORC at the CMIP gene downstream of dCas9-binding sites resulted in the absence of MCM complexes and thereby DNA replication initiation (Fig. 5b; Additional file 1: Figure S6b-e).
The MCM complex experiences two phases, loading and activation, when DNA replication is initiated in eukaryotic cells . In the loading step, which occurs during the late M/G1 phase, multiple head-to-head MCM double hexamers are loaded onto chromatin by the ORC at each replication origin . These loaded MCM double hexamers are inactive and encircle the duplex DNA, allowing them to move after loading [29, 30, 35, 36]. In this context, we detected uncoupling of the ORC and MCM complexes in both cycling and G1-arrested human cells (Fig. 4b, d; Additional file 1: Figure S5a), consistent with previous reports [60, 61]. Since activated MCM double hexamers, but not the ORC, determine whether DNA replication initiation loci are functional [32, 62, 63], RNA polymerase II can effectively shape early DNA replication initiation by redistributing MCM double hexamers prior to activation. In addition, active transcription shapes the genome-wide distribution of MCM complexes in Drosophila . Moreover, in vitro studies demonstrated that MCM double-hexamers from Saccharomyces cerevisiae can slide on naked duplex DNA [29, 30], and T7 RNA polymerase can push MCM complexes off of the ends of linear dsDNA  via a process supported by studies of MCM-DNA structure [64, 65]. In the context of chromatin, in vivo ChIP-seq also showed that MCM complexes were relocated from their loading sites to adjacent downstream regions along the transcription direction upon transcription read-through in budding yeast [36, 37]. In this study, we found that genome-wide redistribution of MCM complexes was shaped by RNA polymerase II-dependent transcription in mammalian cells (Fig. 4). Moreover, RNA polymerase II can relocate MCM complexes to dCas9-binding sites within active gene bodies (Fig. 5), supporting the “transcription bulldozing” model. Furthermore, we observed that the abundance of MCM complexes on chromatin was comparable before and after transcription perturbation (Fig. 4a), implying that RNA polymerase II may push MCM complexes along chromatin rather than disassociating them from transcribed regions, in line with the observation that T7 RNA polymerase cannot disassociate budding yeast MCM complexes from circular DNA in vitro . However, we cannot fully exclude the possibility that RNA polymerase II may disassociate MCM complexes from chromatin in the G1 phase, as it has been reported that transcription can disrupt pre-RCs in budding yeast [53,54,55]. The mechanism by which RNA polymerase II redistributes MCM complexes on chromatin remains to be explored.
During the S phase, DNA replication initiation sites are prone to DNA damage and transcription is robust [9, 15, 39, 66]. Therefore, suppression of DNA replication initiation within transcribed regions is important because it allows cells to avoid collisions between DNA and RNA polymerases, thus maintaining genome stability during early DNA replication. Accordingly, we detected an increased level of genome instability when the transcription machinery encountered induced DNA replication initiation at gene bodies (Fig. 6a–c). Moreover, inducing a disordered cell cycle via overexpression of oncogenes (c-MYC or cyclin E) results in DNA replication initiation in gene bodies, which promotes chromosomal translocation and tumorigenesis . In this context, employing RNA polymerase II to exclude inactive MCM double hexamers from active gene bodies is a superbly effective strategy by which cells initiate DNA replication while maintaining robust gene transcription.
Both early DNA replication initiation and transcription occur in active chromatin compartments; however, the mechanisms through which these pivotal cellular processes coordinate remain unclear, especially in mammalian cells. Here, we report our development of the NAIL-seq method to monitor the interplay between early replication initiation and transcription, which allowed us to reveal that transcription actively shapes early DNA replication initiation in mammalian cells. Moreover, we describe the mechanism through which RNA polymerase II actively redistributes replication initiation factor MCM complexes to non-transcribed regions to avoid vulnerable early DNA replication initiation within transcribed genes. We also report that inducing collisions between transcription and DNA replication initiation within gene bodies leads to gross DNA damage. Therefore, we propose a “transcription bulldozing” model to describe the key role of transcription in preserving genomic stability during DNA replication initiation in mammalian cells.
Materials and resource table
See the Additional file 5: Table S4.
Cell culture, cell cycle synchronization, and incorporation of thymidine analogs for NAIL-seq
GM12878 and K562 (3111C0001CCC000039, National Infrastructure of Cell line resource, China) cells were cultured in RPMI1640 media supplied with 15% FBS as described previously [1, 67]. For G1 arrest, GM12878 or K562 cells were incubated with 1 μM or 5 μM palbociclib (SelleckChem) for 36 h, respectively. For EdU and BrdU dual-labeling experiments, G1-arrested GM12878 and K562 cells were washed with pre-warmed RPMI1640 and released into fresh medium for 3 or 2.5 h, respectively. Next, the cells were cultured with 10 μM EdU for 15 min, followed by a wash with pre-warmed RPMI1640, and then cultured with 50 μM BrdU for an additional 15 min. For EdU/HU treatment, G1-synchronized GM12878 or K562 cells were released into fresh medium with 10 μM EdU plus 5 mM HU for 24 h or 10 μM EdU plus 10 mM HU 12 h, respectively. The concentration of EdU was comparable with that of endogenous dTTP . EdU is efficiently incorporated into nascent DNA, and a 2-min period is sufficient for EdU to label nascent Okazaki fragments in OK-seq . Moreover, G1-released K562 cells require an additional 2.5 h to reach the G1/S transition. Therefore, the vast majority of initiation sites are subjected to EdU labeling in the EdU/HU assay. For ⍺-amanitin or DRB mediated-RNA polymerase II inhibition, G1-arrested K562 cells were released into fresh media containing the indicated concentrations of ⍺-amanitin or DRB with HU and EdU for 12 h.
Murine embryonic stem V6.5 (mES-V6.5) cells were grown in 2i media containing 15% FBS, 1 μM PD0325901 (SelleckChem), 3 μM CHIR99021 (SelleckChem), 1000 U/mL mouse LIF (Millipore), and other essential supplements as described previously . mES-V6.5 cells were cultured on 0.2% gelatin (Sigma)-treated plates and synchronized by thymidine and nocodazole as described previously . Mitotically arrested cells were shaken off, washed with pre-warmed 2i media, and cultured in fresh media for 1.5 h, followed by sequential labeling with 10 μM EdU and 50 μM BrdU for 10 min each. To reduce replication fork progression speed, M phase-synchronized mES cells were released into fresh media and incubated with 4 mM HU and 50 μM EdU for 3 h. To induce degradation of RNA polymerase II, homozygous mAID-tagged RNA polymerase II cells were first synchronized as wild-type mES cells. Next, 1 μg/mL doxycycline (Dox) was added to induce production of OsTIR1. In total, 500 μM IAA was added 1 h before the cells were released from the nocodazole treatment to induce degradation of RNA polymerase II. The nocodazole-arrested cells were then cultured in the presence of Dox, IAA, EdU, and HU for 3 h before harvesting.
Wild-type C57BL/6NCr mice (6–8 weeks old) were subjected to splenic B cell purification by following the manufacturer’s instructions included with the EasySep Mouse B cell isolation kit (STEMCELL). Both male and female mice were indiscriminately used in this study. Mice were housed and handled according to the standards set by the Peking University laboratory animal center, and all animal experiments were approved by the institutional animal care and use committee at Peking University. Naïve splenic B cells were isolated and activated with LPS (25 mg/mL; Sigma) and IL-4 (5 ng/mL; Sigma). B cells were treated with 10 mM HU and 10 μM EdU for 28 h beginning from the activation . For ⍺-amanitin treatment, primary B cells were incubated with 10 mM HU, 10 μM EdU, and ⍺-amanitin for 28 h after LPS and IL-4 were added.
Cells were harvested, washed with PBS, and subjected to genomic DNA (gDNA) extraction with Proteinase K digestion as described previously . Purified genomic DNA was subjected to the Click reaction with 100 μM Biotin-PEG-azide, 2 mM CuSO4, 4 mM THPTA, and 10 mM sodium ascorbate at 25 °C for 1 h, followed by DNA precipitation. The copper in the Click reaction hydrolyzes gDNA and fragments it into small fragments (around 200~400 bp in length), so that each fragmented nascent DNA has little chance to contain both EdU and BrdU simultaneously. The dissolved DNA was denatured at 95 °C for 5 min and phosphorylated by T4 PNK at 37 °C for 30 min.
For EdU- and BrdU-labeled samples, the denatured and modified DNA was purified by anti-BrdU (BU1/75) and Protein G beads to isolate the BrdU-labeled portion, and the supernatants were denatured again and incubated with Streptavidin C1 beads (Dynabeads) to isolate the EdU-labeled portion. For EdU/HU-labeled samples, T4 PNK-treated DNA was incubated with C1 beads alone. All of the enrichment steps were performed at room temperature for 4 h.
Isolated EdU- or BrdU-labeled ssDNA samples were ligated onto the beads with two types of bridge adapters (Bridge adapter -1 and -2, see Additional file 5: Table S4 for sequence details) simultaneously at room temperature for at least 4 h. The ligation reaction volume contained 1 μM of each Bridge adapter, 1,600 U T4 DNA ligase and 15% PEG-8000. Next, the ligated DNA was carefully washed and tagged with Illumina P5-I5 and P7-I7 sequences as described elsewhere . All libraries were subjected to 2 × 150 bp Hiseq sequencing.
ORC2 and MCM5 ChIP-seq
For the ORC2 ChIP-seq, ORC2 tagged with Avi-tag  at its N-terminal was inserted into the pMAX-GFP (Lonza) backbone to replace GFP. Bacterial BirA  was inserted into pX330 to replace Cas9. Ten million K562 cells transfected with 30 μg BirA and 30 μg Avi-ORC2 plasmids were immediately arrested at the G1 phase by 5 μM palbociclib for 36 h and then incubated with or without 10 μg/mL ⍺-amanitin for 12 h in the presence of palbociclib. For MCM5, wild-type K562 cells were treated with palbociclib for 36 h and then cultured in the presence or absence of 10 μg/mL ⍺-amanitin for 12 h with palbociclib.
Cells were fixed by 1% formaldehyde (F1635, Sigma) for 10 min at room temperature and quenched by 125 mM glycine for 5 min. Ten million fixed cells were washed twice with PBS and lysed with ice-cold NP40 lysis buffer (10 mM Tris-HCl, pH 7.5; 150 mM NaCl; 0.05% NP40) on ice for 15 min. The cell lysate was separated by 24% (w/v) sucrose in NP40 lysis buffer. The nuclei pellet was washed with 1 mM EDTA/PBS, resuspended in glycerol buffer (20 mM Tris-HCl, pH 8.0; 75 mM NaCl; 0.5 mM EDTA; 0.85 mM DTT; 50% glycerol), and lysed using nuclei lysis buffer (10 mM HEPES, pH 7.6; 1 mM DTT; 7.5 mM MgCl2; 0.2 mM EDTA; 0.3 M NaCl; 1 M urea; 1% NP40). Chromatin pellets were washed twice with 1 mM EDTA/PBS and resuspended in sonication buffer (20 mM Tris-HCl, pH 8.0; 150 mM NaCl; 2 mM EDTA; 0.1% SDS; 1% Triton X-100; 4 mM CaCl2) containing 40 U MNase for 15 min at 37 °C. The nuclease was quenched by adding 5 mM EDTA and 5 mM EGTA. The chromatin in the samples was further fragmented by sonication and separated from the soluble fraction.
For ORC2 ChIP-seq, the supernatants were pre-cleared using 40 μL protein G dynabeads for 1 h and then incubated with 40 μL Streptavidin T1 beads (Dynabeads) for 6 h at 4 °C. The beads were carefully washed with freshly made 2% SDS, high-salt buffer (500 mM NaCl; 0.1% SDS; 1% Triton X-100; 2 mM EDTA; 20 mM Tris-HCl, pH 8.0), LiCl buffer (0.25 M LiCl; 1% NP-40; 10 mM Tris-HCl; pH 8.0, 1 mM EDTA), and TE buffer. The bead-bound chromatin was eluted with elution buffer (20 mM Tris-HCl, pH 8.0; 10 mM EDTA; 1% SDS) at 65 °C overnight. For MCM5 ChIP-seq, the soluble fractions were directly subjected to anti-MCM5 (1:50) enrichment at 4 °C overnight. Next, protein G dynabeads were added and the mixture was incubated for 1 h. The beads were washed using the same procedure used for the ORC2/5 ChIP, except that the wash with 2% SDS was not performed. The chromatin on the beads was eluted twice with elution buffer at 65 °C for 15 min.
The eluate fraction was digested by RNase A and incubated with Proteinase K overnight. DNA was purified, end-repaired, and subjected to poly-dC tailing as described previously . After poly-dC tailing, DNA was tagged with a biotinylated primer, captured by Streptavidin C1 beads, and ligated with a selected bridge adapter (see Additional file 5: Table S4) overnight at room temperature. The ligated products were then washed and tagged with Illumina P5-I5 and P7-I7 sequences via PCR for 13 cycles. Sequencing was performed on the Illumina Hiseq platform (2 × 150 bp).
dCas9-mediated transcription elongation blockade
Four gRNAs (see Additional file 5: Table S4 for sequence details) binding to the non-template strand of the fourth exon of CMIP were designed in a 134-bp region and inserted into a single plasmid expressing nuclease-dead SpCas9 (dCas9) following the procedures for CRISPRi . One million K562 cells transfected with 5 μg plasmids were arrested at the G1 phase by 5 μM palbociclib for 36 h. The G1-arrested cells were either released into the early S phase or grown through the G1 phase with a further incubation with HU and EdU for 12 h. The cells were harvested, after which early DNA replication initiation and the transcription level of the CMIP gene were measured.
For MCM5 ChIP-qPCR, 10 million K562 cells were transfected with the indicated plasmids, incubated with 5 μM palbociclib immediately, and cultured for 48 h. The G1-arrested cells were harvested and ChIP-ed as described in the ChIP-seq section, except that the MNase treatment was not performed. The ChIP eluate fractions were subjected to RNA removal, decrosslinking, and purification. Next, the resuspended DNA was subjected to qPCR detection with the indicated primer sets (see Additional file 5: Table S4 for sequence details).
To block the elongation of RNA polymerase II on GALNT10, eight gRNAs (see Additional file 5: Table S4 for sequence details) targeting the second exon of GALNT10 were designed in a 160-bp region and cloned into two plasmids expressing dCas9. The two plasmids were transfected into K562 cells, simultaneously, and transfected cells were treated as described above. Next, the cells were harvested, after which early DNA replication initiation and the transcription level of GALNT10 were measured.
Flow cytometry analysis
GM12878 or K562 cells treated with or without palbociclib were labeled with 50 μM BrdU for 30 min and fixed by 4% PFA for 1 h at 4 °C. The BrdU pulse-labeled cells were denatured by 3 M HCl, incubated with anti-BrdU (× 100, BD) for 40 min, and then stained with 7-AAD (× 250, BD) for 20 min. For EdU/HU-treated samples, EdU-labeled cells were subjected to the Click reaction following the manufacturer’s instructions (Click-iT EdU Alexa Flour 488 Flow Cytometry Assay Kit). Samples were acquired on a BD FACSVerse and analyzed with FlowJo 10.4.
Cells were labeled with EdU or BrdU as described in the aforementioned sections and fixed with 4% PFA for 15 min, followed by a wash in PBS. Cells were denatured using 3 M HCl and neutralized in 0.1 M Na2B4O7. To detect the cross-reaction between anti-BrdU and EdU (Additional file 1: Figure S1d), denatured cells were Click-reacted with TAMRA-azide-biotin or DMSO, blocked with BSA, and incubated with anti-BrdU (BU1/75), followed by incubation with the indicated secondary antibody. For cell cycle progression detection (Additional file 1: Figure S1b), denatured samples were blocked with 5% BSA, incubated with an antibody against BrdU (BU1/75) for 1 h followed by a thorough wash, and then incubated with the indicated secondary antibody. For γ-H2AX detection (Additional file 1: Figure S4c), G1-arrested K562 cells were released into fresh medium containing 10 mM HU plus 10 μM EdU with or without amanitin treatment for 12 h. The fixed cells were Click-reacted with Azide-Alexa-555 and then incubated with antibodies against γ-H2AX and RNA polymerase II simultaneously. For the early S phase samples used in the DRB removal assay, G1-arrested K562 cells were released into fresh medium supplied with 10 mM HU, 50 μM BrdU, and 30 μM DRB for 12 h. Next, the cells were washed and released into fresh medium containing 10 mM HU for 3.5 h. For the G1 phase samples used in the DRB removal assay, G1-arrested K562 cells were cultured in the presence of 30 μM DRB for 12 h and released into fresh medium containing 5 μM palbociclib for an additional 3.5 h. Next, the fixed cells were denatured and incubated with antibodies against BrdU or γ-H2AX. For the RNase H1 overexpression assay, cells were transfected with plasmids containing GFP or GFP-RNase H1 before cell arrest was induced via palbociclib. Images were acquired using an LSM 710 NLO & DuoScan System (× 63, 1.4 NA or × 40, 1.2 NA) and analyzed using ImageJ according to the instructions included with the software.
Quantification and statistical analysis
Sequencing data analysis
NAIL-seq data processing
R1 and R2 reads were stitched using PEAR  and the bridge adapter sequences were trimmed with custom scripts. The stitched reads were aligned to the human genome (assembly hg19) or mouse genome (assembly mm10) by BWA-MEM (-k 20). PCR duplicates were distinguished by random molecular barcodes (RMBs) and only one was kept. Reads with mapping quality (MAPQ) ≥ 30 were kept for further analysis. We used samtools to convert .sam files into .bam files.
ChIP-seq data processing
First, we trimmed the adapter sequences in the reads using cutadapt (-g NGGGGGGGGG -g GGGGG --times=2 --minimum-length=70 --error-rate=0.2 -O 3) . Next, we aligned reads to the human genome (assembly hg19) using BWA-MEM (-k 20). For MCM ChIP-seq, we first generated RPKM to normalize the total reads by converting alignment files in .bam to .BigWig format using bamCoverage from deeptools  with parameters of -bs 1000 and --normalizeUsing RPKM. To compare with the control, the fold change of the RPKM signal between the MCM ChIP-seq and the input control was generated for further analysis. For ORC2 ChIP-seq, peaks were identified by MACS 1.4.2 using -p 1e-5 --nolambda --nomodel --keep-dup=all . The peaks within 500 bp were merged using BEDTools 2.27.0 . Peaks within 500 bp from two biological replicates without ⍺-amanitin treatment were defined as reproducible ORC2 peaks and used for further analysis. Genome-wide ORC2 tracks in .BigWig format were generated by bamCoverage from deeptools with -bs 50 and –normalizeUsing RPKM.
Identifying early replication initiation zones (ERIZs) from NAIL-seq
We developed a pipeline named “RepFind” for NAIL-seq analysis, which is described in the following section.
Defining E-B peaks from EdU/BrdU-labeled NAIL-seq samples
Both EdU and BrdU reads were normalized by reads per million (RPM) in 5-kb bins tiling the whole genome and marked as E and B, respectively. ΔEB was defined as: ΔEB=E-B. Bins, with ΔEB > 0.3, were taken as seeds and the neighbor bins with ΔEB > 0 were merged with the seeds to mark EdU-rich regions. The E-B peaks were the EdU-rich regions whose lengths were > 20 kb and fell in EdU peaks (defined from EdU reads using SICER with options -w 5000 -g 3). Peaks that occurred in at least two biological replicates were kept as bona fide E-B peaks.
Defining EdU/HU peaks
The peak calling procedure was modified as reported previously . Peaks were called by MACS (version 1.4.2) (-p 1e-5 --nolambda --nomodel --keep-dup=all). Next, we filtered the peaks with ≥ 400-fold enrichment against a random Poisson distribution generated from MACS14. Neighbor peaks were merged if their interval was shorter than 10 kb for K562 or 20 kb for GM12878. The merged peaks with < 10 kb size were discarded. Peaks in ChrY, mitochondria DNA, and blacklist regions were omitted.
The ERIZs were defined as EdU/HU peaks that overlapped with E-B peaks, while non-ERIZs were defined as EdU/HU peaks that did not overlap with E-B peaks. The replication timing value was obtained from wavelet-smoothed signals from six fractions of the ENCODE Repli-seq profile. The mean value of replication timing for each region was calculated by bwtool summary. The distribution was plotted by ggplot2 histogram. Domains with a replication timing value > 0.5 were defined as early replication domains, and other domains were classified as late replication domains (Additional file 1: Figure S2c).
Comparing OK-seq and SNS-seq with NAIL-seq
OK-seq data from K562 cells were analyzed using the BAMscale pipeline with default settings (https://github.com/ncbi/BAMscale/wiki/Detailed-Use:-OKseq-RFD-(Replication-Fork-Directionality)-Track-Generation). Replication initiation zones (IZs) were identified by using OKseq_switches.R with default settings (https://github.com/ncbi/BAMscale/wiki/Finding-OK-seq-strand-switched-from-the-RFD-track). Replication origins identified by SNS-seq were downloaded from GSE46189 (GSE46189_Ori-Peak.bed.gz). The overlaps between ERIZs/OK-seq and ERIZs/SNS-seq were calculated using intersect in Bedtools.
Visualization of sequencing data
To generate genome tracks, we converted the .bam files into .BigWig files. For NAIL-seq and ORC2 ChIP-seq, files were generated by bamCoverage (in deepTools) with the parameter --binSize 1000 --normalizeUsing RPKM. For MCM5 ChIP-seq, the fold change value of the ChIP-ed samples over the input was generated by bamCompare with parameter --binSize 1000, --normalizeUsing RPKM –operation ratio. The data were presented by IGV 2.4.14 .
All heatmaps were generated by the R package “EnrichedHeatmap” . For the non-TR heatmaps (see the “GRO-seq analysis” section for region definition), the 20–100-kb non-TRs were centered at the midpoint and ordered by increasing width, and the ±100-kb regions from the midpoint were used for display. Each region was binned with a 1-kb window filled with the mean value of each signal. For early replication or ORC2 ChIP-seq, the signal was generated by RPKM. Specially for MCM5 ChIP-seq, the z-score-transformed fold change value was displayed. For the transcribed gene-related heatmaps, ERIZ-flanked transcribed regions (defined by TR in the “GRO-seq” analysis) larger than 50 kb were ranked by gene width, with the smallest at the top. For display, all transcribed regions were scaled to the same width and aligned at both TSS and TTS. Additional 50-kb regions upstream or downstream of the TSS or TTS, respectively, were shown. The signal processing was performed in the same manner as that performed for the non-TR-related heatmaps.
A/B compartment analysis
GM12878 and K562 HiC data (.hic file) were downloaded from GSE63525 . Next, eigenvector from juicer  (KR BP 100000) was used to identify A/B compartments in GM12878 or K562 cells with 100-kb bins. The arbitrarily assigned “+” or “–” from eigenvector was corrected by defining “+” as an eigenvector positively correlated with an active histone marker (such as H3K27ac). The A compartments were defined as regions with a positive corrected value (“+”) for the eigenvector.
Prediction of ERIZ locations using epigenetic features
To predict ERIZ locations using epigenetic features within A compartments, we built a logistic regression model based on the peaks of the ENCODE ChIP-seq data (Fig. 2b; Additional file 1: Figure S3a). First, we converted the location of each ERIZ into a binary format using ChromHMM with 50-kb windows within an A compartment, where the window overlapping with the ERIZ is designated as “1” (“ERIZ-1”), and all others are designated as “0” (“ERIZ-0”) . Second, to generate the epigenetic predictor matrix for ERIZs, we defined the enrichment score of each 50-kb window using the log odds ratio of the observed and expected occurrence of ChIP-seq peaks, where the observed occurrence was defined as the proportion of 1-kb bins that overlapped with ChIP-seq peaks. The expected occurrence was estimated as the mean occurrence of the random peaks shuffled within the same chromosome. Third, we subsampled bins from the “ERIZ-0” bins to balance the sample size of the “ERIZ-1” bins. Finally, we built the logistic regression model using “glm” in R. The coefficients were used to evaluate the contribution to the prediction of ERIZ location.
K562 and GM12878 GRO-seq data from a previous report were acquired for re-analysis . For activated mouse B cells and mES V6.5 cells, GRO-seq data were acquired from published reports [15, 69]. Adapter-trimmed reads were aligned using bowtie (-l 25 -v 1 -k 1 -m 1 -S -q --best) . Reads mapped to rRNA were omitted. For K562 and GM12878 cells, the read density (reads per kilobase, RPK) normalized to 10 million sequencing reads in the gene bodies or promoter regions was calculated via analyzeRepeats.pl from Homer (analyzeRepeats.pl rna hg19 -count genes -condenseGenes -strand + -norm 1e7) . For activated mouse B cells and mES V6.5 cells, the read density was obtained from Homer (analyzeRepeats.pl rna mm10 -count genes -condenseGenes -strand + -norm 1e6).
Defining active genes
Active genes in K562 or GM12878 cells were defined as those with read density > 0 at the promoter region and RPK > 4 at the gene body from the GRO-seq data. Silent genes were defined as those with no read at the promoter region and RPK ≤ 1 at the gene body. Mouse active genes in activated B cells and mES V6.5 cells were defined as those with RPK > 1. Silent genes were genes with zero RPK from the GRO-seq data.
Defining non-transcribed regions (non-TRs) and transcribed regions (TRs)
The non-transcribed regions used in heatmaps were defined as interval regions between two active genes within A compartments via bedtools subtract. Non-transcribed regions with a size of 20–100 kb that overlapped with ERIZs were used for the non-TR-related heatmaps. The active genes flanked by ERIZ-occupied non-TRs were defined as transcribed regions (TRs) for display in heatmaps (Fig. 4f–h).
Comparing signal enrichment in non-transcribed regions and transcribed regions
The non-transcribed regions were further trimmed by 200 bp at the head and tail to exclude the TSS and TTS regions (Figs. 3b, d, f and 4c, e; Additional file 1: Figure S4f). The flanked transcribed regions were also trimmed as active gene body regions (TSS + 200 bp, TTS − 200 bp). The read density was calculated as follows: for each non-transcribed region, reads in the region were counted and then normalized by region length. Read counts in transcribed regions upstream or downstream of non-transcribed regions were summed and normalized by the summed length of the transcribed regions. In order to compare the signal enrichment of non-transcribed and transcribed regions based on NAIL-seq and ORC2 ChIP-seq, the log2 fold change of the ratio of the read density in the non-transcribed region to that of the transcribed region was calculated to allow quantification of differences. For MCM ChIP-seq, the fold enrichment of the read density relative to that of the input control was calculated and used to calculate the log2 fold change in the non-transcribed regions and transcribed regions.
Analysis of dCas9 as a transcription barrier
To compare the read density adjacent to the dCas9 binding sites, the RPK ratio was defined as follows: the read count per kilobase (RPK) was calculated in each bin and then divided by the mean RPK of the B compartments at the same chromosome to normalize the background of each sample. For the CMIP locus, the ± 5-kb region centered on the dCas9 binding site was tiled in 1-kb bins (with a 200-bp sliding) for calculation of the RPK ratio (Fig. 5b, d). A 2-kb bin upstream or downstream of the dCas9 binding site was used to calculate the RPK ratios shown in Fig. 5c, e. For the GALNT10 locus, the ± 6-kb region centered on the dCas9 binding sites was tiled in 3-kb bins (with a 300-bp sliding) (Additional file 1: Figure S6h).
We performed Student’s t test to test the statistical significance between RPK ratios for dCas9 with gRNA-CMIP and dCas9 with scrambled gRNA. The t-test was applied to assess the difference between the dCas9-mediated transcription blockade and control groups in the MCM5 ChIP-qPCR assay. The t-test was also used to analyze the distribution of γ-H2AX intensity with or without DRB treatment. The Wilcoxon rank-sum test was performed to analyze the differences between the log fold change NAIL-seq signal and ORC2 or MCM5 ChIP-seq signals in non-transcribed regions in comparison with transcribed regions among the groups subjected to different treatments with transcription inhibitors.
Availability of data and materials
All sequencing data presented in the present study were deposited at NODE (National Omics Data Encyclopedia) under the accession number of OEP000658  and were also deposited in NCBI Gene Expression Omnibus (GEO) database under accession GSE174680 . All the other data generated in the present study can be found in the manuscript and additional files. All public data used in this study have been listed and well referenced in Additional file 5: Table S4, including Repli-seq of K562: GSM923448 , SNS-seq of K562: GSE46189 , ORC2 ChIP-seq from K562: GSE70165 , Repli-seq of GM12878: GSM923451 , MCM3 ChIP-seq from K562: GSE92216 , MCM5 ChIP-seq from K562: GSE127403 , HiC from K562 and GM12878: GSE63525 , GRO-seq for GM12878: GSM1480326 , GRO-seq for K562: GSM1480325 , G4 ChIP-seq from K562: GSM2876090 , EdU-seq-HU in mouse splenic B cells: GSE116321 , nsRNA-seq for mouse splenic B cells: GSM3227964 , and GRO-seq for mESC: GSM1665566  in GEO database and OK-seq of K562: PRJEB25180  in European Nucleotide Archive database. Scripts of RepFind pipeline for NAIL-seq in this study are available in GitHub repository under GNU general public license Version 3  and in Zenodo . All cell lines used in this study have been authenticated and are available upon request.
Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159(7):1665–80. https://doi.org/10.1016/j.cell.2014.11.021.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485(7398):376–80. https://doi.org/10.1038/nature11082.
Hnisz D, Day DS, Young RA. Insulated neighborhoods: structural and functional units of mammalian gene control. Cell. 2016;167(5):1188–200. https://doi.org/10.1016/j.cell.2016.10.024.
Kieffer-Kwon KR, Tang Z, Mathe E, Qian J, Sung MH, Li G, et al. Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation. Cell. 2013;155(7):1507–20. https://doi.org/10.1016/j.cell.2013.11.039.
Pope BD, Ryba T, Dileep V, Yue F, Wu W, Denas O, et al. Topologically associating domains are stable units of replication-timing regulation. Nature. 2014;515(7527):402–5. https://doi.org/10.1038/nature13986.
Fragkos M, Ganier O, Coulombe P, Méchali M. DNA replication origin activation in space and time. Nat Rev Mol Cell Biol. 2015;16(6):360–74. https://doi.org/10.1038/nrm4002.
Marchal C, Sima J, Gilbert DM. Control of DNA replication timing in the 3D genome. Nat Rev Mol Cell Biol. 2019;20(12):721–37. https://doi.org/10.1038/s41580-019-0162-y.
Garcia-Muse T, Aguilera A. Transcription-replication conflicts: how they occur and how they are resolved. Nat Rev Mol Cell Biol. 2016;17(9):553–63. https://doi.org/10.1038/nrm.2016.88.
Hamperl S, Bocek MJ, Saldivar JC, Swigut T, Cimprich KA. Transcription-replication conflict orientation modulates r-loop levels and activates distinct DNA damage responses. Cell. 2017;170(4):774–86 e719. https://doi.org/10.1016/j.cell.2017.07.043.
Petryk N, Kahli M, d'Aubenton-Carafa Y, Jaszczyszyn Y, Shen Y, Silvain M, et al. Replication landscape of the human genome. Nat Commun. 2016;7(1):10208. https://doi.org/10.1038/ncomms10208.
Petryk N, Dalby M, Wenger A, Stromme CB, Strandsby A, Andersson R, et al. MCM2 promotes symmetric inheritance of modified histones during DNA replication. Science. 2018;361(6409):1389–92. https://doi.org/10.1126/science.aau0294.
Chen YH, Keegan S, Kahli M, Tonzi P, Fenyo D, Huang TT, et al. Transcription shapes DNA replication initiation and termination in human cells. Nat Struct Mol Biol. 2019;26(1):67–77. https://doi.org/10.1038/s41594-018-0171-0.
Macheret M, Halazonetis TD. Intragenic origins due to short G1 phases underlie oncogene-induced DNA replication stress. Nature. 2018;555(7694):112–6. https://doi.org/10.1038/nature25507.
Blin M, Le Tallec B, Nähse V, Schmidt M, Brossas C, Millot GA, et al. Transcription-dependent regulation of replication dynamics modulates genome stability. Nat Struct Mol Biol. 2019;26(1):58–66. https://doi.org/10.1038/s41594-018-0170-1.
Tubbs A, Sridharan S, van Wietmarschen N, Maman Y, Callen E, Stanlie A, et al. Dual roles of poly(dA:dT) tracts in replication initiation and fork collapse. Cell. 2018;174(5):1127–42 e1119. https://doi.org/10.1016/j.cell.2018.07.011.
Ekundayo B, Bleichert F. Origins of DNA replication. PLoS Genet. 2019;15(9):e1008320. https://doi.org/10.1371/journal.pgen.1008320.
Prioleau MN, MacAlpine DM. DNA replication origins-where do we begin? Genes Dev. 2016;30(15):1683–97. https://doi.org/10.1101/gad.285114.116.
MacAlpine DM. ORChestrating the human DNA replication program. Proc Natl Acad Sci U S A. 2016;113(33):9136–8. https://doi.org/10.1073/pnas.1610336113.
Hansen RS, Thomas S, Sandstrom R, Canfield TK, Thurman RE, Weaver M, et al. Sequencing newly replicated DNA reveals widespread plasticity in human replication timing. Proc Natl Acad Sci U S A. 2010;107(1):139–44. https://doi.org/10.1073/pnas.0912402107.
Zhao PA, Sasaki T, Gilbert DM. High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells. Genome Biol. 2020;21(1):76. https://doi.org/10.1186/s13059-020-01983-8.
Cayrou C, Coulombe P, Vigneron A, Stanojcic S, Ganier O, Peiffer I, et al. Genome-scale analysis of metazoan replication origins reveals their organization in specific but flexible sites defined by conserved features. Genome Res. 2011;21(9):1438–49. https://doi.org/10.1101/gr.121830.111.
Besnard E, Babled A, Lapasset L, Milhavet O, Parrinello H, Dantec C, et al. Unraveling cell type-specific and reprogrammable human replication origin signatures associated with G-quadruplex consensus motifs. Nat Struct Mol Biol. 2012;19(8):837–44. https://doi.org/10.1038/nsmb.2339.
Picard F, Cadoret JC, Audit B, Arneodo A, Alberti A, Battail C, et al. The spatiotemporal program of DNA replication is associated with specific combinations of chromatin marks in human cells. PLoS Genet. 2014;10(5):e1004282. https://doi.org/10.1371/journal.pgen.1004282.
Wu X, Kabalane H, Kahli M, Petryk N, Laperrousaz B, Jaszczyszyn Y, et al. Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions. Nucleic Acids Res. 2018;46(19):10157–72. https://doi.org/10.1093/nar/gky797.
O'Donnell M, Langston L, Stillman B. Principles and concepts of DNA replication in bacteria, archaea, and eukarya. Cold Spring Harb Perspect Biol. 2013;5(7). https://doi.org/10.1101/cshperspect.a010108.
Douglas ME, Ali FA, Costa A, Diffley JFX. The mechanism of eukaryotic CMG helicase activation. Nature. 2018;555(7695):265–8. https://doi.org/10.1038/nature25787.
Yeeles JT, Deegan TD, Janska A, Early A, Diffley JF. Regulated eukaryotic DNA replication origin firing with purified proteins. Nature. 2015;519(7544):431–5. https://doi.org/10.1038/nature14285.
Ilves I, Petojevic T, Pesavento JJ, Botchan MR. Activation of the MCM2-7 helicase by association with Cdc45 and GINS proteins. Mol Cell. 2010;37(2):247–58. https://doi.org/10.1016/j.molcel.2009.12.030.
Evrin C, Clarke P, Zech J, Lurz R, Sun J, Uhle S, et al. A double-hexameric MCM2-7 complex is loaded onto origin DNA during licensing of eukaryotic DNA replication. Proc Natl Acad Sci U S A. 2009;106(48):20240–5. https://doi.org/10.1073/pnas.0911500106.
Remus D, Beuron F, Tolun G, Griffith JD, Morris EP, Diffley JF. Concerted loading of Mcm2-7 double hexamers around DNA during DNA replication origin licensing. Cell. 2009;139(4):719–30. https://doi.org/10.1016/j.cell.2009.10.015.
Miotto B, Ji Z, Struhl K. Selectivity of ORC binding sites and the relation to replication timing, fragile sites, and deletions in cancers. Proc Natl Acad Sci U S A. 2016;113(33):E4810–9. https://doi.org/10.1073/pnas.1609060113.
Hyrien O. How MCM loading and spreading specify eukaryotic DNA replication initiation sites. F1000Res. 2016;5:2063. https://doi.org/10.12688/f1000research.9008.1.
Yankulov K, Todorov I, Romanowski P, Licatalosi D, Cilli K, McCracken S, et al. MCM proteins are associated with RNA polymerase II holoenzyme. Mol Cell Biol. 1999;19(9):6154–63. https://doi.org/10.1128/MCB.19.9.6154.
Snyder M, Huang XY, Zhang JJ. The minichromosome maintenance proteins 2-7 (MCM2-7) are necessary for RNA polymerase II (Pol II)-mediated transcription. J Biol Chem. 2009;284(20):13466–72. https://doi.org/10.1074/jbc.M809471200.
Powell SK, MacAlpine HK, Prinz JA, Li Y, Belsky JA, MacAlpine DM. Dynamic loading and redistribution of the Mcm2-7 helicase complex through the cell cycle. EMBO J. 2015;34(4):531–43. https://doi.org/10.15252/embj.201488307.
Gros J, Kumar C, Lynch G, Yadav T, Whitehouse I, Remus D. Post-licensing specification of eukaryotic replication origins by facilitated Mcm2-7 sliding along DNA. Mol Cell. 2015;60(5):797–807. https://doi.org/10.1016/j.molcel.2015.10.022.
Foss EJ, Gatbonton-Schwager T, Thiesen AH, Taylor E, Soriano R, Lao U, et al. Sir2 suppresses transcription-mediated displacement of Mcm2-7 replicative helicases at the ribosomal DNA repeats. PLoS Genet. 2019;15(5):e1008138. https://doi.org/10.1371/journal.pgen.1008138.
Kumar V, Alt FW, Frock RL. PAXX and XLF DNA repair factors are functionally redundant in joining DNA breaks in a G1-arrested progenitor B-cell line. Proc Natl Acad Sci U S A. 2016;113(38):10619–24. https://doi.org/10.1073/pnas.1611882113.
Barlow JH, Faryabi RB, Callen E, Wong N, Malhowski A, Chen HT, et al. Identification of early replicating fragile sites that contribute to genome instability. Cell. 2013;152(3):620–32. https://doi.org/10.1016/j.cell.2013.01.006.
Ge XQ, Jackson DA, Blow JJ. Dormant origins licensed by excess Mcm2-7 are required for human cells to survive replicative stress. Genes Dev. 2007;21(24):3331–41. https://doi.org/10.1101/gad.457807.
Hoffman EA, McCulley A, Haarer B, Arnak R, Feng W. Break-seq reveals hydroxyurea-induced chromosome fragility as a result of unscheduled conflict between DNA replication and transcription. Genome Res. 2015;25(3):402–12. https://doi.org/10.1101/gr.180497.114.
Hu J, Sun L, Shen F, Chen Y, Hua Y, Liu Y, et al. The intra-S phase checkpoint targets Dna2 to prevent stalled replication forks from reversing. Cell. 2012;149(6):1221–32. https://doi.org/10.1016/j.cell.2012.04.030.
Ryba T, Hiratani I, Lu J, Itoh M, Kulik M, Zhang J, et al. Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types. Genome Res. 2010;20(6):761–70. https://doi.org/10.1101/gr.099655.109.
Long H, Zhang L, Lv M, Wen Z, Zhang W, Chen X, et al. H2A.Z facilitates licensing and activation of early replication origins. Nature. 2020;577(7791):576–81. https://doi.org/10.1038/s41586-019-1877-9.
Core LJ, Martins AL, Danko CG, Waters CT, Siepel A, Lis JT. Analysis of nascent RNA identifies a unified architecture of initiation regions at mammalian promoters and enhancers. Nat Genet. 2014;46(12):1311–20. https://doi.org/10.1038/ng.3142.
Ballabeni A, Park IH, Zhao R, Wang W, Lerou PH, Daley GQ, et al. Cell cycle adaptations of embryonic stem cells. Proc Natl Acad Sci U S A. 2011;108(48):19252–7. https://doi.org/10.1073/pnas.1116794108.
Bushnell DA, Cramer P, Kornberg RD. Structural basis of transcription: alpha-amanitin-RNA polymerase II cocrystal at 2.8 A resolution. Proc Natl Acad Sci U S A. 2002;99(3):1218–22. https://doi.org/10.1073/pnas.251664698.
Langley AR, Graf S, Smith JC, Krude T. Genome-wide identification and characterisation of human DNA replication origins by initiation site sequencing (ini-seq). Nucleic Acids Res. 2016;44(21):10230–47. https://doi.org/10.1093/nar/gkw760.
Jiang Y, Huang J, Lun K, Li B, Zheng H, Li Y, et al. Genome-wide analyses of chromatin interactions after the loss of Pol I, Pol II, and Pol III. Genome Biol. 2020;21(1):158. https://doi.org/10.1186/s13059-020-02067-3.
Shen W, Sun H, De Hoyos CL, Bailey JK, Liang XH, Crooke ST. Dynamic nucleoplasmic and nucleolar localization of mammalian RNase H1 in response to RNAP I transcriptional R-loops. Nucleic Acids Res. 2017;45(18):10672–92. https://doi.org/10.1093/nar/gkx710.
Anglana M, Apiou F, Bensimon A, Debatisse M. Dynamics of DNA replication in mammalian somatic cells: nucleotide pool modulates origin choice and interorigin spacing. Cell. 2003;114(3):385–94. https://doi.org/10.1016/S0092-8674(03)00569-5.
Pourkarimi E, Bellush JM, Whitehouse I. Spatiotemporal coupling and decoupling of gene transcription with DNA replication origins during embryogenesis in C. elegans. Elife. 2016;5:e21728. https://doi.org/10.7554/eLife.21728.
Lõoke M, Reimand J, Sedman T, Sedman J, Järvinen L, Värv S, et al. Relicensing of transcriptionally inactivated replication origins in budding yeast. J Biol Chem. 2010;285(51):40004–11. https://doi.org/10.1074/jbc.M110.148924.
Mori S, Shirahige K. Perturbation of the activity of replication origin by meiosis-specific transcription. J Biol Chem. 2007;282(7):4447–52. https://doi.org/10.1074/jbc.M609671200.
Tanaka S, Halter D, Livingstone-Zatchej M, Reszel B, Thoma F. Transcription through the yeast origin of replication ARS1 ends at the ABFI binding site and affects extrachromosomal maintenance of minichromosomes. Nucleic Acids Res. 1994;22(19):3904–10. https://doi.org/10.1093/nar/22.19.3904.
Andrade-Lima LC, Veloso A, Ljungman M. Transcription Blockage Leads to New Beginnings. Biomolecules. 2015;5(3):1600–17. https://doi.org/10.3390/biom5031600.
Prorok P, Artufel M, Aze A, Coulombe P, Peiffer I, Lacroix L, et al. Involvement of G-quadruplex regions in mammalian replication origin activity. Nat Commun. 2019;10(1):3274. https://doi.org/10.1038/s41467-019-11104-0.
Valton AL, Hassan-Zadeh V, Lema I, Boggetto N, Alberti P, Saintome C, et al. G4 motifs affect origin positioning and efficiency in two vertebrate replicators. EMBO J. 2014;33(7):732–46. https://doi.org/10.1002/embj.201387506.
Bell SP, Labib K. Chromosome duplication in Saccharomyces cerevisiae. Genetics. 2016;203(3):1027–67. https://doi.org/10.1534/genetics.115.186452.
Kirstein N, Buschle A, Wu X, Krebs S, Blum H, Hammerschmidt W, et al. Active transcription regulates ORC/MCM distribution whereas replication timing correlates with ORC density in human cells. bioRxiv. 2019:778423. https://doi.org/10.1101/778423.
Takahashi TS, Wigley DB, Walter JC. Pumps, paradoxes and ploughshares: mechanism of the MCM2-7 DNA helicase. Trends Biochem Sci. 2005;30(8):437–44. https://doi.org/10.1016/j.tibs.2005.06.007.
Hua XH, Newport J. Identification of a preinitiation step in DNA replication that is independent of origin recognition complex and cdc6, but dependent on cdk2. J Cell Biol. 1998;140(2):271–81. https://doi.org/10.1083/jcb.140.2.271.
Rowles A, Tada S, Blow JJ. Changes in association of the Xenopus origin recognition complex with chromatin on licensing of replication origins. J Cell Sci. 1999;112(Pt 12):2011–8. https://doi.org/10.1242/jcs.112.12.2011.
Noguchi Y, Yuan Z, Bai L, Schneider S, Zhao G, Stillman B, et al. Cryo-EM structure of Mcm2-7 double hexamer on DNA suggests a lagging-strand DNA extrusion model. Proc Natl Acad Sci U S A. 2017;114(45):E9529–38. https://doi.org/10.1073/pnas.1712537114.
Li N, Zhai Y, Zhang Y, Li W, Yang M, Lei J, et al. Structure of the eukaryotic MCM complex at 3.8 Å. Nature. 2015;524(7564):186–91. https://doi.org/10.1038/nature14685.
Liu Y, Chen S, Wang S, Soares F, Fischer M, Meng F, et al. Transcriptional landscape of the human cell cycle. Proc Natl Acad Sci U S A. 2017;114(13):3473–8. https://doi.org/10.1073/pnas.1617636114.
Chiarle R, Zhang Y, Frock RL, Lewis SM, Molinie B, Ho YJ, et al. Genome-wide translocation sequencing reveals mechanisms of chromosome breaks and rearrangements in B cells. Cell. 2011;147(1):107–19. https://doi.org/10.1016/j.cell.2011.07.049.
Chen P, Liu Z, Liu S, Xie Z, Aimiuwu J, Pang J, et al. A LC-MS/MS method for the analysis of intracellular nucleoside triphosphate levels. Pharm Res. 2009;26(6):1504–15. https://doi.org/10.1007/s11095-009-9863-9.
Sigova AA, Abraham BJ, Ji X, Molinie B, Hannett NM, Guo YE, et al. Transcription factor trapping by RNA in gene regulatory elements. Science. 2015;350(6263):978–81. https://doi.org/10.1126/science.aad3346.
Yin J, Liu M, Liu Y, Wu J, Gan T, Zhang W, et al. Optimizing genome editing strategy by primer-extension-mediated sequencing. Cell Discov. 2019;5(1):18. https://doi.org/10.1038/s41421-019-0088-8.
Fairhead M, Howarth M. Site-specific biotinylation of purified proteins using BirA. Methods Mol Biol. 2015;1266:171–84. https://doi.org/10.1007/978-1-4939-2272-7_12.
Peng X, Wu J, Brunmeir R, Kim SY, Zhang Q, Ding C, et al. TELP, a sensitive and versatile library construction method for next-generation sequencing. Nucleic Acids Res. 2015;43(6):e35. https://doi.org/10.1093/nar/gku818.
Qi LS, Larson MH, Gilbert LA, Doudna JA, Weissman JS, Arkin AP, et al. Repurposing CRISPR as an RNA-guided platform for sequence-specific control of gene expression. Cell. 2013;152(5):1173–83. https://doi.org/10.1016/j.cell.2013.02.022.
Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics. 2014;30(5):614–20. https://doi.org/10.1093/bioinformatics/btt593.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17:10–2.
Ramírez F, Dündar F, Diehl S, Grüning BA. Manke T: deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42(W1):W187–91. https://doi.org/10.1093/nar/gku365.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based Analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137. https://doi.org/10.1186/gb-2008-9-9-r137.
Quinlan AR. BEDTools: the Swiss-army tool for genome feature analysis. Curr Protoc Bioinformatics. 2014;47:11.12.11–34 editoral board, Andreas D Baxevanis [et al].
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6. https://doi.org/10.1038/nbt.1754.
Gu Z, Eils R, Schlesner M, Ishaque N. EnrichedHeatmap: an R/Bioconductor package for comprehensive visualization of genomic signal associations. BMC Genomics. 2018;19(1):234. https://doi.org/10.1186/s12864-018-4625-x.
Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Methods. 2012;9(3):215–6. https://doi.org/10.1038/nmeth.1906.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25. https://doi.org/10.1186/gb-2009-10-3-r25.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Molecular Cell. 2010;38(4):576–89. https://doi.org/10.1016/j.molcel.2010.05.004.
Ai C: Transcription orchestrates early DNA replication initiation in mammalian cells. OEP000658. National Omics Data Encyclopedia (NODE). 2019. https://www.biosino.org/node/project/detail/OEP000658.
Liu Y, Ai C, Gan T, Hu J: Transcription shapes DNA replication initiation to preserve genome integrity. GSE174680. Gene Expression Omnibus (GEO). 2021. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE174680.
Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74. https://doi.org/10.1038/nature11247.
Mao SQ, Ghanbarian AT, Spiegel J, Martínez Cuesta S, Beraldi D, Di Antonio M, et al. DNA G-quadruplex structures mold the DNA methylome. Nat Struct Mol Biol. 2018;25(10):951–7. https://doi.org/10.1038/s41594-018-0131-8.
Liu Y, Ai C, Gan T, Hu J: Transcription shapes DNA replication initiation to preserve genome integrity. Github. 2021. https://github.com/JiazhiHuLab/RepFind_for_Nailseq.
Liu Y, Ai C, Gan T, Hu J: Transcription shapes DNA replication initiation to preserve genome integrity. Zenodo. 2021. https://zenodo.org/record/4784707#.YKzVuahKjIU.
We acknowledge the members of the Hu laboratory for their helpful discussions. The GFP-RNaseH1 plasmid was a kind gift from Dr. Sijie Liu and Dr. Daochun Kong. We thank the Flow Cytometry Core at the National Center for Protein Sciences, Peking University, particularly Liying Du, for her technical help. We also thank Dr. Siying Qin at the Core Facilities of the School of Life Sciences, Peking University for assistance with optical/confocal imaging.
The review history is available as Additional file 6.
Peer review information
Anahita Bishop was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
This work was supported by the National Key R&D Program of China (2017YFA0506700 to J.H., 2017YFA0506600 to X.J., 2019YFA0508900 and 2016YFA0500700 to N.G.) and the NSFC grant (31771485 to J.H., 31871309 to X.J., 31830048 to Q.L., 31725007, and 31630087 to N.G.). N.G., Q.L., X.J., and J.H. are investigators at the PKU-TSU Center for Life Sciences. J.H. is also a Bayer Investigator and Y.L. is supported by the Boehringer Ingelheim-Peking University Postdoctoral Program.
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Experimental procedure for NAIL-seq. Figure S2. Comparison among NAIL-seq, Repli-seq, OK-seq, SNS-seq, and EdU/HU-seq in K562 or mouse primary splenic B cells. Figure S3. The landscape of ERIZs. Figure S4. Transcription disturbance relocates ERIZs in K562 cells, mSBCs, and mESCs. Figure S5. Transcription inhibition does not lead to ORC2 relocation. Figure S6. dCas9-mediated transcription blockade induces MCM5 accumulation and replication initiation. Figure S7. Collision of transcription and DNA replication initiation induces genome instability. Figure S8. The full uncropped pictures for Western blots.
The locations of ERIZs in GM12878 cells, aligned to hg19.
The locations of ERIZs in K562 cells, aligned to hg19.
The locations of ERIZs in mESCs, aligned to mm10.
Key resource table.
About this article
Cite this article
Liu, Y., Ai, C., Gan, T. et al. Transcription shapes DNA replication initiation to preserve genome integrity. Genome Biol 22, 176 (2021). https://doi.org/10.1186/s13059-021-02390-3
- DNA replication initiation
- MCM redistribution
- Transcription-replication initiation collision
- Genome instability
- DNA damage