Prediction of Saccharomyces cerevisiae replication origins

Defining the replication program through the chromatin landscape

DNA replication is an essential cell cycle event required for the accurate and timely duplication of the chromosomes. It is essential that the genome is replicated accurately and completely within the confines of S-phase. Failure to completely copy the genome has the potential to result in catastrophic genomic instability. Replication initiates in a coordinated manner from multiple locations, termed origins of replication, distributed across each of the chromosomes. The selection of these origins of replication is a dynamic process responding to both developmental and tissue-specific signals. In this review, we explore the role of the local chromatin environment in regulating the DNA replication program at the level of origin selection and activation. Finally, there is increasing molecular evidence that the DNA replication program itself affects the chromatin landscape, suggesting that DNA replication is critical for both genetic and epigenetic inheritance.

Genome-scale analysis of metazoan replication origins reveals their organization in specific but flexible sites defined by conserved features

In metazoans, thousands of DNA replication origins (Oris) are activated at each cell cycle. Their genomic organization and their genetic nature remain elusive. Here, we characterized Oris by nascent strand (NS) purification and a genome-wide analysis in Drosophila and mouse cells. We show that in both species most CpG islands (CGI) contain Oris, although methylation is nearly absent in Drosophila, indicating that this epigenetic mark is not crucial for defining the activated origin. Initiation of DNA synthesis starts at the borders of CGI, resulting in a striking bimodal distribution of NS, suggestive of a dual initiation event. Oris contain a unique nucleotide skew around NS peaks, characterized by G/T and C/A overrepresentation at the 5' and 3' of Ori sites, respectively. Repeated GC-rich elements were detected, which are good predictors of Oris, suggesting that common sequence features are part of metazoan Oris. In the heterochromatic chromosome 4 of Drosophila, Oris correlated with HP1 binding sites. At the chromosome level, regions rich in Oris are early replicating, whereas Ori-poor regions are late replicating. The genome-wide analysis was coupled with a DNA combing analysis to unravel the organization of Oris. The results indicate that Oris are in a large excess, but their activation does not occur at random. They are organized in groups of site-specific but flexible origins that define replicons, where a single origin is activated in each replicon. This organization provides both site specificity and Ori firing flexibility in each replicon, allowing possible adaptation to environmental cues and cell fates.

Preferential localization of human origins of DNA replication at the 5'-ends of expressed genes and at evolutionarily conserved DNA sequences

BACKGROUND

Replication of mammalian genomes requires the activation of thousands of origins which are both spatially and temporally regulated by as yet unknown mechanisms. At the most fundamental level, our knowledge about the distribution pattern of origins in each of the chromosomes, among different cell types, and whether the physiological state of the cells alters this distribution is at present very limited.

METHODOLOGY/PRINCIPAL FINDINGS

We have used standard λ-exonuclease resistant nascent DNA preparations in the size range of 0.7-1.5 kb obtained from the breast cancer cell line MCF-7 hybridized to a custom tiling array containing 50-60 nt probes evenly distributed among genic and non-genic regions covering about 1% of the human genome. A similar DNA preparation was used for high-throughput DNA sequencing. Array experiments were also performed with DNA obtained from BT-474 and H520 cell lines. By determining the sites showing nascent DNA enrichment, we have localized several thousand origins of DNA replication. Our major findings are: (a) both array and DNA sequencing assay methods produced essentially the same origin distribution profile; (b) origin distribution is largely conserved (>70%) in all cell lines tested; (c) origins are enriched at the 5'ends of expressed genes and at evolutionarily conserved intergenic sequences; and (d) ChIP on chip experiments in MCF-7 showed an enrichment of H3K4Me3 and RNA Polymerase II chromatin binding sites at origins of DNA replication.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that the program for origin activation is largely conserved among different cell types. Also, our work supports recent studies connecting transcription initiation with replication, and in addition suggests that evolutionarily conserved intergenic sequences have the potential to participate in origin selection. Overall, our observations suggest that replication origin selection is a stochastic process significantly dependent upon local accessibility to replication factors.

High-resolution analysis of four efficient yeast replication origins reveals new insights into the ORC and putative MCM binding elements

In budding yeast, the eukaryotic initiator protein ORC (origin recognition complex) binds to a bipartite sequence consisting of an 11 bp ACS element and an adjacent B1 element. However, the genome contains many more matches to this consensus than actually bind ORC or function as origins in vivo. Although ORC-dependent loading of the replicative MCM helicase at origins is enhanced by a distal B2 element, less is known about this element. Here, we analyzed four highly active origins (ARS309, ARS319, ARS606 and ARS607) by linker scanning mutagenesis and found that sequences adjacent to the ACS contributed substantially to origin activity and ORC binding. Using the sequences of four additional B2 elements we generated a B2 multiple sequence alignment and identified a shared, degenerate 8 bp sequence that was enriched within 228 known origins. In addition, our high-resolution analysis revealed that not all origins exist within nucleosome free regions: a class of Sir2-regulated origins has a stably positioned nucleosome overlapping or near B2. This study illustrates the conserved yet flexible nature of yeast origin architecture to promote ORC binding and origin activity, and helps explain why a strong match to the ORC binding site is insufficient to identify origins within the genome.

Nuclear mitochondrial DNA activates replication in Saccharomyces cerevisiae

The nuclear genome of eukaryotes is colonized by DNA fragments of mitochondrial origin, called NUMTs. These insertions have been associated with a variety of germ-line diseases in humans. The significance of this uptake of potentially dangerous sequences into the nuclear genome is unclear. Here we provide functional evidence that sequences of mitochondrial origin promote nuclear DNA replication in Saccharomyces cerevisiae. We show that NUMTs are rich in key autonomously replicating sequence (ARS) consensus motifs, whose mutation results in the reduction or loss of DNA replication activity. Furthermore, 2D-gel analysis of the mrc1 mutant exposed to hydroxyurea shows that several NUMTs function as late chromosomal origins. We also show that NUMTs located close to or within ARS provide key sequence elements for replication. Thus NUMTs can act as independent origins, when inserted in an appropriate genomic context or affect the efficiency of pre-existing origins. These findings show that migratory mitochondrial DNAs can impact on the replication of the nuclear region they are inserted in.

Chromatin signatures of the Drosophila replication program

DNA replication initiates from thousands of start sites throughout the Drosophila genome and must be coordinated with other ongoing nuclear processes such as transcription to ensure genetic and epigenetic inheritance. Considerable progress has been made toward understanding how chromatin modifications regulate the transcription program; in contrast, we know relatively little about the role of the chromatin landscape in defining how start sites of DNA replication are selected and regulated. Here, we describe the Drosophila replication program in the context of the chromatin and transcription landscape for multiple cell lines using data generated by the modENCODE consortium. We find that while the cell lines exhibit similar replication programs, there are numerous cell line-specific differences that correlate with changes in the chromatin architecture. We identify chromatin features that are associated with replication timing, early origin usage, and ORC binding. Primary sequence, activating chromatin marks, and DNA-binding proteins (including chromatin remodelers) contribute in an additive manner to specify ORC-binding sites. We also generate accurate and predictive models from the chromatin data to describe origin usage and strength between cell lines. Multiple activating chromatin modifications contribute to the function and relative strength of replication origins, suggesting that the chromatin environment does not regulate origins of replication as a simple binary switch, but rather acts as a tunable rheostat to regulate replication initiation events.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.

A comprehensive genome-wide map of autonomously replicating sequences in a naive genome

Eukaryotic chromosomes initiate DNA synthesis from multiple replication origins. The machinery that initiates DNA synthesis is highly conserved, but the sites where the replication initiation proteins bind have diverged significantly. Functional comparative genomics is an obvious approach to study the evolution of replication origins. However, to date, the Saccharomyces cerevisiae replication origin map is the only genome map available. Using an iterative approach that combines computational prediction and functional validation, we have generated a high-resolution genome-wide map of DNA replication origins in Kluyveromyces lactis. Unlike other yeasts or metazoans, K. lactis autonomously replicating sequences (KlARSs) contain a 50 bp consensus motif suggestive of a dimeric structure. This motif is necessary and largely sufficient for initiation and was used to dependably identify 145 of the up to 156 non-repetitive intergenic ARSs projected for the K. lactis genome. Though similar in genome sizes, K. lactis has half as many ARSs as its distant relative S. cerevisiae. Comparative genomic analysis shows that ARSs in K. lactis and S. cerevisiae preferentially localize to non-syntenic intergenic regions, linking ARSs with loci of accelerated evolutionary change.

Conserved nucleosome positioning defines replication origins

The origin recognition complex (ORC) specifies replication origin location. The Saccharomyces cerevisiae ORC recognizes the ARS (autonomously replicating sequence) consensus sequence (ACS), but only a subset of potential genomic sites are bound, suggesting other chromosomal features influence ORC binding. Using high-throughput sequencing to map ORC binding and nucleosome positioning, we show that yeast origins are characterized by an asymmetric pattern of positioned nucleosomes flanking the ACS. The origin sequences are sufficient to maintain a nucleosome-free origin; however, ORC is required for the precise positioning of nucleosomes flanking the origin. These findings identify local nucleosomes as an important determinant for origin selection and function.

Structural diversity and dynamics of genomic replication origins in Schizosaccharomyces pombe

DNA replication origins (ORI) in Schizosaccharomyces pombe colocalize with adenine and thymine (A+T)-rich regions, and earlier analyses have established a size from 0.5 to over 3 kb for a DNA fragment to drive replication in plasmid assays. We have asked what are the requirements for ORI function in the chromosomal context. By designing artificial ORIs, we have found that A+T-rich fragments as short as 100 bp without homology to S. pombe DNA are able to initiate replication in the genome. On the other hand, functional dissection of endogenous ORIs has revealed that some of them span a few kilobases and include several modules that may be as short as 25-30 contiguous A+Ts capable of initiating replication from ectopic chromosome positions. The search for elements with these characteristics across the genome has uncovered an earlier unnoticed class of low-efficiency ORIs that fire late during S phase. These results indicate that ORI specification and dynamics varies widely in S. pombe, ranging from very short elements to large regions reminiscent of replication initiation zones in mammals.

The right half of the Escherichia coli replication origin is not essential for viability, but facilitates multi-forked replication

Replication initiation is a key event in the cell cycle of all organisms and oriC, the replication origin in Escherichia coli, serves as the prototypical model for this process. The minimal sequence required for oriC function was originally determined entirely from plasmid studies using cloned origin fragments, which have previously been shown to differ dramatically in sequence requirement from the chromosome. Using an in vivo recombineering strategy to exchange wt oriCs for mutated ones regardless of whether they are functional origins or not, we have determined the minimal origin sequence that will support chromosome replication. Nearly the entire right half of oriC could be deleted without loss of origin function, demanding a reassessment of existing models for initiation. Cells carrying the new DnaA box-depleted 163 bp minimal oriC exhibited little or no loss of fitness under slow-growth conditions, but were sensitive to rich medium, suggesting that the dense packing of initiator binding sites that is a hallmark of prokaryotic origins, has likely evolved to support the increased demands of multi-forked replication.

An ensemble model of competitive multi-factor binding of the genome

Hundreds of different factors adorn the eukaryotic genome, binding to it in large number. These DNA binding factors (DBFs) include nucleosomes, transcription factors (TFs), and other proteins and protein complexes, such as the origin recognition complex (ORC). DBFs compete with one another for binding along the genome, yet many current models of genome binding do not consider different types of DBFs together simultaneously. Additionally, binding is a stochastic process that results in a continuum of binding probabilities at any position along the genome, but many current models tend to consider positions as being either binding sites or not. Here, we present a model that allows a multitude of DBFs, each at different concentrations, to compete with one another for binding sites along the genome. The result is an "occupancy profile," a probabilistic description of the DNA occupancy of each factor at each position. We implement our model efficiently as the software package COMPETE. We demonstrate genome-wide and at specific loci how modeling nucleosome binding alters TF binding, and vice versa, and illustrate how factor concentration influences binding occupancy. Binding cooperativity between nearby TFs arises implicitly via mutual competition with nucleosomes. Our method applies not only to TFs, but also recapitulates known occupancy profiles of a well-studied replication origin with and without ORC binding. Importantly, the sequence preferences our model takes as input are derived from in vitro experiments. This ensures that the calculated occupancy profiles are the result of the forces of competition represented explicitly in our model and the inherent sequence affinities of the constituent DBFs.

Subtelomeric ACS-containing proto-silencers act as antisilencers in replication factors mutants in Saccharomyces cerevisiae

Subtelomeric genes are either fully active or completely repressed and can switch their state about once per 20 generations. This meta-stable telomeric position effect is mediated by strong repression signals emitted by the telomere and relayed/enhanced by weaker repressor elements called proto-silencers. In addition, subtelomeric regions contain sequences with chromatin partitioning and antisilencing activities referred to as subtelomeric antisilencing regions. Using extensive mutational analysis of subtelomeric elements, we show that ARS consensus sequence (ACS)-containing proto-silencers convert to antisilencers in several replication factor mutants. We point out the significance of the B1 auxiliary sequence next to ACS in mediating these effects. In contrast, an origin-derived ACS does not convert to antisilencer in mutants and its B1 element has little bearing on silencing. These results are specific for the analyzed ACS and in addition to the effects of each mutation (relative to wild type) on global silencing. Another line of experiments shows that Mcm5p possesses antisilencing activity and is recruited to telomeres in an ACS-dependent manner. Mcm5p persists at this location at the late stages of S phase. We propose that telomeric ACS are not static proto-silencers but conduct finely tuned silencing and antisilencing activities mediated by ACS-bound factors.

Computational detection of significant variation in binding affinity across two sets of sequences with application to the analysis of replication origins in yeast

Background

In analyzing the stability of DNA replication origins in Saccharomyces cerevisiae we faced the question whether one set of sequences is significantly enriched in the number and/or the quality of the matches of a particular position weight matrix relative to another set.

Results

We present SADMAMA, a computational solution to a address this problem. SADMAMA implements two types of statistical tests to answer this question: one type is based on simplified models, while the other relies on bootstrapping, and as such might be preferable to users who are averse to such models. The bootstrap approach incorporates a novel "site-protected" resampling procedure which solves a problem we identify with naive resampling.

Conclusion

SADMAMA's utility is demonstrated here by offering a plausible explanation to the differential ARS activity observed in our previous mcm1-1 mutant experiments [1], by suggesting the relevance of multiple weak ACS matches to efficient replication origin function in Saccharomyces cerevisiae, and by suggesting an explanation to the observed negative effect FKH2 has on chromatin silencing [2]. SADMAMA is available for download from .

Multiple mechanisms contribute to Schizosaccharomyces pombe origin recognition complex-DNA interactions

Eukaryotic DNA replication requires the assembly of multiprotein pre-replication complexes (pre-RCs) at chromosomal origins of DNA replication. Here we describe the interactions of highly purified Schizosaccharomyces pombe pre-RC components, SpORC, SpCdc18, and SpCdt1, with each other and with ars1 origin DNA. We show that SpORC binds DNA in at least two steps. The first step likely involves electrostatic interactions between the AT-hook motifs of SpOrc4 and AT tracts in ars1 DNA and results in the formation of a salt-sensitive complex. In the second step, the salt-sensitive complex is slowly converted to a salt-stable complex that involves additional interactions between SpORC and DNA. Binding of SpORC to ars1 DNA is facilitated by negative supercoiling and is accompanied by changes in DNA topology, suggesting that SpORC-DNA complexes contain underwound or negatively writhed DNA. Purified human origin recognition complex (ORC) induces similar topological changes in origin DNA, indicating that this property of ORC is conserved in eukaryotic evolution and plays an important role in ORC function. We also show that SpCdc18 and SpCdt1 form a binary complex that has greater affinity for DNA than either protein alone. In addition, both proteins contribute significantly to the stability of the initial SpORC-DNA complex and enhance the SpORC-dependent topology changes in origin DNA. Thus, the formation of stable protein-DNA complexes at S. pombe origins of replication involves binary interactions among all three proteins, as well as interactions of both SpORC and SpCdt1-SpCdc18 with origin DNA. These findings demonstrate that SpORC is not the sole determinant of origin recognition.

Analysis of chromosome III replicators reveals an unusual structure for the ARS318 silencer origin and a conserved WTW sequence within the origin recognition complex binding site

Saccharomyces cerevisiae chromosome III encodes 11 autonomously replicating sequence (ARS) elements that function as chromosomal replicators. The essential 11-bp ARS consensus sequence (ACS) that binds the origin recognition complex (ORC) has been experimentally defined for most of these replicators but not for ARS318 (HMR-I), which is one of the HMR silencers. In this study, we performed a comprehensive linker scan analysis of ARS318. Unexpectedly, this replicator depends on a 9/11-bp match to the ACS that positions the ORC binding site only 6 bp away from an Abf1p binding site. Although a largely inactive replicator on the chromosome, ARS318 becomes active if the nearby HMR-E silencer is deleted. We also performed a multiple sequence alignment of confirmed replicators on chromosomes III, VI, and VII. This analysis revealed a highly conserved WTW motif 17 to 19 bp from the ACS that is functionally important and is apparent in the 228 phylogenetically conserved ARS elements among the six sensu stricto Saccharomyces species.

Global regulation of genome duplication in eukaryotes: an overview from the epifluorescence microscope

In eukaryotes, DNA replication is initiated along each chromosome at multiple sites called replication origins. Locally, each replication origin is "licensed" or specified at the end of the M and the beginning of the G1 phases of the cell cycle. During the S phase when DNA synthesis takes place, origins are activated in stages corresponding to early and late-replicating domains. The staged and progressive activation of replication origins reflects the need to maintain a strict balance between the number of active replication forks and the rate at which DNA synthesis proceeds. This suggests that origin densities (frequency of initiation) and replication fork movement (rates of elongation) must be coregulated to guarantee the efficient and complete duplication of each subchromosomal domain. Emerging evidence supports this proposal and suggests that the ATM/ATR intra-S phase checkpoint plays an important role in the coregulation of initiation frequencies and rates of elongation. In this paper, we review recent results concerning the mechanisms governing the global regulation of DNA replication and discuss the roles these mechanisms play in maintaining genome stability during both a normal and perturbed S phase.

Analysis of genetic interactions on a genome-wide scale in budding yeast: diploid-based synthetic lethality analysis by microarray

Comprehensive collections of open reading frame (ORF) deletion mutant strains exist for the budding yeast Saccharomyces cerevisiae. With great prescience, these strains were designed with short molecular bar codes or TAGs that uniquely mark each deletion allele, flanked by shared priming sequences. These features have enabled researchers to handle yeast mutant collections as complex pools of approximately 6000 strains. The presence of any individual mutant within a pool can be assessed indirectly by measuring the relative abundance of its corresponding TAG(s) in genomic DNA prepared from the pool. This is readily accomplished by wholesale polymerase chain reaction (PCR) amplification of the TAGs using fluorescent oligonucleotide primers that recognize the common flanking sequences, followed by hybridization of the labeled PCR products to a TAG oligonucleotide microarray. Here we describe a method-diploid-based synthetic lethality analysis by microarray (dSLAM)-whereby such pools can be manipulated to rapidly construct and assess the fitness of 6000 double-mutant strains in a single experiment. Analysis of double-mutant strains is of growing importance in defining the spectrum of essential cellular functionalities and in understanding how these functionalities interrelate.

AT excursion: a new approach to predict replication origins in viral genomes by locating AT-rich regions

Background

Replication origins are considered important sites for understanding the molecular mechanisms involved in DNA replication. Many computational methods have been developed for predicting their locations in archaeal, bacterial and eukaryotic genomes. However, a prediction method designed for a particular kind of genomes might not work well for another. In this paper, we propose the AT excursion method, which is a score-based approach, to quantify local AT abundance in genomic sequences and use the identified high scoring segments for predicting replication origins. This method has the advantages of requiring no preset window size and having rigorous criteria to evaluate statistical significance of high scoring segments.

Results

We have evaluated the AT excursion method by checking its predictions against known replication origins in herpesviruses and comparing its performance with an existing base weighted score method (BWS₁). Out of 43 known origins, 39 are predicted by either one or the other method and 26 origins are predicted by both. The excursion method identifies six origins not predicted by BWS₁, showing that the AT excursion method is a valuable complement to BWS₁. We have also applied the AT excursion method to two other families of double stranded DNA viruses, the poxviruses and iridoviruses, of which very few replication origins are documented in the public domain. The prediction results are made available as supplementary materials at [1]. Preliminary investigation shows that the proposed method works well on some larger genomes too.

Conclusion

The AT excursion method will be a useful computational tool for identifying replication origins in a variety of genomic sequences.

Genome-wide mapping of ORC and Mcm2p binding sites on tiling arrays and identification of essential ARS consensus sequences in S. cerevisiae

BACKGROUND

Eukaryotic replication origins exhibit different initiation efficiencies and activation times within S-phase. Although local chromatin structure and function influences origin activity, the exact mechanisms remain poorly understood. A key to understanding the exact features of chromatin that impinge on replication origin function is to define the precise locations of the DNA sequences that control origin function. In S. cerevisiae, Autonomously Replicating Sequences (ARSs) contain a consensus sequence (ACS) that binds the Origin Recognition Complex (ORC) and is essential for origin function. However, an ACS is not sufficient for origin function and the majority of ACS matches do not function as ORC binding sites, complicating the specific identification of these sites.

RESULTS

To identify essential origin sequences genome-wide, we utilized a tiled oligonucleotide array (NimbleGen) to map the ORC and Mcm2p binding sites at high resolution. These binding sites define a set of potential Autonomously Replicating Sequences (ARSs), which we term nimARSs. The nimARS set comprises 529 ORC and/or Mcm2p binding sites, which includes 95% of known ARSs, and experimental verification demonstrates that 94% are functional. The resolution of the analysis facilitated identification of potential ACSs (nimACSs) within 370 nimARSs. Cross-validation shows that the nimACS predictions include 58% of known ACSs, and experimental verification indicates that 82% are essential for ARS activity.

CONCLUSION

These findings provide the most comprehensive, accurate, and detailed mapping of ORC binding sites to date, adding to the emerging picture of the chromatin organization of the budding yeast genome.

The spatial arrangement of ORC binding modules determines the functionality of replication origins in budding yeast

In the quest to define autonomously replicating sequences (ARSs) in eukaryotic cells, an ARS consensus sequence (ACS) has emerged for budding yeast. This ACS is recognized by the replication initiator, the origin recognition complex (ORC). However, not every match to the ACS constitutes a replication origin. Here, we investigated the requirements for ORC binding to origins that carry multiple, redundant ACSs, such as ARS603. Previous studies raised the possibility that these ACSs function as individual ORC binding sites. Detailed mutational analysis of the two ACSs in ARS603 revealed that they function in concert and give rise to an initiation pattern compatible with a single bipartite ORC binding site. Consistent with this notion, deletion of one base pair between the ACS matches abolished ORC binding at ARS603. Importantly, loss of ORC binding in vitro correlated with the loss of ARS activity in vivo. Our results argue that replication origins in yeast are in general comprised of bipartite ORC binding sites that cannot function in random alignment but must conform to a configuration that permits ORC binding. These requirements help to explain why only a limited number of ACS matches in the yeast genome qualify as ORC binding sites.

Genome-wide identification of replication origins in yeast by comparative genomics

We discovered that sequences essential for replication origin function are frequently conserved in sensu stricto Saccharomyces species. Here we use analysis of phylogenetic conservation to identify replication origin sequences throughout the Saccharomyces cerevisiae genome at base pair resolution. Origin activity was confirmed for each of 228 predicted sites--representing 86% of apparent origin regions. This is the first study to determine the genome-wide location of replication origins at a resolution sufficient to identify the sequence elements bound by replication proteins. Our results demonstrate that phylogenetic conservation can be used to identify the origin sequences responsible for replicating a eukaryotic genome.

Asymmetric positioning of nucleosomes and directional establishment of transcriptionally silent chromatin by Saccharomyces cerevisiae silencers

In Saccharomyces cerevisiae, silencers flanking the HML and HMR loci consist of various combinations of binding sites for Abf1p, Rap1p, and the origin recognition complex (ORC) that serve to recruit the Sir silencing complex, thereby initiating the establishment of transcriptionally silent chromatin. There have been seemingly conflicting reports concerning whether silencers function in an orientation-dependent or -independent manner, and what determines the directionality of a silencer has not been explored. We demonstrate that chromatin plays a key role in determining the potency and directionality of silencers. We show that nucleosomes are asymmetrically distributed around the HML-I or HMR-E silencer so that a nucleosome is positioned close to the Abf1p side but not the ORC side of the silencer. This coincides with preferential association of Sir proteins and transcriptional silencing on the Abf1p side of the silencer. Elimination of the asymmetry in nucleosome positioning at a silencer leads to comparable silencing on both sides. Asymmetric nucleosome positioning in the immediate vicinity of a silencer is independent of its orientation and genomic context, indicating that it is the inherent property of the silencer. Moreover, it is also independent of the Sir complex and thus precedes the formation of silent chromatin. Finally, we demonstrate that asymmetric positioning of nucleosomes and directional silencing by a silencer depend on ORC and Abf1p. We conclude that the HML-I and HMR-E silencers promote asymmetric positioning of nucleosomes, leading to unequal potentials of transcriptional silencing on their sides and, hence, directional silencing.

Genome-wide hierarchy of replication origin usage in Saccharomyces cerevisiae

Replication origins in a genome are inherently different in their base sequence and in their response to temporal and cell cycle regulation signals for DNA replication. To investigate the chromosomal determinants that influence the efficiency of initiation of DNA replication genome-wide, we made use of a reverse strategy originally used for the isolation of replication initiation mutants in Saccharomyces cerevisiae. In yeast, replication origins isolated from chromosomes support the autonomous replication of plasmids. These replication origins, whether in the context of a chromosome or a plasmid, will initiate efficiently in wild-type cells but show a dramatically contrasted efficiency of activation in mutants defective in the early steps of replication initiation. Serial passages of a genomic library of autonomously replicating sequences (ARSs) in such a mutant allowed us to select for constitutively active ARSs. We found a hierarchy of preferential initiation of ARSs that correlates with local transcription patterns. This preferential usage is enhanced in mutants defective in the assembly of the prereplication complex (pre-RC) but not in mutants defective in the activation of the pre-RC. Our findings are consistent with an interference of local transcription with the assembly of the pre-RC at a majority of replication origins.

The length of S phase regulated by the rate of DNA synthesis varies dramatically during the development of metazoans. Key to this regulation is the number of replication origins utilized in different developmental stages. A fundamental question is whether there is a hierarchy in the usage of replication origins under different conditions and if so, what are the determinants for preferential usage. In Saccharomyces cerevisiae, replication origins isolated in DNA fragments are known as autonomously replicating sequences (ARSs). To gain insight into the determinants that regulate replication origin usage, genomic ARSs that are preferentially used under adverse conditions for replication initiation were identified. One of the determinants appears to be the local transcription pattern. Transcriptional activity directed towards an ARS correlates with reduced efficiency of replication initiation of that ARS. This transcriptional interference appears to be targeted at the assembly of the prereplication complex. These results are consistent with the deregulated initiation patterns observed in early developing Xenopus embryos that are devoid of transcription. Other yet-to-be-identified factors are also important in determining the efficiency of replication origin usage.

Genome-wide analysis of re-replication reveals inhibitory controls that target multiple stages of replication initiation

DNA replication must be tightly controlled during each cell cycle to prevent unscheduled replication and ensure proper genome maintenance. The currently known controls that prevent re-replication act redundantly to inhibit pre-replicative complex (pre-RC) assembly outside of the G1-phase of the cell cycle. The yeast Saccharomyces cerevisiae has been a useful model organism to study how eukaryotic cells prevent replication origins from reinitiating during a single cell cycle. Using a re-replication-sensitive strain and DNA microarrays, we map sites across the S. cerevisiae genome that are re-replicated as well as sites of pre-RC formation during re-replication. Only a fraction of the genome is re-replicated by a subset of origins, some of which are capable of multiple reinitiation events. Translocation experiments demonstrate that origin-proximal sequences are sufficient to predispose an origin to re-replication. Origins that reinitiate are largely limited to those that can recruit Mcm2-7 under re-replicating conditions; however, the formation of a pre-RC is not sufficient for reinitiation. Our findings allow us to categorize origins with respect to their propensity to reinitiate and demonstrate that pre-RC formation is not the only target for the mechanisms that prevent genomic re-replication.

Positioning and the specific sequence of each 13-mer motif are critical for activity of the plasmid RK2 replication origin

The minimal replication origin of the broad-host-range plasmid RK2, oriV, contains five iterons which are binding sites for the plasmid-encoded replication initiation protein TrfA, four DnaA boxes, which bind the host DnaA protein, and an AT-rich region containing four 13-mer sequences. In this study, 26 mutants with altered sequence and/or spacing of 13-mer motifs have been constructed and analysed for replication activity in vivo and in vitro. The data show that the replacement of oriV 13-mers by similar but not identical 13-mer sequences from Escherichia coli oriC inactivates the origin. In addition, interchanging the positions of the oriV 13-mers results in greatly reduced activity. Mutants with T/A substitutions are also inactive. Furthermore, introduction of single-nucleotide substitutions demonstrates very restricted sequence requirements depending on the 13-mer position. Only two of the mutants are host specific, functional in Pseudomonas aeruginosa but not in E. coli. Our experiments demonstrate considerable complexity in the plasmid AT-rich region architecture required for functionality. It is evident that low internal stability of this region is not the only feature contributing to origin activity. Our studies suggest a requirement for sequence-specific protein interactions within the 13-mers during assembly of replication complexes at the plasmid origin.

Scoring schemes of palindrome clusters for more sensitive prediction of replication origins in herpesviruses

Many empirical studies show that there are unusual clusters of palindromes, closely spaced direct and inverted repeats around the replication origins of herpesviruses. In this paper, we introduce two new scoring schemes to quantify the spatial abundance of palindromes in a genomic sequence. Based on these scoring schemes, a computational method to predict the locations of replication origins is developed. When our predictions are compared with 39 known or annotated replication origins in 19 herpesviruses, close to 80% of the replication origins are located within 2% of the genome length. A list of predicted locations of replication origins in all the known herpesviruses with complete genome sequences is reported.

Control of replication initiation and heterochromatin formation in Saccharomyces cerevisiae by a regulator of meiotic gene expression

Heterochromatinization at the silent mating-type loci HMR and HML in Saccharomyces cerevisiae is achieved by targeting the Sir complex to these regions via a set of anchor proteins that bind to the silencers. Here, we have identified a novel heterochromatin-targeting factor for HML, the protein Sum1, a repressor of meiotic genes during vegetative growth. Sum1 bound both in vitro and in vivo to HML via a functional element within the HML-E silencer, and sum1Delta caused HML derepression. Significantly, Sum1 was also required for origin activity of HML-E, demonstrating a role of Sum1 in replication initiation. In a genome-wide search for Sum1-regulated origins, we identified a set of autonomous replicative sequences (ARS elements) that bound both the origin recognition complex and Sum1. Full initiation activity of these origins required Sum1, and their origin activity was decreased upon removal of the Sum1-binding site. Thus, Sum1 constitutes a novel global regulator of replication initiation in yeast.

Susceptibility to superhelically driven DNA duplex destabilization: a highly conserved property of yeast replication origins

Strand separation is obligatory for several DNA functions, including replication. However, local DNA properties such as A+T content or thermodynamic stability alone do not determine the susceptibility to this transition in vivo. Rather, superhelical stresses provide long-range coupling among the transition behaviors of all base pairs within a topologically constrained domain. We have developed methods to analyze superhelically induced duplex destabilization (SIDD) in genomic DNA that take into account both this long-range stress-induced coupling and sequence-dependent local thermodynamic stability. Here we apply this approach to examine the SIDD properties of 39 experimentally well-characterized autonomously replicating DNA sequences (ARS elements), which function as replication origins in the yeast Saccharomyces cerevisiae. We find that these ARS elements have a strikingly increased susceptibility to SIDD relative to their surrounding sequences. On average, these ARS elements require 4.78 kcal/mol less free energy to separate than do their immediately surrounding sequences, making them more than 2,000 times easier to open. Statistical analysis shows that the probability of this strong an association between SIDD sites and ARS elements arising by chance is approximately 4 × 10⁻¹⁰. This local enhancement of the propensity to separate to single strands under superhelical stress has obvious implications for origin function. SIDD properties also could be used, in conjunction with other known origin attributes, to identify putative replication origins in yeast, and possibly in other metazoan genomes.

Several DNA functions require the two strands of the DNA duplex to transiently separate. Examples include the initiation of gene expression and of DNA replication. Here the authors examine the strand separation properties of the DNA duplex at autonomously replicating sequences (ARS elements), which are the potential replication origins in yeast.

In vivo, susceptibility to strand separation does not depend only on local DNA properties such as adenine plus thymine content or thermodynamic stability. Rather, stresses imposed on the DNA in vivo couple together the strand-opening behaviors of all base pairs that experience them. The authors use computational methods for analyzing stress-driven strand separation to examine the susceptibility to opening of 39 experimentally well-characterized ARS elements. They show that these ARS elements have strikingly increased susceptibilities to stress-induced separation relative to the surrounding sequences. On average, these ARS elements require 4.78 kcal/mol less free energy to separate than do surrounding sequences, making them more than 2,000 times easier to open. This enhanced susceptibility to stress-driven strand separation has obvious implications for the mechanisms that begin the process of replication. This property is also shared by bacterial and viral replication start points, suggesting that it may be a general attribute of replication origins.

Identification of cis-acting elements that mediate the replication and maintenance of human papillomavirus type 16 genomes in Saccharomyces cerevisiae

Papillomaviruses contain small double-stranded DNA genomes that are maintained in persistently infected mammalian host epithelia as nuclear plasmids and rely upon the host replication machinery for replication. Papillomaviruses encode a DNA helicase, E1, which can specifically bind to the viral genome and support DNA synthesis. Under some conditions in mammalian cells, E1 is not required for viral DNA synthesis, leading to the hypothesis that papillomavirus DNA can be replicated solely by the host replication machinery. This machinery is highly conserved among eukaryotes. We and others found that papillomavirus DNA could replicate in a simple eukaryote, Saccharomyces cerevisiae. Specifically, papillomavirus DNA could substitute for the function of the autonomously replicating sequence (ARS) and centromere (CEN) elements that are normally both required for the stable replication of extrachromosomal DNAs in yeast. Furthermore, this form of replication in yeast was E1 independent. In this study, we map the elements in the human papillomavirus type 16 (HPV16) genome that can substitute for yeast ARS and CEN elements. A single element, termed rep, was identified that can substitute for ARS, and multiple elements, termed mtc, could substitute for CEN. The location of one of these mtc elements overlaps the location of rep, and this approximately 1,000-bp region of HPV16 was sufficient to support stable replication of a bacterial-yeast shuttle plasmid deleted of both ARS and CEN elements.

Independence of replisomes in Escherichia coli chromosomal replication

In Escherichia coli DNA replication is carried out by the coordinated action of the proteins within a replisome. After replication initiation, the two bidirectionally oriented replisomes from a single origin are colocalized into higher-order structures termed replication factories. The factory model postulated that the two replisomes are also functionally coupled. We tested this hypothesis by using DNA combing and whole-genome microarrays. Nascent DNA surrounding oriC in single, combed chromosomes showed instead that one replisome, usually the leftward one, was significantly ahead of the other 70% of the time. We next used microarrays to follow replication throughout the genome by measuring DNA copy number. We found in multiple E. coli strains that the replisomes are independent, with the leftward replisome ahead of the rightward one. The size of the bias was strain-specific, varying from 50 to 130 kb in the array results. When we artificially blocked one replisome, the other continued unabated, again demonstrating independence. We suggest an improved version of the factory model that retains the advantages of threading DNA through colocalized replisomes at about equal rates, but allows the cell flexibility to overcome obstacles encountered during elongation.

Flexibility and governance in eukaryotic DNA replication

Eukaryotic DNA replication begins at numerous but often poorly characterized sequences called origins, which are distributed fairly regularly along chromosomes. The elusive and idiosyncratic nature of origins in higher eukaryotes is now understood as resulting from a strong epigenetic influence on their specification, which provides flexibility in origin selection and allows for tailoring the dynamics of chromosome replication to the specific needs of cells. By contrast, the factors that assemble in trans to make these origins competent for replication and the kinases that trigger initiation are well conserved. Genome-wide and single-molecule approaches are being developed to elucidate the dynamics of chromosome replication. The notion that a well-coordinated progression of replication forks is crucial for many aspects of the chromosome cycle besides simply duplication begins to be appreciated.

DNA replication origins in the Schizosaccharomyces pombe genome

Origins of DNA replication in Schizosaccharomyces pombe lack a specific consensus sequence analogous to the Saccharomyces cerevisiae autonomously replicating sequence (ARS) consensus, raising the question of how they are recognized by the replication machinery. Because all well characterized S. pombe origins are located in intergenic regions, we analyzed the sequence properties and biological activity of such regions. The AT content of intergenes is very high ( approximately 70%), and runs of A's or T's occur with a significantly greater frequency than expected. Additionally, the two DNA strands in intergenes display compositional asymmetry that strongly correlates with the direction of transcription of flanking genes. Importantly, the sequence properties of known S. pombe origins of DNA replication are similar to those of intergenes in general. In functional studies, we assayed the in vivo origin activity of 26 intergenes in a 68-kb region of S. pombe chromosome 2. We also assayed the origin activity of sets of randomly chosen intergenes with the same length or AT content. Our data demonstrate that at least half of intergenes have potential origin activity and that the relative ability of an intergene to function as an origin is governed primarily by AT content and length. We propose a stochastic model for initiation of DNA replication in the fission yeast. In this model, the number of AT tracts in a given sequence is the major determinant of its probability of binding SpORC and serving as a replication origin. A similar model may explain some features of origins of DNA replication in metazoans.

In search of the holy replicator

After 40 years of searching for the eukaryotic replicator sequence, it is time to abandon the concept of 'the' replicator as a single genetic entity. Here I propose a 'relaxed replicon model' in which a positive initiator-replicator interaction is facilitated by a combination of several complex features of chromatin. An important question for the future is whether the positions of replication origins are simply a passive result of local chromatin structure or are actively localized to coordinate replication with other chromosomal activities.

Genomic specification and epigenetic regulation of eukaryotic DNA replication origins

Identification of DNA replication origins (ORIs) at a genome-wide level in eukaryotes has proved to be difficult due to the high degree of degeneracy of their sequences. Recent structural and functional approaches, however, have circumvented this limitation and have provided reliable predictions of their genomic distribution in the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, and they have also significantly increased the number of characterized ORIs in animals. This article reviews recent evidence on how ORIs are specified and maintained in these systems and on their regulation and sensitivity to epigenetic signals. It also discusses the possible additional involvement of ORIs in processes other than DNA replication.

Initiation sites for human DNA replication at a putative ribulose-5-phosphate 3-epimerase gene

Replication of the human genome requires the activation of thousands of replicons distributed along each one of the chromosomes. Each replicon contains an initiation, or origin, site, at which DNA synthesis begins. However, very little information is known about the nature and positioning of these initiation sites along human chromosomes. We have recently focused our attention to a 1.1 kb region of human chromosome 2 which functioned as an episomal origin in the yeast Saccharomyces cerevisiae. This region corresponded to the largest exon of a putative ribulose-5-phosphate-3-epimerase gene (RPE). In the present study we have used a real-time PCR-based nascent strand DNA abundance assay to map initiation sites for DNA replication in in vivo human chromosomes around a 13.4 kb region encompassing the putative RPE gene. By applying this analysis to a 1-1.4 kb nascent strand DNA fraction isolated from both normal skin fibroblasts, and the breast cell line MCF10; we have identified five initiation sites within the 13.4 kb region of chromosome 2. The initiation sites appear to map to similar positions in both cell lines and occur outside the coding regions of the putative RPE gene.