Whole-genome sequencing of sake yeast Saccharomyces cerevisiae Kyokai no. 7

Genome-wide analysis of intraspecific transposon diversity in yeast

Background

In the model organism Saccharomyces cerevisiae, the transposable elements (TEs) consist of LTR (Long Terminal Repeat) retrotransposons called Ty elements belonging to five families, Ty1 to Ty5. They take the form of either full-length coding elements or non-coding solo-LTRs corresponding to remnants of former transposition events. Although the biological features of Ty elements have been studied in detail in S. cerevisiae and the Ty content of the reference strain (S288c) was accurately annotated, the Ty-related intra-specific diversity has not been closely investigated so far.

Results

In this study, we investigated the Ty contents of 41 available genomes of isolated S. cerevisiae strains of diverse geographical and ecological origins. The strains were compared in terms of the number of Ty copies, the content of the potential transpositionally active elements and the genomic insertion maps. The strain repertoires were also investigated in the closely related Ty1 and Ty2 families and subfamilies.

Conclusions

This is the first genome-wide analysis of the diversity associated to the Ty elements, carried out for a large set of S. cerevisiae strains. The results of the present analyses suggest that the current Ty-related polymorphism has resulted from multiple causes such as differences between strains, between Ty families and over time, in the recent transpositional activity of Ty elements. Some new Ty1 variants were also identified, and we have established that Ty1 variants have different patterns of distribution among strains, which further contributes to the strain diversity.

Comparative polygenic analysis of maximal ethanol accumulation capacity and tolerance to high ethanol levels of cell proliferation in yeast

The yeast Saccharomyces cerevisiae is able to accumulate ≥17% ethanol (v/v) by fermentation in the absence of cell proliferation. The genetic basis of this unique capacity is unknown. Up to now, all research has focused on tolerance of yeast cell proliferation to high ethanol levels. Comparison of maximal ethanol accumulation capacity and ethanol tolerance of cell proliferation in 68 yeast strains showed a poor correlation, but higher ethanol tolerance of cell proliferation clearly increased the likelihood of superior maximal ethanol accumulation capacity. We have applied pooled-segregant whole-genome sequence analysis to identify the polygenic basis of these two complex traits using segregants from a cross of a haploid derivative of the sake strain CBS1585 and the lab strain BY. From a total of 301 segregants, 22 superior segregants accumulating ≥17% ethanol in small-scale fermentations and 32 superior segregants growing in the presence of 18% ethanol, were separately pooled and sequenced. Plotting SNP variant frequency against chromosomal position revealed eleven and eight Quantitative Trait Loci (QTLs) for the two traits, respectively, and showed that the genetic basis of the two traits is partially different. Fine-mapping and Reciprocal Hemizygosity Analysis identified ADE1, URA3, and KIN3, encoding a protein kinase involved in DNA damage repair, as specific causative genes for maximal ethanol accumulation capacity. These genes, as well as the previously identified MKT1 gene, were not linked in this genetic background to tolerance of cell proliferation to high ethanol levels. The superior KIN3 allele contained two SNPs, which are absent in all yeast strains sequenced up to now. This work provides the first insight in the genetic basis of maximal ethanol accumulation capacity in yeast and reveals for the first time the importance of DNA damage repair in yeast ethanol tolerance.

The yeast Saccharomyces cerevisiae is unique in being the most ethanol tolerant organism known. This property lies at the basis of its ecological competitiveness in sugar-rich ecological niches and its use for the production of alcoholic beverages and bioethanol, both of which involve accumulation of high levels of ethanol. Up to now, all research on yeast ethanol tolerance has focused on tolerance of cell proliferation to high ethanol levels. However, the most ecologically and industrially relevant aspect is the capacity of fermenting yeast cells to accumulate high ethanol levels in the absence of cell proliferation. Using QTL mapping by pooled-segregant whole-genome sequence analysis, we show that maximal ethanol accumulation capacity and tolerance of cell proliferation to high ethanol levels have a partially different genetic basis. We identified three specific genes responsible for high ethanol accumulation capacity, of which one gene encodes a protein kinase involved in DNA damage repair. Our work provides the first insight in the genetic basis of maximal ethanol accumulation capacity, shows that it involves different genetic elements compared to tolerance of cell proliferation to high ethanol levels, and reveals for the first time the importance of DNA damage repair in ethanol tolerance.

Genetic Basis of Ammonium Toxicity Resistance in a Sake Strain of Yeast: A Mendelian Case

High concentrations of ammonium at physiological concentrations of potassium are toxic for the standard laboratory strain of Saccharomyces cerevisiae. In the original description of this metabolic phenotype, we focused on the standard laboratory strains of Saccharomyces. In this study, we screened a large collection of S. cerevisiae natural isolates and identified one strain that is resistant to high concentrations of ammonium. This strain, K12, was isolated in sake breweries. When the K12 strain was crossed to the standard laboratory strain (FY4), the resulting tetrads displayed 2:2 segregation of the resistance phenotype, suggesting a single gene trait. Using a bulk segregant analysis strategy, we mapped this trait to a 150-kb region on chromosome X containing the TRK1 gene. This gene encodes a transporter required for high-affinity potassium transport in S. cerevisiae. Data from reciprocal hemizygosity experiments with TRK1 deletion strains in K12 and BY backgrounds, as well as analysis of the deletion of this gene in the K12 strain, demonstrate that the K12 allele of TRK1 is responsible for ammonium toxicity resistance. Furthermore, we determined the minimal amount of potassium required for both the K12 and laboratory strain needed for growth. These results demonstrate that the gene encoded by the K12 allele of TRK1 has a greater affinity for potassium than the standard allele of TRK1 found in Saccharomyces strains. We hypothesize that this greater-affinity allele of the potassium transporter reduces the flux of ammonium into the yeast cells under conditions of ammonium toxicity. These findings further refine our understanding of ammonium toxicity in yeast and provide an example of using natural variation to understand cellular processes.

The new modern era of yeast genomics: community sequencing and the resulting annotation of multiple Saccharomyces cerevisiae strains at the Saccharomyces Genome Database

The first completed eukaryotic genome sequence was that of the yeast Saccharomyces cerevisiae, and the Saccharomyces Genome Database (SGD; http://www.yeastgenome.org/) is the original model organism database. SGD remains the authoritative community resource for the S. cerevisiae reference genome sequence and its annotation, and continues to provide comprehensive biological information correlated with S. cerevisiae genes and their products. A diverse set of yeast strains have been sequenced to explore commercial and laboratory applications, and a brief history of those strains is provided. The publication of these new genomes has motivated the creation of new tools, and SGD will annotate and provide comparative analyses of these sequences, correlating changes with variations in strain phenotypes and protein function. We are entering a new era at SGD, as we incorporate these new sequences and make them accessible to the scientific community, all in an effort to continue in our mission of educating researchers and facilitating discovery.

Database URL:http://www.yeastgenome.org/

Genetic instability of constitutive acid phosphatase in shochu and sake yeast

Genetic instability of constitutive acid phosphatase (cAPase) activity was observed in a shochu brewer's yeast strain (Ko), which consistently produced 0.3-1% progeny without cAPase when it had been subcultured for a long period of time in barley shochu mash or in conventional complete medium. Genetic analysis showed that the cAPase-negative phenotype was associated with a single mutation in the PHO3 gene and that the Ko strain had heteroallelic PHO3/pho3 genes, while the PHO3⁻ mutants had the homoallelic pho3/pho3 defect. Some sake yeast strains that are cAPase negative, such as K6, K7 and K9, also had the same homoallelic defect, whereas another sake yeast strain K3, with heteroallelic PHO3/pho3 genes, displayed similar genetic instability of cAPase activity. In all cases, the pho3-defective genes were generated by deletion of an approximately 1.9 kb region between the PHO5-PHO3 tandem genes on chromosome II, resulting in chimeric PHO5/3 fusion genes with different fusion points. By integrating a lys2 marker, which is linked with the pho3 allele on the arm of chromosome II in the Ko strain, we demonstrated that the pho3/pho3 defect originated either from a loss of heterozygosity at the heteroallelic PHO3/pho3 locus or from a looping out of the PHO3 region. Although fermentation experiments have not yet indicated any correlation between cAPase activity and alcohol production, the PHO3⁻ mutation itself could prove to be a useful selective marker for yeast strains carrying a number of advantageous mutations for fermentation and which display phenotypic diversity and stability.

The genetic architecture of biofilm formation in a clinical isolate of Saccharomyces cerevisiae

Biofilms are microbial communities that form on surfaces. They are the primary form of microbial growth in nature and can have detrimental impacts on human health. Some strains of the budding yeast Saccharomyces cerevisiae form colony biofilms, and there is substantial variation in colony architecture between biofilm-forming strains. To identify the genetic basis of biofilm variation, we developed a novel version of quantitative trait locus mapping, which leverages cryptic variation in a clinical isolate of S. cerevisiae. We mapped 13 loci linked to heterogeneity in biofilm architecture and identified the gene most closely associated with each locus. Of these candidate genes, six are members of the cyclic AMP-protein kinase A pathway, an evolutionarily conserved cell signaling network. Principal among these is CYR1, which encodes the enzyme that catalyzes production of cAMP. Through a combination of gene expression measurements, cell signaling assays, and gene overexpression, we determined the functional effects of allelic variation at CYR1. We found that increased pathway activity resulting from protein coding and expression variation of CYR1 enhances the formation of colony biofilms. Four other candidate genes encode kinases and transcription factors that are targets of this pathway. The protein products of several of these genes together regulate expression of the sixth candidate, FLO11, which encodes a cell adhesion protein. Our results indicate that epistatic interactions between alleles with both positive and negative effects on cyclic AMP-protein kinase A signaling underlie much of the architectural variation we observe in colony biofilms. They are also among the first to demonstrate genetic variation acting at multiple levels of an integrated signaling and regulatory network. Based on these results, we propose a mechanistic model that relates genetic variation to gene network function and phenotypic outcomes.

QTL mapping of the production of wine aroma compounds by yeast

Background

Wine aroma results from the combination of numerous volatile compounds, some produced by yeast and others produced in the grapes and further metabolized by yeast. However, little is known about the consequences of the genetic variation of yeast on the production of these volatile metabolites, or on the metabolic pathways involved in the metabolism of grape compounds. As a tool to decipher how wine aroma develops, we analyzed, under two experimental conditions, the production of 44 compounds by a population of 30 segregants from a cross between a laboratory strain and an industrial strain genotyped at high density.

Results

We detected eight genomic regions explaining the diversity concerning 15 compounds, some produced de novo by yeast, such as nerolidol, ethyl esters and phenyl ethanol, and others derived from grape compounds such as citronellol, and cis-rose oxide. In three of these eight regions, we identified genes involved in the phenotype. Hemizygote comparison allowed the attribution of differences in the production of nerolidol and 2-phenyl ethanol to the PDR8 and ABZ1 genes, respectively. Deletion of a PLB2 gene confirmed its involvement in the production of ethyl esters. A comparison of allelic variants of PDR8 and ABZ1 in a set of available sequences revealed that both genes present a higher than expected number of non-synonymous mutations indicating possible balancing selection.

Conclusions

This study illustrates the value of QTL analysis for the analysis of metabolic traits, and in particular the production of wine aromas. It also identifies the particular role of the PDR8 gene in the production of farnesyldiphosphate derivatives, of ABZ1 in the production of numerous compounds and of PLB2 in ethyl ester synthesis. This work also provides a basis for elucidating the metabolism of various grape compounds, such as citronellol and cis-rose oxide.

Genome sequencing and genetic breeding of a bioethanol Saccharomyces cerevisiae strain YJS329

Background

Environmental stresses and inhibitors encountered by Saccharomyces cerevisiae strains are the main limiting factors in bioethanol fermentation. Strains with different genetic backgrounds usually show diverse stress tolerance responses. An understanding of the mechanisms underlying these phenotypic diversities within S. cerevisiae populations could guide the construction of strains with desired traits.

Results

We explored the genetic characteristics of the bioethanol S. cerevisiae strain YJS329 and elucidated how genetic variations in its genome were correlated with specified traits compared to similar traits in the S288c-derived strain, BYZ1. Karyotypic electrophoresis combined with array-comparative genomic hybridization indicated that YJS329 was a diploid strain with a relatively constant genome as a result of the fewer Ty elements and lack of structural polymorphisms between homologous chromosomes that it contained. By comparing the sequence with the S288c genome, a total of 64,998 SNPs, 7,093 indels and 11 unique genes were identified in the genome of YJS329-derived haploid strain YJSH1 through whole-genome sequencing. Transcription comparison using RNA-Seq identified which of the differentially expressed genes were the main contributors to the phenotypic differences between YJS329 and BYZ1. By combining the results obtained from the genome sequences and the transcriptions, we predicted how the SNPs, indels and chromosomal copy number variations may affect the mRNA expression profiles and phenotypes of the yeast strains. Furthermore, some genetic breeding strategies to improve the adaptabilities of YJS329 were designed and experimentally verified.

Conclusions

Through comparative functional genomic analysis, we have provided some insights into the mechanisms underlying the specific traits of the bioenthanol strain YJS329. The work reported here has not only enriched the available genetic resources of yeast but has also indicated how functional genomic studies can be used to improve genetic breeding in yeast.

A complete sequence of Saccharomyces paradoxus mitochondrial genome that restores the respiration in S. cerevisiae

We determined the complete sequence of 71 355-bp-long mitochondrial genome from Saccharomyces paradoxus entirely by direct sequencing of purified mitochondrial DNA (mtDNA). This mtDNA possesses the same features as its close relative Saccharomyces cerevisiae - A + T content 85.9%, set of genes coding for the three components of cytochrome oxidase, cytochrome b, three subunits of ATPase, both ribosomal subunits, gene for ribosomal protein, rnpB gene, tRNA package (24) and yeast genetic code. Genes are interrupted by nine group I and group II introns, two of which are in positions unknown in S. cerevisiae, but recognized in Saccharomyces pastorianus. The gene products are related to S. cerevisiae, and the identity of amino acid residues varies from 100% for cox2 to 83% for rps3. The remarkable differences from S. cerevisiae are (1) different gene order (translocation of trnF-trnT1-trnV-cox3-trnfM-rnpb-trnP and transposition of trnW-rns), (2) occurrence of two unusual GI introns, (3) eight active ori elements, and (4) reduced number of GC clusters and divergent intergenic spacers. Despite these facts, the sequenced S. paradoxus mtDNA introduced to S. cerevisiae was able to support the respiratory function to the same extent as the original mtDNAs.

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the ~12-Mb genome of CAT-1, when compared with the reference S228c genome, contains ~36,000 homozygous and ~30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.

Turbidostat culture of Saccharomyces cerevisiae W303-1A under selective pressure elicited by ethanol selects for mutations in SSD1 and UTH1

We investigated the genetic causes of ethanol tolerance by characterizing mutations selected in Saccharomyces cerevisiae W303-1A under the selective pressure of ethanol. W303-1A was subjected to three rounds of turbidostat, in a medium supplemented with increasing amounts of ethanol. By the end of selection, the growth rate of the culture has increased from 0.029 to 0.32 h(-1) . Unlike the progenitor strain, all yeast cells isolated from this population were able to form colonies on medium supplemented with 7% ethanol within 6 days, our definition of ethanol tolerance. Several clones selected from all three stages of selection were able to form dense colonies within 2 days on solid medium supplemented with 9% ethanol. We sequenced the whole genomes of six clones and identified mutations responsible for ethanol tolerance. Thirteen additional clones were tested for the presence of similar mutations. In 15 of 19 tolerant clones, the stop codon in ssd1-d was replaced with an amino acid-encoding codon. Three other clones contained one of two mutations in UTH1, and one clone did not contain mutations in either SSD1 or UTH1. We showed that the mutations in SSD1 and UTH1 increased tolerance of the cell wall to zymolyase and conclude that stability of the cell wall is a major factor in increased tolerance to ethanol.

De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113-7D, a model for modern industrial biotechnology

Saccharomyces cerevisiae CEN.PK 113-7D is widely used for metabolic engineering and systems biology research in industry and academia. We sequenced, assembled, annotated and analyzed its genome. Single-nucleotide variations (SNV), insertions/deletions (indels) and differences in genome organization compared to the reference strain S. cerevisiae S288C were analyzed. In addition to a few large deletions and duplications, nearly 3000 indels were identified in the CEN.PK113-7D genome relative to S288C. These differences were overrepresented in genes whose functions are related to transcriptional regulation and chromatin remodelling. Some of these variations were caused by unstable tandem repeats, suggesting an innate evolvability of the corresponding genes. Besides a previously characterized mutation in adenylate cyclase, the CEN.PK113-7D genome sequence revealed a significant enrichment of non-synonymous mutations in genes encoding for components of the cAMP signalling pathway. Some phenotypic characteristics of the CEN.PK113-7D strains were explained by the presence of additional specific metabolic genes relative to S288C. In particular, the presence of the BIO1 and BIO6 genes correlated with a biotin prototrophy of CEN.PK113-7D. Furthermore, the copy number, chromosomal location and sequences of the MAL loci were resolved. The assembled sequence reveals that CEN.PK113-7D has a mosaic genome that combines characteristics of laboratory strains and wild-industrial strains.

A loss-of-function mutation in the PAS kinase Rim15p is related to defective quiescence entry and high fermentation rates of Saccharomyces cerevisiae sake yeast strains

Sake yeast cells have defective entry into the quiescent state, allowing them to sustain high fermentation rates. To reveal the underlying mechanism, we investigated the PAS kinase Rim15p, which orchestrates initiation of the quiescence program in Saccharomyces cerevisiae. We found that Rim15p is truncated at the carboxyl terminus in modern sake yeast strains as a result of a frameshift mutation. Introduction of this mutation or deletion of the full-length RIM15 gene in a laboratory strain led to a defective stress response, decreased synthesis of the storage carbohydrates trehalose and glycogen, and impaired G(1) arrest, which together closely resemble the characteristic phenotypes of sake yeast. Notably, expression of a functional RIM15 gene in a modern sake strain suppressed all of these phenotypes, demonstrating that dysfunction of Rim15p prevents sake yeast cells from entering quiescence. Moreover, loss of Rim15p or its downstream targets Igo1p and Igo2p remarkably improved the fermentation rate in a laboratory strain. This finding verified that Rim15p-mediated entry into quiescence plays pivotal roles in the inhibition of ethanol fermentation. Taken together, our results suggest that the loss-of-function mutation in the RIM15 gene may be the key genetic determinant of the increased ethanol production rates in modern sake yeast strains.

Association of constitutive hyperphosphorylation of Hsf1p with a defective ethanol stress response in Saccharomyces cerevisiae sake yeast strains

Modern sake yeast strains, which produce high concentrations of ethanol, are unexpectedly sensitive to environmental stress during sake brewing. To reveal the underlying mechanism, we investigated a well-characterized yeast stress response mediated by a heat shock element (HSE) and heat shock transcription factor Hsf1p in Saccharomyces cerevisiae sake yeast. The HSE-lacZ activity of sake yeast during sake fermentation and under acute ethanol stress was severely impaired compared to that of laboratory yeast. Moreover, the Hsf1p of modern sake yeast was highly and constitutively hyperphosphorylated, irrespective of the extracellular stress. Since HSF1 allele replacement did not significantly affect the HSE-mediated ethanol stress response or Hsf1p phosphorylation patterns in either sake or laboratory yeast, the regulatory machinery of Hsf1p is presumed to function differently between these types of yeast. To identify phosphatases whose loss affected the control of Hsf1p, we screened a series of phosphatase gene deletion mutants in a laboratory strain background. Among the 29 mutants, a Δppt1 mutant exhibited constitutive hyperphosphorylation of Hsf1p, similarly to the modern sake yeast strains, which lack the entire PPT1 gene locus. We confirmed that the expression of laboratory yeast-derived functional PPT1 recovered the HSE-mediated stress response of sake yeast. In addition, deletion of PPT1 in laboratory yeast resulted in enhanced fermentation ability. Taken together, these data demonstrate that hyperphosphorylation of Hsf1p caused by loss of the PPT1 gene at least partly accounts for the defective stress response and high ethanol productivity of modern sake yeast strains.