<title>Background</title> <p>The adaptation of genome expression to the chemical environment is a complex but crucial challenge for all living cells. Functional genomics analyses in budding yeast have shown that environmental stress responses may involve rapid changes in the expression of up to 30% of the genome. A common response to all stresses, named ESR (Environmental Stress Response), has been described, which consists in the inhibition of the cytosolic translation apparatus and the activation of the energy storage pathways [##REF##11102521##1##]. However, pathways responding specifically to the parameters of the environment also form a key part of the stress response. These pathways involve specific transcriptional modules that rapidly sense the environment as a series of chemical and physical features (e.g. redox, pH, osmolarity, temperature, etc.) and act together to adapt genome expression to the specific nature of each stress [##REF##17158869##2##]. For instance, at least eight different transcription factors act together to define the first-hour response of yeast cells to the toxic metalloid arsenite [##REF##15575969##3##]. These global and rapid responses are highly dynamic, involving sequential waves of gene activation and repression [##REF##11102521##1##,##REF##17158869##2##,##REF##15713640##4##]. This requires tight temporal coordination between different transcriptional routes, which can be achieved in two complementary ways. First, the transcription factors involved in stress responses, despite responding to different signals, may have overlapping sets of targets [##UREF##0##5##]. Second, cross-regulation between transcription factors may ensure the coordinated activation of different pathways [##REF##16556235##6##]. We focus here on the cross-talks between three transcriptional modules responsible for the oxidative stress response, the ubiquitine-mediated protein degradation and the pleiotropic drug resistance, respectively. These cellular pathways exist in all species, from bacteria to mammals and plants. In <italic>S. cerevisiae</italic>, the oxidative stress response is controlled principally by the Yap1p transcription factor of the AP1-like leucine zipper family. Yap1p acts as a secondary sensor for oxidative molecules, and thus responds to a wide spectra of toxic compounds, such as hydrogen peroxide, metals and metalloids, organic nucleophilic molecules and internal metabolic oxidative stress due to the production of toxic by-products during glycolysis [##REF##15575969##3##,##REF##15713640##4##,##REF##12437921##7##,##REF##15367692##8##]. Yap1p recognizes YRE (Yap1p response elements, 5'-TKACTMA-3') in the promoters of genes involved in redox homeostasis and in xenobiotic export at the plasma membrane. The proteasome is involved in both the degradation of damaged or aggregated proteins and in the post-translational regulation of several biological processes, playing a key role in many stress responses [##REF##15571806##9##]. Expression of the genes involved in proteasome biogenesis and activity, and in ubiquitin-dependent proteolysis, is controlled by the C<sub>2</sub>H<sub>2 </sub>zinc finger protein Rpn4p, which recognizes the PACE (proteasome associated control element, 5'-GGTGGCAAA-3') sequence in the promoter of its target genes [##REF##15575969##3##,##REF##10350051##10##,##REF##11248031##11##]. Pleiotropic drug resistance involves the upregulation of membrane proteins involved in drug efflux. The corresponding genes are controlled principally by two Zn<sub>2</sub>Cys<sub>6 </sub>Gal4p-like transcription factors: Pdr1p and Pdr3p. These two transcription factors have largely overlapping sets of targets and recognize the same DNA motif (named PDRE, 5'-TCCGYGGR-3'), but have different roles and regulatory properties [##REF##17158869##2##,##UREF##0##5##,##REF##11415981##12##, ####REF##12882520##13##, ##REF##10734226##14####10734226##14##]. The Yap1p and Rpn4p pathways are simultaneously involved in the yeast response to arsenate [##REF##15575969##3##]. Yap1p acts together with the Pdr1p/Pdr3p pathway to induce a drug specific response to the antifungal drug benomyl [##REF##15713640##4##]. No transcriptional regulation has been described for <italic>PDR1</italic>, but <italic>YAP1</italic>, <italic>RPN4 </italic>and <italic>PDR3 </italic>are induced by stress [##REF##15575969##3##,##UREF##0##5##,##REF##16556235##6##]. In this work, we used selenite as a model stress to investigate further the interactions between these three transcriptional modules. Selenium is an essential oligoelement that replaces the sulfur atom of some methionine and cysteine in proteins involved in various essential cell functions [##REF##17508906##15##]. Selenium is also a promising agent for cancer therapy and anti-aging treatments [##REF##15387717##16##]. However, high doses of selenium are toxic to eukaryotic cells [##REF##16935405##17##]. In yeast, selenium alters genome stability [##REF##17900630##18##] and is detoxified in the vacuole after reacting with glutathione [##REF##17261587##19##]. Yeast cells have a high level of selenium tolerance, and yeast enriched in selenium have been used in therapeutic trials [##REF##14748935##20##]. We showed in this study that toxic doses of selenite activated various yeast stress response pathways, including the proteasome, oxidative stress, iron homeostasis and general stress pathways. We demonstrated that, in these growth conditions, the expression of <italic>PDR1 </italic>and <italic>RPN4 </italic>was coordinated through a positive transcriptional loop. This loop contributed to the optimal Yap1p-dependent oxidative stress induction of several genes encoding membrane proteins, including <italic>FLR1</italic>, <italic>ATR1 </italic>and <italic>FRM2</italic>. This function seemed to be conserved in the pathogenic yeast species <italic>C. glabrata</italic>. Finally, our data provide evidence for direct transcriptional regulation of the iron homeostasis regulator Aft2p and of the multidrug resistance regulator Yrr1p by Yap1p, indicating a broader role for this factor in coordination of the oxidative stress response.</p>
<title>Methods</title> <title>Yeast strains</title> <p>The <italic>Saccharomyces cerevisiae </italic>strains used were all of the BY4742 (<italic>MATa</italic>; <italic>his3Δ1</italic>;<italic>leu2Δ0; lysΔ0; ura3Δ0</italic>) background. The <italic>rpn4Δ</italic>, <italic>yap1Δ</italic>, <italic>pdr1Δ </italic>and <italic>pdr3Δ </italic>strains were purchased from Euroscarf [##UREF##3##59##]. The Pdr1-myc and the Yap1-myc strains have been described elsewhere [##REF##17158869##2##,##REF##11013218##40##]. The <italic>Candida glabrata </italic>strain was CBS418.</p> <title>Growth conditions and time-course analyses of stress responses</title> <p>Cells were grown at 30°C in YPD (1%(w/v) bacto-yeast extract, 2% (w/v) bacto-peptone, 2% (v/w) glucose) to an OD<sub>600 nm </sub>of 0.5. The cultures were then split in two. Sodium selenite (1 mM), hydrogen peroxide (0.3 mM) or benomyl (20 μg/ml) was added to one of the two half-cultures and water or DMSO (mock treatment) was added to the other. The cells were incubated for an appropriate period of time (see text), and were then either flash-frozen in cold ethanol for RNA extraction or treated with formaldehyde for chromatin immunoprecipitation (see below).</p> <title>RNA extractions</title> <p>Cell culture (15 ml) was flash-frozen in 30 ml of absolute ethanol at -80°C. The cells were harvested by centrifugation (4 minutes at 3000 g). The cell pellets were stored at -80°C. Total RNAs was extracted as previously described [##REF##17158869##2##].</p> <title>Transcriptome and quantitative RT PCR analyses</title> <p>The <italic>S. cerevisiae </italic>microarrays used are fully described in Array express ([##UREF##4##60##]; accession numbers A-MEXP-337, A-MEXP-114 and A-MEXP-1064). The <italic>C. glabrata </italic>microarrays are described in the Gene Expression Omnibus database ([##UREF##5##61##]; accession number: GPL3922). We used 10 μg of total RNA for cDNA synthesis and labeling. The microarray experiments were conducted as previously described [##REF##17158869##2##]. Raw data were normalized using global lowess followed by print-tip median methods, with background removal, as implemented in Goulphar [##REF##17059595##62##]. Experiments with wild-type strains were carried out 4 times, with dye swapping. The statistical significance of the differences in expression observed was determined with the TMEV version of SAM, with a FDR of 5% and the exact number of permutations [##REF##11309499##21##,##REF##16939790##63##,##REF##12613259##64##]. Only genes with less than 25% missing values were considered for the SAM analyses. The remaining missing values were imputed by the KNN input method directly in the TMEV application [##REF##11309499##21##,##REF##16939790##63##,##REF##12613259##64##]. Hierarchical clustering was performed using TMEV, with Euclidean distances and average linkage [##REF##16939790##63##,##REF##12613259##64##]. The complete <italic>S. cerevisiae </italic>transcriptome data are available as additional files ##SUPPL##0##1## and ##SUPPL##4##5##. The raw data can be downloaded from the array express database (accession number: E-TABM-439). The <italic>C. glabrata </italic>data can be downloaded from the Gene Expression Omnibus database (accession number: GSE10244).</p> <p>We used 500 ng of total RNA for quantitative RT-PCR, which was performed as previously described [##REF##17158869##2##]. The oligonucleotides used are described in additional file ##SUPPL##9##10##. <italic>ACT1 </italic>was used as a reference for normalization.</p> <title>Chromatin immunoprecipitation</title> <p>Chromatin immunoprecipitation experiments, followed by intergenic microarray or quantitative PCR analyses, were performed as previously described [##REF##17158869##2##]. The <italic>S. cerevisiae </italic>intergenic arrays are described in the Array express database (accession number: A-MEXP-1065). The array results were normalized using the print-tip median [##REF##17059595##62##]. The statistical significance of the ChIP enrichments was assessed with the TMEV version of SAM with a FDR of 1% and the exact number of permutations [##REF##11309499##21##,##REF##16939790##63##,##REF##12613259##64##]. The complete ChIP-chip results are available as additional file ##SUPPL##10##11##. The raw data can be downloaded from the array express database (accession number: E-TABM-437). The sequence of the oligonucleotides used for quantitative PCR can be found in additional file ##SUPPL##9##10##.</p> <title>Data mining</title> <p>Functional analyses and network mapping of the genome-wide data were carried out with T-profiler [##REF##15980543##22##], SGD GO term finder [##REF##14681421##32##] or Yeastract [##REF##16381908##33##], using the default parameters. Promoter sequence analyses were performed with the DNA pattern search tool from RSA tools [##REF##12824373##65##] and the genome Browser tool from the SGD [##REF##12368253##66##].</p>
<title>Results</title> <title>Gene ontology mapping of the selenite response</title> <p>We analyzed the transcriptome of <italic>S. cerevisiae </italic>cells treated with 1 mM sodium selenite for 2, 5, 10, 20, 40, 60 and 80 minutes. The corresponding cDNAs were competitively hybridized on DNA microarrays, with cDNAs obtained from cells mock-treated for an equivalent period. Each experiment was carried out four times, with independent biological samples. We used SAM to evaluate the significance of variations in expression of each gene [##REF##11309499##21##]. The dose of selenite used was the lowest dose that significantly altered cell survival in our preliminary growth assays (data not shown). This toxic dose of selenite induced large changes in the transcriptome of the cells (see additional file ##SUPPL##0##1##). In our data set, about 30% of the yeast ORFs displayed a significant change in expression level for at least two consecutive time points, with similar numbers of up- and downregulated. The earliest significant effects on gene expression were detected at 10 minutes, with the peak response at 40 minutes (see additional file ##SUPPL##1##2## and additional file ##SUPPL##0##1##). We investigated the functional pathways involved in this large and complex response using the t-profiler software [##REF##15980543##22##], which identified the Gene Ontology (GO) categories significantly up- or downregulated over the whole dataset (figure ##FIG##0##1A##, see additional file ##SUPPL##2##3##). T-profiler analysis indicated that iron homeostasis genes were among the first to be induced (figure ##FIG##0##1A##). Selenite also induced an oxidative stress response characterized by the up-regulation of genes involved in redox homeostasis and in proteasome activity (figure ##FIG##0##1A##). Finally, our selenite treatment, similarly to all stresses with deleterious effects on growth, triggered a large ESR [##REF##11102521##1##]. This response was characterized by the general repression of genes involved in cytosolic translation (ribosome biogenesis, rRNA and tRNA processing, etc.) and by the induction of chaperone proteins and genes involved in energy storage and carbohydrate metabolism (figure ##FIG##0##1A##). This probabilistic model of the selenite response was confirmed by the expression profiles of genes from the GO categories cited above (figure ##FIG##0##1B## and see additional file ##SUPPL##3##4##). We also used T-profiler to predict the transcription factors likely to be responsible for the gene expression patterns observed, based on previous genome-wide chromatin immunoprecipitation experiments and DNA motifs in the promoters of the corresponding genes (figure ##FIG##1##2A##). This led to the identification of 14 transcription factors as possibly positively regulating gene expression in response to selenite (figure ##FIG##1##2A##, see additional file ##SUPPL##2##3##). This list includes regulators involved in the iron homeostasis (Aft1p, Aft2p), two regulators of the oxidative stress response (Yap1p and Skn7p) and certain Yap1p homologues (Yap7p, Cin5p and Cad1p), the regulator of proteasome expression (Rpn4p) and several regulators involved in the general stress response (Hsf1p, Msn2p, Msn4p and Gcn4p). Several transcription factors controlling stress response pathways were themselves found to be upregulated at the mRNA level following selenite treatment (figure ##FIG##1##2B##). For instance, <italic>YAP1 </italic>and <italic>RPN4 </italic>displayed a remarkable pattern of co-induction. <italic>CIN5 </italic>(<italic>YAP4</italic>), which encodes a Yap1p homologue involved in redox homeostasis and salt tolerance [##REF##15565582##23##], was also induced in our experiments, together with <italic>YRR1</italic>, which encode a zinc-finger protein conferring cell resistance to drugs such as 4-NQO or reveromycine and which is upregulated by oxidizing agents such as MMS or benomyl [##REF##15713640##4##,##REF##9767597##24##, ####REF##11909958##25##, ##REF##9990050##26####9990050##26##]. More surprisingly, <italic>PDR1</italic>, which encodes the major regulator of pleiotropic drug resistance in <italic>S. cerevisiae </italic>[##REF##12882520##13##], and <italic>AFT2</italic>, which encodes one of the two main regulators of iron homeostasis [##REF##11448968##27##,##REF##16024809##28##], were both induced by selenite, although no transcriptional regulation of these two genes has ever been reported before. We investigated possible connections between the mechanisms regulating these transcription factors.</p> <title>A transcriptional loop connects RPN4 to PDR1</title> <p>We first examined the basis of the regulation of <italic>PDR1 </italic>and <italic>RPN4</italic>. The induction of these two genes by selenite was confirmed by quantitative RT-PCR (figure ##FIG##2##3##). The induction of <italic>RPN4 </italic>in response to oxidative stress was previously shown to be largely dependent on Yap1p [##REF##15575969##3##,##REF##16556235##6##]. <italic>RPN4 </italic>has also been shown to be activated by gain-of-function alleles of <italic>PDR1 </italic>[##REF##11415981##12##] and the <italic>RPN4 </italic>promoter contains a PDRE which is actually bound by Pdr1p <italic>in vivo </italic>and <italic>in vitro </italic>[##REF##17158869##2##,##REF##11918814##29##]. This PDRE was recently shown to be essential to the full induction of <italic>RPN4 </italic>in response to oxidative stress. This effect was principally attributed to Pdr3p, which also recognizes PDRE [##REF##16556235##6##]. However, <italic>pdr1Δ </italic>cells have normal basal levels of <italic>RPN4 </italic>mRNA and of proteasome activity [##REF##11918814##29##] and <italic>RPN4 </italic>is not induced by drugs which efficiently trigger the Pdr1p/Pdr3p multidrug resistance response [##REF##17158869##2##,##REF##15713640##4##,##UREF##0##5##]. Finally, the impact of Pdr1p on the <italic>RPN4 </italic>expression in response to stress still had to be established. We therefore examined the expression of <italic>RPN4 </italic>in the presence of selenite in <italic>pdr1Δ </italic>cells (figure ##FIG##2##3##). We observed a slight decrease of the <italic>RPN4 </italic>induction from 60 minutes after selenite exposure. This result proved that Pdr1p has a modest but significant role in the induction of <italic>RPN4 </italic>in response to oxidative stress.</p> <p>In searches of SGD [##REF##15343339##30##, ####REF##12748633##31##, ##REF##14681421##32####14681421##32##] and Yeastract [##REF##16381908##33##] databases for transcriptional regulators that might account for the selenite-dependent induction of <italic>PDR1</italic>, we identified a conserved Rpn4p recognition element in the promoters of the <italic>PDR1 </italic>orthologues in <italic>Saccharomyces sensu stricto </italic>species. Moreover, Rpn4p was found to bind to the <italic>PDR1 </italic>promoter in a global study of the genomic locations of yeast transcription factors binding sites [##REF##15343339##30##]. We therefore monitored <italic>PDR1 </italic>expression in an <italic>rpn4Δ </italic>strain. The Rpn4p had no apparent role in <italic>PDR1 </italic>basal expression but the inactivation of <italic>RPN4 </italic>severely reduced the sensitivity of <italic>PDR1 </italic>to selenite (figure ##FIG##2##3##). We conclude from these results that a positive transcriptional loop connects <italic>RPN4 </italic>and <italic>PDR1 </italic>in response to selenite.</p> <title>The PDR1/RPN4 loop optimizes the Yap1p dependent oxidative stress response</title> <p>We investigated the role of the <italic>PDR1</italic>/<italic>RPN4 </italic>loop, by carrying out a genome-wide analyses of the contributions to the selenite response of Rpn4p, Pdr1p, Yap1p, which controls <italic>RPN4 </italic>expression in response to metalloids [##REF##15575969##3##], and Pdr3p, the functional homologue of Pdr1p. Fluorescent cDNAs from <italic>yap1Δ</italic>, <italic>rpn4Δ</italic>, <italic>pdr1Δ </italic>or <italic>pdr3Δ </italic>cells treated for 0, 40, 60 or 80 minutes with 1 mM of selenite, were hybridized to microarrays together with cDNA from wild type cells subjected to identical treatment (see additional file ##SUPPL##4##5##). We then focused on 175 genes displaying levels of induction lower than those observed for the wild type in at least one of the mutants for at least two consecutive time points (see additional file ##SUPPL##5##6##). These genes were classified as a function of their gene expression patterns using hierarchical clustering (see additional file ##SUPPL##5##6##). Based on the clustering tree obtained, we defined five groups corresponding to five different patterns of expression (figure ##FIG##3##4A##). We used Yeastract [##REF##16381908##33##] and the SGD GO term finder [##REF##15297299##34##] to address the main transcriptional regulators and functional annotations available for these genes (figure ##FIG##3##4B##). Yap1p was the factor regulating the largest number of these genes (about 80%), which were, as expected, mostly involved in redox homeostasis and the chemical stress response (clusters 1, 2, 4 and 5). Many of these genes had already been annotated as Yap1p targets and their selenite induction was abolished in the absence of <italic>YAP1 </italic>(clusters 1 and 2). Yap1p also had a slight effect on some genes involved in proteolysis (part of cluster 3), possibly indirectly due to the significant influence of Yap1p on <italic>RPN4 </italic>expression (see additional file ##SUPPL##5##6##). <italic>RPN4 </italic>deletion had an effect on about 60% of the genes shown in figure ##FIG##3##4A##. Most of these genes encoded proteins involved in proteolysis (cluster 3), as expected, and one third of them had already been identified as Rpn4p targets (figure ##FIG##3##4B##). More surprisingly, the deletion of <italic>RPN4 </italic>slightly but reproducibly decreased the selenite induction of several Yap1p targets (clusters 1 and 4). The inactivation of <italic>PDR1 </italic>or <italic>PDR3 </italic>altered the expression of few genes (less than 20% of the genes in figure ##FIG##3##4A##), all of which were Yap1p targets (upper parts of clusters 1 and 4) and encoded proteins involved in chemical stress response and xenobiotic export (<italic>FLR1</italic>, <italic>ATR1</italic>, <italic>FRM2</italic>, etc.). Remarkably, these effects were similar to those of the <italic>RPN4 </italic>deletion on these genes (figure ##FIG##3##4C##). We carried out chromatin immunoprecipitation (ChIP) analyses, together with real-time quantitative PCR, to investigate the direct binding of Pdr1p to these Pdr1p-dependent genes, using a myc-tagged version of Pdr1p [##REF##17158869##2##]. We used the <italic>PDR5 </italic>promoter as a positive control, as this sequence constitutively binds Pdr1p [##REF##17158869##2##]. We first focused on the <italic>FLR1 </italic>promoter, which contains a PDRE and is partially controlled by Pdr3p in response to oxidative stress [##REF##10572257##35##]. We detected no significant binding of Pdr1p to the promoter of <italic>FLR1</italic>, despite using several pairs of oligonucleotides to scan the entire <italic>FLR1 </italic>promoter region (figure ##FIG##4##5##). Similarly, the promoters of <italic>ATR1</italic>, <italic>FRM2 </italic>and <italic>YCR102c </italic>did not bind by Pdr1p in the presence of selenite (data not shown). <italic>RPN4 </italic>was the only gene tested that displayed Pdr1p-dependent selenite induction, its promoter binding directly to Pdr1p (Figure ##FIG##4##5##).</p> <title>Conservation of the YAP1/RPN4/PDR1 co-regulation in C. glabrata</title> <p>We assessed the conservation of <italic>RPN4</italic>, <italic>PDR1 </italic>and <italic>YAP1 </italic>co-regulation in the pathogenic yeast <italic>C. glabrata</italic>. This species does not belong to the <italic>Saccharomyces sensu stricto </italic>group but is much more closely related to <italic>S. cerevisiae </italic>than the other pathogenic <italic>Candida </italic>species [##REF##15229592##36##]. Especially, clear <italic>C. glabrata </italic>orthologs of <italic>PDR1 </italic>(called <italic>CgPDR1</italic>), <italic>YAP1 </italic>(<italic>CgAP1</italic>) and <italic>RPN4 </italic>(<italic>CgRPN4</italic>) could be defined, whereas there seems to be no ortholog of <italic>PDR3 </italic>in this species [##REF##16381905##37##]. As in <italic>S. cerevisiae</italic>, <italic>CgRPN4 </italic>is upregulated in strains harboring gain-of-function alleles of <italic>CgPDR1 </italic>[##REF##16803598##38##]. We used <italic>C. glabrata </italic>microarrays to investigate the expression patterns of <italic>CgAP1</italic>, <italic>CgRPN4 </italic>and <italic>CgPDR1 </italic>in response to the oxidative stress caused by the antifungal drug, benomyl (figure ##FIG##5##6A##). A complete comparison of the oxidative stress responses of <italic>C. glabrata </italic>and <italic>S. cerevisiae </italic>is presented elsewhere (Lelandais et al., manuscript in preparation). We observed that benomyl induced the <italic>C. glabrata </italic>homologues of <italic>YAP1</italic>, <italic>RPN4 </italic>and <italic>PDR1 </italic>(figure ##FIG##5##6A##). We therefore analyzed the DNA regions upstream from the <italic>CgRPN4, CgAP1 </italic>and <italic>CgPDR1 </italic>ORFs (figure ##FIG##5##6B##). As in <italic>S. cerevisiae</italic>, two PDRE and one YRE were found in the <italic>CgRPN4 </italic>promoter, whereas one PACE was found in the <italic>CgPDR1 </italic>and <italic>CgAP1 </italic>promoters (figure ##FIG##5##6B##). These data strongly suggest that the cross-regulation of <italic>PDR1 </italic>and <italic>RPN4 </italic>on one hand, and of <italic>YAP1 </italic>and <italic>RPN4 </italic>on the other, is conserved from <italic>S. cerevisiae </italic>to <italic>C. glabrata</italic>. These findings are consistent with this network playing a role in the fine-tuning of oxidative stress responses in yeasts.</p> <title>Expanding the Yap1 network: connections with the Yrr1 and Aft2 regulons</title> <p>The <italic>YRR1 </italic>and <italic>AFT2 </italic>genes were both induced by selenite (figure ##FIG##1##2##). <italic>YRR1 </italic>was shown to be regulated by oxidizing agents such as benomyl or mancozeb in a Yap1p-dependent manner [##REF##15713640##4##,##REF##18086556##39##]. No transcriptional regulation of <italic>AFT2 </italic>has ever before been reported. The selenite-driven expression of <italic>AFT2 </italic>was weaker in the <italic>yap1Δ </italic>strain than in the wild type (see additional file ##SUPPL##5##6## and cluster 5 of figure ##FIG##3##4A##) and one YRE (TTAGTCA) conserved in <italic>Saccharomyces sensu stricto </italic>species, was found in the promoter of this gene 147 base pairs upstream from the ATG codon [##REF##12748633##31##,##REF##14681421##32##]. We conducted genome-wide location analyses of a myc-tagged version of Yap1p [##REF##11013218##40##] using the ChIP-chip technique [##REF##15343339##30##]. Previous experiments indicated that Yap1p could discriminate between different sources of oxidative stress but that its DNA-binding properties were independent of the oxidizing agents used [##REF##15713640##4##,##REF##14556853##41##]. We therefore increased the relevance of our ChIP-chip data, by carrying out analyses of three different sources of oxidative stress: selenite (1 mM), hydrogen peroxide (0.3 mM) and benomyl (20 μg/ml). The doses of hydrogen peroxide and benomyl and the time of exposure (20 minutes) used were chosen based on previous experiments [##REF##15713640##4##,##REF##12437921##7##,##REF##11013218##40##]. Each experiment was carried out independently four times. We used the SAM method [##REF##11309499##21##] to obtain lists of significantly and reproducibly bound targets in the three sets of conditions. We considered as potential targets only those sequences bound in all the oxidative conditions tested. We controlled for false positive by conducting similar experiments in control conditions, on cells treated with DMSO. In these conditions, Yap1p was efficiently excluded from the nucleus [##REF##9130715##42##]. We then deleted the DNA sequences which were significantly enriched in the DMSO experiments, to obtain the final list of Yap1p direct targets presented in the additional file ##SUPPL##6##7##. We identified 310 DNA sequences, corresponding to about 300 promoter regions, that were reproducibly and significantly enriched in all conditions. According to the SGD GO term finder [##REF##15297299##34##], the genes regulated by these promoter regions were mostly involved in the cell response to chemical stimuli (p = 2.21 × 10<sup>-18</sup>), oxidative stress (p = 1.16 × 10<sup>-9</sup>), and drugs (p = 1.35 × 10<sup>-3</sup>); and in sulfur metabolism (p = 4.1 × 10<sup>-4</sup>). This list contained many genes previously known to be regulated by Yap1p (e.g. <italic>TRX2</italic>, <italic>FRM2</italic>, <italic>TRR1</italic>, <italic>TSA2</italic>, <italic>SOD1</italic>, etc.), including <italic>CIN5 </italic>and <italic>RPN4 </italic>(see additional file ##SUPPL##6##7##). Remarkably, 75% of these 310 sequences contained a YRE (see additional file ##SUPPL##6##7##) but only 23% of them were induced by selenite in a Yap1p-dependent manner, which confirmed that Yap1p binding is necessary but not sufficient for transcriptional regulation [##REF##15713640##4##,##REF##15343339##30##,##REF##14556853##41##]. Our ChIP-chip results overlapped with the clusters defined in figure ##FIG##3##4##. The groups of genes which were highly dependent on Yap1p for selenite induction (clusters 1 and 2), contained 72% and 86%, respectively, of sequences binding Yap1p in our experiments. Clusters 4 and 5 included only 50% and 33% of sequences binding Yap1p. Only 5% of the genes in cluster 3, which was mainly composed of targets of Rpn4p, were found to bind Yap1p, confirming that the mild effect of Yap1p on the expression of about half of this group is indirect. Finally, the sequences corresponding to the <italic>YRR1 </italic>and <italic>AFT2 </italic>promoters were significantly and specifically enriched in all the oxidative conditions tested (Figure ##FIG##6##7##). Based on these data, we conclude that the selenite-dependent induction of <italic>YRR1 </italic>and <italic>AFT2 </italic>is directly controlled by Yap1p.</p>
<title>Discussion</title> <title>Network mapping for the selenite response in budding yeast</title> <p>We used selenite as a model stress to decipher part of the transcriptional network controlling the adaptation of the genome expression to toxic environmental conditions in yeast. Our data suggested that one of the earliest effects of selenite on gene expression was the induction of an iron starvation-like response, known to be controlled by the Aft1p and Aft2p transcription factors [##REF##16024809##28##,##REF##14668481##43##]. Selenite may affect iron homeostasis at two, non exclusive, levels. First, selenium can interact with iron, with a high affinity, in the culture medium [##UREF##1##44##,##UREF##2##45##], rendering this metal unavailable to the cell. Second, it may interfere with iron homeostasis by replacing sulfur in iron-sulfur cluster protein biosynthesis in the mitochondria [##REF##15649888##46##]. Interestingly, the sets of yeast mutant strains which accumulating selenium or iron are very similar, suggesting that these two elements are metabolized through similar cellular routes and affect similar cellular processes [##REF##16168084##47##]. Like arsenite [##REF##15575969##3##], selenite induces a strong oxidative stress response, triggering redox homeostasis pathways and proteasome activity. Selenite may affect redox homeostasis in several ways. First, each selenite molecule contains three atoms of oxygen, which may generate reactive oxygen species (ROS) during selenite reduction. Second, selenite is metabolized through interaction with thiol derivatives, including gluthatione [##REF##17261587##19##], probably leading to imbalance in redox homeostasis. Third, as mentioned above, selenite may interfere with iron homeostasis and iron-sulfur cluster protein biosynthesis, potentially affecting mitochondrial activity and redox homeostasis [##REF##17659286##48##]. Another feature of the oxidative stress response generated by metals and metalloids is the upregulation of genes involved in methionine and sulfur metabolism [##REF##15575969##3##,##REF##11983164##49##]. This effect was not detected in our transcriptome analyses of the yeast response to selenite. This was certainly due to our experiments being carried out in rich media, in which the abundance of sulfur amino acids efficiently switched off the MET gene transcription [##REF##12150908##50##], whereas the studies mentioned above were conducted in minimal media. Despite the simplicity of its molecular structure, selenite induces many different stress response pathways. Moreover, the dose of selenite used clearly compromised the cells ability to respond to stress efficiently. We were therefore able to observe the regulation of mRNA levels for many transcriptional regulators involved in stress responses. This made it possible to reveal new connections between these stress response pathways.</p> <title>Transcriptional loops connect proteasome to pleiotropic drug resistance and oxidative stress response</title> <p>The proteasome functions in many cellular processes, some of which are essential for cell survival as cell cycle progression or the adaptation to environmental changes [##REF##15571806##9##,##REF##16861887##51##]. Its activity is therefore tightly regulated. One level of regulation is the expression of its subunits, which is controlled by the Rpn4p transcription factor in yeast [##REF##17418826##52##]. <italic>RPN4 </italic>is itself positively regulated by a complex array of transcriptional controls connected to environmental stress. These include the heat shock factor Hsf1p, the multidrug resistance regulators Pdr1p and Pdr3p and the oxidative stress response major regulator, Yap1p [##REF##15575969##3##,##REF##16556235##6##]. In this work, we confirmed that Pdr1p and Yap1p were required for the full selenite-driven induction of <italic>RPN4 </italic>(figure ##FIG##2##3## and see additional file ##SUPPL##5##6##). Recent data have suggested that <italic>RPN4 </italic>not only is a target for stress response pathways but Rpn4p also has a direct impact on the expression of some of its regulators. Indeed, the response of <italic>YAP1 </italic>to selenite and arsenite was diminished in the absence of <italic>RPN4 </italic>(see additional file ##SUPPL##5##6##,[##REF##15575969##3##]), as Rpn4p binds to a PACE present in the <italic>YAP1 </italic>promoter [##REF##15343339##30##], providing strong evidence in favor of a transcriptional loop connecting <italic>RPN4 </italic>and <italic>YAP1</italic>. Similarly, one of the major findings of this study was that Rpn4p controls the expression of <italic>PDR1 </italic>in response to selenite (figure ##FIG##2##3##). Therefore, <italic>RPN4 </italic>establishes positive feedback loops with both the oxidative stress response and the pleiotropic drug resistance network. These loops seem to optimize part of the oxidative stress response, as the deletion of <italic>RPN4</italic>, <italic>PDR1 </italic>or <italic>PDR3 </italic>decreases the positive effect of Yap1p on some of its target genes by 50% (figure ##FIG##3##4C##). These effects on gene expression have apparently no impact on cell survival in laboratory conditions (see additional file ##SUPPL##7##8##), but evidence for their biological significance is provided by the observation of <italic>RPN4</italic>, <italic>YAP1 </italic>and <italic>PDR1 </italic>co-regulation in response to oxidative stress in the yeast <italic>C. glabrata </italic>(figure ##FIG##5##6A##). Our results are also consistent with those of a recent study showing that the Rpn4p, Yap1p, Pdr3p and Yrr1p transcription factors collaborate in the upregulation of <italic>FLR1 </italic>in response to oxidative stress [##REF##18086556##39##]. Positive feedback loops may induce bistability in biological systems [##REF##17398098##53##]. Bistability has obvious advantages in the responses of microbial cells to environmental changes. After transient exposure to stress, bistability allows some cells to maintain a particular pattern of gene expression long after the stimulus has ended. This provides a mechanism for the anticipation of future environmental changes based on past environmental conditions and can potentially generate heterogeneity in isogenic cell populations by generating bimodal population responses [##REF##17398098##53##,##REF##17880264##54##]. For Rpn4p and Yap1p, the transcriptional positive feedback loop is counterbalanced by a post-translational negative feedback loop, with the levels of Rpn4p and Yap1p being negatively regulated by the proteasome [##REF##11248031##11##,##REF##16485023##55##].</p> <title>A third mode of functioning for the PDR pathway</title> <p>Two modes of functioning of the PDR pathway have been extensively described. First, Pdr1p and Pdr3p are both involved in the basal expression and drug-dependent up-regulation of multidrug resistance transporter genes, such as <italic>PDR5</italic>, <italic>SNQ2 </italic>and <italic>YOR1 </italic>[##REF##12882520##13##]. Second, Pdr3p is involved in a retrograde response negatively connecting the expression of drug resistance transporters to the mitochondrial respiratory activity [##REF##15896930##56##,##REF##11943188##57##]. The tight connection described above between <italic>RPN4 </italic>and <italic>PDR1 </italic>reflects a third mode of functioning for the PDR pathway. In response to oxidative stress, Pdr1p and Pdr3p control the induction of <italic>RPN4 </italic>(figure ##FIG##2##3##, [##REF##16556235##6##]), which in turn affects the expression of Yap1p target genes such as <italic>FLR1 </italic>(figure ##FIG##3##4##, [##REF##10572257##35##]). Remarkably, <italic>FLR1 </italic>and <italic>RPN4 </italic>were not induced by drugs that efficiently trigger the "classic" PDR response, such as progesterone or fluphenazine [##REF##17158869##2##,##UREF##0##5##], or by mitochondrial defects that stimulate Pdr3p activity [##REF##11943188##57##]. <italic>PDR1 </italic>and <italic>PDR3 </italic>also had original patterns of regulation in response to selenite. As described above, <italic>PDR1 </italic>is under the control of Rpn4p (figure ##FIG##2##3##). By contrast, <italic>PDR3 </italic>has been reported to be controlled by Hsf1p in response to oxidative stress [##REF##16556235##6##]. The induction of <italic>PDR3 </italic>by selenite (see additional file ##SUPPL##8##9##) together with the other targets of Hsf1p (see additional file ##SUPPL##3##4##), suggests that the regulation of <italic>PDR3 </italic>by Hsf1p is effective in the presence of selenite (figure ##FIG##7##8##).</p> <title>Yap1p: a central node in the oxidative stress network</title> <p>Yap1p is the major regulator of the genes encoding proteins involved in redox homeostasis in response to various sources of oxidative stress [##REF##11102521##1##,##REF##10347154##58##]. Yap1p has been shown to influence other regulatory pathways, for instance by regulating <italic>RPN4 </italic>and <italic>CIN5 </italic>[##REF##15575969##3##,##REF##15565582##23##]. A previous global analysis of DNA binding sites for most of the transcription factors of yeast provided the first evidence to suggest that Yap1p may regulate the expression of transcription factors [##REF##15343339##30##]. We have shown that the induction of <italic>YRR1 </italic>[##REF##15713640##4##] and <italic>AFT2 </italic>(figure ##FIG##3##4##) in response to stress is dependent on <italic>YAP1 </italic>and that Yap1p actually binds to the promoters of these genes in the presence of various sources of oxidative stress (figure ##FIG##6##7##). In response to both benomyl and selenite, <italic>YRR1 </italic>was induced later than other Yap1p targets (figure ##FIG##1##2##, [##REF##15713640##4##]). Remarkably, several Yrr1p targets (e.g. <italic>FLR1</italic>, <italic>SNG1</italic>, <italic>SNQ2</italic>) [##REF##11909958##25##] were found to be also directly regulated by Yap1p in response to selenite (see additional files ##SUPPL##5##6## and ##SUPPL##6##7##). The role of Yrr1p therefore seems to be to support Yap1p in the long-term regulation of these genes. The physiological connection between iron metabolism and redox homeostasis has been established before. In particular, it has been shown that the activity of the Aft1p transcription factor, which senses iron through iron-sulfur protein biogenesis status [##REF##15649888##46##], is influenced by the glutathione biosynthesis pathway and the Grx3 and Grx4 glutaredoxins, which are involved in the thiol redox system [##REF##17659286##48##]. Moreover, the double deletion of <italic>AFT1 </italic>and <italic>AFT2 </italic>induces cell hypersensitivity to oxidative stress [##REF##11448968##27##]. Our finding that Yap1p controls <italic>AFT2 </italic>expression provides a direct transcriptional connection between the two pathways. It may seem paradoxical that Yap1p positively controls a system that responds to iron starvation, given that iron uptake is likely to cause oxidative stress. Noteworthy, Aft2p specifically controls the expression of genes involved in the transport of iron from the cytosol to the vacuole and mitochondria [##REF##16024809##28##]. By contrast, Aft1p, which actually controls iron uptake from the environment, does not seem to be positively regulated by Yap1p.</p> <p>In conclusion, Yap1p is a central node in the oxidative stress response network, coordinating the expression of at least four transcription factors involved in various stress response pathways (figure ##FIG##7##8##).</p> <title>Structure and dynamics of the regulatory networks driving cell adaptation to environmental changes</title> <p>We analyzed the structure and dynamics of the transcriptional regulatory network which controls the adaptation of yeast transcriptome to toxic doses of selenite (figure ##FIG##7##8##). Our findings hinted several important features of the regulatory networks involved in chemical stress responses. First, these networks are highly interconnected. The cross-regulation of different regulators makes it possible to transmit information between the different transcriptional routes, resulting in the tight coordination of the various cellular pathways required for cell survival. Second, these networks have versatile and dynamic structures and properties. This plasticity is based on the combination of different transcription factors responsive to different chemical and physical parameters, but also on the fact that the same transcription factor can change its protein partners and/or its DNA binding properties to adapt its activity to the physiological conditions (e.g. the three modes of functioning of the Pdr1p/Pdr3p combination). Third, there is a clear hierarchy in these networks, as illustrated by unidirectional regulations (e.g. Hsf1p on <italic>RPN4</italic>) and the unequal weightings of different relationships (e.g. <italic>RPN4</italic>/<italic>YAP1 </italic>and <italic>RPN4</italic>/<italic>PDR1 </italic>loops), which may also be a function of time and physiological status of the cell.</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>Stress responses provide valuable models for deciphering the transcriptional networks controlling the adaptation of the cell to its environment. We analyzed the transcriptome response of yeast to toxic concentrations of selenite. We used gene network mapping tools to identify functional pathways and transcription factors involved in this response. We then used chromatin immunoprecipitation and knock-out experiments to investigate the role of some of these regulators and the regulatory connections between them.</p> <title>Results</title> <p>Selenite rapidly activates a battery of transcriptional circuits, including iron deprivation, oxidative stress and protein degradation responses. The mRNA levels of several transcriptional regulators are themselves regulated. We demonstrate the existence of a positive transcriptional loop connecting the regulator of proteasome expression, Rpn4p, to the pleiotropic drug response factor, Pdr1p. We also provide evidence for the involvement of this regulatory module in the oxidative stress response controlled by the Yap1p transcription factor and its conservation in the pathogenic yeast <italic>C. glabrata</italic>. In addition, we show that the drug resistance regulator gene <italic>YRR1 </italic>and the iron homeostasis regulator gene <italic>AFT2 </italic>are both directly regulated by Yap1p.</p> <title>Conclusion</title> <p>This work depicted a highly interconnected and complex transcriptional network involved in the adaptation of yeast genome expression to the presence of selenite in its chemical environment. It revealed the transcriptional regulation of <italic>PDR1 </italic>by Rpn4p, proposed a new role for the pleiotropic drug resistance network in stress response and demonstrated a direct regulatory connection between oxidative stress response and iron homeostasis.</p>
<title>Abbreviations</title> <p>ESR: environmental stress response; YRE: Yap1p response element; PACE: proteasome associated control element; PDRE: pleiotropic drug response element; ORF: open reading frame; GO: gene ontology; MMS: methyl methane sulfonate; 4-NQO: nitroquinoline oxide; DMSO: dimethyl sulfoxyde; SGD: saccharomyces genome database; ChIP: chromatine immunoprecipitation; ChIP-chip: ChIP associated to microarrays; SAM: significance analysis of microarrays; ROS: reactive oxygen species; MET: methionine biosynthesis pathway; PDR: pleiotropic drug resistance.</p> <title>Authors' contributions</title> <p>HS, VF and VT carried out the transcriptome and Yap1p DNA binding experiments. EP performed the ChIP analyses of Pdr1p DNA binding. GL and SL contributed to the bioinformatic analyses of transcriptome and ChIP-chip results. CJ contributed to the writing of the manuscript. FD directed the work, carried out the transcriptome analyses of mutant strains, analyzed the data and wrote the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>VF was the recipient of a fellowship from the French ministery of research (MESR). HS had a post-doctoral grant from the "Toxicologie nucléaire" French national program. EP was funded by the Italian Istituto Pasteur-Fondazione Cenci bolognetti. We are grateful to the Service de Génomique Fonctionnelle, CEA-Evry, and in particular to Frank Amyot and Xavier Gidrol, who kindly provided the yeast intergenic microarrays. This work was supported by the "Toxicologie nucléaire" program.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Global Gene Ontology mapping of the selenite response in <italic>S. cerevisiae</italic></bold>. (A): The gene ontology categories with levels of expression significantly changed by selenite were identified by T-profiler [##REF##15980543##22##]. The graph represents only the most significant functional categories, as a function of time. More complete T-profiler results are available in additional file ##SUPPL##2##3##. The color code is as follows: black: E-value &gt; 0.05 (non significant variation); red: E-value &lt; 0.05 and t-value &gt; 0 (positive significant variation); green: E-value &lt; 0.05 and t-value &lt; 0 (negative significant variation). (B): Gene expression patterns for the functional categories discussed further in the text. Wild-type cells were treated with selenite and gene expression was evaluated by microarray analyses at different time points, using untreated cells as a reference. Note that the lists of genes represented here are not exhaustive and represent a sample of the genes present in these GO categories. More complete information can be found in additional file ##SUPPL##0##1##.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Transcriptional regulation and selenite response</bold>. (A): T-profiler prediction of the transcription factors involved in the positive response of yeast to selenite. The color code is the same as for figure 1A. These predictions were based on previous ChIP-chip results [##REF##15343339##30##]. (B): Expression patterns of genes encoding transcription factors involved in stress response pathways. Full results are presented in additional file ##SUPPL##0##1##. Wild-type cells were treated with selenite and gene expression was evaluated by microarray analyses at different time points, using untreated cells as a reference.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold><italic>RPN4 </italic>and <italic>PDR1 </italic>are dependent on each other for regulation by selenite</bold>. The levels of expression of <italic>PDR1 </italic>(upper panel) and <italic>RPN4 </italic>(lower panel) after selenite or mock treatment were quantified in wild-type and <italic>rpn4Δ </italic>(upper panel) or <italic>pdr1Δ </italic>(lower panel) cells, using real-time quantitative PCR. The expression values were normalized using the gene encoding actin (<italic>ACT1</italic>, see methods). The values represented here are the ratios of the normalized expression levels of <italic>PDR1 </italic>(upper panel) or <italic>RPN4 </italic>(lower panel) in the presence of selenite, versus the normalized expression levels of these genes in mock experiments. Each measurement was repeated three times on independent samples. The standard errors are indicated.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Respective weights of <italic>RPN4</italic>, <italic>PDR1</italic>, <italic>PDR3 </italic>or <italic>YAP1 </italic>in the selenite response</bold>. (A) Hierarchical clustering of the genes dependent on one or several of the four transcription factors studied for selenite induction. DNA microarrays were used to compare gene expression levels between wild-type and <italic>rpn4Δ</italic>, <italic>pdr1Δ</italic>, <italic>pdr3Δ </italic>or <italic>yap1Δ </italic>exposed to selenite (times 40, 60 and 80 minutes) or mock-treated (time 0). The 175 genes displaying an alteration of selenite induction in at least one mutant strain were clustered into 5 groups (see additional files ##SUPPL##4##5## and ##SUPPL##5##6##). (B): Schematic representation of the importance of each transcription factor in the regulation of the five clusters defined in (A). The arrows symbolize the positive regulation of each cluster by the transcription factor: the width of the arrows indicates the importance of the regulation (large arrows: strong effect, thin arrows: weak effect). Dashed arrows indicate that the transcription factor controls the expression of only some of the genes present in the cluster. Solid arrows mean that all the genes present in the cluster are regulated. The relevance of the Gene Ontology categories and regulatory relationships in each of the clusters was investigated with the SGD GO term finder and Yeastract tools [##REF##16381908##33##,##REF##15297299##34##]. Only the main GO category with a p-value &lt; 0.0001 and the main transcription factor known to regulate the genes in one cluster were indicated. (C): Enlargement of the part of cluster 1 containing the genes most sensitive to the deletion of <italic>PDR1 </italic>and <italic>PDR3</italic>. Gene names are indicated. The presence of a PDRE in the promoters of these genes is indicated by black dots.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Pdr1p binds to the promoter of <italic>RPN4 </italic>but not to that of <italic>FLR1 </italic>in the presence of selenite</bold>. The binding of Pdr1p to promoters was assessed by combining chromatin immunoprecipitation with real-time quantitative PCR, using strains harboring a chromosomal tagged version of Pdr1p (myc-Pdr1p). Sequence enrichment in the ChIP (i.e. ChIP/whole cell extract ratio) was normalized using the <italic>ACT1 </italic>ORF as a reference. Similar experiments were conducted on cells with the untagged Pdr1p as a negative control ("mock ChIP"). The <italic>PDR5 </italic>promoter was used as a positive control for Pdr1p binding. The results shown for <italic>FLR1 </italic>were obtained using a pair of oligonucleotides spanning the PDRE motif present in the FLR1 promoter. The cells were exposed to 1 mM of selenite (+) or mock-treated (-) for 60 minutes before the beginning of the ChIP procedure.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold><italic>PDR1</italic>, <italic>RPN4 </italic>and <italic>YAP1 </italic>are induced by oxidative stress in <italic>C. glabrata</italic></bold>. (A) DNA microarrays were used to analyze the transcriptome response of <italic>C. glabrata </italic>cells to the oxidizing agent benomyl. The data obtained were analyzed and discussed fully in another manuscript (Lelandais et al., in preparation). (B) We searched for consensus DNA binding sites for Rpn4p, Yap1p and Pdr1p in the <italic>CgRPN4, CgAP1 </italic>and <italic>CgPDR1 </italic>promoters using RSA tools software [##REF##12824373##65##]. The positions of the motifs identified are indicated relative to the start codon of the corresponding ORF. The <italic>Saccharomyces sensu stricto </italic>species used here are <italic>S. cerevisiae</italic>, <italic>S. paradoxus</italic>, <italic>S. bayanus</italic>, <italic>S. kuudriazvezii </italic>and <italic>S. mikatae</italic>. The positions of the motifs indicated for these species refer to the <italic>S. cerevisiae </italic>genes.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Yap1p binds to the promoters of <italic>YRR1 </italic>and <italic>AFT2 </italic>in oxidative conditions</bold>. The global DNA binding pattern of Yap1p was determined by chromatin immunoprecipitation using strains harboring a tagged version of Yap1p (myc-Yap1p). The immunoprecipitate DNA was hybridized with intergenic DNA microarrays, together with genomic DNA (see methods). The cells were previously treated with benomyl, hydrogen peroxide or selenite (oxidative conditions); or with DMSO (control condition). Enrichment ratios were normalized using the print-tip median method (see methods). Each experiment was carried out four (oxidative conditions) or two (control condition) times. The histograms indicate the enrichment of the sequence in the immunoprecipitate, normalized with respect to the value obtained for whole cell extract. The results obtained for the mitochondrial ORF <italic>COX1 </italic>are shown as a negative control. The standard errors are indicated. The full results are shown in additional file ##SUPPL##10##11##.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>A highly interconnected transcriptional network is involved in the response to selenite</bold>. This diagram indicates the regulatory interactions involved in the response to selenite, as identified through this study. All the arrows represent positive transcriptional regulations of transcription factor encoding genes or of groups of target genes (symbolized then by functional categories) by a transcription factor. Red arrows indicate interactions demonstrated for the first time by this work. The arrows in bold indicate high-weighted interactions and the thin arrows symbolize low-weighted interactions, according to the results of this study (in the case of Rpn4p, Yap1p, Pdr1p, Pdr3p) or measurements published by Hahn et al (in the case of Hsf1p) [##REF##16556235##6##].</p></caption></fig>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Transcriptome analyses of the response to selenite in wild-type <italic>S. cerevisiae </italic>cells</bold>. The ORF names are indicated, according to the SGD annotation. For each ORF, the numbers are the median Log2 values of the ratio of expression between treated and untreated cells at various time points, calculated from four independent experiments. The SAM results are indicated as follows: NS: non significant; NA: not applicable (less than 3 independent measurements available); S: significant. At 2 and 5 minutes, no variation in gene expression was considered significant by SAM. We therefore indicate no SAM results for these time points.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Number of genes significantly up-or downregulated for two consecutive time points in the selenite response.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>T-profiler results</bold>. Upper part: The gene ontology categories displaying significant changes in expression in response to selenite were identified by t-profiler analysis [##REF##15980543##22##]. The t-values are indicated. Only T-values in bold were significant (E-value &lt; 0.05). Only GO categories significant for at least two consecutive time points are presented in this table. Lower part: The transcription factors with ChIP-chip target genes [##REF##15343339##30##] displaying a significant change in expression during the selenite response were identified by t-profiler [##REF##15980543##22##]. The t-values are indicated. Only T-values in bold were significant (E-value &lt; 0.05). Only transcription factors significant for at least two consecutive time points are presented in this table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p><bold>Gene expression patterns for Gene Ontology categories identified by t-profiler, following exposure to selenite</bold>. Wild-type cells were treated with selenite and gene expression levels were evaluated by microarray analysis, using untreated cells as a reference. Note that the lists of genes given is not exhaustive and corresponds to a sample of the genes present in these GO categories. More complete information can be found in additional file ##SUPPL##0##1##.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p><bold>Transcriptome analyses of the response to selenite in mutant cells</bold>. The ORF names are indicated, according to the SGD annotation. For each ORF, the numbers are the median Log2 values of the ratio of expression between wild-type and mutant cells at different time points, calculated from two independent experiments. NA: not available.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Expression values associated with the eisengram in Figure </bold>##FIG##3##4A##. The gene and PRF names are indicated, according to the SGD annotation. Genes are ordered according to the hierarchical cluster of figure ##FIG##3##4A##. For each gene, the numbers are the median Log2 values of the ratio of expression between wild-type and mutant cells at different time points, calculated from two independent experiments.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p><bold>Myc-Yap1p ChIP-chip results</bold>. The global DNA binding pattern of Yap1p was determined by chromatin immunoprecipitation using strains harboring a tagged version of Yap1p (myc-Yap1p). The immunoprecipitate obtained was hybridized with intergenic DNA microarrays, together with whole cell genomic DNA (see methods). The cells were previously treated with benomyl, hydrogen peroxide or selenite (oxidative conditions), or with DMSO (reducing condition). SAM was used to sequences significantly enriched in each set of conditions. This table lists the features significantly enriched in all oxidative conditions but not in the reducing condition (negative control). The genes which are potentially under the transcriptional control of the corresponding DNA regions are indicated, based on a previous annotation [##REF##15343339##30##]. The presence of YRE in the corresponding sequences and the presence of the regulated genes in one of the five clusters defined in figure ##FIG##3##4## are indicated.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p><bold>Phenotype of strains deleted for <italic>RPN4</italic>, <italic>PDR1</italic>, <italic>PDR3 </italic>or <italic>YAP1 </italic>in presence of selenite</bold>. The cells were grown in YPD to an OD<sub>600 nm </sub>of 0.25 (early exponential phase). They were then treated by 0, 0.5, 0.75 or 1 mM of sodium selenite. The graphs represent the OD<sub>600 nm </sub>(Y axis) as a function of the time of exposure to selenite (X-axis). The wild type, <italic>pdr1Δ </italic>and <italic>pdr3Δ </italic>strains exhibited the same sensitivity to selenite in these conditions. The <italic>rpn4Δ </italic>and <italic>yap1Δ </italic>strains were more sensitive than the wild type and this defect is more severe in the case of <italic>yap1D </italic>cells.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p><bold>Expression profile for <italic>PDR3 </italic>in response to selenite in wild-type yeast cells</bold>. Levels of <italic>PDR3 </italic>expression were quantified in wild-type cells, by real-time quantitative PCR. Expression values were normalized, using the gene encoding actin (<italic>ACT1</italic>, see methods). The values shown here are the ratios of normalized levels of <italic>PDR3 </italic>expression in the presence of selenite to normalized levels of expression of this gene in mock experiments. Each measurement was repeated three times, on independent samples. The standard errors are indicated.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional file 10</title><p>Oligonucleotides used for quantitative PCR analyses of gene expression or transcription factor DNA binding.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S11"><caption><title>Additional file 11</title><p><bold>Genome-wide chromatin immunoprecipitation of Yap1p in different conditions</bold>. The ORF names are indicated, according to the SGD annotation. The intergenic regions are annotated as follows: iORFX is the intergenic region 5' to ORF X, when ORF X is on the Crick strand and 3' when ORF X is on the Watson strand. For each feature, the numbers are the median Log2 values of the ratio of expression between the immunoprecipitate and the whole cell extract, calculated from four independent experiments. Standard deviations (SD) are indicated. NA: value not available.</p></caption></supplementary-material>
<graphic xlink:href="1471-2164-9-333-1"/> <graphic xlink:href="1471-2164-9-333-2"/> <graphic xlink:href="1471-2164-9-333-3"/> <graphic xlink:href="1471-2164-9-333-4"/> <graphic xlink:href="1471-2164-9-333-5"/> <graphic xlink:href="1471-2164-9-333-6"/> <graphic xlink:href="1471-2164-9-333-7"/> <graphic xlink:href="1471-2164-9-333-8"/>
<media xlink:href="1471-2164-9-333-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S2.tiff" mimetype="image" mime-subtype="tiff"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S4.tiff" mimetype="image" mime-subtype="tiff"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S5.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S6.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S7.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S8.tiff" mimetype="image" mime-subtype="tiff"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S9.tiff" mimetype="image" mime-subtype="tiff"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S10.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-333-S11.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Banerjee", "Lelandais", "Shukla", "Mukhopadhyay", "Jacq", "Devaux", "Prasad"], "given-names": ["D", "G", "S", "G", "C", "F", "R"], "article-title": ["The steroid responses of pathogenic and non-pathogenic yeast species enlighten the functioning and evolution of multidrug resistance transcriptional networks"], "source": ["Eukaryot Cell"], "year": ["2007"]}, {"surname": ["Su", "Suarez"], "given-names": ["C", "DL"], "article-title": ["Selenate and selenite sorption on iron oxides: an infrared and electrophoretic study"], "source": ["Soil Sci Soc Am J"], "year": ["2000"], "volume": ["64"], "fpage": ["101"], "lpage": ["111"]}, {"surname": ["Seby", "Potin-Gautier", "Giffaut", "Donard"], "given-names": ["F", "M", "E", "OFX"], "article-title": ["Assessing the speciation and the biogeochemical processes affecting the mobility of selenium from a geological repository of radioactive wastes to the biosphere."], "source": ["Analysis"], "year": ["1998"], "volume": ["26"], "fpage": ["193"], "lpage": ["198"]}, {"article-title": ["EUROSCARF"]}, {"article-title": ["ArrayExpress"]}, {"article-title": ["GEO"]}]
BMC Genomics. 2008 Jul 15; 9:333