# Isolation and analysis of differentially expressed genes in Penicillium glabrum subjected to thermal stress

Laurent Nevarez, Valérie Vasseur, Gwenola Le Drean, Arnaud Tanguy, Isabelle Guisle-Marsollier, Rémi Houlgatte, Georges Barbier

► To cite this version:

Laurent Nevarez, Valérie Vasseur, Gwenola Le Drean, Arnaud Tanguy, Isabelle Guisle-Marsollier, et al. Isolation and analysis of differentially expressed genes in Penicillium glabrum subjected to thermal stress. Microbiology, 2008, 154 (154), pp.3752-3765. 10.1099/mic.0.2008/021386-0. hal-00551256

HAL Id: hal-00551256

https://hal.univ-brest.fr/hal-00551256

Submitted on 3 Jan 2011

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Isolation and analysis of differentially expressed genes in
Penicillium glabrum submitted to thermal stress

L. Nevarez<sup>a*</sup>, V. Vasseur<sup>a</sup>, G. Le Dréan<sup>a</sup>, A. Tanguy<sup>b</sup>, I. Guisle-Marsollier<sup>c</sup>,
R. Houlgatte<sup>c</sup>, G. Barbier<sup>a</sup>

<sup>a</sup> Laboratoire Universitaire de Biodiversité et Ecologie Microbienne, Université Européenne de Bretagne, Ecole Supérieure de
Microbiologie et Sécurité Alimentaire de Brest, Technopôle Brest-Iroise, 28280 Plouzané, France

<sup>b</sup> Evolution et Génétique des Populations Marines, UMR CNRS 7144, Université Pierre et Marie Curie, Station Biologique de Roscoff,
Place Georges Teissier, 29682 Roscoff Cedex, France

<sup>c</sup> Plate-forme Transcriptomique Ouest-Génopôle, Institut du Thorax INSERM U533, 1 Rue Gaston Veil, BP 53508, 44035 Nantes

Cedex 1, France

<footnote><sup>*</sup> Corresponding author. Tel: +33 2 98 05 61 26; Fax : +33 2 98 05 61 01</footnote>

<footnote>E-mail address: laurent.nevarez@univ-brest.fr</footnote>
Summary

Penicillium glabrum is a filamentous fungus frequently involved in food contamination. Numerous environmental factors (temperature, humidity, atmosphere composition, etc.) or food characteristics (water activity, pH, preservatives, etc.) could represent potential sources of stress for micro-organisms. These factors can directly affect physiology of these spoilage micro-organisms: growth, conidiation, synthesis of secondary metabolites, etc. This study investigates the transcriptional response to temperature in P. glabrum, since this factor is one of the most important for fungal growth. Gene expression was first analysed by using suppression subtractive hybridization to generate two libraries containing 445 different up- and down-regulated Expressed Sequenced Tags (ESTs). Expression of these ESTs was then assessed for different thermal stress conditions, with cDNA microarrays, resulting in the identification of 35 and 49 significantly up- and down-regulated ESTs, respectively. These ESTs encode heat shock proteins, ribosomal proteins, superoxide dismutase, trehalose-6-phosphate synthase and a large variety of identified or unknown proteins. Some of them may be potential molecular markers for thermal stress response in P. glabrum. To our knowledge, this work represents the first study of the transcriptional response of a food spoilage filamentous fungus under thermal stress conditions.

**Keywords:** Penicillium glabrum, filamentous fungi, food spoilage, thermal stress, suppressive subtractive hybridization, cDNA microarrays, qRT-PCR, transcriptional study, gene regulation
Fungi are ubiquitous micro-organisms often associated with spoilage and
biodeterioration of a large variety of foods and feedstuffs. Penicillium is one of the most
widespread fungal genera isolated from food products. In addition to the economic losses
caused by these contaminants, many Penicillium species can produce mycotoxins that
represent a potential health risk for humans and animals (Pitt & Hocking, 1997; Samson et al.,
2004). Due to its ability to disperse a large number of spores in the environment, P. glabrum
is very frequently encountered in the food manufacturing industry (Pitt & Hocking, 1997).
This fungus is responsible for spoilage of many different products, such as cheese (Northolt et
al., 1980), stored cereals (Borjesson et al., 1992), maize (Mislivec & Tuite, 1970), rice
(Kurata et al., 1968), jam (Udagawa et al., 1977), commercially marketed chestnuts (Overy et
al., 2003), soda (Ancasi et al., 2006) and mineral water (Cabral & Fernandez Pinto, 2002).
Numerous intrinsic factors (water activity, pH, preservatives, etc.) in foods, provide ideal
conditions for micro-organisms to develop. Combined with extrinsic factors (temperature,
humidity, atmosphere composition, etc.), modification of these conditions are potential
sources of stress that can directly affect the physiology (growth, sporulation, synthesis of
secondary metabolites, etc.) of spoilage micro-organisms. In vitro data obtained with different
Penicillium strains, have shown significant variation in the production of secondary
metabolites among substrate (Chang et al., 1991; Kokkonen et al., 2005). Environmental
conditions may also modify mycotoxin production as reported in Fusarium species (Magan et
al., 2002). Some modification, such as osmotic stress have been shown to influence the
synthesis of secondary metabolites in P. notatum (Fiedurek, 1997).
Temperature and water activity are the most important factors that determine the ability of moulds to grow (Dantigny *et al.*, 2005). Fungi and more generally, every organism in each kingdom, have developed strategies to respond rapidly to environmental fluctuations. Among them, elevated temperatures are one of the foremost and best-known stress inducers at the cellular level. The general response to thermal stress is characterised by acquired thermotolerance, which increases survival capacity during subsequent exposure to high temperatures (Piper, 1993). Response to thermal stress has been described as a tolerance mechanism against adverse growth conditions and associated with synthesis of heat shock proteins (HSPs) (Lindquist, 1986). This mechanism aims to stabilize cell metabolism and prevents the accumulation of denatured proteins which can be toxic for the cell. Heat shock response appears to be widespread, spanning many taxonomic ranks and is found in drosophila (Ashburner & Bonner, 1979), rabbit (Cosgrove & Brown, 1983), *E. coli* (Yamamori & Yura, 1982), *Saccharomyces cerevisae* (McAlister & Finkelstein, 1980), *Aspergillus nidulans* (Newbury & Peberdy, 1996), etc. To our knowledge, no data are available on the heat shock transcriptional response of food-contaminating moulds, such as *P. glabrum*. From a food industry perspective, understanding the effects of various factors on the growth and physiology of a spoilage mould, is crucial to better manage the risk of alteration and toxicity of food products. It would therefore be useful to gain more insight on the mechanisms involved in thermal stress and to identify the genes that are differentially regulated in response to thermal stress.

In the present study, the effects of temperature on the transcriptional response of *P. glabrum* was investigated by using two complementary techniques: suppression subtractive hybridization (SSH) and cDNA microarrays. Transcriptional response to heat shock was compared between thermal stress conditions (40°C, 120 min) and the optimal growth temperature of 25°C. SSH is a powerful method based on suppression PCR where subtracted
cDNA libraries can be constructed to identify differentially expressed genes in response to an
experimental variation in physiological, environmental or other factors (Diatchenko et al.,
1996). This method has been previously used in fungi to isolate differentially expressed genes
associated with different morphologies or different growth conditions. It has been applied in
*Aspergillus nidulans* to characterise genes specifically expressed in conidiating cultures and
mature conidia (Osherov et al., 2002). SSH has been also used to investigate the differential
gene expression between mycelium and yeast forms of the human pathogens *Penicillium
marneffei* (Liu et al., 2007) and *Paracoccidioides brasiliensis* (Marques et al., 2004) with
respect to temperature (25° or 37°C).

Expressed Sequenced Tags (ESTs) isolated by SSH, were then investigated, in
different thermal stress conditions, by assessing their expression with the cDNA microarray
technique and confirmed with real-time reverse transcription PCR (qRT-PCR) on few selected
ESTs. The advent of cDNA microarrays has provided the possibility of analysing the global
changes in gene expression that occur in a cell, under different conditions (Shalon et al.,
1996). Transcriptional studies using whole-genome microarrays have been conducted on a
wide range of environmental conditions (e.g. thermal, osmotic and oxidative stresses) in the
yeasts S. cerevisiae (Gasch et al., 2000; Causton et al., 2001; Sakaki et al., 2003) and
*Schizosaccharomyces pombe* (Chen et al., 2003). As the whole genome of P. glabrum is not
available, we constructed a custom-made cDNA microarray and we report here the results
obtained. This work represents, to our knowledge, the first study investigating the
transcriptional response of a food spoilage fungus under thermal stress.
Biological material and culture conditions

The fungal strain used in this study was isolated from contaminated aromatized mineral
water and precisely characterised as *Penicillium glabrum* (Wehmer) Westling according to the
reference method used for classifying *Penicillium* species (Pitt, 2000). It was registered as
LMSA 1.01.421 in the fungi collection ‘’Souchothèque de Bretagne’’ (Brest, France). For
conservation and collecting spores, the fungus was cultured in tubes of potato dextrose agar
medium (PDA) (Difco Laboratories) at 25°C. Spores were collected from seven-day-old
mycelium by adding, 2 mL of sterile water containing 0,01 % of Tween 80 (Sigma-Aldrich)
to each tube and flooding with agitation at 250 rpm. Mycelia were cultured from a suspension
of 5·10⁶ spores of *P. glabrum* inoculated in 250 ml Erlenmeyer flasks containing 50 ml of
potato dextrose browth medium (PDB) (Difco Laboratories). For suppression subtractive
hybridization, *P. glabrum* was exposed briefly to two different conditions: stress (40°C, 120
min) and control (25 °C), hereafter referred as S and C, respectively. In both case, mycelia
were first grown in 50 ml of PDB medium shaken at 25°C, 120 rpm for 48 h, then sterile-
filtered, transferred into 50 ml of sterile PDB media and grown for 120 min at 25°C or 40°C
shaking at 120 rpm. Mycelia were then aseptically harvested by filtration, washed thoroughly
with sterile water, quickly frozen in liquid nitrogen and stored at –80°C. For microarray
experiments, five different thermal stress conditions were studied, as defined in a reduced
central composite design (CCD) (Box et al., 1978): S1 (33 °C, 42 min), S2 (33 °C, 198 min),
S3 (40 °C, 120 min), S4 (47°C, 42 min) and S5 (47°C, 198 min) (Fig. 1). These five culture
conditions were defined from a previous CCD used for physiological investigations of
thermal stress in *P.glabrum* (unpublished data). For each condition, mycelia were grown 48 h
in triplicate (a, b and c), heat shocked according to the different conditions and harvested
following the method previously described above. Mycelia referred as ‘‘Reference’’(R) and
‘‘Control’’ (C) were also grown in triplicate, at 25°C, in order to constitute respectively the
unique reference pool (pool R) or the Ca, Cb, Cc samples used for statistical comparison with
results obtained under stress conditions.

Suppressive Subtractive Hybridization

For each sample, the frozen mycelium was reduced to a fine powder by grinding in a
mortar containing liquid nitrogen and total RNA was extracted with Trizol (Invitrogen,
Carlsbad, CA, USA) according to the manufacturer's instructions. Total RNA was quantified
using Nanodrop spectrophotometer (Labtech) and their integrity was verified with
Bioanalyser (Agilent Technologies) and RNA Nano LabChip (Agilent Technologies). Poly
(A)⁺ RNAs were then isolated with PolyATtract kit (Promega) according to the
manufacturer's instructions. The Poly (A)⁺ RNAs were then precipitated overnight at -20°C
with 0.1 volume of 3 M sodium acetate (pH 5.2) and 1 volume of isopropanol. Poly (A)⁺
RNAs were then washed twice in 75 % (v/v) ethanol and finally re-suspended with 12µL of
RNase-free water. The Poly (A)⁺ RNAs extracted were then quantified using Nanodrop
spectrophotometer (Labtech).

Two reciprocal subtractive cDNA libraries were constructed by SSH on poly (A)⁺ RNAs
corresponding to C (25°C) and S (40°C, 120 min) conditions using the PCR-Select cDNA
subtraction kit (BD Biosciences-Clontech) according to the manufacturer's instructions. After
subtraction, isolated cDNA fragments were amplified by PCR with Advantage cDNA PCR kit
(BD Biosciences-Clontech), purified with QIAquick PCR purification kit (Qiagen) and
cloned into pGEM-T easy vector (Promega) in JM 109 Bacteria (Promega) according to the
manufacturer's instructions. Bacterial clones of both forward and reverse libraries were
cultured for 24 h at 37°C, in Luria-Bertani (LB) medium (supplemented with 100 mg
ampicillin 1'). Plasmids were then extracted from these bacterial clones using an alkaline lysis
plasmid minipreparation (Birnboim & Doly, 1979). Positive clones containing a single insert
were then selected after EcoRI digestion and 1% agarose gel electrophoresis.

DNA manipulations and sequence analysis

The cDNAs isolated from positive clones, were sequenced with BigDye Terminator
chemistry on an AB 3130 xl sequencer (Applied Biosystems). All sequencing was performed
with the Sp6 primer (5'-GATTTAGGTGACACTATAG-3', Tm = 56°C) at a concentration of
5 µM for each reaction. The obtained raw sequences were trimmed using Phred software
(Ewing et al., 1998), plasmid sequences were eliminated with SeqClean software (TIGR
software tools, http://www.tigr.org/tdb/tgi/software). The cleaned sequences were then
compared using BLASTx program (Altschul et al., 1990) with non-redundant GenBank and
Swiss-Prot protein sequence databases. Sequences of cDNA contained in both subtracted
libraries were annotated and the corresponding ESTs were classified based on the function of
their putative encoded proteins. Functional classification used in this work comes from a
similar study on *Paracoccidioides brasiliensis* (Felipe et al., 2005) and is originally based on
MIPS clustering (Ruepp et al., 2004).
## cDNA microarray construction

In order to normalise the amount of cDNA used for microarray, the bacterial clones used for the subtracted libraries, were cultured in LB medium (with 100 mg ampicillin $1^{-1}$) for 6 h at 37°C, at 140 rpm. Twenty microliters of each culture were diluted in 80 µL of sterile water, then incubated at 95°C for 8 min and 4°C for 10 min. Two microliters of each lysed bacterial clone solution, were used in duplicate as cDNA template for PCR amplification in 50 µL of final reaction volume with an initial denaturation at 94 °C, 3 min, 40 cycles of 94°C for 20 s, 56°C for 30 s, 68°C for 2 min and a final elongation step at 68°C for 5 min. The reaction was performed with Taq DNA polymerase Kit (Amersham-Pharmacia Biotech): 3 U Taq polymerase per reaction, 0.25 mM dNTP, 0.01 µM T7 primer (5'-TAATACGACTCACTATAGGG-3', Tm = 56°C) and 0.01 µM Sp6 primer (5'-GATTTAGGTGACACTATAG-3', Tm = 56°C). PCR product duplicates were then pooled. Amplified cDNAs were purified using Multiscreen plate MSNU 03010 (Millipore) and resuspended in 3X saline sodium citrate (SSC) (Invitrogen) containing 1.5 M betaine (Sigma-Aldrich). A semi-quantitative estimation was then performed on agarose gel electrophoresis with Smart Ladder SF (Eurogentec) in order to validate for each clone, the presence of single amplified cDNA band and its normalised intensity in each clone. Purified PCR products, were then spotted in triplicate onto GAPS II coated slides (Corning) with Eurogrid Spotter (Eurogentec).

## Target preparation and hybridization

For target preparation, total RNAs were extracted from biological triplicates (a, b, c) of each mycelium sample grown in control (C and R) and stress conditions (S1, S2, S3, S4, S5),
as previously described. Total RNA from each sample was extracted with Trizol (Invitrogen) and genomic DNA leftovers were eliminated by DNase I treatment with RNase-free DNase set (Qiagen) and RNeasy MinElute Cleanup kit (Qiagen) according to the manufacturer's instructions. Purified total RNAs were then analysed with Nanodrop spectrophotometer (Labtech) and Bioanalyser (Agilent Technologies). Total RNAs extracted from the 3 distinct samples Ra, Rb, Rc were then pooled in order to form a single reference pool (pool R). This pooled sample was used as a unique reference for hybridizations against samples corresponding to different conditions (C, S1, S2, S3, S4, S5). Five hundred nanograms of each RNA sample were used to prepare cDNAs by Reverse Transcription (RT). Those cDNAs were then 'in vitro' transcribed and labelled in amplified RNA (aRNA) targets with Amino-Allyl MessageAmp II kit (Ambion) and Cy3 / Cy5 dye (Amersham-Pharmacia Biotech), according to the manufacturer's instructions.

Before target hybridization, the cDNA microarray slides were first hydrated (10 s), dried (10 s) and fixed for 20 s at 60 000 µJoules with Spectrolinker UV Crosslinker (Spectroline). Microarrays were then agitated 15 min at 100 rpm in a solution containing 130 mL 1-methyl-2-pyrrolidone (Sigma-Aldrich), 0.22 M succinic anhydride (Sigma-Aldrich) and 10 mL 1 M borate buffer, pH 8. The arrayed cDNA fragments were then denatured by placing the slides in milliQ water at 95°C for 2 min, washed in 95% (v/v) ethanol (Prolabo) and immediately dried by centrifugation at 700 rpm for 3 min at room temperature. Slides were then pre-hybridized by incubating at 42°C for 1 h in a solution containing 3X SSC, 0.3% SDS, 1% Bovine Serum Albumine (BSA) sterilized by 0.22 µm filtration. Microarrays were then washed in 5 different milliQ water baths and then immediately dried by centrifugation at 700 rpm for 3 min at room temperature.
Competitive hybridizations were then realised by first mixing 500 ng of Cy3-labeled aRNA targets (S) and 500 ng of Cy5-labeled aRNA targets (C), denaturing 2 min at 100°C and incubating 30 min at 37°C. Mixed targets were then transferred onto the microarray surface, a cover slip was applied on top and microarrays were then incubated for 12 h at 42°C. After hybridization, microarrays were rinsed for 2 min at 100 rpm in four successive baths (respectively 2 X SSC, 0.1 % SDS; 1 X SSC, 0.2 X SSC twice) and immediately dried by centrifugation at 700 rpm for 3 min at 25°C.

## Microarray data analysis

Hybridized slides were scanned at 10-µm resolution with ScanArray 4000 (Perkin Elmer) and ScanArray Express software (Perkin-Elmer). Pixel fluorescence intensity and differential expression ratios were determined from the picture files (.tiff) for both channel (Cy3 and Cy5) by using GenePix Pro 6.0 software (Molecular Devices Corporation). Lowess fitness regression (Yang *et al.*, 2002) was applied for global normalisation of raw expression ratios. Normalised data for the 913 cDNAs were log₂ transformed and where an EST was represented by more than one cDNA, the median value was recorded. A statistical analysis of the results obtained for each experimental condition, was carried out with Microsoft Excel. A global ANOVA on results obtained for C, S1 and S2 conditions followed by a LSD test, were applied in order to isolate ESTs which were, significantly (*p*≤0.05) differentially expressed between S1-C and S2-C. Among these ESTs, only those with S:C expression ratios, corresponding with differences of at least 0.7 log₂ (i.e. 1.62 fold change) were selected to make the up- and down-regulated EST lists.
qRT-PCR experiments

Reverse transcription (RT) was performed on 50 µg of poly (A)⁺ RNA purified from total RNA samples used for preparing microarray labeled targets. ImProm-II Reverse Transcription System (Promega) was used according to the manufacturer's instructions, Real-time PCR (qPCR) experiments were performed with a Mini-Opticon thermocycler (Bio-Rad). Experiments were performed on 1 µl of two twice-diluted cDNA samples with 0.5 M each primer and 7.5 µl 2X iQ SYBR Green Supermix (Bio-Rad) in a final reaction volume of 15 µl. Cycling conditions were 2 min at 50°C, 2 min at 95 °C, 40 cycles (15 s at 95 °C, 1 min at 60 °C) and 15 s at 95 °C.

In the first step of the qRT-PCR experiments, the internal normalisation control was determined. A set of eight genes was chosen and composed of β-actin (frequently used as a housekeeping gene) and seven other genes presenting relative thermal stability according to microarray results. This group contained seven heterogeneous genes encoding adenosine kinase, 14-3-3 protein homolog, FAD-dependent sulfhydryl oxidase, copper-fist DNA binding domain protein, ATP synthase gamma chain, pre-mRNA splicing factor ini1 and S-phase kinase-associated protein 1A. The sequences of these genes *in r. P. glabrum* were unavailable in public databases. Nevertheless, for β-actin gene, partial sequences were obtained from a previous study. Moreover, the seven other genes were previously cloned and partially sequenced for SSH experiments.

Another set of ten ESTs was also selected from the microarray analysis in order to assess their expression by qRT-PCR and validate our results. These ESTs were selected because their expression appeared to vary significantly with thermal stress. These ESTs encode for
heat shock protein 70, heat shock protein 98, trehalose-6-phosphate synthase, import
translocase tim 50, NADP alcohol dehydrogenase, DMP phosphatase, hypothetical protein
AN5480.2, hypothetical protein EAL 90164.1, AAA family ATPase, transcriptional activator
HAA1.

For each gene, a pair of specific primers was designed using Primer 3 software
(http://frodo.wi.mit.edu/cgi-bin/primer3/primer3.cgi) (Rozen & Skaletsky, 2000) in order to
produce 120-180 bp amplicons (Table 1a and 1b).

Each primer pair was first tested with qRT-PCR in order to check that its PCR
amplification efficiency was superior than 0.90. For each gene, the melting curve profile was
first inspected in order to verify amplification specificity. Standard curves were then
determined by plotting the threshold cycle (Ct) values of serial 10-fold dilutions (5, 50, 500,
5000) of a calibrator cDNA template as a linear function of the log of the dilution factor. The
calibrator cDNA was obtained by reverse transcription on a mix of different stressed or
control RNAs. Expression values were then calculated with the threshold cycle (Ct) of the
fluorescence measures obtained. Normalised data were calculated with the 2<sup>-ΔΔCt</sup> method
(Livak & Schmittgen, 2001) using the expression of an internal control gene. Microarray and
qRT-PCR results obtained for the ten selected genes, were then compared using a Bland-
Altman plot (Bland & Altman, 1996) as implimented by others (Steenman *et al.*, 2005).
Subtracted cDNA libraries

More than 1300 clones were randomly selected from the forward and reverse cDNA
subtracted libraries constructed for Penicillium glabrum after comparison between C (25 °C)
and S (40 °C, 120 min). Among them, 1056 positive clones were sequenced and finally 913
clean and useful cDNA sequences were obtained. Using BLASTx on GenBank and Swiss-
Prot protein databases, 80 % of these sequences matched with ORFs of other fungi:
Aspergillus nidulans, Aspergillus fumigatus, Emericella nidulans, Neurospora crassa,
Saccharomyces cerevisiae, Schizosaccharomyces pombe (Fig 2a). Analysis of these data
allowed isolation of 445 non-redundant Expressed Sequence Tags (ESTs) with E value lower
than 10⁻⁵. Only eight of the 445 ESTs were represented by more than eight different cDNA
clones. The two reciprocal subtracted libraries contained 197 and 275 potentially up- and
down-regulated ESTs, respectively; 27 of them were shared by both libraries. The results
obtained were available from our website at http://pagesperso.univ-
brest.fr/~vasseur/Nevarez_et_al_2008/. The 197 potentially up-regulated ESTs were
composed by 108 identified ESTs, 48 ESTs with hypothetical functions and 41 ESTs that
appeared unknown in the databases. From the 275 ESTs contained in the down-regulated
library, 195 matched with annotated sequences, 47 had hypothetical functions and 33
appeared to be unknown. Of all the isolated ESTs, 64 % showed similarity with annotated
ORFs, 20 % encoded for hypothetical proteins and 16 % did not match any sequence in the
public databases (Fig. 2b).
The analysis of both subtracted libraries showed a wide variety of ESTs encoding for proteins of different functional classes (Fig. 2c and 2d). Several ESTs coding for proteins often associated with cellular stress (heat shock proteins, trehalose-6-phosphate synthase, superoxide dismutase) appeared up-regulated after thermal stress (Fig. 2c). ESTs encoding other categories of proteins (metabolism, protein fate, transport proteins, etc.) showed relatively similar representation in both libraries. Some ESTs appeared down-regulated in *P. glabrum* after heat shock (Fig. 2d), as those encoding proteins involved in transcription (RNA polymerase, initiation factor, splicing factors, etc.) or implicated in protein synthesis (ribosomal proteins, translation initiation factors, elongation factors, etc.). Other ESTs isolated in each library coded for hypothetical or unknown proteins

## cDNA microarray construction and RNA target preparation

The differential expression of ESTs isolated from SSH was analysed using custom-made cDNA microarrays. Amplified cDNAs isolated from up- and down-regulated subtracted libraries were spotted on cDNA microarrays. Each slide contained, in triplicate, the 913 subtracted cDNA corresponding to 445 different ESTs isolated from SSH. Many ESTs were represented by one to four different cDNA sequences and some were even represented by up to eight cDNAs. This redundancy was useful for evaluating the reproducibility of the expression results obtained for each EST.

Before target preparation, the quality of total RNAs extracted was first verified. Analysis of total RNA triplicates from C conditions revealed a convenient profile with three different and clearly defined peaks, corresponding to migration marker, 18S and 28S rRNAs. Similar results were obtained for total RNA triplicates extracted from S1, S2, S3 and S4 (data not shown). Analysis of total RNA samples from the more drastic S5 condition (47°C, 42
min) showed a migration marker peak and a large amount of very small RNA fragments indicating substantial total RNA degradation. We thus excluded the S5 stress condition from the central composite design.

## Determination of housekeeping genes for validation of gene expression using RT-q-PCR

The final validation of the microarray results by qRT-PCR, requires previous selection of internal reference genes for qRT-PCR normalisation. Expression of the eight potential housekeeping genes selected (β-actin, adenosine kinase, 14-3-3 protein homolog, FAD dependent sulfhydryl oxidase, copper-fist DNA binding domain protein, ATP synthase gamma chain, Pre-mRNA splicing factor ini1 and S-phase kinase-associated protein 1A) was assessed under the different experimental conditions (C, S1, S2, S3, S4) with qRT-PCR. For each primer pair tested, PCR specificity was confirmed by melting curve analysis (data not shown). Expression results obtained for the eight genes tested were comparable, where Ct values were stable under C, S1, S2 conditions but substantially increased under S3 and S4 conditions (Fig. 3). Given these results, normalisation of qRT-PCR results could not be applied to S3 and S4 conditions since stable expression of internal control genes is required for each tested condition. As microarray data could not be validated for these thermal stress conditions, analysis was only done on S1 and S2 conditions. Beta-actin was selected as housekeeping gene for qRT-PCR normalisation, owing to its stable expression under C, S1, S2 conditions and its frequent use in the literature. ATP synthase gamma chain and adenosine kinase were also useful genes and were used to further confirm qRT-PCR normalisation (data not shown).
Microarray results for S1 and S2 thermal stress conditions

Based on the 913 cDNA sequences spotted in triplicate representing 445 different ESTs, the statistical analysis of results with a global ANOVA and LSD test indicated 35 and 49 significantly up- and down-regulated ESTs. Their S-to-C expression ratios corresponded to at least a +/-.7 log2 (1.62-fold change). Expression of most of the isolated ESTs represented a ca. +/-2-fold change, but a few ESTs showed up to +3.5- and even a -5.0-fold change for up- and down-regulated ESTs, respectively.

Analysis of the microarray results revealed 49 highly significant down-regulated ESTs
under thermal stress conditions S1 and S2 (Table 2). From this set, 26 ESTs encoded
ribosomal proteins and accounted for 53% of the results. This large group also contained
ESTs encoding proteins implicated in transcription (bZip transcription factor CpscA,
transcriptional activator HAA1, histone H4, etc.), RNA splicing (U6 small nuclear RNA
associated protein, mago nashi protein homolog) and protein synthesis (translation elongation
factors G2). A few proteases (glutamate carboxypeptidase-like protein, candidapepsin 2
precursor (aspartate protease)) and proteins involved in general metabolism (xanthine
phophoribosyltransferase 1 and ornithine decarboxylase) were also observed. Few down-
regulated ESTs also encoded miscellaneous proteins and few hypothetical or unknown
proteins.

From the microarray results, 35 highly significant up-regulated ESTs were identified
for thermal stress conditions S1 and S2 (Table 3). From this set, six and four ESTs encoded
heat shock proteins (HSP 30, HSP 60, HSP 70, mitochondrial HSP 70, HSP 90, HSP 98/104)
and proteins with antioxidant functions (superoxide dismutase, glutathione S-transferase 3,
cytochrome c oxidase polypeptide VIa, flavohemoprotein), respectively. The up-regulated EST set also included seven other ESTs coding for metabolism proteins (L-xylulose reductase, mannitol-1-phosphate-5-dehydrogenase, 2,3-diketo-5-methylthio-1-phosphopentane phosphatase, carboxylesterase, NADP-dependent alcohol dehydrogenase, trehalose-6-phosphate synthase). Up-regulated ESTs were also completed by few ESTs encoding proteases (aspergillopepsin A and vacuolar protease A) and other miscellaneous proteins. About 37% of identified up-regulated ESTs coded for hypothetical or unknown proteins.

## Validation of microarray results with RT-q-PCR

The microarray results obtained for S1 and S2 conditions were validated by assessing the expression level of a set of 10 differentially expressed ESTs selected from the microarray analysis using qRT-PCR. Relative expression ratios of the 10 selected ESTs were normalised against β-actin gene expression and results obtained were comparable to results from the microarray (Table 4). The Bland-Altman plot (Fig. 4) carried out with values obtained for S1 and S2 conditions, showed that all of the log₂-transformed microarray data and qRT-PCR data lay within the ‘95 % confidence intervals’, delimited by a difference of 1. Moreover, more than 65% of the comparisons differed by less than 0.5. This similarity between both sets of data validate the microarray results for S1 and S2 conditions. Reliability of the data was further verified as similar results were also obtained by normalising qRT-PCR results with expression of the two other housekeeping genes selected: adenosine kinase and ATP synthase (data not shown).
In the first part of the study, SSH was applied on *Penicillium glabrum* to isolate 913 cDNA sequences corresponding to 445 differentially expressed ESTs in response to thermal stress (i.e. 40°C, 120 min). SSH is frequently employed to isolate differentially expressed genes in response to varying nutritional or environmental conditions. SSH overcomes the problem of transcript abundance differences in the mRNA population by incorporating a normalisation step which equalises the relative quantity of the cDNA. It enhances the probability of identifying the increased expression of low-abundance transcripts (Diatchenko *et al.*, 1996; Gurskaya *et al.*, 1996).

Many of the isolated ESTs encoded proteins with a wide variety of functions while some of them corresponded to hypothetical or unknown proteins. The large number of ESTs is remarkable compared to results obtained from similar studies on other fungi. For example, SSH was used to isolate specific transcripts in either yeast or mycelial forms of *Ophiostoma piceae* or *P. marneffei* and only revealed 50 and 43 genes, respectively (Dogra & Breuil, 2004; Liu *et al.*, 2007). Other studies conducted on the effects of different nutritional growth conditions, isolated 37 genes implicated in growth with penicillin-repressing or non-repressing carbon sources in *P. chrysogenum* (Castillo *et al.*, 2006). The large number of ESTs isolated here in *P. glabrum* under heat shock conditions (i.e. 40°C, 120 min) and their wide functional classification, could be explained by the fact that temperature is involved in general cellular mechanisms and influences many functional categories. It has been shown that two-thirds of the *Saccharomyces cerevisiae* genome is involved in response to environmental changes such as temperature, oxidation, nutriments, pH and osmolarity (Causton *et al.*, 2001).
Expression of ESTs contained in both subtracted libraries was further assessed, using microarrays to confirm their differential expression under thermal stress conditions. Custom cDNA microarray is a very convenient tool for studying subtracted libraries in details and are more appropriate than Northern blots or qRT-PCR methods for assessing the expression of a large number of ESTs. The combination of SSH and custom cDNA microarrays has been previously employed in *Gibberella zeae* to identify genes involved in sexual reproduction (mediated by ascospores) and implicated in disease occurring in cereal crops (Lee *et al.*, 2006). In our study, the custom microarray used in triplicates of the 913 subtracted cDNA sequences isolated corresponded to 445 different ESTs. Beyond just confirming the subtracted library approach, microarrays were also employed to analyse expression of these ESTs under five different thermal stress conditions (S1, S2, S3, S4 and S5).

Degradation total RNA observed in S5 conditions (47°C, 198min) leads us to eliminate this experiment from the study. This result may be correlated with a decrease of mycelium growth observed under these conditions: fungal dry weight decreased about 35% compared to dry weight in optimal growth conditions (unpublished data). Microarray analysis for the four other stress conditions, required final validation by qRT-PCR and prior selection of an internal reference gene for normalisation. As no housekeeping genes have yet been defined for *P. glabrum* or closely related fungal species, comparative studies and microarrays were used. Other than the seven potential housekeeping genes selected after microarray analysis, the β-actin gene was also chosen among other well-known internal reference genes such as β-tubulin, glyceraldehyde-3-phosphate dehydrogenase and 18S ribosomal RNA (Bustin, 2000; Livak & Schmittgen, 2001). Selection of β-actin was based on its frequent use in qRT-PCR normalisation in *Candida albicans* (Nailis *et al.*, 2006), *Phytophthora infestans* (Armstrong *et al.*, 2005) and *Neurospora crassa* (Mohsenzadeh *et al.*, 1998). The systematic increase in Ct values observed using qRT-PCR for every gene tested in S3 and S4 conditions
compared to C, S1 and S2 conditions in spite of the fact that mRNA quantity and quality was
validated before RT. Variation in cDNA abundance can be partially explained by mRNA
instability after exposure to high temperatures (Bond, 2006). For example, substantial decay
of Hsp70 mRNA was reported in *S. cerevisiae* after a 3-h exposure to 48°C (Kapoor *et al*.,
1995). Messenger RNA decay primarily involves deadenylation or decapping (Tucker &
Parker, 2000; Wilusz *et al*., 2001). These processes probably affect mRNA stability and
reverse transcription (Fleige & Pfaffl, 2006) without any change in the global quantity of
mRNA. Given the increase in Ct values, housekeeping genes were not stable and microarray
analyses were not possible for S3 and S4 conditions. This result was surprising since a
preliminary study in *P. glabrum* only showed a slight decrease in growth under similar
conditions (40°C, 120 min; unpublished data). Analysis of the quality of the total RNAs
extracted from S3 and S4 conditions indicated that RNA integrity appeared conserved under
both conditions. In addition, other thermal stress studies on fungi have been successfully
conducted at temperatures close to 40°C: 37°C for *S. cerevisiae* (Gasch *et al*., 2000; Causton
*et al*., 2001; Sakaki *et al*., 2003), 39° for *Schizosaccharomyces pombe* (Chen *et al*., 2003), 42
°C for *N. crassa* (Mohsenzadeh *et al*., 1998), 45 °C for *C. albicans* (Zeuthen & Howard,
1989) and 50°C for *Aspergillus* (Bai *et al*., 2003). Nevertheless, results obtained for S3 and
S4 conditions clearly showed that defining a standard “housekeeping” gene can be difficult
for experiments conducted under diverse conditions as reported by some authors (Schmittgen
& Zakrajsek, 2000; Yan & Liou, 2006).

Given these results, both S3 and S4 conditions were excluded from the microarray
analysis and results obtained for S1 and S2 were then confirmed using qRT-PCR on 10 genes
as observed in the Bland-Altman plot (see Fig. 4). Statistical treatment of the microarray data,
identified 35 significantly up-regulated ESTs and 49 significantly down-regulated. Some of
the down-regulated ESTs encoded proteases or general metabolism proteins. To our knowledge none of these genes have ever been described as being down-regulated in previous thermal stress studies performed on fungi. Another group of down-regulated ESTs encoded proteins implicated in transcription, RNA splicing and protein synthesis. Several studies carried out on *S. cerevisiae* exposed to thermal stress (25°C to 37°C) have also reported down-regulation of genes coding for histone H4 (Gasch *et al.*, 2000; Causton *et al.*, 2001; Sakaki *et al.*, 2003), translation elongation factor G2 (Sakaki *et al.*, 2003) and U6 small nuclear RNA associated protein (Gasch *et al.*, 2000). The group of ESTs involved in translation included a large set of 27 ESTs encoding ribosomal proteins. Down-regulation of these ESTs and repression of the corresponding protein synthesis have been frequently reported in the literature. In *S. cerevisiae*, it has been shown that ribosomal protein synthesis is repressed after exposure to 33°C for 90 min (Gorenstein & Warner, 1976), seemingly linked to a decrease of their corresponding mRNAs (Hereford & Rosbash, 1977) caused by inhibition of their transcription (Kim & Warner, 1983). Microarray studies at 25° to 37°C in *S. cerevisiae* and at 30° to 39°C in *Sc. pombe* (Gasch *et al.*, 2000; Causton *et al.*, 2001; Chen *et al.*, 2003; Sakaki *et al.*, 2003) also confirmed down-regulation of this genes in response to thermal stress. For example, in *S. cerevisiae* mRNAs have been reported to represent an important part of total poly(A)⁺ mRNAs; repression of these genes may therefore allow energy to be dedicated to increased expression of genes involved in protective response (Jelinsky & Samson, 1999). This last group of ESTs encoding for ribosomal proteins, may contain some potential candidates for molecular markers of thermal stress in *P. glabrum*.

In addition, we identified 35 up-regulated ESTs that encoded proteases and some proteins implicated in metabolism. As for down-regulated ESTs, and to our knowledge, none of these genes have been previously reported in the literature on fungi, except trehalose-6-
phosphate synthase (tpsA), frequently described as being induced under heat shock conditions
in yeasts (Gasch et al., 2000; Causton et al., 2001; Chen et al., 2003; Sakaki et al., 2003).
Trehalose has been described to act as a protectant against various environmental stresses and
has been showed to stabilize cellular structures under stress conditions (Tereshina, 2005).
This disaccharide has been previously described in *S. cerevisae* to play a role in the
acquisition of stress tolerance (De Virgilio et al., 1994; Hottiger et al., 1994; Singer &
Lindquist, 1998). For *P. glabrum*, although it has never been reported in the literature, an
accumulation of trehalose concurs with our preliminary results (unpublished data).

A group of up-regulated ESTs encoding HSPs was also identified in our study. These
highly conserved proteins are certainly the most investigated element of the thermal stress
response and increase in expression of the corresponding genes has been frequently observed
in thermal stress studies (Lindquist, 1986). These proteins are classified into several families
according to their molecular mass and, despite their different roles and proprieties, most HSPs
have been shown to generally act as molecular chaperones. They are mainly involved in
folding, assembly, regulation, degradation of other proteins and play an important role in
thermotolerance (Craig, 1985; Feder & Hofmann, 1999).

Exposure to elevated temperatures has been previously reported to generate reactive
oxygen species (ROS) responsible for oxidative damage in yeasts and filamentous fungi
(Davidson et al., 1996; Noventa-Jordao et al., 1999; Bai et al., 2003). As observed in our
results, antioxidant systems appear to be extremely important for detoxifying cells. The
synthesis of enzymes such as superoxide dismutase, glutathione S transferase, cytochrome c
oxidase VIa, may serve to reduce the level of ROS in cells (Choi et al., 1998; Marques et al.,
2004). In *S. cerevisiae*, it was shown that the synthesis of superoxide dismutase is enhanced in
thermal stress conditions (Piper, 1995) and plays a role in resistance to heat shock from 30°C
to 50°C (Davidson et al., 1996). Microarray studies on the yeasts *S. cerevisiae* and *Sc. pombe*
have also confirmed induction of genes encoding superoxide dismutase, cytosolic glutathione S-transferase and cytochrome c oxidase popypeptide VIb (which has similarity to cytochrome c oxidase popypeptide VIa gene isolated in this study) (Gasch et al., 2000; Causton et al., 2001; Chen et al., 2003). These ESTs coding for antioxidant proteins also appeared to be induced in heat shocked *P. glabrum* and could be potential molecular markers of thermal stress.

Considering these results, the present study provides a preliminary basis for the genetic regulation of thermal stress response in *P. glabrum*. Temperature is one of the most important factors for fungal growth and as others may represent potential sources of stress that can directly affect the physiology of fungi (growth, conidiation, synthesis of secondary metabolites, etc.). The transcriptional analysis performed here combined efficient and complementary molecular techniques making it possible to isolate 84 heat-regulated ESTs having a wide variety of functions and to assess their expression under thermal stress conditions. Some of these ESTs have already been associated with different stress-related responses as they encode heat shock proteins, trehalose-6-phosphate synthase, superoxide dismutase, ribosomal proteins or proteins involved in transcription/translation. In *P. glabrum*, they may be “common environmental stress response” genes and respond in a similar manner to many different environmental changes as those identified in studies on yeasts (Gasch et al., 2000; Chen et al., 2003). On the other hand, a large number of isolated ESTs encoded general metabolism proteins, proteases, miscellaneous proteins and hypothetical or unknown proteins. To our knowledge, a large majority of them have never been associated with fungal stress before and may represent new potential stress markers to characterise the physiological state of *P. glabrum* (stressed or unstressed). Nevertheless as reported by some authors, caution should be taken when analysing microarray results as mRNA stability and decay induced by
stress treatments may vary and have some influence gene expression profiling (Monje-Casas
et al., 2004; Michan et al., 2005). Further investigations may study more precisely the
stability of each mRNA and investigate the response of the different potential markers
isolated.

Acknowledgements

This work benefited from financial support from the Brittany Regional Council and
the invaluable help of the sequencing and microarray facilities of Ouest-Genopole located in
Roscoff and Nantes, respectively. The authors also thank Dr. Jean-Marie Heslan (Inserm,
Nantes) and Dr. David Mazurais (Ifremer, Brest) for their helpful advice and technical support
for qRT-PCR experiments. Authors would also like to thank Emeric Dubois (Inserm, Nantes)
for his assistance in the microarray analysis.
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. (1990). Basic local alignment search tool. *J Mol Biol* **215**, 403-410.

Ancasi, E. G., Carrillo, L., Benitez Ahrendts, M. R. (2006). Moulds and yeasts in bottled water and soft drinks. *Rev Argent Microbiol* **38**, 93-96.

Armstrong, M. R., Whisson, S. C., Pritchard, L., Bos, J. I., Venter, E., Avrova, A. O., Rehmany, A. P., Bohme, U., Brooks, K. & other authors (2005). An ancestral oomycete locus contains late blight avirulence gene Avr3a, encoding a protein that is recognized in the host cytoplasm. *Proc Natl Acad Sci U S A* **102**, 7766-7771.

Ashburner, M., Bonner, J. J. (1979). The induction of gene activity in drosophila by heat shock. *Cell* **17**, 241-254.

Bai, Z., Harvey, L. M., McNeil, B. (2003). Oxidative stress in submerged cultures of fungi. *Crit Rev Biotechnol* **23**, 267-302.

Birnboim, H. C., Doly, J. (1979). A rapid alkaline extraction procedure for screening recombinant plasmid DNA. *Nucleic Acids Res* **7**, 1513-1523.

Bland, J. M., Altman, D. G. (1996). Measurement error. *BMJ* **313**, 744.

Bond, U. (2006). Stressed out! Effects of environmental stress on mRNA metabolism. *FEMS Yeast Res* **6**, 160-170.

Borjesson, T., Stollman, U., Schnurer, J. (1992). Volatile metabolites produced by six fungal species compared with other indicators of fungal growth on cereal grains. *Appl Environ Microbiol* **58**, 2599-2605.

Box, G. E. P., Hunter, W. G., Hunter, J. S. (1978). Statistics for experimenters. An introduction to design data analysis and model building. Wiley, New York.

Bustin, S. A. (2000). Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. *J Mol Endocrinol* **25**, 169-193.

Cabral, D., Fernandez Pinto, V. E. (2002). Fungal spoilage of bottled mineral water. *Int J Food Microbiol* **72**, 73-76.

Castillo, N. I., Fierro, F., Gutierrez, S., Martin, J. F. (2006). Genome-wide analysis of differentially expressed genes from *Penicillium chrysogenum* grown with a repressing or a non-repressing carbon source. *Curr Genet* **49**, 85-96.

Causton, H. C., Ren, B., Koh, S. S., Harbison, C. T., Kanin, E., Jennings, E. G., Lee, T. I., True, H. L., Lander, E. S. & other authors (2001). Remodeling of yeast genome expression in response to environmental changes. *Mol Biol Cell* **12**, 323-337.
Chang, S. C., Wei, Y. H., Wei, D. L., Chen, Y. Y., Jong, S. C. (1991). Factors affecting the production of eremofortin C and PR toxin in *Penicillium roqueforti*. *Appl Environ Microbiol* **57**, 2581-2585.

Chen, D., Toone, W. M., Mata, J., Lyne, R., Burns, G., Kivinen, K., Brazma, A., Jones, N., Bahler, J. (2003). Global transcriptional responses of fission yeast to environmental stress. *Mol Biol Cell* **14**, 214-229.

Choi, J. H., Lou, W., Vancura, A. (1998). A novel membrane-bound glutathione S-transferase functions in the stationary phase of the yeast *Saccharomyces cerevisiae*. *J Biol Chem* **273**, 29915-29922.

Cosgrove, J. W., Brown, I. R. (1983). Heat shock protein in mammalian brain and other organs after a physiologically relevant increase in body temperature induced by D-lysergic acid diethylamide. *Proc Natl Acad Sci U S A* **80**, 569-573.

Craig, E. A. (1985). The heat shock response. *CRC Crit Rev Biochem* **18**, 239-280.

Dantigny, P., Guilmart, A., Bensoussan, M. (2005). Basis of predictive mycology. *Int J Food Microbiol* **100**, 187-196.

Davidson, J. F., Whyte, B., Bissinger, P. H., Schiestl, R. H. (1996). Oxidative stress is involved in heat-induced cell death in *Saccharomyces cerevisiae*. *Proc Natl Acad Sci U S A* **93**, 5116-21.

De Virgilio, C., Hottiger, T., Dominguez, J., Boller, T., Wiemken, A. (1994). The role of trehalose synthesis for the acquisition of thermotolerance in yeast. I. Genetic evidence that trehalose is a thermoprotectant. *Eur J Biochem* **219**, 179-86.

Diatchenko, L., Lau, Y. F., Campbell, A. P., Chenchik, A., Moqadam, F., Huang, B., Lukyanov, S., Lukyanov, K., Gurskaya, N. & other authors (1996). Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. *Proc Natl Acad Sci U S A* **93**, 6025-30.

Dogra, N., Breuil, C. (2004). Suppressive subtractive hybridization and differential screening identified genes differentially expressed in yeast and mycelial forms of *Ophiostoma piceae*. *FEMS Microbiol Lett* **238**, 175-81.

Ewing, B., Hillier, L., Wendl, M. C., Green, P. (1998). Base-calling of automated sequencer traces using phred. I. Accuracy assessment. *Genome Res* **8**, 175-185.

Feder, M. E., Hofmann, G. E. (1999). Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. *Annu Rev Physiol* **61**, 243-82.

Felipe, M. S., Andrade, R. V., Arraes, F. B., Nicola, A. M., Maranhao, A. Q., Torres, F. A., Silva-Pereira, I., Pocas-Fonseca, M. J., Campos, E. G. & other authors (2005). Transcriptional profiles of the human pathogenic fungus *Paracoccidioides brasiliensis* in mycelium and yeast cells. *J Biol Chem* **280**, 24706-24714.

Fiedurek, J. (1997). Enhancement of â-galactosidase production and secretion by high osmotic stress in *Penicillium notatum*. *Microbiol Res* 65-69.
**Fleige, S., Pfaffl, M. W. (2006).** RNA integrity and the effect on the real-time qRT-PCR performance. *Mol Aspects Med* **27**, 126-139.

**Gasch, A. P., Spellman, P. T., Kao, C. M., Carmel-Harel, O., Eisen, M. B., Storz, G., Botstein, D., Brown, P. O. (2000).** Genomic expression programs in the response of yeast cells to environmental changes. *Mol Biol Cell* **11**, 4241-4257.

**Gorenstein, C., Warner, J. R. (1976).** Coordinate regulation of the synthesis of eukaryotic ribosomal proteins. *Proc Natl Acad Sci U S A* **73**, 1547-1551.

**Gurskaya, N. G., Diatchenko, L., Chenchik, A., Siebert, P. D., Khaspekov, G. L., Lukyanov, K. A., Vagner, L. L., Ermolaeva, O. D., Lukyanov, S. A. & other authors (1996).** Equalizing cDNA subtraction based on selective suppression of polymerase chain reaction: cloning of Jurkat cell transcripts induced by phytohemaglutinin and phorbol 12-myristate 13-acetate. *Anal Biochem* **240**, 90-7.

**Hereford, L. M., Rosbash, M. (1977).** Regulation of a set of abundant mRNA sequences. *Cell* **10**, 463-467.

**Hottiger, T., De Virgilio, C., Hall, M. N., Boller, T., Wiemken, A. (1994).** The role of trehalose synthesis for the acquisition of thermotolerance in yeast. II. Physiological concentrations of trehalose increase the thermal stability of proteins in vitro. *Eur J Biochem* **219**, 187-193.

**Jelinsky, S. A., Samson, L. D. (1999).** Global response of *Saccharomyces cerevisiae* to an alkylating agent. *Proc Natl Acad Sci U S A* **96**, 1486-1491.

**Kapoor, M., Curle, C. A., Runham, C. (1995).** The hsp70 gene family of Neurospora crassa: cloning, sequence analysis, expression, and genetic mapping of the major stress-inducible member. *J Bacteriol* **177**, 212-221.

**Kim, C. H., Warner, J. R. (1983).** Mild temperature shock alters the transcription of a discrete class of *Saccharomyces cerevisiae* genes. *Mol Cell Biol* **3**, 457-465.

**Kokkonen, M., Jestoi, M., Rizzo, A. (2005).** The effect of substrate on mycotoxin production of selected *Penicillium* strains. *Int J Food Microbiol* **99**, 207-214.

**Kurata, H., Sakabe, F., Udagawa, S., Ichinoe, M., Suzuki, M. (1968).** [A mycological examination for the presence of mycotoxin-producers on the 1954-1967's stored rice grains]. *Eisei Shikenjo Hokoku* **86**, 183-188.

**Lee, S. H., Lee, S., Choi, D., Lee, Y. W., Yun, S. H. (2006).** Identification of the down-regulated genes in a mat1-2-deleted strain of *Gibberella zeae*, using cDNA subtraction and microarray analysis. *Fungal Genet Biol* **43**, 295-310.

**Lindquist, S. (1986).** The heat-shock response. *Annu Rev Biochem* **55**, 1151-91.

**Liu, H., Xi, L., Zhang, J., Li, X., Liu, X., Lu, C., Sun, J. (2007).** Identifying differentially expressed genes in the dimorphic fungus *Penicillium marneffei* by suppression subtractive hybridization. *FEMS Microbiol Lett* **270**, 97-103.
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. *Methods* **25**, 402-408.

Magan, N., Hope, R., Colleate, A., Baxter, E. S. (2002). Relationship between growth and mycotoxin production by *Fusarium* species, biocides and environment. *Eur J Plant Pathol* **108**, 685-690.

Marques, E. R., Ferreira, M. E., Drummond, R. D., Felix, J. M., Menossi, M., Savoldi, M., Travassos, L. R., Puccia, R., Batista, W. L. & other authors (2004). Identification of genes preferentially expressed in the pathogenic yeast phase of *Paracoccidioides brasiliensis*, using suppression subtraction hybridization and differential macroarray analysis. *Mol Genet Genomics* **271**, 667-677.

McAlister, L., Finkelstein, D. B. (1980). Heat shock proteins and thermal resistance in yeast. *Biochem Biophys Res Commun* **93**, 819-824.

Michan, C., Monje-Casas, F., Pueyo, C. (2005). Transcript copy number of genes for DNA repair and translesion synthesis in yeast: contribution of transcription rate and mRNA stability to the steady-state level of each mRNA along with growth in glucose-fermentative medium. *DNA Repair* (Amst) **4**, 469-478.

Mislivec, P. B., Tuite, J. (1970). Species of *Penicillium* occurring in freshly-harvested and in stored dent corn kernels. *Mycologia* **62**, 67-74.

Mohsenzadeh, S., Saupe-Thies, W., Steier, G., Schroeder, T., Fracella, F., Ruoff, P., Rensing, L. (1998). Temperature adaptation of house keeping and heat shock gene expression in *Neurospora crassa*. *Fungal Genet Biol* **25**, 31-43.

Monje-Casas, F., Michan, C., Pueyo, C. (2004). Absolute transcript levels of thioredoxin- and glutathione-dependent redox systems in *Saccharomyces cerevisiae*: response to stress and modulation with growth. *Biochem J* **383**, 139-147.

Nailis, H., Coenye, T., Van Nieuwerburgh, F., Deforce, D., Nelis, H. J. (2006). Development and evaluation of different normalization strategies for gene expression studies in *Candida albicans* biofilms by real-time PCR. *BMC Mol Biol* **7**, 25.

Newbury, J., Peberdy, J. F. (1996). Characterization of the heat shock response in protoplasts of *Aspergillus nidulans*. *Mycol Res* **100**, 1325-1332.

Northolt, M. D., van Egmond, H. P., Soentoro, P., Deijll, E. (1980). Fungal growth and the presence of sterigmatocystin in hard cheese. *J Assoc Off Anal Chem* **63**, 115-119.

Noventa-Jordao, M. A., Couto, R. M., Goldman, M. H., Aguirre, J., Iyer, S., Caplan, A., Terenzi, H. F., Goldman, G. H. (1999). Catalase activity is necessary for heat-shock recovery in *Aspergillus nidulans* germlings. *Microbiology* **145** (Pt 11), 3229-34.

Osherov, N., Mathew, J., Romans, A., May, G. S. (2002). Identification of conidial-enriched transcripts in *Aspergillus nidulans* using suppression subtractive hybridization. *Fungal Genet Biol* **37**, 197-204.

Overy, D. P., Seifert, K. A., Savard, M. A., Frisvad, J. C. (2003). Spoilage fungi and their mycotoxins in commercially marketed chestnuts. *Int J Food Microbiol* **69-77**.
**Piper, P. W. (1993).** Molecular events associated with acquisition of heat tolerance by the yeast *Saccharomyces cerevisiae*. *FEMS Microbiol Rev* **11**, 339-355.

**Piper, P. W. (1995).** The heat shock and ethanol stress response of yeast exhibit extensive similarity and functional overlap. *FEMS Microbiol Lett* **134**, 121-127.

**Pitt, J. I. (2000).** A laboratory guide to common *Penicillium* species, Third edition. Food Science Australia.

**Pitt, J. I., Hocking, A. D. (1997).** Fungi and food spoilage. Blackie Academic and Professional, London.

**Rozen, S., Skaletsky, H. (2000).** Primer3 on the WWW for general users and for biologist programmers. *Methods Mol Biol* **132**, 365-386.

**Ruepp, A., Zollner, A., Maier, D., Albermann, K., Hani, J., Mokrejs, M., Tetko, I., Guldener, U., Mannhaupt, G. & other authors (2004).** The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. *Nucleic Acids Res* **32**, 5539-5545.

**Sakaki, K., Tashiro, K., Kuhara, S., Mihara, K. (2003).** Response of genes associated with mitochondrial function to mild heat stress in yeast *Saccharomyces cerevisiae*. *J Biochem* (*Tokyo*) **134**, 373-384.

**Samson, R. A., Hoekstra, E. S., Frisvad, J. C., Filtenborg, O. (2004).** Introdustion to food and airborne Fungi. Centraalbureau voor Schimmelcultures (CBS), Utrecht, Netherlands.

**Schmittgen, T. D., Zakrajsek, B. A. (2000).** Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. *J Biochem Biophys Methods* **46**, 69-81.

**Shalon, D., Smith, S. J., Brown, P. O. (1996).** A DNA microarray system for analyzing complex DNA samples using two-color fluorescent probe hybridization. *Genome Res* **6**, 639-645.

**Singer, M. A., Lindquist, S. (1998).** Thermotolerance in *Saccharomyces cerevisiae*: the Yin and Yang of trehalose. *Trends Biotechnol* **16**, 460-468.

**Steenman, M., Lamirault, G., Le Meur, N., Le Cunff, M., Escande, D., Leger, J. J. (2005).** Distinct molecular portraits of human failing hearts identified by dedicated cDNA microarrays. *Eur J Heart Fail* **7**, 157-165.

**Tereshina, V. M. (2005).** Thermotolerance in fungi: the role of heat shock proteins and trehalose. *Microbiology* **74**, 293-304.

**Tucker, M., Parker, R. (2000).** Mechanisms and control of mRNA decapping in *Saccharomyces cerevisiae*. *Annu Rev Biochem* **69**, 571-595.

**Udagawa, S. I., Kobatake, M., Kurata, H. (1977).** [Re--estimation of preservation effectiveness of potassium sorbate (food additive) in jams and marmalade (author's transl)]. *Eisei Shikenjo Hokoku* **88-92**.
Wilusz, C. J., Gao, M., Jones, C. L., Wilusz, J., Peltz, S. W. (2001). Poly(A)-binding proteins regulate both mRNA deadenylation and decapping in yeast cytoplasmic extracts. *RNA* 7, 1416-1424.

Yamamori, T., Yura, T. (1982). Genetic control of heat-shock protein synthesis and its bearing on growth and thermal resistance in *Escherichia coli* K-12. *Proc Natl Acad Sci U S A* 79, 860-864.

Yan, H. Z., Liou, R. F. (2006). Selection of internal control genes for real-time quantitative RT-PCR assays in the oomycete plant pathogen *Phytophthora parasitica*. *Fungal Genet Biol* 43, 430-438.

Yang, Y. H., Dudoit, S., Luu, P., Lin, D. M., Peng, V., Ngai, J., Speed, T. P. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. *Nucleic Acids Res* 30, 15.

Zeuthen, M. L., Howard, D. H. (1989). Thermotolerance and the heat-shock response in *Candida albicans*. *J Gen Microbiol* 35, 2509-2518.
Tables :

Table 1:

(a) Primer sequences for the set of eight potential housekeeping genes. (b) Primer sequences for the set of 10 genes selected for validating the microarray results with qRT-PCR

a

<table><thead><tr><th>EST</th><th>Primer</th><th>Primer sequence 5'-&gt;3'</th></tr></thead><tbody><tr><td rowspan="2">β-actin</td><td>β-actin-F</td><td>ACGTTGTCCCCATCTACGAA</td></tr><tr><td>β-actin-R</td><td>GCTCAGCCAGGATCTTCATC</td></tr><tr><td rowspan="2">Adenosine Kinase</td><td>AK-F</td><td>CTGGAGGACTCCGTTGACAT</td></tr><tr><td>AK-R</td><td>TCCCAGAGTGGAAAATCCTG</td></tr><tr><td rowspan="2">14-3-3 protein homolog</td><td>14-3-3-F</td><td>CCGAGCGCTATGAAGAGATG</td></tr><tr><td>14-3-3-R</td><td>ACTCCTCCTTCTGCTCGATG</td></tr><tr><td rowspan="2">FAD dependent sulfhydryl oxidase</td><td>FAD-F</td><td>AACCTCATTATGCGCCTCAC</td></tr><tr><td>FAD-R</td><td>GGGCTTCCTGAAGCTTTTCT</td></tr><tr><td rowspan="2">Copper fist DNA binding domain protein</td><td>Copper fist-F</td><td>AGGTAGTCCGACCAATGACG</td></tr><tr><td>Copper fist-R</td><td>GGCGATATTTCCGAAGTTGA</td></tr><tr><td rowspan="2">ATP synthase gamma chain</td><td>ATP synth-F</td><td>AACCTCATTATGCGCCTCAC</td></tr><tr><td>ATP synth-R</td><td>GGGCTTCCTGAAGCTTTTCT</td></tr><tr><td rowspan="2">Pre-mRNA splicing factor ini1</td><td>Pre-mRNA-F</td><td>TACCCCTCCACTCCACAAAG</td></tr><tr><td>Pre-mRNA-R</td><td>GGTACACGCCTTGAGAAGGA</td></tr><tr><td rowspan="2">S-phase kinase-associated protein 1A</td><td>S-Phase-F</td><td>ATGATGACAATCGTCGCAAG</td></tr><tr><td>S-Phase-R</td><td>GGGGATTTGCCCTTAATCAT</td></tr></tbody></table>
b

<table><thead><tr><th>EST</th><th>Primer</th><th>Primer sequence 5'-&gt;3'</th></tr></thead><tbody><tr><td rowspan="2">Heat shock protein 70</td><td>Hsp 70-F</td><td>CCCGTCATTGAGGTTGAGTT</td></tr><tr><td>Hsp 70-R</td><td>TGGGAGTCGTTGAAGTAGGC</td></tr><tr><td rowspan="2">Heat shock protein 98/104</td><td>Hsp 98-F</td><td>GAATGATTGGTGCACCTCCT</td></tr><tr><td>Hsp 98-R</td><td>CCATCAGCTGCAAAAGAACA</td></tr><tr><td rowspan="2">Trehalose-6-phosphate synthase</td><td>TpsA-F</td><td>TGGGTGTCTTGCATTGTGAT</td></tr><tr><td>TpsA-R</td><td>CTTTTCAGGGTCGATTCCAA</td></tr><tr><td rowspan="2">Import inner membrane translocase tim50</td><td>Tim50-F</td><td>AGGGTAGCCAGATCCCTGTT</td></tr><tr><td>Tim50-R</td><td>AGACCTCGCCATCCATAGTG</td></tr><tr><td rowspan="2">NADP-dependant alcohol dehydrogenase</td><td>NADP Alcohol DH-F</td><td>TGACAGGCTTGATTTGGTGA</td></tr><tr><td>NADP Alcohol DH-F</td><td>GGTGACCTGAAGGTTCCGTA</td></tr><tr><td rowspan="2">2,3-diketo-5-methylthio-1-phosphopentane phosphatase</td><td>2,3 phosphatase-F</td><td>TCAAGCCAGCACTGATAACG</td></tr><tr><td>2,3 phosphatase-R</td><td>CGGTCTTTCCGTGAAGTCTC</td></tr><tr><td rowspan="2">Conserved hypothetical protein AN5480.2</td><td>AN5480.2-F</td><td>AGAATCTTCCAGCCGTTCAA</td></tr><tr><td>AN5480.2-R</td><td>ATCATGACTCGCGGAGAAAG</td></tr><tr><td rowspan="2">Hypothetical protein EAL90164.1</td><td>EAL90164.1-F</td><td>CCTCCCAACCTCAAACTCAA</td></tr><tr><td>EAL90164.1-R</td><td>GTCATCATCGGACTGGGACT</td></tr><tr><td rowspan="2">AAA family ATPase</td><td>AAA ATPase-F</td><td>GCGACACTCTCAGCAGTCAG</td></tr><tr><td>AAA ATPase-R</td><td>CAAGGGTCAGGAGGACTTCA</td></tr><tr><td rowspan="2">Transcriptional activator HAA1</td><td>HAA1-F</td><td>TCGCATAAAGCGAAAAGGTC</td></tr><tr><td>HAA1-R</td><td>TGTGAATTTGCGGTTGTGTT</td></tr></tbody></table>
Table 2:

List of ESTs isolated using microarrays as significantly down-regulated in thermal stress conditions S1 and S2.

<table><thead><tr><th rowspan="2">Gene homolog</th><th colspan="2">Accession number</th><th colspan="2">Fold change</th><th rowspan="2">Significance</th></tr><tr><th>Swiss-prot</th><th>Genbank</th><th>(S1-C)</th><th>(S2-C)</th></tr></thead><tbody><tr><td colspan="6"><strong>Ribosomal proteins:</strong></td></tr><tr><td>40S ribosomal protein S2</td><td>P25443</td><td></td><td>-1,99</td><td>-1,87</td><td>***</td></tr><tr><td>40S ribosomal protein S5 a / b</td><td>O14277 / Q9P3T6</td><td></td><td>-2,20</td><td>-1,77</td><td>***</td></tr><tr><td>40S ribosomal protein S8 a / b</td><td>O14049 / Q9P7B2</td><td></td><td>-2,21</td><td>-1,69</td><td>***</td></tr><tr><td>40S ribosomal protein S9</td><td>P52810</td><td></td><td>-2,44</td><td>-1,87</td><td>***</td></tr><tr><td>40S ribosomal protein S11</td><td>P79013</td><td></td><td>-2,05</td><td>-1,65</td><td>***</td></tr><tr><td>40S ribosomal protein S12</td><td>O59936</td><td></td><td>-2,85</td><td>-2,20</td><td>***</td></tr><tr><td>40S ribosomal protein S16</td><td>Q42340</td><td></td><td>-2,39</td><td>-1,85</td><td>***</td></tr><tr><td>40S ribosomal protein S17</td><td>P27770</td><td></td><td>-2,12</td><td>-1,53</td><td>***</td></tr><tr><td>40S ribosomal protein S23</td><td>Q873W8</td><td></td><td>-2,56</td><td>-2,19</td><td>***</td></tr><tr><td>40S ribosomal protein S25</td><td>Q7SC06</td><td></td><td>-2,03</td><td>-1,68</td><td>***</td></tr><tr><td>60S ribosomal protein L2</td><td>Q75AP7</td><td></td><td>-1,68</td><td>-1,40</td><td>***</td></tr><tr><td>60S ribosomal protein L10a</td><td>Q7RZS0</td><td></td><td>-3,02</td><td>-2,42</td><td>***</td></tr><tr><td>60S ribosomal protein L11</td><td>Q758S7</td><td></td><td>-1,75</td><td>-1,76</td><td>***</td></tr><tr><td>60S ribosomal protein L14 a</td><td>P36105</td><td></td><td>-2,49</td><td>-1,83</td><td>***</td></tr><tr><td>60S ribosomal protein L15</td><td>O13418</td><td></td><td>-2,92</td><td>-2,21</td><td>***</td></tr><tr><td>60S ribosomal protein L18b</td><td>Q8TFH1</td><td></td><td>-2,17</td><td>-1,90</td><td>***</td></tr><tr><td>60S ribosomal protein L23</td><td>P04451</td><td></td><td>-2,50</td><td>-1,85</td><td>***</td></tr><tr><td>60S ribosomal protein L26b</td><td>Q9FJX2</td><td></td><td>-2,38</td><td>-1,92</td><td>***</td></tr><tr><td>60S ribosomal protein L30</td><td>Q7S7F1</td><td></td><td>-4,02</td><td>-2,27</td><td>***</td></tr><tr><td>60S ribosomal protein L34 b</td><td>P40525</td><td></td><td>-3,76</td><td>-2,76</td><td>***</td></tr><tr><td>60S ribosomal protein L35</td><td>Q8L805</td><td></td><td>-2,31</td><td>-1,65</td><td>***</td></tr><tr><td>60S ribosomal protein L38</td><td>Q9C2B9</td><td></td><td>-2,90</td><td>-2,15</td><td>***</td></tr><tr><td>60S ribosomal protein L44</td><td>P52809</td><td></td><td>-2,83</td><td>-2,25</td><td>***</td></tr><tr><td>60S ribosomal protein L5</td><td>O59953</td><td></td><td>-1,83</td><td>-1,46</td><td>***</td></tr><tr><td>60S ribosomal protein L8</td><td>O13672</td><td></td><td>-1,90</td><td>-1,55</td><td>***</td></tr><tr><td>Ribosome biogenesis protein RPF2</td><td>P36160</td><td></td><td>-1,72</td><td>-1,46</td><td>**</td></tr><tr><td colspan="6"><strong>Transcription and translation proteins:</strong></td></tr><tr><td>Histone H4</td><td>P04914</td><td></td><td>-1,82</td><td>-1,34</td><td>***</td></tr></tbody></table>
<table><tr><td>Translation elongation factor G2</td><td></td><td>EAL91882.1</td><td>-1,85</td><td>-1,50</td><td>***</td></tr><tr><td>Translationally controlled tumor protein homolog</td><td>P35691</td><td></td><td>-2,06</td><td>-1,58</td><td>***</td></tr><tr><td>Transcriptional activator HAA1</td><td>Q12753</td><td></td><td>-2,77</td><td>-2,90</td><td>**</td></tr><tr><td>bZIP transcription factor CpcA</td><td></td><td>EAL89546.1</td><td>-1,93</td><td>-1,96</td><td>**</td></tr><tr><td>U6 snRNA-associated protein LSm3 (splicing)</td><td>Q9Y7M4</td><td></td><td>-1,92</td><td>-1,57</td><td>***</td></tr><tr><td>Mago nashi protein homolog (splicing)</td><td>P49030</td><td></td><td>-1,95</td><td>-1,64</td><td>*</td></tr></table>

**Metabolism proteins:**

<table><tr><td>Xanthine phosphoribosyltransferase 1</td><td>P47165</td><td></td><td>-1,85</td><td>-1,54</td><td>**</td></tr><tr><td>Ornithine decarboxylase</td><td>P27121</td><td></td><td>-2,58</td><td>-2,14</td><td>*</td></tr></table>

**Proteases:**

<table><tr><td>Candidapepsin 2 precursor (Aspartate protease)</td><td>P28871</td><td></td><td>-2,01</td><td>-1,97</td><td>***</td></tr><tr><td>Glutamate carboxypeptidase-like protein</td><td>Q9P6I2</td><td></td><td>-1,63</td><td>-1,37</td><td>*</td></tr></table>

**Miscellaneous proteins:**

<table><tr><td>AAA family ATPase</td><td></td><td>EAL86805.1</td><td>-1,73</td><td>-2,19</td><td>***</td></tr><tr><td>Guanine nucleotide-binding protein (protein G)</td><td>Q01369</td><td></td><td>-2,01</td><td>-1,53</td><td>***</td></tr><tr><td>Prefoldin subunit 6 (actin chaperone protein)</td><td></td><td>EAL92646.1</td><td>-2,18</td><td>-1,58</td><td>***</td></tr></table>

**Hypothetical proteins / unknown proteins:**

<table><tr><td>Hypothetical protein Afu1g07570</td><td></td><td>EAL88465.1</td><td>-1,64</td><td>-1,47</td><td>***</td></tr><tr><td>Hypothetical protein AN0409.2</td><td></td><td>XP_658013.1</td><td>-2,13</td><td>-1,81</td><td>***</td></tr><tr><td>Hypothetical protein AN5226.2</td><td></td><td>XP_662830.1</td><td>-2,21</td><td>-1,64</td><td>**</td></tr><tr><td>Hypothetical protein MG04362.4</td><td></td><td>XP_361917.1</td><td>-2,21</td><td>-2,07</td><td>*</td></tr><tr><td>Unknown protein F17</td><td></td><td></td><td>-1,94</td><td>-1,41</td><td>***</td></tr><tr><td>Unknown protein F26</td><td></td><td></td><td>-1,74</td><td>-1,49</td><td>***</td></tr><tr><td>Unknown protein F6</td><td></td><td></td><td>-2,84</td><td>-2,60</td><td>***</td></tr><tr><td>Unknown protein F9</td><td></td><td></td><td>-2,08</td><td>-1,77</td><td>***</td></tr><tr><td>Unknown protein R16</td><td></td><td></td><td>-3,70</td><td>-2,71</td><td>***</td></tr></table>

<footnote>* Significant</footnote>

<footnote>** Highly significant</footnote>

<footnote>*** Very highly significant</footnote>
**Table 3:**

List of ESTs isolated using microarrays as significantly up-regulated in thermal stress conditions S1 and S2.

<table><thead><tr><th rowspan="2">Gene homolog</th><th colspan="2">Accession number</th><th colspan="2">Fold change</th><th rowspan="2">Significance</th></tr><tr><th>Swiss-prot</th><th>Genbank</th><th>(S1-C)</th><th>(S2-C)</th></tr></thead><tbody><tr><td colspan="6"><strong>Chaperones:</strong></td></tr><tr><td>Heat shock protein 30</td><td>P40920</td><td></td><td>+2,67</td><td>+2,51</td><td>***</td></tr><tr><td>Heat shock protein 60</td><td>O60008</td><td></td><td>+1,74</td><td>+1,43</td><td>***</td></tr><tr><td>Heat shock protein 70</td><td>Q92260</td><td></td><td>+2,98</td><td>+1,93</td><td>***</td></tr><tr><td>Heat shock protein 90</td><td>P40292</td><td></td><td>+3,55</td><td>+2,06</td><td>***</td></tr><tr><td>Heat shock protein 98 / 104</td><td>P31540</td><td></td><td>+2,28</td><td>+2,06</td><td>***</td></tr><tr><td>Mitochondrial Hsp70 chaperone (Ssc70)</td><td></td><td>EAL93290.1</td><td>+1,66</td><td>+1,35</td><td>**</td></tr><tr><td colspan="6"><strong>Antioxidant proteins:</strong></td></tr><tr><td>Cytochrome c oxidase polypeptide Vla</td><td>P32799</td><td></td><td>+1,65</td><td>+1,52</td><td>**</td></tr><tr><td>Superoxide dismutase</td><td>Q9UQX0</td><td></td><td>+1,65</td><td>+1,62</td><td>***</td></tr><tr><td>Flavohemoprotein (Hemoglobin-like protein)</td><td></td><td>EAL84490.1</td><td>+1,69</td><td>+1,26</td><td>***</td></tr><tr><td>Glutathione S-transferase 3</td><td></td><td>AAX07321</td><td>+1,87</td><td>+1,60</td><td>***</td></tr><tr><td colspan="6"><strong>Metabolism proteins:</strong></td></tr><tr><td>Trehalose-6-phosphate synthase tpsA</td><td>O59921</td><td></td><td>+2,36</td><td>+1,67</td><td>***</td></tr><tr><td>2,3-diketo-5-methylthio-1-phosphopentane phosphatase</td><td></td><td>EAL86676.1</td><td>+4,34</td><td>+2,33</td><td>***</td></tr><tr><td>Carboxylesterase</td><td></td><td>EAL92659.1</td><td>+2,14</td><td>+1,60</td><td>***</td></tr><tr><td>L-xylulose reductase</td><td>Q8NK50</td><td></td><td>+3,00</td><td>+1,92</td><td>***</td></tr><tr><td>Mannitol-1-phosphate 5-dehydrogenase</td><td>Q6AHH7</td><td></td><td>+3,22</td><td>+2,18</td><td>***</td></tr><tr><td>Plasma membrane ATPase (proton pump)</td><td>P24545</td><td></td><td>+1,69</td><td>+1,42</td><td>***</td></tr><tr><td>NADP-dependent alcohol dehydrogenase</td><td></td><td>EAL84214.1</td><td>+2,64</td><td>+1,77</td><td>***</td></tr><tr><td colspan="6"><strong>Proteases:</strong></td></tr><tr><td>Aspergillopepsin A</td><td>Q06902</td><td></td><td>+1,68</td><td>+1,63</td><td>**</td></tr><tr><td>Vacuolar protease A</td><td>Q01294</td><td></td><td>+1,78</td><td>+1,63</td><td>***</td></tr><tr><td colspan="6"><strong>Miscellaneous proteins:</strong></td></tr><tr><td>Metalloreductase, putative</td><td></td><td>EAL85757.1</td><td>+1,71</td><td>+1,60</td><td>**</td></tr></tbody></table>
<table><tr><td>Import inner membrane translocase subunit tim50</td><td>Q4WI16</td><td>+1,66</td><td>+1,66</td><td>**</td></tr><tr><td>Protein bli-3 (blue light-inducible protein)</td><td>EAL90736.1</td><td>+2,30</td><td>+2,05</td><td>**</td></tr></table>

**Hypothetical proteins / unknown proteins:**

<table><tr><td>Integral membrane protein</td><td>EAL87419.1</td><td>+2,51</td><td>+2,02</td><td>***</td></tr><tr><td>Predicted protein</td><td>XP_322336.1</td><td>+1,71</td><td>+1,86</td><td>***</td></tr><tr><td>Conserved hypothetical protein</td><td>EAL90164.1</td><td>+3,58</td><td>+2,26</td><td>***</td></tr><tr><td>Hypothetical protein AN0677.2</td><td>XP_658281.1</td><td>+1,76</td><td>+1,46</td><td>**</td></tr><tr><td>Hypothetical protein AN5056.2</td><td>XP_662660.1</td><td>+2,42</td><td>+2,00</td><td>**</td></tr><tr><td>Hypothetical protein AN5480.2</td><td>XP_663084.1</td><td>+3,23</td><td>+1,99</td><td>***</td></tr><tr><td>Hypothetical protein AN5537.2</td><td>XP_663141.1</td><td>+2,12</td><td>+1,68</td><td>***</td></tr><tr><td>Hypothetical protein AN9515.2</td><td>XP_868897.1</td><td>+3,01</td><td>+1,87</td><td>***</td></tr><tr><td>Unknown protein F18</td><td></td><td>+1,77</td><td>+1,66</td><td>***</td></tr><tr><td>Unknown protein F23</td><td></td><td>+1,70</td><td>+1,78</td><td>**</td></tr><tr><td>Unknown protein R10</td><td></td><td>+2,57</td><td>+1,71</td><td>***</td></tr><tr><td>Unknown protein R29</td><td></td><td>+1,96</td><td>+1,60</td><td>**</td></tr><tr><td>Unknown protein R36</td><td></td><td>+2,93</td><td>+1,91</td><td>**</td></tr></table>

* Significant

** Highly significant

*** Very highly significant
Table 4:

Microarray expression ratios and relative expression results obtained with qRT-PCR in S1 and S2 conditions for the 10 selected ESTs. β-actin gene expression was used for normalisation of qRT-PCR results.

<table><thead><tr><th rowspan="3">EST</th><th rowspan="3">Accession number</th><th colspan="2">S1-C</th><th colspan="2">S2-C</th></tr><tr><th colspan="2">Relative expression</th><th colspan="2">Relative expression</th></tr><tr><th>Microarrays</th><th>qRT-PCR</th><th>Microarrays</th><th>qRT-PCR</th></tr></thead><tbody><tr><td>Heat shock protein 70</td><td>Q92260</td><td>+2,98</td><td>+4,41</td><td>+1,93</td><td>+2,71</td></tr><tr><td>Heat shock protein 98 / 104</td><td>P31540</td><td>+2,28</td><td>+3,34</td><td>+2,06</td><td>+3,68</td></tr><tr><td>Trehalose-6-phosphate synthase</td><td>O59921</td><td>+2,36</td><td>+4,22</td><td>+1,67</td><td>+2,07</td></tr><tr><td>Import inner membrane translocase subunit tim50</td><td>Q4WI16</td><td>+1,66</td><td>+1,87</td><td>+1,66</td><td>+1,98</td></tr><tr><td>NADP-dependent alcohol deshydrogenase</td><td>EAL84214.1</td><td>+2,64</td><td>+4,17</td><td>+1,77</td><td>+2,49</td></tr><tr><td>2,3-diketo-5-methylthio-1-phosphopentane phosphatase</td><td>EAL86676.1</td><td>+4,34</td><td>+3,80</td><td>+2,33</td><td>+2,17</td></tr><tr><td>Hypothetical protein AN5480.2</td><td>XP_663084.1</td><td>+3,23</td><td>+2,19</td><td>+1,99</td><td>+1,24</td></tr><tr><td>Conserved hypothetical protein EAL90164.1</td><td>EAL90164.1</td><td>+3,58</td><td>+3,84</td><td>+2,26</td><td>+1,80</td></tr><tr><td>AAA family ATPase</td><td>EAL86805.1</td><td>-1,73</td><td>-2,25</td><td>-2,19</td><td>-2,81</td></tr><tr><td>Transcriptional activator HAA1</td><td>Q12753</td><td>-2,77</td><td>-2,67</td><td>-2,90</td><td>-2,91</td></tr></tbody></table>
**Figure legends:**

Fig. 1. The central composite design used to investigate transcriptional effect of thermal stress in *Penicillium glabrum*.

Fig. 2. Transcriptional analysis of *Penicillium glabrum* under thermal stress (40°C, 120 min) with SSH. (a) Global distribution of BlastX best hits of up- and down-regulated subtracted libraries among organisms; (b) Global view of proteins encoded by ESTs isolated in both up- and down-regulated subtracted libraries. (c) Functional classification of proteins encoded by ESTs isolated in the up-regulated library and (d) the down-regulated subtracted library.

Fig. 3. Results obtained with qRT-PCR (Ct) for potential housekeeping genes in control (C) and thermal stress conditions (S1, S2, S3 and S4). □ ATP synthase gamma chain, ▲ β-actin, ♦ adenosine kinase, ○ S-phase kinase-associated protein 1A, ■ 14-3-3 protein homolog, △ copper-fist DNA binding domain protein, ● FAD dependent sulfhydryl oxidase, ◇ pre-mRNA splicing factor ini1.

Fig. 4. Bland–Altman plot showing the agreement between expression ratios obtained in S1 and S2 conditions, with qRT-PCR and microarray analysis for the 10 selected genes.