- Research article
- Open Access
Systems-level organization of yeast methylotrophic lifestyle
BMC Biologyvolume 13, Article number: 80 (2015)
Some yeasts have evolved a methylotrophic lifestyle enabling them to utilize the single carbon compound methanol as a carbon and energy source. Among them, Pichia pastoris (syn. Komagataella sp.) is frequently used for the production of heterologous proteins and also serves as a model organism for organelle research. Our current knowledge of methylotrophic lifestyle mainly derives from sophisticated biochemical studies which identified many key methanol utilization enzymes such as alcohol oxidase and dihydroxyacetone synthase and their localization to the peroxisomes. C1 assimilation is supposed to involve the pentose phosphate pathway, but details of these reactions are not known to date.
In this work we analyzed the regulation patterns of 5,354 genes, 575 proteins, 141 metabolites, and fluxes through 39 reactions of P. pastoris comparing growth on glucose and on a methanol/glycerol mixed medium, respectively. Contrary to previous assumptions, we found that the entire methanol assimilation pathway is localized to peroxisomes rather than employing part of the cytosolic pentose phosphate pathway for xylulose-5-phosphate regeneration. For this purpose, P. pastoris (and presumably also other methylotrophic yeasts) have evolved a duplicated methanol inducible enzyme set targeted to peroxisomes. This compartmentalized cyclic C1 assimilation process termed xylose-monophosphate cycle resembles the principle of the Calvin cycle and uses sedoheptulose-1,7-bisphosphate as intermediate. The strong induction of alcohol oxidase, dihydroxyacetone synthase, formaldehyde and formate dehydrogenase, and catalase leads to high demand of their cofactors riboflavin, thiamine, nicotinamide, and heme, respectively, which is reflected in strong up-regulation of the respective synthesis pathways on methanol. Methanol-grown cells have a higher protein but lower free amino acid content, which can be attributed to the high drain towards methanol metabolic enzymes and their cofactors. In context with up-regulation of many amino acid biosynthesis genes or proteins, this visualizes an increased flux towards amino acid and protein synthesis which is reflected also in increased levels of transcripts and/or proteins related to ribosome biogenesis and translation.
Taken together, our work illustrates how concerted interpretation of multiple levels of systems biology data can contribute to elucidation of yet unknown cellular pathways and revolutionize our understanding of cellular biology.
Methylotrophic yeasts accept a broad range of carbon sources. Multicarbon sources, such as sugars and sugar alcohols like glucose, glycerol, or mannitol, are utilized at similar efficiency as reduced C1-compounds like methanol . Besides the proper equipment of the cells with enzymes necessary for substrate metabolism, their coordinated expression is a prerequisite for the successful utilization of different carbon and energy sources. The methylotrophic yeast Pichia pastoris (syn. Komagataella sp.) is widely used for recombinant protein production with several biopharmaceuticals on the market  and an expanding portfolio of industrial enzymes produced . Recently, the application of P. pastoris as a model system for peroxisome and secretory organelle proliferation has also expanded [4, 5]. The methylotrophic lifestyle has been the main driving force for this development, as it involves strong and regulated promoters used for expression of recombinant genes , as well as specialized organelles, the peroxisomes. Peroxisomes are defined as intracellular compartments accommodating hydrogen peroxide (H2O2) forming oxidases together with the H2O2 detoxifying enzyme catalase. Also the fatty acid beta-oxidation pathway of P. pastoris is located in these organelles . Yeast peroxisomal oxidases are predominantly involved in the metabolism of various unusual carbon and nitrogen sources (e.g. alcohols, fatty acids, D-amino acids, or primary amines) . In methylotrophic yeasts, peroxisomes, which harbor the initial steps of the methanol utilization pathway, are highly abundant in methanol-grown cells but become heavily decreased in both number and volume upon catabolite repression . When grown on glucose, Hansenula polymorpha, another methylotrophic yeast, harbors only a single, small peroxisome which can serve as a source for proliferation by fission when induction is triggered by shifting the cells to methanol [10, 11]. In addition to genes encoding structural peroxisomal proteins, the expression of methanol utilization related genes is strongly induced on methanol. The first steps of methanol assimilation involve an alcohol oxidase (AOX) to convert methanol to formaldehyde, and a special transketolase named dihydroxyacetone synthase (DAS) to form a C-C bond with the C1 molecule formaldehyde. The reactions of these two enzymes and their localization to peroxisomes are well described [12, 13]. The further reaction cycle of methanol assimilation is supposed to involve pentose phosphate reactions, but the details are not fully clarified to date.
While there are several studies analyzing cellular reactions of P. pastoris to methanol induction in context of recombinant protein production [14–18], the response of non-recombinant strains to the different carbon sources is largely unknown. Thus, we decided to investigate the cellular responses of P. pastoris cells not producing a recombinant protein to methanol and glucose, respectively, which are the two most widely used substrates for cultivation. To enable the same chemostat-controlled constant specific growth rates for direct comparability the methanol cultures were co-fed with glycerol. Availability of whole genome sequences made a number of transcriptome regulation studies of P. pastoris, analyzing the implications of growth rate , unfolded protein response (UPR) induction , oxygen availability , osmotic stress , or heterologous protein production [16, 23], become feasible. Analyses of the host proteome gave further insights into characteristics of P. pastoris grown at different temperatures , osmolarity , UPR induction , and oxygen supply . More recently, P. pastoris strains producing an insulin precursor were analyzed for changes in the cellular proteome as adaptation response to methanol induction during fed batch cultivation using 2D-DIGE and subsequent mass spectrometry identification of differentially abundant proteins. High abundance of enzymes from the dissimilatory methanol metabolism and induction of the UPR were observed . Regulation of cellular enzyme concentrations will cause changes in metabolic fluxes, eventually also leading to changes in free metabolite concentrations. Quantitative determination of intracellular fluxes is the key to a better understanding of metabolic networks. First genome-scale metabolic network models of P. pastoris [26, 27] and flux distributions of central carbon metabolism [28–30] indicate growth rate-related methanol (co-)assimilation with proposed implications for the pentose phosphate pathway .
The work at hand incorporates transcriptomics, proteomics, metabolomics, and fluxomics analyses of non-producing P. pastoris in steady-state cultures at a uniform specific growth rate comprising the carbon source as the investigated variable. This integrated systems level analysis allowed to reveal cellular processes that are co-regulated with methanol metabolism, such as vitamin biosynthesis and amino acid metabolism. Furthermore, these co-regulation patterns were the pre-requisite to elucidate the thus far unidentified steps of sugar phosphate rearrangements recycling xylulose-5-phosphate for methanol fixation. We propose, herein, a new model for the assimilation of methanol as a separate strictly regulated pathway, originating from duplication of the involved genes.
Results and discussion
Growth parameters of P. pastoris differ significantly on different substrates
P. pastoris CBS7435 was cultivated in chemostat cultivations at a fixed specific growth rate of 0.1 h−1, corresponding to approximately 60 % of μmax on glucose . Constant growth is a prerequisite to avoid growth rate-dependent effects during genome-scale analyses. As the maximum specific growth rate on pure methanol as a carbon source would be significantly lower, and intracellular carbon fluxes could not be analyzed on methanol alone, a mixed feed strategy applying glycerol-methanol co-feeding was employed. A methanol-glycerol mix of 8.5 g/L methanol and 49.0 g/L glycerol was employed based on experiments with P. pastoris proving that total methanol utilization and full induction of the methylotrophic pathway were realized under these conditions. Chemostats were run in three biological replicates per condition and samples for transcriptomics, proteomics, and metabolite analyses were taken in steady state after seven residence times as described in the Methods section. For metabolic flux analysis, separate chemostat cultivations employing 13C-labelled substrates were performed. Substrate limitation of all cultures, i.e. no residual glucose or methanol/glycerol, respectively, was confirmed by HPLC. The growth parameters derived from these cultures are summarized in Table 1. The CO2 exchange rate of cells grown on methanol/glycerol was 13 % lower compared to those grown on glucose while their oxygen uptake rate was 30 % higher. The higher oxygen uptake rate of methanol/glycerol-grown P. pastoris can be explained by the higher degree of reduction of methanol and glycerol compared to glucose. As methanol oxidation to formaldehyde by AOX is an exothermic oxygen consuming reaction, an equimolar amount of oxygen is needed only to pass methanol into cellular metabolism. The biomass yield was slightly higher for cells grown on methanol/glycerol compared to glucose, which is in good agreement with data from the literature [28, 30]. Transcriptional regulation was analyzed using P. pastoris-specific DNA microarrays [20, 32], liquid chromatography-tandem mass spectrometry (LC/MS-MS) was used for differential proteomics and quantification of metabolites. Additionally, distribution of specific lipid classes were analyzed. Flux ratios were calculated from 13C labelling patterns in proteinogenic amino acids. The numerical results of these genome scale analyses can be found in Additional file 1.
Transcriptome and proteome are significantly co-regulated
At the transcriptional level, 406 of 5,354 genes were significantly differentially expressed on methanol/glycerol and glucose. As protein abundance, however, does not necessarily directly correlate with transcription , we also measured differential proteome regulation using 2D-LC-MS of Tandem Mass Tag labelled total protein samples and obtained quantitative data for 575 cellular proteins. In agreement with the literature [33–35], where a positive correlation between protein concentration and the abundance of the transcript has been described, we could mainly quantify proteins with higher transcript levels. Figure 1 shows 575 data pairs with mean log2 fold changes in transcript and protein levels of the methanol/glycerol experiments compared to the glucose experiments. Proteins and their transcripts were significantly co-regulated (r = 0.78, r2 = 0.61). Lu et al.  have shown that, in Saccharomyces cerevisiae, protein levels are determined to 73 % by transcription. Similarly, we observed that transcriptional control determined the regulation of protein abundance by 61 %.
Of the 575 genes with available transcriptomics and proteomics data, 130 (23 %) were differentially regulated at the transcript and/or protein level, the largest group being up-regulated at both levels. Based on the differential changes at both protein and transcript levels, data have been allocated to seven groups and analyzed for overrepresentation of functional groups (Fig. 1 and Additional file 1). As expected, during growth on methanol/glycerol there are strongly increased levels of transcripts and proteins involved in methanol metabolism and peroxisome formation, while two proteins needed mainly on glucose, hexokinase and high affinity glucose transporter, had lower abundance on both transcript and protein level. Higher levels of proteins of the translation machinery and cytoskeleton organization in cells grown on methanol/glycerol were not met by higher transcript levels, thus indicating a post-transcriptional regulation while the significant down-regulation of transcripts for lower glycolysis and fatty acid beta-oxidation was not reflected in protein levels (Fig. 1). These processes will be described in more detail below.
To confirm that the gene regulations attributed to methanol cultivation in this work are truly due to methanol, and not to glycerol as a co-substrate, we compared these regulation patterns to transcript regulation obtained in fed batch cultivations using methanol, glycerol, or glucose (Additional file 2). Thereby, we could confirm that all genes discussed to be regulated by methanol utilization in the present study are actually induced by methanol.
Peroxisome proliferation is strongly up-regulated on methanol
Methanol-induced cells show an up-regulation in genes encoding proteins essential for peroxisome biogenesis and proliferation. Fourteen PEX-genes were markedly up-regulated on the transcriptome level, three of them also on the proteome level (Additional file 1). Genes of the peroxisomal import machinery encoding the docking complex (PEX13, PEX14, PEX17) and RING-finger complex (PEX2, PEX10, PEX12) are up-regulated at almost equal levels. Among peroxins required for the import of peroxisomal matrix proteins , PEX11 and its isoform PEX11C were amongst the highest up-regulated genes (3- to 6-fold higher levels), while members of the Pex23-family were not induced. This distinction seems to be specific for P. pastoris, as both gene groups were up-regulated upon methanol induction in H. polymorpha . A similar regulation of P. pastoris PEX genes has also been described recently by Prielhofer et al. , who observed that genes encoding the Pex7/Pex20-mediated import machinery and the Pex23-family were only up-regulated in conditions inducing expression of beta-oxidation genes but not upon methanol induction. Most key players of the peroxisomal methanol utilization pathway, such as AOX, catalase, and DAS, rely on the Pex5-mediated PTS1 import pathway ; thus, up-regulation of receptors, which recognize a peroxisomal targeting signal sequence (PTS), was restricted to PTS1-specific PEX5 . Conversely, PEX7, encoding a signal receptor for the signal sequence PTS2, was slightly down regulated. This finding is in agreement with previous studies [39, 40] demonstrating that H. polymorpha and P. pastoris do not require the PTS2 import pathway for growth on methanol. Unchanged expression levels of Pex7 and Pex20 encoding its accessory protein also correlated with unaltered protein levels of beta-oxidation enzymes. Genes encoding auxiliary functions in matrix protein import and quality control were induced: the putative peroxisomal Lon-protease Pim1-2 and the peroxisomal ATP importer Pmp47 (PAS_chr3_0099) required for import of DAS were strongly up-regulated (log2FC +3.78). Further, up-regulation at both transcriptome and protein level was detected for the glutathione peroxidase Pmp20, a peroxisomal protein which might also be involved in the detoxification of H2O2 in the peroxisome of methanol growing cells.
The xylulose-monophosphate cycle of methanol assimilation utilizes a duplicated methanol inducible enzyme set and is entirely localized to peroxisomes
The key players of the methanol utilization pathway have been identified during the last 30 years [12, 13]; however, major steps of the assimilation pathway still remain to be resolved. Briefly, methanol is oxidized to formaldehyde by AOX (Aox1 and Aox2 in P. pastoris) within the peroxisomes, thereby generating stoichiometric amounts of H2O2. Formaldehyde is further converted in two possible routes, either dissimilatory by glutathione-dependent formaldehyde dehydrogenase, S-formyl glutathione hydrolase, and formate dehydrogenase yielding NADH and CO2, or assimilatory by the action of DAS (Das1 and Das2 in P. pastoris). DAS catalyzes the fusion of formaldehyde to xylulose-5-phosphate (XYL5P), thereby generating dihydroxyacetone and glyceraldehyde-3-phosphate (GAP). These intermediates are further converted by dihydroxyacetone kinase (DAK), fructose-1,6-bisphosphate aldolase, and fructose-1,6-bisphosphatase, to finally yield one molecule of GAP per three molecules of methanol, which is then used for the generation of biomass and energy. It is generally assumed that XYL5P gets recycled through rearrangements in the pentose phosphate pathway (PPP), although the detailed mechanism of these rearrangements as well as the interplay of PPP and peroxisomes is still unknown.
All known enzymes of the methanol utilization pathway have significantly higher transcript and protein levels when methanol is present (Fig. 1, upper right quadrant). Interestingly, our analysis revealed that P. pastoris does not only have a second isoform of fructose-1,6-bisphosphate aldolase (designated as Fba1-2) as reported by Küberl et al. , but also isoforms of the PPP enzymes transaldolase (Tal1-2), ribose-5-phosphate ketol-isomerase (Rki1-2), and ribulose-5-phosphate 3-epimerase (Rpe1-2). All these isoforms were found among the group of up-regulated gene-protein pairs, except for Rpe1-2 which was not identified at the proteomic level. Sequence analysis predicted that Fba1-2, Tal1-2, Rki1-2, and Rpe1-2 each contain a PTS1 peroxisomal targeting signal [42, 43], indicating their potential involvement in a separate peroxisomal methanol assimilation pathway. On the contrary, their respective cytosolic or mitochondrial isoforms (Fba1-1, Tal1-1, Rki1-1, and Rpe1-1) were not differentially regulated and do not contain a peroxisomal targeting sequence. The same regulation pattern was observed comparing cultures grown on methanol alone to those grown on glycerol or glucose (Additional file 2).
Subsequently, cellular fractions enriched of highly pure peroxisomes were isolated from methanol- or glucose-grown P. pastoris according to the protocol established by Wriessnegger et al. , and subjected to proteomics analyses. In this way, we demonstrated that all relevant enzymes for methanol assimilation were present only in methanol-derived peroxisomal fractions, but not in glucose-derived peroxisome fractions (Table 2). The relative enrichment of proteins in the peroxisomal fractions compared to total cell homogenates was quantified as average weighted ratios of the peak areas of respective peptides. This was consequently only possible with methanol-derived samples where peptides of the proteins of interest had been identified. Table 3 shows MASCOT scores as indicators of identification, and the average ratios of protein abundance in peroxisomal vs. homogenate samples, normalized to Aox1. The methanol assimilation pathway enzymes discussed above were enriched at the same level or higher than Aox1 in the peroxisomal fractions, just as several selected peroxisomal proteins, while cytosolic proteins of glycolysis, PPP, and methanol dissimilation were either not identified at all or markedly depleted in peroxisomal preparations. Only DAK of the methanol assimilation pathway was rather depleted compared to Aox1. Luers et al.  described that DAK localizes to the cytosol despite having a PTS1 signal. We could, however, quantify DAK also in peroxisomal fractions of methanol-grown cells, indicating that this enzyme can localize in more than one compartment.
Additionally, one of the unidentified ORFs present among the up-regulated gene/protein pairs was identified to be the homolog of S. cerevisiae, YKR043C, which was recently reported to encode sedoheptulose-1,7-bisphosphatase (Shb17) . Further, this protein was found to be enriched in the peroxisomes in methanol-grown P. pastoris (Table 2). Shb17 was shown to hydrolyze sedoheptulose-1,7-bisphosphate (S1,7BP) to sedoheptulose-7-phosphate in a thermodynamically driven pathway for the synthesis of pentose-5-phosphates alternative to PPP . S1,7BP was not among the quantified metabolites in our initial metabolomics analyses due to the lack of a commercially available standard. After receiving purified S1,7BP from Amy Caudy (University of Toronto, CA), a previously unidentified substance with differential abundance could be unambiguously assigned as S1,7BP. The signal-to-noise ratios of all samples with glucose-treatment were below 5. In the methanol-grown samples, the signal-to-noise ratios were 84, 111, and 116, clearly indicating the presence of S1,7BP in methanol-grown P. pastoris cells, contrary to glucose-grown cells (Fig. 2). This result prompted us to reconsider the pentose phosphate rearrangements leading to the formation of XYL5P for methanol assimilation.
Methanol assimilation employs an alternative xylulose-5-phosphate forming pathway via sedoheptulose-1,7-bisphosphate
While it would be stoichiometrically possible that XYL5P is regenerated through the canonical non-oxidative branch of the pentose phosphate pathway, our genomic, transcriptomic, and proteomic data point to another direction of C1 assimilation. It appears most likely that P. pastoris and other methylotrophic yeasts evolved a specialized set of enzymes for sugar phosphate rearrangements which is specifically induced by growth on methanol and localizes to peroxisomes. Figure 3 shows our proposed pathway for the rearrangement reactions, with 1 GAP molecule per 3 molecules of methanol as the net result. Thereby, F6P (generated from GAP and dihydroxyacetone phosphate (DHAP) by the action of Fba1 and fructose-1,6-bisphosphatase) and another GAP are interconverted to erythrose-4-phopsphate and XYL5P in a transketolase reaction. Erythrose-4-phopsphate is then condensed with DHAP to form S1,7BP, a reaction shown to be catalyzed by the aldolase Fba1 in yeast and plants . We propose that peroxisomal Fba1-2 or Tal1-2 might be the responsible enzyme for this reaction in P. pastoris. Shb17 catalyzes the dephosphorylation of S1,7BP to sedoheptulose-7-phosphate, which is finally converted to two XYL5P by transketolase, Rki1-2, and Rpe1-2. As P. pastoris Tkl1 is cytosolic and not induced in the presence of methanol, we propose that Das1 and/or Das2, both homologs of Tkl1, catalyze this reaction. Overall, in a process driven by the net loss of one high-energy phosphate bond, Sbh17, together with transaldolase and transketolase, convert five moles of triose-phosphate into the three moles of XYL5P required for fixation of three moles of formaldehyde by DAS. Localization of this entire pathway in the same compartment makes import of XYL5P into peroxisomes obsolete, which was proposed to be necessary by Douma et al.  according to the classical model of methanol assimilation. Thus, the net peroxisomal flux of carbon would be one mole DHAP or GAP out of peroxisomes per three moles of methanol.
Based on our data, we propose here a novel carbon assimilation pathway (Fig. 3, left) that shares the concept of compartmentalization with plants  and cyanobacteria . According to this model, DAS is responsible for C-C bond formation similar to the mechanism of RuBisCO, followed by a cyclic pathway (the equivalent to the Calvin cycle; Fig. 3, right) for regeneration of the pentose phosphate substrate of the carboxylation reaction. Shb17 has been shown to drive the flux from erythrose-4-phosphate and DHAP toward ribose-5-phosphate during riboneogenesis in S. cerevisiae in a reaction similar to the Calvin cycle . A similar mechanism driving the flux towards XYL5P is proposed here in methanol-induced P. pastoris and probably also in other methylotrophic yeasts. Supporting this hypothesis we also found PTS1 containing isoforms of Fba1 and Tal1 in the H. polymorpha genome sequence by BLAST analysis.
Tandem gene duplication occurs with high frequency and has been reported to be a major contributor of new genetic material . Several models for the occurrence of gene duplications have been proposed (reviewed in [50, 51]). Among them unequal crossing over can lead to tandem duplication, as it is observed here. Duplicated genes have a high probability of being lost again unless they acquire a new function . Byun-McKay and Geeta  have proposed that subcellular relocalization of duplicate gene products may play an important role in stabilizing duplications and acquiring new functions. While they extend their idea only to N-terminal mutations modifying targeting sequences to the endoplasmic reticulum, mitochondria, or chloroplasts, it may well be that C-terminal mutations may have enabled peroxisomal relocation of duplicate gene products in an ancestor of methylotrophic yeasts. One may envisage that compartmentalized xylulose-5-monophosphate pathway enzymes would constitute novel functions which underwent positive selection, leading to a highly regulated peroxisomal pathway as observed herein, while leaving the PPP unaffected.
The central carbon metabolism is reverted to gluconeogenesis
Growth on non-carbohydrate carbon sources necessitates the synthesis of hexoses and pentoses for the biosynthesis of macromolecules, which is accomplished by reverting the carbon flux to gluconeogenesis. Glycolysis and gluconeogenesis share several enzymes, while the irreversible, highly exergonic steps of glycolysis are bypassed. Therefore, exactly these reactions are the control steps of flux direction, and their regulation indicates the activity of gluconeogenesis. Both methanol and glycerol enter the central carbon metabolism at the level of C3-molecules (DHAP and GAP). One of the key regulatory enzymes of the upper part of glycolysis/gluconeogenesis is fructose-1,6-bisphosphatase, which we found to be up-regulated in methanol/glycerol-grown cells at the transcriptomic and proteomic level (Fig. 4). The other key regulatory enzymes of the lower part, pyruvate carboxylase and phosphoenolpyruvate carboxykinase, showed no differential regulation comparing both conditions, which is consistent with the fact that carbon flux from methanol and glycerol enters the central carbon metabolism at the point of glyceraldehyde-3-phosphate. Correspondingly, there was only a minor difference in the calculated fluxes from glyceraldehyde-3-phosphate towards pyruvate while the upper glycolytic flux was reverted on methanol/glycerol towards glucose-6-phosphate (Fig. 5). Additionally, the absence of extracellular glucose rendered low and high affinity glucose transporters HXT1 and GTH1 as well as hexokinase HXK1 obsolete which all had significantly lower transcript (and protein) levels on methanol/glycerol (Fig. 1, lower left).
Flux through the lower branch of glycolysis was lower on methanol/glycerol (Fig. 5) which fits to the observed lower TCA-cycle flux and the accumulation of three glycolytic intermediates (2-phosphoglycerate, 3-phosphoglycerate, and phosphoenolpyruvate). Transcript and protein levels of most glycolytic enzymes did not change (Fig. 4). Similarly, in S. cerevisiae, a poor correlation between fluxes and transcript levels of genes of this pathway was observed . The low TCA-cycle flux indicates that methanol dissimilation is a major source for NADH and energy production in cells grown on methanol/glycerol.
Pentose phosphate pathway flux is increased on methanol/glycerol
The PPP serves for the generation of NADPH for reductive assimilatory processes and for the generation of ribose-5-phosphate as a precursor for nucleic acids. Other PPP intermediates are used as precursors for other metabolic pathways like synthesis of histidine, nucleotides, and riboflavin. Methanol assimilation, using sugar phosphate intermediates in a cyclic fashion, has to be regarded separately from PPP, as outlined above.
The high gluconeogenic flux on methanol/glycerol was accompanied by a high flux through the PPP (Fig. 5), enabling a high specific production rate of reduced NADPH of appr. 3 mmol g−1 h−1. About 10 % of this higher NADPH production is needed for amino acid synthesis for the higher protein content of cells grown on methanol/glycerol. A higher PPP flux also provides for more ribose for nucleotide and riboflavin synthesis (see below). Cytosolic pentose phosphate pathway genes, however, were not differentially regulated at the different media (Table 2), contrary to the methanol-induced peroxisomal isoforms suggested to be employed in methanol assimilation in this compartment.
Up-regulation of the malate-aspartate shuttle serves for mitochondrial import of NADH generated by methanol dissimilation and correlates with decreased TCA-cycle flux in methanol/glycerol-grown P. pastoris
The TCA-cycle as a central hub for cellular metabolism is dedicated to energy production and supply of precursors for several other metabolic pathways. Again, the most striking C-source-dependent differences were observed on the flux level, being 3.2 times lower on methanol/glycerol, mainly being controlled by the lower influx of acetyl-CoA into the TCA-cycle (Fig. 5). The low TCA-cycle flux on methanol is mainly diverted to glutamate, thus contributing only marginally to energy production. No significant changes in transcript or protein levels of genes connected to the TCA-cycle were observed. Nevertheless, we found marked differences in TCA-cycle metabolites. Citrate levels in methanol/glycerol-grown cells were higher than in glucose-grown cells, whereas levels for isocitrate, fumarate, and malate were lower (Fig. 4). Taken together, these data indicate that, on methanol/glycerol, the TCA-cycle reactions are mainly employed for production of metabolic precursors for biomass formation rather than producing energy through the respiratory chain. Methanol utilization has a major impact on the energy state of the cells, as two moles of NADH are produced via dissimilation of one mole of methanol to CO2. Intracellular flux calculation showed that about half of the methanol was dissimilated to CO2 and therefore additional NADH was produced which may consequently lead to down-regulation of TCA-cycle flux.
Dissimilatory oxidation of formaldehyde takes place in the cytosol. Therefore, the produced NADH has to be transported via the inner mitochondrial membrane to drive the generation of ATP. For the transport of electrons via the mitochondrial membrane several shuttle systems exist, most importantly the malate-aspartate shuttle. The homologs of the malate-α-ketoglutarate transporter Odc1  and the glutamate-aspartate transporter Agc1  were both highly up-regulated at transcript and protein levels, indicating the relevance of this NADH shuttle for methylotrophic ATP generation in mitochondria.
The glyoxylate cycle is active in methanol/glycerol-grown P. pastoris and mainly localizes to the peroxisomes
The glyoxylate cycle is necessary for the utilization of non-fermentable carbon sources because of its ability to convert acetyl-CoA into C4 compounds that can be used for ATP generation in the mitochondria . Isocitrate lyase converts isocitrate to glyoxylate and succinate, the former intermediate is then condensed with acetyl-CoA to form malate by malate synthase. Additionally, malate dehydrogenase, citrate synthase, and aconitase are required. This process is assumed to take place in the peroxisomes in non-Saccharomyces yeasts . Indeed, we found most of the enzymes to be present in the peroxisomal fraction in methanol-grown cells (Table 2), only isocitrate lyase was found solely on methanol but predominantly in the cytosolic fraction. Furthermore, methanol-grown cells had increased transcript and/or protein levels of both putative peroxisomal malate synthase (PAS_chr1-4_0459, which we named Mls2) as well as of the cytosolic malate synthase and malate dehydrogenase. In agreement with this data, we also found increased glyoxylate cycle fluxes in cells grown on methanol/glycerol. We propose that the generated C4 compounds are mainly used as precursors for the biosynthesis of TCA-cycle-derived amino acids when using methanol/glycerol as substrate, rather than being shuttled to gluconeogenesis. In this line, no up-regulation of phosphoenolpyruvate carboxykinase, a key enzyme in lower gluconeogenesis, could be seen.
Methanol-grown P. pastoris cells have a higher protein but lower free amino acid content
Protein is the largest macromolecular component of cells, creating the highest demand for energy, reduction equivalents, and carbon flux. A change in substrate forces the cell to adapt, e.g. by varying the total protein content. Methanol/glycerol-grown cells had a 35 % higher protein content, with 0.54 g(protein) g(cell dry weight (CDW))−1 compared to glucose-grown cells with 0.40 g(protein) g(CDW)−1 (Table 4). Consequently, the levels of protein-bound amino acids were generally higher in cells grown on methanol/glycerol, while the level of free intracellular amino acids was approximately 20 % lower in methanol/glycerol-grown cells. The higher specific protein synthesis rate on methanol/glycerol (0.054 g(protein) g(CDW)−1 h−1 vs 0.040 g(protein) g(CDW)−1 h−1 on glucose) creates a higher drain of amino acids towards protein synthesis, which explains the generally low levels of intracellular free amino acids. The higher demand for amino acids creates a metabolic pull for the respective synthesis pathways. The higher demand for amino acids was well supported by the transcriptional and/or post-transcriptional up-regulation of genes involved in the biosynthesis of all twenty amino acids (Fig. 6). For histidine which derives from one intermediate of the pentose phosphate pathway, we saw regulation of HIS1, while for the amino acids which derive from glycolytic intermediates, regulation of Ser2 (serine), AGX1 and GLY1 (glycine), TRP5-2 (tryptophan), Aro7 and Aro8 (tyrosine and phenylalanine), Ilv2 and Ilv5 (leucine, valine, isoleucine), and Alt1 (alanine) was observed. If we consider amino acids which derive from intermediates of the TCA-cycle, we saw regulation of LYS20, LYS21 and Lys2 (lysine), Gdh2 and Gdh3 (glutamate), AAT1 (aspartate), ASP1 (asparagine), Hom2 (precursor of threonine, methionine, and cysteine), THR1 (threonine), Met17 (precursor of methionine and cysteine), and Cys3 (cysteine).
While methanol utilization enzymes were up-regulated at the transcriptional level and thus increased at the protein level, the second major class of more abundant proteins (ribosomal proteins) were not transcriptionally regulated. This finding indicates an efficient post-transcriptional regulation mechanism, as described for S. cerevisiae , and a higher steady state demand of translational capacity on methanol. This last observation is supported also by a higher total protein content of biomass grown on methanol/glycerol and by the up-regulation of amino acid synthesis pathways. It remains to be elucidated whether this higher translational capacity on methanol is related to the observed higher recombinant protein production capacity of methanol-based expression strains. The general upregulation of protein synthesis in methanol-induced cultures did not coincide with a higher abundance of enzymes in the protein folding machinery. An accumulation of misfolded proteins as a result of heterologous gene expression has been observed many times in recombinant P. pastoris, leading to UPR activation (reviewed by Puxbaum et al. ) but appears to be absent in non-recombinant P. pastoris cultivated on methanol.
Protein folding, secretion, and degradation pathways are not affected by methanol as substrate
While we and others have observed a transient up-regulation of the UPR immediately following methanol induction [59, 60], no such regulation pattern was noticed in the methanol-adapted cells in steady state in the present study, thus ruling out the possibility that a permanently induced UPR positively influences recombinant protein production in methanol-grown cells. Contrary to Liang et al. , who detected up-regulation of endoplasmic reticulum protein processing, N-glycan biosynthesis, and protein export pathways when comparing recombinant protein secreting P. pastoris in chemostats with methanol as substrate, we did not see any changes in protein folding, secretory pathway, N-glycosylation, or proteasome both at the proteome and transcriptome level in the non-expressing strains in our study (Additional file 1). Protein synthesis, however, was obviously up-regulated on methanol, as described above.
High levels of methanol utilization enzymes require overproduction of vitamins and cofactors
Alcohol oxidase requires high riboflavin synthesis
AOX, which catalyzes the first reaction of methanol utilization, is a homooctamer with flavin adenine dinucleotide (FAD) as non-covalently bound prosthetic group. When methylotrophic yeast cells grow on methanol, AOX can account for up to 30 % of total cellular protein , and the FAD content of AOX alone amounts to 1.7 mg/g biomass. AOX predominantly oligomerizes in the peroxisomal matrix of methylotrophic yeasts [62, 63]. Experiments with H. polymorpha and P. pastoris revealed that insertion of FAD is an essential step prior to the assembly of AOX .
Almost the entire pathway leading to FAD is transcriptionally up-regulated when methanol is present (Fig. 7 and Additional files 1 and 2): RIB1 and RIB3 encode the first steps of the riboflavin biosynthesis pathway (with GTP and ribulose-5-phosphate as precursors, respectively), while RIB4 and RIB5 code for the last enzymes in the pathway. The induction of the riboflavin pathway during growth on methanol has been previously observed , but was not linked to AOX biosynthesis. Via the up-regulated FMN1, riboflavin is converted to flavin mononucleotide (FMN), a strong oxidizing cofactor of mitochondrial NADH-dehydrogenases (which, however, are not regulated). The generation of FAD from FMN is catalyzed by Fad1, which is strongly transcriptionally up-regulated (~8-fold). The bulk of FAD apparently goes into AOX as other cellular flavoproteins  are rather unaffected during growth on methanol with the FAD-requiring Gut2, Hem14 (both up-regulated), Pox1, Fmo1-1, and Fre2 (all three down-regulated) as exceptions. We observed no changes in free riboflavin, indicating that flux to riboflavin is up-regulated upon its high demand while synthesis is tightly regulated by its free intracellular concentration as described by Marx et al. .
Thiamine synthesis is strongly up-regulated due to demand of peroxisomal transketolases
There is also significant up-regulation of genes involved in thiamine (vitamin B1) and thiamine pyrophosphate (TPP) biosynthesis (THI20, THI6, THI80, THI4, THI13, THI73, THI21; Fig. 7 and Additional files 1 and 2) in cells grown on methanol. We have shown before that severe thiamine limitation is required for the induction of THI13 , indicating that induction of the methanol utilization pathway leads to intracellular thiamine deficiency.
TPP, the active derivative of thiamine, is the co-factor of decarboxylases, transketolases, and phosphoketolases. The homodimeric enzymes bind one Mg2+ ion and one TPP per subunit. In P. pastoris cultures grown on methanol, upregulation of TPP biosynthesis coincided with the high abundance of the TPP-containing enzymes Das1 and Das2, catalyzing the fixation of formaldehyde to XYL5P in the peroxisome. In this study, Das1 and Das2 had very high changes in transcript levels (~68-fold and 46-fold up-regulation, respectively) and one of the highest changes in protein level (~6.5-fold and 4.7-fold, respectively). On the contrary, the level of cytosolic transketolase (Tkl1) was unaffected. Thus, we conclude that the strong induction of DAS1/2 led to a limitation of thiamine availability which was compensated by induction of the thiamine synthesis pathway.
High nicotinamide levels are required for formaldehyde detoxification
In the methanol dissimilation pathway, formaldehyde is oxidized to carbon dioxide by two consecutive reactions catalyzed by formaldehyde dehydrogenase and formate dehydrogenase. On methanol, both enzymes are strongly increased both on transcript (5-fold and 19-fold, respectively) and protein (formaldehyde dehydrogenase not quantified, formate dehydrogenase 5-fold) levels. The two enzymes, which are mainly located to the cytosol, are required for detoxification of formaldehyde and formate and both use nicotinamide adenine dinucleotide as cofactor. The generated NADH provides energy for growth on methanol. In this respect, the total amount of nicotinamide in cells grown in the presence of methanol was nearly 10-fold higher than in the glucose-grown cells, and total NAD content is approx. 50 % higher. Expression of NMA1, encoding nicotinic acid mononucleotide adenylyltransferase, which is involved in the de novo biosynthesis of NAD as well as in the NAD salvage pathway , was up-regulated 2.4-fold.
Heme synthesis is up-regulated upon catalase demand for peroxide detoxification
Toxic H2O2 and formaldehyde are generated in the first step of methanol metabolism. The peroxisomal enzyme catalase, which is involved in the detoxification of H2O2, is transcriptionally up-regulated when methanol is present. Properly folded catalase incorporates a heme cofactor with an iron ion in the center, and needs to tetramerize to become active . CTA1 expression is up-regulated on methanol and, consequently, expression of almost all heme biosynthesis genes was up-regulated, including the rate-limiting steps HEM2 and HEM3 (log2FC +0.43 and +0.54). Pet18, a heme oxygenase-like protein, was also up-regulated at both the transcript and protein levels. Heme oxygenases catalyze the degradation of heme and produce iron. Down-regulation (0.6-fold) of a low-affinity Fe(II) transporter (FET4-2) and up-regulation (1.65-fold) of FTH1-1, a putative high affinity iron transporter involved in intravacuolar iron storage, points towards low iron levels in the presence of methanol.
On methanol, the general lipid metabolism is altered to allow peroxisome formation at the expense of lipid droplets
Environmental conditions and nutritional modifications often have dramatic effects on the composition of cellular membranes, which becomes apparent in lipid composition and regulation of lipid metabolism. Major key enzymes of lipid-related pathways were apparently not affected when culture conditions varied between the supply of glucose or methanol/glycerol. Important components of biological membranes, sterols and phospholipids, were elevated only slightly in methanol/glycerol-grown cells (Table 5). The observed increase of building blocks for membranes can most likely be explained by the enhanced occurrence of internal membranes. The total amount of peroxisomes was strongly increased in methanol/glycerol-grown cells which caused a weak effect on the total amounts of phospholipids and sterols of internal membranes. Wriessnegger et al. [44, 68] already showed in previous work that utilization of glucose or methanol as the sole carbon source does not lead to major differences in the distribution of phospholipids, although the culture conditions and sampling points were not the same as in the present study. The slight increase in the total amount of phospholipids observed here was not matched by any significant regulation of lipid biosynthetic genes involved in the complex pathways of phospholipid formation, except for INO1 and OPI3, which were both down-regulated (log2FC of −1.38 and −0.89). The pattern of fatty acids from methanol-grown cells as well mostly resembled glucose-grown cells, although some minor changes were detected. A decrease in oleic acid (C18:1) by roughly 20 % was accompanied by an increase in palmitic acid (C16:0), palmitoleic acid (C16:1), and linolenic acid (C18:3). Again, the influence of the intracellularly predominant peroxisomal membranes most likely was the reason for the observed changes of bulk membrane fatty acid composition.
The strongest effect on lipid classes resulting from cultivation on different carbon sources was on triacylglycerols (TAG), the major non-polar lipid of P. pastoris. Both TAG synthases, DGA1 and LRO1, were transcriptionally down-regulated on methanol/glycerol (log2FC −1.07 and −0.43). As a direct result, TAG were reduced in methanol/glycerol-grown cells by more than 50 % (Table 5). The significant decrease of TAG was accompanied by a severe reduction of lipid droplets in P. pastoris cultivated on methanol, which was observed by electron microscopy (Fig. 8). While the amount of TAG was severely reduced, precursors of TAG (diacylglycerols and free fatty acids) were increased by approximately 40 %. Upon mobilization of TAG by TAG lipases, activated fatty acids could serve as substrates either for β-oxidation or as building blocks for membrane formation. In comparison to glucose we observed on methanol/glycerol a down-regulation of transcripts encoding β-oxidation relevant genes as well as TAG forming enzymes, which was not followed at the protein level. Notably, it has been previously shown that genes involved in fatty acid utilization are differentially regulated upon using glycerol or glucose as the carbon source, and depend on substrate availability  (Additional file 2). The utilization of methanol enables P. pastoris cells for proper growth based on energy supply by alcohol oxidation, but apparently does not provide excess carbons to be incorporated in storage material. Therefore, non-polar lipid synthesizing enzymes are down-regulated. As a direct consequence, no alternative supply of fatty acids may be available and β-oxidation relevant enzymes are shut down as well because of the limited substrate available.
ERG20, encoding farnesyl pyrophosphate synthetase, is the only lipid biosynthetic gene which was found to be up-regulated when comparing glucose to methanol-grown cells. Erg20 is part of the sterol biosynthetic pathway, which is composed of more than 20 enzymes. However, all other sterol biosynthetic genes remained transcriptionally unaffected. Erg20 is located at an important branching point of this biosynthetic pathway. The product of the Erg20 catalyzed reaction, farnesyl pyrophosphate, cannot only serve as a substrate for the formation of structural lipid compounds which is one of the major routes, but can be directed towards several other pathways, among them heme biosynthesis. As the formation of heme was found to be transcriptionally up-regulated to serve as a prosthetic group of catalase, we anticipated that ERG20 was up-regulated predominantly to provide sufficient substrate for the de-novo formation of heme.
Methylotrophy is a unique ability of microorganisms to live on C1 molecules that requires efficient pathways to form C-C bonds and to oxidize C1 compounds via toxic intermediates. This systems level investigation provides comprehensive insight into regulatory and metabolic specificities of the methylotrophic yeast P. pastoris. Co-regulation of enzymes with AOX and DAS at the transcript and protein level allowed us to identify in detail the putative pathway for XYL5P regeneration during methanol assimilation. We revealed that the xylulose-monophosphate cycle is employing a specialized set of methanol-induced enzymes located in the peroxisome, rather than the PPP proteins, which are essentially not transcriptionally or translationally regulated in this study. For this purpose, P. pastoris has acquired a second copy of the relevant genes, each adjacent to the canonical PPP gene. During growth on methanol, the peroxisomes also harbor an active glyoxylate cycle, while the TCA cycle flux is reduced, indicating that methanol dissimilation is a major source for NADH and energy production on this substrate. Furthermore, growth on methanol/glycerol leads to a higher amino acid synthesis rate and a higher translational capacity which is reflected by a higher total protein content and may indicate a higher capacity for production of heterologous proteins as well. The observed changes in lipid metabolism can be explained by the high abundance of peroxisomes and the absence of lipid droplets in methanol-grown P. pastoris. The methylotrophic lifestyle reflects a low energy status, thus impeding lipid storage. During growth on methanol, the methanol utilization enzymes are produced in high amounts. Consequently, the biosynthetic pathways for the corresponding prosthetic groups and co-enzymes are also strongly up-regulated. Up-regulation of the pathways to riboflavin, thiamine, nicotinamide, and heme clearly indicates their high steady state demand in methanol-grown cells.
This work provides a unique data set on the methylotrophic metabolism of P. pastoris, and enables the redefinition of the methanol assimilation pathway. These findings will also have major impact on the understanding and evolution of methylotrophy in other yeasts.
Strains & chemostat cultivation
The chemostat cultivations were performed in a 1.4-L bioreactor (DASGIP Parallel Bioreactor System, Germany) with a working volume of 400 mL.
Briefly, 100 mL pre-culture medium (per liter: 10 g yeast extract, 20 g peptone, 10 g glycerol) were inoculated with 750 μL cryostock of P. pastoris CBS7435 and grown at 28 °C and 150 rpm overnight. This culture was used for inoculation of the bioreactor at an optical density (OD600) of 1.0. After a batch phase of approximately 24 h, the cells were grown in carbon-limited chemostats with a dilution rate of 0.1 h−1 for at least seven residence times before taking the samples. For each condition, three independent chemostat cultivations were performed. Temperature, pH, and dissolved oxygen were maintained at 25 °C, 5.0 (with 8 M KOH) and 20 % (by controlling the stirrer speed and inlet air), respectively.
Batch medium contained per liter: 39.9 g glycerol, 1.8 g citric acid, 12.6 g (NH4)2HPO4, 0.022 g CaCl2·2H2O, 0.9 g KCl, 0.5 g MgSO4·7H2O, 2 mL biotin (0.2 g L−1), 4.6 mL trace salts stock solution. The pH was set to 5.0 with 32 % (w/w) HCl.
Chemostat medium (Glucose) contained per liter: 55 g glucose·H2O, 2.3 g citric acid, 21.75 g (NH4)2HPO4, 0.04 g CaCl2·2H2O, 2.5 g KCl, 1.0 g MgSO4·7H2O, 2 g biotin (0.2 g L−1), and 2.43 g trace salts stock solution. The pH was set to 5.0 with 32 % (w/w) HCl.
Chemostat medium (methanol/glycerol) contained per liter: 57 g glycerol (86 %), 8.5 g methanol (100 %), 2.3 g citric acid, 21.75 g (NH4)2HPO4, 0.04 g CaCl2·2H2O, 2.5 g KCl, 1.0 g MgSO4·7H2O, 2 g biotin (0.2 g L−1), and 2.43 g trace salts stock solution. The pH was set to 5.0 with 32 % (w/w) HCl.
Trace salts stock solution contained per liter: 6.0 g CuSO4·5H2O, 0.08 g NaI, 3.0 g MnSO4·H2O, 0.2 g Na2MoO4·2H2O, 0.02 g H3BO3, 0.5 g CoCl2, 20.0 g ZnCl2, 5.0 g FeSO4·7H2O, and 5.0 mL H2SO4 (95–98 % w/w).
Sampling and quenching
For transcriptomics, 9 mL of culture were added to 4.5 mL of freshly prepared pre-chilled (−20 °C) fixing solution (5 % v/v phenol in ethanol abs.), mixed, and 1.5 mL were aliquoted into ribolyzer tubes and centrifuged at 13,000 rpm for 1 min at 4 °C. The supernatant was discarded and the tubes containing the fixed cell pellets were immediately stored at −80 °C. For protein analysis, 2 mL of culture were centrifuged and the cell pellet was stored at −80 °C. The supernatant was also stored at −80 °C for analysis of extracellular metabolites.
Samples for analysis of intracellular metabolites were taken immediately by using a pump. Approximately 50 mL fermentation broth were quenched in 200 mL of 60 % (v/v) methanol at −27 °C. After quenching, 2 mL of quenched cells (corresponding to approximately 10 mg biomass per filter) were filtered with cellulose acetate filter (0.45 μm, Satorius Biolab Products) using a vacuum pump.
The cells were washed once with cold 60 % (v/v) methanol and then the filter was kept on dry ice. Using two filtration units (Polycarbonat Filter Holders, Satorius Lab Technologies Product), 6 samples per chemostat cultivation were taken.
Biomass was determined by drying duplicates of 2 mL chemostat culture to constant weight at 105 °C in pre-weight beakers.
Total protein determination
Cell pellets from 2 mL chemostat culture were washed with 0.9 % NaCl and resuspended in 1 mL of PBS (pH 7.0). The protein extraction was done accordingly to Verduyn et al. , by addition of NaOH and incubation at 95 °C. After incubation, 0.8 M HCl were added and cell debris were collected via centrifugation. The supernatant was used for the determination of the total protein content using Bradford. The total protein content was related to the yeast dry mass (%).
Lipid extraction from P. pastoris chemostat samples was performed as described by Folch et al. . For quantitative determination of non-polar lipids (TG, SE), free fatty acids and free ergosterol, lipid extracts were loaded to Silica Gel 60 plates (Macherey-Nagel, Düren, Germany) and chromatograms were developed in an ascending manner by using the solvent system light petroleum/diethyl ether/acetic acid (70:30:2; per vol.) for approximately the first third of the distance. Subsequently, plates were briefly dried and further developed using the solvent system light petroleum/diethyl ether (49:1; per vol.) until the solvent front reached the top of the plate. Unesterified sterols and steryl esters were quantified densitrometrically using a TLC scanner (Camag TLC Scanner 3) at 275 nm using ergosterol as standard. Other lipids were irreversibly stained by dipping the TLC plates into a charring solution (0.63 g MnCl2·4H2O, 60 mL water, 60 mL methanol, and 4 mL concentrated sulfuric acid) and heated at 100 °C for 30 min. Densitrometric scanning was performed at a wavelength of 400 nm, and lipids were quantified with ergosterol, oleic acid, or triolein as standard.
For estimation of total amounts of glycerophospholipids separate bands from non-polar lipid analysis (see above) were used. Glycerophospholipids were visualized on plates by reversible staining with iodine vapor, scraped off, and subjected to quantification by the method of Broekhyuse .
Analysis of total fatty acids was achieved by conversion to methyl esters by methanolysis using 2.5 % sulfuric acid in methanol and heating at 85 °C for 90 min. FAMEs (fatty acid methyl ester) were extracted twice in a mixture of light petroleum and water (3:1; v/v.) and subjected to gas liquid chromatography (Hewlett-Packard 6890 Gas-chromatograph) using an HP-INNOWax capillary column (15 m × 0.25 mm i.d. × 0.50 μm film thickness) with helium as carrier gas. Fatty acids were identified by comparison to commercially available fatty acid methyl ester standard mix GLC-68B (NuCheck, Inc., Elysian, MN, USA) and quantified by using pentadecanoic acid (Sigma) as an internal standard.
Microarrays and data analysis
The RNA was isolated from chemostat sample cells using TRI reagent according to the supplier’s instructions (Ambion, USA). RNA integrity was analyzed using RNA nano chips (Agilent). In-house designed P. pastoris-specific oligonucleotide arrays (AMAD-ID: 034821, 8x15K custom arrays, Agilent) were used [20, 32]. cRNA synthesis, hybridization, and scanning were performed according to the Agilent protocol for two-color expression arrays. Each sample was hybridized against an RNA reference pool sample in dye swap. The microarray data was not background normalized. Normalization steps and statistical analysis of microarray data included removal of color bias using locally weighted MA-scatterplot smoothing (LOESS) followed by between array normalization using the “Aquantile” method. The P values associated with the differential expression values were calculated using a linear model fit (limma R package), subsequently they were adjusted for multiple testing using the method of Benjamini and Yekutieli  using the BY method of limma R package. To identify differentially expressed genes, the following criteria were applied: fold change cut-off of at least 1.5 > FC >1/1.5 and adjusted P value <0.05. All steps were performed using the R software package , and the limma package. Transcriptomics data were deposited at Gene Expression Omnibus with the accession number GSE67690. Data can be accessed with following link http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=stopswyszunfkf&acc=GSE67690.
Cell lysis and sample preparation
Cells were lysed with glass beads as described by Dragosits et al.  in 100 mM triethylammoniumbicarbonate (TEAB) buffer, containing 30 mM tris(2-carboxyethyl)phosphine hydrochloride and 2 % SDS. After incubation for 45 min at 56 °C (to reduce cysteine bridges) cellular proteins were extracted with chloroform/methanol, dried, dissolved in TEAB buffer, and digested with trypsin. Tandem Mass Tag (Thermo Scientific) labelling was performed as described by Pichler et al.  following the manufacturer’s protocol.
2D-LC and MS analysis
Samples were separated by high pH C18 HPLC applying an elution gradient of 12.5–80 % acetonitrile at pH 10 (200 mM ammonium formiate). Eighteen fractions were collected, partially pooled and applied to a C18 nano-column on a Bruker maxis 4G ETD QTOF LC-MS instrument, and separated with a 5–32 % acetonitrile gradient with 0.1 % formic acid (followed by a 32–80 % gradient to elute large peptides). The mass spectrometer was equipped with the captive spray source (1350 V capillary voltage, 3 L/min dry gas). Mass spectrometry scans were recorded in DDA mode (range: 150–2200 Da) and the 10 highest peaks were selected for fragmentation. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD002036.
The software Mascot was used for the identification of peptides and proteins by matching the observed spectra with a database containing unique P. pastoris protein sequences. Mascot uses the MOWSE (MOlecular Weight SEarch) score: the more matches, the higher the peptide score. Protein scores are the sum of the peptide scores. Protein identification requires the match of at least two independent peptides with a score of >25.
For quantitative analysis of the proteomics data, the software Isobar Version 1.7.5 was used . Mascot identification and quantification data were normalized using Isobar's default normalization method, which corrects for differences in reporter channel median intensities. Intensity measurement noise was corrected with a noise model comparing identical samples in multiple channels. For obtaining the protein ratios, Isobar calculates a weighted average of the peptide spectra after eliminating outliers. Comparing different distributions showed that a t-distribution fitted the random protein ratio distribution of our data best, and was selected for P value calculation of differentially expressed proteins.
Three biological replicates, with two technical replicates each, had been analyzed leading to six replicate data sets of both growth conditions. For every identified peptide, Isobar calculated the log2 of the ratio between the methanol/glycerol samples and the glucose samples (log2 FC) and the P value. Peptides with ion intensity values smaller then 300 and protein ratios deriving from single peptide spectra were excluded from the analysis, as well as proteins that were identified only in one or two replicates; 1,066 proteins fulfilled those criteria.
Proteins meeting the following criteria were defined as significantly changed between growth conditions: |mean FC| >1.5 and P values <0.1 or 1.3 < |mean FC| ≤1.5 and P value <0.05. Proteins with |mean FC| ≤1.3 and P value >0.05 were defined as not changed between the two growth conditions. To further increase stringency of evaluation we defined that >50 % of the replicates in which a given protein could be identified must have the same regulatory characteristics. Proteins that did not fulfill these criteria were not further considered. From the 1,066 proteins identified, 575 could be quantitatively evaluated.
Extraction and measurement of intracellular metabolites
For the measurement of intracellular concentrations of free metabolites quenched cells on cellulose acetate filters were used. Prior to the extraction, uniformly labelled 13C internal standard was added to the samples. Free intracellular metabolites were extracted by addition of 4 mL boiling HPLC grade ethanol (82 %; v/v; tempered at 85 °C). After addition of the boiling ethanol the quenched cells were immediately suspended by vortexing for 30 s. Suspended cells were heated for 3 min in total at 85 °C using a water bath. After 1.5 min of extraction samples were vortexed for 10 s and put back to the water bath at 85 °C. After 3 min of extractions extracted cells were immediately cooled down on dry ice. The cooled sample was then centrifuged to remove cell debris (10 min, −20 °C, 4000 g). The ethanolic extract was decanted into a fresh cooled 15 mL tube and kept on dry ice until sample preparation for LC-MS/MS and GC-MS/MS analysis. LC-MS/MS analysis of free intracellular metabolites was performed according to Klavins et al. , whereas GC-MS/MS analysis of sugar phosphates was performed after automated derivatization via ethoximation followed by trimethylsilylation. Both methods employed quantification by external calibration utilizing a uniformly 13C-labeled ethanolic extract of P. pastoris for internal standardization .
Detection of sedoheptulose-1,7-bisphosphate in cell extracts of P. pastoris
Acetonitrile, water, and formic acid (all LC-MS grade) were purchased at Sigma-Aldrich. P. pastoris cells were grown in glucose- or methanol-limited conditions in bioreactors. A set of three samples from glucose-fed cellular extracts was compared to a set of three methanol-grown samples. Each sample was derived from a separate biological replicate. After extraction (see above) the samples were stored at −80 °C until analysis. 500 μL of the sample were evaporated to dryness using a Savant SPD 121P SpeedVac Concentrator (Thermo Scientific). The residues were reconstituted in 100 μL water and directly analyzed.
Liquid chromatography separation was performed on a Hypercarb 150 × 2.1 mm, 3 μm particle size column (Thermo Scientific) with a Hypercarb guard cartridge (10 × 2.1 mm, 3 μm) using a 1260 BinPumpSL (Agilent Technologies) combined with a CTC Pal autosampler (CTC Analytics AG). The flow rate was 250 μL min−1 and the column oven was set to 40 °C. Sample injection volume was 5 μL. Mobile phase A was 100 % water, whereas mobile phase B contained 80 % acetonitrile, 10 % water, and 10 % formic acid. A gradient was applied as follows: starting conditions of 1 % B were held for 2.5 min and then increased to 40 % within 14 min. This composition was held for 1 min, before returning to 1 % B in 0.1 min for re-equilibration. The total analysis time was 20 min.
An Agilent 6220 LC-TOFMS system equipped with a dual-ESI-Source was used for the LC-MS analysis. Source parameters for negative mode were set as follows: 350 °C gas temperature, 10 L min-1 drying gas flow, 25 psig nebulizer gas pressure, 3500 V capillary voltage, 140 V fragmentor voltage, and 60 V skimmer voltage. The mass spectrometer was operated in the 2 GHz mode (extended dynamic range) recording the mass range from 50 to 1000 m/z with an acquisition rate of 1.03 spectra s-1 (9644 transients per spectrum). Data evaluation was performed using the Agilent MassHunter Qualitative Analysis B.07.00.
The identification of sedoheptulose-1,7-bisphosphate was confirmed by comparing the signals obtained in the samples to a standard which was provided by Amy A. Caudy (University of Toronto, Canada). The difference between the measured accurate mass and the calculated exact mass was below 2 ppm for all samples where sedoheptulose-1,7-bisphosphate was detected.
13C-Metabolic flux analysis
13C-labelling experiments were performed as described in Baumann et al. . The cells grew in a chemostat at a constant growth rate of 0.1 h−1 on a mixture of 20 % fully 13C-labelled substrate and 80 % naturally labelled substrate, either glucose or methanol/glycerol. The labelling pattern of protein-bound amino acids was determined via GC-MS. The GC-MS spectra were used for the calculation of mass distribution vectors of the protein bound amino acids . The Matlab-based software package Openflux was used for 13C-Metabolic flux analysis. For the calculation standard settings were applied . The flux calculation was performed with a stoichiometric model of P. pastoris central carbon metabolism. The model is analogous to the model already published by Jorda et al. . As a constraint, the labelling pattern of protein-bound amino acids and the calculated uptake and segregation rates of extracellular metabolites were used.
Isolation and proteome characterization of the peroxisomes
P. pastoris cells were cultivated on glucose (YPD) or methanol (YPM) until they reached the late logarithmic growth phase. Cellular fractions enriched of highly pure peroxisomes were isolated following the procedure which had previously been established for P. pastoris by Wriessnegger et al. . Isolated peroxisomes from methanol- and glucose-grown cells and the respective homogenates were evaluated for specific marker protein enrichment by Western blots (Additional file 3) and subjected to proteomics identification.
Samples were analyzed with a nano LC system as described above in 2D-LC and MS analysis. A standard 180 min gradient, using 0.1 % formic acid and 80 % acetonitrile as solvents, was applied. Data interpretation was performed manually (quantification) using DataAnalysis 4.0 and the files were converted to XML files for protein identification. XML files are suitable for performing a MS/MS ion search with ProteinScape (Bruker software, MASCOT embedded). At least two peptides and a MASCOT score of 30 were minimum thresholds for a positive hit. For quantification, the extracted ion chromatograms of the most intense peptides of each protein were integrated and peak areas were calculated. The sum of peak areas of each protein quantified of peroxisomal preparations was set in relation to the sum of peak areas of homogenate samples. Thus, a value higher than 1 reflects a relative enrichment in comparison to Aox1 in the peroxisomal preparation and a value lower than 1 shows a lower abundance in the peroxisomal fraction in comparison to Aox1. Peptide and protein hits, and the peak areas of peptides used for quantification are provided in Additional file 4. The mass spectrometry proteomics data were deposited in the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD002831.
Cells were cultivated at 25 °C with shaking at 150 rpm in baffled flasks using YPD until reaching the late exponential phase. Washed cells were fixed for 5 min in a 1 % aqueous solution of KMnO4 at room temperature, washed with double distilled water, and fixed again in a 1 % aqueous solution of KMnO4 for 20 min. Fixed cells were washed four times in distilled water and incubated in 0.5 % aqueous uranyl acetate overnight at 4 °C. Samples were then dehydrated for 20 min, in a graded series of 50 %, 70 %, 90 %, and 100 % ethanol, each. Pure ethanol was then changed to propylene oxide and specimens were gradually infiltrated with increasing concentrations (30 %, 50 %, 70 %s and 100 %) of Agar 100 epoxy resin mixed with propylene oxide for a minimum of 3 h per step. Samples were embedded in pure, fresh Agar 100 epoxy resin and polymerized at 60 °C for 48 h. Ultrathin 80-nm sections were stained for 3 min with lead citrate and viewed with a Philips CM 10 transmission electron microscope.
Cell dry weight
Flavin adenine dinucleotide
- H2O2 :
Pentose phosphate pathway
Peroxisomal targeting signal sequence
Unfolded protein response
Anthony C. The biochemistry of methylotrophs. London: Academic; 1982.
Meehl MA, Stadheim TA. Biopharmaceutical discovery and production in yeast. Curr Opin Biotechnol. 2014;30:120–7.
Liu L, Yang H, Shin HD, Chen RR, Li J, Du G, et al. How to achieve high-level expression of microbial enzymes: strategies and perspectives. Bioengineered. 2013;4:212–23.
Ma C, Agrawal G, Subramani S. Peroxisome assembly: matrix and membrane protein biogenesis. J Cell Biol. 2011;193:7–16.
Losev E, Reinke CA, Jellen J, Strongin DE, Bevis BJ, Glick BS. Golgi maturation visualized in living yeast. Nature. 2006;441:1002–6.
Vogl T, Glieder A. Regulation of Pichia pastoris promoters and its consequences for protein production. N Biotechnol. 2013;30:385–404.
Kohlwein SD, Veenhuis M, van der Klei IJ. Lipid droplets and peroxisomes: key players in cellular lipid homeostasis or a matter of fat--store ‘em up or burn’em down. Genetics. 2013;193:1–50.
van der Klei IJ, Veenhuis M. Yeast peroxisomes: function and biogenesis of a versatile cell organelle. Trends Microbiol. 1997;5:502–9.
van der Klei IJ, Yurimoto H, Sakai Y, Veenhuis M. The significance of peroxisomes in methanol metabolism in methylotrophic yeast. Biochim Biophys Acta. 1763;2006:1453–62.
van Zutphen T, Baerends RJ, Susanna KA, de Jong A, Kuipers OP, Veenhuis M, et al. Adaptation of Hansenula polymorpha to methanol: a transcriptome analysis. BMC Genomics. 2010;11:1.
Nagotu S, Krikken AM, Otzen M, Kiel JA, Veenhuis M, van der Klei IJ. Peroxisome fission in Hansenula polymorpha requires Mdv1 and Fis1, two proteins also involved in mitochondrial fission. Traffic. 2008;9:1471–84.
Hartner FS, Glieder A. Regulation of methanol utilisation pathway genes in yeasts. Microb Cell Fact. 2006;5:39.
Yurimoto H, Oku M, Sakai Y. Yeast methylotrophy: metabolism, gene regulation and peroxisome homeostasis. Int J Microbiol. 2011;2011:101298.
Vanz AL, Lunsdorf H, Adnan A, Nimtz M, Gurramkonda C, Khanna N, et al. Physiological response of Pichia pastoris GS115 to methanol-induced high level production of the Hepatitis B surface antigen: catabolic adaptation, stress responses, and autophagic processes. Microb Cell Fact. 2012;11:103.
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:23.
Liang S, Wang B, Pan L, Ye Y, He M, Han S, et al. Comprehensive structural annotation of Pichia pastoris transcriptome and the response to various carbon sources using deep paired-end RNA sequencing. BMC Genomics. 2012;13:738.
Sauer M, Branduardi P, Gasser B, Valli M, Maurer M, Porro D, et al. Differential gene expression in recombinant Pichia pastoris analysed by heterologous DNA microarray hybridisation. Microb Cell Fact. 2004;3:17.
Jorda J, Rojas HC, Carnicer M, Wahl A, Ferrer P, Albiol J. Quantitative metabolomics and instationary 13C-metabolic flux analysis reveals impact of recombinant protein production on Trehalose and energy metabolism in Pichia pastoris. Metabolites. 2014;4:281–99.
Rebnegger C, Graf AB, Valli M, Steiger MG, Gasser B, Maurer M, et al. In Pichia pastoris, growth rate regulates protein synthesis and secretion, mating and stress response. Biotechnol J. 2014;9:511–25.
Graf A, Gasser B, Dragosits M, Sauer M, Leparc G, Tuechler T, et al. Novel insights into the unfolded protein response using Pichia pastoris specific DNA microarrays. BMC Genomics. 2008;9:390.
Baumann K, Carnicer M, Dragosits M, Graf AB, Stadlmann J, Jouhten P, et al. A multi-level study of recombinant Pichia pastoris in different oxygen conditions. BMC Syst Biol. 2010;4:141.
Dragosits M, Stadlmann J, Graf A, Gasser B, Maurer M, Sauer M, et al. The response to unfolded protein is involved in osmotolerance of Pichia pastoris. BMC Genomics. 2010;11:207.
Hesketh AR, Castrillo JI, Sawyer T, Archer DB, Oliver SG. Investigating the physiological response of Pichia (Komagataella) pastoris GS115 to the heterologous expression of misfolded proteins using chemostat cultures. Appl Microbiol Biotechnol. 2013;97:9747–62.
Dragosits M, Stadlmann J, Albiol J, Baumann K, Maurer M, Gasser B, et al. The effect of temperature on the proteome of recombinant Pichia pastoris. J Proteome Res. 2009;8:1380–92.
Lin XQ, Liang SL, Han SY, Zheng SP, Ye YR, Lin Y. Quantitative iTRAQ LC-MS/MS proteomics reveals the cellular response to heterologous protein overexpression and the regulation of HAC1 in Pichia pastoris. J Proteomics. 2013;91:58–72.
Chung BK, Selvarasu S, Andrea C, Ryu J, Lee H, Ahn J, et al. Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement. Microb Cell Fact. 2010;9:50.
Sohn SB, Graf AB, Kim TY, Gasser B, Maurer M, Ferrer P, et al. Genome-scale metabolic model of methylotrophic yeast Pichia pastoris and its use for in silico analysis of heterologous protein production. Biotechnol J. 2010;5:705–15.
Jorda J, Jouhten P, Camara E, Maaheimo H, Albiol J, Ferrer P. Metabolic flux profiling of recombinant protein secreting Pichia pastoris growing on glucose:methanol mixtures. Microb Cell Fact. 2012;11:57.
Jorda J, Suarez C, Carnicer M, ten Pierick A, Heijnen JJ, van Gulik W, et al. Glucose-methanol co-utilization in Pichia pastoris studied by metabolomics and in stationary 13C flux analysis. BMC Syst Biol. 2013;7:17.
Jorda J, de Jesus SS, Peltier S, Ferrer P, Albiol J. Metabolic flux analysis of recombinant Pichia pastoris growing on different glycerol/methanol mixtures by iterative fitting of NMR-derived 13C-labelling data from proteinogenic amino acids. N Biotechnol. 2014;31:120–32.
Solà A, Jouhten P, Maaheimo H, Sánchez-Ferrando F, Szyperski T, Ferrer P. Metabolic flux profiling of Pichia pastoris grown on glycerol/methanol mixtures in chemostat cultures at low and high dilution rates. Microbiology. 2007;153:281–90.
Prielhofer R, Maurer M, Klein J, Wenger J, Kiziak C, Gasser B, et al. Induction without methanol: novel regulated promoters enable high-level expression in Pichia pastoris. Microb Cell Fact. 2013;12(1):5.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13:227–32.
Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007;25:117–24.
Lee MV, Topper SE, Hubler SL, Hose J, Wenger CD, Coon JJ, et al. A dynamic model of proteome changes reveals new roles for transcript alteration in yeast. Mol Syst Biol. 2011;7:514.
Kiel JA, Veenhuis M, van der Klei IJ. PEX genes in fungal genomes: common, rare or redundant. Traffic. 2006;7:1291–303.
Prielhofer R, Cartwright SP, Graf AB, Valli M, Bill RM, Mattanovich D, et al. Pichia pastoris regulates its gene-specific response to different carbon sources at the transcriptional, rather than the translational, level. BMC Genomics. 2015;16:167.
Terlecky SR, Nuttley WM, McCollum D, Sock E, Subramani S. The Pichia pastoris peroxisomal protein PAS8p is the receptor for the C-terminal tripeptide peroxisomal targeting signal. EMBO J. 1995;14:3627–34.
Faber KN, Haima P, Gietl C, Harder W, Ab G, Veenhuis M. The methylotrophic yeast Hansenula polymorpha contains an inducible import pathway for peroxisomal matrix proteins with an N-terminal targeting signal (PTS2 proteins). Proc Natl Acad Sci U S A. 1994;91:12985–9.
Elgersma Y, Elgersma-Hooisma M, Wenzel T, McCaffery JM, Farquhar MG, Subramani S. A mobile PTS2 receptor for peroxisomal protein import in Pichia pastoris. J Cell Biol. 1998;140:807–20.
Küberl A, Schneider J, Thallinger GG, Anderl I, Wibberg D, Hajek T, et al. High-quality genome sequence of Pichia pastoris CBS7435. J Biotechnol. 2011;154:312–20.
Neuberger G, Maurer-Stroh S, Eisenhaber B, Hartig A, Eisenhaber F. Prediction of peroxisomal targeting signal 1 containing proteins from amino acid sequence. J Mol Biol. 2003;328:581–92.
PTS1 Predictor. http://mendel.imp.ac.at/mendeljsp/sat/pts1/PTS1predictor.jsp. Last access April 28, 2015.
Wriessnegger T, Gübitz G, Leitner E, Ingolic E, Cregg J, de la Cruz B, et al. Lipid composition of peroxisomes from the yeast Pichia pastoris grown on different carbon sources. Biochim Biophys Acta. 1771;2007:455–61.
Luers GH, Advani R, Wenzel T, Subramani S. The Pichia pastoris dihydroxyacetone kinase is a PTS1-containing, but cytosolic, protein that is essential for growth on methanol. Yeast. 1998;14:759–71.
Clasquin MF, Melamud E, Singer A, Gooding JR, Xu X, Dong A, et al. Riboneogenesis in yeast. Cell. 2011;145:969–80.
Douma AC, Veenhuis M, de Koning W, Evers M, Harder W. Dihydroxyacetone synthase is localized in the peroxisomal matrix of methanol-grown Hansenula polymorpha. Arch Microbiol. 1985;143:237–43.
Anderson LE, Carol AA. Enzyme co-localization with rubisco in pea leaf chloroplasts. Photosynth Res. 2004;82:49–58.
Rae BD, Long BM, Whitehead LF, Forster B, Badger MR, Price GD. Cyanobacterial carboxysomes: microcompartments that facilitate CO2 fixation. J Mol Microbiol Biotechnol. 2013;23:300–7.
Hahn MW. Distinguishing among evolutionary models for the maintenance of gene duplicates. J Hered. 2009;100:605–17.
Zhang J. Evolution by gene duplication: an update. Trends Ecol Evol. 2003;18:292–8.
Byun-McKay SA, Geeta R. Protein subcellular relocalization: a new perspective on the origin of novel genes. Trends Ecol Evol. 2007;22:338–44.
Daran-Lapujade P, Rossell S, van Gulik WM, Luttik MA, de Groot MJ, Slijper M, et al. The fluxes through glycolytic enzymes in Saccharomyces cerevisiae are predominantly regulated at posttranscriptional levels. Proc Natl Acad Sci U S A. 2007;104:15753–8.
Tibbetts AS, Sun Y, Lyon NA, Ghrist AC, Trotter PJ. Yeast mitochondrial oxodicarboxylate transporters are important for growth on oleic acid. Arch Biochem Biophys. 2002;406:96–104.
Cavero S, Vozza A, del Arco A, Palmieri L, Villa A, Blanco E, et al. Identification and metabolic role of the mitochondrial aspartate-glutamate transporter in Saccharomyces cerevisiae. Mol Microbiol. 2003;50:1257–69.
Kunze M, Pracharoenwattana I, Smith SM, Hartig A. A central role for the peroxisomal membrane in glyoxylate cycle function. Biochim Biophys Acta. 1763;2006:1441–52.
Warner JR, Mitra G, Schwindinger WF, Studeny M, Fried HM. Saccharomyces cerevisiae coordinates accumulation of yeast ribosomal proteins by modulating mRNA splicing, translational initiation, and protein turnover. Mol Cell Biol. 1985;5:1512–21.
Puxbaum V, Mattanovich D, Gasser B. Quo vadis? The challenges of recombinant protein folding and secretion in Pichia pastoris. Appl Microbiol Biotechnol. 2015;99:2925–38.
Hohenblum H, Gasser B, Maurer M, Borth N, Mattanovich D. Effects of gene dosage, promoters, and substrates on unfolded protein stress of recombinant Pichia pastoris. Biotechnol Bioeng. 2004;85:367–75.
Resina D, Bollók M, Khatri N, Valero F, Neubauer P, Ferrer P. Transcriptional response of P. pastoris in fed-batch cultivations to Rhizopus oryzae lipase production reveals UPR induction. Microb Cell Fact. 2007;6:21.
Ozimek P, van Dijk R, Latchev K, Gancedo C, Wang DY, van der Klei IJ, et al. Pyruvate carboxylase is an essential protein in the assembly of yeast peroxisomal oligomeric alcohol oxidase. Mol Biol Cell. 2003;14:786–97.
Stewart MQ, Esposito RD, Gowani J, Goodman JM. Alcohol oxidase and dihydroxyacetone synthase, the abundant peroxisomal proteins of methylotrophic yeasts, assemble in different cellular compartments. J Cell Sci. 2001;114:2863–8.
Waterham HR, Russell KA, Vries Y, Cregg JM. Peroxisomal targeting, import, and assembly of alcohol oxidase in Pichia pastoris. J Cell Biol. 1997;139:1419–31.
Gudipati V, Koch K, Lienhart WD, Macheroux P. The flavoproteome of the yeast Saccharomyces cerevisiae. Biochim Biophys Acta. 1844;2014:535–44.
Marx H, Mattanovich D, Sauer M. Overexpression of the riboflavin biosynthetic pathway in Pichia pastoris. Microb Cell Fact. 2008;7:23.
Stadlmayr G, Mecklenbrauker A, Rothmuller M, Maurer M, Sauer M, Mattanovich D, et al. Identification and characterisation of novel Pichia pastoris promoters for heterologous protein production. J Biotechnol. 2010;150:519–29.
Sporty J, Lin SJ, Kato M, Ognibene T, Stewart B, Turteltaub K, et al. Quantitation of NAD+ biosynthesis from the salvage pathway in Saccharomyces cerevisiae. Yeast. 2009;26:363–9.
Wriessnegger T, Sunga AJ, Cregg JM, Daum G. Identification of phosphatidylserine decarboxylases 1 and 2 from Pichia pastoris. FEMS Yeast Res. 2009;9:911–22.
Verduyn C, Postma E, Scheffers WA, van Dijken JP. Physiology of Saccharomyces cerevisiae in anaerobic glucose-limited chemostat cultures. J Gen Microbiol. 1990;136:395–403.
Folch J, Lees M, Sloane Stanley GH. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957;226:497–509.
Broekhuyse RM. Phospholipids in tissues of the eye. I. Isolation, characterization and quantitative analysis by two-dimensional thin-layer chromatography of diacyl and vinyl-ether phospholipids. Biochim Biophys Acta. 1968;152:307–15.
Benjamini Y, Heller R, Yekutieli D. Selective inference in complex research. Philos Trans A Math Phys Eng Sci. 2009;367:4255–71.
R project. http://www.r-project.org. Last access January 31, 2012.
Pichler P, Kocher T, Holzmann J, Mazanek M, Taus T, Ammerer G, et al. Peptide labeling with isobaric tags yields higher identification rates using iTRAQ 4-plex compared to TMT 6-plex and iTRAQ 8-plex on LTQ Orbitrap. Anal Chem. 2010;82:6549–58.
Breitwieser FP, Muller A, Dayon L, Kocher T, Hainard A, Pichler P, et al. General statistical modeling of data from protein relative expression isobaric tags. J Proteome Res. 2011;10:2758–66.
Klavins K, Neubauer S, Al Chalabi A, Sonntag D, Haberhauer-Troyer C, Russmayer H, et al. Interlaboratory comparison for quantitative primary metabolite profiling in Pichia pastoris. Anal Bioanal Chem. 2013;405:5159–69.
Neubauer S, Haberhauer-Troyer C, Klavins K, Russmayer H, Steiger MG, Gasser B, et al. U13C cell extract of Pichia pastoris--a powerful tool for evaluation of sample preparation in metabolomics. J Sep Sci. 2012;35:3091–105.
Zamboni N, Fischer E, Sauer U. FiatFlux--a software for metabolic flux analysis from 13C-glucose experiments. BMC Bioinformatics. 2005;6:209.
Quek LE, Wittmann C, Nielsen LK, Kromer JO. OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis. Microb Cell Fact. 2009;8:25.
This work has been supported by the Federal Ministry of Science, Research and Economy (BMWFW), the Federal Ministry of Traffic, Innovation and Technology (bmvit), the Styrian Business Promotion Agency SFG, the Standortagentur Tirol, the Government of Lower Austria and ZIT – Technology Agency of the City of Vienna through the COMET-Funding Program managed by the Austrian Research Promotion Agency FFG. Further support by Biomin Research Center, Boehringer-Ingelheim RCV, Lonza AG, Biocrates Life Sciences AG, VTU Technology GmbH, and Sandoz GmbH is acknowledged. EQ BOKU VIBT GmbH is acknowledged for providing mass spectrometry and bioreactor instrumentation. The authors are grateful to Amy A. Caudy (University of Toronto, Canada) for kindly providing a standard of seduheptulose bisphosphate and thank Jonas Hohlweg for excellent support in transcriptomics analyses.
The authors declare that they have no competing interests.
HR and MB performed all cultivations, sampling, and transcriptomics analyses. CG measured quantitative proteomes, supervised by FA. MV coordinated data management and evaluation, and analyzed proteomics data together with GM and CG. KHG performed lipid analysis and purification of peroxisomes. RG, KK, SN, and AC measured the metabolomes under the guidance of GK, SH, and DS. CT measured 13C labelling patterns for metabolic flux analysis which was calculated by HR. HD analyzed S1,7BP with support by SH. GZ provided electron microscopy. ABG supervised and participated in evaluation of transcriptomics and proteomics data. DM and BG designed this study with support of MS, MSt, and GD. Systems level data interpretation was performed by HR, MB, MV, KHG, GD, MSt, MS, DM, and BG. HR, MB, MV, and KHG drafted the manuscript, which was revised by BG and DM. All authors read and approved the final manuscript.
Hannes Rußmayer and Markus Buchetics contributed equally to this work.
Transcriptomic, proteomic, and metabolomic regulation of P. pastoris during methylotrophic growth. Containing the following eight sheets: Summary Omics Data: number of significantly regulated genes, proteins or metabolites (e.g. “up” refers to up-regulation in methanol/glycerol compared to glucose). Transcriptomics and proteomics: Average fold changes and P values of transcriptomics and proteomics comparing P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat. Average values derive from three biological replicates per condition. Metabolomics: Average fold changes and P values of metabolomics measurements comparing P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat cultivations. Average values derive from three biological replicates per condition. Co-regulation (related to Fig. 1 in the text): Regulation of the 575 gene-protein pairs with transcriptomics and proteomics data available and assignment to regulatory groups. Central carbon metabolism (related to Fig. 4 in the text): Average fold changes and P values of transcriptomics, proteomics, and metabolomics measurement depicted in Fig. 4. Amino acid metabolism (related to Fig. 6 in the text): Average fold changes and P values of transcriptomics, proteomics, and metabolomics measurement depicted in Fig. 6. Vitamin biosynthesis (related to Fig. 7 in the text): Average fold changes and P values of transcriptomics, proteomics, and metabolomics measurement depicted in Fig. 7. Peroxisomal gene regulation: Average fold changes and P values of transcriptomics and proteomics for all mentioned peroxisomal genes. Average values derive from three biological replicates per condition. (XLSX 2348 kb)
Comparison of gene regulation in P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat to transcriptomics data obtained in fed batch cultivation with methanol glycerol or glucose as carbon source. Average fold changes and P values of regulation patterns for methanol vs glycerol and methanol vs glucose 1 h our after starting the methanol feed are shown. Average values derive from three biological replicates per condition. Containing the following four sheets: Description: Experimental setup of the fed batch cultivations, description, and discussion of observed similarities and differences in gene regulations. Upregulated_methanol: Comparison of all genes upregulated in P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat to transcriptomics data obtained by analyzing methanol vs glycerol and methanol vs glucose fed batches. Downregulated_methanol: Comparison of all genes downregulated in P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat to transcriptomics data obtained by analyzing methanol vs glycerol and methanol vs glucose fed batches. All_data: Transcriptomics data of Additional file 1 (P. pastoris cultivated with methanol/glycerol or glucose as carbon source in chemostat) compared to transcriptomics data obtained by analyzing methanol vs glycerol and methanol vs glucose fed batches. (XLSX 1178 kb)
Enrichment of the peroxisomal marker protein Pex3p in the peroxisomal fraction. (PDF 271 kb)
Proteomic identification and quantification of methanol metabolic enzymes and control proteins in peroxisomal fractions and homogenates of P. pastoris cells grown on methanol. Containing the following three sheets: Protein hits: contains all identified proteins that met the threshold in at least one sample, with their respective MASCOT scores, number of peptides, and percent sequence coverage. Peptide hits: list of all identified peptides, their MASCOT scores, mass and charge values, and intensities. Peptides used for quant + areas: lists all peptides of the proteins in Table 3 that were used for quantification, and their respective peak areas in the different samples. (XLSX 879 kb)