QTL analysis reveals genomic variants linked to high-temperature fermentation performance in the industrial yeast

Genes controlling hydrolysate toxin tolerance identified by QTL analysis of the natural Saccharomyces cerevisiae BCC39850

Lignocellulosic material can be converted to valorized products such as fuels. Pretreatment is an essential step in conversion, which is needed to increase the digestibility of the raw material for microbial fermentation. However, pretreatment generates by-products (hydrolysate toxins) that are detrimental to microbial growth. In this study, natural Saccharomyces strains isolated from habitats in Thailand were screened for their tolerance to synthetic hydrolysate toxins (synHTs). The Saccharomyces cerevisiae natural strain BCC39850 (toxin-tolerant) was crossed with the laboratory strain CEN.PK2-1C (toxin-sensitive), and quantitative trait locus (QTL) analysis was performed on the segregants using phenotypic scores of growth (OD600) and glucose consumption. VMS1, DET1, KCS1, MRH1, YOS9, SYO1, and YDR042C were identified from QTLs as candidate genes associated with the tolerance trait. CEN.PK2-1C knockouts of the VMS1, YOS9, KCS1, and MRH1 genes exhibited significantly greater hydrolysate toxin sensitivity to growth, whereas CEN.PK2-1C knock-ins with replacement of VMS1 and MRH1 genes from the BCC39850 alleles showed significant increased ethanol production titers compared with the CEN.PK2-1C parental strain in the presence of synHTs. The discovery of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin tolerance in S. cerevisiae indicates the roles of the endoplasmic-reticulum-associated protein degradation pathway, plasma membrane protein association, and the phosphatidylinositol signaling system in this trait. KEY POINTS: • QTL analysis was conducted using a hydrolysate toxin-tolerant S. cerevisiae natural strain • Deletion of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin-sensitivity • Replacement of VMS1 and MRH1 with natural strain alleles increased ethanol production titers in the presence of hydrolysate toxins.

Identification of genes associated with the high-temperature fermentation trait in the Saccharomyces cerevisiae natural isolate BCC39850

The fermentative model yeast Saccharomyces cerevisiae has been extensively used to study the genetic basis of stress response and homeostasis. In this study, we performed quantitative trait loci (QTL) analysis of the high-temperature fermentation trait of the progeny from the mating of the S. cerevisiae natural isolate BCC39850 (haploid#17) and the laboratory strain CEN.PK2-1C. A single QTL on chromosome X was identified, encompassing six candidate genes (GEA1, PTK2, NTA1, NPA3, IRT1, and IML1). The functions of these candidates were tested by reverse genetic experiments. Deletion mutants of PTK2, NTA1, and IML1 showed growth defects at 42 °C. The PTK2 knock-out mutant also showed significantly reduced ethanol production and plasma membrane H+ ATPase activity and increased sensitivity to acetic acid, ethanol, amphotericin B (AMB), and β-1,3-glucanase treatment. The CRISPR-Cas9 system was used to construct knock-in mutants by replacement of PTK2, NTA1, IML1, and NPA3 genes with BCC39850 alleles. The PTK2 and NTA1 knock-in mutants showed increased growth and ethanol production titers at 42 °C. These findings suggest an important role for the PTK2 serine/threonine protein kinase in regulating plasma membrane H+ ATPase activity and the NTA1 N-terminal amidase in protein degradation via the ubiquitin-proteasome system machinery, which affects tolerance to heat stress in S. cerevisiae.

Crucial aspects of metabolism and cell biology relating to industrial production and processing of Saccharomyces biomass

The multitude of applications to which Saccharomyces spp. are put makes these yeasts the most prolific of industrial microorganisms. This review considers biological aspects pertaining to the manufacture of industrial yeast biomass. It is proposed that the production of yeast biomass can be considered in two distinct but interdependent phases. Firstly, there is a cell replication phase that involves reproduction of cells by their transitions through multiple budding and metabolic cycles. Secondly, there needs to be a cell conditioning phase that enables the accrued biomass to withstand the physicochemical challenges associated with downstream processing and storage. The production of yeast biomass is not simply a case of providing sugar, nutrients, and other growth conditions to enable multiple budding cycles to occur. In the latter stages of culturing, it is important that all cells are induced to complete their current budding cycle and subsequently enter into a quiescent state engendering robustness. Both the cell replication and conditioning phases need to be optimized and considered in concert to ensure good biomass production economics, and optimum performance of industrial yeasts in food and fermentation applications. Key features of metabolism and cell biology affecting replication and conditioning of industrial Saccharomyces are presented. Alternatives for growth substrates are discussed, along with the challenges and prospects associated with defining the genetic bases of industrially important phenotypes, and the generation of new yeast strains."I must be cruel only to be kind: Thus bad begins, and worse remains behind." William Shakespeare: Hamlet, Act 3, Scene 4.

Genotypic and Phenotypic Diversity of Kluyveromyces marxianus Isolates Obtained from the Elaboration Process of Two Traditional Mexican Alcoholic Beverages Derived from Agave: Pulque and Henequen (Agave fourcroydes) Mezcal

Seven Kluyveromyces marxianus isolates from the elaboration process of pulque and henequen mezcal were characterized. The isolates were identified based on the sequences of the D1/D2 domain of the 26S rRNA gene and the internal transcribed spacer (ITS-5.8S) region. Genetic differences were found between pulque and henequen mezcal isolates and within henequen mezcal isolates, as shown by different branching patterns in the ITS-5.8S phylogenetic tree and (GTG)5 microsatellite profiles, suggesting that the substrate and process selective conditions may give rise to different K. marxianus populations. All the isolates fermented and assimilated inulin and lactose and some henequen isolates could also assimilate xylose and cellobiose. Henequen isolates were more thermotolerant than pulque ones, which, in contrast, presented more tolerance to the cell wall-disturbing agent calcofluor white (CFW), suggesting that they had different cell wall structures. Additionally, depending on their origin, the isolates presented different maximum specific growth rate (µmax) patterns at different temperatures. Concerning tolerance to stress factors relevant for lignocellulosic hydrolysates fermentation, their tolerance limits were lower at 42 than 30 °C, except for glucose and furfural. Pulque isolates were less tolerant to ethanol, NaCl, and Cd. Finally, all the isolates could produce ethanol by simultaneous saccharification and fermentation (SSF) of a corncob hydrolysate under laboratory conditions at 42 °C.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

Multiple epistatic DNA variants in a single gene affect gene expression in trans

DNA variants that alter gene expression in trans are important sources of phenotypic variation. Nevertheless, the identity of trans-acting variants remains poorly understood. Single causal variants in several genes have been reported to affect the expression of numerous distant genes in trans. Whether these simple molecular architectures are representative of trans-acting variation is unknown. Here, we studied the large RAS signaling regulator gene IRA2, which contains variants with extensive trans-acting effects on gene expression in the yeast Saccharomyces cerevisiae. We used systematic CRISPR-based genome engineering and a sensitive phenotyping strategy to dissect causal variants to the nucleotide level. In contrast to the simple molecular architectures known so far, IRA2 contained at least seven causal nonsynonymous variants. The effects of these variants were modulated by nonadditive, epistatic interactions. Two variants at the 5'-end affected gene expression and growth only when combined with a third variant that also had no effect in isolation. Our findings indicate that the molecular basis of trans-acting genetic variation may be considerably more complex than previously appreciated.

QTL mapping of a Brazilian bioethanol strain links the cell wall protein-encoding gene GAS1 to low pH tolerance in S. cerevisiae

BACKGROUND

Saccharomyces cerevisiae is largely applied in many biotechnological processes, from traditional food and beverage industries to modern biofuel and biochemicals factories. During the fermentation process, yeast cells are usually challenged in different harsh conditions, which often impact productivity. Regarding bioethanol production, cell exposure to acidic environments is related to productivity loss on both first- and second-generation ethanol. In this scenario, indigenous strains traditionally used in fermentation stand out as a source of complex genetic architecture, mainly due to their highly robust background-including low pH tolerance.

RESULTS

In this work, we pioneer the use of QTL mapping to uncover the genetic basis that confers to the industrial strain Pedra-2 (PE-2) acidic tolerance during growth at low pH. First, we developed a fluorescence-based high-throughput approach to collect a large number of haploid cells using flow cytometry. Then, we were able to apply a bulk segregant analysis to solve the genetic basis of low pH resistance in PE-2, which uncovered a region in chromosome X as the major QTL associated with the evaluated phenotype. A reciprocal hemizygosity analysis revealed the allele GAS1, encoding a β-1,3-glucanosyltransferase, as the casual variant in this region. The GAS1 sequence alignment of distinct S. cerevisiae strains pointed out a non-synonymous mutation (A631G) prevalence in wild-type isolates, which is absent in laboratory strains. We further showcase that GAS1 allele swap between PE-2 and a low pH-susceptible strain can improve cell viability on the latter of up to 12% after a sulfuric acid wash process.

CONCLUSION

This work revealed GAS1 as one of the main causative genes associated with tolerance to growth at low pH in PE-2. We also showcase how GAS1PE-2 can improve acid resistance of a susceptible strain, suggesting that these findings can be a powerful foundation for the development of more robust and acid-tolerant strains. Our results collectively show the importance of tailored industrial isolated strains in discovering the genetic architecture of relevant traits and its implications over productivity.

Rapid Enabling of Gluconobacter oxydans Resistance to High D-Sorbitol Concentration and High Temperature by Microdroplet-Aided Adaptive Evolution

Gluconobacter oxydans is important in the conversion of D-sorbitol into l-sorbose, which is an essential intermediate for industrial-scale production of vitamin C. In a previous study, the strain G. oxydans WSH-004 could directly produce 2-keto-l-gulonic acid (2-KLG). However, its D-sorbitol tolerance was poor compared with that of other common industrial G. oxydans strains, which grew well in the presence of more than 200 g/L of D-sorbitol. This study aimed to use the microbial microdroplet culture (MMC) system for the adaptive evolution of G. oxydans WSH-004 so as to improve its tolerance to high substrate concentration and high temperature. A series of adaptively evolved strains, G. oxydans MMC1-MMC10, were obtained within 90 days. The results showed that the best strain MMC10 grew in a 300 g/L of D-sorbitol medium at 40°C. The comparative genomic analysis revealed that genetic changes related to increased tolerance were mainly in protein translation genes. Compared with the traditional adaptive evolution method, the application of microdroplet-aided adaptive evolution could improve the efficiency in terms of reducing time and simplifying the procedure for strain evolution. This research indicated that the microdroplet-aided adaptive evolution was an effective tool for improving the phenotypes with undemonstrated genotypes in a short time.

Stress tolerance enhancement via SPT15 base editing in Saccharomyces cerevisiae

BACKGROUND

Saccharomyces cerevisiae is widely used in traditional brewing and modern fermentation industries to produce biofuels, chemicals and other bioproducts, but challenged by various harsh industrial conditions, such as hyperosmotic, thermal and ethanol stresses. Thus, its stress tolerance enhancement has been attracting broad interests. Recently, CRISPR/Cas-based genome editing technology offers unprecedented tools to explore genetic modifications and performance improvement of S. cerevisiae.

RESULTS

Here, we presented that the Target-AID (activation-induced cytidine deaminase) base editor of enabling C-to-T substitutions could be harnessed to generate in situ nucleotide changes on the S. cerevisiae genome, thereby introducing protein point mutations in cells. The general transcription factor gene SPT15 was targeted, and total 36 mutants with diversified stress tolerances were obtained. Among them, the 18 tolerant mutants against hyperosmotic, thermal and ethanol stresses showed more than 1.5-fold increases of fermentation capacities. These mutations were mainly enriched at the N-terminal region and the convex surface of the saddle-shaped structure of Spt15. Comparative transcriptome analysis of three most stress-tolerant (A140G, P169A and R238K) and two most stress-sensitive (S118L and L214V) mutants revealed common and distinctive impacted global transcription reprogramming and transcriptional regulatory hubs in response to stresses, and these five amino acid changes had different effects on the interactions of Spt15 with DNA and other proteins in the RNA Polymerase II transcription machinery according to protein structure alignment analysis.

CONCLUSIONS

Taken together, our results demonstrated that the Target-AID base editor provided a powerful tool for targeted in situ mutagenesis in S. cerevisiae and more potential targets of Spt15 residues for enhancing yeast stress tolerance.

Random Transfer of Ogataea polymorpha Genes into Saccharomyces cerevisiae Reveals a Complex Background of Heat Tolerance

Horizontal gene transfer, a process through which an organism acquires genes from other organisms, is a rare evolutionary event in yeasts. Artificial random gene transfer can emerge as a valuable tool in yeast bioengineering to investigate the background of complex phenotypes, such as heat tolerance. In this study, a cDNA library was constructed from the mRNA of a methylotrophic yeast, Ogataea polymorpha, and then introduced into Saccharomyces cerevisiae. Ogataea polymorpha was selected because it is one of the most heat-tolerant species among yeasts. Screening of S. cerevisiae populations expressing O. polymorpha genes at high temperatures identified 59 O. polymorpha genes that contribute to heat tolerance. Gene enrichment analysis indicated that certain S. cerevisiae functions, including protein synthesis, were highly temperature-sensitive. Additionally, the results confirmed that heat tolerance in yeast is a complex phenotype dependent on multiple quantitative loci. Random gene transfer would be a useful tool for future bioengineering studies on yeasts.

Rapid selection response to ethanol in Saccharomyces eubayanus emulates the domestication process under brewing conditions

Although the typical genomic and phenotypic changes that characterize the evolution of organisms under the human domestication syndrome represent textbook examples of rapid evolution, the molecular processes that underpin such changes are still poorly understood. Domesticated yeasts for brewing, where short generation times and large phenotypic and genomic plasticity were attained in a few generations under selection, are prime examples. To experimentally emulate the lager yeast domestication process, we created a genetically complex (panmictic) artificial population of multiple Saccharomyces eubayanus genotypes, one of the parents of lager yeast. Then, we imposed a constant selection regime under a high ethanol concentration in 10 replicated populations during 260 generations (6 months) and compared them with propagated controls exposed solely to glucose. Propagated populations exhibited a selection differential of 60% in growth rate in ethanol, mostly explained by the proliferation of a single lineage (CL248.1) that competitively displaced all other clones. Interestingly, the outcome does not require the entire time‐course of adaptation, as four lineages monopolized the culture at generation 120. Sequencing demonstrated that de novo genetic variants were produced in all propagated lines, including SNPs, aneuploidies, INDELs and translocations. In addition, the different propagated populations showed correlated responses resembling the domestication syndrome: genomic rearrangements, faster fermentation rates, lower production of phenolic off‐flavours and lower volatile compound complexity. Expression profiling in beer wort revealed altered expression levels of genes related to methionine metabolism, flocculation, stress tolerance and diauxic shift, likely contributing to higher ethanol and fermentation stress tolerance in the evolved populations. Our study shows that experimental evolution can rebuild the brewing domestication process in ‘fast motion’ in wild yeast, and also provides a powerful tool for studying the genetics of the adaptation process in complex populations.

Identification of novel genes involved in neutral lipid storage by quantitative trait loci analysis of Saccharomyces cerevisiae

BACKGROUND

The accumulation of intracellular fat depots is a polygenic trait. Therefore, the extent of lipid storage in the individuals of a species covers a broad range and is determined by many genetic factors. Quantitative trait loci analysis can be used to identify those genetic differences between two strains of the same species that are responsible for the differences in a given phenotype. We used this method and complementary approaches to identify genes in the yeast Saccharomyces cerevisiae that are involved in neutral lipid storage.

RESULTS

We selected two yeast strains, the laboratory strain BY4741 and the wine yeast AWRI1631, with a more than two-fold difference in neutral lipid content. After crossing, sporulation and germination, we used fluorescence activated cell sorting to isolate a subpopulation of cells with the highest neutral lipid content from the pool of segregants. Whole genome sequencing of this subpopulation and of the unsorted pool of segregants implicated several loci that are involved in lipid accumulation. Three of the identified genes, PIG1, PHO23 and RML2, were investigated in more detail. Deletions of these genes and the exchange of the alleles between the two parental strains confirmed that the encoded proteins contribute to neutral lipid storage in S. cerevisiae and that PIG1, PHO23 and RML2 are the major causative genes. Backcrossing of one of the segregants with the parental strains for seven generations revealed additional regions in the genomes of both strains with potential causative genes for the high lipid accumulation phenotype.

CONCLUSIONS

We identified several genes that contribute to the phenotype of lipid accumulation in an allele-specific manner. Surprisingly, no allelic variations of genes with known functions in lipid metabolism were found, indicating that the level of storage lipid accumulation is determined by many cellular processes that are not directly related to lipid metabolism.

Differential proteomic analysis by SWATH-MS unravels the most dominant mechanisms underlying yeast adaptation to non-optimal temperatures under anaerobic conditions

Elucidation of temperature tolerance mechanisms in yeast is essential for enhancing cellular robustness of strains, providing more economically and sustainable processes. We investigated the differential responses of three distinct Saccharomyces cerevisiae strains, an industrial wine strain, ADY5, a laboratory strain, CEN.PK113-7D and an industrial bioethanol strain, Ethanol Red, grown at sub- and supra-optimal temperatures under chemostat conditions. We employed anaerobic conditions, mimicking the industrial processes. The proteomic profile of these strains in all conditions was performed by sequential window acquisition of all theoretical spectra-mass spectrometry (SWATH-MS), allowing the quantification of 997 proteins, data available via ProteomeXchange (PXD016567). Our analysis demonstrated that temperature responses differ between the strains; however, we also found some common responsive proteins, revealing that the response to temperature involves general stress and specific mechanisms. Overall, sub-optimal temperature conditions involved a higher remodeling of the proteome. The proteomic data evidenced that the cold response involves strong repression of translation-related proteins as well as induction of amino acid metabolism, together with components related to protein folding and degradation while, the high temperature response mainly recruits amino acid metabolism. Our study provides a global and thorough insight into how growth temperature affects the yeast proteome, which can be a step forward in the comprehension and improvement of yeast thermotolerance.

Polygenic analysis of very high acetic acid tolerance in the yeast Saccharomyces cerevisiae reveals a complex genetic background and several new causative alleles

BACKGROUND

High acetic acid tolerance is of major importance in industrial yeast strains used for second-generation bioethanol production, because of the high acetic acid content of lignocellulose hydrolysates. It is also important in first-generation starch hydrolysates and in sourdoughs containing significant acetic acid levels. We have previously identified snf4 ^E269* as a causative allele in strain MS164 obtained after whole-genome (WG) transformation and selection for improved acetic acid tolerance.

RESULTS

We have now performed polygenic analysis with the same WG transformant MS164 to identify novel causative alleles interacting with snf4 ^E269* to further enhance acetic acid tolerance, from a range of 0.8-1.2% acetic acid at pH 4.7, to previously unmatched levels for Saccharomyces cerevisiae. For that purpose, we crossed the WG transformant with strain 16D, a previously identified strain displaying very high acetic acid tolerance. Quantitative trait locus (QTL) mapping with pooled-segregant whole-genome sequence analysis identified four major and two minor QTLs. In addition to confirmation of snf4 ^E269* in QTL1, we identified six other genes linked to very high acetic acid tolerance, TRT2, MET4, IRA2 and RTG1 and a combination of MSH2 and HAL9, some of which have never been connected previously to acetic acid tolerance. Several of these genes appear to be wild-type alleles that complement defective alleles present in the other parent strain.

CONCLUSIONS

The presence of several novel causative genes highlights the distinct genetic basis and the strong genetic background dependency of very high acetic acid tolerance. Our results suggest that elimination of inferior mutant alleles might be equally important for reaching very high acetic acid tolerance as introduction of rare superior alleles. The superior alleles of MET4 and RTG1 might be useful for further improvement of acetic acid tolerance in specific industrial yeast strains.

Selection and subsequent physiological characterization of industrial Saccharomyces cerevisiae strains during continuous growth at sub- and- supra optimal temperatures

A phenotypic screening of 12 industrial yeast strains and the well-studied laboratory strain CEN.PK113-7D at cultivation temperatures between 12 °C and 40 °C revealed significant differences in maximum growth rates and temperature tolerance. From those 12, two strains, one performing best at 12 °C and the other at 40 °C, plus the laboratory strain, were selected for further physiological characterization in well-controlled bioreactors. The strains were grown in anaerobic chemostats, at a fixed specific growth rate of 0.03 h-1 and sequential batch cultures at 12 °C, 30 °C, and 39 °C. We observed significant differences in biomass and ethanol yields on glucose, biomass protein and storage carbohydrate contents, and biomass yields on ATP between strains and cultivation temperatures. Increased temperature tolerance coincided with higher energetic efficiency of cell growth, indicating that temperature intolerance is a result of energy wasting processes, such as increased turnover of cellular components (e.g. proteins) due to temperature induced damage.

Quantitative trait loci (QTL) underlying phenotypic variation in bioethanol-related processes in Neurospora crassa

Bioethanol production from lignocellulosic biomass has received increasing attention over the past decade. Many attempts have been made to reduce the cost of bioethanol production by combining the separate steps of the process into a single-step process known as consolidated bioprocessing. This requires identification of organisms that can efficiently decompose lignocellulose to simple sugars and ferment the pentose and hexose sugars liberated to ethanol. There have been many attempts in engineering laboratory strains by adding new genes or modifying genes to expand the capacity of an industrial microorganism. There has been less attention in improving bioethanol-related processes utilizing natural variation existing in the natural ecotypes. In this study, we sought to identify genomic loci contributing to variation in saccharification of cellulose and fermentation of glucose in the fermenting cellulolytic fungus Neurospora crassa through quantitative trait loci (QTL) analysis. We identified one major QTL contributing to fermentation of glucose and multiple putative QTL's underlying saccharification. Understanding the natural variation of the major QTL gene would provide new insights in developing industrial microbes for bioethanol production.

Multimodal Microorganism Development: Integrating Top-Down Biological Engineering with Bottom-Up Rational Design

Biological engineering has unprecedented potential to solve society's most pressing challenges. Engineering approaches must consider complex technical, economic, and social factors. This requires methods that confer gene/pathway-level functionality and organism-level robustness in rapid and cost-effective ways. This article compares foundational engineering approaches - bottom-up, gene-targeted engineering, and top-down, whole-genome engineering - and identifies significant complementarity between them. Cases drawn from engineering Saccharomyces cerevisiae exemplify the synergy of a combined approach. Indeed, multimodal engineering streamlines strain development by leveraging the complementarity of whole-genome and gene-targeted engineering to overcome the gap in design knowledge that restricts rational design. As biological engineers target more complex systems, this dual-track approach is poised to become an increasingly important tool to realize the promise of synthetic biology.

Correction to: QTL analysis reveals genomic variants linked to high-temperature fermentation performance in the industrial yeast

[This corrects the article DOI: 10.1186/s13068-019-1398-7.].