Skip to main content
Download PDF

A unique nucleosome arrangement, maintained actively by chromatin remodelers facilitates transcription of yeast tRNA genes

BMC Genomics volume 14, Article number: 402 (2013) Cite this article

Abstract

Background

RNA polymerase (pol) III transcribes a unique class of genes with intra-genic promoters and high transcriptional activity. The major contributors to the pol III transcriptome, tRNAs genes are found scattered on all chromosomes of yeast. A prototype tDNA of <150 bp length, is generally considered nucleosome-free while some pol III-transcribed genes have been shown to have nucleosome-positioning properties.

Results

Using high resolution ChIP-chip and ChIP-seq methods, we found several unique features associated with nucleosome profiles on all tRNA genes of budding yeast, not seen on nucleosome-dense counterparts in fission yeast and resting human CD4+ T cells. The nucleosome-free region (NFR) on all but three yeast tDNAs is found bordered by an upstream (US) nucleosome strongly positioned at −140 bp position and a downstream (DS) nucleosome at variable positions with respect to the gene terminator. Perturbation in this nucleosomal arrangement interferes with the tRNA production. Three different chromatin remodelers generate and maintain the NFR by targeting different gene regions. Isw1 localizes to the gene body and makes it nucleosome-depleted, Isw2 maintains periodicity in the upstream nucleosomal array, while RSC targets the downstream nucleosome. Direct communication of pol III with RSC serves as a stress-sensory mechanism for these genes. In its absence, the downstream nucleosome moves towards the gene terminator. Levels of tRNAs from different families are found to vary considerably as different pol III levels are seen even on isogenes within a family. Pol III levels show negative correlation with the nucleosome occupancies on different genes.

Conclusions

Budding yeast tRNA genes maintain an open chromatin structure, which is not due to sequence-directed nucleosome positioning or high transcription activity of genes. Unlike 5′ NFR on pol II-transcribed genes, the tDNA NFR, which facilitates tDNA transcription, results from action of chromatin remodeler Isw1, aided by Isw2 and RSC. The RSC-regulated nucleosome dynamics at the 3′ gene-end serves as a novel regulatory mechanism for pol III transcription in vivo, probably by controlling terminator-dependent facilitated recycling of pol III. Salient features of yeast tDNA chromatin structure reported in this study can explain the basis of the novel non-transcriptional roles ascribed to tDNAs.

Background

Positioning of the nucleosomes on specific locations with respect to factor binding sites can control the accessibility of underlying DNA and allow the transcription machinery to work effectively in a chromatin environment [1]. Genome-wide maps of nucleosomes in a number of eukaryotes have established that RNA polymerase II promoters are nucleosome- depleted but coding regions may have a regular array of statistically arranged nucleosomes [2–8]. The nucleosome-free region (NFR) just upstream of transcription start site (TSS), results from sequence preferences as well as actions of chromatin remodelers [9, 10].

RNA polymerase (pol) III transcribes a heterogeneous set of small non-coding RNA genes with majority constituted of tRNA genes (70–130 bp). Earlier, active transcription of a tRNA gene was reported to exclude nucleosome from the gene in vivo[11]. Pol III-transcribed genes are found in comparatively nucleosome-depleted intergenic regions [4, 6], which gives the general impression that chromatin may not have much to do with their transcription. Activation of pol III transcription from several sites after global nucleosome depletion in yeast cells, suggests otherwise [12]. Presence of chromatin remodelers RSC and Isw2 on pol III-transcribed genes has also been reported [13, 14]. RSC shows different remodeling activity on pol II and pol III-transcribed genes [15] and regulates pol III-transcribed SNR6 and SUP4 genes [16, 17]. Recent studies with human pol III-transcribed genes have reported that the presence of pol II and pol II-specific chromatin marks near pol III-transcribed genes generates an active chromatin milieu for them [18–20].

This study shows that nucleosome depletion mechanisms at NFR on tDNAs and 5’ NFR on pol II-transcribed genes are fundamentally different. A stronger upstream (US) nucleosome positioning, persistence of NFR under repression on most of the genes, sequence-independence of NFR and requirement of remodeler activities to maintain the chromatin structure in active state, are some of the observations which suggest tDNA NFRs are different from the 5’ NFR on pol II-transcribed genes. An array of Isw2-dependent, regularly spaced nucleosomes especially in gene upstream region reflects stronger US nucleosome barrier as compared to downstream of the gene terminator. We identify Isw1 as the new remodeler for pol III genes, which along with RSC maintains the tDNA NFR. Finally, pol III communication with RSC regulates the downstream nucleosome mobility. This downstream nucleosome dynamics can regulate the terminator-dependent pol III recycling and hence the rate of transcription.

Results

Chromatin structure of tDNAs in active and repressed states

We used high-resolution tiling microarray to map nucleosome positions around budding yeast tDNAs in a 2 kb window. A large number of positioned as well as fuzzy or delocalized nucleosomes (Figure 1A) were seen, which give similar profiles on alignment of the whole data with respect to TSS or TTS (transcription termination site) of all the tRNA genes, probably due to short transcribed length. Conforming to earlier reports [5, 6], all tDNAs except three reside within a nucleosome-depleted region from −70 to +180 bp with respect to TSS (Figure 1B). However, some unique and interesting features associated with the gene-flanking nucleosomes can be noticed. Most tDNAs (82%) show a strongly positioned upstream (−1, US) nucleosome, centered at −140 bp position with respect to TSS while position of the downstream (+1, DS) nucleosome varies across the genes (Figure 1A, B; Additional file 1: Figure S1A). Nucleosomes in the further upstream regions are more regularly spaced compared to those in downstream of the genes. Similar NFR is seen on two other pol III-transcribed genes, RNA170 and ZOD1 (Figure 1C). We had earlier reported a similar arrangement on SNR6 with NFR spanning over its TATA box and A box regions [17]. Such a common scheme of arrangement on rest of the non-tRNA genes or any special arrangement with reference to upstream TATA box or introns found in some of tDNAs is not seen.

Figure 1
figure 1

Nucleosome arrangement around yeast tRNA genes. A) Arrangements of nucleosomes around 273 yeast tRNA genes with respect to TSS (bent arrow) show NFR on gene regions. Cartoon B) gives the summary and typical arrangement of nucleosomes around a tDNA. C) Heat maps showing nucleosome arrangements around two non-tRNA genes RNA170 and ZOD1 and four initiator tRNAMet (IMT) genes, three of which have fuzzy nucleosomes on the gene body (Nuc-plus). D) Average nucleosome profiles are compared between normal growth and starved condition for 52 genes that showed changes. E) In vivo nucleosome profile (this study) is compared with the sequence-based prediction of nucleosome profile [9].

Full size image

Interestingly, nucleosomal arrangement is different on tRNAMet family members. While EMT (elongator) genes show the general NFR-based arrangement of nucleosomes, three genes belonging to IMT (initiator) family show a fuzzy nucleosome covering the 5’ gene ends (Nuc-plus genes, Figure 1C, lower panel). This may be of functional significance, with a specialized role of chromatin in their regulation, as discussed later. We analyzed nucleosome positions on tDNAs in fission yeast and human CD4+ cells also from published datasets [2, 7]. As depicted in Additional file 1: Figures S1B and C, tDNA chromatin in these organisms differs from the budding yeast. With a higher nucleosome density, ~50% of the genes in both species are embedded in nucleosomes and do not show clear NFR, strong US nucleosome or phased array (Additional file 1: Figure S1B,C). The tDNAs in human and mouse show tissue specific expression [20, 21]. Non-uniform nucleosomal arrangement across all the genes could give differential expression of the genes in a tissue-specific manner.

Pol III transcription is repressed under stress conditions [22, 23]. In order to see relevance of the tDNA chromatin structure to transcription, we used the same micro-array to find nucleosomal changes, if any near tDNAs, after prolonged nutrient deprivation. After 4 hr of starvation, nucleosomal changes are seen on 52 genes, within a window of 100 bp upstream and downstream of the gene body (Additional file 1: Figure S1D). An increase in average occupancy of both US and DS nucleosomes along with a small decrease in occupancy of −2 and +2 nucleosomes on both the sides can be noticed (Figure 1D). The average US nucleosome position is not disturbed while the gene terminator is encroached by the DS nucleosome, resulting in shrinkage of average NFR. This observation suggests that the DS nucleosome on the tDNAs either regulates their transcription or high transcription in active state keeps it away from the genes. However, NFR on rest of the genes is not disturbed under repression; suggesting nucleosome-depletion on active tDNAs is not a result of their high transcription rate.

tDNA NFR is not sequence-directed

The 5’ NFR at pol II-transcribed genes is established with a significant contribution of sequence-preferences of nucleosome positions [9]. We asked whether NFR seen at tDNAs is directed because of sequence preferences of flanking nucleosomes. Sequence-based average nucleosome occupancy profile, as predicted (trained on real data) by Kaplan et al. around tDNAs [9], is opposite of the in vivo nucleosome profile found in this study (Figure 1E), indicating that nucleosome depletion on tDNAs is not because of underlying DNA sequences. Nucleosome density data from genome wide in vitro nucleosome assembly [24] also supports this observation and negates the possibility that variable positions of DS nucleosome are because of positioning signals (Additional file 1: Figure S1E,F). These observations suggest that nucleosome arrangement around tDNAs in vivo is maintained by active mechanisms, such as ATP-dependent chromatin remodeling, which help over-ride the sequence-directed nucleosome positioning information.

Isw1 targets tDNAs and maintains tDNA chromatin in collaboration with Isw2

Isw2 chromatin remodeler is recruited in upstream regions of tRNA genes by Bdp1 subunit of transcription factor IIIB [25]. We mapped nucleosomes in isw1∆2∆ cells and used the published data [26, 27] on remodelers as well as nucleosomes occupancy on tRNA genes in single deletion mutants (Figure 2; Additional file 2: Figure S2) to find their effect on these genes. A total of 112 genes show nucleosome gain (77 genes) or loss (35 genes) within 100 bp around the gene ends (cf. Figure 2A,B). Moreover, most of the genes show perturbation in upstream nucleosome array in isw1∆2∆ cells (Figure 2B,C; Additional file 2: Figure S2A), probably because of Isw2 absence, since Isw2 mutation has been reported to disrupt chromatin structure in upstream regions of few tRNA genes [25] and in isw2∆ cells, a disturbance in nucleosome array on tRNA genes is seen only in the upstream region (Additional file 2: Figure S2B). We found that majority of Isw2-occupied tDNAs are also enriched for Isw1 (Additional file 2: Figure S2C) but on different gene regions (Figure 2D). Only ~60% of 260 Isw2-occupied genes show nucleosome gain or loss in isw1∆2∆ (Figure 2C; Additional file 2: Figure S2D). In contrast, most (249 out of 262) of the Isw1-occupied genes show nucleosome gain over gene body in isw1∆ cells, without disruption of flanking nucleosomal periodicity (Figure 2E; Additional file 2: Figure S2E). This observation indicates that Isw1 may be the major remodeler responsible for nucleosome-free status of tDNAs in wild type cells. It is intriguiging that when both Isw1 and Isw2 are deleted, the double deletion mutant shows nucleosome gain on far less number of genes (Figure 2E). This suggests that individual effects of deleting the remodelers are nullified on many genes probably due to the opposite directionality of remodelers’ actions. Isw2 aligns and pushes nucleosome from upstream towards TSS while Isw1 opposes its action by pushing the nucleosome away from the gene body. As a result of combined but antagonistic actions of Isw1 and Isw2, a nucleosome free region is created on the gene body along with strong positioning of the US nucleosome.

Figure 2
figure 2

ISWI maintains tDNA chromatin. Heat maps in A) and B) compare the nucleosomes on 112 genes, which show nucleosome changes in isw1∆2∆ cells compared to wild type cells. C) Comparison of average nucleosome occupancy profiles under two conditions as in panels A and B. D) Isw1 and Isw2 occupy different regions on yeast tDNAs [13, 26]. E) Venn intersections among Isw1 target genes and the genes that show nucleosome gain on gene body either in isw1∆ or isw1∆2∆ cells.

Full size image

RSC-pol III communication maintains the active state chromatin structure

RSC is enriched on many tRNA genes [14] and compromising its catalytic activity was shown to cause increase in nucleosome density on these genes [15]. A mutation in the Rsc4 subunit of RSC is reported to abolish its interaction with pol III [28]. We had earlier shown that single nucleosomes near SNR6 and SUP4 are affected in localized manner in rsc4-∆4 mutant [16, 17]. Nucleosome mapping in this mutant revealed a large number of nucleosomal shifts (120 genes) on tDNAs (Figure 3A). Nucleosome in the mutant covers even the gene body in several cases (Figure 3A, right panel), as seen in the example of tA(AGC)D, validated by the ChIP-Real Time PCR measurement (Figure 3B). On the average, the US nucleosome position is not disturbed and most of the shifts are seen in the downstream region, where DS nucleosome encroaches the gene region while NFR shrinks (Figure 3C). Unlike isw1∆2∆ mutant, periodicity in flanking regions is largely unaffected but nucleosomes show lower average occupancy on both sides of the NFR (Figures 3C). A similar gain in density of nucleosomes near pol III-transcribed genes upon loss of Sth1, the catalytic subunit of RSC was previously reported [15]. Pol III is not lost from the genes in this mutant as exemplified by three of the affected genes (Figure 3D). Similarly, the catalytic subunit Sth1 continues to occupy the genes in the rsc4-∆4 mutant (Figure 3E), as reported even earlier [28]. This data indicates that the dynamics of the DS nucleosome seen in this mutant is guided by a cross-talk between RSC and pol III while bound to the genes.

Figure 3
figure 3

NFR shrinkage and DS nucleosome dynamics in rsc4-∆4 mutant. A) Nucleosome heat maps of 120 genes showing nucleosome changes in rsc4-∆4 mutant as compared to the wild type cells. B) Real Time ChIP-PCR validation of a nucleosome encroaching the otherwise NFR, tA(AGC)D gene body. C) Comparison of the averaged nucleosome occupancy for 120 genes shown in the panel A. D) RPC160-myc (pol III) occupancy does not differ significantly in the two strains. E) RSC (Sth1-myc) occupancy shows somewhat higher levels in the mutant.

Full size image

tDNA chromatin is required for optimum transcription

We enquired whether nucleosomal perturbations seen in different conditions affect transcription by measuring cellular tRNA levels using Real Time qPCR analysis. Based on anti-codon sequence, budding yeast tRNA genes are grouped into 42 families, with highly similar sequences of multiple isogenes. Quantification of tRNA levels for 33 families in wild type cells (Additional file 3: Figure S3A; Additional file 4: Table S1) shows different intra-cellular tRNA levels across the families, which match with the codon usage of yeast cells [29], giving tRNAGlu at the highest level, followed by tRNAGly (Additional file 3: Figure S3A). As compared to tRNA levels for these families (in a descending order) for wild type cells, the corresponding levels show gross perturbations under repressed or remodeler mutation conditions, like isw1∆2∆ and rsc4-∆4 cells (Figure 4A). Some families show higher levels in mutants probably because both gain and loss of nucleosomes on genes are observed in these mutants. Transcription of tDNAs is not affected due to Isw2 loss [13] probably because Isw2 on tDNAs maintains mainly upstream nucleosome periodicity (Figure 2). This indicates that the tRNA changes seen in isw1∆2∆ cells are because of Isw1 absence, which targets the gene body. On the other hand, transcriptional changes in rsc4-∆4 mutant could be due to changes in downstream nucleosome dynamics. This data shows nucleosome dynamics has profound influence on tRNA production.

Figure 4
figure 4

Differential Pol III Occupancy and tRNA gene expression. A) Heat map of tRNA levels showing perturbation of tRNA levels in wild type starved and remodeler mutants as compared to wild type active state levels. B) Heat maps show pol III tag counts binned around TSSs of all tDNAs in active and starved conditions. C) Scatter plot for nucleosome occupancy (−60 to TSS region) and pol III tag counts (−100 to +100 region) for isogenes of tP(UGG) family shows negative correlation. D) Pol III occupancy compared to gene-flanking nucleosome occupancies is shown for the three members of the tRNAGly, tG(UCC) family. Genes denoted by G, N and O are color coded and US, DS nucleosomes are marked. Highest level of pol III on tG(UCC)N is complemented by the lowest nucleosome density and longest NFR.

Full size image

Pol III shows differential enrichment on isogenes

For the intra-cellular levels of tRNA species to be according to the copy number of the coding tRNA gene [29]; every gene copy within a family should be transcribed and occupied by pol III similarly. Different strengths of A and B-boxes may also result in different pol III levels on tDNAs. We used ChIP-seq method to measure the pol III level on different genes. Despite the sequence redundancy within and across the families, ChIP-seq data analysis with or without reallocation of non-unique reads [30] does not give significantly different results (Additional file 3: Figure S3B), probably because the average ChIP fragment size (~180 bp) is longer than the tRNA gene length in our experiment. No pol III on pseudo tRNA genes and high but differential pol III enrichment on all tRNA genes (Additional file 3: Figure S3C; Additional file 4: Table S2) is found. Pol III factors, TFIIIB and TFIIIC persist on genes after starvation [22] but a drastic loss of pol III is seen on all the genes (Figure 4B). Surprisingly, total pol III levels do not match with the relative RNA levels from the gene pools covered by each primer set (Additional file 3: Figure S3D). Moreover, different pol III levels on identical gene copies within a family (Additional file 3: Figure S3E,F) indicate different transcriptional activity on different isogenes. From these data, it appears that tRNA abundance (corresponding to copy number), is probably achieved by high transcription from few isogenes and not by equal transcription from all isogenes of a family.

Nucleosome and pol III occupancy levels on tDNAs are inversely correlated

Pol III initiation factor TFIIIB binds at −30 bp from TSS and recruits pol III [31]. Differential pol III levels on isogenes could be because of different 5’ flanking sequences which affect TFIIIB binding [32–35]. Therefore, high nucleosome occupancy in this region would interfere with TFIIIB binding and as a result, pol III recruitment could be affected. Conversely, pol III presence on the genes could exclude nucleosomes. Hence, we compared nucleosome occupancy in upstream region (−60 to TSS) with pol III enrichment over gene region for tDNAs within each family. An inverse relationship (Pearson correlation between −0.5 to −0.9) for 16 families covering 135 genes (barring 1–2 outliers in each) was found. For example, the scatter plot for one of these families, tProUGG with 10 isogenes present on different chromosomes and Pearson correlation of −0.6544, shows greater enrichment of pol III near TSS of isogenes and lower nucleosome occupancy upstream (Figure 4C). In another example of the tGlyUCC family with only three members and Pearson correlation of −0.93, highest pol III levels are found on tG(UCC)N. The gene has lowest US but similar DS nucleosome occupancy when compared to tG(UCC)G and tG(UCC)O, which show comparatively lower pol III levels (Figure 4D). Thus, pol III occupancy (hence transcription) on most of the target genes has an inverse relationship with nucleosome occupancy in immediate upstream, which is not visible if average plots for individual familes are compared. It gets revealed better in studies of individual examples.

DS nucleosome interferes with transcription from individual genes

Nucleosome mapping results show that under different conditions, the US nucleosome is generally not disturbed. In comparison, the DS nucleosome shows greater dynamism (Figure 1D, 3C), suggesting it may have regulatory influence on gene transcription. Five more examples were studied, for which we could design unique primer sets for RNA estimation (Figure 5; Additional file 4: Table S1). Five of them are found on different chromosomes and belong to different families. RNA yield from tP(UGG)O3 is lowest of the five (Figure 5A). Among the rest four, tR(CCG)L and tS(CGA)C RNA levels are higher than tR(CCU)J and tT(CGU)K RNA levels. The US nucleosome occupancy on all of them is found to be similar but pol III levels still show differences. The DS nucleosome occupancy on tS(CGA)C giving higher RNA yield, is lower than on tR(CUU)J with lower RNA yield (Figure 5B). Similarly, tR(CCG)L yields higher RNA and has lower DS nucleosome occupancy than tT(CGU)K (Figure 5C). In comparison to these, on tP(UGG)O3, pol III levels are comparatively low with both US and DS nucleosomes showing high occupancy (Figure 5D), which may be the reason for lowest RNA yield from this gene. A direct negative correlation of the DS nucleosome occupancy and pol III occupancy/RNA yield from other gene cases could not be made due to sequence redundancy problems associated with isogenes as mentioned above, but the differences in relative RNA levels from genes with mostly similar pol III occupancy can most probably be correlated with the DS nucleosome occupancy, as discussed later. It may be noticed in this context that on many genes, the DS nucleosome shows mobility towards the gene under repressed state (Figure 1C).

Figure 5
figure 5

DS nucleosome is inhibitory for transcription. A) Relative RNA levels estimated by using five of the unique primer sets which could be designed for the five examples shown. Panels B-D) Comparison of nucleosome occupancy profiles with similar pol III levels on different genes. US and DS nucleosomes are marked. B) Comparison of nucleosome occupancy (shaded areas) profiles for tR(CCU)J and tS(CGA)C which have similar pol III levels (solid lines). Higher nucleosome density in downstream region and lower RNA levels of tR(CCU)J suggest that higher nucleosome occupancy in downstream region could be inhibitory for transcription. C) Comparison of nucleosome occupancy (shaded areas) profiles for tR(CCG)L and tT(CGU)K which have similar pol III levels. Higher nucleosome density in downstream region and lower RNA levels of tT(CGU)K suggest that higher nucleosome occupancy in downstream region could be inhibitory for transcription. The tT(CGU)K DNA has a shorter NFR and higher nucleosome occupancy even in the upstream region. D) On tP(UGG)O3 gene, the high nucleosome occupancy on both the gene ends may be inhibitory for transcription, despite similar pol III levels. E) Pol III occupancy (fold enrichment over mock sample, from ChIP-seq data) on initiator tRNAMet (IMT) genes is shown. Average with scatter from three independent experiments is plotted.

Full size image

Different nucleosomal arrangement on initiator and elongator methionyl tRNA genes presents another interesting example of pol III gene regulation by nucleosome dynamics. Initiator tRNAMet (IMT, Met1) with four identified isogenes is found at levels higher than the elongator counterpart (Additional file 3: Figure S3A; Met1 vs. Met2). Being initiator, high but regulated levels of this tRNA may be required by the cell for normal protein synthesis. It is interesting that out of four initiator tRNA (IMT) genes, those with a nucleosome covering the gene region show lower pol III enrichment compared to two other sequence identical copies (Figure 5E). Nucleosomes are generally repressive for transcription. Accordingly, lower pol III levels seen on the two initiator family members (tM(CAU)J3 and tM(CAU)O1), both of which have a fuzzy nucleosome covering their 5’ end and a very short NFR (Figure 1C); suggest either no or very low transcription from the two isogenes. The other two family members, tM(CAU)E and tM(CAU)P; which have higher pol III levels (Figure 5E) and more clear gene regions (Figure 1C) can probably meet the normal cellular requirements. At the same time, losing the nucleosome from the repressed copies provides a possibility of enhanced yield required under conditions of higher demand.

Above examples imply that the nucleosome occupancy and dynamics influence the pol III level and transcription from the tRNA gene loci. Different chromatin remodelers target different gene regions. Isw 2 maintains nucleosome arrangements in upstream region of tDNAs and Isw1 keeps the gene body nucleosome-depleted. RSC action is seen mostly in the downstream region, although the role of other remodelers can not be excluded at this stage. The DS nucleosome mobility towards the gene may inhibit the transcription by restricting the access of pol III to the gene terminator. Therefore, RSC activity keeps the DS nucleosome away from the TTS of the gene.

Discussion

Sequence-directed nucleosome positions can be further modulated by DNA binding proteins and chromatin remodelers [36]. Nucleosome arrangement at tDNAs follows the same principles. Short pol III genes with intra-genic factor binding sites are mostly a size of single nucleosomal DNA. Transcription machinery of pol III has to compete with nucleosomes for access to the tRNA gene sequence, which favors nucleosome formation [9]. Protection of complete gene region by its transcription factors and high transcription rate have been implied as cause of nucleosome depletion at tDNAs. Persistence of tDNA NFR even after prolonged repression suggests that high transcription rate is not required for nucleosome depletion. Nucleosomes are rather actively excluded from the gene body by employing chromatin remodelers. The NFR at tDNAs is an elegant example where nucleosome-favoring sequences are depleted of nucleosomes in vivo by employing chromatin remodelers to facilitate transcription.

In yeast, several ATP-dependent chromatin remodelers are shown to direct the positions of most of the nucleosomes [26]. This study has shown the importance and novelty of the roles of two chromatin remodelers in tDNA regulation. Nevertheless, the nucleosome-depleted status of several gene bodies even in the ISWI double deletion mutant suggests the involvement of some other remodelers as well. We extracted and analysed the data on Chd1 and nucleosome occupancy on tDNAs from published studies [26]. Chd1 showed strong association with tDNAs but nucleosome profiles in wild type and chd1∆ cells did not show discernible differences making the analysis results inconclusive. Isw2 activity in upstream region of tDNAs has been reported earlier [25, 37], where it is found enriched [13]. Our study with double deletion mutant helped us find that Isw1 and Isw2 work in opposite manner on tDNA chromatin and it is the Isw1 action, which makes tDNAs nucleosome-depleted. The opposing actions of Isw1 and Isw2 position the US nucleosome strongly in most cases, which strongly supports the Isw2-dependent nucleosome periodicity in the further upstream region. The strong positioning of US nucleosome restricts the 3’ boundary of the upstream nucleosomal array to −70 bp, which could keep the −30 bp position nucleosome-depleted. This would facilitate TFIIIB binding and promote transcription in turn. RSC activity on the other hand, is seen in the gene downstream region. The earlier reported effects of RSC on tDNAs as judged by transcription of few genes, have been found in a condition where its catalytic subunit is lost [15]. We have used a RSC mutant in which RSC and its catalytic activity remain intact yet remodeling activity is lost due to its deficiency in interaction with pol III. Nucleosome positions and transcription output change across most of the genes in this mutant. This links RSC directly with pol III transcription and gives evidence in support of the mechanism by which RSC controls the tDNA activity and chromatin structure via modulating downstream nucleosome dynamics.

Regulation via nucleosome dynamics is apt for genes that are transcribed at very high rate since it can allow a quick regulatory response in adverse conditions and transcription can resume as soon as repressive condition is over. The interaction of pol III with RSC plays an important role in this by negotiating with the downstream nucleosome. Occluding the whole gene body and DS nucleosome dynamics offer novel regulatory mechanisms, adopted by a majority of pol III-transcribed genes in yeast. Transcription terminator facilitates the recycling of pol III on the same template giving faster re-initiation and hence transcription rate in vitro [38]. Blocking the gene terminator with a nucleosome can reduce the transcription by inhibiting the terminator-dependent recycling of pol III. In the absence of any known classical activation route available for these genes, facilitated re-initiation could be a very important mode of enhancing transcript yield in vivo. DS nucleosome mobility was shown to control the accessibility of TTS and transcription on the SUP4 gene in vivo [16]. Thus, a differential dynamics of DS nucleosome on different tDNAs and their isogenes can lead to different transcript yields from and within a family, as observed in this study. Such a mechanism of gene regulation from downstream is not known for other classes of genes and could be a mechanism specially adapted by and evolved at pol III-transcribed loci. A DS nucleosome covering almost half of the gene body for most of tRNA genes in S. pombe may be significant in this connection (Additional file 1: Figure S1C).

The presence of a strong upstream nucleosome and NFR on tDNAs, which are found scattered across the genome may play a significant role in barrier functions attributed to tRNA genes [39] as nucleosome-free regions have been shown to serve as boundary elements [40]. Replication forks are reported to stall with high frequency near tRNA genes [41], probably due to the well-positioned US nucleosome on tDNAs seen in this study. The study enhances our perception about how nucleosome dynamics participates in pol III transcription regulation and how their arrangements on different gene classes have evolved to match the requirements.

Conclusions

It is fully evident from this data that a strongly positioned upstream (US) nucleosome, the NFR on tDNAs and nucleosome downstream of the gene are not sequence-directed, rather the result of a fine balance between specific activities of the chromatin remodelers RSC and ISWI. Disturbance in this arrangement interferes with transcript yield but NFR on most of the genes is not lost under repression, suggesting the cause of nucleosome depletion is not the high transcriptional activity of tRNA genes. Collaboration between polymerase and ATP-dependent chromatin remodelers is important for optimum transcription from pol III-transcribed genes too. Involvement of other chromatin modifiers and histone modifications remains implicit. The strongly positioned US nucleosome and NFR on tDNAs can form the basis of the barrier functions reported for tRNA genes.

Methods

Strains and growth conditions

Details of yeast strains used in this study are given in Additional file 4: Table S3 [28, 42–44]. Cells were grown to mid log phase (A600 ~0.8) in 100 ml 1X YEPD medium at 30°C. For setting up the repression, cells were peletted down in mid-log phase and shifted to pre- warmed 0.15X YEP medium without any carbon source and incubated for 4 hours [22].

Tiling Array Design

A high density tiling array (244 K Chip) of Agilent platform containing 60-mer probes was designed. Probes covering all pol III-transcribed genes (70–120 bp length) with 1 kb upstream and 1 kb downstream gene regions were tiled at a resolution of 4 bases, giving a 56 bases overlap with the next probe. Non-unique probes that fell within the gene regions were not filtered out to prevent the possible loss of signal in all tRNA gene regions with highly similar sequences. As a corrective measure, the probes in gene and 200 bp immediate upstream or downstream regions were tiled with 59 base overlap instead of 56 bases, generating even higher resolution of 1 bp on and around the genes. This could allow us to identify the false nucleosomal signal, since any signal that covers only the gene region and not the flanks could be because of cross hybridization. A true nucleosomal signal would cover either upstream or downstream flank of the gene or both.

Nucleosomal DNA preparation and array hybridization

Mono-nucleosomal and size-matched naked DNA samples were hybridized to the high resolution, custom-designed micro-array in triplicates. Mono-nucleosomal DNA was prepared by digesting formaldehyde-fixed chromatin with MNase according to previously described method [45] with minor changes. Briefly, spheroplasts were made from fixed cells in mid log phase, subjected to controlled MNase digestion and digested DNA was purified and checked on 1.25% agarose gels (Additional file 5: Figure S4A). To serve as a control, naked genomic DNA was digested with MNase such that it gave a fragment distribution ranging from 100–300 bp (Additional file 5: Figure S4B). The band corresponding to mono-nucleosomal DNA was cut out and purified from the gel. Equal amounts of mono-nucleosomal DNA and digested genomic DNA (2 μg) were labeled with Cy5 and Cy3 respectively, mixed and hybridized to array using manufacturer’s protocol. Post hybridization, slides were scanned with Agilent scanner and data extracted with Agilent Feature extraction software.

Analysis of microarray data

Raw data from feature extraction files was imported into Agilent Chip Analytics software (Agilent Genomic Workbench Lite Edition 6.5). Signal intensities in both channels were normalized against the median of intensities at blank spots followed by intra-array intensity dependent Lowess normalization. Biological replicates showed good reproducibility as indicated by Pearson correlation and box-plots for log2 intensity ratios (Additional file 5: Figure S4C). Log2 intensity ratios from replicates were averaged for each probe and imported to MATLAB for peak calling. A previously used hidden Markov model [45] was applied for detection of nucleosomes through MLM package [46] in MATLAB, which suited well for our experiment and design as it is faster and allows nucleosome detection from incontiguous regions tiled in our array. Signals were smoothed by 3 point smoothing window and nucleosome peaks were called in different datasets. Nucleosomal peaks were then binned around TSS and terminator of the genes in question in 20 bp bins. Matrix files resulting from binning were used to make heat maps and clustering in Multiple Experiment Viewer [47]. Out of 275 tRNA genes, probes in microarray for two of them, tK(CUU)C and tM(CAU)C failed. These genes are located close together and we found aberrant hybridization signals from this region in a previous genome wide study too [5], suggesting the DNA sequence in this region may be problematic (data not shown). The bin-wise averages of log2 intensity ratios (nucleosomal DNA/genomic DNA) in the given window from microarray data were plotted against the distance from the reference point (TSS or TTS).

Validation of Nucleosome shifts

Nucleosomal and ChIP DNA from independent preparations was quantified by Real Time quantitative PCR (ABI Prism) for validation of nucleosome shifts in different conditions. Genomic DNA sample served as a common calibrator. Primers were designed for selected genes (Additional file 4: Table S4) and all samples were diluted to same concentration (5 ng/ml). Data was normalized to a region in telomere at the right arm of chromosome 6 (TELVIR), as this region does not show any significant change in nucleosome occupancy under the conditions in question.

Pol III ChIP and deep sequencing

Samples were prepared from 100 ml cultures of RPC128-FLAG tagged strain grown to A600 0.8-0.9 for pol III ChIP as described earlier [16]. Sepharose CL4B beads served as no-antibody control (mock) while, Anti-FLAG M2-Agarose was used to immuno-precipitate the DNA cross linked to FLAG tagged RPC128. Quality of ChIP samples prepared was checked with quantitative real time PCR on known target genes. Libraries for single end sequencing were made according to standard protocol from Illumina for mock and IP samples. Two biological replicates for IP and one for mock sample in each condition were sequenced using Illumina GA-II platform.

Pol III ChIP-Seq data analysis

Raw data on 36 bp reads from single end sequencing were preprocessed with FASTX-Toolkit [48]. Reads having a quality score of at least 20, on at least 80% of the bases, were exported for further analysis. Reads were aligned to the yeast genome (sacCer3 assembly) with BOWTIE [49], allowing only unique and best alignments for each read. Reads were extended to average ChIP-fragment size (~180) and coverage on each base in the genome was calculated. Read coverage across the genome was scaled to one million and normalized to the total number of reads in each dataset so that comparisons between the samples can be done. Pol III peaks were called in HOMER package [50]. Reproducibility of replicates was checked with USeq package [51] using a scanning window of 500 bp with a step size of 250 bp and a minimum window score of 5.

Considering the very high degree of sequence similarity among tRNA genes, we aligned the reads to yeast genome with BWA [52] allowing 3 mismatches and multiple alignments for non-unique reads (multi-reads), up to 20 places in the genome. Next, multi-reads were reallocated to the most probable target positions according to a previously described algorithm [30] followed by calculation of coverage across the genome and peak calling by HOMER. This algorithm assigns non-unique reads to one of the candidate mapped locations, based on which location gets higher number of uniquely mapped reads. Comparing results from two strategies of alignment, we did not find any significant difference in the results which is attributable to a very small size of pol III target genes in comparison to the typical ChIP fragment size used.

Nucleosome maps in fission yeast and human CD4 cells

Raw sequence reads from MNase-seq experiments [2, 7] were quality filtered as described above. Filtered reads were aligned to respective genomes allowing 2 mismatches and only uniquely aligning reads were considered. Aligned reads were extended to 150 bp and read coverage for each base in the genome was calculated. Nucleosomal peaks were called with GeneTrack [53] (sigma: 20; exclusion: 147). Heat maps in Additional file 1: Figure S1B and C show nucleosomal peaks binned around tRNA genes. Sequences and annotations used for S. pombe were downloaded from PomBase [54] as on 01/09/2012. For human dataset, “hg18” assembly from UCSC was used.

Nucleosome changes in Isw1 and Isw2 mutants

For comparison of nucleosome density in wild type and isw1∆ cells we used Raw MNAse-seq reads from published data [27]. The dataset had reads from S. cerevisiae and S. paradoxus. We aligned the sequencing reads competitively to both species by indexing both the genomes together using Bowtie. We made sure that a read was aligned only to one of the two genomes. Alignments to budding yeast were taken and genome wide coverage was calculated for each base as described above. Data for isw2∆ cells was obtained from the authors’ website [37] and bin-wise averages of normalized log2 intensity ratios (nucleosomal DNA/genomic DNA) for tRNA genes were plotted.

Analysis of Isw1 ChIP-chip data

2-color ChIP-chip data for Isw1 ChIP was taken from a previous study [26]. Data from three replicates were quantile normalized and log transformed using biotoolbox [55]. Normalized log2 ratios were binned around the tRNA genes TSS in appropriately sized bins and data matrix used for generating heat maps or average intensity plots.

Quantification of tRNAs by real time PCR

Total RNA was extracted from 100 ml mid log phase cultures by acidic hot phenol method followed by DNase I treatment. RNA was then subjected to poly-A tailing with poly-A polymerase (Epicenter) using manufacturer’s protocol. First strand cDNA synthesis was performed by AMVRT primed with oligo-dT (18-mer) oligo followed by degradation of template RNA with RNAse A. The cDNA was extracted with Phenol:Chloroform:Isoamyl- alcohol, ethanol precipitated and quantified by Real Time qPCR using ∆∆Ct method. An equivalent amount of poly-A tailed RNA was kept as a control of reverse transcription reaction to make sure DNA contamination did not affect the quantification. U4 snRNA, a pol II transcript was used as internal control [17]. A common genomic DNA sample (5 ng/mL) was used as calibrator so that tRNA expression values could be compared between the gene families, without a complication due to copy number variation. Close to 100 specific primers for tRNA genes were designed with PRISE [56] covering all 42 families and tested for specificity and efficiency in real-time q-PCR using serial dilutions of cDNA templates. Finally 33 primer sets could be validated for use (Additional file 4: Table S1). Most of the primer sets are not unique to isogenes, and measure the transcript and mature tRNA levels together.

References

  1. Bell p, Tiwari VK, Thoma NH, Schubeler D: Determinants and dynamics of genome accessibility. Nat Rev Genet. 2011, 12 (8): 554-564. 10.1038/nrg3017.

    Article  CAS  PubMed  Google Scholar 

  2. Givens RM, Lai WK, Rizzo JM, Bard JE, Mieczkowski PA, Leatherwood J, Huberman JA, Buck MJ: Chromatin architectures at fission yeast transcriptional promoters and replication origins. Nucleic Acids Res. 2012, 40 (15): 7176-7189. 10.1093/nar/gks351.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Jansen A, Verstrepen KJ: Nucleosome positioning in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2011, 75 (2): 301-320. 10.1128/MMBR.00046-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat Genet. 2004, 36 (8): 900-905. 10.1038/ng1400.

    Article  CAS  PubMed  Google Scholar 

  5. Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: A high-resolution atlas of nucleosome occupancy in yeast. Nat Genet. 2007, 39 (10): 1235-1244. 10.1038/ng2117.

    Article  CAS  PubMed  Google Scholar 

  6. Mavrich TN, Ioshikhes IP, Venters BJ, Jiang C, Tomsho LP, Qi J, Schuster SC, Albert I, Pugh BF: A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome. Genome Res. 2008, 18 (7): 1073-1083. 10.1101/gr.078261.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132 (5): 887-898. 10.1016/j.cell.2008.02.022.

    Article  CAS  PubMed  Google Scholar 

  8. Shivaswamy S, Bhinge A, Zhao Y, Jones S, Hirst M, Iyer VR: Dynamic remodeling of individual nucleosomes across a eukaryotic genome in response to transcriptional perturbation. PLoS Biol. 2008, 6 (3): e65-10.1371/journal.pbio.0060065.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, et al: The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009, 458 (7236): 362-366. 10.1038/nature07667.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Yen K, Vinayachandran V, Batta K, Koerber RT, Pugh BF: Genome-wide nucleosome specificity and directionality of chromatin remodelers. Cell. 2012, 149 (7): 1461-1473. 10.1016/j.cell.2012.04.036.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Morse RH, Roth SY, Simpson RT: A transcriptionally active tRNA gene interferes with nucleosome positioning in vivo. Mol Cell Biol. 1992, 12 (9): 4015-4025.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Guffanti E, Percudani R, Harismendy O, Soutourina J, Werner M, Iacovella MG, Negri R, Dieci G: Nucleosome depletion activates poised RNA polymerase III at unconventional transcription sites in Saccharomyces cerevisiae. J Biol Chem. 2006, 281 (39): 29155-29164. 10.1074/jbc.M600387200.

    Article  CAS  PubMed  Google Scholar 

  13. Gelbart ME, Bachman N, Delrow J, Boeke JD, Tsukiyama T: Genome-wide identification of Isw2 chromatin-remodeling targets by localization of a catalytically inactive mutant. Genes Dev. 2005, 19 (8): 942-954. 10.1101/gad.1298905.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Ng HH, Robert F, Young RA, Struhl K: Genome-wide location and regulated recruitment of the RSC nucleosome-remodeling complex. Genes Dev. 2002, 16 (7): 806-819. 10.1101/gad.978902.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Parnell TJ, Huff JT, Cairns BR: RSC regulates nucleosome positioning at Pol II genes and density at Pol III genes. EMBO J. 2008, 27 (1): 100-110. 10.1038/sj.emboj.7601946.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Mahapatra S, Dewari PS, Bhardwaj A, Bhargava P: Yeast H2A.Z, FACT complex and RSC regulate transcription of tRNA gene through differential dynamics of flanking nucleosomes. Nucleic Acids Res. 2011, 39 (10): 4023-4034. 10.1093/nar/gkq1286.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Arimbasseri AG, Bhargava P: Chromatin structure and expression of a gene transcribed by RNA polymerase III are independent of H2A.Z deposition. Mol Cell Biol. 2008, 28 (8): 2598-2607. 10.1128/MCB.01953-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Barski A, Chepelev I, Liko D, Cuddapah S, Fleming AB, Birch J, Cui K, White RJ, Zhao K: Pol II and its associated epigenetic marks are present at Pol III-transcribed noncoding RNA genes. Nat Struct Mol Biol. 2010, 17 (5): 629-634. 10.1038/nsmb.1806.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Moqtaderi Z, Wang J, Raha D, White RJ, Snyder M, Weng Z, Struhl K: Genomic binding profiles of functionally distinct RNA polymerase III transcription complexes in human cells. Nat Struct Mol Biol. 2010, 17 (5): 635-640. 10.1038/nsmb.1794.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Oler AJ, Alla RK, Roberts DN, Wong A, Hollenhorst PC, Chandler KJ, Cassiday PA, Nelson CA, Hagedorn CH, Graves BJ, et al: Human RNA polymerase III transcriptomes and relationships to Pol II promoter chromatin and enhancer-binding factors. Nat Struct Mol Biol. 2010, 17 (5): 620-628. 10.1038/nsmb.1801.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Canella D, Bernasconi D, Gilardi F, Lemartelot G, Migliavacca E, Praz V, Cousin P, Delorenzi M, Hernandez N: A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver. Genome Res. 2012, 22 (4): 666-680. 10.1101/gr.130286.111.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Roberts DN, Stewart AJ, Huff JT, Cairns BR: The RNA polymerase III transcriptome revealed by genome-wide localization and activity-occupancy relationships. Proc Natl Acad Sci U S A. 2003, 100 (25): 14695-14700. 10.1073/pnas.2435566100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Moir RD, Willis IM: Regulation of pol III transcription by nutrient and stress signaling pathways. Biochim Biophys Acta. 2013, 1829: 361-375. 10.1016/j.bbagrm.2012.11.001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Zhang Y, Moqtaderi Z, Rattner BP, Euskirchen G, Snyder M, Kadonaga JT, Liu XS, Struhl K: Intrinsic histone-DNA interactions are not the major determinant of nucleosome positions in vivo. Nat Struct Mol Biol. 2009, 16 (8): 847-852. 10.1038/nsmb.1636.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Bachman N, Gelbart ME, Tsukiyama T, Boeke JD: TFIIIB subunit Bdp1p is required for periodic integration of the Ty1 retrotransposon and targeting of Isw2p to S. cerevisiae tDNAs. Genes Dev. 2005, 19 (8): 955-964. 10.1101/gad.1299105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Gkikopoulos T, Schofield P, Singh V, Pinskaya M, Mellor J, Smolle M, Workman JL, Barton GJ, Owen-Hughes T: A role for Snf2-related nucleosome-spacing enzymes in genome-wide nucleosome organization. Science. 2011, 333 (6050): 1758-1760. 10.1126/science.1206097.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Tirosh I, Sigal N, Barkai N: Divergence of nucleosome positioning between two closely related yeast species: genetic basis and functional consequences. Mol Syst Biol. 2010, 6: 365-

    PubMed Central  PubMed  Google Scholar 

  28. Soutourina J, Bordas-Le Floch V, Gendrel G, Flores A, Ducrot C, Dumay-Odelot H, Soularue P, Navarro F, Cairns BR, Lefebvre O, et al: Rsc4 connects the chromatin remodeler RSC to RNA polymerases. Mol Cell Biol. 2006, 26 (13): 4920-4933. 10.1128/MCB.00415-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Percudani R, Pavesi A, Ottonello S: Transfer RNA gene redundancy and translational selection in Saccharomyces cerevisiae. J Mol Biol. 1997, 268 (2): 322-330. 10.1006/jmbi.1997.0942.

    Article  CAS  PubMed  Google Scholar 

  30. Kutter C, Brown GD, Goncalves A, Wilson MD, Watt S, Brazma A, White RJ, Odom DT: Pol III binding in six mammals shows conservation among amino acid isotypes despite divergence among tRNA genes. Nat Genet. 2011, 43 (10): 948-955. 10.1038/ng.906.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Kassavetis GA, Geiduschek EP: Transcription factor TFIIIB and transcription by RNA polymerase III. Biochem Soc Trans. 2006, 34 (Pt 6): 1082-1087.

    Article  CAS  PubMed  Google Scholar 

  32. Raymond KC, Raymond GJ, Johnson JD: In vivo modulation of yeast tRNA gene expression by 5'-flanking sequences. EMBO J. 1985, 4 (10): 2649-2656.

    PubMed Central  CAS  PubMed  Google Scholar 

  33. Shaw KJ, Olson MV: Effects of altered 5'-flanking sequences on the in vivo expression of a Saccharomyces cerevisiae tRNATyr gene. Mol Cell Biol. 1984, 4 (4): 657-665.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Giuliodori S, Percudani R, Braglia P, Ferrari R, Guffanti E, Ottonello S, Dieci G: A Composite Upstream Sequence Motif Potentiates tRNA Gene Transcription in Yeast. J Mol Biol. 2003, 333: 1-20. 10.1016/j.jmb.2003.08.016.

    Article  CAS  PubMed  Google Scholar 

  35. Zhang G, Lukoszek R, Mueller-Roeber B, Ignatova Z: Different sequence signatures in the upstream regions of plant and animal tRNA genes shape distinct modes of regulation. Nucleic Acids Res. 2011, 39 (8): 3331-3339. 10.1093/nar/gkq1257.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Vinayachandran V, Pusarla RH, Bhargava P: Multiple sequence-directed possibilities provide a pool of nucleosome position choices in different states of activity of a gene. Epigenetics Chromatin. 2009, 2 (1): 4-10.1186/1756-8935-2-4.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Whitehouse I, Rando OJ, Delrow J, Tsukiyama T: Chromatin remodelling at promoters suppresses antisense transcription. Nature. 2007, 450 (7172): 1031-1035. 10.1038/nature06391.

    Article  CAS  PubMed  Google Scholar 

  38. Dieci G, Sentenac A: Facilitated recycling pathway for RNA polymerase III. Cell. 1996, 84 (2): 245-252. 10.1016/S0092-8674(00)80979-4.

    Article  CAS  PubMed  Google Scholar 

  39. Haldar D, Kamakaka RT: tRNA genes as chromatin barriers. Nat Struct Mol Biol. 2006, 13 (3): 192-193. 10.1038/nsmb0306-192.

    Article  CAS  PubMed  Google Scholar 

  40. Bi X, Yu Q, Sandmeier JJ, Zou Y: Formation of boundaries of transcriptionally silent chromatin by nucleosome-excluding structures. Mol Cell Biol. 2004, 24 (5): 2118-2131. 10.1128/MCB.24.5.2118-2131.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Donze D: Extra-transcriptional functions of RNA Polymerase III complexes: TFIIIC as a potential global chromatin bookmark. Gene. 2012, 493 (2): 169-175. 10.1016/j.gene.2011.09.018.

    Article  CAS  PubMed  Google Scholar 

  42. Ng H, Robert F, Young R, Struhl K: Genome-wide location and regulated recruitment of the RSC nucleosome-remodeling complex. Genes & development. 2002, 16 (7): 806-819. 10.1101/gad.978902.

    Article  CAS  Google Scholar 

  43. Schwabish MA, Struhl K: The Swi/Snf complex is important for histone eviction during transcriptional activation and RNA polymerase II elongation in vivo. Mol Cell Biol. 2007, 27 (20): 6987-6995. 10.1128/MCB.00717-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. James P, Whelen S, Hall BD: The RET1 gene of yeast encodes the second-largest subunit of RNA polymerase III. Structural analysis of the wild-type and ret1-1 mutant alleles. J Biol Chem. 1991, 266 (9): 5616-5624.

    CAS  PubMed  Google Scholar 

  45. Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJ: Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005, 309 (5734): 626-630. 10.1126/science.1112178.

    Article  CAS  PubMed  Google Scholar 

  46. Di Gesu V, Lo Bosco G, Pinello L, Yuan GC, Corona DF: A multi-layer method to study genome-scale positions of nucleosomes. Genomics. 2009, 93 (2): 140-145. 10.1016/j.ygeno.2008.09.012.

    Article  CAS  PubMed  Google Scholar 

  47. Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Methods Enzymol. 2006, 411: 134-193.

    Article  CAS  PubMed  Google Scholar 

  48. FastX toolkit: FASTQ/A short-reads pre-processing tools. http://hannonlab.cshl.edu/fastx_toolkit/index.html,

  49. Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK: Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010, 38 (4): 576-589. 10.1016/j.molcel.2010.05.004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Malone BM, Tan F, Bridges SM, Peng Z: Comparison of four ChIP-Seq analytical algorithms using rice endosperm H3K27 trimethylation profiling data. PLoS One. 2011, 6 (9): e25260-10.1371/journal.pone.0025260.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Albert I, Wachi S, Jiang C, Pugh BF: GeneTrack–a genomic data processing and visualization framework. Bioinformatics. 2008, 24 (10): 1305-1306. 10.1093/bioinformatics/btn119.

    Article  CAS  PubMed  Google Scholar 

  54. Wood V, Harris MA, McDowall MD, Rutherford K, Vaughan BW, Staines DM, Aslett M, Lock A, Bahler J, Kersey PJ, Oliver SG: PomBase: a comprehensive online resource for fission yeast. Nucleic Acids Res. 2012, 40: D695-699. 10.1093/nar/gkr853.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. biotoolbox: Tools for querying and analysis of genomic data. [http://code.google.com/p/biotoolbox/]

  56. Fu Q, Ruegger P, Bent E, Chrobak M, Borneman J: PRISE (PRImer SElector): software for designing sequence-selective PCR primers. J Microbiol Methods. 2008, 72 (3): 263-267. 10.1016/j.mimet.2007.12.004.

    Article  CAS  PubMed  Google Scholar 

  57. Pavesi A, Conterio F, Bolchi A, Dieci G, Ottonello S: Identification of new eukaryotic tRNA genes in genomic DNA databases by a multistep weight matrix analysis of transcriptional control regions. Nucleic Acids Res. 1994, 22 (7): 1247-1256. 10.1093/nar/22.7.1247.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to Akshay Bhinge, Chris Benner and Rohit Singh for helpful tips in data analysis; to Gordon Brown for sharing the code for probabilistic reallocation of non-unique reads; to Kevin Struhl and Bradley Cairns for sharing strains. This work was supported by the grants to PB from DBT, Govt. of India and IFCPAR. YK is recipient of Senior Research Fellowship from Council of Scientific and Industrial Research (CSIR).

Accession codes

All datasets from this study have been submitted to NCBI Geo for public access. Microarray data is linked to accession ID-GSE39802 and pol III ChIP-seq data to accession ID-GSE39566.

Author information

Authors and Affiliations

  1. Centre for Cellular and Molecular Biology, Council of Scientific and Industrial Research, Uppal Road, Hyderabad, 500007, India

    Yatendra Kumar & Purnima Bhargava

Authors
  1. Yatendra Kumar
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Purnima Bhargava
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Purnima Bhargava.

Additional information

Competing interests

The authors declare that they have no competing financial interests.

Authors’ contributions

PB conceived the study and wrote the manuscript, PB and YK designed the experiments; YK performed the experiments and analyzed data. Both authors read and approved the final manuscript.

Electronic supplementary material

12864_2013_5111_MOESM1_ESM.pdf

Additional file 1: Figure S1: A prominent NFR at gene region is a characteristic feature of yeast tRNA genes. A) Heat maps show nucleosomes aligned according to TSS of 270 tDNAs, plotted by arranging the genes according to first, the DS nucleosome distance from TTS and then according to the families. One nucleosome can be always found in a narrow window upstream of the TSS (position marked by a bent arrow) whereas downstream nucleosome does not align to the same position according to TSS or TTS. Variations if any, in positions of US/DS nucleosomes and lengths of NFRs at gene region are better visible in this panel. A small number of genes, with DS nucleosomes farthest from TTS, are marked as “Select” by a short bar on right hand side of the panel. B) Heat maps for nucleosomes on 538 out of total 631 tRNA genes in human resting CD4+ cells were extracted from Schones et al. [7] and occupancy scores are given as color gradients. C) Heat maps for nucleosomes on 125 out of total 170 tRNA genes in S. pombe were extracted from Givens et al. [2]. D) Heat map of nucleosome profile after 4 hrs of nutrient starvation. Genes are arranged similar to panel A. Nucleosome dynamics mainly at 3’ end of several genes becomes evident on comparison of the two panels. E) and F) Bin-wise log2 intensity ratios for nucleosome signals on both sides for the tRNA genes according to TSS or TTS are plotted. The data was extracted and analyzed from the Zhang et al. [24]. Averaged data for all the 270 genes or a set of 30 genes (Select, panel A) are plotted together. The plots in both the panels show the presence of a nucleosome on the gene bodies of all 270 genes. An additional downstream nucleosome seen on the select genes suggests the variable position of the DS nucleosome is also not sequence-directed. (PDF 941 KB)

12864_2013_5111_MOESM2_ESM.pdf

Additional file 2: Figure S2: Nucleosome profile over tRNA genes in isw1∆2∆ mutant. A) Normalized log2 [nucleosomal DNA/genomic DNA] intensity ratios for nucleosomes were binned around the TSS and TTS and averaged for all the 270 genes in wild type or isw1∆2∆ cells. Perturbation of the upstream nucleosomal array is evident in case of the mutant. B) Figure shows normalized log2 ratios [nucleosomal DNA/genomic DNA] binned around the TSS of tRNA genes in wild type and isw2∆ strain from Whitehouse et al. [37]. In this dataset, a sharp signal over tRNA gene body is also seen, which most likely is an artifact of array design as it is also obsereved in another dataset by Lee et al. [5] which used the same microarray. C) Venn intersections of the genes occupied by Isw1 or Isw2 [13, 26]. D) Venn intersection of genes occupied by Isw2p and 112 genes showing nucleosome changes near them in isw1∆2∆ mutant cells. E) Comparison of nucleosome profiles on 273 tDNAs in wild type and isw1∆ cells according to available data [27] shows NFR on tDNA is occupied by a nucleosome in the mutant, which could be same as sequence-directed nucleosome on tDNAs (Figure 1E, Additional file 1: Figure S1E-F) [9, 24]. (PDF 531 KB)

12864_2013_5111_MOESM3_ESM.pdf

Additional file 3: Figure S3: Pol III and nucleosome occupancy influence transcription. A) RNA levels normalized for copy number of genes by using genomic DNA as calibrator for Real Time PCR quantification of the cDNA for 33 tDNA families, relative to U4 snRNA. Family name descriptions and primer sequences are given in the Additional file 4: Table S3. B) Data analysis using either unique reads only or multi-reads reallocated to the loci, gives similar results. Normalized tag counts from pol III ChIP-seq experiment were calculated in two ways. All unique reads represent tag coverage binned around TSS and averaged for all genes, using only uniquely aligned reads. All reallocated reads represent tag coverage taking unique reads along with non-unique reads that map up to 20 places in genome and reallocating the non- unique reads probabilistically using a previously reported algorithm [30]. Both the plots show high overlap as tRNA genes although repetitive in sequence, are very small in comparison to the ChIP fragment size used in typical ChIP experiments. C) Normalized and averaged tag coverage profiles of all tDNAs and pseudo-tRNA (Additional file 4: Table S2) for RPC128-FLAG IP and mock samples. D) Family wise total pol III tag counts for all the genes that are covered in panel A and Figure 4A. Pol III tag count in a region −100 to +100 bp of every gene was averaged. Then for each primer set, this average pol III count is summed up for all genes that are covered by the primer set. E) Pol III fold-enrichment over mock from ChIP-Seq data for selected genes showing very high to very low pol III levels. F) The data in panel E is validated by ChIP-Real time PCR quantification for pol III on the same genes. Pol III occupancy data for both IP and mock samples are compared with the same on TELVIR region. (PDF 454 KB)

12864_2013_5111_MOESM4_ESM.xls

Additional file 4: Tables listing strains and primers used in this study. Table S1. List of all primers and their sequences, used for tRNA expression analysis by Real Time PCR quantification method. Primers unique for five tRNAs are highlighted. Table S2. Details of tDNA like elements (pseudo-tDNAs) in yeast genome. Data is extracted from “pol3- scan” data [57]. Table S3. Details of yeast strains used in this study. Table S4. Details of primers used for quantifications in ChIP assays. (XLS 26 KB)

12864_2013_5111_MOESM5_ESM.pdf

Additional file 5: Figure S4: Quality control of microarray samples and data. (A) Quality check for the mono-nucleosomal DNA preparation. Lane M shows 100 bp DNA ladder used as molecular size marker. Purity of DNA before (lane 1) and after (lane 2) MNase digestion is shown. B) Genomic DNA control for nucleosome mappings after MNase digestion was checked. M shows the 100 bp DNA ladder used as size marker. C) Box plots showing near identical distribution of normalized Log2 ratios among three biological replicates. Pearson correlations between all possible replicate-pairs were 0.97 to 0.99, which confirmed the reproducibility of biological replicates. (PDF 218 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Kumar, Y., Bhargava, P. A unique nucleosome arrangement, maintained actively by chromatin remodelers facilitates transcription of yeast tRNA genes. BMC Genomics 14, 402 (2013). https://doi.org/10.1186/1471-2164-14-402

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-402

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us