Whole-Genome Sequence and Variant Analysis of W303, a Widely-Used Strain of Saccharomyces cerevisiae

Genetically Encoded, Noise-Tolerant, Auxin Biosensors in Yeast

Auxins are crucial signaling molecules that regulate the growth, metabolism, and behavior of various organisms, most notably plants but also bacteria, fungi, and animals. Many microbes synthesize and perceive auxins, primarily indole-3-acetic acid (IAA, referred to as auxin herein), the most prevalent natural auxin, which influences their ability to colonize plants and animals. Understanding auxin biosynthesis and signaling in fungi may allow us to better control interkingdom relationships and microbiomes from agricultural soils to the human gut. Despite this importance, a biological tool for measuring auxin with high spatial and temporal resolution has not been engineered in fungi. In this study, we present a suite of genetically encoded, ratiometric, protein-based auxin biosensors designed for the model yeast Saccharomyces cerevisiae. Inspired by auxin signaling in plants, the ratiometric nature of these biosensors enhances the precision of auxin concentration measurements by minimizing clonal and growth phase variation. We used these biosensors to measure auxin production across diverse growth conditions and phases in yeast cultures and calibrated their responses to physiologically relevant levels of auxin. Future work will aim to improve the fold change and reversibility of these biosensors. These genetically encoded auxin biosensors are valuable tools for investigating auxin biosynthesis and signaling in S. cerevisiae and potentially other yeast and fungi and will also advance quantitative functional studies of the plant auxin perception machinery, from which they are built.

An extrinsic motor directs chromatin loop formation by cohesin

The ring-shaped cohesin complex topologically entraps two DNA molecules to establish sister chromatid cohesion. Cohesin also shapes the interphase chromatin landscape with wide-ranging implications for gene regulation, and cohesin is thought to achieve this by actively extruding DNA loops without topologically entrapping DNA. The 'loop extrusion' hypothesis finds motivation from in vitro observations-whether this process underlies in vivo chromatin loop formation remains untested. Here, using the budding yeast S. cerevisiae, we generate cohesin variants that have lost their ability to extrude DNA loops but retain their ability to topologically entrap DNA. Analysis of these variants suggests that in vivo chromatin loops form independently of loop extrusion. Instead, we find that transcription promotes loop formation, and acts as an extrinsic motor that expands these loops and defines their ultimate positions. Our results necessitate a re-evaluation of the loop extrusion hypothesis. We propose that cohesin, akin to sister chromatid cohesion establishment at replication forks, forms chromatin loops by DNA-DNA capture at places of transcription, thus unifying cohesin's two roles in chromosome segregation and interphase genome organisation.

Dbf4-dependent kinase promotes cell cycle controlled resection of DNA double-strand breaks and repair by homologous recombination

DNA double-strand breaks (DSBs) can be repaired by several pathways. In eukaryotes, DSB repair pathway choice occurs at the level of DNA end resection and is controlled by the cell cycle. Upon cell cycle-dependent activation, cyclin-dependent kinases (CDKs) phosphorylate resection proteins and thereby stimulate end resection and repair by homologous recombination (HR). However, inability of CDK phospho-mimetic mutants to bypass this cell cycle regulation, suggests that additional cell cycle regulators may be important. Here, we identify Dbf4-dependent kinase (DDK) as a second major cell cycle regulator of DNA end resection. Using inducible genetic and chemical inhibition of DDK in budding yeast and human cells, we show that end resection and HR require activation by DDK. Mechanistically, DDK phosphorylates at least two resection nucleases in budding yeast: the Mre11 activator Sae2, which promotes resection initiation, as well as the Dna2 nuclease, which promotes resection elongation. Notably, synthetic activation of DDK allows limited resection and HR in G1 cells, suggesting that DDK is a key component of DSB repair pathway selection.

Quantifying turgor pressure in budding and fission yeasts based upon osmotic properties

Walled cells, such as plants, fungi, and bacteria cells, possess a high internal hydrostatic pressure, termed turgor pressure, that drives volume growth and contributes to cell shape determination. Rigorous measurement of turgor pressure, however, remains challenging, and reliable quantitative measurements, even in budding yeast are still lacking. Here, we present a simple and robust experimental approach to access turgor pressure in yeasts based upon the determination of isotonic concentration using protoplasts as osmometers. We propose three methods to identify the isotonic condition - three-dimensional cell volume, cytoplasmic fluorophore intensity, and mobility of a cytGEMs nano-rheology probe - that all yield consistent values. Our results provide turgor pressure estimates of 1.0 ± 0.1 MPa for Schizosaccharomyces pombe, 0.49 ± 0.01 MPa for Schizosaccharomyces japonicus, 0.5 ± 0.1 MPa for Saccharomyces cerevisiae W303a and 0.31 ± 0.03 MPa for Saccharomyces cerevisiae BY4741. Large differences in turgor pressure and nano-rheology measurements between the Saccharomyces cerevisiae strains demonstrate how fundamental biophysical parameters can vary even among wild-type strains of the same species. These side-by-side measurements of turgor pressure in multiple yeast species provide critical values for quantitative studies on cellular mechanics and comparative evolution.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease activity. Foldback priming at DNA double-strand breaks (DSBs) is one mechanism proposed for the generation of inverted duplications. However, the other pathway steps have not been fully elucidated. Here, we show that a DSB induced near natural inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity, but not Rad1 nuclease, trims the heterologous flaps formed after foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric inversion chromosome. Finally, we show that stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences or by deletion of one centromere.

Quantifying turgor pressure in budding and fission yeasts based upon osmotic properties

Walled cells, such as plants, fungi, and bacteria cells, possess a high internal hydrostatic pressure, termed turgor pressure, that drives volume growth and contributes to cell shape determination. Rigorous measurement of turgor pressure, however, remains challenging, and reliable quantitative measurements, even in budding yeast are still lacking. Here, we present a simple and robust experimental approach to access turgor pressure in yeasts based upon the determination of isotonic concentration using protoplasts as osmometers. We propose three methods to identify the isotonic condition - 3D cell volume, cytoplasmic fluorophore intensity, and mobility of a cytGEMs nano-rheology probe - that all yield consistent values. Our results provide turgor pressure estimates of 1.0 ± 0.1 MPa for S. pombe, 0.49 ± 0.01 MPa for S. japonicus, 0.5 ± 0.1 MPa for S. cerevisiae W303a and 0.31 ± 0.03 MPa for S. cerevisiae BY4741. Large differences in turgor pressure and nano-rheology measurements between the S. cerevisiae strains demonstrate how fundamental biophysical parameters can vary even among wildtype strains of the same species. These side-by-side measurements of turgor pressure in multiple yeast species provide critical values for quantitative studies on cellular mechanics and comparative evolution.

Sequence and chromatin features guide DNA double-strand break resection initiation

DNA double-strand breaks (DSBs) are cytotoxic genome lesions that must be accurately and efficiently repaired to ensure genome integrity. In yeast, the Mre11-Rad50-Xrs2 (MRX) complex nicks 5'-terminated DSB ends to initiate nucleolytic processing of DSBs for repair by homologous recombination. How MRX-DNA interactions support 5' strand-specific nicking and how nicking is influenced by the chromatin context have remained elusive. Using a deep sequencing-based assay, we mapped MRX nicks at single-nucleotide resolution next to multiple DSBs in the yeast genome. We observed that the DNA end-binding Ku70-Ku80 complex directed DSB-proximal nicks and that repetitive MRX cleavage extended the length of resection tracts. We identified a sequence motif and a DNA meltability profile that is preferentially nicked by MRX. Furthermore, we found that nucleosomes as well as transcription impeded MRX incisions. Our findings suggest that local DNA sequence and chromatin features shape the activity of this central DSB repair complex.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ and Rad51-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease. Foldback priming at naturally occurring inverted repeats is one mechanism proposed for the generation of inverted duplications. However, the initiating lesion for these events and the mechanism by which they form has not been fully elucidated. Here, we show that a DNA double-strand break (DSB) induced near natural short, inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity acts non-redundantly with Rad1 nuclease to remove heterologous tails formed during foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric isochromosome. Stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences and requires Rad51.

Long-Term Adaptation to Galactose as a Sole Carbon Source Selects for Mutations Outside the Canonical GAL Pathway

Galactose is a secondary fermentable sugar that requires specific regulatory and structural genes for its assimilation, which are under catabolite repression by glucose. When glucose is absent, the catabolic repression is attenuated, and the structural GAL genes are fully activated. In Saccharomyces cerevisiae, the GAL pathway is under selection in environments where galactose is present. However, it is unclear the adaptive strategies in response to long-term propagation in galactose as a sole carbon source in laboratory evolution experiments. Here, we performed a 4,000-generation evolution experiment using 48 diploid Saccharomyces cerevisiae populations to study adaptation in galactose. We show that fitness gains were greater in the galactose-evolved population than in identically evolved populations with glucose as a sole carbon source. Whole-genome sequencing of 96 evolved clones revealed recurrent de novo single nucleotide mutations in candidate targets of selection, copy number variations, and ploidy changes. We find that most mutations that improve fitness in galactose lie outside of the canonical GAL pathway. Reconstruction of specific evolved alleles in candidate target of selection, SEC23 and IRA1, showed a significant increase in fitness in galactose compared to glucose. In addition, most of our evolved populations (28/46; 61%) fixed aneuploidies on Chromosome VIII, suggesting a parallel adaptive amplification. Finally, we show greater loss of extrachromosomal elements in our glucose-evolved lineages compared with previous glucose evolution. Broadly, these data further our understanding of the evolutionary pressures that drive adaptation to less-preferred carbon sources.

Prediction of strain engineerings that amplify recombinant protein secretion through the machine learning approach MaLPHAS

This article presents a discussion of the process of precision fermentation (PF), describing the history of the space, the expected 70% growth over the next 5 years, various applications of precision fermented products, and the markets available to be disrupted by the technology. A range of prokaryotic and eukaryotic host organisms used for PF are described, with the advantages, disadvantages and applications of each. The process of setting up PF and strain engineering is described, as well as various ways that computational analysis and design techniques can be employed to assist PF engineering. The article then describes the design and implementation of a machine learning method, machine learning predictions having amplified secretion (MaLPHAS) to predict strain engineerings, which optimise the secretion of a recombinant protein. This approach showed an in silico cross-validated R 2 accuracy on the training data of up to 46.6% and in an in vitro test on a Komagataella phaffii strain, identified one gene engineering out of five predicted, which was shown to double the secretion of a heterologous protein and outperform three of the best-known edits from the literature for improving secretion in K. phaffii.

Nab3 nuclear granule accumulation is driven by respiratory capacity

Numerous biological processes involve proteins capable of transiently assembling into subcellular compartments necessary for cellular functions. One process is the RNA polymerase II transcription cycle which involves initiation, elongation, co-transcriptional modification of nascent RNA, and termination. The essential yeast transcription termination factor Nab3 is required for termination of small non-coding RNAs and accumulates into a compact nuclear granule upon glucose removal. Nab3 nuclear granule accumulation varies in penetrance across yeast strains and a higher Nab3 granule accumulation phenotype is associated with petite strains, suggesting a possible ATP-dependent mechanism for granule disassembly. Here, we demonstrate the uncoupling of mitochondrial oxidative phosphorylation by drug treatment or deletions of nuclear-encoded ATP synthase subunit genes were sufficient to increase Nab3 granule accumulation and led to an inability to proliferate during prolonged glucose deprivation, which requires respiration. Additionally, by enriching for respiration competent cells from a petite-prone strain, we generated a low granule-accumulating strain from a relatively high one, providing another link between respiratory competency and Nab3 granules. Consistent with the resulting idea that ATP is involved in granule accumulation, the addition of extracellular ATP to semi-permeabilized cells was sufficient to reduce Nab3 granule accumulation. Deleting the SKY1 gene, which encodes a kinase that phosphorylates nuclear SR repeat-containing proteins and is involved in efficient stress granule disassembly, also resulted in increased granule accumulation. This observation implicates Sky1 in Nab3 granule biogenesis. Taken together, these findings suggest there is normally an equilibrium between termination factor granule assembly and disassembly mediated by ATP-requiring nuclear machinery.

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

Sir2 and Reb1 antagonistically regulate nucleosome occupancy in subtelomeric X-elements and repress TERRAs by distinct mechanisms

Telomere chromatin structure is pivotal for maintaining genome stability by regulating the binding of telomere-associated proteins and inhibiting the DNA damage response. In Saccharomyces cerevisiae, silent information regulator (Sir) proteins bind to terminal repeats and to subtelomeric X-elements, resulting in transcriptional silencing. Herein, we show that sir2 mutant strains display a specific loss of a nucleosome residing in the X-elements and that this deficiency is remarkably consistent between different telomeres. The X-elements contain several binding sites for the transcription factor Reb1 and we found that Sir2 and Reb1 compete for stabilizing/destabilizing this nucleosome, i.e. inactivation of Reb1 in a sir2 background reinstated the lost nucleosome. The telomeric-repeat-containing RNAs (TERRAs) originate from subtelomeric regions and extend into the terminal repeats. Both Sir2 and Reb1 repress TERRAs and in a sir2 reb1 double mutant, TERRA levels increased synergistically, showing that Sir2 and Reb1 act in different pathways for repressing TERRAs. We present evidence that Reb1 restricts TERRAs by terminating transcription. Mapping the 5′-ends of TERRAs from several telomeres revealed that the Sir2-stabilized nucleosome is the first nucleosome downstream from the transcriptional start site for TERRAs. Finally, moving an X-element to a euchromatic locus changed nucleosome occupancy and positioning, demonstrating that X-element nucleosome structure is dependent on the local telomere environment.

Telomeres are specialized structures at the end of linear chromosomes that protect the genetic material from degradation and mistaken recognition as sites of damage. Telomere dysfunction has been linked to several diseases and senescence. The telomeres contain repetitive DNA sequences bound by specialized proteins. Here, we describe two such proteins, Sir2 and Reb1, which regulate the formation of nucleosomes at a repetitive sequence known as the X-element. Sir2 has very important roles in regulating the accessibility of telomeres to the cellular machinery that reads and transcribes the genetic material. Reb1 had not been previously implicated in telomere biology, but is rather known as a general regulator of transcription. We explored the effects of removing either or both of these factors on telomeric features and their relationship in regulating the structure and accessibility of the telomeres in budding yeast. We show that Sir2 and Reb1 have opposing roles in stabilizing and de-stabilizing a nucleosome at the telomeres, but that both inhibit the accumulation of a non-coding RNA molecule transcribed from the telomeres.

Overproduction of Membrane-Associated, and Integrated, Proteins Using Saccharomyces cerevisiae

Structural and functional eukaryotic membrane protein research continues to grow at an increasing rate, placing greater significance on leveraging productive protein expression pipelines to feed downstream studies. Bacterial expression systems (e.g., E. coli) are often the preferred system due to their simple growth conditions, relative simplicity in experimental workflow, low overall cost per liter of cell growth, and ease of genetic manipulation. However, overproduction success of eukaryotic membrane proteins in bacterial systems is hindered by the limited native processing ability of bacterial systems for important protein folding interactions (e.g., disulfide bonds), post-translational modifications (e.g., glycosylation), and inherent disadvantages in protein trafficking and folding machinery compared to other expression systems.In contrast, Saccharomyces cerevisiae expression systems combine positive benefits of simpler bacterial systems with those of more complex eukaryotic systems (e.g., mammalian cells). Benefits include inexpensive growth, robust DNA repair and recombination machinery, amenability to high density growths in bioreactors, efficient transformation, and robust post-translational modification machinery. These characteristics make S. cerevisiae a viable first-alternative when bacterial overproduction is insufficient. Thus, this chapter provides a framework, using methods that have proven successful in prior efforts, for overproducing membrane anchored or membrane integrated proteins in S. cerevisiae. The framework is designed to improve yields for all levels of overexpression expertise, providing optimization insights for the variety of processes involved in heterologous protein expression.

A viral histone-like protein exploits antagonism between linker histones and HMGB proteins to obstruct the cell cycle

Virus infection necessarily requires redirecting cellular resources toward viral progeny production. Adenovirus encodes the histone-like protein VII, which causes catastrophic global reorganization of host chromatin to promote virus infection. Protein VII recruits the family of high mobility group box (HMGB) proteins to chromatin along with the histone chaperone SET. As a consequence of this recruitment, we find that protein VII causes chromatin depletion of several linker histone H1 isoforms. The relationship between linker histone H1 and the functionally opposite HMGB proteins is critical for higher-order chromatin structure. However, the physiological consequences of perturbing this relationship are largely unknown. Here, we employ complementary systems in Saccharomyces cerevisiae and human cells to demonstrate that adenovirus protein VII disrupts the H1-HMGB balance to obstruct the cell cycle. We find that protein VII causes an accumulation of G2/M cells both in yeast and human systems, underscoring the high conservation of this chromatin vulnerability. In contrast, adenovirus E1A and E1B proteins are well established to override cell cycle regulation and promote transformation of human cells. Strikingly, we find that protein VII obstructs the cell cycle, even in the presence of E1A and E1B. We further show that, in a protein-VII-deleted infection, several cell cycle markers are regulated differently compared to wild-type infection, supporting our model that protein VII plays an integral role in hijacking cell cycle regulation during infection. Together, our results demonstrate that protein VII targets H1-HMGB1 antagonism to obstruct cell cycle progression, revealing an unexpected chromatin vulnerability exploited for viral benefit.

Chromosome-specific telomere lengths and the minimal functional telomere revealed by nanopore sequencing

We developed a method to tag telomeres and measure telomere length by nanopore sequencing in the yeast S. cerevisiae. Nanopore allows long-read sequencing through the telomere, through the subtelomere, and into unique chromosomal sequence, enabling assignment of telomere length to a specific chromosome end. We observed chromosome end–specific telomere lengths that were stable over 120 cell divisions. These stable chromosome-specific telomere lengths may be explained by slow clonal variation or may represent a new biological mechanism that maintains equilibrium unique to each chromosome end. We examined the role of RIF1 and TEL1 in telomere length regulation and found that TEL1 is epistatic to RIF1 at most telomeres, consistent with the literature. However, at telomeres that lack subtelomeric Y′ sequences, tel1Δ rif1Δ double mutants had a very small, but significant, increase in telomere length compared with the tel1Δ single mutant, suggesting an influence of Y′ elements on telomere length regulation. We sequenced telomeres in a telomerase-null mutant (est2Δ) and found the minimal telomere length to be ∼75 bp. In these est2Δ mutants, there were apparent telomere recombination events at individual telomeres before the generation of survivors, and these events were significantly reduced in est2Δ rad52Δ double mutants. The rate of telomere shortening in the absence of telomerase was similar across all chromosome ends at ∼5 bp per generation. This new method gives quantitative, high-resolution telomere length measurement at each individual chromosome end and suggests possible new biological mechanisms regulating telomere length.

Overdominant Mutations Restrict Adaptive Loss of Heterozygosity at Linked Loci

Loss of heterozygosity is a common mode of adaptation in asexual diploid populations. Because mitotic recombination frequently extends the full length of a chromosome arm, the selective benefit of loss of heterozygosity may be constrained by linked heterozygous mutations. In a previous laboratory evolution experiment with diploid yeast, we frequently observed homozygous mutations in the WHI2 gene on the right arm of Chromosome XV. However, when heterozygous mutations arose in the STE4 gene, another common target on Chromosome XV, loss of heterozygosity at WHI2 was not observed. Here, we show that mutations at WHI2 are partially dominant and that mutations at STE4 are overdominant. We test whether beneficial heterozygous mutations at these two loci interfere with one another by measuring loss of heterozygosity at WHI2 over 1,000 generations for ∼300 populations that differed initially only at STE4 and WHI2. We show that the presence of an overdominant mutation in STE4 reduces, but does not eliminate, loss of heterozygosity at WHI2. By sequencing 40 evolved clones, we show that populations with linked overdominant and partially dominant mutations show less parallelism at the gene level, more varied evolutionary outcomes, and increased rates of aneuploidy. Our results show that the degree of dominance and the phasing of heterozygous beneficial mutations can constrain loss of heterozygosity along a chromosome arm, and that conflicts between partially dominant and overdominant mutations can affect evolutionary outcomes.

Potassium and Sodium Salt Stress Characterization in the Yeasts Saccharomyces cerevisiae, Kluyveromyces marxianus, and Rhodotorula toruloides

Biotechnology requires efficient microbial cell factories. The budding yeast Saccharomyces cerevisiae is a vital cell factory, but more diverse cell factories are essential for the sustainable use of natural resources. Here, we benchmarked nonconventional yeasts Kluyveromyces marxianus and Rhodotorula toruloides against S. cerevisiae strains CEN.PK and W303 for their responses to potassium and sodium salt stress. We found an inverse relationship between the maximum growth rate and the median cell volume that was responsive to salt stress. The supplementation of K⁺ to CEN.PK cultures reduced Na⁺ toxicity and increased the specific growth rate 4-fold. The higher K⁺ and Na⁺ concentrations impaired ethanol and acetate metabolism in CEN.PK and acetate metabolism in W303. In R. toruloides cultures, these salt supplementations induced a trade-off between glucose utilization and cellular aggregate formation. Their combined use increased the beta-carotene yield by 60% compared with that of the reference. Neural network-based image analysis of exponential-phase cultures showed that the vacuole-to-cell volume ratio increased with increased cell volume for W303 and K. marxianus but not for CEN.PK and R. toruloides in response to salt stress. Our results provide insights into common salt stress responses in yeasts and will help design efficient bioprocesses. IMPORTANCE Characterization of microbial cell factories under industrially relevant conditions is crucial for designing efficient bioprocesses. Salt stress, typical in industrial bioprocesses, impinges upon cell volume and affects productivity. This study presents an open-source neural network-based analysis method to evaluate volumetric changes using yeast optical microscopy images. It allows quantification of cell and vacuole volumes relevant to cellular physiology. On applying salt stress in yeasts, we found that the combined use of K⁺ and Na⁺ improves the cellular fitness of Saccharomyces cerevisiae strain CEN.PK and increases the beta-carotene productivity in Rhodotorula toruloides, a commercially important antioxidant and a valuable additive in foods.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

Chromosome-level genome assembly and structural variant analysis of two laboratory yeast strains from the Peterhof Genetic Collection lineage

Thousands of yeast genomes have been sequenced with both traditional and long-read technologies, and multiple observations about modes of genome evolution for both wild and laboratory strains have been drawn from these sequences. In our study, we applied Oxford Nanopore and Illumina technologies to assemble complete genomes of two widely used members of a distinct laboratory yeast lineage, the Peterhof Genetic Collection (PGC), and investigate the structural features of these genomes including transposable element content, copy number alterations, and structural rearrangements. We identified numerous notable structural differences between genomes of PGC strains and the reference S288C strain. We discovered a substantial enrichment of mid-length insertions and deletions within repetitive coding sequences, such as in the SCH9 gene or the NUP100 gene, with possible impact of these variants on protein amyloidogenicity. High contiguity of the final assemblies allowed us to trace back the history of reciprocal unbalanced translocations between chromosomes I, VIII, IX, XI, and XVI of the PGC strains. We show that formation of hybrid alleles of the FLO genes during such chromosomal rearrangements is likely responsible for the lack of invasive growth of yeast strains. Taken together, our results highlight important features of laboratory yeast strain evolution using the power of long-read sequencing.

Exploring a Local Genetic Interaction Network Using Evolutionary Replay Experiments

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using “evolutionary replay” experiments, we identified additional mutations that have positive genetic interactions with the kel1-P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1-P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1-P344T background. We reconstructed mutations in the ancestral and kel1-P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1, we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network.

Strand-specific ChIP-seq at DNA breaks distinguishes ssDNA versus dsDNA binding and refutes single-stranded nucleosomes

In a first step of DNA double-strand break (DSB) repair by homologous recombination, DNA ends are resected such that single-stranded DNA (ssDNA) overhangs are generated. ssDNA is specifically bound by RPA and other factors, which constitutes a ssDNA-domain on damaged chromatin. The molecular organization of this ssDNA and the adjacent dsDNA domain is crucial during DSB signaling and repair. However, data regarding the presence of nucleosomes, the most basic chromatin components, in the ssDNA domain have been contradictory. Here, we use site-specific induction of DSBs and chromatin immunoprecipitation followed by strand-specific sequencing to analyze in vivo binding of key DSB repair and signaling proteins to either the ssDNA or dsDNA domain. In the case of nucleosomes, we show that recently proposed ssDNA nucleosomes are not a major, persistent species, but that nucleosome eviction and DNA end resection are intrinsically coupled. These results support a model of separated dsDNA-nucleosome and ssDNA-RPA domains during DSB repair.

Mutational analysis of the alpha subunit of eIF2B provides insights into the role of eIF2B bodies in translational control and VWM disease

Eukaryotic initiation factor 2B (eIF2B) serves as a vital control point within protein synthesis and regulates translation initiation in response to cellular stress. Mutations within eIF2B result in the fatal disease, leukoencephalopathy with vanishing white matter (VWM). Previous biochemical studies on VWM mutations have illustrated that changes in the activity of eIF2B poorly correlate with disease severity. This suggests that there may be additional characteristics of eIF2B contributing to VWM pathogenesis. Here, we investigated whether the localization of eIF2B to eIF2B bodies was integral for function and whether this localization could provide insight into the pathogenesis of VWM. We demonstrate that the regulatory subunit, eIF2Bα, is required for the assembly of eIF2B bodies in yeast and that loss of eIF2B bodies correlates with an inability of cells to regulate eIF2B activity. Mutational analysis of eIF2Bα showed that missense mutations that disrupt the regulation of eIF2B similarly disrupt the assembly of eIF2B bodies. In contrast, when eIF2Bα mutations that impact the catalytic activity of eIF2B were analyzed, eIF2B bodies were absent and instead eIF2B localized to small foci, termed microfoci. Fluorescence recovery after photobleaching analysis highlighted that within these microfoci, eIF2 shuttles more slowly indicating that formation of eIF2B bodies correlates with full eIF2B activity. When eIF2Bα VWM mutations were analyzed, a diverse impact on localization was observed, which did not seem to correlate with eIF2B activity. These findings provide key insights into how the eIF2B body assembles and suggest that the body is a fundamental part of the translational regulation via eIF2α phosphorylation.

Biochemical and genetic analysis of Ecm14, a conserved fungal pseudopeptidase

Background

Like most major enzyme families, the M14 family of metallocarboxypeptidases (MCPs) contains a number of pseudoenzymes predicted to lack enzyme activity and with poorly characterized molecular function. The genome of the yeast Saccharomyces cerevisiae encodes one member of the M14 MCP family, a pseudoenzyme named Ecm14 proposed to function in the extracellular matrix. In order to better understand the function of such pseudoenzymes, we studied the structure and function of Ecm14 in S. cerevisiae.

Results

A phylogenetic analysis of Ecm14 in fungi found it to be conserved throughout the ascomycete phylum, with a group of related pseudoenzymes found in basidiomycetes. To investigate the structure and function of this conserved protein, His6-tagged Ecm14 was overexpressed in Sf9 cells and purified. The prodomain of Ecm14 was cleaved in vivo and in vitro by endopeptidases, suggesting an activation mechanism; however, no activity was detectable using standard carboxypeptidase substrates. In order to determine the function of Ecm14 using an unbiased screen, we undertook a synthetic lethal assay. Upon screening approximately 27,000 yeast colonies, twenty-two putative synthetic lethal clones were identified. Further analysis showed many to be synthetic lethal with auxotrophic marker genes and requiring multiple mutations, suggesting that there are few, if any, single S. cerevisiae genes that present synthetic lethal interactions with ecm14Δ.

Conclusions

We show in this study that Ecm14, although lacking detectable enzyme activity, is a conserved carboxypeptidase-like protein that is secreted from cells and is processed to a mature form by the action of an endopeptidase. Our study and datasets from other recent large-scale screens suggest a role for Ecm14 in processes such as vesicle-mediated transport and aggregate invasion, a fungal process that has been selected against in modern laboratory strains of S. cerevisiae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-020-00330-w.

Mitochondrial-nuclear coadaptation revealed through mtDNA replacements in Saccharomyces cerevisiae

BACKGROUND

Mitochondrial function requires numerous genetic interactions between mitochondrial- and nuclear- encoded genes. While selection for optimal mitonuclear interactions should result in coevolution between both genomes, evidence for mitonuclear coadaptation is challenging to document. Genetic models where mitonuclear interactions can be explored are needed.

RESULTS

We systematically exchanged mtDNAs between 15 Saccharomyces cerevisiae isolates from a variety of ecological niches to create 225 unique mitochondrial-nuclear genotypes. Analysis of phenotypic profiles confirmed that environmentally-sensitive interactions between mitochondrial and nuclear genotype contributed to growth differences. Exchanges of mtDNAs between strains of the same or different clades were just as likely to demonstrate mitonuclear epistasis although epistatic effect sizes increased with genetic distances. Strains with their original mtDNAs were more fit than strains with synthetic mitonuclear combinations when grown in media that resembled isolation habitats.

CONCLUSIONS

This study shows that natural variation in mitonuclear interactions contributes to fitness landscapes. Multiple examples of coadapted mitochondrial-nuclear genotypes suggest that selection for mitonuclear interactions may play a role in helping yeasts adapt to novel environments and promote coevolution.

Effects of a high-cultivation temperature on the physiology of three different Yarrowia lipolytica strains

Despite the increasing relevance, ranging from academic research to industrial applications, only a limited number of non-conventional, oleaginous Yarrowia lipolytica strains are characterized in detail. Therefore, we analyzed three strains in regard to their metabolic and physiological properties, especially with respect to important characteristics of a production strain. By investigating different cultivation conditions and media compositions, similarities and differences between the distinct strain backgrounds could be derived. Especially sugar alcohol production, as well as an agglomeration of cells were found to be connected with growth at high temperatures. In addition, sugar alcohol production was independent of high substrate concentrations under these conditions. To investigate the genotypic basis of particular traits, including growth characteristics and metabolite concentrations, genomic analysis were performed. We found sequence variations for one third of the annotated proteins but no obvious link to all phenotypic features.

Lager-brewing yeasts in the era of modern genetics

The yeast Saccharomyces pastorianus is responsible for the annual worldwide production of almost 200 billion liters of lager-type beer. S. pastorianus is a hybrid of Saccharomyces cerevisiae and Saccharomyces eubayanus that has been studied for well over a century. Scientific interest in S. pastorianus intensified upon the discovery, in 2011, of its S. eubayanus ancestor. Moreover, advances in whole-genome sequencing and genome editing now enable deeper exploration of the complex hybrid and aneuploid genome architectures of S. pastorianus strains. These developments not only provide novel insights into the emergence and domestication of S. pastorianus but also generate new opportunities for its industrial application. This review paper combines historical, technical and socioeconomic perspectives to analyze the evolutionary origin and genetics of S. pastorianus. In addition, it provides an overview of available methods for industrial strain improvement and an outlook on future industrial application of lager-brewing yeasts. Particular attention is given to the ongoing debate on whether current S. pastorianus originates from a single or multiple hybridization events and to the potential role of genome editing in developing industrial brewing yeast strains.

The COP9 signalosome mediates the Spt23 regulated fatty acid desaturation and ergosterol biosynthesis

The COP9 signalosome (CSN) is a conserved eukaryotic complex, essential for vitality in all multicellular organisms and critical for the turnover of key cellular proteins through catalytic and non-catalytic activities. Saccharomyces cerevisiae is a powerful model organism for studying fundamental aspects of the CSN complex, since it includes a conserved enzymatic core but lacks non-catalytic activities, probably explaining its non-essentiality for life. A previous transcriptomic analysis of an S. cerevisiae strain deleted in the CSN5/RRI1 gene, encoding to the CSN catalytic subunit, revealed a downregulation of genes involved in lipid metabolism. We now show that the S. cerevisiae CSN holocomplex is essential for cellular lipid homeostasis. Defects in CSN assembly or activity lead to decreased quantities of ergosterol and unsaturated fatty acids (UFA); vacuole defects; diminished lipid droplets (LDs) size; and to accumulation of endoplasmic reticulum (ER) stress. The molecular mechanism behind these findings depends on CSN involvement in upregulating mRNA expression of SPT23. Spt23 is a novel activator of lipid desaturation and ergosterol biosynthesis. Our data reveal for the first time a functional link between the CSN holocomplex and Spt23. Moreover, CSN-dependent upregulation of SPT23 transcription is necessary for the fine-tuning of lipid homeostasis and for cellular health.

Molecular signatures of aneuploidy-driven adaptive evolution

Alteration of normal ploidy (aneuploidy) can have a number of opposing effects, such as unbalancing protein abundances and inhibiting cell growth but also accelerating genetic diversification and rapid adaptation. The interplay of these detrimental and beneficial effects remains puzzling. Here, to understand how cells develop tolerance to aneuploidy, we subject disomic (i.e. with an extra chromosome copy) strains of yeast to long-term experimental evolution under strong selection, by forcing disomy maintenance and daily population dilution. We characterize mutations, karyotype alterations and gene expression changes, and dissect the associated molecular strategies. Cells with different extra chromosomes accumulated mutations at distinct rates and displayed diverse adaptive events. They tended to evolve towards normal ploidy through chromosomal DNA loss and gene expression changes. We identify genes with recurrent mutations and altered expression in multiple lines, revealing a variant that improves growth under genotoxic stresses. These findings support rapid evolvability of disomic strains that can be used to characterize fitness effects of mutations under different stress conditions.

Aneuploidy (abnormal chromosome number) can enable rapid adaptation to stress conditions, but it also entails fitness costs from gene imbalance. Here, the authors experimentally evolve yeast while forcing maintenance of aneuploidy to identify the mechanisms that promote tolerance of aneuploidy.

Dysfunctional CAF-I reveals its role in cell cycle progression and differential regulation of gene silencing

Chromatin Assembly Factor I (CAF-I) plays a central role in the reassembly of H3/H4 histones during DNA replication. In S. cerevisiae CAF-I is not essential and its loss is associated with reduced gene silencing at telomeres and increased sensitivity to DNA damage. Two kinases, Cyclin Dependent Kinase (CDK) and Dbf4-Dependent Kinase (DDK), are known to phosphorylate the Cac1p subunit of CAF-I, but their role in the regulation of CAF-I activity is not well understood. In this study we systematically mutated the phosphorylation target sites of these kinases. We show that concomitant mutations of the CDK and DDK target sites of Cac1p lead to growth retardation and significant cell cycle defects, altered cell morphology and increased sensitivity to DNA damage. Surprisingly, some mutations also produced flocculation, a phenotype that is lost in most laboratory strains, and displayed elevated expression of FLO genes. None of these effects is observed upon the destruction of CAF-I. In contrast, the mutations that caused flocculation did not affect gene silencing at the mating type and subtelomeric loci. We conclude that dysfunctional CAF-I produces severe phenotypes, which reveal a possible role of CAF-I in the coordination of DNA replication, chromatin reassembly and cell cycle progression. Our study highlights the role of phosphorylation of Cac1p by CDK and a putative role for DDK in the transmission and re-assembly of chromatin during DNA replication.

Histone chaperones and the Rrm3p helicase regulate flocculation in S. cerevisiae

Background

Biofilm formation or flocculation is a major phenotype in wild type budding yeasts but rarely seen in laboratory yeast strains. Here, we analysed flocculation phenotypes and the expression of FLO genes in laboratory strains with various genetic backgrounds.

Results

We show that mutations in histone chaperones, the helicase RRM3 and the Histone Deacetylase HDA1 de-repress the FLO genes and partially reconstitute flocculation. We demonstrate that the loss of repression correlates to elevated expression of several FLO genes, to increased acetylation of histones at the promoter of FLO1 and to variegated expression of FLO11. We show that these effects are related to the activity of CAF-1 at the replication forks. We also demonstrate that nitrogen starvation or inhibition of histone deacetylases do not produce flocculation in W303 and BY4742 strains but do so in strains compromised for chromatin maintenance. Finally, we correlate the de-repression of FLO genes to the loss of silencing at the subtelomeric and mating type gene loci.

Conclusions

We conclude that the deregulation of chromatin maintenance and transmission is sufficient to reconstitute flocculation in laboratory yeast strains. Consequently, we propose that a gain in epigenetic silencing is a major contributing factor for the loss of flocculation phenotypes in these strains. We suggest that flocculation in yeasts provides an excellent model for addressing the challenging issue of how epigenetic mechanisms contribute to evolution.

Engineered mitochondrial production of monoterpenes in Saccharomyces cerevisiae

Monoterpene indole alkaloids (MIAs) from plants encompass a broad class of structurally complex and medicinally valuable natural products. MIAs are biologically derived from the universal precursor strictosidine. Although the strictosidine biosynthetic pathway has been identified and reconstituted, extensive work is required to optimize production of strictosidine and its precursors in yeast. In this study, we engineered a fully integrated and plasmid-free yeast strain with enhanced production of the monoterpene precursor geraniol. The geraniol biosynthetic pathway was targeted to the mitochondria to protect the GPP pool from consumption by the cytosolic ergosterol pathway. The mitochondrial geraniol producer showed a 6-fold increase in geraniol production compared to cytosolic producing strains. We further engineered the monoterpene-producing strain to synthesize the next intermediates in the strictosidine pathway: 8-hydroxygeraniol and nepetalactol. Integration of geraniol hydroxylase (G8H) from Catharanthus roseus led to essentially quantitative conversion of geraniol to 8-hydroxygeraniol at a titer of 227 mg/L in a fed-batch fermentation. Further introduction of geraniol oxidoreductase (GOR) and iridoid synthase (ISY) from C. roseus and tuning of the relative expression levels resulted in the first de novo nepetalactol production. The strategies developed in this work can facilitate future strain engineering for yeast production of later intermediates in the strictosidine biosynthetic pathway.

Integrative Meta-Assembly Pipeline (IMAP): Chromosome-level genome assembler combining multiple de novo assemblies

BACKGROUND

Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data.

METHODS AND RESULTS

In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1.

CONCLUSIONS

We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at http://github.com/jkimlab/IMAP including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest.

SSD1 suppresses phenotypes induced by the lack of Elongator-dependent tRNA modifications

The Elongator complex promotes formation of 5-methoxycarbonylmethyl (mcm5) and 5-carbamoylmethyl (ncm5) side-chains on uridines at the wobble position of cytosolic eukaryotic tRNAs. In all eukaryotic organisms tested to date, the inactivation of Elongator not only leads to the lack of mcm5/ncm5 groups in tRNAs, but also a wide variety of additional phenotypes. Although the phenotypes are most likely caused by a translational defect induced by reduced functionality of the hypomodified tRNAs, the mechanism(s) underlying individual phenotypes are poorly understood. In this study, we show that the genetic background modulates the phenotypes induced by the lack of mcm5/ncm5 groups in Saccharomyces cerevisiae. We show that the stress-induced growth defects of Elongator mutants are stronger in the W303 than in the closely related S288C genetic background and that the phenotypic differences are caused by the known polymorphism at the locus for the mRNA binding protein Ssd1. Moreover, the mutant ssd1 allele found in W303 cells is required for the reported histone H3 acetylation and telomeric gene silencing defects of Elongator mutants. The difference at the SSD1 locus also partially explains why the simultaneous lack of mcm5 and 2-thio groups at wobble uridines is lethal in the W303 but not in the S288C background. Collectively, our results demonstrate that the SSD1 locus modulates phenotypes induced by the lack of Elongator-dependent tRNA modifications.

Yeast Expression Systems: Overview and Recent Advances

Yeasts are outstanding hosts for the production of functional recombinant proteins with industrial or medical applications. Great attention has been emerged on yeast due to the inherent advantages and new developments in this host cell. For the production of each specific product, the most appropriate expression system should be identified and optimized both on the genetic and fermentation levels, considering the features of the host, vector and expression strategies. Currently, several new systems are commercially available; some of them are private and need licensing. The potential for secretory expression of heterologous proteins in yeast proposed this system as a candidate for the production of complex eukaryotic proteins. The common yeast expression hosts used for recombinant proteins' expression include Saccharomyces cerevisiae, Pichia pastoris, Hansenula polymorpha, Yarrowia lipolytica, Arxula adeninivorans, Kluyveromyces lactis, and Schizosaccharomyces pombe. This review is dedicated to discuss on significant characteristics of the most common methylotrophic and non-methylotrophic yeast expression systems with an emphasis on their advantages and new developments.

The ATP-dependent chromatin remodelling enzyme Uls1 prevents Topoisomerase II poisoning

Topoisomerase II (Top2) is an essential enzyme that decatenates DNA via a transient Top2-DNA covalent intermediate. This intermediate can be stabilized by a class of drugs termed Top2 poisons, resulting in massive DNA damage. Thus, Top2 activity is a double-edged sword that needs to be carefully controlled to maintain genome stability. We show that Uls1, an adenosine triphosphate (ATP)-dependent chromatin remodelling (Snf2) enzyme, can alter Top2 chromatin binding and prevent Top2 poisoning in yeast. Deletion mutants of ULS1 are hypersensitive to the Top2 poison acriflavine (ACF), activating the DNA damage checkpoint. We map Uls1's Top2 interaction domain and show that this, together with its ATPase activity, is essential for Uls1 function. By performing ChIP-seq, we show that ACF leads to a general increase in Top2 binding across the genome. We map Uls1 binding sites and identify tRNA genes as key regions where Uls1 associates after ACF treatment. Importantly, the presence of Uls1 at these sites prevents ACF-dependent Top2 accumulation. Our data reveal the effect of Top2 poisons on the global Top2 binding landscape and highlights the role of Uls1 in antagonizing Top2 function. Remodelling Top2 binding is thus an important new means by which Snf2 enzymes promote genome stability.

Detecting genetic interactions using parallel evolution in experimental populations

Eukaryotic genomes contain thousands of genes organized into complex and interconnected genetic interaction networks. Most of our understanding of how genetic variation affects these networks comes from quantitative-trait loci mapping and from the systematic analysis of double-deletion (or knockdown) mutants, primarily in the yeast Saccharomyces cerevisiae. Evolve and re-sequence experiments are an alternative approach for identifying novel functional variants and genetic interactions, particularly between non-loss-of-function mutations. These experiments leverage natural selection to obtain genotypes with functionally important variants and positive genetic interactions. However, no systematic methods for detecting genetic interactions in these data are yet available. Here, we introduce a computational method based on the idea that variants in genes that interact will co-occur in evolved genotypes more often than expected by chance. We apply this method to a previously published yeast experimental evolution dataset. We find that genetic targets of selection are distributed non-uniformly among evolved genotypes, indicating that genetic interactions had a significant effect on evolutionary trajectories. We identify individual gene pairs with a statistically significant genetic interaction score. The strongest interaction is between genes TRK1 and PHO84, genes that have not been reported to interact in previous systematic studies. Our work demonstrates that leveraging parallelism in experimental evolution is useful for identifying genetic interactions that have escaped detection by other methods. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Exploring the xylose paradox in Saccharomyces cerevisiae through in vivo sugar signalomics of targeted deletants

BACKGROUND

There have been many successful strategies to implement xylose metabolism in Saccharomyces cerevisiae, but no effort has so far enabled xylose utilization at rates comparable to that of glucose (the preferred sugar of this yeast). Many studies have pointed towards the engineered yeast not sensing that xylose is a fermentable carbon source despite growing and fermenting on it, which is paradoxical. We have previously used fluorescent biosensor strains to in vivo monitor the sugar signalome in yeast engineered with xylose reductase and xylitol dehydrogenase (XR/XDH) and have established that S. cerevisiae senses high concentrations of xylose with the same signal as low concentration of glucose, which may explain the poor utilization.

RESULTS

In the present study, we evaluated the effects of three deletions (ira2∆, isu1∆ and hog1∆) that have recently been shown to display epistatic effects on a xylose isomerase (XI) strain. Through aerobic and anaerobic characterization, we showed that the proposed effects in XI strains were for the most part also applicable in the XR/XDH background. The ira2∆isu1∆ double deletion led to strains with the highest specific xylose consumption- and ethanol production rates but also the lowest biomass titre. The signalling response revealed that ira2∆isu1∆ changed the low glucose-signal in the background strain to a simultaneous signalling of high and low glucose, suggesting that engineering of the signalome can improve xylose utilization.

CONCLUSIONS

The study was able to correlate the previously proposed beneficial effects of ira2∆, isu1∆ and hog1∆ on S. cerevisiae xylose uptake, with a change in the sugar signalome. This is in line with our previous hypothesis that the key to resolve the xylose paradox lies in the sugar sensing and signalling networks. These results indicate that the future engineering targets for improved xylose utilization should probably be sought not in the metabolic networks, but in the signalling ones.

The Genetic Basis of Mutation Rate Variation in Yeast

Mutations are the root source of genetic variation and underlie the process of evolution. Although the rates at which mutations occur vary considerably between species, little is known about differences within species, or the genetic and molecular basis of these differences. Here, we leveraged the power of the yeast Saccharomyces cerevisiae as a model system to uncover natural genetic variants that underlie variation in mutation rate. We developed a high-throughput fluctuation assay and used it to quantify mutation rates in seven natural yeast isolates and in 1040 segregant progeny from a cross between BY, a laboratory strain, and RM, a wine strain. We observed that mutation rate varies among yeast strains and is heritable (H² = 0.49). We performed linkage mapping in the segregants and identified four quantitative trait loci underlying mutation rate variation in the cross. We fine-mapped two quantitative trait loci to the underlying causal genes, RAD5 and MKT1, that contribute to mutation rate variation. These genes also underlie sensitivity to the DNA-damaging agents 4NQO and MMS, suggesting a connection between spontaneous mutation rate and mutagen sensitivity.

Altered access to beneficial mutations slows adaptation and biases fixed mutations in diploids

Ploidy varies considerably in nature. However, our understanding of the impact of ploidy on adaptation is incomplete. Many microbial evolution experiments characterize adaptation in haploid organisms, but few focus on diploid organisms. Here, we perform a 4,000-generation evolution experiment using diploid strains of the yeast Saccharomyces cerevisiae. We show that the rate of adaptation and spectrum of beneficial mutations are influenced by ploidy. Haldane's sieve effectively alters access to recessive beneficial mutations in diploid populations, leading to a slower rate of adaptation and a spectrum of beneficial mutations that is shifted towards dominant mutations. Genomic position also has an important role, as the prevalence of homozygous mutations is largely dependent on their proximity to a recombination hotspot. Our results demonstrate key aspects of diploid adaptation that have previously been understudied and provide support for several proposed theories.

Nanopore sequencing enables near-complete de novo assembly of Saccharomyces cerevisiae reference strain CEN.PK113-7D

The haploid Saccharomyces cerevisiae strain CEN.PK113-7D is a popular model system for metabolic engineering and systems biology research. Current genome assemblies are based on short-read sequencing data scaffolded based on homology to strain S288C. However, these assemblies contain large sequence gaps, particularly in subtelomeric regions, and the assumption of perfect homology to S288C for scaffolding introduces bias. In this study, we obtained a near-complete genome assembly of CEN.PK113-7D using only Oxford Nanopore Technology's MinION sequencing platform. Fifteen of the 16 chromosomes, the mitochondrial genome and the 2-μm plasmid are assembled in single contigs and all but one chromosome starts or ends in a telomere repeat. This improved genome assembly contains 770 Kbp of added sequence containing 248 gene annotations in comparison to the previous assembly of CEN.PK113-7D. Many of these genes encode functions determining fitness in specific growth conditions and are therefore highly relevant for various industrial applications. Furthermore, we discovered a translocation between chromosomes III and VIII that caused misidentification of a MAL locus in the previous CEN.PK113-7D assembly. This study demonstrates the power of long-read sequencing by providing a high-quality reference assembly and annotation of CEN.PK113-7D and places a caveat on assumed genome stability of microorganisms.