Budding Yeast Strains Research Articles

Reliable yeast replicative lifespan determination with a single-channel microfluidic chip.

Saccharomyces cerevisiae is a powerful model for aging research due to its short lifespan and genetic malleability. Microfluidic devices offer an attractive approach enabling rapid monitoring of hundreds of cells during their entire replicative lifespan (RLS). Yet, key operational issues such as contaminations, cell loss, and cell aggregates-dependent flow obstruction can hinder RLS experiments. We report the development of a microfluidic device configuration that effectively prevents flow blockage. We conducted comprehensive performance characterization, evaluating trapping efficiency, cell retention, budding orientation, and cell aggregate formation. The optimized device successfully supported long-term culturing and reliable RLS measurements of budding yeast strains. For accurate lifespan determination, a detailed workflow is provided that includes device fabrication, live microscopy setup, and characterization of cell age distribution. This work describes an accessible and reliable microfluidic device for yeast RLS studies, promoting further exploration in aging research.

Specificity of Membrane-Associated J-Domain Protein, Caj1, in Amphotericin B Tolerance in Budding Yeast.

Hsp70:J-domain protein (JDP) machineries play pivotal roles in maintaining cellular proteostasis and governing various aspects of fungal physiology. While Hsp70 is known for its involvement in conferring tolerance to diverse antifungal drugs, the specific contribution of JDPs remains unclear. In this study, we examined the sensitivity of cytosolic JDP deletion strains of budding yeast to amphotericin B (AmB), a polyene antifungal agent widely utilized in fungal disease treatment due to its ability to disrupt the fungal plasma membrane (PM). Deleting Caj1, a PM-associated class II JDP, heightened susceptibility to AmB, and the protection conferred by Caj1 against AmB necessitated both its N-terminal J-domain and C-terminal lipid binding domain. Moreover, Caj1 deficiency compromised PM integrity as evidenced by increased phosphate efflux and exacerbated AmB sensitivity, particularly at elevated temperatures. Notably, phytosphingosine (PHS) addition as well as overexpression of PMP3, a positive PM integrity regulator, significantly rescued AmB sensitivity of caj1Δ cells. Our results align with the notion that Caj1 associates with the PM and cooperates with Hsp70 to regulate PM proteostasis, thereby influencing PM integrity in budding yeast. Loss of Caj1 function at the PM compromises PM protein quality control, thereby rendering yeast cells more susceptible to AmB.

Spindle architecture constrains karyotype evolution

The eukaryotic cell division machinery must rapidly and reproducibly duplicate and partition the cell’s chromosomes in a carefully coordinated process. However, chromosome numbers vary dramatically between genomes, even on short evolutionary timescales. We sought to understand how the mitotic machinery senses and responds to karyotypic changes by using a series of budding yeast strains in which the native chromosomes have been successively fused. Using a combination of cell biological profiling, genetic engineering and experimental evolution, we show that chromosome fusions are well tolerated up until a critical point. Cells with fewer than five centromeres lack the necessary number of kinetochore-microtubule attachments needed to counter outward forces in the metaphase spindle, triggering the spindle assembly checkpoint and prolonging metaphase. Our findings demonstrate that spindle architecture is a constraining factor for karyotype evolution.

Plasticity of the mitotic spindle in response to karyotype variation

Karyotypes, composed of chromosomes, must be accurately partitioned by the mitotic spindle for optimal cell health. However, it is unknown how underlying characteristics of karyotypes, such as chromosome number and size, govern the scaling of the mitotic spindle to ensure accurate chromosome segregation and cell proliferation. We utilize budding yeast strains engineered with fewer chromosomes, including just two "mega chromosomes," to study how spindle size and function are responsive to, and scaled by, karyotype. We determined that deletion and overexpression of spindle-related genes are detrimental to the growth of strains with two chromosomes, suggesting that mega chromosomes exert altered demands on the spindle. Using confocal microscopy, we demonstrate that cells with fewer but longer chromosomes have smaller spindle pole bodies, fewer microtubules, and longer spindles. Moreover, using electron tomography and confocal imaging, we observe elongated, bent anaphase spindles with fewer core microtubules in strains with mega chromosomes. Cells harboring mega chromosomes grow more slowly, are delayed in mitosis, and a subset struggle to complete chromosome segregation. We propose that the karyotype of the cell dictates the microtubule number, type, spindle pole body size, and spindle length, subsequently influencing the dynamics of mitosis, such as the rate of spindle elongation and the velocity of pole separation. Taken together, our results suggest that mitotic spindles are highly plastic ultrastructures that can accommodate and adjust to a variety of karyotypes, even within a species.

Spindle architecture constrains karyotype in budding yeast.

The eukaryotic cell division machinery must rapidly and reproducibly duplicate and partition the cell's chromosomes in a carefully coordinated process. However, chromosome number varies dramatically between genomes, even on short evolutionary timescales. We sought to understand how the mitotic machinery senses and responds to karyotypic changes by using a series of budding yeast strains in which the native chromosomes have been successively fused. Using a combination of cell biological profiling, genetic engineering, and experimental evolution, we show that chromosome fusions are well tolerated up until a critical point. Cells with fewer than five centromeres lack the necessary number of kinetochore-microtubule attachments needed to counter outward forces in the metaphase spindle, triggering the spindle assembly checkpoint and prolonging metaphase. Our findings demonstrate that spindle architecture is a constraining factor for karyotype evolution.

Nonhomologous tails direct heteroduplex rejection and mismatch correction during single-strand annealing in Saccharomyces cerevisiae.

Single-strand annealing (SSA) is initiated when a double strand break (DSB) occurs between two flanking repeated sequences, resulting in a deletion that leaves a single copy of the repeat. We studied budding yeast strains carrying two 200-bp URA3 sequences separated by 2.6 kb of spacer DNA (phage lambda) in which a site-specific DSB can be created by HO or Cas9 endonucleases. Repeat-mediated deletion requires removal of long 3'-ended single-stranded tails (flaps) by Rad1-Rad10 with the assistance of Msh2-Msh3, Saw1 and Slx4. A natural 3% divergence of unequally spaced heterologies between these repeats (designated F and A) causes a significant reduction in the frequency of SSA repair. This decrease is caused by heteroduplex rejection in which mismatches (MMs) in the annealed intermediate are recognized by the MutS (Msh2 and Msh6) components of the MM repair (MMR) pathway coupled to unwinding of the duplex by the Sgs1-Rmi1-Top3 helicase. MutL homologs, Mlh1-Pms1 (MutL), are not required for rejection but play their expected role in mismatch correction. Remarkably, heteroduplex rejection is very low in strains where the divergent repeats were immediately adjacent (Tailless strains) and the DSB was induced by Cas9. These results suggest that the presence of nonhomologous tails strongly stimulates heteroduplex rejection in SSA. DNA sequencing analysis of SSA products from the FA Tailed strain showed a gradient of correction favoring the sequence opposite each 3' end of the annealed strand. Mismatches located in the center of the repair intermediate were corrected by Msh2-Msh6 mediated mismatch correction, while correction of MMs at the extremity of the SSA intermediate often appears to use a different mechanism, possibly by 3' nonhomologous tail removal that includes part of the homologous sequence. In contrast, in FA Tailless strains there was a uniform repair of the MMs across the repeat. A distinctive pattern of correction was found in the absence of MSH2, in both Tailed and Tailless strains, different from the spectrum seen in a msh3Δ msh6Δ double mutant. Previous work has shown that SSA is Rad51-independent but dependent on the strand annealing activity of Rad52. However Rad52 becomes dispensable in a Tailless construct where the DSB is induced by Cas9 or in transformation of a plasmid where SSA occurs in the absence of nonhomologous tails.

Fitness effects of a demography-dispersal trade-off in expandingSaccharomyces cerevisiaemats.

Fungi expand in space and time to form complex multicellular communities. The mechanisms by which they do so can vary dramatically and determine the life-history and dispersal traits of expanding populations. These traits influence deterministic and stochastic components of evolution, resulting in complex eco-evolutionary dynamics during colony expansion. We perform experiments on budding yeast strains genetically engineered to display rough-surface and smooth-surface phenotypes in colony-like structures called 'mats'. Previously, it was shown that the rough-surface strain has a competitive advantage over the smooth-surface strain when grown on semi-solid media. We experimentally observe the emergence and expansion of segments with a distinct smooth-surface phenotype during rough-surface mat development. We propose a trade-off between dispersal and local carrying capacity to explain the relative fitness of these two phenotypes. Using a modified stepping-stone model, we demonstrate that this trade-off gives the high-dispersing, rough-surface phenotype a competitive advantage from standing variation, but that it inhibits this phenotype's ability to invade a resident smooth-surface population via mutation. However, the trade-off improves the ability of the smooth-surface phenotype to invade in rough-surface mats, replicating the frequent emergence of smooth-surface segments in experiments. Together, these computational and experimental findings advance our understanding of the complex eco-evolutionary dynamics of fungal mat expansion.

Genetic effects on molecular network states explain complex traits

The complexity of many cellular and organismal traits results from the integration of genetic and environmental factors via molecular networks. Network structure and effect propagation are best understood at the level of functional modules, but so far, no concept has been established to include the global network state. Here, we show when and how genetic perturbations lead to molecular changes that are confined to small parts of a network versus when they lead to modulation of network states. Integrating multi‐omics profiling of genetically heterogeneous budding and fission yeast strains with an array of cellular traits identified a central state transition of the yeast molecular network that is related to PKA and TOR (PT) signaling. Genetic variants affecting this PT state globally shifted the molecular network along a single‐dimensional axis, thereby modulating processes including energy and amino acid metabolism, transcription, translation, cell cycle control, and cellular stress response. We propose that genetic effects can propagate through large parts of molecular networks because of the functional requirement to centrally coordinate the activity of fundamental cellular processes.

Effects of Desiccation and Freezing on Microbial Ionizing Radiation Survivability: Considerations for Mars Sample Return.

Increasingly, national space agencies are expanding their goals to include Mars exploration with sample return. To better protect Earth and its biosphere from potential extraterrestrial sources of contamination, as set forth in the Outer Space Treaty of 1967, international efforts to develop planetary protection measures strive to understand the danger of cross-contamination processes in Mars sample return missions. We aim to better understand the impact of the martian surface on microbial dormancy and survivability. Radiation resistance of microbes is a key parameter in considering survivability of microbes over geologic times on the frigid, arid surface of Mars that is bombarded by solar and galactic cosmic radiation. We tested the influence of desiccation and freezing on the ionizing radiation survival of six model microorganisms: vegetative cells of two bacteria (Deinococcus radiodurans, Escherichia coli) and a strain of budding yeast (Saccharomyces cerevisiae); and vegetative cells and endospores of three Bacillus bacteria (B. subtilis, B. megaterium, B. thuringiensis). Desiccation and freezing greatly increased radiation survival of vegetative polyploid microorganisms when applied separately, and when combined, desiccation and freezing increased radiation survival even more so. Thus, the radiation survival threshold of polyploid D. radiodurans cells can be extended from the already high value of 25 kGy in liquid culture to an astonishing 140 kGy when the cells are both desiccated and frozen. However, such synergistic radioprotective effects of desiccation and freezing were not observed in monogenomic or digenomic Bacillus cells and endospores, which are generally sterilized by 12 kGy. This difference is associated with a critical requirement for survivability under radiation, that is, repair of genome damage caused by radiation. Deinococcus radiodurans and S. cerevisiae accumulate similarly high levels of the Mn antioxidants that are required for extreme radiation resistance, as do endospores, though they greatly exceed spores in radioresistance because they contain multiple identical genome copies, which in D. radiodurans are joined by persistent Holliday junctions. We estimate ionizing radiation survival limits of polyploid DNA-based life-forms to be hundreds of millions of years of background radiation while buried in the martian subsurface. Our findings imply that forward contamination of Mars will essentially be permanent, and backward contamination is a possibility if life ever existed on Mars.

Drug-dependent growth curve reshaping reveals mechanisms of antifungal resistance in Saccharomyces cerevisiae

Microbial drug resistance is an emerging global challenge. Current drug resistance assays tend to be simplistic, ignoring complexities of resistance manifestations and mechanisms, such as multicellularity. Here, we characterize multicellular and molecular sources of drug resistance upon deleting the AMN1 gene responsible for clumping multicellularity in a budding yeast strain, causing it to become unicellular. Computational analysis of growth curve changes upon drug treatment indicates that the unicellular strain is more sensitive to four common antifungals. Quantitative models uncover entwined multicellular and molecular processes underlying these differences in sensitivity and suggest AMN1 as an antifungal target in clumping pathogenic yeasts. Similar experimental and mathematical modeling pipelines could reveal multicellular and molecular drug resistance mechanisms, leading to more effective treatments against various microbial infections and possibly even cancers.

Electromagnetic field exposure increases the protein misfolding in vivo

Extremely low frequency electromagnetic field (ELF-EMF) exposure can elevate the risk of neurodegenerative diseases, which are closely related with protein misfolding. We used yeast prion as model to analyze whether the exposure to ELF-EMF has effect on protein folding/misfolding. [URE3] is a yeast prion phenotype, in which its prion determinant protein, Ure2, is misfolded and aggregated. We used two budding yeast strains, NT64C and SB34, both of which contained [URE3] prion form but their genetic background was quite different.

Purine Auxotrophic Starvation Evokes Phenotype Similar to Stationary Phase Cells in Budding Yeast.

Purine auxotrophy is an abundant trait among eukaryotic parasites and a typical marker for many budding yeast strains. Supplementation with an additional purine source (such as adenine) is necessary to cultivate these strains. If not supplied in adequate amounts, purine starvation sets in. We explored purine starvation effects in a model organism, a budding yeast Saccharomyces cerevisiae ade8 knockout, at the level of cellular morphology, central carbon metabolism, and global transcriptome. We observed that purine-starved cells stopped their cycle in G1/G0 state and accumulated trehalose, and the intracellular concentration of AXP decreased, but adenylate charge remained stable. Cells became tolerant to severe environmental stresses. Intracellular RNA concentration decreased, and massive downregulation of ribosomal biosynthesis genes occurred. We proved that the expression of new proteins during purine starvation is critical for cells to attain stress tolerance phenotype Msn2/4p targets are upregulated in purine-starved cells when compared to cells cultivated in purine-rich media. The overall transcriptomic response to purine starvation resembles that of stationary phase cells. Our results demonstrate that the induction of a strong stress resistance phenotype in budding yeast can be caused not only by natural starvation, but also starvation for metabolic intermediates, such as purines.

Isolation, identification, and characterization of wild budding yeasts from rose flowers in Fukuyama city, Hiroshima, Japan, and their application in bread and wine production.

In this study, we isolated 741 wild budding yeast strains from the flowers of 45 rose cultivars growing in Fukuyama city, Hiroshima, Japan. Of these 741 strains, 21 were found to have high fermentation abilities in yeast extract-peptone-dextrose (YPD) medium. Four of the 21 strains were able to ferment bread dough to make bread. These yeasts were identified as Saccharomyces cerevisiae, Lachancea fermentati, Lachancea kluyveri, and a Torulaspora sp. based on DNA sequences from the 26S rDNA D1/D2 regions. The CO2 production profiles of the bread dough generated by the rose yeasts were evaluated using a Fermograph. Saccharomyces cerevisiae FRY2915 exhibited the highest fermentation capacity. Furthermore, FRY2915 was able to ferment grape juice to produce wine, yielding an alcohol concentration of more than 12%. The four rose yeasts isolated during this study have the potential to produce various types of unique fermented foods, thus enhancing the value of the microbiota associated with rose flowers.

SUMO orchestrates multiple alternative DNA-protein crosslink repair pathways.

Endogenous metabolites, environmental agents, and therapeutic drugs promote formation of covalent DNA-protein crosslinks (DPCs). Persistent DPCs compromise genome integrity and are eliminated by multiple repair pathways. Aberrant Top1-DNA crosslinks, or Top1ccs, are processed by Tdp1 and Wss1 functioning in parallel pathways in Saccharomyces cerevisiae. It remains obscure how cells choose between diverse mechanisms of DPC repair. Here, we show that several SUMO biogenesis factors (Ulp1, Siz2, Slx5, and Slx8) control repair of Top1cc or an analogous DPC lesion. Genetic analysis reveals that SUMO promotes Top1cc processing in the absence of Tdp1 but has an inhibitory role if cells additionally lack Wss1. In the tdp1Δ wss1Δ mutant, the E3 SUMO ligase Siz2 stimulates sumoylation in the vicinity of the DPC, but not SUMO conjugation to Top1. This Siz2-dependent sumoylation inhibits alternative DPC repair mechanisms, including Ddi1. Our findings suggest that SUMO tunes available repair pathways to facilitate faithful DPC repair.

Comparison between instrumental variable and mediation-based methods for reconstructing causal gene networks in yeast.

Causal gene networks model the flow of information within a cell. Reconstructing causal networks from omics data is challenging because correlation does not imply causation. When genomics and transcriptomics data from a segregating population are combined, genomic variants can be used to orient the direction of causality between gene expression traits. Instrumental variable methods use a local expression quantitative trait locus (eQTL) as a randomized instrument for a gene's expression level, and assign target genes based on distal eQTL associations. Mediation-based methods additionally require that distal eQTL associations are mediated by the source gene. A detailed comparison between these methods has not yet been conducted, due to the lack of a standardized implementation of different methods, the limited sample size of most multi-omics datasets, and the absence of ground-truth networks for most organisms. Here we used Findr, a software package providing uniform implementations of instrumental variable, mediation, and coexpression-based methods, a recent dataset of 1012 segregants from a cross between two budding yeast strains, and the Yeastract database of known transcriptional interactions to compare causal gene network inference methods. We found that causal inference methods result in a significant overlap with the ground-truth, whereas coexpression did not perform better than random. A subsampling analysis revealed that the performance of mediation saturates at large sample sizes, due to a loss of sensitivity when residual correlations become significant. Instrumental variable methods on the other hand contain false positive predictions, due to genomic linkage between eQTL instruments. Instrumental variable and mediation-based methods also have complementary roles for identifying causal genes underlying transcriptional hotspots. Instrumental variable methods correctly predicted STB5 targets for a hotspot centred on the transcription factor STB5, whereas mediation failed due to Stb5p auto-regulating its own expression. Mediation suggests a new candidate gene, DNM1, for a hotspot on Chr XII, whereas instrumental variable methods could not distinguish between multiple genes located within the hotspot. In conclusion, causal inference from genomics and transcriptomics data is a powerful approach for reconstructing causal gene networks, which could be further improved by the development of methods to control for residual correlations in mediation analyses, and for genomic linkage and pleiotropic effects from transcriptional hotspots in instrumental variable analyses.

A C6/C5 co‐fermentingSaccharomyces cerevisiaestrain with the alleviation of antagonism between xylose utilization and robustness

AbstractDuring second‐generation bioethanol production from lignocellulosic biomass, the desired traits for fermenting microorganisms, such asSaccharomyces cerevisiae, are high xylose utilization and high robustness to inhibitors in lignocellulosic hydrolysates. However, as observed previously, these two traits easily showed the antagonism, one rising and the other falling, in the C6/C5 co‐fermentingS. cerevisiaestrain. In this study, LF1 obtained in our previous study is an engineered budding yeast strain with a superior co‐fermentation capacity of glucose and xylose, and was then mutated by atmospheric and room temperature plasma (ARTP) mutagenesis to improve its robustness. The ARTP‐treated cells were grown in 50% (v/v) leachate from lignocellulose pretreatment with high inhibitors content for adaptive evolution. After 30 days, the generated mutant LF1‐6 showed significantly enhanced tolerance, with a six‐fold increase in cell density in the above leachate. Unfortunately, its xylose utilization dropped markedly, indicating the recurrence of the negative correlation between xylose utilization and robustness. To alleviate this antagonism, LF1‐6 cells were iteratively mutated with ARTP mutagenesis and then anaerobically grown using xylose as the sole carbon source, and xylose utilization was restored in the resulting strain 6M‐15. 6M‐15 also exhibited increased co‐fermentation performance of xylose and glucose with the highest ethanol productivity reported to date (0.525 g g−1 h−1) in high‐level mixed sugars (80 g L−1glucose and 40 g L−1xylose) with no inhibitors. Meanwhile, its fermentation time was shortened by 8 h compared to that of LF1. During the fermentation of non‐detoxified lignocellulosic hydrolysate with high inhibitor concentrations at pH ~3.5, 6M‐15 can efficiently convert glucose and xylose with an ethanol yield of 0.43 g g−1. 6M‐15 is also regarded as a potential chassis cell for further design of a customized strain suitable for production of second‐generation bioethanol or other high value‐added products from lignocellulosic biomass.

Getting to grips with circular chromosomes.

A strain of budding yeast that contains one large chromosome reveals how the telomere capping complex CST maintains linear but not circular chromosomes.

Functional Interactions of Ribosomal Intersubunit Bridges in Saccharomyces cerevisiae.

Ribosomes of Archaea and Eukarya share higher homology with each other than with bacterial ribosomes. For example, there is a set of 35 r-proteins that are specific only for archaeal and eukaryotic ribosomes. Three of these proteins-eL19, eL24, and eL41-participate in interactions between ribosomal subunits. The eukaryote-specific extensions of r-proteins eL19 and eL24 form two intersubunit bridges eB12 and eB13, which are present only in eukaryotic ribosomes. The third r-protein, eL41, forms bridge eB14. Notably, eL41 is found in all eukaryotes but only in some Archaea. It has been shown that bridges eB12 and eB13 are needed for efficient translation, while r-protein eL41 plays a minor role in ribosome function. Here, the functional interactions between intersubunit bridges were studied using budding yeast strains lacking different combinations of the abovementioned bridges/proteins. The growth phenotypes, levels of in vivo translation, ribosome-polysome profiles, and in vitro association of ribosomal subunits were analyzed. The results show a genetic interaction between r-protein eL41 and the eB12 bridge-forming region of eL19, and between r-proteins eL41 and eL24. It was possible to construct viable yeast strains with Archaea-like ribosomes lacking two or three eukaryote-specific bridges. These strains display slow growth and a poor translation phenotype. In addition, bridges eB12 and eB13 appear to cooperate during ribosome subunit association. These results indicate that nonessential structural elements of r-proteins become highly important in the context of disturbed subunit interactions. Therefore, eukaryote-specific bridges may contribute to the evolutionary success of eukaryotic translation machinery.

Paralog buffering contributes to the variable essentiality of genes in cancer cell lines.

What makes a gene essential for cellular survival? In model organisms, such as budding yeast, systematic gene deletion studies have revealed that paralog genes are less likely to be essential than singleton genes and that this can partially be attributed to the ability of paralogs to buffer each other's loss. However, the essentiality of a gene is not a fixed property and can vary significantly across different genetic backgrounds. It is unclear to what extent paralogs contribute to this variation, as most studies have analyzed genes identified as essential in a single genetic background. Here, using gene essentiality profiles of 558 genetically heterogeneous tumor cell lines, we analyze the contribution of paralogy to variable essentiality. We find that, compared to singleton genes, paralogs are less frequently essential and that this is more evident when considering genes with multiple paralogs or with highly sequence-similar paralogs. In addition, we find that paralogs derived from whole genome duplication exhibit more variable essentiality than those derived from small-scale duplications. We provide evidence that in 13-17% of cases the variable essentiality of paralogs can be attributed to buffering relationships between paralog pairs, as evidenced by synthetic lethality. Paralog pairs derived from whole genome duplication and pairs that function in protein complexes are significantly more likely to display such synthetic lethal relationships. Overall we find that many of the observations made using a single strain of budding yeast can be extended to understand patterns of essentiality in genetically heterogeneous cancer cell lines.

Comprehensive Analysis of Replication Origins in Saccharomyces cerevisiae Genomes.

DNA replication initiates from multiple replication origins (ORIs) in eukaryotes. Discovery and characterization of replication origins are essential for a better understanding of the molecular mechanism of DNA replication. In this study, the features of autonomously replicating sequences (ARSs) in Saccharomyces cerevisiae have been comprehensively analyzed as follows. Firstly, we carried out the analysis of the ARSs available in S. cerevisiae S288C. By evaluating the sequence similarity of experimentally established ARSs, we found that 94.32% of ARSs are unique across the whole genome of S. cerevisiae S288C and those with high sequence similarity are prone to locate in subtelomeres. Subsequently, we built a non-redundant dataset with a total of 520 ARSs, which are based on ARSs annotation of S. cerevisiae S288C from SGD and then supplemented with those from OriDB and DeOri databases. We conducted a large-scale comparison of ORIs among the diverse budding yeast strains from a population genomics perspective. We found that 82.7% of ARSs are not only conserved in genomic sequence but also relatively conserved in chromosomal position. The non-conserved ARSs tend to distribute in the subtelomeric regions. We also conducted a pan-genome analysis of ARSs among the S. cerevisiae strains, and a total of 183 core ARSs existing in all yeast strains were determined. We extracted the genes adjacent to replication origins among the 104 yeast strains to examine whether there are differences in their gene functions. The result showed that the genes involved in the initiation of DNA replication, such as orc3, mcm2, mcm4, mcm6, and cdc45, are conservatively located adjacent to the replication origins. Furthermore, we found the genes adjacent to conserved ARSs are significantly enriched in DNA binding, enzyme activity, transportation, and energy, whereas for the genes adjacent to non-conserved ARSs are significantly enriched in response to environmental stress, metabolites biosynthetic process and biosynthesis of antibiotics. In general, we characterized the replication origins from the genome-wide and population genomics perspectives, which would provide new insights into the replication mechanism of S. cerevisiae and facilitate the design of algorithms to identify genome-wide replication origins in yeast.