- 1State Key Laboratory of Genetic Engineering, School of Life Sciences, Fudan University, Shanghai, China
- 2Shanghai Engineering Research Center of Industrial Microorganisms, Shanghai, China
Soy leghemoglobin, when bound to heme, imparts a meat-like color and flavor and can serve as a substitute for animal-derived proteins. Enhancing cellular heme synthesis improves the recombinant expression of leghemoglobin in yeast. To achieve high-level expression of leghemoglobin A (LBA) in Kluyveromyces marxianus, a food-safe yeast, large-scale heme synthesis modules were transferred into K. marxianus using yeast artificial chromosomes (KmYACs). These modules contained up to 8 native and heterologous genes to promote the supply of heme precursors and downstream synthesis. Next, eight genes inhibiting heme or LBA synthesis were individually or combinatorially deleted, with the lsc1Δssn3Δ mutant yielding the best results. Subsequently, heme synthesis modules were combined with the lsc1Δssn3Δ mutant. In the resulting strains, the module genes were all actively expressed. Among these module genes, heterologous S. cerevisiae genes in the downstream heme synthesis pathway significantly enhanced the expression of their counterparts in K. marxianus, resulting in high heme content and LBA yield. After optimizing the medium recipe by adjusting the concentrations of glucose, glycine, and FeSO4·7H2O, a heme content of 66.32 mg/L and an intracellular LBA titer of 7.27 g/L were achieved in the engineered strain in a 5 L fermentor. This represents the highest intracellular expression of leghemoglobin in microorganisms to date. The leghemoglobin produced by K. marxianus can be utilized as a safe ingredient for plant-based protein products.
Introduction
Leghemoglobin was discovered in Glycine max in 1939 (Hoy and Hargrove, 2008). It is synthesized collaboratively by symbiotic rhizobia bacteria and the host plant in the root nodules, where it plays a crucial role in nitrogen fixation (Jiang et al., 2021). Leghemoglobin has a molecular weight of approximately 150–170 kDa and consists of four major components in the root nodules: leghemoglobin A (LBA), leghemoglobin C1 (LBC1), leghemoglobin C2 (LBC2), and leghemoglobin C3 (LBC3) (Singh and Varma, 2017). Leghemoglobin is composed of a single polypeptide chain and a heme group, similar to the structure of myoglobin found in mammals. The bound heme residues (porphyrin ring and iron ion complex) regulate the oxygen concentration and participate in the nitrogen fixation process of the rhizobia bacteria (Jiang et al., 2021). Due to the structural and functional similarity between leghemoglobin and myoglobin, leghemoglobin can produce meat-like color and flavor, making it a potential substitute for animal-derived proteins (Ismail, Hwang, and Joo, 2020).
Using microbial cell factories for leghemoglobin production offers several advantages, including industrial scalability, low cost, and short production cycles. As a result, it proves more appealing than extracting leghemoglobin from plant root nodules. Currently, secretory expression of leghemoglobin (LBC2) has only been achieved in Pichia pastoris, with a titer of 3.5 g/L (Shao et al., 2022). Intracellular expression of leghemoglobin (LBC2) in Saccharomyces cerevisiae achieved a level of 108 mg/L (Xue et al., 2023), and in Escherichia coli, intracellular expression of leghemoglobin II (LbII) reached 0.23 mg/g DCW (Arredondo-Peter et al., 1997). Different strategies have been applied to enhance hemoglobin expression in microorganisms. For instance, the deletion of VPS10 and PEP4, which are involved in the degradation of misfolded proteins in the vacuole, has been employed to reduce hemoglobin degradation, thereby increasing the expression of human hemoglobin (Ishchuk et al., 2021). By using a strong promoter and increasing the copy number of the hemoglobin expression cassette, the expression of leghemoglobin was significantly improved (Shao et al., 2022). Another successful strategy to promote hemoglobin production is to enhance cellular heme synthesis. Heme binds to the nascent globin chains to promote the co-translational folding of the growing polypeptide (Komar et al., 1997). The supplementation of heme significantly contributed to the solubility and half-life of recombinant human hemoglobin expressed in E. coli (Weickert and Curry, 1997). In S. cerevisiae, a five-fold increase in bound heme levels was achieved by eliminating spatial barriers in the heme synthesis process, enhancing the synthesis of heme precursors, and moderately improving the rate-limiting steps of heme synthesis. Consequently, the engineered strains exhibited a 46% increase in the expression of soybean hemoglobin (Liu et al., 2014). In P. pastoris, the overexpression of all heme biosynthetic pathway genes resulted in a nine-fold increase in the secretion of active soybean hemoglobin H (Shao et al., 2022).
Heme is widely synthesized in animals, plants, and microorganisms. The precursor for heme synthesis is 5-aminolevulinic acid (5-ALA). Depending on the species, the synthesis of 5-ALA can occur via the C5 pathway or the C4 pathway (Su et al., 2023). Most bacteria use the C5 pathway, which involves the catalysis of glutamyl-tRNA synthetase (GltX), glutamyl-tRNA reductase (HemA), and glutamate-1-semialdehyde aminotransferase (HemL) to convert L-glutamate into 5-ALA (Yang et al., 2016). In eukaryotes, the C4 pathway is utilized, where 5-ALA is synthesized from glycine and succinyl-CoA by ALA synthase (Hem1). The conversion of 5-ALA to heme requires seven downstream enzymes: porphobilinogen synthase (Hem2), porphobilinogen deaminase (Hem3), uroporphyrinogen III synthase (Hem4), uroporphyrinogen III decarboxylase (Hem12), coproporphyrinogen III oxidase (Hem13), protoporphyrinogen oxidase (Hem14), and ferrochelatase (Hem15) (Kornitzer and Roy, 2020). Different metabolic engineering strategies have been applied to improve the amount of cellular heme. For example, improving the supply of precursors for heme synthesis increases heme production. In Corynebacterium glutamicum, the metabolic flux of heme precursor synthesis was enhanced by overexpressing Rhodobacter capsulatus alaS in the C4 pathway, Salmonella typhimurium hemAM and E. coli hemL in the C5 pathway. The synthesis of 5-ALA was increased by at least 5.6-fold, resulting in a 4.49-fold increase in heme content (Ko et al., 2021). Supplementation of L-glutamate, a precursor in the C5 pathway, resulted in an increase in the total heme content from 115.5 mg/L to 239.3 mg/L (Zhao, Choi, and Lee, 2018). Heme production can also be improved by alleviating the repression of key heme enzymes in cells. In S. cerevisiae, the transcription factor Rox1 inhibits the expression of HEM13 under high heme conditions (Keng, 1992), and HAP1 activates the transcription of ROX1 (Martínez et al., 2015). The protein complex encoded by SSN6 and TUP1 inhibits the repressive regulatory genes of ROX1 (Amillet, Buisson, and Labbe-Bois, 1995). Deletion of ROX1, HAP1, SSN6 or TUP1 significantly reduces the inhibition of HEM13 in aerobic cells, thereby increasing cellular heme synthesis (Kastaniotis and Zitomer, 2000). Combined strategies were applied to achieve the highest heme content in several microbes. For instance, in E. coli, the overexpression of enzymes in the downstream pathway, along with the heme transporter CcmABC, and the elimination of heme degradation and competitive pathway genes, yielded a total heme titer of 115.5 mg/L and an extracellular heme titer of 73.4 mg/L (Zhao, Choi, and Lee, 2018). In C. glutamicum, the combination of C4 and C5 pathways increased the carbon flux toward heme and porphyrin intermediate biosynthesis, resulting in a total heme titer of 309 mg/L and an extracellular heme titer of 242 mg/L (Ko et al., 2021). In S. cerevisiae, eleven genes involved in heme biosynthesis, Fe(II) transportation, serine and glycine metabolism were modified, resulting in an intracellular heme titer of 53.5 mg/L, which was 70 times higher than the titer in the control strain (Ishchuk et al., 2022).
Kluyveromyces marxianus (K. marxianus) is a food-grade yeast certified by the U.S. Food and Drug Administration and the European Union Food Safety Authority (Lane and Morrissey, 2010). K. marxianus possesses various advantageous properties for industrial production, including thermotolerance (Raimondi et al., 2013), high biomass production, fast growth rate, and a broad spectrum of utilizable carbon sources, which make it a promising microbial cell factory for producing heterologous proteins, especially food-grade proteins. In this study, K. marxianus was applied in the production of leghemoglobin (LBA) for the first time. Yeast artificial chromosomes (KM-YACs) were used to introduce designed heme synthesis modules into K. marxianus, and genes affecting heme accumulation were deleted, resulting in the remodeling of the heme synthesis pathway. The crosstalk between heterologous and native heme synthesis genes was investigated. In an engineered strain, KMLSC4D, the relative level of heme and leghemoglobin increased by 6.57 times and 113%, respectively. After optimization of the media recipe, the content of heme in KMLSC4D reached 66.32 mg/L, and the titer of LBA reached 7.27 g/L in a 5 L fermentor, which was the highest intracellular expression of leghemoglobin in microorganisms so far. The LBA expressed by food-grade K. marxianus provides a safe ingredient for substitutes of animal-based proteins.
Materials and methods
Plasmids and strains
Plasmids used in this study were listed in Supplementary Table S1. To construct a LBA-expressing plasmid, ORF of LBA (gene ID:100527427) was inserted between the Not I and Spe I sites of pUKDN115 (Duan et al., 2019). The resulting plasmid was named pLBA. A His6-tag was fused to the C-terminus of LBA in pLBA to obtain pLBA-His6. To construct a plasmid carrying a single-gene cassette of HEM1Km, ORF of HEM1 (gene ID: 34715801), promoter and terminator of ADH1 (gene ID:34716095) were amplified from the genome of FIM1ΔU, and ligated with pMD18-T vector (6,011, Takara, China) by Seamless Assembly Mix (RK21020, ABclonal, China). The resulting plasmid was named LHZ1081. Similarly, plasmids carrying other single-gene cassettes were constructed. ORF used in the construction included HEM2Km (gene ID: 34714765), HEM3Sc (gene ID: 851322), HEM4Km (gene ID: 34713986), HEM12Sc (gene ID: 851617), HEM13Km (gene ID: 34716108), HEM14Sc (gene ID: 856733), HEM15Sc (gene ID: 854347), GLTXKm (gene ID: 34715977), HEMASty (gene ID: NC_003197.2), and HEMLEco (gene ID: 946892). Promoter and terminator of PDC1 (gene ID: 34717303), PGK1 (gene ID: 34714197), ENO2 (gene ID: 34714125), INU1 (gene ID: 34714176), AFT1 (gene ID: 8197809), TEF (gene ID: 4619833), OM45 (gene ID: 34716359), FBA1 (gene ID: 34716374) and HXT4 (gene ID: 34716967) were used in the construction. The resulting plasmids were named LHZ1110∼LHZ1119. To construct a plasmid carrying a multiple-gene cassette of HEM2Km-HEM3Sc, HEM2Km cassette was released from LHZ1110 by Sac II digestion, and HEM3Sc cassette was released from LHZ1111 by Swa I and Sac I digestion. Two released cassettes were ligated together by Seamless Assembly Mix and then ligated with pMD18-T vector. The resulting plasmid was named LHZ1083. Similarly, a plasmid carrying cassette of HEMASty-HEMLEco-GLTXKm-HEM1Km was constructed and named LHZ1082. A plasmid carrying cassette of HEM4Km-HEM12Sc-HEM13Km-HEM14Sc-HEM15Sc, was constructed and named LHZ1084. Sequences of single-gene and multiple-gene cassettes are listed in Supplementary Table S2. To construct CRISPR plasmids (LHZ1128∼1135) used in deleting genes, primers containing 20 bp target sequences were annealed in pairs and inserted into Sap I or Aar I sites of LHZ531 (Shi et al., 2021). Primers used in the construction are listed in Supplementary Table S3.
The strains used in this study are listed in Table 1. All K. marxianus strains originated from FIM-1 (China General Microbiological Culture Collection Center, CGMCC No.10621) (Zhou et al., 2018). To delete LSC1 (gene ID: 34715035), 500 bp sequences upstream and downstream of LSC1 ORF were amplified and ligated together as the donor sequence. CRISPR plasmid LHZ1128 was co-transformed with the donor sequence into T1a by a lithium acetate method to obtain a LSC1 deletion mutant (KML1) (Iborra, 1993). Using the same strategy, LSC2 (gene ID: 34714005), ROX1 (gene ID: 34718061), HAP1 (gene ID: 34716284), SSN3 (gene ID: 34717610), TUP1 (gene ID: 34714724), VPS10 (gene ID: 34714525) and PEP4 (gene ID: 34717789) were individually deleted in T1a to obtain KML2, KMR1, KMH1, KMS3, KMT1, KMV10, and KMP4, respectively. SSN3 was deleted in KML1 to obtain KMLS. PEP4 was deleted in KML1 to obtain KMLP. ROX1 was deleted in KMLS to obtain KMLSR. TUP1 was deleted in KMLS to obtain KMLST. VPS10 was deleted in KMLP to obtain KMLPV. ROX1 was deleted in KMLPV to obtain KMLPVR. SSN3 was deleted in KMLPV to obtain KMLPVS. TUP1 was deleted in KMLPVS to obtain KMLPVST. HIS3 (gene ID:34716654) and TRP1 (gene ID:34717357) were deleted in FIM1ΔU to obtain KMHW. HIS3 and TRP1 were deleted in KMLS to obtain KMLSHW. Construction of strains containing KmYACs was described below.
Construction of KmYAC carrying designed heme synthesis modules
KmYAC vector (LHZ1015) was composed of HIS3, ARS1/CEN5, HphMX4, TRP1 and KmTEL-filler-KmTEL cassette, in which KmTEL referred to 25 tandem repeats of K. marxianus telomeric sequence (Wu et al., 2023). To construct the KmYAC [C4+Down], LHZ1015 was digested with Not Ⅰ and BamH Ⅰ to release two arms and the filler sequence. After dephosphorylation, the digested product was ligated with HEM1Km cassette released from LHZ1081, HEM2Km-HEM3Sc cassette released from LHZ111, HEM4Km-HEM12Sc-HEM13Km-HEM14Sc-HEM15Sc cassette released from LHZ1084, by Seamless Assembly Mix. The ligation product, containing KmYAC [C4+Down], was transformed into KMΔHW and KMLSHW, to obtain KMC4D and KMLSC4D, respectively. To construct the KmYAC [C4+C5+Down], the digested and dephosphorylated product of LHZ1015 was ligated with GLTXKm-HEMASty-HEMLEco-HEM1Km cassette released from LHZ1082, HEM2Km-HEM3Sc cassette released from LHZ111, HEM4Km-HEM12Sc-HEM13Km-HEM14Sc-HEM15Sc cassette released from LHZ1084. The ligation product, containing KmYAC [C4+C5+Down], was transformed into KMΔHW and KMLSHW, to obtain KMC45D and KMLSC45D, respectively. Sequences of [C4+Down] and [C4+C5+Down] modules were listed in Supplementary Table S2.
qPCR
Cells were grown in 50 mL YP medium (40 g/L Glucose, 20 g/L Yeast extract) for 48 h or 72 h. RNA was extracted from cells using a ZR Fungal/Bacterial Miniprep kit (R2014, Zymo Research, United States) and was reverse transcribed using HiScript III All-in-one RT SuperMix Perfect (R333, Vazyme, China). The qPCR was performed using ChamQ Universal SYBR qPCR TB Master Mix (Q711-02, Vazyme, China) in LightCyclery480 (Roche, United States). Primers used in qPCR are listed in Supplementary Table S3.
Stability of KmYAC and pLBA
KMLSC4D-LBA were grown in YD medium overnight and culture was inoculated into fresh 50 mL YD medium to start at an OD600 of 0.2. Cells were grown for 24 h (about 7 generations). This process was repeated for 10 consecutive days (approximately 70 generations). The diluted cultures collected after indicated generations were spread on YPD (10 g/L yeast extract, 20 g/L hipolypepton, 20 g/L glucose, 20 g agar) plates, synthetic dropout media without histidine and tryptophan (SC-His-Trp) plates, and synthetic dropout media without histidine, tryptophan and uracil (SC-His-Trp-Ura) plates (Amberg, Burke, and Strathern, 2005). The plates were incubated at 30°C for 1–2 days. The percentage of cells containing KmYAC [C4+Down] was determined by calculating the ratio of colony numbers on SC-His-Trp plates to YPD plates. The percentage of cells containing both KmYAC [C4+Down] and pLBA was determined by calculating the ratio of colony numbers on SC-His-Trp-Ura plates to YPD plates.
Fed-batch fermentation
Fed-batch fermentations were carried out in a 5 L bioreactor (BXBIO, Shanghai, China). A seed culture was prepared by growing cells in 150 mL SM medium (0.5% (NH4)2SO4, 1% glucose, 0.05% MgSO4·7H2O, 0.3% KH2PO4) for 16 h. The seed culture was inoculated into 1.5 L SM medium supplemented with 0.1 g/L glycine, 6% glucose, and 100 µM FeSO4·7H2O. A total of 1,300 g glucose was continuously added, and the fermentation temperature was maintained at 32°C. The pH was automatically controlled at 5.5 using ammonia, and sterile air was supplied at a rate of 3 L/min for oxygenation. Dissolved oxygen and agitation speed were controlled in a cascade manner within the range of 200–850 rpm.
Measuring heme content
Cells were grown in 50 mL YD medium for 72 h. The total content of heme and free porphyrin in the cells was measured by a fluorescence-based method as described before (Liu et al., 2014). To accurately measure the heme content in cells grown in a feed-batch fermentor, a modified HPLC method was applied (Zhao, Choi, and Lee, 2018). A total of 1 mL of culture was collected at the indicated time. Cells were pelleted and mixed with 1 mL of acidic acetone (Ko et al., 2021). After centrifugation, the supernatant was mixed with 1 mL of methanol. A hemin chloride standard solution was prepared using acidic acetone and diluted with methanol. The samples were separated on a ZORBAX Eclipse Plus C18 column (250 × 4.6 mm, 5 μm) using a 1,260 Infinity II LC System (Agilent, United States). The column temperature was set at 40°C, and the flow rate was 0.6 mL/min. The mobile phase consisted of 0.05% trifluoroacetic acid (Solution A) and HPLC-grade methanol (Solution B). The gradient program was as follows: 0–1 min, 30% solution B; 1–19 min, a transition from solution A to solution B; 19–30 min, 100% solution B; 30–31 min, a transition from solution B to solution A until 30% solution B was reached; 31–35 min, 30% solution B. The absorbance at 400 nm was measured, and a standard curve was plotted using the peak heights and the concentration of the hemin chloride standard. The heme content in the samples was quantified based on the standard curve.
Comparison and quantification of LBA titer
Cells were grown in 50 mL YD medium for 72 h. Proteins were extracted by a modified post-alkaline method (Kushnirov, 2000). In brief, 100 μL culture was pelleted and cells were washed by H2O. Cells were resuspended in 200 μL NaOH and kept at room temperature for 5 min. Cells were pelleted again and resuspended in 125 μL SDS-PAGE sample buffer (50 mM Tris-HCl (pH 6.8), 2% SDS, 0.1% Bromophenol blue, 10% Glycerol and 1% 2-Mercaptoethanol). Samples were boiled and centrifuged. The supernatant was subjected to SDS-PAGE analysis. LBA-His6 was detected by Western blot analysis by using anti-His antibody (AE003, Abclonal, China). Bands in SDS-PAGE were scanned using the Gene Gnome system (Syngene, United Kingdom). The intensity of bands at the position corresponding to LBA in the samples not expressing LBA was considered as the background signal, which was then subtracted from the intensities of bands in the samples expressing LBA. The adjusted intensities of LBA bands were then used to compare the relative yield of LBA in different samples. To quantify the yield of LBA in cells grown in a feed-batch fermentor, cells were collected at the indicated time and proteins were extracted as described above. Samples containing varying concentrations of lactoglobulin (L3908, Sigma-Aldrich, United States) were subjected to SDS-PAGE alongside the LBA samples. Lactoglobulin and LBA are both globular proteins and have a similar size. Lactoglobulin displays a uniform band on SDS-PAGE, making it suitable for the quantitative analysis. The intensities of lactoglobulin bands were quantified, and a standard curve correlating band intensity to protein concentration was established. The LBA concentrations in the culture were then calculated based on this curve.
Results
Transferring designed heme synthesis modules by KmYACs to improve heme synthesis and leghemoglobin production
The bioinformatics analysis revealed the presence of a complete heme synthesis pathway in K. marxianus. Similar to S. cerevisiae, K. marxianus possesses the C4 pathway for the synthesis of 5-ALA and the key gene GLTX involved in the first step of the C5 pathway. The expressions of heme synthesis genes were investigated in T1a strain. T1a originated from a wild-type FIM1ΔU strain and displayed an improved yield of heterologous proteins by attenuating autophagy (Liu et al., 2018). After T1a was cultured in YD medium for 48 h, HEM1, HEM13, HEM15, and HEM2 were highly expressed, compared to a housekeeping gene SWC4, while the mRNA levels of HEM3, HEM4, HEM12, and HEM14 were comparable to that of SWC4 (Figure 1A). A LBA-expressing plasmid (pLBA) was transformed into T1a to obtain T1a-LBA and LBA was detected in the lysate of T1a-LBA (Figure 1B). To confirm the expression of LBA, a His6 tag was fused to the C-terminus of LBA and transformed into T1a to obtain T1a-LBA-His6. A clear band was detected by Western blot using anti-His6 antibody in the cell lysate of T1a-LBA-His6, but not in that of T1a, indicating that the band corresponded to LBA-His6 (Figure 1C). The theoretical MW of LBA-His6 is 16.2 KDa, but the LBA-His6 band detected by Western was still lower than 15 KDa. Previous reports showed that bovine hemoglobin, deer mouse hemoglobin, and cowpea leghemoglobin II exhibited similar faster mobilities in SDS-PAGE (Arredondo-Peter et al., 1997; Natarajan et al., 2011). A recent study observed LBA expressed intracellularly in C. glutamicum exhibiting a migration pattern on SDS-PAGE also below 15 kDa (Wang et al., 2023). The results indicate that LBA exhibits a faster mobility in SDS-PAGE, which might be related to the conformation of LBA. Compared with T1a, the mRNA levels of HEM2, HEM3 and HEM15 increased by 87%, 103%, and 115%, respectively in T1a-LBA (Figure 1A). Therefore, expression of heterologous LBA in T1a significantly increased the transcription levels of certain heme synthesis genes. After culturing for 72 h, the total content of heme and free porphyrin in T1a-LBA was 92.8% higher than the content in T1a (Figure 1D). This increase might be related to the elevated expression of HEM2, HEM3, and HEM15 in T1a. Leghemoglobin needs to combine with heme to form the stable form of leghemoglobin-heme. The consumption of heme might alleviate the feedback inhibition on the enzyme(s) within the heme biosynthesis pathway, thereby promoting heme synthesis. Therefore, the synthesis of hemoglobin and heme can mutually promote each other (Ishchuk et al., 2021; Shao et al., 2022).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g001.jpg)
FIGURE 1. Correlation between heme synthesis and LBA production in K. marxianus. (A) Relative mRNA levels of the heme synthesis genes pathway in T1a and T1a-LBA. Cells were collected after culturing in YD medium for 48 h. The mRNA level of the heme synthesis gene was calculated relative to that of SWC4. Values represent the mean ± S.D. (n = 3). (B) Expression of LBA in T1a. pLBA was transformed into T1a, resulting in T1a-LBA. Cells were collected after culturing in YD medium for 72 h. The cell lysate was subjected to SDS-PAGE. Four independent clones (1–4) of T1a-LBA were analyzed. (C) Western of cell lysate of T1a and T1a-LBA-His6. pLBA-His6 was transformed into T1a, resulting in T1a-LBA-His6. Cells were collected after culturing in YD medium for 72 h. The cell lysate was subjected to SDS-PAGE and Western blot. (D) Relative Content of heme and free porphyrin in T1a and T1a-LBA. Cells were cultured in YD for 72 h. The content in T1a was designated as 1. Values represent the mean ± S.D. (n = 3). ***, p < 0.001.
Therefore, we aimed to increase the production of LBA in K. marxianus by enhancing heme synthesis. To design a novel heme synthesis module, we selected heme synthesis genes from different species, including three genes from the C5 pathway, the native K. marxianus HEM1 from the C4 pathway, and seven genes from the heme downstream synthesis pathway (Figure 2A). Each gene’s expression was driven by a strong promoter and terminator of K. marxianus in a single-gene cassette (Figure 2A) (Wu et al., 2022). Single-gene cassettes were ligated to construct multiple-gene cassettes. Different single-gene and multiple-gene cassettes were then combined to build two designed heme synthesis modules, [C4+Down] and [C4+C5+Down]. [C4+down] module included native HEM1Km from the C4 pathway, three native genes (HEM2Km, HEM4Km, HEM13Km) and four S. cerevisiae genes (HEM3Sc, HEM12Sc, HEM14Sc, HEM15Sc) from downstream pathways, while the [C4+C5+down] module contained an additional native GTLXKm and two bacterial genes (HEMASty, HEMLEco) from C5 pathways. Two modules were ligated with a yeast artificial chromosome (KmYAC), respectively. The resultant KmYAC [C4+Down] and KmYAC [C4+C5+Down] were transformed into KMΔHW to obtain KMC4D and KMC45D, respectively (Figure 2B). After 72 h of culturing, the total heme/free porphyrin content in KMC4D and KMC45D increased by 23% and 28%, respectively, compared to that in KMΔHW (Figure 2C). pLBA was transformed into KMΔHW to obtain KMΔHW-LBA, KMC4D-LBA and KMC45D-LBA, respectively. Compared with KMΔHW-LBA, the levels of LBA in KMC4D-LBA and KMC45D-LBA increased by 10.1% and 14.0%, respectively, (Figure 2D). These results demonstrated that the designed heme synthesis modules, transferred by KmYACs, significantly improved heme synthesis and the yield of LBA in K. marxianus.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g002.jpg)
FIGURE 2. Remodeling of K. marxianus heme synthesis pathway by YAC. (A) Genes included in the heme synthesis module. (B) Flow chart of constructing KM-YAC that carried designed heme synthesis modules. KmYAC [C4+Down] and KmYAC [C4+C5+Down] were transformed into KMDHW to obtain KMC4D and KMC45D, respectively. (C) Relative content of heme and free porphyrin in KMΔHW, KMC4D, and KMC45D. Cells were cultured in YD medium for 72 h. The content in KMΔHW was designated as 1. Values represent the mean ± S.D. (n = 3). **, p < 0.01. ns, not significant. (D) Relative levels of LBA in KMΔHW-LBA, KMC4D-LBA and KMC45D-LBA. pLBA was transformed into KMΔHW, KMC4D, and KMC45D, resulting in KMΔHW-LBA, KMC4D-LBA and KMC45D-LBA, respectively. The transformants were then cultured in YD medium for 72 h. The level of LBA in KMΔHW-LBA was designated as 1. Values represent the mean ± S.D. (n = 3). *, p < 0.05.
Improving heme synthesis and LBA production by deleting inhibitory genes
To improve cellular heme content and the yield of LBA, several genes exhibiting inhibitory effects on heme synthesis or LBA production were selected for deletion (Figure 3A). In C. glutamicum, sucC and sucD encode two subunits of succinyl-CoA ligase, responsible for converting succinyl-CoA into succinate. Deletion of sucC and sucD abolishes the flux to succinate and instead enhances the synthesis of 5-ALA from succinyl-CoA and glycine (Yang et al., 2016). Deletion of ROX1, HAP1, SSN6, or TUP1 reduces the inhibition of HEM13, which encodes a key enzyme in the downstream synthesis pathway of heme (Keng, 1992; Kastaniotis and Zitomer, 2000; Martínez et al., 2015). The deletion of VPS10 and PEP4 reduces the degradation of hemoglobin in the vacuole (Ishchuk et al., 2021). Based on sequence homology analysis, 8 genes were identified as homologues of the above-mentioned genes in K. marxianus. They were individually deleted in T1a, resulting in eight single-gene deletion mutants. The biomass, the total content of heme and free porphyrin, and LBA production of T1a and deletion mutants were measured and compared side by side (Figure 3B). Deletion of ROX1 or VPS10 significantly reduced biomass, while other deletions did not affect biomass. Deletion of LSC1, LSC2, HAP1, SSN3, TUP1 or PEP4 significantly improves the total content of heme and free porphyrin.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g003.jpg)
FIGURE 3. Improving heme synthesis and LBA production by deleting inhibitory genes. (A) Schematic diagram of effects of deletions. Deletions were expected to abolish the synthesis of succinate (lsc1Δ, lsc2Δ), repression of HEM13 (rox1Δ, hap1Δ, tup1Δ, ssn3Δ), or degradation of misfolded proteins (vps10Δ, pep4Δ). (B,C) Effect of deleting a single gene (B) or multiple genes (C) on the biomass, heme and free porphyrin content, and LBA production. Predicted effects of deletion(s) on enhancing 5-ALA synthesis, reducing repression of HEM13, and reducing protein degradation were checked. Genes were deleted in T1a. OD600 and the total content of heme and free porphyrin in T1a and deletion mutant were compared after culturing in YD medium for 72 h. To compare the production of LBA, pLBA was transformed into T1a and mutants, respectively. The resulting transformants were cultured in YD medium for 72 h. Heme content and level of LBA in T1a were designated as 1 when compared to mutants. Values represent the mean ± S.D. (n = 3). *, p < 0.05, **, p < 0.01, ***, p < 0.001.
Among the eight mutants, the lsc1Δ mutant demonstrated the greatest improvement in heme and leghemoglobin production. Specifically, compared with T1a, the lsc1Δ mutant exhibited a 35% increase in the total content of heme and free porphyrin, as well as a 28% increase in LBA production. Therefore, eight multiple-gene deletion mutants were constructed in the lsc1Δ background (Figure 3C). All eight multiple-gene deletion mutants showed decreased biomass after 72 h of culturing, indicating a burden on cell growth due to multiple deletions. Meanwhile, four mutants showed simultaneously increased heme content and LBA production. Among them, the lsc1Δssn3Δ mutant performed the best, displaying an 87% increase in the total content of heme and free porphyrin, as well as a 97% increase in LBA production. This mutant was designated as KMLS and selected for further engineering.
Remodeling the heme synthesis pathway by combining designed heme synthesis modules and deletions of inhibitory genes
Our results indicate that the designed heme synthesis modules transferred by KmYACs, as well as the deletions of LSC1 and SSN3, significantly improve heme synthesis and the yield of LBA in K. marxianus. Therefore, we aimed to combine these two strategies to achieve further improvement. HIS3 and TRP1 were deleted in KMLS to create auxotrophic markers to accept KmYAC and the resultant strain named KMLSHW. KmYAC [C4+Down] and KmYAC [C4+C5+Down] were transformed into KMLSHW to obtain KMLSC4D and KMLSC45D, respectively.
The expressions of heme synthesis genes in KMLSHW, KMLSC4D and KMLSC45D were investigated (Figure 4). In KMLSC4D, the expression of HEM1Km, HEM2Km, HEM4Km, or HEM13Km was derived from one native gene in the genome and an additional copy transferred by KmYAC (Figure 4B, gray column). The relative mRNA levels of HEM1KM, HEM2Km, HEM4Km, and HEM13Km in KMLSC4D were 32, 136, 43, and 82 times higher than the corresponding genes in KMLSHW, where the expression of HEM1Km, HEM2Km, HEM4Km, and HEM13Km only originated from the native genes (Figure 4A, black column). The results indicate transferring K. marxianus heme synthesis gene by KmYAC significantly upregulates the overall expression of these genes. Four heme synthesis genes from S. cerevisiae, including HEM3Sc, HEM12Sc, HEM14Sc, and HEM15Sc were all actively expressed in KMLSC4D (Figure 4B, white columns). Notably, the presence of these heterologous genes significantly enhanced the expression of their orthologues (Figure 4B, black columns), suggesting a mechanism to promote the co-upregulation of orthologous genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g004.jpg)
FIGURE 4. Expressions of heme synthesis genes in strains containing designed heme synthesis modules and deletions of inhibitory genes. (A–C) The relative mRNA levels of the heme synthesis genes pathway in KMLSHW (A), KMLSC4D (B), and KMLSC45D (C). KmYAC [C4+Down] and KmYAC [C4+C5+Down] were transformed into KMLSHW to obtain KMLSC4D and KMLSC45D, respectively. Cells were collected after culturing in YD medium for 72 h. The mRNA levels of heme synthesis genes were calculated relative to that of SWC4. The relative mRNA levels of genomic K. marxianus genes were represented by black columns, while those of the heterologous genes introduced by KmYAC were represented by white columns. The relative mRNA levels of K. marxianus genes expressed from one genomic copy and one additional copy introduced by KmYAC were represented by grey columns. The values represent the mean ± S.D. (n = 3).
In KMLSC45D, the expression of GLTXKm, HEMKm, HEM2Km, HEM4Km, or HEM13Km was a combined expression from a native genomic copy and a KmYAC-transferred copy (Figure 4C, gray columns). Similar to the situation in KMLSC4D, the relative mRNA levels of GLTXKm, HEM1Km, HEM2Km, HEM4Km, and HEM13Km in KMLSC45D were 12, 21, 25, 6, and 26 times higher respectively than the corresponding genes in KMLSHW (Figure 4A, black column). In KMLSC4D, the heterologous HEMASty, HEMLEco, HEM3sc, HEM12Sc, HEM14Sc, and HEM15Sc introduced by KmYAC were all actively expressed (Figure 4C, white columns). However, the expressions of endogenous HEM3Km, HEM14Km, and HEM15Km were not significantly different from those in KMLSHW, while the expression of HEM12Km in KMLSC45D was significantly reduced compared to that in KMLSHW (Figures 4A, C, black columns). The result suggests that the transfer of heterologous genes in C5 pathways might interfere with the co-upregulation of orthologous genes in the downstream pathway.
There was no difference in the expressions of HEM3sc, HEM12Sc, HEM14Sc, and HEM15Sc in KMLSC4D and KMLSC45D (Figures 4B, C, white columns), indicating the expressions of S. cerevisiae genes were stable in different environments. However, the expressions of K. marxianus genes in the downstream pathway were significantly higher in KMLSC4D compared to the KMLSC45D strain (Figures 4B, C, gray and black columns).
As shown in Figure 5, the total contents of heme and free porphyrin in KMLSC4D and KMLSC45D were 7.57 times and 5.34 times that of KMLSHW, respectively (Figure 5A). pLBA were transformed into KMLSHW, KMLSC4D and KMLSC45D, resulting in KMLSHW-LBA, KMLSC4D-LBA and KMLSC45D-LBA, respectively. The levels of LBA produced in KMLSC4D-LBA and KMLSC45D-LBA increased by 113% and 76.9% respectively, compared to that in KMLSHW-LBA (Figure 5B). Compared to KMLSC45D, higher heme accumulation and LBA production in KMLSC4D might be attributed to the relatively high expressions of K. marxianus genes in the downstream pathway.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g005.jpg)
FIGURE 5. The heme content and LBA production in strains containing designed heme synthesis modules and deletion of inhibitory genes. (A) The relative contents of heme and free porphyrin in KMLSC4D and KMLSC45D. Cells were cultured in YD medium for 72 h. The heme content in KMLSHW was designated as 1. The values represent the mean ± S.D. (n = 3). **, p < 0.01, ***, p < 0.001. (B) The relative levels of LBA produced in KMLSC4D and KMLSC45D. pLBA was transformed into KMLSHW, KMLSC4D and KMLSC45D, resulting in KMLSHW-LBA, KMLSC4D-LBA and KMLSC45D-LBA, respectively. The transformants were then cultured in YD medium for 72 h. The level of LBA produced in KMLSHW-LBA was designated as 1.
Optimization of the medium recipe and feed-batch fermentation for improved yield of LBA
The improved concentration of cellular heme and the yield of LBA were dependent on the presence of KmYAC [C4+Down] and pLBA. Therefore, the stabilities of KmYAC [C4+Down] and pLBA in KMLSC4D-LBA became a concern during industrial production and were thus investigated. KMLSC4D-LBA was cultured in YD medium without selective pressure for 70 generations. After 70 generations, 84.09% of cells contained KmYAC [C4+Down], and 82.46% of cells contained both KmYAC [C4+Down] and pLBA (Figure 6). The result indicate the robust and stable propagation of both KmYAC [C4+Down] and pLBA during mitosis, which will facilitate their applications in industrial production.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g006.jpg)
FIGURE 6. Stability of KmYAC and LBA-expressing plasmid in KMLSC4D. KMLSC4D-LBA were grown for 70 generations without selective pressure. The percentage of cells containing KmYAC [C4+Down] and pLBA was calculated after growing for 7, 14, 21, 42, 56, and 70 generations. The values represent the mean ± S.D. (n = 3). *, p < 0.05. **, p < 0.01. ns, not significant.
To further enhance LBA yield, we conducted an orthogonal experiment to optimize the recipe of synthetic medium (SM). We varied three concentrations each of glucose (2%, 4%, and 6%), glycine (0.1 g/L, 10 g/L, and 20 g/L), and FeSO4·7H2O (20 μM, 100 μM, 200 μM), resulting in nine distinct media. KMLSC4D-LBA was cultured in these media, and titers of LBA were quantified by the grayscale intensities of the bands in SDS-PAGE (Figure 7A). The analysis of variance revealed that the concentrations of glucose, glycine, and FeSO4·7H2O significantly influenced LBA production. According to the range analysis, the impact of these factors on LBA titers was determined as follows: glycine concentration > glucose concentration > iron ion concentration. The recipe that achieved the highest yield of LBA was the SM supplemented with 0.1 g/L of glycine, 6% glucose, and 100 µM FeSO4·7H2O.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1329016/fbioe-11-1329016-HTML/image_m/fbioe-11-1329016-g007.jpg)
FIGURE 7. Optimization of the medium recipe and feed-batch fermentation for improved yield of LBA. (A) Optimization of the medium recipe to improve the yield of LBA. KMLSC4D-LBA was grown in a medium with the indicated concentrations of glucose, glycine, and FeSO4·7H2O (Fe2+) for 72 h. The cell lysate was subjected to SDS-PAGE. The bands of LBA were scanned in grayscale, and the intensities of the bands were quantified. The values represent the mean ± S.D. (n = 3). (B) Growth curves of KMLSC4D-LBA and T1a-LBA in a 5L fermentor. pLBA was transformed into T1a, resulting in T1a-LBA, which served as a control. Cells were grown in the optimized medium in a feed-batch 5L fermentor for 72 h. (C,D) Comparison of heme concentrations in KMLSC4D-LBA and T1a-LBA during feed-batch fermentation. Cells were collected at the indicated times during feed-batch fermentation described in (B). The concentrations of heme in the cell lysate at each time point were quantified by HPLC and compared (C). A chromatogram map of 36-h samples was shown in (D). (E) Image of cells in the fermentor after 36 h. (F) SDS-PAGE of cell lysate of KMLSC4D-LBA and T1a-LBA during feed-batch fermentation. Cells were collected at the indicated times during feed-batch fermentation described in (B). The cell lysate was diluted and subjected to SDS-PAGE. The amount of LBA was semi-quantified using a β-lactoglobulin standard. (G) Comparison of LBA titer in KMLSC4D-LBA and T1a-LBA during feed-batch fermentation. The values represent the mean ± S.D. (n = 3). ***, p < 0.001.
The optimized recipe was applied in the SM medium for the feed-batch fermentation of KMLSC4D-LBA. As a control, LBA was transformed into T1a, resulting in T1a-LBA. The fermentation of both strains was carried out in a 5L fermentor at 32°C for 72 h. The growth rate of KMLSC4D-LBA was slightly lower than that of T1a-LBA during the exponential stage, but the biomass of both strains reached the same level after 72 h (Figure 7B). The heme content in the cell lysate of KMLSC4D-LBA and T1a-LBA were analyzed by quantitative HPLC and compared. The highest heme contents were observed after 36 h, with the content in KMLSC4D-LBA reaching 66.32 mg/L, which was 23.74 times higher than that in T1a-LBA (Figures 7C, D). The higher heme content in KMLSC4D-LBA resulted in a noticeably redder culture compared to that of T1a-LBA (Figure 7E). In KMLSC4D-LBA, accumulated heme might exceed its consumption by LBA and other cellular processes. Excess heme could potentially induce toxicity within the cells (Anzaldi and Skaar, 2010), leading to cell lysis or heme degradation, subsequently causing a noticeable decrease in heme titers after 36 h (Figure 7C). The titer of LBA in KMLSC4D-LBA was quantified by a lactoglobulin standard (Figure 7F). Throughout the fermentation, the titer of LBA in KMLSC4D-LBA consistently surpassed that in T1a-LBA at every time point (Figure 7G). Additionally, the titer of LBA in KMLSC4D-LBA exhibited gradual increments during fermentation. In contrast, the yield of LBA in T1a-LBA declined, possibly due to the instability of LBA resulting from insufficient heme available for binding. After 72 h, the titer of LBA in KMLSC4D-LBA reached 7.27 g/L, which was 80% higher than the titer in T1a-LBA (4.04 g/L) and represented the highest titer of intracellular leghemoglobin reported to date (Figure 7G).
Discussion
Microbial cell factories offer a sustainable, high-yield, and cost-effective source of leghemoglobin to enhance the sensory qualities of plant-based meat. Secretory expression of leghemoglobin simplifies purification procedures and is favored in industrial production. Currently, the highest secretory expression of LBC2 has been achieved in P. pastoris, with a titer of 3.5 g/L, making P. pastoris a favorable host for leghemoglobin production. However, while some proteins expressed by P. pastoris are considered food-safe, P. pastoris itself has not yet obtained Generally Recognized as Safe (GRAS) status (Barone et al., 2023). During the construction of P. pastoris strains for leghemoglobin production, the high-copy integration of the target gene into the genome is achieved through the selection of the G418 antibiotic marker. Within a certain range, the yield of the recombinant protein is positively correlated with the degree of G418 resistance (Lin-Cereghino et al., 2008). Embedding antibiotic resistance of P. pastoris strains might bring potential risks (Baghban et al., 2019). Additionally, P. pastoris is a methylotrophic yeast, and methanol induction is required to induce the expression of recombinant proteins, raising safety concerns related to the final product. Studies have reported no allergenicity or toxicity associated with the recombinant LBC2 expressed in P. pastoris, as well as with the 17 P. pastoris host proteins co-purified with LBC2 (Jin et al., 2018). Nevertheless, the presence of host proteins from a species lacking a GRAS certificate still raises certain safety concerns.
Four types of leghemoglobin are induced at different stages in root nodules. The predominant type of leghemoglobin expressed by P. pastoris is LBC2. In this study, we chose to work with LBA because it lacks predicted glycosylation sites (N-X-S/T) and has demonstrated a relatively higher yield in K. marxianus. Considering that K. marxianus has attained GRAS status (Leonel et al., 2021), there are fewer safety concerns associated with the LBA product produced by K. marxianus. We investigated both intracellular and secretory expressions of LBA in K. marxianus. While achieving a high yield of intracellular expression of LBA in this study, we observed a low yield in secretory expression. Currently, we are optimizing the secretory expression method.
K. marxianus is a crabtree-negative strain known for its rapid growth and high glucose consumption under aerobic conditions. The results of the orthogonal experiment indicate glucose is a pivotal factor influencing LBA production. Consequently, during the feed-batch fermentation, the glucose feeding rate was controlled based on dissolved oxygen levels, typically maintained within the range of 5%–10%. This approach ensured both rapid cell growth and efficient LBA production. The highest LBA production was achieved in 72 h of fermentation, ensuring optimal efficiency for industrial-scale production.
The heme synthesis pathway of KMLSC4D was remodeled by deleting LSC1 and SSN3, which inhibited heme synthesis, and by introducing a designed C4+down module through KmYAC. The expression of native heme synthesis genes was significantly upregulated in KMLSC4D (Figure 4B, black columns). The introduction of additional heme synthesis modules might accelerate the consumption of intermediate metabolites of heme and promote the cycling of the heme metabolism. As a result, the total heme and porphyrin content in KMLSC4D was 6.57 times higher compared to the control. In comparison with KMLSC4D, KMLSC45D contained additional heme synthesis genes from the C5 pathway. However, the enhancement in total heme and porphyrin content in KMLSC45D was reduced to 5.34 times that of the control. The expressions of S. cerevisiae heme synthesis genes did not differ significantly between KMLSC4D and KMLSC45D, while the expressions of K. marxianus heme synthesis genes in KMLSC4D were significantly higher than those in KMLSC45D. These results suggest that heterologous modules transferred by KmYAC, represented by S. cerevisiae heme synthesis genes, exhibit a certain degree of orthogonality and are more robust in countering disturbances caused by the introduction of additional genes. Therefore, the introduction of heterologous modules, rather than native K. marxianus modules, tends to result in the expected and orthogonal expressions in K. marxianus, thereby facilitating the engineering of strains for improved yield.
In summary, in food-grade K. marxianus, we intentionally deleted genes that inhibit heme synthesis and strategically combined them with heme synthesis modules carried by KmYACs. Compared to the wild-type strain, the engineered strains achieved enhanced heme content, increasing by 23.74 times to reach 66.32 mg/L. LBA produced in the engineered strain reached 7.2 g/L, marking the highest intracellular expression of recombinant leghemoglobin to date. The utilization of recombinant LBA expressed in K. marxianus as a substitute ingredient for animal meat will enhance food safety and reduce reliance on animal resources, making a significant contribution to sustainability and environmental protection.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
TT: Conceptualization, Data curation, Formal Analysis, Investigation, Writing–original draft, Writing–review and editing. XW: Data curation, Writing–review and editing. PW: Data curation, Methodology, Writing–review and editing. XL: Data curation, Writing–review and editing. QW: Data curation, Writing–review and editing. YL: Data curation, Writing–review and editing. CL: Data curation, Writing–review and editing. JZ: Data curation, Methodology, Writing–review and editing. YY: Conceptualization, Supervision, Writing–review and editing, Writing–original draft. HL: Conceptualization, Funding acquisition, Supervision, Writing–original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Key Research and Development Program of China 2021YFA0910601 and 2021YFC2100203, Shanghai Municipal Education Commission (2021-03-52) and Science and Technology Research Program of Shanghai 19DZ2282100.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2023.1329016/full#supplementary-material
References
Amberg, D. C., Burke, D. J., and Strathern, J. N. (2005). Methods in yeast genetics: a cold spring harbor laboratory course manual. Cold Spring Harbor, New York, United States: Cold Spring Harbor Laboratory Press.
Amillet, J. M., Buisson, N., and Labbe-Bois, R. (1995). Positive and negative elements involved in the differential regulation by heme and oxygen of the HEM13 gene (coproporphyrinogen oxidase) in Saccharomyces cerevisiae. Curr. Genet. 28 (6), 503–511. doi:10.1007/bf00518161
Anzaldi, L. L., and Skaar, E. P. (2010). Overcoming the heme paradox: heme toxicity and tolerance in bacterial pathogens. Infect. Immun. 78 (12), 4977–4989. doi:10.1128/iai.00613-10
Arredondo-Peter, R., Moran, J. F., Sarath, G., Luan, P., and Klucas, R. V. (1997). Molecular cloning of the cowpea leghemoglobin II gene and expression of its cDNA in Escherichia coli. Purification and characterization of the recombinant protein. Plant Physiol. 114 (2), 493–500. doi:10.1104/pp.114.2.493
Baghban, R., Farajnia, S., Rajabibazl, M., Ghasemi, Y., Mafi, A., Hoseinpoor, R., et al. (2019). Yeast expression systems: overview and recent advances. Mol. Biotechnol. 61 (5), 365–384. doi:10.1007/s12033-019-00164-8
Barone, G. D., Emmerstorfer-Augustin, A., Biundo, A., Pisano, I., Coccetti, P., Mapelli, V., et al. (2023). Industrial production of proteins with Pichia pastoris-komagataella phaffii. Biomolecules 13 (3), 441. doi:10.3390/biom13030441
Duan, J., Yang, D., Chen, L., Yu, Y., Zhou, J., and Lu, H. (2019). Efficient production of porcine circovirus virus-like particles using the nonconventional yeast Kluyveromyces marxianus. Appl. Microbiol. Biotechnol. 103 (2), 833–842. doi:10.1007/s00253-018-9487-2
Hoy, J. A., and Hargrove, M. S. (2008). The structure and function of plant hemoglobins. Plant Physiol. Biochem. 46 (3), 371–379. doi:10.1016/j.plaphy.2007.12.016
Iborra, F. (1993). High efficiency transformation of Kluyveromyces marxianus by a replicative plasmid. Curr. Genet. 24 (1-2), 181–183. doi:10.1007/bf00324685
Ishchuk, O. P., Domenzain, I., Sánchez, B. J., Muñiz-Paredes, F., Martínez, J. L., Nielsen, J., et al. (2022). Genome-scale modeling drives 70-fold improvement of intracellular heme production in Saccharomyces cerevisiae. Proc. Natl. Acad. Sci. 119 (30), e2108245119. doi:10.1073/pnas.2108245119
Ishchuk, O. P., Frost, A. T., Muñiz-Paredes, F., Matsumoto, S., Laforge, N., Eriksson, N. L., et al. (2021). Improved production of human hemoglobin in yeast by engineering hemoglobin degradation. Metab. Eng. 66, 259–267. doi:10.1016/j.ymben.2021.05.002
Ismail, I., Hwang, Y. H., and Joo, S. T. (2020). Meat analog as future food: a review. J. Anim. Sci. Technol. 62 (2), 111–120. doi:10.5187/jast.2020.62.2.111
Jiang, S., Jardinaud, M. F., Gao, J., Pecrix, Y., Wen, J., Mysore, K., et al. (2021). NIN-like protein transcription factors regulate leghemoglobin genes in legume nodules. Science 374 (6567), 625–628. doi:10.1126/science.abg5945
Jin, Y., He, X., Andoh-Kumi, K., Fraser, R. Z., Lu, M., and Goodman, R. E. (2018). Evaluating potential risks of food allergy and toxicity of soy leghemoglobin expressed in Pichia pastoris. Mol. Nutr. Food Res. 62 (1), 1700297. doi:10.1002/mnfr.201700297
Kastaniotis, A. J., and Zitomer, R. S. (2000). Rox1 mediated repression. Oxygen dependent repression in yeast. Adv. Exp. Med. Biol. 475, 185–195. doi:10.1007/0-306-46825-5_18
Keng, T. (1992). HAP1 and ROX1 form a regulatory pathway in the repression of HEM13 transcription in Saccharomyces cerevisiae. Mol. Cell Biol. 12 (6), 2616–2623. doi:10.1128/mcb.12.6.2616-2623.1992
Ko, Y. J., Kim, M., You, S. K., Shin, S. K., Chang, J., Choi, H. J., et al. (2021). Animal-free heme production for artificial meat in Corynebacterium glutamicum via systems metabolic and membrane engineering. Metab. Eng. 66, 217–228. doi:10.1016/j.ymben.2021.04.013
Komar, A. A., Kommer, A., Krasheninnikov, I. A., and Spirin, A. S. (1997). Cotranslational folding of globin. J. Biol. Chem. 272 (16), 10646–10651. doi:10.1074/jbc.272.16.10646
Kornitzer, D., and Roy, U. (2020). Pathways of heme utilization in fungi. Biochimica Biophysica Acta (BBA) - Mol. Cell Res. 1867 (11), 118817. doi:10.1016/j.bbamcr.2020.118817
Kushnirov, V. V. (2000). Rapid and reliable protein extraction from yeast. Yeast 16 (9), 857–860. doi:10.1002/1097-0061(20000630)16:9<857::AID-YEA561>3.0.CO;2-B
Lane, M. M., and Morrissey, J. P. (2010). Kluyveromyces marxianus: a yeast emerging from its sister's shadow. Fungal Biol. Rev. 24 (1), 17–26. doi:10.1016/j.fbr.2010.01.001
Leonel, L. V., Arruda, P. V., Chandel, A. K., Felipe, M. G. A., and Sene, L. (2021). Kluyveromyces marxianus: a potential biocatalyst of renewable chemicals and lignocellulosic ethanol production. Crit. Rev. Biotechnol. 41 (8), 1131–1152. doi:10.1080/07388551.2021.1917505
Lin-Cereghino, J., Hashimoto, M. D., Moy, A., Castelo, J., Orazem, C. C., Kuo, P., et al. (2008). Direct selection of Pichia pastoris expression strains using new G418 resistance vectors. Yeast 25 (4), 293–299. doi:10.1002/yea.1587
Liu, L., Martínez, J. L., Liu, Z., Petranovic, D., and Nielsen, J. (2014). Balanced globin protein expression and heme biosynthesis improve production of human hemoglobin in Saccharomyces cerevisiae. Metab. Eng. 21, 9–16. doi:10.1016/j.ymben.2013.10.010
Liu, Y., Mo, W. J., Shi, T. F., Wang, M. Z., Zhou, J. G., Yu, Y., et al. (2018). Mutational Mtc6p attenuates autophagy and improves secretory expression of heterologous proteins in Kluyveromyces marxianus. Microb. Cell Fact. 17 (1), 144. doi:10.1186/s12934-018-0993-9
Martínez, J. L., Liu, L., Petranovic, D., and Nielsen, J. (2015). Engineering the oxygen sensing regulation results in an enhanced recombinant human hemoglobin production by Saccharomyces cerevisiae. Biotechnol. Bioeng. 112 (1), 181–188. doi:10.1002/bit.25347
Natarajan, C., Jiang, X., Fago, A., Weber, R. E., Moriyama, H., and Storz, J. F. (2011). Expression and purification of recombinant hemoglobin in Escherichia coli. PLoS One 6 (5), e20176. doi:10.1371/journal.pone.0020176
Raimondi, S., Zanni, E., Amaretti, A., Palleschi, C., Uccelletti, D., and Rossi, M. (2013). Thermal adaptability of Kluyveromyces marxianus in recombinant protein production. Microb. Cell Fact. 12, 34. doi:10.1186/1475-2859-12-34
Shao, Y., Xue, C., Liu, W., Zuo, S., Wei, P., Huang, L., et al. (2022). High-level secretory production of leghemoglobin in Pichia pastoris through enhanced globin expression and heme biosynthesis. Bioresour. Technol. 363, 127884. doi:10.1016/j.biortech.2022.127884
Shi, T., Zhou, J., Xue, A., Lu, H., He, Y., and Yu, Y. (2021). Characterization and modulation of endoplasmic reticulum stress response target genes in Kluyveromyces marxianus to improve secretory expressions of heterologous proteins. Biotechnol. Biofuels 14 (1), 236. doi:10.1186/s13068-021-02086-7
Singh, S., and Varma, A. (2017). “Structure, function, and estimation of leghemoglobin,” in Rhizobium biology and Biotechnology. Editors A. P. Hansen, D. K. Choudhary, P. K. Agrawal, and A. Varma (Cham, Germany: Springer International Publishing), 309–330.
Su, H., Chen, X., Chen, S., Guo, M., and Liu, H. (2023). Applications of the whole-cell system in the efficient biosynthesis of heme. Int. J. Mol. Sci. 24 (9), 8384. doi:10.3390/ijms24098384
Wang, M., Shi, Z., Gao, N., Zhou, Y., Ni, X., Chen, J., et al. (2023). Sustainable and high-level microbial production of plant hemoglobin in Corynebacterium glutamicum. Biotechnol. Biofuels Bioprod. 16 (1), 80. doi:10.1186/s13068-023-02337-9
Weickert, M. J., and Curry, S. R. (1997). Turnover of recombinant human hemoglobin in Escherichia coli occurs rapidly for insoluble and slowly for soluble globin. Arch. Biochem. Biophys. 348 (2), 337–346. doi:10.1006/abbi.1997.0410
Wu, P., Mo, W., Tian, T., Song, K., Lyu, Y., Ren, H., et al. (2023). Transfer of disulfide bond formation modules via yeast artificial chromosomes promotes the expression of heterologous proteins in Kluyveromyces marxianus. Biorxiv. doi:10.1101/2023.11.30.569359
Wu, P., Zhou, J., Yu, Y., and Lu, H. (2022). Characterization of essential elements for improved episomal expressions in Kluyveromyces marxianus. Biotechnol. J. 17 (4), e2100382. doi:10.1002/biot.202100382
Xue, J., Zhou, J., Li, J., Du, G., Chen, J., Wang, M., et al. (2023). Systematic engineering of Saccharomyces cerevisiae for efficient synthesis of hemoglobins and myoglobins. Bioresour. Technol. 370, 128556. doi:10.1016/j.biortech.2022.128556
Yang, P., Liu, W., Cheng, X., Wang, J., Wang, Q., and Qi, Q. (2016). A new strategy for production of 5-aminolevulinic acid in recombinant corynebacterium glutamicum with high yield. Appl. Environ. Microbiol. 82 (9), 2709–2717. doi:10.1128/aem.00224-16
Zhao, X. R., Choi, K. R., and Sang, Y. (2018). Metabolic engineering of Escherichia coli for secretory production of free haem. Nat. Catal. 1 (9), 720–728. doi:10.1038/s41929-018-0126-1
Keywords: Kluyveromyces marxianus, leghemoglobin, heme biosynthesis, yeast artificial chromosomes, high-level expression
Citation: Tian T, Wu X, Wu P, Lu X, Wang Q, Lin Y, Liu C, Zhou J, Yu Y and Lu H (2024) High-level expression of leghemoglobin in Kluyveromyces marxianus by remodeling the heme metabolism pathway. Front. Bioeng. Biotechnol. 11:1329016. doi: 10.3389/fbioe.2023.1329016
Received: 27 October 2023; Accepted: 19 December 2023;
Published: 08 January 2024.
Edited by:
Haifeng Zhao, South China University of Technology, ChinaReviewed by:
Liu-Jing Wei, East China University of Science and Technology, ChinaRodrigo Ledesma-Amaro, Imperial College London, United Kingdom
Mangyao Lu, Imperial College London, in collaboration with reviewer RL-A
Copyright © 2024 Tian, Wu, Wu, Lu, Wang, Lin, Liu, Zhou, Yu and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yao Yu, yaoyu@fudan.edu.cn; Hong Lu, honglv@fudan.edu.cn