Novel features of ARS selection in budding yeast Lachancea kluyveri

Various short autonomously replicating sequences from the yeast Kluyveromyces marxianus seemingly without canonical consensus

Eukaryotic autonomously replicating sequences (ARSs) are composed of three domains, A, B, and C. Domain A is comprised of an ARS consensus sequence (ACS), while the B domain has the DNA unwinding element and the C domain is important for DNA-protein interactions. In Saccharomyces cerevisiae and Kluyveromyces lactis ARS101, the ACS is commonly composed of 11 bp, 5ˊ-(A/T)AAA(C/T)ATAAA(A/T)-3ˊ. This core sequence is essential for S. cerevisiae and K. lactis ARS activity. In this study, we identified ARS-containing sequences from genomic libraries of the yeast Kluyveromyces marxianus DMKU3-1042 and validated their replication activities. The identified K. marxianus DMKU3-1042 ARSs (KmARSs) have very effective replication ability but their sequences are divergent and share no common consensus. We have carried out point mutations, deletions, and base pairs substitutions within the sequences of some of the KmARSs to identify the sequence(s) that influence the replication activity. Consensus sequences same as the 11 bp ACS of S. cerevisiae and K. lactis were not found in all minimum functional KmARSs reported here except KmARS7. Moreover, partial sequences from different KmARSs are interchangeable among each other to retain the ARS activity. We have also specifically identified the essential nucleotides, which are indispensable for replication, within some of the KmARSs. Our deletions analysis revealed that only 21 bp in KmARS18 could retain the ARS activity. The identified KmARSs in this study are unique compared to other yeasts’ ARSs, do not share common ACS, and are interchangeable.

•

Identification of minimal autonomously replicating sequences (ARSs) from the yeast Kluyveromyces marxianus DMKU3-1042.

•

The identities of the isolated ARSs are divergent and have no common consensus with the ARSs of other yeasts.

•

A short ARS sequence of twenty-one nucleotides functions as an effective replicator in K. marxianus DMKU3-1042.

•

Segments of ARSs from the yeast K. marxianus are interchangeable among each other.

•

Functional ARSs are found in both the intergenic and coding sequences of the strain DMKU3-1042.

Conservation of a DNA Replication Motif among Phylogenetically Distant Budding Yeast Species

Eukaryotic DNA replication begins at genomic loci termed origins, which are bound by the origin recognition complex (ORC). Although ORC is conserved across species, the sequence composition of origins is more varied. In the budding yeast Saccharomyces cerevisiae, the ORC-binding motif consists of an A/T-rich 17 bp “extended ACS” sequence adjacent to a B1 element composed of two 3-bp motifs. Similar sequences occur at origins in closely related species, but it is not clear when this type of replication origin arose and whether it predated a whole-genome duplication that occurred around 100 Ma in the budding yeast lineage. To address these questions, we identified the ORC-binding sequences in the nonduplicated species Torulaspora delbrueckii. We used chromatin immunoprecipitation followed by sequencing and identified 190 ORC-binding sites distributed across the eight T. delbrueckii chromosomes. Using these sites, we identified an ORC-binding motif that is nearly identical to the known motif in S. cerevisiae. We also found that the T. delbrueckii ORC-binding sites function as origins in T. delbrueckii when cloned onto a plasmid and that the motif is required for plasmid replication. Finally, we compared an S. cerevisiae origin with two T. delbrueckii ORC-binding sites and found that they conferred similar stabilities to a plasmid. These results reveal that the ORC-binding motif arose prior to the whole-genome duplication and has been maintained for over 100 Myr.

The evolution of the temporal program of genome replication

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

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

DNA replication timing influences gene expression level

Eukaryotic genomes are replicated in a reproducible temporal order whose physiological significance is poorly understood. Müller and Nieduszynski compare the temporal order of genome replication in phylogenetically diverse yeast species and identify genes for which conserved replication timing contributes to maximal expression.

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.

Context based computational analysis and characterization of ARS consensus sequences (ACS) of Saccharomyces cerevisiae genome

Genome-wide experimental studies in Saccharomyces cerevisiae reveal that autonomous replicating sequence (ARS) requires an essential consensus sequence (ACS) for replication activity. Computational studies identified thousands of ACS like patterns in the genome. However, only a few hundreds of these sites act as replicating sites and the rest are considered as dormant or evolving sites. In a bid to understand the sequence makeup of replication sites, a content and context-based analysis was performed on a set of replicating ACS sequences that binds to origin-recognition complex (ORC) denoted as ORC-ACS and non-replicating ACS sequences (nrACS), that are not bound by ORC. In this study, DNA properties such as base composition, correlation, sequence dependent thermodynamic and DNA structural profiles, and their positions have been considered for characterizing ORC-ACS and nrACS. Analysis reveals that ORC-ACS depict marked differences in nucleotide composition and context features in its vicinity compared to nrACS. Interestingly, an A-rich motif was also discovered in ORC-ACS sequences within its nucleosome-free region. Profound changes in the conformational features, such as DNA helical twist, inclination angle and stacking energy between ORC-ACS and nrACS were observed. Distribution of ACS motifs in the non-coding segments points to the locations of ORC-ACS which are found far away from the adjacent gene start position compared to nrACS thereby enabling an accessible environment for ORC-proteins. Our attempt is novel in considering the contextual view of ACS and its flanking region along with nucleosome positioning in the S. cerevisiae genome and may be useful for any computational prediction scheme.

Genome-wide replication landscape of Candida glabrata

Background

The opportunistic pathogen Candida glabrata is a member of the Saccharomycetaceae yeasts. Like its close relative Saccharomyces cerevisiae, it underwent a whole-genome duplication followed by an extensive loss of genes. Its genome contains a large number of very long tandem repeats, called megasatellites. In order to determine the whole replication program of the C. glabrata genome and its general chromosomal organization, we used deep-sequencing and chromosome conformation capture experiments.

Results

We identified 253 replication fork origins, genome wide. Centromeres, HML and HMR loci, and most histone genes are replicated early, whereas natural chromosomal breakpoints are located in late-replicating regions. In addition, 275 autonomously replicating sequences (ARS) were identified during ARS-capture experiments, and their relative fitness was determined during growth competition. Analysis of ARSs allowed us to identify a 17-bp consensus, similar to the S. cerevisiae ARS consensus sequence but slightly more constrained. Megasatellites are not in close proximity to replication origins or termini. Using chromosome conformation capture, we also show that early origins tend to cluster whereas non-subtelomeric megasatellites do not cluster in the yeast nucleus.

Conclusions

Despite a shorter cell cycle, the C. glabrata replication program shares unexpected striking similarities to S. cerevisiae, in spite of their large evolutionary distance and the presence of highly repetitive large tandem repeats in C. glabrata. No correlation could be found between the replication program and megasatellites, suggesting that their formation and propagation might not be directly caused by replication fork initiation or termination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0177-6) contains supplementary material, which is available to authorized users.

Origin replication complex binding, nucleosome depletion patterns, and a primary sequence motif can predict origins of replication in a genome with epigenetic centromeres

Origins of DNA replication are key genetic elements, yet their identification remains elusive in most organisms. In previous work, we found that centromeres contain origins of replication (ORIs) that are determined epigenetically in the pathogenic yeast Candida albicans. In this study, we used origin recognition complex (ORC) binding and nucleosome occupancy patterns in Saccharomyces cerevisiae and Kluyveromyces lactis to train a machine learning algorithm to predict the position of active arm (noncentromeric) origins in the C. albicans genome. The model identified bona fide active origins as determined by the presence of replication intermediates on nondenaturing two-dimensional (2D) gels. Importantly, these origins function at their native chromosomal loci and also as autonomously replicating sequences (ARSs) on a linear plasmid. A “mini-ARS screen” identified at least one and often two ARS regions of ≥100 bp within each bona fide origin. Furthermore, a 15-bp AC-rich consensus motif was associated with the predicted origins and conferred autonomous replicating activity to the mini-ARSs. Thus, while centromeres and the origins associated with them are epigenetic, arm origins are dependent upon critical DNA features, such as a binding site for ORC and a propensity for nucleosome exclusion.

DNA replication machinery is highly conserved, yet the definition of exactly what specifies a replication origin differs in different species. Here, we utilized computational genomics to predict origin locations in Candida albicans by combining locations of binding sites for the conserved origin replication complex, necessary for replication initiation, together with chromatin organization patterns. We identified predicted sequences that exhibited bona fide origin function and developed a linear plasmid assay to delimit the DNA fragments necessary for origin function. Additionally, we found that a short AC-rich motif, which is enriched in predicted origins, is required for origin function. Thus, we demonstrated a new machine learning paradigm for identification of potential origins from a genome with no prior information. Furthermore, this work suggests that C. albicans has two different types of origins: “hard-wired” arm origins that rely upon specific sequence motifs and “epigenetic” centromeric origins that are recruited to kinetochores in a sequence-independent manner.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

An autonomously replicating sequence for use in a wide range of budding yeasts

The initiation of DNA replication at replication origins is essential for the duplication of genomes. In yeast, the autonomously replicating sequence (ARS) property of replication origins is necessary for the stable maintenance of episomal plasmids. However, because the sequence determinants of ARS function differ among yeast species, current ARS modules are limited for use to a subset of yeasts. Here, we describe a short ARS sequence that functions in at least 10 diverse species of budding yeast. These include, but are not limited to members of the Saccharomyces, Lachancea, Kluyveromyces, and Pichia (Komagataella) genera spanning over 500 million years of evolution. In addition to its wide species range, this ARS and an optimized derivative confer improved plasmid stability relative to other currently used ARS modules.

Specification of DNA replication origins and genomic base composition in fission yeasts

In the "Replicon Theory", Jacob, Brenner and Cuzin proposed the existence of replicators and initiators as the two major actors in DNA replication. Over the years, many protein components of initiators have been shown to be conserved in different organisms during evolution. By contrast, replicator DNA sequences (often referred to as replication origins) have diverged beyond possible comparison between eukaryotic genomes. Replication origins in the fission yeast Schizosaccharomyces pombe are made up of A+T-rich sequences that do not share any consensus elements. The information encoded in these replicators is interpreted by the Orc4 subunit of the ORC (origin recognition complex), which is unique among eukaryotes in that it contains a large domain harboring nine AT-hook subdomains that target ORC to a great variety of A+T-rich sequences along the chromosomes. Recently, the genomes of other Schizosaccharomyces species have been sequenced and the regions encompassing their replication origins have been identified. DNA sequence analysis and comparison of the organization of their Orc4 proteins have revealed species-specific differences that contribute to our understanding of how the specification of replication origins has evolved during the phylogenetic divergence of fission yeasts.

DNA replication origins

The onset of genomic DNA synthesis requires precise interactions of specialized initiator proteins with DNA at sites where the replication machinery can be loaded. These sites, defined as replication origins, are found at a few unique locations in all of the prokaryotic chromosomes examined so far. However, replication origins are dispersed among tens of thousands of loci in metazoan chromosomes, thereby raising questions regarding the role of specific nucleotide sequences and chromatin environment in origin selection and the mechanisms used by initiators to recognize replication origins. Close examination of bacterial and archaeal replication origins reveals an array of DNA sequence motifs that position individual initiator protein molecules and promote initiator oligomerization on origin DNA. Conversely, the need for specific recognition sequences in eukaryotic replication origins is relaxed. In fact, the primary rule for origin selection appears to be flexibility, a feature that is modulated either by structural elements or by epigenetic mechanisms at least partly linked to the organization of the genome for gene expression.

Genetic and epigenetic determinants of DNA replication origins, position and activation

In the genome of eukaryotic cells, DNA synthesis is initiated at multiple sites called origins of DNA replication. Origins must fire only once per cell cycle and how this is achieved is now well understood. However, little is known about the mechanisms that determine when and where replication initiates in a given cell. A large body of evidence indicates that origins are not equal in terms of efficiency and timing of activation. Origin usage also changes concomitantly with the different cell differentiation programs. As DNA replication occurs in the context of chromatin, initiation could be influenced by multiple parameters, such as nucleosome positioning, histone modifications, and three-dimensional (3D) organization of the nucleus. This view is supported by recent genome-wide studies showing that DNA replication profiles are shaped by genetic and epigenetic processes that act both at the local and global levels to regulate origin function in eukaryotic cells.

The spatiotemporal program of replication in the genome of Lachancea kluyveri

We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition. This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication. The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

High-resolution mapping, characterization, and optimization of autonomously replicating sequences in yeast

DNA replication origins are necessary for the duplication of genomes. In addition, plasmid-based expression systems require DNA replication origins to maintain plasmids efficiently. The yeast autonomously replicating sequence (ARS) assay has been a valuable tool in dissecting replication origin structure and function. However, the dearth of information on origins in diverse yeasts limits the availability of efficient replication origin modules to only a handful of species and restricts our understanding of origin function and evolution. To enable rapid study of origins, we have developed a sequencing-based suite of methods for comprehensively mapping and characterizing ARSs within a yeast genome. Our approach finely maps genomic inserts capable of supporting plasmid replication and uses massively parallel deep mutational scanning to define molecular determinants of ARS function with single-nucleotide resolution. In addition to providing unprecedented detail into origin structure, our data have allowed us to design short, synthetic DNA sequences that retain maximal ARS function. These methods can be readily applied to understand and modulate ARS function in diverse systems.

Conservation of replication timing reveals global and local regulation of replication origin activity

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.

Maintaining replication origins in the face of genomic change

Origins of replication present a paradox to evolutionary biologists. As a collection, they are absolutely essential genomic features, but individually are highly redundant and nonessential. It is therefore difficult to predict to what extent and in what regard origins are conserved over evolutionary time. Here, through a comparative genomic analysis of replication origins and chromosomal replication patterns in the budding yeasts Saccharomyces cerevisiae and Lachancea waltii, we assess to what extent replication origins survived genomic change produced from 150 million years of evolution. We find that L. waltii origins exhibit a core consensus sequence and nucleosome occupancy pattern highly similar to those of S. cerevisiae origins. We further observe that the overall progression of chromosomal replication is similar between L. waltii and S. cerevisiae. Nevertheless, few origins show evidence of being conserved in location between the two species. Among the conserved origins are those surrounding centromeres and adjacent to histone genes, suggesting that proximity to an origin may be important for their regulation. We conclude that, over evolutionary time, origins maintain sequence, structure, and regulation, but are continually being created and destroyed, with the result that their locations are generally not conserved.