A unique nucleosome arrangement, maintained actively by chromatin remodelers facilitates transcription of yeast tRNA genes
© Kumar and Bhargava; licensee BioMed Central Ltd. 2013
Received: 14 February 2013
Accepted: 6 June 2013
Published: 17 June 2013
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.
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.
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.
KeywordsChromatin Isw1 Isw2 Nucleosomes RNA polymerase III RSC Transcription Yeast tRNA
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 . 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. 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 . 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  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.
Chromatin structure of tDNAs in active and repressed states
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 . 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 , 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  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
RSC-pol III communication maintains the active state chromatin structure
tDNA chromatin is required for optimum transcription
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 ; 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  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  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 . 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
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.
Sequence-directed nucleosome positions can be further modulated by DNA binding proteins and chromatin remodelers . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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  as nucleosome-free regions have been shown to serve as boundary elements . Replication forks are reported to stall with high frequency near tRNA genes , 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.
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.
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 .
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  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  was applied for detection of nucleosomes through MLM package  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 . 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 , 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 . 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 . 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 , 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 . Reproducibility of replicates was checked with USeq package  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  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  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  (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  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 . 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  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 . Data from three replicates were quantile normalized and log transformed using biotoolbox . 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 . 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  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.
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).
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Jansen A, Verstrepen KJ: Nucleosome positioning in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2011, 75 (2): 301-320. 10.1128/MMBR.00046-10.PubMed CentralView 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 (8): 900-905. 10.1038/ng1400.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Kassavetis GA, Geiduschek EP: Transcription factor TFIIIB and transcription by RNA polymerase III. Biochem Soc Trans. 2006, 34 (Pt 6): 1082-1087.View ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitehouse I, Rando OJ, Delrow J, Tsukiyama T: Chromatin remodelling at promoters suppresses antisense transcription. Nature. 2007, 450 (7172): 1031-1035. 10.1038/nature06391.View ArticlePubMedGoogle Scholar
- Dieci G, Sentenac A: Facilitated recycling pathway for RNA polymerase III. Cell. 1996, 84 (2): 245-252. 10.1016/S0092-8674(00)80979-4.View ArticlePubMedGoogle Scholar
- Haldar D, Kamakaka RT: tRNA genes as chromatin barriers. Nat Struct Mol Biol. 2006, 13 (3): 192-193. 10.1038/nsmb0306-192.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- FastX toolkit: FASTQ/A short-reads pre-processing tools. http://hannonlab.cshl.edu/fastx_toolkit/index.html,
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- biotoolbox: Tools for querying and analysis of genomic data. [http://code.google.com/p/biotoolbox/]
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar