Deconstructing the genetic basis of spent sulphite liquor tolerance using deep sequencing of genome-shuffled yeast

How adaptive laboratory evolution can boost yeast tolerance to lignocellulosic hydrolyses

The yeast Saccharomyces cerevisiae is an excellent candidate for establishing cell factories to convert lignocellulosic biomass into chemicals and fuels. To enable this technology, yeast robustness must be improved to withstand the fermentation inhibitors (e.g., weak organic acids, phenols, and furan aldehydes) resulting from biomass pretreatment and hydrolysis. Here, we discuss how evolution experiments performed in the lab, a method commonly known as adaptive laboratory evolution (ALE), may contribute to lifting yeast tolerance against the inhibitors of lignocellulosic hydrolysates (LCHs). The key is that, through the combination of whole-genome sequencing and reverse engineering, ALE provides a robust platform for discovering and testing adaptive alleles, allowing to explore the genetic underpinnings of yeast responses to LCHs. We review the insights gained from past evolution experiments with S. cerevisiae in LCH inhibitors and propose experimental designs to optimise the discovery of genetic variants adaptive to biomass toxicity. The knowledge gathered through ALE projects is envisaged as a roadmap to engineer superior yeast strains for biomass-based bioprocesses.

Rational and evolutionary engineering of Saccharomyces cerevisiae for production of dicarboxylic acids from lignocellulosic biomass and exploring genetic mechanisms of the yeast tolerance to the biomass hydrolysate

BACKGROUND

Lignosulfonates are significant wood chemicals with a $700 million market, produced by sulfite pulping of wood. During the pulping process, spent sulfite liquor (SSL) is generated, which in addition to lignosulfonates contains hemicellulose-derived sugars-in case of hardwoods primarily the pentose sugar xylose. The pentoses are currently underutilized. If they could be converted into value-added chemicals, overall economic profitability of the process would increase. SSLs are typically very inhibitory to microorganisms, which presents a challenge for a biotechnological process. The aim of the present work was to develop a robust yeast strain able to convert xylose in SSL to carboxylic acids.

RESULTS

The industrial strain Ethanol Red of the yeast Saccharomyces cerevisiae was engineered for efficient utilization of xylose in a Eucalyptus globulus lignosulfonate stream at low pH using CRISPR/Cas genome editing and adaptive laboratory evolution. The engineered strain grew in synthetic medium with xylose as sole carbon source with maximum specific growth rate (µ_max) of 0.28 1/h. Selected evolved strains utilized all carbon sources in the SSL at pH 3.5 and grew with µ_max between 0.05 and 0.1 1/h depending on a nitrogen source supplement. Putative genetic determinants of the increased tolerance to the SSL were revealed by whole genome sequencing of the evolved strains. In particular, four top-candidate genes (SNG1, FIT3, FZF1 and CBP3) were identified along with other gene candidates with predicted important roles, based on the type and distribution of the mutations across different strains and especially the best performing ones. The developed strains were further engineered for production of dicarboxylic acids (succinic and malic acid) via overexpression of the reductive branch of the tricarboxylic acid cycle (TCA). The production strain produced 0.2 mol and 0.12 mol of malic acid and succinic acid, respectively, per mol of xylose present in the SSL.

CONCLUSIONS

The combined metabolic engineering and adaptive evolution approach provided a robust SSL-tolerant industrial strain that converts fermentable carbon content of the SSL feedstock into malic and succinic acids at low pH.in production yields reaching 0.1 mol and 0.065 mol per mol of total consumed carbon sources.. Moreover, our work suggests potential genetic background of the tolerance to the SSL stream pointing out potential gene targets for improving the tolerance to inhibitory industrial feedstocks.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

Proteome response of two natural strains of Saccharomyces cerevisiae with divergent lignocellulosic inhibitor stress tolerance

Strains of Saccharomyces cerevisiae with improved tolerance to plant hydrolysates are of utmost importance for the cost-competitive production of value-added chemicals and fuels. However, engineering strategies are constrained by a lack of understanding of the yeast response to complex inhibitor mixtures. Natural S. cerevisiae isolates display niche-specific phenotypic and metabolic diversity, encoded in their DNA, which has evolved to overcome external stresses, utilise available resources and ultimately thrive in their challenging environments. Industrial and laboratory strains, however, lack these adaptations due to domestication. Natural strains can serve as a valuable resource to mitigate engineering constraints by studying the molecular mechanisms involved in phenotypic variance and instruct future industrial strain improvement to lignocellulosic hydrolysates. We, therefore, investigated the proteomic changes between two natural S. cerevisiae isolates when exposed to a lignocellulosic inhibitor mixture. Comparative shotgun proteomics revealed that isolates respond by regulating a similar core set of proteins in response to inhibitor stress. Furthermore, superior tolerance was linked to NAD(P)/H and energy homeostasis, concurrent with inhibitor and reactive oxygen species detoxification processes. We present several candidate proteins within the redox homeostasis and energy management cellular processes as possible targets for future modification and study. Data are available via ProteomeXchange with identifier PXD010868.

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

BACKGROUND

High-temperature fermentation is desirable for the industrial production of ethanol, which requires thermotolerant yeast strains. However, yeast thermotolerance is a complicated quantitative trait. The understanding of genetic basis behind high-temperature fermentation performance is still limited. Quantitative trait locus (QTL) mapping by pooled-segregant whole genome sequencing has been proved to be a powerful and reliable approach to identify the loci, genes and single nucleotide polymorphism (SNP) variants linked to quantitative traits of yeast.

RESULTS

One superior thermotolerant industrial strain and one inferior thermosensitive natural strain with distinct high-temperature fermentation performances were screened from 124 Saccharomyces cerevisiae strains as parent strains for crossing and segregant isolation. Based on QTL mapping by pooled-segregant whole genome sequencing as well as the subsequent reciprocal hemizygosity analysis (RHA) and allele replacement analysis, we identified and validated total eight causative genes in four QTLs that linked to high-temperature fermentation of yeast. Interestingly, loss of heterozygosity in five of the eight causative genes including RXT2, ECM24, CSC1, IRA2 and AVO1 exhibited positive effects on high-temperature fermentation. Principal component analysis (PCA) of high-temperature fermentation data from all the RHA and allele replacement strains of those eight genes distinguished three superior parent alleles including VPS34, VID24 and DAP1 to be greatly beneficial to high-temperature fermentation in contrast to their inferior parent alleles. Strikingly, physiological impacts of the superior parent alleles of VPS34, VID24 and DAP1 converged on cell membrane by increasing trehalose accumulation or reducing membrane fluidity.

CONCLUSIONS

This work revealed eight novel causative genes and SNP variants closely associated with high-temperature fermentation performance. Among these genes, VPS34 and DAP1 would be good targets for improving high-temperature fermentation of the industrial yeast. It also showed that loss of heterozygosity of causative genes could contribute to the improvement of high-temperature fermentation capacities. Our findings would provide guides to develop more robust and thermotolerant strains for the industrial production of ethanol.

Development and genomic elucidation of hybrid yeast with improved glucose-xylose co-fermentation at high temperature

Enhanced capability of co-fermenting glucose and xylose at high temperature is highly desirable for yeast application in second-generation bioethanol production. Here, we obtained hybrid strains with improved glucose-xylose co-fermentation properties at high temperature by combining genome shuffling and adaptive evolution. Genome resequencing of these strains suggested predominantly inherited genetic information from one parental strain Spathaspora passalidarum SP rather than the other parental strain Saccharomyces cerevisiae ScY01, possibly due to that the CUG codon system of S. passalidarum might have systematically eliminated most of the functional proteins from S. cerevisiae through misfolding. Compared to SP, one-copy loss of a 146-kb fragment was found in the hybrid strain and regained after being evolved for a while, whereas one-copy loss of an 11-kb fragment was only found after being evolved for a longer time. Besides, the genes affected by nonsynonymous variants were also identified, especially the mutation S540F in the endoplasmic reticulum chaperon Kar2. Structural prediction indicated that S540F might change the substrate binding activity of Kar2, and thus play a role in preventing protein aggregation in yeast at high temperature. Our results illustrated genomic alterations during this process and revealed some genomic factors that might be involved to determine yeast thermotolerance.

Determinants of selection in yeast evolved by genome shuffling

BACKGROUND

Genome shuffling (GS) is a widely adopted methodology for the evolutionary engineering of desirable traits in industrially relevant microorganisms. We have previously used genome shuffling to generate a strain of Saccharomyces cerevisiae that is tolerant to the growth inhibitors found in a lignocellulosic hydrolysate. In this study, we expand on previous work by performing a population-wide genomic survey of our genome shuffling experiment and dissecting the molecular determinants of the evolved phenotype.

RESULTS

Whole population whole-genome sequencing was used to survey mutations selected during the experiment and extract allele frequency time series. Using growth curve assays on single point mutants and backcrossed derivatives, we explored the genetic architecture of the selected phenotype and detected examples of epistasis. Our results reveal cohorts of strongly correlated mutations, suggesting prevalent genetic hitchhiking and the presence of pre-existing founder mutations. From the patterns of apparent selection and the results of direct phenotypic assays, our results identify key driver mutations and deleterious hitchhikers.

CONCLUSIONS

We use these data to propose a model of inhibitor tolerance in our GS mutants. Our results also suggest a role for compensatory evolution and epistasis in our genome shuffling experiment and illustrate the impact of historical contingency on the outcomes of evolutionary engineering.

Improvement of Candida parapsilosis by genome shuffling for the efficient production of arabitol from l-arabinose

Arabitol is used in the food industry as a low-calorie sweetener. It is produced by yeasts during the biotransformation process of l-arabinose. Genome shuffling was performed in Candida parapsilosis DSM 70125, an efficient producer of arabitol, to obtain fusants with improved arabitol production ability. Four mutants from the parental library were used for the first round of genome shuffling. The best fusants, GSI-1 and GSI-10A, were subjected to a second round of genome shuffling. Finally, two fusants, GSII-3 and GSII-16, produced concentrations of arabitol that were 50% higher than that of the wild-type strain during selection culture. Under the optimal conditions established for C. parapsilosis, the two fusants produced 11.83 and 11.75 g/L of arabitol and were approximately 15-16% more efficient than the wild-type strain. Flow cytometry analysis showed that the ploidy of the new strains did not change.

Saccharomyces cerevisiae strains for second-generation ethanol production: from academic exploration to industrial implementation

The recent start-up of several full-scale 'second generation' ethanol plants marks a major milestone in the development of Saccharomyces cerevisiae strains for fermentation of lignocellulosic hydrolysates of agricultural residues and energy crops. After a discussion of the challenges that these novel industrial contexts impose on yeast strains, this minireview describes key metabolic engineering strategies that have been developed to address these challenges. Additionally, it outlines how proof-of-concept studies, often developed in academic settings, can be used for the development of robust strain platforms that meet the requirements for industrial application. Fermentation performance of current engineered industrial S. cerevisiae strains is no longer a bottleneck in efforts to achieve the projected outputs of the first large-scale second-generation ethanol plants. Academic and industrial yeast research will continue to strengthen the economic value position of second-generation ethanol production by further improving fermentation kinetics, product yield and cellular robustness under process conditions.

Quantitative Trait Loci (QTL)-Guided Metabolic Engineering of a Complex Trait

Engineering complex phenotypes for industrial and synthetic biology applications is difficult and often confounds rational design. Bioethanol production from lignocellulosic feedstocks is a complex trait that requires multiple host systems to utilize, detoxify, and metabolize a mixture of sugars and inhibitors present in plant hydrolysates. Here, we demonstrate an integrated approach to discovering and optimizing host factors that impact fitness of Saccharomyces cerevisiae during fermentation of a Miscanthus x giganteus plant hydrolysate. We first used high-resolution Quantitative Trait Loci (QTL) mapping and systematic bulk Reciprocal Hemizygosity Analysis (bRHA) to discover 17 loci that differentiate hydrolysate tolerance between an industrially related (JAY291) and a laboratory (S288C) strain. We then used this data to identify a subset of favorable allelic loci that were most amenable for strain engineering. Guided by this "genetic blueprint", and using a dual-guide Cas9-based method to efficiently perform multikilobase locus replacements, we engineered an S288C-derived strain with superior hydrolysate tolerance than JAY291. Our methods should be generalizable to engineering any complex trait in S. cerevisiae, as well as other organisms.

Seamless site-directed mutagenesis of the Saccharomyces cerevisiae genome using CRISPR-Cas9

CRISPR assisted homology directed repair enables the introduction of virtually any modification to the Saccharomyces cerevisiae genome. Of obvious interest is the marker-free and seamless introduction of point mutations. To fulfill this promise, a strategy that effects single nucleotide changes while preventing repeated recognition and cutting by the gRNA/Cas9 complex is needed. We demonstrate a two-step method to introduce point mutations at 17 positions in the S. cerevisiae genome. We show the general applicability of the method, enabling the seamless introduction of single nucleotide changes at any location, including essential genes and non-coding regions. We also show a quantifiable phenotype for a point mutation introduced in gene GSH1. The ease and wide applicability of this general method, combined with the demonstration of its feasibility will enable genome editing at an unprecedented level of detail in yeast and other organisms.

Genomic structural variation contributes to phenotypic change of industrial bioethanol yeast Saccharomyces cerevisiae

Genomic structural variation (GSV) is a ubiquitous phenomenon observed in the genomes of Saccharomyces cerevisiae strains with different genetic backgrounds; however, the physiological and phenotypic effects of GSV are not well understood. Here, we first revealed the genetic characteristics of a widely used industrial S. cerevisiae strain, ZTW1, by whole genome sequencing. ZTW1 was identified as an aneuploidy strain and a large-scale GSV was observed in the ZTW1 genome compared with the genome of a diploid strain YJS329. These GSV events led to copy number variations (CNVs) in many chromosomal segments as well as one whole chromosome in the ZTW1 genome. Changes in the DNA dosage of certain functional genes directly affected their expression levels and the resultant ZTW1 phenotypes. Moreover, CNVs of large chromosomal regions triggered an aneuploidy stress in ZTW1. This stress decreased the proliferation ability and tolerance of ZTW1 to various stresses, while aneuploidy response stress may also provide some benefits to the fermentation performance of the yeast, including increased fermentation rates and decreased byproduct generation. This work reveals genomic characters of the bioethanol S. cerevisiae strain ZTW1 and suggests that GSV is an important kind of mutation that changes the traits of industrial S. cerevisiae strains.

CRISPR/Cas9 advances engineering of microbial cell factories

One of the key drivers for successful metabolic engineering in microbes is the efficacy by which genomes can be edited. As such there are many methods to choose from when aiming to modify genomes, especially those of model organisms like yeast and bacteria. In recent years, clustered regularly interspaced palindromic repeats (CRISPR) and its associated proteins (Cas) have become the method of choice for precision genome engineering in many organisms due to their orthogonality, versatility and efficacy. Here we review the strategies adopted for implementation of RNA-guided CRISPR/Cas9 genome editing with special emphasis on their application for metabolic engineering of yeast and bacteria. Also, examples of how nuclease-deficient Cas9 has been applied for RNA-guided transcriptional regulation of target genes will be reviewed, as well as tools available for computer-aided design of guide-RNAs will be highlighted. Finally, this review will provide a perspective on the immediate challenges and opportunities foreseen by the use of CRISPR/Cas9 genome engineering and regulation in the context of metabolic engineering.

Modifying Yeast Tolerance to Inhibitory Conditions of Ethanol Production Processes

Saccharomyces cerevisiae strains having a broad range of substrate utilization, rapid substrate consumption, and conversion to ethanol, as well as good tolerance to inhibitory conditions are ideal for cost-competitive ethanol production from lignocellulose. A major drawback to directly design S. cerevisiae tolerance to inhibitory conditions of lignocellulosic ethanol production processes is the lack of knowledge about basic aspects of its cellular signaling network in response to stress. Here, we highlight the inhibitory conditions found in ethanol production processes, the targeted cellular functions, the key contributions of integrated -omics analysis to reveal cellular stress responses according to these inhibitors, and current status on design-based engineering of tolerant and efficient S. cerevisiae strains for ethanol production from lignocellulose.

Determinants of tolerance to inhibitors in hardwood spent sulfite liquor in genome shuffled Pachysolen tannophilus strains

Genome shuffling was used to obtain Pachysolen tannophilus mutants with improved tolerance to inhibitors in hardwood spent sulfite liquor (HW SSL). Genome shuffled strains (GHW301, GHW302 and GHW303) grew at higher concentrations of HW SSL (80 % v/v) compared to the HW SSL UV mutant (70 % v/v) and the wild-type (WT) strain (50 % v/v). In defined media containing acetic acid (0.70-0.90 % w/v), GHW301, GHW302 and GHW303 exhibited a shorter lag compared to the acetic acid UV mutant, while the WT did not grow. Genome shuffled strains produced more ethanol than the WT at higher concentrations of HW SSL and an aspen hydrolysate. To identify the genetic basis of inhibitor tolerance, whole genome sequencing was carried out on GHW301, GHW302 and GHW303 and compared to the WT strain. Sixty single nucleotide variations were identified that were common to all three genome shuffled strains. Of these, 40 were in gene sequences and 20 were within 5 bp-1 kb either up or downstream of protein encoding genes. Based on the mutated gene products, mutations were grouped into functional categories and affected a variety of cellular functions, demonstrating the complexity of inhibitor tolerance in yeast. Sequence analysis of UV mutants (UAA302 and UHW303) from which GHW301, GHW302 and GHW303 were derived, confirmed the success of our cross-mating based genome shuffling strategy. Whole-genome sequencing analysis allowed identification of potential gene targets for tolerance to inhibitors in lignocellulosic hydrolysates.

Cell periphery-related proteins as major genomic targets behind the adaptive evolution of an industrial Saccharomyces cerevisiae strain to combined heat and hydrolysate stress

BACKGROUND

Laboratory evolution is an important tool for developing robust yeast strains for bioethanol production since the biological basis behind combined tolerance requires complex alterations whose proper regulation is difficult to achieve by rational metabolic engineering. Previously, we reported on the evolved industrial Saccharomyces cerevisiae strain ISO12 that had acquired improved tolerance to grow and ferment in the presence of lignocellulose-derived inhibitors at high temperature (39 °C). In the current study, we used comparative genomics to uncover the extent of the genomic alterations that occurred during the evolution process and investigated possible associations between the mutations and the phenotypic traits in ISO12.

RESULTS

Through whole-genome sequencing and variant calling we identified a high number of strain-unique SNPs and INDELs in both ISO12 and the parental strain Ethanol Red. The variants were predicted to have 760 non-synonymous effects in both strains combined and were significantly enriched in Gene Ontology terms related to cell periphery, membranes and cell wall. Eleven genes, including MTL1, FLO9/FLO11, and CYC3 were found to be under positive selection in ISO12. Additionally, the FLO genes exhibited changes in copy number, and the alterations to this gene family were correlated with experimental results of multicellularity and invasive growth in the adapted strain. An independent lipidomic analysis revealed further differences between the strains in the content of nine lipid species. Finally, ISO12 displayed improved viability in undiluted spruce hydrolysate that was unrelated to reduction of inhibitors and changes in cell wall integrity, as shown by HPLC and lyticase assays.

CONCLUSIONS

Together, the results of the sequence comparison and the physiological characterisations indicate that cell-periphery proteins (e.g. extracellular sensors such as MTL1) and peripheral lipids/membranes are important evolutionary targets in the process of adaptation to the combined stresses. The capacity of ISO12 to develop complex colony formation also revealed multicellularity as a possible evolutionary strategy to improve competitiveness and tolerance to environmental stresses (also reflected by the FLO genes). Although a panel of altered genes with high relevance to the novel phenotype was detected, this study also demonstrates that the observed long-term molecular effects of thermal and inhibitor stress have polygenetic basis.