- Open Access
Spatial reorganization of telomeres in long-lived quiescent cells
Genome Biologyvolume 16, Article number: 206 (2015)
The spatiotemporal behavior of chromatin is an important control mechanism of genomic function. Studies in Saccharomyces cerevisiae have broadly contributed to demonstrate the functional importance of nuclear organization. Although in the wild yeast survival depends on their ability to withstand adverse conditions, most of these studies were conducted on cells undergoing exponential growth. In these conditions, as in most eukaryotic cells, silent chromatin that is mainly found at the 32 telomeres accumulates at the nuclear envelope, forming three to five foci.
Here, combining live microscopy, DNA FISH and chromosome conformation capture (HiC) techniques, we report that chromosomes adopt distinct organizations according to the metabolic status of the cell. In particular, following carbon source exhaustion the genome of long-lived quiescent cells undergoes a major spatial re-organization driven by the grouping of telomeres into a unique focus or hypercluster localized in the center of the nucleus. This change in genome conformation is specific to quiescent cells able to sustain long-term viability. We further show that reactive oxygen species produced by mitochondrial activity during respiration commit the cell to form a hypercluster upon starvation. Importantly, deleting the gene encoding telomere associated silencing factor SIR3 abolishes telomere grouping and decreases longevity, a defect that is rescued by expressing a silencing defective SIR3 allele competent for hypercluster formation.
Our data show that mitochondrial activity primes cells to group their telomeres into a hypercluster upon starvation, reshaping the genome architecture into a conformation that may contribute to maintain longevity of quiescent cells.
The spatiotemporal behavior of genomes and their regulatory proteins is an important control mechanism of genomic function. One of the most pervasive features of nuclear organization is the existence of subnuclear compartments, which are thought to create microenvironments that favor or impede specific DNA- or RNA-related processes . Deciphering how the dynamics of this subnuclear compartmentalization are regulated in relation to changes in genome activity is a key step in understanding how nuclear organization participates in nuclear function.
Well-characterized examples of subnuclear compartments include clusters of specific genes or repetitive DNA sequences , such as telomeric repeats (in budding yeast) or centromeric satellites (in fission yeast, fly and mammals) and retrotransposons (in fission yeast, Tn2/Ku70-mediated clustering) . These repetitive sequences generally nucleate patterns of histone modifications that are recognized by histone-binding repressors, and their clustering results in the sequestration of these general repressors into subcompartments. Besides its role in concentrating silencing factors, this evolutionarily conserved phenomenon has a dominant impact on chromosome folding and positioning. In metazoans, a cell type-specific nuclear distribution of heterochromatin is established upon cell differentiation, and is often compromised in cancer cells . In budding yeast, the clustering of silent chromatin provides an excellent model of a subnuclear compartment.
Most Saccharomyces cerevisiae functional and structural studies have been conducted on exponentially growing cell cultures. In these conditions, silent chromatin is mainly found at telomeres and at the cryptic mating type loci (HM loci), where it is generated by the recruitment of the SIR complex comprising Sir2, Sir3, and Sir4. At telomeres, this nucleation event is achieved by the transcription factor Rap1, which binds the telomere TG repeats and interacts with Sir3 and Sir4. Sir4 heterodimerizes with the NAD + −dependent histone deacetylase Sir2, which deacetylates H4 histone tails from neighboring nucleosomes, thus generating binding sites for Sir3. The SIR complex thus spreads over a 2–3-kb subtelomeric region leading to the transcriptional repression of subtelomeric regions.
The clustering of telomeres into perinuclear foci generates a zone that favors SIR-mediated repression at the nuclear periphery [5, 6] and ensures that SIR proteins do not bind promiscuously to repress other sites in the genome [7, 8]. Furthermore, telomere anchorage in S phase contributes to proper telomerase control and suppresses recombination among telomere repeats [9, 10].
The average large-scale organization of budding yeast chromosomes during exponential growth has been described through genome-wide capture of chromosome conformation (3C) experiments . This analysis unveiled a polarized configuration with chromosome arms extending away from the centromeres that are held by the spindle-pole body, in agreement with microscopy data . This so called Rabl organization — initially observed by Carl Rabl in rapidly dividing nuclei of salamanders  — can be mimicked to some extent by polymer models using a limited number of assumptions [11, 14–16]. However, it remains unclear how specific biological processes could affect this robust average organization.
As mentioned above, most of the studies characterizing genome organization and its functional consequences in budding yeast have been conducted in nutrient-replete conditions with cells undergoing exponential growth. However, yeast cells rarely experience such a lush environment and their survival in the wild depends on their ability to withstand adverse conditions.
It is well known that yeast cells finely tune their growth and behavior to their environment, adapting to nutritional depletion or stresses by engaging specific developmental programs . When grown in rich media containing glucose, they progress through distinct metabolic programs (Fig. 1a), with each transition being accompanied by widespread transcriptional reprogramming [18, 19]. In the first phase (exponential phase), cells metabolize glucose predominantly by glycolysis, releasing ethanol in the medium. When glucose becomes limiting, the culture enters diauxic shift, a transition characterized by a decreased growth rate and a metabolic switch from glycolysis to aerobic utilization of ethanol. Finally, when no other carbon source is available cells enter stationary phase (SP). During that stage most cells are in quiescence, a non-proliferative state that maintains the ability to resume growth following restoration of missing nutrients.
Recent studies in different species demonstrated that a hostile environment (i.e., caloric restriction or the presence of mild oxygen stresses) can trigger a “vaccination-like” adaptive response leading to the acquisition of anti-aging functions . Following the same principle, budding yeast can reach different quiescent states depending on the conditions that induce the cell cycle exit, each of them leading to different outcome in terms of chronological lifespan (CLS) . Deciphering the key features that differentiate each metabolic state is essential to understand mechanisms that extend lifespan in yeast.
Here we show that, following carbon source exhaustion, the silencing factor Sir3 drives the telomeres of quiescent cells to group together, forming a discrete, large cluster (hypercluster) at the center of the nucleus. This organization is specific to quiescent cells able to sustain long-term viability. Our data strongly support a model in which mitochondrial activity, through the production of reactive oxygen species (ROS) during cell respiration, commits cells to form a telomere hypercluster upon starvation. Importantly, sir3∆ cultures, which are defective in forming telomere hyperclusters in SP, show reduced CLS. Furthermore, expressing a silencing-defective SIR3 allele rescues both telomere distribution and the CLS of a sir3 null strain, strongly arguing that telomere clustering directly contributes to cell survival during quiescence.
Massive telomere reorganization upon carbon source exhaustion
To investigate telomere organization in live cells, we monitored the subnuclear distribution of the telomeric protein Rap1 fused to green fluorescent protein (GFP)  at different stages of a liquid culture, from glycolysis to late respiration to SP. We observed dramatic changes in the distribution of the Rap1-GFP signal during this time course (Fig. 1a, b). In agreement with previous reports [6, 22], Rap1-GFP formed three to five foci during the logarithmic phase, quantified using our custom-made software (Fig. 1c; adapted from ). In cells undergoing respiration (after 2 days in culture), Rap1-GFP foci were fewer and brighter, with 50 % of the cells showing two foci and 30 % of the cells having only one focus (versus 23 % and 10 %, respectively, during fermentation). In SP 62 % of the cells exhibited a unique focus with a median intensity that was fivefold higher than in the exponential phase. Moreover, we noticed that when the number of foci per cell decreases, the intensities of the remaining foci increase (Fig. 1c), suggesting that smaller foci group into larger ones. Importantly, we verified that the brightness of the Rap1-GFP clusters observed in SP was not due to an overall increase in Rap1-GFP levels (Additional file 1: Figure S1a). Furthermore, we observed a similar clustering with SIR complex proteins fused to GFP (Sir2/3/4; Additional file 1: Figure S1b). We confirmed that Rap1-GFP foci coincided with the Y’ telomeric clusters and Sir3 foci in SP cells by combined immunostaining and fluorescence in situ hybridization (immuno-FISH; Fig. 1d) and in vivo imaging (Additional file 1: Figure S1c). Thus, telomere-associated silent chromatin groups into “hyperclusters” in SP cells.
As in exponentially growing cells, telomere hyperclustering requires SIR3 and SIR4 in SP cells (Fig. 1e; Additional file 1: Figure S1d). Although the brightest Rap1-GFP focus was mainly found adjacent to the nuclear envelope in exponentially growing cells [6, 22], telomere hyperclusters were overwhelmingly found in the innermost area in SP cells (>90 % in zone 3; Fig. 1f).
We next evaluated whether other nuclear landmarks were also altered in SP. In agreement with previous reports, we found that the nuclear diameter (data not shown, inferred from experiment Fig. 1f) was smaller and the nucleolus more compact in cells after the diauxic shift (Additional file 1: Figure S1e) . Consistent with , we noticed that kinetochore proteins formed a “bundle” in a subpopulation of cells; however, this structure did not correlate with telomere hyperclusters (Additional file 1: Figure S1f). Furthermore, we did not observe major changes in the distribution of the centromere-associated protein Cse4 in SP cells containing telomere hyperclusters (Additional file 1: Figure S1g). Thus, a specific SIR-dependent re-organization of telomeres occurs in a subpopulation of SP cells.
Hyperclustering of telomeres occurs only in the long-lived fraction of SP cells
As previously reported , SP cultures consist of different types of cells. Equilibrium density-gradient centrifugation enables the separation of a dense fraction mainly composed of small unbudded daughter cells that are able to sustain long-term viability, and a lighter fraction that includes both budded and unbudded cells that rapidly lose the ability to perpetuate over time. Calcofluor staining revealed that cells with hyperclusters (defined as cells containing one or two foci and at least one Rap1-GFP focus with intensity levels above 95 % of foci in exponentially growing cells) are essentially small unbudded cells (Fig. 2a). Sorting SP cells by density gradient enriched the population of cells showing hyperclusters from 69 % to 84 % in the densest fraction (HD) while most cells from the less dense fraction (LD) showed a distribution of Rap1-GFP foci similar to the post-diauxic shift cells (Figs. 1b, c and 2b, c). Moreover, we confirmed that the viability is significantly lower for the lighter cells than for the denser ones that show hyperclusters (37 % versus 99 %, respectively). We thus conclude that telomere hyperclustering occurs specifically in quiescent SP cells.
The global chromosome organization in long-lived SP cells is constrained by centromere and telomere clustering
To decipher the three-dimensional (3D) organization of the entire genome in long-lived SP cells, we turned to 3C . We used an untagged strain to avoid any possible artifact related to the expression of tagged telomere proteins. As cells from the dense fraction of SP are small unbudded cells (Fig. 2a), we compared the genomic contact maps of these cells with G1 daughter cells elutriated from an exponential culture to avoid the contribution of the cell cycle in this latter case. In order to facilitate the interpretation of the contact map [Fig. 3a(ii)], the matrix was converted into a 3D map in which the distance between each pair of genome segments was optimized to reach the inverse of their measured contact frequency [Fig. 3a(i); Additional file 2) . This 3D reconstruction of the entire contact map provided a remarkable overview of the average yeast genome organization in a population of G1 cells, with the rDNA clearly isolated from the rest of the genome, a dense centromere cluster, and a tendency for subtelomeres to co-localize, consistent with the well-documented perinuclear clustering of telomeres .
In excellent agreement with our microscopy data, contacts between telomeres became prominent in quiescent cells, generating a remarkable hypercluster. The influence of chromosome arm length on the subtelomere contacts — which in exponentially growing cells discriminates two groups of telomeres exhibiting preferential contacts with each other’s — is alleviated by the formation of the hypercluster, suggesting the formation of a grid-like/homogeneous disposition of telomeres (Fig. 3b, c) [11, 28]. In addition, regions closer to the telomeres exhibited an increased number of contacts in SP, whereas the number of contacts between centromeres decreased slightly (Fig. 3d; Additional file 3: Figure S2a). Thus, the frequency of contacts increases specifically between telomeres, imposing a general constraint on the whole genome organization, with each chromosome arm now being tethered at two points of the nuclear space (Fig. 3a(iii); Additional file 4). As a result, the average contacts between chromosome arms, which are primarily constrained by their sizes and centromere clustering in G1 (Fig. 3a(i), c), appear distorted due to subtelomere interactions in G0 (Fig. 3a(iii), d). Importantly, these observations were confirmed in two different genetic backgrounds (BY and W303; Additional file 3: Figure S2c, d).
To test whether this reorganization is driven by increased telomere–telomere interactions, we compared the genomic contact map of cells in which SIR3 had been deleted and wild-type cells from the dense fraction of a SP culture. In agreement with our microscopy data (Fig. 1e) we observed that sir3∆ cells were not able to generate a hypercluster upon entry into quiescence (Fig. 3e, f; Additional file 3: Figure S2b) and that the general organization of chromosomes in sir3∆ quiescent cells resembles the organization of wild-type G1 cells, with similar levels of contacts between subtelomeric regions (Fig. 3d, f; Additional file 3: Figure S2b). We thus conclude that the main changes in chromosome organization that occur as cells enter quiescence are driven by an increase in Sir3-dependent telomere clustering.
Telomeres form hyperclusters specifically in conditions inducing long-lived quiescent cells
To test whether telomere hyperclusters were a general feature of quiescence we compared telomere subnuclear distribution in quiescent cells induced by different means. As mentioned above, although quiescent cells are by definition viable, their CLS properties depend on the method/metabolic changes used to induce the cell cycle exit  (Fig. 4a). At day 7 of CLS 61 % of quiescent cells arising from progressive carbon source exhaustion (SP) had formed telomere hyperclusters and these retained >90 % viability (Fig. 4a, b). In contrast, quiescent cells induced by nitrogen starvation formed hyperclusters at a much lower rate (18 % had done so) and lost viability more rapidly, as previously reported . Hence, the grouping of telomere foci into hyperclusters is not a consequence of cell cycle arrest but rather a specific feature of long-lived quiescent cells induced by carbon source exhaustion.
The ability to form telomere hyperclusters upon starvation is acquired during respiration
Interestingly, when abruptly starved from carbon source, cells respond differently depending on their initial metabolic status: few cells previously undergoing glucose fermentation formed telomere hyperclusters upon starvation (7 %) and showed a strong decrease in viability at day 7 (≈40 %), in agreement with previous reports [21, 29]. In contrast, 73 % of cells previously undergoing respiration (post-diauxic shift) formed telomere hyperclusters upon starvation and these retained ≈ 90 % of viability at day 7. Thus, only cells that experienced respiration before entering quiescence had a long CLS (>90 % viability after 1 week of starvation) and formed telomere hyperclusters at rates of more than 60 % (Fig. 4a, b). These characteristics could be attributed either to their metabolic activity or to their growth rates, as cells undergoing respiration divide slower and slow growth confers resistance to various stresses . However, slow growth was not sufficient to prime cells to form a hypercluster upon starvation, as cells grown slowly in glucose at 25 °C and starved after fermentation did not form hyperclusters (Additional file 5). To determine if respiration was an obligatory step to induce telomere hyperclustering upon starvation, we monitored telomere clustering in respiratory deficient cells (rho-) after glucose exhaustion (Fig. 4b) or upon abrupt starvation (data not shown). These conditions led to a very low rate of cells with bright Rap1-GFP foci (3 %; Fig. 4b) indicating that respiration, or at least mitochondrial metabolism, favors the formation of telomere hyperclusters upon abrupt starvation. It is noteworthy that rho- cells show very short chronological lifespan in SP (Fig. 4a), consistent with our observation that telomere hyperclusters are a feature of long-lived quiescent cells. These data indicate that the ability to form hyperclusters is favored by mitochondrial activity.
Hormetic ROS during exponential phase prime cells to form hyperclusters upon starvation and to sustain long-term viability
We reasoned that ROS, as byproducts of the respiration process, could prime cells to form telomere hyperclusters upon starvation. Indeed, studies in model organisms show that a mild increase in ROS levels can positively influence health and lifespan, a process defined as mitochondrial hormesis or mitohormesis [20, 31]. Since hydrogen peroxide (H2O2) has emerged as a ROS signaling molecule able to induce an adaptive response , we tested the effect of increasing intracellular H2O2 on telomere hypercluster formation. This was achieved either by deleting the gene encoding the cytoplasmic catalase Ctt1, which scavenges H2O2 , or by overexpressing the superoxide dismutase Sod2, which converts O2- into H2O2 (Fig. 5a, b). In agreement with our hypothesis, we observed that telomere hyperclusters formed more efficiently in SP of ctt1∆ cells, and appeared earlier in cells overexpressing SOD2, compared with wild-type cells (Fig. 5a, b). Importantly, these strains deleted for CTT1 or overexpressing SOD2 both show extended lifespan [33, 34].
We next tested whether increasing ROS levels in fermenting cells by treating them with H2O2 would bypass the requirement for the respiration phase and promote hypercluster formation upon starvation. As expected, untreated cells were unable to form telomere hyperclusters after starvation (Fig. 5c) and had a short CLS (Fig. 5d). In contrast, H2O2 pre-treated cells contained brighter and fewer Rap1-GFP foci (Fig. 5c). Importantly, like SP HD cells, H2O2 pre-treated cells had >90 % viability at day 7 of CLS (Fig. 5d). Combined, these data strongly suggest that ROS exposure prior to starvation promotes telomere grouping and long-term viability during starvation.
Sir3-dependent telomere clustering favors long term survival during quiescence
We previously demonstrated that telomere grouping in exponentially growing cells is dependent on Sir3 protein amount but independent of silencing . We found that telomere hyperclustering in wild-type quiescent cells is not driven by an increase in Sir3 protein levels as revealed by western blot analysis (Additional file 6: Figure S4a). Furthermore, monitoring Sir3 occupancy genome-wide by chromatin immunoprecipitation (ChIP) revealed no significant changes in Sir3 spreading between exponentially growing cells and SP cells showing telomere hyperclusters (Additional file 6: Figure S4b).
To evaluate whether the silencing function of Sir3 is required for telomere hyperclustering and for longevity in quiescent cells, we transformed sir3Δ cells (defective for telomere clustering) with either a wild-type or a silencing dead copy of SIR3 (sir3-A2Q)  and assessed their CLS. We found that the insertion of either SIR3 or sir3-A2Q rescued the telomere hyperclustering in quiescent cells (Fig. 6a). We noticed that Rap1-GFP foci in the sir3-A2Q mutant were dimmer than in the SIR3 strain, probably due to a lower stability of this mutant form of Sir3 in SP (Fig. 6b). Nevertheless, this establishes that the silencing function of Sir3 is not required for telomere clustering in quiescence.
The sir3∆ strain had viability similar to wild-type cells at days 4 and 7 (Fig. 6c and not shown), arguing that this mutant enters properly into quiescence upon carbon source exhaustion. In agreement with this, we observed that the sir3∆ strain generates dense cells following the diauxic shift. Furthermore, these cells are as thermo-tolerant as their wild-type counterpart (Additional file 6: Figure S4c–e). In contrast, the sir3∆ strain shows a modest decrease in viability after day 10 compared with wild type, suggesting that while Sir3 is dispensable to enter into the quiescent state, it contributes to the maintenance of this specific cellular state. Importantly, expressing the sir3-A2Q mutant rescued the viability defect observed in the sir3∆ strain (Fig. 6c). Thus, Sir3-mediated telomere clustering but not silencing is required for the maintenance of the quiescent state.
We report that the organization of the budding yeast genome changes drastically depending on the metabolic status of the cell. In particular, quiescent cells that sustain long-term viability or increased CLS form a discrete subcompartment of telomeric silent chromatin in the most internal part of the nucleus (Fig. 6d).
Dynamics of nuclear organization upon carbon source exhaustion
We describe the dynamics of nuclear organization upon two major metabolic transitions: from fermentation to respiration and from respiration to SP. First, we show that telomere clusters, which are known to form three to five foci at the nuclear periphery in cells undergoing fermentation, form brighter and fewer foci after the diauxic shift. Furthermore, when cells exhaust the carbon source after respiration and enter the SP, these foci further group into a hypercluster located in the center of the nucleus in SP cells able to sustain long-term viability.
SIR-mediated telomere clustering drives chromosome conformation in long-lived quiescent cells
Genomic 3C analyses reveal that long-lived SP cells display increased constraints in their nuclear architecture, which appears to be driven by the clustering of telomeres. Because S. cerevisiae chromosomes exhibit such differences in size, mechanical constraints are likely to play significant roles on the organization of chromosomes tethered at both their centromere and telomeric regions. The positioning of the cluster in the middle of the nuclear space may actually reflect this physical constraint imposed by the smallest chromosome arms. As SIR3-deleted cells are unable to form telomere hyperclusters in quiescence and show a global organization that is similar to that of G1 cells, we conclude that SIR-mediated telomere clustering drives the global reorganization of chromosomes in long-lived quiescent cells. Although both Sir3 and Sir4 are required for telomere hyperclustering, gene silencing is not necessary for this event, as demonstrated by expressing a silencing defective version of Sir3 . Furthermore, telomere hyperclustering in quiescent cells is not driven by an increase in Sir3 protein or an increase of Sir3 spreading. As Sir3 may bind nucleosomes in more than one conformation , it is possible that telomere clustering after the diauxic shift is driven by specific post-translational modifications that increase Sir3 clustering function.
Mitochondrial ROS commit cells to form telomere hyperclusters upon starvation
Importantly, we show that increased telomere clustering is not a general feature of cell cycle arrest, as it is observed only in quiescent cells able to sustain long-term viability. Furthermore, the ability to form telomere hyperclusters required mitochondrial activity and is acquired post-diauxic shift in the quiescent fraction of cells shown to have a six-fold higher respiration rate compared with the non-quiescent fraction of cells . ROS, and more specifically H2O2 produced during respiration, are obvious candidates to trigger the commitment to form hyperclusters upon starvation . Indeed, we show that mutants known to increase the cellular level of H2O2 form hyperclusters with a higher rate and faster kinetics than wild-type cells. Furthermore, treating pre-diauxic shift cells with a sub-lethal dose of H2O2 is sufficient to commit cells to form telomere hyperclusters upon starvation and to sustain long-term viability. This commitment could be mediated by the checkpoint kinase Rad53, which is activated at these levels of H2O2 , thus allowing crosstalk between mitochondria and the nucleus [38, 39].
Potential benefits of telomere hyperclustering for CLS
Although alterations of nuclear architecture have been reported upon differentiation  and in quiescent metazoan cells , the function of this reorganization remains elusive. Interestingly, dramatic changes in the distribution of silent chromatin are observed in mammalian senescent cells with the formation of senescence-associated heterochromatin foci, which are thought to contribute to the stability of the cell cycle arrest . Another striking example of genome reorganization comes from rod photoreceptor cells of nocturnal rodents. In these cells, the nuclei exhibit an “inverted organization” — that is, reminiscent to the hypercluster observed in long live yeast cells — probably as an adaptation to limited light .
The large reorganization of budding yeast telomeres into a hypercluster concomitant with an important metabolic adaptation most likely provides a survival advantage in the long-term. Accordingly, sir3∆ strains, which cannot form telomere clusters, show a modest reduction in longevity compared with wild-type strains, when SP cultures (after 3 days in rich medium) were shifted to water. This is consistent with the findings of . However, quiescent cells purified from 7-day cultures of prototroph W303 strains showed no difference in the lifespan of sir3∆ or sir4∆ and wild-type cells (Linda Breeden, personal communication), possibly due to strain or experimental procedure variations. Importantly, the viability defect that we observed is rescued by expressing a SIR3 allele that is competent for telomere clustering but defective for silencing (sir3-A2Q mutant ), indicating that telomere clustering in quiescence has a positive effect on CLS independent of gene silencing under our conditions.
We propose that telomere hyperclusters could influence survival by protecting telomeres from degradation, fusion, and/or ectopic recombination events. Alternatively, telomere hyperclustering in quiescence could also be a way to sequester multifunctional factors that could have deleterious effects if localized to nuclear subcompartments where they are not needed. Such a factor could be the sirtuin Sir2, since it plays a pro-aging role by regulating cytoplasmic enzymes involved in carbon metabolism [44, 45].
By establishing that the nuclear organization of quiescent cells significantly differs from the well-described organization of cells grown in nutrient-replete conditions, our study sets the ground to (re)interpret studies on nuclear processes in the context of quiescence and aging. Moreover, our results unravel a novel connection between nuclear organization and aging, paving the way for future experiments analyzing the importance of nuclear organization for chronological lifespan.
Materials and methods
Media and growth conditions
All yeast strains used in this work are listed in Additional file 7 and are from the W303 background  except for the strains used for the HiC experiment (BY4741 ). Gene deletions and gene tagging were performed by PCR-based gene targeting [46, 47].
Yeast cells were grown in rich medium (YPD, yeast extract-peptone-dextrose) at 30 °C.
Induction of quiescence by carbon source exhaustion was performed as follows. Yeast cells were inoculated in YPD and grown overnight. The following day, cultures were diluted to an optical density of 0.2 (OD600nm) and grown at 30 °C in agitation for 5–6 h (fermentation), 24–48 h (respiration) or more than 7 days (SP). Levels of glucose in the medium were determined by using the D-Glucose HK assay kit (Megazyme). Induction of quiescence by carbon source starvation was performed by growing the cells in YPD at 30 °C (before or after glucose exhaustion) and then transferring them to exhausted YPD or sterile water for at least 16 h. For nitrogen starvation experiments, cells were grown to an OD600nm of 1 and transferred to a synthetic medium containing 0.17 % yeast nitrogen base (MP Biomedical) and 2 % glucose.
Density gradient fractionation
For density gradient fractionation, a solution of Percoll (Sigma-Aldrich) with a final NaCl concentration of 167 mM was added to a 30 ml Corex tube and centrifuged at 13,000 rpm for 20 min.
Approximately 2 × 109 cells were harvested, resuspended in 1 ml Tris buffer, added to the preformed gradient and centrifuged at 400 gav for 60 min at 20 °C. Density gradient tubes were imaged, and fractions collected, washed once in water, and used directly for assays or split into aliquots, pelleted, and frozen in liquid nitrogen. Cell number was determined for each fraction.
Viability (colony forming unit) assay
To test quiescent cells’ colony forming ability, cultures were grown as indicated. After 24 h of quiescence induction (day 1 CLS), 50 μl of each culture was collected, diluted 1:1.2 × 106 and plated in YPD plates. Culture tubes were agitated at 30 °C for 7 days and plated. Colonies were counted after 3 days at 30 °C. Day 7 CLS was normalized to day 1 CLS. Plots represent the mean value obtained for at least three independent experiments; error bars correspond to standard error of the mean.
To test whether direct addition of ROS in the medium of cultures undergoing fermentation could commit nuclei to form telomere hyperclusters during starvation, cells grown overnight were diluted to 0.002 OD600nm/ml in fresh YPD containing no drugs or H2O2 1 mM, grown until they reached 1 OD600nm/ml, and then starved in water for at least 24 h.
For protein isolation, 200 μl of trichloroacetic acid (TCA) 20 %, 200 μl of TCA buffer [20 mM Tris–HCl pH 8, 50 mM ammonium acetate, 2 mM EDTA, 1 mM phenylmethylsulfonyl fluoride (PMSF)], 1 μl of Protease inhibitor cocktail (Sigma-Aldrich), and 400 μl of acid-washed glass beads (710–1180 μm; Sigma-Aldrich) were added to 1 × 108 pelleted cells. Cells were then disrupted by vigorous vortexing (1 min, two times). Resulting extracts were centrifuged for 30 min at 4 °C at 14,000 rpm, and pellets were resuspended in 200 μl of TCA-Laemmli loading buffer (120 mM Tris base, 3.5 % sodium dodecyl sulfate (SDS), 8 mM EDTA, 5 % β-mercaptoethanol, 1 mM PMSF, 15 % glycerol, 0.01 % bromophenol blue). Samples were boiled for 10 min and centrifuged at 14,000 rpm for 10 min. Aliquots were immediately loaded or frozen. For immunoblotting, we used custom-made polyclonal antibodies against Rap1 (Agrobio, raised against Rap1[358–828] recombinant protein (a generous gift from M.H. LeDu, CEA Saclay) and Sir3 at 1:5000 . Loading was normalized according to H2A at 1:5000 (Abcam).
Immuno-FISH experiments were performed as in  with minor modifications. For quiescent cells, spheroplasting time was increased (20 min instead of 10 min).
Sets of images from any given figure panel were acquired the same day using identical acquisition parameters, except for time course experiments where the same culture was imaged at different time points, using identical acquisition parameters and using a wild-type growing culture as control. Details are provided in Additional file 8.
Quantification of Rap1 foci
A dedicated tool has been designed to find and quantify the telomere cluster in the 3D images acquired with fluorescence microscopy. Details are provided in Additional file 8.
Construction of 3C libraries and sequencing
S. cerevisiae G1 daughter cells (strain BY4741) were recovered from an exponentially growing population through an elutriation procedure . Long-lived quiescent cells were recovered as described above. 3C libraries were generated as described  with minor changes in the protocol. Briefly, the cells were cross-linked for 20 minutes with fresh formaldehyde (3 % final concentration), pooled as aliquots of 3 × 109 cells, and stored at −80 °C until use. Aliquots were thawed on ice and resuspended in 6 ml 1× DpnII buffer (NEB). The cells were then split into four tubes and lysed using a Precellys grinder (3 cycles of 6500 rpm, 30 s ON/60 s OFF) and VK05 beads. The cells were incubated for 3 h with 50 units of restriction enzyme under agitation (DpnII; NEB). The digestion mix was then diluted into ligation buffer and a ligation was performed at 16 °C for 4 h followed by a decrosslinking step consisting of an overnight incubation at 65 °C in the presence of 250 μg/ml proteinase K in 6.2 mM EDTA. DNA was then precipitated, resuspended in TE buffer, and treated with RNAse.
The resulting 3C libraries were sheared and processed into Illumina libraries using custom-made versions of the Illumina paired-end adapters (Paired-End DNA Sample Prep Kit, Illumina PE-930-1001). Fragments of sizes between 400 and 800 bp were purified using a PippinPrep apparatus (SAGE Science), PCR amplified, and paired-end sequenced on an Illumina platform (HiSeq2000; 2 × 100 bp).
Processing of paired-end reads
The raw data from each 3C experiment were processed as follows First, PCR duplicates were collapsed using the six Ns present on each of the custom-made adapters. Reads were then aligned using Bowtie 2 in its most sensitive mode against the S. cerevisiae reference genome . Paired-end reads were aligned as follows: for each read the length of the sequence mapped was increased gradually from 20 bp until the mapping became unambiguous (mapping quality >40). Paired reads were aligned independently.
Generation of contact maps
Each mapped read was assigned to a restriction fragment. Genome-wide contact matrices were built by binning the genome into units of 20 restriction fragments, resulting in 1797 × 1797 contact maps. The contact maps were subsequently filtered and normalized using the sequential component normalization procedure described in . This procedure ensures that the sum over the column and lines of the matrix equals 1 and reduces the influence of biases inherent to the protocol. Full resolution contact maps binned at ten restriction fragments are available in the supplemental material section (Additional files 9, 10, 11 and 12). The 3D structures were directly computed from the normalized contact maps using ShRec3D . The algorithm first computes a distance matrix from the contact map by assuming that the distance between each pair of beads is equal to the shortest path on the weighted graph associated with the inverse of the contact frequency between the two corresponding nodes. Multi-dimensional scaling is then applied to recover the optimal 3D coordinates from this distance matrix. To allow direct comparison between the structures obtained in different conditions we first re-scaled them to equalize the volume occupied by their associated convex hull. We then scaled the distances in each structure to account for the measured difference in size between nuclei in G0 and G1 daughter cells (1.5 and 1.7 μm, respectively; data not shown and ). Telomere pair distances were then directly computed from the structures to assess telomere re-organization.
The sequences of the chromosome conformation capture experiments reported in this paper have been deposited in BioProject with accession number PRJNA291473 . Microarray data are available from the Gene Expression Omnibus (GEO) under the accession number [GEO:GSE71273]. Microscopy data are available from Figshare .
capture of chromosome conformation
fluorescence in situ hybridization
green fluorescent protein
reactive oxygen species
Van Bortle K, Corces VG. Nuclear organization and genome function. Annu Rev Cell Dev Biol. 2012;28:163–87.
Meister P, Taddei A. Building silent compartments at the nuclear periphery: a recurrent theme. Curr Opin Genet Dev. 2013;23:96–103.
Tanaka A, Tanizawa H, Sriswasdi S, Iwasaki O, Chatterjee AG, Speicher DW, et al. Epigenetic regulation of condensin-mediated genome organization during the cell cycle and upon DNA damage through histone H3 lysine 56 acetylation. Mol Cell. 2012;48:532–46.
Carone DM, Lawrence JB. Heterochromatin instability in cancer: from the Barr body to satellites and the nuclear periphery. Semin Cancer Biol. 2012;23:99.
Andrulis ED, Neiman AM, Zappulla DC, Sternglanz R. Perinuclear localization of chromatin facilitates transcriptional silencing. Nature. 1998;394:592–5.
Gotta M, Laroche T, Formenton A, Maillet L, Scherthan H, Gasser SM. The clustering of telomeres and colocalization with Rap1, Sir3, and Sir4 proteins in wild-type Saccharomyces cerevisiae. J Cell Biol. 1996;134:1349–63.
Maillet L, Boscheron C, Gotta M, Marcand S, Gilson E, Gasser SM. Evidence for silencing compartments within the yeast nucleus: a role for telomere proximity and Sir protein concentration in silencer-mediated repression. Genes Dev. 1996;10:1796–811.
Taddei A, Van Houwe G, Nagai S, Erb I, van Nimwegen E, Gasser SM. The functional importance of telomere clustering: Global changes in gene expression result from SIR factor dispersion. Genome Res. 2009;19:611–25.
Ferreira HC, Luke B, Schober H, Kalck V, Lingner J, Gasser SM. The PIAS homologue Siz2 regulates perinuclear telomere position and telomerase activity in budding yeast. Nat Cell Biol. 2011;13:867–74.
Schober H, Ferreira H, Kalck V, Gehlen LR, Gasser SM. Yeast telomerase and the SUN domain protein Mps3 anchor telomeres and repress subtelomeric recombination. Genes Dev. 2009;23:928–38.
Duan Z, Andronescu M, Schutz K, McIlwain S, Kim YJ, Lee C, et al. A three-dimensional model of the yeast genome. Nature. 2010;465:363–7.
Taddei A, Gasser SM. Structure and function in the budding yeast nucleus. Genetics. 2012;192(1):107-29.
Rabl C. über Zellteilung. Morphol Jahrbuch. 1885;10:214–330.
Tjong H, Gong K, Chen L, Alber F. Physical tethering and volume exclusion determine higher-order genome organization in budding yeast. Genome Res. 2012;22:1295–305.
Tokuda N, Terada TP, Sasai M. Dynamical modeling of three-dimensional genome organization in interphase budding yeast. Biophys J. 2012;102:296–304.
Wong H, Marie-Nelly H, Herbert S, Carrivain P, Blanc H, Koszul R, et al. A predictive computational model of the dynamic 3D interphase yeast nucleus. Curr Biol. 2012;22:1881–90.
Broach JR. Nutritional control of growth and development in yeast. Genetics. 2012;192:73–105.
Martinez MJ, Roy S, Archuletta AB, Wentzell PD, Anna-Arriola SS, Rodriguez AL, et al. Genomic analysis of stationary-phase and exit in Saccharomyces cerevisiae: gene expression and identification of novel essential genes. Mol Biol Cell. 2004;15:5295–305.
Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, et al. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000;11:4241–57.
Ristow M, Schmeisser S. Extending life span by increasing oxidative stress. Free Radic Biol Med. 2011;51:327–36.
Klosinska MM, Crutchfield CA, Bradley PH, Rabinowitz JD, Broach JR. Yeast cells can access distinct quiescent states. Genes Dev. 2011;25:336–49.
Ruault M, De Meyer A, Loiodice I, Taddei A. Clustering heterochromatin: Sir3 promotes telomere clustering independently of silencing in yeast. J Cell Biol. 2011;192:417–31.
Fuge EK, Braun EL, Werner-Washburne M. Protein synthesis in long-term stationary-phase cultures of Saccharomyces cerevisiae. J Bacteriol. 1994;176:5802–13.
Laporte D, Courtout F, Salin B, Ceschin J, Sagot I. An array of nuclear microtubules reorganizes the budding yeast nucleus during quiescence. J Cell Biol. 2013;203:585–94.
Allen C, Buttner S, Aragon AD, Thomas JA, Meirelles O, Jaetao JE, et al. Isolation of quiescent and nonquiescent cells from yeast stationary-phase cultures. J Cell Biol. 2006;174:89–100.
Dekker J, Rippe K, Dekker M, Kleckner N. Capturing chromosome conformation. Science. 2002;295:1306–11.
Lesne A, Riposo J, Cournac A, Mozziconacci J. A fast algorithm to reconstruct 3D genome structures from chromosomal contact maps. Nat Methods. 2014;11:1141–3.
Therizols P, Duong T, Dujon B, Zimmer C, Fabre E. Chromosome arm length and nuclear constraints determine the dynamic relationship of yeast subtelomeres. Proc Natl Acad Sci U S A. 2010;107:2025–30.
Li L, Miles S, Melville Z, Prasad A, Bradley G, Breeden LL. Key events during the transition from rapid growth to quiescence in budding yeast require posttranscriptional regulators. Mol Biol Cell. 2013;24:3697–709.
Lu C, Brauer MJ, Botstein D. Slow growth induces heat-shock resistance in normal and respiratory-deficient yeast. Mol Biol Cell. 2009;20:891–903.
Yun J, Finkel T. Mitohormesis. Cell Metab. 2014;19:757–66.
Veal EA, Day AM, Morgan BA. Hydrogen peroxide sensing and signaling. Mol Cell. 2007;26:1–14.
Mesquita A, Weinberger M, Silva A, Sampaio-Marques B, Almeida B, Leao C, et al. Caloric restriction or catalase inactivation extends yeast chronological lifespan by inducing H2O2 and superoxide dismutase activity. Proc Natl Acad Sci U S A. 2010;107:15123–8.
Fabrizio P, Liou LL, Moy VN, Diaspro A, Valentine JS, Gralla EB, et al. SOD2 functions downstream of Sch9 to extend longevity in yeast. Genetics. 2003;163:35–46.
Oppikofer M, Kueng S, Gasser SM. SIR-nucleosome interactions: structure-function relationships in yeast silent chromatin. Gene. 2013;527:10–25.
Davidson GS, Joe RM, Roy S, Meirelles O, Allen CP, Wilson MR, et al. The proteomics of quiescent and nonquiescent cell differentiation in yeast stationary-phase cultures. Mol Biol Cell. 2011;22:988–98.
Leroy C, Mann C, Marsolier MC. Silent repair accounts for cell cycle specificity in the signaling of oxidative DNA lesions. EMBO J. 2001;20:2896–906.
Schroeder EA, Raimundo N, Shadel GS. Epigenetic silencing mediates mitochondria stress-induced longevity. Cell Metab. 2013;17:954–64.
Tsang CK, Liu Y, Thomas J, Zhang Y, Zheng XF. Superoxide dismutase 1 acts as a nuclear transcription factor to regulate oxidative stress resistance. Nat Commun. 2014;5:3446.
Dixon JR, Jung I, Selvaraj S, Shen Y, Antosiewicz-Bourget JE, Lee AY, et al. Chromatin architecture reorganization during stem cell differentiation. Nature. 2015;518:331–6.
Bridger JM, Boyle S, Kill IR, Bickmore WA. Re-modelling of nuclear architecture in quiescent and senescent human fibroblasts. Curr Biol. 2000;10:149–52.
Corpet A, Stucki M. Chromatin maintenance and dynamics in senescence: a spotlight on SAHF formation and the epigenome of senescent cells. Chromosoma. 2014;123:423.
Solovei I, Kreysing M, Lanctot C, Kosem S, Peichl L, Cremer T, et al. Nuclear architecture of rod photoreceptor cells adapts to vision in mammalian evolution. Cell. 2009;137:356–68.
Fabrizio P, Gattazzo C, Battistella L, Wei M, Cheng C, McGrew K, et al. Sir2 blocks extreme life-span extension. Cell. 2005;123:655–67.
Casatta N, Porro A, Orlandi I, Brambilla L, Vai M. Lack of Sir2 increases acetate consumption and decreases extracellular pro-aging factors. Biochim Biophys Acta. 2013;1833:593–601.
Thomas BJ, Rothstein R. Elevated recombination rates in transcriptionally active DNA. Cell. 1989;56:619–30.
Longtine MS, McKenzie 3rd A, Demarini DJ, Shah NG, Wach A, Brachat A, et al. Additional modules for versatile and economical PCR-based gene deletion and modification in Saccharomyces cerevisiae. Yeast. 1998;14:953–61.
Marbouty M, Ermont C, Dujon B, Richard GF, Koszul R. Purification of G1 daughter cells from different Saccharomycetes species through an optimized centrifugal elutriation procedure. Yeast. 2014;31:159–66.
Marie-Nelly H, Marbouty M, Cournac A, Liti G, Fischer G, Zimmer C, et al. Filling annotation gaps in yeast genomes using genome-wide contact maps. Bioinformatics. 2014;30:2105–13.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Cournac A, Marie-Nelly H, Marbouty M, Koszul R, Mozziconacci J. Normalization of a chromosomal contact map. BMC Genomics. 2012;13:436.
Heun P, Laroche T, Shimada K, Furrer P, Gasser SM. Chromosome dynamics in the yeast interphase nucleus. Science. 2001;294:2181–6.
BioProject: Saccharomyces cerevisiae strain:BY4741/W303 (baker’s yeast). http://www.ncbi.nlm.nih.gov/bioproject/291473.
figshare: Guidi et al. 2015. http://dx.doi.org/10.6084/m9.figshare.1505174.
We thank Patricia Le Baccon for assistance with microscopy, M.H. LeDu for providing recombinant Rap1, and Sheldon Decombe for his help on Fig. 6b. We also thank Geneviève Almouzni for critical reading of the manuscript and the members of the Taddei laboratory. AT and MG thank Isabelle Sagot for discussion and sharing unpublished results. RK and MM thank Alain Jacquier and Gwenael Badis-Breard for their constructive suggestions and comments and their help in the preparation of BY quiescent strains. M. Guidi was supported by a doctoral fellowship from the Association pour la Recherche sur le Cancer and the Institut Curie. M.M. is the recipient of an Association pour la Recherche sur le Cancer fellowship 20100600373. The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Program (FP7/2007 2013/European Research Council grant agreement 210508 to AT and 260822 to RK) and the investissement d’avenir (Labex DEEP) Agence Nationale de la Recherche to A.T.
All the authors declare that they have no competing interests.
MG generated strains and performed microscopy, lifespan, western blots and ROS treatment experiments. MR generated strains and contributed to microscopy experiments and writing/revising the manuscript. MM performed Hi-C experiments. IL contributed to lifespan experiments. AC conducted the bioinformatics analysis of the contact data. CB developed image analysis tools. AH performed and analyzed ChIP on ChIP experiments. JM built the 3D structures and performed bioinformatics analysis. AT, MG, MR and RK contributed to the design and interpretation of the study, drafting the figures and writing/revising the manuscript. All authors read and approved the final manuscript.
Myriam Ruault and Martial Marbouty contributed equally to this work.
Characterization of the SP silent chromatin hypercluster. a Western blot against Rap1 on crude extracts from exponential, respiratory, or stationary cultures of a WT strain (yAT1684). H2A antibody was used for the loading control. b Representative fluorescent images of wild-type (WT) strains tagged with Rap1-GFP “yAT 1684”, GFP-Sir2 “yAT405”, Sir3-GFP “yAT779” and GFP-Sir4 “yAT431” strains. Overnight liquid cultures were diluted to 0.2 OD600nm/ml and images were acquired after 5 h (1 OD600nm/ml, fermentation phase) and 7 days (40 OD600nm/ml, stationary phase). c Representative fluorescent image of a Rap1-GFP Sir3-mCherry-tagged strain “yAT194” from stationary phase cultures. We note that Sir3 associates with both telomeres and the rDNA in stationary phase cells. d Representative fluorescent images of Rap1-GFP in stationary cultures of WT “yAT1684” and sir4∆ “yAT2092” strains. e Representative fluorescent images of the nucleolar protein Sik1 tagged with mCherry during fermentation, respiration, and stationary phase (“yAT340”). f Representative fluorescent image of Rap1-GFP Dad2-mRFP (Duo1 And Dam1 interacting, an essential component of the microtubule–kinetochore interface) tagged stationary phase cells (“yAT2279”). g Representative fluorescent image of Sir3-mCherry Cse4-GFP-tagged strain “yAT2280” from stationary phase. Scale bar is 1 μm. (PDF 1343 kb)
Related to Fig. 3. Animated 3D reconstruction of the entire contact map of G1 cells. Each chromosome is represented as a chain of beads (1 bead = 20 kb), and the color code reflects the chromosome arm length, from blue for short arms to red for long arms. Each chromosome carries a black bead that corresponds to the centromere position. Yellow beads = subtelomeric regions; black beads = centromeres; purple beads = boundaries of the rDNA cluster on chromosome XII (in pink/red). (GIF 15597 kb)
SIR-mediated telomere clustering drives chromosome conformation in the dense fraction of SP cells. a Mean contacts frequencies between 100-kb centromeres windows in G1 (blue) and G0 quiescent cells (red). Black and green curves: contacts between 100-kb segments randomly sampled in both conditions, to illustrate the absence of coverage biases after normalization. b Chromosome organization of WT and sir3∆ quiescent cells (the cryptic mating type locus HML was deleted to prevent pseudo-diploid effect). ii) Normalized contact matrix obtained for hml∆* (left) and hml∆ sir3∆ (right) cells. Color scale: contact frequencies from rare (white) to frequent (dark blue). Red arrowheads: centromeres contacts; green and yellow arrowheads: telomere–telomere contacts in hml∆ and hml∆ sir3∆ G0 cells, respectively. The 3D representations of the hml∆ and hml∆ sir3∆ matrices are represented next to the contact maps. Each chromosome is represented as a chain of beads (1 bead = 20 kb), with color code reflecting the chromosome arm lengths, from short (blue) to long (red) arms. Yellow beads: subtelomeric regions; black beads: centromeres; purple beads: boundaries of the rDNA cluster. c Contact maps of W303 strain during exponentially growth (EXPO, left) and quiescence (G0, right). Red arrowheads: centromere clustering; green and yellow arrowheads: telomere–telomere contacts of two chromosomes (XIII and XV) in expo and G0 cells, respectively. Because of the low sequencing coverage and quality, the signal is not as strong as for data in Fig. 3 and the bins are larger (1 vector: 80 DpnII RFs). d Quantification of colocalization of 30-kb telomeric regions (red dots) compared with the distribution of the colocalization scores (box plot, two standard deviations) computed for 1000 random sets of 32 windows of 30 kb in the genome (excluding centromeric regions). The colocalization score is normalized by the sequencing depth for each dataset.
Telomere hyperclustering is not due to slow growth. a Representative fluorescent image of Rap1-GFP tagged strain grown either at 30 °C or 25 °C in exponential phase (top) and then starved for 16 h in water before imaging (bottom). b Calcofluor staining of LD and HD fractions of a post DS culture after gradient separation. c Heat shock (HS) assay on the LD and HD fractions used in b. (PDF 11591 kb)
Mechanism driving telomere clustering in long-lived SP cells. a Western blot against Sir3 and H2A on crude extracts from exponential, respiratory, or stationary cultures of a wild-type (WT) strain (yAT1684). b Sir3 spreading at yeast subtelomeres in cells from an exponentially growing culture (Fermentation) or in cells isolated from the dense fraction of a SP culture (Stationary HD). ChIP-chip profiles (Sir3 enrichment Z score) correspond to the mean of two independent experiments. Pearson correlation between conditions is 0.95. Sir3 spreading at TELVIR was confirmed in independent experiments by ChIP-quantitative PCR for both conditions (not shown). Each panel spans the first 30 kb from each telomere and the heading color for each panel indicates the middle repeat element content of the corresponding telomere: Y’ XCR XCS (beige), Y’ XCS (green), XCS (red), or XCR XCS (blue). Each dot represents a data point and lines are drawn for visual purposes. c Quiescent sir3∆ cells are as thermotolerant as quiescent WT cells to heat shock (HS). Dilution assays are shown (starting at DO600nm = 5 and diluted 1/5 each time). Left: growth control of exponential cells or 24 h LD cells. Middle: sensitivity to HS of WT exponential cells, 24 h LD cells or 24 h HD cells. Right: sensitivity to HS of WT or sir3∆ LD cells. d Stationary WT, sir3∆, sir3-A2Q cells are resistant to HS like WT cells. Dilution assays are shown (starting at DO600nm = 1 and diluted 1/5 each time). Left: growth control. Middle: 30 min 52 °C HS. Left: 1 h 52 °C HS. e Stationary WT and sir3∆ cells that spent 14 days in water after glucose exhaustion show the same extent of thermotolerance to a 1 h 52 °C HS. Dilution assays are shown (starting at DO600nm = 1 and diluted 1/5 each time). (PDF 12384 kb)
Strains used in this study. (DOCX 22 kb)
Additional experimental procedures. (DOCX 32 kb).
Full resolution contact map binned at 10 RF of G1 population presented in Fig. 3 . (DAT 24.7 mb)
Full resolution contact map binned at 10 RF of G0 population presented in Fig. 3 . (DAT 16.3 mb)
Full resolution contact map binned at 10 RF of G0 WT population presented in Additional file 3 . (DAT 24.6 mb)