Nucleosome occupancy as a novel chromatin parameter for replication origin functions

Sir2 and Fun30 regulate ribosomal DNA replication timing via Mcm helicase positioning and nucleosome occupancy

The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well-known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). We have previously reported that in the absence of SIR2, a derepressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the non-displaced MCMs. Furthermore, we found that both, activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Investigation of the acetic acid stress response in Saccharomyces cerevisiae with mutated H3 residues

Enhanced levels of acetic acid reduce the activity of yeast strains employed for industrial fermentation-based applications. Therefore, unraveling the genetic factors underlying the regulation of the tolerance and sensitivity of yeast towards acetic acid is imperative for optimising various industrial processes. In this communication, we have attempted to decipher the acetic acid stress response of the previously reported acetic acid-sensitive histone mutants. Revalidation using spot-test assays and growth curves revealed that five of these mutants, viz., H3K18Q, H3S28A, H3K42Q, H3Q68A, and H3F104A, are most sensitive towards the tested acetic acid concentrations. These mutants demonstrated enhanced acetic acid stress response as evidenced by the increased expression levels of AIF1, reactive oxygen species (ROS) generation, chromatin fragmentation, and aggregated actin cytoskeleton. Additionally, the mutants exhibited active cell wall damage response upon acetic acid treatment, as demonstrated by increased Slt2-phosphorylation and expression of cell wall integrity genes. Interestingly, the mutants demonstrated increased sensitivity to cell wall stress-causing agents. Finally, screening of histone H3 N-terminal tail truncation mutants revealed that the tail truncations exhibit general sensitivity to acetic acid stress. Some of these N-terminal tail truncation mutants viz., H3 [del 1-24], H3 [del 1-28], H3 [del 9-24], and H3 [del 25-36] are also sensitive to cell wall stress agents such as Congo red and caffeine suggesting that their enhanced acetic acid sensitivity may be due to cell wall stress induced by acetic acid.

Genome-wide binding sites of Plasmodium falciparum mini chromosome maintenance protein MCM6 show new insights into parasite DNA replication

Multiple rounds of DNA replication take place in various stages of the life cycle in the human malaria parasite Plasmodium falciparum. Previous bioinformatics analysis has shown the presence of putative Autonomously Replicating Sequence (ARS) like sequences in the Plasmodium genome. However, the actual sites and frequency of replication origins in the P. falciparum genome based on experimental data still remain elusive. Minichromosome maintenance (MCM) proteins are recruited by the Origin recognition complex (ORC) to the origins of replication in eukaryotes including P. falciparum. We used PfMCM6 for chromatin immunoprecipitation followed by sequencing (ChIP-seq) in the quest for identification of putative replication origins in the parasite. PfMCM6 DNA binding sites annotation revealed high enrichment at exon regions. This is contrary to higher eukaryotes that show an inclination of origin sites towards transcriptional start sites. ChIP-seq results were further validated by ChIP-qPCR results as well as nascent strand abundance assay at the selected PfMCM6 enriched sites that also showed preferential binding of PfORC1 suggesting potential of these sites as origin sites. Further, PfMCM6 ChIP-seq data showed a positive correlation with previously published histone H4K8Ac genome-wide binding sites but not with H3K9Ac sites suggesting epigenetic control of replication initiation sites in the parasites. Overall, our data show the genome-wide distribution of PfMCM6 binding sites with their potential as replication origins in this deadly human pathogen that not only broadens our knowledge of parasite DNA replication and its unique biology, it may help to find new avenues for intervention processes.

Impact of Chromosomal Context on Origin Selection and the Replication Program

Eukaryotic DNA replication is regulated by conserved mechanisms that bring about a spatial and temporal organization in which distinct genomic domains are copied at characteristic times during S phase. Although this replication program has been closely linked with genome architecture, we still do not understand key aspects of how chromosomal context modulates the activity of replication origins. To address this question, we have exploited models that combine engineered genomic rearrangements with the unique replication programs of post-quiescence and pre-meiotic S phases. Our results demonstrate that large-scale inversions surprisingly do not affect cell proliferation and meiotic progression, despite inducing a restructuring of replication domains on each rearranged chromosome. Remarkably, these alterations in the organization of DNA replication are entirely due to changes in the positions of existing origins along the chromosome, as their efficiencies remain virtually unaffected genome wide. However, we identified striking alterations in origin firing proximal to the fusion points of each inversion, suggesting that the immediate chromosomal neighborhood of an origin is a crucial determinant of its activity. Interestingly, the impact of genome reorganization on replication initiation is highly comparable in the post-quiescent and pre-meiotic S phases, despite the differences in DNA metabolism in these two physiological states. Our findings therefore shed new light on how origin selection and the replication program are governed by chromosomal architecture.

Histone-Net: a multi-paradigm computational framework for histone occupancy and modification prediction

Deep exploration of histone occupancy and covalent post-translational modifications (e.g., acetylation, methylation) is essential to decode gene expression regulation, chromosome packaging, DNA damage, and transcriptional activation. Existing computational approaches are unable to precisely predict histone occupancy and modifications mainly due to the use of sub-optimal statistical representation of histone sequences. For the establishment of an improved histone occupancy and modification landscape for multiple histone markers, the paper in hand presents an end-to-end computational multi-paradigm framework “Histone-Net”. To learn local and global residue context aware sequence representation, Histone-Net generates unsupervised higher order residue embeddings (DNA2Vec) and presents a different application of language modelling, where it encapsulates histone occupancy and modification information while generating higher order residue embeddings (SuperDNA2Vec) in a supervised manner. We perform an intrinsic and extrinsic evaluation of both presented distributed representation learning schemes. A comprehensive empirical evaluation of Histone-Net over ten benchmark histone markers data sets for three different histone sequence analysis tasks indicates that SuperDNA2Vec sequence representation and softmax classifier-based approach outperforms state-of-the-art approach by an average accuracy of 7%. To eliminate the overhead of training separate binary classifiers for all ten histone markers, Histone-Net is evaluated in multi-label classification paradigm, where it produces decent performance for simultaneous prediction of histone occupancy, acetylation, and methylation.

Differential condensation of sister chromatids acts with Cdc6 to ensure asynchronous S-phase entry in Drosophila male germline stem cell lineage

During Drosophila melanogaster male germline stem cell (GSC) asymmetric division, preexisting old versus newly synthesized histones H3 and H4 are asymmetrically inherited. However, the biological outcomes of this phenomenon have remained unclear. Here, we tracked old and new histones throughout the GSC cell cycle through the use of high spatial and temporal resolution microscopy. We found unique features that differ between old and new histone-enriched sister chromatids, including differences in nucleosome density, chromosomal condensation, and H3 Ser10 phosphorylation. These distinct chromosomal features lead to their differential association with Cdc6, a pre-replication complex component, and subsequent asynchronous DNA replication initiation in the resulting daughter cells. Disruption of asymmetric histone inheritance abolishes differential Cdc6 association and asynchronous S-phase entry, demonstrating that histone asymmetry acts upstream of these critical cell-cycle progression events. Furthermore, disruption of these GSC-specific chromatin features leads to GSC defects, indicating a connection between histone inheritance, cell-cycle progression, and cell fate determination.

Nucleosome positioning based on DNA sequence embedding and deep learning

Background

Nucleosome positioning is the precise determination of the location of nucleosomes on DNA sequence. With the continuous advancement of biotechnology and computer technology, biological data is showing explosive growth. It is of practical significance to develop an efficient nucleosome positioning algorithm. Indeed, convolutional neural networks (CNN) can capture local features in DNA sequences, but ignore the order of bases. While the bidirectional recurrent neural network can make up for CNN's shortcomings in this regard and extract the long-term dependent features of DNA sequence.

Results

In this work, we use word vectors to represent DNA sequences and propose three new deep learning models for nucleosome positioning, and the integrative model NP_CBiR reaches a better prediction performance. The overall accuracies of NP_CBiR on H. sapiens, C. elegans, and D. melanogaster datasets are 86.18%, 89.39%, and 85.55% respectively.

Conclusions

Benefited by different network structures, NP_CBiR can effectively extract local features and bases order features of DNA sequences, thus can be considered as a complementary tool for nucleosome positioning.

Preventing excess replication origin activation to ensure genome stability

Cells activate distinctive regulatory pathways that prevent excessive initiation of DNA replication to achieve timely and accurate genome duplication. Excess DNA synthesis is constrained by protein-DNA interactions that inhibit initiation at dormant origins. In parallel, specific modifications of pre-replication complexes prohibit post-replicative origin relicensing. Replication stress ensues when the controls that prevent excess replication are missing in cancer cells, which often harbor extrachromosomal DNA that can be further amplified by recombination-mediated processes to generate chromosomal translocations. The genomic instability that accompanies excess replication origin activation can provide a promising target for therapeutic intervention. Here we review molecular pathways that modulate replication origin dormancy, prevent excess origin activation, and detect, encapsulate, and eliminate persistent excess DNA.

Chemical map-based prediction of nucleosome positioning using the Bioconductor package nuCpos

Background

Assessing the nucleosome-forming potential of specific DNA sequences is important for understanding complex chromatin organization. Methods for predicting nucleosome positioning include bioinformatics and biophysical approaches. An advantage of bioinformatics methods, which are based on in vivo nucleosome maps, is the use of natural sequences that may contain previously unknown elements involved in nucleosome positioning in vivo. The accuracy of such prediction attempts reflects the genomic coordinate resolution of the nucleosome maps applied. Nucleosome maps are constructed using micrococcal nuclease digestion followed by high-throughput sequencing (MNase-seq). However, as MNase has a strong preference for A/T-rich sequences, MNase-seq may not be appropriate for this purpose. In addition to MNase-seq-based maps, base pair-resolution chemical maps of in vivo nucleosomes from three different species (budding and fission yeasts, and mice) are currently available. However, these chemical maps have yet to be integrated into publicly available computational methods.

Results

We developed a Bioconductor package (named nuCpos) to demonstrate the superiority of chemical maps in predicting nucleosome positioning. The accuracy of chemical map-based prediction in rotational settings was higher than that of the previously developed MNase-seq-based approach. With our method, predicted nucleosome occupancy reasonably matched in vivo observations and was not affected by A/T nucleotide frequency. Effects of genetic alterations on nucleosome positioning that had been observed in living yeast cells could also be predicted. nuCpos calculates individual histone binding affinity (HBA) scores for given 147-bp sequences to examine their suitability for nucleosome formation. We also established local HBA as a new parameter to predict nucleosome formation, which was calculated for 13 overlapping nucleosomal DNA subsequences. HBA and local HBA scores for various sequences agreed well with previous in vitro and in vivo studies. Furthermore, our results suggest that nucleosomal subsegments that are disfavored in different rotational settings contribute to the defined positioning of nucleosomes.

Conclusions

Our results demonstrate that chemical map-based statistical models are beneficial for studying nucleosomal DNA features. Studies employing nuCpos software can enhance understanding of chromatin regulation and the interpretation of genetic alterations and facilitate the design of artificial sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04240-2.

Comparative analysis and prediction of nucleosome positioning using integrative feature representation and machine learning algorithms

Background

Nucleosome plays an important role in the process of genome expression, DNA replication, DNA repair and transcription. Therefore, the research of nucleosome positioning has invariably received extensive attention. Considering the diversity of DNA sequence representation methods, we tried to integrate multiple features to analyze its effect in the process of nucleosome positioning analysis. This process can also deepen our understanding of the theoretical analysis of nucleosome positioning.

Results

Here, we not only used frequency chaos game representation (FCGR) to construct DNA sequence features, but also integrated it with other features and adopted the principal component analysis (PCA) algorithm. Simultaneously, support vector machine (SVM), extreme learning machine (ELM), extreme gradient boosting (XGBoost), multilayer perceptron (MLP) and convolutional neural networks (CNN) are used as predictors for nucleosome positioning prediction analysis, respectively. The integrated feature vector prediction quality is significantly superior to a single feature. After using principal component analysis (PCA) to reduce the feature dimension, the prediction quality of H. sapiens dataset has been significantly improved.

Conclusions

Comparative analysis and prediction on H. sapiens, C. elegans, D. melanogaster and S. cerevisiae datasets, demonstrate that the application of FCGR to nucleosome positioning is feasible, and we also found that integrative feature representation would be better.

INO80C Remodeler Maintains Genomic Stability by Preventing Promiscuous Transcription at Replication Origins

The proper coordination of transcription with DNA replication and repair is central for genomic stability. We investigate how the INO80C chromatin remodeling enzyme might coordinate these genomic processes. We find that INO80C co-localizes with the origin recognition complex (ORC) at yeast replication origins and is bound to replication initiation sites in mouse embryonic stem cells (mESCs). In yeast· INO80C recruitment requires origin sequences but does not require ORC· suggesting that recruitment is independent of pre-replication complex assembly. In both yeast and ESCs· INO80C co-localizes at origins with Mot1 and NC2 transcription factors· and genetic studies suggest that they function together to promote genome stability. Interestingly· nascent transcript sequencing demonstrates that INO80C and Mot1 prevent pervasive transcription through origin sequences· and absence of these factors leads to formation of new DNA double-strand breaks. We propose that INO80C and Mot1/NC2 function through distinct pathways to limit origin transcription· maintaining genomic stability.

The INO80C chromatin remodeler is known to regulate transcription and genomic stability. Topal et al. find that INO80C functions with Mot1 repressor to prevent pervasive transcription at replication origins in yeast and mESCs. Increased ncRNAs lead to new DNA double-strand breaks at origins, linking transcriptional regulation to genomic stability.

Chromatin structure restricts origin utilization when quiescent cells re-enter the cell cycle

Quiescent cells reside in G0 phase, which is characterized by the absence of cell growth and proliferation. These cells remain viable and re-enter the cell cycle when prompted by appropriate signals. Using a budding yeast model of cellular quiescence, we investigated the program that initiated DNA replication when these G0 cells resumed growth. Quiescent cells contained very low levels of replication initiation factors, and their entry into S phase was delayed until these factors were re-synthesized. A longer S phase in these cells correlated with the activation of fewer origins of replication compared to G1 cells. The chromatin structure around inactive origins in G0 cells showed increased H3 occupancy and decreased nucleosome positioning compared to the same origins in G1 cells, inhibiting the origin binding of the Mcm4 subunit of the MCM licensing factor. Thus, quiescent yeast cells are under-licensed during their re-entry into S phase.

Molecular mechanisms of eukaryotic origin initiation, replication fork progression, and chromatin maintenance

Eukaryotic DNA replication is a highly dynamic and tightly regulated process. Replication involves several dozens of replication proteins, including the initiators ORC and Cdc6, replicative CMG helicase, DNA polymerase α-primase, leading-strand DNA polymerase ε, and lagging-strand DNA polymerase δ. These proteins work together in a spatially and temporally controlled manner to synthesize new DNA from the parental DNA templates. During DNA replication, epigenetic information imprinted on DNA and histone proteins is also copied to the daughter DNA to maintain the chromatin status. DNA methyltransferase 1 is primarily responsible for copying the parental DNA methylation pattern into the nascent DNA. Epigenetic information encoded in histones is transferred via a more complex and less well-understood process termed replication-couple nucleosome assembly. Here, we summarize the most recent structural and biochemical insights into DNA replication initiation, replication fork elongation, chromatin assembly and maintenance, and related regulatory mechanisms.

Arabidopsis DNA Replication Initiates in Intergenic, AT-Rich Open Chromatin

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2'-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

ZCMM: A Novel Method Using Z-Curve Theory- Based and Position Weight Matrix for Predicting Nucleosome Positioning

Nucleosomes are the basic units of eukaryotes. The accurate positioning of nucleosomes plays a significant role in understanding many biological processes such as transcriptional regulation mechanisms and DNA replication and repair. Here, we describe the development of a novel method, termed ZCMM, based on Z-curve theory and position weight matrix (PWM). The ZCMM was trained and tested using the nucleosomal and linker sequences determined by support vector machine (SVM) in Saccharomyces cerevisiae (S. cerevisiae), and experimental results showed that the sensitivity (Sn), specificity (Sp), accuracy (Acc), and Matthews correlation coefficient (MCC) values for ZCMM were 91.40%, 96.56%, 96.75%, and 0.88, respectively, and the average area under the receiver operating characteristic curve (AUC) value was 0.972. A ZCMM predictor was developed to predict nucleosome positioning in Homo sapiens (H. sapiens), Caenorhabditis elegans (C. elegans), and Drosophila melanogaster (D. melanogaster) genomes, and the accuracy (Acc) values were 77.72%, 85.34%, and 93.62%, respectively. The maximum AUC values of the four species were 0.982, 0.861, 0.912 and 0.911, respectively. Another independent dataset for S. cerevisiae was used to predict nucleosome positioning. Compared with the results of Wu's method, it was found that the Sn, Sp, Acc, and MCC of ZCMM results for S. cerevisiae were all higher, reaching 96.72%, 96.54%, 94.10%, and 0.88. Compared with the Guo's method 'iNuc-PseKNC', the results of ZCMM for D. melanogaster were better. Meanwhile, the ZCMM was compared with some experimental data in vitro and in vivo for S. cerevisiae, and the results showed that the nucleosomes predicted by ZCMM were highly consistent with those confirmed by these experiments. Therefore, it was further confirmed that the ZCMM method has good accuracy and reliability in predicting nucleosome positioning.

Pervasive transcription fine-tunes replication origin activity

RNA polymerase (RNAPII) transcription occurs pervasively, raising the important question of its functional impact on other DNA-associated processes, including replication. In budding yeast, replication originates from Autonomously Replicating Sequences (ARSs), generally located in intergenic regions. The influence of transcription on ARSs function has been studied for decades, but these earlier studies have neglected the role of non-annotated transcription. We studied the relationships between pervasive transcription and replication origin activity using high-resolution transcription maps. We show that ARSs alter the pervasive transcription landscape by pausing and terminating neighboring RNAPII transcription, thus limiting the occurrence of pervasive transcription within origins. We propose that quasi-symmetrical binding of the ORC complex to ARS borders and/or pre-RC formation are responsible for pausing and termination. We show that low, physiological levels of pervasive transcription impact the function of replication origins. Overall, our results have important implications for understanding the impact of genomic location on origin function.

Noncoding transcription influences the replication initiation program through chromatin regulation

In eukaryotic organisms, replication initiation follows a temporal program. Among the parameters that regulate this program in Saccharomyces cerevisiae, chromatin structure has been at the center of attention without considering the contribution of transcription. Here, we revisit the replication initiation program in the light of widespread genomic noncoding transcription. We find that noncoding RNA transcription termination in the vicinity of autonomously replicating sequences (ARSs) shields replication initiation from transcriptional readthrough. Consistently, high natural nascent transcription correlates with low ARS efficiency and late replication timing. High readthrough transcription is also linked to increased nucleosome occupancy and high levels of H3K36me3. Moreover, forcing ARS readthrough transcription promotes these chromatin features. Finally, replication initiation defects induced by increased transcriptional readthrough are partially rescued in the absence of H3K36 methylation. Altogether, these observations indicate that natural noncoding transcription into ARSs influences replication initiation through chromatin regulation.

The Protective Role of Dormant Origins in Response to Replicative Stress

Genome stability requires tight regulation of DNA replication to ensure that the entire genome of the cell is duplicated once and only once per cell cycle. In mammalian cells, origin activation is controlled in space and time by a cell-specific and robust program called replication timing. About 100,000 potential replication origins form on the chromatin in the gap 1 (G1) phase but only 20⁻30% of them are active during the DNA replication of a given cell in the synthesis (S) phase. When the progress of replication forks is slowed by exogenous or endogenous impediments, the cell must activate some of the inactive or "dormant" origins to complete replication on time. Thus, the many origins that may be activated are probably key to protect the genome against replication stress. This review aims to discuss the role of these dormant origins as safeguards of the human genome during replicative stress.

Recent advances in understanding DNA replication: cell type-specific adaptation of the DNA replication program

DNA replication is an essential process occurring prior to cell division. Cell division coupled to proliferation ensures the growth and renewal of a large variety of specialized cell types generated during embryonic development. Changes in the DNA replication program occur during development. Embryonic undifferentiated cells show a high replication rate and fast proliferation, whereas more differentiated cells are characterized by reduced DNA synthesis and a low proliferation rate. Hence, the DNA replication program must adapt to the specific features of cells committed to different fates. Recent findings on DNA synthesis regulation in different cell types open new perspectives for developing efficient and more adapted therapies to treat various diseases such as genetic diseases and cancer. This review will put the emphasis on recent progress made in this field.

Application of MNase-Seq in the Global Mapping of Nucleosome Positioning in Plants

The precise positioning of nucleosomes along the underlying DNA is critical for a variety of biological processes, especially in regulating transcription. The interplay between nucleosomes and transcription factors for accessing the underlying DNA sequences is one of the key determinants that affect transcriptional regulation. Moreover, nucleosomes with various packing statuses confer distinct functions in regulating gene expressions in response to various internal or external signals. Therefore, global mapping of nucleosome positions is one informative way to elucidate the relationship between patterns of nucleosome positioning/occupancy and transcriptional regulations. MNase digestion coupled with high-throughput sequencing (MNase-seq) has been utilized widely for global mapping of nucleosome positioning in eukaryotes that have a sequenced genome. We have developed a robust MNase-seq procedure in plants. It mainly includes plant nuclei isolation, treatment of purified nuclei with MNase, gel recovery of MNase-trimmed mononucleosomal DNA with an approximate size of 150 bp, MNase-seq library preparation followed by Illumina sequencing, and data analysis. MNase-seq has already been successfully applied to identify genome-wide nucleosome positioning in model plants, rice, and Arabidopsis thaliana.

The evolution of the temporal program of genome replication

Genome replication is highly regulated in time and space, but the rules governing the remodeling of these programs during evolution remain largely unknown. We generated genome-wide replication timing profiles for ten Lachancea yeasts, covering a continuous evolutionary range from closely related to more divergent species. We show that replication programs primarily evolve through a highly dynamic evolutionary renewal of the cohort of active replication origins. We found that gained origins appear with low activity yet become more efficient and fire earlier as they evolutionarily age. By contrast, origins that are lost comprise the complete range of firing strength. Additionally, they preferentially occur in close vicinity to strong origins. Interestingly, despite high evolutionary turnover, active replication origins remain regularly spaced along chromosomes in all species, suggesting that origin distribution is optimized to limit large inter-origin intervals. We propose a model on the evolutionary birth, death, and conservation of active replication origins.

Temporal programs of genome replication show different levels of conservation between closely or distantly related species. Here, the authors generate genome-wide replication timing profiles for ten yeast species, and analyze their evolutionary dynamics.

Yeast heterochromatin regulators Sir2 and Sir3 act directly at euchromatic DNA replication origins

Most active DNA replication origins are found within euchromatin, while origins within heterochromatin are often inactive or inhibited. In yeast, origin activity within heterochromatin is negatively controlled by the histone H4K16 deacetylase, Sir2, and at some heterochromatic loci also by the nucleosome binding protein, Sir3. The prevailing view has been that direct functions of Sir2 and Sir3 are confined to heterochromatin. However, growth defects in yeast mutants compromised for loading the MCM helicase, such as cdc6-4, are suppressed by deletion of either SIR2 or SIR3. While these and other observations indicate that SIR2,3 can have a negative impact on at least some euchromatic origins, the genomic scale of this effect was unknown. It was also unknown whether this suppression resulted from direct functions of Sir2,3 within euchromatin, or was an indirect effect of their previously established roles within heterochromatin. Using MCM ChIP-Seq, we show that a SIR2 deletion rescued MCM complex loading at ~80% of euchromatic origins in cdc6-4 cells. Therefore, Sir2 exhibited a pervasive effect at the majority of euchromatic origins. Using MNase-H4K16ac ChIP-Seq, we show that origin-adjacent nucleosomes were depleted for H4K16 acetylation in a SIR2-dependent manner in wild type (i.e. CDC6) cells. In addition, we present evidence that both Sir2 and Sir3 bound to nucleosomes adjacent to euchromatic origins. The relative levels of each of these molecular hallmarks of yeast heterochromatin–SIR2-dependent H4K16 hypoacetylation, Sir2, and Sir3 –correlated with how strongly a SIR2 deletion suppressed the MCM loading defect in cdc6-4 cells. Finally, a screen for histone H3 and H4 mutants that could suppress the cdc6-4 growth defect identified amino acids that map to a surface of the nucleosome important for Sir3 binding. We conclude that heterochromatin proteins directly modify the local chromatin environment of euchromatic DNA replication origins.

When a cell divides, it must copy or “replicate” its DNA. DNA replication starts at chromosomal regions called origins when a collection of replication proteins gains local access to unwind the two DNA strands. Chromosomal DNA is packaged into a protein-DNA complex called chromatin and there are two major structurally and functionally distinct types. Euchromatin allows DNA replication proteins to access origin DNA, while heterochromatin inhibits their access. The prevalent view has been that the heterochromatin proteins required to inhibit origins are confined to heterochromatin. In this study, the conserved heterochromatin proteins, Sir2 and Sir3, were shown to both physically and functionally associate with the majority of origins in euchromatin. This observation raises important questions about the chromosomal targets of heterochromatin proteins, and how and why the majority of origins exist within a potentially repressive chromatin structure.