Steady-state and transient-state analyses of aerobic fermentation in Saccharomyces kluyveri

Metabolic engineering of Kluyveromyces marxianus for biomass-based applications

Kluyveromyces marxianus ATCC 26,548 was cultivated in aerobic chemostats with [1-13C] and [U-13C] glucose as carbon source under three different growth conditions (0.10, 0.25, and 0.5 h-1) to evaluate metabolic fluxes. Carbon balances closed always within 97-102%. Growth was carbon limited, and the cell yield on glucose was the same. The extracellular side-product formation was very low, totaling 0.0008 C-mol C-mol-1 substrate at 0.5 h-1. The intracellular flux ratios did not show significant variation for metabolic flux analysis from labelling and biomass composition and metabolic flux ratio analysis from labelling. The observed strictly oxidative metabolism and the stability of the metabolism in terms of fluxes even at high growth rates, without triggering out the synthesis of by-products, is an extremely desired condition that underlines the potential of K. marxianus for biotechnological biomass-related applications and the comprehension of the metabolic pools and pathways is an important step to engineering this organism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03324-x.

Reconstruction and analysis of genome-scale metabolic model of weak Crabtree positive yeast Lachancea kluyveri

Lachancea kluyveri, a weak Crabtree positive yeast, has been extensively studied for its unique URC pyrimidine catabolism pathway. It produces more biomass than Saccharomyces cerevisiae due to the underlying weak Crabtree effect and resorts to fermentation only in oxygen limiting conditions that renders it as a suitable industrial host. The yeast also produces ethyl acetate as a major overflow metabolite in aerobic conditions. Here, we report the first genome-scale metabolic model, iPN730, of L. kluyveri comprising of 1235 reactions, 1179 metabolites, and 730 genes distributed in 8 compartments. The in silico viability in different media conditions and the growth characteristics in various carbon sources show good agreement with experimental data. Dynamic flux balance analysis describes the growth dynamics, substrate utilization and product formation kinetics in various oxygen-limited conditions. We have also demonstrated the effect of switching carbon sources on the production of ethyl acetate under varying oxygen uptake rates. A phenotypic phase plane analysis described the energetic cost penalty of ethyl acetate and ethanol production on the specific growth rate of L. kluyveri. We generated the context specific models of L. kluyveri growing on uracil or ammonium salts as the sole nitrogen source. Differential flux calculated using flux variability analysis helped us in highlighting pathways like purine, histidine, riboflavin and pyrimidine metabolism associated with uracil degradation. The genome-scale metabolic construction of L. kluyveri will provide a better understanding of metabolism behind ethyl acetate production as well as uracil catabolism (pyrimidine degradation) pathway. iPN730 is an addition to genome-scale metabolic models of non-conventional yeasts that will facilitate system-wide omics analysis to understand fungal metabolic diversity.

Comparative genomics of protoploid Saccharomycetaceae

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

Increasing NADH oxidation reduces overflow metabolism in Saccharomyces cerevisiae

Respiratory metabolism plays an important role in energy production in the form of ATP in all aerobically growing cells. However, a limitation in respiratory capacity results in overflow metabolism, leading to the formation of byproducts, a phenomenon known as "overflow metabolism" or "the Crabtree effect." The yeast Saccharomyces cerevisiae has served as an important model organism for studying the Crabtree effect. When subjected to increasing glycolytic fluxes under aerobic conditions, there is a threshold value of the glucose uptake rate at which the metabolism shifts from purely respiratory to mixed respiratory and fermentative. It is well known that glucose repression of respiratory pathways occurs at high glycolytic fluxes, resulting in a decrease in respiratory capacity. Despite many years of detailed studies on this subject, it is not known whether the onset of the Crabtree effect is due to limited respiratory capacity or is caused by glucose-mediated repression of respiration. When respiration in S. cerevisiae was increased by introducing a heterologous alternative oxidase, we observed reduced aerobic ethanol formation. In contrast, increasing nonrespiratory NADH oxidation by overexpression of a water-forming NADH oxidase reduced aerobic glycerol formation. The metabolic response to elevated alternative oxidase occurred predominantly in the mitochondria, whereas NADH oxidase affected genes that catalyze cytosolic reactions. Moreover, NADH oxidase restored the deficiency of cytosolic NADH dehydrogenases in S. cerevisiae. These results indicate that NADH oxidase localizes in the cytosol, whereas alternative oxidase is directed to the mitochondria.

Improvement of galactose uptake in Saccharomyces cerevisiae through overexpression of phosphoglucomutase: example of transcript analysis as a tool in inverse metabolic engineering

Through genome-wide transcript analysis of a reference strain and two recombinant Saccharomyces cerevisiae strains with different rates of galactose uptake, we obtained information about the global transcriptional response to metabolic engineering of the GAL gene regulatory network. One of the recombinant strains overexpressed the gene encoding the transcriptional activator Gal4, and in the other strain the genes encoding Gal80, Gal6, and Mig1, which are negative regulators of the GAL system, were deleted. Even though the galactose uptake rates were significantly different in the three strains, we surprisingly did not find any significant changes in the expression of the genes encoding the enzymes catalyzing the first steps of the pathway (i.e., the genes encoding Gal2, Gal1, Gal7, and Gal10). We did, however, find that PGM2, encoding the major isoenzyme of phosphoglucomutase, was slightly up-regulated in the two recombinant strains with higher galactose uptake rates. This indicated that PGM2 is a target for overexpression in terms of increasing the flux through the Leloir pathway, and through overexpression of PGM2 the galactose uptake rate could be increased by 70% compared to that of the reference strain. Based on our findings, we concluded that phosphoglucomutase plays a key role in controlling the flux through the Leloir pathway, probably due to increased conversion of glucose-1-phosphate to glucose-6-phosphate. This conclusion was supported by measurements of sugar phosphates, which showed that there were increased concentrations of glucose-6-phosphate, galactose-6-phosphate, and fructose-6-phosphate in the strain construct overexpressing PGM2.

Pyruvate decarboxylases from the petite-negative yeast Saccharomyces kluyveri

Saccharomyces kluyveri is a petite-negative yeast, which is less prone to form ethanol under aerobic conditions than is S. cerevisiae. The first reaction on the route from pyruvate to ethanol is catalysed by pyruvate decarboxylase, and the differences observed between S. kluyveri and S. cerevisiae with respect to ethanol formation under aerobic conditions could be caused by differences in the regulation of this enzyme activity. We have identified and cloned three genes encoding functional pyruvate decarboxylase enzymes (PDCgenes) from the type strain of S. kluyveri (Sk- PDC11, Sk- PDC12 and Sk- PDC13). The regulation of pyruvate decarboxylase in S. kluyveri was studied by measuring the total level of Sk- PDC mRNA and the overall enzyme activity under various growth conditions. It was found that the level of Sk- PDC mRNA was enhanced by glucose and oxygen limitation, and that the level of enzyme activity was controlled by variations in the amount of mRNA. The mRNA level and the pyruvate decarboxylase activity responded to anaerobiosis and growth on different carbon sources in essentially the same fashion as in S. cerevisiae. This indicates that the difference in ethanol formation between these two yeasts is not due to differences in the regulation of pyruvate decarboxylase(s), but rather to differences in the regulation of the TCA cycle and the respiratory machinery. However, the PDC genes of Saccharomyces/ Kluyveromyces yeasts differ in their genetic organization and phylogenetic origin. While S. cerevisiae and S. kluyveri each have three PDC genes, these have apparently arisen by independent duplications and specializations in each of the two yeast lineages.