Pichia pastoris is a powerful and broadly used host for recombinant protein production (RPP), where past bioprocess performance has often been directed with the methanol regulated AOX1 promoter (PAOX1), and the constitutive GAP promoter (PGAP). Since promoters play a crucial role in an expression system and the bioprocess efficiency, innovative alternatives are constantly developed and implemented. Here, a thorough comparative kinetic characterization of two expression systems based on the commercial PDF and UPP promoters (PPDF, PUPP) was first conducted in chemostat cultures. Most promising conditions were subsequently tested in fed-batch cultivations. These new alternatives were compared with the classical strong promoter PGAP, using the Candida antarctica lipase B (CalB) as model protein for expression system performance.
Both the PPDF and PUPP-based expression systems outperformed similar PGAP-based expression in chemostat cultivations, reaching ninefold higher specific production rates (qp). CALB transcription levels were drastically higher when employing the novel expression systems. This higher expression was also correlated with a marked upregulation of unfolded protein response (UPR) related genes, likely from an increased protein burden in the endoplasmic reticulum (ER). Based on the chemostat results obtained, best culture strategies for both PPDF and PUPP expression systems were also successfully implemented in 15 L fed-batch cultivations where qp and product to biomass yield (YP/X*) values were similar than those obtained in chemostat cultivations.
As an outcome of the macrokinetic characterization presented, the novel PPDF and PUPP were observed to offer much higher efficiency for CalB production than the widely used PGAP-based methanol-free alternative. Thus, both systems arise as highly productive alternatives for P. pastoris-based RPP bioprocesses. Furthermore, the different expression regulation patterns observed indicate the level of gene expression can be adjusted, or tuned, which is interesting when using Pichia pastoris as a cell factory for different products of interest.
The non-conventional yeast Komagataella phaffii, widely known under the former name Pichia pastoris, is a distinguished host for recombinant protein production (RPP) [1,2,3,4,5,6,7] and metabolite production . Among the many positive features that make P. pastoris a good choice for RPP, and historically one of the most relevant, is strong and tightly regulated expression based on the alcohol oxidase 1 promoter (PAOX1) [9,10,11,12]. When using the PAOX1 promoter, induction occurs in the presence of methanol, whereas glycerol or glucose fully repress expression . De-repression is not sufficient for significant gene expression; therefore, a simple recombinant protein production process is typically divided into two phases. First, a glucose/glycerol-based batch phase, where a relatively high amount of biomass is generated without recombinant protein production. Subsequently, the methanol feeding phase triggers strong PAOX1-driven protein production. However, such tightly controlled induction and strong expression levels by PAOX1 causes operational drawbacks due to the use of methanol, including high oxygen requirements and heat production, as well as increased costs derived from methanol storage and handling [14, 15]. To address these challenges, mutated promoter variants  or co-substrate feeding strategies had been employed .
In order to open new opportunities, innovative alternatives are constantly developed, evaluated and implemented. In terms of promoter strength, the other methanol inducible promoters such as DAS1 and DAS2 (PDAS1; PDAS2), demonstrate similar strength . In addition to the numerous attempts that had been made to modify PAOX1 regulation by mutagenesis or synthetic fusions , the co-expression of transcription factors was demonstrated as an interesting alternative to induction by methanol. In addition, numerous methanol-independent expression systems have been developed and tested with promising results such as PGTH1, PCAT1, PTHI11 PHpFMD or PTEF, among others [20,21,22,23,24,25].
Historically, RPP improvements have been mainly obtained through strain and promoter system engineering [16, 26,27,28,29]. Multiple clones with different expression cassettes or random integration variants with the same expression cassette are tested and compared in parallel in shake flasks or microtiter plates. This approach is considered fast and cost-effective; however, most testing platforms do not allow control of key bioprocess parameters such as dissolved oxygen, pH, as well as growth and feed rates. Since these parameters affect target protein expression, selection of a “best performing” clone might not always be optimal. Accordingly, the performance of the production clones candidates should be compared in cultivation platforms such as bench-top bioreactors [22, 30, 31], and/or alternative systems that allow controlled substrate delivery. Using bioreactors, production processes can be carried out applying optimal ranges of the key bioprocess parameters [32, 33]. Chemostat systems, where cultures are maintained at non-dynamic, steady-state conditions, have become a valuable tool for bioprocess characterization and further optimization . In this way, a full kinetic characterization of the candidate cell factories can be performed. Furthermore, interestingly, sampling for ‘omics’ analyses can be reliably carried out on cells from chemostat, where cultures have constant key process parameters, and the cell population is highly homogeneous .
Studies including precise strain characterization by chemostat cultivations, have revealed how the specific growth rate (µ) significantly affects recombinant protein production (RPP) rates [6, 15, 35,36,37,38,39,40]. Importantly, the relationship between µ and the specific production rate (qp), also called production kinetics, is dependent on both the expression system used and the recombinant protein expressed. In previous studies, García-Ortega et al. , and Nieto-Taype et al.  described a linear µ-qp relationship when producing a human Fab, as well as the Candida rugose lipase 1 (Crl1), respectively, both under the control of the constitutive PGAP. The same trend was observed for the production of Lipase B from Candida antarctica using the constitutive PGK promoter . These authors concluded that since the constitutive PGAP has a pivotal role in the growth-associated glycolysis, therefore, one should expect the RPP to be growth-coupled. On the other hand, curved/non-linear µ-qp trends were observed [21, 39, 42], suggesting non-coupled transcriptional regulation, or bottlenecks in the protein processing pathway. In particular, Garrigós-Martínez et al.  remarked that the µ-qp bell-shaped trend observed in the PAOX1 regulated production of Crl1 was probably caused by an alternative transcriptional regulation. This conclusion was based on the determination that at different µ, target protein production profiles and the relative transcripts did not present the usual linearity of growth-coupled expression systems.
In this work, the performance characterization of two novel expression systems for RPP with P. pastoris are based on: (1) The new PDF promoter (PPDF, a commercial variant of the Hansenula polymorpha FMD promoter [25, 43], which drives strong transcription by simple methanol-free de-repression and can be also further induced with methanol), and (2) UPP promoter (PUPP, a constitutive commercial variant of a Pichia promoter called GCW14, ). Both promoter systems have been thoroughly studied and compared with PGAP, the most well characterized constitutive promoter, considered a reference standard for methanol-free expression systems. Expression strains for the lipase B from Candida antarctica (CalB) were constructed, all with the same parental strain, gene dosage and identical vectors, except for the promoter sequence driving CALB gene expression. To compare these expression systems, a set of chemostat cultivations designed to assess the effect of different µ values on the production kinetics was performed. Furthermore, CALB transcript levels were determined and compared to the expression levels for each condition tested in chemostat. Finally, working at the µ ranges that generated the best results in chemostat mode, the same selected clones were cultivated in 15 L fed-batch processes to evaluate their performance in this operational mode.
Results and discussion
Strain generation, screening and gene dosage
Isogenic clones were generated to compare the performance of the promoters PUPP, PPDF and PGAP for the CalB expression as a model recombinant protein. Considering the potential clonal variability often observed in Pichia clone generation methods, care was taken to select a clone for each expression system with a single expression cassette integrated into the genome [29, 44]. Subsequently, around 90 individual transformants were analyzed in a high-throughput screen based on deep well plate (DWP) system, to develop a “landscape” of expression data for clone characterization according to Weis et al. . Putative single-copy integration transformants for each of the different promoter constructs were picked from the majority of transformants which showed very similar lipase activity in the supernatant after cultivation and induction in 96-DWP and an initial screen measuring CalB activity of secreted reporter enzyme. Among the discarded clones, secreting higher amounts of CalB were suspected to be associated with multicopy or random integration events; while lower activity observed in some clones might be from with detrimental effects exerted by ectopic integration . Subsequently, a second round of DWP screening among the potential single copy integration candidates was performed, which were tested at least the biological triplicates (data presented on Additional file 1: Figure S1A, B, C). Therefore, for the candidate clones, gene dosage was determined by droplet digital PCR (ddPCR, data presented on Additional file 2: Tables S2A, B) to confirm each construct only contained a single-copy of the respective expression vector in the Pichia genome. Confirmed single copy clones for each expression system presenting an average CalB production were therefore selected to start the expression systems characterization, thus ensuring that production differences are only a result of the effect of each promoter’s specific influence on CalB recombinant expression, and not due to a different gene dosage.
Physiological state comparison of the P. pastoris clones harboring different expression system
Chemostat cultivations were performed with one selected CalB production clone for each of the three different expression systems (GAP-C, PDF-C and UPP-C). This comparison test was performed at three different dilution rates (D): 0.05 h−1, 0.10 h−1 and 0.15 h−1. This characterization allowed to determine the range of dilution rate to significantly improve CalB production in subsequent fed-batch (FB) cultivations.
PUPP and PPDF clones had significantly higher expression levels than those based on PGAP. Furthermore, high levels of recombinant protein expression have been shown to cause a burden on the protein secretion machinery likely due to an overload of the processing capacity [47, 48]. Therefore, an impact of the three expression systems on the physiological state was tested in chemostat cultivations by analyzing glycerol and O2 consumption rates, and CO2 production rates (Fig. 1).
In this regard, no significant differences were evident across the D tested. The specific glycerol consumption rate (qs) and overall biomass-to-substrate yield (YX/S*) were rather similar. As expected, qs increases linearly over the D, whereas YX/S* values were constant (only slight differences could be observed at the highest D). All the clones presented similar specific CO2 production rates (qCO2) and specific O2 consumption rates (qO2) and followed standard linear trends. Consequently, similar respiratory quotient values were exhibited by all the clones studied (RQ, i.e. about 0.62, see Fig. 1b). Based on the analysis at macrokinetic level, it can be stated that the higher CalB production provided by the new generation expression systems based on the promoters PPDF or PUPP did not alter any of the studied physiological parameters compared to the GAP-C, which presents lower CALB expression levels.
Novel expression systems outperform PGAP-based CalB production
Compared to PGAP, both of the new promoters resulted in notably higher qp values, between 4 and ninefold higher at any D (Table 1). UPP-C also had qp values significantly higher than PDF-C at the lowest and middle D. At the highest D, UPP-C was similar to PDF-C, with only slightly higher qp.
Different production kinetic profiles, qp at different D, were obtained for all three expression systems compared (Fig. 2a). UPP-C presented a bell-shape profile with a maximum at mid D, 0.10 h−1. On the other hand, PDF-C generated a clearly saturated profile. And, GAP-C had a relatively linear pattern, with a slight saturation trend at higher D values. This result differs from other examples reported using PGAP in which qp clearly increases linearly with D [37, 38]. Thus, these results indicate that production kinetics, in most cases, are protein dependent.
Another important key parameter to be considered is the overall product to biomass yield (YP/X*); it determines the overall capacity of cells to produce recombinant protein under certain conditions. The PUPP and PPDF expression systems are similar to other systems where increasing D is detrimental for YP/X* [38, 39], as shown in Fig. 3. As observed with qp, the biggest difference between UPP-C and PDF-C YP/X* values was demonstrated at the lowest D, yet similar at higher D. Importantly, the highest YP/X* values for UPP-C and PDF-C were notably higher than those obtained with the GAP-C (i.e. up to 8.9-fold higher with UPP-C at 0.05 h−1) (Fig. 3).
Based on CalB production-related parameters qp and YP/X*, both novel expression systems should be considered good candidate promoters to produce recombinant proteins. Different conditions could be optimal depending on the objective: to reach higher product titer, or maximal productivity. If the objective is the highest protein titer, the lowest D should be selected [38, 39, 49], especially in the case of UPP-C, where YP/X* reductions with D were more pronounced than with PDF-C. On the other hand, in order to maximize qp, the conditions recommended for high production should be at a D of 0.10 h−1, observed most markedly with PDF-C, where D has a bigger impact (Fig. 2a).
The new promoters enable increased tunability of recombinant protein expression processes in Pichia pastoris
mRNA levels are not always directly correlated with the level of recombinant protein production obtained [22, 38, 39], especially since high CALB transcript levels can cause enhanced cellular stress . In the present work, variable CALB transcript levels were observed in all three expression systems compared: GAP-C, UPP-C, and PDF-C. As shown in Fig. 2b, a linear profile of CALB relative transcript levels (RTL) was observed across D for GAP-C, confirming the widely reported constitutive and growth-coupled regulation pattern with PGAP. For this GAP-C, even both RTL-D and qp-D present similar profiles, a slight saturation trend of qp can be observed at high D (Fig. 2a), likely because the low production rates observed with PGAP are not likely to overload of the processing and secretory capacity. For UPP-C, only slight differences in CALB RTL were observed among the different D tested. Therefore, the regulation of CALB expression under PUPP control should be considered growth independent. Strikingly, the CALB mRNA expression patterns do not correlate with the bell-shaped qp profile described in the previous section (highest at a 0.10 D, Fig. 2a). For the PDF-C, RTL presents a bell-shape trend, while the qp-D profile presents a saturation pattern. Therefore, according to the RTL results (Fig. 2b), the PPDF-based expression system exhibits a growth-rate dependent regulation, which thus can be considered a system with a promising tunable expression pattern.
As presented in Fig. 4, the comparison of CALB expression regulated by PUPP and PPDF, relative to the PGAP, illustrates an interesting contrasting behavior. The weaker, growth-coupled, PGAP-based expression system performs better at higher µ; demonstrated here with the single-copy expression strain. In this case, the highest target transcription levels may result in a recombinant protein “burden” can still be sorted, or processed, properly in the ER. On the other hand, with both novel promoters, CALB transcription levels and specially qp ratios, decrease over D, thus indicating that high specific growth rates are detrimental for these more productive systems. Both UPP-C and PDF-C generated CALB transcripts, as well as secreted protein at significantly higher levels than GAP-C. The higher transcript levels may be overwhelming the secretory pathway, triggering the unfolded protein response (UPR). Therefore, UPP-C and PDF-C might be better at low and medium µ, when most of transcription can generate protein, as is demonstrated by the higher productivity rates. Consequently, suitable conditions to balance both growth and protein production are needed to improve production protocols.
UPR influence on CalB production
In order to assess potential endoplasmic reticulum (ER) stress derived from the excessive heterologous protein production, the expression of key UPR genes were analyzed. The reporters selected were two well-known ER-resident chaperones, KAR2 and ERO1, and a gene product generally considered to be an UPR master regulator, HAC1 [2, 51, 52]. The relative transcription levels of these three UPR related genes were measured in UPP-C and PDF-C and compared to the GAP-C. Data is presented as a log2 fold change relative to the GAP-C levels (Fig. 5).
Overall, the expression of all three of the UPR reporter genes, HAC1, KAR2 and ERO1, had similar expression patterns across all the D rates tested. Their expression was the highest when the clones were cultivated at intermediate D of 0.10 h−1, whereas only moderate expression increases were found at D 0.15 h−1. These results are in-line with other published work reporting that not only recombinant protein production, but also µ, having a significant impact on UPR induction [36, 53]. The cited studies are also similar to the data presented here, showing unremarkable UPR levels at lower µ, whereas µ increases led to UPR upregulation.
Strikingly, the growth-coupled expression of CALB by GAP-C was the only case in which the CALB RTLs presented a similar pattern to the qp. The UPP-C qp values mimic the UPR sensor gene expression profiles across the D tested, regardless of CALB transcription rates, which were rather similar for all the D rates. Therefore, UPR might have an influence in subsequent steps of CalB processing and secretion. The UPR impact on CALB expression is demonstrated by comparing qp and CALB RTL values at D 0.10 h−1. In the Fig. 5, UPR-associated gene expression is higher for UPP-C (Fig. 5b) than PDF-C (Fig. 5a) and may explain why UPP-C qp is higher than PDF-C qp at this D despite UPP-C presenting 84% less CALB RTL. In this sense, it has been also described that the co-expression of protein disulfide isomerase, which is also upregulated at higher UPR, enhances active lipase production by P. pastoris . Lastly, the comparison between two D conditions for PDF-C in continuous cultivations supports this hypothesis. At both D 0.05 h−1 and D 0.15 h−1, the CALB RTL levels are rather similar. However, the UPR-related gene expression is growth coupled, enhanced at higher D. The increased UPR could be contributing the 50% higher qp observed at the highest D, even though the target gene RTL levels are rather similar.
Together, these analyses indicate that qp for CalB in P. pastoris is influenced by several factors: heterologous gene transcription rates, recombinant protein-associated UPR, and D-associated UPR. Comparing the different alternatives for methanol-free expression presented in this work, it could be observed that the new generation constructs, based on the PDF and UPP promoters, allowed to achieve CALB transcription levels of up to eightfold higher than with PGAP-regulated expression, for all the D tested. However, at the high target protein expression levels, a direct correlation between CALB RTL and qp was not observed. In contrast, higher qp values were usually observed at D conditions with enhanced expression of UPR-related genes, suggesting a relevant impact of UPR on CalB production.
Fed-batch cultures for further scalable bioprocess development
Through chemostat cultivations, a systematic characterization of the three expression systems studied at different dilution rates was carried out, generating information both at macrokinetic, stoichiometric and transcriptional level. From these results, the range of µ that improves CalB production were determined in order to achieve higher product yields and/or productivities. Ideally, the best µ values found in chemostat cultivations should be implemented to fed batch cultures (FB), which is currently the most widely used scalable operational mode for industrial production of recombinant proteins. However, production kinetics may present relevant differences between the different operational modes [38, 39]. Therefore, fed-batch cultivations were also conducted to confirm CalB production kinetics pattern determined in chemostat cultivations. Carbon-limited fed-batch cultivations, the culture strategy usually considered as most efficient with P. pastoris methanol free processes, were performed with the UPP-C and PDF-C to obtain biomass and CalB production profiles (Fig. 6). Accordingly, based on the results obtained in chemostat cultures in which strong expression systems perform better at mid-low µ, the µ of 0.15 h−1 was discarded for further fed-batch implementation.
As expected, biomass production of all the cultures presented the targeted exponential profile, reaching a maximum between 90 and 100 g L−1 of dry cell weight (DCW), which is considered a standard endpoint for Pichia high-cell density fed-batch (Fig. 6a). CalB production, expressed in activity units (kAU), also increased pseudo-exponentially over time (Fig. 6b). Product titers obtained at the lowest µ tested (~ 0.05 h−1) were substantially higher than with the intermediate µ (~ 0.10 h−1), being 38% higher with PDF-C and 67% with UPP-C. Therefore, YP/X* values were also markedly higher at the low µ cultures (Table 1). Comparing product-related parameters between chemostat and FB cultivations, PDF-C performed better in FB mode (24% higher, on average qp, and 41% on average higher YP/X* at the equivalent µ, Table 1). On the other hand, UPP-C presented smaller differences of qp and YP/X* values at low µ, on average below 10%. Strikingly, for the range µ 0.10 h−1, UPP-C presented performance parameters were significantly worse in FB cultures relative to the chemostat cultures (Table 1). With respect to product quality and purity, it is worth to mention that the supernatant impurity fingerprints were similar between all four fed-batch cultivations; all had very low levels of Pichia native proteins secreted to the culture broth. Additionally, as observed in the SDS-PAGE presented, no signs of CalB degradation is observed, indicating protein quality is similar between the different expression systems and fed-batch cultivations (Additional file 3: Figure S3).
At the end of the batch phase, a twofold higher CalB titer was obtained with UPP-C compared to PDF-C. This points to a difference in the expression systems regulation in presence of excess carbon, a situation that only occurs during the batch phase. Whereas, PUPP presents a constitutive and strong expression in excess of C-source, PPDF-regulated expression was repressed and increase under C-limiting conditions, such as at the end of the batch phase, and during the fed-batch phase. Thus, PPDF generates an interesting tunable expression system, allowing simple uncoupling of biomass growth and achieving different transcription levels without altering the carbon source.
In this study focused on developing methanol-free alternatives for RPP with P. pastoris, a macrokinetic characterization of two promising expression systems were conducted in chemostat cultures, PUPP and PPDF, benchmarking their performance with PGAP. In terms of substrate and respiration-related parameters, all three expression systems behaved similarly, suggesting that the potential differences in CalB production does not significantly alter the yeast homeostasis in chemostat cultivations. Overall, the CalB production kinetics with the two novel expression systems generated significantly higher levels of recombinant protein than the reference, GAP-C, up to ninefold higher in terms of qp. The differences in product-related parameters were primarily attributed to the significantly higher CALB transcription levels. Interestingly, under carbon-limiting conditions, the PPDF-based expression system showed a D-dependent tunable expression, while PUPP-regulated expression was more constant, independent of the growth rates tested. Furthermore, an UPR up-regulation was noted most markedly with the UPP-C at D = 0.1 h−1. At this dilution rate, the three UPR reporter genes monitored were at their highest level. Notable is that the highest qp was also at this D condition.
The chemostat results were used to design strategies based on the µ ranges that provide the best production results for its further implementation in fed-batch cultivations. Thus, both novel expression systems, based on PUPP, PPDF, were also tested in this operational mode. The difference in regulation patterns was reproduced in a fed-batch mode as UPP-C had around twofold higher CalB production than PDF-C at the end of batch phase, illustrating a strong constitutive CALB expression under PUPP regulation. On the other hand, highest CALB expression in the PDF-C was obtained under C-limiting conditions, in which the expression is derepressed, thus presenting an interesting tunable expression pattern. Concerning the production kinetics, UPP-C showed much better performance at low µ in fed-batch, as the qp at this µ outperformed the levels obtained at mid µ by 11%. PDF-C expression, on the other hand, was enhanced at mid µ, as qp was significantly higher than under low µ conditions. In all cases, no significant difference in the product quality was observed among the different fed-batch cultivations performed, presenting for all cases low levels of protein impurity fingerprints and no signs of CalB degradation.
This work, testing alternative promoter designs, validate the approach that confirms the transferability from small scale screenings to its characterization in chemostats. Most importantly, outcomes obtained during the characterization should be considered as highly relevant to establish further strategies to be finally industrially implemented at large scale. Furthermore, these production conditions should also be taken into consideration for further improvements of the host at a molecular level, closing the duty cycle of synthetic biology based on iterative Design-Build-Test-Learn steps (DBTL). The successful results obtained in this work, are expected to make this approach transferable to other processes based on other expression systems and also different microbial cell factories towards rationally improve the efficiency of bioprocesses.
Materials and methods
Clone construction, selection and expression testing
The parental strain, P. pastoris BSYBG11(aox1-/MutS), which is a BioGrammatics (Carlsbad, CA) Pichia pastoris BG11 strain, deposited at Bisy in Austria, was transformed with each of the expression vectors, only differing in the respective promoters. This strain is derived from the BioGrammatics K. phaffii (P. pastoris) strain BG10 but with a slower methanol metabolization phenotype, MutS . The transformation method to express Candida antarctica lipase B (CalB) under the regulation of PGAP, PUPP and PPDF has been described elsewhere . Recombinant vectors with the selected promoters were based on pPpT4_Alpha_S vector . To avoid biasing of the results by transformant variability, low amounts of the linearized plasmid DNA (< 1 µg of DNA) were used to avoid multi copy expression cassette integration as described elsewhere [29, 44].
Candidate clones expression screening was done at microscale cultures in deep well plates (DWP´s), as described by Krainer et al.  with minor modifications. During the induction phase, methanol from BMM2 and BMM10 was replaced with glycerol as the carbon source with 1% and 5% glycerol (w/v), respectively. All media were buffered at pH 7.0.
An initial screen in DWP’s allowed to select the producer candidates from the clonal variability after transformation. Furthermore, a second screening employing DWP’s with seven replicates per clone was used to validate previous results with biological replicates towards the selection of producer clones that integrated only one copy of the desired expression cassettes. These final clones were used in the chemostat and fed-batch cultivations.
Gene dosage determination
The gene dosage/copy number was determined for each selected producer clone using droplet digital PCR (ddPCR) as described elsewhere [38, 39]. ddPCR was performed with primers to amplify the CALB gene present in the expression cassette, as well as with primers for amplification of the Actin 1 gene (ACT1) as a reference. The ACT1 gene had been demonstrated before, and thus to be a reference for single copy gene of the haploid P. pastoris. The list of primers used is presented as Additional file 2: Table S2A.
Total RNA extraction, cDNA synthesis and transcript level determination
RNA was isolated from 1 mL culture samples taken from the chemostat culture broth as defined by Landes et al. ; RNA was prepared using the SV Total RNA Isolation System (Promega, Madison, Wisconsin, US) following the manufacturer’s instructions.
RNA integrity was checked by agarose electrophoresis and RNA concentration was measured with Nanodrop 2000 (Thermo Scientifc™, Waltham, MA, US).
cDNA was synthetized using iScript™ cDNA Synthesis kit (Bio-Rad, Hercules, CA, USA), following the manufacturer’s instructions. Primers were designed to analyze the relative transcript levels (RTL) of the following target genes by qPCR: CALB, KAR2 and ERO1 (two ER-resident chaperones), and HAC (UPR master regulator). MTH1 gene was used as housekeeping gene for the transcriptional analysis. The qPCR procedure, including equipment, qPCR master mix solution and housekeeping gene is detailed in prior work . The annealing-extension temperature was adjusted to 59 °C. The list of primers used is presented as Additional file 2: Table S2A.
Chemostat cultivations were performed in duplicate, all in 2 L Biostat B plus Bioreactors (Sartorius Stedim, Goettingen, Germany) according to García Ortega et al. . Batch and chemostat media compositions are stoichiometrically identical to the detailed in the reference ; however, the concentrations were reduced by half.
Cultivation conditions were monitored and controlled at the following set points: pH, 5.0 with addition of 15% (v/v) ammonium hydroxide; temperature, 25 °C; stirring rate, 700 rpm; air flow, 0.8 vvm and pO2 values were variable depending on the dilution rate. pO2 values were above 20% in all conditions tested. An exhaust gas condenser with cooling water at 4 °C minimized mass loses by water evaporation and other possible volatile compounds.
A broad range of dilution rates were covered for the three expression systems tested. Taking into consideration that 0.19 h−1 was the P. pastoris µmax of GAP-C at 25 °C (data not shown), the following dilution rates were used: 0.05 h−1, 0.10 h−1 and 0.15 h−1 (dilution rates were tested as low, middle and high growth rate conditions, respectively). In order to ensure that the steady state was reached, the stability of the parameters of interest were monitored from the third residence time until the fifth one, where samples were taken.
Fed-batch cultivations were performed in duplicate, all in New Brunswick BioFlo 510 bioreactor (Eppendorf, Germany), connected to BioCommand Control software. Batch media composition was the same as that used in the chemostat runs, stoichiometrically identical to that used in prior work , except glycerol, instead of glucose, was used in the fed-batch feeding. The cultures were grown at 25 °C under overpressure (0.2 bar) and had a 7.5 L starting volume, including 1 L seed. The pH was kept at 5.0 by the automated addition of 12.5% NH4OH. Dissolved oxygen (DO) was maintained above 30% of air saturation with the automatic control of stirrer speed (400–700 rpm), constant airflow 10 L min−1, and enriched with O2 when needed (0–50% of total inlet flow rate).
Glycerol feeding was started upon depletion of batch medium glycerol. An exponential pre-programming feeding rate was performed to maintain the specific growth rate (0.05 h−1 or 0.10 h−1) constant at the selected set-point. All cultivations were grown under carbon-limiting conditions. The procedure is described in detail elsewhere .
Biomass concentrations were measured in triplicate as DCW values, as described elsewhere . The relative standard deviation (RSD) of the measurements was about 4%.
Quantification of the carbon source and byproducts
The concentration of both glycerol and the potential fermentation byproducts, were measured by HPLC. The column and the program used are described elsewhere . RSD was invariably less than 1%.
A BlueInOne Cell gas analyzer (BlueSens, Herten, Germany) was used for monitoring the CO2 and O2 molar fraction of the chemostat cultivations off-gas. Off-gas pressure and humidity measurements were used to calculate the oxygen uptake rate (OUR), carbon dioxide evolution rate (CER), and their corresponding specific rates (qO2 and qCO2) and respiratory quotient (RQ), as previously described. RSD was less than 5% in all cases.
CalB activity was determined by an esterase activity assay, similar to the assay described by Krainer et al. . Briefly, 100 μL of culture supernatant was mixed with 900 μL of fresh assay solution containing 4 mM p-NPB in 300 mM Tris–HCl, pH 7.4, 1% acetone. The absorbance increase at 405 nm was monitored at 30 °C for 2 min (Specord 200 Plus spectrophotometer from Analytic Jena Germany). One activity unit was defined as the amount of enzyme needed to release 1 μmol of p-nitrophenol per minute under assay conditions. RSD was less than 4%.
SDS-PAGE analysis were performed with the culture supernatants collected during the bioreactor cultivations, which were diluted into water prior to loading into precast 4–20% kD Criterion TGX Stain-Free Gel (Biorad, Hercules, CA, USA). Fifteen µL of the diluted samples were mixed with 5 µL of 4 × loading buffer (20% glycerol, 4% SDS, 0.3 mM bromophenol blue (Merck), 10% β-mercaptoethanol, 0.1 M Tris, pH 6.8), and incubated at 95 °C for 5 min. Samples were cooled, centrifuged quickly, and 19 µL loaded into 4–20% kD Criterion TGX Stain-Free Gel together with Precision Plus Protein molecular weight marker (Biorad, Hercules, CA, USA). The visualization of the gel was performed using the Gel Doc EZ (Biorad, Hercules, CA, USA).
Process parameter determination, consistency test and data reconciliation
Mass balance and stoichiometric equations
All equations used to calculate yields and rates are based on mass balances at continuous and fed-batch operation and can be found elsewhere [11, 37]. The P. pastoris elemental composition grown on glycerol as the sole C-source was determined as previously reported . Carbon and electron balances were checked and less than 5% of deviation was observed prior to reconciliation.
Consistency test and data reconciliation
Measurement consistency was checked by using the standard test with carbon and electron balances as constraints. Both online and offline measurements enabled the calculation of five key specific rates in the black-box process model: biomass generation (μ), glucose uptake rate (qs), product generation rate (qp), oxygen uptake rate (qO2) and carbon dioxide production rate (qCO2). The method used for this purpose is described in detail elsewhere .
Availability of data and materials
The datasets and materials used and/or analysed during the current study are available from the corresponding author on reasonable request.
Cereghino GPL, Cereghino JL, Ilgen C, Cregg JM. Production of recombinant proteins in fermenter cultures of the yeast Pichia pastoris. Curr Opin Biotech. 2002;13:329–32.
Barrigón JM, Montesinos JL, Valero F. Searching the best operational strategies for Rhizopus oryzae lipase production in Pichia pastoris Mut+ phenotype: Methanol limited or methanol non-limited fed-batch cultures? Biochem Eng J. 2013;75:47–54.
Barrigon JM, Valero F, Montesinos JL. A macrokinetic model-based comparative meta-analysis of recombinant protein production by Pichia pastoris under AOX1 promoter. Biotechnol Bioeng. 2015;112:1132–45.
Ponte X, Montesinos-Seguí JL, Valero F. Bioprocess efficiency in Rhizopus oryzae lipase production by Pichia pastoris under the control of PAOX1 is oxygen tension dependent. Process Biochem. 2016;51:1954–63.
Ponte X, Barrigón JM, Maurer M, Mattanovich D, Valero F, Montesinos-Seguí JL. Towards optimal substrate feeding for heterologous protein production in Pichia pastoris (Komagataella spp) fed-batch processes under PAOX1 control: a modeling aided approach. J Chem Technol Biotechnol. 2018;93:3208–18.
Heyland J, Fu J, Blank LM, Schmid A. Quantitative physiology of Pichia pastoris during glucose-limited high-cell density fed-batch cultivation for recombinant protein production. Biotechnol Bioeng. 2010;107:357–68.
Çalık P, Ata Ö, Güneş H, Massahi A, Boy E, Keskin A, Örtürk S, Zerze GH, Özdamar TH. Recombinant protein production in Pichia pastoris under glyceraldehyde-3-phosphate dehydrogenase promoter: from carbon source metabolism to bioreactor operation parameters. Biochem Eng J. 2015;95:20–36.
Berrios J, Flores MO, Díaz-Barrera A, Altamirano C, Martínez I, Cabrera Z. A comparative study of glycerol and sorbitol as co-substrates in methanol-induced cultures of Pichia pastoris: temperature effect and scale-up simulation. J Ind Microbiol Biotechnol. 2017;44:407–11.
Prielhofer R, Reichinger M, Wagner N, Claes K, Kiziak C, Gasser B, Mattanovich D. Superior protein titers in half the fermentation time: Promoter and process engineering for the glucose-regulated GTH1 promoter of Pichia pastoris. Biotechnol Bioeng. 2018;115:2479–88.
Landes N, Gasser B, Vorauer-Uhl K, Lhota G, Mattanovich D, Maurer M. The vitamin-sensitive promoter PTHI11 enables pre-defined autonomous induction of recombinant protein production in Pichia pastoris. Biotechnol Bioeng. 2016;113:2633–43.
Wang J, Zhang T, Li Y, Li L, Wang Y, Yang B, Wang Y. High-level expression of Thermomyces dupontii thermo-alkaline lipase in Pichia pastoris under the control of different promoters. Biotech. 2019;9:33.
Vogl T, Fischer JE, Hyden P, Wasmayer R, Sturmberger L, Glieder A. Orthologous promoters from related methylotrophic yeasts surpass expression of endogenous promoters of Pichia pastoris. AMB Express. 2020;10:38.
Liu J, Han Q, Cheng Q, Chen Y, Wang R, Li X, Liu Y, Yan D. Efficient expression of human lysozyme through the increased gene dosage and co-expression of transcription factor Hac1p in Pichia pastoris. Curr Microbiol. 2020;77:846–54.
Yang J, Lu Z, Chen J, Chu P, Cheng Q, Liu J, Ming F, Huang C, Xiao A, Cai H, Zhang L. Effect of cooperation of chaperones and gene dosage on the expression of porcine PGLYRP-1 in Pichia pastoris. Appl Microbiol Biotechnol. 2016;100:5453–65.
Vogl T, Sturmberger L, Kickenweiz T, Wasmayer R, Schmid C, Hatzl AM, Gerstmann MA, Pitzer J, Wagner M, Thallinger GG, Geier M, Glieder A. A toolbox of diverse promoters related to methanol utilization: functionally verified parts for heterologous pathway expression in Pichia pastoris. ACS Synth Biol. 2016;5:172–86.
Cámara E, Landes N, Albiol J, Gasser B, Mattanovich D, Ferrer P. Increased dosage of AOX1 promoter-regulated expression cassettes leads to transcription attenuation of the methanol metabolism in Pichia pastoris. Sci Rep. 2017;7:44302.
Cámara E, Monforte S, Albiol J, Ferrer P. Deregulation of methanol metabolism reverts transcriptional limitations of recombinant Pichia pastoris (Komagataella spp) with multiple expression cassettes under control of the AOX1 promoter. Biotechnol Bioeng. 2019;116:1710–20.
Krause M, Ukkonen K, Haataja T, Ruottinen M, Glumoff T, Neubauer A, Neubauer P, Vasala A. A novel fed-batch based cultivation method provides high cell-density and improves yield of soluble recombinant proteins in shaken cultures. Microb Cell Fact. 2010;9:11.
Panula-Perälä J, Šiurkus J, Vasala A, Wilmanowski R, Casteleijn MG, Neubauer P. Enzyme controlled glucose auto-delivery for high cell density cultivations in microplates and shake flasks. Microb Cell Fact. 2008;7:31.
Nieto-Taype MA, Garcia-Ortega X, Albiol J, Montesinos-Seguí JL, Valero F. Continuous cultivation as a tool toward the rational bioprocess development with Pichia Pastoris cell factory. Front Bioeng Biotechnol. 2020;8:632.
Rebnegger C, Graf AB, Valli M, Steiger MG, Gasser B, Maurer M, Mattanovich D. In Pichia pastoris, growth rate regulates protein synthesis and secretion, mating and stress response. Biotechnol J. 2014;9:511–25.
Garcia-Ortega X, Adelantado N, Ferrer P, Montesinos JL, Valero F. A step forward to improve recombinant protein production in Pichia pastoris: from specific growth rate effect on protein secretion to carbon-starving conditions as advanced strategy. Process Biochem. 2016;51:681–91.
Nieto-Taype MA, Garrigós-Martínez J, Sánchez-Farrando M, Valero F, Garcia-Ortega X, Montesinos-Seguí JL. Rationale-based selection of optimal operating strategies and gene dosage impact on recombinant protein production in Komagataella phaffii (Pichia pastoris). Microb Biotechnol. 2020;13:315–27.
García-Ortega X, Cámara E, Ferrer P, Albiol J, Montesinos-Seguí JL, Valero F. Rational development of bioprocess engineering strategies for recombinant protein production in Pichia pastoris (Komagataella phaffii) using the methanol-free GAP promoter Where do we stand? New Biotechnol. 2019;53:24–34.
De Macedo RJ, Garcia-Ortega X, Montesinos-Seguí JL, Guimaraes Freire DM, Valero F. Continuous operation, a realistic alternative to fed-batch fermentation for the production of recombinant lipase B from Candida antarctica under the constitutive promoter PGK in Pichia pastoris. Biochem Eng J. 2019;147:39–47.
Canales C, Altamirano C, Berrios J. Effect of dilution rate and methanol-glycerol mixed feeding on heterologous Rhizopus oryzae lipase production with Pichia pastoris Mut+ phenotype in continuous culture. Biotechnol Prog. 2015;31:707–14.
Delic M, Göngrich R, Mattanovich D, Gasser B. Engineering of protein folding and secretion-strategies to overcome bottlenecks for efficient production of recombinant proteins. Antioxid Redox Signal. 2014;21:414–37.
Garcia-Ortega X, Ferrer P, Montesinos JL, Valero F. Fed-batch operational strategies for recombinant Fab production with Pichia pastoris using the constitutive GAP promoter. Biochem Eng J. 2013;79:172–81.
Zhu T, Guo M, Zhuang Y, Chu J, Zhang S. Understanding the effect of foreign gene dosage on the physiology of Pichia pastoris by transcriptional analysis of key genes. Appl Microbiol Biotechnol. 2011;89:1127–35.
Raschmanová H, Zamora I, Borčinová M, Meier P, Weninger A, Mächler D, Glieder A, Melzoch K, Knejzlík Z, Kovar K. Single-cell approach to monitor the unfolded protein response during biotechnological processes with Pichia pastoris. Front Microbiol. 2019;10:335.
Sha C, Yu XW, Lin NX, Zhang M, Xu Y. Enhancement of lipase r27RCL production in Pichia pastoris by regulating gene dosage and co-expression with chaperone protein disulfide isomerase. Enzyme Microb Technol. 2013;53:438–43.
Lin-Cereghino J, Wong WW, Xiong S, Giang W, Luong LT, Vu J, Johnson SD, Lin-Cereghino GP. Condensed protocol for competent cell preparation and transformation of the methylotrophic yeast Pichia pastoris. Biotechniques. 2005;38:44.
Cos O, Serrano A, Montesinos JL, Ferrer P, Cregg JM, Valero F. Combined effect of the methanol utilization (Mut) phenotype and gene dosage on recombinant protein production in Pichia pastoris fed-batch cultures. J Biotechnol. 2005;116:321–35.
Garcia-Ortega X, Valero F, Montesinos-Seguí JL. Physiological state as transferable operating criterion to improve recombinant protein production in Pichia pastoris through oxygen limitation. J Chem Technol Biotechnol. 2017;92:2573–82.
Krainer FW, Dietzsch C, Hajek T, Herwig C, Spadiut O, Glieder A. Recombinant protein expression in Pichia pastoris strains with an engineered methanol utilization pathway. Microb Cell Fact. 2012;11:22.
Tomàs-Gamisans M, Ferrer P, Albiol J. Fine-tuning the P pastoris iMT1026 genome-scale metabolic model for improved prediction of growth on methanol or glycerol as sole carbon sources. Microb Biotechnol. 2018;11:224–37.
The Spanish group is member of 2017-SGR-1462 and the Reference Network in Biotechnology (XRB, Generalitat de Catalunya). JGM acknowledges a PIF scholarship from the Universitat Autònoma de Barcelona. MAN acknowledges award by the National Council of Science Technology and Technological Innovation (CONCYTEC) through its executing unit and the National Fund for Scientific, Technological and Technological Innovation Development (FONDECYT). This work was performed within the EU project IBISBA1.0 including the two IBISBA partners UAB and VTT in partnership with the bisy GmbH. EU-IBISBA obtained an ESFRI status in 2018, with the aim to accelerate biotechnology and synthetic biology activities in Europe through service infrastructure and know-how. EU-IBISBA is gratefully acknowledged for facilitating part of the work in this study through its trans-national access that received funding from the European Union’s Horizon 2020 research and innovation program (IBISBA 1.0 project) under grant N° 730,976.
This work was funded by the EU-IBISBA through its trans-national access funded by the European Union’s Horizon 2020 research and innovation program (IBISBA 1.0 project) under grant N° 730,976.
Authors and Affiliations
Department of Chemical, Biological and Environmental Engineering, Universitat Autònoma de Barcelona, 08193, Bellaterra (Cerdanyola del Vallès), Spain
Javier Garrigós-Martínez, Miguel Angel Nieto-Taype, José Luis Montesinos-Seguí, Francisco Valero & Xavier Garcia-Ortega
Solutions for Natural Resources and Environment, VTT Technical Research Centre of Finland Ltd, Tietotie 2, Espoo, Finland
JG-M: Chemostat cultivations, tech transfer between UAB and VTT, data treatment for UAB and VTT fermentations, transcriptional analysis, manuscript drafting; KV: design and execution of fed-batch fermentations at VTT, enzymatic assays, interpretation of data and contributed to the drafting of manuscript; MAN-T: Strain construction, clone screening and transcriptional analysis; JT: tech transfer and fed-batch fermentations at VTT, data analysis; JU: participated in the design of experiments and assisted drafting manuscript; PT: Coordination of partners and commented manuscript draft; CS: participated in the strain engineering of CalB producer clones from Bisy; IT: Promoter and strain design and qPCR results interpretation; KM: qPCR results interpretation and manuscript preparation; MP: participated in the interpretation of results and assisted drafting manuscript; JLM-S: Assistance on data treatment, data interpretation and manuscript drafting; FV: Data interpretation and assistance on manuscript drafting; AG: Conceived the study, participated in the interpretation of results and assisted manuscript drafting; XG-O: Conceived the study, participated in the interpretation of results and assisted manuscript drafting. All authors read and approved the final manuscript.
The authors of Bisy GmbH and BioGrammatics declare that they have interests to commercialize the promoters used for this recombinant protein production study. The authors of UAB and VTT have no any competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
S1. Figures that present as landscapes the second round of clone screening results performed for the expression systems studied. Mean values for each and standard deviations are plotted and expressed in terms of activity units (AU) normalized per biomass concentration (OD600). S1A–PGAP-based clone GAP-C; S1B–PPDF-based clone PDF-C; S1C–PUPP-based clone UPP-C.
A, table listing the primer pairs used for gene dosage analyses and relative transcription levels (RTL) determination by means of ddPCR and qPCR, respectively. B, table that presents the gene dosage determination of CalB producer clones by digital droplet PCR (ddPCR). Analyses were performed by triplicates, using Actin gene (ACT1) as reference. Two positive controls which contain 3 and 5 copies of the expression cassette for Candida rugosa lipase 1 (CRL1) were also analyzed as controls.
Figure that presents the analysis of the product quality by SDS-PAGE. Samples from fed-batch cultivations—PDF-C and UPP-C—run at A: µsp = 0.10 h−1 and B: µsp = 0.05 h−1 were analyzed. Different samples obtained at different feeding time (FT) supernatants were loaded on SDS-PAGE. BSA standards at different concentrations were also loaded in SDS-PAGE as reference (lanes 1–4). A. µsp = 0.10 h−1.PDF 7: 17h FT; PDF 8: 18h FT; PDF 9: 20h FT; PDF 10: 21h FT; UPP 8: 19h FT; UPP 9: 21h FT; UPP 10: 22h FT. B. µsp = 0.05 h−1. PDF 1: Batch end, PDF 5: 19h FT; PDF 7: 24h FT; PDF 8: 27h FT; PDF 9: 36h FT; PDF 10: 41h FT; PDF 11: 43h FT; UPP 1: Batch end, UPP 7: 27h FT; UPP 8: 37h FT; UPP 9: 40h FT; UPP 10: 42h FT; UPP 11:44h FT.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Garrigós-Martínez, J., Vuoristo, K., Nieto-Taype, M.A. et al. Bioprocess performance analysis of novel methanol-independent promoters for recombinant protein production with Pichia pastoris.
Microb Cell Fact20, 74 (2021). https://doi.org/10.1186/s12934-021-01564-9