H3 Lysine 4 Is Acetylated at Active Gene Promoters and Is Regulated by H3 Lysine 4 Methylation

Methylation of histone H3 lysine 4 (H3K4me) is an evolutionarily conserved modification whose role in the regulation of gene expression has been extensively studied. In contrast, the function of H3K4 acetylation (H3K4ac) has received little attention because of a lack of tools to separate its function from that of H3K4me. Here we show that, in addition to being methylated, H3K4 is also acetylated in budding yeast. Genetic studies reveal that the histone acetyltransferases (HATs) Gcn5 and Rtt109 contribute to H3K4 acetylation in vivo. Whilst removal of H3K4ac from euchromatin mainly requires the histone deacetylase (HDAC) Hst1, Sir2 is needed for H3K4 deacetylation in heterochomatin. Using genome-wide chromatin immunoprecipitation (ChIP), we show that H3K4ac is enriched at promoters of actively transcribed genes and located just upstream of H3K4 tri-methylation (H3K4me3), a pattern that has been conserved in human cells. We find that the Set1-containing complex (COMPASS), which promotes H3K4me2 and -me3, also serves to limit the abundance of H3K4ac at gene promoters. In addition, we identify a group of genes that have high levels of H3K4ac in their promoters and are inadequately expressed in H3-K4R, but not in set1Δ mutant strains, suggesting that H3K4ac plays a positive role in transcription. Our results reveal a novel regulatory feature of promoter-proximal chromatin, involving mutually exclusive histone modifications of the same histone residue (H3K4ac and H3K4me).


Introduction
Histones are covalently modified on many lysine residues. The outcome of these modifications is either to modify chromatin structure directly or provide docking sites for the binding of nonhistone proteins [1]. The same chemical modification (e.g. methylation) on different residues can lead to distinct outcomes. Moreover, some histone modifications function in a combinatorial fashion to generate different functional outcomes [2][3][4]. This led to the notion that histone modifications may represent an epigenetic code that influences gene expression and serves as a ''memory'' of cell identity during development of cell lineages [5][6]. The recent development of high resolution mass spectrometry has enabled the identification of a great number of new histone modifications [7][8][9]. Elucidation of the functions of these new modifications is greatly facilitated in model organisms where, in contrast to vertebrate cells, histone gene mutations that abolish specific modifications can be readily introduced.
A close link exists between H3K4me2 and -me3 and acetylation of the H3 tail. H3K4me3 is often found together with acetylation of other residues (i.e. K9, K14, K18, K23 and K27) on the same H3 molecules [45][46][47][48][49]. In a subset of yeast genes, H3K4me3 directly binds to the PHD finger domain of Yng1, a subunit of the NuA3 histone acetyltransferase (HAT) complex that modifies H3K14, which couples the acetylation and methylation of H3 on different residues [33]. In contrast, through the recruitment of the SET3 complex, H3K4me2 in coding regions promotes deacetylation of H3 in the wake of RNA polymerase II (RNApol II) [50]. These results suggest that there is a highly dynamic and coordinated interplay between histone H3K4 methylation and the enzymes that control H3 acetylation during transcription.
Despite extensive studies of histone H3K4 methylation, the functional implications of other modifications that occur on the same residue have not been investigated. Here, we identified H3K4 acetylation (H3K4ac) in S. cerevisiae using mass spectrometry and a highly specific antibody that we developed. We found that GCN5 and, to a lesser extent, RTT109 are needed for both H3K4ac and H3K9ac in vivo. H3K4 deacetylation in euchromatin is mainly dependent upon HST1, and partly on SIR2. In contrast, removal of H3K4ac from heterochromatin requires SIR2 only. Genome-wide ChIP experiments revealed that H3K4ac is generally found upstream of H3K4me3 in active gene promoters, a pattern which has been conserved at many human CD4 + T-cell promoters [51]. We further demonstrate that H3K4me2 and -me3 mediated by the COMPASS complex limits global levels of H3K4ac at promoters and prevents it from spreading into the 59-ends of coding regions. Using a genetic approach to separate the functions of H3K4ac and H3K4me, we identified a subset of S. cerevisiae genes whose expression depends upon H3K4ac, but not H3K4me. Altogether, our results strongly support a positive role for H3K4ac in gene transcription and identify a novel interplay between two modifications of the same histone residue.

Histone H3 lysine 4 is acetylated in Saccharomyces cerevisiae
One of the most studied and conserved histone modifications is H3K4 methylation (H3K4me), which is tightly linked to transcription in many species [10]. Because of this, and the fact that methylation and acetylation of the same lysine residue are mutually exclusive, we investigated whether H3K4ac was present in S. cerevisiae, a powerful model organism for genetic analysis. For this purpose, we raised an antibody against H3K4ac and confirmed its specificity with peptide binding assays and immunoblotting competition assays ( Figure 1A, 1B). Previous studies have detected the presence of H3K4ac in human cells and an antibody was made commercially available, however the specificity of this antibody was questioned as it showed crossreactivity against the acetylated H4 N-terminal tail and H3K9ac peptides [51][52]. In contrast, our antibody does not show any strong cross-reactivity against H3K9ac or tetra-acetylated H4 peptides ( Figure 1C).
We readily detected a single H3K4ac band in immunoblots of extracts prepared from wild-type (WT) S. cerevisiae cells ( Figure 1D). Importantly, there was no H3K4ac signal in extracts from cells that express either H3-K4R or H3-K4Q as the only source of H3 ( Figure 1D). We also confirmed the presence of H3K4ac by mass spectrometry (MS). Histones were isolated from exponentially growing yeast cells ( Figure 1E). After trypsin digestion, histone peptides were analysed by nano-LC MS/MS. In histones purified from wild-type cells, we identified a doubly charged precursor ion with the expected mass-to-charge ratio (m/z) for the tryptic peptide 3-TKacQTAR-8. The m/z ratio for the monoisotopic form of this precursor ion was determined with high accuracy as 373.7113 ( Figure 1F). The experimental mass of the non-protonated peptide corresponding to this precursor ion (745.4070 Da) was in excellent agreement (21.6 ppm) with the theoretical mass expected for a K4-acetylated peptide (745.4082 Da) and clearly different from the predicted mass of a peptide containing trimethylated K4 (745.4592 Da). Fragmentation of this precursor peptide by collision-activated dissociation resulted in band y-type ion series ( Figure 1G) with diagnostic fragment ions (e.g. the mass difference between the y3 and y5 fragments and the mass of the b2 fragment) that unambiguously proved the presence of an acetyl group on H3K4. However, technical limitations precluded us from determining the abundance of H3K4ac by MS. Using synthetic peptides, we found that the tryptic peptide containing H3K4ac was extremely poorly retained during reverse phase HPLC. Because of this we were not able to accurately assess the stoichiometry of H3K4 acetylation.
To investigate the presence of H3K4ac in chromatin, we performed chromatin immunoprecipitation (ChIP) in strains expressing either H3-WT or H3-K4R. Our initial screen of various genomic regions by qPCR revealed a strong enrichment of H3K4ac (relative to ChIP of non-modified H3) at some gene promoters (MAE1, PHO85, RPL25) a modest enrichment at other genes (BNA5, MCM1, BAP2, BAP3) and a very low enrichment at a region ,1 Kb from telomere III-R (TEL03R) ( Figure 1H). These results are consistent with our genome-wide ChIP-chip data showing an enrichment of H3K4ac in euchromatin at the promoters of highly transcribed genes (see below). We performed control experiments to confirm the target lysine specificity of our antibody in ChIP assays. The ChIP signal obtained with the H3K4ac antibody dropped significantly at all regions tested in the H3-K4R strain ( Figure 1H). As an average for all the promoters tested, the signal obtained by ChIP with our H3K4ac antibody in the WT strain is about 7-fold higher than the non-specific signal (noise) derived from the H3-K4R strain. The average signal-to-noise ratio increases to 9-fold when considering only the four promoters where the signal is most abundant in WT cells. Certain promoters show a higher background, such as the RPL25 promoter (3-fold signal-tonoise), which may reflect some degree of non-specific binding of

Author Summary
In the nucleus of mammals and yeast, DNA is packaged by forming complexes with histone proteins in a structure called the nucleosome, the basic building block of chromatin. The tails of the histones protrude from the nucleosome and can be marked on many amino acid residues by chemical modifications such as methylation and acetylation. A highly studied modification, which is robustly associated with active gene promoters, is histone H3 lysine 4 methylation. We describe here a novel modification, histone H3 lysine 4 acetylation (H3K4ac), which can occur on the same lysine of the histone H3 tail (but not at the same time as methylation). We have identified the enzymes responsible for depositing and removing this mark and mapped its distribution throughout the yeast genome. We found that H3K4ac is present on active genes and is important for the full expression of a subset of them. Strikingly, H3K4 methylation was found in the same promoters as H3K4ac and contributes to regulate the abundance and localisation of H3K4ac. This example of cross-talk between two different modifications of the same residue has fundamental implications for understanding how genes are activated and how their packaging in the nucleus controls this process. Figure 1. H3K4 Is Acetylated in S. cerevisiae. (A) Synthetic peptides derived from the histone H3 N-terminal tail (residues 1 to 10) with modifications at different residues were applied to a nitrocellulose membrane. The membrane was incubated with the H3K4ac antibody, a horseradish peroxidase-conjugated anti-rabbit IgG, and signals visualised by chemiluminescence (left). Ponceau S staining showing that all the peptides bound to nitrocellulose (right). (B) Immunoblots carried out with nuclear extracts from mouse thymocytes. The membrane was cut into strips and incubated with H3K4ac antibody that was pre-incubated with the indicated competitor peptides (1 mg/ml), resulting in a molar ratio of our H3K4ac to other acetylated lysine residues present in proteins in the vicinity of the RPL25 promoter. Nevertheless, our H3K4ac-specific enrichment values are comparable with previously published ChIP data using acetyl-lysine specific antibodies, where the signal-to-noise ratio ranged from 3 to 10-fold [53]. These results validate the specificity of our antibody in ChIP assays and confirm the presence of H3K4ac in gene promoters.

H3K4 acetylation depends upon GCN5 and RTT109
To screen for potential HATs that act on H3K4, we monitored H3K4ac by immunoblotting extracts derived from strains carrying gene disruptions of previously identified HATs. Using this approach, we found that GCN5 and RTT109 were required for H3K4ac in vivo (Figure 2A, lanes 4 and 7). Consistent with a recent report [54], gcn5D and rtt109D single mutants also exhibited lower levels of H3K9ac than WT cells (Figure 2A). A gcn5D rtt109D double mutant showed levels of H3K4ac comparable to the background observed in histones from H3-K4R cells ( Figure 2B, 2C), suggesting that, as is the case for H3K9ac, GCN5 and RTT109 are responsible for essentially all H3K4ac in yeast.

H3K4 deacetylation is dependent upon HST1 and SIR2
To identify potential HDACs responsible for H3K4 deacetylation we repeated our immunoblotting screen with extracts from different HDAC deletion mutants. This revealed that H3K4 deacetylation is, at least in part, dependent upon HST1 and HDA1 ( Figure 2D, lanes 2 and 7). Deletion of HDA1 also led to an increase in H3K9ac, whereas the hst1D mutant showed an increase in H3K4ac but not H3K9ac ( Figure 2D, bottom panel). Because of this, we focussed our experiments on Hst1. S. cerevisiae Hst1 is a member of a family of five NAD + -dependent HDACs, known as sirtuins because of their homology to Sir2. These enzymes are readily inhibited by nicotinamide, one of the products of the deacetylation reaction. Consistent with this, exposure of WT cells to nicotinamide caused an increase in H3K4ac ( Figure S1). Treatment of hst1D cells with nicotinamide led to a further increase in H3K4ac ( Figure S1), which suggested that other sirtuins also contribute to H3K4 deacetylation. We analysed H3K4ac in a collection of sirtuin double mutants and found that deletion of SIR2 in cells lacking Hst1 increased H3K4ac above the levels observed in WT cells and hst1D single mutants ( Figure 2E, 2F, lanes 1-3). This suggests that SIR2 can compensate for the loss of HST1 for global deacetylation of H3K4. In S. cerevisiae, Sir2 promotes heterochromatin formation through histone deacetylation in specific chromosomal regions, namely the silent mating type loci (HMRa and HMLa ), sub-telomeric regions and a subset of ribosomal DNA (rDNA) repeats [55]. By immunoblotting whole-cell extracts, we found no global increase in H3K4ac in a sir2D single mutant ( Figure 2D, lanes 1 and 11). We hypothesised that, in WT cells, Sir2 might deacetylate H3K4 mainly at heterochromatic loci, which represent a relatively small fraction of the genome. To test this hypothesis, we performed ChIP in WT, hst1D, sir2D and hst1D sir2D double mutant cells at different euchromatic and heterochromatic regions known to be deacetylated by Hst1 and Sir2 respectively: the promoters of BNA5 and BTN2, targeted by Hst1 [56], a sub-telomeric region (TEL03R) and the silent mating type locus HMLa, targeted by Sir2 [55]. We also analysed an origin of replication (ARS1233) and the BSC1 gene promoter, which are not known to be targeted by either Sir2 or Hst1. The sir2D single mutation elevated H3K4ac at the TEL03R and HMLa heterochromatic loci to levels comparable to those observed in the hst1D sir2D double mutant ( Figure 2G). The deletion of HST1 had virtually no effect at these regions. We therefore conclude that H3K4 deacetylation in heterochromatic regions is solely dependent upon SIR2. Conversely, deletion of HST1 increased H3K4ac at the promoters of BNA5 and BTN2, but not at other euchromatic regions ( Figure 2G). The SIR2 deletion affected H3K4ac only slightly at those loci whereas the hst1D sir2D double mutant showed higher levels of H3K4ac than in any of the single mutants ( Figure 2G). This is consistent with a previous study [57] and suggests that Sir2 can also contribute to deacetylate H3K4 at Hst1-targeted loci.

Genome-wide presence of H3K4ac in euchromatin and active promoters
To localise H3K4ac throughout the genome, DNA recovered from ChIPs was analysed with high-resolution tiling arrays (Affymetrix 1.0R). Nucleosome density was controlled by normalising H3K4ac to the results of a ChIP experiment performed from the same cell lysates using an antibody against the non-modified C-terminal tail of H3 (H3-C).
The data showed that H3K4ac peaks are distributed along all chromosomes and relatively low levels of H3K4ac are observed at most sub-telomeric regions ( Figure S2). This is most apparent on chromosome III, where H3K4ac is markedly depleted in large sub-telomeric regions that encompass HMLa and HMRa ( Figure  S2), where Sir2 was previously shown to promote histone deacetylation. Closer manual inspection of the results using the UCSC genome browser [58] indicated that the strongest peaks of H3K4ac are located at highly transcribed genes, such as ribosomal protein genes. Two examples (RPL37A and RPS31) are highlighted in Figure 3A. We performed a systematic analysis aimed at determining whether H3K4ac was generally enriched at promoters. For this purpose, we aligned H3K4ac data derived from genes with upstream intergenic regions (IGRs) that contain divergent promoters. Based on this criterion, these IGRs are devoid of DNA sequences corresponding to 39-untranslated regions (39-UTRs) and transcriptional termination sequences. This analysis revealed that H3K4ac peaks slightly upstream of the 59-ends of genes ( Figure 3B). Furthermore, the degree of H3K4ac enrichment correlated with transcriptional activity (as measured in [59]). For the most highly transcribed genes (50+ mRNA/hr), H3K4ac peaks upstream and appears to spread into coding regions ( Figure 3B). In less frequently transcribed genes, H3K4ac also peaks just upstream of 59ends, but sharply drops in coding regions ( Figure 3B). We performed a similar analysis for genes that have converging 39-ends. The corresponding IGRs contain approximately 1000:1 (peptide:antibody). (C) Dot blots of synthetic peptides (100ng each) containing H3K4ac, H3K9ac, and tetra-acetylated histone H4, (K5, 8,12,16)ac, were incubated with H3K4ac and H3K9ac antibodies. Antibody binding was detected as described above. (D) Whole-cell lysates from strains expressing either WT or mutant versions of H3 were analysed by immunoblotting for H3K4ac or a non-modified C-terminal peptide of H3. (E) Coomassie stained gel of histones purified from exponentially growing yeast cells. (F) Doubly charged precursor ion with the expected m/z ratio for tryptic peptide 3-TKacQTAR-8. The experimental and theoretical masses of the non-charged peptides are indicated in the inset along with the mass difference between the empirically determined and the predicted mass of the K4-acetylated peptide. (G) MS/MS spectrum of the doubly charged precursor ion with m/z 373.7113. The peptide sequence is shown from its C-terminus to its N-terminus above the spectrum. (H) Chromatin immunoprecipitation (ChIP) assays were analysed by quantitative PCR to determine the abundance of H3K4ac relative to non-modified H3 ChIP at different loci in WT and H3-K4R cells. doi:10.1371/journal.pgen.1001354.g001 39-UTRs and transcriptional termination sequences, but are devoid of promoters. These chromosomal regions revealed only a modest enrichment for H3K4ac in the most highly transcribed gene group ( Figure 3C, coordinates 0.0 and beyond). This effect probably results from the aforementioned spreading of H3K4ac peaks over the whole body of highly transcribed genes (compare (D-E) Whole-cell lysates obtained from strains with deletions of known HDAC encoding genes were analysed by immunoblotting with the indicated antibodies. (F) Quantification of the immunoblots shown in (E), calculated as described above. (G) Chromatin immunoprecipitation (ChIP) was performed by quantitative PCR to determine the abundance of H3K4ac at different loci in WT, hst1D, sir2D or hst1D sir2D strains. The six target regions were divided into three groups: regions targeted by Sir2, gene promoters targeted by Hst1 and other gene regions not targeted by Sir2 or Hst1. The TEL03R primers amplify a portion of ARS319, a unique sequence region approximately 1 Kb from the telomeric repeats on the right arm of chromosome III. Results are expressed as a ratio of ChIP signals obtained from the same extracts with antibodies against H3K4ac and a non-modified H3 C-terminal peptide. doi:10.1371/journal.pgen.1001354.g002 Figure 3B and 3C). From this data, it is clear that H3K4ac peaks at gene promoters and the degree of enrichment is proportional to the rate of transcription.

H3K4ac is upstream of H3K4me3 at gene promoters
We found that H3K4ac was enriched at transcriptionally active promoters ( Figure 3B) and it has been reported that H3K4me3 is also present at this location in active genes [20]. To explore the relationship between H3K4ac and H3K4me, we compared our genome-wide H3K4ac localisation data with the previously published high-resolution data for H3K4me1, -me2 and -me3 [60]. Strong peaks of H3K4ac were generally located upstream of H3K4me3 in gene promoters, followed by the typical pattern of H3K4me2 and H3K4me1 enrichment along coding regions ( Figure 4A and more examples in Figure S3). There is a considerable overlap between H3K4ac with H3K4me3, although their localisation patterns are not identical. We repeated the H3K4me3 ChIP-chip under our conditions (as described for H3K4ac) and our results (data not shown) produced an enrichment pattern nearly identical to the Kirmizis et al. dataset [60]. Systematic analysis showed that H3K4me3 spreads further into coding regions than H3K4ac in moderately and highly transcribed genes (from 4 to 50+ mRNA/hr) ( Figure 4B). The peak intensity for the H3K4ac signals are slightly upstream of H3K4me3 ( Figure 4B), which is consistent with the UCSC track observations ( Figure 4A).
To investigate whether this pattern was evolutionarily conserved, we examined the localisation of H3K4 modifications in human CD4 + lymphocytes from previously published ChIP-seq data [19,51]. Manual inspection of the data with the UCSC genome browser revealed H3K4ac and H3K4me3 overlapping at many transcriptional start sites, but H3K4ac was slightly upstream of H3K4me3 (representative examples are shown in Figure S4). Interestingly, H3K4ac peaks also overlap with RNApol II more closely than H3K4me3 at these genes ( Figure S4), suggesting a function for H3K4ac in transcriptional regulation in human cells. From these data, we conclude that H3K4ac is present at transcribed gene promoters, and is adjacent to a unidirectional gradient of H3K4me3, -me2 and -me1 that proceeds from promoters into coding regions. This pattern of H3K4 modifications has been conserved during evolution.

Genome-wide overlap of H3K4ac and H3K4me3 at promoters
To investigate and compare the proportion of promoters that are enriched in H3K4ac or H3K4me3, we transformed the logratio data into Z-scores (see Materials and Methods). This was performed to normalise datasets with different scales, thus allowing a direct comparison of the H3K4ac and H3K4me3 datasets. We then established a threshold value above which a promoter was judged enriched in H3K4ac or H3K4me3 relative to H3. When the cut-off was set at the 90 th percentile (only promoters with enrichment values scoring in the top 10% of all values), H3K4ac was found to be enriched at 38% of all gene promoters, whereas H3K4me3 was enriched at 46% of promoters ( Figure 4C and 4F). Of the promoters enriched in H3K4ac, 74% were also enriched in H3K4me3. This value increased to 93% when the cut-off was set to the top 80 th percentile ( Figure 4D and 4F) and 96% at the top 70 th percentile ( Figure 4E and 4F), indicating that most promoters enriched in H3K4ac also have a significant enrichment in H3K4me3. However, a smaller proportion of H3K4me3-enriched promoters also contain H3K4ac: 63%, 69% and 79% for the top 90 th , 80 th and 70 th percentile cut-offs respectively ( Figure 4F). From this we conclude that, genome-wide, H3K4me3 is enriched at a greater number of promoters than H3K4ac, and that most promoters which harbour H3K4ac are also enriched in H3K4me3.

Set1-mediated H3K4me2 and me3 limit the accumulation of H3K4ac
Since H3K4ac and H3K4me3 overlap considerably at the promoters of many genes, we sought to determine if there is a competition for the lysine 4 substrate. Based on mass spectrometry, we found that a SET1 deletion caused a global increase in H3K4 acetylation ( Figure 5A). This result was corroborated by immunoblotting data showing that H3K4ac, but not H3K9ac was elevated in set1D cells ( Figure 5B, 5C). The immunoblotting result truly reflected an increase in H3K4ac in set1D cells, rather than cross-reaction to other sites of acetylation, because no signal was observed in set1D H3-K4R mutant cells ( Figure 5D). In addition, an H3K4ac peptide, but not H3K9ac or non-modified K4 peptides, could compete the signal observed with the K4ac antibody in WT and set1D cells ( Figure 5E). Set1 is part of COMPASS, a complex that has multiple subunits which contribute to the different degrees of H3K4 methylation (me1, me2 or me3) [61][62][63]. To determine which subunits of COMPASS were important to limit the accumulation of H3K4ac, we performed immunoblots using strains with single gene deletions of all the COMPASS subunits. Cells lacking SET1, BRE2, SDC1, SWD3 or SWD1 have in common a complete absence of H3K4me2 and me3, but mutations of these genes differentially affect H3K4me1 ( Figure 5F). For instance, bre2D and sdc1D mutants have WT levels of H3K4me1 ( Figure 5F, lanes 1, 3 and 5), whereas set1D, swd3D and swd1D mutants lack all forms of H3K4me (lanes 2, 6 and 8). Nonetheless, relative to WT cells, all these mutants increased H3K4ac to similar degrees ( Figure 5F, compare lane 1 with lanes 2, 3, 5, 6 and 8). An spp1D mutant with nearly normal levels of H3K4me2, but reduced H3K4me3 showed a moderate increase of H3K4ac ( Figure 5F, lanes 1 and 7). In striking contrast, an shg1D mutant that retained normal levels of H3K4me2 and H3K4me3 did not show a significant increase in H3K4ac compared with WT cells (Figure 5F, lanes 1 and 4). These results argue that the increase in H3K4ac observed in several mutants of the COMPASS complex occurs mainly because of the absence of H3K4me2 and me3, rather than a lack of H3K4me1.
Since a global reduction in H3K4me2 and me3 led to an increase in H3K4ac, we wished to investigate if the opposite was also true. Immunoblots of cell lysates prepared from gcn5D and gcn5D rtt109D strains, with strongly reduced H3K4ac, showed no significant increase in any form of H3K4me ( Figure 2B, 2C). Conversely, we observed no apparent decrease in H3K4me1, me2 or me3 in cells with elevated levels of H3K4ac (hst1D or hst1D sir2D, (Figure 2E, 2F). Therefore, although H3K4ac increases in mutant cells lacking H3K4me2 and me3, changes in H3K4ac do not influence global levels of H3K4 methylation. These results suggest that the relative abundance of H3K4ac does not merely reflect a simple competition for the H3K4 substrate.

Set1 regulates H3K4ac at promoters and the 59 end of coding regions
To see if the deletion of SET1 led to a change in the genomewide distribution of H3K4ac, we performed ChIP-chip experiments to compare the localisation of H3K4ac in set1D and WT cells. Manual inspection of the data in the UCSC Genome Browser revealed changes in the distribution and the relative abundance of H3K4ac in set1D cells, at several gene promoters. For example, at the FKH1 and SUR7 genes, H3K4ac is increased at the promoter region and spreads towards the coding region in the set1D strain ( Figure 5G, 5H). Systematic alignment of divergent promoter regions showed that the deletion of SET1 led to a global increase in the relative abundance of H3K4ac at gene promoters ( Figure 5I, blue versus black dots). We also observed a slight shift of H3K4ac towards the 59 ends of the coding region in set1D cells ( Figure 5I). To better assess the changes in localisation of H3K4ac in the set1D strain, the H3K4ac ChIP-chip data from WT cells were subtracted from those of set1D cells. This analysis clearly showed that the strongest difference in H3K4ac occurs at position 0 relative to gene start ( Figure 5I, green dots), which corresponds to the 59 ends of coding regions. The fact that the shoulder of the H3K4ac peak in the set1D strain is skewed over downstream sequences argues that this increase is not merely due to increased peak signal in set1D cells, since this would lead to a symmetric increase in each shoulder. Instead, this shows that, on average, H3K4ac spreads further into the 59 ends of coding regions in set1D cells. We also observed a global decrease in H3K4ac at the 39 ends of coding regions in set1D mutants ( Figure 5I). The significance of this decrease is unknown. Nevertheless, we can conclude that the global increase in H3K4ac observed in set1D mutant cells by MS and immunoblotting occurs more prominently at the promoters and 59 ends of coding regions relative to other regions.
The global average increase of H3K4ac in set1D cells likely occurs at a subset of genes and not all genes. Manual inspection of the datasets clearly shows that the H3K4ac pattern remains unaffected at certain genes, for example GAL80 and YML053c ( Figure 5H). To provide further evidence that the increase in H3K4ac observed in the set1D strain did not occur equally at all genes, we performed ChIP-qPCR experiments in WT and set1D strains at three genes: CPA2, ARG7 and FKH1. These genes were selected because they respectively showed little, moderate and The H3K4ac / H3 (black) ChIP-chip data were aligned with the H3K4me3 (burgundy), -me2 (purple) and -me1 (blue) data from [60] and displayed in the UCSC genome browser (http://genome.ucsc.edu/). Each vertical line represents one probe in the dataset and the intensity of the color is proportional to the enrichment of the modification. For clarity, only the probes with a log 2 ratio above 0 are shown for each dataset. (B) Systematic alignement of H3K4ac / H3 (red shades) and H3K4me3 / H3 (blue shades) ChIP-chip data at intergenic regions (IGRs) with divergent promoters (as in Figure 3B). The raw ChIP-chip data were converted into Z-scores in order to be compared on the same scale (see Material and Methods). (C-E) Venn diagrams illustrating the overlapping number of promoters that are enriched for H3K4ac or H3K4me3. The cut-off indicates the H3K4(mod) / H3-C Z-score threshold at which a given promoter was judged to be enriched in either H3K4ac or H3K4me3. For example, with a cut-off at the 90 th percentile (C), the red circle contains the number of promoters that have H3K4ac / H3 signal ratios in the top 10% of the genome. (F) A table displaying the data from the Venn diagrams as percentages. doi:10.1371/journal.pgen.1001354.g004 large changes in H3K4ac in set1D strains according to our initial observations in the UCSC Genome Browser. Our ChIP-qPCR analyses confirmed the ChIP-chip results for H3K4ac in the set1D strain ( Figure S5A-S5C). Therefore, although the absence of H3K4me2 and me3 leads to a global increase in H3K4ac ( Figure 5A, 5B), ChIP assays demonstrate that the increase in H3K4ac is more prominent at some genes than others.

H3K4ac and gene expression
Since H3K4ac is enriched at highly transcribed gene promoters, we explored the possibility that H3K4ac might play a role in gene transcription. Genetic analysis of the function of H3K4ac is complicated by the fact that mutation of H3K4 abolishes both acetylation and methylation. To circumvent this problem, we compared the global mRNA profile of a set1D strain, which lacks In (E), a peptide competition assay was performed as described in Figure 1B. (G-H) Examples of the localisation pattern of H3K4ac in either WT (purple) or set1D (green) strains at two genomic regions displayed using the UCSC genome browser (genome.ucsc.edu) [58]. The data represents the normalised (Z-score) log 2 ratio of H3K4ac over H3-C (see Materials and Methods). (I) Systematic analyses of H3K4ac in WT versus set1D cells in all intergenic regions (IGRs) that contain divergent promoters. The raw H3K4ac ChIP-chip data in WT or set1D cells were normalised either against their respective H3-C ChIP (black: WT, blue: set1D), converted into Z-scores (see Material and Methods) and aligned relative to the 59-end of the corresponding ORF (normalised as described in Figure 3). The difference between the set1D and WT datasets (set1D -WT, green) is also shown. doi:10.1371/journal.pgen.1001354.g005 H3K4me but retains H3K4ac, with that of an H3-K4R mutant where both H3K4 acetylation and methylation are lost. In order to get a global view of the expression profiles in both mutants, we ranked the genes by log 2 fold change in mutant strains relative to WT cells ( Figure 6A). At one end of the spectrum, we found many genes that are up-regulated in both the H3-K4R and the set1D mutant strains ( Figure 6A, left part). Many of these genes are upregulated only in H3-K4R mutant cells ( Figure 6A, 6B). One possibility is that expression of some of these genes may be affected simply because of the lysine-to-arginine mutations, rather than a lack of modification. However, there is considerable overlap between the genes up-regulated in H3-K4R and set1D mutant strains ( Figure 6B). We conclude that the repressed state of these genes in WT cells is likely to be dependent on H3K4me, rather than H3K4ac. This is consistent with previous reports showing a repressive function for SET1 [64][65][66].
At the other end of the spectrum, we found a group of genes which were down-regulated in the H3-K4R strain, but were unaffected in the set1D strain ( Figure 6A, right part), suggesting a positive role for H3K4ac in the expression of these genes. We identified 37 genes that were down-regulated 2-fold or more in the H3-K4R strain and 44 in the set1D strain. Strikingly, there was no overlap between these two groups of genes ( Figure 6C), suggesting a function of H3K4 in gene expression that is independent of methylation. We next compared the H3K4ac abundance at promoter regions in WT (from our ChIP-chip dataset) with the expression change of the gene in the H3-K4R strain. The genes that are down-regulated in the H3-K4R strain have high levels of H3K4ac at their promoters in WT cells ( Figure 6D), which is consistent with a direct role of H3K4ac in gene expression. Gene ontology term enrichment analysis of the genes down-regulated 2fold or more in the H3-K4R strain, but whose expression was not affected in a set1D strain, identified genes encoding amino acid transporters (p = 1.42610 25 ): BAP2/YBR068C, TAT1/YBR069C, BAP3/YDR046C, GNP1/YDR508C, TAT2/YOL020W. Inspection of our ChIP-chip data for these genes confirmed that they all had high levels of H3K4ac at their promoters in WT cells ( Figure 6E-6H). We confirmed by RT-qPCR that these genes are down-regulated in an H3-K4R strain, but unaffected in a set1D strain ( Figure 6I-6M). In contrast to the other RNAs that we monitored, the expression of the TAT2 RNA was unexpectedly restored in a set1D H3-K4R double mutant strain, suggesting an indirect effect of the SET1 deletion on this gene ( Figure 6M). Taken together, our data suggest that H3K4ac has a positive and direct role in the transcription of a subset of yeast genes.

Discussion
Histone H3K4 methylation is an extensively studied modification in a wide range of organisms, and particularly in the budding yeast Saccharomyces cerevisiae. Using MS and highly specific antibodies we demonstrated that H3K4 is also acetylated in budding yeast. An earlier study of histone modifications had shown that both H3K4ac and H3K4me exist in human, mouse and ciliates [52] and, very recently, in Schizosaccharomyces pombe [67]. The versatility of genetic tools available in S. cerevisiae, and the ability to mutate histone H3K4, has enabled us to uncover a positive function for H3K4ac in the expression of a subset of genes that carry H3K4ac in their promoters. The presence of H3K4ac at active gene promoters has been conserved in human CD4 + cells [51], which suggests a similar role for H3K4ac in gene expression in mammals. We also discovered that the COMPASS complex and H3K4me2 and me3 appear to globally limit the abundance of H3K4ac at promoters and its spreading into the 59-coding regions of many genes.

Acetylation of the H3 N-terminal tail by specific HATs
Our targeted screen of known HAT-encoding genes has allowed us to identify GCN5 and RTT109 as responsible for H3K4ac in vivo. Gcn5 is the catalytic subunit of the ADA, SAGA, SLIK/ SALSA transcriptional co-activator complexes [68] and has been highly conserved during evolution [69]. Our results are consistent with the fact that Gcn5 has been shown to acetylate multiple lysine residues within the H3 N-terminal tail [70][71] and is localised to active gene promoters dispersed throughout the genome [56]. Interestingly, RTT109 is also important for the acetylation of H3K4. Furthermore, the gcn5D rtt109D double mutant has markedly less H3K4ac than the single mutants. Collectively, our results indicate that GCN5 and RTT109 promote H3K4ac via different pathways and/or in different contexts. If they were fully redundant, there should be no decrease in H3K4ac in the single mutants, which is not the case. However, we cannot exclude the possibility that they may also have partially redundant functions in mediating H3K4ac. RTT109 encodes a HAT that mediates H3K56ac [72][73][74][75], which is involved in nucleosome assembly during DNA replication and the DNA damage response [76][77][78]. In addition, Rtt109 has recently been implicated in acetylating H3K9 and H3K27 in vivo and in vitro [54,[79][80]. As is the case for H3K4ac, H3K9ac and H3K27ac are virtually undetectable by immunoblotting and mass spectrometry in gcn5D rtt109D cells [79]. Therefore, the same HATs are responsible for the bulk of both H3K4 and H3K9 acetylation.

H3 N-terminal tail acetylation and gene transcription
The genome-wide localisation pattern of H3K4ac showed strong peaks at the promoters of active genes. This enrichment has also been observed in the genome-wide patterns of other acetylated lysines in the H3 N-terminal tail, such as H3K9ac, H3K14ac and H3K18ac [20][21], suggesting that these modifications often occur together on chromatin. Our microarray analysis identified a relatively small group of genes (37 genes downregulated by 2-fold or more, Figure 6C) that require H3K4ac, but not H3K4me, for normal expression. However, based on our ChIP-chip data, many other promoters contain high levels of H3K4ac. It is tempting to speculate that acetylation of H3K4 functions in a partially redundant manner with other sites of acetylation at many promoters to increase accessibility to the transcriptional initiation machinery and thereby facilitate gene expression. Interestingly, our gene ontology term analyses revealed a group of amino acid transporter genes (BAP2, BAP3, TAT1, TAT2, GNP1) that are co-regulated by levels of extracellular amino acids [81][82][83][84]. An interesting possibility is that the transcription factors involved in their expression require H3K4ac either alone versus WT cells (derived from the mRNA expression microarray data). A sliding window average was applied to the data on both axes (window = 50, step = 1). (E-H) Localisation of H3K4ac near genes that belong to the ontology group of ''amino acid transporters'' and are poorly expressed in H3-K4R mutants compared with WT cells: (BAP2/YBR068C, TAT1/YBR069C, BAP3/YDR046C, GNP1/YDR508C, TAT2/YOL020W). The start sites of these genes are indicated by arrows and vertical dashed lines indicate the peaks of H3K4ac in their upstream regions. The display is from the UCSC genome browser (http://genome.ucsc.edu/). (I-M) RT-qPCR analyses of RNA extracted from the indicated strains (WT, H3-K4R, set1D, H3-K4R set1D). The values were normalised to ACT1. The data represent an average of three biological replicates. doi:10.1371/journal.pgen. 1001354.g006 or in combination with other acetylated residues to specifically promote initiation of these genes.

Regional and functional specificity of HST1 and SIR2
We identified a role for HST1 in global deacetylation of H3K4 in euchromatin. HST1 encodes a class III HDAC found in at least two complexes. Hst1 is part of a complex with Sum1 and Rfm1 [85][86]. Interestingly, the deletion of SUM1 did not result in a significant global increase in H3K4ac based on immunoblotting (data not shown). Hst1 has also been found as a sub-stoichiometric subunit of the SET3 complex (SET3C) [87], which represses middlesporulation genes during mitosis [88]. However, neither deletion of HOS2 (encoding the main HDAC in the complex) ( Figure 2C), nor that of any other components of SET3C led to a global increase in H3K4ac (data not shown). Our results suggest that Hst1 might act in a non-targeted manner to promote global deacetylation of H3K4. Alternatively, Hst1 might be part of another, as yet unidentified complex involved in H3K4 deacetylation.
Interestingly, we also observed some redundancy between HST1 and SIR2 for deacetylation of H3K4 in euchromatic regions. In specific conditions, Hst1 and Sir2 have been shown to act outside of the chromosomal regions where they ''normally'' exert their functions. For example, overexpression of HST1, but not HST2 or HST3, can specifically compensate for the loss of silencing at HMRa that occurs in a sir2D mutant [89]. Our results are also consistent with studies showing a strong connection between these two genes. Both genes share 71% identity in their sequence over 84% of the length of HST1 [90]. SIR2 was previously shown to affect the expression of euchromatic genes when HTZ1 and SET1 were deleted [91] and to interfere with the firing of a number of replication origins [92]. Moreover, studies of chimaeric proteins demonstrated that the N-terminal domains of Hst1 and Sir2 control the chromosomal location and function of both proteins. In contrast, other domains of Hst1 and Sir2 could be swapped, supporting the notion that these enzymes have similar substrate specificity [93]. We propose here that H3K4 is an important target for deacetylation by both Hst1 and Sir2. HST1 and SIR2 are both part of the same phylogenetic branch (subgroup Ia) of NAD +dependent HDACs (class III) that also includes human SIRT1 and Drosophila Sir2 [94], which are, therefore likely candidates for deacetylation of H3K4 in these organisms.

Patterns of mutually exclusive histone modifications
The relationship between H3K4ac and H3K4me differs from that of other modifications that occur on the same lysine residue. The best studied case is that of H3K9, which is either acetylated or methylated. H3K9ac and H3K9me are present at different genomic regions [51] and have opposing roles in transcriptional regulation. H3K9ac is found in promoter regions and stimulates transcription through the recruitment of TFIID [95], whereas H3K9me2 and me3 contribute to gene repression and pericentric heterochromatin structure by binding HP1 proteins through their chromo-domains [96][97]. In this case, deacetylation of H3K9 is essential for the establishment of H3K9 methylation [98][99][100]. Another example is H3K36ac and H3K36me3, which are respectively present in promoters and coding regions of transcriptionally active genes [101]. In contrast to these modifications, H3K4ac overlaps considerably with H3K4me3 in promoter regions. Furthermore, H3K4me2/me3 globally limit the levels of H3K4ac, which suggests the existence of mechanisms to establish the correct patterns of localisation of these mutually exclusive modifications. For example, H3K4me2/me3 could prevent H3K4ac from occurring in gene coding regions by simple competition for the H3K4 substrate. Another possibility is that H3K4 methylation might be a prerequisite step in a negative feedback loop that enables Hst1 to fine tune the levels of H3K4 acetylation in promoters and 59-coding regions. Mechanisms that confine H3K4ac to some promoters may contribute to prevent spurious transcriptional initiation events [102] or attenuation of transcriptional initiation [66]. The potential to restrict the spreading of histone H3K4ac suggests a novel function for H3K4 methylation and reveals a previously unrecognised layer of chromatin regulation linked to the regulation of transcription in vivo.

H3K4ac antibody production and validation
Immunisation procedures and animal handling were carried out by Eurogentec. Sera from ten animals were screened prior to immunisation. Rabbits were immunised with the following peptide: ART(acetyl-K)QTARKSC. The efficiency of antibody production was monitored using ELISA. The antibody was affinity-purified using the same peptide coupled to a Sulfolink column (Thermo) and its specificity tested by dot blots with synthetic peptides and experiments where 1 mg/ml of the peptides were incubated with the antibody prior to immunoblotting.

Yeast strains and media
For most experiments, yeast strains were grown at 30uC in YPD (1% yeast extract, 2% peptone, 2% dextrose). A complete list of yeast strains is available in Table S1.

Histone sample preparation and mass spectrometry
Histones were acid extracted from set1D hst1D cells essentially as described previously [104], except that we added 100 mM sodium butyrate and 100 mM nicotinamide in the nuclear isolation buffer and wash buffers. H3 was purified and digested as previously described [105]. Intact core histones (approximately 10 mg of total protein acid extract) were separated using an Agilent 1200 HPLC system equipped with a fraction collector. Separations were performed using an ACE C 8 column (5 mm, 300 Å ), 15064.6 mm i.d., with a solvent system consisting of 0.1% trifluoroacetic acid (TFA) in water (v/v) (A) and 0.1% TFA in acetonitrile (v/v) (B). Gradient elution was performed from 5-70% B in 60 minutes at 0.7 ml/min. Fractions were collected in conical tubes in 60 second time slices. Individual histone peaks that eluted over multiple fractions were pooled together. Histone fractions were evaporated in a Speed-Vac. Dried fractions were resuspended in 0.1 M ammonium bicarbonate (without pH adjustment) and digested overnight at 37uC using 1 mg of trypsin. Samples were acidified with 5% TFA in water (v/v) prior to LC-MS/MS analysis. Tryptic digests of histone extracts were injected onto a Thermo Electron LTQ-Orbitrap XL mass spectrometer equipped with an Eksigent nanoLC separation module. Samples were loaded onto a C 18 trapping column at 10 mL/min using 0.2% formic acid (v/v) in water. Peptides were eluted onto a C 18 analytical column (10 cm6150 mm i.d.) at 0.6 ml/min. Gradient elution was performed using a solvent system of 0.2% formic acid in water (A), and 0.2% formic acid in acetonitrile (B). Peptides were separated by gradient elution from 0 to 70% B in 60 minutes. For mass calibration, we used an internal lock mass [protonated (Si(CH 3 ) 2 O)) 6 ; m/z 445.12057]. In each data collection cycle, one full MS scan (m/z 300-2000) was acquired in the Orbitrap (60000 resolution setting, automatic gain control (AGC) target of 10 6 ), followed by 3 datadependent MS/MS scans in the LTQ (AGC target 10000; threshold 5000) for the 3 most abundant ions and collision induced dissociation (CID) for fragmentation.

Chromatin immunoprecipitation (ChIP)
Extracts containing fragments for chromatin imunoprecipitation (ChIP) were prepared as follows. 200 ml of exponentially growing cells (O.D. 600 of 0.6 to 0.8) in YPD were treated with 1% formaldehyde at room temperature for 10 minutes. Formaldehyde was quenched with 0.125 M of glycine and cells were washed three times in cold water. Pellets (approximately 1.2610 9 cells) were resuspended in 1.2 ml of lysis buffer (50 mM HEPES-KOH pH 7.5, 140 mM NaCl, 1% Triton X-100, 0.1% Na-deoxycholate, 1 mM EDTA, 1 mM PMSF, 30 mM sodium butyrate, complete EDTA-free protease inhibitor cocktail, Roche), and split into three tubes. After adding one volume of glass beads to each sample, cells were disrupted using a bead beater (Disruptor Genie, Scientific Industries) at maximal speed for 2h at 4uC. After removing the glass beads, the lysates were transferred to fresh tubes and sonicated for 15 minutes (30 seconds ON, 30 seconds OFF) at high intensity in a Bioruptor (Diagenode) connected to a water cooler at 4uC. Lysates were clarified by centrifugation and pooled into a 2 ml tube. After addition of 800 ml of fresh lysis buffer and mixing, lysates were clarified again by centrifugation. 4 ml of lysate (1% of the whole cell extract, WCE) was kept aside as an input control. For immunoprecipitation, aliquots of 400 ml were prepared for each ChIP and mixed with the appropriate antibody: either 5 ml of anti-H3K4ac, 5 ml of anti-H3-C (Abcam ab1791), 3 ml of anti-H3K4me3 (Abcam ab8580) for at least 1h at 4uC. Before use, Dynabeads M-280 coupled to sheep anti-rabbit IgG (InVitrogen) (30 ml per ChIP) were washed in lysis buffer containing 5 mg/ml bovine serum albumin. The beads were added to each lysate and incubated at 4uC for at least 1 h. Beads were washed twice with lysis buffer, twice with lysis buffer containing 500 mM NaCl, twice in wash buffer (10 mM Tris-HCl pH 8.0, 250 mM LiCl, 0.5% Nonidet-P40, 0.5% Na-deoxycholate, 1 mM EDTA) and once in TE (10 mM Tris-HCl pH 8.0, 1 mM EDTA). To elute the DNA, 100 ml of a 10% Chelex 100 (Biorad) suspension in water was added to the beads (as well as to the 1% WCE) and incubated at 100uC for 10 minutes. Tubes were then cooled and treated with 0.2 mg/ml ribonuclease A at 50uC for 30 min, followed by 0.2 mg/ ml proteinase K at 50uC for 30 min. Samples were incubated again at 100uC for 10 min to inactivate the proteinase K and then cooled down and spun down to get rid of debris. 80 ml of each supernatant was kept at 220uC for qPCR or ChIP-chip analysis.

ChIP-chip
All ChIP-on-chip datasets have been deposited in NCBI's Gene Expression Omnibus [106] and are accessible through GEO Series accession number GSE27307 (http://www.ncbi.nlm. nih.gov/geo/query/acc.cgi?acc = GSE27307). Samples were processed according to the Affymetrix protocol for amplification, fragmentation and labelling of DNA for hybridisation to a GeneChip S. cerevisiae Tiling 1.0R Array (www.affymetrix.com). All ChIP-chips were performed in triplicate from independent colonies. Pre-processing of the data was carried out with the Affymetrix Tiling Analysis Software with default settings at 200bp bandwidth and expressed as a log 2 ratio of the ChIP signal obtained for modified H3 (K4ac or K4me3) over that obtained from non-modified C-terminal H3 ChIP. The datasets used in Figure 4C and Figure 5F, 5G were normalised using the standard score (Z-score) equation: Z = (x2m)/s, where x is the [log 2 (H3K4(mod) / H3-C)] raw value for each probe; m is the mean raw value of all probes and s is the standard deviation of the raw values for all probes. The resulting datasets have a normal distribution with a mean value of 0 and a standard deviation value of 1, and thus can be compared directly on the same scales. Data was displayed on the UCSC genome browser with the October 2003 S.cerevisiae genome assembly (http://genome.ucsc. edu/). Systematic alignment of ChIP-chip data to gene promoters and terminators were performed using the database management tools provided in the web application Galaxy (http://main.g2.bx. psu.edu/) and graphics were created with R (http://www. r-project.org/).

RNA extraction and expression analysis using microarrays
Total RNA was extracted from three independent colonies using the hot phenol method as described previously [107]. A total of 7 mg of RNA were processed for microarray analysis using the One-Cycle Target Labeling package from Affymetrix and hybridised to a GeneChip Yeast Genome 2.0 Array following the manufacturer's instructions (http://www.affymetrix.com/). Probe intensity data was pre-processed using the mas5 program of the R-Bioconductor software (http://www.bioconductor.org/ index.html) [108]. The data from either H3-K4R or set1D mutants were normalised to their WT counterpart and the genes were sorted according to log 2 fold change. The log 2 fold change values for all genes are presented in Table S2 and represented graphically in Figure S6. For RT-qPCR analyses, 1 mg of total RNA was treated with DNase I, which was then inactivated following the manufacturer's instructions (DNA-free kit, Ambion). RT was performed using the First-Strand cDNA Synthesis Kit with M-MLV RT (Invitrogen) and oligo-dT primers, following the manufacturer's instructions.
Quantitative PCR 1 ml of the ChIP sample was mixed with 50 nM of regionspecific primers and SYBR Green JumpStart Taq ReadyMix (Sigma) in a total volme of 20 ml and analysed on an Opticon real-time PCR machine (MJ Research/Bio-Rad). Relative binding values were extrapolated by subtracting the Ct of nonmodified H3-C to that of modified H3 (K4ac or K4me3), and converting the DCt with the following formula: x = 2 (2DCt) , where x represents the relative binding value of modified H3 on nonmodified H3-C. Figure S1 Nicotinamide Treatment Increases H3K4ac in WT and hst1D cells. Whole-cell lysates from cells treated with 0, 5 or 25 mM nicotinamide for the indicated times were analysed by immunoblotting. Found at: doi:10.1371/journal.pgen.1001354.s001 (0.40 MB TIF) Figure S2 Genome-Wide Location of H3K4ac along the 16 S. cerevisiae Chromosomes. The log 2 of the H3K4ac / H3 derived from ChIP-chip data was aligned to the chromosomes in the UCSC genome browser (http://genome.ucsc.edu). The graphs are drawn to scale for chromosome length, except for the chromosome XII, on which the repeated ribosomal DNA locus is represented only once on the array. The y-axis is scaled according to maximal and minimal values for each individual chromosome. In general, regions that are enriched (peaks) or depleted (troughs) for H3K4ac are interspersed and evenly distributed across all chromosomes. One notable exception is the ends of chromosome III, where large regions that extend from telomeres to beyond the silent mating type loci (HMLa and HMRa; highlighted) are depleted in H3K4ac. Found at: doi:10.1371/journal.pgen.1001354.s002 (2.00 MB TIF) Figure S3 Localisation of H3K4ac, H3K4me3, -me2 and me1 at several genomic regions in yeast. The log 2 of the H3K4ac / H3 derived from ChIP-chip data was aligned to the chromosomes in the UCSC genome browser and aligned with the H3K4me1, me2 and me3 data published previously [59]. Each vertical line represents one probe in the dataset and the intensity of the color is proportional to the enrichment of the modification. For clarity, only the probes with a log 2 ratio above 0 are shown for each dataset. Found at: doi:10.1371/journal.pgen.1001354.s003 (1.05 MB TIF) Figure S4 Localisation of H3K4 Modifications (H3K4ac, H3K4me3, H3K4me2, H3K4me1) and RNA polymerase II (Pol II) in Human CD4 + T Cells. ChIP-seq data from human CD4 + lymphocytes previously published by the Zhao laboratory [19,51] were displayed on the UCSC genome browser (http://genome. ucsc.edu/). Genes are illustrated at the bottom of each panel with an arrow indicating the transcriptional start site and direction. Found at: doi:10.1371/journal.pgen.1001354.s004 (2.50 MB TIF) Figure S5 ChIP-qPCR confirmation of H3K4ac regulation by SET1 at specific gene promoters. We performed ChIP against either H3K4ac (top row), H3K4me2 (middle row) or H3K4me3 (bottom row) in WT and set1D cells. We analysed the enrichment over three genomic regions by qPCR: CPA2 (A), ARG7 (B) and FKH1 (C) which showed low, moderate and high increase in H3K4ac in set1D cells compared to WT in the coding region in our ChIP-chip analyses. Modification-specific ChIPs were normalised with ChIPs against non-modified H3-C using the same whole cell extract and experiments were done in biological triplicates. Found at: doi:10.1371/journal.pgen.1001354.s005 (0.51 MB TIF) Figure S6 Gene Expression Changes in H3-K4R and set1D Mutants Relative to WT Cells. All genes were color-coded and ranked according to log 2 fold change in the H3-K4R strain or the set1D strain versus WT cells. On the right, the top and bottom part of the graph were cropped to reveal the general overlap between the genes that are up-regulated (yellow), but not down-regulated (blue). The right part of the figure contains the list of all genes with a fold-change equal or greater than two in the H3-K4R strain. A complete list of all the genes and their fold-change values is presented in Table S2.