Modulation of the transcription regulatory program in yeast cells committed to sporulation
© Friedlander et al.; licensee BioMed Central Ltd. 2006
Received: 12 October 2005
Accepted: 9 February 2006
Published: 8 March 2006
Meiosis in budding yeast is coupled to the process of sporulation, where the four haploid nuclei are packaged into a gamete. This differentiation process is characterized by a point of transition, termed commitment, when it becomes independent of the environment. Not much is known about the mechanisms underlying commitment, but it is often assumed that positive feedback loops stabilize the underlying gene-expression cascade.
We describe the gene-expression program of committed cells. Sporulating cells were transferred back to growth medium at different stages of the process, and their transcription response was characterized. Most sporulation-induced genes were immediately downregulated upon transfer, even in committed cells that continued to sporulate. Focusing on the metabolic-related transcription response, we observed that pre-committed cells, as well as mature spores, responded to the transfer to growth medium in essentially the same way that vegetative cells responded to glucose. In contrast, committed cells elicited a dramatically different response.
Our results suggest that cells ensure commitment to sporulation not by stabilizing the process, but by modulating their gene-expression program in an active manner. This unique transcriptional program may optimize sporulation in an environment-specific manner.
Meiosis is a specialized cell division by which haploid gametes are generated from diploid cells. The principal features of meiosis are common to all eukaryotic organisms and include a single round of DNA replication ('premeiotic' replication) followed by two consecutive nuclear divisions, meiosis I and meiosis II. In the first meiotic division homologous chromosomes segregate to opposite poles, whereas in the second division the two sister chromatids separate from each other. Meiosis is characterized by a high frequency of recombination events, occurring during a prolonged prophase that separates DNA replication from the first meiotic division. This genetic exchange between homologous chromosomes ensures that they segregate properly and that the offspring differ genetically from their parents and from each other.
Diploid yeast cells undergoing early stages of meiosis and sporulation may return to the mitotic cell division if provided with nutrients, especially glucose and a nitrogen source. These return-to-growth (RTG) cells complete some meiotic processes, such as high-frequency recombination [12, 13], but switch to the mitotic mode of chromosome segregation and produce diploid cells, following a single division. It is not yet clear how RTG cells resolve the high frequency of double-strand breaks to ensure chromosomal fidelity, and how the cells reinitiate the mitotic cell-cycle program after DNA replication but before chromosomal separation.
RTG may occur only up to a certain stage in the sporulation process. Beyond this stage, cells will continue with sporulation events even if challenged with rich nutritional conditions. This stage of irreversibility was termed "commitment to meiosis" ; it occurs after premeiotic DNA replication and recombination have taken place, but before meiosis I. This commitment is probably essential for ensuring cell viability, as the sporulation process involves drastic changes in morphology, chromosomal state and cell-wall composition, which may be hazardous if abandoned before completion.
The molecular mechanism underlying commitment is not understood. It was proposed that commitment to the meiotic process involves the separation of spindle pole bodies (SPBs) . Indeed, SPB separation was shown to correlate with commitment to mitosis . Conditions that impair late sporulation events may also impact on commitment . In particular, cells mutated in the SPO14 gene, which codes for phospholipase D and is required for late sporulation events, are able to return to growth even after completing meiosis I . The role of phospholipase D-dependent signaling in the progression of the meiotic program is not yet clear.
To better understand the process of commitment to meiosis, we examined the genome-wide transcription response triggered by the transfer of sporulating cells to rich growth medium at different stages of the process. At all stages, the transfer initiated large-scale changes in the gene-expression pattern. The majority of genes that were induced during sporulation were immediately repressed upon the transfer to growth medium, even in committed cells that continued to sporulate. At the same time, committed cells displayed a unique response to nutrients, which was dramatically different from that observed in all other cell types examined (pre-committed cells, spores or vegetative cells exposed to glucose). This unique response consisted of metabolic-related genes, as well as genes involved in competing developmental processes such as pseudo-hyphal growth. Our findings suggest that commitment to meiosis is achieved not by stabilizing the transcription cascade, but rather by an active modulation of the gene-expression program, the detailed nature of which is only starting to be unraveled.
To examine the transcription program associated with commitment to sporulation, we employed the RTG experimental paradigm [12, 18]. The sporulation process was initiated by transferring cells to sporulation medium (SPM). Cells were allowed to progress through the process for varying lengths of times, and were then transferred back to nutrient-containing media (Figure 1c). Cells that were transferred early, before the initiation of the first meiotic division, grew buds and resumed mitotic divisions. In contrast, most cells transferred at later stages continued the sporulation process to produce spores [12, 14, 18] (Figure 1d-e).
To define the time of commitment, we first followed meiotic landmarks events (Figure 1d). Previous studies associated commitment with a time before the completion of meiosis I, which in our conditions occurred at around 5-6 hours in SPM. Indeed, cells that were transferred to yeast extract/peptone/dextrose medium (YPD) after five hours in SPM continued through the second meiotic division also in YPD (Figure 1e). As a more direct assay of commitment, we transferred cells from SPM to a solution containing 4% glucose (see Materials and methods). Glucose is a potent inhibitor of sporulation initiation, and indeed, cells that were transferred early (less than five hours in SPM) abandoned the sporulation program and arrested the cell cycle. In contrast, most cells transferred at a later stage, after being more than five hours in SPM, continued the sporulation process and generated mature spores (Figure 1f). Taken together, we conclude that commitment occurs at around 5-6 hours in SPM. Significantly, mature spores began to appear only at about 8 hours in SPM, reaching the maximum percentage at about 12 hours in SPM (Figure 1d). Thus, the sporulation process continued in rich medium for a considerable period of time before cells became mature spores.
The transfer of sporulating cells to YPD led to a large-scale change in gene expression. Following the transfer, about 1,000 genes were induced by at least twofold. The identity of the induced genes differed between early or late sporulating cells (see below). The number of repressed genes varied somewhat between early and late sporulating cells: around 1,200 genes were repressed over twofold in cells that were transferred early, whereas only around 480 genes were repressed when cells were transferred later in the sporulation process. Evidently, sporulating cells sense and respond to the growth medium at all stages of the sporulation process, even after they have become committed to its completion.
The return to growth (RTG) transcription program
To try and infer the stage through which the cells re-enter the mitotic cell cycle, we examined the response of genes known to be regulated during different phases of the mitotic cell cycle [20, 21]. The G1 cyclin gene CLN3 was immediately induced on the addition of YPD (Figure 3c). Indeed, Cln3 promotes entry to the cell cycle  and serves as a negative regulator of meiotic initiation [23, 24]. Its induction during the return to growth may thus assist in switching from meiosis to mitotic growth. None of the other cyclins or cyclin-dependent kinases was similarly induced. The overall pattern of gene induction, however, was hard to interpret as it did not resemble any of the mitotic cell-cycle phases. For example, although a significant portion of the genes that are upregulated during either G1 or the M phase of the mitotic cycle were induced also during RTG, others, with similar expression patterns during the mitotic cycle, were not (see Additional data file 3). A likely explanation is that individual cells induce different cell-cycle genes depending on which stage of meiosis they are coming from, and this mixture of cell-cycle genes reflects the incomplete synchronization of our culture. In support of that interpretation, genes that are induced during the S phase of the mitotic cell cycle were induced during RTG only when the transfer was done early enough (two hours in SPM) but not later (see Additional data file 3).
In addition to reprogramming their gene expression, RTG cells need to resolve meiosis-specific events and structures. For example, the meiosis-specific alignment of homologous chromosomes during the first meiotic division [25–27] does not exist during mitosis. Indeed, the synaptonemal complex (SC) structures, associated with meiotic pairing, disappear rapidly upon RTG . Moreover, meiosis is characterized by a high level of post-replication double-strand breaks (DSBs), which are required for initiating recombination events. Single-gene studies indicated that in addition to the meiotic repair pathways, additional mitotic DNA repair pathways are used to resolve those breaks during RTG . Not much is known about the processes that are involved in the RTG program, however.
To characterize RTG pathways we analyzed the expression pattern of genes associated with DNA-related processes . Interestingly, most such genes were induced either during sporulation or during the return to the mitotic cell cycle, with only a few induced during both processes (see, for example, Figure 3d). It is likely that genes that are induced during RTG, such as HAM1, RAD55, and MSH1, participate in this process. A comprehensive classification of genes according to their time of induction is provided in Additional data file 3. The classification of homologous genes is also given in Additional data file 3.
Downregulation of sporulation-specific genes in committed cells
Large-scale phenotypic studies have shown that only a portion of the genes that are induced during sporulation are also required for the process [9, 29]. One possibility is that the sporulation-induced genes that were repressed upon the transfer to YPD are not in fact required for sporulation. To test this possibility, we examined in detail the properties of the repressed genes. This set included numerous cell-cycle genes that are required for the two consecutive meiotic divisions [30–32], such as genes coding for most components of the anaphase-promoting complex (APC), its activator Cdc20, the B-type cyclins Clb1,3-6, the polo-like kinase Cdc5, the securin Pds1, the gene CDC15 and the kinesin-like coding gene KAR3 (see Additional data file 3). Several genes involved in spore-wall formation were also identified. Thus, many of the genes that were immediately repressed upon the addition of YPD have a well-established role in sporulation.
Moreover, the major meiotic regulator genes such as IME2 and SMK1, were also repressed upon transfer to rich medium (see Additional data file 2). In particular, NDT80, the master regulator of the mid-sporulation genes, was downregulated as well (Figure 4b). This indicates that the capacity of Ndt80 to autoactivate its gene expression is not sufficient for stabilizing its expression on exposure to rich growth medium. Interestingly, neither glucose alone nor acetate- and nitrogen-containing growth medium (yeast extract/peptone/potassium acetate (YPA)) were sufficient to cause this downregulation, indicating a combinatorial requirement for both factors (Figure 4a).
Downregulation of sporulation-induced genes in committed cells is independent of NDT80 or SUM1
The general response to glucose is dramatically modified in committed cells
Taken together, our results indicated that the transfer of sporulating cells to glucose-containing growth medium altered the expression of most sporulation-specific genes, even in committed cells that proceed to become viable spores in rich medium. We next asked whether the transfer of sporulating cells to YPD altered the expression of additional genes, which are not part of the normal sporulation cascade. Such a response could reflect general effects of the growth medium. Alternatively, it could also indicate specific regulation that is important for ensuring the continuation of meiotic progression in rich medium.
Glucose is also a potent regulator of gene expression in vegetative cells. The transcription response of yeast cells to glucose was characterized in a recent study , showing a dramatic modification of the transcription program. Altered expression was observed for numerous genes involved in metabolic processes, as well as in protein synthesis . We asked whether the addition of glucose-containing growth medium to sporulating cells invokes a similar metabolic response. To eliminate differences resulting from the response of sporulation-specific genes, we focused on 936 genes whose expression was altered by the addition of glucose to vegetative cells  and which showed a significant response also in our experiment.
To examine whether this distinct transcription response to glucose is a property of cells maturing in the sporulation process, we also examined the response of fully developed spores, which have been in sporulation medium for three days. In those cells, the addition of rich medium initiates the process of germination. Interestingly, spores responded to rich medium in essentially the same manner as vegetative or pre-committed cells (Figure 6a). We conclude that the modified transcription response, observed in committed cells, characterizes the sporulation process itself, and not the mature spore state.
Gene expression during sporulation in YPD
Next we asked which genes are induced specifically in committed cells on transfer to YPD. Such genes could be classified into two groups. First, a 'process-linked' group was defined. Genes in this group were induced during late meiosis in SPM (at around 8-9 hours), and were also induced in YPD, at variable times that appeared to be linked to the progression of sporulation (Figure 6e). Second, a 'commitment-specific' group was defined, which included genes that were induced by YPD specifically in committed cells (Figure 6f). Genes in this group were not induced during normal meiosis in SPM, and were also insensitive to YPD during early sporulation.
About 60 genes display a 'process-linked' expression pattern (see Figure 6e and Additional data file 3). Among these we identified two sporulation-specific genes, which are involved in spore-wall formation (DIT1 and SPS100). Surprisingly, however, most genes in this group are not in fact classified as sporulation genes, and do not have a recognized role in the process. Rather, most genes in this group are associated with the general environmental stress response (ESR), which is invoked in response to a variety of different environmental stresses [36, 37]. Examples include genes coding for heat-shock proteins (such as HSP12, GRE1, and SIP18), for the stress-induced transcription factor Xbp1, and for the Tor1 kinase. Accordingly, the promoter regions of most of those genes contained the stress-response element CCCCT. The induction of this gene group may reflect some stress signal associated with the progression of sporulation.
The 'commitment-specific' gene group included around 65 genes that were induced by YPD specifically after commitment. In parallel, about 50 genes were repressed by YPD specifically after commitment (see Figure 6f-g and Additional data file 3). Those genes may be part of a modified transcription program that assists the continuation of the sporulation process in rich medium. Interestingly, the repressed group includes two genes that function in the pseudohyphal growth pathway (MUC1 and MSB2). It is likely that repression of these genes assists in the inhibition of this alternative developmental route.
Differentiation processes proceed through a sequence of events involving coordinate modulations of morphology and gene expression. In many cases, once the process has passed a certain stage it becomes determined and loses its dependence on environmental signals. This transition is termed commitment. Mechanisms for achieving commitment can be classified into two broad classes. First, a passive mechanism might be imagined, whereby the process becomes effectively insulated from the external signals. Alternatively, commitment might require an active modulation of the signaling apparatus, which senses the external signal, but interprets it in a specific manner that optimizes the continuation of the process in changing conditions.
Differentiating cells could be insulated from the environment through the inhibition of essential sensory receptors, such as the receptors for glucose or nitrogen in the case of yeast sporulation. Effective insulation of the differentiation process could also be achieved through positive feedback loops that would render the process self-sustaining. Such positive feedbacks are ubiquitous in differentiation cascades. In the case of the Xenopus oocyte at least this positive feedback loop has been shown to be important for commitment . Also in the case of yeast sporulation, the capacity of Ndt80 to autoactivate its own gene expression could in principle implement such a feedback.
Our gene-expression analysis ruled out the possibility that commitment to sporulation is achieved through a passive mechanism. Cells responded to nutrients at all stages of the sporulation process (see Additional data file 5 for a matlab program that enables the reader to view our data). In fact, within 5 minutes of the addition of rich medium, we observed a change in the expression pattern of over 1,000 genes (around 15% of the yeast genome). Importantly, this large-scale transcription response was dramatically different in committed versus non-committed cells, and was highly specific to the type of medium added. In particular, the response seen on addition of YPD (containing both glucose and nitrogen) was distinct from that observed on the addition of YPA (containing nitrogen and acetate, but not glucose) or on the addition of glucose alone. This indicates that the reduction in mRNA is a result of combinatorial regulation. The machinery of committed cells responds to the combination of glucose and nitrogen signals. This behavior suggests that active modulation of the response to external signals takes place. This modulation may play a central role in enabling the continuation of the sporulation process after the addition of nutrients.
How do the cells achieve commitment to meiosis? The continuation of sporulation in rich media poses two complementary challenges. First, cells need to overcome mitogenic signals, which in vegetative cells facilitate the mitotic program. Second, cells need to ensure the continuation of the sporulation process itself. Our data suggest that these two aspects are goverened by complementary molecular mechanisms.
The mitogenic effects of glucose are mediated to a large extent by the PKA pathway [42, 43]. This pathway represents a central junction in the choice between the meiotic and the mitotic programs: Entry into the mitotic cell cycle requires high activity of this pathway, whereas entry into meiosis requires low PKA activity [40, 44–48]. Interestingly, our data suggest that committed cells respond to glucose not by activating the PKA pathway, but rather by inhibiting it through the induction of BCY1, the gene coding for the negative regulatory unit of the PKAs. This inhibition is probably required in order to overcome mitogenic signals, to ensure the continuation of the meiotic program upon the addition of nutrients.
The addition of either YPD or glucose to committed cells failed to induce the typical set of glucose-responsive genes. In contrast, the addition of YPD (but not of glucose alone) had a dramatic effect on the sporulation cascade itself. In fact, the vast majority of genes that are induced during sporulation in SPM were immediately repressed. This repression did not reflect the lack of functional requirement, as many of the repressed genes had a well defined role in the two meiotic divisions or in the formation of the spore wall. It may be that the proteins encoded by these genes are already synthesized at the time of commitment in sufficient amounts to complete the meiotic division. Alternatively, it is possible that the downregulation of at least some of the mid-sporulation genes is in fact required to prevent a shift back to the mitotic cell cycle upon the addition of nutrients, thus augmenting the cells' commitment. Further work is needed in order to better understand the nature of this downregulation and its potential contribution to the commitment process.
Whereas the transcription program characterizing sporulation in sporulation medium (SPM) was practically eliminated in YPD, a group of around 60 late-sporulation genes was also induced in YPD, in a temporal manner that appears to be linked to the progression of the sporulation process. Surprisingly, the vast majority of these genes do not have a known role in sporulation. Rather, this group included a majority of stress-related genes, which are also induced in response to a variety of environmental stresses . The timely induction of these genes may indicate the generation of some stress-related signal. Such stress could be initiated, for example, by the onset of spore-wall deposition, similarly to the stress signal that ensues after the formation of a mating projection (a shmoo) during the mating process . In the case of mating, this stress signal is propagated through the protein kinase C (PKC) pathway to induce the second wave of gene expression. It is tempting to propose that here also, an analogous stress signal may be associated with the commitment, by marking the initiation of a particular developmental stage (for example, spore-wall formation) and triggering subsequent processes required for the completion of sporulation. It is of interest that the formation of the pro-spore membrane is initiated on the SPBs that were implicated in commitment . Such a proposal may also be consistent with a potential role for SPO14 in both the formation of a spore membrane and the commitment to sporulation . This proposition awaits further experimental validation.
In conclusion, previous theoretical and experimental work has shown that bistability, generated through positive feedback loops, can render a process irreversible by effectively isolating it from external signals [41, 51–53]. In contrast, sporulating cells modulate their transcription program in response to changing conditions, proceeding through different alternative paths. This alternative strategy is likely to be beneficial when a differentiating cell is challenged not only with the removal of the initiating signal but also with competing signals that could potentially interfere with structural or morphological processes, or could signal alternative fates. The experimental design presented here could be applied to distinguish between these alternatives. Further studies are required to examine which of the two strategies is more prominent during cellular differentiation.
Materials and methods
Yeast strains and plasmids
All strains used for this study are of the SK1 genetic background and are listed in Additional data file 3. The pSPS4-3HA-NDT80 strain (GF18) was constructed by a one-step PCR-based replacement method  using the plasmid pFA6a-pSPS4-3HA-KanMX6. The latter was constructed by replacing the GAL1 promoter fragment in pFA6a-pGAL-3HA-KanMX6 with a 1 kb PCR fragment of the SPS4 promoter. Transformed haploids NKY1712 were mated with strain NE30. Integration to the correct site was verified by PCR. Diploids were selected on minimally supplemented yeast synthetic dextrose (SD) plates and were checked for sporulation efficiency on sporulation plates.
Sporulation, return to growth, and germination experiments
The compositions of all media are given in Additional data file 3. Cells (NKY1551, JPY214, or GF18) were grown to saturation in YPD for 24 hours, diluted into YPA and grown overnight. The cells were then washed twice in sterile water and resuspended in SPM to initiate the sporulation process. At different stages of the sporulation process, cells were centrifuged and resuspended in rich medium. At each time point, 5 × 108 cells were frozen in liquid nitrogen. To monitor the meiotic divisions, 2 ml of cells were kept for staining with 4,6-diamidino-2-phenylindole (DAPI). To monitor recombination frequencies at the HIS4 locus, cells were plated on both YPD and -His plates, and colonies were counted after three days incubation at 30°C. All liquid cultures were grown at 30°C with constant shaking. Sporulation efficiency was determined by counting asci in the overnight SPM culture. For germination experiments, cells (DS1) were grown to saturation in YPD for 24 hours and plated on sporulation plates. Three-day-old asci were harvested and suspended in YPD at 30°C to initiate spore germination.
Cells (2 ml) were centrifuged for 30 seconds, fixed in 0.25 mM Tris, 70% ethanol and kept at 4°C. Before staining, cells were washed in 1× PBS solution (137 mM NaCl, 2.7 mM KCl, 8 mM Na2HPO4, 2 mM KH2PO4), sonicated for 5 seconds and centrifuged. One microliter DAPI (at a concentration of 25 μg/ml) was added to the pellet and incubated in the dark for 10 minutes at room temperature. The number of nuclei per cell was counted using a fluorescence microscope. Around 300 cells were analyzed for each sample.
Defining the timing of commitment
To define the time of commitment, we transferred sporulating cells from SPM to a solution containing 4% glucose at various times (Figure 1f). Glucose is a potent inhibitor of sporulation initiation, but glucose alone does not support growth or germination. For each time point, we calculated the fraction of cells that became spores in over night glucose cultures. We used glucose rather then YPD in this experiment because glucose does not support growth of vegetative cells or germination of mature spores. This allows for a precise quantification of the fraction of spores after overnight incubation. In contrast, on transfer to YPD, the uncommitted cells initiate growth and cell division, and quickly take over the population. Moreover, as the culture is not fully synchronized, some spores may start to germinate early, while others have not yet become spores. We normalized for this effect when looking at the completion of meiotic division II among DAPI-stained cells (Figure 1e), but such normalization is difficult for the longer times corresponding to spore maturation. Thus, although we clearly saw the formation of mature spores in cells transferred to YPD after five hours in SPM, it was difficult to use this fraction as a quantitative measure for commitment.
RNA extraction and labeling
Total RNA was extracted using the RNeasy Midi Kit (Qiagen, Valencia, USA) and reverse transcribed using M-MLV reverse transcriptase RNase H Minus (Promega, Madison, USA). cDNA products were labeled with Cy3 and Cy5 by the indirect amino-allyl method  with minor modifications. Total RNA (20 μg) was combined with 8 μg of oligo(dT) 12- to 18-mer (Amersham Pharmacia Biotech, Little Chalfont, UK) in a final volume of 26 μl and heated to 70°C for 10 minutes to denature RNA tertiary structures. The mixture was then transferred to 45°C to anneal the primers. A preheated (50°C) 12 μl mixture containing reaction buffer (Promega), 2.5 mM MgCl2, 0.5 mM dATP, dCTP and dGTP, 0.1 mM dTTP, 0.4 mM amino-allyl dUTP (Ambion, Austin, USA), RNase H Minus and 40 Units RNasin (Promega) was added to the RNA, and the reaction mixture was incubated for 2 hours at 45°C with 400 units of Superscript II (Gibco, Carlsbad, USA). An additional 400 units of the enzyme was applied to the reaction mixture, which was incubated for another 2 hours. After incubation, RNA was hydrolyzed by adding 10 μl of 1 M NaOH and incubating at 65°C for 30 minutes. The reaction mixture was neutralized by titration with 1 M Tris-HCl pH 7.5, and DNA was precipitated with ethanol, resuspended in 10 μl 0.1 M sodium carbonate buffer pH 9.3. Cy3 and Cy5 (Amersham) were dissolved in 10 μl DMSO. Cy5 or Cy3 (1.25 μl) was added to the sample and incubated at 25°C for 1 hour. Labeled cDNA was purified by Strataprep PCR purification kit (Stratagene, La Jolla, USA) according to the manufacturer's instructions. Dye incorporation was measured using a spectrophotometer. Reference RNA for all microarrays in this work was from YPA-grown cells, just before the transfer to SPM.
For each hybridization, cDNA samples were labeled with Cy3 and Cy5 and combined with blockers: 5 μg herring sperm (Promega), 5 μg tRNA (Gibco) and 17.5 μg poly(A) (poly(A) oligonucleotides were synthesized at mixed lengths of 40, 50, and 60 adenine residues). The labeled cDNAs were concentrated to 40 μl using Microcon (Millipore, Bedford, USA) and 40 μl of 2× hybridization solution (10× SSC, 50% formamide, 0.2% SDS) was added. Microarrays containing all yeast open reading frames (ORFs) were prehybridized by incubating at 42°C for 45 minutes in a solution containing 1% BSA, 25% formamide, 5× SSC and 0.1% SDS. The slides were washed in sterile water and dried by centrifugation (3 minutes, 2,000 rpm). The labeled samples were boiled for 5 minutes, centrifuged for 1 minute, hybridized on the slide and placed in a hybridization chamber (Corning, Corning, USA) for overnight incubation at 42°C. The slides were then washed for 5 minutes at 42°C with a solution containing 2× SSC and 0.1% SDS. An additional wash was performed at room temperature with a solution containing 0.1× SSC and 0.1% SDS, followed by three additional washes at room temperature in 0.1× SSC solution. For most experiments, we used arrays purchased from the University Health Network Microarray Center, Toronto, Canada, where PCR products are printed. For the YPA experiment (35 arrays) we used slides obtained from F. Holstege (UMC Utrecht, the Netherlands), where genes are represented by 70-base oligonucleotides . Control experiments have been performed to check that differences in the data were not caused by differences in the arrays (these data have been deposited in the Gene Expression Omnibus (GEO) database run by the National Center for Biotechnology Information (NCBI) [57–59] and are accessible through GEO Series accession number GSE3820).
Scanning and quantification
Images were obtained using ScanArray 4000 (Packard BioScience, Meriden, USA). Image analysis was performed using QuantArray version 3 software (PerkinElmer Life Sciences, Boston, USA), with spots quantified using the adaptive method. Low-quality spots were discarded following detailed visual inspection. The data were then transformed into log2 ratios, and normalized by subtracting the median. Values of replicate spots on the slides were averaged (see Additional data file 1 for MIAME data). The data discussed in this publication have been deposited in GEO  and are accessible through GEO Series accession numbers: GSE3814, GSE3815, GSE3816, GSE3817, GSE3818, GSE3819 and GSE3820 [57, 58] (see also Additional data file 4).
Calculating the correlation coefficients with results of previous experiments
To compare our sporulation time course with previous experiments, the different time courses were first aligned such that similar time points corresponded to each other (progression into sporulation) [60, 61]. In each pairwise comparison, one experiment (A) was held constant and the other (B) was transformed to best match the time points of A. The transformation was done in two steps: first, smoothing splines were used to add intermediate time points to B (by creating a continuous representation of B, and re-sampling it at a higher frequency) ; second, we searched for the time points of B that maximizes the correlations between each time point in A and the corresponding time point in B over all genes. We only considered the time points of A, or transformations of the kind T' = T*(1 + k1) + k2, where T are the time points of A, k1 is a 'stretching' parameter, and k2 is a 'shifting' parameter. After transforming B, we calculated the average correlation between corresponding time points in A and B, over all genes (the same measure that was used to determine the optimal transformation). For each pairwise comparison of A and B this procedure was done twice, once when starting with A and once when starting with B. The average correlation is given.
Comparison of the responses of sporulating cells, vegetative cells and spores to glucose
We included in this analysis genes whose expression was altered by the addition of glucose to vegetative cells. To reduce noise, we filtered out genes whose expression was not altered by YPD in our experiment (transfer of sporulating cells to YPD). We ended with 936 genes (see Additional data file 3). For the vegetative cell data, we calculated the average expression profile over the three given time points (20, 40, and 60 minutes), and considered a given gene to be down- or upregulated if its average expression change exceeded twofold. For the sporulation data, we considered the YPD expression relative to that observed during sporulation, and calculated the average expression change observed at 20 and 40 minutes after the transfer. For spores that completed the sporulation process, we calculated the average expression response observed at 15 and 30 minutes following the addition of YPD. Correlations were calculated between these average profiles. The correlation matrix of the sporulation data was clustered  (separately for glucose-induced and glucose-repressed genes). Only genes belonging to homogeneous clusters were kept.
Identification of gene groups
A few genes with a common expression pattern were chosen as template genes. To find more genes with the same expression pattern we calculated the correlation of each gene in the genome with each of the template genes. The genes were ordered according to their correlation with the template gene. Genes with the highest correlation were defined as genes in the group. For each template gene the threshold was determined by inspection of the expression pattern of the genes obtained. The top genes with the same expression pattern were chosen. The qualitative results of our study do not depend on these thresholds.
Utilization of Gene Ontology annotations for identification of polarized gene expression during sporulation and RTG
For each of the chosen Gene Ontology (GO)  groups we filtered out genes in the group which their expression did not change (we kept only genes whose expression changed more than twofold in one array at least). We then obtained a matrix of the expression of the selected genes in the following conditions: the sporulation time points and transfers to YPD at early times (2, 3, and 4 hours). The correlation between any two genes in this matrix was computed to give the gene-gene correlation matrix. We then clustered  this gene-gene correlation matrix and chose homogeneous clusters.
Additional data files
The following additional data files are included with this article: Additional data file 1 includes minimum information about a microarray experiment (MIAME). Additional data file 2 includes supplementary figures S1 to S3. Figure S1 shows polarized expression during sporulation and RTG. Figure S2 shows the expression pattern of IME2 and SMK1. Figure S3 shows a matrix of pairwise correlations that exemplifies the general metabolic response to glucose. Additional data file 3 includes supplemental tables S1 to S17. Table S1 includes a list of early sporulation genes. Table S2 includes a complete list of genes used in Figure 3d of the paper. Table S3 lists examples of genes that according to their expression during RTG are regulated by the mitotic cell cycle. Table S4 lists genes that are induced during RTG (20 minutes after the transfer). Table S5 lists genes that are induced during RTG (immediately after the transfer). Table S6 lists homologous genes which are induced during sporulation or RTG. Table S7 lists middle sporulation genes that are repressed upon transfer to YPD. Table S8 lists insulated middle sporulation genes. Table S9 lists the 936 genes used for Figure 6a. Table S10 lists rRNA-processing genes. Table S11 lists gluconeogenesis genes. Table S12 lists genes that encode ribosomal proteins. Table S13 lists genes induced in a time-dependent manner. Table S14 lists genes that are induced in response to YPD in committed cells. Table S15 lists genes that are repressed in response to YPD in committed cells. Table S16 includes a list of the yeast strains used in the present study. Table S17 includes the composition of the media used in the present study. Additional data file 4 includes the normalized data of the present study (in log2 ratios). Addtional file 5 is a zip file containing a matlab program that enables to view the expression data discussed in this article. It contains also a help file.
We thank members of our groups for discussions. We are grateful to Itay Tirosh for help in comparing experiments. This work was supported by the NIH Grant No. A150562, and the Israel Ministry of Science and Technology.
- Kassir Y, Granot D, Simchen G: IME1, a positive regulator gene of meiosis in S. cerevisiae. Cell. 1988, 52: 853-862. 10.1016/0092-8674(88)90427-8.PubMedView ArticleGoogle Scholar
- Shefer-Vaida M, Sherman A, Ashkenazi T, Robzyk K, Kassir Y: Positive and negative feedback loops affect the transcription of IME1, a positive regulator of meiosis in Saccharomyces cerevisiae. Dev Genet. 1995, 16: 219-228. 10.1002/dvg.1020160302.PubMedView ArticleGoogle Scholar
- Rubin-Bejerano I, Mandel S, Robzyk K, Kassir Y: Induction of meiosis in Saccharomyces cerevisiae depends on conversion of the transcriptional represssor Ume6 to a positive regulator by its regulated association with the transcriptional activator Ime1. Mol Cell Biol. 1996, 16: 2518-2526.PubMedPubMed CentralView ArticleGoogle Scholar
- Covitz PA, Mitchell AP: Repression by the yeast meiotic inhibitor RME1. Genes Dev. 1993, 7: 1598-1608.PubMedView ArticleGoogle Scholar
- Sagee S, Sherman A, Shenhar G, Robzyk K, Ben-Doy N, Simchen G, Kassir Y: Multiple and distinct activation and repression sequences mediate the regulated transcription of IME1, a transcriptional activator of meiosis-specific genes in Saccharomyces cerevisiae. Mol Cell Biol. 1998, 18: 1985-1995.PubMedPubMed CentralView ArticleGoogle Scholar
- Chu S, Herskowitz I: Gametogenesis in yeast is regulated by a transcriptional cascade dependent on Ndt80. Mol Cell. 1998, 1: 685-696. 10.1016/S1097-2765(00)80068-4.PubMedView ArticleGoogle Scholar
- Chu S, DeRisi J, Eisen M, Mulholland J, Botstein D, Brown PO, Herskowitz I: The transcriptional program of sporulation in budding yeast. Science. 1998, 282: 699-705. 10.1126/science.282.5389.699.PubMedView ArticleGoogle Scholar
- Primig M, Williams RM, Winzeler EA, Tevzadze GG, Conway AR, Hwang SY, Davis RW, Esposito RE: The core meiotic transcriptome in budding yeasts. Nat Genet. 2000, 26: 415-423. 10.1038/82539.PubMedView ArticleGoogle Scholar
- Enyenihi AH, Saunders WS: Large-scale functional genomic analysis of sporulation and meiosis in Saccharomyces cerevisiae. Genetics. 2003, 163: 47-54.PubMedPubMed CentralGoogle Scholar
- Deutschbauer AM, Williams RM, Chu AM, Davis RW: Parallel phenotypic analysis of sporulation and postgermination growth in Saccharomyces cerevisiae. Proc Natl Acad Sci USA. 2002, 99: 15530-15535. 10.1073/pnas.202604399.PubMedPubMed CentralView ArticleGoogle Scholar
- Briza P, Bogengruber E, Thur A, Rutzler M, Munsterkotter M, Dawes IW, Breitenbach M: Systematic analysis of sporulation phenotypes in 624 non-lethal homozygous deletion strains of Saccharomyces cerevisiae. Yeast. 2002, 19: 403-422. 10.1002/yea.843.PubMedView ArticleGoogle Scholar
- Sherman F, Roman H: Evidence for two types of allelic recombination in yeast. Genetics. 1963, 48: 255-261.PubMedPubMed CentralGoogle Scholar
- Zenvirth D, Loidl J, Klein S, Arbel A, Shemesh R, Simchen G: Switching yeast from meiosis to mitosis: double-strand break repair, recombination and synaptonemal complex. Genes Cells. 1997, 2: 487-498. 10.1046/j.1365-2443.1997.1370335.x.PubMedView ArticleGoogle Scholar
- Simchen G, Pinon R, Salts Y: Sporulation in Saccharomyces cerevisiae: premeiotic DNA synthesis, readiness and commitment. Exp Cell Res. 1972, 75: 207-218. 10.1016/0014-4827(72)90538-1.PubMedView ArticleGoogle Scholar
- Horesh O, Simchen G, Friedmann A: Morphogenesis of the synapton during yeast meiosis. Chromosoma. 1979, 75: 101-115. 10.1007/BF00330628.PubMedView ArticleGoogle Scholar
- Hirschberg J, Simchen G: Commitment to the mitotic cell cycle in yeast in relation to meiosis. Exp Cell Res. 1977, 105: 245-252. 10.1016/0014-4827(77)90122-7.PubMedView ArticleGoogle Scholar
- Honigberg SM, Esposito RE: Reversal of cell determination in yeast meiosis: postcommitment arrest allows return to mitotic growth. Proc Natl Acad Sci USA. 1994, 91: 6559-6563.PubMedPubMed CentralView ArticleGoogle Scholar
- Ganesan AT, Holter H, Roberts C: Some observations on sporulation in saccharomyces. CR Trav Lab Carlsberg [Chim]. 1958, 13: 1-6.Google Scholar
- Smith HE, Mitchell AP: A transcriptional cascade governs entry into meiosis in Saccharomyces cerevisiae. Mol Cell Biol. 1989, 9: 2142-2152.PubMedPubMed CentralView ArticleGoogle Scholar
- Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell. 1998, 9: 3273-3297.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho RJ, Campbell MJ, Winzeler EA, Steinmetz L, Conway A, Wodicka L, Wolfsberg TG, Gabrielian AE, Landsman D, Lockhart DJ, et al: A genome-wide transcriptional analysis of the mitotic cell cycle. Mol Cell. 1998, 2: 65-73. 10.1016/S1097-2765(00)80114-8.PubMedView ArticleGoogle Scholar
- Wijnen H, Landman A, Futcher B: The G(1) cyclin Cln3 promotes cell cycle entry via the transcription factor Swi6. Mol Cell Biol. 2002, 22: 4402-4418. 10.1128/MCB.22.12.4402-4418.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Purnapatre K, Piccirillo S, Schneider BL, Honigberg SM: The CLN3/SWI6/CLN2 pathway and SNF1 act sequentially to regulate meiotic initiation in Saccharomyces cerevisiae. Genes Cells. 2002, 7: 675-691. 10.1046/j.1365-2443.2002.00551.x.PubMedView ArticleGoogle Scholar
- Colomina N, Gari E, Gallego C, Herrero E, Aldea M: G1 cyclins block the Ime1 pathway to make mitosis and meiosis incompatible in budding yeast. EMBO J. 1999, 18: 320-329. 10.1093/emboj/18.2.320.PubMedPubMed CentralView ArticleGoogle Scholar
- Martini E, Keeney S: Sex and the single (double-strand) break. Mol Cell. 2002, 9: 700-702. 10.1016/S1097-2765(02)00512-9.PubMedView ArticleGoogle Scholar
- Kleckner N: Meiosis: how could it work?. Proc Natl Acad Sci USA. 1996, 93: 8167-8174. 10.1073/pnas.93.16.8167.PubMedPubMed CentralView ArticleGoogle Scholar
- Roeder GS: Meiotic chromosomes: it takes two to tango. Genes Dev. 1997, 11: 2600-2621.PubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
- Rabitsch KP, Toth A, Galova M, Schleiffer A, Schaffner G, Aigner E, Rupp C, Penkner AM, Moreno-Borchart AC, Primig M, et al: A screen for genes required for meiosis and spore formation based on whole-genome expression. Curr Biol. 2001, 11: 1001-1009. 10.1016/S0960-9822(01)00274-3.PubMedView ArticleGoogle Scholar
- Shanks RM, Bascom-Slack C, Dawson DS: Analysis of the kar3 meiotic arrest in Saccharomyces cerevisiae. Cell Cycle. 2004, 3: 363-371.PubMedView ArticleGoogle Scholar
- Cohen-Fix O: Meiosis: polo, FEAR and the art of dividing reductionally. Curr Biol. 2003, 13: R603-R605. 10.1016/S0960-9822(03)00526-8.PubMedView ArticleGoogle Scholar
- Lee BH, Amon A: Role of Polo-like kinase CDC5 in programming meiosis I chromosome segregation. Science. 2003, 300: 482-486. 10.1126/science.1081846.PubMedView ArticleGoogle Scholar
- Xie J, Pierce M, Gailus-Durner V, Wagner M, Winter E, Vershon AK: Sum1 and Hst1 repress middle sporulation-specific gene expression during mitosis in Saccharomyces cerevisiae. EMBO J. 1999, 18: 6448-6454. 10.1093/emboj/18.22.6448.PubMedPubMed CentralView ArticleGoogle Scholar
- Pak J, Segall J: Regulation of the premiddle and middle phases of expression of the NDT80 gene during sporulation of Saccharomyces cerevisiae. Mol Cell Biol. 2002, 22: 6417-6429. 10.1128/MCB.22.18.6417-6429.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Y, Pierce M, Schneper L, Guldal CG, Zhang X, Tavazoie S, Broach JR: Ras and gpa2 mediate one branch of a redundant glucose signaling pathway in yeast. PLoS Biol. 2004, 2: E128-10.1371/journal.pbio.0020128.PubMedPubMed CentralView ArticleGoogle Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.PubMedPubMed CentralView ArticleGoogle Scholar
- Causton HC, Ren B, Koh SS, Harbison CT, Kanin E, Jennings EG, Lee TI, True HL, Lander ES, Young RA: Remodeling of yeast genome expression in response to environmental changes. Mol Biol Cell. 2001, 12: 323-337.PubMedPubMed CentralView ArticleGoogle Scholar
- Thevelein JM: Signal transduction in yeast. Yeast. 1994, 10: 1753-1790. 10.1002/yea.320101308.PubMedView ArticleGoogle Scholar
- Broach JR, Deschenes RJ: The function of ras genes in Saccharomyces cerevisiae. Adv Cancer Res. 1990, 54: 79-139.PubMedView ArticleGoogle Scholar
- Matsumoto K, Uno I, Ishikawa T: Initiation of meiosis in yeast mutants defective in adenylate cyclase and cyclic AMP-dependent protein kinase. Cell. 1983, 32: 417-423. 10.1016/0092-8674(83)90461-0.PubMedView ArticleGoogle Scholar
- Xiong W, Ferrell JE Jr: A positive-feedback-based bistable 'memory module' that governs a cell fate decision. Nature. 2003, 426: 460-465. 10.1038/nature02089.PubMedView ArticleGoogle Scholar
- Thevelein JM, de Winde JH: Novel sensing mechanisms and targets for the cAMP-protein kinase A pathway in the yeast Saccharomyces cerevisiae. Mol Microbiol. 1999, 33: 904-918. 10.1046/j.1365-2958.1999.01538.x.PubMedView ArticleGoogle Scholar
- Broach JR: RAS genes in Saccharomyces cerevisiae: signal transduction in search of a pathway. Trends Genet. 1991, 7: 28-33. 10.1016/0168-9525(91)90018-L.PubMedView ArticleGoogle Scholar
- Gancedo JM: Control of pseudohyphae formation in Saccharomyces cerevisiae. FEMS Microbiol Rev. 2001, 25: 107-123.PubMedView ArticleGoogle Scholar
- Shenhar G, Kassir Y: A positive regulator of mitosis, Sok2, functions as a negative regulator of meiosis in Saccharomyces cerevisiae. Mol Cell Biol. 2001, 21: 1603-1612. 10.1128/MCB.21.5.1603-1612.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Matsuura A, Treinin M, Mitsuzawa H, Kassir Y, Uno I, Simchen G: The adenylate cyclase/protein kinase cascade regulates entry into meiosis in Saccharomyces cerevisiae through the gene IME1. EMBO J. 1990, 9: 3225-3232.PubMedPubMed CentralGoogle Scholar
- Shilo V, Simchen G, Shilo B: Initiation of meiosis in cell cycle initiation mutants of Saccharomyces cerevisiae. Exp Cell Res. 1978, 112: 241-248. 10.1016/0014-4827(78)90206-9.PubMedView ArticleGoogle Scholar
- Rubin-Bejerano I, Sagee S, Friedman O, Pnueli L, Kassir Y: The in vivo activity of Ime1, the key transcriptional activator of meiosis-specific genes in Saccharomyces cerevisiae, is inhibited by the cyclic AMP/protein kinase A signal pathway through the glycogen synthase kinase 3-beta homolog Rim11. Mol Cell Biol. 2004, 24: 6967-6979. 10.1128/MCB.24.16.6967-6979.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Roberts CJ, Nelson B, Marton MJ, Stoughton R, Meyer MR, Bennett HA, He YD, Dai H, Walker WL, Hughes TR, et al: Signaling and circuitry of multiple MAPK pathways revealed by a matrix of global gene expression profiles. Science. 2000, 287: 873-880. 10.1126/science.287.5454.873.PubMedView ArticleGoogle Scholar
- Honigberg SM, Conicella C, Espositio RE: Commitment to meiosis in Saccharomyces cerevisiae: involvement of the SPO14 gene. Genetics. 1992, 130: 703-716.PubMedPubMed CentralGoogle Scholar
- Jacob F, Monod J: Genetic regulatory mechanisms in the synthesis of proteins. J Mol Biol. 1961, 3: 318-356.PubMedView ArticleGoogle Scholar
- Lisman JE: A mechanism for memory storage insensitive to molecular turnover: a bistable autophosphorylating kinase. Proc Natl Acad Sci USA. 1985, 82: 3055-3057.PubMedPubMed CentralView ArticleGoogle Scholar
- Hofer T, Nathansen H, Lohning M, Radbruch A, Heinrich R: GATA-3 transcriptional imprinting in Th2 lymphocytes: a mathematical model. Proc Natl Acad Sci USA. 2002, 99: 9364-9368. 10.1073/pnas.142284699.PubMedPubMed CentralView ArticleGoogle Scholar
- Longtine MS, McKenzie A, Demarini DJ, Shah NG, Wach A, Brachat A, Philippsen P, Pringle JR: Additional modules for versatile and economical PCR-based gene deletion and modification in Saccharomyces cerevisiae. Yeast. 1998, 14: 953-961. 10.1002/(SICI)1097-0061(199807)14:10<953::AID-YEA293>3.0.CO;2-U.PubMedView ArticleGoogle Scholar
- Guterman I, Shalit M, Menda N, Piestun D, Dafny-Yelin M, Shalev G, Bar E, Davydov O, Ovadis M, Emanuel M, et al: Rose scent: genomics approach to discovering novel floral fragrance-related genes. Plant Cell. 2002, 14: 2325-2338. 10.1105/tpc.005207.PubMedPubMed CentralView ArticleGoogle Scholar
- van de Peppel J, Kemmeren P, van Bakel H, Radonjic M, van Leenen D, Holstege FC: Monitoring global messenger RNA changes in externally controlled microarray experiments. EMBO Rep. 2003, 4: 387-393. 10.1038/sj.embor.embor798.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, Rudnev D, Lash AE, Fujibuchi W, Edgar R: NCBI GEO: mining millions of expression profiles - database and tools. Nucleic Acids Res. 2005, D562-D566.Google Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMedPubMed CentralView ArticleGoogle Scholar
- Gene Expression Omnibus (GEO). [http://www.ncbi.nlm.nih.gov/geo/]
- Bar-Joseph Z, Gerber GK, Gifford DK, Jaakkola TS, Simon I: Continuous representations of time-series gene expression data. J Comput Biol. 2003, 10: 341-356. 10.1089/10665270360688057.PubMedView ArticleGoogle Scholar
- Aach J, Church GM: Aligning gene expression time series with time warping algorithms. Bioinformatics. 2001, 17: 495-508. 10.1093/bioinformatics/17.6.495.PubMedView ArticleGoogle Scholar
- Alon U, Barkai N, Notterman DA, Gish K, Ybarra S, Mack D, Levine AJ: Broad patterns of gene expression revealed by clustering analysis of tumor and normal colon tissues probed by oligonucleotide arrays. Proc Natl Acad Sci USA. 1999, 96: 6745-6750. 10.1073/pnas.96.12.6745.PubMedPubMed CentralView ArticleGoogle Scholar
- Matlab Central, Proportional Venn Diagrams. [http://www.mathworks.com/atlabcentrafileexchangloadFile.do?objectId=6116&objectType=file]
- Ihmels J, Bergmann S, Barkai N: Defining transcription modules using large-scale gene expression data. Bioinformatics. 2004, 20: 1993-2003. 10.1093/bioinformatics/bth166.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.