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.

Improving Industrially Relevant Phenotypic Traits by Engineering Chromosome Copy Number in Saccharomyces pastorianus

The lager-brewing yeast Saccharomyces pastorianus is a hybrid between S. cerevisiae and S. eubayanus with an exceptional degree of aneuploidy. While chromosome copy number variation (CCNV) is present in many industrial Saccharomyces strains and has been linked to various industrially-relevant traits, its impact on the brewing performance of S. pastorianus remains elusive. Here we attempt to delete single copies of chromosomes which are relevant for the production of off-flavor compound diacetyl by centromere silencing. However, the engineered strains display CNV of multiple non-targeted chromosomes. We attribute this unintended CCNV to inherent instability and to a mutagenic effect of electroporation and of centromere-silencing. Regardless, the resulting strains displayed large phenotypic diversity. By growing centromere-silenced cells in repeated sequential batches in medium containing 10% ethanol, mutants with increased ethanol tolerance were obtained. By using CCNV mutagenesis by exposure to the mitotic inhibitor MBC, selection in the same set-up yielded even more tolerant mutants that would not classify as genetically modified organisms. These results show that CCNV of alloaneuploid S. pastorianus genomes is highly unstable, and that CCNV mutagenesis can generate broad diversity. Coupled to effective selection or screening, CCNV mutagenesis presents a potent tool for strain improvement.

Chromosome level assembly and comparative genome analysis confirm lager-brewing yeasts originated from a single hybridization

BACKGROUND

The lager brewing yeast, S. pastorianus, is a hybrid between S. cerevisiae and S. eubayanus with extensive chromosome aneuploidy. S. pastorianus is subdivided into Group 1 and Group 2 strains, where Group 2 strains have higher copy number and a larger degree of heterozygosity for S. cerevisiae chromosomes. As a result, Group 2 strains were hypothesized to have emerged from a hybridization event distinct from Group 1 strains. Current genome assemblies of S. pastorianus strains are incomplete and highly fragmented, limiting our ability to investigate their evolutionary history.

RESULTS

To fill this gap, we generated a chromosome-level genome assembly of the S. pastorianus strain CBS 1483 from Oxford Nanopore MinION DNA sequencing data and analysed the newly assembled subtelomeric regions and chromosome heterozygosity. To analyse the evolutionary history of S. pastorianus strains, we developed Alpaca: a method to compute sequence similarity between genomes without assuming linear evolution. Alpaca revealed high similarities between the S. cerevisiae subgenomes of Group 1 and 2 strains, and marked differences from sequenced S. cerevisiae strains.

CONCLUSIONS

Our findings suggest that Group 1 and Group 2 strains originated from a single hybridization involving a heterozygous S. cerevisiae strain, followed by different evolutionary trajectories. The clear differences between both groups may originate from a severe population bottleneck caused by the isolation of the first pure cultures. Alpaca provides a computationally inexpensive method to analyse evolutionary relationships while considering non-linear evolution such as horizontal gene transfer and sexual reproduction, providing a complementary viewpoint beyond traditional phylogenetic approaches.

Himalayan Saccharomyces eubayanus Genome Sequences Reveal Genetic Markers Explaining Heterotic Maltotriose Consumption by Saccharomyces pastorianus Hybrids

Saccharomyces pastorianus strains are hybrids of Saccharomyces cerevisiae and Saccharomyces eubayanus that have been domesticated for centuries in lager beer brewing environments. As sequences and structures of S. pastorianus genomes are being resolved, molecular mechanisms and evolutionary origins of several industrially relevant phenotypes remain unknown. This study investigates how maltotriose metabolism, a key feature in brewing, may have arisen in early S. eubayanus × S. cerevisiae hybrids. To address this question, we generated a nearly complete genome assembly of Himalayan S. eubayanus strains of the Holarctic subclade. This group of strains has been proposed to be the S. eubayanus subgenome origin of current S. pastorianus strains. The Himalayan S. eubayanus genomes harbored several copies of an S. eubayanusAGT1 (SeAGT1) α-oligoglucoside transporter gene with high sequence identity to genes encountered in S. pastorianus Although Himalayan S. eubayanus strains cannot grow on maltose and maltotriose, their maltose-hydrolase and SeMALT1 and SeAGT1 maltose transporter genes complemented the corresponding null mutants of S. cerevisiae Expression, in Himalayan S. eubayanus of a functional S. cerevisiae maltose metabolism regulator gene (MALx3) enabled growth on oligoglucosides. The hypothesis that the maltotriose-positive phenotype in S. pastorianus is a result of heterosis was experimentally tested by constructing an S. cerevisiae × S. eubayanus laboratory hybrid with a complement of maltose metabolism genes that resembles that of current S. pastorianus strains. The ability of this hybrid to consume maltotriose in brewer's wort demonstrated regulatory cross talk between subgenomes and thereby validated this hypothesis. These results support experimentally the new postulated hypothesis on the evolutionary origin of an essential phenotype of lager brewing strains and valuable knowledge for industrial exploitation of laboratory-made S. pastorianus-like hybrids.IMPORTANCES. pastorianus, an S. cerevisiae × S. eubayanus hybrid, is used for production of lager beer, the most produced alcoholic beverage worldwide. It emerged by spontaneous hybridization and colonized early lager brewing processes. Despite accumulation and analysis of genome sequencing data of S. pastorianus parental genomes, the genetic blueprint of industrially relevant phenotypes remains unresolved. Assimilation of maltotriose, an abundant sugar in wort, has been postulated to be inherited from the S. cerevisiae parent. Here, we demonstrate that although Asian S. eubayanus isolates harbor a functional maltotriose transporter SeAGT1 gene, they are unable to grow on α-oligoglucosides, but expression of S. cerevisiae regulator MAL13 (ScMAL13) was sufficient to restore growth on trisaccharides. We hypothesized that the S. pastorianus maltotriose phenotype results from regulatory interaction between S. cerevisiae maltose transcription activator and the promoter of SeAGT1 We experimentally confirmed the heterotic nature of the phenotype, and thus these results provide experimental evidence of the evolutionary origin of an essential phenotype of lager brewing strains.

Allele-specific genome editing using CRISPR-Cas9 is associated with loss of heterozygosity in diploid yeast

Targeted DNA double-strand breaks (DSBs) with CRISPR–Cas9 have revolutionized genetic modification by enabling efficient genome editing in a broad range of eukaryotic systems. Accurate gene editing is possible with near-perfect efficiency in haploid or (predominantly) homozygous genomes. However, genomes exhibiting polyploidy and/or high degrees of heterozygosity are less amenable to genetic modification. Here, we report an up to 99-fold lower gene editing efficiency when editing individual heterozygous loci in the yeast genome. Moreover, Cas9-mediated introduction of a DSB resulted in large scale loss of heterozygosity affecting DNA regions up to 360 kb and up to 1700 heterozygous nucleotides, due to replacement of sequences on the targeted chromosome by corresponding sequences from its non-targeted homolog. The observed patterns of loss of heterozygosity were consistent with homology directed repair. The extent and frequency of loss of heterozygosity represent a novel mutagenic side-effect of Cas9-mediated genome editing, which would have to be taken into account in eukaryotic gene editing. In addition to contributing to the limited genetic amenability of heterozygous yeasts, Cas9-mediated loss of heterozygosity could be particularly deleterious for human gene therapy, as loss of heterozygous functional copies of anti-proliferative and pro-apoptotic genes is a known path to cancer.

Phenotype-Independent Isolation of Interspecies Saccharomyces Hybrids by Dual-Dye Fluorescent Staining and Fluorescence-Activated Cell Sorting

Interspecies hybrids of Saccharomyces species are found in a variety of industrial environments and often outperform their parental strains in industrial fermentation processes. Interspecies hybridization is therefore increasingly considered as an approach for improvement and diversification of yeast strains for industrial application. However, current hybridization methods are limited by their reliance on pre-existing or introduced selectable phenotypes. This study presents a high-throughput phenotype-independent method for isolation of interspecies Saccharomyces hybrids based on dual dye-staining and subsequent mating of two strains, followed by enrichment of double-stained hybrid cells from a mating population by fluorescence-activated cell sorting (FACS). Pilot experiments on intra-species mating of heterothallic haploid S. cerevisiae strains showed that 80% of sorted double-stained cells were hybrids. The protocol was further optimized by mating an S. cerevisiae haploid with homothallic S. eubayanus spores with complementary selectable phenotypes. In crosses without selectable phenotype, using S. cerevisiae and S. eubayanus haploids derived from laboratory as well as industrial strains, 10 to 15% of double-stained cells isolated by FACS were hybrids. When applied to rare mating, sorting of double-stained cells consistently resulted in about 600-fold enrichment of hybrid cells. Mating of dual-stained cells and FACS-based selection allows efficient enrichment of interspecies Saccharomyces hybrids within a matter of days and without requiring selectable hybrid phenotypes, both for homothallic and heterothallic strains. This strategy should accelerate the isolation of laboratory-made hybrids, facilitate research into hybrid heterosis and offer new opportunities for non-GM industrial strain improvement and diversification.

In vivo recombination of Saccharomyces eubayanus maltose-transporter genes yields a chimeric transporter that enables maltotriose fermentation

Saccharomyces eubayanus is the non-S. cerevisiae parent of the lager-brewing hybrid S. pastorianus. In contrast to most S. cerevisiae and Frohberg-type S. pastorianus strains, S. eubayanus cannot utilize the α-tri-glucoside maltotriose, a major carbohydrate in brewer’s wort. In Saccharomyces yeasts, utilization of maltotriose is encoded by the subtelomeric MAL gene family, and requires transporters for maltotriose uptake. While S. eubayanus strain CBS 12357^T harbors four SeMALT genes which enable uptake of the α-di-glucoside maltose, it lacks maltotriose transporter genes. In S. cerevisiae, sequence identity indicates that maltotriose and maltose transporters likely evolved from a shared ancestral gene. To study the evolvability of maltotriose utilization in S. eubayanus CBS 12357^T, maltotriose-assimilating mutants obtained after UV mutagenesis were subjected to laboratory evolution in carbon-limited chemostat cultures on maltotriose-enriched wort. An evolved strain showed improved maltose and maltotriose fermentation in 7 L fermenter experiments on industrial wort. Whole-genome sequencing revealed a novel mosaic SeMALT413 gene, resulting from repeated gene introgressions by non-reciprocal translocation of at least three SeMALT genes. The predicted tertiary structure of SeMalT413 was comparable to the original SeMalT transporters, but overexpression of SeMALT413 sufficed to enable growth on maltotriose, indicating gene neofunctionalization had occurred. The mosaic structure of SeMALT413 resembles the structure of S. pastorianus maltotriose-transporter gene SpMTY1, which has high sequences identity to alternatingly S. cerevisiae MALx1, S. paradoxus MALx1 and S. eubayanus SeMALT3. Evolution of the maltotriose transporter landscape in hybrid S. pastorianus lager-brewing strains is therefore likely to have involved mechanisms similar to those observed in the present study.

Fermentation of the wort sugar maltotriose is critical for the flavor profile obtained during beer brewing. The recently discovered yeast Saccharomyces eubayanus is gaining popularity as an alternative to S. pastorianus and S. cerevisiae for brewing, however it is unable to utilize maltotriose. Here, a combination of non-GMO mutagenesis and laboratory evolution of the S. eubayanus type strain CBS 12357^T was used to enable maltotriose fermentation and improve brewing performance. The improved strain expressed a novel transporter gene, SeMALT413, which was formed by recombination between three different SeMALT maltose-transporter genes. Overexpression of SeMALT413 in CBS 12357^T confirmed its neofunctionalization as a maltotriose transporter. As the S. pastorianus maltotriose transporter SpMty1 has a mosaic structure similar to SeMalT413, maltotriose utilization likely involved similar recombination events during the domestication of current lager brewing strains. Based on a posteriori sequence analysis, the emergence of gene functions has been attributed to gene neofunctionalization in a broad range of organisms. The real-time observation of neofunctionalization during laboratory evolution constitutes an important validation of the relevance and importance of this mechanism for Darwinian evolution.

Laboratory Evolution of a Saccharomyces cerevisiae × S. eubayanus Hybrid Under Simulated Lager-Brewing Conditions

Saccharomyces pastorianus lager-brewing yeasts are domesticated hybrids of S. cerevisiae x S. eubayanus that display extensive inter-strain chromosome copy number variation and chromosomal recombinations. It is unclear to what extent such genome rearrangements are intrinsic to the domestication of hybrid brewing yeasts and whether they contribute to their industrial performance. Here, an allodiploid laboratory hybrid of S. cerevisiae and S. eubayanus was evolved for up to 418 generations on wort under simulated lager-brewing conditions in six independent sequential batch bioreactors. Characterization of 55 single-cell isolates from the evolved cultures showed large phenotypic diversity and whole-genome sequencing revealed a large array of mutations. Frequent loss of heterozygosity involved diverse, strain-specific chromosomal translocations, which differed from those observed in domesticated, aneuploid S. pastorianus brewing strains. In contrast to the extensive aneuploidy of domesticated S. pastorianus strains, the evolved isolates only showed limited (segmental) aneuploidy. Specific mutations could be linked to calcium-dependent flocculation, loss of maltotriose utilization and loss of mitochondrial activity, three industrially relevant traits that also occur in domesticated S. pastorianus strains. This study indicates that fast acquisition of extensive aneuploidy is not required for genetic adaptation of S. cerevisiae × S. eubayanus hybrids to brewing environments. In addition, this work demonstrates that, consistent with the diversity of brewing strains for maltotriose utilization, domestication under brewing conditions can result in loss of this industrially relevant trait. These observations have important implications for the design of strategies to improve industrial performance of novel laboratory-made hybrids.

Genome editing in Kluyveromyces and Ogataea yeasts using a broad-host-range Cas9/gRNA co-expression plasmid

While CRISPR-Cas9-mediated genome editing has transformed yeast research, current plasmids and cassettes for Cas9 and guide-RNA expression are species specific. CRISPR tools that function in multiple yeast species could contribute to the intensifying research on non-conventional yeasts. A plasmid carrying a pangenomic origin of replication and two constitutive expression cassettes for Cas9 and ribozyme-flanked gRNAs was constructed. Its functionality was tested by analyzing inactivation of the ADE2 gene in four yeast species. In two Kluyveromyces species, near-perfect targeting (≥96%) and homologous repair (HR) were observed in at least 24% of transformants. In two Ogataea species, Ade- mutants were not observed directly after transformation, but prolonged incubation of transformed cells resulted in targeting efficiencies of 9% to 63% mediated by non-homologous end joining (NHEJ). In an Ogataea parapolymorpha ku80 mutant, deletion of OpADE2 mediated by HR was achieved, albeit at low efficiencies (<1%). Furthermore the expression of a dual polycistronic gRNA array enabled simultaneous interruption of OpADE2 and OpYNR1 demonstrating flexibility of ribozyme-flanked gRNA design for multiplexing. While prevalence of NHEJ prevented HR-mediated editing in Ogataea, such targeted editing was possible in Kluyveromyces. This broad-host-range CRISPR/gRNA system may contribute to exploration of Cas9-mediated genome editing in other Saccharomycotina yeasts.

CRISPR-Cas9 mediated gene deletions in lager yeast Saccharomyces pastorianus

BACKGROUND

The ease of use of CRISPR-Cas9 reprogramming, its high efficacy, and its multiplexing capabilities have brought this technology at the forefront of genome editing techniques. Saccharomyces pastorianus is an aneuploid interspecific hybrid of Saccharomyces cerevisiae and Saccharomyces eubayanus that has been domesticated for centuries and is used for the industrial fermentation of lager beer. For yet uncharacterised reasons, this hybrid yeast is far more resilient to genetic alteration than its ancestor S. cerevisiae.

RESULTS

This study reports a new CRISPR-Cas9 method for accurate gene deletion in S. pastorianus. This method combined the Streptococcus pyogenes cas9 gene expressed from either a chromosomal locus or from a mobile genetic element in combination with a plasmid-borne gRNA expression cassette. While the well-established gRNA expression system using the RNA polymerase III dependent SNR52 promoter failed, expression of a gRNA flanked with Hammerhead and Hepatitis Delta Virus ribozymes using the RNA polymerase II dependent TDH3 promoter successfully led to accurate deletion of all four alleles of the SeILV6 gene in strain CBS1483. Furthermore the expression of two ribozyme-flanked gRNAs separated by a 10-bp linker in a polycistronic array successfully led to the simultaneous deletion of SeATF1 and SeATF2, genes located on two separate chromosomes. The expression of this array resulted in the precise deletion of all five and four alleles mediated by homologous recombination in the strains CBS1483 and Weihenstephan 34/70 respectively, demonstrating the multiplexing abilities of this gRNA expression design.

CONCLUSIONS

These results firmly established that CRISPR-Cas9 significantly facilitates and accelerates genome editing in S. pastorianus.

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.

Accurate and predictive antibody repertoire profiling by molecular amplification fingerprinting

High-throughput antibody repertoire sequencing (Ig-seq) provides quantitative molecular information on humoral immunity. However, Ig-seq is compromised by biases and errors introduced during library preparation and sequencing. By using synthetic antibody spike-in genes, we determined that primer bias from multiplex polymerase chain reaction (PCR) library preparation resulted in antibody frequencies with only 42 to 62% accuracy. Additionally, Ig-seq errors resulted in antibody diversity measurements being overestimated by up to 5000-fold. To rectify this, we developed molecular amplification fingerprinting (MAF), which uses unique molecular identifier (UID) tagging before and during multiplex PCR amplification, which enabled tagging of transcripts while accounting for PCR efficiency. Combined with a bioinformatic pipeline, MAF bias correction led to measurements of antibody frequencies with up to 99% accuracy. We also used MAF to correct PCR and sequencing errors, resulting in enhanced accuracy of full-length antibody diversity measurements, achieving 98 to 100% error correction. Using murine MAF-corrected data, we established a quantitative metric of recent clonal expansion-the intraclonal diversity index-which measures the number of unique transcripts associated with an antibody clone. We used this intraclonal diversity index along with antibody frequencies and somatic hypermutation to build a logistic regression model for prediction of the immunological status of clones. The model was able to predict clonal status with high confidence but only when using MAF error and bias corrected Ig-seq data. Improved accuracy by MAF provides the potential to greatly advance Ig-seq and its utility in immunology and biotechnology.

A new laboratory evolution approach to select for constitutive acetic acid tolerance in Saccharomyces cerevisiae and identification of causal mutations

BACKGROUND

Acetic acid, released during hydrolysis of lignocellulosic feedstocks for second generation bioethanol production, inhibits yeast growth and alcoholic fermentation. Yeast biomass generated in a propagation step that precedes ethanol production should therefore express a high and constitutive level of acetic acid tolerance before introduction into lignocellulosic hydrolysates. However, earlier laboratory evolution strategies for increasing acetic acid tolerance of Saccharomyces cerevisiae, based on prolonged cultivation in the presence of acetic acid, selected for inducible rather than constitutive tolerance to this inhibitor.

RESULTS

Preadaptation in the presence of acetic acid was shown to strongly increase the fraction of yeast cells that could initiate growth in the presence of this inhibitor. Serial microaerobic batch cultivation, with alternating transfers to fresh medium with and without acetic acid, yielded evolved S. cerevisiae cultures with constitutive acetic acid tolerance. Single-cell lines isolated from five such evolution experiments after 50-55 transfers were selected for further study. An additional constitutively acetic acid tolerant mutant was selected after UV-mutagenesis. All six mutants showed an increased fraction of growing cells upon a transfer from a non-stressed condition to a medium containing acetic acid. Whole-genome sequencing identified six genes that contained (different) mutations in multiple acetic acid-tolerant mutants. Haploid segregation studies and expression of the mutant alleles in the unevolved ancestor strain identified causal mutations for the acquired acetic acid tolerance in four genes (ASG1, ADH3, SKS1 and GIS4). Effects of the mutations in ASG1, ADH3 and SKS1 on acetic acid tolerance were additive.

CONCLUSIONS

A novel laboratory evolution strategy based on alternating cultivation cycles in the presence and absence of acetic acid conferred a selective advantage to constitutively acetic acid-tolerant mutants and may be applicable for selection of constitutive tolerance to other stressors. Mutations in four genes (ASG1, ADH3, SKS1 and GIS4) were identified as causative for acetic acid tolerance. The laboratory evolution strategy as well as the identified mutations can contribute to improving acetic acid tolerance in industrial yeast strains.