- Research article
- Open Access
Genomewide identification of pheromone-targeted transcription in fission yeast
BMC Genomics volume 7, Article number: 303 (2006)
Fission yeast cells undergo sexual differentiation in response to nitrogen starvation. In this process haploid M and P cells first mate to form diploid zygotes, which then enter meiosis and sporulate. Prior to mating, M and P cells communicate with diffusible mating pheromones that activate a signal transduction pathway in the opposite cell type. The pheromone signalling orchestrates mating and is also required for entry into meiosis.
Here we use DNA microarrays to identify genes that are induced by M-factor in P cells and by P-factor in M-cells. The use of a cyr1 genetic background allowed us to study pheromone signalling independently of nitrogen starvation. We identified a total of 163 genes that were consistently induced more than two-fold by pheromone stimulation. Gene disruption experiments demonstrated the involvement of newly discovered pheromone-induced genes in the differentiation process. We have mapped Gene Ontology (GO) categories specifically associated with pheromone induction. A direct comparison of the M- and P-factor induced expression pattern allowed us to identify cell-type specific transcripts, including three new M-specific genes and one new P-specific gene.
We found that the pheromone response was very similar in M and P cells. Surprisingly, pheromone control extended to genes fulfilling their function well beyond the point of entry into meiosis, including numerous genes required for meiotic recombination. Our results suggest that the Ste11 transcription factor is responsible for the majority of pheromone-induced transcription. Finally, most cell-type specific genes now appear to be identified in fission yeast.
The fission yeast Schizosaccharomyces pombe provides an attractive model system for studying cellular differentiation processes in response to environmental cues. It is a unicellular, haploid organism with a conventionally organized cell cycle , and it is easily subjected to genetic analysis . Furthermore, since its entire genome sequence is known , genome-wide analysis tools can be readily applied.
In response to nitrogen starvation fission yeast cells undergo sexual differentiation. In this process haploid P and M cells first mate to form dipoid zygotes , which then undergo meiosis and sporulation . Nitrogen starvation leads to the establishment of a pheromone signalling system by which M and P cells communicate with each other (see Figure 1A), causing growth arrest in the G1 phase of the cell cycle. Pairs of G1 arrested M- and P- cells subsequently agglutinate, mate, and undergo karyogamy, after which the resulting diploid zygotes enter meiosis and finally sporulate.
The HMG-box transcription factor Ste11 is a master regulator of the sexual differentiation programme . Upon starvation, Ste11 activates numerous genes required for mating and meiosis, including the mat genes specifying cellular mating type, many genes encoding components of the pheromone communication system, and mei2 that encodes the ultimate activator of meiosis. Ste11 binds to one or several Ste11-sites (AACAAAGAAA) present in the promoter regions of the genes it directly controls.
The two pheromones are small, secreted peptides that bind to and activate cell-type specific 7TM receptors on the cell surface (Mam2 in M cells and Map3 in P cells), thereby activating the receptor-coupled G-protein, Gpa1. This stimulates a MAP kinase cascade, consisting of Byr2 (a MAP3K), Byr1 (a MAP2K) and Spk1 (a MAPK). The pheromone response, which is shared by the two cell types, activates genes required at various steps throughout the differentiation process, with entry into meiosis being the most downstream event reported to rely on the signal .
Several lines of evidence suggest that Ste11 is also the transcription factor activated by pheromone signalling. Pheromone-controlled genes often have Ste11-boxes in their promoter regions, and hyper-activation of the MAP2K Byr2 causes constitutive activation of Ste11 . Spk1 can in fact directly phosphorylate Ste11 on threonines 305 and 317 in vitro, and mutational analysis of these residues in vivo is consistent with a direct activation mechanism . However, since Ste11 is absolutely required for the expression of many components in the response pathway, it is difficult to evaluate whether pheromone signalling can stimulate transcription in the absence of Ste11.
The analysis of pheromone-controlled transcription is further complicated by the nutritional control over the differentiation pathway. In order to make them respond to mating pheromones, it is necessary to starve cells for nitrogen, which in itself causes dramatic changes in the expression of many genes [9, 28]. Furthermore, the transcriptional responses to nitrogen starvation and pheromone stimulation are partially overlapping, and Ste11-controlled genes appear to be induced to varying degrees by both these environmental cues (Figure 1B). Some genes, such as ste6  are already strongly induced by nitrogen starvation and become further enhanced by pheromone (Figure 1B, panel 1). Expression of other genes like fus1  has an absolute requirement for both nitrogen starvation and pheromone signalling (Figure 1B, panel 2). Genetic analysis has shown that byr2, byr1 and ste11 are required for induction of pheromone controlled genes by nitrogen starvation , suggesting that this signal is also conveyed by the pheromone-responding MAPK pathway.
An additional layer of regulation is conferred by cellular mating type. Whilst most pheromone-controlled genes appear to be equally expressed in the two cell types, a small subset is either M-specific (Figure 1B, panel 3) or P-specific (Figure 1B, panel 4). These genes are involved in synthesis and processing of the cell-type specific extracellular components of the pheromone-communication system, i.e. the pheromones and receptors. Their expression is controlled by the allele present at the mat1 locus of the cell . In M cells the mat1-Mc gene causes activation of M-specific genes, while the mat1-Pc gene activates P-specific genes in P cells. Mat1-Mc is an HMG-box protein that appears to recruit Ste11 to a weaker version of its binding site (AACAAAGA) present in the promoter region of M-specific genes . The mechanism by which Mat1-Pc activates the P-specific genes is less clear. Genetic analysis suggests that it cooperates with the Map1 protein, a putative transcription factor of the MADS-box family [15, 16].
In the present study we use genomic microarrays to identify pheromone-induced genes in fission yeast. In order to discriminate pheromone-induction from the effect of nitrogen starvation, we utilize the fact that cyr1 strains defective in cAMP production are derepressed for mating activities, and hence respond to pheromone without being starved . We have identified classes of genes induced by M-factor in P cells and genes induced by P-factor in M-cells. As expected, these two groups of genes were highly overlapping. Moreover, a direct comparison of M and P cells allowed us to identify several new genes that are either M- or P-specific. Our results support the notion that Ste11 is the major transcription factor bringing about pheromone-stimulated transcription.
Identification of genes induced by M- and P-factor
In order to characterize the changes in transcription profiles caused by mating pheromone stimulation, we performed two different experiments; one where we added M-factor to h+cyr1 cells, and a second adding P-factor to h-cyr1 sxa2 cells. The cyr1 mutation was included in both strains in order to enable pheromone stimulation without first subjecting the cells to nitrogen starvation . Furthermore, the h- strain was inactivated in the sxa2 gene, which encodes a very potent P-factor degrading protease .
We took samples for RNA preparation after 2, 4 and 6 hours from both cultures. Total RNA from each time point was hybridised to genomic arrays using the 0 hour time point from the same mating type as a reference. The complete data set is available from the Gene Expression Omnibus database (see Materials and Methods). Analysis of the data showed that pheromone treatment caused dramatic induction of a relatively large number of genes. K-means cluster analysis of these genes revealed that most genes were already induced after 2 hours of pheromone stimulation, and that their expression remained high for the duration of the experiment (data not shown), suggesting that pheromone induction occurred rapidly. For each pheromone, we therefore identified genes that were consistently induced more than twofold at all three time-points after pheromone addition. Using these criteria, we identified 120 M-factor induced genes and 93 P-factor induced genes [see Additional file 1]. A direct comparison revealed that most pheromone-induced genes were expressed at similar levels in untreated h+ and h- cells [see Additional file 1].
Comparison of the sets of M-factor and P-factor induced gene lists showed that they were highly overlapping, containing many genes in common. This is consistent with the fact that the pheromone response pathway is common to both mating types, and most of the proteins required in the differentiation pathway perform similar functions in the two cell types. To analyse the similarity of the responses further we compared the M-factor and P-factor induced changes for both sets of induced genes. Figure 2 shows that the fold change values for M- and P- treatment are highly correlated for most genes. Most of the minority of genes that are less well correlated still show a tendency to induction under both conditions.
The list identified by our array experiments contained genes such as bqt1, fus1 and rep1 previously demonstrated to be pheromone controlled, but also many new genes (see Table 1). To investigate the function of the pheromone induced genes further, we identified Gene Ontolology (GO) categories that were significantly associated with M-factor and P-factor induced genes (Figure 3).
As expected, genes involved in conjugation were enriched in both M-factor and P-factor induced genes. In addition, we found genes such as spk1,ste4 and ste11 encoding proteins directly involved in response to pheromone.
Other functions, specifically involved in the morphological alterations observed during mating, were also enriched among the pheromone-induced genes. These included glucan catabolism, cytoskeleton-dependent intracellular transport, and microtubule-mediated nuclear migration and nucleus localization. Presumably, the induction by pheromone of a specific subset of genes involved in these processes redirects them to accommodate their specific behaviour during differentiation.
The mat1-Pm and mat1-Mm genes required for entry into meiosis are the most downstream functions so far reported to be pheromone controlled . It was therefore of interest to find gene categories performing specific meiotic functions that lie further downstream (regulation of meiotic transcription, meiosis I and meiotic chromosome segregation) among the up-regulated genes. Genes in these categories include the meiotic master regulator mei2 and the meiotic forkhead transcription factor mei4 but our results also identify new potential regulators such as the putative forkhead transcription factor (SPAP14E8.02). Finally, a relatively large group of genes belonging to the categories meiotic recombination and meiotic gene conversion were found to be pheromone induced. These include dmc1, eme1, rad25, rad32, rec10, rec12, rec25, rec27, rti1 and swi5.
Validation and functional characterization
To confirm that the newly identified genes were actually pheromone controlled, we performed real-time RCR analysis on several of them (Figure 4 and data not shown). In order to do this, we did the pheromone stimulation with cyr1+ cells that were also starved for nitrogen. This allowed us not only to validate the results in a different genetic background, but also to evaluate the inducibility by nitrogen starvation of the pheromone-induced genes. As an example, Figure 4A shows the analysis of the gene SPBC4.01, which was strongly induced by both pheromones in the array experiment. The real-time PCR data showed that this gene was already induced by nitrogen starvation, and confirmed that its expression was further elevated by pheromone addition in both mating types, leading to a more than 50 fold induction in total.
SPBC4.01 encodes a sequence orphan with no homologs in other organisms. In order to examine whether the gene was required for sexual differentiation, we disrupted the gene. As judged by iodine starvation, elimination of SPBC4.01 indeed caused a dramatic defect in the ability of cells to undergo sporulation (Figure 4B), as both the rate of entry into differentiation and the final level of sporulation were strongly reduced. Microscopic examination revealed that the disrupted strain was defective in both cell fusion and spore formation (Figure 4C). In contrast, loss of SPBC4.01 did not seem to affect mitotic cell growth severely (data not shown). In conclusion, the array data allowed us to identify novel pheromone-induced genes involved in sexual differentiation. Furthermore, the real-time PCR confirmed their regulation by pheromone, and demonstrated that such genes could also be induced by nitrogen starvation alone.
Pheromone controlled genes have an elevated occurrence of putative Ste11 binding sites
As outlined above, several lines of evidence suggest that Ste11 is the transcription factor activated in response to pheromone signalling. We therefore asked whether the pheromone-controlled genes identified by our array study had Ste11 binding sites in their promoter regions. The number of putative Ste11 binding sites as well as the sequence and position of the best putative site are provided in Additional file 1. Figure 5 shows that Ste11 binding sites within 500 bp upstream of the ATG start codon are about four-fold more frequent in pheromone-regulated genes than in other genes. Ste11 sites are also found further upstream but their frequency gradually drops towards a background level of about 2%. This background level is also seen for the remaining genes in the genome that are not included in our pheromone-induced set. For these genes the background level is fairly constant showing only a slight increase in the -500 to -1 region. Thus the increase in the frequency of Ste11 sites in the pheromone-induced genes appears to be specific. We conclude that the frequency of Ste11 binding sites is highly enriched in the promoters of pheromone-regulated genes, consistent with the notion that Ste11 is the major transcriptional regulator involved in their induction.
Identification of cell-type specific genes
As mentioned in the introduction, the M- and P-specific genes controlled by the mat1 locus are also induced by pheromone. We therefore sought to identify cell-type specific genes by looking for transcripts that were differentially induced by the two mating pheromones. We hybridised the two RNA preparations obtained after 6 hours of M-factor and P-factor stimulation to the same array, and selected transcripts that had a more than 2 fold difference in expression levels. By these criteria we identified 8- M-specific and 4 P-specific genes (see below). We confirmed that these genes were selectively induced by either M-factor or P-factor in the array time course described above (data not shown). In Table 2 we have summarized current knowledge on functions and expression patterns of the cell-type specific genes in fission yeast.
The M-specific transcripts included 5 genes previously described to encode M-specific mating functions (mat1-Mc, sxa2, mam1, mam2 and mam4); the remaining known M-specific genes (mfm1- 3 and mat1-Mm) are all very short and are not present on the arrays we used. Three new genes (mam3, cwp1 and 11H11.03) were identified, and we confirmed all three to be M-specific by real-time PCR on RNA from cyr1+ control cells (Figure 6A).
The first gene, mam3, is annotated to encode the M cell agglutinin , and hence, its expression was expected to be M-cell specific. The second gene, cwp1 encodes the α-subunit of the enzyme geranyl-gernyl transferase I (GGTase I) , and its induction by M-factor presumably reflects its role in farnesylation of M-factor . Analogous to the cwg2 gene encoding the β-subunit of GGTase I , we presume that cwp1 is an essential gene, and we have not evaluated its involvement in M-factor production. Consistent with a role of cwp1 in other cellular processes, its expression does not have an absolute requirement for M-specific pheromone stimulation (Figure 6A).
The third gene, 11H11.03, encodes a hypothetical protein that is conserved in eukaryotes. Interestingly, this gene is divergently transcribed from the same promoter region as the M-specific gene mam2, and we speculate that the two genes are co-regulated. This gene is also expressed at a relatively high level in P-cells (Figure 6A), and we have not analysed its involvement in mating activities of M-cells.
With these findings, the total number of known M-specific genes is 12. In Table 3 we have aligned their promoter sequences upstream of their Ste11-box. Except for mam4 all the M-specific genes contain imperfect Ste11-boxes (AACAAAG) to which Ste11 only binds in the presence of Mat1-Mc . mam4 contains a perfect Ste11-box (AACAAAGAAA) to which Ste11 can bind on its own and consistent with this mam4 is not strictly M-specific, as it is expressed in both cell types, albeit 2–3 fold higher in M cells [, Additional file 1]. Similar to the strict M-specific genes, mam4 also harbours an M-box (a Mat1-Mc binding site) 14–26 bp upstream of the Ste11-box (Table 3), which presumably boosts its expression in M cells.
Apart from the two mating-type genes mat1-Pc and mat1-Pm (which are not present on our arrays) only two P-specific genes have been described, namely map2 and map3, encoding, respectively P-factor and the receptor for M-factor. Both these genes were selectively induced by M-factor in our arrays (Figure 6B, Table 2 and data not shown). The map4 gene encoding P-cell specific agglutinin  was not identified by our array approach, but real-time PCR analysis showed that its expression indeed was P-specific (Figure 6B). The two additional P-specific genes identified in the arrays (SPAC1565.03 and SPAC3F10.09), were analysed further. SPAC1565.03 was confirmed to be P-specific by real-time PCR analysis (Figure 6B), and hence represents a new bona fide P-specific gene. However, disruption of SPAC1565.03 had no obvious effect on mating of P-cells (data not shown), leaving the reason for its P-specific regulation unexplained.
The other P-specific gene, SPAC3F10.09, represented a puzzle. Its sequence suggests that it encodes an enzyme homologous to S. cerevisiae HIS6, and hence it is likely to be involved in histidine biosynthesis; yet real-time PCR analysis confirmed that its expression was highly induced by pheromone in a P-cell specific manner (Figure 7A). We therefore sought alternative explanations for its apparent P-cell specific expression pattern. We noticed that this gene is partially overlapping with the P-specific gene map3 (see Figure 7B). The two genes are convergently transcribed, and their reading frames overlap by 27 codons. Hence, we speculated that the array spots for SPAC3F10.09 might have picked up the map3-specific signal from the opposite strand. In order to test this hypothesis, we did strand-specific RT-PCR with the same two SPAC3F10.09 primers that were used to generate the PCR product constituting the array spots (Figure 7C). Indeed, the P-specific, pheromone-induced transcript originates from the map3-specific strand. The SPAC3F10.09 specific transcript from the opposite strand gives rise to a somewhat smaller PCR product due to the splicing out of the 66 bases long third intron in the gene (see Figure 7C). Curiously, this transcript seemed to be down-regulated when map3 was induced by pheromone, suggesting that strong transcription in the opposite direction can negatively affect gene expression. Consistent with the array data, we mapped the 3'-end of the map3 mRNA to the sequence 5'-ACATCCGAGA-3' in the second intron of SPAC3F10.09 (see Figure 7B).
In fission yeast, nitrogen starvation causes induction of gene expression by at least three different signalling pathways. First, it leads to a modest reduction in cellular cAMP, which in turn down regulates the glucose sensing Pka1 protein kinase that directly represses transcription of genes such as fbp1 and ste11 . Secondly, nitrogen starvation activates the Spc1/Sty1 stress response pathway , which causes induction of the atf1 transcription factor that also controls ste11 . Finally, genetic analysis suggests that the pheromone-responding MAPK pathway conveys the signal of nitrogen starvation to pheromone controlled promoters .
In the present study we have used a cyr1 genetic background to separate pheromone-induced transcription from the general effects of nitrogen starvation. Consequently, we found that pheromone stimulation of cyr1 cells caused a much less dramatic change in expression profiles than nitrogen starvation of wild-type cells [9, 28]. By this approach we have identified a total of 163 genes that were more than two-fold induced by pheromone stimulation, including several genes previously reported to be pheromone controlled.
Among the pheromone-induced genes we found the cAMP degrading phosphodiesterase Cgs2, suggesting that pheromone signalling enforces the reduction in cAMP caused by nitrogen starvation . This gene was also shown to be induced by activation of the Spc1/Sty1 stress response pathway , suggesting that cgs2 constitutes a junction for cross talk between different starvation sensing pathways. When cAMP levels drop, the APC/C E3 ubiquitin ligase becomes activated , and this causes G1 arrest via degradation of the S-phase activating cyclin Cig2 . Consistent with an increased APC/C activity during differentiation, we found apc2 encoding a subunit of APC/C among the pheromone-induced genes.
Our identification of components in the pheromone signal transduction pathway (spk1,ste4 and ste11) suggests that initial pheromone stimulation triggers a positive feedback loop enforcing the response. Furthermore, the whi5 gene, whose counterpart in S. cerevisiae encodes an inhibitor of transcriptional activities required for vegetative S-phase entry [33, 34], was among the induced genes, suggesting that pheromone signalling inhibits vegetative growth. Taken together these observations underscore the autocatalytic and irreversible nature of the developmental switch.
On the other hand, functions that dampen the pheromone response were also identified. Thus, the msa1 gene encoding an RNA-binding protein that negatively regulates differentiation specific transcripts , and rgs1 coding for a GTPase that down-regulates the receptor-coupled G-protein Gpa1  were both pheromone induced, stressing the importance of a tight modulation of the response during differentiation.
The results presented here are consistent with Ste11 being the transcription factor activated by pheromone signalling. The pheromone-regulated genes in Figure 5 contain perfect matches to a matrix constructed from known sites and they together account for about 30% of the pheromone-induced transcripts. The remaining pheromone-regulated genes almost invariably contain more diverged derivatives of the Ste11-binding site, often in multiple copies. It is thus likely that Ste11 regulates the majority of pheromone-regulated genes, perhaps in some cases via co-operative binding to a number of relatively weak binding sites as previously reported . It is, however, likely that some pheromone-regulated genes are indirectly induced by other transcription factors, either alone or in combination with Ste11. Indeed, several transcription factors such as rep1 and mei4 are found within the group of pheromone-induced genes.
The isp6 gene, which we found to be pheromone induced, was actually reported to be induced by nitrogen starvation independently of ste11 , and it does not have an obvious Ste11-box in its promoter region [see Additional file 1]. The cAMP pathway negatively regulates isp6 transcription , and presumably the pheromone effect on this gene is caused indirectly by cross talk with other signalling pathways as discussed above.
The cig2 gene, encoding an S-phase activating cyclin, is specifically expressed in the G1 phase of the cell cycle , and hence its induction by pheromone could be an indirect consequence of pheromone-mediated G1 arrest However, the cig2 promoter contains a perfect Ste11-box [see Additional file 1], suggesting that Ste11 and the pheromone pathway may also directly control this gene. The regulation of cig2 represents a paradox: Transcription of the gene is on the one hand induced by pheromone, yet the same signal also mediates proteolysis of the Cig2 protein . Presumably this regulation reflects the fact that Cig2 needs to be down-regulated in order to activate the differentiation pathway, but then Cig2 is needed again soon after for activation of pre-meiotic S-phase .
By comparing the effects of M- and P-factor stimulation we were able to identify P-and M-specific genes. Given the extensive hunt over many years for cell-type specific mutants, it was no surprise that we only found very few new genes, and an important conclusion of this work is that the list in Table 2 can be regarded, as close to being complete. The larger number of M-specific genes largely reflects the more elaborate synthesis of M-factor relative to P-factor.
Our results confirm the proposed mechanism of M-specific gene activation , where the Mat1-Mc protein both recruits Ste11 to a weak binding site and also binds to an M-box sequence situated approximately 1.5 or 2.5 helical turns upstream of the Ste11-box (Table 3). The control of P-specific transcription is still an enigma that awaits functional promoter dissection studies to be resolved.
Our observation that the P-specific transcript map3 extends into the neighbouring gene SPAC3F10.09 demonstrates the importance of mapping the genome-wide location of transcripts and the use of this information in the design of future DNA microarrays. Using a different array system, SPAC3F10.09 was reported to be induced in meiosis , presumably also reflecting an erroneous detection of map3 transcripts transcribed from the opposite strand. Convergent overlapping reading frames are relatively rare in the fission yeast genome, but our finding that the map3 transcript extends more than 500 bases downstream of the reading frame suggests that convergent transcripts may overlap even when the reading frames do not.
We also found 92 genes that were more than twofold down regulated by both pheromones (data nor shown). GO categories representing nucleotide and amino acid metabolism were significantly over-represented among these genes, suggesting that specific mechanism exist either to repress these genes or degrade their transcripts during differentiation.
Pheromone communication in fungi was originally considered to be a mechanism by which cells of opposite mating type identified each other, but subsequent studies have demonstrated that it controls many other functions than mere cell recognition. In fission yeast pheromone signalling brings about G1 cell-cycle arrest prior to mating, cell elongation towards the partner cell, and is also involved in cell fusion and karyogamy. Finally, it has been known for a decade that entry into meiosis also requires activation of the pheromone pathway. This is because the two meiotic subfunctions mat1-Pm and mat1-Mm are pheromone controlled; if they are expressed from a heterologous promoter, meiosis is readily activated in the absence of pheromone signalling . Our results in the present study show that – although dispensable for completion of meiosis – pheromone control extends beyond the point of meiotic entry. Thus, several genes required for specific meiotic sub-functions in particular for meiotic recombination turned out being pheromone induced. Consistent with this, expression of the dmc1 gene encoding a RecA like helicase was previously shown to be dependent on Ste11 . However, we note that many of these genes do not have optimal Ste11 binding sites in their promoters [see Additional file 1], suggesting that their induction by pheromone may be indirect. Nevertheless, it is evident from our array results that pheromone stimulation in itself can bring about a robust induction of these genes.
Previously, the bqt1 gene mediating telomere clustering prior to meiotic horsetail movement was shown to be pheromone induced . These nuclear movements are already initiated in haploid cells challenged with pheromone . We found that dhc1 and dic1, encoding dynein heavy and intermediate chains respectively were both pheromone induced. Dynein is required for the horsetail movements, and a null mutant in dhc1 significantly reduces meiotic recombination . Furthermore, the dynein-associated Ssm4 protein is also pheromone regulated  and this study). Taken together, these observations suggest that pheromone signalling orchestrates the entire process of pre-meiotic chromosome movement and associated recombination. It will be interesting to learn whether oscillations in pheromone signalling play a dynamic role in bringing about these processes.
We have identified 163 pheromone induced genes in fission yeast. Most genes are induced by pheromone in both cell types, and analysis of GO annotations associated with pheromone induction showed that the genes were involved in many aspects of sexual differentiation, including conjugation, and morphological processes associated with mating. Surprisingly, many meiotic gene functions were also pheromone controlled. Our results confirmed that the Ste11 transcription factor controls many pheromone-induced genes. Finally, most cell-type specific genes now appear to be identified in fission yeast.
Strains, growth conditions and RNA preparation
Standard genetic procedures were carried out as described previously . For array experiments we used the strains Eg794 (genotype: h+ Δmat2,::LEU2 cyr1::ura4+ ura4-D18, ) and Eg796 (genotype: h- Δmat2,3::LEU2 cyr1::ura4+ ura4-D18 sxa2, this study). Cells were grown at 30°C in MSL  to a density of 2 × 106 cells/ml, and chemically synthesized M-factor (made by Schafer-N, Copenhagen, Denmark) or P-factor (made by Alto Bioscience, Birmingham, UK) was added to a final concentration of 1 μg/ml. Samples were taken for RNA preparation after 2, 4 and 6 hours. For validation of the results in cyr1+ genetic background cells, we grew the strains Eg544 (genotype: h- Δmat2,3::LEU2) and Eg545 (genotype: h+ Δmat2,3::LEU2) at 30°C in MSL to a density of 2 × 106 cells/ml. Cells were collected and resuspended in, respectively MSL, MSL – nitrogen and MSL – nitrogen + the appropriate pheromone (same concentration as above). After 5 hours of incubation, cells were harvested for preparation of RNA.
Total RNA for microarray analysis was prepared as described in . Microarray probe labeling and hybridization was performed as previously described , using arrays from Eurogentic SA containing PCR-generated probes in duplicate for 4927 nuclear genes. Slides were scanned using a BIO-RAD Vers-Array ChipReader and quantified with Imagene 4.2 software. Data were normalized by Lowess normalization using GeneSpring software (Agilent Technologies). Temporal expression patterns of genes in response to pheromone treatment were studied using K-means clustering using GeneSpring software (Agilent Technologies). Larger clusters of genes manifesting patterns that were robustly identified independently of cluster number where regarded as significant. Pheromone regulated genes were selected by a variety of criteria including consistent regulation of different samples and replicates by = 2-fold as well as the magnitude of mean fold changes, their significance (one-samole T-test) and variation (SD).
GoMiner  was used to identify gene ontology (GO) categories that were significantly (p ≤ 0.05) enriched in sets of both M-factor and P-factor regulated genes in relation to the frequency of their occurrence in the whole S. pombe genome. Microarray data generated in this study are available at Gene Expression Omnibus  under accession number GSE5323.
Prediction of DNA binding sites
The Patser algorithm, available via the Regulatory Site Analysis Tools package  was used to search for putative binding sites for the Ste11 transcription factor upstream of the coding regions of genes. The frequency matrix for nucleotides at each position (Table 4) was constructed using sequences found in known Ste11 target genes .
The significance of the elevated Ste11-site frequency in pheromone-induced genes was tested in relation to the mean frequency observed in 10 random gene sets of equal size. A one-sample t-test was used to determine whether the mean Ste11-site frequency in the random sets was significantly different from that observed in the sets of pheromone-induced genes. The number of putative Ste11 sites in intergenic regions was also estimated by searching for strings of nucleotides matching the optimal Ste11 site (the inverse sequence of AACAAAGAAA) allowing for 0, 1 or 2 mismatches.
Real-time and RT PCR analysis
Total RNA was extracted using hot acidic phenol as described above. To remove DNA the RNA samples were digested with DNaseI prior to RT-PCR and real-time PCR analyses. The gene-specific primer sets used for PCR have been described previously . Quantitative real-time PCR reactions were carried out in 96-wells optical reaction plates; using 50 ng DNaseI treated RNA and the QuantiTect SYBR-green RT-PCR kit (Qiagen) supplemented with 10 nM fluorescein. Real-time RT-PCR was developed using the iCycler (BIO-RAD) by single-step amplification. The thermal cycling conditions comprised reverse transcription at 50°C for 30 min and an initial denaturation at 95°C for 15 min, followed by 40 cycles of 95°C for 15 s, 50°C for 30 s and 72°C for 60 s. Each sample was assayed in triplicate. To minimize errors arising from variations in the amount of starting RNA among samples, amplification of act1 mRNA served as an internal reference. Standard curves for each relevant gene and for act1 were produced by performing serial two-fold dilutions of control RNA (a 1:1 mixture of RNA from pheromone-stimulated cells of opposite mating types). For each experimental sample, the amount of specific transcript was normalized to the amount of act1 transcript measured in the same sample.
RT-PCR was carried out using 50 ng DNaseI-treated RNA and the Qiagen OneStep RT-PCR kit. The thermal cycling conditions comprised reverse transcription at 50°C for 30 min and an initial denaturation at 95°C for 15 min, followed by 20–25 cycles of 95°C for 15 s, 50°C for 30 s and 72°C for 60 s. For strand-specific RT-PCR, the reactions were split in two and one of the primers (either sense or anti-sense) was not added until after inactivation of reverse transcriptase.
Functional characterisation of genes
The open reading frame of SPBC4.01 was replaced by the kmx-cassette by means of the PCR-based one-step gene disruption procedure, using the oligos OBP421 and ONP422 . The strain deleted for SPAC1565.03 was obtained from Bioneer Corporation, Korea. Sporulation was evaluated by iodine staining as previously described .
The 3' end of the map3 mRNA was determined by 3' RACE (Rapid Amplification of cDNA Ends). In brief, RNA extracted from pheromone stimulated P cells was reverse transcribed using the oligonucleotide SKP62. Next, map3 was amplified by PCR using the oligonucleotides SKP63 and ONP441, and the 3' end was identified by sequencing.
ONP421: TTATAATAGACTTACATTCCATAATAGTTCTTCTACTCCCTAAACTAATT TTGCTAAAGT AAACAATTTAAAATTACGGTTTAAGACGGATCCCCGGTTAATTAA
ONP422: 5' GGTTTATGATTCATTAAAAGCTAATTCTCAATACATTTTTGCTCAAATTAGTAATATGGTTTACAAAATCTTATTTATGTTTGTGAATTCGAGCTCGTTTAAAC
Forsburg SL, Nurse P: Cell cycle regulation in the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. Annu Rev Cell Biol. 1991, 7: 227-256. 10.1146/annurev.cb.07.110191.001303.
Moreno S, Klar A, Nurse P: Molecular genetic analysis of fission yeast Schizosaccharomyces pombe. Methods Enzymol. 1991, 194: 795-823.
Wood V, Gwilliam R, Rajandream MA, Lyne M, Lyne R, Stewart A, Sgouros J, Peat N, Hayles J, Baker S, et al: The genome sequence of Schizosaccharomyces pombe . Nature. 2002, 415: 871-880. 10.1038/nature724.
Nielsen O: Mating-type control and differentiation. The Molecular Biology of Schizosaccharomyces pombe. 2004, Egel R: Springer, 281-296.
Yamamoto M: Initiation of meiosis. The Molecular Biology of Schizosaccharomyces pombe. 2004, Egel R: Springer, 297-309.
Sugimoto A, Iino Y, Maeda T, Watanabe Y, Yamamoto M: Schizosaccharomyces pombe ste11+ encodes a transcription factor with an HMG motif that is a critical regulator of sexual development. Genes Dev. 1991, 5: 1990-1999.
Willer M, Hoffmann L, Styrkarsdottir U, Egel R, Davey J, Nielsen O: Two-step activation of meiosis by the mat1 locus in Schizosaccharomyces pombe. Mol Cell Biol. 1995, 15: 4964-4970.
Kjaerulff S, Lautrup-Larsen I, Truelsen S, Pedersen M, Nielsen O: Constitutive activation of the fission yeast pheromone-responsive pathway induces ectopic meiosis and reveals ste11 as a mitogen-activated protein kinase target. Mol Cell Biol. 2005, 25: 2045-2059. 10.1128/MCB.25.5.2045-2059.2005.
Mata J, Lyne R, Burns G, Bahler J: The transcriptional program of meiosis and sporulation in fission yeast. Nat Genet. 2002, 32: 143-147. 10.1038/ng951.
Hughes DA, Yabana N, Yamamoto M: Transcriptional regulation of a Ras nucleotide-exchange factor gene by extracellular signals in fission yeast. J Cell Sci. 1994, 107: 3635-3642.
Petersen J, Weilguny D, Egel R, Nielsen O: Characterization of fus1 of Schizosaccharomyces pombe: a developmentally controlled function needed for conjugation. Mol Cell Biol. 1995, 15: 3697-3707.
Kjaerulff S, Davey J, Nielsen O: Analysis of the structural genes encoding M-factor in the fission yeast Schizosaccharomyces pombe: identification of a third gene, mfm3. Mol Cell Biol. 1994, 14: 3895-3905.
Kelly M, Burke J, Smith M, Klar A, Beach D: Four mating-type genes control sexual differentiation in the fission yeast. EMBO J. 1988, 7: 1537-1547.
Kjaerulff S, Dooijes D, Clevers H, Nielsen O: Cell differentiation by interaction of two HMG-box proteins: Mat1-Mc activates M cell-specific genes in S.pombe by recruiting the ubiquitous transcription factor Ste11 to weak binding sites. EMBO J. 1997, 16: 4021-4033. 10.1093/emboj/16.13.4021.
Nielsen O, Friis T, Kjaerulff S: The Schizosaccharomyces pombe map1 gene encodes an SRF/MCM1-related protein required for P-cell specific gene expression. Mol Gen Genet. 1996, 253: 387-392.
Yabana N, Yamamoto M: Schizosaccharomyces pombe map1+ encodes a MADS-box-family protein required for cell-type-specific gene expression. Mol Cell Biol. 1996, 16: 3420-3428.
Maeda T, Mochizuki N, Yamamoto M: Adenylyl cyclase is dispensable for vegetative cell growth in the fission yeast Schizosaccharomyces pombe . Proc Natl Acad Sci USA. 1990, 87: 7814-7818. 10.1073/pnas.87.20.7814.
Davey J, Nielsen O: Mutations in cyr1 and pat1 reveal pheromone-induced G1 arrest in the fission yeast Schizosaccharomyces pombe. Curr Genet. 1994, 26: 105-112. 10.1007/BF00313796.
Imai Y, Yamamoto M: Schizosaccharomyces pombe sxa1+ and sxa2+ encode putative proteases involved in the mating response. Mol Cell Biol. 1992, 12: 1827-1834.
Arellano M, Coll PM, Yang W, Duran A, Tamanoi F, Perez P: Characterization of the geranylgeranyl transferase type I from Schizosaccharomyces pombe. Mol Microbiol. 1998, 29: 1357-1367. 10.1046/j.1365-2958.1998.01009.x.
Davey J: Mating pheromones of the fission yeast Schizosaccharomyces pombe: purification and structural characterization of M-factor and isolation and analysis of two genes encoding the pheromone. EMBO J. 1992, 11: 951-960.
Imai Y, Davey J, Kawagishi-Kobayashi M, Yamamoto M: Genes encoding farnesyl cysteine carboxyl methyltransferase in Schizosaccharomyces pombe and Xenopus laevis. Mol Cell Biol. 1997, 17: 1543-1551.
Sharifmoghadam MR, Bustos-Sanmamed P, Valdivieso MH: The fission yeast Map4 protein is a novel adhesin required for mating. FEBS Lett. 2006, 580: 4457-4462. 10.1016/j.febslet.2006.07.016.
Higuchi T, Watanabe Y, Yamamoto M: Protein kinase A regulates sexual development and gluconeogenesis through phosphorylation of the Zn finger transcriptional activator Rst2p in fission yeast. Mol Cell Biol. 2002, 22: 1-11. 10.1128/MCB.22.1.1-11.2002.
Shiozaki K, Russell P: Conjugation, meiosis, and the osmotic stress response are regulated by Spc1 kinase through Atf1 transcription factor in fission yeast. Genes Dev. 1996, 10: 2276-2288.
Wilkinson MG, Samuels M, Takeda T, Toone WM, Shieh JC, Toda T, Millar JB, Jones N: The Atf1 transcription factor is a target for the Sty1 stress-activated MAP kinase pathway in fission yeast. Genes Dev. 1996, 10: 2289-2301.
Mata J, Bahler J: Global roles of Ste11p, cell type, and pheromone in the control of gene expression during early sexual differentiation in fission yeast. Proc Natl Acad Sci USA. 2006, 103: 15517-15522. 10.1073/pnas.0603403103.
Mochizuki N, Yamamoto M: Reduction in the intracellular cAMP level triggers initiation of sexual development in fission yeast. Mol Gen Genet. 1992, 233: 17-24. 10.1007/BF00587556.
Davidson MK, Shandilya HK, Hirota K, Ohta K, Wahls WP: Atf1-Pcr1-M26 complex links stress-activated MAPK and cAMP-dependent protein kinase pathways via chromatin remodeling of cgs2+. J Biol Chem. 2004, 279: 50857-50863. 10.1074/jbc.M409079200.
Yamashita YM, Nakaseko Y, Samejima I, Kumada K, Yamada H, Michaelson D, Yanagida M: 20S cyclosome complex formation and proteolytic activity inhibited by the cAMP/PKA pathway. Nature. 1996, 384: 276-279. 10.1038/384276a0.
Stern B, Nurse P: Cyclin B proteolysis and the cyclin-dependent kinase inhibitor rum1p are required for pheromone-induced G1 arrest in fission yeast. Mol Biol Cell. 1998, 9: 1309-1321.
Costanzo M, Nishikawa JL, Tang X, Millman JS, Schub O, Breitkreuz K, Dewar D, Rupes I, Andrews B, Tyers M: CDK activity antagonizes Whi5, an inhibitor of G1/S transcription in yeast. Cell. 2004, 117: 899-913. 10.1016/j.cell.2004.05.024.
de Bruin RA, McDonald WH, Kalashnikova TI, Yates J, Wittenberg C: Cln3 activates G1-specific transcription via phosphorylation of the SBF bound repressor Whi5. Cell. 2004, 117: 887-898. 10.1016/j.cell.2004.05.025.
Jeong HT, Ozoe F, Tanaka K, Nakagawa T, Matsuda H, Kawamukai M: A novel gene, msa1, inhibits sexual differentiation in Schizosaccharomyces pombe. Genetics. 2004, 167: 77-91. 10.1534/genetics.167.1.77.
Pereira PS, Jones NC: The RGS domain-containing fission yeast protein, Rgs1p, regulates pheromone signalling and is required for mating. Genes Cells. 2001, 6: 789-802. 10.1046/j.1365-2443.2001.00465.x.
Nakashima A, Ueno M, Ushimaru T, Uritani M: Involvement of a CCAAT-binding complex in the expression of a nitrogen-starvation-specific gene, isp6+, in Schizosaccharomyces pombe. Biosci Biotechnol Biochem. 2002, 66: 2224-2227. 10.1271/bbb.66.2224.
Yamano H, Kitamura K, Kominami K, Lehmann A, Katayama S, Hunt T, Toda T: The spike of S phase cyclin Cig2 expression at the G1-S border in fission yeast requires both APC and SCF ubiquitin ligases. Mol Cell. 2000, 6: 1377-1387. 10.1016/S1097-2765(00)00135-0.
Borgne A, Murakami H, Ayte J, Nurse P: The G1/S cyclin Cig2p during meiosis in fission yeast. Mol Biol Cell. 2002, 13: 2080-2090.
Fukushima K, Tanaka Y, Nabeshima K, Yoneki T, Tougan T, Tanaka S, Nojima H: Dmc1 of Schizosaccharomyces pombe plays a role in meiotic recombination. Nucleic Acids Res. 2000, 28: 2709-2716. 10.1093/nar/28.14.2709.
Chikashige Y, Tsutsumi C, Yamane M, Okamasa K, Haraguchi T, Hiraoka Y: Meiotic proteins bqt1 and bqt2 tether telomeres to form the bouquet arrangement of chromosomes. Cell. 2006, 125: 59-69. 10.1016/j.cell.2006.01.048.
Chikashige Y, Ding DQ, Imai Y, Yamamoto M, Haraguchi T, Hiraoka Y: Meiotic nuclear reorganization: switching the position of centromeres and telomeres in the fission yeast Schizosaccharomyces pombe. EMBO J. 1997, 16: 193-202. 10.1093/emboj/16.1.193.
Miki F, Okazaki K, Shimanuki M, Yamamoto A, Hiraoka Y, Niwa O: The 14-kDa Dynein Light Chain-Family Protein Dlc1 Is Required for Regular Oscillatory Nuclear Movement and Efficient Recombination during Meiotic Prophase in Fission Yeast. Mol Biol Cell. 2002, 13: 930-946. 10.1091/mbc.01-11-0543.
Niccoli T, Yamashita A, Nurse P, Yamamoto M: The p150-Glued Ssm4p regulates microtubular dynamics and nuclear movement in fission yeast. J Cell Sci. 2004, 117: 5543-5556. 10.1242/jcs.01475.
Egel R, Willer M, Kjaerulff S, Davey J, Nielsen O: Assessment of pheromone production and response in fission yeast by a halo test of induced sporulation. Yeast. 1994, 10: 1347-1354. 10.1002/yea.320101012.
Lyne R, Burns G, Mata J, Penkett CJ, Rustici G, Chen D, Langford C, Vetrie D, Bahler J: Whole-genome microarrays of fission yeast: characteristics, accuracy, reproducibility, and processing of array data. BMC Genomics. 2003, 4: 27-10.1186/1471-2164-4-27.
Xue Y, Haas SA, Brino L, Gusnanto A, Reimers M, Talibi D, Vingron M, Ekwall K, Wright AP: A DNA microarray for fission yeast: minimal changes in global gene expression after temperature shift. Yeast. 2004, 21: 25-39. 10.1002/yea.1053.
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/projects/geo]
Regulatory Site Analysis Tools package. [http://rsat.ulb.ac.be/rsat/]
Bahler J, Wu JQ, Longtine MS, Shah NG, McKenzie A, Steever AB, Wach A, Philippsen P, Pringle JR: Heterologous modules for efficient and versatile PCR-based gene targeting in Schizosaccharomyces pombe. Yeast. 1998, 14: 943-951. 10.1002/(SICI)1097-0061(199807)14:10<943::AID-YEA292>3.0.CO;2-Y.
Stern B, Nurse P: Fission yeast pheromone blocks S-phase by inhibiting the G1 cyclin B-p34cdc2 kinase. EMBO J. 1997, 16: 534-544. 10.1093/emboj/16.3.534.
Sugiyama A, Tanaka K, Okazaki K, Nojima H, Okayama H: A zinc finger protein controls the onset of premeiotic DNA synthesis of fission yeast in a Mei2-independent cascade. EMBO J. 1994, 13: 1881-1887.
Christensen PU, Davey J, Nielsen O: The Schizosaccharomyces pombe mam1 gene encodes an ABC transporter mediating secretion of M-factor. Mol Gen Genet. 1997, 255: 226-236. 10.1007/s004380050493.
Kitamura K, Shimoda C: The Schizosaccharomyces pombe mam2 gene encodes a putative pheromone receptor which has a significant homology with the Saccharomyces cerevisiae Ste2 protein. EMBO J. 1991, 10: 3743-3751.
Nielsen O, Davey J, Egel R: The ras1 function of Schizosaccharomyces pombe mediates pheromone-induced transcription. EMBO J. 1992, 11: 1391-1395.
Imai Y, Yamamoto M: The fission yeast mating pheromone P-factor: its molecular structure, gene structure, and ability to induce gene expression and G1 arrest in the mating partner. Genes Dev. 1994, 8: 328-338.
Tanaka K, Davey J, Imai Y, Yamamoto M: Schizosaccharomyces pombe map3 + encodes the putative M-factor receptor. Mol Cell Biol. 1993, 13: 80-88.
We thank Karin Holm for expert technical assistance and Richard Egel for comments on the manuscript. This work was supported by the Danish Natural Science Research Council (ON), the Swedish Research Council (AW), the European Commission (QLRT-1999-30174, AW) and the Nordic Academy for Advanced Study (AW and ON). AW is a senior investigator supported by the Swedish Research Council.
Y-XF did all microarray experiments including normalization, data processing, identification of enriched GO categories and contributed to the writing process. SK mapped the 3'end of SPAC3F10.09 and analysed the mating behaviour of ΔSPAC1565.03. CH characterised the ΔSPBC4.01 phenotype. AW did data mining, analysed the relationship between P-factor and M-factor induced genes, quantiated the Ste11 binding sites and contributed to the writing process. ON did all strain constructions, RNA preps and realtime/RT PCR analyses, and also prepared the manuscript draft. All authors have read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Xue-Franzén, Y., Kjærulff, S., Holmberg, C. et al. Genomewide identification of pheromone-targeted transcription in fission yeast. BMC Genomics 7, 303 (2006). https://doi.org/10.1186/1471-2164-7-303
- Fission Yeast
- Meiotic Recombination
- Nitrogen Starvation
- Opposite Mating Type
- Pheromone Response