Cell-Cycle–Dependent Chromatin Dynamics at Replication Origins

Origins of DNA replication are specified by the ordered recruitment of replication factors in a cell-cycle–dependent manner. The assembly of the pre-replicative complex in G1 and the pre-initiation complex prior to activation in S phase are well characterized; however, the interplay between the assembly of these complexes and the local chromatin environment is less well understood. To investigate the dynamic changes in chromatin organization at and surrounding replication origins, we used micrococcal nuclease (MNase) to generate genome-wide chromatin occupancy profiles of nucleosomes, transcription factors, and replication proteins through consecutive cell cycles in Saccharomyces cerevisiae. During each G1 phase of two consecutive cell cycles, we observed the downstream repositioning of the origin-proximal +1 nucleosome and an increase in protected DNA fragments spanning the ARS consensus sequence (ACS) indicative of pre-RC assembly. We also found that the strongest correlation between chromatin occupancy at the ACS and origin efficiency occurred in early S phase, consistent with the rate-limiting formation of the Cdc45–Mcm2-7–GINS (CMG) complex being a determinant of origin activity. Finally, we observed nucleosome disruption and disorganization emanating from replication origins and traveling with the elongating replication forks across the genome in S phase, likely reflecting the disassembly and assembly of chromatin ahead of and behind the replication fork, respectively. These results provide insights into cell-cycle–regulated chromatin dynamics and how they relate to the regulation of origin activity.


Introduction
Duplication of a cell's genetic information occurs every cell cycle in S phase. While DNA replication is restricted to S phase, the DNA replication program is established earlier in the cell cycle with the licensing of DNA replication origins in G1 [1]. These licensed origins are then activated with an inherent efficiency during S phase [2]. While many studies have carefully examined the kinetics of DNA replication progression through S phase [3][4][5][6][7], few have examined the chromatin dynamics of replication origins as cells progress through consecutive cell cycles. Instead, most chromatin-based studies have interrogated static snapshots of chromatin from discrete cell-cycle phases (e.g., G1) or from an asynchronous population of cells [8][9][10][11]. Understanding how chromatin structure and organization change at DNA replication origins as they progress through consecutive unperturbed cell cycles will provide important insights into the chromatin features that modulate origin usage and efficiency.
The selection and activation of replication origins involve the recruitment of a series of replication factors in an ordered manner through multiple phases of the cell cycle. In S. cerevisiae, potential origins (which are defined by an autonomous replicating sequence (ARS) and contain a conserved T-rich ARS consensus sequence (ACS)) are recognized and bound by the heterohexameric origin recognition complex (ORC) [12,13]. In G1, the Mcm2-7 replicative helicase is loaded in a Cdc6-and Cdt1-dependent manner to form the pre-replicative complex (pre-RC) and license origins for activation [14]. Additional initiation factors, including Cdc45 and GINS, are recruited in S phase to assemble the active helicase known as the Cdc45-Mcm2-7-GINS (CMG) complex that commits origins for activation [15]. The concentrations of Cdc45 and components of GINS are significantly lower than those of the pre-RC and serve as a rate-limiting step in the activation of individual origins [16,17].
Only a small subset of ACS motif matches in the yeast genome are bona fide ORC binding sites, and chromatin architecture is believed to be an important factor in defining origins [8,18,19]. Origins of replication in the budding yeast have a stereotypical chromatin architecture with well-positioned nucleosomes surrounding a nucleosome-free region (NFR) containing the ACS [8,20]. The chromatin architecture with well-positioned nucleosomes appears to be a conserved feature of eukaryotic origins [21][22][23]. In addition, the inherent initiation efficiency [2] and the time of activation [3,24] of each origin are also thought to be regulated, in part, by the local chromatin environment and the levels of chromatin-associated ORC [20,25] and Mcm2-7 [26,27].
The selection and activation of DNA replication origins are tightly coupled with the cell cycle. Concomitant with the recruitment of factors to the origin in G1 to assemble and load the helicase, the NFR expands in a Cdc6-dependent manner to accommodate pre-RC assembly [20,28]. Further, as cells enter S phase in the absence of primase activity, the disorder or entropy of the origin-flanking chromatin increases markedly, presumably due to helicase activation and nucleosome eviction [29]. Finally, the stability and occupancy of ORC on the DNA throughout the cell cycle are origin-dependent and predictive of origin efficiency [25]. The majority of origins exhibit a protected footprint representing ORC and/or the pre-RC in the NFR in both G1 and G2 phases of the cell cycle; however, select origins only exhibit a defined footprint in G1, suggesting a more dynamic or transient interaction with ORC which may be stabilized by pre-RC assembly [20]. The stability of ORC and/or pre-RC components on the DNA is predictive of origin efficiency, with more efficient origins having a protected footprint in both G1 and G2.
To better understand cell-cycle-regulated chromatin dynamics with high spatiotemporal resolution, we generated genome-wide chromatin occupancy profiles [20,30] of chromatin sampled at multiple points throughout two consecutive cell cycles. This approach provides factor-agnostic occupancy profiles of DNA-binding proteins, including nucleosomes, replication factors, and transcription factors at nucleotide resolution. We comprehensively profiled the dynamics of protected fragments at the ACS and the organization of ACS-flanking nucleosomes throughout two cell cycles. Our study describes the chromatin architecture at individual origins in a synchronized cell population and associates cell-cycle-dependent chromatin features with origin efficiency, thus providing mechanistic insight into the dynamic interplay between chromatin architecture and origin function.

Yeast Culture and Cell-Cycle Time Courses
In this study, we used the W303 yeast strain with the genotype MATa, leu2-3,112, BAR1::TRP, can1-100, URA3::BrdU-Inc, ade2-1, his3-11,15. Yeast cells were grown in YPD (1% yeast extract, 2% peptone, 2% dextrose) at 30 • C to an OD 600 of~0.3 and arrested in G1 with α-factor (GenWay) at a final concentration of 50 ng/mL for 2 h. Samples were taken right before release and labeled as the "α-factor" time point. Cells were then washed twice in sterile water, resuspended in fresh YPD medium, and samples were collected every 10 min until 150 min post-release. For each time point, 40 mL of culture was crosslinked with a final concentration of 1% formaldehyde at room temperature for 30 min, quenched with 0.125 M glycine for 5 min, washed, and flash frozen. In parallel, 1 mL of culture was resuspended in 70% ethanol and fixed overnight at 4 • C for flow cytometry. Independent biological duplicates were performed.

Chromatin Digestion with MNase and Sequencing Library Preparation
MNase digestion of chromatin and sequencing library preparation were performed as previously described [20,30] with the following modifications: 2 µg of digested DNA was used as input; NEBNext multiplex oligos for Illumina Kit (New England Biolabs, Ipswich, MA, USA) were used in adapter ligation and PCR amplification steps; PCR reactions were performed with 12 cycles; and libraries were purified using Agencourt AMPure XP beads (Beckman, Brea, CA, USA). Libraries were sequenced on NextSeq 500 High-Output 25 bp PE platform (Illumina, San Diego, CA, USA).

Flow Cytometry
Fixed yeast cells were washed with water, briefly sonicated, and incubated in 50 mM sodium citrate (pH 7.4) with 0.3 mg/mL RNase A for 2 h at 50 • C. Then, 0.6 mg/mL Proteinase K (Worthington, Columbus, OH, USA) was added and incubated for an additional 2 h at 50 • C. Finally, cell pellets were resuspended in 50 mM sodium citrate with 1:5000 SYTOX green (Invitrogen, Waltham, MA, USA) and incubated for 1 h at room temperature. Flow cytometry was performed on a BD FACSCanto analyzer, and 30,000 cells were recorded for each sample.

Sequencing Data Processing and Analysis
All reads were aligned to the sacCer3/R64 version of the S. cerevisiae genome using Bowtie 0.12.7 [31]. MNase-seq reads were mapped in paired-end mode with the following Bowtie parameters: -n 2 -l 20 -phred33-quals -m 1 -best -strata -y. Data analysis was performed in R version 3.2.0. All genomic data are publicly available at the NCBI GEO repository with the accession number GSE168699.
Because the position of each MNase-seq fragment could be determined by the coordinate of one end and the fragment length, only reads mapped on the forward strand were kept. MNase-seq data from the biological duplicates were randomly subsampled and merged to reduce bias from MNase digestion, library preparation, and sequencing depth before downstream analysis. For each replicate over all time points, the fewest number of reads for each fragment size (from 20 to 250 bp) was identified and used as the subsampling depth. The MNase-seq data for each time point were then subsampled to the above depth per fragment size so as to assign an equal number of reads for each fragment length among all time points. Reads mapped to mitochondrial DNA (chrM) or the ribosomal DNA locus (chrXII:451,575-489,469) were excluded. After subsampling, the total number of reads for each time point was~17 million and~21 million for Replicates 1 and 2, respectively. The matched time points between duplicates were merged for downstream analysis.

Quantification of Nucleosome Occupancy
For each time point, a pileup matrix of fragment size by fragment midpoint position was calculated for the aggregate MNase-seq signal of 8632 unique nucleosome positions on chrIV which were mapped by a sensitive chemical mapping method [32]. This matrix represents the approximate size and coverage distribution of MNase-seq reads centered at a canonical well-positioned nucleosome. A two-dimensional kernel was then derived using a bivariate Gaussian distribution parameterized by the marginal means and variances of the matrix [33]. The variance of fragment size dimension (y-axis) was set to 1/16 of the original marginal variance and the variance of midpoint position dimension (x-axis) was set to 1/4 of the original marginal variance. To quantify the occupancy signal of a nucleosome at a given chromosomal location, a cross-correlation score was computed between the local MNase-seq signal matrix and the model nucleosome kernel.
To correct for replication-dependent DNA copy number variation throughout the cell cycle for a given chromosomal location, the RPKM of all MNase-seq reads for a 1001 bp window centered at the given position was calculated for each time point, and the ratio of this RPKM to the RPKM of the α-factor time point (G1) was considered to be the copy number. For each time point, the nucleosome score for any chromosomal location was normalized by its copy number.

Quantification of Small Fragment Occupancy
For each chromosome, the midpoint density of fragments smaller than 120 bp was estimated using a Gaussian kernel with a bandwidth of 50 bp. The small fragment occupancy of a given chromosomal location was calculated as the product of density and chromosome length and normalized by the copy number of the given position. To adjust for variations in MNase digestion among samples, the average signal of the aggregate small fragment occupancy within ±100 bp around 151 Abf1p binding sites identified by Mac-Isaac et al. [34] was calculated for each time point and the reciprocal of this average was used as a scale factor.

Quantification of Nucleosome Disorganization by Shannon Entropy
For a region X of size n bp, the probability of nucleosome positioning at location i was defined as where nuc(i) is the nucleosome score at position i. The disorganization of nucleosome positioning for region X was measured using Shannon entropy:

Chromatin Occupancy Profiling of Replication Origins Throughout the Cell Cycle
We sought to profile the cell-cycle-dependent changes in chromatin organization surrounding replication origins throughout the yeast genome. Cells were synchronized in late G1 using α-factor. Cells were then released from α-factor arrest and samples were collected every 10 min for approximately two complete cell cycles (150 min) ( Figure 1A). Biological replicates were performed, and progression through the cell cycle was monitored by flow cytometry (Supplementary Figure S1). Samples were aligned by cell-cycle progression and data from corresponding samples were merged for downstream analysis.
To comprehensively interrogate chromatin dynamics at nucleotide resolution, we generated genome-wide chromatin occupancy profiles by digesting chromatin with MNase followed by paired-end sequencing [20,30,35]. DNA fragments protected by DNA-binding factors (nucleosomes, transcription factors, ORC, etc.) were recovered and subjected to next-generation sequencing. The length and location of the mapped reads provide an unbiased view of chromatin occupancy throughout the genome. For example, nucleosomes protect DNA fragments of~150 bp, and smaller DNA-binding factors (e.g., transcription factors and ORC) protect DNA fragments less than 120 bp. The mode of the distribution of fragment sizes recovered was 166 bp (Supplementary Figure S2), consistent with the majority of DNA being packaged into nucleosomes throughout the genome. To visualize chromatin architecture at individual replication origins, we plotted the length of each fragment as a function of the chromosomal position of its midpoint. At ARS1623, an efficient and early origin, we observed an array of well-positioned nucleosomes flanking the origin, which appeared as clusters of fragment midpoints centered at~160 bp. We also observed an accumulation of smaller fragments at the ACS which represents an ORCdependent footprint [20] ( Figure 1B). Similar chromatin organization was observed at ARS228.5, an inefficient and late origin, albeit with "fuzzier" ACS-proximal nucleosomes, a narrower nucleosome-free region (NFR), and a substantially weaker small fragment footprint at the ACS ( Figure 1C). We observed dynamic changes in chromatin organization as cells proceeded through the cell cycle. As shown for select time points representing late G1 (α-factor), early S (20 min), and M (60 min) phases, we observed fluctuations in the dyad positions of the ACS-proximal nucleosomes and in the occupancy of ACS-bound small fragments ( Figure 1B,C). A score for nucleosome occupancy and position was calculated using a two-dimensional nucleosome kernel modeled on the MNase fragments associated with nucleosomes mapped by an orthogonal chemical cleavage method [32,33] (Supplementary Figure S3). A score for small factor occupancy (e.g., ORC, pre-RC, and pre-IC) at origins was generated by calculating the density of small protected fragments less than 120 bp at each origin of DNA replication. To adjust for variation in sample-specific MNase digestion, we normalized the occupancy of the small fragment footprint at the ACS by the occupancy of the footprint at Abf1p binding sites (Supplementary Figure S4). At the efficient ARS1623, the downstream (+1) nucleosome of the ACS was displaced further from the ACS at the α-factor time point compared with other time points in the first cell cycle ( Figure 1D). The movement of the nucleosome away from the ACS in G1 is Cdc6-dependent and presumably facilitates pre-RC assembly [20]. In contrast, the +1 nucleosome of the inefficient ARS228.5 appeared to be more static ( Figure 1E). While the occupancy of ACS-bound small fragments is significantly stronger at ARS1623, both origins showed fluctuating small fragment occupancy throughout the cell cycle ( Figure 1D,E), likely reflecting the cell-cycle-coupled dynamics of helicase loading in G1 and the subsequent activation of the CMG holohelicase complex and its movement away from the origin as the cell progresses into and through S phase.

Cell-Cycle-Dependent Changes in Replication Initiation Factor Occupancy at Replication Origins
The loading of the helicase to form the pre-RC and the recruitment of replication initiation factors at the origin are tightly coupled to the progression of the cell cycle. Each replication origin in the genome has an inherent efficiency which is thought to be determined, at least in part, by ORC binding, Mcm2-7 loading, the recruitment of activation factors, and the local chromatin environment [2,20,26,27,36,37]. We had previously used genome-wide chromatin occupancy profiling to identify ORC-dependent small fragment occupancy footprints at replication origins in two discrete phases of the cell cycle, G1 and G2 [20]. In that study, we identified two classes of origins: the first consisted of 264 origins that exhibited an ORC-dependent footprint in both G1 and G2 while a second, less efficient class of 128 origins exhibited an ORC-dependent footprint only in G1. To precisely quantify the dynamic changes in small fragment footprints at the ACS throughout the cell cycle, we calculated the small fragment density at the ACS of each origin with a bandwidth of 50 bp at each time point. We observed a periodic protection footprint at the ACS for both classes of origins (Figure 2A,B) which peaked in late G1/early S phase in consecutive cell cycles and likely represents the recruitment of Cdc45 and GINS to the pre-RC to form the CMG holohelicase complex at the most efficient origins. After peaking in late G1/early S phase, the density of small fragments at the ACS gradually declines to a nadir near mitosis. This decline in signal throughout S phase likely reflects the activation and the bi-directional movement of the CMG holohelicase complex with the replication fork away from the origin and/or the disassembly of the pre-RC at any passively replicated origins. By late G2 and through mitosis, it is likely that any remaining small fragments at the origins are due to ORC.
The signal in the origins exhibiting an ORC-dependent footprint in "G1 only" was significantly dampened relative to the "G1 & G2" class of origins ( Figure 2C). Although we still observe an oscillatory pattern, the small fragment occupancy signals in G2 and M are barely at the detectable limit. Given the decreased signal observed throughout the cell cycle, this suggests a defect in either the recruitment or the stability of ORC on the DNA at these origins which ultimately leads to a stochastic defect in downstream helicase loading.
We previously found that stable ORC binding in both G1 and G2 was a determinant of efficient origins [20]. However, with the more comprehensive view of small fragment occupancy at the ACS throughout the cell cycle that our data provide, we reasoned that a better predictor of origin efficiency might be the level of protected fragments at each origin as cells progress from G1 into S phase. To test this hypothesis, we calculated the correlation between ACS-bound small fragment occupancy and origin efficiency at each time point throughout the cell cycle ( Figure 2D; Supplementary Figure S5). We found a cyclic pattern, with the correlation between origin efficiency and small fragment occupancy peaking in early S phase and reaching its lowest point near mitosis. These results are consistent with recruitment of origin activation factors like Cdc45 and GINS to form the CMG holohelicase complex in late G1 and early S being a key determinant of origin activation. Heatmaps of aggregate small fragment (<120 bp) occupancy at 264 origins with a previously described ORC-dependent footprint in both "G1 & G2" (A) and 128 origins with an ORC-dependent footprint in "G1 only" (B) plotted throughout the cell cycle [20]. All origins are oriented by the T-rich ACS strand. (C) Average small fragment occupancy ±100 bp surrounding the peak of the aggregate ORC-dependent footprint for each class of origins. (D) Spearman correlation between log2-transformed ACS-bound small fragment footprint density and activation efficiency [2] for 371 origins exhibiting an ORC-dependent footprint at each time point.

Cell-Cycle-Regulated Nucleosome Occupancy Dynamics Around Replication Origins
The position of nucleosomes relative to the ACS impacts origin function [8,9,[38][39][40]. We analyzed the aggregate nucleosome signal of each origin class to profile the cell-cycledependent dynamics of origin-flanking nucleosomes. For both classes, the conserved chromatin structure of replication origins-an NFR at the ACS surrounded by two wellpositioned nucleosomes-was maintained throughout the cell cycle ( Figure 3A). We first analyzed the positioning of the NFR-proximal -1 and +1 nucleosomes throughout the cell cycle. We defined the dyad of either nucleosome as the position with maximal nucleosome score upstream or downstream of the ACS. The NFRs of origins with a footprint in both G1 and G2 are wider throughout the cell cycle (average dyad-to-dyad distance: 257.8 bp versus 220.9 bp). The increased width of the NFR was due to the downstream localization of the +1 nucleosome while the position of the -1 nucleosomes remained static ( Figure 3B). Notably, the +1 nucleosome of origins with a G1 and G2 footprint exhibited periodic repositioning throughout the cell cycle, coincident with the helicase loading and recruitment of factors to form the Cdc45-Mcm2-7-GINS (CMG) complex. For example, the farthest +1 nucleosome positioning from the ACS was observed in the first G1/S transition (10 min), indicating recruitment of the CMG complex and pre-IC assembly. After origin activation, we observed that +1 nucleosomes gradually started to move back towards the ACS with cell-cycle progression, likely reflecting the subsequent dissociation of replication factors from origins. The extent of the nucleosome shift for the second cell cycle became weaker due to the gradual loss of cell synchrony. The +1 nucleosome positioning of origins with a "G1 only" footprint did not fluctuate in a cell-cycle-dependent manner which may be attributed to reduced formation of the rate-limiting CMG complex at less efficient origins [17,40].

Replication-Coupled Nucleosome Disorganization
The progression of the replication fork results in eviction of histones and disruption of nucleosome positioning [19]. We utilized Shannon entropy to assess the level of nucleosome organization [33] throughout the cell cycle. Well-positioned nucleosomes have a low entropy while disorganized nucleosomes exhibit a high entropy. We recently reported that helicase activation in the absence of primase activity resulted in the disorganization of origin-proximal nucleosomes [29]. We reasoned that nucleosome disruption by helicaseinduced unwinding would be a feature of sequences at the active replication fork. We calculated the entropy of nucleosome signals in 1 kb windows for 30 kb surrounding the 69 most efficient or active origins ( Figure 4A) and the 69 least efficient origins which are passively replicated ( Figure 4B). The mean entropy scores for each window were then ordered by their distance from the origin (rows) as a function of progression through the cell cycle (columns). For both the active and passive origins, we observed increased entropy in S phase, consistent with disruption of chromatin by the passage of the replication fork. For the sequences proximal to active origins, we observed a distance-dependent temporal shift in peak entropy from early S (20 min) to mid S phase (30 min), consistent with earlier timing of activation for efficient origins.  [2]. Each row represents a genomic window and rows are ordered by the distance of that window from the nearest origin. Nucleosome entropy is standardized into z-scores across each row.

Discussion
We interrogated chromatin dynamics at replication origins over the course of two consecutive cell cycles. As cells progressed through each cell cycle, we observed specific changes in chromatin occupancy that correspond to key steps in the selection and activation of start sites of DNA replication. We found that for each cell cycle, the downstream +1 nucleosome relative to the T-rich ACS moved further downstream in late G1 and early S phase before returning back to its original position by G2. Similarly, we observed a cyclic increase in small protected fragments spanning the ACS that likely represent pre-RC assembly and the recruitment of origin activation factors in late G1/early S. The accumulation of protected fragments at the ACS in early S phase was more predictive of active and efficient origins than at other times in the cell cycle. Finally, we observed genome-wide disruption of nucleosomes traveling with active replication forks during S phase, reflecting helicase-induced DNA unwinding at the fork. Together, these data provide insights into the dynamic interplay between origin function and local chromatin environment with high spatiotemporal resolution.
In contrast to chromatin immunoprecipitation, genome-wide chromatin occupancy profiling by MNase digestion offers a holistic view of DNA-protein interactions that neither requires factor-specific antibodies nor is encumbered by epitope accessibility [5]. A consequence of genome-wide chromatin occupancy being factor-agnostic is that it only reports if a DNA sequence is protected or occupied. The factor responsible for the occupancy has to be inferred from other sources such as sequence context (e.g., an ACS motif) or existing factor-specific localization experiments. Thus, we need to infer that the subtle cell-cycle-dependent changes in chromatin occupancy at the ACS reflect distinct steps in licensing and activation of replication origins. Importantly, in support of these inferences, our prior work using static G1-and G2-arrested samples demonstrated that chromatin occupancy at the ACS was ORC-dependent and that the increase in occupancy in G1 was dependent on origin licensing or pre-RC assembly [20].
The determinants of origin efficiency are poorly understood at the chromatin level. While it is clear that rate-limiting activation factors establish the temporal order of origin activation [16,17], what is less clear is why and how individual origins are more or less sensitive to these factors. Local histone modifications have been shown to modulate replication timing, with early work demonstrating that origin activity could be repressed by moving an active early origin to a heterochromatic region of the genome [41]. Global perturbation of histone acetylation levels also impacts genome-wide replication timing [10,42], but it is unclear if this effect is direct or indirect through modulation of the accessibility of the multicopy rDNA locus and consequent sequestration of origin activation factors [11]. The association of ORC with DNA can be driven by sequence or chromatin elements, with the chromatin-dependent ORC class being correlated with early replicating efficient origins [25]. The relative amount of Mcm2-7 loading in G1 is also predictive of origin function [27]. Our work identifies an increase in protected fragment occupancy at the ACS in early S phase that is predictive of efficient origin activation, consistent with the formation of the CMG holohelicase complex at licensed origins. While this does not necessarily preclude earlier steps like ORC binding or Mcm2-7 loading from being deterministic, it is notable that the accumulation of protected fragments at the ACS during G2 (with ORC alone) or during G1 were less predictive of origin efficiency.
It has been proposed that efficient loading of the Mcm2-7 double-hexamer is established by a quasi-symmetrical loading mechanism that requires the binding of two ORC molecules to two ORC binding sites at yeast origins [43]. We examined the chromatin architecture at two example origins explored by Coster et al. [43], namely ARS606 and ARS1216. During the first cell cycle, we noted two regions protected by small fragments near the ACSs, likely reflecting multiple sites of ORC binding (Supplementary Figure S6). However, these origins appeared to be the exception rather than the rule, as we did not detect multiple protected regions in the vast majority of replication origins. A more recent model posits that a single ORC molecule may relocate from the ACS to the downstream B2 element in order to load the Mcm2-7 double-hexamer [44]. Our data do not discount this model, but we lack the temporal resolution to detect the relocation of ORC to the B2 element during the transient process of helicase loading.
A consequence of the passage of the replication fork is disruption of chromatin ahead of the fork and the restoration of chromatin in the wake of the fork. Our chromatin occupancy profiling was able to capture the dynamic and transient disruption of chromatin associated with the replication fork. Specifically, we observed a transient and S-phasespecific increase in entropy that was temporally linked to the distance from the nearest replication origin. The reassembly of nucleosomes and the re-establishment of the chromatin landscape behind the fork is critical for epigenetic inheritance, and it is known that factors that impair or delay assembly can significantly impact gene regulation [45][46][47] and differentiation [48]. Our ability to discern differences in entropy associated with replication-coupled nucleosome assembly will enable future studies, with increased temporal resolution, to identify and characterize locus-specific differences in chromatin assembly that may be governed in part by chromosome position, the local chromatin environment, and transcription.

Conclusions
Chromatin surrounding DNA replication origins is dynamically organized throughout the cell cycle. During G1 and S phase, the downstream +1 nucleosome is re-positioned further away from the origin to accommodate helicase loading and activation. Similarly, we detect cell-cycle-dependent increases in small protected DNA fragments at the origin that likely represent helicase loading in G1 and the formation of the CMG holohelicase complex in S phase. The protection of small DNA fragments at origins in early S phase was most predictive of origin efficiency, consistent with the formation of the CMG holohelicase complex being a rate-limiting step for origin activation. Finally, we observed the S-phasespecific disruption of chromatin at active DNA replication forks.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/genes12121998/s1, Figure S1: Flow cytometry profiles monitoring cell-cycle progression for each biological duplicate. Figure S2: Distribution of MNase fragment lengths after subsampling and merging replicates. Figure S3: Density plot of the two-dimensional kernel modeled from an orthogonal set of nucleosome dyads. Figure S4: Small fragment footprint occupancy at Abf1p binding sites throughout the cell cycle. Figure S5: Correlation between small factor occupancy at the ACS and origin efficiency. Figure S6: Chromatin profiles at ARS606 and ARS1216 for select time points during the first cell cycle.