Comparative Genomic and Transcriptomic Analysis Reveals Specific Features of Gene Regulation in Kluyveromyces marxianus

Engineering a carbon source-responsive promoter for improved biosynthesis in the non-conventional yeast Kluyveromyces marxianus

Many desired biobased chemicals exhibit a range of toxicity to microbial cell factories, making industry-level biomanufacturing more challenging. Separating microbial growth and production phases is known to be beneficial for improving production of toxic products. Here, we developed a novel synthetic carbon-responsive promoter for use in the rapidly growing, stress-tolerant yeast Kluyveromyces marxianus, by fusing carbon-source responsive elements of the native ICL1 promoter to the strong S. cerevisiae TDH3 or native NC1 promoter cores. Two hybrids, P IT350 and P IN450 , were validated via EGFP fluorescence and demonstrated exceptional strength, partial repression during growth, and late phase activation in glucose- and lactose-based medium, respectively. Expressing the Gerbera hybrida 2-pyrone synthase (2-PS) for synthesis of the polyketide triacetic acid lactone (TAL) under the control of P IN450 increased TAL more than 50% relative to the native NC1 promoter, and additional promoter engineering further increased TAL titer to 1.39 g/L in tube culture. Expression of the Penicillium griseofulvum 6-methylsalicylic acid synthase (6-MSAS) under the control of P IN450 resulted in a 6.6-fold increase in 6-MSA titer to 1.09 g/L and a simultaneous 1.5-fold increase in cell growth. Finally, we used P IN450 to express the Pseudomonas savastanoi IaaM and IaaH proteins and the Salvia pomifera sabinene synthase protein to improve production of the auxin hormone indole-3-acetic acid and the monoterpene sabinene, respectively, both extremely toxic to yeast. The development of carbon-responsive promoters adds to the synthetic biology toolbox and available metabolic engineering strategies for K. marxianus, allowing greater control over heterologous protein expression and improved production of toxic metabolites.

Coupling thermotolerance and high production of recombinant protein by CYR1N1546K mutation via cAMP signaling cascades

In recombinant protein-producing yeast strains, cells experience high production-related stresses similar to high temperatures. It is possible to increase recombinant protein production by enhancing thermotolerance, but few studies have focused on this topic. Here we aim to identify cellular regulators that can simultaneously activate thermotolerance and high yield of recombinant protein. Through screening at 46 °C, a heat-resistant Kluyveromyces marxianus (K. marxianus) strain FDHY23 is isolated. It also exhibits enhanced recombinant protein productivity at both 30 °C and high temperatures. The CYR1N1546K mutation is identified as responsible for FDHY23's improved phenotype, characterized by weakened adenylate cyclase activity and reduced cAMP production. Introducing this mutation into the wild-type strain greatly enhances both thermotolerance and recombinant protein yields. RNA-seq analysis reveals that under high temperature and recombinant protein production conditions, CYR1 mutation-induced reduction in cAMP levels can stimulate cells to improve its energy supply system and optimize material synthesis, meanwhile enhance stress resistance, based on the altered cAMP signaling cascades. Our study provides CYR1 mutation as a novel target to overcome the bottleneck in achieving high production of recombinant proteins under high temperature conditions, and also offers a convenient approach for high-throughput screening of recombinant proteins with high yields.

Transcriptome analysis of Kluyveromyces marxianus under succinic acid stress and development of robust strains

Kluyveromyces marxianus has become an attractive non-conventional yeast cell factory due to its advantageous properties such as high thermal tolerance and rapid growth. Succinic acid (SA) is an important platform molecule that has been applied in various industries such as food, material, cosmetics, and pharmaceuticals. SA bioproduction may be compromised by its toxicity. Besides, metabolite-responsive promoters are known to be important for dynamic control of gene transcription. Therefore, studies on global gene transcription under various SA concentrations are of great importance. Here, comparative transcriptome changes of K. marxianus exposed to various concentrations of SA were analyzed. Enrichment and analysis of gene clusters revealed repression of the tricarboxylic acid cycle and glyoxylate cycle, also activation of the glycolysis pathway and genes related to ergosterol synthesis. Based on the analyses, potential SA-responsive promoters were investigated, among which the promoter strength of IMTCP2 and KLMA_50231 increased 43.4% and 154.7% in response to 15 g/L SA. In addition, overexpression of the transcription factors Gcr1, Upc2, and Ndt80 significantly increased growth under SA stress. Our results benefit understanding SA toxicity mechanisms and the development of robust yeast for organic acid production. KEY POINTS: • Global gene transcription of K. marxianus is changed by succinic acid (SA) • Promoter activities of IMTCP2 and KLMA_50123 are regulated by SA • Overexpression of Gcr1, Upc2, and Ndt80 enhanced SA tolerance.

Characterization and optimization of 5´ untranslated region containing poly-adenine tracts in Kluyveromyces marxianus using machine-learning model

BACKGROUND

The 5´ untranslated region (5´ UTR) plays a key role in regulating translation efficiency and mRNA stability, making it a favored target in genetic engineering and synthetic biology. A common feature found in the 5´ UTR is the poly-adenine (poly(A)) tract. However, the effect of 5´ UTR poly(A) on protein production remains controversial. Machine-learning models are powerful tools for explaining the complex contributions of features, but models incorporating features of 5´ UTR poly(A) are currently lacking. Thus, our goal is to construct such a model, using natural 5´ UTRs from Kluyveromyces marxianus, a promising cell factory for producing heterologous proteins.

RESULTS

We constructed a mini-library consisting of 207 5´ UTRs harboring poly(A) and 34 5´ UTRs without poly(A) from K. marxianus. The effects of each 5´ UTR on the production of a GFP reporter were evaluated individually in vivo, and the resulting protein abundance spanned an approximately 450-fold range throughout. The data were used to train a multi-layer perceptron neural network (MLP-NN) model that incorporated the length and position of poly(A) as features. The model exhibited good performance in predicting protein abundance (average R2 = 0.7290). The model suggests that the length of poly(A) is negatively correlated with protein production, whereas poly(A) located between 10 and 30 nt upstream of the start codon (AUG) exhibits a weak positive effect on protein abundance. Using the model as guidance, the deletion or reduction of poly(A) upstream of 30 nt preceding AUG tended to improve the production of GFP and a feruloyl esterase. Deletions of poly(A) showed inconsistent effects on mRNA levels, suggesting that poly(A) represses protein production either with or without reducing mRNA levels.

CONCLUSION

The effects of poly(A) on protein production depend on its length and position. Integrating poly(A) features into machine-learning models improves simulation accuracy. Deleting or reducing poly(A) upstream of 30 nt preceding AUG tends to enhance protein production. This optimization strategy can be applied to enhance the yield of K. marxianus and other microbial cell factories.

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.

Fluorescence lifetime imaging of NAD(P)H upon oxidative stress in Kluyveromyces marxianus

K. marxianus is a promising cell factory for producing heterologous proteins. Oxidative stresses were raised during overexpression of heterologous proteins, leading to the shift of the redox state. How to measure the redox state of live K. marxianus cells without perturbing their growth remains a big challenge. Here, a fluorescence lifetime imaging (FLIM)-based method was developed in live K. marxianus cells. During the early exponential growth, K. marxianus cells exhibited an increased mean fluorescence lifetime (τ-mean) of NAD(P)H compared with Saccharomyces cerevisiae cells, which was consistent with the preference for respiration in K. marxianus cells and that for fermentation in S. cerevisiae cells. Upon oxidative stresses induced by high temperature or H2O2, K. marxianus cells exhibited an increased τ-mean in company with decreased intracellular NAD(P)H/NAD(P)+, suggesting a correlation between an increased τ-mean and a more oxidized redox state. The relationship between τ-mean and the expression level of a heterologous protein was investigated. There was no difference between the τ-means of K. marxianus strains which were not producing a heterologous protein. The τ-mean of a strain yielding a high level of a heterologous protein was higher than that of a low-yielding strain. The results suggested the potential application of FLIM in the non-invasive screen of high-yielding cells.

Identification of a novel gene required for competitive growth at high temperature in the thermotolerant yeast Kluyveromyces marxianus

It is important to understand the basis of thermotolerance in yeasts to broaden their application in industrial biotechnology. The capacity to run bioprocesses at temperatures above 40 °C is of great interest but this is beyond the growth range of most of the commonly used yeast species. In contrast, some industrial yeasts such as Kluyveromyces marxianus can grow at temperatures of 45 °C or higher. Such species are valuable for direct use in industrial biotechnology and as a vehicle to study the genetic and physiological basis of yeast thermotolerance. In previous work, we reported that evolutionarily young genes disproportionately changed expression when yeast were growing under stressful conditions and postulated that such genes could be important for long-term adaptation to stress. Here, we tested this hypothesis in K. marxianus by identifying and studying species-specific genes that showed increased expression during high-temperature growth. Twelve such genes were identified and 11 were successfully inactivated using CRISPR-mediated mutagenesis. One gene, KLMX_70384, is required for competitive growth at high temperature, supporting the hypothesis that evolutionary young genes could play roles in adaptation to harsh environments. KLMX_70384 is predicted to encode an 83 aa peptide, and RNA sequencing and ribo-sequencing were used to confirm transcription and translation of the gene. The precise function of KLMX_70384 remains unknown but some features are suggestive of RNA-binding activity. The gene is located in what was previously considered an intergenic region of the genome, which lacks homologues in other yeasts or in databases. Overall, the data support the hypothesis that genes that arose de novo in K. marxianus after the speciation event that separated K. marxianus and K. lactis contribute to some of its unique traits.

Meiosis-Based Laboratory Evolution of the Thermal Tolerance in Kluyveromyces marxianus

Kluyveromyces marxianus is the fastest-growing eukaryote and a promising host for producing bioethanol and heterologous proteins. To perform a laboratory evolution of thermal tolerance in K. marxianus, diploid, triploid and tetraploid strains were constructed, respectively. Considering the genetic diversity caused by genetic recombination in meiosis, we established an iterative cycle of “diploid/polyploid - meiosis - selection of spores at high temperature” to screen thermotolerant strains. Results showed that the evolution of thermal tolerance in diploid strain was more efficient than that in triploid and tetraploid strains. The thermal tolerance of the progenies of diploid and triploid strains after a two-round screen was significantly improved than that after a one-round screen, while the thermal tolerance of the progenies after the one-round screen was better than that of the initial strain. After a two-round screen, the maximum tolerable temperature of Dip2-8, a progeny of diploid strain, was 3°C higher than that of the original strain. Whole-genome sequencing revealed nonsense mutations of PSR1 and PDE2 in the thermotolerant progenies. Deletion of either PSR1 or PDE2 in the original strain improved thermotolerance and two deletions displayed additive effects, suggesting PSR1 and PDE2 negatively regulated the thermotolerance of K. marxianus in parallel pathways. Therefore, the iterative cycle of “meiosis - spore screening” developed in this study provides an efficient way to perform the laboratory evolution of heat resistance in yeast.

Characterization and modulation of endoplasmic reticulum stress response target genes in Kluyveromyces marxianus to improve secretory expressions of heterologous proteins

BACKGROUND

Kluyveromyces marxianus is a promising cell factory for producing bioethanol and that raised a demand for a high yield of heterologous proteins in this species. Expressions of heterologous proteins usually lead to the accumulation of misfolded or unfolded proteins in the lumen of the endoplasmic reticulum (ER) and then cause ER stress. To cope with this problem, a group of ER stress response target genes (ESRTs) are induced, mainly through a signaling network called unfolded protein response (UPR). Characterization and modulation of ESRTs direct the optimization of heterologous expressions. However, ESRTs in K. marxianus have not been identified so far.

RESULTS

In this study, we characterized the ER stress response in K. marxianus for the first time, by using two ER stress-inducing reagents, dithiothreitol (DTT) and tunicamycin (TM). Results showed that the Kar2-Ire1-Hac1 pathway of UPR is well conserved in K. marxianus. About 15% and 6% of genes were upregulated during treatment of DTT and TM, respectively. A total of 115 upregulated genes were characterized as ESRTs, among which 97 genes were identified as UPR target genes and 37 UPR target genes contained UPR elements in their promoters. Genes related to carbohydrate metabolic process and actin filament organization were identified as new types of UPR target genes. A total of 102 ESRTs were overexpressed separately in plasmids and their effects on productions of two different lignocellulolytic enzymes were systematically evaluated. Overexpressing genes involved in carbohydrate metabolism, including PDC1, PGK and VID28, overexpressing a chaperone gene CAJ1 or overexpressing a reductase gene MET13 substantially improved secretion expressions of heterologous proteins. Meanwhile, overexpressing a novel gene, KLMA_50479 (named ESR1), as well as overexpressing genes involved in ER-associated protein degradation (ERAD), including HRD3, USA1 andYET3, reduced the secretory expressions. ESR1 and the aforementioned ERAD genes were deleted from the genome. Resultant mutants, except the yet3Δ mutant, substantially improved secretions of three different heterologous proteins. During the fed-batch fermentation, extracellular activities of an endoxylanase and a glucanase in hrd3Δ cells improved by 43% and 28%, respectively, compared to those in wild-type cells.

CONCLUSIONS

Our results unveil the transcriptional scope of the ER stress response in K. marxianus and suggest efficient ways to improve productions of heterologous proteins by manipulating expressions of ESRTs.

Characterization and optimization of the LAC4 upstream region for low-leakage expression in Kluyveromyces marxianus

Kluyveromyces marxianus is a promising host for the production of heterologous proteins, chemicals, and bioethanol. One superior feature of this species is its capacity to assimilate lactose, which is rendered by the LAC12-LAC4 gene pair encoding a lactose permease and a β-galactosidase enzyme. Little is known about the regulation of LAC4 in K. marxianus. In this study, we showed the presence of weak glucose repression in the regulation of LAC4 and that might contribute to the leaky expression of LAC4 in the glucose medium. In a mutagenesis screen of 1000-bp LAC4 upstream region, one mutant region, named H1, drove low-leakage expression of a URA3 reporter gene in glucose medium. Two mutations inside a polyadenosine stretch (poly(A)) of 5' UTR were major contributors to the low-leakage phenotype of H1. H1 directed low-leakage expression of GFP on a plasmid and that of LAC4 in situ in the glucose medium, which was not due to the reduction of mRNA levels. Meanwhile, H1 did not affect the induction of GFP or LAC4 by lactose. Cre recombinase expressed by H1 caused lower toxicity in the repressive condition and achieved higher yield after induction, compared with that expressed by a wild-type LAC4 upstream region or a strong INU1 promoter. Our study suggested that poly(A) inside 5' UTR played a role in regulating the expression of LAC4 in the repressive condition. Meanwhile, H1 provided a base for the development of a strict inducible system for expressing industrial proteins, especially toxic proteins.