Previous Article | Next Article ![]()
Eukaryotic Cell, August 2006, p. 1337-1346, Vol. 5, No. 8
1535-9778/06/$08.00+0 doi:10.1128/EC.00101-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
School of Life Sciences, Södertörns Högskola, and Department of Biosciences and Nutrition, Karolinska Institute, SE-141 89 Huddinge, Sweden
Received 12 April 2006/ Accepted 18 June 2006
|
|
|---|
|
|
|---|
Genes differ in their requirement for coactivators, such as Gcn5, during their activation. For example, Gcn5 is required for activation of the Saccharomyces cerevisiae PHO8 gene but not the PHO5 gene, even though both are strongly activated by the Pho4 activator protein during phosphate starvation (17). The ability of activator proteins to recruit different repertoires of coactivator proteins to different target genes may be a characteristic feature of their activation domains. The activation domains of many activators show a low propensity for secondary and tertiary structure formation but adopt a folded conformation during interaction with target proteins, such as coactivators (11, 16). The structural flexibility of activator proteins provides potential for flexible recruitment of coactivators, and thus, consistent with experimental evidence, the identity of activator proteins cannot be expected to predict the coactivators they recruit to different target genes.
Various approaches have been taken to understand why some genes require the SAGA complex for activation while others do not. Evidence from experiments using fluorescence resonance energy transfer has shown that the SAGA complex is an important target for the S. cerevisiae Gal4 activator protein in vivo (3). DNA microarray studies in S. cerevisiae suggest that the SAGA complex is redundant with the TFIID coactivator complex (22). More recently, it has been suggested that the SAGA complex is important for the activation of highly induced stress-regulated genes while the TFIID complex plays a role in the expression of housekeeping genes whose expression level does not change greatly under different growth conditions (20). Another recent study shows that the Gcn5 component of the SAGA complex is recruited to essentially all expressed genes with an efficiency corresponding to their expression level (27).
Fission yeast exhibits several chromatin features that distinguish it from S. cerevisiae and which are similar to chromatin in mammals. Most notably, the fission yeast genome contains larger regions of heterochromatin that are established via binding of the HP1-like Swi6 protein to methylated histone H3 (lysine 9) and a mechanism involving the RNA interference machinery (8). In Schizosaccharomyces pombe, Gcn5 is involved in histone acetylation and chromatin remodeling events that are important for mitotic and meiotic recombination (32), but its role in gene regulation has not been studied previously.
S. cerevisiae contains several mitogen-activated protein (MAP) kinase signal transduction pathways that together contribute to the specificity of stress responses. S. pombe uses only one MAP kinase pathway and selectivity of the stress response depends on Tup1-like corepressors that mediate stress-specific changes in chromatin structure (18). Indeed, gene duplication has lead to the acquisition of two Tup1-like proteins in S. pombe (Tup11 and Tup12), one of which plays a specific role in the stress response (10). Given the chromatin differences between S. cerevisiae and S. pombe, the apparent key role of Gcn5 in the stress response of S. cerevisiae, and the differences in the mechanisms regulating stress in the two yeasts, S. pombe provides an interesting comparative model system for studies of Gcn5. Here we show that Gcn5 is not important for the growth of S. pombe under normal conditions but that it does play a critical role in the response of cells to specific stress conditions. We describe the use of the physiological adaptation to one of these stresses (KCl) as a model to understand how Gcn5 contributes to reprogramming genome expression in response to external signals.
|
|
|---|
For growth assays strains were cultivated in YEA medium (0.5% yeast extract, 3% glucose, 0.2% Casamino Acids with 100 mg/liter of adenine, uracil, and leucine). For cell plating assays, YEA medium alone or supplemented with KCl (1 M), CaCl2 (0.1 M), sorbitol (1.2 M), NaCl (0.18 M), MnCl2 (0.25 mM), or FeCl3 (0.1 mM) was used. For microarray experiments cells were cultivated in YES medium (31) and YES medium supplemented with 1 M KCl. The parental strain Fy368 was used as a wild-type control for all experiments.
Microarray analysis. For expression profiling of the gcn5 strain under nonstress growth conditions three biological replicates were used. For KCl stress experiments, at least two biological replicates were used. Cells were grown at 25°C to a cell density of 1 x 107 to 2 x 107 cells/ml with the cultures split in two. To one half, YES medium containing 2 M KCl was added to a final concentration of 1 M KCl, and the cells were frozen in liquid nitrogen after 15 or 60 min. The other half of the cultures were used as controls and immediately frozen in liquid nitrogen. In the case of the gcn5 strain, the cells were subjected to KCl for 60 min. RNA was collected as previously described (31). For each experiment, about 25 µg of total RNA was subjected to reversed transcription (catalog no. 11904-018; Invitrogen) and labeled with Cy3- or Cy5-labeled dCTP (Amersham). Reverse transcription primers (Eurogentec SA,) include an S. pombe-specific primer mix plus oligo(dT)18-21 and T20VN anchor primers. Labeled probes were hybridized to DNA microarrays containing S. pombe gene probes spotted in duplicate (31) (Eurogentec SA, Seraing, Belgium). Hybridization was performed as described previously (31) except that the hybridization volume was 72 µl. Slides were scanned using a Scan Express laser scanner and quantified with Imagene 4.2 software. Data were analyzed and normalized by the Lowess per-spot per-chip method using Genespring software (Silicon Genetics). Genes with changed expression were identified by selecting genes for which the mean ratio (sample to reference) reproducibly exceeded a chosen threshold value. Reproducibility was assessed using Student's t test to assess the significance of the difference of the mean ratio from a ratio of 1, expected for nonregulated genes. GeneSpring software was also used to calculate overlaps between gene groups using Fisher's exact test, for which P values were corrected for multiple testing (Bonferroni). Overlap between two groups of genes with a P value higher than 0.05 was not considered significant.
Promoter sequence analysis. RSAT software (http://rsat.ulb.ac.be/) was used to analyze promoter sequences. We used the Pattern Discovery tool to look for overrepresented sequences of six to eight nucleotides in the 800-bp upstream region of Gcn5-dependent KCl-induced genes. The Pattern Matching tool was then used to look at the occurrences of these sequence motifs as well as known DNA sequence motifs in different groups of genes selected according to their change in expression in response to KCl stress in the wild-type and gcn5 mutant strains. As controls, nine groups of 30 genes were selected at random; in addition, the 800-bp upstream sequence of each S. pombe gene was also used.
Gene ontology terms.
GoMiner (http://discover.nci.nih.gov/gominer/) was used to find enriched gene ontology (GO) terms of groups of genes based on their expression profiles during the KCl stress response. The different groups of genes were compared to all genes with known GO terms. Only GO terms that were significantly enriched (P
0.05) were included in the results.
Microarray data accession number. All microarray data are available at Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) under accession number GSE5227 [NCBI GEO] .
|
|
|---|
![]() View larger version (51K): [in a new window] |
FIG. 1. Gcn5 is not required for normal growth. wt, wild type. (A) Fivefold dilutions of wild-type and gcn5 strains were spotted on rich medium (YEA) and the plates incubated at 30°C for 3 days. (B) Growth rate of log-phase wild-type (filled points) and gcn5 (clear points) cells in liquid culture (mean relative densities of three cultures ± standard deviations are plotted). (C) There is no morphological difference between wild-type and gcn5 cells in log phase. Cells were stained with Hoechst stain to visualize nuclear DNA. (D) gcn5 deletion causes small effects on gene expression. A scatter plot shows the mean signal intensity values of three independent microarray experiments. The flanking diagonal lines indicate a threshold value of a twofold change in expression (see Table S1 in the supplemental material for gene names, mean changes in induction [n-fold], and P values of the 260 genes with a reproducible change of >1.5-fold).
|
|
View this table: [in a new window] |
TABLE 1. Genes with changed expression in the gcn5 mutant
|
![]() View larger version (74K): [in a new window] |
FIG. 2. Gcn5 is required for growth under specific environmental stress conditions. The gcn5 mutant spotted on YEA plates at 36°C or YEA plates supplemented with 1 M KCl or 100 mM CaCl2 displayed reduced growth relative to the wild type (wt), whereas gcn5 mutant cells spotted on YEA plates supplemented with 1.2 M sorbitol showed no difference in growth from wild-type cells. Plates were incubated for 3 to 5 days at 25°C unless otherwise indicated.
|
![]() View larger version (27K): [in a new window] |
FIG. 3. Characterization of gene expression changes during the response to 1 M KCl. (A) Gene expression changes after 15-min KCl treatment are characterized by robust down-regulation of genes. A scatter plot of DNA microarray signal intensities after 15 min of KCl stress plotted against signal intensities of nonstressed control cells is shown. The flanking diagonal lines indicate a threefold change in expression. (B) Gene expression changes after 60 min of KCl treatment are characterized by robust up- and down-regulation of genes. A scatter plot of DNA microarray signal intensities after 60 min of KCl stress plotted against signal intensities of nonstressed control cells is shown. Annotations are as shown in panel A.
|
|
View this table: [in a new window] |
TABLE 2. Number of genes with altered expression during KCl adaptation
|
|
View this table: [in a new window] |
TABLE 3. Biological processes specifically associated with KCl response genes
|
![]() View larger version (27K): [in a new window] |
FIG. 4. Gcn5 is required for regulation of a subset of genes involved in the KCl stress response. (A) Scatter plot showing an increase in the number of Gcn5-dependent genes and the extent of their dependence during KCl stress. n, number of genes with a mean change in induction of 2-fold (P 0.05). The mean change in induction (n-fold) for all genes in each class is indicated (dark line). (B) There is a significant overlap between genes involved in adaptation to KCl and genes with changed expression in the gcn5 mutant during KCl stress. A P of 0.00317 is the probability with which the observed overlap of 54 genes would be expected by chance. (C) KCl-induced genes show a significant overlap (P = 6.32 x 1012) with genes that are dependent on Gcn5 during KCl adaptation. The overlaps (*) between genes with increased expression in the gcn5 mutant and KCl-repressed genes and between Gcn5-dependent genes and genes that are not involved in KCl adaptation were not considered significant (P > 0.05).
|
In order to get an overall view of Gcn5 dependence during KCl treatment in relation to the transcriptional response to KCl, we compared the change in expression of KCl adaptation genes (60 min) with the effect of gcn5 deletion on expression (KCl, 60 min) for genes that were either induced or repressed by KCl stress. For both sets of genes, the relationship between expression change in response to KCl and the degree of Gcn5 dependence (positive or negative) was essentially random (data not shown). Thus, there is no correlation between the extent of KCl induction and Gcn5 dependence. We conclude that Gcn5 dependence is not correlated to the extent to which genes respond to KCl-mediated stress. Therefore, there must be an alternative explanation that accounts for why some genes require Gcn5 for their regulation while others do not.
Gcn5 dependence is associated with a subset of KCl adaptation processes. To further investigate the reason for the Gcn5 dependence of some genes, we used gene ontology biological process terms to determine whether Gcn5-dependent genes during KCl adaptation are associated with any of the biological functions that change in response to KCl treatment (Table 3). Table 4 shows that genes with reduced expression in the gcn5 strain are significantly enriched and associated with some of the same functions of enriched genes in the set that was induced by KCl after 60 min. Most significantly, these functions include carbohydrate metabolism and nicotinamide metabolism, whose genes were further enriched in the set of Gcn5-dependent genes, even in relation to the enrichment level observed during KCl adaptation. However, other terms that are closely related in the gene ontology, such as vitamin metabolism, generation of precursor metabolites and energy, alcohol metabolism, main pathways of carbohydrate metabolism, and glycolysis, had genes that were also significantly enriched. While several of these terms are also associated with KCl adaptation genes, many classes of adaptation genes are not enriched in the Gcn5-dependent group (e.g., cation transport, cell communication, M phase, meiosis, response to stimulus, and signal transduction). Furthermore, the genes associated with ontology terms were significantly enriched in the gene set with increased expression in the gcn5 mutant (Table 4) and do not overlap with those for the genes that are down-regulated after KCl treatment for 60 min (Table 3). This is consistent with the low significance of the overlap between these sets of genes (Fig. 4). We conclude that a subset of functions associated with KCl-induced genes may have evolved together with Gcn5 dependence, while other adaptation responses may have evolved independently of Gcn5.
|
View this table: [in a new window] |
TABLE 4. Biological processes specifically associated with Gcn5-dependent genes during KCl adaptation
|
![]() View larger version (18K): [in a new window] |
FIG. 5. Identification of DNA sequence motifs that are specifically overrepresented in Gcn5-dependent KCl-induced genes. The bar charts show the frequency of four indicated DNA sequence motifs upstream (800 bp from ATG) of genes within different classes relative to the genomic frequency of the motifs. The classes contain (i) 30 genes that are induced 2-fold after KCl treatment for 60 min, (ii) 52 genes that are induced 2-fold by KCl at 60 min but are not Gcn5 dependent (gcn5 NC; change in induction [n-fold] 1.1-fold), (iii) 27 genes that show reduced expression in the gcn5 mutant after 60 min of KCl treatment ( 2-fold) but are not KCl induced (KCl NC; change in induction [n-fold] 1.3-fold), and (iv) 9 x 30 randomly selected genes. The proportion (%) of genes in each class that contain each motif is indicated above the bars.
|
|
|
|---|
It has been suggested that, in S. cerevisiae, Gcn5 function is characterized by the highly inducible nature of the genes it regulates such that a large proportion of stress-induced genes are associated with Gcn5 function (20). In S. pombe, Gcn5 is required for adaptation to some stress conditions but not others. Furthermore, during the KCl stress response, only a minority of the induced genes are Gcn5 dependent, and we have not observed any correlation between the extent of gene induction and Gcn5 dependence. Consistently, it has been suggested that very small changes in nucleosome positioning in individual promoters can have a large effect on the requirement for histone acetylation during gene activation (24).
Interestingly, we have recently shown that the stress sensitivity profile of gcn5 mutants (KCl and CaCl2) is shared by strains that lack the Tup12 corepressor (10). Tup12, together with its homologue Tup11, are involved in the regulation of KCl-induced genes, such as cta3 (14), and they have recently been shown to function as specificity factors for the stress response in S. pombe (18). In strains lacking both Tup11 and Tup12, the cta3 gene, which is normally specifically derepressed by KCl, is derepressed by a broadened range of stress conditions (18). Tup11 and Tup12 have also been shown to repress chromatin remodeling of the ade6-M26 locus (19), where Gcn5 is required for H3 acetylation and chromatin remodeling during meiosis (32). In addition, the S. cerevisiae GAL1 gene, which is repressed by Tup1-Ssn6, requires the SAGA complex during induction (26). Thus, one possibility that we will investigate in future work is that Gcn5 has a specific role in the derepression of at least some of the genes that are repressed by Tup11 and Tup12 corepressors.
Over half of the genes that are induced in a Gcn5-dependent fashion by KCl after 60 min have been identified previously as CESR genes that are induced by many environmental stresses (6). The overlapping group of genes includes zym1 (encoding metallothionein), mex6 (encoding an mRNA export factor), gut2 (encoding glycerol-3-phosphate dehydrogenase), and pyp2 (encoding a serine/threonine protein phosphatase). Metallothionein is a zinc-binding protein that is thought to transfer zinc to newly synthesized zinc finger transcription factors, and thus, its function, together with those of Mex6 and glycerol-3-phosphate dehydrogenase, can easily be understood in the context of a cell adapting its proteome and metabolism to external stress. Many CESR genes have been shown to respond to the Sty1/Atf1 stress response system (6). Phosphorylation by the Sty1 MAP kinase is thought to activate the Atf1 transcription factor, which subsequently activates many CESR genes. Pyp2 has been shown to inactivate Sty1 in response to stress (25) and may thereby constitute a negative feedback mechanism that regulates the stress response.
As expected from the overlap with CESR genes, we found that the Atf1-binding site (TKACGT) is significantly overrepresented in the upstream region Gcn5-dependent KCl-induced genes. Interestingly, the same site is not significantly overrepresented in KCl-induced genes that are Gcn5 independent or in Gcn5-dependent genes that are not induced by KCl. We found several other sequence motifs that are significantly overrepresented, specifically in the upstream regions of Gcn5-dependent KCl-induced genes. The TCTAGAA motif is contained within the binding site for the S. cerevisiae heat shock factor Hsf1 (33) and is thus a likely binding site for the S. pombe heat shock factor. The occurrence of this motif in CESR genes might be expected, since heat shock is one of the stresses by which they are induced. Indeed, five out of six Gcn5-dependent KCl-induced genes containing the motif have been shown to be induced by heat shock (
2-fold) (6). The other motifs are more obscure. The central hexamer within the ATATATAT motif is a high-affinity binding site for the TBP (15). However, the first and last bases of the motif are not part of the sequence required for TBP binding, and therefore, these sites are not like ordinary TATA boxes. Additionally, ordinary TATA boxes would not be expected to be specifically associated with this set of genes. Since the distance between the TATA box and transcription site is around 30 bp in S. pombe (23), it is possible to calculate the predicted length of the 5' untranscribed regions (UTRs) for genes containing the ATATATAT motif if the motif were to represent a specialized type of TATA box. In this case, the mean 5' UTR length would be about 290 nucleotides, which is significantly longer than the mean length of 170 nucleotides for the 390 5' UTRs that have been mapped (P = 0.056). Taken together, our data suggest that the ATATATAT motifs are not TATA boxes. Proteins that might bind to the GTCATA motif have not to our knowledge been identified.
We have found that a set of related biological-process gene ontology functions are specifically enriched in the group of Gcn5-dependent KCl-induced genes. As shown in Fig. 6, the ontology terms are associated with carbohydrate metabolism, glycolysis, and nicotinamide metabolism. The latter term includes metabolism of nicotinamide adenine dinucleotide, which is an essential cofactor in glycolysis. The only enriched ontology term that is not closely connected to the others is protein folding. Interestingly, the group of related gene ontology terms overlaps with a similar group of genes that are enriched by genes that are induced by KCl treatment for 60 min. As shown in Fig. 6, some terms are enriched by both groups, and two terms are significantly further enriched by the Gcn5-dependent group in relation to their enrichment level by the KCl-induced group. Thus, the group of carbohydrate metabolism-associated terms that is associated with Gcn5-dependent genes is independently associated with KCl-induced genes. Genes induced by KCl after 60 min are also associated with other gene ontology terms, such as cation transport, mitosis, meiosis, and cell communication. Therefore, only a subset of related functions associated with KCl induction are dependent upon Gcn5.
![]() View larger version (36K): [in a new window] |
FIG. 6. Gcn5-dependent KCl-induced genes are involved in similar biological processes. Relationships between GO processes enriched in the groups of KCl-induced genes (gray) and Gcn5-dependent genes (black) (Tables 3 and 4) are shown. Two GO processes are further enriched by Gcn5-dependent, KCl-induced genes (gray and black) (Table 3).
|
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»