Genome-Wide Characterization of the Phosphate Starvation Response in Schizosaccharomyces pombe

Background Inorganic phosphate is an essential nutrient required by organisms for growth. During phosphate starvation, Saccharomyces cerevisiae activates the phosphate signal transduction (PHO) pathway, leading to expression of the secreted acid phosphatase, PHO5. The fission yeast, Schizosaccharomyces pombe, regulates expression of the ScPHO5 homolog (pho1+) via a non-orthologous PHO pathway involving genetically identified positive (pho7+) and negative (csk1+) regulators. The genes induced by phosphate limitation and the molecular mechanism by which pho7+ and csk1+ function are unknown. Here we use a combination of molecular biology, expression microarrays, and chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) to characterize the role of pho7+ and csk1+ in the PHO response. Results We define the set of genes that comprise the initial response to phosphate starvation in S. pombe. We identify a conserved PHO response that contains the ScPHO5 (pho1+), ScPHO84 (SPBC8E4.01c), and ScGIT1 (SPBC1271.09) orthologs. We identify members of the Pho7 regulon and characterize Pho7 binding in response to phosphate-limitation and Csk1 activity. We demonstrate that activation of pho1+ requires Pho7 binding to a UAS in the pho1+ promoter and that Csk1 repression does not regulate Pho7 enrichment. Further, we find that Pho7-dependent activation is not limited to phosphate-starvation, as additional environmental stress response pathways require pho7+ for maximal induction. Conclusions We provide a global analysis of the transcriptional response to phosphate limitation in S. pombe. Our results elucidate the conserved core regulon induced in response to phosphate starvation in this ascomycete distantly related to S. cerevisiae and provide a better understanding of flexibility in environmental stress response networks.

We provide a global analysis of the transcriptional response to phosphate limitation in S. pombe. Our results elucidate the conserved core regulon induced in response to phosphate starvation in this ascomycete distantly related to S. cerevisiae and provide a better understanding of flexibility in environmental stress response networks.

KEYWORDS
phosphate starvation pathway, pho7 + , csk1 + , S. pombe ChIP-Seq, gene expression BACKGROUND Inorganic phosphate (Pi) is an essential nutrient required for signal transduction, energy metabolism, and biochemistry in all organisms. Maintaining a constant, stable concentration of internal inorganic phosphate is a major challenge for biological systems. Because external concentrations of inorganic phosphate fluctuate unpredictably, microorganisms have evolved strategies to sense external phosphate levels [1][2][3], communicate this information to the nucleus [4,5], and induce transcription to respond to phosphate flux [6][7][8]. The phosphate signal transduction (PHO) pathway in the budding yeast, Saccharomyces cerevisiae, is the most thoroughly studied example of phosphate homeostasis in eukaryotes [9][10][11].
The transcription factors Pho4 and Pho2 play a key role in the phosphate starvation response in S. cerevisiae. When cells are grown in conditions where inorganic phosphate is plentiful, Pho4 is multiply phosphorylated by the cyclindependent kinase-cyclin (CDK-cyclin) complex, Pho85-Pho80 [12]. When Pho4 is phosphorylated, it is localized to the cytoplasm [13,14], does not interact with Pho2 [15], and the PHO regulon is not expressed. During phosphate starvation, the CDK inhibitor Pho81 binds to the secondary metabolite myo-D-inositol heptakisphosphate (IP 7 ) and inhibits the Pho85-Pho80 complex [16,17]. Inhibition of Pho85-Pho80 allows Pho4 to be dephosphorylated, enter the nucleus [18], co-operate with Pho2 [7], and induce a set of genes responsible for harvesting inorganic phosphate from the environment [6]. Pho4 function can be conveniently monitored by measuring the activity of the secreted acid phosphatase, Pho5, which is one of the most highly induced members of the PHO response [19,20]. The genes that comprise the PHO regulon have been well characterized and the precise sites within the genome where Pho4 binds during phosphate starvation are known [7,21]. Pho4 regulation occurs in response to changes in external phosphate levels and Pho4 activity is not thought to be regulated by other stress responses.
In this study we ask the following: is the PHO transcriptional response observed in S. cerevisiae conserved in the distantly related ascomycete, Schizosaccharomyces pombe? S. pombe presents an interesting opportunity for addressing this question because: (1) S. pombe did not experience a recent whole-genome duplication eventthought to contribute to specialization [22] -possibly preventing the PHO response from developing a dedicated regulatory network; (2) the orthologs for the PHO pathway either do not exist (PHO81, PHO2, PHO4) or are not involved in the PHO response (PHO80, PHO85) in S. pombe [23]; and (3) recent work utilizing a deletion collection in S. pombe has outlined a basic regulatory structure for the Pi-inducible, secreted acid phosphatase pho1 + (the ortholog to PHO5) creating an opportunity for comparison with the S. cerevisiae PHO response [24]. During phosphate starvation S. pombe Pho7 -a putative transcription factor containing a Zn 2 Cys 6 binuclear cluster [25] -activates pho1 + expression. Csk1 -a CDK-activating kinase-activating kinase (CAKAK) [26] -represses pho1 + expression in high-Pi conditions. Epistasis analysis indicates that Pho7 acts downstream of Csk1.
In this study, we explore how these factors affect transcriptional output by characterizing the PHO transcriptional response in S. pombe. We analyze this response as a function of phosphate, Pho7, and Csk1 availability using DNA microarrays. We delineate a core PHO transcriptional response comprising the genes pho1 + , SPBC8E4.01c (an S. cerevisiae PHO84 ortholog), and SPBC1271.09 (an S. cerevisiae GIT1 ortholog), whose induction in response to phosphate starvation is conserved between S. cerevisiae and S. pombe. Interestingly, while these three genes share a functionally analogous regulatory pathway (i.e. activation through a transcription factor that is normally repressed by a kinase) we find that the mechanism for regulation differs widely between species. Our analysis of the Pho7-regulated transcriptional outputcoupled with a global profile of Pho7 binding to promoters of stress responsive genesleads us to the conclusion that, unlike Pho4, Pho7 plays a role in multiple stress response pathways. We conclude that while there is a core PHO transcriptional response shared between these two ascomycetes, the systems logic and specialization of PHO components varies widely.

pho7 + and csk1 + Regulate a Core Subset of the PHO Response in S. pombe
The kinetics and maximal output of transcription vary widely between different environmental stress response pathways [27][28][29]. To outline the specific PHO response in S. pombe for subsequent analysis, and to avoid indirect activation of non-phosphate starvation regulated genes, we performed a single time-dependent, genome-wide expression analysis of wild-type S. pombe cells in medium lacking inorganic phosphate.
The starvation time-course revealed two distinct responses to phosphate starvation. The rapid response contained 63 genes that exhibited an increase in expression at 120 minutes post-starvation (red lines in Additional File 1, genes listed in Additional File 6; see materials and methods for gene selection criteria). This class contains the secreted acid-phosphatase, pho1 + (orthologous to ScPHO5), SPBC8E4.01c (orthologous to ScPHO84), and SPBC1271.09 (orthologous to ScGIT1).
As pho1 + and SPBC8E4.01c induction has been previously observed in response to Pi starvation [24,30], we believe that this set accurately reflects the genes that respond rapidly to changes in external Pi. In contrast, the slower response (86 additional genes, blue lines in Additional File 1, genes listed in Additional File 6; see materials and methods for gene selection criteria) was significantly enriched for genes previously implicated in a generalized stress response [27]. We focused our attention on the fast responding genes to avoid indirect effects caused by persistent stress in cells.
Previous work indicated that pho7 + and csk1 + are important regulators of pho1 + expression [24]; we expected that they would also play a significant role in regulating additional components of the PHO response. To test our hypothesis, we probed the transcriptional profiles of wild-type, pho7Δ, csk1Δ, and pho7Δcsk1Δ strains in high-Pi and no-Pi conditions at 120 minutes post-starvation using DNA hybridization microarrays.
We identified 22 genes induced in response to Pi limitation ( Figure 1A, first column, applying a p-value cutoff of 0.10; see materials and methods for gene selection criteria). If expression of these genes is dependent on the activity of pho7 + , then their transcript abundance should decrease in a pho7Δ strain when compared to a wild-type background (both starved of Pi). Induction of 31.8% (7/22) of these genes in Pi starvation is pho7 + dependent ( Figure 1A, third column). If the repressor (csk1 + ) prevents this induction by pho7 + in high-Pi conditions, then transcript abundance in a csk1Δ strain compared to the csk1 + background in high-Pi conditions should increase.
Three genes display this response ( Figure 1A, fourth column). A complete listing of all genes regulated by Pi, pho7 + , and/or csk1 + , along with their orthologs in S. cerevisiae can be found in Additional File 7. Finally, comparisons of the pho7Δ, csk1Δ, and pho7Δcsk1Δ double deletion strains confirm the previously described epistatic relationship between pho7 + and csk1 + (i.e. a loss of pho7 + in a csk1Δ switches constitutive expression to uninducible expression) ( Figure 1A, last two columns).
When phosphate is limiting, S. cerevisiae Pho4, along with Pho2, induces the transcription of genes required for phosphate acquisition [7]. The orthologs of Pho4, Pho2, Pho81, and Pho80 are not found in S. pombe or involved in the PHO response, raising the question: does a functionally analogous signaling pathway involving pho7 + and csk1 + regulate the PHO transcriptional response in S. pombe? Comparing the Pho4, Pi-starvation induced genes in S. cerevisiae [7] with the pho7 + , Pi-starvation induced genes in S. pombe reveals an overlap of only three orthologs ( Figure 1B). For pho1 + , SPBC8E4.01c, and SPBC1271.09 a similar system of transcription factor activation, with repression by a kinase in high-Pi conditions, occurs. Unlike Pho85-Pho80 regulation of Pho4 in S. cerevisiae, most of the pho7 + -mediated response is independent of csk1 + regulation. Further, a large segment (82%) of the observed S. pombe PHO response is not conserved between S. pombe and S. cerevisiae. Therefore, the primary role for pho7 + and csk1 + in the PHO pathway is regulating the core set of phosphate harvesting and transport genes (pho1 + , SPBC8E4.01c, and SPBC1271.09).

Pho7 is Enriched at the PHO Core Promoters During Pi Starvation
Pho7 is classified as a putative transcription factor because it possesses a Zn 2 Cys 6 binuclear cluster (ZC), a DNA binding domain for a number of transcription factors [31,32]. To test if Pho7 binds to pho7 + -regulated promoters, cells containing a functional, epitope tagged version of Pho7 (Pho7-TAP) were grown in high-Pi or no-Pi medium and purified DNA associated with Pho7 was processed via high-throughput sequencing (ChIP-Seq, see materials and methods).
Surprisingly, there is widespread Pho7 binding even in high-Pi conditions (1676 peaks out of 4208 passed peak-ID thresholds; see materials and methods) (Additional File 2). During Pi starvation 367 Pho7-bound sites exhibit an increase in Pho7 enrichment. The highest levels of enrichment were observed in the promoters of Pho7regulated genes identified in the microarray analysis (for a complete list of the identified peaks, see Additional File 8). Further, there is a distinct overlap between genes whose expression levels are regulated by Pi and/or pho7 + and those whose promoters display enrichment in Pho7 binding. 13 of the 22 Pi dependent genes (p-value = 8.1e-8), 16 of the 29 pho7 + dependent genes (p-value = 9.5e-9), and 6 of the 7 genes regulated by both Pi and pho7 + (p-value = 1.3e-5) have promoters that are enriched for Pho7 binding in no-Pi medium (p-values were determined using a hypergeometric test, see materials and methods for details). These results are very different from the global binding profile of Pho4 in S. cerevisiae. In that system, Pho4 is only recruited to the promoters of PHO regulated genes during phosphate starvation and, even then, only to relatively few locations (115) within the genome [7].
Pho7 binding is significantly enriched in the promoters of pho1 + , SPBC8E4.01c, SPBC1271.09, SPAC9E9.09C, SPBPB2B2.05, and SPCC330.06c in no-Pi conditions ( Figure 2). In the case of pho1 + , SPBC8E4.01c, SPAC9E9.09C, and SPBPB2B2.05 there is enrichment of Pho7 even in high-Pi conditions (when compared to mock treated samples). This explains the previously noted basal expression of pho1 + in S. pombe [30,33] -it appears that Pho7 is bound and activating transcription even in the absence of stress ( Figure 1A, second column). Based on our microarray results, we note that increased Pho7 enrichment during Pi starvation is not due to an increase in pho7 + transcript abundance (0.05 log 2 fold-change in no-Pi vs. high-Pi conditions). Moreover, we observe no significant difference in Pho7-TAP protein levels during Pi starvation (Additional File 3). Together these observations suggest that the enrichment of Pho7 in the promoters of Pho7-regulated genes is attributable to increased affinity for the promoter and not due to an increase in Pho7 abundance. Further, the number of distinct Pho7 binding events varies between different promoters (Figure 2), and we do not observe a clear correlation between the number of binding events and transcriptional up-regulation. Our attempts to identify a DNA binding motif for Pho7 based solely on our ChIP-Seq data were unsuccessful (data not shown). strain (pho1 + and SPBC1271.09), which is still well below the enrichment seen during Pi starvation ( Figure 3). As we observed during Pi starvation, the loss of Csk1 does not influence either pho7 + transcript abundance (-0.12 log 2 fold-change, [csk1Δ vs. csk1 + ] +Pi) or Pho7-TAP protein levels (Additional File 3). We draw two conclusions from these data: (1) the level of Pho7 bound in high-Pi conditions would be sufficient to induce high levels of transcription if not for the repressive action of Csk1; and (2) Csk1 does not repress Pho7 activity by preventing Pho7 from binding to the promoters of responsive genes.

Module Control pho1 + Expression
Based on our ChIP-Seq results, we know Pho7 binds between nucleotides -280 and -180 in the pho1 + promoter. To determine whether the sequences in this region are necessary and/or sufficient for Pho7-dependent, Pi-starvation induced expression, we utilized an in vivo strategy for confirming Pho7-promoter interactions using exogenous expression plasmids. Briefly, differing lengths of the pho1 + promoter driving the expression of yellow fluorescent protein (yfp + ) were constructed in a vector and transformed into pho7 + , pho7Δ, and csk1Δ backgrounds ( Figure 4A). yfp + expression was measured using FACS with mean YFP intensity serving as a proxy for promoter activity (see materials and methods). The 2 kb segment of the pho1 + promoter activates yfp + expression during Pi-starvation -it exhibits a ~7-fold increase in YFP intensity upon starvation ( Figure 4B). This induction is dependent on Pho7, as it is abolished in a pho7Δ background ( Figure 4B). Surprisingly, we find that loss of Csk1 does not result in high-levels of yfp + expression with the 2 kb reporter in csk1Δ cells grown in either high-Pi or no-Pi growth media. While the levels of expression are above those seen in pho7 + cells grown in high-Pi conditions and pho7Δ cells grown in any condition, they are significantly lower than the levels observed in a pho7 + background in no-Pi conditions ( Figure 4B). Based on our evidence that trimming the pho1 + promoter beyond 1 kb results in high level, constitutive expression of yfp + in csk1Δ cells ( Figure 4D, discussed below), our results with the 2 kb fragment in the csk1Δ background might be complicated by: (1) an additional repressor element located between 2 kb and 280 bp in the pho1 + promoter; and/or (2) differences between the chromatin structure/promoter architecture of the endogenous pho1 + locus and the exogenous 2 kb pho1 + pr-yfp + vector that are influencing transcription in the csk1Δ cells. Therefore, we focused on the behavior of the shorter construct (280 bp pho1 + pr-yfp + ; Figure 4D) as a proxy for the interaction of Csk1 at the region bound by Pho7-TAP in the pho1 + promoter.
If the pho1 + promoter sequence bound by Pho7-TAP in the ChIP-Seq experiment is necessary for activation of pho1 + transcription during Pi-starvation, then deletion of this region should result in a loss of yfp + expression during Pi-starvation. To test this hypothesis, we generated a construct in which the 20 bp centered under the Pho7-TAP ChIP-Seq signal were deleted ( Figure 4A, UASΔ pho1 + pr-yfp + ). In high-Pi growth conditions, the loss of the Pho7 bound region results in a slight increase in yfp + expression ( Figure 4C, compare blue columns to the blue columns in Figure 4B), and in Pi-starvation this UASΔ pho1 + pr-yfp + construct no longer fully activates yfp + expression ( Figure 4C, compare column 2 to column 2 in Figure 4B). The loss of pho7 + results in a further decrease in expression from this UASΔ pho1 + pr-yfp + construct ( Figure 4C, lane 4). It is possible that Pho7 recognizes additional segments of the promoter, though such contributions to activation in Pi-limiting conditions are modest. Together these results demonstrate that the Pho7-TAP bound promoter element is necessary for Pho7dependent transcriptional activation during Pi-limitation. We have termed this region the Pho7 upstream activating sequence (UAS).
To test whether the Pho7 UAS is sufficient for Pho7-dependent, Pi-limitation induced transcriptional activation we deleted all but the first 280 bp of the pho1 + promoter and assayed in vivo yfp + expression ( Figure 4A, 280 bp pho1 + pr-yfp + ). yfp + expression from the 280 bp pho1 + pr-yfp + construct is elevated in high-Pi conditions, and is only marginally activated during Pi-starvation ( Figure 4D, pho7 + , 1st column).
Expression from this construct is reduced in a pho7Δ background and is unaffected by Pi limitation -thus, expression in high Pi conditions and the modest expression increase in Pi-limiting conditions are dependent on Pho7 ( Figure 4D, pho7Δ, 2 nd column). In each background tested, the mean YFP intensities from the 280 bp pho1 + pr-yfp + construct vary by less than 1.5 fold between high-Pi and no-Pi conditions ( Figure 4D, left panel). This is in contrast to both yfp + expression from the 2 kb pho1 + pryfp + construct and endogenous expression of pho1 + during Pi-starvation which exhibit >10-fold induction in Pi limitation. Thus, the Pho7 UAS is necessary but not sufficient for Pho7-dependent transcriptional activation during Pi-starvation.
However, the 280 bp pho1 + pr-yfp + construct is not capable of relieving Csk1 repression in Pi-starvation conditions -expression is not induced in response to Pi limitation.
Trimming beyond the Pho7 UAS results in transcriptionally inactive promoters in all backgrounds tested ( Figure 4D, columns 4-6, 180 bp pho1 + pr-yfp + ). We conclude that there must be promoter elements present in the region between -2 kb and -280 bp in the pho1 + promoter that act as a Pi-sensor: (1) preventing partial Pho7-dependent activation in high-Pi conditions; and (2) are important for Csk1 de-repression during Pistarvation.
Our FACS and ChIP-Seq results lead us to the following model for Pho7 and Csk1 regulation at the pho1 + promoter. In high-Pi conditions some Pho7 is bound to the UAS in the pho1 + promoter. Pho7 in this state drives basal expression of pho1 + . Csk1, through an interaction with either Pho7 or elements near the UAS (directly or indirectly), prevents the full activation of pho1 + expression. The upstream Pi-sensor in the promoter ensures that Csk1 remains repressive in these conditions through an as yet unspecified mechanism. During Pi starvation, Csk1 repression is relieved and additional Pho7 is recruited to the pho1 + promoter, driving maximal expression ( Figure 4E). Investigating the promoter elements and transcription factors that comprise the Pi-sensor -as well as the use of this promoter structure at additional Pho7-dependent and Pi-starvation inducible promoters -is an exciting area for future research.

Pho7 Regulates Gene Expression in Response to Multiple Stress Conditions
During our expression analysis we noticed a set of genes with decreased expression in high phosphate conditions in a pho7Δ background that are not induced during phosphate starvation (Additional File 7). Additionally, Pho7 is bound to the promoters of a number of these genes in the ChIP-Seq analysis (Additional File 5).
These observations raise the following question: is Pho7 dedicated solely to the phosphate starvation pathway, like Pho4, or does it play a broader role in the stress response?
To answer this question, we identified the genes from our microarray analysis that display an increase in expression between the wild-type and pho7Δ backgrounds (in either +Pi or -Pi conditions; thresholds described in materials and methods) and asked whether their promoters contained a significant peak of bound Pho7 within 800 bp of the start codon -there are 63 genes that meet these requirements. This gene set was then processed through the Gene Ontology Tools: Term Enrichment algorithm [34], allowing us to look for biological functions over-represented in our gene set. For a GO term to be considered enriched, we required at least three unique gene products be included in the term family and the enrichment must have a p-value ≤ 0.01. 32 of the original 63 genes meet these requirements and are classified by the highest parent term available ( Figure 5A). 15 of these 32 genes are involved in transmembrane transport, including 3 involved in phosphate ion transport. Characterizing the transmembrane category in further depth we find genes involved in iron, copper, and zinc transport ( Figure 5B). A complete listing of GO Terms and the enriched gene set can be found in Additional File 9.
From the 32 identified genes, we utilized gpd1 + , hxk2 + , fio1 + , and ctr4 + expression as proxies for Pho7-mediated transcriptional induction in various stress conditions. Gpd1 is a glycerol-3-phosphate dehydrogenase that synthesizes glycerol and is essential for survival during osmotic stress [35]. During osmotic stress, glycerol pools increase, protecting the cell. Hxk2 is a hexokinase that plays a role in regulating alternative carbon utilization when glucose sources are limited [36]. It is maximally induced in response to a switch from glucose to glycerol as a carbon source. Fio1, in conjunction with Fip1, comprises the oxidase-permease iron transport system responsible for harvesting iron in depleted conditions [37]. During iron repletion fio1 + is repressed by the activity of Fep1, an iron sensing transcription factor [38]. During iron starvation fio1 + is de-repressed and induced ~70-fold. Finally, Ctr4 is a high-affinity copper transporter that is induced in copper depletion conditions by Cuf1, a copper sensing transcription factor [39]. We designed RT-qPCR primer sets for each of these genes and measured their expression as a function of osmotic, iron, copper, and carbon utilization stress in both pho7 + and pho7Δ backgrounds (see materials and methods).
As previously demonstrated [24], loss of Pho7 completely abrogates induction of pho1 + in no-Pi medium ( Figure 5C, first panel). For each of the additional stresses tested, the loss of Pho7 causes a significant decrease (p-value ≤ 0.05) in the maximal induction of the target gene ( Figure 5C). The Pho7 dependence of these genes varies, with some (fio1 + , ctr4 + ) showing a relatively minor Pho7 component, while others (gpd1 + , hxk2 + ) appear fully dependent on Pho7 to reach a peak level of induction in stress. None of the Pho7-regulated responses were as dramatic as that observed for Pi starvation, which may indicate that Pho7 plays a more subtle role in coordinating expression at Pi-independent loci. Given this subtle response, the only recent availability of a deletion collection for S. pombe [40], and the fact that Pho7 function was unverified until recently, it is not surprising that the more general role for Pho7 has not been previously observed.
It remains possible that the stress response effects we observe are artifacts limited to only the few genes we studied using RT-qPCR. Using microarray analysis with RNA collected from pho7 + or pho7Δ cells grown in non-stress and stress conditions, we examined the stress response mediated by pho7 + in the above conditions. We find that of 274 genes induced in the stress conditions studied, 44 genes are pho7 + dependent (7 genes are induced in multiple stress conditions) ( Figure 5D and Additional File 10). pho7 + is responsible for coordinating between 11.7-23.5% of the total stress response in each condition (for comparison, pho7 + coordinates 21.7% of the Pi starvation response using these thresholds). In each stress we find enrichment of distinct GO terms. The set of iron-responsive, pho7 + -regulated genes is significantly enriched for the biological process of iron assimilation (fio1 + , fip1 + , str3 + , sib2 + , p-value = 5.2e-09), the set of copper-responsive, pho7 + -regulated genes is enriched for copper ion transport (ctr4 + and ctr5 + , p-value = 5.6e-05), and the osmotic shock-responsive, pho7 + -regulated gene set is enriched for metal ion transport (zrt1 + and ctr4 + , p-value = 7.8e-04). With the carbon switch, pho7 + -regulated response we see an enrichment of genes responsible for small molecule metabolism as well as conjugation. The reasoning for the conjugation process enrichment or why pho7 + would be involved is unclear.
Overall, the biological processes are generally linked with transmembrane transport, suggesting that pho7 + is responsible for coordinating the correct transport of nutrients required for each stress response. Unlike the system in S. cerevisiae, where the central Pi-starvation regulator is tightly linked with the PHO response, the pho7 + based system in S. pombe functions differentially in a number of stress response networks.

CONCLUSIONS
In this study we have defined and characterized the gene regulatory network in S. pombe responsible for coordinating the response to inorganic phosphate starvation.
There are two distinct temporal responses in the PHO pathway in S. pombe: a fast response concerned with immediately harvesting inorganic phosphate from the environment and transporting it into the cell, and a slower one associated with a general stress response. Within the fast response we define a core PHO regulon comprised of the pho1 + , SPBC8E4.01c, and SPBC1271.09 genes whose induction in response to phosphate starvation, and regulatory behavior, has been conserved between S. pombe and S. cerevisiae.
Unlike the PHO response in S. cerevisiae, however, the positive regulator in S.
pombe is bound to the pho1 + and SPBC8E4.01c promoters irrespective of external phosphate availability. Using our Pho7-TAP ChIP-Seq dataset, we identified a single Pho7 binding region in the pho1 + promoter. This region -located between nucleotides -265 and -245 -serves as a Pho7-dependent UAS, which is necessary for transcriptional activation of pho1 + . The interaction of Pho7 with the UAS leads to a basal expression of the secreted acid phosphatase in high-Pi conditions. The 1.8 kb region preceding the UAS in the pho1 + promoter is required for Pi sensing, coordinating the activation by Pho7 and repression by Csk1 based on Pi availability. Csk1 prevents full activation of Pho7 during phosphate replete conditions and repression is maintained even in the minimal UAS construct (280 bp pho1 + pr-yfp + ). During phosphate starvation this inhibition is relieved (through an unknown mechanism) and additional Pho7 is recruited to a number of sites throughout the genome, causing further induction of pho7 + -dependent genes (see Figure 6 for a comparison of the S. cerevisiae and S. pombe PHO pathways). In previous work, Csk1 was shown to regulate transcription by activating the positive transcription elongation factor b (P-TEFb) ortholog, Cdk9 [41].
Cdk9 coordinates transcript elongation and processing, and its full activation by Csk1 leads to an increase in CTD kinase activity [42]. How this generally positive regulatory network is switched to an inhibitory role in the PHO system remains an open question.
We were also surprised to find that Pho7 was bound throughout the genome in both high-Pi and no-Pi conditions. We had thought based on previous evidence that Pho7, like Pho4, would be specific to the PHO response. Instead we demonstrate that Pho7 binds within the promoters of additional stress responsive genes and plays a role in iron, copper, osmotic, and alternative carbon utilization stress. Each stressor elicits a different pho7 + -dependent transcriptional response, though it appears that the main regulatory role of pho7 + is coordinating stress-specific transmembrane transport. There must exist some mechanism to either direct Pho7 to the proper location for inducing the correct genes or activate Pho7 at only the appropriate locations (or some mixture of both). In S. cerevisiae, the osmotic, oxidative, and glucose limitation stress responses are mediated by the transcription factor Msn2 [43]. In normal conditions, Msn2 is phosphorylated and its entry into the nucleus is limited [44]. Different stresses elicit distinct dynamics of nuclear transport, leading to a different transcriptional outputs [29].
Given that Pho7 is bound to the genome constitutively, we do not expect that nuclear exclusion will play as large a role as it does with Msn2 regulation, but it remains possible that differential post-translational modifications are responsible for this combinatorial activation by Pho7. Pho7 may be playing a more passive role in regulation, with additional factors determining Pho7 genomic localization.
Nonetheless, we have demonstrated that within the evolutionary parallel signal transduction networks that comprise the PHO pathway there exists a core PHO transcriptional regulon. The specific mechanisms involved in regulating the PHO response in S. cerevisiae and S. pombe show remarkable flexibility. An interesting area for future research centers on the environmental factors that contributed to the development of these two parallel networks. Why is the PHO response in S. pombe under the control of a general stress transcription factor, Pho7, while S. cerevisiae has developed the phosphate starvation specific pathway for Pho4 activation? What are the environmental pressures that favor a "leaky" response in S. pombe and a tightly controlled one in S. cerevisiae? Broadly speaking, our study provides a framework for determining the fundamental requirements for regulating phosphate homeostasis in Ascomycota and the specific points in the signal transduction pathway that can be altered as conditions merit. for inorganic phosphate starvation, all starvation experiments were conducted with cells incubated in a 90%SD-10%EMM media, which has been previously described [24].

Microarray Analysis and Data Processing
Strains were grown in 90%SD-10%EMM medium containing 10 mM KH 2 PO 4 (high-Pi) at 30°C until they reached early-log phase (OD 600 =0.15-0.25). Cells were collected via filtration, washed twice, transferred to fresh media lacking Pi (no-Pi), and grown at 30°C for up to 4 hours. Immediately prior to starvation (t=0), 20 mL of cells were added to 30 mL of methanol kept at -65°C to prevent further transcription or RNA degradation. At 30, 60, 120, and 240 minutes post-starvation this process was repeated. Cells were left in methanol for 10 minutes, pelleted, washed quickly in autoclaved water, re-suspended in 750 uL of RNAlater (Ambion), and snap-frozen in liquid nitrogen. RNA was extracted using the RNeasy Mini kit (Qiagen). cDNA was generated in a reverse transcriptase reaction using 10 µg total RNA with a 1:1 mixture of oligo-dT and random hexamer primers (Operon) and a 2:3 ratio of amino-allyl-dUTP:dTTP (Sigma). Superscript II RT (Invitrogen) was added and the reaction mixture was incubated at 42°C for 2.5 hours. cDNA was purified using a PCR purification kit To determine the extent of pho7 + and csk1 + regulation within the PHO response we grew the relevant strains (DP1, DP81, DP106, DP113) as described above, with the exception that cells were split into either high-Pi or no-Pi media and grown for 2 hours prior to RNA collection. Two independent biological replicates were performed for each of the conditions tested except for the pho7 + csk1Δ/pho7Δcsk1Δ comparison in no-Pi media. For each of these arrays the two probes used to detect each ORF were averaged and treated as single data points with p-values determined using a student's ttest with a one-tailed distribution against the null hypothesis in the MATLAB software.
Thresholds were set at ≥ 1.8 log 2 fold-change to facilitate comparison with the previously characterized S. cerevisiae data set [7]. Genes passing the induction threshold also had to pass a p-value threshold of ≤ 0.10. Clustering analysis was completed using k-means clustering in the Cluster 3.0 program [47] after empirically determining the optimal number of clusters using the MATLAB bioinformatics toolbox.

Chromatin Immunoprecipitation of Pho7-TAP with High-Throughput Sequencing (ChIP-Seq)
ChIP-Seq was performed on the DP1, DP94, and DP115 strains as previously described [7,[48][49][50]. Cells were grown to early log-phase (OD 600~0 .18) in high-Pi media at 30°C and split into either 200 mL of high-Pi or no-Pi media and grown for 2 hours. To determine which of our enriched regions were actually attributable to a Pho7-TAP binding event we used a modified method from [48]. For each condition analyzed we set a lower threshold for peak discovery equal to the genomic average of reads per base. We set the upper threshold equal to the highest observed read count within the given sample. Using 380 equal increments between these two thresholds we defined peaks that were larger than 100 nucleotides and separated by at least 20 nucleotides.
Peaks were compiled at the highest threshold at which they met those standards and full peaks were required to be at least 150 nucleotides distant from the nearest neighbor. Statistical analysis comparing sample enrichment to mock enrichment was performed in MATLAB using previously described methods [48]. Peaks used for subsequent analysis had a ≥ 2-fold enrichment over the genome average and a p-value ≤ 0.005 compared to the mock sample.
To determine the likelihood that the genes determined from microarray analysis to be regulated by Pi and/or pho7 would also have promoters that contain at least one Pho7 binding site we utilized a hypergeometric test. In this case the p-value is given by: where N is the total number of genes probed by the microarray (5046), M is the number of peaks within -800 and 0 bp of any start codon in no-Pi conditions (570), n is the full set of regulated genes (22 for Pi, 29 for pho7 + , and 7 for both), and k is the set of regulated genes with at least one Pho7 peak in the promoter (13 for Pi, 16 for pho7 + , and 6 for both).
Western blotting with Pho7-TAP DP1, DP94, or DP114 Cells were grown to log-phase in high-Pi media at 30°C.
Following collection of a high-Pi sample, cells were washed three times in no-Pi media, transferred to no-Pi media, and grown for either 60 or 120 minutes. Cells were lysed by beadbeating in urea lysis buffer (20 mM Tris-HCl, pH 8.0, 50 mM Na 2 HPO 4 , 8M urea, and 1 mM PMSF). Total protein was quantified on a Nanodrop 2000 (Thermo Scientific) using a BSA standard (Bio-Rad). Equal amounts of total protein for each sample (30 ug) were subjected to separation by SDS-PAGE and transferred to nitrocellulose.
Immunoblotting was performed with rabbit IgG (1:1000, Jackson ImmunoResearch) followed by incubation with goat anti-rabbit-HRP (1:5000, Thermo Scientific). Blots were developed using the SuperSignal West Femto Chemiluminescent Substrate (Thermo Scientific) and analyzed on an Alpha Innotech Gel Imagining System.

pho1 + Promoter Deletion Analysis
Segments of the pho1 + promoter were amplified using PCR and cloned into a yfp + plasmid using homologous recombination [51] containing the selectable ura + marker, creating a pho1 + pr-yfp + fusion. Plasmids were transformed into DP18, DP109, or DP111 backgrounds using lithium acetate and polyethylene glycol 8000. Cells were selected based on their ability to grow in EMM-ura media. Cells containing the various plasmids were grown to early log-phase (OD 600~0 .18) in high-Pi media (lacking uracil) at 30°C, collected, washed twice in sterile water, and split into either high-Pi or no-Pi media (both lacking uracil). Cells were grown for 4 hours at 30°C and 100 uL of 10% buffered formalin (Sigma) was added to 900 uL culture. Fixation proceeded for 5 minutes at room temperature prior to washing: once with 0.1M potassium phosphate buffer (pH 8.5) and once with 1.2M sorbitol in 0.1M potassium phosphate buffer (pH 8.5). Cells were resuspended in 1.2 M sorbitol, 0.1M potassium phosphate buffer (pH 8.5) and incubated overnight at 4°C.
FACS counting with each sample was performed using a LSR II Analyzer (BD Biosciences). DP18, DP109, and DP111 cells lacking the yfp + expression system were used to normalize forward-scatter, side-scatter, and autofluorescence for each experiment. 50,000 cells were counted for each experimental condition tested; cells with forward-and side-scatter values between 50,000-150,000 and YFP expression ≥ mean autofluorescence were subject to further analysis. Three biological replicates were performed and the average YFP intensity for the replicates is reported ± SE.

pho7 + Regulation in Additional Stress Response Pathways
For each individual stress response, initial cultures of DP1 and DP81 were grown in 90%SD-10%EMM media containing 10 mM KH 2 PO 4 to early-log phase at 30°C. Cells Shift -1.2 M NaCl instead of 0.1M NaCl, and Carbon Switch -2% glycerol/1% ethanol (GE) instead of 2% glucose (G). Cells were grown for 2 hours and harvested as described above. Recovered RNA was converted into cDNA using the iScript cDNA synthesis kit (Bio-Rad) and subjected to RT-qPCR. Amplification of the gpd1 + , fio1 + , ctr4 + , hxk2 + , and pho1 + transcripts were measured for the three independent replicates and transcript abundance was normalized to act1 + . Shown is the average ± SE. Primers used in the RT-qPCR analysis can be found in Additional File 11.
Extracted RNA was also subjected to microarray analysis as detailed above.
Expression from pho7 + cells in replete conditions was compared to that in stress conditions for each individual stressor to determine the base set of genes that respond in each given stress. The dependence upon pho7 + was determined by comparing the levels of induction in pho7 + cells in stress to induction in pho7Δ cells in stress. Based on previous reports [27] the osmotic shift conditions were assayed 20 minutes post-shift to provide a more accurate measure of genes directly induced by osmotic pressure.
Extraction of Cy3-Cy5 fluorescence intensity was performed using the GenePix 5.1 software and normalization was completed using the MATLAB bioinformatics toolbox. At least two independent biological replicates were performed for each of the conditions tested. For each of these arrays the two probes used to detect each ORF were averaged and treated as single data points with p-values determined using a student's ttest with a one-tailed distribution against the null hypothesis in the MATLAB software.
Thresholds were set at ≥ 1.

COMPETING INTERESTS
The authors declare that they have no competing interests.

AUTHOR'S CONTRIBUTIONS
The manuscript was written by ICO. Microarrays, ChIP-Seq, and FACS quantization were performed by ICO. DDW provided S. pombe strains. MP prepared the pho1 + pryfp + constructs. Computational analysis was performed by ICO using modified methods from previously published sources [7,46,48]. EKO and DDW provided supervision throughout and EKO and DDW are the corresponding authors.

DESCRIPTION OF ADDITIONAL DATA FILES
The following data files are available with the online version of this manuscript.     Coordinates are relative to the initial ATG, with gene products identified by black bars.
Reads were normalized as previously described.   Shown are ChIP-Seq profiles for the genes identified in Figure 6. As previously described, wild-type cells containing Pho7-TAP were grown in either high-Pi (blue) or no-Pi (red) conditions and ChIP-Seq libraries were prepared from purified DNA. For comparison, the ChIP-Seq signal from mock (black) cells grown in no-Pi and csk1Δ (cyan) cells incubated in high-Pi is included. The gene product of interest is plotted based on transcript direction with the plus (+) strand above and the minus (-) strand below. Reads were normalized as described in the text and the location within the genome is plotted on the x-axis.
Table1.xlsx Genes induced following two hours (fast response, bold) or four hours (slow response) of Pi depletion were tabulated and their normalized expression levels (to t=0 min) are shown. Induction thresholds for the fast response were set at ≥ 2σ + median log 2 fold change for each time point (1.00 log 2 fold change at 120 minutes, 1.24 log 2 fold change at 240 minutes). Genes above threshold at both 120 minutes and 240 minutes poststarvation were classified as the rapid response. Genes above only the 240-minute threshold were classified as the slow response.
Table2.xlsx Genes possessing a pho7 + -dependency independent of Pi availability were crossreferenced with the Pho7-TAP peak list to identify pho7 + -regulated genes with promoters enriched by Pho7-TAP. Thresholds for induction and peak enrichment are described in materials and methods. The gene set was processed through the AmiGO toolkit [34,[52][53][54][55] and GO terms containing at least 3 genes with p-values ≤ 0.01 are displayed.
All primers utilized in this manuscript are listed with their primer ID, nucleotide sequence, and experimental purpose.