|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: HP MC BL. Performed the experiments: HP MD MC. Analyzed the data: BL MC SA LK RR HP MD. Contributed reagents/materials/analysis tools: HP MD MC. Wrote the paper: BL.
Using a genome-wide screening approach, we have established the genetic requirements for proper telomere structure in Saccharomyces cerevisiae. We uncovered 112 genes, many of which have not previously been implicated in telomere function, that are required to form a fold-back structure at chromosome ends. Among other biological processes, lysine deacetylation, through the Rpd3L, Rpd3S, and Hda1 complexes, emerged as being a critical regulator of telomere structure. The telomeric-bound protein, Rif2, was also found to promote a telomere fold-back through the recruitment of Rpd3L to telomeres. In the absence of Rpd3 function, telomeres have an increased susceptibility to nucleolytic degradation, telomere loss, and the initiation of premature senescence, suggesting that an Rpd3-mediated structure may have protective functions. Together these data reveal that multiple genetic pathways may directly or indirectly impinge on telomere structure, thus broadening the potential targets available to manipulate telomere function.
Impaired telomere elongation eventually results in telomere dysfunction and can lead to diseases such as dyskeratosis congenita, which is associated with bone-marrow failure and pulmonary fibrosis. Cancer cells require continuous telomere maintenance to ensure continued cellular proliferation. Therefore the regulation of telomere function, both positively (in the case of dyskeratosis congenita) and negatively (for cancer), may be of therapeutic benefit. In this study we have used yeast to determine which genetic factors are important for a certain telomeric structure (the loop structure), which may help to maintain chromosome ends in a protected state. We found that multiple genetic factors and pathways affect telomere structure, ranging from metabolic signaling to specific telomere-binding proteins. We found that proper chromatin structure at the telomere is essential to maintain a telomere fold-back structure. Importantly, there was a strong correlation between telomere structure and function, as the mutants found in our screen (looping defective) were often associated with rapid senescence and telomere dysfunction phenotypes. We believe that, through the regulation of the various genetic pathways uncovered in our screen, one may be able to both positively and negatively influence telomere function.
The physical ends of linear chromosomes resemble double-strand breaks (DSBs) in many respects with the exception that DSBs result in the activation of the DNA damage response and are eventually subject to repair; activities to which telomeres are refractory . This essential quality of telomeres exists as a result of their repetitive sequence that is bound by specific proteins (shelterin and CST complexes), which in turn inhibit DNA damage checkpoints, DNA repair activities and exonuclease-mediated degradation –. In yeast, the CST (Cdc13-Stn1-Ten1) complex is essential for viability and prevents the accessibility of 5′ exonucleases (primarily Exo1) to the telomere –. Upon inactivation of CST with temperature-sensitive alleles of CDC13 and STN1, cells undergo a DNA damage-mediated checkpoint arrest due to the accumulation of single-stranded (ss) telomeric DNA –. In parallel, the Rap1, Rif1 and Rif2 complex, also contribute to telomere end protection by limiting telomeric ssDNA accumulation and subsequent checkpoint activation –.
In most human somatic cells, telomeres shorten during each cell division due, in part, to the end-replication problem –. Eventually, the loss of telomeric DNA leads to telomere dysfunction, checkpoint activation and cellular senescence. Some cell types as well as most cancer cells avoid telomere attrition-induced senescence by expressing the specialized reverse transcriptase, telomerase. Telomerase elongates telomeres through the iterative addition of short sequence repeats to the 3′ ends of telomeres, compensating for the end-replication problem . Wild type S. cerevisiae constitutively express telomerase, however the cellular senescence phenotype can be induced following its inactivation/deletion .
In yeast, reporter genes become silenced when placed in the vicinity of telomeres . This telomere-induced silencing is dependent on the Sir2/3/4 lysine deacetylation (KDAC) complex, which is recruited to chromosome ends via the telomere binding protein, Rap1 . Apart from the Sir2/3/4 complex, other KDACs also contribute to the heterochromatic constitution of telomeres and sub-telomeres. The class I KDAC, Rpd3 (the yeast ortholog of human KDAC1), consisting of two sub-complexes, Rpd3L and Rpd3S , also localizes to telomeres and is important to establish the euchromatin/heterochromatin boundary in the sub-telomeric regions , as well as to prevent hyper-silencing . The class II KDAC, Hda1, also contributes to chromatin regulation at yeast telomeres . The relationships between heterochromatin and telomere structure/function remain unclear.
It has been postulated that telomere protection may stem, in part, from a higher-order chromatin structure. Analysis of telomeric DNA from human and mouse cells has revealed that the telomere terminus can be hidden in a lariat-like structure termed a t-loop , –. T-loops, thought to form through the strand invasion of the 3′ telomeric overhang into the double-stranded region of the telomere, have also been found in chickens, worms, plants, and protozoa –. Via electron microscopy, telomeric loops have also been observed in yeast (K. lactis) with over-elongated telomeres , and the telomere associated S. pombe protein, Taz1, has been shown to re-model model DNA substrates into t-loops . However, due to the small size of yeast telomeres it has been difficult to both prepare and analyze wild type length yeast telomeres via electron microscopy . In the budding yeast, S. cerevisiae, both genetic and chromatin immunoprecipitation-based experiments have revealed that wild type telomeres do fold-back onto themselves and into the subtelomeric region , –, suggesting that loops or fold-back structures are indeed important for telomere function in yeast. Apart from the Sir2/3/4 deacetylase complex in budding yeast being important for this fold-back , and the shelterin component, TRF2 in human cells being required for t-loop formation , , the regulation of such telomeric structures remains poorly understood. In this study we have taken an unbiased genome-wide screening approach in yeast to better understand how telomere structure/fold-back is regulated in vivo. We demonstrate that multiple biological processes influence telomere structure, including the state of the subtelomeric heterochromatin as dictated by multiple lysine deacetylases. Furthermore, we find that there are direct correlations between the inability of a telomere to fold-back and telomere dysfunction, implying that the loop structure may make important contributions to telomere protection.
By placing a TATA-less galactose-inducible UAS (upstream activating sequence) downstream of the URA3 gene (from hereon referred to as construct 2), URA3 transcription is only achieved when the UAS loops back and comes into proximity with the URA3 promoter – (Figure 1A). Fold-back-induced transcription only takes place when this construct is integrated at the telomere and does not occur when it is integrated at an internal chromosomal locus . Transcription of URA3 results in lethality on media containing the drug 5-fluoroortic acid (5-FOA), providing a robust readout (cell death on 5-FOA) for successful telomere looping.
To better understand how the telomere fold-back structure in yeast is regulated, we introduced construct 2 into the yeast haploid deletion collection using the synthetic genetic array (SGA) procedure , resulting in the construction of ~4800 haploid deletion mutants harboring construct 2 (Figure 1B). Robotic pinning of these strains in quadruplicate onto galactose media in the presence and absence of 5-FOA revealed potential looping defective mutants that grew on 5-FOA (Figure 1B, bottom panel example of looping defective mutant). All positively scoring mutants were independently re-constructed and spotted as serial dilutions onto galactose +/− 5-FOA media in duplicate. We confirmed 112 yeast mutants that were defective for telomere looping and subsequently ranked them qualitatively for growth on 5-FOA (Table 1, Figure S1A). Using the Cytoscape BinGO plugin , the statistically over-represented GO (gene ontology) categories were determined for our positive scoring candidates (Figure 1C). The confirmed mutants formed the “positive hit set” whereas the “reference set” consisted of all 4800 genes screened. The analysis used a hypergeometric test and significance was tested at 5% (p<0.05) after applying Benjamini & Hochberg False Discovery Rate (FDR) correction for multiple testing. This protocol revealed histone deacetylation as a significantly enriched GO term in the GO-Cellular Component and the GO-Biological Process ontologies. Moreover, the Rpd3L, Rpd3S and Hda1 KDAC complexes were specifically over-represented (Figure 1C, Table 1). We introduced the looping construct into deletion mutants of all members of the Rpd3L/S and Hda1 complexes (including those that did not score positive in the screen) and determined that all complex members tested were important for wild type-like telomere structure (Figure 1D). Importantly, we replicated the 5-FOA plates onto media lacking uracil to ensure that the FOA resistance observed was not due to inactivation of the URA3 gene (Figure 1D). FOA resistant strains maintain the ability to grow on media lacking uracil due to low basal levels of the URA3 transcript (see Figure S1D). To confirm that there was not an inherent problem of inducing transcription within the subtelomere of these mutants, we generated and introduced construct 4 (Figure 1E) into Rpd3L, Rpd3S and Hda1 mutants where the TATA-less UAS was placed upstream of URA3 and found that (unlike with construct 2) upon galactose induction all mutants were dead on 5-FOA containing media (Figure 1E). Figure 1E demonstrates that in the mutants of the Rpd3L/S and Hda1 complexes, if the UAS in construct 2 were able to loop back to the URA3 promoter it would be able to induce transcription to an extent that would result in cell death on FOA, as is the case with wild type cells. We excluded that telomere length variation may affect the looping read-out in the rpd3Δ and hda1Δ mutants, as no significant changes in telomere length were detectable when comparing the KDAC mutants to isogenic wild type cells (Figure S1B, S1C). Finally, we demonstrated that the plate read-out effects that we have observed with construct 2 on 5-FOA are due to changes in levels of the URA3 transcript (Figure S1D) as has previously been reported  and not an unrelated artifact of 5-FOA.
In conclusion, an unbiased genome-wide screen has implicated lysine deacetylation through the Rpd3L, Rpd3S and Hda1 complexes in promoting a structural change (likely a fold-back) at budding yeast telomeres.
To demonstrate that the looping defect we observe is not specific to modified telomere 7L (construct 2), we employed a previously established chromatin immunoprecipitation (ChIP) technique where it has been shown at a natural telomere (telomere 6R) that α-Rap1 antibodies are able to precipitate subtelomeric DNA greater than 2 kb away from the start of the telomeric tract, despite the fact that chromatin is sheared into fragments of 0.5 kb . From this study, it was concluded that the subtelomeric ChIP signal from cross-linked Rap1 extracts was a result of the telomere looping back into the subtelomeric region (Figure 2A, top). We predicted that the subtelomeric signal would be lost in mutants identified in our above-described screen (Figure 2A, bottom). In agreement with previous reports, we could detect cross-linked Rap1 at a position 0.5 kb and to a lesser extent 1 kb away from the subtelomere/telomere transition point in wild type cells, indicative of a telomeric loop-back structure (Figure 2B). Unlike the previous report , we did not detect reproducible differences at positions farther than 1.5 kb from the telomere (not shown), which is likely due to the smaller chromatin fragment size used in our ChIP protocol. Strikingly, the Rap1 signal was diminished in the subtelomere in hda1Δ mutants and lost to a greater extent in sin3Δ cells (Rpd3L/S common subunit), consistent with the 5-FOA assay using construct 2 (Figure 2B, Figure 1D). sir4Δ cells were used as a looping defective positive control for the Rap1 ChIP assay (Figure 2B)  and indeed displayed the greatest loss of Rap1 signal in the subtelomeric region. Importantly, the bulk of our chromatin was sheared to 0.3 kbp fragments or less, excluding the possibility that our subtelomeric signals come from inefficient sonication (Figure S2A). Furthermore, Rap1 protein levels were not affected in any of the above-mentioned mutant backgrounds (Figure S2B).
In order to rule out the unlikely possibility that Rap1 spreading into the subtelomere may account for a portion of the ChIP signal in the assay described above (Figure 2A), we repeated the ChIP experiments using an epitope-tagged Cdc13-TAP (Tandem Affinity Purification) allele. Cdc13 associates with the 3′ ssDNA telomeric overhang, and therefore is not prone to spread into the subtelomere. Furthermore we reasoned that by using Cdc13-TAP we would be able to more easily reconcile differences at the -1 kb position due to its distal positioning at the 3′end. Consistently we were also able to detect Cdc13-TAP associated with subtelomeric DNA 1000 bp upstream of the telomeric tract (indicative of a fold-back), and the signal was reduced to background levels in the sin3Δ mutant (Figure 2C). The sin3Δ mutation did not result in decreased Cdc13-TAP protein levels, which could have potentially accounted for the reduced ChIP signal (Figure S2C). It is important to note that we did not enrich significant amounts of subtelomeric DNA at the -6 and -500 bp positions following the Cdc13-TAP ChIP (Figure 2C) although Cdc13 has been previously shown by similar methods to localize to subtelomeres . We interpret this to indicate that our sonication was extremely efficient and fragments above 400 bp (the approximate length of the wild type telomere) were extremely rare and did not give a signal significantly above background (untagged control). In order verify this notion we deleted telomerase (EST2) in Cdc13-TAP cells and let telomeres shorten over 25 and 50 generations (Figure S2D and S2E). We predicted that upon telomere shortening we would be able to increasingly detect a signal at the -6 position as telomeres would be shorter than our sheared chromatin fragments (see Figure 2C for visualization). Indeed, we found that as telomere length decreased to under 300 bp (our average chromatin fragment size) we were able to detect a robust Cdc13-TAP ChIP signal at the -6 position at natural telomere 6R (Figure S2F).
In summary we have confirmed that the looping defect observed in hda1Δ and sin3Δ (Rpd3L/S subunit) mutants using construct 2 and 5-FOA as a read-out (Figure 1D) can be recapitulated using an independent method (ChIP) at natural telomere 6R (Figure 2B, 2C).
Among the list of looping defective mutants we were intrigued that RIF1, a regulator of telomere length, was also implicated in promoting a fold-back structure (Table 1). Since Rif2 works in parallel with Rif1 to regulate telomere length we introduced construct 2 into rif2Δ cells and found that like rif1Δ cells, rif2Δ mutants also displayed a looping defect (Figure 3A). Using the Rap1 ChIP assay (Figure 2A) it was also evident that both rif1Δ and rif2Δ mutants had structural defects in terms of folding back into the subtelomere (Figure 3B). As with the above described Rap1 ChIP, we confirmed that Rap1 protein levels were not altered in rif1Δ and rif2Δ cells which may have accounted for observed differences (Figure S3A). From here on we have performed further analysis only with the rif2Δ mutant rather that rif1Δ cells due to the fact that apart from telomere length regulation, Rif1 also plays an important role in telomere capping , checkpoint regulation – as well as telomere localization , which greatly complicated the interpretation of rif1Δ cells and their genetic interactions. Ongoing studies are directed at better understand the contributions of Rif1 in promoting the telomere fold-back structure.
We constructed double mutants between rif2Δ and mutants of the Rpd3L, Rpd3S and Hda1 complexes harboring construct 2 in order to assess potential genetic interactions between different pathways involved in telomere looping. Whereas the Rpd3S specific mutants (eaf3Δ and rco1Δ) displayed a slight additive growth advantage on 5-FOA when combined with rif2Δ mutants compared to the respective single mutants (Figure 3C, bottom panels) there were no additive effects with rif2Δ mutants and the Rpd3L complex (sap30Δ and rxt2Δ) (Figure 3C, top panel). sin3Δ rif2Δ double mutants were not additive in comparison to the respective single mutants (Figure 3C, middle panel), as Sin3 belongs to both L and S complexes. As expected, rxt2Δ (Rpd3L) and eaf3Δ (Rpd3S) double mutants had a slight additive looping defect in comparison to the single mutants and further deletion of RIF2 did not exacerbate this defect (Figure 3D). The looping defect of hda1Δ and hda2Δ mutants was also additive with rif2Δ mutants (Figure 3E). The results of this genetic epistasis analysis suggest that Rif2 and Rpd3L may function together in a common pathway to promote a telomere fold-back.
To mechanistically understand the genetic relationships between the Rpd3L complex and Rif2, we performed ChIP experiments to determine if the rif2Δ mutation had an effect on Rpd3L (Rxt2-TAP) localization at telomeres. Subtelomeric DNA was enriched above non-tagged wild type control cells (background) with an epitope-tagged Rxt2 allele both close to (-6 bp) (Figure 3F) and up to 2000 base pairs away from the telomeric tracts (Figure 3F) as previously described . This enrichment was decreased to near background levels in a rif2Δ mutant (Figure 3F). The loss of ChIP signal is not due to altered expression levels of Rxt2-TAP in rif2Δ mutants as confirmed by western blot analysis (Figure S3B). Unlike Rpd3L, the ability to cross-link the Rpd3S complex (Rco1-TAP) to subtelomeric regions was not altered in rif2Δ cells (Figure 3G and Figure S3C). Together these data suggest that the Rif2 promotes a structural alteration at telomeres through the recruitment of the Rpd3L KDAC complex. Moreover, the Hda1 KDAC as well as the Rpd3S complex promote the same fold-back structure, but independent of the Rif2/Rpd3L pathway.
To better understand the function of the loop structure and whether or not it may have a protective role at the telomere we impaired looping (deletion of SIN3) in various genetic backgrounds where telomere function was compromised. In cdc13-1 sin3Δ double mutants we observed a temperature-dependent synthetic lethality in the double mutant, compared to the respective single mutants (Figure 4A), indicating that partially uncapped telomeres (cdc13-1) may require a Rpd3-mediated structure for viability. The negative genetic interaction was suppressed by the further deletion of EXO1, the nuclease responsible for the majority of telomere resection at dysfunctional telomeres (Figure 4A). To ensure that this interaction was a direct consequence of telomere dysfunction, we assayed the accumulation of telomeric ssDNA following the shift of nocodazole-arrested cdc13-1 and cdc13-1 sin3Δ cells from 23°C to the semi-permissive temperature of 26°C (Figure 4B). In agreement with the negative genetic interaction (Figure 4A), we observed an Exo1 dependent increase in telomeric ssDNA in the double mutant above that seen in the cdc13-1 single mutant (Figure 4B). Importantly, we did not detect an increase in ssDNA at telomeres in sin3Δ single mutants when compared to isogenic wild type control cells (Figure S4A).
To determine if Rpd3 dependent telomere structure may have an influence on the rate of cellular senescence, we compared senescence onset in both est2Δ rad52Δ (where HR and telomerase-mediated telomere elongation are impaired) and est2Δ rad52Δ sin3Δ mutants. The deletion of SIN3 resulted in a dramatic increase in the rate of cellular senescence in the absence of telomerase and homologous recombination (Figure 4C). The accelerated loss of viability associated with the sin3Δ mutation is specifically related to the absence of telomerase as rad52Δ cells maintain viability to a similar extent as rad52Δ sin3Δ cells (Figure 4C). To better understand the cause of premature senescence in sin3Δ mutants, genomic DNA was prepared from the senescence curves (Figure 4C) and both single stranded telomeric DNA accumulation and telomere length were analyzed. As is the case when combined with the cdc13-1 mutant, we found that sin3Δ est2Δ rad52Δ mutants had increased levels of ssDNA at telomeres compared to isogenic est2Δ rad52Δ cells (Figure 4D). The telomere shortening rate is unaffected between sin3Δ est2Δ rad52Δ and est2Δ rad52Δ strains (Figure S4B, S4C), however there is evidence of early rapid telomere loss events at some telomeres in the triple mutant compared to the isogenic double (Figure S4B).
Taken together these results imply that the Rpd3 lysine deacetylase is essential to prevent excessive nuclease-mediated resection specifically at uncapped telomeres. In the absence of a telomere lengthening mechanism this resection may lead to excessive telomere shortening. We propose that this protective function may involve the formation of a fold-back structure at the chromosome ends.
We have demonstrated that multiple biological processes influence the ability of telomeres to form a higher-order fold-back structure. High-throughput screening coupled to stringent bioinformatic analysis, has revealed that class I (Rpd3) and class II (Hda1) KDAC activities were among the most significantly enriched biological processes required to promote the formation/maintenance of a telomere loop. In addition, the Rap1 binding proteins, Rif1 and Rif2, which localize directly to telomeres, were implicated in telomere fold-back establishment. Through genetic epistasis analysis, we found that rif2Δ and mutants of the Rpd3L complex did not have additive structural defects at telomeres, suggesting that they may function together in a single pathway. This genetic interaction was confirmed by demonstrating that Rpd3L was no longer able to localize to telomeres in rif2Δ cells. Furthermore, we have shown that the telomeres in mutants with looping defects are more susceptible to uncapping, nucleolytic degradation, telomere loss and promote accelerated rates of cellular senescence in the absence of telomere maintenance. A recent study has reported that the Rpd3 complex is required to prevent chromosome end fusions in Drosophila melanogaster . Together, these data suggest that the regulation of such telomere fold-back structures may be conserved, and furthermore indicate that multiple cellular processes, apart from those that directly impinge on DNA metabolism, may have effects on telomere structure/function.
Although we have focused on the telomere dysfunction phenotypes associated rpd3 mutants, it is of interest that many of the other mutants recovered in our telomere structure screen have been previously identified to have negative synthetic interactions with cdc13-1 (rrd1Δ, swr1Δ, arp6Δ, spe1Δ, spe3Δ, oca1Δ, oca2Δ, oca5Δ, elp4Δ, pep8Δ, rif1Δ, yme1Δ, htd2Δ, ski3Δ, rim101Δ, ede1Δ) , – as well as an increased rate of replicative senescence (hda1Δ, sin3Δ, sec28Δ, med1Δ, asf1Δ, rif1Δ, rif2Δ, arp6Δ, elp4Δ) – (see Table 1 for a complete overview of overlaps between our screen and other selected telomere function screens). Moreover, the Sir2/3/4 complex, which is required to form a telomere fold-back structure  also prevents premature senescence . This overlap between our screen and previously published data suggests that telomere structure may make significant contributions towards preserving telomere integrity when telomere function is compromised (Figure 5). Consistently, sin3Δ single mutants do not exhibit increased ssDNA accumulation at telomeres nor do they have any changes in telomere length in comparison to isogenic wild type strains. This would suggest that in the absence of a telomere fold-back, the CST complex and other capping factors are sufficient to maintain a protected state (Figure 5). However, when CST function is compromised (e.g. cdc13-1) in combination with an inability to fold-back, resection becomes accelerated (Figure 4B, ,5).5). In terms of telomere-induced cellular senescence in the absence of telomere maintenance, the increased resection in non-looped mutants would also lead to increased telomere shortening (Figure 5). It will be of interest to perform an extensive genetic epistasis of all mutants isolated in the loop screen in order to determine if the negative interactions with cdc13-1 are epistatic (i.e. due to a fold-back defect).
There was also a large overlap between mutants found in our screen and mutants that have been implicated in both the positive and negative regulation of telomere length (vps28Δ, trk1Δ, ctk1Δ, mrt4Δ, rif1Δ, sur4Δ, nut1Δ, siw14Δ, leo1Δ, hit1Δ, pcp1Δ, pdx3Δ) –. This would suggest that either telomere looping has a direct effect on telomere length regulation, or conversely, may indicate that telomere length changes impinge on the ability to form a fold-back. Our results suggest that telomere looping does not affect telomere length homeostasis directly, as many of the mutants recovered in our screen, even those with “strong” looping defects, have wild type telomere length. On the other hand, telomere length changes could indeed have a drastic effect on telomere structure in terms of the chromatin alterations that occur with respect to length changes. Telomere shortening, for example, results in de-silencing in the subtelomeric region  due to the decreased capacity of shortened telomeres to recruit the Sir2/3/4 histone deacetylase complex, which is required for loop formation in yeast . Long telomeres, in contrast, promote a hyper-silenced state in the subtelomere , much like what occurs in rif1Δ and rif2Δ mutants, where, in the case of the latter mutant, Rpd3L fails to localize to telomeres. Indeed, Rpd3 mutants (L and S) are hypersilenced in the subtelomeric zone. One possibility would be that long telomeres (as seen in rif2Δ) mutants fail to properly localize Rpd3L, which leads to a subsequent fold-back defect. Although our screen implicates proper length regulation as a key regulator of telomere looping, they remain correlative and require further investigation in order to draw concrete conclusions.
Whereas both the Rpd3S and Hda1 complexes were additive with rif2Δ mutants in terms of a looping defect, the Rpd3L complex was epistatic. This relationship was confirmed mechanistically as we noticed that Rpd3L is not able to properly localize to telomeres in rif2Δ cells. These epistasis analyses revealed that multiple KDACs contribute to telomere looping (Figure 5). Rpd3L, Rpd3S and Hda1 all promote telomere looping in parallel pathways whereas the relationship between the Sir2/3/4 complex and the other KDACs in terms of telomere structure remains enigmatic. The KDACs are best known for their deacetylation of histones and they are known to contribute significantly to silencing at subtelomeric loci . Consistent with a connection between chromatin modification and the telomere fold-back structure, we also found that many members of the Swr1 chromatin remodeling complex as well as the histone chaperone Asf1 (Table 1) were important for telomere folding. A future challenge will be to determine the targets of the KDACs. Although subtelomeric histones are prime candidates, telomeric proteins themselves may be targets. Furthermore, it will be important to characterize how the other mutants revealed in the screen contribute to telomere looping and to understand if these mutants are epistatic with the KDACs. Indeed, it has been difficult to understand why mutations that affect such diverse biological pathways may have synthetic growth defects with cdc13-1 or in some cases, senesce rapidly –, . Since many signaling pathways activate effectors via activation/repression of target genes through chromatin remodeling and/or histone acetylation/deacetylation, we propose that activation or repression of these pathways may influence the ability of the KDACS (Rpd3/Hda1 and Sir2) to act at telomeres and in turn directly or indirectly influence telomere structure.
This work has uncovered multiple regulators of the telomere fold-back structure, including lysine deactylation and chromatin remodeling. The results of our screen correlate well with screens that have been performed to elucidate genes implicated in telomere function and cellular senescence suggesting that the fold-back structure may be important for chromosome end protection. Previous models have speculated that the fold-back in yeast may be important to establish silent chromatin within the telomeric/subtelomeric loci –. As an alternative, it was also suggested that silent telomeric chromatin may be required to establish a particular architecture that contributes to chromosome end protection . Our results indicate that silent chromatin can be established in the absence of a telomere fold-back since many of the mutants recovered in our screen are not compromised for silencing (Table 1) or even have slightly enhanced silencing (e.g. rpd3 mutants). Sir2/3/4-mediated silencing as well as other chromatin modifications are however important to establish a telomere loop, which likely promotes end protection. Interestingly, Rpd3 dependent histone deacetylation has been shown to prevent Sir2/3/4 protein spreading towards the centromere , which may potentially deplete SIR protein levels immediately adjacent to the telomere, raising the possibility that SIR2/3/4 disruption and RPD3L/S disruption may be one in the same in terms of a fold-back defect. Further characterization of the yeast fold-back structure and the relationship to silent chromatin will be essential in order to clarify these issues. In summary, by understanding how different biological processes impinge on chromosome end structure, we increase the possibilities to manipulate telomere function, both positively and negatively, which may have important implications for diseases that stem from telomere dysfunction.
For spotting assays, yeast cells were incubated overnight at appropriate temperature in YPD. Cells were diluted to OD600 0.5 and spotted in ten-fold dilutions onto 2% raffinose, 1% galactose plates either with (+FOA) or without (−FOA) 5-FOA. Cells were incubated for 3 days at proper temperature, imaged and then replica plated on SD-URA plates for 2 days at the same temperature before imaging.
Telomere PCR was performed using 100 ng genomic DNA diluted in 1× NEB4 buffer and water. Samples were denatured for 10 min at 96°C and cooled to 4°C. Tailing mix (4 U/µl terminal transferase (NEB), 1× NEB4 buffer, 1 mM CTPs) was added to a final concentration of 10%. Tailing reaction was performed as the follows: 37°C 30 min, 65°C 10 min, 96°C 5 min, arrest at 65°C. 3× volume of preheated PCR-MIX (1 µM oligo dG reverse primer, 1 µM telomere specific forward primer either 1L, 6R, 7L or Y′, 0.267 mM dNTPs, 0.083 U/µl Phusion polymerase (NEB), PCR buffer (89.11 mM Tris-HCl pH 8.8, 21.28 mM (NH4)2SO4, 6.65% glycerol, 0.0133% Tween-20) was added and PCR reaction was performed using: 95°C 3 min, 45 cycles: (95°C 30 s, 68°C 15 s, 72°C 20 s), 68°C 5 min, hold on 12°C.
Samples were mixed with DNA loading buffer and separated on a 1.8% agarose gel for 30 min at 100 V. Bands were detected using LAS-4000 (Fujifilm) and quantified using Multi Gauge Software (Fujifilm). A complete lis of oligonucleotides used in this study can be found in Table S2.
Spore-colonies of dissected heterozygous diploids were suspended in water and diluted in 5 ml YPD medium to a final concentration of OD600 0.01. Cells were incubated for 24 h at 30°C and absorption at 600 nm was measured. Cultures were re-diluted to OD600 0.01 in 5 ml YPD and inoculated for further 24 h at 30°C. Each day cell samples were harvested and genomic DNA was prepared for telomere length analysis (Quiagen genomic DNA prep. Kit). Population doublings (PD) were calculated as log2(OD60024 h/0.01). All PD values refer to PD after the spore colony had been harvested from the dissection plate (about 25 generations). Graphs were made in Prism5 (GraphPad).
Synthetic Genetic Array (SGA) methodology was used (strains R1459 and R1460) to obtain haploid gene deletion mutants containing either construct 1A or construct 2 (Tong and Boone, 2006). Cells were then replica-pinned onto media containing galactose, with and without 5-FOA. Growth on 5-FOA was compared between construct 2-containing deletion mutants and construct 1-containing mutants (comparison 1). Growth of construct 2-containing mutants was also compared with and without the presence of 5-FOA (comparison 2). Construct 2-containing mutants that grew better by either comparison 1 or 2 were selected for validation. Validation was carried out by manually crossing and dissecting tetrads from independent starter strains followed by duplicate spot assays onto media with and without 5-FOA.
Yeast cells were grown over night at 30°C and diluted to OD600 0.2. They were grown until exp. Phase (OD600 0.6–1.0), crosslinked for 8 min (20 min for Rxt2-TAP-ChIP, 10 min for Cdc13-TAP-Chip) with formaldehyde (final conc. 1.2%) and quenched with glycine (360 mM final). After adjusting the volume to the same OD all samples were washed two times with 1× PBS, resuspended in FA Lysis Buffer (50 mM HEPES-KOH pH 7.5, 140 mM NaCl, 1 mM EDTA pH 8, 1% Triton X-100, 0.1% Sodium deoxycholate) and lysed with Matrix C tubes via FastPrep (6.5 M/s, 2°—30 sec with 1 min break). Cell extracts were recovered, centrifuged and the soluble potion of the lysate was discarded. Pellets were resuspended in FA buffer +SDS (2% final) and split up for sonication. Chromatin was sheared 30 sec on/off for 15 min. Supernantant (ChIP extract) was diluted to 1 mg/ml protein concentration in FA buffer and used for immunoprecipitation (IP).
Pre incubated protein G sepharose beads (washed with 1×PBS, FA Buffer and pre-incubated with 5% BSA for 1 h at 4°C) were added to the 1 mg/ml solution to perform an addition precleaning step before the IP (1 h at 4°C). After precleaning anti-Rap1 antibody (Santa Cruz) was added to the solution (1100) and incubated with fresh beads over night at 4°C, rotating.
For Tap-ChIPs (Rxt2 and Cdc13) IgG-Sepharose Beads (washed with 1× PBS and FA-buffer) were added to the 1 mg/ml solution and IP was incubated over night at 4°C.
Sonication efficiency was tested via cleaning 100 µl of the ChIP extract and performing agarose gelelectrophoresis. IP was washed with FA-Lysis buffer, FALysis buffer 500 (FA buffer with 500 mM NaCl), Buffer3 (10 mM Tris-HCl pH 8 1 mM EDTA pH 8, 250 mM LiCl, 1%NP-40, 1% Sodium deoxycholate), and TE (pH 8). For elution buffer B (50 mM Tris-HCl ph 7.5, 1% SDS, 10 mM EDTA pH 8) was added and IP was incubated at 65°C for 8 mins. For reverse-crosslinking proteinase K was added to the IP and INPUT control (ChIP extract, 1 mg/ml solution without IP) and incubated at 65°C, rotating overnight. Samples were cleaned with Quiagen “QIAquick PCR Purification Kit” and qPCR analysis was performed using Roche standard PCR protocol for Sybr-Green detection with 55°C annealing temperature, all oligonucleotides used are listen in Table S2. Measured ct values were corrected to INPUT and normalized to the actin signal using the following formulas.
Cells were grown overnight at 23°C in 10 ml YPD. Saturated cultures were diluted to OD600 0.2 in 150 ml YPD and incubated at 23°C until they reach log phase (0.6–0.8). Nocodazol (20 µg/ml final) was added and cells were incubated for a further 3 h at 23°C, shaking. Cells were checked under the microscope until >90% were largebudded. “Pre” samples were harvested and cells were subsequently shifted to 26°C.
Additional samples were collected for all time points (30 min, 60 min, 90 min and 120 min) after the shift. For ssDNA analysis, dot blotting was performed. DNA was extracted using genomic DNA Kit (Quiagen). Isolated DNA was either denatured using 0.2 M NaOH and 65°C for 15 min or kept on ice for native conditions.
For blotting, 4 µg DNA (native) or 0.5 µg (denatured) were suspended in 200 µl 2×SSC and loaded to the dot blot apparatus using nylon membrane (GE Healthcare Amersham H-bond). After crosslinking (UV Stratalinker 2400, Stratagene) DIG labeling (DIG labeled probe oBL207) and detection was performed as described by the product guidelines (Roche DIG oligonucleotide 3′labeling KIT).
Cells were grown overnight at 30°C in 5 ml SD medium containing 2% Raffinose (S-Raf). Saturated cultures were diluted to OD600 0.2 in 8 ml S-Raf and incubated at 30°C until they reach log phase (0.6–0.8). Cells were split and 2% galactose or 2% glucose (final) were added. Cells were incubated for 2 1/2 h at 30°C, shaking. Cells were centrifuged down and RNA was extracted and Northern Blotting was performed as described previously . URA3 and actin were detected using DiG labeled PCR products (Roche, “DiG High Prime” labeling) gained from a PCR reaction with oBL17, oBL18 for URA3 and oBL292, oBL293 for actin (1. 98°C 30 sec, 2. 98°C 10 sec, 3. 60°C 30 sec, 4. 72°C 1 min, 5. 72°C 5 min, 12°C forever, repeating steps 2–4 for 33 cycles). Quantification was performed using Multi Gauge software (Fujiifilm) and signal was displayed as URA3 over actin.
3 ml of culture (log phase) was centrifuged at 13000 rpm for 2 min. Pellets were resuspended with 150 µl solution 1 (0,97 M β-mercaptoethanol, 1,8 M NaOH) and incubated on ice for 10 min. 150 µl, 50% TCA was added and cells were incubated 10 min on ice, centrifuged at 13000 rpm for 2 min at 4°C and the pellet was resuspended with 1 ml acetone. Solution was centrifuged at 13000 rpm for 2 min at 4°C and the pellet was resuspended in 140 µl UREA buffer (120 mM Tris-HCL pH 6,8, 5% Glycerol final, 8 M Urea final, 143 M 2-mercaptoethanol final, 8% SDS final, a little bit of bromphenol blue indicator). Protein extract was incubated 5 min at 95°C, centrifuged and loaded on a pre-cast gradient Gel (BioRad).
A) Examples of the different qualitative assignments given to looping mutants from the screen where +++=strong defect, ++=moderate and +=weak, the mutant ORFs are indicated. B), C) The indicated mutants were grown for over 100 generations and subject to telomere PCR at telomere 1L and at Y′ telomeres. Telomere lengths are represented with respect to wild type (which is set to 0). rif2Δ and est2Δ mutants serve as controls for mutants with long and short telomeres, respectively. D) Northern blot quantification of URA3 expression in strain harboring construct 2. Strains were either grown in the presents of glucose (uninduced) or galactose (induced). RNA was isolated and northern blot bands were quantified using Multi Gauge software (FujiiFilm). Signal was calculated as URA3 RNA over actin RNA.
Sonication control for the Rap1-ChIP; RAP1 and Cdc13-TAP protein expression. A) DNA-sonication control of the chromatin fragments for Rap1-ChIP (hda1Δ and sin3Δ). All fragments have an average bulk size of 300 bp. B) Western blot for Rap1 total protein level in hda1Δ and sin3Δ or sir4Δ mutants reveals no differences in general RAP1 expression. Rap1-TAP causes a shift in the Rap1 protein band due to the TAP-tag. C) Western blot for Cdc13-TAP protein level in wild type or sin3Δ cells. Tap-tag was detected via peroxidase-anti-peroxidase antibody. D and E) Telomere PCR at telomere 1L revealed a shortening upon EST2 disruption following 25 and 50 generations of growth. F) Enrichment of the -6 position at natural telomere 6R is increased upon telomere shortening following ChIP directed against Cdc13-TAP.
Rap1, Rxt2-TAP- and Rco1-TAP protein expression. A), B) and C) Western blot analysis of Rap1, Rxt2-TAP and Rco1-TAP protein levels. Protein was extracted from cells with the indicated genotypes and western blotting was performed with either an anti-Rap1 antibody (A) or a mouse, anti-actin antibody which through its Fc domain cross reacts with the TAP epitope tag on Rxt2-TAP and Rco1-TAP respectively (B). Expression levels of Rap1, Rxt2-TAP and Rco1-TAP do not differ when comparing a wild type background with a rif2Δ (and rif1Δ) background.
ssDNA analysis of sin3Δ mutants and telomere length analysis for senescence curves on telomere 6R, 7L and Y′telomeres. A) Cells of the indicated genotype were arrested in nocodazole for 2.5 hours before genomic DNA was extracted and dot blotted in native and denaturing conditions. Quantification of a dot blot analysis where the amount of ss telomeric DNA is normalized to total telomeric DNA. B) Telomere PCR for telomere 6R, and Y′ telomeres to follow the rate of telomere shortening from the survivor curves described in Figure 4C.
Yeast strains used in this study. All strains used in this study were derived from BY4741 background (his3-1, leu2-0, ura3-0, met1-0) unless indicated otherwise.
Oligonucleotides used in this study. All oligos used in this study are listed and are in the 5′to 3′ (left to right) direction.
We thank S. Marcand, V. Zakian, A. Lustig, and M. Ptashne for strains and reagents and Sarah-Luke Glaser for critical reading of the manuscript.