H2B ubiquitylation is part of chromatin architecture that marks exon-intron structure in budding yeast
- Grace S Shieh†1,
- Chin-Hua Pan†2,
- Jia-Hong Wu†1,
- Yun-Ju Sun1,
- Chia-Chang Wang1,
- Wei-Chun Hsiao2,
- Chia-Yeh Lin2,
- Luh Tung3,
- Tien-Hsien Chang3,
- Alastair B Fleming5,
- Cory Hillyer7,
- Yi-Chen Lo4,
- Shelley L Berger6,
- Mary Ann Osley7Email author and
- Cheng-Fu Kao2Email author
© Shieh et al; licensee BioMed Central Ltd. 2011
Received: 13 September 2011
Accepted: 22 December 2011
Published: 22 December 2011
The packaging of DNA into chromatin regulates transcription from initiation through 3' end processing. One aspect of transcription in which chromatin plays a poorly understood role is the co-transcriptional splicing of pre-mRNA.
Here we provide evidence that H2B monoubiquitylation (H2BK123ub1) marks introns in Saccharomyces cerevisiae. A genome-wide map of H2BK123ub1 in this organism reveals that this modification is enriched in coding regions and that its levels peak at the transcribed regions of two characteristic subgroups of genes. First, long genes are more likely to have higher levels of H2BK123ub1, correlating with the postulated role of this modification in preventing cryptic transcription initiation in ORFs. Second, genes that are highly transcribed also have high levels of H2BK123ub1, including the ribosomal protein genes, which comprise the majority of intron-containing genes in yeast. H2BK123ub1 is also a feature of introns in the yeast genome, and the disruption of this modification alters the intragenic distribution of H3 trimethylation on lysine 36 (H3K36me3), which functionally correlates with alternative RNA splicing in humans. In addition, the deletion of genes encoding the U2 snRNP subunits, Lea1 or Msl1, in combination with an htb-K123R mutation, leads to synthetic lethality.
These data suggest that H2BK123ub1 facilitates cross talk between chromatin and pre-mRNA splicing by modulating the distribution of intronic and exonic histone modifications.
Genome-wide histone modification maps have now been generated for a number of eukaryotic organisms. These maps have revealed the preferential localization of specific marks to active or silent chromatin and the association of marks of active transcription with different regions of genes. For example, H3 trimethylation on lysine 4 (H3K4me3) is enriched at the 5' ends of actively transcribed genes while H3 trimethylation on lysine 36 (H3K36me3) is localized towards the 3' ends of coding regions. These localization patterns are related to the roles that the marks play in transcription: H3K4me3 regulates the efficiency of transcription initiation and early steps in transcription elongation, and H3K36me3 prevents the utilization of cryptic initiation sites in coding regions and controls aspects of transcription termination and processing [1–4]. Most eukaryotic genes are modular, containing multiple exons interrupted by introns. Genome-wide histone modification maps from C. elegans and human revealed that intron-exon chromatin is also preferentially marked, with transcriptionally active modifications generally excluded from introns and concentrated in exons [5–11]. These studies concluded that this pattern was primarily the consequence of different levels of nucleosome occupancy in these regions because nucleosomes were depleted in introns relative to exons. However, a recent analysis of published human epigenomic data found that 10 histone modifications were enriched in the 5' introns of human genes independently of the level of nucleosome occupancy . It was suggested that the presence of these marks reflects an aspect of the splicing process such as exon definition and could play a direct role in regulating splicing. Thus, the location of intragenic histone modifications and the functional roles associated with different localization patterns remain areas of intense investigation.
One important intragenic histone modification is the monoubiquitylation of H2B (H2BK123ub1). H2B is ubiquitylated co-transcriptionally and in turn regulates the presence of other active chromatin marks during the transcription process, including H3K4, H3K36, and H3K79 methylation [13–19]. The presence of H2BK123ub1 in chromatin has been associated with both nucleosome stabilization and destabilization. H2BK123ub1 and the histone chaperone, Spt16, have been shown to function interdependently during transcription elongation to regulate nucleosome reassembly and preserve chromatin integrity [20–22]. Biochemical evidence and genomic nucleosome occupancy data also indicate that the presence of H2BK123ub1 generally promotes nucleosome stability [23, 24]. However, it was recently shown that synthetic nucleosome arrays containing H2BK123ub1 are less compact and exhibit an increase in inter-nucleosomal distance, as compared to arrays containing unmodified H2B . In addition, two recent reports described a putative role for this modification in mediating chromatin decondensation at DNA damage sites [26, 27]. Thus, H2BK123ub1 may differentially affect chromatin structure in a context dependent manner.
In this report, we generated a genome-wide map of H2BK123ub1 occupancy in budding yeast to determine if the distribution of this modification could be related to additional biological processes. We found that H2BK123ub1 was enriched across gene coding regions and marked both introns and exons of ribosomal protein (RP) genes, and that the level of this mark was further increased at 3' intron-exon boundaries. The presence of H2BK123ub1 in introns of RP genes was separable from nucleosome occupancy, which was generally lower in introns compared to exons. In addition, we noted that disruption of H2B ubiquitylation tended to alter the distribution of H3K36 trimethylation in intragenic regions. H3K36me3 has been functionally linked to pre-mRNA splicing in worms and humans [6, 7, 28]. Furthermore, when an htb-K123R mutation was combined with deletions of LEA1 and MSL1, whose products facilitate U2 snRNA association with pre-mRNA [29, 30], we found a synthetic lethal phenotype. These data suggest that by modulating the distribution of intronic and exonic histone modifications, H2BK123ub1 facilitates cross talk between chromatin and pre-mRNA splicing.
H2B ubiquitylation is enriched in transcribed regions
A recent ChIP-chip study using a polyclonal antibody directed against yeast H2BK123ub1 also found that this mark was enriched at gene coding regions . We compared the two genome-wide data sets and found almost identical patterns of H2BK123ub1 localization at a number of different genomic locations, with somewhat higher levels of H2BK123ub1 detected by the ChDIP method (Additional file 1, Figures S2A, S3A, S4A) Thus, two different approaches to detect H2BK123ub1 lead to the conclusion that this mark is primarily associated with the transcribed regions of genes.
Next, we investigated the relationship between H2BK123ub1 levels and transcription rate. We found that the level of H2BK123ub1 at coding regions was proportional to transcription rate, with the highest levels occurring at genes that have a transcription rate of > 16 copies of mRNA per hour (Figure 1C; Additional file 1, Figure S1). Genes in this category represent 8.2% (389 genes) of the 4802 genes whose transcription is detectable in rich media (Additional file 1, Tables S1 and S2) . The strong correlation between H2BK123ub1 levels in coding regions and transcriptional activity suggests that the majority of H2BK123ub1 in the yeast genome is established during transcription elongation through the association of the H2B ubiquitylation machinery with elongating RNA polymerase II .
Intriguingly, the level of H2BK123ub1 also showed a correlation with gene length, with higher levels of H2BK123ub1 present at genes with longer coding regions (Figure 1D). Long genes have a higher probability of containing cryptic transcription initiation sites in their coding regions, and increased levels of H3K36 methylation at these genes have been correlated with the suppressed utilization of these sites [1, 36]. In the case of H3K36 methylation, increased levels of the modification are found predominantly at long genes that are infrequently transcribed . However, H2BK123ub1 is present at long genes with both low and high rates of transcription (Figure 1E, F; Additional file 1, Table S2). We have previously shown that H2BK123ub1 plays a role in restoring chromatin structure displaced by elongating RNAPII . Thus, at long genes, the increased levels of H2K123ub1 could be important for preventing the utilization of cryptic transcription initiation sites by ensuring nucleosome reassembly. Conversely, at shorter genes, where suppression of cryptic transcription is not a general issue, the enrichment of H2BK123ub1 at coding regions could serve additional functions, such as regulating the levels of other transcription-related histone modifications or other events in the transcription process [4–18, 33].
Comparison of the genome-wide patterns of H2B ubiquitylation and H3K4me3, H3K79me2/3, and H3K36me3
Histone H3 is modified by methylation at lysines 4 and 36 at its N terminus and at lysine 79 in the globular domain of H3. These three H3 methylations have been functionally linked to the levels of H2BK123ub1, with H3K4 and H3K79 di- and tri-methylation levels reduced, and H3K36 tri-methylation levels enhanced, in the absence of H2BK123ub1 [15–19, 37]. To investigate the links between H2BK123ub1 and these H3 methylations at the genomic level, we compared the map of H2BK123ub1 occupancy to published genome-wide maps of H3K4me3, K36me3 and K79me2/me3 methylation (Additional file 1, Figure S2, S3). By superimposing the genomic profiles of H2BK123ub1 and H3K4me3  over the average gene and at individual chromosomal loci, we found that the two histone marks occupy distinct regions at yeast genes: H3K4me3 peaks at promoters and 5' coding regions, whereas H2BK123ub1 levels are very low at promoters and peak over coding regions (Additional file 1, Figure S2A, B).
Both H3K36me3 and H2BK123ub1 have been connected to the suppression of internal transcription initiation in gene coding regions [20, 39, 40], and the results discussed above indicate that the two marks are enriched at the coding regions of long genes. The genome-wide and regional distribution patterns of H2BK123ub1 and H3K36me3 are also similar, with each mark absent from transcription start sites and concentrated along coding regions (Additional file 1, Figure S2C, D). However, H3K36me3 is more enriched towards the 3' end of coding regions, consistent with its additional roles in 3' end processing . Finally, a comparison of our H2BK123ub1 profile with a recent analysis of H3K79me2/me3 genome-wide occupancy  confirmed the high degree of co-localization between H2BK123ub1 and H3K79me3 over gene coding regions (Additional file 1, Figure S2E, F; Figure S3B; Figure S4C) and the general anti-correlation between the distribution of H2BK123ub1 and H3K79me2, which is enriched at intergenic regions (Additional file 1, Figure S4B). However, unlike H2BK123ub1, neither H3K79 methylation state correlates closely with transcription rate [8, 38].
H2BK123ub1 is present in introns and exons of ribosomal protein genes
To determine whether chromatin structure contributed to the distribution of H2BK123ub1 at intron-containing genes, we utilized published data sets  to measure nucleosome occupancy in introns and exons of the these genes (Figure 2 D, E, F). Consistent with observations in C. elegans and human [10, 43], we found that nucleosome occupancy in all intron-containing genes was reduced in the intron compared to flanking exons (Figure 2D). In addition, the overall nucleosome occupancy in both the introns and exons of RP genes (Figure 2E) was generally lower than in the introns of non-RP genes (Figure 2F). This is consistent with the observation that highly transcribed genes have lower histone occupancy . We validated this nucleosomal pattern at several intron-containing RP genes using a highly accurate, probe-based quantitative PCR protocol. The results showed that both H2B and H3 levels were reduced in the introns of these genes compared to exons, indicating that overall nucleosome occupancy is lower in introns (Figure 2 G, H).
Next, we addressed whether the sharp peak of H2BK123ub1 observed at the 3' intron-exon boundary of RP genes (Figure 2B) was a special feature of the chromatin structure at these regions. We noted that the average nucleosome occupancy in RP gene introns appeared to be lowest just before the 3' splice site, gradually increasing towards the 3' intron-exon boundary and extending into the 5' region of exon 2 (Figure 2E). By comparing the profiles of H2BK123ub1 and nucleosome occupancy at RP genes (Figure 2B, E), we conclude that the peak of H2BK123ub1 at the 3' intron-exon boundary is probably not the consequence of increased nucleosome occupancy in this region. Thus, our analysis suggests that the 3' intron-exon boundary of RP genes could contain a signal that facilitates H2BK123ub1 deposition or prevents its turnover.
Loss of H2B ubiquitylation alters the distribution of H3K36 trimethylation in introns
The regulation of H3K4, H3K36, and H3K79 methylation by H2BK123ub1 raised the question of whether these marks would also be present in introns as well as exons. The published genome-wide datasets of histone modifications in yeast do not have the resolution to map the distributions of H3K4, H3K36, and H3K79 methylation at intron-exon boundaries . Thus, to address this question, we used conventional ChIP to measure the levels of H3K36me3, H3K4me3, and H3K79me3 at the intron and 3' exon of the same RP genes used to analyze H2BK123ub1 (Figure 3). When compared to the levels of H3 at these same locations, it is clear that these three H3 marks are also enriched in introns (Figure 3B, C, D), similar to the enrichment of H2BK123ub1 (Figure 3A). The analysis of H3K36me3 levels in human and C. elegans suggested that H3K36me3 was enriched in exons and depleted in introns ; however, this difference vanished after H3K36me3 levels were normalized to nucleosome occupancy . This emphasizes that it is essential to consider histone occupancy when comparing histone modification profiles across various genomic regions. Interestingly, in the absence of H2BK123ub1, H3K36me3 levels were significantly reduced in several RP gene introns (RPL16A, RPL16B and RPL22A) and increased in exons (RPL16B, RPL22A, and RPS27B) (Figure 3B). Because both H3K36me3 and H2BK123ub1 mark gene-coding regions, it is tempting to speculate that a dynamic balance between the distributions of the two histone marks contributes to the expression of these intron-containing genes. This notion is consistent with previous findings that altered levels of H2B ubiquitylation affect the levels of H3K36 methylation at the GAL1 gene [16, 18]. Although H3K4me3 and H3K79me3 were also enriched in both introns and exons, we found that H3K4me3 was present at higher levels in introns compared to H3K79me3 (Figure 3C and 3D), perhaps reflecting the 5' bias in the distribution of H3K4me3 coupled with the typical 5' location of RP introns. Together, the data suggest that H2BK123ub1 and its downstream marks of H3K36me3, H3K4me3, and H3K79me3 are part of a chromatin architecture that defines RP gene introns.
H2B ubiquitylation genetically interacts with RNA processing factors
To further address whether H2B ubiquitylation was responsible for the observed bre1Δ interactions with mutations in genes encoding RNA processing factors, we combined an htb-K123R mutation, which abolished H2B ubiquitylation, with the following mutations: a deletion of MUD2, which encodes a component of the pre-mRNA-U1 snRNP [48–50]; a deletion of SAC3 or EDC2, two genes encoding factors with roles in mRNA export and mRNA decapping [51, 52]; and a deletion of LSM1, which encodes a protein involved in degradation of cytoplasmic mRNAs [48–50]. We found that the absence of H2BK123ub1 had no effect on the growth of mud2, edc2Δ, sac3Δ and lsm1Δ mutants (Additional file 1, Figure S5). However, when htb-K123R was combined with deletions of LEA1 and MSL1, whose products facilitate U2 snRNA association with pre-mRNA [29, 30], we found a synthetic lethal phenotype (Figure 4B). A likely interpretation of these genetic interactions is that Bre1-mediated H2B ubiquitylation is functionally linked to U2 snRNP assembly.
The positioning of histone modifications at genes has been associated with numerous co-transcriptional processes ranging from initiation to elongation to 3' end processing. In this report, we investigated the genome-wide distribution of H2B monoubiquitylation in budding yeast. We found that H2BK123ub1, a mark of active transcription, is predominantly localized across gene coding regions, consistent with its postulated roles in transcription elongation [20, 22, 53]. The mark was also proportional to transcription rate, a likely consequence of the association of the H2B ubiquitylation machinery with elongating RNA polymerase II [32, 54, 55]. Both of these features are similar to the distribution of H2B monoubiquitylation in the human genome, supporting the view that both the regulation of this histone modification and its functional roles in transcription are evolutionarily conserved [34, 56]. We also found a correlation between the levels of H2BK123ub1 and gene length, with the mark enriched on long genes. Because long genes frequently contain cryptic transcription initiation sites in their coding regions, this correlation supports the view that the presence of H2BK123ub1 contributes to the restoration of chromatin structure during transcription elongation to suppress the utilization of these sites .
A novel finding of our studies was the discovery that H2BK123ub1 is present in both the introns and exons of the yeast ribosomal protein genes which accounts for the majority of splicing events in yeast cells. Bioinformatic analyses of genome-wide data sets from worms and humans revealed that exons are enriched with active histone modifications while introns are generally depleted of these same marks [5, 7, 9]. This bias has been attributed primarily to chromatin architecture, with nucleosome occupancy higher in exons compared to introns [5–7, 9, 10]. Using published data sets on nucleosome occupancy in yeast , we found that yeast introns also have reduced nucleosome levels in introns compared to exons. However, this difference does not correlate with a reduction in H2BK123ub1 levels in the introns of the RP genes, which account for ~90% of the splicing activity of intron-containing genes . The finding that H2BK123ub1 is a feature of yeast introns is similar to a recent finding, again from bioinformatic analysis of chromatin modification data sets, that H2B ubiquitylation is one of 10 marks of 5' introns in humans . Most yeast genes with introns contain a single intron that is located close to the 5' end of the coding region [57–59]. Thus, the presence of H2B ubiquitylation in the 5' introns of human genes may reflect a similarity in the chromatin architecture of promoter-proximal introns between yeast and human intron-containing genes. We also found that H2BK123ub1 peaked at the 3' intron-exon boundary, particularly in the RP genes. We speculate that the chromatin structure of 3' intron-exon boundaries in these genes could carry a signal that enhances the accessibility of the enzymatic machinery that mediates H2B ubiquitylation (Rad6-Bre1) [27, 28]. Alternatively, enzymes that target H2B for de-ubiquitylation (Ubp8 and Ubp10) [15, 29] could be preferentially prevented from associating with these regions. It is currently unclear what feature of chromatin architecture at these boundaries regulates either the deposition or removal of H2B ubiquitylation. Likewise, it is not known if this distinct chromatin structure plays a role in transcriptional regulation, including pre-mRNA splicing.
H2B ubiquitylation controls the methylation of three lysine residues in histone H3 in trans. A comparison of the genome-wide distributions of H2BK123ub1, H3K4me3, H3K79me2/me3, and H3K36me3 showed that the four modifications occupy distinct regions in genes. As previously reported, H3K123ub1 and H3K79me3 co-localize across coding regions, while H2BK123ub1 and H3K79me2 show an anti-correlation at intergenic regions . H3K4me3 is localized predominantly at 5' gene regions, while H3K36me3, like H2BK123ub1, spreads across coding regions, but with a more pronounced enrichment at the 3' end of genes. How H2BK123ub1 controls these particular distribution patterns remains an area of intense investigation [16, 18, 24, 60–63]. H3K4me3 and H3K79me3, like H2BK123ub1, were also present in the introns of five ribosomal protein genes that were analyzed, and all of these marks were separable from nucleosome occupancy. Together, the results are similar to the reported presence of H2B ubiquitylation and H3K79me3 in the 5' introns of humans .
Unlike the situation in humans and worms [6, 7, 12], H3K36me3 is present in both the introns and exons of yeast genes. The restricted presence of H3K36me3 in exons in higher eukaryotic genomes has been correlated with the regulation of alternative mRNA splicing [6, 7, 28]. These observations suggest two possible roles for H3K36me3 in pre-mRNA splicing, and specifically in the regulation of pre-RNA splicing. First, H3K36me3 could mark exons as a part of a gene structure and along with cis-splicing elements facilitate the decision of whether to include a specific exon [6, 12]. Alternatively, H3K36me3 could act as an anchor site for recruiting splicing factors that regulate alternative mRNA splicing [7, 28]. However, in the yeast genome, the majority of intron-containing genes contain a single intron, and the regulation of splicing efficiency is thus more important than alternative mRNA splicing. It has been suggested that splicing efficiency is a function of the rate of RNA polymerase II elongation. This scenario is supported by a recent finding that RNAP II pauses transiently around the 3'end of introns and that this pause coincides with splicing factor recruitment . Thus, the presence of histone marks in both introns and exons might promote splicing efficiency by controlling RNA polymerase II elongation. The finding that H2BK123ub1 levels are enhanced at 3' intron-exon boundaries could provide a mechanism to couple RNAPII elongation to the recruitment of splicing factors. Moreover, the dynamic relationship between the levels of H2BK123ub1 and H3K36me3 in introns and exons supports a redundant mechanism to ensure optimal RNAP II elongation, in turn promoting efficient pre-mRNA splicing.
Because the loss of H2B ubiquitylation, H3K4/K79 methylation, or H3K36 methylation does not compromise cell viability, the histone modifications cannot play an essential role in splicing. We suggest that the presence of these marks in introns, together with reduced nucleosome occupancy in these regions, are part of a chromatin architecture that facilitates the recognition of exons and introns by splicing regulators. Further support for this mechanism comes from the observation that a synthetic lethal phenotype resulted from combining an htb-K123R mutation with deletions of genes with roles in pre-mRNA splicing, specifically in U2 splicesome assembly. For example, the histone modifications might serve as binding sites for proteins that in turn interact with splicing factors. Such a scenario has been proposed for H3K4me3 and the Chd1 protein, which contains a chromodomain that recognizes the methyl mark and interacts with the U2 snRNP complex in both humans and yeast to promote efficient splicing [65–67].
In summary, the co-transcriptional formation of H2BK123ub1 leads to its spread across gene coding regions. This mark in turn defines the distribution of H3K4 tri-methylation at the 5' end of coding regions and H3K79 tri-methylation toward the center of coding regions during the processes of transcription initiation and elongation. H2BK123ub1 also marks introns, and its presence in these regions leads to the presence of H3K4me3 and H3K79me3 at the introns of several ribosomal protein genes examined. The presence of H2BK123ub1 also influences the distribution of H3K36me3 across coding regions, particularly in introns. As a consequence, we suggest that this coordinated pattern of histone marks along transcribed regions facilitates the recognition of exons and introns and allows for efficient co-transcriptional pre-mRNA processing.
Yeast strains and growth conditions
Strains used in this study are listed in Additional file 1, Table S3. Yeast cultures were grown to mid-log phase in YPD medium (2% yeast extract; 2% peptone; 2% glucose) at 30°C for all ChIP-chip and growth studies.
ChIP was performed as described previously (Kao et al., 2004) with minor modifications. Chromatin pellets from formaldehyde fixed cells were lysed by glass bead vortexing for 30 min at 4°C. For all ChIP experiments analyzing individual genes, cell lysates were digested with 160 U of Micrococcal nuclease (MNase) per 100 ml of cells for 15 min (Nuclease S7, Roche Applied Science, Taiwan). The following antibodies were used with the equivalence of 10 OD units of cell lysate: FLAG, 20 μl (M2-F3165; Sigma-Aldrich, St. Louis, MO); HA, 5 μl (3F10; Roche Applied Science, Taiwan); H3, 4 μl (ab1791; Abcam, Cambridge, England); H3K4me3, 2 μl (ab8580; Abcam); H3K79me3, 10 μl (ab 2621, Abcam); and H3K36me3, 4 μl (ab9050, Abcam). Antibodies were prebound to protein A or protein G sepharose or dynabeads®. Purified DNAs were analyzed by probe-based real-time quantitative PCR on a Roche LightCycler® 480 Real-Time PCR System. Probe numbers are listed in Additional file 1, Table S4 and probes were ordered from Roche Universal Probe Library, Roche Applied Science, Taiwan. The IP/Input ratios of H2BK123ub1, H3K4me3, H3K36me3, and H3K79me3 were normalized to the IP/Input ratio of INT-V sequences as described previously .
ChDIP-chip hybridization to oligonucleotide arrays
Chromatin pellets from formaldehyde fixed cells were lysed by glass bead vortexing for 30 min at 4°C, and the cell lysates were sonicated on a Branson Sonifier. ChDIP to measure H2BK123ub1 was performed as described previously [20, 31] using a 2 step sequential ChIP. Briefly, the first ChIP was performed using anti-FLAG antibody, and the immune complexes were eluted with FA-lysis buffer containing 200 μg/ml of 3 × FLAG peptide. One-tenth of the eluate from the first ChIP was reserved for "input", and anti-HA antibody (12C5A; Roche) was added to a final concentration of 15 μg/ml to the remaining eluate. The immune complexes were then eluted with 1% SDS/50 mM Tris (pH 8.0) and represent "IP". After reversal of the crosslinks and purification, the IP and input DNA were amplified according to the Affymetrix protocol. IP and input samples from two biological replicates were hybridized to an Affymetrix 1.0R S. cerevisiae microarray, which comprised over 3.2 million probes covering the entire genome at 5 bp resolution.
ChDIP-chip data analysis
Preprocessing and normalization
Two biological replicates for each paired ChIP (HA ubiquitin modified Flag-H2B; denoted by H2Bub1) and control (the level of H2B genomewide) were included in the experiments. Signals for ChIP-chip were extracted by the TileMap algorithm , in which perfect match (PM) only intensities were used. For each probe, the log ratio of intensities of ChIP and control were normalized using the quantile normalization method .
Averaged fold enrichments of all genes, intergenic regions and other regions of interest (in log2 scale)
These averaged log ratios for all probes (including 4802 genes with assigned transcription rates) were allocated, via a code, to a genomic locus (ORF) and a 5' or 3' intergenic region, which were mapped evenly to the 20th-60th, 1th-20th, 60th -80th bins on the x-axis, respectively. In each bin, we averaged all probes of each gene, then averaging over all genes to result in the averaged intensity, which was plotted on top of the 5' intergenic, 'Average Gene' and the 3' intergenic regions to depict their fold enrichments. The averaged fold enrichments of whole genome, coding regions (ORF), intergenic regions and silent chromatin regions: HM mating type cassettes, telomeres and rDNA, were also computed. Similarly, the fold enrichment of all genes grouped by transcription rate, by gene length, and certain subgroups of interest were calculated, then smoothed by a moving window with size 5 (the averaged intensity of every five adjacent bins) , and plotted.
Analysis of intron-containing genes
Averaged fold enrichment of 5'exon-intron-3'exon (in log2 scale)
For intron-containing genes with defined transcription rates, each of their 5' exon (denoted as Exon 1), intron (Intron), and 3' exon (Exon 2) sequences were allocated to the 1st -10th , 10th -40th and 40 th -80th bins, and averaged fold enrichments were calculated over genes. For the few multiple-intron-containing genes, the multiple copies of Exon 1-Intron-Exon 2 sequences were aligned by each Intron and allocated accordingly. Averaged fold enrichments (allocated to the middle point of each bin such as 2.5) were then calculated over genes and smoothed by a moving window with size 5, except for the first two and the last two bins of Exon 1, Intron, and Exon 2.
Averaged nucleosome occupancy
Nucleosome occupancy was defined as the percentage of DNA bound to a nucleosome, following the analysis in . Scaled occupancy levels of nucleosomes in Exon 1 (Intron or Exon 2) ranging from 0 to 100% were weight averaged over 211-intron containing genes, proportional to the length of each Exon 1 (Intron or Exon 2), to yield the averaged nucleosome occupancy.
Raw sequencing data are available at the NCBI Gene Expression Omnibus (GEO)
(Accession number: GSE34325)
Supported by NIH grants GM40118 (M.A.O.) and GM55360 (S.L.B.) and by Leukemia and Lymphoma Society grant LLS 3213-06 (C.F.K.). Research in the laboratories of C.F.K. and G.S.S. was funded by the National Science Council (NSC 98-2320-B-001 -015 -MY3 and NSC 100-2118-M-001-003-MY2) and Academia Sinica, Taiwan.
- Li B, Carey M, Workman JL: The role of chromatin during transcription. Cell. 2007, 128: 707-719. 10.1016/j.cell.2007.01.015.View ArticlePubMedGoogle Scholar
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.View ArticlePubMedGoogle Scholar
- Lee JS, Smith E, Shilatifard A: The language of histone crosstalk. Cell. 2010, 142: 682-685. 10.1016/j.cell.2010.08.011.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith E, Shilatifard A: The chromatin signaling pathway: diverse mechanisms of recruitment of histone-modifying enzymes and varied biological outcomes. Mol Cell. 2010, 40: 689-701. 10.1016/j.molcel.2010.11.031.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersson R, Enroth S, Rada-Iglesias A, Wadelius C, Komorowski J: Nucleosomes are well positioned in exons and carry characteristic histone modifications. Genome Res. 2009, 19: 1732-1741. 10.1101/gr.092353.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Hon G, Wang W, Ren B: Discovery and annotation of functional chromatin signatures in the human genome. PLoS Comput Biol. 2009, 5: e1000566-10.1371/journal.pcbi.1000566.View ArticlePubMedPubMed CentralGoogle Scholar
- Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41: 376-381. 10.1038/ng.322.View ArticlePubMedPubMed CentralGoogle Scholar
- Schulze JM, Jackson J, Nakanishi S, Gardner JM, Hentrich T, Haug J, Johnston M, Jaspersen SL, Kobor MS, Shilatifard A: Linking cell cycle to histone modifications: SBF and H2B monoubiquitination machinery and cell-cycle regulation of H3K79 dimethylation. Mol Cell. 2009, 35: 626-641. 10.1016/j.molcel.2009.07.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Spies N, Nielsen CB, Padgett RA, Burge CB: Biased chromatin signatures around polyadenylation sites and exons. Mol Cell. 2009, 36: 245-254. 10.1016/j.molcel.2009.10.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Tilgner H, Nikolaou C, Althammer S, Sammeth M, Beato M, Valcarcel J, Guigo R: Nucleosome positioning as a determinant of exon recognition. Nat Struct Mol Biol. 2009, 16: 996-1001. 10.1038/nsmb.1658.View ArticlePubMedGoogle Scholar
- Wilhelm BT, Marguerat S, Aligianni S, Codlin S, Watt S, Bahler J: Differential patterns of intronic and exonic DNA regions with respect to RNA polymerase II occupancy, nucleosome density and H3K36me3 marking in fission yeast. Genome Biol. 2011, 12: R82-10.1186/gb-2011-12-8-r82.View ArticlePubMedPubMed CentralGoogle Scholar
- Huff JT, Plocik AM, Guthrie C, Yamamoto KR: Reciprocal intronic and exonic histone modification regions in humans. Nat Struct Mol Biol. 2010, 17: 1495-1499. 10.1038/nsmb.1924.View ArticlePubMedPubMed CentralGoogle Scholar
- Osley MA: Regulation of histone H2A and H2B ubiquitylation. Brief Funct Genomic Proteomic. 2006, 5: 179-189. 10.1093/bfgp/ell022.View ArticlePubMedGoogle Scholar
- Weake VM, Workman JL: Histone ubiquitination: triggering gene activity. Mol Cell. 2008, 29: 653-663. 10.1016/j.molcel.2008.02.014.View ArticlePubMedGoogle Scholar
- Briggs SD, Xiao T, Sun ZW, Caldwell JA, Shabanowitz J, Hunt DF, Allis CD, Strahl BD: Gene silencing: trans-histone regulatory pathway in chromatin. Nature. 2002, 418: 498-10.1038/nature00970.View ArticlePubMedGoogle Scholar
- Henry KW, Wyce A, Lo WS, Duggan LJ, Emre NC, Kao CF, Pillus L, Shilatifard A, Osley MA, Berger SL: Transcriptional activation via sequential histone H2B ubiquitylation and deubiquitylation, mediated by SAGA-associated Ubp8. Genes Dev. 2003, 17: 2648-2663. 10.1101/gad.1144003.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun ZW, Allis CD: Ubiquitination of histone H2B regulates H3 methylation and gene silencing in yeast. Nature. 2002, 418: 104-108. 10.1038/nature00883.View ArticlePubMedGoogle Scholar
- Wyce A, Xiao T, Whelan KA, Kosman C, Walter W, Eick D, Hughes TR, Krogan NJ, Strahl BD, Berger SL: H2B ubiquitylation acts as a barrier to Ctk1 nucleosomal recruitment prior to removal by Ubp8 within a SAGA-related complex. Mol Cell. 2007, 27: 275-288. 10.1016/j.molcel.2007.01.035.View ArticlePubMedGoogle Scholar
- Dover J, Schneider J, Tawiah-Boateng MA, Wood A, Dean K, Johnston M, Shilatifard A: Methylation of histone H3 by COMPASS requires ubiquitination of histone H2B by Rad6. J Biol Chem. 2002, 277: 28368-28371. 10.1074/jbc.C200348200.View ArticlePubMedGoogle Scholar
- Fleming AB, Kao CF, Hillyer C, Pikaart M, Osley MA: H2B ubiquitylation plays a role in nucleosome dynamics during transcription elongation. Mol Cell. 2008, 31: 57-66. 10.1016/j.molcel.2008.04.025.View ArticlePubMedGoogle Scholar
- Pavri R, Zhu B, Li G, Trojer P, Mandal S, Shilatifard A, Reinberg D: Histone H2B monoubiquitination functions cooperatively with FACT to regulate elongation by RNA polymerase II. Cell. 2006, 125: 703-717. 10.1016/j.cell.2006.04.029.View ArticlePubMedGoogle Scholar
- Tanny JC, Erdjument-Bromage H, Tempst P, Allis CD: Ubiquitylation of histone H2B controls RNA polymerase II transcription elongation independently of histone H3 methylation. Genes Dev. 2007, 21: 835-847. 10.1101/gad.1516207.View ArticlePubMedPubMed CentralGoogle Scholar
- Chandrasekharan MB, Huang F, Sun ZW: Ubiquitination of histone H2B regulates chromatin dynamics by enhancing nucleosome stability. Proc Natl Acad Sci USA. 2009, 106: 16686-16691. 10.1073/pnas.0907862106.View ArticlePubMedPubMed CentralGoogle Scholar
- Batta K, Zhang Z, Yen K, Goffman DB, Pugh BF: Genome-wide function of H2B ubiquitylation in promoter and genic regions. Genes Dev. 2011, 25: 2254-2265. 10.1101/gad.177238.111.View ArticlePubMedPubMed CentralGoogle Scholar
- Fierz B, Chatterjee C, McGinty RK, Bar-Dagan M, Raleigh DP, Muir TW: Histone H2B ubiquitylation disrupts local and higher-order chromatin compaction. Nat Chem Biol. 2011, 7: 113-119. 10.1038/nchembio.501.View ArticlePubMedPubMed CentralGoogle Scholar
- Moyal L, Lerenthal Y, Gana-Weisz M, Mass G, So S, Wang SY, Eppink B, Chung YM, Shalev G, Shema E, et al: Requirement of ATM-dependent monoubiquitylation of histone H2B for timely repair of DNA double-strand breaks. Mol Cell. 2011, 41: 529-542. 10.1016/j.molcel.2011.02.015.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakamura K, Kato A, Kobayashi J, Yanagihara H, Sakamoto S, Oliveira DV, Shimada M, Tauchi H, Suzuki H, Tashiro S, et al: Regulation of homologous recombination by RNF20-dependent H2B ubiquitination. Mol Cell. 2011, 41: 515-528. 10.1016/j.molcel.2011.02.002.View ArticlePubMedGoogle Scholar
- Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T: Regulation of alternative splicing by histone modifications. Science. 2010, 327: 996-1000. 10.1126/science.1184208.View ArticlePubMedPubMed CentralGoogle Scholar
- Caspary F, Seraphin B: The yeast U2A'/U2B complex is required for pre-spliceosome formation. EMBO J. 1998, 17: 6348-6358. 10.1093/emboj/17.21.6348.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang J, Abovich N, Rosbash M: Identification and characterization of a yeast gene encoding the U2 small nuclear ribonucleoprotein particle B» protein. Mol Cell Biol. 1996, 16: 2787-2795.View ArticlePubMedPubMed CentralGoogle Scholar
- Kao CF, Hillyer C, Tsukuda T, Henry K, Berger S, Osley MA: Rad6 plays a role in transcriptional activation through ubiquitylation of histone H2B. Genes Dev. 2004, 18: 184-195. 10.1101/gad.1149604.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiao T, Kao CF, Krogan NJ, Sun ZW, Greenblatt JF, Osley MA, Strahl BD: Histone H2B ubiquitylation is associated with elongating RNA polymerase II. Mol Cell Biol. 2005, 25: 637-651. 10.1128/MCB.25.2.637-651.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Emre NC, Ingvarsdottir K, Wyce A, Wood A, Krogan NJ, Henry KW, Li K, Marmorstein R, Greenblatt JF, Shilatifard A, Berger SL: Maintenance of low histone ubiquitylation by Ubp10 correlates with telomere-proximal Sir2 association and gene silencing. Mol Cell. 2005, 17: 585-594.View ArticlePubMedGoogle Scholar
- Minsky N, Shema E, Field Y, Schuster M, Segal E, Oren M: Monoubiquitinated H2B is associated with the transcribed region of highly expressed genes in human cells. Nat Cell Biol. 2008, 10: 483-488. 10.1038/ncb1712.View ArticlePubMedGoogle Scholar
- Holstege FC, Jennings EG, Wyrick JJ, Lee TI, Hengartner CJ, Green MR, Golub TR, Lander ES, Young RA: Dissecting the regulatory circuitry of a eukaryotic genome. Cell. 1998, 95: 717-728. 10.1016/S0092-8674(00)81641-4.View ArticlePubMedGoogle Scholar
- Li B, Gogol M, Carey M, Pattenden SG, Seidel C, Workman JL: Infrequently transcribed long genes depend on the Set2/Rpd3S pathway for accurate transcription. Genes Dev. 2007, 21: 1422-1430. 10.1101/gad.1539307.View ArticlePubMedPubMed CentralGoogle Scholar
- Shahbazian MD, Zhang K, Grunstein M: Histone H2B ubiquitylation controls processive methylation but not monomethylation by Dot1 and Set1. Mol Cell. 2005, 19: 271-277. 10.1016/j.molcel.2005.06.010.View ArticlePubMedGoogle Scholar
- Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, et al: Genome-wide Map of Nucleosome Acetylation and Methylation in Yeast. Cell. 2005, 122: 517-527. 10.1016/j.cell.2005.06.026.View ArticlePubMedGoogle Scholar
- Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, Shia WJ, Anderson S, Yates J, Washburn MP, Workman JL: Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005, 123: 581-592. 10.1016/j.cell.2005.10.023.View ArticlePubMedGoogle Scholar
- Keogh MC, Kurdistani SK, Morris SA, Ahn SH, Podolny V, Collins SR, Schuldiner M, Chin K, Punna T, Thompson NJ, et al: Cotranscriptional set2 methylation of histone H3 lysine 36 recruits a repressive Rpd3 complex. Cell. 2005, 123: 593-605. 10.1016/j.cell.2005.10.025.View ArticlePubMedGoogle Scholar
- Warner JR: The economics of ribosome biosynthesis in yeast. Trends Biochem Sci. 1999, 24: 437-440. 10.1016/S0968-0004(99)01460-7.View ArticlePubMedGoogle Scholar
- Jiang C, Pugh BF: A compiled and systematic reference map of nucleosome positions across the Saccharomyces cerevisiae genome. Genome Biol. 2009, 10: R109-10.1186/gb-2009-10-10-r109.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwartz S, Meshorer E, Ast G: Chromatin organization marks exon-intron structure. Nat Struct Mol Biol. 2009, 16: 990-995. 10.1038/nsmb.1659.View ArticlePubMedGoogle Scholar
- Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat Genet. 2004, 36: 900-905. 10.1038/ng1400.View ArticlePubMedGoogle Scholar
- Hwang WW, Venkatasubrahmanyam S, Ianculescu AG, Tong A, Boone C, Madhani HD: A conserved RING finger protein required for histone H2B monoubiquitination and cell size control. Mol Cell. 2003, 11: 261-266. 10.1016/S1097-2765(02)00826-2.View ArticlePubMedGoogle Scholar
- Wood A, Krogan NJ, Dover J, Schneider J, Heidt J, Boateng MA, Dean K, Golshani A, Zhang Y, Greenblatt JF, et al: Bre1, an E3 ubiquitin ligase required for recruitment and substrate selection of Rad6 at a promoter. Mol Cell. 2003, 11: 267-274. 10.1016/S1097-2765(02)00802-X.View ArticlePubMedGoogle Scholar
- Breitkreutz BJ, Stark C, Reguly T, Boucher L, Breitkreutz A, Livstone M, Oughtred R, Lackner DH, Bahler J, Wood V, et al: The BioGRID Interaction Database: 2008 update. Nucleic Acids Res. 2008, 36: D637-640.View ArticlePubMedGoogle Scholar
- Abovich N, Liao XC, Rosbash M: The yeast MUD2 protein: an interaction with PRP11 defines a bridge between commitment complexes and U2 snRNP addition. Genes Dev. 1994, 8: 843-854. 10.1101/gad.8.7.843.View ArticlePubMedGoogle Scholar
- He W, Parker R: Functions of Lsm proteins in mRNA degradation and splicing. Curr Opin Cell Biol. 2000, 12: 346-350. 10.1016/S0955-0674(00)00098-3.View ArticlePubMedGoogle Scholar
- Kress TL, Krogan NJ, Guthrie C: A single SR-like protein, Npl3, promotes pre-mRNA splicing in budding yeast. Mol Cell. 2008, 32: 727-734. 10.1016/j.molcel.2008.11.013.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunckley T, Tucker M, Parker R: Two related proteins, Edc1p and Edc2p, stimulate mRNA decapping in Saccharomyces cerevisiae. Genetics. 2001, 157: 27-37.PubMedPubMed CentralGoogle Scholar
- Fischer T, Strasser K, Racz A, Rodriguez-Navarro S, Oppizzi M, Ihrig P, Lechner J, Hurt E: The mRNA export machinery requires the novel Sac3p-Thp1p complex to dock at the nucleoplasmic entrance of the nuclear pores. EMBO J. 2002, 21: 5843-5852. 10.1093/emboj/cdf590.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim J, Guermah M, McGinty RK, Lee JS, Tang Z, Milne TA, Shilatifard A, Muir TW, Roeder RG: RAD6-Mediated transcription-coupled H2B ubiquitylation directly stimulates H3K4 methylation in human cells. Cell. 2009, 137: 459-471. 10.1016/j.cell.2009.02.027.View ArticlePubMedPubMed CentralGoogle Scholar
- Song YH, Ahn SH: A Bre1-associated protein, large 1 (Lge1), promotes H2B ubiquitylation during the early stages of transcription elongation. J Biol Chem. 2010, 285: 2361-2367. 10.1074/jbc.M109.039255.View ArticlePubMedGoogle Scholar
- Ng HH, Dole S, Struhl K: The Rtf1 component of the Paf1 transcriptional elongation complex is required for ubiquitination of histone H2B. J Biol Chem. 2003, 278: 33625-33628. 10.1074/jbc.C300270200.View ArticlePubMedGoogle Scholar
- Shema E, Tirosh I, Aylon Y, Huang J, Ye C, Moskovits N, Raver-Shapira N, Minsky N, Pirngruber J, Tarcic G, et al: The histone H2B-specific ubiquitin ligase RNF20/hBRE1 acts as a putative tumor suppressor through selective regulation of gene expression. Genes Dev. 2008, 22: 2664-2676. 10.1101/gad.1703008.View ArticlePubMedPubMed CentralGoogle Scholar
- Spingola M, Grate L, Haussler D, Ares M: Genome-wide bioinformatic and molecular analysis of introns in Saccharomyces cerevisiae. RNA. 1999, 5: 221-234. 10.1017/S1355838299981682.View ArticlePubMedPubMed CentralGoogle Scholar
- Davis CA, Grate L, Spingola M, Ares M: Test of intron predictions reveals novel splice sites, alternatively spliced mRNAs and new introns in meiotically regulated genes of yeast. Nucleic Acids Res. 2000, 28: 1700-1706. 10.1093/nar/28.8.1700.View ArticlePubMedPubMed CentralGoogle Scholar
- Clark TA, Sugnet CW, Ares M: Genomewide analysis of mRNA processing in yeast using splicing-specific microarrays. Science. 2002, 296: 907-910. 10.1126/science.1069415.View ArticlePubMedGoogle Scholar
- Vitaliano-Prunier A, Menant A, Hobeika M, Geli V, Gwizdek C, Dargemont C: Ubiquitylation of the COMPASS component Swd2 links H2B ubiquitylation to H3K4 trimethylation. Nat Cell Biol. 2008, 10: 1365-1371. 10.1038/ncb1796.View ArticlePubMedGoogle Scholar
- McGinty RK, Kim J, Chatterjee C, Roeder RG, Muir TW: Chemically ubiquitylated histone H2B stimulates hDot1L-mediated intranucleosomal methylation. Nature. 2008, 453: 812-816. 10.1038/nature06906.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee JS, Shukla A, Schneider J, Swanson SK, Washburn MP, Florens L, Bhaumik SR, Shilatifard A: Histone crosstalk between H2B monoubiquitination and H3 methylation mediated by COMPASS. Cell. 2007, 131: 1084-1096. 10.1016/j.cell.2007.09.046.View ArticlePubMedGoogle Scholar
- Schulze JM, Hentrich T, Nakanishi S, Gupta A, Emberly E, Shilatifard A, Kobor MS: Splitting the task: Ubp8 and Ubp10 deubiquitinate different cellular pools of H2BK123. Genes Dev. 2011, 25: 2242-2247. 10.1101/gad.177220.111.View ArticlePubMedPubMed CentralGoogle Scholar
- Alexander RD, Innocente SA, Barrass JD, Beggs JD: Splicing-dependent RNA polymerase pausing in yeast. Mol Cell. 2010, 40: 582-593. 10.1016/j.molcel.2010.11.005.View ArticlePubMedPubMed CentralGoogle Scholar
- Pray-Grant MG, Daniel JA, Schieltz D, Yates JR, Grant PA: Chd1 chromodomain links histone H3 methylation with SAGA- and SLIK-dependent acetylation. Nature. 2005, 433: 434-438. 10.1038/nature03242.View ArticlePubMedGoogle Scholar
- Sims RJ, Millhouse S, Chen CF, Lewis BA, Erdjument-Bromage H, Tempst P, Manley JL, Reinberg D: Recognition of trimethylated histone H3 lysine 4 facilitates the recruitment of transcription postinitiation factors and pre-mRNA splicing. Mol Cell. 2007, 28: 665-676. 10.1016/j.molcel.2007.11.010.View ArticlePubMedPubMed CentralGoogle Scholar
- Gunderson FQ, Johnson TL: Acetylation by the transcriptional coactivator Gcn5 plays a novel role in co-transcriptional spliceosome assembly. PLoS Genet. 2009, 5: e1000682-10.1371/journal.pgen.1000682.View ArticlePubMedPubMed CentralGoogle Scholar
- Ji H, Wong WH: TileMap: create chromosomal map of tiling array hybridizations. Bioinformatics. 2005, 21: 3629-3636. 10.1093/bioinformatics/bti593.View ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Åstrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Gonzalez RC, Woods RE: Digital Image Processing. 2002, Upper Saddle River, NJ: Prentice Hall, 2Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.