Search

Schizosaccharomyces japonicus belongs to the single-genus class Schizosaccharomycetes, otherwise known as "fission yeasts." As part of a composite model system with its widely studied S. pombe sister species, S. japonicus has provided critical insights into the workings and the evolution of cell biological mechanisms. Furthermore, its divergent biology makes S. japonicus a valuable model organism in its own right. However, the currently available genome assembly contains gaps and has been unable to resolve centromeres and other repeat-rich chromosomal regions. Here we present a telomere-to-telomere long-read genome assembly of the S. japonicus genome. This includes the three megabase-length chromosomes, with centromeres hundreds of kilobases long, rich in 5S ribosomal RNA genes, transfer RNA genes, long terminal repeats, and short repeats. We identify a gene-sparse region on chromosome 2 that resembles a 331 kb centromeric duplication. We revise the genome size of S. japonicus to at least 16.6 Mb and possibly up to 18.12 Mb, at least 30% larger than previous estimates. Our whole genome assembly will support the growing S. japonicus research community and facilitate research in new directions, including centromere and DNA repeat evolution, and yeast comparative genomics., (© 2024 The Authors. Yeast published by John Wiley & Sons Ltd.)

The fission yeast species Schizosaccharomyces japonicus is currently divided into two varieties-S. japonicus var. japonicus and S. japonicus var. versatilis. Here we examine the var. versatilis isolate CBS5679. The CBS5679 genome shows 88% identity to the reference genome of S. japonicus var. japonicus at the coding sequence level, with phylogenetic analyses suggesting that it has split from the S. japonicus lineage 25 million years ago. The CBS5679 genome contains a reciprocal translocation between chromosomes 1 and 2, together with several large inversions. The products of genes linked to the major translocation are associated with 'metabolism' and 'cellular assembly' ontology terms. We further show that CBS5679 does not generate viable progeny with the reference strain of S. japonicus. Although CBS5679 shares closer similarity to the 'type' strain of var. versatilis as compared to S. japonicus, it is not identical to the type strain, suggesting population structure within var. versatilis. We recommend that the taxonomic status of S. japonicus var. versatilis is raised, with it being treated as a separate species, Schizosaccharomyces versatilis., (© 2023 The Authors. Yeast published by John Wiley & Sons Ltd.)

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs., Competing Interests: Declaration of interests J.D.B. is a Founder and Director of CDI Labs, Inc.; a Founder of and consultant to Neochromosome, Inc.; a Founder, SAB member of, and consultant to ReOpen Diagnostics, LLC; and serves or served on the Scientific Advisory Board of the following: Logomix, Inc., Modern Meadow, Inc., Rome Therapeutics, Inc., Sample6, Inc., Sangamo, Inc., Tessera Therapeutics, Inc., and the Wyss Institute., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 ( CPAR2_104610 ) and ARR3 ( CPAR2_601050 ), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3 , which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3 , 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700 , coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene ( RTA3 or ARR3 ), and they have occurred multiple times in multiple isolates. The copy number of RTA3 , a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Cell plating, the spreading out of a liquid suspension of cells on a surface followed by colony growth, is a common laboratory procedure in microbiology. Despite this, the exact impact of its parameters on colony growth has not been extensively studied. A common protocol involves the shaking of glass beads within a Petri dish containing solid growth media. We investigated the effects of multiple parameters in this protocol: the number of beads, the shape of movement, and the number of movements. Standard suspensions of Escherichia coli were spread while varying these parameters to assess their impact on colony growth. Results were assessed by a variety of metrics: the number of colonies, the mean distance between closest colonies, and the variability and uniformity of their spatial distribution. Finally, we devised a mathematical model of shifting billiard to explain the heterogeneities in the observed spatial patterns. Exploring the parameters that affect the most fundamental techniques in microbiology allows us to better understand their function, giving us the ability to precisely control their outputs for our exact needs.

Saccharomyces cerevisiae exhibits gene expression homeostasis, which is defined as the buffering of transcription levels against changes in DNA copy number during the S phase of the cell cycle. It has been suggested that S. cerevisiae employs an active mechanism to maintain gene expression homeostasis through Rtt109-Asf1-dependent acetylation of histone H3 on lysine 56 (H3K56). Here, we show that gene expression homeostasis can be achieved independently of H3K56 acetylation by Tos4 (Target of Swi6-4). Using Nanostring technology, we establish that Tos4-dependent gene expression homeostasis depends on its forkhead-associated (FHA) domain, which is a phosphopeptide recognition domain required to bind histone deacetylases (HDACs). We demonstrate that the mechanism of Tos4-dependent gene expression homeostasis requires its interaction with the Rpd3L HDAC complex. However, this is independent of Rpd3's well-established roles in both histone deacetylation and controlling the DNA replication timing program, as established by deep sequencing of Fluorescence-Activated Cell Sorted (FACS) S and G2 phase populations. Overall, our data reveals that Tos4 mediates gene expression homeostasis through its FHA domain-dependent interaction with the Rpd3L complex, which is independent of H3K56ac., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

A eukaryotic chromosome relies on the function of multiple spatially distributed DNA replication origins for its stable inheritance. The spatial location of an origin is determined by the chromosomal position of an MCM complex, the inactive form of the DNA replicative helicase that is assembled onto DNA in G1-phase (also known as origin licensing). While the biochemistry of origin licensing is understood, the mechanisms that promote an adequate spatial distribution of MCM complexes across chromosomes are not. We have elucidated a role for the Sir2 histone deacetylase in establishing the normal distribution of MCM complexes across Saccharomyces cerevisiae chromosomes. In the absence of Sir2, MCM complexes accumulated within both early-replicating euchromatin and telomeric heterochromatin, and replication activity within these regions was enhanced. Concomitantly, the duplication of several regions of late-replicating euchromatin were delayed. Thus, Sir2-mediated attenuation of origin licensing within both euchromatin and telomeric heterochromatin established the normal spatial distribution of origins across yeast chromosomes important for normal genome duplication., Competing Interests: The authors declare no competing interest.

Biological systems are made up of components that change their actions (and interactions) over time and coordinate with other components nearby. Together with a large state space, the complexity of this behaviour can make it difficult to create concise mathematical models that can be easily extended or modified. This paper introduces the Beacon Calculus, a process algebra designed to simplify the task of modelling interacting biological components. Its breadth is demonstrated by creating models of DNA replication dynamics, the gene expression dynamics in response to DNA methylation damage, and a multisite phosphorylation switch. The flexibility of these models is shown by adapting the DNA replication model to further include two topics of interest from the literature: cooperative origin firing and replication fork barriers. The Beacon Calculus is supported with the open-source simulator bcs (https://github.com/MBoemo/bcs.git) to allow users to develop and simulate their own models., Competing Interests: The authors have declared that no competing interests exist.

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

Recent studies have accomplished the extraordinary feat of measuring the exact status of DNA replication in individual cells. We outline how these studies have revealed surprising uniformity in how cells replicate their DNA, and consider the implications of this remarkable technological advance.

Despite the conserved essential function of centromeres, centromeric DNA itself is not conserved. The histone-H3 variant, CENP-A, is the epigenetic mark that specifies centromere identity. Paradoxically, CENP-A normally assembles on particular sequences at specific genomic locations. To gain insight into the specification of complex centromeres, here we take an evolutionary approach, fully assembling genomes and centromeres of related fission yeasts. Centromere domain organization, but not sequence, is conserved between Schizosaccharomyces pombe, S. octosporus and S. cryophilus with a central CENP-A Cnp1 domain flanked by heterochromatic outer-repeat regions. Conserved syntenic clusters of tRNA genes and 5S rRNA genes occur across the centromeres of S. octosporus and S. cryophilus, suggesting conserved function. Interestingly, nonhomologous centromere central-core sequences from S. octosporus and S. cryophilus are recognized in S. pombe, resulting in cross-species establishment of CENP-A Cnp1 chromatin and functional kinetochores. Therefore, despite the lack of sequence conservation, Schizosaccharomyces centromere DNA possesses intrinsic conserved properties that promote assembly of CENP-A chromatin.

Replication of eukaryotic genomes is highly stochastic, making it difficult to determine the replication dynamics of individual molecules with existing methods. We report a sequencing method for the measurement of replication fork movement on single molecules by detecting nucleotide analog signal currents on extremely long nanopore traces (D-NAscent). Using this method, we detect 5-bromodeoxyuridine (BrdU) incorporated by Saccharomyces cerevisiae to reveal, at a genomic scale and on single molecules, the DNA sequences replicated during a pulse-labeling period. Under conditions of limiting BrdU concentration, D-NAscent detects the differences in BrdU incorporation frequency across individual molecules to reveal the location of active replication origins, fork direction, termination sites, and fork pausing/stalling events. We used sequencing reads of 20-160 kilobases to generate a whole-genome single-molecule map of DNA replication dynamics and discover a class of low-frequency stochastic origins in budding yeast. The D-NAscent software is available at https://github.com/MBoemo/DNAscent.git .

DNA replication is a stochastic process with replication forks emanating from multiple replication origins. The origins must be licenced in G1, and the replisome activated at licenced origins in order to generate bi-directional replication forks in S-phase. Differential firing times lead to origin interference, where a replication fork from an origin can replicate through and inactivate neighbouring origins (origin obscuring). We developed a Bayesian algorithm to characterize origin firing statistics from Okazaki fragment (OF) sequencing data. Our algorithm infers the distributions of firing times and the licencing probabilities for three consecutive origins. We demonstrate that our algorithm can distinguish partial origin licencing and origin obscuring in OF sequencing data from Saccharomyces cerevisiae and human cell types. We used our method to analyse the decreased origin efficiency under loss of Rat1 activity in S. cerevisiae, demonstrating that both reduced licencing and increased obscuring contribute. Moreover, we show that robust analysis is possible using only local data (across three neighbouring origins), and analysis of the whole chromosome is not required. Our algorithm utilizes an approximate likelihood and a reversible jump sampling technique, a methodology that can be extended to analysis of other mechanistic processes measurable through Next Generation Sequencing data., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

3D genome organization is strongly predictive of DNA replication timing in mammalian cells. This work tested the extent to which loop-based genome architecture acts as a regulatory unit of replication timing by using an auxin-inducible system for acute cohesin ablation. Cohesin ablation in a population of cells in asynchronous culture was shown not to disrupt patterns of replication timing as assayed by replication sequencing (RepliSeq) or BrdU-focus microscopy. Furthermore, cohesin ablation prior to S phase entry in synchronized cells was similarly shown to not impact replication timing patterns. These results suggest that cohesin-mediated genome architecture is not required for the execution of replication timing patterns in S phase, nor for the establishment of replication timing domains in G1., Competing Interests: The authors declare no conflict of interest.

Genomes are replicated in a reproducible temporal pattern. Current methods for assaying allele replication timing are time consuming and/or expensive. These include high-throughput sequencing which can be used to measure DNA copy number as a proxy for allele replication timing. Here, we use droplet digital PCR to study DNA replication timing at multiple loci in budding yeast and human cells. We establish that the method has temporal and spatial resolutions comparable to the high-throughput sequencing approaches, while being faster than alternative locus-specific methods. Furthermore, the approach is capable of allele discrimination. We apply this method to determine relative replication timing across timing transition zones in cultured human cells. Finally, multiple samples can be analysed in parallel, allowing us to rapidly screen kinetochore mutants for perturbation to centromere replication timing. Therefore, this approach is well suited to the study of locus-specific replication and the screening of cis- and trans-acting mutants to identify mechanisms that regulate local genome replication timing.

The common ancestry of archaea and eukaryotes is evident in their genome architecture. All eukaryotic and several archaeal genomes consist of multiple chromosomes, each replicated from multiple origins. Three scenarios have been proposed for the evolution of this genome architecture: 1) mutational diversification of a multi-copy chromosome; 2) capture of a new chromosome by horizontal transfer; 3) acquisition of new origins and splitting into two replication-competent chromosomes. We report an example of the third scenario: the multi-origin chromosome of the archaeon Haloferax volcanii has split into two elements via homologous recombination. The newly generated elements are bona fide chromosomes, because each bears "chromosomal" replication origins, rRNA loci, and essential genes. The new chromosomes were stable during routine growth but additional genetic manipulation, which involves selective bottlenecks, provoked further rearrangements. To the best of our knowledge, rearrangement of a naturally evolved prokaryotic genome to generate two new chromosomes has not been described previously.

The Rif1 protein negatively regulates telomeric TG repeat length in the budding yeast Saccharomyces cerevisiae, but how it prevents telomere over-extension is unknown. Rif1 was recently shown to control DNA replication by acting as a Protein Phosphatase 1 (PP1)-targeting subunit. Therefore, we investigated whether Rif1 controls telomere length by targeting PP1 activity. We find that a Rif1 mutant defective for PP1 interaction causes a long-telomere phenotype, similar to that of rif1Δ cells. Tethering PP1 at a specific telomere partially substitutes for Rif1 in limiting TG repeat length, confirming the importance of PP1 in telomere length control. Ablating Rif1-PP1 interaction is known to cause precocious activation of telomere-proximal replication origins and aberrantly early telomere replication. However, we find that Rif1 still limits telomere length even if late replication is forced through deletion of nearby replication origins, indicating that Rif1 can control telomere length independent of replication timing. Moreover we find that, even at a de novo telomere created after DNA synthesis during a mitotic block, Rif1-PP1 interaction is required to suppress telomere lengthening and prevent inappropriate recruitment of Tel1 kinase. Overall, our results show that Rif1 controls telomere length by recruiting PP1 to directly suppress telomerase-mediated TG repeat lengthening.

Background: Understanding DNA replication initiation is essential to understand the mis-regulation of replication seen in cancer and other human disorders. DNA replication initiates from DNA replication origins. In eukaryotes, replication is dependent on cell cycle kinases which function during S phase. Dbf4-dependent kinase (DDK) and cyclin-dependent kinase (CDK) act to phosphorylate the DNA helicase (composed of mini chromosome maintenance proteins: Mcm2-7) and firing factors to activate replication origins. It has recently been found that Rif1 can oppose DDK phosphorylation. Rif1 can recruit protein phosphatase 1 (PP1) to dephosphorylate MCM and restricts origin firing. In this study, we investigate a potential role for another phosphatase, protein phosphatase 2A (PP2A), in regulating DNA replication initiation. The PP2A regulatory subunit Rts1 was previously identified in a large-scale genomic screen to have a genetic interaction with ORC2 (a DNA replication licensing factor). Deletion of RTS1 synthetically rescued the temperature-sensitive (ts-) phenotype of ORC2 mutants. Methods: We deleted RTS1 in multiple ts-replication factor Saccharomyces cerevisiae strains, including ORC2 .  Dilution series assays were carried out to compare qualitatively the growth of double mutant ∆rts1 ts-replication factor strains relative to the respective single mutant strains.   Results: No synthetic rescue of temperature-sensitivity was observed. Instead we found an additive phenotype, indicating gene products function in separate biological processes. These findings are in agreement with a recent genomic screen which found that RTS1 deletion in several ts-replication factor strains led to increased temperature-sensitivity. Conclusions: We find no evidence that Rts1 is involved in the dephosphorylation of DNA replication initiation factors., Competing Interests: No competing interests were disclosed.

Eukaryotic genomes are replicated in a reproducible temporal order; however, the physiological significance is poorly understood. We compared replication timing in divergent yeast species and identified genomic features with conserved replication times. Histone genes were among the earliest replicating loci in all species. We specifically delayed the replication of HTA1 - HTB1 and discovered that this halved the expression of these histone genes. Finally, we showed that histone and cell cycle genes in general are exempt from Rtt109-dependent dosage compensation, suggesting the existence of pathways excluding specific loci from dosage compensation mechanisms. Thus, we have uncovered one of the first physiological requirements for regulated replication time and demonstrated a direct link between replication timing and gene expression., (© 2017 Müller and Nieduszynski.)

Here, we report the successful design, construction, and characterization of a 770-kilobase synthetic yeast chromosome II (synII). Our study incorporates characterization at multiple levels-including phenomics, transcriptomics, proteomics, chromosome segregation, and replication analysis-to provide a thorough and comprehensive analysis of a synthetic chromosome. Our Trans-Omics analyses reveal a modest but potentially relevant pervasive up-regulation of translational machinery observed in synII, mainly caused by the deletion of 13 transfer RNAs. By both complementation assays and SCRaMbLE (synthetic chromosome rearrangement and modification by loxP -mediated evolution), we targeted and debugged the origin of a growth defect at 37°C in glycerol medium, which is related to misregulation of the high-osmolarity glycerol response. Despite the subtle differences, the synII strain shows highly consistent biological processes comparable to the native strain., (Copyright © 2017, American Association for the Advancement of Science.)

Three eukaryotic DNA polymerases are essential for genome replication. Polymerase (Pol) α-primase initiates each synthesis event and is rapidly replaced by processive DNA polymerases: Polɛ replicates the leading strand, whereas Polδ performs lagging-strand synthesis. However, it is not known whether this division of labor is maintained across the whole genome or how uniform it is within single replicons. Using Schizosaccharomyces pombe, we have developed a polymerase usage sequencing (Pu-seq) strategy to map polymerase usage genome wide. Pu-seq provides direct replication-origin location and efficiency data and indirect estimates of replication timing. We confirm that the division of labor is broadly maintained across an entire genome. However, our data suggest a subtle variability in the usage of the two polymerases within individual replicons. We propose that this results from occasional leading-strand initiation by Polδ followed by exchange for Polɛ.

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

DNA replication initiates at defined sites called origins, which serve as binding sites for initiator proteins that recruit the replicative machinery. Origins differ in number and structure across the three domains of life and their properties determine the dynamics of chromosome replication. Bacteria and some archaea replicate from single origins, whereas most archaea and all eukaryotes replicate using multiple origins. Initiation mechanisms that rely on homologous recombination operate in some viruses. Here we show that such mechanisms also operate in archaea. We use deep sequencing to study replication in Haloferax volcanii and identify four chromosomal origins of differing activity. Deletion of individual origins results in perturbed replication dynamics and reduced growth. However, a strain lacking all origins has no apparent defects and grows significantly faster than wild type. Origin-less cells initiate replication at dispersed sites rather than at discrete origins and have an absolute requirement for the recombinase RadA, unlike strains lacking individual origins. Our results demonstrate that homologous recombination alone can efficiently initiate the replication of an entire cellular genome. This raises the question of what purpose replication origins serve and why they have evolved.

Eukaryotic genome replication is stochastic, and each cell uses a different cohort of replication origins. We demonstrate that interpreting high-resolution Saccharomyces cerevisiae genome replication data with a mathematical model allows quantification of the stochastic nature of genome replication, including the efficiency of each origin and the distribution of termination events. Single-cell measurements support the inferred values for stochastic origin activation time. A strain, in which three origins were inactivated, confirmed that the distribution of termination events is primarily dictated by the stochastic activation time of origins. Cell-to-cell variability in origin activity ensures that termination events are widely distributed across virtually the whole genome. We propose that the heterogeneity in origin usage contributes to genome stability by limiting potentially deleterious events from accumulating at particular loci., (Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.)

During S phase, the entire genome must be precisely duplicated, with no sections of DNA left unreplicated. Here, we develop a simple mathematical model to describe the probability of replication failing due to the irreversible stalling of replication forks. We show that the probability of complete genome replication is maximized if replication origins are evenly spaced, the largest inter-origin distances are minimized, and the end-most origins are positioned close to chromosome ends. We show that origin positions in the yeast Saccharomyces cerevisiae genome conform to all three predictions thereby maximizing the probability of complete replication if replication forks stall. Origin positions in four other yeasts-Kluyveromyces lactis, Lachancea kluyveri, Lachancea waltii and Schizosaccharomyces pombe-also conform to these predictions. Equating failure rates at chromosome ends with those in chromosome interiors gives a mean per nucleotide fork stall rate of ∼5 × 10(-8), which is consistent with experimental estimates. Using this value in our theoretical predictions gives replication failure rates that are consistent with data from replication origin knockout experiments. Our theory also predicts that significantly larger genomes, such as those of mammals, will experience a much greater probability of replication failure genome-wide, and therefore will likely require additional compensatory mechanisms.

Inside the nucleus, DNA replication is organized at discrete sites called replication factories, consisting of DNA polymerases and other replication proteins. Replication factories play important roles in coordinating replication and in responding to replication stress. However, it remains unknown how replicons are organized for processing at each replication factory. Here we address this question using budding yeast. We analyze how individual replicons dynamically organized a replication factory using live-cell imaging and investigate how replication factories were structured using super-resolution microscopy. Surprisingly, we show that the grouping of replicons within factories is highly variable from cell to cell. Once associated, however, replicons stay together relatively stably to maintain replication factories. We derive a coherent genome-wide mathematical model showing how neighboring replicons became associated stochastically to form replication factories, which was validated by independent microscopy-based analyses. This study not only reveals the fundamental principles promoting replication factory organization in budding yeast, but also provides insight into general mechanisms by which chromosomes organize sub-nuclear structures.

Chromosome duplication normally initiates through the assembly of replication fork complexes at defined origins. DNA synthesis by any one fork is thought to cease when it meets another travelling in the opposite direction, at which stage the replication machinery may simply dissociate before the nascent strands are finally ligated. But what actually happens is not clear. Here we present evidence consistent with the idea that every fork collision has the potential to threaten genomic integrity. In Escherichia coli this threat is kept at bay by RecG DNA translocase and by single-strand DNA exonucleases. Without RecG, replication initiates where forks meet through a replisome assembly mechanism normally associated with fork repair, replication restart and recombination, establishing new forks with the potential to sustain cell growth and division without an active origin. This potential is realized when roadblocks to fork progression are reduced or eliminated. It relies on the chromosome being circular, reinforcing the idea that replication initiation is triggered repeatedly by fork collision. The results reported raise the question of whether replication fork collisions have pathogenic potential for organisms that exploit several origins to replicate each chromosome.

Centromeres play several important roles in ensuring proper chromosome segregation. Not only do they promote kinetochore assembly for microtubule attachment, but they also support robust sister chromatid cohesion at pericentromeres and facilitate replication of centromeric DNA early in S phase. However, it is still elusive how centromeres orchestrate all these functions at the same site. Here, we show that the budding yeast Dbf4-dependent kinase (DDK) accumulates at kinetochores in telophase, facilitated by the Ctf19 kinetochore complex. This promptly recruits Sld3-Sld7 replication initiator proteins to pericentromeric replication origins so that they initiate replication early in S phase. Furthermore, DDK at kinetochores independently recruits the Scc2-Scc4 cohesin loader to centromeres in G1 phase. This enhances cohesin loading and facilitates robust pericentromeric cohesion in S phase. Thus, we have found the central mechanism by which kinetochores orchestrate early S phase DNA replication and robust sister chromatid cohesion at microtubule attachment sites., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Background: Comparative genomics is a formidable tool to identify functional elements throughout a genome. In the past ten years, studies in the budding yeast Saccharomyces cerevisiae and a set of closely related species have been instrumental in showing the benefit of analyzing patterns of sequence conservation. Increasing the number of closely related genome sequences makes the comparative genomics approach more powerful and accurate., Results: Here, we report the genome sequence and analysis of Saccharomyces arboricolus, a yeast species recently isolated in China, that is closely related to S. cerevisiae. We obtained high quality de novo sequence and assemblies using a combination of next generation sequencing technologies, established the phylogenetic position of this species and considered its phenotypic profile under multiple environmental conditions in the light of its gene content and phylogeny., Conclusions: We suggest that the genome of S. arboricolus will be useful in future comparative genomics analysis of the Saccharomyces sensu stricto yeasts.

Eukaryotic DNA replication origins are selected in G1-phase when the origin recognition complex (ORC) binds chromosomal positions and triggers molecular events culminating in the initiation of DNA replication (a.k.a. origin firing) during S-phase. Each chromosome uses multiple origins for its duplication, and each origin fires at a characteristic time during S-phase, creating a cell-type specific genome replication pattern relevant to differentiation and genome stability. It is unclear whether ORC-origin interactions are relevant to origin activation time. We applied a novel genome-wide strategy to classify origins in the model eukaryote Saccharomyces cerevisiae based on the types of molecular interactions used for ORC-origin binding. Specifically, origins were classified as DNA-dependent when the strength of ORC-origin binding in vivo could be explained by the affinity of ORC for origin DNA in vitro, and, conversely, as 'chromatin-dependent' when the ORC-DNA interaction in vitro was insufficient to explain the strength of ORC-origin binding in vivo. These two origin classes differed in terms of nucleosome architecture and dependence on origin-flanking sequences in plasmid replication assays, consistent with local features of chromatin promoting ORC binding at 'chromatin-dependent' origins. Finally, the 'chromatin-dependent' class was enriched for origins that fire early in S-phase, while the DNA-dependent class was enriched for later firing origins. Conversely, the latest firing origins showed a positive association with the ORC-origin DNA paradigm for normal levels of ORC binding, whereas the earliest firing origins did not. These data reveal a novel association between ORC-origin binding mechanisms and the regulation of origin activation time., Competing Interests: The authors have declared that no competing interests exist.

DNA replication initiates from defined locations called replication origins; some origins are highly active, whereas others are dormant and rarely used. Origins also differ in their activation time, resulting in particular genomic regions replicating at characteristic times and in a defined temporal order. Here we report the comparison of genome replication in four budding yeast species: Saccharomyces cerevisiae, S. paradoxus, S. arboricolus, and S. bayanus. First, we find that the locations of active origins are predominantly conserved between species, whereas dormant origins are poorly conserved. Second, we generated genome-wide replication profiles for each of these species and discovered that the temporal order of genome replication is highly conserved. Therefore, active origins are not only conserved in location, but also in activation time. Only a minority of these conserved origins show differences in activation time between these species. To gain insight as to the mechanisms by which origin activation time is regulated we generated replication profiles for a S. cerevisiae/S. bayanus hybrid strain and find that there are both local and global regulators of origin function.

Eukaryotic DNA replication is initiated from multiple sites on the chromosome, but little is known about the global and local regulation of replication. We present a mathematical model for the spatial dynamics of DNA replication, which offers insight into the kinetics of replication in different types of organisms. Most biological experiments involve average quantities over large cell populations (typically >10(7) cells) and therefore can mask the cell-to-cell variability present in the system. Although the model is formulated in terms of a population of cells, using mathematical analysis we show that one can obtain signatures of stochasticity in individual cells from averaged quantities. This work generalizes the result by Retkute et al. [Phys. Rev. Lett. 107, 068103 (2011)] to a broader set of parameter regimes.

Cell-cycle progression is governed by a series of essential regulatory proteins. Two major regulators are cell-division cycle protein 20 (CDC20) and its homologue, CDC20 homologue 1 (CDH1), which activate the anaphase-promoting complex/cyclosome (APC/C) in mitosis, and facilitate degradation of mitotic APC/C substrates. The malaria parasite, Plasmodium, is a haploid organism which, during its life-cycle undergoes two stages of mitosis; one associated with asexual multiplication and the other with male gametogenesis. Cell-cycle regulation and DNA replication in Plasmodium was recently shown to be dependent on the activity of a number of protein kinases. However, the function of cell division cycle proteins that are also involved in this process, such as CDC20 and CDH1 is totally unknown. Here we examine the role of a putative CDC20/CDH1 in the rodent malaria Plasmodium berghei (Pb) using reverse genetics. Phylogenetic analysis identified a single putative Plasmodium CDC20/CDH1 homologue (termed CDC20 for simplicity) suggesting that Plasmodium APC/C has only one regulator. In our genetic approach to delete the endogenous cdc20 gene of P. berghei, we demonstrate that PbCDC20 plays a vital role in male gametogenesis, but is not essential for mitosis in the asexual blood stage. Furthermore, qRT-PCR analysis in parasite lines with deletions of two kinase genes involved in male sexual development (map2 and cdpk4), showed a significant increase in cdc20 transcription in activated gametocytes. DNA replication and ultra structural analyses of cdc20 and map2 mutants showed similar blockage of nuclear division at the nuclear spindle/kinetochore stage. CDC20 was phosphorylated in asexual and sexual stages, but the level of modification was higher in activated gametocytes and ookinetes. Changes in global protein phosphorylation patterns in the Δcdc20 mutant parasites were largely different from those observed in the Δmap2 mutant. This suggests that CDC20 and MAP2 are both likely to play independent but vital roles in male gametogenesis.

OriDB (http://www.oridb.org/) is a database containing collated genome-wide mapping studies of confirmed and predicted replication origin sites. The original database collated and curated Saccharomyces cerevisiae origin mapping studies. Here, we report that the OriDB database and web site have been revamped to improve user accessibility to curated data sets, to greatly increase the number of curated origin mapping studies, and to include the collation of replication origin sites in the fission yeast Schizosaccharomyces pombe. The revised database structure underlies these improvements and will facilitate further expansion in the future. The updated OriDB for S. cerevisiae is available at http://cerevisiae.oridb.org/ and for S. pombe at http://pombe.oridb.org/.

Background: Natural variation offers a powerful approach for assigning function to DNA sequence-a pressing challenge in the age of high throughput sequencing technologies., Scope of Review: Here we review comparative genomic approaches that are bridging the sequence-function and genotype-phenotype gaps. Reverse genomic approaches aim to analyse sequence to assign function, whereas forward genomic approaches start from a phenotype and aim to identify the underlying genotype responsible., Major Conclusions: Comparative genomic approaches, pioneered in budding yeasts, have resulted in dramatic improvements in our understanding of the function of both genes and regulatory sequences. Analogous studies in other systems, including humans, demonstrate the ubiquity of comparative genomic approaches. Recently, forward genomic approaches, exploiting natural variation within yeast populations, have started to offer powerful insights into how genotype influences phenotype and even the ability to predict phenotypes., General Significance: Comparative genomic experiments are defining the fundamental rules that govern complex traits in natural populations from yeast to humans. This article is part of a Special Issue entitled Systems Biology of Microorganisms., (Copyright © 2011 Elsevier B.V. All rights reserved.)

We present a mathematical model for the spatial dynamics of DNA replication. Using this model we determine the probability distribution for the time at which each chromosomal position is replicated. From this we show, contrary to previous reports, that mean replication time curves cannot be used to directly determine origin parameters. We demonstrate that the stochastic nature of replication dynamics leaves a clear signature in experimentally measured population average data, and we show that the width of the activation time probability distribution can be inferred from this data. Our results compare favorably with experimental measurements in Saccharomyces cerevisae.

The fission yeast clade--comprising Schizosaccharomyces pombe, S. octosporus, S. cryophilus, and S. japonicus--occupies the basal branch of Ascomycete fungi and is an important model of eukaryote biology. A comparative annotation of these genomes identified a near extinction of transposons and the associated innovation of transposon-free centromeres. Expression analysis established that meiotic genes are subject to antisense transcription during vegetative growth, which suggests a mechanism for their tight regulation. In addition, trans-acting regulators control new genes within the context of expanded functional modules for meiosis and stress response. Differences in gene content and regulation also explain why, unlike the budding yeast of Saccharomycotina, fission yeasts cannot use ethanol as a primary carbon source. These analyses elucidate the genome structure and gene regulation of fission yeast and provide tools for investigation across the Schizosaccharomyces clade.

All chromosomes must be completely replicated prior to cell division, a requirement that demands the activation of a sufficient number of appropriately distributed DNA replication origins. Here we investigate how the activity of multiple origins on each chromosome is coordinated to ensure successful replication. We present a stochastic model for whole chromosome replication where the dynamics are based upon the parameters of individual origins. Using this model we demonstrate that mean replication time at any given chromosome position is determined collectively by the parameters of all origins. Combining parameter estimation with extensive simulations we show that there is a range of model parameters consistent with mean replication data, emphasising the need for caution in interpreting such data. In contrast, the replicated-fraction at time points through S phase contains more information than mean replication time data and allowed us to use our model to uniquely estimate many origin parameters. These estimated parameters enable us to make a number of predictions that showed agreement with independent experimental data, confirming that our model has predictive power. In summary, we demonstrate that a stochastic model can recapitulate experimental observations, including those that might be interpreted as deterministic such as ordered origin activation times.

The origin recognition complex (ORC) marks chromosomal sites as replication origins and is essential for replication initiation. In yeast, ORC also binds to DNA elements called silencers, where its primary function is to recruit silent information regulator (SIR) proteins to establish transcriptional silencing. Indeed, silencers function poorly as chromosomal origins. Several genetic, molecular, and biochemical studies of HMR-E have led to a model proposing that when ORC becomes limiting in the cell (such as in the orc2-1 mutant) only sites that bind ORC tightly (such as HMR-E) remain fully occupied by ORC, while lower affinity sites, including many origins, lose ORC occupancy. Since HMR-E possessed a unique non-replication function, we reasoned that other tight sites might reveal novel functions for ORC on chromosomes. Therefore, we comprehensively determined ORC "affinity" genome-wide by performing an ORC ChIP-on-chip in ORC2 and orc2-1 strains. Here we describe a novel group of orc2-1-resistant ORC-interacting chromosomal sites (ORF-ORC sites) that did not function as replication origins or silencers. Instead, ORF-ORC sites were comprised of protein-coding regions of highly transcribed metabolic genes. In contrast to the ORC-silencer paradigm, transcriptional activation promoted ORC association with these genes. Remarkably, ORF-ORC genes were enriched in proximity to origins of replication and, in several instances, were transcriptionally regulated by these origins. Taken together, these results suggest a surprising connection among ORC, replication origins, and cellular metabolism., Competing Interests: The authors have declared that no competing interests exist.

Effective experimental techniques are available to identify replication origin regions in eukaryotic cells. Genome-wide identification of the precise sequence elements that direct origin activity is however still not straightforward, even in the yeast Saccharomyces cerevisiae which has the best characterised eukaryotic replication origins. The availability of genome sequences for a series of closely related (sensu stricto) budding yeasts has allowed us to take a 'comparative genomics' approach to this problem. Since they represent functional protein-binding sites, origin sequences are conserved better than the surrounding intergenic sequence within the genomes of closely related yeasts. We describe here how phylogenetic comparison data can be used to identify candidate replication origin sequences in the S. cerevisiae genome, and how large numbers of such candidate sites can simultaneously be assayed for ability to initiate replication. Similar approaches could potentially be used to identify protein-binding sequence elements having other functions, as well as replication origin sites in other species.

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.

Replication of eukaryotic chromosomes initiates at multiple sites called replication origins. Replication origins are best understood in the budding yeast Saccharomyces cerevisiae, where several complementary studies have mapped their locations genome-wide. We have collated these datasets, taking account of the resolution of each study, to generate a single list of distinct origin sites. OriDB provides a web-based catalogue of these confirmed and predicted S.cerevisiae DNA replication origin sites. Each proposed or confirmed origin site appears as a record in OriDB, with each record comprising seven pages. These pages provide, in text and graphical formats, the following information: genomic location and chromosome context of the origin site; time of origin replication; DNA sequence of proposed or experimentally confirmed origin elements; free energy required to open the DNA duplex (stress-induced DNA duplex destabilization or SIDD); and phylogenetic conservation of sequence elements. In addition, OriDB encourages community submission of additional information for each origin site through a User Notes facility. Origin sites are linked to several external resources, including the Saccharomyces Genome Database (SGD) and relevant publications at PubMed. Finally, a Chromosome Viewer utility allows users to interactively generate graphical representations of DNA replication data genome-wide. OriDB is available at www.oridb.org.

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.

The mini-chromosome maintenance proteins Mcm2-7 are essential for DNA replication. They are loaded onto replication origins during G1 phase of the cell cycle to form a pre-replication complex (pre-RC) that licenses each origin for subsequent initiation. We have investigated the DNA elements that determine the dependence of yeast replication origins on Mcm2-7 activity, i.e. the sensitivity of an origin to mcm mutations. Using chimaeric constructs from mcm sensitive and mcm insensitive origins, we have identified two main elements affecting the requirement for Mcm2-7 function. First, transcription into an origin increases its dependence on Mcm2-7 function, revealing a conflict between pre-RC assembly and transcription. Second, sequence elements within the minimal origin influence its mcm sensitivity. Replication origins show similar differences in sensitivity to mutations in other pre-RC proteins (such as Origin Recognition Complex and Cdc6), but not to mutations in initiation and elongation factors, demonstrating that the mcm sensitivity of an origin is determined by its ability to establish a pre-RC. We propose that there is a hierarchy of replication origins with respect to the range of pre-RC protein concentrations under which they will function. This hierarchy is both 'hard-wired' by the minimal origin sequences and 'soft-wired' by local transcriptional context.

We have investigated whether the Ku complex is involved in regulating DNA replication in the yeast Saccharomyces cerevisiae. We find that Ku proteins control the replication time of telomeric regions; replication origins located close to telomeres or within subtelomeric repeat sequences normally initiate late, but are activated much earlier in mutants lacking Ku function. In contrast, origins distant from telomeres initiate replication at the normal time. Ku is one of the first components identified as important for replication timing, and specification of the replication time of chromosome ends by Ku is consistent with its role in maintaining telomere localization.

Background: Multiple A- and B-type cyclins have been identified in animals, but their study is complicated by varying degrees of functional redundancy. A non-essential phenotype may reflect redundancy with a known or as yet unknown gene. Complete sequencing of several animal genomes has allowed us to determine the size of the mitotic cyclin gene family and therefore to start to address this issue., Results: We analyzed the Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens genomes to identify known and novel A- and B-type cyclin genes and distinguish them from related pseudogenes. We find only a single functional A-type cyclin gene in invertebrates but two in vertebrates. In addition to the single functional cyclin A gene, the C. elegans genome contains numerous cyclin A pseudogenes. In contrast, the number and relationship of B-type cyclins varies considerably between organisms but all contain at least one cyclin B1-like gene and a cyclin B3 gene., Conclusions: There are three conserved families of mitotic cyclins in animals: A-, B3- and B-type. The precise number of genes within the A- and B-type families varies in different organisms, possibly as an adaptation to their distinct developmental strategies.