<title>Background</title> <p>Transcriptomics, the parallel quantification of many, or all transcripts of an organism in given conditions, has become a favorite tool for basic research [##REF##11714516##1##]. Messenger-RNA regulation patterns of model organisms under many different conditions have become available during the last years. However, these methods are still not applicable for many industrially important organisms, mainly due to the lack of DNA microarrays targeting these organisms. A typical example is the yeast <italic>Pichia pastoris</italic>, which is widely applied for the production of recombinant proteins. Several approaches have been taken to derive transcriptomic data without specific microarrays. Sauer et al. [##REF##15610561##2##] have applied heterologous hybridization of <italic>P. pastoris </italic>samples to <italic>Saccharomyces cerevisiae </italic>microarrays. Alternative methodological concepts like Transcript Analysis with the Aid of Affinity Capture (TRAC) [##REF##16198435##3##] may be applied preferentially to subsets of the transcriptome [##REF##17578563##4##], provided that genome sequence data are available. If this is not the case, total cDNA may be utilized as a source of probes, either by applying expressed sequence tags to microarrays [##REF##14732266##5##] or employing RNA fingerprinting like cDNA-amplified fragment length polymorphism (cDNA-AFLP) [##REF##8653120##6##], which has recently been applied to <italic>Trichoderma reesei </italic>[##REF##16504068##7##]. These unannotated methods bear of course the disadvantage that specific hits may only be identified after sequencing their respective probes.</p> <p>Therefore oligonucleotide microarrays have become the method of choice for many applications, although their design depends on the availability of a genomic sequence with good gene identification and annotation. The genome sequence of <italic>P. pastoris </italic>is not published yet. The data available from Integrated Genomics (IG, Chicago, IL, USA; [##UREF##0##8##]) contain a partial gene identification and annotation, so that additional effort in this direction was a first step necessary towards development of comprehensive DNA microarrays for this yeast species. There is a wide choice of computational gene finders available at the moment which can be classified into intrinsic and extrinsic prediction programs. Intrinsic or <italic>de novo </italic>gene finder only use information from the sequences to be studied, building statistical models to distinguish between coding and non-coding regions of the genome on the basis of biological sequence patterns [##REF##15980510##9##, ####REF##15145805##10##, ##REF##12364589##11####12364589##11##]. Extrinsic gene finder utilize homology search to determine where protein coding regions are in the genome. Their applicability is therefore limited to organisms that have homologs in current databases that are correctly annotated. Because of this limitation it is common to integrate homology search with <italic>de novo </italic>prediction [##REF##15691859##12##]. Most state of the art gene finders use a form of Hidden Markov Model (HMM) differing in the implementation and complexity of the model as well as the ease in which users can adapt the application to their needs [##REF##12538242##13##].</p> <p>It is well known that cross-hybridization can confound microarray results rendering good probe design an essential requirement for accurate microarray analyses. The specificity of oligonucleotides is determined by the Gibbs free energy (ΔG) of the hybridization reaction between potential binding partners. Highly specific probes will bind their target transcript much more strongly than any other transcript. Considering that microarray experiments are non-equilibrium measurements, it is desirable that microarray probes exhibit uniform thermodynamic properties, which many probe design tools aim to achieve by demanding a narrow distribution of the probe-target melting temperature T<sub>m</sub>. Ideally, probes should have a uniform binding free energy at the hybridization temperature <italic>T</italic><sub><italic>hyb </italic></sub>[##REF##16938547##14##].</p> <p>Previous studies have demonstrated that industrial production strains may behave quite differently to laboratory strains and model organisms [##REF##12748050##15##], which emphasizes the importance of analytical tools for industrially relevant strains and species. As an example, the unfolded protein response (UPR), a regulation circuit of high relevance for heterologous protein production in eukaryotic cells [##REF##12676684##16##], has been shown to be differentially regulated in <italic>P. pastoris </italic>[##REF##17578563##4##] compared to <italic>S. cerevisiae </italic>[##REF##10847680##17##], which is the typical model species for hemiascomycete yeasts. The development of specific microarrays for <italic>P. pastoris </italic>was intended to allow a detailed analysis of UPR regulation in <italic>P. pastoris</italic>. As in previous transcriptomics work with <italic>S. cerevisiae </italic>the induction of UPR was either accomplished by addition of dithiothreitol (DTT) or tunicamycin, this work aimed at a comparison of DTT induced gene regulation in <italic>P. pastoris </italic>to that in <italic>S. cerevisiae </italic>published by Travers et al. [##REF##10847680##17##]. Finally we aimed at the comparison of DTT induced regulation to the regulatory response to overexpression of <italic>HAC1</italic>, the transcription factor controlling the UPR. Transcriptional regulation of <italic>HAC1 </italic>overexpression has not been studied for yeasts so far, so that we expected valuable data to better define the core UPR regulated transcriptome.</p>
<title>Methods</title> <title>Gene Prediction and Sequence Selection</title> <p>Gene prediction and the selection of sequences for oligonucleotide probes were based on sequenced contigs of the <italic>P. pastoris </italic>genome including predictions of protein coding genes, available through Integrated Genomics [##UREF##0##8##]. The number of predicted genes was 5,425 of which 3,680 had an assigned function. The ORFs were made up of experimentally identified genes, as well as ORFs predicted by a proprietary gene finder [##REF##12519973##38##].</p> <p>To validate and possibly improve these predictions, <italic>de novo </italic>gene finding was conducted. First three <italic>de novo </italic>gene finder (GeneMark, Glimmer3, GlimmerHMM) were tested on the genome sequence of <italic>S. cerevisiae </italic>(data from BioMart, [##UREF##5##39##]) to evaluate their performance on yeast genomes. As described in Results and Discussion, GeneMark [##UREF##6##40##] was selected for further gene prediction on the <italic>P. pastoris </italic>genome sequence. To run the gene prediction it was necessary to train GeneMark on <italic>S. cerevisiae </italic>by building a matrix with transition probabilities for coding and non-coding regions used by the Hidden Markov Model (HMM) of the program. With the amount of data available we were able to generate a matrix of the 7<sup>th </sup>order. The genes of <italic>P. pastoris </italic>were predicted using the <italic>S. cerevisiae </italic>matrix and the lowest possible probability score cut-off (t = 0.05). In the initial stage of the microarray design the aim was to predict as many putative ORFs as possible. In this context a higher false positives rate was accepted in order to keep the false negatives rate as low as possible.</p> <p>The predicted sequences were merged with data from IG and clustered by running cd-hit [##REF##16731699##20##] with a similarity cut-off of 90%. For all of the resulting sequences a BLASTX search was done against <italic>S. cerevisiae </italic>using WU-BLAST [##UREF##1##19##]. Blast data was further filtered for length (cutoff 55 bp) and low prediction score. Clusters comprised of more than one gene were represented by the longest sequence, or curated manually, if appropriate.</p> <p>From this first gene list (PpaV1) microarrays were analyzed as described below. Spots with a positive signal were determined using the mean plus one standard deviation of the negative control probes as a cut-off. Sequences were selected if they were positive in at least 8 out of 12 arrays. This criterion was chosen to fill the array capacity. Additionally all sequences with a probability score higher than 0.5 or having an annotation were kept for the second set of sequences (PpaV2).</p> <title>Annotation</title> <p>For the PpaV2 sequence set the program cd-hit-est [##REF##16731699##20##] was used to find all ORFs that had a global identity of &gt; 80% with <italic>S. cerevisiae</italic>. WU-BLASTX and WU-TBLASTN searches were conducted against <italic>S. cerevisiae</italic>, using a low complexity filter and E &lt; 10<sup>-7</sup>. For all the sequences that did not have a match with <italic>S. cerevisiae </italic>under these conditions the two BLAST searches were repeated against the SwissProt/TrEMBL [##UREF##7##41##] database. A perl script was developed to summarize and compare the BLAST results.</p> <title>Oligo Design and Array platform</title> <p>Oligos for the PpaV1 sequences were designed with the Program OligoArray 2.1 [##REF##12799432##21##] to match the melting temperature distribution of Agilent's <italic>S. cerevisiae </italic>oligos on the Yeast Oligo Microarray (V2), design number 013384.</p> <p>The oligo-set for the PpaV2 sequence set was designed using the thermodynamic model-based oligoset optimizer 'TherMODO'. This tool incorporates advanced quantitative models for probe-target binding region accessibility and position-dependent target labelling efficiency, and replaces the common greedy search algorithm by a global set optimization step, achieving high discrimination power for particularly uniform probe sets [##UREF##2##23##]. Probes for Agilent arrays are limited to a maximum length of 60 nucleotides by the manufacturing process. For increased flexibility in the probe design, the oligoset design optimization considered probes ranging in length from 57 to 60 nucleotides.</p> <p>These arrays were produced on Agilent 60 mer oligonucleotide high density arrays 4 × 44 K (with 42,034 available features) for PpaV1 and 8 × 15 K (with 15,208 available features) for PpaV2.</p> <title>Experimental Design</title> <p>For the first batch of arrays a same-same design was used, employing six replicates each of Pool 1 and of Pool 2. The aim of this experiment was to determine which of the probes hybridize to <italic>P. pastoris </italic>targets. For the second batch of arrays a two-state comparison set up was chosen with 6 replicates for each experiment of which 3 were dye swapped.</p> <title>Strains und Cultures</title> <p>For the first batch of arrays the aim was to determine which of the predicted probes hybridize with targets from <italic>P. pastoris</italic>. To make sure that many genes were active it was important to pool samples from various conditions of the cells. Samples were taken from two different <italic>P. pastoris </italic>strains, X-33 and CBS2612, grown on different media and taken at both exponential and stationary growth phase. The media were YP Medium (1% yeast extract, 2% peptone and either 2% glucose, 2% glycerol or 0.5% methanol as carbon source), Buffered Minimal Medium (1.34% yeast nitrogen base, 4 × 10<sup>-5</sup>% biotin, 100 mM potassium phosphate pH 6.0 and either 2% glucose, 2% glycerol or 0.5% methanol as carbon source), and Buffered Minimal Medium described above supplemented with amino acids (0.005% of L-glutamic acid, L-methionine, L-lysine, L-leucine and L-isoleucine). The samples were combined into two pools with Pool 1 containing 18 samples from the exponential growth phase and Pool 2 containing 18 samples from the stationary phase. Both pools additionally contained seven chemostat samples of the strain X-33 3H6Fab, grown as in [##REF##18078287##42##].</p> <p>For the UPR experiments, strains GS115 HAC1, constitutively overproducing the activated form of <italic>S. cerevisiae </italic>Hac1, as described in Gasser et al. [##REF##16570317##33##,##REF##17578563##4##], as well as GS115 transformed with the empty vector pGAPHIS (a histidine prototrophic isogenic strain of GS115) were cultivated in YPD (YP as above with glucose) at 28°C. After growing the cultures to an OD<sub>600 </sub>= 5.7, dithiothreitol (2.5 mM) was added where appropriate. After 1 more hour of cultivation, 1 ml culture was added to 0.5 ml precooled phenol solution (5% in absolute ethanol) and centrifuged immediately for 30 sec at 13.000 rpm. After discarding the supernatants the pellets were frozen at -80°C.</p> <title>RNA Isolation</title> <p>All samples were resuspended with 1 mL TRI Reagent (Sigma). Cells were disrupted after addition of 500 μL glass-beads with a Thermo Savant Fastprep FP120 Ribolyzer by treatments of 2 × 20 sec at 6.5 ms<sup>-1</sup>. RNA was extracted with chloroform, precipitated with isopropanol, washed with 75% ethanol and dissolved with diethylpyrocarbonate treated water. The extracted RNAs were quantified via absorption at 260 and 280 nm. The quality of the RNA samples was verified with the Agilent Bioanalyzer 2100 and RNA 6000 Nano Assay kit (Agilent Technologies, California).</p> <title>Labeling and Hybridization</title> <p>Hybridization targets for <italic>P. pastoris </italic>microarrays were prepared according to Agilent's Two-Color Microarray-Based Gene Expression Analysis protocol (Version 5.5, February 2007). Purification of the labelled and amplified RNA was conducted using RNeasy mini spin columns (Qiagen). The quality of labelled cRNA was evaluated on the Agilent Bioanalyzer 2100 and quantified using a ND-1000 NanoDrop Spectrophotometer. Fragmented cRNA samples were applied to the individual arrays. The slides were placed into Agilent hybridization oven and hybridized for 17 h, at 65°C and 10 rpm.</p> <title>Microarray Analysis</title> <p>Slides were scanned with an Agilent MicroArray Scanner and intensities were extracted using Agilent's Feature Extraction software (version 9.1). The resulting data was imported into R where data pre-processing and normalization was performed. In the pre-processing step all outliers and saturated spots were given the weight zero. After plotting the data we decided to refrain from background correction since it has the tendency to add more noise to the data [##REF##17472750##43##]. The data were normalized using locally weighted MA-scatterplot smoothing (LOESS) followed by a between array scale normalization. Both functions are available within the limma package of R [##REF##14597310##44##]. For the selection of differentially expressed genes linear models were fitted to the log-ratios of the expression data separately for each gene. An empirical Bayes approach was used to shrink the probe-wise sample variances towards a common value yielding a moderated <italic>t</italic>-statistic per gene [##REF##16646809##45##]. <italic>P</italic>-values were corrected for multiple testing using Holm's method [##UREF##8##46##]. Features were defined as differentially expressed if they had a <italic>p</italic>-value &lt; 0.05. For the identification of stronger regulatory effects an additional cut-off for the fold change (FC) of 1.5 &gt; FC &gt; 1/1.5 was applied. Description of the platform, array, raw data as well as processed data were deposited at ArrayExpress [##UREF##9##47##] under the accession numbers A-MEXP-1157.</p> <p>All annotated <italic>P. pastoris </italic>genes were categorized into GO biological process terms using the SGD GO slim tool [##UREF##4##28##], whereby <italic>P. pastoris </italic>specific genes were included into the term 'other'. The significance of a deviation of the number of up- or downregulated genes in each group from the average was verified with a Fisher test (Additional file ##SUPPL##3##4##).</p>
<title>Results and Discussion</title> <title>Gene prediction and Oligo Design</title> <p>To evaluate available gene finders for their performance on yeast genomes, three <italic>de novo </italic>gene finders (GeneMark, Glimmer3, GlimmerHMM) were tested on the genome sequence of <italic>S. cerevisiae</italic>. GeneMark and Glimmer3 work with a prokaryotic Hidden Markov Model (HMM) whereas GlimmerHMM employs a eukaryotic gene model. GeneMark was trained with coding and non-coding sequences of <italic>S. cerevisiae</italic>, building an HMM transition probability matrix of the 7<sup>th </sup>order. Glimmer3 and GlimmerHMM could be trained directly on the genome in question without specifying coding and non-coding regions. In Lomsadze et al. [##REF##16314312##18##] and Besemer and Borodovsky [##REF##15980510##9##] the difficulty of eukaryotic gene finders in the prediction of genes for organisms with few introns is discussed and linked to a lack of data for representative exon – intron models. Our results confirmed that a gene finder written for eukaryotes (GlimmerHMM) could not be trained well on yeast genomes, introducing far too many introns into the predicted genes. Both prokaryotic versions performed much better, with GeneMark predicting less false negatives but more false positives than Glimmer3 (Table ##TAB##0##1##). Even though the positive prediction value was somewhat lower with GeneMark it was more important not to miss true positives than to achieve a lower rate of false positives. A further improvement could be achieved by a GeneMark model for lower eukaryotes, in which the prokaryotic algorithm is modified to use Kozak start sites instead of prokaryotic ribosomal binding sites. <italic>P. pastoris </italic>genes were predicted using this version of GeneMark with the lowest possible threshold (probability score t = 0.05) so that filter conditions could be better controlled at a later state. The prediction yielded a total of 26,471 putative genes for the genome of <italic>P. pastoris</italic>.</p> <p>In a WU-BLASTN search against <italic>S. cerevisiae</italic>, 6,374 sequences that were predicted by GeneMark, and 3,964 of the IG predictions produced hits with <italic>S. cerevisiae </italic>using an <italic>E </italic>value (Expectation value, [##UREF##1##19##]) of &lt; 10<sup>-4</sup>, a hit length &gt; 100 nucleotides and an identity of &gt;50%. To reduce the redundancy within the data set the predicted genes were clustered into groups sharing more than 90% similarity using cd-hit [##REF##16731699##20##]. From a total of 31,896 candidate sequences (GeneMark and IG predictions), 22,020 cd-hit groups were obtained. From the cluster file it was clear that some of the clusters had to be analyzed further before selecting target sequences for the oligo design. After the removal of all sequences that had a short length and a low prediction value, complex clusters were defined as clusters for which the minimum relative length of all sequences was smaller than 0.9. A total of 2,612 clusters fell into this category and were excluded at a first design stage.</p> <p>Finally 19,508 predicted target sequences remained to be tested in the first microarray experiments. OligoArray 2.1 [##REF##12799432##21##] was able to design oligonucleotide probes for 17,161 sequences ranging in length from 57 to 60 nucleotides.</p> <title>Validation arrays for the first list of predicted transcript sequences (Same-Same experiment)</title> <p>With these probes 4 × 44 K slides were produced on the Agilent microarray platform and employed for an initial validation of the predicted transcript sequences by hybridization with the Pool samples of <italic>P. pastoris </italic>(for preparation of Pool samples see Material and Methods). One slide had to be discarded because of quality issues. For the remaining 12 arrays the number of probes showing a signal varied between 10,708 and 15,598. Of these, 7,980 had a signal on all 12 arrays, and only 951 probes showed no hybridization on all 12 arrays.</p> <title>Second, curated list of predicted target sequences and second oligo design</title> <p>The results of the initial validation arrays were utilized to adapt the list of predicted genes, keeping all predictions for which a hybridization signal could be observed for all arrays plus all predictions with significant sequence similarity to annotated genes as well as all sequences with an average gene prediction score &gt; 0.5. This approach allows for the fact that not all genes will have been actively expressed in the target samples. Additionally, predicted transcripts resulting from a subsequent analysis of the complex clusters were included at this stage. Of the 2,612 complex cluster that were not included in the design for the first batch of arrays, only 223 contained more than 2 sequences and for a further 14 no subsequence match of at least 60 nucleotides could be found within the last 1000 bases at the 3'-end. These 237 clusters were manually curated while the rest could be automatically reduced to one sequence. To make full use of the 15,208 features available on the Agilent microarray platform, it was decided to also include predicted sequences with somewhat lower gene prediction score that showed a hybridization signal in at least 8 of the 12 arrays. Finally, a selected set of 15,253 predicted transcript sequences was used as targets for probe design of a comprehensive <italic>P. pastoris </italic>microarray. While it is obvious that this list is larger than the expected number of open reading frames (6,000–7,000), as judged in comparison to other yeast species [##REF##17868481##22##], we intentionally included more putative transcript sequences, as false positives with a distinct sequence will not negatively affect microarray design or experiments, in contrast to the damage of falsely excluding a potential transcript target.</p> <p>Oligonucleotide probes were designed using a probe design tool developed in-house, a thermodynamic model-based oligoset optimizer ('TherMODO', [##UREF##2##23##]). TherMODO designed probes for 15,035 sequences, of which only 665 were predicted as having cross-hybridization potential. The TherMODO design was compared to probe design with eArray [##UREF##3##24##]. The distributions of ΔG and T<sub>m </sub>of both designs are shown in additional file ##SUPPL##0##1##. Clearly the TherMODO designed probes are more uniform in respect to the Gibbs free energy ΔG, indicating a superior hybridization performance [##REF##16938547##14##].</p> <p>The final probe design was manufactured on 8 × 15 K slides by Agilent, and evaluated for reproducibility and biological meaningfulness. Pool samples were applied to 2 arrays on 2 slides each, including dye swap. The scatterplots show uniformly high correlations &gt; 97% both within and between arrays, both on same and different slides, indicating high reproducibility of hybridization signals between identical samples. Exemplarily, a scatterplot of signal intensities derived from the same samples (wild type strain untreated) is shown in Figure ##FIG##0##1##. For the final gene list the annotation was improved in addition to the annotation provided by IG. This resulted in 3954 annotated ORFs, of which 2989 had an IG annotation. 965 newly annotated ORFs were found, and the annotation of 288 hypothetical proteins was confirmed. All annotated genes are listed in Additional file ##SUPPL##1##2##.</p> <title>Biological evaluation of the new microarrays</title> <p>The performance of the new arrays was examined by a hybridization experiment using samples, for which transcript regulation data have been obtained before [##REF##17578563##4##]. The biological question evaluated was the regulatory response of <italic>P. pastoris </italic>to constitutive overexpression of the active form of <italic>S. cerevisiae HAC1</italic>, the transcription factor controlling UPR target genes. By this approach, the regulation of 52 genes which have been studied before using TRAC [##REF##16198435##3##] could be verified, with 80% of these genes showing the same regulation pattern for both methods (genes highlighted in bold in Additional file ##SUPPL##1##2##). This correlation is statistically significant based on calculating the regression (p = 8.8 · 10<sup>-6</sup>).</p> <p>The similarities and differences of UPR induction and reaction to DTT stress have been discussed before [##REF##17578563##4##,##REF##17561995##25##,##REF##15870366##26##]. To achieve further insight into this technologically relevant issue, we compared the gene regulation patterns of a <italic>HAC1 </italic>overexpressing strain <italic>vs </italic>wildtype control with the regulation pattern of the wildtype treated with DTT for 60 min <italic>vs </italic>the untreated control. Genes were qualified as significantly regulated with a <italic>p</italic>-value &lt; 0.05 (adjusted for multiple testing). 11,262 of all features on the microarrays appeared as differentially regulated either upon DTT treatment or <italic>HAC1 </italic>overexpression, or both. 8,480 reacted to <italic>HAC1</italic>, and 6,870 to DTT, with an overlap of 4,088. Considering only the 3,954 annotated genes, a similar pattern is observed with roughly half of the regulated genes overlapping between DTT and <italic>HAC1</italic>, and another half being typical only for either of the treatments (Figure ##FIG##1##2##). Accordingly, the correlation of log fold changes of the two treatments is apparent but rather weak (Figure ##FIG##2##3##). While DTT treatment is widely accepted as a standard inducer of UPR, these observations indicate that the gene regulation pattern triggered by the UPR transcription factor Hac1 differs to a significant extent from that exerted by DTT.</p> <p>As previous research on transcriptome regulation upon UPR induction usually employs a fold change (FC) cut-off to highlight the strongly regulated genes, we decided to introduce FC &gt; 1.5 as a second criterion to identify more strongly regulated genes for further detailed analysis (Volcano plots visualizing the two criteria are provided in Additional file ##SUPPL##2##3##). Although the introduction of a FC cut-off alters the absolute number of regulated genes, it does not alter the relative distribution of regulated genes categorized into functional groups (GO slim biological process), as can be seen in Figure ##FIG##3##4## and Additional file ##SUPPL##3##4##.</p> <title>Comparison of UPR induction by DTT in <italic>P. pastoris </italic>and <italic>S. cerevisiae</italic></title> <p>In order to compare the effects of DTT treatment in <italic>S. cerevisiae </italic>with those in <italic>P. pastoris</italic>, the data published by Travers et al. [##REF##10847680##17##] for 60 min treatment of <italic>S. cerevisiae </italic>with DTT were evaluated alongside with our results for <italic>P. pastoris</italic>. All genes of <italic>S. cerevisiae </italic>which were listed in [##REF##10847680##17##] and for which homologs in <italic>P. pastoris </italic>were identified were classified as upregulated, downregulated or unregulated. In order to compare the two data sets, a cutoff of 1.5 fold differential expression was set in both to define regulated genes. A significance threshold on <italic>p</italic>-values could not be employed, as these data were not provided for <italic>S. cerevisiae</italic>. 48% of these genes defined as regulated or unregulated reacted in <italic>P. pastoris </italic>just as in <italic>S. cerevisiae</italic>.</p> <p>A closer evaluation revealed that certain GO groups were regulated very similarly in both yeast species, while others showed only a low degree of similarity (Table ##TAB##1##2##). Fisher's exact test was performed to evaluate the significance of groups with low similarity. Especially the GO groups 'translocation', 'protein folding', 'protein degradation', and to some extent 'glycosylation' and 'transport' showed high degrees of similarity. In some GO groups, only some subgroups reacted similarly while others behaved differently in the two yeasts. Of the 'glycosylation' group, core oligosaccharide synthesis and glycosyltransferase genes behaved very similarly, while glycoprotein processing, GPI anchoring and O-glycosylation related genes were regulated significantly different (p &lt; 0.05). In the 'protein degradation' group, more similarity was observed for ERAD genes than for ubiquitin/proteasome related genes. Among the 'transport' gene group, budding, fusion and retrieval of ER to Golgi showed a high degree of similar regulation, contrary to the subgroup distal secretion. Low similarities were observed for 'lipid metabolism', 'vacuolar protein sorting' and 'cell wall biogenesis' genes. It becomes obvious that core UPR genes related to protein translocation, folding and ER transport, as well as core N-glycosylation react similarly to DTT treatment in <italic>P. pastoris </italic>as compared to <italic>S. cerevisiae</italic>, while genes involved in processes which are more distal from ER protein folding behave more differently, indicating that those processes (like functions in the Golgi, [##REF##10189369##27##]) differ significantly between the two yeasts.</p> <title>Overexpression of Hac1 triggers a different regulation pattern compared to DTT treatment</title> <p>In most previous studies of the UPR in lower eukaryotic cells, treatment with DTT or tunicamycin, or heterologous protein expression has been employed to trigger the UPR. This study clearly indicates that the set of regulatory events triggered by DTT analysis only partially overlaps with the reactions to constitutive expression of the activated form of the UPR transcription factor Hac1 (see Figures ##FIG##1##2## and ##FIG##2##3##). Interestingly, both treatments resulted in the same amount of genes being down-regulated as being up-regulated, a fact that has been neglected to some extent in the existing literature.</p> <p>Those genes appearing beyond the threshold (<italic>p</italic>-value &lt; 0.05 and FC &gt;1.5) were subjected to a more detailed comparison between the effects of DTT treatment and Hac1 induced regulation. The relative numbers of up- and downregulated genes in each GO biological process term based on the SGD GO slim tool [##UREF##4##28##] are depicted in Figure ##FIG##3##4##.</p> <p>A pattern common to both treatments is the down-regulation of major metabolic processes like carbohydrate, amino acid and lipid metabolism, as well as that of vitamins, cofactors and aromatic and heterocyclic compounds. This makes it obvious that the UPR has a major impact on decreasing both catabolic and anabolic processes. On the other side, both treatments lead to up-regulation of protein folding and vesicular transport. These effects are in line with the published literature, indicating the cellular reaction towards alleviation of the UPR [##REF##17578563##4##,##REF##17561995##25##,##REF##15870366##26##,##REF##10847680##17##].</p> <p>As expected, the genes coding for classical UPR targets are induced both in Hac1 overproducing and in DTT stressed cells, and genes underlined in the following paragraphs have been identified as UPR targets in previous studies. Especially the ER folding catalysts <italic><underline>PDI1</underline></italic> and <italic><underline>ERO1</underline></italic>, the DnaJ homologs <italic><underline>JEM1</underline></italic> and <italic><underline>SCJ1</underline></italic>, the ER resident chaperones <italic>CNE1 </italic>(calnexin), <underline><italic>KAR2</italic>/BiP</underline> and <italic><underline>LHS1</underline></italic> and the mitochondrial chaperones <italic>HSP60 </italic>and <italic>SSC1 </italic>are significantly up-regulated in both conditions. Among the functional group of 'protein modification' the majority of up-regulated genes belong to the core oligosaccharide synthesis (<italic><underline>DPM1</underline></italic>, <italic>DIE2</italic>), oligosaccharyltransferase complex (<italic>OST1</italic>, <italic><underline>OST2</underline></italic>, <italic><underline>OST3</underline></italic>, <italic><underline>SWP1</underline></italic>, <italic>STT3</italic>, <italic><underline>WBP1</underline></italic>), glycoprotein processing (<italic>ALG2, <underline>ALG7</underline></italic>, <italic>SEC53</italic>), GPI anchor biosynthesis (<italic><underline>GPI2</underline></italic>, <italic>GPI14, PSA1</italic>) and Golgi/O-linked glycosylation (<italic><underline>PMT1</underline></italic>, <italic><underline>PMT2</underline></italic>, <italic>PMT4, PMT6</italic>). Besides these, several genes coding for the translocon pore complex (<italic><underline>SEC61</underline></italic>, <italic><underline>SEC62</underline></italic>, <italic>SEC63</italic>, <italic><underline>SEC72</underline></italic>, <italic>SSS1</italic>), which aid the translocation of nascent polypeptides into the ER, are induced. Higashio and Kohno [##REF##12176018##29##] describe the stimulation of ER-to-Golgi transport through the UPR by inducing COPII vesicle formation. In this context, we see <italic>SEC23</italic>, <italic><underline>SEC24</underline></italic>, <italic><underline>SFB2</underline></italic>, <italic>YIP3</italic>, and <italic>ERV2 </italic>upregulated. However, also proteins building the COPI coatomer, which are required for retrograde Golgi-to-ER transport, show increased transcription levels upon ER stress in our experiments (<italic>COP1</italic>, <italic><underline>RET2</underline></italic>, <italic>SEC21</italic>, <italic><underline>SEC27</underline></italic>).</p> <p>While we cannot give any information on ERAD regulation, as <italic><underline>HRD1</underline></italic> is the only annotated gene of this protein degradation process (up-regulated in the Hac1 strain), we observed the down-regulation of some components involved in the assembly of the 20 S core of the 26 S proteasome (<italic>ADD66</italic>, <italic>PRE1</italic>, <italic>PRE4</italic>, <italic>SCL1</italic>) and ubiquitin <italic>UBI4 </italic>upon constitutive UPR activation. In this context, Shaffer et al. [##REF##15345222##30##] describe reduced degradation of newly synthesized proteins in XBP1-overexpressing human Raji cells.</p> <p>Induction of genes encoding cytosolic chaperones (Cns1, Jjj3, Hsp82, Ssa1, Ssa2, Sse1, Ydj1, Zuo1) can only be seen in the Hac1-overproducing strain. Additionally, the ER-resident Pdi homolog <underline>Mpd1</underline> and two members of the PPIases (<italic>FPR4 </italic>and <italic>CPR6</italic>) are only up-regulated in the engineered strain, but not upon DTT addition.</p> <p>One of the most striking patterns is the significant up-regulation of a large number of genes with functions in ribosomal biogenesis (233 genes assigned to the GO-categories 'ribosome biogenesis and assembly' and 'RNA metabolic process'). Most of these genes are contributing to rRNA processing (RRP family) and ribosome subunit nuclear export and assembly, while the ribosomal proteins (RPS and RPL families) themselves are not among the regulated genes for <italic>P. pastoris </italic>(see Additional file ##SUPPL##1##2##). No genes with a function in mRNA decay show increased transcription levels. The induction of the above functional categories came as a surprise, as translational down-regulation of proteins involved in ribosomal biogenesis was recently reported when <italic>S. cerevisiae </italic>cell were treated with DTT [##REF##18206654##31##]. In contrast, the transcription levels of 9 out of the 16 mRNAs listed by these authors are enhanced in our study. Transcriptional down-regulation of ribosomal proteins during ER stress conditions was also revealed when reanalysing the raw data provided by Travers et al. [##REF##10847680##17##]. However, Shaffer et al. [##REF##15345222##30##] describe an increase in total protein synthesis as well as in the number of assembled ribosomes upon the overexpression of the mammalian Hac1 homolog XBP1 in Raji cells, but did not observe upregulation of genes related to ribosome biogenesis. A similar effect was observed after XBP1 overexpression in CHO-K1 cells [##REF##16635796##32##]. These results may be an indication that the positive effect of overexpression of the UPR transcription factor on heterologous protein production [##REF##16570317##33##,##REF##12676684##16##,##REF##14660339##35##] results not just from stimulation of folding and secretion of proteins but also their synthesis. The induction of protein folding related genes upon Hac1 overexpression is in line with the literature on UPR effects, while an impact on organelle biosynthesis other than ER and Golgi has so far only been described for mammalian cells.</p> <p>The stimulatory effects of XBP1 induction on ribosomes and organelle synthesis in mammalian cells like lymphocytes have been attributed to their function as dedicated protein factories. On the other hand the UPR in lower eukaryotes should rather serve to alleviate the load of unfolded, aggregation prone protein. It will be of interest in the future to investigate whether Hac1 stimulates ribosome biogenesis in other yeasts and fungi as well, and whether this leads to increased translation.</p> <p>In this context, it is worthwhile to mention the induction of two pathways leading to the unusual post-translationally modified amino acid derivatives diphthamide and hypusine which are exclusively found in eukaryotic translation elongation factors 2 (eEF2) and 5 (eEF5), respectively [##REF##8406038##36##,##REF##17476569##37##]. As these biosynthetic pathways are rather complex, and outstanding in the otherwise downregulated group of 'amino acid biosynthesis', this induction underlines the increased demand for protein synthesis.</p> <p>Furthermore, we observe that ER stress leads to increased transcription of genes coding for the large and small subunits of the mitochondrial ribosomes (<italic>MRPS</italic>, <italic>RSM </italic>and <italic>MRPL </italic>families), mitochondrial translation initiation and elongation factors (<italic>IFM1</italic>, <italic>MEF1</italic>, <italic>MEF2</italic>) and mitochondrial DNA polymerase (<italic>MIP1</italic>). Several essential constituents of the mitochondrial inner membrane presequence translocase (<italic>TIM </italic>family) are also up-regulated, indicating increased necessity for protein import into the mitochondria. Similarly, XBP1 was shown to increase mitochondrial mass and function in two types of mammalian cells [##REF##15345222##30##].</p> <p>While previous studies analysing UPR regulation mainly focus on up-regulated genes [##REF##10847680##17##], more than half of the genes identified in our study to be regulated are strongly down-regulated (at least 1.5 fold). As can be seen in Figure ##FIG##3##4##, anabolic processes such as vitamin production, amino acid and aromatic compound biosynthesis, heterocycle metabolic processes, carbohydrate, lipid and cofactor metabolism are among the most prominent repressed classes in both DTT-treated as well as Hac1-overproducing cells. The down-regulation of energy consuming biosynthetic pathways emerges as a general picture during ER stress conditions. However, it becomes obvious that the response to the folding perturbation agent DTT strongly differs from constitutive UPR induction by Hac1-overproduction. Especially the prominent down-regulation of genes belonging to 'electron transport' and 'cellular respiration' can easily be explained by the strong reducing capacities of DTT. Prominent members of the mitochondrial inner membrane electron transport chain such as subunits of the cytochrome c oxidase (<italic>COX4, COX4, COX5A, COX13</italic>) and the ubiquinol cytochrome-c reductase complex (<italic>COR1, QRC6, QRC7, QRC9, RIP1</italic>) are significantly repressed upon DTT treatment. Additionally, cytochrome c (<italic>CYC1</italic>), cytochrome c1 (<italic>CYT1</italic>) and cytochrome c heme lyase (<italic>CYC3</italic>) are only under DTT-dependent repression (GO: 'generation of precursor metabolites and energy'). The reducing features of DTT are most probably also the reason for the up-regulation of genes involved in the upkeeping of 'cellular homeostasis' and clearly, addition of DTT is provoking a 'response to a chemical stimulus'.</p> <p>Down-regulated genes appearing in both Hac1 and DTT in the 'protein modification' group focus on protein kinases (<italic>CDC5, CDH1, DBF2</italic>) and components of the ubiquitinylation complex (<italic>BUL1, CUL3</italic>) involved in cell cycle regulation driving the cells towards mitotic exit (<italic>CDC5, CDH1</italic>, <italic>MOB1</italic>). These effects are even more pronounced in the Hac1-strain, where several more histone modifying enzymes as well as cycline-dependent protein kinases and components of the protein kinase C signalling pathway show reduced transcription levels compared to the wild type. Unlike reported for the filamentous fungi <italic>T. reesei </italic>[##REF##16504068##7##] and <italic>A. nidulans </italic>[##REF##15870366##26##], genes encoding the histones H2A, H2B, H3 and H4 appear to be down-regulated upon secretion stress in <italic>P. pastoris</italic>.</p> <p>No clear picture emerges regarding the regulation of 'lipid metabolism': While sterol and ergosterol biosynthesis tend to be inhibited, the production of sphingolipid precursor substances is enhanced. On the other hand, a down-regulation of the major cell wall constituents (β-1,3 glucanases <italic>BGL2 </italic>and <italic>EXG1</italic>, cell wall mannoproteins <italic>CCW12</italic>, <italic>CWP2 </italic>and <italic>TPI1</italic>, GPI-glycoproteins <italic>GAS1 </italic>and <italic>SED1</italic>, <italic>PST1</italic>) and genes coding for proteins required for the transport of cell wall components to the cell surface (<italic>SBE22</italic>) is manifest. Taken together, these results indicate a significant remodelling process regarding the <italic>P. pastoris </italic>cell envelope during ER stress conditions.</p> <p>Interestingly, the major groups of metabolic genes were down-regulated upon Hac1 overexpression, indicating a decrease of the supply of metabolites. However, it should be noted that no reduction of the specific growth rate was observed as compared to the wild type strain (μ = 0.37 and 0.39 h<sup>-1</sup>, respectively). A reduction of metabolic processes, and amino acid synthesis in particular, is contradictory to translation stimulation. Further research will be needed to elucidate the overall regulatory pattern of UPR in respect to protein synthesis.</p>
<title>Results and Discussion</title> <title>Gene prediction and Oligo Design</title> <p>To evaluate available gene finders for their performance on yeast genomes, three <italic>de novo </italic>gene finders (GeneMark, Glimmer3, GlimmerHMM) were tested on the genome sequence of <italic>S. cerevisiae</italic>. GeneMark and Glimmer3 work with a prokaryotic Hidden Markov Model (HMM) whereas GlimmerHMM employs a eukaryotic gene model. GeneMark was trained with coding and non-coding sequences of <italic>S. cerevisiae</italic>, building an HMM transition probability matrix of the 7<sup>th </sup>order. Glimmer3 and GlimmerHMM could be trained directly on the genome in question without specifying coding and non-coding regions. In Lomsadze et al. [##REF##16314312##18##] and Besemer and Borodovsky [##REF##15980510##9##] the difficulty of eukaryotic gene finders in the prediction of genes for organisms with few introns is discussed and linked to a lack of data for representative exon – intron models. Our results confirmed that a gene finder written for eukaryotes (GlimmerHMM) could not be trained well on yeast genomes, introducing far too many introns into the predicted genes. Both prokaryotic versions performed much better, with GeneMark predicting less false negatives but more false positives than Glimmer3 (Table ##TAB##0##1##). Even though the positive prediction value was somewhat lower with GeneMark it was more important not to miss true positives than to achieve a lower rate of false positives. A further improvement could be achieved by a GeneMark model for lower eukaryotes, in which the prokaryotic algorithm is modified to use Kozak start sites instead of prokaryotic ribosomal binding sites. <italic>P. pastoris </italic>genes were predicted using this version of GeneMark with the lowest possible threshold (probability score t = 0.05) so that filter conditions could be better controlled at a later state. The prediction yielded a total of 26,471 putative genes for the genome of <italic>P. pastoris</italic>.</p> <p>In a WU-BLASTN search against <italic>S. cerevisiae</italic>, 6,374 sequences that were predicted by GeneMark, and 3,964 of the IG predictions produced hits with <italic>S. cerevisiae </italic>using an <italic>E </italic>value (Expectation value, [##UREF##1##19##]) of &lt; 10<sup>-4</sup>, a hit length &gt; 100 nucleotides and an identity of &gt;50%. To reduce the redundancy within the data set the predicted genes were clustered into groups sharing more than 90% similarity using cd-hit [##REF##16731699##20##]. From a total of 31,896 candidate sequences (GeneMark and IG predictions), 22,020 cd-hit groups were obtained. From the cluster file it was clear that some of the clusters had to be analyzed further before selecting target sequences for the oligo design. After the removal of all sequences that had a short length and a low prediction value, complex clusters were defined as clusters for which the minimum relative length of all sequences was smaller than 0.9. A total of 2,612 clusters fell into this category and were excluded at a first design stage.</p> <p>Finally 19,508 predicted target sequences remained to be tested in the first microarray experiments. OligoArray 2.1 [##REF##12799432##21##] was able to design oligonucleotide probes for 17,161 sequences ranging in length from 57 to 60 nucleotides.</p> <title>Validation arrays for the first list of predicted transcript sequences (Same-Same experiment)</title> <p>With these probes 4 × 44 K slides were produced on the Agilent microarray platform and employed for an initial validation of the predicted transcript sequences by hybridization with the Pool samples of <italic>P. pastoris </italic>(for preparation of Pool samples see Material and Methods). One slide had to be discarded because of quality issues. For the remaining 12 arrays the number of probes showing a signal varied between 10,708 and 15,598. Of these, 7,980 had a signal on all 12 arrays, and only 951 probes showed no hybridization on all 12 arrays.</p> <title>Second, curated list of predicted target sequences and second oligo design</title> <p>The results of the initial validation arrays were utilized to adapt the list of predicted genes, keeping all predictions for which a hybridization signal could be observed for all arrays plus all predictions with significant sequence similarity to annotated genes as well as all sequences with an average gene prediction score &gt; 0.5. This approach allows for the fact that not all genes will have been actively expressed in the target samples. Additionally, predicted transcripts resulting from a subsequent analysis of the complex clusters were included at this stage. Of the 2,612 complex cluster that were not included in the design for the first batch of arrays, only 223 contained more than 2 sequences and for a further 14 no subsequence match of at least 60 nucleotides could be found within the last 1000 bases at the 3'-end. These 237 clusters were manually curated while the rest could be automatically reduced to one sequence. To make full use of the 15,208 features available on the Agilent microarray platform, it was decided to also include predicted sequences with somewhat lower gene prediction score that showed a hybridization signal in at least 8 of the 12 arrays. Finally, a selected set of 15,253 predicted transcript sequences was used as targets for probe design of a comprehensive <italic>P. pastoris </italic>microarray. While it is obvious that this list is larger than the expected number of open reading frames (6,000–7,000), as judged in comparison to other yeast species [##REF##17868481##22##], we intentionally included more putative transcript sequences, as false positives with a distinct sequence will not negatively affect microarray design or experiments, in contrast to the damage of falsely excluding a potential transcript target.</p> <p>Oligonucleotide probes were designed using a probe design tool developed in-house, a thermodynamic model-based oligoset optimizer ('TherMODO', [##UREF##2##23##]). TherMODO designed probes for 15,035 sequences, of which only 665 were predicted as having cross-hybridization potential. The TherMODO design was compared to probe design with eArray [##UREF##3##24##]. The distributions of ΔG and T<sub>m </sub>of both designs are shown in additional file ##SUPPL##0##1##. Clearly the TherMODO designed probes are more uniform in respect to the Gibbs free energy ΔG, indicating a superior hybridization performance [##REF##16938547##14##].</p> <p>The final probe design was manufactured on 8 × 15 K slides by Agilent, and evaluated for reproducibility and biological meaningfulness. Pool samples were applied to 2 arrays on 2 slides each, including dye swap. The scatterplots show uniformly high correlations &gt; 97% both within and between arrays, both on same and different slides, indicating high reproducibility of hybridization signals between identical samples. Exemplarily, a scatterplot of signal intensities derived from the same samples (wild type strain untreated) is shown in Figure ##FIG##0##1##. For the final gene list the annotation was improved in addition to the annotation provided by IG. This resulted in 3954 annotated ORFs, of which 2989 had an IG annotation. 965 newly annotated ORFs were found, and the annotation of 288 hypothetical proteins was confirmed. All annotated genes are listed in Additional file ##SUPPL##1##2##.</p> <title>Biological evaluation of the new microarrays</title> <p>The performance of the new arrays was examined by a hybridization experiment using samples, for which transcript regulation data have been obtained before [##REF##17578563##4##]. The biological question evaluated was the regulatory response of <italic>P. pastoris </italic>to constitutive overexpression of the active form of <italic>S. cerevisiae HAC1</italic>, the transcription factor controlling UPR target genes. By this approach, the regulation of 52 genes which have been studied before using TRAC [##REF##16198435##3##] could be verified, with 80% of these genes showing the same regulation pattern for both methods (genes highlighted in bold in Additional file ##SUPPL##1##2##). This correlation is statistically significant based on calculating the regression (p = 8.8 · 10<sup>-6</sup>).</p> <p>The similarities and differences of UPR induction and reaction to DTT stress have been discussed before [##REF##17578563##4##,##REF##17561995##25##,##REF##15870366##26##]. To achieve further insight into this technologically relevant issue, we compared the gene regulation patterns of a <italic>HAC1 </italic>overexpressing strain <italic>vs </italic>wildtype control with the regulation pattern of the wildtype treated with DTT for 60 min <italic>vs </italic>the untreated control. Genes were qualified as significantly regulated with a <italic>p</italic>-value &lt; 0.05 (adjusted for multiple testing). 11,262 of all features on the microarrays appeared as differentially regulated either upon DTT treatment or <italic>HAC1 </italic>overexpression, or both. 8,480 reacted to <italic>HAC1</italic>, and 6,870 to DTT, with an overlap of 4,088. Considering only the 3,954 annotated genes, a similar pattern is observed with roughly half of the regulated genes overlapping between DTT and <italic>HAC1</italic>, and another half being typical only for either of the treatments (Figure ##FIG##1##2##). Accordingly, the correlation of log fold changes of the two treatments is apparent but rather weak (Figure ##FIG##2##3##). While DTT treatment is widely accepted as a standard inducer of UPR, these observations indicate that the gene regulation pattern triggered by the UPR transcription factor Hac1 differs to a significant extent from that exerted by DTT.</p> <p>As previous research on transcriptome regulation upon UPR induction usually employs a fold change (FC) cut-off to highlight the strongly regulated genes, we decided to introduce FC &gt; 1.5 as a second criterion to identify more strongly regulated genes for further detailed analysis (Volcano plots visualizing the two criteria are provided in Additional file ##SUPPL##2##3##). Although the introduction of a FC cut-off alters the absolute number of regulated genes, it does not alter the relative distribution of regulated genes categorized into functional groups (GO slim biological process), as can be seen in Figure ##FIG##3##4## and Additional file ##SUPPL##3##4##.</p> <title>Comparison of UPR induction by DTT in <italic>P. pastoris </italic>and <italic>S. cerevisiae</italic></title> <p>In order to compare the effects of DTT treatment in <italic>S. cerevisiae </italic>with those in <italic>P. pastoris</italic>, the data published by Travers et al. [##REF##10847680##17##] for 60 min treatment of <italic>S. cerevisiae </italic>with DTT were evaluated alongside with our results for <italic>P. pastoris</italic>. All genes of <italic>S. cerevisiae </italic>which were listed in [##REF##10847680##17##] and for which homologs in <italic>P. pastoris </italic>were identified were classified as upregulated, downregulated or unregulated. In order to compare the two data sets, a cutoff of 1.5 fold differential expression was set in both to define regulated genes. A significance threshold on <italic>p</italic>-values could not be employed, as these data were not provided for <italic>S. cerevisiae</italic>. 48% of these genes defined as regulated or unregulated reacted in <italic>P. pastoris </italic>just as in <italic>S. cerevisiae</italic>.</p> <p>A closer evaluation revealed that certain GO groups were regulated very similarly in both yeast species, while others showed only a low degree of similarity (Table ##TAB##1##2##). Fisher's exact test was performed to evaluate the significance of groups with low similarity. Especially the GO groups 'translocation', 'protein folding', 'protein degradation', and to some extent 'glycosylation' and 'transport' showed high degrees of similarity. In some GO groups, only some subgroups reacted similarly while others behaved differently in the two yeasts. Of the 'glycosylation' group, core oligosaccharide synthesis and glycosyltransferase genes behaved very similarly, while glycoprotein processing, GPI anchoring and O-glycosylation related genes were regulated significantly different (p &lt; 0.05). In the 'protein degradation' group, more similarity was observed for ERAD genes than for ubiquitin/proteasome related genes. Among the 'transport' gene group, budding, fusion and retrieval of ER to Golgi showed a high degree of similar regulation, contrary to the subgroup distal secretion. Low similarities were observed for 'lipid metabolism', 'vacuolar protein sorting' and 'cell wall biogenesis' genes. It becomes obvious that core UPR genes related to protein translocation, folding and ER transport, as well as core N-glycosylation react similarly to DTT treatment in <italic>P. pastoris </italic>as compared to <italic>S. cerevisiae</italic>, while genes involved in processes which are more distal from ER protein folding behave more differently, indicating that those processes (like functions in the Golgi, [##REF##10189369##27##]) differ significantly between the two yeasts.</p> <title>Overexpression of Hac1 triggers a different regulation pattern compared to DTT treatment</title> <p>In most previous studies of the UPR in lower eukaryotic cells, treatment with DTT or tunicamycin, or heterologous protein expression has been employed to trigger the UPR. This study clearly indicates that the set of regulatory events triggered by DTT analysis only partially overlaps with the reactions to constitutive expression of the activated form of the UPR transcription factor Hac1 (see Figures ##FIG##1##2## and ##FIG##2##3##). Interestingly, both treatments resulted in the same amount of genes being down-regulated as being up-regulated, a fact that has been neglected to some extent in the existing literature.</p> <p>Those genes appearing beyond the threshold (<italic>p</italic>-value &lt; 0.05 and FC &gt;1.5) were subjected to a more detailed comparison between the effects of DTT treatment and Hac1 induced regulation. The relative numbers of up- and downregulated genes in each GO biological process term based on the SGD GO slim tool [##UREF##4##28##] are depicted in Figure ##FIG##3##4##.</p> <p>A pattern common to both treatments is the down-regulation of major metabolic processes like carbohydrate, amino acid and lipid metabolism, as well as that of vitamins, cofactors and aromatic and heterocyclic compounds. This makes it obvious that the UPR has a major impact on decreasing both catabolic and anabolic processes. On the other side, both treatments lead to up-regulation of protein folding and vesicular transport. These effects are in line with the published literature, indicating the cellular reaction towards alleviation of the UPR [##REF##17578563##4##,##REF##17561995##25##,##REF##15870366##26##,##REF##10847680##17##].</p> <p>As expected, the genes coding for classical UPR targets are induced both in Hac1 overproducing and in DTT stressed cells, and genes underlined in the following paragraphs have been identified as UPR targets in previous studies. Especially the ER folding catalysts <italic><underline>PDI1</underline></italic> and <italic><underline>ERO1</underline></italic>, the DnaJ homologs <italic><underline>JEM1</underline></italic> and <italic><underline>SCJ1</underline></italic>, the ER resident chaperones <italic>CNE1 </italic>(calnexin), <underline><italic>KAR2</italic>/BiP</underline> and <italic><underline>LHS1</underline></italic> and the mitochondrial chaperones <italic>HSP60 </italic>and <italic>SSC1 </italic>are significantly up-regulated in both conditions. Among the functional group of 'protein modification' the majority of up-regulated genes belong to the core oligosaccharide synthesis (<italic><underline>DPM1</underline></italic>, <italic>DIE2</italic>), oligosaccharyltransferase complex (<italic>OST1</italic>, <italic><underline>OST2</underline></italic>, <italic><underline>OST3</underline></italic>, <italic><underline>SWP1</underline></italic>, <italic>STT3</italic>, <italic><underline>WBP1</underline></italic>), glycoprotein processing (<italic>ALG2, <underline>ALG7</underline></italic>, <italic>SEC53</italic>), GPI anchor biosynthesis (<italic><underline>GPI2</underline></italic>, <italic>GPI14, PSA1</italic>) and Golgi/O-linked glycosylation (<italic><underline>PMT1</underline></italic>, <italic><underline>PMT2</underline></italic>, <italic>PMT4, PMT6</italic>). Besides these, several genes coding for the translocon pore complex (<italic><underline>SEC61</underline></italic>, <italic><underline>SEC62</underline></italic>, <italic>SEC63</italic>, <italic><underline>SEC72</underline></italic>, <italic>SSS1</italic>), which aid the translocation of nascent polypeptides into the ER, are induced. Higashio and Kohno [##REF##12176018##29##] describe the stimulation of ER-to-Golgi transport through the UPR by inducing COPII vesicle formation. In this context, we see <italic>SEC23</italic>, <italic><underline>SEC24</underline></italic>, <italic><underline>SFB2</underline></italic>, <italic>YIP3</italic>, and <italic>ERV2 </italic>upregulated. However, also proteins building the COPI coatomer, which are required for retrograde Golgi-to-ER transport, show increased transcription levels upon ER stress in our experiments (<italic>COP1</italic>, <italic><underline>RET2</underline></italic>, <italic>SEC21</italic>, <italic><underline>SEC27</underline></italic>).</p> <p>While we cannot give any information on ERAD regulation, as <italic><underline>HRD1</underline></italic> is the only annotated gene of this protein degradation process (up-regulated in the Hac1 strain), we observed the down-regulation of some components involved in the assembly of the 20 S core of the 26 S proteasome (<italic>ADD66</italic>, <italic>PRE1</italic>, <italic>PRE4</italic>, <italic>SCL1</italic>) and ubiquitin <italic>UBI4 </italic>upon constitutive UPR activation. In this context, Shaffer et al. [##REF##15345222##30##] describe reduced degradation of newly synthesized proteins in XBP1-overexpressing human Raji cells.</p> <p>Induction of genes encoding cytosolic chaperones (Cns1, Jjj3, Hsp82, Ssa1, Ssa2, Sse1, Ydj1, Zuo1) can only be seen in the Hac1-overproducing strain. Additionally, the ER-resident Pdi homolog <underline>Mpd1</underline> and two members of the PPIases (<italic>FPR4 </italic>and <italic>CPR6</italic>) are only up-regulated in the engineered strain, but not upon DTT addition.</p> <p>One of the most striking patterns is the significant up-regulation of a large number of genes with functions in ribosomal biogenesis (233 genes assigned to the GO-categories 'ribosome biogenesis and assembly' and 'RNA metabolic process'). Most of these genes are contributing to rRNA processing (RRP family) and ribosome subunit nuclear export and assembly, while the ribosomal proteins (RPS and RPL families) themselves are not among the regulated genes for <italic>P. pastoris </italic>(see Additional file ##SUPPL##1##2##). No genes with a function in mRNA decay show increased transcription levels. The induction of the above functional categories came as a surprise, as translational down-regulation of proteins involved in ribosomal biogenesis was recently reported when <italic>S. cerevisiae </italic>cell were treated with DTT [##REF##18206654##31##]. In contrast, the transcription levels of 9 out of the 16 mRNAs listed by these authors are enhanced in our study. Transcriptional down-regulation of ribosomal proteins during ER stress conditions was also revealed when reanalysing the raw data provided by Travers et al. [##REF##10847680##17##]. However, Shaffer et al. [##REF##15345222##30##] describe an increase in total protein synthesis as well as in the number of assembled ribosomes upon the overexpression of the mammalian Hac1 homolog XBP1 in Raji cells, but did not observe upregulation of genes related to ribosome biogenesis. A similar effect was observed after XBP1 overexpression in CHO-K1 cells [##REF##16635796##32##]. These results may be an indication that the positive effect of overexpression of the UPR transcription factor on heterologous protein production [##REF##16570317##33##,##REF##12676684##16##,##REF##14660339##35##] results not just from stimulation of folding and secretion of proteins but also their synthesis. The induction of protein folding related genes upon Hac1 overexpression is in line with the literature on UPR effects, while an impact on organelle biosynthesis other than ER and Golgi has so far only been described for mammalian cells.</p> <p>The stimulatory effects of XBP1 induction on ribosomes and organelle synthesis in mammalian cells like lymphocytes have been attributed to their function as dedicated protein factories. On the other hand the UPR in lower eukaryotes should rather serve to alleviate the load of unfolded, aggregation prone protein. It will be of interest in the future to investigate whether Hac1 stimulates ribosome biogenesis in other yeasts and fungi as well, and whether this leads to increased translation.</p> <p>In this context, it is worthwhile to mention the induction of two pathways leading to the unusual post-translationally modified amino acid derivatives diphthamide and hypusine which are exclusively found in eukaryotic translation elongation factors 2 (eEF2) and 5 (eEF5), respectively [##REF##8406038##36##,##REF##17476569##37##]. As these biosynthetic pathways are rather complex, and outstanding in the otherwise downregulated group of 'amino acid biosynthesis', this induction underlines the increased demand for protein synthesis.</p> <p>Furthermore, we observe that ER stress leads to increased transcription of genes coding for the large and small subunits of the mitochondrial ribosomes (<italic>MRPS</italic>, <italic>RSM </italic>and <italic>MRPL </italic>families), mitochondrial translation initiation and elongation factors (<italic>IFM1</italic>, <italic>MEF1</italic>, <italic>MEF2</italic>) and mitochondrial DNA polymerase (<italic>MIP1</italic>). Several essential constituents of the mitochondrial inner membrane presequence translocase (<italic>TIM </italic>family) are also up-regulated, indicating increased necessity for protein import into the mitochondria. Similarly, XBP1 was shown to increase mitochondrial mass and function in two types of mammalian cells [##REF##15345222##30##].</p> <p>While previous studies analysing UPR regulation mainly focus on up-regulated genes [##REF##10847680##17##], more than half of the genes identified in our study to be regulated are strongly down-regulated (at least 1.5 fold). As can be seen in Figure ##FIG##3##4##, anabolic processes such as vitamin production, amino acid and aromatic compound biosynthesis, heterocycle metabolic processes, carbohydrate, lipid and cofactor metabolism are among the most prominent repressed classes in both DTT-treated as well as Hac1-overproducing cells. The down-regulation of energy consuming biosynthetic pathways emerges as a general picture during ER stress conditions. However, it becomes obvious that the response to the folding perturbation agent DTT strongly differs from constitutive UPR induction by Hac1-overproduction. Especially the prominent down-regulation of genes belonging to 'electron transport' and 'cellular respiration' can easily be explained by the strong reducing capacities of DTT. Prominent members of the mitochondrial inner membrane electron transport chain such as subunits of the cytochrome c oxidase (<italic>COX4, COX4, COX5A, COX13</italic>) and the ubiquinol cytochrome-c reductase complex (<italic>COR1, QRC6, QRC7, QRC9, RIP1</italic>) are significantly repressed upon DTT treatment. Additionally, cytochrome c (<italic>CYC1</italic>), cytochrome c1 (<italic>CYT1</italic>) and cytochrome c heme lyase (<italic>CYC3</italic>) are only under DTT-dependent repression (GO: 'generation of precursor metabolites and energy'). The reducing features of DTT are most probably also the reason for the up-regulation of genes involved in the upkeeping of 'cellular homeostasis' and clearly, addition of DTT is provoking a 'response to a chemical stimulus'.</p> <p>Down-regulated genes appearing in both Hac1 and DTT in the 'protein modification' group focus on protein kinases (<italic>CDC5, CDH1, DBF2</italic>) and components of the ubiquitinylation complex (<italic>BUL1, CUL3</italic>) involved in cell cycle regulation driving the cells towards mitotic exit (<italic>CDC5, CDH1</italic>, <italic>MOB1</italic>). These effects are even more pronounced in the Hac1-strain, where several more histone modifying enzymes as well as cycline-dependent protein kinases and components of the protein kinase C signalling pathway show reduced transcription levels compared to the wild type. Unlike reported for the filamentous fungi <italic>T. reesei </italic>[##REF##16504068##7##] and <italic>A. nidulans </italic>[##REF##15870366##26##], genes encoding the histones H2A, H2B, H3 and H4 appear to be down-regulated upon secretion stress in <italic>P. pastoris</italic>.</p> <p>No clear picture emerges regarding the regulation of 'lipid metabolism': While sterol and ergosterol biosynthesis tend to be inhibited, the production of sphingolipid precursor substances is enhanced. On the other hand, a down-regulation of the major cell wall constituents (β-1,3 glucanases <italic>BGL2 </italic>and <italic>EXG1</italic>, cell wall mannoproteins <italic>CCW12</italic>, <italic>CWP2 </italic>and <italic>TPI1</italic>, GPI-glycoproteins <italic>GAS1 </italic>and <italic>SED1</italic>, <italic>PST1</italic>) and genes coding for proteins required for the transport of cell wall components to the cell surface (<italic>SBE22</italic>) is manifest. Taken together, these results indicate a significant remodelling process regarding the <italic>P. pastoris </italic>cell envelope during ER stress conditions.</p> <p>Interestingly, the major groups of metabolic genes were down-regulated upon Hac1 overexpression, indicating a decrease of the supply of metabolites. However, it should be noted that no reduction of the specific growth rate was observed as compared to the wild type strain (μ = 0.37 and 0.39 h<sup>-1</sup>, respectively). A reduction of metabolic processes, and amino acid synthesis in particular, is contradictory to translation stimulation. Further research will be needed to elucidate the overall regulatory pattern of UPR in respect to protein synthesis.</p>
<title>Conclusion</title> <p>Additional gene finding and annotation added to the available data for <italic>P. pastoris </italic>lead to a list of approximately 4,000 genes with a putative identification of their function, and 11,000 more potential open reading frames. An oligonucleotide probe set was designed, the hybridization results were evaluated for reproducibility, and results from a biologically relevant analysis were tested for meaningfulness. In a direct comparison to <italic>S. cerevisiae </italic>employing DTT treatment for UPR induction, 45 out of 93 genes reacted similarly. The differences thus observed between <italic>P. pastoris </italic>and <italic>S. cerevisiae </italic>underline the importance of DNA microarrays for industrial production strains. <italic>HAC1 </italic>overexpression in <italic>P. pastoris </italic>obviously leads to induction of many genes involved in translation: most genes of ribosome biogenesis, as well as many related to RNA metabolism and translation were up-regulated, an effect that has never been observed in yeasts and filamentous fungi so far.</p> <p>The upregulation of ribosomal biogenesis, RNA metabolism, translation, and organelle biosynthesis is specific for <italic>HAC1 </italic>overexpression and not observed with DTT treatment, while the latter leads specifically to the upregulation of genes related to chemical stimulus, and the downregulation in the groups electron transport and respiration, so that these reactions have to be regarded as specific for the treatment with a reducing agent rather than UPR regulated.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>DNA Microarrays are regarded as a valuable tool for basic and applied research in microbiology. However, for many industrially important microorganisms the lack of commercially available microarrays still hampers physiological research. Exemplarily, our understanding of protein folding and secretion in the yeast <italic>Pichia pastoris </italic>is presently widely dependent on conclusions drawn from analogies to <italic>Saccharomyces cerevisiae</italic>. To close this gap for a yeast species employed for its high capacity to produce heterologous proteins, we developed full genome DNA microarrays for <italic>P. pastoris </italic>and analyzed the unfolded protein response (UPR) in this yeast species, as compared to <italic>S. cerevisiae</italic>.</p> <title>Results</title> <p>By combining the partially annotated gene list of <italic>P. pastoris </italic>with <italic>de novo </italic>gene finding a list of putative open reading frames was generated for which an oligonucleotide probe set was designed using the probe design tool TherMODO (a thermodynamic model-based oligoset design optimizer). To evaluate the performance of the novel array design, microarrays carrying the oligo set were hybridized with samples from treatments with dithiothreitol (DTT) or a strain overexpressing the UPR transcription factor <italic>HAC1</italic>, both compared with a wild type strain in normal medium as untreated control. DTT treatment was compared with literature data for <italic>S. cerevisiae</italic>, and revealed similarities, but also important differences between the two yeast species. Overexpression of <italic>HAC1</italic>, the most direct control for UPR genes, resulted in significant new understanding of this important regulatory pathway in <italic>P. pastoris</italic>, and generally in yeasts.</p> <title>Conclusion</title> <p>The differences observed between <italic>P. pastoris </italic>and <italic>S. cerevisiae </italic>underline the importance of DNA microarrays for industrial production strains. <italic>P. pastoris </italic>reacts to DTT treatment mainly by the regulation of genes related to chemical stimulus, electron transport and respiration, while the overexpression of <italic>HAC1 </italic>induced many genes involved in translation, ribosome biogenesis, and organelle biosynthesis, indicating that the regulatory events triggered by DTT treatment only partially overlap with the reactions to overexpression of <italic>HAC1</italic>. The high reproducibility of the results achieved with two different oligo sets is a good indication for their robustness, and underlines the importance of less stringent selection of regulated features, in order to avoid a large number of false negative results.</p>
<title>Authors' contributions</title> <p>AG performed gene finding and annotation, statistical data analysis, supported data evaluation, and drafted part of the manuscript. BG performed data evaluation, supported annotation, study design, array design and drafted part of the manuscript. MD performed the cultivations and hybridizations. MS contributed to study design, annotation and array design. GGL developed the employed probe design tool and supported the array design. TT developed the quantitative model for the position-dependent target labelling efficiency and adapted it for the relevant end-primed labelling protocol. DPK supervised gene identification and annotation, supervised and contributed to the development of the employed probe design tool, and contributed to the manuscript. DM conceived of the study, and participated in data evaluation and manuscript drafting. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported by the Austrian Science Fund (project No. I37-B03), the European Science Foundation (programme EuroSCOPE), and the Austrian Research Promotion Agency (programme FHplus).</p> <p>The Vienna Science Chair of Bioinformatics gratefully acknowledges support by the Vienna Science and Technology Fund (WWTF), Baxter AG, Austrian Research Centres (ARC) Seibersdorf, and Austrian Centre of Biopharmaceutical Technology (ACBT). TT acknowledges partial funding by the GenAu BIN-II PhD programme.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Correlation of signal intensities</bold>. Scatterplots of untreated wild type strain samples on (A) different arrays of the same slide; (B) different arrays on different slides. Red line: linear regression of the data; blue line: theoretical perfect correlation.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Venn diagrams of differentially expressed genes upon DTT treatment or <italic>HAC1</italic> overexpression</bold>. (A, B) Regulated hits with annotation; (C, D) all regulated features; (A, C) cut-off adjusted <italic>p</italic>-value &lt; 0.05; (B, D) cut-off adjusted <italic>p</italic>-value &lt; 0.05 and FC &gt; 1.5.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Comparison of expression changes induced by DTT treatment and <italic>HAC1</italic> overexpression, respectively</bold>. Log<sub>2 </sub>values of expression changes (log<sub>2 </sub>FC) caused by DTT (DTT treated wildtype <italic>vs </italic>untreated wildtype) and by Hac1 (<italic>HAC1 </italic>overexpression <italic>vs </italic>wildtype) are compared. The correlation coefficient r<sup>2 </sup>is indicated. Red line: linear regression of the data; blue line: theoretical perfect correlation.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Fractions of up- and downregulated genes in functional groups</bold>. Relative numbers of upregulated (red), downregulated (blue) and unregulated (yellow) genes categorized in GO biological process terms upon <italic>HAC1 </italic>overexpression (left panel) and DTT treatment (right panel). Shaded in black: regulated in both treatments. Upper panels: cut-off <italic>p</italic>-value &lt; 0.05, lower panel cut-off <italic>p</italic>-value &lt; 0.05 and FC &gt; 1.5. The results of significance testing (Fisher's exact test) are given in additional file ##SUPPL##3##4##.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Comparison of gene finder performance on yeast genomic sequence data</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene finder</td><td align="left">True positives</td><td align="left">Partly</td><td align="left">False positives</td><td align="left">False negatives</td><td align="left">Sensitivity (%)</td><td align="left">Positive prediction value (%)</td></tr></thead><tbody><tr><td align="left">Glimmer3</td><td align="left">75</td><td align="left">13</td><td align="left">21</td><td align="left">31</td><td align="left">73.9</td><td align="left">68.8</td></tr><tr><td align="left">GlimmerHMM</td><td align="left">1</td><td align="left">3</td><td align="left">68(234)</td><td align="left">115</td><td align="left">3.2</td><td align="left">1.4</td></tr><tr><td align="left">GeneMark</td><td align="left">81</td><td align="left">16</td><td align="left">32</td><td align="left">22</td><td align="left">81.5</td><td align="left">62.7</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Similarity of gene regulation between <italic>P. pastoris</italic> and <italic>S. cerevisiae</italic> upon DTT treatment</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Function</td><td align="left">Subfunction</td><td align="left">No. of similarly regulated/total</td><td align="center">% similar regulation</td></tr></thead><tbody><tr><td align="left">Translocation</td><td align="left">total</td><td align="left">4/6</td><td align="center">67</td></tr><tr><td align="left">Glycosylation</td><td align="left">total</td><td align="left">11/22</td><td align="center">50</td></tr><tr><td/><td align="left">Core oligosaccharide synthesis</td><td align="left">3/4</td><td align="center">75</td></tr><tr><td/><td align="left">Oligosaccharyltransferase</td><td align="left">4/4</td><td align="center">100</td></tr><tr><td/><td align="left">Glycoprotein processing</td><td align="left">1/5</td><td align="center">20</td></tr><tr><td/><td align="left">GPI anchoring</td><td align="left">1/4</td><td align="center">25</td></tr><tr><td/><td align="left">Golgi/O-linked</td><td align="left">2/5</td><td align="center">40</td></tr><tr><td align="left">Protein Folding</td><td align="left">total</td><td align="left">5/8</td><td align="center">63</td></tr><tr><td/><td align="left">Chaperones</td><td align="left">3/5</td><td align="center">60</td></tr><tr><td/><td align="left">Disulfide bond formation</td><td align="left">2/3</td><td align="center">67</td></tr><tr><td align="left">Protein Degradation</td><td align="left">total</td><td align="left">4/5</td><td align="center">80</td></tr><tr><td/><td align="left">ERAD</td><td align="left">3/3</td><td align="center">100</td></tr><tr><td/><td align="left">Ubiquitin/Proteasome</td><td align="left">1/2</td><td align="center">50</td></tr><tr><td align="left">Transport</td><td align="left">total</td><td align="left">11/20</td><td align="center">55</td></tr><tr><td/><td align="left">Budding (ER-Golgi)</td><td align="left">4/7</td><td align="center">57</td></tr><tr><td/><td align="left">Fusion (ER-Golgi)</td><td align="left">1/1</td><td align="center">100</td></tr><tr><td/><td align="left">Retrieval (ER-Golgi)</td><td align="left">4/5</td><td align="center">80</td></tr><tr><td/><td align="left">Distal secretion</td><td align="left">2/7</td><td align="center">29</td></tr><tr><td align="left">Lipid Metabolism</td><td align="left">total</td><td align="left">5/18</td><td align="center">28</td></tr><tr><td/><td align="left">Fatty acid metabolism</td><td align="left">0/4</td><td align="center">0</td></tr><tr><td/><td align="left">Heme biosynthesis</td><td align="left">2/5</td><td align="center">40</td></tr><tr><td/><td align="left">Phospholipid biosynthesis</td><td align="left">2/6</td><td align="center">33</td></tr><tr><td/><td align="left">Sphingolipid biosynthesis</td><td align="left">0/1</td><td align="center">0</td></tr><tr><td/><td align="left">Sterol metabolism</td><td align="left">1/2</td><td align="center">50</td></tr><tr><td align="left">Vacuolar Protein Sorting</td><td align="left">total</td><td align="left">1/4</td><td align="center">25</td></tr><tr><td align="left">Cell Wall Biogenesis</td><td align="left">total</td><td align="left">4/10</td><td align="center">40</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional File 1</title><p><bold>Thermodynamic properties of the TherMODO probe design compared to probes designed through Agilent's eArray</bold>. Distribution of Gibbs free energy ΔG (A) and the probe-target melting temperature T<sub>m </sub>(B) of the oligo sets. The upper row (1) shows the oligos designed through eArray and the lower row (2) the oligos designed with TherMODO. PpaV2 is the name of the second set of sequences as described in the Materials and Methods section.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional File 2</title><p><bold>Differential expression values of all annotated genes upon DTT treatment and <italic>HAC1 </italic>overexpression</bold>. Differential expression values and adjusted <italic>p</italic>-values of all annotated genes of <italic>P. pastoris</italic>, denominated with the gene name of their respective <italic>S. cerevisiae </italic>homolog. Genes that were tested with TRAC as a different method for transcript quantification are highlighted in bold letters. Legend of headers: id – internal unique identifier of sequence; sequ_id – ERGO identifier (RPPA.) or gene finder identifier (orf.) respectively; DTT_logFC – log<sub>2</sub> fold change of DTT treatment compared to control; <italic>HAC1 </italic> logFC – log<sub>2</sub> fold change of <italic>HAC1 </italic> overexpression compared to control; Gene name – Standard gene name or if missing systematic ORF name according to <italic>S. cerevisiae </italic>nomenclature; GO – Gene Ontology term (for descriptions see additional file ##SUPPL##3##4##). If a gene is present in a certain GO group it has a 1 in the respective column, if not it has a 0.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional File 3</title><p><bold>Volcano plots of fold change vs. adjusted <italic>p</italic>-values</bold>. (A) DTT treatment; (B) <italic>HAC1 </italic>overexpression. Blue line: <italic>p</italic>-value cut-off <italic>p </italic>&gt; 0.05; red lines: optional fold change cut-off FC &gt; 1.5.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional File 4</title><p><bold>Fisher's exact test of the up/down regulated gene groups upon DTT treatment and <italic>HAC1 </italic>overexpression</bold>. Fisher's exact test was applied to test significance of the up- and downregulated gene groups displayed in figure ##FIG##3##4##. <italic>p</italic><sub><italic>adj </italic></sub>values are given for each GO group. Legend of headers: group – Gene Ontology term; Description – Gene Ontology description; odds.ratio – measure of independence between variables; adj.p – Holm adjusted <italic>p</italic>-value; <italic>HAC1 </italic> up/down – up/down regulated in <italic>HAC1 </italic> overexpression experiment; DTT up/down – up/down regulated in DTT experiment. The first work sheet represents results using only a <italic>p</italic>-value cut-off <italic>p </italic>&gt; 0.05, the second work sheet represents results using a <italic>p</italic>-value cut-off <italic>p </italic>&gt; 0.05 and a fold change cut-off FC &gt; 1.5.</p></caption></supplementary-material>
<table-wrap-foot><p>Three different gene finders were tested on the genome sequence of <italic>S. cerevisiae </italic>chromosome 1 to evaluate the quality of gene prediction. Sensitivity = TP/(TP + FN), positive prediction value = TP/(TP + FP); For Glimmer HMM the column False Positives contains the number of genes and in brackets the number of exons.</p></table-wrap-foot> <table-wrap-foot><p>All genes that were indicated in [##REF##10847680##17##] as core UPR genes in <italic>S. cerevisiae </italic>and having an annotation in <italic>P. pastoris </italic>were grouped by their GO process functions. Similar regulation of a gene means upregulated, downregulated or below cut-off, respectively, in both yeasts.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-390-1"/> <graphic xlink:href="1471-2164-9-390-2"/> <graphic xlink:href="1471-2164-9-390-3"/> <graphic xlink:href="1471-2164-9-390-4"/>
<media xlink:href="1471-2164-9-390-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-390-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-390-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-390-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Genomics"], "given-names": ["I"], "article-title": ["ERGO bioinformatics suite"]}, {"surname": ["University"], "given-names": ["W"], "article-title": ["WU-BLAST"]}, {"surname": ["Leparc", "Tuechler", "Striedner", "Bayer", "Sykacek", "Hofacker", "Kreil"], "given-names": ["GG", "T", "G", "K", "P", "I", "DP"], "article-title": ["Model based probe set optimization for high-performance microarrays."], "source": ["Nucleic Acids Research"], "year": ["2008"], "volume": ["submitted"]}, {"collab": ["Agilent"], "article-title": ["eArray"]}, {"collab": ["SGD"], "article-title": ["SGD Gene Ontology Slim Mapper"]}, {"surname": ["(EBI)", "(CSHL)"], "given-names": ["EBI", "CSHL"], "article-title": ["BioMart"]}, {"surname": ["Technology"], "given-names": ["GI"], "article-title": ["GeneMark"]}, {"surname": ["(SIB)"], "given-names": ["SIB"], "article-title": ["Swiss-Prot/TrEMBL"]}, {"surname": ["Holm"], "given-names": ["S"], "article-title": ["A Simple Sequentially Rejective Bonferroni Test"], "source": ["Scandinavian Journal of Statistics"], "year": ["1979"], "fpage": ["65 "], "lpage": ["670"]}, {"collab": ["EMBL-EBI"], "article-title": ["ArrayExpress"]}]
BMC Genomics. 2008 Aug 19; 9:390