The spatiotemporal program of replication in the genome of Lachancea kluyveri

Absence of chromosome axis protein recruitment prevents meiotic recombination chromosome-wide in the budding yeast Lachancea kluyveri

Meiotic recombination shows broad variations across species and along chromosomes and is often suppressed at and around genomic regions determining sexual compatibility such as mating type loci in fungi. Here, we show that the absence of Spo11-DSBs and meiotic recombination on Lakl0C-left, the chromosome arm containing the sex locus of the Lachancea kluyveri budding yeast, results from the absence of recruitment of the two chromosome axis proteins Red1 and Hop1, essential for proper Spo11-DSBs formation. Furthermore, cytological observation of spread pachytene meiotic chromosomes reveals that Lakl0C-left does not undergo synapsis. However, we show that the behavior of Lakl0C-left is independent of its particularly early replication timing and is not accompanied by any peculiar chromosome structure as detectable by Hi-C in this yet poorly studied yeast. Finally, we observed an accumulation of heterozygous mutations on Lakl0C-left and a sexual dimorphism of the haploid meiotic offspring, supporting a direct effect of this absence of meiotic recombination on L. kluyveri genome evolution and fitness. Because suppression of meiotic recombination on sex chromosomes is widely observed across eukaryotes, the mechanism for recombination suppression described here may apply to other species, with the potential to impact sex chromosome evolution.

An evolutionary model identifies the main evolutionary biases for the evolution of genome-replication profiles

Recent results comparing the temporal program of genome replication of yeast species belonging to the Lachancea clade support the scenario that the evolution of the replication timing program could be mainly driven by correlated acquisition and loss events of active replication origins. Using these results as a benchmark, we develop an evolutionary model defined as birth-death process for replication origins and use it to identify the evolutionary biases that shape the replication timing profiles. Comparing different evolutionary models with data, we find that replication origin birth and death events are mainly driven by two evolutionary pressures, the first imposes that events leading to higher double-stall probability of replication forks are penalized, while the second makes less efficient origins more prone to evolutionary loss. This analysis provides an empirically grounded predictive framework for quantitative evolutionary studies of the replication timing program.

Lachancea fermentati Strains Isolated From Kombucha: Fundamental Insights, and Practical Application in Low Alcohol Beer Brewing

With a growing interest in non-alcoholic and low alcohol beer (NABLAB), researchers are looking into non-conventional yeasts to harness their special metabolic traits for their production. One of the investigated species is Lachancea fermentati, which possesses the uncommon ability to produce significant amounts of lactic acid during alcoholic fermentation, resulting in the accumulation of lactic acid while exhibiting reduced ethanol production. In this study, four Lachancea fermentati strains isolated from individual kombucha cultures were investigated. Whole genome sequencing was performed, and the strains were characterized for important brewing characteristics (e.g., sugar utilization) and sensitivities toward stress factors. A screening in wort extract was performed to elucidate strain-dependent differences, followed by fermentation optimization to enhance lactic acid production. Finally, a low alcohol beer was produced at 60 L pilot-scale. The genomes of the kombucha isolates were diverse and could be separated into two phylogenetic groups, which were related to their geographical origin. Compared to a Saccharomyces cerevisiae brewers' yeast, the strains' sensitivities to alcohol and acidic conditions were low, while their sensitivities toward osmotic stress were higher. In the screening, lactic acid production showed significant, strain-dependent differences. Fermentation optimization by means of response surface methodology (RSM) revealed an increased lactic acid production at a low pitching rate, high fermentation temperature, and high extract content. It was shown that a high initial glucose concentration led to the highest lactic acid production (max. 18.0 mM). The data indicated that simultaneous lactic acid production and ethanol production occurred as long as glucose was present. When glucose was depleted and/or lactic acid concentrations were high, the production shifted toward the ethanol pathway as the sole pathway. A low alcohol beer (<1.3% ABV) was produced at 60 L pilot-scale by means of stopped fermentation. The beer exhibited a balanced ratio of sweetness from residual sugars and acidity from the lactic acid produced (13.6 mM). However, due to the stopped fermentation, high levels of diacetyl were present, which could necessitate further process intervention to reduce concentrations to acceptable levels.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

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.

Discovery and Evolution of New Domains in Yeast Heterochromatin Factor Sir4 and Its Partner Esc1

Sir4 is a core component of heterochromatin found in yeasts of the Saccharomycetaceae family, whose general hallmark is to harbor a three-loci mating-type system with two silent loci. However, a large part of the Sir4 amino acid sequences has remained unexplored, belonging to the dark proteome. Here, we analyzed the phylogenetic profile of yet undescribed foldable regions present in Sir4 as well as in Esc1, an Sir4-interacting perinuclear anchoring protein. Within Sir4, we identified a new conserved motif (TOC) adjacent to the N-terminal KU-binding motif. We also found that the Esc1-interacting region of Sir4 is a Dbf4-related H-BRCT domain, only present in species possessing the HO endonuclease and in Kluveryomyces lactis. In addition, we found new motifs within Esc1 including a motif (Esc1-F) that is unique to species where Sir4 possesses an H-BRCT domain. Mutagenesis of conserved amino acids of the Sir4 H-BRCT domain, known to play a critical role in the Dbf4 function, shows that the function of this domain is separable from the essential role of Sir4 in transcriptional silencing and the protection from HO-induced cutting in Saccharomyces cerevisiae. In the more distant methylotrophic clade of yeasts, which often harbor a two-loci mating-type system with one silent locus, we also found a yet undescribed H-BRCT domain in a distinct protein, the ISWI2 chromatin-remodeling factor subunit Itc1. This study provides new insights on yeast heterochromatin evolution and emphasizes the interest of using sensitive methods of sequence analysis for identifying hitherto ignored functional regions within the dark proteome.

André Goffeau's imprinting on second generation yeast "genomologists"

All authors of the present paper have worked in labs that participated to the sequencing effort of the Saccharomyces cerevisiae reference genome, and we owe to this the fact that we have all chosen to work on genomics of yeasts. S. cerevisiae has been a popular model species for genetics since the 20th century as well as being a model for general eukaryotic cellular processes. Although it has also been used empirically in fermentation for millennia, there was until recently, a lack of knowledge about the natural and evolutionary history of this yeast. The achievement of the international effort to sequence its genome was the foundation for understanding many eukaryotic biological processes but also represented the first step towards the study of the genome and ecological diversity of yeast populations worldwide. We will describe recent advances in yeast comparative and population genomics that find their origins in the S. cerevisiae genome project initiated and pursued by André Goffeau.

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.

Cell-to-cell variability and robustness in S-phase duration from genome replication kinetics

Genome replication, a key process for a cell, relies on stochastic initiation by replication origins, causing a variability of replication timing from cell to cell. While stochastic models of eukaryotic replication are widely available, the link between the key parameters and overall replication timing has not been addressed systematically. We use a combined analytical and computational approach to calculate how positions and strength of many origins lead to a given cell-to-cell variability of total duration of the replication of a large region, a chromosome or the entire genome. Specifically, the total replication timing can be framed as an extreme-value problem, since it is due to the last region that replicates in each cell. Our calculations identify two regimes based on the spread between characteristic completion times of all inter-origin regions of a genome. For widely different completion times, timing is set by the single specific region that is typically the last to replicate in all cells. Conversely, when the completion time of all regions are comparable, an extreme-value estimate shows that the cell-to-cell variability of genome replication timing has universal properties. Comparison with available data shows that the replication program of three yeast species falls in this extreme-value regime.

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.

Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.

A checkpoint control orchestrates the replication of the two chromosomes of Vibrio cholerae

Bacteria with multiple chromosomes represent up to 10% of all bacterial species. Unlike eukaryotes, these bacteria use chromosome-specific initiators for their replication. In all cases investigated, the machineries for secondary chromosome replication initiation are of plasmid origin. One of the important differences between plasmids and chromosomes is that the latter replicate during a defined period of the cell cycle, ensuring a single round of replication per cell. Vibrio cholerae carries two circular chromosomes, Chr1 and Chr2, which are replicated in a well-orchestrated manner with the cell cycle and coordinated in such a way that replication termination occurs at the same time. However, the mechanism coordinating this synchrony remains speculative. We investigated this mechanism and revealed that initiation of Chr2 replication is triggered by the replication of a 150-bp locus positioned on Chr1, called crtS. This crtS replication-mediated Chr2 replication initiation mechanism explains how the two chromosomes communicate to coordinate their replication. Our study reveals a new checkpoint control mechanism in bacteria, and highlights possible functional interactions mediated by contacts between two chromosomes, an unprecedented observation in bacteria.

A Versatile Procedure to Generate Genome-Wide Spatiotemporal Program of Replication in Yeast Species

Here, we describe a complete protocol, comprising both the experimental and the analytical procedures, that allows to generate genome-wide spatiotemporal program of replication and to find the location of chromosomally active replication origins in yeast. The first step consists on synchronizing a cell population by physical discrimination of G1 cells according to their sedimentation coefficient. G1 cells are then synchronously released into S-phase and time-point samples are regularly taken until they reach the G2 phase. Progression through the cell cycle is monitored by measuring DNA content variation by flow cytometry. DNA samples, covering the entire S-phase, are then extracted and analyzed using deep sequencing. The gradual change of DNA copy number is measured to determine the mean replication time along the genome. A simple method of peak calling allows to infer from the replication profile the location of replication origins along the chromosomes. Our protocol is versatile enough to be applied to virtually any yeast species of interest and generate its replication profile.

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.

Total synthesis of a functional designer eukaryotic chromosome

Rapid advances in DNA synthesis techniques have made it possible to engineer viruses, biochemical pathways and assemble bacterial genomes. Here, we report the synthesis of a functional 272,871-base pair designer eukaryotic chromosome, synIII, which is based on the 316,617-base pair native Saccharomyces cerevisiae chromosome III. Changes to synIII include TAG/TAA stop-codon replacements, deletion of subtelomeric regions, introns, transfer RNAs, transposons, and silent mating loci as well as insertion of loxPsym sites to enable genome scrambling. SynIII is functional in S. cerevisiae. Scrambling of the chromosome in a heterozygous diploid reveals a large increase in a-mater derivatives resulting from loss of the MATα allele on synIII. The complete design and synthesis of synIII establishes S. cerevisiae as the basis for designer eukaryotic genome biology.

Purification of G1 daughter cells from different Saccharomycetes species through an optimized centrifugal elutriation procedure

Centrifugal elutriation discriminates cells according to their sedimentation coefficients, generating homogeneous samples well suited for genomic comparative approaches. It can, for instance, isolate G1 daughter cells from a Saccharomyces cerevisiae unsynchronized population, alleviating ageing and cell-cycle biases when conducting genome-wide/single-cell studies. The present report describes a straightforward and robust procedure to determine whether a cell population of virtually any yeast species can be efficiently elutriated, while offering solutions to optimize success. This approach was used to characterize elutriation parameters and S-phase progression of four yeast species (S. cerevisiae, Candida glabrata, Lachancea kluyveri and Pichia sorbitophila) and could theoretically be applied to any culture of single, individual cells.

The dynamics of genome replication using deep sequencing

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.

Replisome stall events have shaped the distribution of replication origins in the genomes of yeasts

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.

Global expression analysis of the yeast Lachancea (Saccharomyces) kluyveri reveals new URC genes involved in pyrimidine catabolism

Pyrimidines are important nucleic acid precursors which are constantly synthesized, degraded, and rebuilt in the cell. Four degradation pathways, two of which are found in eukaryotes, have been described. One of them, the URC pathway, has been initially discovered in our laboratory in the yeast Lachancea kluyveri. Here, we present the global changes in gene expression in L. kluyveri in response to different nitrogen sources, including uracil, uridine, dihydrouracil, and ammonia. The expression pattern of the known URC genes, URC1-6, helped to identify nine putative novel URC genes with a similar expression pattern. The microarray analysis provided evidence that both the URC and PYD genes are under nitrogen catabolite repression in L. kluyveri and are induced by uracil or dihydrouracil, respectively. We determined the function of URC8, which was found to catalyze the reduction of malonate semialdehyde to 3-hydroxypropionate, the final degradation product of the pathway. The other eight genes studied were all putative permeases. Our analysis of double deletion strains showed that the L. kluyveri Fui1p protein transported uridine, just like its homolog in Saccharomyces cerevisiae, but we demonstrated that is was not the only uridine transporter in L. kluyveri. We also showed that the L. kluyveri homologs of DUR3 and FUR4 do not have the same function that they have in S. cerevisiae, where they transport urea and uracil, respectively. In L. kluyveri, both of these deletion strains grew normally on uracil and urea.