Adaptation of Hansenula polymorpha to methanol: a transcriptome analysis
© van Zutphen et al. 2010
Received: 15 September 2009
Accepted: 4 January 2010
Published: 4 January 2010
Skip to main content
© van Zutphen et al. 2010
Received: 15 September 2009
Accepted: 4 January 2010
Published: 4 January 2010
Methylotrophic yeast species (e.g. Hansenula polymorpha, Pichia pastoris) can grow on methanol as sole source of carbon and energy. These organisms are important cell factories for the production of recombinant proteins, but are also used in fundamental research as model organisms to study peroxisome biology. During exponential growth on glucose, cells of H. polymorpha typically contain a single, small peroxisome that is redundant for growth while on methanol multiple, enlarged peroxisomes are present. These organelles are crucial to support growth on methanol, as they contain key enzymes of methanol metabolism.
In this study, changes in the transcriptional profiles during adaptation of H. polymorpha cells from glucose- to methanol-containing media were investigated using DNA-microarray analyses.
Two hours after the shift of cells from glucose to methanol nearly 20% (1184 genes) of the approximately 6000 annotated H. polymorpha genes were significantly upregulated with at least a two-fold differential expression. Highest upregulation (> 300-fold) was observed for the genes encoding the transcription factor Mpp1 and formate dehydrogenase, an enzyme of the methanol dissimilation pathway. Upregulated genes also included genes encoding other enzymes of methanol metabolism as well as of peroxisomal β-oxidation.
A moderate increase in transcriptional levels (up to 4-fold) was observed for several PEX genes, which are involved in peroxisome biogenesis. Only PEX11 and PEX32 were higher upregulated. In addition, an increase was observed in expression of the several ATG genes, which encode proteins involved in autophagy and autophagy processes. The strongest upregulation was observed for ATG8 and ATG11.
Approximately 20% (1246 genes) of the genes were downregulated. These included glycolytic genes as well as genes involved in transcription and translation.
Transcriptional profiling of H. polymorpha cells shifted from glucose to methanol showed the expected downregulation of glycolytic genes together with upregulation of the methanol utilisation pathway. This serves as a confirmation and validation of the array data obtained. Consistent with this, also various PEX genes were upregulated. The strong upregulation of ATG genes is possibly due to induction of autophagy processes related to remodeling of the cell architecture required to support growth on methanol. These processes may also be responsible for the enhanced peroxisomal β-oxidation, as autophagy leads to recycling of membrane lipids. The prominent downregulation of transcription and translation may be explained by the reduced growth rate on methanol (td glucose 1 h vs td methanol 4.5 h).
Hansenula polymorpha is an important cell factory for the production of pharmaceutical proteins . Moreover, it is extensively used in fundamental research aiming at understanding the molecular principles of peroxisome biology .
As cell factory, H. polymorpha has several important advantages. First, it contains very strong and inducible promoters derived from the methanol metabolism pathway. Also, the organism is thermotolerant (it can grow at high temperatures up to 48°C, ) and tolerates various environmental stresses. H. polymorpha does not hyperglycosylate secreted proteins, which often is a problem in heterologous protein production in S. cerevisiae.
In H. polymorpha peroxisomes massively develop during growth on methanol as sole source of carbon and energy. Methanol is oxidized by the enzyme alcohol oxidase (AOX), which is localized in peroxisomes together with catalase and dihydroxyacetone synthase (DHAS), the first enzyme of the assimilation pathway. Peroxisomes are not required for primary metabolism when cells are grown on glucose. Moreover, glucose represses the key enzymes of methanol metabolism AOX and DHAS. Therefore, during growth on glucose H. polymorpha cells contain only a single, small peroxisome. Upon a shift to methanol media, the enzymes of methanol metabolism are induced paralleled by an increase in peroxisome size and abundance. The initial small peroxisome serves as the target organelle for the enzymes of methanol metabolism and proliferates by fission . Ultimately, in exponentially growing cells, each cell contains several enlarged peroxisomes .
A wealth of information is now available of individual genes encoding enzymes of the methanol metabolism as well as on genes involved in peroxisome formation (PEX genes). However, genomics approaches in H. polymorpha are still rare.
We used genome-wide transcriptional profiling to dissect the initial events accompanying the adaptation of glucose grown H. polymorpha cells to methanol metabolism. This will gain information on the induction and repression of metabolic genes as well as on non-metabolic genes, including PEX genes.
All experiments described in this paper were performed in batch cultures. We chose not to grow the cells in carbon-limited chemostats, as glucose-limitation results in derepression of genes involved in methanol metabolism .
Replicates were obtained by growing 4 independent cultures on glucose that were independently transferred to fresh media containing methanol. Of each culture, mRNA isolated from the glucose and the methanol sample was labeled (and dye-swapped) and hybridized on two arrays per culture. In addition, as a control AOX transcript levels of these samples were determined by RT-PCR, confirming the absence of transcript in the glucose-grown pre-cultures and the presence of AOX transcript after 2 hours incubation (data not shown).
The DNA microarray analyses data were analyzed to generate the ratio between the transcripts on methanol and glucose for each gene to identify any differential expression and to determine the p-value to assess the significance and the A-value to check the intensity of the signals. [Additional file 1 : Supplemental table S1] presents an overview of the array results. Of the nearly 6000 annotated H. polymorpha genes that are listed, approximately 20% (1184 genes) are upregulated, while another 20% (1246 genes) are downregulated with at least a two-fold differential expression, meeting the significance and signal intensity criteria.
Of the upregulated genes, 13 are more than 100 times upregulated, 192 genes show a 10-100-fold upregulation, 156 genes increase between 5 and 10-fold and the remaining 823 genes are less than 5-fold upregulated. Highest upregulated are the central methanol metabolism regulator MPP1 (394-fold, Hp27g360, see below) and the gene encoding formate dehydrogenase (347-fold), required in methanol catabolism. Also the other components of the methanol metabolic pathway are highly upregulated. Moreover, CRC1 is highly induced, encoding a mitochondrial inner membrane carnitine transporter, which is required for the transport of acetyl-CoA from peroxisomes to mitochondria during fatty acid beta-oxidation (111-fold). In line with CRC1, also the genes involved in the beta-oxidation of fatty acids are overrepresented among the highly upregulated genes (for details see below). Furthermore, approximately 13% of the more than 10-fold upregulated genes is involved in transport. The upregulation of hexose transporters may be important for the uptake of the residual glucose that was present in the inoculum. Of the downregulated genes, the highest fold downregulation (65-fold) is observed for Hp24g956, encoding a protein with strong similarity to Sik1p of S. cerevisiae, which is involved in pre-rRNA processing. This predicted function is consistent with the general trend among the downregulated genes, since of the 179 genes that are over 10-fold downregulated, nearly 50% code for products that function in either transcription or translation processes. Of the other downregulated genes, 269 show a 5- to 10-fold downregulation. Of these, nearly 40% encode proteins involved in transcription and translation. The other 789 genes are less than 5-fold reduced in transcript levels.
As expected, genes involved in metabolic pathways strongly contribute to both the up- and downregulated genes (20% and 15.5%, respectively), reflecting the large-scale adaptations accompanying the shift from glucose to methanol. However, metabolism is a large group also in the total genome and the contribution in percentages does not reflect the nature of the metabolic pathways that are regulated (see below).
In contrast to metabolism, some other functional categories display a difference in contribution to the upregulated compared to the downregulated genes. One such functional category is the group of protein synthesis genes, which is almost absent among upregulated genes (0.25%), while it composes a large portion of the downregulated genes (12.1%). Of the total genome of H. polymorpha only approximately 4% is involved in protein synthesis, reflecting the considerable effect of a shift to methanol on protein synthesis. In addition, also the group of genes involved in transcription is more predominant among downregulated genes (12.7% versus 6.3% of the upregulated genes and an intermediate 9% of the total genome). In concurrence with the trend of genes related to transcription and translation, also genes related to nucleotide biosynthesis are mostly downregulated (42 of 51 genes), yet genes involved in amino acid biosynthesis show a less clear trend (30 down-, 9 upregulated, 54 not differentially expressed). The observed downregulation of major anabolic processes most likely is related to the reduction in doubling time (td of cells on methanol relative to growth on glucose (td methanol = 4.5 h, td glucose = 1 h) and may reflect the accompanying decrease in DNA replication, RNA transcription and translation.
Stress response genes form only relatively small categories among both the upregulated genes relative to the downregulated genes (2.9% vs 1.6%). Based on the Functional Catalogue, only 185 of the nearly 6000 annotated H. polymorpha genes are indicated as stress response genes. However, based on several studies by Gasch in Saccharomyces cerevisiae [9–11], many more genes could contribute to the cellular stress response. Hence, most likely also genes classified in other groups play a role in coping with stress accompanying a shift in cultivation conditions. Thus, the contribution of the stress response to the total differential expression in H. polymorpha upon transfer to methanol medium is probably larger than the observed 2.9% upregulation and 1.6% downregulation.
A last remarkable group in the Functional Catalogue diagram is the category of unclassified proteins, showing that 33% of the upregulated genes and 17.9% of the downregulated genes are thus far not experimentally characterized, relative to 25% of the genes of the total genome of H. polymorpha. This observation suggests that our current knowledge on adaptation to methanol is far from complete.
Promoter studies in Candida boidinii using phosphatase as a reporter enzyme revealed that upon a shift to methanol medium FMD was induced first, followed by DHAS and even later AOX . The early induction of FMD (347-fold up 2 hours after the shift to methanol medium) relative to AOX and DHAS (17.3-, 19-fold up respectively) suggests that a similar induction pattern may exist in H. polymorpha. The differences in induction of the genes 2 hours after the shift to methanol medium does not reflect the ultimate protein levels in the cells, as AOX and DHAS are the major protein bands in extracts prepared from methanol grown H. polymorpha cells [12, 13].
Expression changes of PEX genes
H. polymorpha PEX25 was upregulated 2.2-fold, whereas PEX11C, whose function is still unknown, showed a 1.6-fold downregulation .
Expression changes of glycolysis and gluconeogenesis genes
Phosphofructokinase alpha subunit
Phosphofructokinase beta subunit
Fructose 1,6-bisphosphate aldolase
Phosphopyruvate hydratase (enolase)
Pyruvate dehydrogenase, alpha subunit
Pyruvate dehydrogenase, beta subunit
Fructose 1,6 bisphosphatase
Accompanying the changes in expression of many metabolic genes, also changes in the underlying regulatory networks are expected. In addition to global changes, the DNA microarray data reflect the initial adaptation to methanol and thus enable the investigation of more specific changes resulting in activation of methanol-dependent genes or in repression of glucose-dependent genes.
Among the upregulated genes, 49 are involved in regulation of transcription. Most of these encode general transcription factors or transcription factors which have not been linked to a specific process. Regulators involved in stress response (6), glucose sensing/derepression (4), and peroxisome-related pathways (4) are overrepresented, as is expected due to the change in carbon source.
The expression of genes coding for peroxisome assembly and function is controlled by a transcriptional regulatory network, which has been studied extensively in S. cerevisiae in response to oleate [17–19]. The Oaf1-Pip2 complex plays a prominent role, although the putative H. polymorpha homolog of PIP2 is only moderately upregulated during adaptation to methanol. However the homolog of ADR1 (23.7-fold), a second activator of oleate-responsive genes is strongly upregulated in H. polymorpha and is also involved in regulation of the response to both oleate and methanol in P. pastoris (named MXR1; . Virtually all known targets of Adr1 and its co-regulator Cat8 were indeed also upregulated in H. polymorpha, suggesting an important role in regulation of methanol metabolism, while most probably also activation by Snf1 is initiated after the shift, since this global regulator is crucial for growth on non-fermentable carbon sources [21, 22]. Mpp1, another transcriptional regulator of methanol metabolism in H. polymorpha, is encoded by the strongest upregulated gene of this study (Hp27g360, 394-fold up), thus suggesting it is a master regulator of methanol-responsive genes .
Swi1 and Snf2 also belong to a regulatory complex involved in gene expression of methanol-related genes as well as peroxisomal assembly, however their encoding genes were not induced in the early adaptation phase to methanol .
Among the downregulated genes, the decreased transcription of RAP1 (Hp16g154, -3.1) is remarkable. This transcriptional regulator is known to activate transcription of genes encoding ribosomal proteins [25, 26] and its downregulation is consistent with the observed massive decrease in transcripts for ribosomal proteins. Interestingly, this gene is also shown to be repressed during the environmental stress response in S. cerevisiae as described by Gasch et al., .
Expression changes of ATG genes
Ultrastructural analysis of cells at different time-points after the shift indeed showed strong morphological evidence for increased autophagy, reflected in the massive uptake of cytoplasmic components in the vacuole (figure 4).
Recent findings showed the importance of autophagy during methanol adaptation of P. pastoris, not only for cell remodeling, but also to provide amino acids . Consistent with these findings, we also observed that H. polymorpha atg mutants showed a slight delay in methanol adaptation (data not shown).
Expression changes of genes related to fatty acid β-oxidation
Peroxisomal ABC-transporter sub-unit 1
Peroxisomal ABC-transporter sub-unit 2
Adenine nucleotide transporter
Fatty-acyl coenzyme A oxidase
Carnitine acetyl-CoA transferase
Carnitine acetyltransferase, YAT1
Carnitine acetyltransferase, YAT2
Mitochondrial carnitine/acyl carnitine carrier
Although mitochondria were long considered as the main source of reactive oxygen species (ROS), also peroxisomes actually defined as organelles that harbour H2O2-producing enzymes as well as catalase are now recognized as a significant contributor to ROS production [2, 35, 36]. Besides catalase (42.8-fold up), peroxisomes also contain the peroxiredoxin Pmp20, that is 11-fold upregulated, both indicative of an increase in peroxisomal ROS production [37, 38].
We also observed an increase of several pivotal genes involved in cytosolic and mitochondrial ROS detoxification; like the superoxide dismutase SOD1 (4.3-fold), the glutathione transferase GTO2 (2.6-fold), the glutathione reductase GPX3 (8-fold), the glutaredoxin GRX5 (2.8-fold) and the thioredoxin peroxidase TSA1 (47.2-fold). The remaining members of the glutathione- and TRX-based system are not differentially expressed (8 other genes found), except for the thioredoxin reductase TRR1 (10.4-fold). Induction of these key enzymes involved in sustaining the redox balance of the cytosol, suggests that methanol-metabolism also results in enhanced cytosolic ROS levels, which may originate from peroxisomal metabolism.
Back in 1996, the proof of principle for the use of DNA microarray technology to investigate transcriptional changes was shown for S. cerevisiae . Since then, DNA microarray analysis has become a regular, well-established tool, facilitated by the now commercial available array slides. For many other yeast species however, thusfar the usage of DNA microarray analysis is still rather limited. Only recently, species-specific DNA microarray studies have been presented for e.g. Debaryomyces hansenii  and P. pastoris . For H. polymorpha, Oh et al.,  presented a partial, DNA microarray study, followed by a full DNA microarray analysis describing the transcriptional response of H. polymorpha to exposure to cadmium . However, this study was not related to the metabolic and architectural alterations associated with a change in cellular metabolism.
Smith et al. [15, 18] published two studies in which bakers' yeast was shifted to oleate. The first study focused mainly on transcriptional changes of genes encoding peroxisomal proteins and peroxins, the second one on the regulatory network coordinating the response to oleate.
The current study is the first in which arrays were used to study the shift from glucose to methanol in the yeast H. polymorpha. Sauer et al.  was the first to publish transcriptional profiling of the methylotrophic yeast P. pastoris, also upon a shift to methanol. However, in that study, yeast cells were transferred from glycerol to methanol and heterologous hybridisation onto S. cerevisiae DNA-microarrays was performed. Still, the same trend in regulated functional groups was observed. Also, similar effects were reported on carbohydrate metabolism and the regulation of ribosomal genes. Although some of the data obtained for H. polymorpha are comparable with those obtained for P. pastoris, a significant drawback of the latter study is the absence of P. pastoris specific or methylotrophic yeast specific genes on the DNA-microarray slides. There is indeed a substantial amount of methanol-responsive H. polymorpha genes found in our study which were not observed in P. pastoris using S. cerevisiae microarrays (± 450 upregulated, ± 350 downregulated) [Additional file 2 : Supplemental Table S2].
The current DNA microarray study revealed the expected repression of genes involved in glucose metabolism concomitant with the induction of the genes of methanol metabolism in response of a shift of H. polymorpha cells from glucose to methanol. Also genes involved in peroxisome biogenesis (PEX genes) are upregulated, with PEX32 being the most strongly upregulated PEX gene. A surprising finding was the upregulation of autophagy- and of beta-oxidation-genes. The first can likely be explained by the need for cellular reorganisations and is confirmed by electron microscopy studies showing active uptake of cytoplasmic components in the vacuoles of the cells. The induction of beta-oxidation is thought to be a consequence of the cellular reorganisations and thereby the high turnover of lipids, serving as substrates for beta-oxidation. A final interesting but not yet uncovered group consists of the 33% of the upregulated genes that have no known function. These genes reflect a large potential of H. polymorpha or methylotrophic yeast specific genes with a specialized role in adaptations to growth on methanol as the sole source of carbon and can form interesting targets for future research.
H. polymorpha open reading frame sequences were collected from the H. polymorpha genome database (Rhein Biotech, unpublished;  and NCBI [45–47]. For the genes on contigs 1 - 48, the annotation was based on the information from RheinBiotech and Ramezani-Rad et al. . The additional NCBI sequences were listed as Hp50 and Hp51 numbers. The annotation of these genes was as described on NCBI (Hp50s) or by manual blast search (Hp51s). All Hp sequences were applied for design of oligonucleotide probes using OLIGOARRAY v2.1  with the following oligonucleotide primer design parameters: a length of 58-60 nucleotides, a melting temperature (Tm) of at least 80°C, a G/C-content of 34-52% and a maximum Tm of secondary structures and cross-hybridisations of 68°C. Oligo's were designed within the 3'-regions of the ORFs (setting: maximum distance between the 5' end of the oligo and the 3' end of the input sequence: 600-nt) to minimise including intron-sequences, since these were not discarded in the input ORF sequences. Paralogous sequences were removed during the final analysis of the design using blastN. Of the 6,248 oligo-probes, 6,002 are from the annotated genes of H. polymorpha and 23 probes are from heterologous genes of specific research interest (data not shown). The remaining 223 probes include positive and negative controls. Subsequently, the oligo-set was printed twice in each of the 8 arrays per slide (8-plex format) using Agilent's SurePrint technology (in situ synthesis; via eArray 4.0-website; Agilent Technologies Netherlands B.V., Amstelveen, the Netherlands).
H. polymorpha strain NCYC 495 leu- was grown in batch cultures on mineral medium containing either 0.5% glucose or 0.5% methanol as carbon source and 0.25% ammonium sulphate as nitrogen source [49, 50]. For transcriptome analysis, cells were extensively pre-cultivated in batch cultures on glucose at 37°C prior to a shift to fresh media containing methanol as sole carbon source. Four independent glucose-grown cultures were used to inoculate fresh medium containing methanol as well as for taking samples for RNA isolation. The methanol cultures were grown for two hours, followed by RNA isolation.
Samples were harvested by freezing them directly in liquid nitrogen, followed by thawing on ice, centrifugation (5.000 G, 4 min, 0°C) and taken up in AE-buffer (50 mM sodium-acetate 10 mM EDTA pH 5.0). 1 volume acid-phenol chloroform isoamylalcohol (125:24:1 pH 4.5, Ambion, Austin USA) and 0.5% SDS was added, and incubated at 65°C for 5 min followed by 15 min at -80°C. After centrifugation (15 min 13.000 G), the upper phase was mixed with 0.5 volume acid phenol/chloroform, centrifuged (4 min 13.000 G) and mixed with 0.5 volume chloroform. The upper phase was used for RNA isolation using column purification according to the manufacturers' instructions (Nucleospin RNA II, Macherey-Nagel, Düren Germany). The Agilent Bioanalyzer 2100 with RNA 600 Nano chips was used to analyze the quality and integrity of the RNA samples.
Transcript levels of methanol-related genes (AOX and DHAS), using those of actin as control, were determined by semi-quantitative RT-PCR, using actin as loading control (Ready-to-go RT-PCR beads, GE Healthcare, Little Chalfont UK).
Using the Low RNA input linear Amplification kit from Agilent, cDNA was generated based on 500 ng of each isolated mRNA sample. Next, cRNA was made using Cy3-CTP or Cy5-CTP incorporation for labeling purposes. For each original mRNA, a portion of Cy3 and Cy5 labeled cRNA was generated. The concentration and incorporation of the cRNA and the dyes are measured using Nanodrop. For each biological replicate, 300 ng Cy3 labeled glucose culture originating cRNA was used for hybridisation against 300 ng Cy5 labeled methanol culture originating cRNA and vice versa for the dye-swap, resulting in 8 hybridisations in total. Hybridisation, washing and scanning were performed according to the Agilent 'Two-color Microarray-based gene expression analysis protocol' (version 5.5, February 2007) by ServiceXS (Leiden, The Netherlands).
To extract the data from the scanned DNA-microarray slides, the feature extraction software version 9.5, Protocol GE2-v5_95_Feb07 from Agilent was used. For the background subtraction the option 'No background subtraction and spatial detrend' was used. For each spot, the ratio between the green and red processed signals was calculated, reflecting the ratio of gene expression on methanol overexpression on glucose. Next, the average ratio per gene was calculated based on the data of 16 spots (8 hybridisations, 2 spots per hybridisation). For reasons of clarity, genes with a ratio of <1 were expressed as -(1 divided by the ratio), thus reflecting the fold downregulation (e.g. -2 instead of 0.5). As a cut-off for differential gene expression, a threshold of (more than) 2-fold up- or downregulation was used, so >2 or <-2. To assess the significance of the data, p-values were computed using the paired-data program at the CyberT interface [53, 54]. Genes were considered to be significantly regulated if they had a p-value below 0.05. In addition, average A-values for each gene were calculated as an indication for the intensity of the signals using A = 0.5*(log2Cy3+log2Cy5). An A-value of 6 was used as a lower limit for trustworthy signal intensity.
All data has been deposited to the NCBI Gene expression omnibus and is accessible under accession number GSE19036.
To show the main represented functions among up- and downregulated genes, H. polymorpha genes were ordered according to their Functional Catalogue (FunCat) as assigned by RheinBiotech. In the diagrams, the main groups of the hierarchical structure are shown as well as the subgroups 'transcription' and 'protein synthesis'. The group 'subcellular localisation' was omitted, while 'control of cellular organisation', which is not in the FunCat structure, was included under 'developmental processes'. Genes can be present in more than one group.
Changes in expression of metabolic pathway genes were investigated using the omics viewer at the Ecocyc website . Since information on H. polymorpha is not included in this database, the genome of Saccharomyces cerevisiae S288C was used as a reference.
Ultrathin sections of KMnO4-fixed cells were used for ultrastructural analysis as described . Analysis of mitochondria was performed using confocal microscopy (Zeiss LSM510). Mitochondria were visualized using MitoTracker Orange (CMTMRos, Invitrogen) and visualization with excitation by a 543 nm Neon-laser (Lasos) and detection using a 565-615 band-pass emission filter.
TvZ is supported by a grant from the Netherlands Organisation for Scientific Research NWO (Earth and Life Sciences). This project was carried out within the research programme of the Kluyver Centre for Genomics of Industrial Fermentation, which is part of the Netherlands Genomics Initiative/Netherlands Organisation for Scientific Research.
We thank Rhein Biotech GmbH, Düsseldorf, Germany for access to the H. polymorpha genome database.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.