|
|
|
|
Vol. 17, Issue 11, 4736-4747, November 2006
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

,
,
,

*Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada M5S 1A8;
Banting and Best Department of Medical Research, University of Toronto, Toronto, Ontario, Canada M5G 1L6; and
Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada M5S 3E1
Submitted May 1, 2006;
Revised August 4, 2006;
Accepted August 16, 2006
Monitoring Editor: Trisha Davis
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The morphological landmarks of the cell cycle stages in budding yeast, most notably the size of the bud relative to the size of the mother cell, allows the identification of mutants blocked at specific stages of the cell cycle and thereby forms the basis for the classic cell cycle screens (Hartwell et al., 1970
; Culotti and Hartwell, 1971
; Hartwell, 1971a
, 1971b
, 1973
; Moir et al., 1982
). These genetic screens, using conditional temperature-sensitive mutants, identified more than 50 genes that are required for specific stages in the cell division cycle and so were termed CDC genes. On average 4.6 alleles were identified for each CDC gene in the original screens, suggesting the number of cdc mutants that could be identified by this approach had reached a plateau (Hartwell et al., 1973
). However, several lines of evidence suggested that additional genes with cell cycle stage-specific functions remained to be identified (Hartwell et al., 1973
; Hartwell, 1974
; Pringle and Hartwell, 1981
). Indeed, cell cycle screens in other model organisms identified additional CDC genes among the essential genes (Nurse, 1975
; Nurse and Bisset, 1976
; Nasmyth and Nurse, 1981
), as have more recent screens using alternative strategies (Prendergast et al., 1990
; Stevenson et al., 2001
; Kanemaki et al., 2003
). Consistent with their roles in a biological process of fundamental importance, most of the CDC genes are conserved and appear to have conserved function in most eukaryotes, including humans.
Several approaches exist for studying the biological function of essential genes. Temperature-sensitive mutations have been used extensively in analysis of cell cycle genes, and many temperature-sensitive mutations lead to rapid depletion of the gene product being analyzed. However, temperature-sensitive mutations are difficult to construct in a systematic manner, and the molecular basis for temperature sensitivity is, in most cases, unknown. Induced proteolysis through the creation of gene fusions to sequences encoding the "N-degron," a temperature-inducible proteolytic degradation signal (Dohmen et al., 1994
), has also been used to identify cell cycle genes (Kanemaki et al., 2003
). The N-degron can be applied systematically, and in many cases leads to rapid gene product depletion. However, in a recent study in which 104 essential genes were fused to the N-degron sequences, nearly 40% of essential genes fused to the N-degron did not result in inviability at the nonpermissive temperature (Kanemaki et al., 2003
), indicating that rapid protein depletion by the N-degron is not uniform across the proteome. Despite this limitation, functional information could be derived from strains in which depletion was incomplete (Aparicio, 2003
; Kanemaki et al., 2003
). Essential gene function can also be studied systematically by gene product depletion using a repressible promoter. Promoter replacement alleles allow the systematic analysis of essential genes, although, like N-degron fusions, the degree of gene product depletion varies from gene to gene, depending on both mRNA and protein half-life. With promoter replacement alleles each open reading frame (ORF) remains intact, and repression conditions with minimal effects on cell physiology can be chosen. In a previous study (Mnaimneh et al., 2004
), we described the construction of tetracycline-regulatable promoter (TetO7 promoter) alleles of
600 essential genes in S. cerevisiae. The TetO7 promoter collection has been used to probe essential gene function in cell size control, cell morphology, mitochondrial morphogenesis, and for gene expression and synthetic genetic interaction profiling (Mnaimneh et al., 2004
; Altmann and Westermann, 2005
; Davierwala et al., 2005
). In the present study, we have expanded the TetO7 promoter collection to encompass 773 essential genes, almost 75% of the essential gene set, making it the most complete resource for systematic analysis of essential gene function in yeast.
Many genes involved in cell growth, cell division, and cell cycle progression are indispensable for these processes and are therefore included among the
1050 yeast genes essential for viability (out of a total of
5800 yeast genes). Here, we utilized the TetO7 promoter collection to analyze essential gene function in cell division and cell cycle progression, using flow cytometric analysis to measure cellular DNA content after promoter shut-off. More than 65% of the strains displayed an altered flow cytometry profile after promoter shut-off, allowing categorization of essential genes on the basis of cell cycle profile. Our systematic analysis is not only useful in illustrating widespread contribution of essential genes to individual stages of the cell cycle, but is also useful for elucidating the functions of uncharacterized genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
1 leu2
0 met15
0. Primers for the PCR amplification for the DNA fragments containing the Tet-off promoter, KanR gene and partial sequences of the target genes and for confirmation of correct promoter replacement are available at http://hugheslab.med.utoronto.ca/Mnaimneh/data/Primers/Hughes_tet_ promoter_primers_Supp.xls. Two strategies were used for PCR confirmation of integrations. The first utilized primers to genomic DNA outside the intended integration site, in which a correct integration yields a band of
2.4 kb, an incorrect integration
400 bases. The second confirmation strategy utilized one primer complimentary to the cassette and one to the genomic DNA, in which a band of
700 bases represents an integration at the correct site. The smc4-1 strain (MATa smc4-1::kanMX his3
1 leu2
0 ura3
0 met15
0) contains the smc4-1 allele (Freeman et al., 2000
Flow Cytometry
For screening the collection of TetO7 promoter alleles, cells were grown at 30°C in the absence or presence of 10 µg/ml doxycycline for 15 h, and
1 x 107 cells for each strain were harvested and fixed with 70% ethanol. Cells were typically harvested at an OD600 of <1.0 and processed for flow cytometry as described (Davierwala et al., 2005
; Figure 1). The complete set of flow cytometry histograms, ordered by position on the array plates of the TetO7 promoter collection, is available at: http://biochemistry.utoronto.ca/brown/data.html.
|
Computational Analysis of Flow Cytometry Profiles
Normalization.
As a preprocessing step, we consecutively aligned the local maximum of all samples in each DNA content region. The alignment of flow cytometry data were done as follows: 1) select the rightmost channel left of the channel representing the DNA content region to align, such that most of the samples do not show DNA content at this channel (i.e., this channel is a sparse region of the flow cytometry matrix), and set this channel as the cutting point; 2) locate the local maximum in a fixed range (±50 channels) around the channel corresponding to the selected DNA content; 3) shift the vector of values starting at the cutting point per each sample so that its local maximum corresponds to the channel representing the selected DNA content; and 4) in case of left shifting discard any value left of the cutting point; otherwise, fill the created spaces with the value at the cutting point.
Classification. Profiles were assigned to six different abnormal categories according to the following criteria (thresholds were manually obtained based on examination of the data, i.e., the computational categorization was tuned to make the same decisions as manual inspection):
Functional Analysis.
For analysis of the functional enrichment of each class, we used the annotations from the Biological Process hierarchy of the Gene Ontology (GO) in the YEAST annotation package 1.10 of Bioconductor (Gentleman et al., 2004
). Before the statistical analysis, GO annotations were up-propagated. GO Categories with a hypergeometric p value smaller than 0.001 were considered to be significantly overrepresented in a class, because this resulted in zero or one false-positive categories on average for each of 10,000 sets of randomly chosen genes of the same sizes as those analyzed here. For visualization, similar significant GO categories were joined together under their most specific common ancestor.
G1 Intragroup Functional Analysis. Genes in the G1 class were first ordered according to the ratio of 1C accumulation peak to 2C accumulation peak and then separated in four equally sized sets. The functional enrichment of each set of genes was determined as described in Functional Analysis above. p values for lower quartiles were calculated excluding genes in upper quartiles.
Clustering. Genes in each of the six DNA content classes shown in the heatmaps (Figures 2A and 3A) were separately clustered using Pearson correlation as a similarity measure and hierarchical clustering with the complete linkage method. Clustering and heatmaps were generated in R 2.2.0. The genes in the cell size heatmap (Figure 3A) are shown in the same order as in the DNA content figure (Figure 2A).
|
|
Raw and normalized flow cytometry data, as well as data used to generate the graphical displays in Figures 2, 3, and 4, and the complete functional enrichment data are available at: http://biochemistry.utoronto.ca/brown/data.html.
|
0.5. Cultures corresponding to 3 OD of cells were then collected, and the cell pellets were resupended in 0.5 ml fresh YPD. Protein extracts were prepared by NaOH lysis and TCA precipitation, and western blotting with monoclonal anti-CPY or anti-ALP antibodies (Molecular Probes, Eugene, OR), or polyclonal anti-Gas1p antibody (kind gift from H. Riezman), as described (Davierwala et al., 2005
Chromatin Immunoprecipitation
Strains were grown in the presence of doxycycline for 5 h and synchronized with
-factor. G1-arrested cells were released into fresh media and
20 OD of cells were harvested every 12 min. Cells were fixed and chromatin immunoprecipitates were prepared essentially as described (Kuras and Struhl, 1999
), using IgG-agarose to immunoprecipitate TAP-tagged Mcm2 from wild-type and TetO7-SMC4 strains. Immunoprecipitated DNA was analyzed by PCR using four primer pairs for specific regions at and around ARS1 as described (Tanaka et al., 1997
). All PCR, gel electrophoresis, and visualization conditions were as described (Tanaka et al., 1997
), and the gel was quantified using ImageQuant.
| RESULTS |
|---|
|
|
|---|
A Cell Cycle Screen of the Titratable Promoter Alleles of the Essential Genes
We used the collection of TetO7 promoter alleles to systematically identify genes that cause a change in cell cycle profile after promoter shut-off. Because the nuclear DNA content of a cell reveals its position in the cell cycle, we measured the DNA contents of populations of cells using flow cytometry, both before and after turning off transcription of each essential gene (Figure 1A). The resulting histograms can be used to determine the fraction of cells in a population that are at each cell cycle stage (G1, S, and G2), and accumulation of cells with a specific DNA content indicates a defect in progression through a particular cell cycle stage. In initial experiments with four of the TetO7 promoter replacement strains, samples were analyzed by flow cytometry at different times after the addition of doxycycline (Supplementary Figure 1). We found that for the strains analyzed, alterations in the flow cytometry histograms were first evident at varying times, from 2 h in the case of CDC27, to 8 h in the case of SLD5. Importantly, the same category of cell cycle accumulation was still evident at the latest time point, 15 h, and in the case of the S phase gene SLD5 the accumulation became increasingly obvious and therefore easier to score at the later points. Thus, in order to detect changes in as many strains as possible, we chose to analyze the promoter replacement collection after 15 h of growth in the presence of doxycycline.
Scoring the Flow Cytometry Profiles
We scored the flow cytometry profiles using both manual inspection and computational analysis (see Materials and Methods). In both cases, the flow cytometry histograms were categorized according to the distribution of DNA contents (Figure 1B). The categories were as follows: strains that accumulate with 1C DNA content (and so are arrested or delayed in G1); strains that accumulate with DNA content between 1C and 2C (and so are delayed in S phase); strains that accumulate with 2C DNA content (and so are arrested or delayed in G2); strains that accumulate with 3C and/or 4C DNA content (an unexpected cell cycle profile, likely relating to defects in cell separation); and strains that were diploid; and strains exhibiting a normal cell cycle profile. In the manual analysis some profiles did not exhibit the characteristics of any of the categories, yet were clearly abnormal. These were grouped together as "other." In the computational analysis the "other" category is made up of strains that showed accumulation of cells with less than 1C DNA contents. Strains that displayed accumulations at more than one DNA content were placed in more than one category, and all strains were placed in at least one category.
Manual inspection indicated that of 838 TetO7 promoter mutants screened, 625 (or 75%) displayed an altered flow cytometry profile. Computational analysis (see Materials and Methods for details) indicated that 579 displayed an altered profile, and 563 mutants (or 67%) were scored as altered by both methods. On average,
80% of assessments made by manual inspection were supported by computational analysis and vice versa. The fraction of genes in each category (supported by both manual and computational analysis) is shown in Figure 1B. Examples of the strongest cell cycle phenotypes exhibited by TetO7 promoter mutants are shown in Figure 1C. However, we scored even modest alterations in flow cytometry profile. For example, in the 1C accumulation category, modeling of the flow cytometry histograms indicated that the fraction of cells with a 1C DNA content ranged from 31 to 97% for mutants scored as 1C, compared with 28% for the wild-type control (see below).
In TetO7 promoter mutant strains the onset of growth cessation varies from gene to gene because of the time required for depletion of the transcript and protein. Even though we cannot be certain of the timing or degree of depletion of the different gene products after promoter shut-off for 15 h, we saw clear and specific effects on cell cycle progression with the majority of the essential genes. The accumulation of specific DNA contents suggests an intimate connection between essential gene function and cell cycle progression. The altered flow cytometry profiles did not appear to reflect cell debris, with the possible exception of some profiles classified as "other," indicating that our results were not artifacts due to simple slow growth or cell death. The flow cytometry profiles of the eight wild-type colonies we analyzed were scored as normal, and most of the TetO7 promoter strains with little or no growth defect on plates also had no perturbation of DNA content (112 of 190 tet promoter mutants with normal growth were scored by both methods as having normal flow cytometry profiles), indicating that our scoring of abnormal flow cytometry profiles was not overly liberal. Strains that did not show a growth defect in the presence of doxycycline likely represent strains in which gene product depletion is insufficient to cause a phenotype, perhaps because of long mRNA or protein half-life. However, 51 strains displayed normal growth and yet had altered cell cycle profiles in both manual and computational scoring (Supplementary Figure 2), indicating that in some instances the screen was sensitive enough to identify cell cycle delays in the absence of a growth defect. Conversely, of the TetO7 promoter strains with normal cell cycle profiles, almost half (82 of 194) exhibited growth defects in the presence of doxycycline. Therefore not all essential genes are linked to cell cycle processes; these 82 strains presumably cease growth without cell cycle stage specificity.
As a measure of the efficiency of the screen in identifying genes with roles in cell division cycle progression, we determined the fraction of known CDC genes that were scored as positive (Supplementary Figure 3). There are 67 CDC genes in the Saccharomyces Genome Database (SGD). Of these, 6 have not been mapped, and 16 are not present in the TetO7 promoter collection. Of the remaining 45, 89% displayed an altered cell cycle profile that was scored manually, computationally, or by both methods. The cell cycle arrest point for 38 of these CDC genes is known or can reasonably be inferred. In our screen, the majority (28 of 38; 74%) showed the expected cell cycle stage-specific accumulation. For example, CDC6, CDC8, CDC9, CDC17, CDC21, CDC45, CDC47, CDC101, CDC102, and CDC105 all displayed accumulation of S phase cells, in agreement with their known functions in DNA replication (Hartwell, 1971a
, 1973
, 1976
; Johnston and Nasmyth, 1978
; Piatti et al., 1995
; Zou et al., 1997
; Burgers, 1998
; Kanemaki et al., 2003
). In addition, CDC4, CDC25, and CDC64 are all required for G1 transit and accumulated cells with 1C DNA contents (Bedard et al., 1981
; Tripp and Pinon, 1986
; Goh and Surana, 1999
). We conclude that neither the extended period of gene product depletion nor the stringency of the categorization cutoffs precluded accurate categorization of the mutants and that the number of false negatives in the screen was low, because cell cycle defects were detected in most TetO7 promoter mutants of known CDC genes.
Overall, our analysis indicated that most essential genes were linked to a specific stage of cell cycle progression. We next asked whether there are relationships between categories of DNA content and other gene properties, and if so, whether they are illuminating with regard to either the nature of the cell cycle or the functions of the individual genes.
G1 Phase
For visualization, flow cytometry histograms from the overlap between manual and computational scoring methods were grouped and subjected to hierarchical clustering (Figure 2A). We then asked if any GO functional annotations were enriched in any of the flow cytometry categories. These data are displayed as a "GO-gram" (Figure 2B), in which the annotations of genes within the categories are displayed; red bars indicate that a given GO annotation is significantly enriched in a given flow cytometry category. The 1C accumulation category, which at 270 genes was the largest, showed a dramatic enrichment for genes involved in protein synthesis (Figure 2B). This included genes involved in ribosome biogenesis, translation, tRNA metabolism, tRNA charging, and RNA polymerase III transcription (GO:0042254, ribosome biogenesis and assembly, p value 2.79e-23; GO:0043037, translation, 1.47 e-08; GO:0006399, tRNA metabolism, 6.70e-07; GO:0043038, amino acid activation, 3.60e-06; and GO:0006383, transcription from RNA polymerase III promoter, 0.0002). Although ribosome biogenesis is not cell cycle regulated in yeast (Bernstein and Baserga, 2004
), there is some prior evidence that overexpression or depletion of ribosomal proteins and factors involved in ribosome biogenesis causes cell cycle defects. For example, depletion of SSU processome proteins leads to G1 arrest due to lack of ribosomes (Bernstein and Baserga, 2004
), as does depletion of the 20S pre-rRNA maturation factor Rio1 (Angermayr et al., 2002
). Similarly, tRNA charging has long been linked to G1 transit as a number of mutants in tRNA synthetases arrest in G1 (Unger and Hartwell, 1976
; Bedard et al., 1981
). Consistent with existing links between protein synthesis and G1 transit, our data indicate that the most dramatic restriction on a cell's ability to progress past the cell cycle commitment point in G1 (known as Start in yeast) is the capacity for protein synthesis, as inhibiting protein synthesis by a range of different means caused accumulation of cells in G1. This observation is reminiscent of previous analyses of cell size, which indicated that cells mutated in ribosome biogenesis factors tend to be small and unbudded (Jorgensen et al., 2002
; Mnaimneh et al., 2004
). However, there does not appear to be a direct relationship between cell size and our categorized flow cytometry profiles. When we grouped the cell size data from Mnaimneh et al. (2004)
according to the assigned flow cytometry categories (Figure 3), we found that although strains in the 1C category tended to be slightly smaller than normal strains, there was not a strong predictive correspondence between the flow cytometry data and the cell size data. Thus, cell cycle accumulation in G1 was far from an absolute predictor of cell size (and vice versa).
Eleven unannotated genes gave rise to G1 arrest after promoter shut-off, suggesting they might be involved in some aspect of ribosome biogenesis or protein biosynthesis. Indeed, high-throughput data suggest several of these genes perform such functions: YPR169W encodes a nucleolar protein of unknown function, which interacts with Tpt1p, a protein involved in tRNA splicing (Hazbun et al., 2003
; Huh et al., 2003
); Ynl313cp is localized to both cytoplasm and nucleus and interacts with proteins involved in RNA processing, ribosomal biogenesis import, and protein metabolism (Gavin et al., 2002
; Huh et al., 2003
; Krogan et al., 2004
); Ynl310cp encodes a protein that interacts with proteins involved in RNA processing (Gavin et al., 2002
; Ho et al., 2002
; Hazbun et al., 2003
); and Yhr020wp is inferred to have tRNA aminoacylation and ligation activity and also affinity-precipitates with Utp13, which is part of the SSU processome involved in pre-18S rRNA processing (Tatusov et al., 2000
; Ho et al., 2002
). Finally, Ynr046wp has been recently identified as subunit of tRNA methyltransferase (Purushothaman et al., 2005
), and Ynr054cp, which encodes a nucleolar protein that interacts with proteins involved in RNA processing, has been recently characterized as a component of the SSU processome (Gavin et al., 2002
; Ho et al., 2002
; Huh et al., 2003
; Hoang et al., 2005
). Thus at least 6 of the 11 unannotated ORFs in the G1 arrest category are likely to function in ribosome biogenesis, tRNA metabolism, or protein metabolism.
Histograms for mutants in the 1C accumulation category revealed that the category encompassed a range of degrees of 1C DNA content accumulation, from 31 to 97% 1C (compared with an average wild-type value of 28% 1C). Although this threshold was not exceeded by any wild-type isolates analyzed (unpublished data), we nonetheless considered that our threshold might have been too lenient, in which case there should be little relationship between growth effect or functional category of the mutated gene for those with relatively low 1C accumulation. We therefore ordered the mutants in the 1C category according to the relative 1C accumulation and asked whether severity of phenotype (in this case 1C accumulation) was related to particular gene functions and growth defect (Figure 4). We found that in the top quartile there was a significant enrichment for annotations related to tRNA synthesis and activation, and translation, whereas the second and third quartiles were enriched for functions related to ribosome biogenesis. The fourth quartile was significantly populated by genes in the broader RNA metabolism category. Additionally, there was a tendency for the mutants with highest 1C accumulation to display a more prominent growth defect. This suggests that given additional phenotypic data, we may be able to understand why some of the tet-promoter mutants display greater growth phenotypes than others; we had previously not been able to draw strong conclusions in this regard (Mnaimneh et al., 2004
). In summary, functional information could be derived from a single flow cytometry parameter (1C:2C ratio) in addition to the information contained in the different categories. Moreover, our choice of threshold was not overly liberal, because there was functional information contained even in the lowest quartile.
S Phase
The S phase accumulation class, which comprises only 27 genes (Figure 5), showed a dramatic enrichment for genes involved in DNA replication (GO:0006260, DNA replication, p value 4.02e-17; GO:0045005, maintenance of fidelity during DNA-dependent DNA replication, 9.97e-05; GO:0006281, and DNA repair, 1.20e-07; Figure 2B). These include CDC6, SLD5, PSF2, POL1, POL2, PRI2, POL30, POL12, DPB11, DPB2, PSF1, and CDC45, which all have direct roles in the initiation and/or the elongation of DNA replication. Several of the genes in this class, although not directly involved in replication, have known roles in nucleotide metabolism (CDC8, CDC21, DFR1, FOL 2, and DUT1). We explored the possibility that CSE1 and SMC4 could play roles in regulating DNA replication. CSE1 is responsible for the nuclear shuttling of the nuclear transporter importin
(Hood and Silver, 1998
; Kunzler and Hurt, 1998
; Solsbacher et al., 1998
) and a role for CSE1 in mitosis (but not in S phase) has been described (Xiao et al., 1993
; Schroeder et al., 1999
). We synchronized cells in G1, released them into the cell cycle after shut-off of TetO7-CSE1, and measured the DNA contents of the cells by flow cytometry (Figure 6A). Cells depleted of CSE1 failed to undergo DNA replication, with 70% of cells still containing 1C DNA contents 90 min after release from G1. We examined whether these cells had passed Start by measuring the fraction of cells that had elicited buds, a marker of cell cycle progression through the G1/S boundary that is independent of the initiation of DNA replication (Figure 6C). At 90 min almost 90% of the TetO7-CSE1 cells had budded (51% had a large bud, and 38% had a small bud), indicating that most cells had passed Start. These data indicate that categorization of TetO7-CSE1 in the S phase accumulation class in the initial screen accurately predicted a role for CSE1 in S phase, either in the initiation of DNA replication or after initiation but before significant DNA synthesis has occurred. Additionally, the depletion of Cse1p was robust enough to observe a clear phenotype in the first cell cycle. Thus, CSE1 is required for DNA replication, suggesting that an essential DNA replication protein is an importin
cargo.
|
|
-factor displayed a shmoo morphology vs. 98% of wild-type cells), suggesting that the 33% of TetO7-SMC4 cells that did not enter S phase during the course of the experiment were not cells that had failed to arrest in G1. Finally, we observed a similar phenotype with the temperature-sensitive smc4-1 strain, in which 35% of cells failed to enter S phase (Figure 6B). As with CSE1, categorization of SMC4 in the initial cell cycle screen accurately predicted a role for SMC4 in entry into S phase.
We explored the role of SMC4 in DNA replication by observing the association of the Mcm2 protein with replication origins and replication forks during S phase, using chromatin immunoprecipitation (ChIP; Figure 7, A and B). In wild-type cells, Mcm2p was preferentially associated with the replication origin ARS1 in G1. As cells progressed into S-phase the association of Mcm2 with the origin decreased and the association of Mcm2p with origin-distal DNA fragments increased. This pattern implies association with the replication fork, as has been previously described (Aparicio et al., 1997
; Tanaka et al., 1997
). We found that when TetO7-SMC4 was shut off, Mcm2p binding to ARS1 decreased as cells progressed into S phase, but Mcm2p failed to then associate with origin-distal regions. Thus, our data indicate that Smc4p depletion affected the binding of Mcm2p to the regions flanking ARS1 during S phase and suggest a role for Smc4p in regulating DNA replication, perhaps in the transition from initiation to elongation. Previous work has hinted at a role for condensin in S phase, including the hydroxyurea sensitivity of condensin mutants in fission yeast (Aono et al., 2002
) and the presence of condensin at the rDNA locus during S phase (Johzuka et al., 2006
). Additionally, genome-wide protein-DNA interaction data suggests that Smc4p is excluded from replication origins (Wang et al., 2005
), although this analysis was not performed during synchronous progression through S phase and therefore may not have detected the relevant localization. Perhaps Smc4p exerts its effect on DNA replication by creating a chromatin environment that facilitates progression of replication forks away from replication origins. We were unable to detect S phase defects for the other condensin subunit genes in the TetO7-promoter collection, SMC2, YCG1, and YCS4. It will be of great interest to determine whether this S phase function is a property of Smc4p or of the condensin complex.
|
3C/4C
The 3C/4C accumulation flow cytometry category showed a clear enrichment for genes associated with protein transport, protein secretion, and glycoprotein metabolism (GO:0016192, vesicle-mediated transport, P value 3.38e-08; GO:0046903, secretion, 1.91e-10; and GO:0009100, glycoprotein metabolism, 4.04e-08; Figure 2B and unpublished data). An enrichment of similar functions was also evident in the category of mutants exhibiting multiple buds in our earlier morphological analysis of the TetO7 promoter collection (Mnaimneh et al., 2004
), although far fewer genes were scored in this category. Four uncharacterized genes, YJL097W/PHS1, YDL193W/NUS1, YNL149C, and YML125C, were included in the 3C/4C category, suggesting they might be involved in some aspect of protein trafficking. The proteins encoded by these four genes all localize to the ER or the ER and vacuole (Huh et al., 2003
) and the YJL097W/PHS1 gene was recently found to be a regulator of sphingosine levels (Schuldiner et al., 2005
). We examined the cellular morphology of these strains after promoter shut-off and found that all exhibited a multiple-budded or large-budded phenotype with normal postmitotic nuclei (Figure 8A), indicating a defect in cytokinesis or in cell separation. Normal nuclear division accompanied by a cytokinesis or cell separation defect likely resulted in the appearance of events with apparent 3C and 4C DNA contents in the flow cytometry analysis. We also examined the TetO7 alleles of YJL097W/PHS1, YDL193W/NUS1, YNL149C, and YML125C for defects in the processing and trafficking of three different glycosylated proteins (Figure 8B). TetO7-YJL097W did not show detectable defects in carboxypeptidase Y (CPY) processing, indicating that it is not essential for ER-to-Golgi transport, N-linked glycosylation, or vacuolar signal peptide cleavage; however, it displayed abnormal processing of alkaline phosphatase (ALP), accumulating the soluble form of ALP. The presence of soluble form ALP is reminiscent of apl5
mutants (Cowles et al., 1997
) and of depletion of the recently described PGA1 and suggests a defect in vacuolar transport. Additionally, TetO7-YJL097W showed defects in the accumulation of mature Gas1p. TetO7-YDL193W accumulated hypoglycosylated forms of CPY, much like TetO7-SEC59, suggesting a role in glycosylation in the ER. TetO7-YDL193W also accumulated immature ER forms of Gas1p. TetO7-YNL149C and TetO7-YML125C both displayed normal CPY processing, but accumulated increased levels of abnormal soluble ALP and slightly decreased levels of mature Gas1p. Together these results indicate roles for YJL097W, YDL193W, YNL149C, and YML125C in protein trafficking. On the basis of these observations we have named the YNL149C and YML125C genes PGA2 and PGA3, for processing of Gas1p and ALP.
|
| DISCUSSION |
|---|
|
|
|---|
It is important to note that the present screen differs from the classic cell cycle screens in that we did not require arrest of a strain at a particular cell cycle stage, but instead scored any significant cell cycle stage-specific accumulation. In this respect our screen is conceptually similar to a recent cell cycle screen in Drosophila cells using gene product depletion by RNA-mediated interference (RNAi; Bjorklund et al., 2006
). One reasonable explanation for the observation that in some strains the cell cycle phenotype is a delay rather than an arrest and that in some strains <100% of cells accumulate at the relevant cell cycle stage is that promoter shut-off will not uniformly result in protein depletion within the time-frame of the experiment. Although this might be perceived as a limitation, we were able to detect cell cycle defects in 89% of known CDC genes as TetO7 alleles and correctly categorize 74% of these, suggesting that the false negative rate in our screen was low. Additionally, the functional enrichment within the 1C and S phase categories was consistent with published work: the 1C category displayed a massive enrichment for genes involved in tRNA and ribosome biogenesis, and the S phase category displayed a massive enrichment for genes involved in DNA replication. Finally, one of the hallmarks of cell cycle arrest is continued cell growth in the absence of cell cycle progression. Strains in the S-phase, 2C, and 3C/4C categories all exhibited larger than normal cell size (Figure 3), consistent with cell growth despite cell cycle arrest. By contrast, strains in the normal cell cycle category exhibited wild-type cell sizes, despite almost half of the strains in the category having clear growth defects after promoter shut-off. Thus we conclude that measuring DNA contents after promoter shut-off of essential genes represents a sensitive and accurate means of identifying cell cycle mutants.
A second major conclusion is that down-regulation of genes with diverse functions in regulating and maintaining the protein biosynthetic capacity of the cell causes arrest in G1. The intimate relationship between protein biosynthesis and cell cycle progression through G1 that we observed is consistent with previous observations linking protein synthesis with G1 transit (Unger and Hartwell, 1976
; Bedard et al., 1981
; Moreno and Nurse, 1994
), with connections between ribosome biogenesis and cell size regulation at Start (Jorgensen et al., 2002
, 2004
), and with links between translation rate and G1 (Polymenis and Schmidt, 1997
) and provides a global view of cellular processes that contribute to passage through the cell cycle commitment point in G1, Start. There is concordance between our G1 data and that derived from cell cycle screening of Drosophila cells after gene product depletion by RNAi (Bjorklund et al., 2006
) in that both screens show clear involvement of protein biosynthesis pathways in G1 transit. However, in the Drosophila screen the genes identified were largely ribosomal subunit genes, whereas our data also illustrate the roles of ribosome assembly factors, and tRNA synthesis, processing, and charging proteins in regulating G1. Given that down-regulation of ribosome biogenesis can cause cell cycle arrest in G1 in mouse and human cells (Strezoska et al., 2000
; Volarevic et al., 2000
; Pestov et al., 2001
; Oliver et al., 2004
), our identification of essential genes that regulate G1 in yeast could be useful in dissecting connections between ribosome biogenesis and cell cycle progression in mammalian cells.
A third conclusion is that sorting genes into different categories on the basis of the flow cytometry profiles allowed prediction of gene function. We derived information about gene function within the flow cytometry categories by identifying biological processes that are enriched within each category. Although in the present screen the cell cycle accumulation stage may not necessarily reflect the execution point or the point of action of the genes, we identified significant functional enrichments within each cell cycle category and so were able to make predictions about the biological function of uncharacterized essential genes. In particular, we found that accumulation of 3C and 4C DNA content reliably identified components of protein trafficking, modification, and secretion pathways. Four uncharacterized genes in this category, NUS1, PHS1, PGA2, and PGA3, had detectable changes in the modifications of one or more of the markers of protein secretion and modification analyzed, indicating roles in protein trafficking. The 1C category predicted function in protein biosynthesis: 6 of the 11 unannotated genes in the G1 accumulation category have connections to protein biosythesis suggested by other functional genomic data, including proteinprotein interaction and protein localization data. Additionally, we found that depletion of the condensin subunit Smc4p or the importin
export receptor Cse1p resulted in defective DNA replication and therefore that accumulation of S phase DNA contents accurately predicted function in DNA synthesis. Thus, despite the possibility that some of the observed cell cycle delays could have been due to indirect effects of protein depletion, this comprehensive data set was useful for prediction of essential gene function.
Our study complements several decades of cell cycle experimentation, and provides a systematic survey of essential gene function in the cell division cycle. Regulation of proliferation is of primary importance in the study of human development and in the study of many human diseases. We anticipate that many of the connections between essential cellular processes and cell division cycle progression illustrated here will be conserved in humans. Indeed, a number of pathways that we find are important for cell cycle progression in yeast are also important for cell cycle progression in Drosophila (Bjorklund et al., 2006
), encouraging the view that this global analysis of the connections between essential genes and the cell cycle will be useful in dissecting the network of processes that impinge on progression through the cell cycle in metazoans.
| ACKNOWLEDGMENTS |
|---|
| Footnotes |
|---|
This article was published online ahead of print in MBC in Press (http://www.molbiolcell.org/cgi/doi/10.1091/mbc.E06-04-0368) on August 30, 2006.
Address correspondence to: Grant W. Brown (grant.brown{at}utoronto.ca) or Timothy R. Hughes (t.hughes{at}utoronto)
| REFERENCES |
|---|
|
|
|---|
Angermayr, M., Roidl, A., Bandlow, W. (2002). Yeast Rio1p is the founding member of a novel subfamily of protein serine kinases involved in the control of cell cycle progression. Mol. Microbiol 44, 309324.[CrossRef][Medline]
Aono, N., Sutani, T., Tomonaga, T., Mochida, S., Yanagida, M. (2002). Cnd2 has dual roles in mitotic condensation and interphase. Nature 417, 197202.[CrossRef][Medline]
Aparicio, O. M. (2003). Tackling an essential problem in functional proteomics of Saccharomyces cerevisiae. Genome Biol 4, 230.[CrossRef][Medline]
Aparicio, O. M., Weinstein, D. M., Bell, S. P. (1997). Components and dynamics of DNA replication complexes in S. cerevisiae: redistribution of MCM proteins and Cdc45p during S phase. Cell 91, 5969.[CrossRef][Medline]
Bedard, D. P., Johnston, G. C., Singer, R. A. (1981). New mutations in the yeast Saccharomyces cerevisiae affecting completion of "Start.". Curr. Genet 4, 204214.
Bernstein, K. A. and Baserga, S. J. (2004). The small subunit processome is required for cell cycle progression at G1. Mol. Biol. Cell 15, 50385046.
Bjorklund, M., Taipale, M., Varjosalo, M., Saharinen, J., Lahdenpera, J., Taipale, J. (2006). Identification of pathways regulating cell size and cell-cycle progression by RNAi. Nature 439, 10091013.[CrossRef][Medline]
Brachmann, C. B., Davies, A., Cost, G. J., Caputo, E., Li, J., Hieter, P., Boeke, J. D. (1998). Designer deletion strains derived from Saccharomyces cerevisiae S288C: a useful set of strains and plasmids for PCR-mediated gene disruption and other applications. Yeast 14, 115132.[CrossRef][Medline]
Burgers, P. M. (1998). Eukaryotic DNA polymerases in DNA replication and DNA repair. Chromosoma 107, 218227.[CrossRef][Medline]
Cowles, C. R., Odorizzi, G., Payne, G. S., Emr, S. D. (1997). The AP-3 adaptor complex is essential for cargo-selective transport to the yeast vacuole. Cell 91, 109118.[CrossRef][Medline]
Culotti, J. and Hartwell, L. H. (1971). Genetic control of the cell division cycle in yeast. III. Seven genes controlling nuclear division. Exp. Cell Res 67, 389401.[CrossRef][Medline]
Davierwala, A. P., et al. (2005). The synthetic genetic interaction spectrum of essential genes. Nat. Genet 37, 11471152.[CrossRef][Medline]
Dohmen, R. J., Wu, P., Varshavsky, A. (1994). Heat-inducible degron: a method for constructing temperature-sensitive mutants. Science 263, 12731276.
Freeman, L., Aragon-Alcaide, L., Strunnikov, A. (2000). The condensin complex governs chromosome condensation and mitotic transmission of rDNA. J. Cell Biol 149, 811824.
Gavin, A. C., et al. (2002). Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 415, 141147.[CrossRef][Medline]
Gentleman, R. C., et al. (2004). Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5, R80.[CrossRef][Medline]
Goh, P. Y. and Surana, U. (1999). Cdc4, a protein required for the onset of S phase, serves an essential function during G(2)/M transition in Saccharomyces cerevisiae. Mol. Cell. Biol 19, 55125522.
Hartwell, L. H. (1971a). Genetic control of the cell division cycle in yeast. II. Genes controlling DNA replication and its initiation. J. Mol. Biol 59, 183194.[CrossRef][Medline]
Hartwell, L. H. (1971b). Genetic control of the cell division cycle in yeast. IV. Genes controlling bud emergence and cytokinesis. Exp. Cell Res 69, 265276.[CrossRef][Medline]
Hartwell, L. H. (1973). Three additional genes required for deoxyribonucleic acid synthesis in Saccharomyces cerevisiae. J. Bacteriol 115, 966974.
Hartwell, L. H. (1974). Saccharomyces cerevisiae cell cycle. Bacteriol. Rev 38, 164198.
Hartwell, L. H. (1976). Sequential function of gene products relative to DNA synthesis in the yeast cell cycle. J. Mol. Biol 104, 803817.[CrossRef][Medline]
Hartwell, L. H., Culotti, J., Reid, B. (1970). Genetic control of the cell-division cycle in yeast. I. Detection of mutants. Proc. Natl. Acad. Sci. USA 66, 352359.
Hartwell, L. H., Mortimer, R. K., Culotti, J., Culotti, M. (1973). Genetic control of the cell division cycle in yeast. V. Genetic analysis of cdc mutants. Genetics 74, 267286.
Hazbun, T. R., et al. (2003). Assigning function to yeast proteins by integration of technologies. Mol. Cell 12, 13531365.[CrossRef][Medline]
Ho, Y., et al. (2002). Systematic identification of protein complexes in Saccharomyces cerevisiae by mass spectrometry. Nature 415, 180183.[CrossRef][Medline]
Hoang, T., Peng, W. T., Vanrobays, E., Krogan, N., Hiley, S., Beyer, A. L., Osheim, Y. N., Greenblatt, J., Hughes, T. R., Lafontaine, D. L. (2005). Esf2p, a U3-associated factor required for small-subunit processome assembly and compaction. Mol. Cell. Biol 25, 55235534.
Hood, J. K. and Silver, P. A. (1998). Cse1p is required for export of Srp1p/importin-alpha from the nucleus in Saccharomyces cerevisiae. J. Biol. Chem 273, 3514235146.