Most growth differences between
met31Δ and
met32Δ cells were identified when Met30 function was abrogated, either permanently (in
met30Δ cells) or transiently (upon exposure of cells to cadmium or temperature shift of the temperature-sensitive mutant
met30-6 strain;
Patton et al., 2000 
;
Barbey et al., 2005 
;
Su et al., 2005 
). To identify transcriptional differences in the roles of Met31 and Met32 when Met30 is inactive, we expressed
MET4 in the absence of Met30 (
met30Δ) in strains lacking either
MET31 or
MET32 (
met31Δ or
met32Δ) and analyzed the strains using genome-wide transcriptional microarrays. Genomic alterations allowed
MET4 to be expressed from the galactose-inducible
GAL1 promoter at its endogenous locus (
met4::GAL1-MET4). Short-term expression of Met4 in the absence of Met30 resembles the physiological response of cells exposed to cadmium without eliciting cadmium-induced transcriptional programs that are independent of Met4 (
Fauchon et al., 2002 
;
Lee et al., 2010 
; Supplemental Figure S1). In addition, this experimental design allows investigation of differences between
met31Δ and
met32Δ cells without the need to factor in a multitude of Met30-mediated effects on Met4 and its cofactors (
Kaiser et al., 2000 
;
Rouillon et al., 2000 
;
Kuras et al., 2002 
;
Flick et al., 2004 
,
2006 
;
Chandrasekaran et al., 2006 
;
Menant et al., 2006 
;
Ouni et al., 2010 
).
Although
MET32 deletion bypasses
met30Δ lethality (
Patton et al., 2000 
), it was unknown whether
MET32 deletion would bypass lethality instigated by the combination of Met30 loss and high
MET4 expression. Deletion of
MET32 alone (and deletion of both
MET31 and
MET32) bypassed lethality caused by this condition, but deletion of
MET31 alone was not sufficient for bypass (). Consistent with the difference in bypassing lethality, we identified clear
met31Δ and
met32Δ transcriptional differences upon galactose induction of
met4::GAL1-MET4 met30Δ cells. Scatterplot comparisons of transcriptional profiles from
met31Δ cells and its wild-type counterpart after 90 min of galactose induction showed that Met31 loss has little effect on the Met4 transcriptional signature (). In contrast, approximately one-third of the Met4 core regulon genes, denoted by red diamonds in the scatterplot, failed to induce in the
met32Δ background ().
Of the 274 statistically significant transcripts shared among all profiles (p < 0.05 following Loess and Bayesian normalization), 93 genes had induction Z-scores of 2 or higher after 90 min of galactose induction of
met4::GAL1-MET4 met30Δ cells. Sixty-three of these 93 transcripts failed to reach induction Z-scores of 2 in
met4::GAL1-MET4 met30Δ met32Δ cells. Induction of 37 of these transcripts depended on Met32 but not on Met31 (with Z-scores of ≥2 in
met4::GAL1-MET4 met30Δ met31Δ cells). Functional specification (Funspec) analysis (
Robinson et al., 2002 
) of transcripts that solely depend on Met32 were enriched for genes involved in sulfate assimilation, ion transport, and allantoin/allantoate transport (Supplemental Figure S2). Although 27 of the 93 transcripts failed to reach induction Z-scores of 2 in the absence of Met31, only one statistically significant target solely depended on Met31 (and not Met32). This finding is consistent with the strong similarity of the
met31Δ transcriptional profile to that of its wild-type counterpart (). Twenty-nine transcripts remained induced (with a Z-score of ≥2) in both the
met31Δ and
met32Δ backgrounds. If these 29 genes are Met4 targets, Met31 and Met32 would be redundant for mediating their activation. These targets were enriched for most of the sulfur metabolic pathways (cysteine biosynthesis, methionine metabolism, homocysteine biosynthesis, AdoMet-homocysteine cycle, sulfate assimilation; Supplemental Figure S2).
Multiple Em for Motif Elicitation (MEME) analysis (
Bailey and Elkan, 1994 
) of the promoter regions of the different gene groups revealed slightly different DNA consensus sequences for the Met31/Met32 motif (). The Met31/Met32 motif determined from the Met32-only–dependent promoters closely resembled the core regulon consensus Met31/Met32 sequence (grayscale motif;
Lee et al., 2010 
). The redundant (Met31/Met32 independent) motif showed less emphasis on the cytosines that flank the central TGTGG sequence (). Both Met32-only–dependent and redundant gene sets contained CACGTGA Cbf1 motifs (). We previously defined the core regulon of Met4 and categorized it into three classes based on Cbf1/Met28 dependence (
Lee et al., 2010 
). Class 1 consists of genes whose transcription depends on Cbf1 and Met28. Class 2 genes require Cbf1 and Met28 under the conditions of sulfur limitation but not upon
MET4 expression in the absence of Met30. Class 3 genes are induced independent of either Cbf1 or Met28. Investigation of
met31Δ and
met32Δ microarray profiles with respect to the core regulon showed a correlation between Met32-dependent transcripts and Cbf1/Met28 dependence (Pearson chi-squared p value = 0.0018), with the most-Met32-dependent transcripts found in the Cbf1/Met28-dependent class gene groups (, classes 1 and 2) and the least-Met32-dependent transcripts found in the Cbf1/Met28-independent class (, class 3). This correlation may be due to decreased
MET28 levels in
met32Δ cells.
MET28 appears dependent on Met32 but not on Met31 (, class 2). Hierarchical clustering on the
met31Δ and
met32Δ microarray profiles with respect to sulfur metabolism genes identified transcripts that failed to induce in the absence of Met32 (, red). Consistent with their correlation with Cbf1/Met28 dependency, many Met32-dependent transcripts lie on the upper half of the sulfur metabolic pathway (, red). Clustering also identified targets that exhibited weak or delayed induction in
met32Δ cells compared with its wild-type and
met31Δ counterparts (, blue).
To further examine differences in Met32-dependent transcription, we performed Northern analysis on eight Met4 targets—four involved in sulfate assimilation (
MET3,
MET14,
MET22,
MET16) and four involved in the management of organic sulfur compounds (
MET17,
MMP1,
MUP1,
GSH1; ). As predicted, the absolute transcript levels for all examined sulfur metabolism genes were low before galactose induction in all strains. Northern analysis detected increases in the sulfate assimilation transcripts in
met32Δ cells upon active Met4 expression, but these increases were significantly weaker or delayed compared with its wild-type and
met31Δ counterparts. For
MET22, the increase was less than twofold. Increases in other targets may have failed microarray detection in
met32Δ cells due to low transcript levels. Consistent with microarray analyses,
MET17,
MMP1, and
MUP1 exhibited weak or delayed transcript induction.
GSH1, the critical gene for glutathione synthesis, was induced less than twofold in
met32Δ cells by both Northern and microarray analysis, consistent with the cadmium sensitivity of
met32Δ cells (
Barbey et al., 2005 
).
On the basis of the transcriptional differences between
met31Δ and
met32Δ cells upon the induction of active Met4, we tested whether
met31Δ and
met32Δ cells exhibited different abilities in using sulfur compounds. We conducted yeast spot assays on
met31Δ,
met32Δ,
met4Δ,
met31Δ met32Δ,
cbf1Δ, and
met28Δ cells (lacking
met4::GAL1-MET4 met30Δ alterations) on minimal media plates containing different levels of sulfate, sulfite, homocysteine, cysteine, and glutathione as sole sulfur sources (). Because all strains use methionine as a sole sulfur source, minimal media plates containing 0.1 mM methionine were used as controls (, top). Both
met4Δ and
met31Δ met32Δ cells failed to use sulfate, sulfite, cysteine, and glutathione as sole sulfur sources and were greatly impaired for growth on homocysteine-supplemented media, consistent with Met4 being the sole transcriptional activator of the sulfur metabolic pathway and our previous finding that Met4-activated transcription is lost in the absence of both Met31 and Met32 (
Thomas and Surdin-Kerjan, 1997 
;
Lee et al., 2010 
). In addition, both
cbf1Δ and
met28Δ cells failed to use sulfate and sulfite as sole sulfur sources. These findings are consistent with the key role of Cbf1 and Met28 in sulfate assimilation and previously reported reduced enzyme activities for sulfate assimilation in
cbf1Δ cells (
Thomas et al., 1992 
;
Lee et al., 2010 
). We also noticed that
cbf1Δ and
met28Δ cells grew less well at lower concentrations of homocysteine and glutathione compared with the highest concentration of 0.5 mM. Met31 and Met32 are redundant for use of sulfate, sulfite, homocysteine, cysteine, and glutathione as sole sulfur sources despite the clear transcriptional differences between
met31Δ and
met32Δ cells upon
MET4 expression in the absence of Met30.
In addition to common sulfur compounds, wild-type yeast can use alternative sulfur sources such as sulfonates (
Uria-Nickelsen et al., 1993 
). We previously identified three genomic clusters of Met4-activated genes that are enriched for the metabolism of alternative sulfur sources (
Choi et al., 1998 
;
Hogan et al., 1999 
;
Rouillon et al., 1999 
;
Hall et al., 2005 
;
Lee et al., 2010 
). Hierarchical clustering on these gene groups revealed that genes involved in the transport and metabolism of sulfonates, such as
YIL166C,
JLP1, and
YLL058W, were most affected by Met32 loss (;
Choi et al., 1998 
;
Hogan et al., 1999 
). Therefore, we examined the ability of
met31Δ and
met32Δ cells to use taurine, isethionate, and cysteate as sole sulfur sources and compared their growth with that of
met4Δ,
met31Δ met32Δ,
cbf1Δ, and
met28Δ cells (). Both
met4Δ and
met31Δmet32Δ cells failed to use these compounds as sole sulfur sources. In addition,
cbf1Δ and
met28Δ cells also failed to use sulfonates. Despite their transcriptional differences, both
met31Δ and
met32Δ cells used all different concentrations of taurine, isethionate, and cysteate as sole sulfur sources ().
Given that met32Δ cells were greatly impaired for MET28 induction () and met28Δ cells failed to use sulfate, sulfite, and sulfonates ( and ), the ability of met32Δ cells to use these compounds was unexpected. MET28 may be expressed at such low levels in met32Δ cells upon active Met4 expression that it failed microarray detection, as was the case for many Met32-dependent transcripts ( compared with ). To investigate this possibility, we examined MET28 transcript levels upon galactose-induced expression of Met4 in the absence of Met30 by Northern analysis. Similar to microarray analyses, MET28 was detected in met31Δ cells but not in met32Δ cells (). Although MET28 was not detected in met32Δ cells within 90 min of Met4 expression, MET28 could have been expressed independent of Met32 under the conditions of the spot assay, which was performed on minimal media in cells that contain Met30. To investigate this possibility, we examined MET28 transcript levels in cells with Met30 upon the Met4-activating condition of sulfur limitation in minimal media. Although MET28 induction was impaired in met31Δ met32Δ cells upon sulfur limitation compared with wild-type cells, MET28 transcripts were detected (). These detectable levels of MET28 may reflect different forms of MET28 regulation between galactose induction of met4::GAL-MET4 met30Δ met32Δ cells and sulfur limitation of met31Δ met32Δ cells. This difference could affect recruitment of Met4 and allow strains in the yeast spot assay to express key genes to allow the use of different sulfur sources.
To compare Met31 and Met32 promoter binding, we used genome-altered strains that contained Myc epitopes fused to the carboxyl termini of Met31 and Met32. Met31 and Met32 exhibited slightly different binding patterns. Consistent with our Northern data on transcript dependence, Met32 bound to promoters of the sulfate assimilation genes
MET14 and
MET16 and the glutathione synthesis gene
GSH1, whereas Met31 did not, supporting the more important role of Met32 in the induction of Cbf1-dependent targets (). We previously compared Met31- and Met32-binding patterns to these promoters upon sulfur limitation in minimal media. On sulfur limitation, Met31 and Met32 exhibit similar relative profiles of promoter binding, with Met31 exhibiting strong binding to
MET14,
MET16, and
GSH1 relative to the other promoters analyzed (
Lee et al., 2010 
). Differences in Met31 promoter-binding patterns between these two tested conditions may explain the discrepancy between the lack of sulfate assimilation gene induction upon galactose induction of
met4::GAL-MET4 met30Δ met32Δ cells and the ability of
met32Δ cells to use sulfate and sulfite as sole sulfur sources. Although the relative promoter-binding patterns for Met31 and Met32 were similar upon sulfur limitation, Met32
Myc had ~10-fold higher percentage capture values than Met31
Myc upon sulfur limitation (
Lee et al., 2010 
). We see a similar difference in percentage capture values between Met31 and Met32 upon Met4 expression in the absence of Met30 (, different scales for Met31 and Met32). As a control, Met4 and Rpb3
HA exhibited very similar percentage captures in the two strains, which differ only in the Myc epitope–tagged target (Supplemental Figure S3).
Despite many parameters that differ between these two Met4-activating conditions (galactose induction of met4::GAL-MET4 met30Δ cells in rich media or sulfur limitation in minimal media), we detected consistent differences between Met31 and Met32 expression ( and ). Transcriptional microarrays on met4::GAL-MET4 met30Δ cells showed that MET32 levels increase with active Met4 expression but showed no change in MET31 ( and later in the paper). Northern analysis confirmed this MET32 pattern (), although a high background-to-signal ratio prevented MET31 assessment. The patterns of Met31 and Met32 steady-state protein levels mimicked microarray patterns of MET31 and MET32, respectively (). Met31 was maintained at a constant level, whereas Met32 accumulated with active Met4 expression. These data support differential regulation of MET31 and MET32, with MET32 as a Met4 target gene and MET31 as a nontarget. Likewise, MET32 and MET31 maintained these profiles upon sulfur limitation (). As before, MET32 transcript and Met32 protein accumulated with the active form of Met4 (phosphorylated and deubiquitylated; and ). This induction was lost in met4Δ cells (), consistent with MET32 as a Met4 target gene. MET31 transcripts and Met31 protein remained unchanged upon sulfur limitation in both wild-type and met4Δ cells ( and ), consistent with MET31 not being a target of Met4. Accumulation of Met32 may contribute to its greater impact on the Met4-activated transcription compared with Met31 ( and ). Cbf1 protein levels remained unchanged in both time courses ().
To get a mechanistic understanding of the prominent role of Met32 in Met4-activated transcription when Met30 is absent, we identified and characterized in vivo binding targets of Met4 and Met32 by conducting chromatin immunoprecipitation (ChIP) studies using genomic tiling arrays (ChIP-chip) after 90 min of galactose induction on
met4::GAL-MET4 met30Δ cells. Anti-Met4 antibodies immunoprecipitated 250 genomic targets that exhibited a significant enrichment compared with control samples that lacked antibody (). Microarray data indicated a distinct correlation between Met4 promoter occupancy and transcriptional induction ( and Supplemental Figure S4). As expected, Met4-bound targets are highly enriched for sulfur metabolism genes and included 39 of the 45 core regulon promoters ( and Supplemental Figure S5). The highest Met4 captured targets (, yellow in Met4 ChIP lane) were enriched for both Cbf1 and Met31/Met32 sites as determined by MAST analysis (, black bands in MAST lanes; Pearson chi-squared analyses, p < 0.005;
Bailey and Gribskov, 1998) 
. Consistent with its prominent role in Met4-activated transcription, ChIP-chip studies showed that Met32 bound all identified Met4-bound DNA targets (). Although original ChIP-chip analyses (as outlined in
Materials and Methods) identified three Met4-bound targets not bound by Met32 in a statistically significant manner, closer examination of tiling arrays showed these three targets bound by Met32 (Supplemental Figure S6). Targets of Met4 and Met32 were enriched for sulfur metabolism genes in addition to stress response and various transport genes (Supplemental Figure S7).
MEME analysis of the top 130 Met4- and Met32-bound promoters revealed that Met31/Met32 and Cbf1 motifs closely matched the core regulon consensus sequences ( compared with , gray motifs). Given the large overlap between Met4-bound and Met32-bound targets, sequence differences were minimal. When nonpromoter targets were included in the analysis, MEME occasionally identified a core TGTGGC Met31/Met32-binding sequence in addition to the full Met31/Met32 motif (Supplemental Figure S8). If analysis included 180 bound regions, MEME-CHIP identified Cbf1 motifs that contained a flanking AAT motif in both Met4- and Met32-bound targets that matched part of the RYAAT motif found in Cbf1- and Met28-responsive promoters (Supplemental Figure S8;
Machanick and Bailey, 2011 
;
Siggers et al., 2011 
).
Close examination of tiling arrays showed almost identical binding patterns for Met4 and Met32 at individual targets, further supporting Met32 as the main anchor for Met4 promoter binding (). Most Met4 and Met32 binding occurred at or near promoters (). Sulfur metabolism genes with known Met31/Met32 promoter sites exhibited significant promoter binding by Met32 and Met4 (). Strong Met4 and Met32 promoter binding was also detected at known Met4 targets
ENO1 and
PDC6, which encode glycolytic enzymes of low sulfur content (;
Fauchon et al., 2002 
;
Pereira et al., 2008 
;
Cormier et al., 2010 
). In addition, Met4 and Met32 bound the endogenous
MET4 and
MET30 promoters, consistent with known feedforward and feedback regulation of Met4 (). Of Met4 cofactors, Met4 and Met32 bound genomic regions associated with
MET28,
MET31, and
MET32 but not
CBF1 (). Whereas Met4 and Met32 bound promoters of
MET28 and
MET32, Met4 and Met32 bound the
MET31 open reading frame (ORF; ). Examination of the
MET31 ORF revealed a full Met31/Met32 consensus binding sequence at the region corresponding to the largest ratios of Met4 and Met32 binding; the
MET32 ORF, which is not bound by Met32 or Met4, lacks this full motif (Supplemental Figure S9). In addition to
MET31, 47 other targets were bound by Met4 and/or Met32 at their ORF or 3′ intergenic regions. These targets failed to induce upon induction of active Met4, either suggestive of roles for Met4 and Met32 that are independent of transcriptional activation or an indication of inappropriate binding by Met4 and Met32 to nontargets due to their high protein levels (Supplemental Figure S1C and ). Although Met4 and Met32 binding to these targets may be an artifact, MEME analysis indicated Met31/Met32 motifs in all retrievable ORFs, showing a strong consensus to an abbreviated CTGTGGC motif (Supplemental Figure S10). These atypical targets were enriched for osmosensing and tRNA processing (Supplemental Figure S10).
In addition to Met4-bound targets, ChIP-chip analyses identified 54 targets that registered as Met32 bound but not Met4 bound (). These “Met32-only”–bound targets did not contain matches to either Met31/Met32 or Cbf1 motifs as determined by MAST and were not induced upon active Met4 expression. In fact, several Met32-only–bound targets were ribosomal protein genes that were repressed upon Met4 expression (Supplemental Figure S6 and ). Met32-only targets were also enriched for carbon metabolic pathways (glycolysis, gluconeogenesis, pentose phosphate shunt) and steroid metabolism (Supplemental Figure S11). Closer examination of most Met32-only targets indicated less Met32 binding compared with known Met4 targets and low levels of Met4 at the precise locations bound by Met32 (Supplemental Figure S11). Because most of these Met32-only transcripts were not induced upon Met4 expression, Met32 may perform roles beyond Met4-activated transcription. Alternatively, this binding may be an artifact of high Met32 levels in this growth condition (). If there is a role for Met32 beyond Met4-activated transcription, many Met32-only–bound regions contained multiple genes, confounding identification of relevant targets (Supplemental Figure S11). Whereas MEME analysis failed to identify Met31/Met32 or Cbf1 motifs among these targets, the binding motif for Rap1 was identified, consistent with enrichment for ribosomal protein genes (Supplemental Figure S11). Excess Met32 and Met4 may simply be drawn to these nontargets through interactions with the basal machinery, as Rap1/Gcr1 activation occurs at nuclear pore–associated transcription factories and potentially accounts for >75% of mRNA production in logarithmically growing cells (
Menon et al., 2005 
). In support of this hypothesis, Met4 and Met32 also bound to the highly expressed
GAL promoters upon galactose induction despite the lack of Met31/Met32 sites in these promoters (Supplemental Figure S12). Like Rap1/Gcr1–activated genes, the
GAL genes are transcribed at high levels upon galactose induction and are known to localize to the nuclear periphery, where transcription factories are found (
Berger et al., 2008 
).