- Open Access
Genome-wide analysis of functional sirtuin chromatin targets in yeast
Genome Biology volume 14, Article number: R48 (2013)
The sirtuins are a conserved family of NAD+-dependent histone/protein deacetylases that regulate numerous cellular processes, including heterochromatin formation and transcription. Multiple sirtuins are encoded by each eukaryotic genome, raising the possibility of cooperativity or functional overlap. The scope and variety of chromatin binding sites of the sirtuins in any specific organism remain unclear.
Here we utilize the ChIP-seq technique to identify and functionally characterize the genome-wide targets of the sirtuins, Sir2, Hst1 to Hst4, and the DNA binding partner of Hst1, Sum 1, in Saccharomyces cerevisiae. Unexpectedly, Sir2, Hst1 and Sum1, but not the other sirtuins, exhibit co-enrichment at several classes of chromatin targets. These include telomeric repeat clusters, tRNA genes, and surprisingly, the open reading frames (ORFs) of multiple highly expressed RNA polymerase II-transcribed genes that function in processes such as fermentation, glycolysis, and translation. Repression of these target genes during the diauxic shift is specifically dependent on Sir2/Hst1/Sum1 binding to the ORF and sufficiently high intracellular NAD+ concentrations. Sir2 recruitment to the ORFs is independent of the canonical SIR complex and surprisingly requires Sum1. The shared Sir2/Hst1/Sum1 targets also significantly overlap with condensin and cohesin binding sites, where Sir2, Hst1, and Sum1 were found to be important for condensin and cohesin deposition, suggesting a possible mechanistic link between metabolism and chromatin architecture during the diauxic shift.
This study demonstrates the existence of overlap in sirtuin function, and advances our understanding of conserved sirtuin-regulated functions, including the regulation of glycolytic gene expression and condensin loading.
The sirtuins are a highly conserved family of NAD+-dependent protein deacetylases that regulate a wide range of cellular processes impacted during aging and in age-associated diseases such as type 2 diabetes and cancer (reviewed in ). They utilize NAD+ as a co-substrate during the deacetylation reaction, such that one molecule of NAD+ is converted into nicotinamide and 2'O-acetyl ADP-ribose for every lysine that is deacetylated [2, 3]. As a result, sufficiently high NAD+ concentrations are required to properly regulate cellular processes in which the sirtuins participate [4, 5]. This link to NAD+ gives sirtuins an inherent ability to 'sense' the intracellular energy status, and regulate target proteins via lysine deacetylation.
The genomes of eukaryotic organisms usually encode multiple sirtuin proteins. For example, there are seven sirtuins in mammals known as SIRT1 through SIRT7 , while the budding yeast Saccharomyces cerevisiae encodes five, known as Sir2, and Hst1 through Hst4 (homologs of Sir two) . Sir2 is the founding family member , and was initially characterized as a factor required for transcriptional silencing at the HML and HMR silent mating-type loci, telomeres, and the ribosomal DNA (rDNA) tandem array, each of which have characteristics of heterochromatin in more complex eukaryotes (reviewed in ). Eventually, Sir2 was found to be a histone deacetylase [9, 10], a seminal discovery that instantly provided a mechanistic role for Sir2 in the formation of heterochromatin. Equally exciting were the implications for aging because Sir2 was also characterized as a limiting factor for replicative life span , which is defined as the number of times a yeast mother cell divides before senescing. Deletion of SIR2 shortens replicative life span, while increased SIR2 gene dosage extends both mean and maximum replicative life span . Similarly, SIRT6 knock out mice prematurely age , and male SIRT6 transgenics are long-lived , suggesting that longevity and/or health span support by sirtuins could be one of their conserved features.
Many histone-modifying enzymes are catalytic subunits of large multi-protein complexes, and the nuclear sirtuins appear to follow this trend. At the HM loci and telomeres, yeast Sir2 is associated with the Sir3 and Sir4 proteins in a complex known as SIR [14–16]. Deleting either of the SIR complex subunits results in a loss of transcriptional silencing [17, 18]. At the rDNA locus, Sir2 associates with Net1 and Cdc14 in the nucleolar silencing complex known as RENT, which silences RNA polymerase II transcription from the intergenic spacers [19, 20]. The Sir2 paralog, Hst1, forms a complex with Sum1 and Rfm1 (the Sum1 complex), which represses specific genes through localized histone deacetylation at promoters [21–23]. SIRT1 is considered the mammalian Sir2 ortholog, and like Sir2, is a histone deacetylase that can function in heterochromatin formation . SIRT1 also regulates the expression of numerous genes and deacetylates numerous non-histone protein targets, but it has only recently been considered part of a larger multi-protein co-repressor complex . Thus, sirtuins are generally recruited to multiple targets (chromatin or non-histone) through interactions with partner proteins that provide an added level of binding specificity. Given such diversity, complete inventories of conserved sirtuin targets, both chromatin and non-chromatin, are needed in order to understand how these factors are truly integrated with various cellular processes via NAD+.
In this study we have focused on the identification and characterization of novel chromatin targets of sirtuins in yeast, with the goal of uncovering additional cellular pathways that are impacted by alterations in cellular NAD+, and determining the extent of functional overlap between the various sirtuins. Conserved targets of this nature are also more likely to be mediators of longevity. We have utilized chromatin immunoprecipitation (ChIP), followed by next generation DNA sequencing (ChIP-seq), to obtain high-resolution chromatin association maps for each yeast sirtuin. Work presented in this study functionally characterizes the overlapping roles for Sir2, Hst1, and Sum1 in telomere maintenance, as well as NAD+-dependent repression of specific genes that are downregulated during the diauxic shift, a time in which yeast shift their metabolism from aerobic fermentation to respiration. These include genes involved in glycolysis, fermentation, and translation. Sir2, Hst1, and Sum1 were also required for the efficient recruitment of condensin and cohesin onto multiple shared target sites, including tRNAs and other genes downregulated during the diauxic shift.
From an earlier expression profiling study, we identified thiamine biosynthesis genes as being upregulated when yeast cells were treated with the sirtuin inhibitor nicotinamide, or when intracellular NAD+ levels were reduced by deleting the NAD+ salvage pathway gene NPT1 . Hundreds of other genes were also upregulated, including known Hst1/Sum1 targets such as middle sporulation and NAD+ biosynthesis genes. These results raised the question of whether each sirtuin specifically regulates a distinct set of gene targets? We therefore set out to obtain a genome-wide picture of binding sites for each sirtuin using ChIP-seq. The sirtuins were carboxy-terminally tagged either with 13 copies of the Myc epitope (Sir2, Hst1, and Hst2) or one copy of the TAP-tag (Hst3 and Hst4). Sum1 was also tagged with 13xMyc to test if it consistently tracked with its binding partner Hst1. Cells growing exponentially in rich YPD medium were then subjected to ChIP with anti-Myc or anti-TAP antibodies, and the recovered DNA sequenced. Certain genomic regions were previously shown to be over- or under-represented when genomic DNA was sequenced by this method . Therefore, to control for representation in the ChIP libraries, crosslinked DNA that went into the immunoprecipitation (IP) reactions (the input) was also sequenced.
Sirtuin binding at the silent mating-type loci and the rDNA tandem array
We first analyzed Sir2 at the silent mating-type loci HML and HMR to confirm the ChIP-seq procedure was effective. Sir2 was clearly enriched at the HML-E and -I silencers as expected (Figure 1a). Sequence reads derived from repetitive regions are discarded by the mapping program, so there is a lack of binding information for the α1 and α2 genes within HML because they are duplicated at MATα. However, Sir2 was still noticeably enriched between the two silencers (Figure 1a), which supported the results of an earlier study suggesting unidirectional spreading from each silencer . At HMR we observed a strong peak of Sir2 at the HMR-E silencer that extended rightward across the body of HMR as expected, but there was no corresponding peak at the HMR-I silencer (Figure 1b). Instead, a very strong peak of Sir2 was located at an adjacent tRNAThr gene, which had already been well characterized as a silencing boundary element . Earlier ChIP results showed that Sir2 does not associate with HMR-I when the E silencer is defective, suggesting that Sir2 is not independently recruited to the I silencer, but rather spreads all the way across HMR until it reaches the boundary . Our data suggest an extension of this model in which the tRNAThr gene plays a dual role as a boundary element and as a protosilencer where Sir2 is independently recruited. A protosilencer is defined as a cis-acting sequence that cannot establish silencing on its own, but instead cooperates in establishment and maintenance of the silent chromatin . Support for such a protosilencer model comes from a recent study where the tRNAThr boundary element was found to impose the cell cycle progression requirement for establishing silencing at HMR .
There was no visible enrichment of Hst2, Hst3 or Hst4 at the silencers or across HML or HMR, but surprisingly, Hst1 and Sum1 were both highly enriched and co-localized with Sir2 at the tRNAThr boundary element flanking HMR (Figure S1 in Additional file 1), suggesting they could also potentially be involved in the cell cycle progression requirement in silencing linked to this cis-acting element. Given the unexpected co-localization of Sir2, Hst1, and Sum1 at the tRNAThr gene, we next asked whether Sum1 and Hst1 were properly associating with several of their known targets, including the promoters of specific NAD+ biosynthesis and middle sporulation genes [21, 33]. As shown in Figure 1c, Sum1 was clearly bound to such promoters (BNA2 and SPS4) with a single strong peak and without any indication of spreading, consistent with previous low-resolution ChIP assays . Hst1 was also highly enriched at the same sites, but its occupancy was extended a short distance (one or two nucleosomes) in both directions (Figure 1c). Importantly, Sir2 does not repress these genes , and we did not detect any binding signal in the ChIP-seq profiles (Figure 1c). Taken together, these results indicated that the ChIP-seq data were highly specific, and the co-association of Sir2, Hst1, and Sum1 at the tRNAThr gene and other genomic sites (see below) is likely meaningful.
Sir2-dependent transcriptional silencing also occurs at the rDNA tandem array [36, 37], which consists of approximately 150 to 200 identical rDNA gene copies on the right arm of chromosome XII. Only the leftmost and rightmost repeats are annotated, so ChIP-seq reads that map to these two repeats are a compiled average for the entire array (Figure 1d). The mapping program parameters were adjusted in this case to allow for repeated sequences, so the read counts for Sir2 peaks in the rDNA approached 200,000, compared to approximately 2,000 for HML and HMR. Consistent with previous silencing and ChIP results [38, 39], Sir2 strongly associated with the rDNA intergenic spacers (IGS1 and IGS2), and its extended enrichment to the left of IGS2 occurs in the same direction as Pol I transcription (Figure 1d). The strongest Sir2 peak in IGS1 corresponded exactly to a replication fork block site (chromosome XII coordinates: 460530-460570), where Fob1 recruits the RENT and cohibin complexes . While there may be some low-level association of Hst1, Sum1, Hst2, Hst3, and Hst4 across the rDNA, it was minimal compared to the overwhelming Sir2 signal (Figure S1 in Additional file), which is presumably derived from the RENT complex.
Telomere length maintenance by Sir2, Hst1, and Sum1
The third known type of silencing in yeast occurs at telomeres. The SIR complex is recruited to telomeres through physical interaction with the Rap1 protein, which directly binds to telomeric TG1-3 repeats [41, 42]. As with the HM loci, sequence read coverage of many telomeric regions was sparse. When there was sufficient coverage, Sir2 association with the terminal TG1-3 repeats was observed (Figure 2a), but its relative distribution among different telomeres was quite variable (Table S1 in Additional file 1). Surprisingly, strong peaks of Hst1 and Sum1 binding were observed directly overlapping with Sir2 at the TG1-3 repeats, regardless of whether they were terminal or internal clusters (Figure 2a,b; Table S1 in Additional file 1). From binding sites identified using BayesPeak , we calculated the enrichment of Hst1, Sum1, or Sir2 for telomeric repeats was at least 30-fold greater than expected by chance (Table S2 in Additional file 1). Sir2 enrichment at terminal and internal TG1-3 clusters typically extended inward up to approximately 1 kb, whereas Hst1 and Sum1 were more restricted to the actual repeats (Figure 2a,b).
To test whether Hst1 and Sum1 were recruited to telomeric repeats via Rap1, their enrichment was assayed by quantitative ChIP in wild type (WT) and rap1-17 mutant strains . The rap1-17 allele truncates the carboxy-terminal SIR interaction domain of Rap1, resulting in mis-localization of the SIR complex and severe silencing defects at the HM loci and telomeres [42, 45]. The repetitive nature of telomeric TG1-3 sequences prevented direct real-time PCR amplification of these regions, so we instead assayed 70 bp adjacent to the cluster on the right arm of chromosome VI, or 7,000 bp away from the repeats as a control. As expected, strong telomeric enrichment of Sir2 at the 70 bp position was eliminated in the rap-17 mutant (Figure 2c, left panel). There was also a low level of Hst1 and Sum1 enrichment 70 bp from the TG repeats in the WT strain, confirming the poor spreading capability of these two proteins, relative to Sir2. Importantly, even this low level Hst1 and Sum1 occupancy was significantly reduced in the rap1-17 mutant (Figure 2c, right panel), suggesting that Rap1 significantly contributes to their recruitment. However, we cannot rule out the possibility that occupancy was indirectly reduced in the mutant because of the general defect in telomeric heterochromatin formation.
Given the precise association of Sir2, Hst1, and Sum1 with terminal TG repeats, we next tested for functional consequences of deleting the sirtuins or SUM1 on telomere length. As shown in Figure 2d, sir2Δ or hst1Δ mutants had little effect on telomere lengths compared to WT strains, but the sir2Δ hst1Δ double mutant clearly had shorter telomeres on average, indicating significant redundancy between the two paralogs. The sum1Δ mutant also had short telomeres similar to the double mutant length. Compared to Sir2, Hst1, and Sum1, the occupancy of Hst2, Hst3, and Hst4 at the TG1-3 repeats was absent or extremely low (Figures 2a,b), suggesting the slightly longer telomeres in hst2Δ and hst3Δ mutants were likely an indirect effect (Figure 2d). Because Sir2, Hst1, and Sum1 appeared to have more specific co-enriched binding sites and functional redundancy than initially anticipated, we focused subsequent analysis on these three proteins.
Sir2, Hst1, and Sum1 associate with and regulate highly transcribed genes that are downregulated during the diauxic shift
Comparisons of BayesPeak-identified Sir2, Hst1 or Sum1 enrichment peaks to gene annotations revealed that the largest percentage of binding sites for each protein was surprisingly found within ORFs (Figure 3a). To determine if there was any relationship between ORF enrichment and expression levels for the associated genes, composite plots were generated for Sir2, Hst1, Sum1, or Hst2 ChIP-Seq read numbers across all yeast genes, with the expression levels divided into five quintiles (Figure 3b). Sir2, Hst1, and Sum1 enrichment, but not Hst2, was strongest toward the 3' end of ORFs with the highest expression levels (5th quintile). There was also significant overlap between the ORFs bound by each factor (Figure 3c), which was easily visualized in the genome browser for highly expressed genes such as PDC1, ENO2, and CDC19 (Figure 3d). Traditional ChIP assays were performed with several target genes to confirm the ORF enrichment was not a sequencing artifact (Figure 3e; Figure S2a in Additional file 1). Hst2 did not have this pattern of binding in the composite plots or in the genome browser, and was not ORF-enriched in standard ChIP assays compared to an untagged control (Figure 3e), consistent with its reported cytoplasmic localization pattern .
Many of the most highly expressed genes during log phase tend to be strongly repressed when cells enter the diauxic shift. We used publicly available expression profiling data , to ask whether genes with Sir2, Hst1, or Sum1 binding on their ORFs were differentially regulated during the diauxic shift (Figure S2b in Additional file 1). For each individual factor there was a trend for association with downregulated genes, which reached strong statistical significance for Hst1. Coupled with this trend, there was also a highly significant depletion for upregulated genes (Figure S2b in Additional file 1,). Furthermore, genes with overlapping Sir2, Hst1, and Sum1 ORF binding were enriched for Gene Ontology terms related to glycolysis, glucose fermentation, translation, and cell wall biosynthesis (Table S3 in Additional file 1), all processes that are highly active during log phase and then repressed during the diauxic shift . Five of the ORF-targeted genes encode enzymes involved in glucose fermentation (Figure 4a, highlighted in red), so we hypothesized that expression of such genes would be dysregulated when SIR2, HST1, or SUM1 were deleted. No significant differences in PDC1 or ENO2 expression were observed in the deletion mutants during log phase (Figure 4b; Figure S3a in Additional file 1), but they were not properly repressed when entering the diauxic shift (Figure 4c; Figure S3b in Additional file 1). PDC1 and ENO2 in the WT strain were strongly repressed within 6 hours after log phase, while repression in the mutants was significantly delayed. This was especially true in the sir2Δ hst1Δ double mutant, which was unable to fully repress PDC1 even 10 hours after log phase (Figure 4d), suggesting some redundancy between Sir2 and Hst1, as was also observed in telomere length regulation (Figure 2d).
During the diauxic shift, cellular metabolism switches from fermentation to respiration, which causes an increase in the NAD+/NADH ratio. We observed this increased ratio in the WT strain within 3 hours post log phase (Figure 5a). Since the higher NAD+/NADH ratio induced by calorie restriction has been reported to promote Sir2 function , we next tested whether simply reducing the intracellular NAD+ concentration by deleting NPT1 would also delay PDC1 repression, which it did (Figure 5b). Similar results were observed for ENO2 and CDC19 (Figure S3c in Additional file 1). The time course of glucose consumption from the YPD medium was almost identical between the WT strain and the sir2Δ, hst1Δ, or sum1Δ mutants (Figure 5c), but it was still formally possible that the mutants were not properly repressing PDC1 simply because their entry into the diauxic shift was delayed. If this were the case, then other genes down-regulated at the diauxic shift whose ORFs were not bound by Sir2, Hst1, or Sum1 should also show a delay in repression in the mutant strains. PDC5 encodes a minor isoform of pyruvate decarboxylase (Figure 4a), but does not have Sir2, Hst1, or Sum1 enriched on the ORF (Figure 5d). Importantly, PDC5 repression at the diauxic shift occurred normally in the sir2Δ, hst1Δ, or npt1Δ mutants (Figure 5e). To make sure this was not a peculiarity of PDC5, we tested three additional genes from this class (CLB1, RIF1, and RPB3) in the sir2Δ and npt1Δ mutants, and again did not observe any defect in repression across the time course (Figure S4 in Additional file 1). These results suggest a model where Sir2 and Hst1 are associated with specific targeted ORFs in a poised state during log phase, and then become functional for transcriptional repression at onset of the diauxic shift, perhaps aided by the increased NAD+/NADH ratio. This would result in the observed early and rapid gene repression, probably before the overall reduction in the general transcription machinery occurs.
We initially hypothesized that Sir2 and Hst1 function in repression of their ORF-targeted genes by deacetylating histones across the ORF when cells enter the diauxic shift. However, using quantitative ChIP, there was no consistent increase in various histone H3 or H4 acetylation marks on the PDC1 or ENO2 genes in sir2Δ, hst1Δ, sum1Δ, or npt1Δ mutants during log phase or the diauxic shift (data not shown). Thinking that there could be some redundancy involved, sir2Δ sum1Δ and sir2Δ hst1Δ double mutants were also tested, but H3 or H4 acetylation was still not elevated compared to WT (Figure S5 in Additional file 1). Even so, the requirement for high NAD+ concentration in Figure 5b strongly suggested that Sir2 and/or Hst1 catalytic activity was required for the repression. The effect of an H364Y mutation in SIR2 that eliminates deacetylase activity was tested and found to block PDC1 repression at onset of the diauxic shift (Figure 5f), similar to the sir2 deletion. Therefore, the mechanism of repression remains unknown, but current evidence points toward either non-traditional lysine modifications or non-histone deacetylation targets such as RNA polymerase or elongation factors being important (see Discussion).
Evidence for novel Sir2 and Hst1 complexes on targeted ORFs
The apparent inactivity of Sir2 and Hst1 on histones at the targeted ORFs suggested they could be functioning independently of their canonical HDAC complexes at these sites. To check this possibility, we first tested whether Sir2 or Hst1 enrichment on the PDC1 ORF was altered in a sum1Δ mutant. Sir2 enrichment was surprisingly lost, while Hst1 was increased (Figure 6a), strongly suggesting that Sum1 is involved in recruiting Sir2, but not Hst1. We next tested whether Sir2 or Hst1 enrichment was altered in sir3Δ or sir4Δ mutants. As shown in Figure 6b, Sir2 enrichment was maintained in both mutants, and actually increased in the absence of SIR4, reminiscent of the increased nucleolar Sir2 localization and rDNA silencing in a sir4Δ mutant caused by redistribution of the telomeric Sir2 pool [50, 51]. Similarly, increased Hst1 enrichment in the sum1Δ mutant could be due to release of Hst1 from its traditional promoter targets. Since Sum1 appeared to be recruiting Sir2, we thought Hst1 may behave oppositely as well, and be recruited via interactions with Sir4. However, Hst1 enrichment at PDC1 was unaffected in the sir4Δ mutant (Figure 6c). Taken together, these ChIP results suggest a model whereby Sir2 is recruited to its ORF targets via Sum1, while Hst1 is recruited through an unknown bridging factor (Figure 6d). It is important to note that independent mechanisms of Sir2 and Hst1 recruitment are consistent with their observed functional redundancy in gene repression during the diauxic shift (Figure 4).
Sir2, Hst1, and Sum1 promote condensin and cohesin recruitment to multiple chromatin loci, including tRNA and other diauxic shift-repressed genes
We earlier noted significant enrichment of Sir2, Hst1, and Sum1 at the tRNAthr boundary element adjacent to HMR-I (Figure 1b). Subsequent binding site identification on Pol III-transcribed genes using BayesPeak revealed fold enrichment values of between 110- and 169-fold over random expectation (Table S4 in Additional file 1). There was also significant overlap between the Sir2, Hst1, and Sum1 associated genes (Figure S6 in Additional file 1). Example Integrative Genomics Viewer (IGV) screenshots for tQ(UUG)H, tE(UUC)E1, and SNR30 binding are shown in Figure 7a, revealing that these small genes are often the center of a broader peak for all three proteins. Binding to tQ(UUG)H and SNR30 was confirmed by standard ChIP assays (Figure 7b).
tRNA genes and the associated TFIIIC binding sites are well established as being chromatin insulator and boundary elements in yeast and human cells [29, 52]. In yeast, the condensin and cohesin complexes have previously been demonstrated through ChIP-chip experiments to globally associate with tRNA genes [53, 54]. We were curious if there was significant overlap of Sir2, Hst1, or Sum1 sites with the published cohesin or condensin binding sites. Indeed, the peaks for all three tagged proteins showed significant overlap with condensin peaks (approximately 140- to 200-fold over random) and cohesin peaks (approximately 24- to 29-fold over random) (Table S5 in Additional file 1). Sir2 and Hst1 were previously shown to function in cohesin recruitment to the silenced rDNA and HMR loci in yeast [55, 56], but such a relationship with condensin had not been reported. To test whether Sir2, Hst1, and Sum1 are involved in the loading of condensin and/or cohesin onto tRNA genes, the condensin subunits Brn1 and Smc4, and the cohesin subunit Mcd1 were carboxy-terminally tagged with the 13xMyc epitope in WT, sir2Δ, hst1Δ, and sum1Δ strains. Steady state protein expression was equivalent in all the strains, as measured by western blotting with the 9E10 α-Myc antibody (Figure S7a in Additional file 1). ChIP assays for the tagged subunits revealed a strong dependence for Sir2, Hst1, and Sum1 in loading both condensin and cohesin onto tQ(UUG)H and tE(UUC)E1 (Figure 7c).
The links between Sir2, Hst1, and Sum1 with condensin and cohesin were not limited to Pol III-transcribed genes. Pol II-transcribed genes such as PDC1 with Sir2, Hst1, and Sum1 coating the ORF often had a condensin peak located near the 3' end of the gene (Figure 8a). This was intriguing given the trend for Sir2, Hst1, and Sum1 association toward the 3' ends of genes (Figure 3b). Quantitative ChIP demonstrated that condensin (Brn1 and Smc4) and cohesin (Mcd1) were indeed highly enriched in the intergenic region 3' of PDC1 and ENO2 (Figure 8b). But interestingly, significant levels of association were also observed across the ORFs. Regardless of the location tested, or the level of enrichment, the sir2Δ, hst1Δ, and sum1Δ mutations impaired the condensin/cohesin recruitment (Figure 8b), which was also observed for targeted ribosomal protein genes (Figure S7b in Additional file 1). We were unable to detect direct physical interactions between Sir2 and tagged condensin subunits through co-IP assays (data not shown), suggesting the recruitment may not be mediated by a direct physical interaction, but perhaps related to Sir2 catalytic activity. Consistent with that idea, the npt1Δ mutant with reduced NAD+ caused a partial defect in cohesin (Mcd1) or condensin (Smc4) subunit association at the rDNA (Figure 8c). Furthermore, loss of Npt1 reduced enrichment of a cohesin loading factor subunit (Scc2), which may also promote functional condensin association with chromosomes . Since condensin and cohesin mediate long-range chromatin interactions, this suggests that sirtuins could potentially make upstream contributions to this type of chromatin organization. Potential implications are addressed in the discussion.
Overlapping functions of Sir2 and Hst1 mediated by Sum1
One of the unexpected findings from this study was the large number of locations where Sir2, Hst1, and Sum1 were co-enriched. Sir2 and Hst1 are paralogs that have acquired differential functions, but are still similar enough to substitute for one another under specific circumstances . For example, a dominant SUM1-1 mutation suppresses the HM silencing defect of strains deleted for SIR genes by directing Hst1 to the silencers [57, 58]. HST1 overexpression can also suppress the HM silencing defects of a sir2Δ mutant , and Sir2 can partially substitute for Hst1 in an hst1Δ background to suppress middle sporulation genes during vegetative growth . It should be noted that Sum1 has repressive activity at the HMR-E silencer and some meiosis genes that is independent of Hst1 [33, 59]. Perhaps Sum1 interacts with Sir2 at such locations when Hst1 is missing, which could even be mediated by another adaptor protein, similar to the role that Rfm1 plays in facilitating Hst1 interactions with Sum1 . Our unexpected finding that Sum1 is required for Sir2 enrichment at the PDC1 ORF supports this hypothesis, and reveals even more diversity in the mechanism of chromatin targeting for Sir2 than was anticipated. However, the Sir2 recruited to ORFs via Sum1 is likely a minor subset of the overall Sir2 population, as indicated by the relatively lower Sir2-myc ChIP-seq read counts at ORFs compared to telomeres, the rDNA, and HM loci. Furthermore, Sir2 interaction with Sum1 in co-IP experiments has only been observed when HST1 is deleted , implying that most Sir2 is not associated with Sum1. Increased Sir2 enrichment in the sir4Δ mutant is likely due to redistribution from telomeres, and suggests the targeted ORFs are in competition with telomeres and rDNA for limiting Sir2. The mechanism of Sum1 and Hst1 recruitment to specific ORFs remains unknown, and merits future investigation.
We were initially surprised to detect strong Hst1 and Sum1 binding to the telomeric TG1-3 repeats, but genome-wide Southern blotting screens for deletion mutants with altered telomere lengths independently identified sum1Δ as having short telomeres [60, 61]. Our results confirm the role of Sum1 in telomere maintenance, but also implicate Sir2 and Hst1, which appear to compensate for each other when one or the other is deleted. It had been speculated that the effect of the Sum1 complex on telomere length was likely an indirect effect on expression of upstream regulators of telomere length . However, the precise association of Sum1 and Hst1 with the TG repeats strongly suggests that the SIR and Sum1 complexes have more direct functions in protecting the telomere ends. Rap1 recruits the SIR complex to telomere repeats , and we have now shown that Hst1 and Sum1 telomeric recruitment also depends on Rap1. Interestingly, there is a positive genetic interaction (P = 1.7e-06) between RAP1-damp and sum1Δ alleles in the DRYGIN database , supporting the idea of a functional relationship.
tRNAs, sirtuins, and repression of adjacent Pol II-transcribed genes
We previously reported that Hst1 and Sir2 both directly repress basal expression of the thiamine biosynthesis gene THI4 by binding to and deacetylating histone H4 at a region approximately 700 to 800 bp upstream of the transcription start site . That region happens to overlap with the tG(UCC)G gene, which, like many other Pol III-transcribed genes in the dataset, is co-enriched for Sir2, Hst1, and Sum1 binding (Table S4 in Additional file 1, and data not shown). We now hypothesize the tG(UCC)G gene is a cis-acting repressive element in this context, and that other tRNA genes with Sir2, Hst1, and Sum1 co-binding may repress adjacent Pol II-transcribed genes under specific growth conditions. Not all tRNA genes showed strong Sir2/Hst1/Sum1 co-enrichment, so any repressive effect would be limited to a subset of tRNA genes, including those that are within sufficient distance to a susceptible promoter, like THI4. Consistent with this model, an earlier genomic analysis found that Pol II-transcribed genes with a tRNA gene located within approximately 400 bp of the promoter were expressed approximately 3.5-fold less than genes without an adjacent tRNA gene . A more recent study also showed that tRNA genes are often associated with Sir3 and Sir4, and those with strong DNA replication pausing activity can enhance silencing when positioned adjacent to an ADE2 reporter gene flanked by the HMR-E and -I silencers .
The Engelke lab previously described a form of tRNA gene-mediated (TGM) silencing where the SUP4 or SUP53 tRNA genes repress expression of an adjacent HIS3 reporter gene . This silencing is dependent on the clustering of tRNA genes at the nucleolar periphery, but does not appear related to HM, telomeric, or rDNA silencing, because deleting the SIR genes has no effect [54, 66]. Instead, TGM silencing is dependent on the association of condensin and cohesin with tRNA genes [54, 67]. It is therefore intriguing that deletion of SIR2, HST1, or SUM1 causes a reduction in condensin and cohesin association with targeted tRNA genes (Figure 7c), but sir2Δ does not affect TGM silencing when tested with the reporter gene system . The enrichment of Sir2, Hst1, and Sum1 is relatively weak at SUP4 and SUP53 compared to many other tRNA genes (data not shown), which could partially explain why their absence has little impact on TGM silencing, which is typically measured in the context of a plasmid. Alternatively, there could be sufficient redundancy between Hst1, Sum1, and Sir2 in condensin and cohesin recruitment to tRNA genes.
The diauxic shift, glycolysis and sirtuins
Saccharomyces cerevisiae is a facultative anaerobe that ferments glucose even in the presence of oxygen. During log phase, genes involved in glycolysis, fermentation and growth-related processes are highly expressed, whereas genes involved in the tricarboxylic acid cycle and mitochondrial respiration are repressed . As glucose is depleted, the cells undergo a change in gene expression and metabolism called the diauxic shift, during which the cells switch from glycolysis to ethanol catabolism and mitochondrial respiration, and prepare for survival in stationary phase (G0). Gene upregulation during the diauxic shift is mediated by transcriptional activators such as Msn2/Msn4 and Gis1, which turn on specific genes in response to the reduction in nutrients . Specific mechanisms of transcriptional repression during the shift are less well understood, and are generally believed to be caused by downregulation and inactivation of the polymerase machinery , rather than by specific repressors. The results from our study indicate that Sir2, Hst1, and Sum1 directly contribute to the repression of genes to which they are bound across the ORF. The repression does not appear to involve traditional H3 or H4 deacetylation like that observed at the silenced domains (Figure S5 in Additional file 1), although we cannot rule out other less studied acetylation sites or alternative lysine modifications being important . One attractive model is that the polymerase machinery or transcription elongation factors associated with the ORFs during transcription are inactivated by direct deacetylation during the diauxic shift. While it remains unclear if such factors are acetylated in yeast, several conserved RNA polymerase and elongation factor subunits have been identified as acetylated proteins in human cells . In such a model, deacetylation by Sir2 and Hst1 could be partially triggered by the higher NAD+/NADH ratio that occurs at the onset of the diauxic shift (Figure 5a), with additional nutritional signaling inputs also likely involved. An alternative, and not mutually exclusive, model is that a repressive higher-ordered chromatin structure is established at these ORFs during the diauxic shift via the Sir2- and Hst1-mediated recruitment of condensin and cohesin (see below).
Earlier ChIP-chip studies on Sir2, Sir3, or Sir4 binding sites identified several euchromatic targets, but either did not test whether such genes were regulated by the SIR complex, or were unable to detect any expression changes in the SIR mutants [71–73]. We also did not observe changes in target gene expression during log phase growth (Figure 4b), but instead uncovered the role in repression of targeted genes downregulated at the diauxic shift, including glycolytic genes. This regulation of glycolytic genes by Sir2 and Hst1 is reminiscent of Sirt6 regulation of glycolytic genes in the mouse , where Sirt6 functions as a co-repressor for Hif1α, a critical regulator of nutrient stress. Sirt6 knock out mice upregulate genes involved in glycolysis and glucose import, and have reduced mitochondrial respiration activity . During exponential growth in glucose, yeast cells almost exclusively ferment, and it is only at the diauxic shift when they will normally begin to respire. In this sense, yeast cells are similar to tumor cells that obtain much of their energy from aerobic glycolysis (the Warburg effect). Like mouse Sirt6, yeast Sir2 and Hst1 appear to favor the normal shift to respiration by repressing glycolytic genes when needed. A key difference is that Sirt6 appears to deacetylate histone tails at promoters  to repress genes, while Sir2 and Hst1 associate with the ORFs and likely deacetylate non-histone proteins. Since the original submission of this manuscript, Sirt6 has also been found to localize on specific ORFs , making the parallel between yeast and mammals even more compelling.
Biochemically, Sirt6 is a histone H3-K9 deacetylase that is required for proper telomere function and maintenance . More recently, male Sirt6-transgenic mice were shown to have significantly extended lifespans , which does not happen with other sirtuin transgenics, and is consistent with the effect of increased Sir2 expression in yeast replicative life span . From a life span standpoint, Sir2 appears to be functionally more similar to Sirt6 than Sirt1, a position that is now further supported by our findings that Sir2 and Hst1 both regulate glycolytic gene expression and are required for proper telomere maintenance. It is therefore tempting to speculate that these non-traditional functions for Sir2 could be related to its role in promoting longevity.
Sirtuins and condensin/cohesin
Sir2-dependent recruitment of cohesin to the HM loci and rDNA is well established [40, 55, 77]. Mutations in cohesin subunits also result in rDNA silencing and recombination suppression defects , as well as impair the silencing boundary function of tRNA genes [78, 79]. Interestingly, the tRNAThr boundary element next to HMR-I was previously shown to be important for establishing cohesion across the HMR locus, but silencing of HMR in the absence of the tRNA gene was not sufficient for establishment . The finding that Sir2 contributes to efficient recruitment of cohesin to tRNA genes helps explain this observation. Sum1 and Hst1 are required for efficient cohesin recruitment at their shared target locations such as tRNA genes (Figure 7c), and therefore likely contribute to sister chromatid cohesion at those locations. We also found that high NAD+ levels are critical for cohesin deposition, suggesting that sirtuin activity is involved. However, a catalytically inactive Sir2 protein can establish cohesion when it is targeted as a fusion protein adjacent to HMR in place of the tRNAThr gene . Perhaps Hst1 is providing a redundant deacetylation function in this context.
We were compelled to assay for Sir2-mediated recruitment of condensin because of the highly significant overlap of Sir2, Hst1, and Sum1 binding sites with condensin binding sites, especially tRNA genes [53, 54]. The negative effects of deleting SIR2, HST1, or SUM1 on condensin association with tRNA genes, the rDNA, or specific ORF regions were almost identical to the effects on cohesin, suggesting that their recruitment to these locations may involve a common factor. One possibility is the Scc2-Scc4 complex, which loads the cohesin complex onto chromosomes , but has also been suggested to contribute to condensin loading onto tRNA genes . Consistent with this idea, Scc2 enrichment at the rDNA was impaired in an npt1Δ mutant (Figure 8c). A recent Sirt1 study showed that it was involved in the binding of condensin during mitotic chromosome condensation , implying another conserved link between yeast and mammalian sirtuin functions. Taken together, the requirements for Sir2, Hst1, and elevated NAD+ concentrations in cohesin and condensin deposition at tRNA genes and genes downregulated at the diauxic shift suggest they may have significant impact on long-range chromosomal architecture beyond their traditional silenced targets.
The budding yeast S. cerevisiae was used as a model system in identifying functional chromatin binding targets for all five sirtuins in this organism. A significant amount of overlap was surprisingly observed between Sir2 and Hst1, which was not shared with Hst2, Hst3, or Hst4. Importantly, these overlapping binding sites were not the previously described functional targets of Sir2 (silenced domains) or Hst1 (promoters of specific genes), but were instead novel targets. Binding of Sir2 and Hst1 at each new target class was functional, including telomere length maintenance at telomeric repeat clusters, diauxic shift-specific repression of glycolytic and translation factor genes through the binding to ORFs, and promoting condensin and cohesin deposition at overlapping positions. Sirtuins have numerous non-histone targets, but the results from this study highlight the idea that the chromatin targets of this enzyme class are equally broad and more highly conserved between yeast and human than previously thought.
Materials and methods
Yeast strains and media
Yeast strains were grown in YPD medium containing 2% glucose as the carbon source. All yeast growth was performed at 30°C. The strains used in this study were derived from the GRF167/JB740 strain background used for Ty1 and rDNA silencing studies , with the exception being a rap1-17 strain set in the W303 background . All strains are listed in Table S6 in Additional file 1. Deletion strains were created by deleting each ORF and replacing them with either kanMX4 or natMX4, using one-step PCR-mediated gene replacement, and then PCR confirmation. Myc tagged strains were produced by fusing 13 copies of the Myc epitope (EQKLISEEDL) at the end of each gene.
Log-phase cultures (200 ml) in YPD were cross-linked with 1% formaldehyde for 20 minutes at 30°C. Cells were pelleted by centrifugation and washed two times with cold Tris-buffered saline. The cells, in 0.6 ml of FA-140 lysis buffer (50 mM HEPES, 140 mM NaCl, 1% Triton X-100, 1 mM EDTA, 0.1% SDS, 0.1 mM phenylmethylsulfonyl fluoride, 2 mM benzamidine, 1× protease inhibitor cocktail (Sigma; St. Louis, MO, USA) were lysed with glass beads in a Mini-BeadBeater (Biospec Products; Bartlesville, OK, USA). The cell lysate was drawn off the beads, sonicated for 60 cycles (30 s 'on' at high level and 30 s 'off' per cycle) in a Bioruptor (Diagenode; Denville, NJ, USA) and spun for 10 minutes at 16,000x g in a microcentrifuge. Equivalent amounts of lysate (2.5 mg protein) were incubated overnight at 4°C with 5 µg of anti-Myc antibody (9E10) and 20 µl of protein G magnetic beads (Millipore; Billerica, MA, USA). The immunoprecipitated chromatin was then recovered and the DNA purified using a Magna ChIP™ G Chromatin Immunoprecipitation Kit (Millipore). The ChIP-sequencing libraries were made using an Illumina ChIP-Seq DNA Sample Prep Kit (catalogue number IP-102-1001), starting with 0.5 µg of DNA isolated from the immunoprecipitation step. The libraries were sequenced using an Illumina GAII system at the University of Virginia Biomolecular Research Facility. Integrative Genomics Viewer (IGV) was used to visualize the data. The scales were normalized to the Sir2 read count for each panel, based on the total number of mapped reads recovered.
Standard ChIP assays
Standard ChIP assays to confirm the sequencing results were performed as previously described . The preparation of chromatin solution from log-phase or post log-phase yeast culture is the same as described for the ChIP-seq procedure, except anti-myc antibody (9E10) and each specific antibody was combined with protein A/G conjugated sepharose beads blocked with salmon sperm DNA in the IP procedure. After incubation of chromatin solution with antibody and beads at 4°C overnight, the beads containing the immune-complex were washed (wash 1, twice with 1 ml of FA-140; wash 2, twice with 1 ml of FA-500 (the buffer was the same as FA-140 except that the NaCl concentration was increased to 500 mM); and wash 3, twice with 1 ml LiCl solution containing 10 mM Tris-HCl, pH 8.0, 250 mM LiCl, 0.5% NP-40, 0.5% sodium dodecyl sulfate, 1 mM EDTA). DNA was then eluted from the beads 2 times with 75 µl of elution buffer (5× TE plus 1% SDS). The combined DNA solution was incubated at 65°C overnight to reverse the cross-linking. The purified DNA samples were analyzed by quantitative real-time PCR, and the results normalized with the input DNA PCR signal, and indicated by relative IP in the graphs.
Analysis of ChIP-sequencing reads
Sequence reads for Sir2, Sum1, Hst1, Hst2 libraries (ChIP DNA), and Input DNA were uniquely mapped to the SGD genome assembly (sacCer2) using BWA version 0.5.7 . Reads that mapped uniquely were filtered on a phred quality score of 20, and were quantified as Sir2 (20760165), Hst1 (15530594), Sum1 (17613707), Hst2 (2637949), Hst3 (19078082), Hst4 (7210751), and Input (34090115). Repeats were allowed only on the rDNA locus on chromosome XII. Individual browsable enriched read wig files for Sir2, Sum1, Hst1, Hst2, and Input DNA were generated from the mapped read files by summing the number of overlapping reads for every genomic coordinate across the yeast genome.
Significant peaks were determined, relative to the Input DNA, using the BioConductor package, BayesPeak version 1.2.3 . Subsequently, 1,391, 1,210, 1,527, and 146 peaks were highlighted as significantly enriched regions for Sir2, Sum1, Hst1, and Hst2 libraries respectively, after the posterior probability for each peak was required to exceed 0.5. Additional computational and statistical analyses are described in the Supplemental Information section in Additional file 1. The ChIP-seq datasets from this study have been deposited in NCBI's Gene Expression Omnibus, and are accessible through the GEO series accession number [GSE41415].
Genomic DNA (10 µg) was digested with XhoI at 37°C, and then separated on a 0.7% agarose gel, followed by transfer to an Immobilon-Ny+ membrane (Millipore). The membrane was prehybridized in QuickHyb solution (Stratagene; Jolla, LA, USA) at 68°C for 20 minutes. To generate the telomeric repeat probe, a 350 bp EcoRI fragment was isolated from pYPLV , and labeled with [α32P] dCTP (3,000 Ci/mmol; Perkin Elmer; Boston, Massachusetts, USA) by random priming. The labeled probe was mixed with 100 µl of sonicated salmon sperm DNA (10 mg/ml), boiled to denature, and then hybridized to the membrane for 1 hour at 68°C in QuickHyb solution. The membrane was washed two times at room temperature in 2×-SSC + 0.1% SDS and once at 60°C in 0.1×-SSC + 0.1% SDS. The washed membrane was exposed to X-ray film for autoradiography.
Quantitative reverse transcriptase (RT) PCR assays
Synthesis of cDNA from total RNA and PCR reactions with SYBR green PCR master mix was performed as previously described . The oligonucleotide primer sequences are provided in Table S7 in Additional file 1. The test mRNA transcript levels were normalized to either ACT1 or ALD2. As indicated in some figures, to determine the fold induction, gene transcript levels in the mutant strains were also normalized to the levels in the wild-type strain. Results reflect the average fold induction (relative to the induction in the wild-type strain) from three biological replicates. Where indicated in the figures, the standard deviation was calculated.
Glucose concentration measurements
The amount of glucose in the growth medium was assayed with a glucose assay kit (Sigma) using the glucose oxidase system. One milliliter of culture was removed at each time point from the indicated cultures and centrifuged at 10,000x g for 5 minutes to clarify the supernatant. The supernatant was further diluted 250-fold with deionized water and 250 µl of the diluted sample was subjected to the assay. Briefly, the reaction was started by adding 500 µl of assay reagent that contains o-dianisidine and glucose oxidase/peroxidase. After reacting for exactly 30 minutes at 37°C, the reaction was stopped by adding 500 µl of 12 N H2SO4, followed by measuring the absorbance at 540 nM.
Intracellular NAD+ and NADH measurements
To determine the NAD+/NADH ratio, we utilized a fluorescent NAD/NADH detection kit (Cell Technology, Inc; Mountain View, CA). Yeast cells were inoculated into 100 ml YPD medium and grown at 30°C in a shaker. A quantity of 2 × 106 cells was collected for each time point and then washed twice with 2 ml of phosphate-buffered saline. After removal of the final supernatant, the cell pellet was resuspended in 200 µl of the NAD or NADH extraction buffer supplied in the kit, and the rest of the protocol performed according to the manufacturer's instructions.
Integrative Genomics Viewer
open reading frame
Houtkooper RH, Pirinen E, Auwerx J: Sirtuins as regulators of metabolism and healthspan. Nat Rev Mol Cell Biol. 2012, 13: 225-238.
Landry J, Slama JT, Sternglanz R: Role of NAD+ in the deacetylase activity of the SIR2-like proteins. Biochem Biophys Res Commun. 2000, 278: 685-690. 10.1006/bbrc.2000.3854.
Tanny JC, Moazed D: Coupling of histone deacetylation to NAD breakdown by the yeast silencing protein Sir2: Evidence for acetyl transfer from substrate to an NAD breakdown product. Proc Natl Acad Sci USA. 2001, 98: 415-420. 10.1073/pnas.98.2.415.
Lin SJ, Defossez PA, Guarente L: Requirement of NAD and SIR2 for life-span extension by calorie restriction in Saccharomyces cerevisiae. Science. 2000, 289: 2126-2128. 10.1126/science.289.5487.2126.
Smith JS, Brachmann CB, Celic I, Kenna MA, Muhammad S, Starai VJ, Avalos JL, Escalante-Semerena JC, Grubmeyer C, Wolberger C, Boeke JD: A phylogenetically conserved NAD+-dependent protein deacetylase activity in the Sir2 protein family. Proc Natl Acad Sci USA. 2000, 97: 6658-6663. 10.1073/pnas.97.12.6658.
Frye RA: Phylogenetic classification of prokaryotic and eukaryotic Sir2-like proteins. Biochem Biophys Res Commun. 2000, 273: 793-798. 10.1006/bbrc.2000.3000.
Brachmann CB, Sherman JM, Devine SE, Cameron EE, Pillus L, Boeke JD: The SIR2 gene family, conserved from bacteria to humans, functions in silencing, cell cycle progression, and chromosome stability. Genes Dev. 1995, 9: 2888-2902. 10.1101/gad.9.23.2888.
Rusche LN, Kirchmaier AL, Rine J: The establishment, inheritance, and function of silenced chromatin in Saccharomyces cerevisiae. Annu Rev Biochem. 2003, 72: 481-516. 10.1146/annurev.biochem.72.121801.161547.
Imai S, Armstrong CM, Kaeberlein M, Guarente L: Transcriptional silencing and longevity protein Sir2 is an NAD-dependent histone deacetylase. Nature. 2000, 403: 795-800. 10.1038/35001622.
Landry J, Sutton A, Tafrov ST, Heller RC, Stebbins J, Pillus L, Sternglanz R: The silencing protein SIR2 and its homologs are NAD-dependent protein deacetylases. Proc Natl Acad Sci USA. 2000, 97: 5807-5811. 10.1073/pnas.110148297.
Kaeberlein M, McVey M, Guarente L: The SIR2/3/4 complex and SIR2 alone promote longevity in Saccharomyces cerevisiae by two different mechanisms. Genes Dev. 1999, 13: 2570-2580. 10.1101/gad.13.19.2570.
Mostoslavsky R, Chua KF, Lombard DB, Pang WW, Fischer MR, Gellon L, Liu P, Mostoslavsky G, Franco S, Murphy MM, Mills KD, Patel P, Hsu JT, Hong AL, Ford E, Cheng HL, Kennedy C, Nunez N, Bronson R, Frendewey D, Auerbach W, Valenzuela D, Karow M, Hottiger MO, Hursting S, Barrett JC, Guarente L, Mulligan R, Demple B, Yancopoulos GD, Alt FW: Genomic instability and aging-like phenotype in the absence of mammalian SIRT6. Cell. 2006, 124: 315-329. 10.1016/j.cell.2005.11.044.
Kanfi Y, Naiman S, Amir G, Peshti V, Zinman G, Nahum L, Bar-Joseph Z, Cohen HY: The sirtuin SIRT6 regulates lifespan in male mice. Nature. 2012, 483: 218-221. 10.1038/nature10815.
Ghidelli S, Donze D, Dhillon N, Kamakaka RT: Sir2p exists in two nucleosome-binding complexes with distinct deacetylase activities. EMBO J. 2001, 20: 4522-4535. 10.1093/emboj/20.16.4522.
Hoppe GJ, Tanny JC, Rudner AD, Gerber SA, Danaie S, Gygi SP, Moazed D: Steps in assembly of silent chromatin in yeast: Sir3-independent binding of a Sir2/Sir4 complex to silencers and role for Sir2-dependent deacetylation. Mol Cell Biol. 2002, 22: 4167-4180. 10.1128/MCB.22.12.4167-4180.2002.
Strahl-Bolsinger S, Hecht A, Luo K, Grunstein M: SIR2 and SIR4 interactions differ in core and extended telomeric heterochromatin in yeast. Genes Dev. 1997, 11: 83-93. 10.1101/gad.11.1.83.
Aparicio OM, Billington BL, Gottschling DE: Modifiers of position effect are shared between telomeric and silent mating-type loci in S. cerevisiae. Cell. 1991, 66: 1279-1287. 10.1016/0092-8674(91)90049-5.
Rine J, Herskowitz I: Four genes responsible for a position effect on expression from HML and HMR in Saccharomyces cerevisiae. Genetics. 1987, 116: 9-22.
Shou W, Seol JH, Shevchenko A, Baskerville C, Moazed D, Chen ZW, Jang J, Charbonneau H, Deshaies RJ: Exit from mitosis is triggered by Tem1-dependent release of the protein phosphatase Cdc14 from nucleolar RENT complex. Cell. 1999, 97: 233-244. 10.1016/S0092-8674(00)80733-3.
Straight AF, Shou W, Dowd GJ, Turck CW, Deshaies RJ, Johnson AD, Moazed D: Net1, a Sir2-associated nucleolar protein required for rDNA silencing and nucleolar integrity. Cell. 1999, 97: 245-256. 10.1016/S0092-8674(00)80734-5.
Bedalov A, Hirao M, Posakony J, Nelson M, Simon JA: NAD+-dependent deacetylase Hst1p controls biosynthesis and cellular NAD+ levels in Saccharomyces cerevisiae. Mol Cell Biol. 2003, 23: 7044-7054. 10.1128/MCB.23.19.7044-7054.2003.
McCord R, Pierce M, Xie J, Wonkatal S, Mickel C, Vershon AK: Rfm1, a novel tethering factor required to recruit the Hst1 histone deacetylase for repression of middle sporulation genes. Mol Cell Biol. 2003, 23: 2009-2016. 10.1128/MCB.23.6.2009-2016.2003.
Zill OA, Rine J: Interspecies variation reveals a conserved repressor of alpha-specific genes in Saccharomyces yeasts. Genes Dev. 2008, 22: 1704-1716. 10.1101/gad.1640008.
Vaquero A, Scher M, Lee D, Erdjument-Bromage H, Tempst P, Reinberg D: Human SirT1 interacts with histone H1 and promotes formation of facultative heterochromatin. Mol Cell. 2004, 16: 93-105. 10.1016/j.molcel.2004.08.031.
Mulligan P, Yang F, Di Stefano L, Ji JY, Ouyang J, Nishikawa JL, Toiber D, Kulkarni M, Wang Q, Najafi-Shoushtari SH, Mostoslavsky R, Gygi SP, Gill G, Dyson NJ, Näär AM: A SIRT1-LSD1 corepressor complex regulates Notch target gene expression and development. Mol Cell. 2011, 42: 689-699. 10.1016/j.molcel.2011.04.020.
Li M, Petteys BJ, McClure JM, Valsakumar V, Bekiranov S, Frank EL, Smith JS: Thiamine biosynthesis in Saccharomyces cerevisiae is regulated by the NAD+-dependent histone deacetylase Hst1. Mol Cell Biol. 2010, 30: 3329-3341. 10.1128/MCB.01590-09.
Teytelman L, Ozaydin B, Zill O, Lefrancois P, Snyder M, Rine J, Eisen MB: Impact of chromatin structures on DNA processing for genomic analyses. PLoS One. 2009, 4: e6700-10.1371/journal.pone.0006700.
Bi X, Braunstein M, Shei GJ, Broach JR: The yeast HML I silencer defines a heterochromatin domain boundary by directional establishment of silencing. Proc Natl Acad Sci USA. 1999, 96: 11934-11939. 10.1073/pnas.96.21.11934.
Donze D, Kamakaka RT: RNA polymerase III and RNA polymerase II promoter complexes are heterochromatin barriers in Saccharomyces cerevisiae. EMBO J. 2001, 20: 520-531. 10.1093/emboj/20.3.520.
Rusche LN, Kirchmaier AL, Rine J: Ordered nucleation and spreading of silenced chromatin in Saccharomyces cerevisiae. Mol Biol Cell. 2002, 13: 2207-2222. 10.1091/mbc.E02-03-0175.
Fourel G, Lebrun E, Gilson E: Protosilencers as building blocks for heterochromatin. Bioessays. 2002, 24: 828-835. 10.1002/bies.10139.
Lazarus AG, Holmes SG: A cis-acting tRNA gene imposes the cell cycle progression requirement for establishing silencing at the HMR locus in yeast. Genetics. 2011, 187: 425-439. 10.1534/genetics.110.124099.
Xie J, Pierce M, Gailus-Durner V, Wagner M, Winter E, Vershon AK: Sum1 and Hst1 repress middle sporulation-specific gene expression during mitosis in Saccharomyces cerevisiae. EMBO J. 1999, 18: 6448-6454. 10.1093/emboj/18.22.6448.
Lynch PJ, Fraser HB, Sevastopoulos E, Rine J, Rusche LN: Sum1p, the origin recognition complex, and the spreading of a promoter-specific repressor in Saccharomyces cerevisiae. Mol Cell Biol. 2005, 25: 5920-5932. 10.1128/MCB.25.14.5920-5932.2005.
Hickman MA, Rusche LN: Substitution as a mechanism for genetic robustness: the duplicated deacetylases Hst1p and Sir2p in Saccharomyces cerevisiae. PLoS Genet. 2007, 3: e126-10.1371/journal.pgen.0030126.
Bryk M, Banerjee M, Murphy M, Knudsen KE, Garfinkel DJ, Curcio MJ: Transcriptional silencing of Ty1 elements in the RDN1 locus of yeast. Genes Dev. 1997, 11: 255-269. 10.1101/gad.11.2.255.
Smith JS, Boeke JD: An unusual form of transcriptional silencing in yeast ribosomal DNA. Genes Dev. 1997, 11: 241-254. 10.1101/gad.11.2.241.
Buck SW, Sandmeier JJ, Smith JS: RNA polymerase I propagates unidirectional spreading of rDNA silent chromatin. Cell. 2002, 111: 1003-1014. 10.1016/S0092-8674(02)01193-5.
Huang J, Moazed D: Association of the RENT complex with nontranscribed and coding regions of rDNA and a regional requirement for the replication fork block protein Fob1 in rDNA silencing. Genes Dev. 2003, 17: 2162-2176. 10.1101/gad.1108403.
Huang J, Brito IL, Villen J, Gygi SP, Amon A, Moazed D: Inhibition of homologous recombination by a cohesin-associated clamp complex recruited to the rDNA recombination enhancer. Genes Dev. 2006, 20: 2887-2901. 10.1101/gad.1472706.
Hecht A, Laroche T, Strahl-Bolsinger S, Gasser SM, Grunstein M: Histone H3 and H4 N-termini interact with SIR3 and SIR4 proteins: a molecular model for the formation of heterochromatin in yeast. Cell. 1995, 80: 583-592. 10.1016/0092-8674(95)90512-X.
Moretti P, Freeman K, Coodly L, Shore D: Evidence that a complex of SIR proteins interacts with the silencer and telomere-binding protein RAP1. Genes Dev. 1994, 8: 2257-2269. 10.1101/gad.8.19.2257.
Spyrou C, Stark R, Lynch AG, Tavare S: BayesPeak: Bayesian analysis of ChIP-seq data. BMC Bioinformatics. 2009, 10: 299-10.1186/1471-2105-10-299.
Kyrion G, Liu K, Liu C, Lustig AJ: RAP1 and telomere structure regulate telomere position effects in Saccharomyces cerevisiae. Genes Dev. 1993, 7: 1146-1159. 10.1101/gad.7.7a.1146.
Cockell M, Palladino F, Laroche T, Kyrion G, Liu C, Lustig AJ, Gasser SM: The carboxy termini of Sir4 and Rap1 affect Sir3 localization: evidence for a multicomponent complex required for yeast telomeric silencing. J Cell Biol. 1995, 129: 909-924. 10.1083/jcb.129.4.909.
Perrod S, Cockell MM, Laroche T, Renauld H, Ducrest AL, Bonnard C, Gasser SM: A cytosolic NAD-dependent deacetylase, Hst2p, can modulate nucleolar and telomeric silencing in yeast. EMBO J. 2001, 20: 197-209. 10.1093/emboj/20.1.197.
Orzechowski Westholm J, Tronnersjo S, Nordberg N, Olsson I, Komorowski J, Ronne H: Gis1 and Rph1 regulate glycerol and acetate metabolism in glucose depleted yeast cells. PLoS One. 2012, 7: e31577-10.1371/journal.pone.0031577.
DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686. 10.1126/science.278.5338.680.
Lin SJ, Ford E, Haigis M, Liszt G, Guarente L: Calorie restriction extends yeast life span by lowering the level of NADH. Genes Dev. 2004, 18: 12-16. 10.1101/gad.1164804.
Kennedy BK, Gotta M, Sinclair DA, Mills K, McNabb DS, Murthy M, Pak SM, Laroche T, Gasser SM, Guarente L: Redistribution of silencing proteins from telomeres to the nucleolus is associated with extension of life span in S. cerevisiae. Cell. 1997, 89: 381-391. 10.1016/S0092-8674(00)80219-6.
Smith JS, Brachmann CB, Pillus L, Boeke JD: Distribution of a limited Sir2 protein pool regulates the strength of yeast rDNA silencing and is modulated by Sir4p. Genetics. 1998, 149: 1205-1219.
Raab JR, Chiu J, Zhu J, Katzman S, Kurukuti S, Wade PA, Haussler D, Kamakaka RT: Human tRNA genes function as chromatin insulators. EMBO J. 2012, 31: 330-350.
D'Ambrosio C, Schmidt CK, Katou Y, Kelly G, Itoh T, Shirahige K, Uhlmann F: Identification of cis-acting sites for condensin loading onto budding yeast chromosomes. Genes Dev. 2008, 22: 2215-2227. 10.1101/gad.1675708.
Haeusler RA, Pratt-Hyatt M, Good PD, Gipson TA, Engelke DR: Clustering of yeast tRNA genes is mediated by specific association of condensin with tRNA gene transcription complexes. Genes Dev. 2008, 22: 2204-2214. 10.1101/gad.1675908.
Kobayashi T, Horiuchi T, Tongaonkar P, Vu L, Nomura M: SIR2 regulates recombination between different rDNA repeats, but not recombination within individual rRNA genes in yeast. Cell. 2004, 117: 441-453. 10.1016/S0092-8674(04)00414-3.
Wu CS, Chen YF, Gartenberg MR: Targeted sister chromatid cohesion by Sir2. PLoS Genet. 2011, 7: e1002000-10.1371/journal.pgen.1002000.
Rusche LN, Rine J: Conversion of a gene-specific repressor to a regional silencer. Genes Dev. 2001, 15: 955-967. 10.1101/gad.873601.
Sutton A, Heller RC, Landry J, Choy JS, Sirko A, Sternglanz R: A novel form of transcriptional silencing by Sum1-1 requires Hst1 and the origin recognition complex. Mol Cell Biol. 2001, 21: 3514-3522. 10.1128/MCB.21.10.3514-3522.2001.
Irlbacher H, Franke J, Manke T, Vingron M, Ehrenhofer-Murray AE: Control of replication initiation and heterochromatin formation in Saccharomyces cerevisiae by a regulator of meiotic gene expression. Genes Dev. 2005, 19: 1811-1822. 10.1101/gad.334805.
Askree SH, Yehuda T, Smolikov S, Gurevich R, Hawk J, Coker C, Krauskopf A, Kupiec M, McEachern MJ: A genome-wide screen for Saccharomyces cerevisiae deletion mutants that affect telomere length. Proc Natl Acad Sci USA. 2004, 101: 8658-8663. 10.1073/pnas.0401263101.
Gatbonton T, Imbesi M, Nelson M, Akey JM, Ruderfer DM, Kruglyak L, Simon JA, Bedalov A: Telomere length as a quantitative trait: genome-wide survey and genetic mapping of telomere length-control genes in yeast. PLoS Genet. 2006, 2: e35-10.1371/journal.pgen.0020035.
Koh JL, Ding H, Costanzo M, Baryshnikova A, Toufighi K, Bader GD, Myers CL, Andrews BJ, Boone C: DRYGIN: a database of quantitative genetic interaction networks in yeast. Nucleic Acids Res. 2010, 38: D502-507. 10.1093/nar/gkp820.
Bolton EC, Boeke JD: Transcriptional interactions between yeast tRNA genes, flanking genes and Ty elements: a genomic point of view. Genome Res. 2003, 13: 254-263. 10.1101/gr.612203.
Dubarry M, Loiodice I, Chen CL, Thermes C, Taddei A: Tight protein-DNA interactions favor gene silencing. Genes Dev. 2011, 25: 1365-1370. 10.1101/gad.611011.
Hull MW, Erickson J, Johnston M, Engelke DR: tRNA genes as transcriptional repressor elements. Mol Cell Biol. 1994, 14: 1266-1277.
Wang L, Haeusler RA, Good PD, Thompson M, Nagar S, Engelke DR: Silencing near tRNA genes requires nucleolar localization. J Biol Chem. 2005, 280: 8637-8639.
Gard S, Light W, Xiong B, Bose T, McNairn AJ, Harris B, Fleharty B, Seidel C, Brickner JH, Gerton JL: Cohesinopathy mutations disrupt the subnuclear organization of chromatin. J Cell Biol. 2009, 187: 455-462. 10.1083/jcb.200906075.
Galdieri L, Mehrotra S, Yu S, Vancura A: Transcriptional regulation in yeast during diauxic shift and stationary phase. OMICS. 2010, 14: 629-638. 10.1089/omi.2010.0069.
Xie Z, Dai J, Dai L, Tan M, Cheng Z, Wu Y, Boeke JD, Zhao Y: Lysine succinylation and lysine malonylation in histones. Mol Cell Proteomics. 2012, 11: 100-107. 10.1074/mcp.M111.015875.
Liu Z, Cao J, Gao X, Zhou Y, Wen L, Yang X, Yao X, Ren J, Xue Y: CPLA 1.0: an integrated database of protein lysine acetylation. Nucleic Acids Res. 2011, 39: D1029-1034. 10.1093/nar/gkq939.
Lieb JD, Liu X, Botstein D, Brown PO: Promoter-specific binding of Rap1 revealed by genome-wide maps of protein-DNA association. Nat Genet. 2001, 28: 327-334. 10.1038/ng569.
Radman-Livaja M, Ruben G, Weiner A, Friedman N, Kamakaka R, Rando OJ: Dynamics of Sir3 spreading in budding yeast: secondary recruitment sites and euchromatic localization. EMBO J. 2011, 30: 1012-1026. 10.1038/emboj.2011.30.
Tsankov AM, Brown CR, Yu MC, Win MZ, Silver PA, Casolari JM: Communication between levels of transcriptional control improves robustness and adaptivity. Mol Syst Biol. 2006, 2: 65-
Zhong L, D'Urso A, Toiber D, Sebastian C, Henry RE, Vadysirisack DD, Guimaraes A, Marinelli B, Wikstrom JD, Nir T, Clish CB, Vaitheesvaran B, Iliopoulos O, Kurland I, Dor Y, Weissleder R, Shirihai OS, Ellisen LW, Espinosa JM, Mostoslavsky R: The histone deacetylase Sirt6 regulates glucose homeostasis via Hif1α. Cell. 2010, 140: 280-293. 10.1016/j.cell.2009.12.041.
Sebastian C, Zwaans BM, Silberman DM, Gymrek M, Goren A, Zhong L, Ram O, Truelove J, Guimaraes AR, Toiber D, Cosentino C, Greenson JK, MacDonald AI, McGlynn L, Maxwell F, Edwards J, Giacosa S, Guccione E, Weissleder R, Bernstein BE, Regev A, Shiels PG, Lombard DB, Mostoslavsky R: The histone deacetylase SIRT6 is a tumor suppressor that controls cancer metabolism. Cell. 2012, 151: 1185-1199. 10.1016/j.cell.2012.10.047.
Michishita E, McCord RA, Berber E, Kioi M, Padilla-Nash H, Damian M, Cheung P, Kusumoto R, Kawahara TL, Barrett JC, Chang HY, Bohr VA, Ried T, Gozani O, Chua KF: SIRT6 is a histone H3 lysine 9 deacetylase that modulates telomeric chromatin. Nature. 2008, 452: 492-496. 10.1038/nature06736.
Chang CR, Wu CS, Hom Y, Gartenberg MR: Targeting of cohesin by transcriptionally silent chromatin. Genes Dev. 2005, 19: 3031-3042. 10.1101/gad.1356305.
Biswas M, Maqani N, Rai R, Kumaran SP, Iyer KR, Sendinc E, Smith JS, Laloraya S: Limiting the extent of the RDN1 heterochromatin domain by a silencing barrier and Sir2 protein levels in Saccharomyces cerevisiae. Mol Cell Biol. 2009, 29: 2889-2898. 10.1128/MCB.00728-08.
Donze D, Adams CR, Rine J, Kamakaka RT: The boundaries of the silenced HMR domain in Saccharomyces cerevisiae. Genes Dev. 1999, 13: 698-708. 10.1101/gad.13.6.698.
Dubey RN, Gartenberg MR: A tDNA establishes cohesion of a neighboring silent chromatin domain. Genes Dev. 2007, 21: 2150-2160. 10.1101/gad.1583807.
Ciosk R, Shirayama M, Shevchenko A, Tanaka T, Toth A, Nasmyth K: Cohesin's binding to chromosomes depends on a separate complex consisting of Scc2 and Scc4 proteins. Mol Cell. 2000, 5: 243-254. 10.1016/S1097-2765(00)80420-7.
Fatoba ST, Okorokov AL: Human SIRT1 associates with mitotic chromatin and contributes to chromosomal condensation. Cell Cycle. 2011, 10: 2317-2322. 10.4161/cc.10.14.15913.
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.
Wellinger RJ, Wolf AJ, Zakian VA: Saccharomyces telomeres acquire single-strand TG1-3 tails late in S phase. Cell. 1993, 72: 51-60. 10.1016/0092-8674(93)90049-V.
We would like to thank Marty Mayo, Natalya Baranova, Neerja Karnani, and Cathleen Brdlik for advice with ChIP-sequencing, and Patrick Grant for use of their Bioruptor apparatus. We also thank David Shore and Danesh Moazed for providing strains and plasmids. Lastly, we thank members of the Smith, Bekiranov, and Grant labs for helpful discussions and suggestions, as well as Mitch Smith and Marty Mayo for suggestions on the project and critical comments on the manuscript.
The authors declare that they have no competing interests.
ML contributed to conceptualization of the experiments, carried out most of the molecular biological experiments, and drafting of the manuscript. VV performed bioinformatic and statistical analysis with the ChIP-seq data and contributed to the writing. KP contributed to statistical analysis of the ChIP-seq data and diauxic shift gene expression analysis. SB contributed to conceptualization of the experiments, drafting the manuscript, and data analysis. JS conceived of the study, contributed to the experimental design, data analysis, and drafting of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Supplementary online material. This file contains Supplemental materials and methods, Figures S1 to S7, and Tables S1 to S7. (PDF 2 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Li, M., Valsakumar, V., Poorey, K. et al. Genome-wide analysis of functional sirtuin chromatin targets in yeast. Genome Biol 14, R48 (2013). https://doi.org/10.1186/gb-2013-14-5-r48
- Telomere Length
- tRNA Gene
- NADH Ratio
- Diauxic Shift
- Replicative Life Span