Skip to main content

Synergistic optimisation of expression, folding, and secretion improves E. coli AppA phytase production in Pichia pastoris

Abstract

Background

Pichia pastoris (Komagataella phaffii) is an important platform for heterologous protein production due to its growth to high cell density and outstanding secretory capabilities. Recent developments in synthetic biology have extended the toolbox for genetic engineering of P. pastoris to improve production strains. Yet, overloading the folding and secretion capacity of the cell by over-expression of recombinant proteins is still an issue and rational design of strains is critical to achieve cost-effective industrial manufacture. Several enzymes are commercially produced in P. pastoris, with phytases being one of the biggest on the global market. Phytases are ubiquitously used as a dietary supplement for swine and poultry to increase digestibility of phytic acid, the main form of phosphorous storage in grains.

Results

Potential bottlenecks for expression of E. coli AppA phytase in P. pastoris were explored by applying bidirectional promoters (BDPs) to express AppA together with folding chaperones, disulfide bond isomerases, trafficking proteins and a cytosolic redox metabolism protein. Additionally, transcriptional studies were used to provide insights into the expression profile of BDPs. A flavoprotein encoded by ERV2 that has not been characterised in P. pastoris was used to improve the expression of the phytase, indicating its role as an alternative pathway to ERO1. Subsequent AppA production increased by 2.90-fold compared to the expression from the state of the AOX1 promoter.

Discussion

The microbial production of important industrial enzymes in recombinant systems can be improved by applying newly available molecular tools. Overall, the work presented here on the optimisation of phytase production in P. pastoris contributes to the improved understanding of recombinant protein folding and secretion in this important yeast microbial production host.

Background

Phytases (myo-inositol hexakisphosphate phosphohydrolases) are routinely added to livestock feed for the liberation of phosphate from phytate. The use of phytase reduces the addition of inorganic phosphate in poultry and swine diets and decreases the anti-nutritional effects of phytate [26, 34].

Phytases can be found in animals and plants but microbial phytases have been the most extensively studied and constitute all commercial formulations available in the livestock feed supplement market [16, 38]. Among commercial phytases, the E. coli AppA phytase shows high catalytic efficiency [37]. AppA phytase and derived variants have been produced in Trichoderma reesei, Schizosaccharomyces pombe and P. pastoris [16, 38]. There has been significant commercial work towards process optimisation to increase phytase production yields to achieve technical and economic manufacturing feasibility, but there is limited available information in the academic literature on P. pastoris strain engineering for phytase production. Previous investigations on the production of AppA phytase in P. pastoris include the expression under the strong constitutive promoter of the glyceraldehyde 3-phosphate dehydrogenase gene (PGAP) and methanol inducible alcohol oxidase 1 promoter (PAOX1) promoters [8, 9, 32, 38, 53, 56].

P. pastoris has become an important platform for heterologous protein production due to its growth to high cell density and high secretory capabilities, while secreting low amounts of endogenous proteins thereby simplifying down-stream processing. Significant development of molecular tools, including synthetic promoters for fine-tuning of expression, glyco-engineered strains and CRISPR/Cas9 technology, have driven the generation of highly secreting strains [65, 70].

Common approaches to improve heterologous protein production in P. pastoris include increasing gene copy number, use of strong inducible promoters and co-expression of chaperones to facilitate folding and secretion [13, 46, 63, 65,66,67]. For some enzymes, the presence of multiple disulfide bonds is challenging when striving for high and industrially compatible levels of expression due to the opportunity for incorrect bond formation and misfolding. AppA phytase has four disulfide bonds with one being non-consecutive, which means that the cysteines forming the pair are not consecutive in the protein sequence and have an increased chance of bonding between incorrect cysteine side chains [39]. Co-expressing suitable chaperone/s in a timely manner with the heterologous protein can therefore aid correct disulfide bond formation and rearrangement leading to correctly folded and active enzyme with reduced accumulation of mis-folded proteins.

When the accumulation of misfolded proteins exceeds the ER capabilities, the unfolded protein response (UPR) is activated through the splicing of HAC1 mRNA. HAC1 is the master regulator of the UPR which leads to transcription of a cascade of genes involved in protein folding and ER-associated protein degradation pathway (ERAD) [2].

Formation of disulfide bonds in the ER is possible due to the local oxidising environment. Proteins undergoing folding form disulfide bonds through dithiol-disulfide exchange with the oxidised form of the thioredoxin-like protein disulfide isomerase (PDI) [20]. PDI also catalyses the isomerisation of disulfide bonds [20, 30, 35, 36, 49].

Several engineering strategies can be applied at different points of the pathway to circumvent limitations for efficient recombinant protein production. Co-expression of folding and secretion chaperones has been successful in many cases but optimisation requires the simultaneous expression fine-tuning of multiple proteins to avoid cellular stress and undesirable mutations that reduce production [63, 64, 66, 67]. Co-expression of the gene of interest and folding helper proteins frequently necessitates fine-tuning of the ratio of the proteins to each other, their overall amounts (to avoid overburdening the cellular machinery), and the timing of expression (i.e. it may be beneficial to express a folding helper protein prior to its protein folding target).

Recently developed bidirectional promoters (BDPs) are particularly useful to achieve co-expression of target enzymes and chaperones that assist expression [63, 66, 67]. BDPs can initiate transcription on opposite sides and present diverse expression strengths, ratios, and regulatory profiles useful for optimisation of multiple gene expression. BDPs available for P. pastoris are not limited to methanol induction, but also offer constitute and derepressed regulatory profiles, which can be beneficial for scale up and industrial production. In this work, we used BDPs to co-express several chaperones and to increase production capability of E. coli AppA phytase in P. pastoris, surpassing expression from the commonly used PAOX1 by nearly three-fold.

Methods

Strains and growth conditions

The strain used in this study was the P. pastoris BG11 strain (derivative of P. pastoris BG10 strain, ΔAOX1 (mutS—methanol utilization slow) from ATUM Inc. (Newark, California, USA). α-Select Silver efficiency competent E. coli strain (Bioline, Australia) was used for cloning. Cultivations were conducted in Luria Broth (LB) media for E. coli and yeast cultures were either grown in YPD medium (1% w/v yeast extract, 2% w/v peptone and 2% w/v glucose), buffered minimal dextrose (BMD) medium (1.34% Yeast Nitrogen Base YNB, 4 × 10−5% biotin, 200 mM potassium phosphate buffer pH 6.0 and 2% glucose), buffered minimal methanol (BMM) medium (1.34% YNB, 4 × 10−5% biotin, 200 mM potassium phosphate buffer pH 6.0) with 1% methanol (BMM2) or 5% methanol (BMM10). Antibiotic Zeocin (Invitrogen) was added to the media when required at a final concentration of 25 µg/mL for E. coli or 100 µg/mL for P. pastoris cultivations. Small scale cultivations of P. pastoris were performed using deep-well-plates (DWP) and induction protocols reported previously [69].

Cloning and transformation of Pichia pastoris

The vector used for cloning and integration into P. pastoris genome was pD912 from ATUM Inc., which is based on the pPpT4S vector family [43]. Monodirectional promoter (MDP) sequences PAOX1, glyceraldehydes-3-phosphate dehydrogenase (PGAP), dihydroxyacetone synthase promoter (PDAS1), catalase 1 (PCAT1), and the orthologous formate dehydrogenase promoter from Hansenula polymorpha (PHpFMD) [62] were amplified with 50 bp overlapping sequences at 5′ and 3′ for Gibson assembly cloning into pD912 derived plasmids [11, 27, 29, 68]. BDP sequences PDAS1/2,, PHpFMD-HpMOX and PAOX1-CAT1 were obtained for Bisy e.U. and cloned into pD912 derived plasmids following the strategy outlined by [63, 66, 67]. The AppA E. coli phytase gene (accession number X05471.1) was codon optimised by ATUM Inc. and expressed in the direction of the first promoter region indicated in the BDP name, using the α-mating factor secretion signal from S. cerevisiae [7]. The genes HAC1 (spliced version) [25], PDI (accession number ACF17572.1), ERV2 (accession number ANZ73509.1), MPDI (accession number XP_002489466), ERO1 (accession number ANZ74048.1), KAR2 (accession number ANZ77450.1), GPX1 (accession number BAH57503.1), SEC1 (accession number CAY67514.1), SLY1 (accession number CCA40908.1), EUG1 from S. cerevisiae (accession number XP_446874) were ordered as gBlocks from Integrated DNA Technologies (IDT) with 50 bp overlapping sequences at 5′ and 3′ for Gibson assembly cloning into pD912 derived plasmids. Chaperones were expressed in the direction of the second promoter region indicated in the BDP name (Additional file 1). P. pastoris cells were transformed with SwaI linearised plasmids following a condensed standard protocol [40] with up to 5 µg of total DNA. After the transformation, screening of transformants was performed in DWP as described previously [69].

Phytase expression

AppA phytase expression was conducted in 250 mL baffled shake flasks following standard expression conditions at 28˚C, 250 rpm. The culture was grown in 50 mL of BMD1 for 65 h following methanol induction with BMM10 and consecutive additions of pure methanol 1% final concentration until harvest at 132 h. The protein concentration in the culture supernatant was determined using the Bradford method [5].

Genomic DNA extraction and gene copy number determination

Genomic DNA was extracted using ISOLATE II Genomic DNA kit (Bioline Pty Ltd.; Alexandria, Australia) following the standard protocol with minor modifications as described previously [63, 66, 67]. Gene copy number determination was performed as described previously [6] using Bio-Rad QX200 Droplet Digital PCR (ddPCR). Genomic DNA was digested with SphI restriction enzyme (New England Biolabs) and 0.5 ng of digested DNA of each sample was added to the reaction mixture. For DNA amplification of the AppA phytase gene a set of primers (Fwd 5′-TTGTCCTCAATCCGGTCAAG-3′ and Rev 5′-AGGGTTAAACAACGGATCGG-3′) and a 5′ hexachloro-fluorescein (HEX) dye-labelled probe were designed (5′-HEX-TTGGCTCCAGACTGTGCTATCACTGT- IB®FQ- 3′). The arginosuccinate lyase (ARG4) gene was used as reference gene and a set of primers (Fwd 5′-TGCGGTTGTATGTCAGAGAC-3′ and Rev 5′-GGTTGAGCTCTTTGCAAGTG-3′) and 5′ 6-fluorescein amidite (6-FAM) dye-labelled probe were designed (5′-FAM-TGGCTGACTATCTGAAGCAGTTCATTCA- IB®FQ- 3′) [1]. The probes were also labelled with Iowa Black® FQ at the 3′ terminus. All primers and probes were synthetised by IDT. Each PCR was performed in a 20 µL volume containing 12.5 µL of ddPCR Supermix for Probes (Bio-rad), 900 nM of each primer, 250 nM of each probe, 0.5 ng of digested genomic DNA and the required amount of MilliQ water. Droplets were generated using the Droplet Generator (Bio-Rad) and transferred into a 96-well plate for amplification. Thermal cycle conditions were 95˚C (10 min), 95˚C (30 s) and 58˚C (1 min) for 40 cycles, 98˚C (10 min). Droplet detection was carried out using QuantaSoft software (Bio-Rad) and the gene copy number was calculated for each sample in triplicate.

RNA extraction and transcriptional studies

RNA was extracted using SV Total RNA Isolation System (Promega) following the standard protocol with minor modifications. Cell pellets of 10 OD600 units were resuspended in 1 M sorbitol, 0.1 M EDTA pH 7.4 and 0.1% β-mercaptoethanol. For cell lysis 175 µL of RNA Lysis Buffer was added and the suspension transferred to a microcentrifuge tube filled with RNAase-free glass beads. Cells were then vortex mixed for 20 min, centrifuged at maximum speed and the supernatant transferred to a new microcentrifuge tube. The extraction was continued according to the standard protocol. Absolute quantification was performed using One-Step RT-ddPCR Advanced Kit for Probes (Bio-Rad) for reverse transcription and PCR amplification following the manufacture’s specifications. For amplification of AppA and ARG4 cDNA the same set of primers and probes described above were used. For pdi cDNA a set of primers (Fwd 5′-GCAGGAGTCGAGTCGCTAGTGT-3′ and Rev 5′-TTGCCTCGGCGATTGTGTCG-3′) and a 5′-hexachloro-fluorescein (HEX) dye-labelled probe were designed (5′-HEX-TTCTCTCTGAGAGTACCGGCAACTCCG-3′). The probe was also labelled with Iowa Black® FQ at the 3′ terminus. All primers and probes were synthetised by IDT, Australia. Each PCR was performed in a 20 µL reaction volume containing 5 µL of Supermix (Bio-Rad), 2 µL of reverse transcriptase (Bio-rad), 1 µL of 300 mM DTT, 900 nM of each primer, 250 nM of ARG4-HEX, AppA-FAM probe or 125 nM of PDI-HEX probe, 0.3 ng of RNA and the required amount of MilliQ water. Droplets were generated using the Droplet Generator (Bio-Rad) and transferred into a 96-well plate for retrotranscription and amplification. Thermal cycle conditions were 60˚C (60 min), 95˚C (10 min), 95˚C (30 s) and 58˚C (1 min) for 40 cycles, 98˚C (10 min). Droplet detection was carried out using QuantaSoft software (Bio-Rad) and absolute quantification determination was calculated for each sample in triplicate.

Enzyme activity quantification

Acid phosphatase activity was assayed to measure phytase expression levels using a assay based on the substrate para-nitrophenylphosphate (p-NPP) (Sigma) at an initial concentration of 5 mM in the assay [17]. Briefly, 10 µL of five enzyme dilutions were incubated with 90 µL of p-NPP substrate in 250 mM sodium acetate buffer pH 5.5 for 15 min at 37 °C. The reaction was stopped by the addition of 10 µL of 1 M NaOH. The released para-nitrophenol was measured at 410 nm after 10 min incubation at room temperature. Reactions were conducted in triplicate in 96-well plates.

One unit of acid phosphatase activity was defined as the amount of enzyme catalysing the formation of 1 µmole of para-nitrophenol per minute under the assay conditions.

Results

Transcriptional optimisation of AppA was conducted using methanol inducible, de-repressed, and constitutive MDPs from P. pastoris and H. polymorpha. Several BDPs were then tested for co-expression of the phytase together with folding and secretion chaperones, comparing phytase production yields with the commonly used PAOX1.

The co-expression of AppA with chaperones for disulfide bond formation, vesicle transport between organelles and ROS detoxification was conducted to address potential bottlenecks for AppA secretion in P. pastoris. Transcriptional analysis using droplet digital PCR was applied to follow the expression of AppA phytase and PDI transcripts during production, bringing new insights into the BDP expression profile.

Transcriptional optimization of AppA phytase expression

The AppA gene from E. coli was expressed under the control of the commonly used methanol inducible P. pastoris MDPs PAOX1, PDAS1 and PCAT1, the methanol inducible PHpFMD from H. polymorpha [62] and the constitutive P. pastoris promoter PGAP (Fig. 1). No difference in growth rate was observed between the strains (data not shown). The highest production of AppA phytase was observed for AppA (PHpFMD) strain showing a 1.53 ± 0.01-fold increase after 132 h compared to AppA (PAOX1) (Table 1). Expression from the constitutive promoter PGAP showed very low production of AppA phytase under the conditions tested (Fig. 1).

Fig. 1
figure 1

AppA expression from PHpFMD surpasses P. pastoris endogenous MDPs (PAOX1, PDAS1, PCAT1 and PGAP). Phytase expression was determined by the p-NPP assay. All strains except AppA (PGAP) were induced with methanol at 65 h and consecutive additions until 132 h. Data are represented as mean values ± standard deviation (n = 3)

Table 1 Fold change of AppA phytase production relative to MDP PAOX1PHpFMD or BDP PHpFMD-HpMOX 

AppA phytase co-expression with HAC1 and ER-associated disulfide bond formation chaperones using bidirectional promoters

AppA phytase was then co-expressed with the unfolded protein response transcriptional regulator HAC1 under the methanol inducible BDPs PAOX1-CAT1, PHpFMD-HpMOX and PDAS1/2 [63, 66, 67]. Previous reports have showed that co-expression with HAC1 increases expression of several recombinant proteins [3, 22, 25, 59, 59, 60, 60]. Figure 2a shows a comparison between co-expression of HAC1 with different BDPs. Expression of AppA phytase showed 1.36 ± 0.05-fold increase after 132 h when co-expressing HAC1 compared to the AppA (PFMD), and 2.08 ± 0.01-fold increase when compared to the expression from the MDP PAOX1 (Table 1). In the case of AppA-HAC1 (PAOX1-CAT1) and AppA-HAC1 (PDAS1/2) strains, the co-expression of HAC1 appeared to improve the production of the phytase by a similar relative ratio compared to the corresponding MDP expression strain AppA (PAOX1) and AppA (PDAS1). AppA-HAC1 (PAOX1-CAT1) showed a 1.31 ± 0.01-fold increase compared to AppA (PAOX1) (Table 1), while phytase production in AppA-HAC1 (PDAS1/2) improved 1.38 ± 0.07-fold compared to AppA (PDAS1). The results show that the effect of HAC1 on phytase production is independent of the promoter driving HAC1 gene expression, (PHpMOX, PCAT1 or PDAS2) leading to a ~ 30% increased yield in all cases, suggesting that the threshold of HAC1 regulator required to activate UPR is reached when expressed from any of the promoters.

Fig. 2
figure 2

HAC1 or PDI co-expression from BDPs further increases phytase yields. a AppA co-expressed with HAC1 under PAOX1-CAT1, PDAS1/2 and PHpFMD-HpMOX. b AppA co-expressed with PDI under the BDPs PAOX1-CAT1, PDAS1/2 and PHpFMD-HpMOX. c AppA co-expressed with MPDI and EUG1 under the BDP PHpFMD-HpMOX. The AppA-HAC1 (PHpFMD-HpMOX) and AppA-PDI (PHpFMD-HpMOX) strains were included for comparison. Phytase expression was determined by the pNPP assay after methanol induction. Data are represented as mean values ± standard deviation (n = 3)

One of the disulfide bonds is AppA phytase is a non-consecutive bond (C155-C430) which increases the opportunity for mispairing [39]. In E. coli the correct formation of this bond depends on DsbC, a protein-disulfide isomerase capable of correcting incorrectly formed bonds [4]. For enhancement of AppA phytase production in P. pastoris, the overexpression of a protein-disulfide isomerase might be advantageous. The co-expression of PDI using the PHpFMD-HpMOX BDP improved production of AppA phytase by 2.14 ± 0.03-fold compared to AppA (PAOX1) and 1.40 ± 0.02-fold compared to AppA (PHpFMD), showing very similar results to AppA-HAC1 (PHpFMD-HpMOX) strain (Table 1). Since the PDI gene is regulated by the UPR, the similarity in phytase production between AppA-PDI (PHpFMD-HpMOX) and AppA-HAC1 (PHpFMD-HpMOX) strains could indicate an activation of PDI gene transcription by HAC1 in the AppA-HAC1 (PHpFMD-HpMOX) strain, hinting towards a dominant effect of PDI in the HAC1 co-expression.

The co-expression of recombinant proteins and chaperones requires fine tuning, as both their relative ratios and cumulative amounts are important for the final expression levels [63, 66, 67]. The co-expression of AppA with PDI under BDPs with alternative expression profiles could result advantageous for final phytase yields by balancing protein amounts relative to each other. Based on this idea, strains AppA-PDI (PAOX1-CAT1) and AppA-PDI (PDAS1/2) were also constructed. As shown in Fig. 2b, the co-expression of AppA and PDI from the PAOX1-CAT1 or PDAS1/2 did not improve production compared to the expression from PHpFMD-HpMOX. In fact, both strains showed decreased phytase production. Nevertheless, the co-expression of PDI was still favourable for AppA phytase production when compared with MDP AppA (PAOX1) and AppA (PDAS1) strains. The fold change observed compared to AppA-PDI (PHpFMD-HpMOX) was 1.71 ± 0.07 and 3.93 ± 0.5 for AppA-PDI (PAOX1-CAT1) and AppA-PDI (PDAS1/2), respectively.

From the BDPs tested by co-expression of HAC1 or PDI, PHpFMD-HpMOX showed the highest phytase yields and corresponding strains were constructed using this BDP for co-expression of different chaperones (Fig. 2a, b). These results demonstrate for the first time that orthologous promoters can be fused into MDPs maintaining their beneficial properties related to strength and regulation, thereby extending the applicability of this concept beyond MDPs [62].

The co-expression of PDI has been shown previously to increase the expression of recombinant proteins in P. pastoris, however PDI homologues have not yet been tested for co-expression of heterologous proteins in this yeast. Considering that different homologues are proposed to interact and mediate thiol exchange with specific protein targets, we tested two PDI homologues for co-expression with AppA phytase (Fig. 2c) [57]. MPDI from P. pastoris presents high similarity to PDI while EUG1 from S. cerevisiae does not have a homologue in P. pastoris [54, 55]. Co-expression of MPDI had a marked detrimental effect on phytase production, while EUG1 improved production compared to AppA (PAOX1) but was less than AppA (PHpFMD) (Fig. 2c, Table 1), suggesting a possible overloading of the protein synthesis cellular machinery.

Characterisation of bidirectional promoter P HpFMD-HpMOX expression profile

Transcriptional analysis of AppA and PDI genes were conducted to study the expression profile from the BDP PHpFMD-HpMOX during phytase production and to verify the expression of PDI (Fig. 3). We used the AppA-PDI (PHpFMD-HpMOX) strain co-expressing AppA from PHpFMD, PDI from PHpMOX, because of the high phytase yields obtained. Since the total expression of PDI in the AppA-PDI (PHpFMD-HpMOX) strain is derived from the native PDI gene and the integrated PDI gene, the AppA (PHpFMD) strain was included for comparison where only expression from the native PDI gene was present. Transcriptional expression was analysed during glucose depletion (24, 48 h), immediately before methanol addition (65 and 72 h) and 1 h after methanol addition (66 and 73 h). The AppA phytase and PDI transcript levels relative to the ARG4 reference gene were consistent and reached similar levels after each methanol addition (Fig. 3). Transcription of the AppA gene reached similar levels in the AppA (PHpFMD) and AppA-PDI (PHpFMD-HpMOX) strains showing that the final production levels were associated with the co-expression of PDI. Both AppA and PDI transcription from PHpFMD-HpMOX started during glucose depletion sometime on or before 24 h. Transcription of both genes kept increasing up to 65 h, most likely due to the continuing drop in glucose levels and enhanced de-repression of the promoters. Transcriptional analysis of PDI in the AppA (PHpFMD) strain indicated the level of transcription associated with the native PDI gene alone. The transcription of PDI reached up to 3.8 times that of the ARG4 gene, consistent with the increased transcription of the AppA phytase during methanol addition (Fig. 3b). The cumulative transcription of the isomerase in AppA-PDI (PHpFMD-HpMOX), corresponding to the native and extra copy of PDI, showed an approximate tenfold increase after each methanol addition compared to the expression of the native gene in AppA (PHpFMD).

Fig. 3
figure 3

Transcriptional studies of AppA and PDI genes during phytase production. a Relative transcription of the AppA gene and b the PDI gene at 24, 48, 65, 66, 72 and 73 h. Methanol additions were performed at 65 and 72 h. ARG4 was used as reference gene. Data are represented as mean values ± standard deviation (n = 3)

Assuming the stability of AppA and PDI mRNAs to be similar, a marked difference in the strength of the PHpFMD and PHpMOX promoters was observed when comparing the transcription of AppA and PDI. Phytase expression from PHpFMD reached a relative ratio of ~ 400 (Fig. 3a) while PDI expression from PHpMOX showed a maximum relative ratio of ~ 40 after 1 h methanol additions (66 and 73 h) (Fig. 3b).

Co-expression of AppA phytase with trafficking associated proteins

Protein trafficking pathways can represent a bottleneck for heterologous protein production in P. pastoris [15, 31]. Two strains co-expressing AppA phytase and proteins involved in trafficking were constructed, AppA-SEC1 (PHpFMD-HpMOX) and AppA-SLY1 (PHpFMD-HpMOX). SLY1 regulates ER-Golgi trafficking and SEC1 interacts with vesicle trafficking between the Golgi and cell membrane [28, 41]. No marked improvement of phytase production was observed for either AppA-SEC1 (PHpFMD-HpMOX) or AppA-SLY1 (PHpFMD-HpMOX) strains when compared to the AppA (PHpFMD) strain (Additional file 2: Fig. S1, Table 1). The results suggest that either vesicular transport of the phytase is not a bottleneck, or that the co-overexpression of SEC1 or SLY1 proteins is not enough to alleviate the burden on the secretion pathway. The co-expression of other transport associated proteins or the combination with PDI might be required to improve AppA production.

AppA phytase production is improved when co-expressing with PDI and ERV2, ERO1 or GPX1

Due to the requirement for restoring oxidised PDI and the potential related limitation in the production of phytase in the AppA-PDI (PHpFMD-HpMOX) strain, four new strains co-expressing flavoproteins ERO1 and ERV2 were constructed [20, 50]. Two different promoters (PGAP and PCAT1) were tested for the expression of ERO1 and ERV2 due to potential problems arising from cumulative amounts of expressed proteins. Furthermore, timing of expression may also play an important role for improved yields hence the constitutive PGAP and methanol inducible PCAT1 were tested for expression of the second chaperone. The strains were named AppA-PDI-ERO1 (PGAP) and AppA-PDI-ERV2 (PGAP), AppA-PDI-ERO1 (PCAT1) and AppA-PDI-ERV2 (PCAT1) when ERO1 and ERV2 were expressed under the control of PGAP or PCAT1, respectively. Figure 4a shows and improvement in phytase production for AppA-PDI-ERO1 (PGAP) and AppA-PDI-ERV2 (PGAP), 1.21 ± 0.06 and 1.35 ± 0.04-fold increase from AppA-PDI (PHpFMD-HpMOX), respectively. On the other hand, the strains AppA-PDI-ERO1 (PCAT1) and AppA-PDI-ERV2 (PCAT1) did not show a marked improvement when compared to AppA-PDI (PHpFMD-HpMOX), with a 1.09 ± 0.11 and 1.13 ± 0.05-fold change, respectively. The expression of the flavoproteins ERO1 and ERV2 under the control of the constitutive promoter PGAP appeared to be more suitable compared to the expression from PCAT1. The effect could be related to the earlier expression of the flavoproteins from PGAP before production of PDI and AppA by glucose de-repression and methanol induction [29, 42, 71]. Expression of ERV2 or ERO1 under PGAP could alleviate stress on the cell machinery by production of the proteins during growth and not during methanol induction when PDI and AppA phytase were being produced (Fig. 5).

Fig. 4
figure 4

ERO1, ERV2 and GPX1 co-expression further improves phytase yields while co-expression of HAC1 or KAR2 does not show improvement. a AppA co-expression with PDI and ERO1 or ERV2. b AppA co-expression with ERO1 or ERV2. c AppA co-expression with GPX1. d AppA co-expression with HAC1 or KAR2. Phytase production was determined by the p-NPP assay after methanol induction. The AppA-PDI (PHpFMD-HpMOX) strain was included for comparison. Data are represented as mean values ± standard deviation (n = 3)

Fig. 5
figure 5

Overview of the strain engineering strategies applied in this work and secretory pathway bottlenecks explored. Column graphs represent AppA phytase activity (mU/mL) after 132 h of cultivation. Strains names are abbreviated for simplicity: AppA (PAOX1) as PAOX1, AppA (PDAS1) as PDAS1, AppA (PCAT1) as PCAT1, AppA (PGAP) as PGAP, AppA (PHpFMD) as PHpFMD, AppA-HAC1 (PHpFMD-HpMOX) as HAC1, AppA-PDI (PHpFMD-HpMOX) as PDI, AppA-EUG1 (PHpFMD-HpMOX) as EUG1, AppA-MPDI (PHpFMD-HpMOX) as MPDI, AppA-ERV2 (PHpFMD-HpMOX) as ERV2, AppA-ERO1 (PHpFMD-HpMOX) as ERO1, AppA-SEC1 (PHpFMD-HpMOX) as SEC1, AppA-SLY1 (PHpFMD-HpMOX) as SLY1, AppA-PDI-KAR2 (PGAP) as PDI-KAR2, AppA-PDI-GPX1 (PCAT1) as PDI-GPX1, AppA-PDI-ERV2 (PGAP) as PDI-ERV2 and AppA-PDI-ERO1 (PGAP) as PDI-ERO1

The improvement of phytase production in AppA-PDI-ERO1 (PGAP) and AppA-PDI-ERV2 (PGAP) strains could also be related to the co-expression of the flavoproteins, independently of PDI overexpression. AppA-ERO1 (PHpFMD-HpMOX) and AppA-ERV2 (PHpFMD-HpMOX) strains were constructed to test this hypothesis. As shown in Fig. 4b phytase production in the AppA-ERV2 (PHpFMD-HpMOX) strain was very similar (1.11 ± 0.01-fold change) to AppA-PDI (PHpFMD-HpMOX). On the other hand, the decrease in phytase production from the AppA-ERO1 (PHpFMD-HpMOX) strain compared to AppA-PDI (PHpFMD-HpMOX) (0.37 ± 0.05-fold change) could indicate a futile overload of the production machinery. Similar production of AppA phytase in AppA-ERV2 (PHpFMD-HpMOX) and AppA-PDI (PHpFMD-HpMOX) suggests that ERV2 flavoprotein might play a direct role in disulfide bond formation as previously proposed [24]. The results suggest that the improvement in AppA-PDI-ERO1 (PGAP) and AppA-PDI-ERV2 (PGAP) strains was most likely a synergistic effect of expressing PDI together with ERO1 or ERV2, and regeneration of oxidised PDI for further disulfide bonding in AppA phytase.

Overexpression of heterologous proteins with a high number of disulfide bonds can affect cellular redox homeostasis [14, 57]. To address a potential redox imbalance as a limitation in phytase production, two strains were constructed co-expressing the antioxidant enzyme GPX1. The strains were named AppA-PDI-GPX1 (PGAP) and AppA-PDI-GPX1 (PCAT1) when GPX1 was expressed under the control of the PGAP or PCAT1, respectively. Contrary to the expression of ERO1 and ERV2 flavoproteins, production of AppA phytase was improved when GPX1 was expressed under PCAT1 (1.30 ± 0.14-fold change), and decreased under PGAP expression (0.80 ± 0.01 fold change) compared to AppA-PDI (PHpFMD-HpMOX) strain (Fig. 4c). Since GPX1 is a peroxidase, expression during growth under PGAP might create an unnecessary redox imbalance that affects phytase production. Importantly, co-expression of GPX1 did not affect growth of either strain (data not shown). The expression of GPX1 under PCAT1 appears to work better with the expression of the AppA phytase under methanol induction when the maintenance of redox homeostasis is required and not necessary during the growth phase.

Based on the initial increase in AppA phytase production when co-expressing HAC1 (Fig. 2a), two more strains were constructed to test for further improvement, AppA-PDI-HAC1 (PCAT1) and AppA-HAC1-PDI (PGAP). In the AppA-PDI-HAC1 (PCAT1) strain HAC1 was expressed from the PCAT1 and in AppA-HAC1-PDI (PGAP) strain, HAC1 was expressed from the BDP PHpFMD-HpMOX and PDI is expressed from PGAP. The AppA-HAC1-PDI (PGAP) strain was designed to test if the effect on phytase production was similar when PDI was expressed under constitutive or inducible expression from PGAP or PHpMOX (in BDP PHpFMD-HpMOX), respectively. AppA-PDI-HAC1 (PCAT1) showed a similar yields to AppA-PDI (PHpFMD-HpMOX) (1.06 ± 0.16-fold change) suggesting that the co-expression of HAC1 under PCAT1 did not significantly contribute to AppA production in a PDI expressing strain background (Fig. 4d). These results may indicate that PDI and HAC1 lead to a similar, non-additive, beneficial effect. Furthermore, when HAC1 was expressed under the BDP PHpFMD-HpMOX and PDI under PGAP, no marked improvement in phytase production was observed compared to AppA (PFMD) (1.13 ± 0.15-fold change) and the expression slightly decreased compared to AppA-PDI (PHpFMD-HpMOX) by 0.80 ± 0.11-fold (Fig. 4d). The expression of PDI under the inducible BDP PHpFMD-HpMOX appeared to play an important role in the improvement of AppA production.

Over-expression of ER chaperone KAR2 has been shown to improve secretion of recombinant proteins [12, 23, 52, 72, 74], to test for improvement of AppA production, strain AppA-PDI-KAR2 (PGAP) was constructed. In this case, the expression of the KAR2 chaperone was only studied using PGAP. Results in Fig. 4d show that the co-expression of the chaperone did not further improve phytase production compared to the AppA-PDI (PHpFMD-HpMOX) strain (0.94 ± 0.02-fold change).

Discussion

The development of P. pastoris as a production platform has made available a diverse toolbox for genetic improvement of strains. Using some of these tools, we have explored bottlenecks in folding and secretion of E. coli AppA phytase to gain a deeper understanding of P. pastoris as a production host. AppA phytase is widely used in animal feed to improve the digestibility of phytic acid. Extensive research has been conducted towards the improvement of the enzymatic properties of AppA as well as optimisation of fermentation processes [8, 9, 18, 21, 42, 44, 47, 53, 73]. Yet, fewer studies on strain engineering for enhancement of production yields have been conducted using P. pastoris as production host.

Transcriptional optimisation of AppA expression was conducted using the methanol inducible PAOX1, PDAS1, PCAT1 and PHpFMD, and the constitutive PGAP, obtaining highest phytase yields under the control of PHpFMD. These results confirms the applicability of the orthologous FMD promoter from H. polymorpha in P. pastoris also for secreted proteins, that are difficult to fold, whereas previous reports have focused mostly on intracellular recombinant proteins [62]. An investigation of consecutive checkpoints along the folding and secretion pathway was then conducted using BDPs (PAOX1-CAT1, PDAS1/2 and PHpFMD-HpMOX). Notably, among the BPDs studied, PHpFMD-HpMOX showed the highest AppA production during co-expression with chaperones (Fig. 2a and b, Table 1). This shows that the two orthologous promoters from H. polymorpha (PHpFMD and PHpMOX) can be combined into BDPs whilst maintaining their specific strength and regulation [62].

Synergistic effects of fusing different pairs of MDPs to each other were previously investigated [63, 66, 67]. For a few BDPs involving shortened MDPs (˂500 bp), where putative 5′ insulating ends were removed, the authors showed elevated expression or ‘transcriptional spill over’. However, nearly all full fusions, where the length of each MDP was maintained (> 500 bp), retained similar expression levels as in the monodirectional state [63, 66, 67]. In this work, the PHpFMD and PHpMOX promoters maintained their typical length and were not drastically shortened for the generation of the BDP fusion [62]. Hence, transcriptional spill over would be unlikely to occur between PHpFMD and PHpMOX considering the long 5′ regions in between. The PHpFMD and the PHpMOX promoters in addition to being strongly inducible by methanol also allow heterologous gene expression upon de-repression in the presence of growth limiting concentrations of carbon sources like glucose and glycerol [62]. This feature can be advantageous for many industrial fermentation facilities that are not certified for the use of large volumes of flammable solvents such as methanol.

Starting with HAC1 and PDI, different parts of the secretion machinery were explored including the over-expression of isomerases and flavoproteins involved in disulfide bond formation (Fig. 2). Digital PCR was used to determine that the gene copy number for all strains was one, enabling appropriate comparison of the chaperone effect without the interference of the effect of different gene copy numbers. Digital PCR was also used to follow AppA and PDI transcripts during phytase production and to study the transcriptional profile of the BDP PHpFMD-HpMOX. An approximate tenfold increase in transcription was observed for the PHpFMD compared to PHpMOX evidenced by the transcript levels of AppA and PDI. Based on this result, and previous reports showing that the expression of the gene of interest is in most cases the limiting factor for production and not the folding chaperone [63, 66, 67], we co-expressed AppA phytase under the stronger promoter PHpFMD and the folding chaperone under the weaker promoter PHpMOX. The rationale was to maximise the yield of the secreted protein of interest, whereas the intracellularly produced, mostly ER- targeted folding helpers are constrained by the cellular machinery and excessive amounts of chaperones may impair the secretion apparatus. Using PHpMOX for PDI (or HAC1) co-expression showed very high AppA expression comparable to expression from PAOX1. We hypothesise that switching to PHpFMD would give extremely high and unnecessary chaperone levels whilst reducing the amount of AppA by using the weaker PHpMOX for its expression. Constitutive activation of the UPR through overexpression of HAC1 has been shown to enhance production of heterologous proteins in several systems including S. cerevisiae, P. pastoris, A. niger and mammalian cells [3, 22, 25, 59, 59, 60, 60]. Results from this work show that co-expression of HAC1 also improves production of AppA (Fig. 1b, Table 1). This regulator is most likely promoting secretion of the phytase by alleviating ER stress through the activation of UPR genes.

The co-expression of PDI improved production of AppA phytase in a similar way to HAC1. The isomerase might be important for the correct formation of the four disulfide bonds present in AppA structure, particularly for the non-consecutive bond. Yet, the expression of both HAC1 and PDI together did not have a synergistic effect in production of the phytase (Table 1). A possibility could be that the expression of HAC1 in AppA-HAC1 strain is inducing expression of PDI to similar levels as the expression of PDI in AppA-PDI (PHpFMD-HpMOX) strain. This hypothesis could also explain the inability to increase production of the phytase when co-expressing HAC1 in the AppA-PDI-HAC1 strain and the decrease of production in AppA-HAC1-PDI (Fig. 4d, Table 1). In fact, results obtained for the AppA-HAC1-PDI (PGAP) strain showed that expression of PDI from PGAP was detrimental for phytase production, supporting the idea of fine-tunning and cumulative amounts as key factors determining success of multiple proteins co-expression [63, 66, 67].

The co-expression of the flavoproteins ERO1 and ERV2 was also investigated. ERO1 has been co-expressed with recombinant proteins in yeast improving the expression of human antibody Fab fragment in S. cerevisiae, while having no effect on the expression of trypsinogen P. pastoris [14, 23]. Due to the intrinsic feedback regulation of ERO1 and shuffling of disulfide bonds through cysteine pairs within the ERO1 structure proposed by Sevier et al. the overexpression of ERO1 without over-expression of PDI might not directly contribute to disulfide bond formation [14]. Results obtained for the AppA-ERO1 strain support this idea. Phytase production was not improved by sole co-expression of ERO1, while the co-expression PDI and ERO1 (from PGAP) slightly improved phytase production by 1.21 ± 0.06-fold (Fig. 4a, Table 1). Co-expression of PDI and ERO1 have been shown to increase the yields of lipase in P. pastoris [51].

The ERV2 flavoprotein has been characterised as an alternative pathway for disulfide bond formation in S. cerevisiae. ERV2 transfers oxidising equivalents to PDI by a dithiol–disulphide exchange reaction in a similar way as that of ERO1-dependent pathway [10, 24, 50, 58]. Gerber et al. reported that ERV2 can also directly oxidise substrate proteins due to the flexibility of the C-terminal arm containing a C-G-C motif that is responsible for disulfide exchange [24]. The putative membrane associated ERV2 flavoprotein from P. pastoris, which has 44% identity to ERV2 from S. cerevisiae, has not yet been characterised. Results from this work suggest that ERV2 in P. pastoris may represent an alternative pathway to ERO1 as in S. cerevisiae and may perhaps directly catalyse the formation of disulfide bonds as proposed by Gerber et al. [24]. The AppA-ERV2 (PHpFMD-HpMOX) strain presents similar phytase production to the AppA-PDI (PHpFMD-HpMOX) strain, and co-expression of PDI and ERV2 (under PGAP) further enhances phytase production by 1.35 ± 0.04-fold (Fig. 4a, Table 1). The capacity to re-establish reduced PDI to its oxidised state, either with ERO1 or ERV2, appears to be advantageous when over-expressing AppA to promote re-shuffling of possible incorrectly formed disulfide bonds. Further, the improvement in phytase secretion was dependent on the promotor used for the co-expression of ERO1 or ERV2. Only the strains where the flavoprotein was expressed from PGAP showed improvement in phytase production. This observation could be related to timing the expression of the chaperones during growth and the requirement for flavoprotein availability ahead of phytase and PDI over-expression during glucose de-repression and methanol induction. On the other hand, the effect could be related to the strength of the promoter, where the expression driven by PCAT1 could potentially saturate the expression machinery and result in decreased AppA production when the cells are over-expressing two proteins (AppA phytase and PDI).

The formation of de novo disulfide bonds and isomerisation of incorrect ones can create a redox imbalance and generate reactive oxygen species (ROS) in the ER [14, 57]. The glutathione redox system controls redox homeostasis by regulating the glutathione redox ratio (GSH/GSSG). As part of this system, the cytosolic glutathione peroxidase GPX1 acts as an antioxidant enzyme that detoxifies ROS at the expense of reduced glutathione [14]. The over-expression of this enzyme has shown to improve the production of recombinant proteins in P. pastoris by increasing the oxidising environment of the ER [14]. Results presented here show that over-expression of GPX1 in AppA-PDI-GPX1 (PCAT1) strain slightly improve phytase production (Fig. 4c, Table 1). In this case, inducible expression of GPX1 under PCAT1 shows better results than constitutive expression under PGAP. The over-expression of GPX1 during growth could potentially create a redox imbalance in conditions where there is no folding oxidative stress.

Lastly, we also explored the co-expression of the vesicular transport proteins SLY1 and SEC1. Once heterologous proteins are correctly folded they can often be retained intracellularly resulting in poor secretion [31, 33]. However, the co-expression of SLY1 or SEC1 did not improve production of AppA phytase (Additional file 2: Fig. S1, Table 1). The limited AppA production in the AppA-SEC1 and AppA-SLY1 strains might still be related to incorrect folding and disulfide bond formation in the ER. Assessing phytase production in strains co-expressing PDI and SEC1 or SLY1 could be used to investigate whether a posterior trafficking bottleneck affects secretion of the phytase and limits productivity, however strains were not constructed in this work to test this hypothesis.

Conclusions

The work presented here explores a wide variety of strain engineering possibilities for optimising protein expression in P. pastoris using E. coli AppA phytase as an example of a challenging protein of industrial importance. Multiple strategies to improve phytase production were tested since a single strategy does not exist to improve secretion of recombinant proteins.

This work used the rational design of strains to improve secretion of AppA, building on previous investigations and bringing new insights into secretion bottlenecks. Transcriptional optimisation revealed that PHpFMD was the best promoter for phytase expression, improving yields by ~ 1.5- fold compared to PAOX1. Co-expression with HAC1 or PDI further improved yields by ~ 1.4-fold compared to the expression without chaperone. Finally, co-expression of a second chaperone, ERV2, further increased yields by ~ 2.0-fold, achieving a total synergistic improvement of ~ 3.0-fold compared to the initial yields obtained with PAOX1. Future work can now be performed to understand if the same production profiles would be observed for each strain under industrial fermentation conditions. The use of newly available tools, like BDPs, shows that protein production improvements can readily be achieved compared to the commonly used PAOX1 promoter and that there are a range of possibilities for further development of this important and widely used yeast microbial factory.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. Abad S, Kitz K, Hörmann A, Schreiner U, Hartner FS, Glieder A. Real-time PCR-based determination of gene copy numbers in Pichia pastoris. Biotechnol J Healthc Nutr Technol. 2010;5(4):413–20.

    CAS  Google Scholar 

  2. Barlowe CK, Miller EA. Secretory protein biogenesis and traffic in the early secretory pathway. Genetics. 2013;193(2):383–410.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Becker E, Florin L, Pfizenmaier K, Kaufmann H. An XBP-1 dependent bottle-neck in production of IgG subtype antibodies in chemically defined serum-free Chinese hamster ovary (CHO) fed-batch processes. J Biotechnol. 2008;135(2):217–23.

    Article  CAS  PubMed  Google Scholar 

  4. Berkmen M, Boyd D, Beckwith J. The nonconsecutive disulfide bond of Escherichia coli phytase (AppA) renders it dependent on the protein-disulfide isomerase. DsbC J Biol Chem. 2005;280(12):11387–94.

    Article  CAS  PubMed  Google Scholar 

  5. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72(1–2):248–54.

    Article  CAS  PubMed  Google Scholar 

  6. Cámara E, Albiol J, Ferrer P. Droplet digital PCR-aided screening and characterization of Pichia pastoris multiple gene copy strains. Biotechnol Bioeng. 2016;113(7):1542–51.

    Article  PubMed  CAS  Google Scholar 

  7. Chahal S, Wei P, Moua P, Park SPJ, Kwon J, Patel A, Vu AT, Catolico JA, Tsai YFT, Shaheen N. Structural characterization of the α-mating factor prepro-peptide for secretion of recombinant proteins in Pichia pastoris. Gene. 2017;598:50–62.

    Article  CAS  PubMed  Google Scholar 

  8. Chen C-C, Wu P-H, Huang C-T, Cheng K-J. A Pichia pastoris fermentation strategy for enhancing the heterologous expression of an Escherichia coli phytase. Enzyme Microbial Technol. 2004;35(4):315–20.

    Article  CAS  Google Scholar 

  9. Chen C-C, Wu P-H, Huang C-T, Cheng K-J. A Pichia pastoris fermentation strategy for enhancing the heterologous expression of an Escherichia coli phytase. Enzyme Microb Technol. 2004;35(4):315–20.

    Article  CAS  Google Scholar 

  10. Collet J-F, Bardwell JC. Disulfides out of thin air. Nat Struct Biol. 2002;9(1):2.

    Article  CAS  PubMed  Google Scholar 

  11. Cregg JM, Madden K, Barringer K, Thill G, Stillman C. Functional characterization of the two alcohol oxidase genes from the yeast Pichia pastoris. Mol Cell Biol. 1989;9(3):1316–23.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Damasceno LM, Anderson KA, Ritter G, Cregg JM, Old LJ, Batt CA. Cooverexpression of chaperones for enhanced secretion of a single-chain antibody fragment in Pichia pastoris. Appl Microbiol Biotechnol. 2007;74(2):381–9.

    Article  CAS  PubMed  Google Scholar 

  13. 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(3):414–37.

    Article  CAS  PubMed  Google Scholar 

  14. Delic M, Rebnegger C, Wanka F, Puxbaum V, Haberhauer-Troyer C, Hann S, Köllensperger G, Mattanovich D, Gasser B. Oxidative protein folding and unfolded protein response elicit differing redox regulation in endoplasmic reticulum and cytosol of yeast. Free Radical Biol Med. 2012;52(9):2000–12.

    Article  CAS  Google Scholar 

  15. Delic M, Valli M, Graf AB, Pfeffer M, Mattanovich D, Gasser B. The secretory pathway: exploring yeast diversity. FEMS Microbiol Rev. 2013;37(6):872–914.

    Article  CAS  PubMed  Google Scholar 

  16. Dersjant-Li Y, Awati A, Schulze H, Partridge G. Phytase in non-ruminant animal nutrition: a critical review on phytase activities in the gastrointestinal tract and influencing factors. J Sci Food Agric. 2015;95(5):878–96.

    Article  CAS  PubMed  Google Scholar 

  17. Engelen AJ, Van Der Heeft FC, Randsdorp PH, Smtt EL. Simple and rapid determination of phytase activity. J AOAC Int. 1994;77(3):760–4.

    Article  CAS  PubMed  Google Scholar 

  18. Fei B, Xu H, Cao Y, Ma S, Guo H, Song T, Qiao D, Cao Y. A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase. J Ind Microbiol Biotechnol. 2013;40(5):457–64.

    Article  CAS  PubMed  Google Scholar 

  19. Frand AR, Cuozzo JW, Kaiser CA. Pathways for protein disulphide bond formation. Trends Cell Biol. 2000;10(5):203–10.

    Article  CAS  PubMed  Google Scholar 

  20. Frand AR, Kaiser CA. Ero1p oxidizes protein disulfide isomerase in a pathway for disulfide bond formation in the endoplasmic reticulum. Mol Cell. 1999;4(4):469–77.

    Article  CAS  PubMed  Google Scholar 

  21. Garrett JB, Kretz KA, O’Donoghue E, Kerovuo J, Kim W, Barton NR, Hazlewood GP, Short JM, Robertson DE, Gray KA. Enhancing the thermal tolerance and gastric performance of a microbial phytase for use as a phosphate-mobilizing monogastric-feed supplement. Appl Environ Microbiol. 2004;70(5):3041–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Gasser B, Maurer M, Gach J, Kunert R, Mattanovich D. Engineering of Pichia pastoris for improved production of antibody fragments. Biotechnol Bioeng. 2006;94(2):353–61.

    Article  CAS  PubMed  Google Scholar 

  23. Gasser B, Sauer M, Maurer M, Stadlmayr G, Mattanovich D. Transcriptomics-based identification of novel factors enhancing heterologous protein secretion in yeasts. Appl Environ Microbiol. 2007;73(20):6499–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Gerber J, Mühlenhoff U, Hofhaus G, Lill R, Lisowsky T. Yeast ERV2p is the first microsomal FAD-linked sulfhydryl oxidase of the Erv1p/Alrp protein family. J Biol Chem. 2001;276(26):23486–91.

    Article  CAS  PubMed  Google Scholar 

  25. Guerfal M, Ryckaert S, Jacobs PP, Ameloot P, Van Craenenbroeck K, Derycke R, Callewaert N. The HAC1 gene from Pichia pastoris: characterization and effect of its overexpression on the production of secreted, surface displayed and membrane proteins. Microb Cell Fact. 2010;9(1):49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Harland B (1999) Phytic acid complex in feed ingredient. Phytase Anim Nutri Waste Manag

  27. Hartner FS, Glieder A. Regulation of methanol utilisation pathway genes in yeasts. Microb Cell Fact. 2006;5(1):39.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Hashizume K, Cheng Y-S, Hutton JL, C-h C, Carr CM. Yeast Sec1p functions before and after vesicle docking. Mol Biol Cell. 2009;20(22):4673–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Hollenberg C, Janowicz Z (1987) DNA-molecules coding for FMDH control regions and structural gene for a protein having FMDH-activity and their use thereof. European Patent Application. 299108.

  30. Holst B, Tachibana C, Winther JR. Active site mutations in yeast protein disulfide isomerase cause dithiothreitol sensitivity and a reduced rate of protein folding in the endoplasmic reticulum. J Cell Biol. 1997;138(6):1229–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hou J, Tyo K, Liu Z, Petranovic D, Nielsen J. Engineering of vesicle trafficking improves heterologous protein secretion in Saccharomyces cerevisiae. Metab Eng. 2012;14(2):120–7.

    Article  CAS  PubMed  Google Scholar 

  32. Hui-Ying L, Huang H-q, Ying-Guo B, Ya-Ru W, Pei-Long Y, Kun M, Tie-Zheng Y, Bin Y. Improving phytase expression by increasing the gene copy number of appA-m in Pichia pastoris. Chin J Biotechnol. 2006;22(4):528–33.

    Article  Google Scholar 

  33. Idiris A, Tohda H, Sasaki M, Okada K, Kumagai H, Giga-Hama Y, Takegawa K. Enhanced protein secretion from multiprotease-deficient fission yeast by modification of its vacuolar protein sorting pathway. Appl Microbiol Biotechnol. 2010;85(3):667–77.

    Article  CAS  PubMed  Google Scholar 

  34. Kies AK, De Jonge LH, Kemme PA, Jongbloed AW. Interaction between protein, phytate, and microbial phytase In vitro studies. J Agric Food Chem. 2006;54(5):1753–8.

    Article  CAS  PubMed  Google Scholar 

  35. Laboissiere MC, Sturley SL, Raines RT. The essential function of protein-disulfide isomerase is to unscramble non-native disulfide bonds. J Biol Chem. 1995;270(47):28006–9.

    Article  CAS  PubMed  Google Scholar 

  36. LaMantia M, Lennarz WJ. The essential function of yeast protein disulfide isomerase does not reside in its isomerase activity. Cell. 1993;74(5):899–908.

    Article  CAS  PubMed  Google Scholar 

  37. Lei XG, Weaver JD, Mullaney E, Ullah AH, Azain MJ. Phytase, a new life for an “old” enzyme. Ann Rev Anim Biosci. 2013;1(1):283–309.

    Article  CAS  Google Scholar 

  38. Li C, Lin Y, Zheng X, Pang N, Liao X, Liu X, Huang Y, Liang S. Combined strategies for improving expression of Citrobacter amalonaticus phytase in Pichia pastoris. BMC Biotechnol. 2015;15(1):88.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Lim D, Golovan S, Forsberg CW, Jia Z. Crystal structures of Escherichia coli phytase and its complex with phytate. Nat Struct Biol. 2000;7(2):108–13.

    Article  CAS  PubMed  Google Scholar 

  40. 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(1):44–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Malsam J, Kreye S, Söllner T. Membrane fusion: SNAREs and regulation. Cell Mol Life Sci. 2008;65(18):2814–32.

    Article  CAS  PubMed  Google Scholar 

  42. Mayer A, Hellmuth K, Schlieker H, Lopez-Ulibarri R, Oertel S, Dahlems U, Strasser A, Van Loon A. An expression system matures: A highly efficient and cost-effective process for phytase production by recombinant strains of Hansenula polymorpha. Biotechnol Bioeng. 1999;63(3):373–81.

    Article  CAS  PubMed  Google Scholar 

  43. Näätsaari L, Mistlberger B, Ruth C, Hajek T, Hartner FS, Glieder A. Deletion of the Pichia pastoris KU70 homologue facilitates platform strain generation for gene expression and synthetic biology. PLoS ONE. 2012;7(6):e39720.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Ndayambaje JB, Meenakshisundaram S. Co-feeding strategy to enhance phytase production in Pichia pastoris. Afr J Biotechnol. 2014;13:21.

    Article  CAS  Google Scholar 

  45. Nørgaard P, Westphal V, Tachibana C, Alsøe L, Holst B, Winther JR. Functional differences in yeast protein disulfide isomerases. J Cell Biol. 2001;152(3):553–62.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Peña DA, Gasser B, Zanghellini J, Steiger MG, Mattanovich D. Metabolic engineering of Pichia pastoris. Metab Eng. 2018;50:2–15.

    Article  PubMed  CAS  Google Scholar 

  47. Rodriguez E, Wood ZA, Karplus PA, Lei XG. Site-directed mutagenesis improves catalytic efficiency and thermostability of Escherichia coli pH 2.5 acid phosphatase/phytase expressed in Pichia pastoris. Arch Biochem Biophys. 2000;382(1):105–12.

    Article  CAS  PubMed  Google Scholar 

  48. Ruohonen L, Toikkanen J, Tieaho V, Outola M, SODERLUND H and KERANEN S, . Enhancement of protein secretion in Saccharomyces cerevisiae by overproduction of Sso protein, a late-acting component of the secretory machinery. Yeast. 1997;13(4):337–51.

    Article  CAS  PubMed  Google Scholar 

  49. Scherens B, Dubois E, Messenguy F. Determination of the sequence of the yeast YCL313 gene localized on chromosome III. Homology with the protein disulfide isomerase (PDI gene product) of other organisms. Yeast. 1991;7(2):185–93.

    Article  CAS  PubMed  Google Scholar 

  50. Sevier CS, Cuozzo JW, Vala A, Åslund F, Kaiser CA. A flavoprotein oxidase defines a new endoplasmic reticulum pathway for biosynthetic disulphide bond formation. Nat Cell Biol. 2001;3(10):874.

    Article  CAS  PubMed  Google Scholar 

  51. Sha C, Yu X-W, Zhang M, Xu Y. Efficient secretion of lipase r27RCL in Pichia pastoris by enhancing the disulfide bond formation pathway in the endoplasmic reticulum. J Ind Microbiol Biotechnol. 2013;40(11):1241–9.

    Article  CAS  PubMed  Google Scholar 

  52. Shusta EV, Raines RT, Plückthun A, Wittrup KD. Increasing the secretory capacity of Saccharomyces cerevisiae for production of single-chain antibody fragments. Nat Biotechnol. 1998;16(8):773–7.

    Article  CAS  PubMed  Google Scholar 

  53. Stahl CH, Wilson DB, Lei XG. Comparison of extracellular Escherichia coli AppA phytases expressed in Streptomyces lividans and Pichia pastoris. Biotechnol Lett. 2003;25(10):827–31.

    Article  CAS  PubMed  Google Scholar 

  54. Tachibana C, Stevens T. The yeast EUG1 gene encodes an endoplasmic reticulum protein that is functionally related to protein disulfide isomerase. Mol Cell Biol. 1992;12(10):4601–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Tachikawa H, Takeuchi Y, Funahashi W, Miura T, Gao X-D, Fujimoto D, Mizunaga T, Onodera K. Isolation and characterization of a yeast gene, MPD1, the overexpression of which suppresses inviability caused by protein disulfide isomerase depletion. FEBS Lett. 1995;369(2–3):212–6.

    Article  CAS  PubMed  Google Scholar 

  56. Tai H-M, Yin L-J, Chen W-C, Jiang S-T. Overexpression of Escherichia coli phytase in Pichia pastoris and its biochemical properties. J Agric Food Chem. 2013;61(25):6007–15.

    Article  CAS  PubMed  Google Scholar 

  57. Tyo KE, Liu Z, Petranovic D, Nielsen J. Imbalance of heterologous protein folding and disulfide bond formation rates yields runaway oxidative stress. BMC Biol. 2012;10(1):16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Vala A, Sevier CS, Kaiser CA. Structural determinants of substrate access to the disulfide oxidase Erv2p. J Mol Biol. 2005;354(4):952–66.

    Article  CAS  PubMed  Google Scholar 

  59. Valkonen M, Penttilä M, Saloheimo M. Effects of inactivation and constitutive expression of the unfolded-protein response pathway on protein production in the yeast Saccharomyces cerevisiae. Appl Environ Microbiol. 2003;69(4):2065–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Valkonen M, Ward M, Wang H, Penttilä M, Saloheimo M. Improvement of foreign-protein production in Aspergillus niger var. awamori by constitutive induction of the unfolded-protein response. Appl Environ Microbiol. 2003;69(12):6979–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Vanz AL, Nimtz M, Rinas U. Decrease of UPR-and ERAD-related proteins in Pichia pastoris during methanol-induced secretory insulin precursor production in controlled fed-batch cultures. Microb Cell Fact. 2014;13(1):23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. 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(1):1–9.

    Article  CAS  Google Scholar 

  63. Vogl T, Gebbie L, Palfreyman RW, Speight R. Effect of plasmid design and type of integration event on recombinant protein expression in Pichia pastoris. Appl Environ Microbiol. 2018;84(6):e02712-02717.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Vogl T, Glieder A. Regulation of Pichia pastoris promoters and its consequences for protein production. Nat Biotechnol. 2013;30(4):385–404.

    CAS  Google Scholar 

  65. Vogl T, Hartner FS, Glieder A. New opportunities by synthetic biology for biopharmaceutical production in Pichia pastoris. Curr Opin Biotechnol. 2013;24(6):1094–101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Vogl T, Kickenweiz T, Pitzer J, Sturmberger L, Weninger A, Biggs BW, Köhler E-M, Baumschlager A, Fischer JE, Hyden P. Engineered bidirectional promoters enable rapid multi-gene co-expression optimization. Nat Comm. 2018;9(1):3589.

    Article  CAS  Google Scholar 

  67. Vogl T, Kickenweiz T, Pitzer J, Sturmberger L, Weninger A, Biggs BW, Köhler E-M, Baumschlager A, Fischer JE, Hyden P. Engineered bidirectional promoters enable rapid multi-gene co-expression optimization. Nat Comm. 2018;9(1):1–13.

    CAS  Google Scholar 

  68. Waterham HR, Digan ME, Koutz PJ, Lair SV, Cregg JM. Isolation of the Pichia pastoris glyceraldehyde-3-phosphate dehydrogenase gene and regulation and use of its promoter. Gene. 1997;186(1):37–44.

    Article  CAS  PubMed  Google Scholar 

  69. Weis R, Luiten R, Skranc W, Schwab H, Wubbolts M, Glieder A. Reliable high-throughput screening with Pichia pastoris by limiting yeast cell death phenomena. FEMS Yeast Res. 2004;5(2):179–89.

    Article  CAS  PubMed  Google Scholar 

  70. Weninger A, Hatzl A-M, Schmid C, Vogl T, Glieder A. Combinatorial optimization of CRISPR/Cas9 expression enables precision genome engineering in the methylotrophic yeast Pichia pastoris. J Biotechnol. 2016;235:139–49.

    Article  CAS  PubMed  Google Scholar 

  71. Weydemann U, Keup P, Piontek M, Strasser A, Schweden J, Gellissen G, Janowicz Z. High-level secretion of hirudin by Hansenula polymorpha—authentic processing of three different preprohirudins. Appl Microbiol Biotechnol. 1995;44(3–4):377–85.

    Article  CAS  PubMed  Google Scholar 

  72. Xu P, Raden D, Doyle FJ III, Robinson AS. Analysis of unfolded protein response during single-chain antibody expression in Saccaromyces cerevisiae reveals different roles for BiP and PDI in folding. Metab Eng. 2005;7(4):269–79.

    Article  CAS  PubMed  Google Scholar 

  73. Yao M-Z, Wang X, Wang W, Fu Y-J, Liang A-H. Improving the thermostability of Escherichia coli phytase, appA, by enhancement of glycosylation. Biotechnol Lett. 2013;35(10):1669–76.

    Article  CAS  PubMed  Google Scholar 

  74. Zhang W, Hl Z, Xue C, Xiong Xh, Xq Y, Xy Li, Hp C, Zm L. Enhanced secretion of heterologous proteins in Pichia pastoris following overexpression of Saccharomyces cerevisiae chaperone proteins. Biotechnol Prog. 2006;22(4):1090–5.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to acknowledge the Central Analytical Research Facility, operated by the Institute for Future Environments (QUT). Access to CARF is supported by funding from the Science and Engineering Faculty (QUT). T.V.’s contribution was supported by an Endeavour Research Fellowship by the Australian Government Department of Education and Training. T.V. gratefully acknowledges supported from the Austrian Science Fund (FWF, Erwin Schrödinger fellowship J 4256).

Funding

Funding for this work was provided by an Advance Queensland Innovation Partnership program grant awarded to the Queensland University of Technology and The University of Queensland, along with co-funding from collaborating industry partners Ridley AgriProducts Pty Ltd, Bioproton Pty Ltd and Kennedy Creek Lime Pty Ltd.

Author information

Authors and Affiliations

Authors

Contributions

LN, TV and RS designed the experiments, interpreted the results, and wrote the manuscript. LN and TV conducted the experimental work. RP, JB, XC, CL-F and JH participated in the results interpretation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Laura Navone.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Additional tables.

Additional file 2: Fig. S1.

SEC1 or SLY1 co-expression does not further improve phytase yields. Phytase expression was determined by the p-NPP assay after methanol induction. The AppA (PHpFMD) strain was included for comparison. Data are represented as mean values ± standard deviation (n = 3).

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Navone, L., Vogl, T., Luangthongkam, P. et al. Synergistic optimisation of expression, folding, and secretion improves E. coli AppA phytase production in Pichia pastoris. Microb Cell Fact 20, 8 (2021). https://doi.org/10.1186/s12934-020-01499-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12934-020-01499-7

Keywords