Noncoding RNA-nucleated heterochromatin spreading is intrinsically labile and requires accessory elements for epigenetic stability

  1. R A Greenstein
  2. Stephen K Jones
  3. Eric C Spivey
  4. James R Rybarski
  5. Ilya J Finkelstein
  6. Bassem Al-Sady  Is a corresponding author
  1. George Williams Hooper Foundation, University of California San Francisco, United States
  2. University of California San Francisco, United States
  3. The University of Texas at Austin, United States

Abstract

The heterochromatin spreading reaction is a central contributor to the formation of gene-repressive structures, which are re-established with high positional precision, or fidelity, following replication. How the spreading reaction contributes to this fidelity is not clear. To resolve the origins of stable inheritance of repression, we probed the intrinsic character of spreading events in fission yeast using a system that quantitatively describes the spreading reaction in live single cells. We show that spreading triggered by noncoding RNA-nucleated elements is stochastic, multimodal, and fluctuates dynamically across time. This lack of stability correlates with high histone turnover. At the mating type locus, this unstable behavior is restrained by an accessory cis-acting element REIII, which represses histone turnover. Further, REIII safeguards epigenetic memory against environmental perturbations. Our results suggest that the most prevalent type of spreading, driven by noncoding RNA-nucleators, is epigenetically unstable and requires collaboration with accessory elements to achieve high fidelity.

https://doi.org/10.7554/eLife.32948.001

Introduction

The formation of gene-repressive heterochromatin domains is critical for genome integrity and for the establishment and maintenance of cell identity. Most heterochromatin formation occurs by a sequence-indifferent spreading reaction that propagates heterochromatic marks, structural proteins, and associated effector proteins outwards from nucleation sites. The precise extent of the spreading reaction has critical heritable consequences for cell identity. For example, in early pluripotent precursors, pre-existing heterochromatin domains spread, sometimes over megabases, to repress specifiers of inappropriate fates. Importantly, the final extent of spreading from a locus depends on the lineage pathway, hence it varies across different precursors (Wen et al., 2009; Zhu et al., 2013) and has to be precise to achieve a stable cell fate and avoid disease (Ceol et al., 2011). Similarly, spreading also specifies cell type in yeasts, where the cell type is maintained by repressing the mating cassettes at the mating type loci (Ekwall et al., 1991; Rusche et al., 2003). Despite the centrality of the spreading reaction in shaping cell identity, its native and intrinsic cellular characteristics, as well as mechanisms for its inter-generational propagation, have remained opaque.

We have some understanding of how cells inherit silencing at nucleation sites, which constitute the DNA-sequence driven component of heterochromatin. Recent results in heterochromatin systems signaled by Histone 3 Lysine 9 and Lysine 27 methylation (H3K9me and H3K27me) indicate that several mechanisms act together to ensure intergenerational inheritance: continuous DNA-mediated recruitment of the histone methylase (Audergon et al., 2015; Jia et al., 2004; Laprell et al., 2017; Ragunathan et al., 2015; Wang and Moazed, 2017), low histone turnover (Aygün et al., 2013; Taneja et al., 2017), as well as the positive ‘read-write’ feedback loop for histone methylases (Al-Sady et al., 2013; Zhang et al., 2008). Additionally, studies suggest that either the histone mark (Gaydos et al., 2014) or the histone methylases (Petruk et al., 2012) can persist trans-generationally.

These insights concerning nucleation sites do not necessarily account for how regions of heterochromatin distal to these sites are maintained. Unlike nucleation, which depends on DNA-based enzyme recruitment (Bayne et al., 2010; Verdel et al., 2004), spreading depends on the ability of the system to propagate along the chromosome, independent of the underlying DNA sequence. Such propagation requires the ‘read-write’ positive feedback function of the system (Al-Sady et al., 2013; Margueron et al., 2009; Müller et al., 2016; Noma et al., 2004; Zhang et al., 2008).

To determine how the spreading reaction acts in the maintenance of cell fate, it is central to understand the native behavior of two interconnected but separable phases of spreading: The initial spreading event, and the propagation of the states formed by this initial event through cell divisions. There is evidence that the initial spreading, at least in contexts outside the native chromosomal position, is stochastic, that is only some nucleation events result in a spreading event. This was first demonstrated by observing position effect variegation (PEV) in flies (Elgin and Reuter, 2013; Muller, 1930). Such stochastic behavior would have to be mitigated across cells to achieve a coherent specification outcome.

Intergenerational propagation of spreading is straightforward to conceptualize when epigenetic information is strongly reinforced, but more challenging in situations where modified nucleosomes are less likely to persist. This is the case for H3K9me-signaled heterochromatin in the fission yeast system, which lacks DNA methylation that can reinforce the epigenetic state. Persistence of the modified state is opposed by an anti-silencing protein Epe1 (Ayoub et al., 2003; Zofall and Grewal, 2006), which acts by antagonizing retention of H3K9me histones (Aygün et al., 2013; Ragunathan et al., 2015), and passage through S-phase, which significantly weakens heterochromatin domains (Chen et al., 2008). For fission yeast, there is evidence in favor of both high fidelity and stochastic propagation of the state formed by spreading. In support of a high fidelity model, theoretical work suggests that heterochromatin can display fundamentally bistable behavior, indicating that the ‘ON’ and ‘OFF’ states are intrinsically highly stable (Dodd et al., 2007). Similar bistable behavior has also been experimentally observed in plants (Angel et al., 2011, 2015). Conversely, the telomere position effect (TPE) observed in budding and fission yeast supports a model where intergenerational inheritance is fundamentally stochastic. In TPE the heterochromatic state is switched at high frequencies in the inheriting generations (Gottschling et al., 1990; Nimmo et al., 1994).

To distinguish whether spreading shapes and enables epigenetic maintenance of a cell identity locus via either of those modes, or combinations thereof, we focused on one of the most well understood heterochromatin loci, the fission yeast MAT locus, as a model. This locus remains tightly repressed to avoid simultaneous expression of both mating cassettes (Ekwall et al., 1991; Noma et al., 2001). The MAT locus contains two cis elements that directly recruit H3K9me. (1) cenH, which is related to the dg and dh repeats at the pericentromere and tlh2 at the subtelomere (Grewal and Klar, 1997; Hansen et al., 2006). These sequences nucleate H3K9me by at least two pathways, which depend on transcription of noncoding RNAs (ncRNAs): the RNAi pathway (Hall et al., 2002; Volpe et al., 2002), and at least one separate pathway dependent on nascent RNA polymerase II transcripts, which requires the budding yeast Nrd1 homology Seb1 (Marina et al., 2013) (collectively ‘ncRNA-nucleation’). Separately and unique to the MAT locus, (2) a region downstream of cenH including the REIII element, which recruits the H3K9 histone methylase, HP1 proteins and histone deacetylases (HDACs). This is dependent on REIII-bound transcription factors (Jia et al., 2004; Kim et al., 2004; Yamada et al., 2005), but is independent of RNA processes. Heterochromatin formation within the MAT locus is confined by boundary elements (Noma et al., 2001, 2006).

In this work, we probe heterochromatin spreading nucleated both at the MAT locus as well as ectopically in the genome with a ‘heterochromatin spreading sensor’ (HSS), which enables quantitative examination of spreading separately from nucleation in single S. pombe cells. Using the HSS, we show that ncRNA-dependent elements trigger epigenetically unstable spreading that is stabilized by an accessory RNA-independent cis-ating element. Both elements collaborate to form a high fidelity domain. The strategy we uncover has important implications for how heterochromatin spreading achieves and maintains ‘epigenetic’ character and can safeguard cell identity against environmental perturbations.

Results

A single-cell heterochromatin spreading sensor (HSS) controls for nucleation and cellular noise

To assess the intrinsic behavior of heterochromatin spreading and what shapes its precise re-establishment with respect to position and extent of repression (‘fidelity’), we employed transcriptionally encoded fluorescent reporters to read silencing by heterochromatin at a given locus, as previously reported. Several critical improvements over prior systems enable documentation of the spreading reaction at high sensitivity (Bintu et al., 2016; Hathaway et al., 2012; Obersriebnig et al., 2016; Osborne et al., 2009; Xu et al., 2006). First, our system has high signal to noise and minimized delay from epigenetic changes to fluorescent output. We accomplish this using the weak, well-characterized ade6 gene promoter (ade6p) (Allshire et al., 1994; Kagansky et al., 2009) to drive production of bright, fast-folding fluorescent proteins (XFPs) (Al-Sady et al., 2016). Second, our system provides separate sensors for nucleation, spreading, and cellular noise. We used ade6p-driven recoded super-folder GFP (Pédelacq et al., 2006) (‘green’) and monomeric Kusabira Orange (Sakaue-Sawano et al., 2008) (‘orange’) to report on nucleation and spreading, respectively (Figure 1A). A third XFP, an ade6p-driven triple fusion of E2Crimson (Strack et al., 2009) (‘red’, noise filter), is fully uncoupled from heterochromatin and inserted in a euchromatic locus. Here it reports on intrinsic or extrinsic noise that arises from cell-to-cell variation in the content of specific and general transcription factors and also translational efficiency (Figure 1A). To validate this reporter system, we characterized the non-heterochromatic state, via null mutation of clr4 (Δclr4), encoding the only S. pombe H3K9 methyltransferase. We show that in the absence of heterochromatin, expression of the noise reporter (‘red’) correlates well with that of reporters for both nucleation (‘green’) and spreading (‘orange’) (Figure 1—figure supplement 1A,B), especially when all cells in the population are considered without applying a size gate (Figure 1—figure supplement 1B, ρ ~0.83–0.93). This analysis mode is required when cell number is limiting. When a smaller subset is considered where all the cells are of similar size and stage of the cell cycle, the correlation still provides useful noise filtering (Figure 1—figure supplement 1A), which becomes evident when the normalization is applied to clr4+ cells that fall in the size gate (Figure 1—figure supplement 1C). Thus, cellular noise is mitigated by dividing the signals from the proximal ‘green’ and distal ‘orange’ heterochromatic reporters by the signal of the ‘red’, euchromatic reporter (‘green’/‘red’; ‘orange’/‘red’). Together, these elements constitute our heterochromatin spreading sensor (HSS) (Figure 1A).

Figure 1 with 1 supplement see all
Heterochromatin spreading from ncRNA-nucleated elements is stochastic and produces intermediate states.

(A) Overview of heterochromatin spreading sensor. Three transcriptionally encoded fluorescent proteins are inserted in the genome: The ‘clamp’ site enables isolation of successful nucleation events, the ‘sensor’ reports on spreading events and the ‘noise filter’ normalizes for cell-to-cell noise. (B) Overview of the ura4::dhHSS1-7kb strains. Genes downstream of the ‘green’ nucleation color are annotated. The alg11 gene is essential. (C) Spreading from ura4::dh visualized by the HSS with ‘orange’ inserted at different distances shown in (B). The ‘red’-normalized ‘orange’ fluorescence distribution of ‘green”OFF cells plotted on a histogram. Inset: 2D-density hexbin plot showing red-normalized ‘green’ and ‘orange’ fluorescence within the size gate, with no ‘green’ or ‘orange’ filtering. The ‘green'OFF population is schematically circled. The fluorescence values are normalized to = 1 for the Δclr4 derivate of each strain. (D) TOP: cartoon overview of the FACS experiment for D. and E. ‘green'OFF cells collected from the ura4::dhHSS3kb were separated in three populations (‘Low’, ‘Intermediate’ and ‘High’) as shown schematically based on the ‘orange’ fluorescence. BOTTOM: ‘orange’ RT-qPCR signal for the indicated populations. The y-axis is scaled to = 1 based on the ‘orange’ signal in Δclr4. Error bars indicate standard deviation of two replicate RNA isolations. (E) ChIP for H3K9me2 and H3K4me3 in the same populations as (D). Each ChIP is normalized over input and scaled to = 1 for a positive control locus (dh repeat for H3K9me2 and act1 promoter for H3K4me3). Error bars indicate standard deviation of two technical ChIP replicates. Primer pairs for RT-qPCR and ChIP are indicated by solid and dashed line, respectively, in the C. ura4::dhHSS3kb diagram.

https://doi.org/10.7554/eLife.32948.002

Spreading from ectopic ncRNA nucleators is stochastic and produces intermediate states

We first examined the intrinsic behavior of the heterochromatin spreading reaction in an ectopic context. We constructed the initial ectopic HSS based on a strain where a part of the centromeric ncRNA-driven nucleation element (dh) is inserted proximal to the endogenous ura4 gene (Canzio et al., 2011; Marina et al., 2013). We replaced the ura4+ open-reading frame (ORF) with ‘green’ to track nucleation element-proximal events. Then, to track distal events, we inserted ‘orange’ at one of several sites downstream from ‘green’ (ura4::dhHSS1kb, ura4::dhHSS3kb, ura4::dhHSS5kb ura4::dhHSS7kb, Figure 1B). The noise filter (‘red’) was inserted between SPBC1711.11 and SPBC1711.12, a bona fide euchromatic region (Garcia et al., 2015). All strains were initially constructed in a Δclr4 background, and we initiated heterochromatin formation by crossing in clr4+. We assessed heterochromatin formation after ~80–100 generations by quantifying the production of ‘green’ and ‘orange’. This period is significantly longer than ~25 generation timeframe required for full formation of a heterochromatic domain (Obersriebnig et al., 2016), ensuring that the population is at equilibrium.

To quantitatively assess the products of heterochromatin formation, we performed steady-state flow cytometry on log-phase cells, which were size-gated for small, recently divided cells (~91% G2, Figure 1—figure supplement 1D and supplemental experimental materials) to remove size- and cell cycle-related effects. At this stage, we only normalize the cells by the ‘red’ noise filter and scale the signal in each channel to Δclr4, giving us a broad overview of the possible expression states of ‘green’ and ‘orange’. We observe no cells that express ‘green’ but repress ‘orange’ (insets, Figure 1C), instead, all cells that are fully or partially ‘orange’ repressed are also robustly ‘green’ repressed. This observation, together with our finding that ‘green’ repression kinetically anticipates ‘orange’ repression (Figure 3—figure supplement 1), is consistent with heterochromatin spreading outward from the ura4::dh nucleator. Considering ‘green’ repression thus a proxy for nucleation, we observed that cells populate a wide range of nucleation states rather than a single state, with the distribution of repressed states varying among the HSS distance sensor strains (ura4::dhHSS1-7kb, Figure 1C). To specifically examine cells that have fully nucleated, we applied a computational ‘nucleation clamp’ that isolates cells with a ‘green’ signal that is lower than the median plus two standard deviations of wild-type cells containing no XFPs (see Appendix 1-Supplemental Materials and methods). Using ‘orange’ as a spreading proxy, we found spreading to be stochastic in nucleated cells, with some cells exhibiting full repression, but others partial repression or full de-repression (Δclr4, x = 1) of the ‘orange’ spreading sensor. The proportion of cells that are fully repressed by spreading declines linearly with distance (scheme, Figure 1B; data, Figure 1C). Intriguingly, cells that are not fully repressed mostly exhibit intermediate levels of repression, which are neither at values of full repression or de-repression.

We next assessed the nature of these intermediate states in the 3 kb distance reporter strain, where ~30% of cells had maximal repression at the ‘orange’ locus and the remainder had intermediate states ranging from strongly to weakly repressed. Using Fluorescence Activated Cell Sorting (FACS), we gated for successful nucleation in the ‘green’ channel and then binned the ‘orange’ channel for fully repressed (low), intermediate and de-repressed (high) populations (Figure 1D, cartoon). We queried each bin for molecular events associated with heterochromatin formation, using RT-qPCR to determine the expression levels of ‘orange’, and Chromatin Immunoprecipitation (ChIP) to query the presence of the marks H3K9me2 and H3K4me3. These marks are thought to be mutually exclusive, associating with repressed heterochromatin and active promoters, respectively (Noma et al., 2001). The message level of ‘orange’ is tightly repressed in the ‘low’ population (0.05 of max), partially repressed in the intermediate population (0.3 of max), and nearly fully ‘de-repressed’ (0.8 of max) in the ‘high’ population. Thus, cells with intermediate fluorescence also exhibit partial gene repression, demonstrating that fluorescence accurately reports on gene expression (Figure 1D, RT primers indicated in diagram in 1C, solid line). Histone modification levels also correlated well with the HSS signals (Figure 1E, ChIP primers indicated in diagram in 1C, dashed line). The ‘low’ fluorescence population has high H3K9me2 (0.9 of dh, positive control) and low H3K4me3 (0.09 of actin, positive control); the intermediate population had intermediate H3K9me2 (0.49 of dh) and H3K4me3 (0.23 of actin), and the high population had low H3K9me2 (0.2 of dh) and higher H3K4me3 (0.44 of actin). Hence, successfully nucleated cells with intermediate fluorescence also exhibit intermediate amounts of the mRNA for ‘orange’ and histone marks reflecting heterochromatin (H3K9me2) and transcriptional activity (H3K4me3). These results support the notion that intermediate states of repression observed by cytometry represent intermediate states of spreading.

These observations are not due to the particularities of the ectopic site chosen or the behavior of the XFPs, as our results are recapitulated at the his1 locus (his1::dhHSS3kb, Figure 1—figure supplement 1E), which contains only one gene (rec10) in the ‘spreading zone’, rather than several transcriptional units. Additionally, switching the nucleation and spreading reporter fluorophores produced similar results (Figure 1—figure supplement 1E). These results suggest that ncRNA-driven heterochromatin spreading at ectopic sites is intrinsically stochastic and multimodal, producing intermediate states of repression.

Distinct forms of heterochromatin spreading at MAT

We next examined spreading behavior at the endogenous mating type locus (MAT), which is tightly repressed (Grewal and Klar, 1997; Thon et al., 2002) and a bona fide high-fidelity locus, as it can behave in a bistable manner with stable epigenetic inheritance even when disrupted (Grewal and Klar, 1996). The MAT locus has two known elements shown to recruit the H3K9 methylase Clr4: the cenH element, homologous to the ncRNA-nucleated dh fragment we inserted at ura4 and his1, and the RNA-independent element termed REIII (Jia et al., 2004; Thon et al., 1999). At REIII, two stress-responsive transcription factors, Atf1 and Pcr1, which form a heterodimer (Wahls and Smith, 1994), recognize two DNA-binding sites within REIII, directly recruit Clr4, Swi6/HP1 and histone deacetylases (HDACs) (Jia et al., 2004; Kim et al., 2004) and are required for heterochromatin formation at MAT when cenH is compromised (Noma et al., 2004). We validated that MAT retains its well-documented tight repression following insertion of the HSS, placing the ‘green’ reporter within the cenH nucleator, and the ‘orange’ reporter proximal to the REIII nucleator. Both colors were fully repressed in the large majority of cells (Figure 2B), which is reproduced when the color orientations are reversed (Figure 2—figure supplement 1A). However, for both reporter configurations, the REIII proximal color showed a small proportion of cells that are slightly de-repressed compared to the cenH internal color, consistent with previous findings (Thon and Friis, 1997). We conclude that the HSS can be used to dissect spreading at the MAT locus.

Figure 2 with 2 supplements see all
ncRNA-dependent and independent nucleation yields qualitatively different spreading reactions in the MAT locus.

(A) Diagram of the reporters within MATHSS and ΔREIIIHSS. WT and m for REIII indicate the presence or deletion of the Atf1/Pcr1 binding sites, respectively. (B) 2D-density hexbin plot showing the ‘red’-normalized ‘green’ and ‘orange’ fluorescence for wild-type MATHSS cells. Scale bar shows every other bin cutoff as a fraction of the bin with the most cells. Inset: histogram of the ‘red’-normalized ‘orange’ fluorescence distribution of ‘green'OFF cells. (C) 2D-density hexbin plot and inset as above for ΔREIIIHSS, which contains two 7 bp Atf1/Pcr1-binding site deletions (m) within the REIII element. (D) ChIP for H3K9me2 (red) and H3K9me3 (grey) for amplicons indicated in (A). normalized to dh. WT, wild-type MATHSS, m, ΔREIIIHSS. (E) TOP: diagram of the reporters within ΔKHSS. The cenH nucleator and additional 5’ sequence is deleted and replaced by ‘orange’. ‘green’ is located directly proximal to REIII and serves as the nucleation clamp. ChIP amplicons are indicated as black bars. BOTTOM: 2D- density hexbin plot and inset as above. LEFT: ChIP for H3K9me2 (red) and H3K9me3 (grey) for ‘green’ and ‘orange’ in isolated ΔKHSS-ON or ΔKHSS-OFF alleles. In hexbin plots, the Δclr4 derivative of each strain was used to normalize the X- and Y-axes to = 1. Error bars indicate standard deviation of technical replicates.

https://doi.org/10.7554/eLife.32948.004

We then examined spreading in cells nucleated solely by the cenH element. The REIII nucleator was inactivated by deleting the critical cis-acting Atf1/Pcr1-binding sites, to create a strain designated ΔREIIIHSS (Figure 2C). To our surprise, the high fidelity that the MAT locus exhibits in the repressed state (Grewal and Klar, 1996) disappeared. Instead, cenH nucleated spreading in the ΔREIII strain behaved similarly to spreading from the ectopic ncRNA-nucleated strains, showing high stochasticity and predominantly intermediate repression states (Figure 2C). We wanted to address if this stochastic silencing is reflected in weakened heterochromatin assembly. We preformed ChIP for H3K9me2 and H3K9me3, marks signaling heterochromatin assembly (Nakayama et al., 2001) and repression or spreading (Al-Sady et al., 2013; Jih et al., 2017; Zhang et al., 2008), respectively. We found that these marks decline progressively towards the distal ‘orange’ reporter in ΔREIIIHSS (Figure 2D), compared to the wild-type (WT) MATHSS. This is consistent with the observed tight repression for WT MATHSS (Figure 2B) and weakened silencing at the distal ‘orange’ in ΔREIIIHSS (Figure 2D). It is possible that this difference in spreading results from an altered heterochromatin structure at cenH in ΔREIIIHSS. However, H3K9me2 and me3 accumulation does not differ between ΔREIIIHSS and WT MATHSS at the cenH nucleator, or the leftward REII locus (Figure 2D). Thus, the observed behavior of ΔREIIIHSS is consistent with stochastic and multimodal spreading, rather than compromised nucleation at cenH.

To examine heterochromatin formation independent of cenH, we used the historical ΔK strain, where the entire cenH nucleation element is deleted and replaced with a ura4+ reporter (Grewal and Klar, 1996). We introduced the HSS into this context (ΔKHSS, Figure 2E), placing the ‘green’ reporter proximal to REIII and the ‘orange’ reporter distally, replacing ura4. We then introduced clr4+ by cross and directly cultured colonies derived from germinated clr4+ spores. We found that although ΔKHSS has very weak nucleation compared to strains with intact ncRNA nucleators, the distribution of cells is sharply bimodal: Cells were either repressed at both reporters (‘OFF’, lower left corner) or de-repressed at both reporters (‘ON’, upper right corner; Figure 2E). We note that isolation of single colonies on nonselective media from original spores of the cross yields mostly ON (ΔKHSS-ON) or OFF (ΔKHSS-OFF) colonies, consistent with each state being metastable (Grewal and Klar, 1996Thon and Friis, 1997). This heterochromatin formation pattern requires REIII, as in 34/34 strains tested, no silencing can be established if Atf1/Pcr1 binding sites are deleted before clr4+ is introduced (Figure 2—figure supplement 2A,B). Additionally, the bimodal behavior does not require the H3K9me-independent gene-repressive REII element (Hansen et al., 2011), as ΔKHSS REII::LEU2, containing a deletion of REII, behaved similarly to ΔKHSS (Figure 2—figure supplement 1C), and is further independent of reporter placement (Figure 2—figure supplement 2C,D). We next characterized the molecular signature of the locus. While in our two color plots cells that were repressed in ‘green’ did not show any de-repression in ‘orange’ (Figure 2E, cells in bottom left corner), we wanted to test if the heterochromatic state at these loci correlated with this silencing pattern. Since we can isolate ΔKHSS-ON and ΔKHSS-OFF alleles by simple plating of ΔKHSS cells, we performed H3K9me2 ChIP on both and H3K9me3 ChIP for ΔKHSS-OFF cells (not detectable for ΔKHSS-ON). We found that methylation correlates with the repression state (Figure 2E) and importantly, does not significantly differ between ‘green’ and ‘orange’. Together, these result indicate that in ΔKHSS-OFF cells heterochromatin spreading is continuous across the locus and does not, unlike cenH-triggered spreading, accumulate any intermediates.

Multi-generational single-cell imaging reveals ncRNA-driven spreading to be unstable

Our measurements thus far cannot reveal the dynamics of transitions between states. This requires long-term imaging of cells over a substantial number of generations (>20), which is difficult with traditional microscopy because of cell crowding effects. To deal with this issue, we used the Fission Yeast Lifespan Micro-dissector (FYLM) microfluidic device (Spivey et al., 2017, 2014), which traps the old pole of a rod shaped S. pombe cell at the bottom of a chamber well for its entire lifetime. Sibling cells generated at the new pole by medial fission eventually exit the chamber. We continuously image the old-pole cell with fluorescence microscopy for up to 60 hr (Figure 3A). We note that unlike Saccharomyces cerevisiae, S. pombe does not execute an aging program but rather dies stochastically (Coelho et al., 2013; Nakaoka and Wakamoto, 2017; Spivey et al., 2017). Thus, imaging S. pombe over long timescales avoids the confounding effects of aging on epigenetic behavior (Guarente, 2000; Li et al., 2017). To capture the long-range dynamics of spreading, we imaged approximately one hundred cells of each strain concurrently (see Figure 3—figure supplement 2B for a summary of cell fates in all experiments). For each cell, we imaged all three channels continuously, and performed similar normalizations as for the flow cytometry data (Appendix 1-Supplemental Materials and methods). We first imaged the HSS distance sensor strain (ectopic ura4::dhHSS3kb). Our ability to observe cells that were initially fully de-repressed allowed us to trace ‘green’ and ‘orange’ repression kinetically. Consistent with linear heterochromatin spread outward of the dh nucleator, we find that ‘orange’ repression is anticipated by repression at ‘green’ (Figure 3—figure supplement 1). While nucleation in this strain is not stable (likely due to ‘green’ being adjacent to, rather than within dh), over time intervals where nucleation does persist, we observed dynamic fluctuations in the distal ‘orange’ color without a fixed temporal pattern (Figure 3—figure supplement 2A and Figure 3—videos 1 and 2), which is not due to the repression state of ‘green’ (Figure 3—figure supplement 2F).

Figure 3 with 7 supplements see all
Single-cell analysis of nucleation and spreading using a Fission Yeast Lifespan Micro-dissector (FYLM).

(A) Overview of the FYLM-based heterochromatin spreading assay. The old-pole cell is trapped at the bottom of one of hundreds of wells in the FYLM microfluidic device and is continuously imaged in brightfield (to enable cell annotation), green, orange and red channels. Hypothetical example traces are shown. (B) Maximum values attained by each nucleated cell for normalized ‘orange’ plotted against normalized ‘green’. Solid horizontal lines correspond to y = 0 and y = 0.5. Dashed line corresponds to an ON cutoff determined by mean less three standard deviations for each strain’s matched Δclr4 strain. Percentage of cells between each line was calculated. (C) FYLM analysis of wild-type MATHSS cells. CELL TRACES: 60 hr of normalized ‘green’ (left) and ‘orange’ (right) fluorescence in cells that maintained nucleation with the same five cells overlaid in different gray line styles in both plots. Gaps indicate loss of focus. HEATMAP: Up to 36 hr of normalized ‘orange’ fluorescence for 30 cells that maintained nucleation is represented from blue (0) to yellow (1). X-Y FLUORESCENCE PLOT: for one representative sample cell, plot of normalized ‘green’ and ‘orange’ fluorescence across its measured lifetime (grayscale). (D) FYLM analysis of ΔREIIIHSS cells as in C. The example cell in the X-Y dot plot is marked with an asterisk(*) on the orange traces (E) FYLM analysis of ΔKHSS-OFF isolate, as in C., D. All cells were normalized to Δclr4 (max, 1).

https://doi.org/10.7554/eLife.32948.007

Next, we analyzed the MAT locus strains and selected cells that maintained nucleation for their entire measured lifespan (Appendix 1-Supplemental Materials and methods). Under this constraint, the three strains exhibit vastly different behaviors (Figure 3B). WT MATHSS cells maintained ‘orange’ repression for the majority of their measured lifespans (Figure 3C, Figure 3—figure supplement 2C and Figure 3—video 3). However, we documented transient loss of ‘orange’ silencing for 20% of the cells. (Figure 3B and C). In contrast, while most cells stay similarly nucleated in ΔREIIIHSS (Figure 3D, Figure 3—figure supplement 2D) 83% of the cells imaged experienced at least half-maximal ‘orange’ de-repression at some time points (Figure 3B). For this strain, 30% of the cells transited through the fully ON state (Figure 3B and D, Figure 3—figure supplement 2D and Figure 3—video 4). In fact, cells sampled a wide range of values from OFF to fully ON, indicating that cells do not occupy ON or OFF states exclusively, but adopt intermediate values across time (Figure 3D). Importantly, ΔREIIIHSS cells, just as ura4::dhHSS3kb cells, fluctuate in their ‘orange’ values, indicating that spreading is unstable and adopts a random walk type behavior. To analyze ΔKHSS cells, which exist predominantly in fully ‘green’ and ‘orange’ ON state (Figure 2C), we analyzed ΔKHSS-OFF cells (see above). ΔKHSS-OFF behaved markedly differently from ΔREIIIHSS: in all the cells analyzed, ‘green’ and ‘orange’ reporters remained OFF throughout the time course (Figure 3B,E and Figure 3—video 5), up to 25 generations, revealing a fundamentally different dynamic behavior between cenH- and REIII-dependent heterochromatin. We note it remains possible that isolation of a ΔKHSS-OFF colony may bias our analysis against potentially more frequent OFF-ON switching events in the primary mixed population derived from continuous propagation of the germinated spore (Figure 2E). However, since the mixed population resolves spontaneously into ON and OFF states once plated, and OFF cells behave similarly in either the mixed ΔKHSS or ΔKHSS-OFF isolated populations (compare Figure 2E and Figure 5C), we believe the stability of ΔKHSS-OFF is intrinsic to the ΔK MAT locus.

Epigenetic stability at MAT is dependent on REIII

To probe memory capacity (i.e. the ability of cells to retain information of an ancestral state established many generations prior), we compared cells containing an intact MAT locus to those lacking either ncRNA- or REIII-dependent heterochromatin. We established two ancestral states (Figure 4A); one with unperturbed heterochromatin, and a second treated with the HDAC inhibitor trichostatin A (TSA), known to fully disrupt the heterochromatin state ([Hall et al., 2002] and Figure 4—figure supplement 1). Following production of the ancestral states, we grew cells either in rich media alone or in a TSA concentration gradient (0–50 µM) for 25 generations and then measured the fraction of fully nucleated cells that effectively silence the ‘orange’ spreading marker (Figure 4A). Cells exhibit memory if the fraction of the population with full spreading (‘orange'OFF) depends on the ancestral state, which would be indicated by separation of the unperturbed (light orange) and perturbed (red) lines. In contrast, no memory is indicated by convergence of the two lines (graphs in Figure 4B–D). We further measure a second parameter we term relative ‘resistance’, which is defined as the TSA concentration at which the fraction of cells with ‘orange'OFF declines to 50% of the no TSA pretreatment value. This value reports on the intrinsic sensitivity to perturbation of the locus formed by spreading.

Figure 4 with 2 supplements see all
ncRNA-nucleated spreading exhibits weak memory and resistance in the absence of REIII.

(A) Experimental schematic for memory and resistance measurements. Cells in log phase were treated with TSA (50 μM) for 10 generations to erase all heterochromatin (de-repressed, yellow) or kept untreated (repressed, gray). Both populations are then grown in a gradient of TSA concentration from 0 to 50 μM for 25 generations. (B) The wild-type MAT locus exhibits memory in silencing ‘orange’ throughout the TSA gradient. The fraction of ‘green'OFF cells that fully silence ‘orange’ normalized to the no TSA pre-treatment, 0 μM TSA point are plotted for each TSA concentration. Red line: cell ancestrally TSA pre-treated; light orange line: cells without pre-treatment. (C) Spreading from cenH exhibits weak memory and low resistance. Cell populations as above. (D) ncRNA-independent spreading exhibits high resistance. The fraction of ‘orange'OFF for all cells is plotted, because in the TSA pre-treatment almost no ‘green'OFF cells can be detected. Dotted lines indicate the half-resistance points: TSA concentration at which 50% of non-pretreated cells fail to form heterochromatin at ‘orange’. Memory is the difference between orange and red lines. One of two full biological repeats of the experiment is shown.( E) Experimental schematic for heat stress and recovery. Cells were grown at either 32 or 38°C for 10 generations and strains subsequently grown continuously for 96 hr at 32°C. (F) The fraction of cells with full spreading (‘green'OFF and ‘orange'OFF) after 38°C exposure and recovery normalized to the fraction of cells with full spreading at 32°C for each strain is plotted over time. For wild-type MATHSS and ΔREIIIHSS strains, we fit a simple sigmoidal dose response curve and determined a t1/2 value. The difference in t1/2 values or Δt1/2 is ~22 hr or ~9–10 generations.

https://doi.org/10.7554/eLife.32948.015

As expected, wild-type MATHSS exhibited clear memory at 25 generations (Figure 4B), which was still weakly evident even at 35 generations (Figure 4—figure supplement 1C). Among fully nucleated (‘green'OFF) cells, those that derived from untreated ancestral cells showed a greater fraction of silencing (‘orange'OFF) than those derived from treated cells throughout the entire TSA gradient, with a half-resistance point of ~2 µM (Figure 4B). Thus, wild-type MATHSS memory is robust in the face of perturbations of the heterochromatic state.

In sharp contrast, when spreading exclusively nucleates from cenH (ΔREIIIHSS strain), memory of silencing (‘orange'OFF) is significantly weaker. Memory collapsed beyond low TSA concentrations (>0.2 µM TSA), with the red and light orange lines coinciding for much of the gradient. Even at 0 μM TSA, history dependence was erased at 35 generations (Figure 4—figure supplement 2C). Interestingly, the half-resistance point was ~0.2 μM, 10-fold lower than that of wild-type MAT (Figure 4C). As cenH-nucleated spreading in ΔREIIIHSS produces little memory capacity and lacks resistance, the memory capacity at MAT does not derive from ncRNA-nucleated spreading. These results are consistent with REIII being required for the memory behavior of WT-MAT.

The ΔKHSS strain at face value had the widest separation in the behavior of the progeny of TSA pretreated and untreated cells. However, ascribing this behavior directly to memory is complicated by the fact that ΔKHSS cells are no longer able to re-nucleate if they were ancestrally TSA treated, consistent with previous findings indicating that RNAi factors are required for heterochromatin establishment at MAT (Hall et al., 2002). However, when examining resistance, that is the behavior of cells not ancestrally TSA pretreated, we observe that the REIII dependent ΔKHSS strain has a half-resistance point of ~3 µM TSA (Figure 4D), similar to the intact locus. This indicates that the increased resistance of the wild-type over ΔREIIIHSS is conferred by REIII. Together these results indicate that REIII is required for epigenetic stability at MAT.

REIII imposes epigenetic behavior under environmental stress conditions

We next studied how REIII contributes to epigenetic stability in the context of a physiological perturbation, such as change in ambient temperature. Consistent with previous reports, we found that ncRNA-nucleated spreading is sensitive to continuous growth at high temperature, likely due to the cytosolic shuttling of RNAi-components (Woolcock et al., 2012; Figure 4—figure supplement 2A). WT MAT behaved in a similarly sensitive manner. In contrast, heterochromatin in ΔKHSS cells was highly resistant to elevated temperature (Figure 4—figure supplement 2A).

We next probed the ability to remember the heterochromatin state after a transient exposure to elevated temperature, by exposing cells to 38˚C for 10 doublings, followed by return to growth at 32˚C (Figure 4E). As expected from our steady-state experiments above, REIII-dependent heterochromatin (ΔKHSS cells) is only minimally affected by the perturbation and regains full spreading rapidly (Figure 4F, Figure 4—figure supplement 2F), whereas WT MAT and ncRNA-nucleated (ΔREIIIHSS) strains lose a significant amount of spreading (Figure 4F, Figure 4—figure supplement 2D,E) and nucleation (Figure insets). Both strains regain nucleation at cenH rapidly (1 day after return to 32˚C; Figure 4—figure supplement 2B,C). However, they are discrepant in their kinetics of restoration to the 32˚C extent of spreading, with WT MAT recovering much more rapidly than the strain nucleated exclusively by ncRNA (∆REIIIHSS) (Figure 4F). Indeed, plot fitting reveals a half-life (t1/2, time to reach 50% of initial state) difference of ~22 hr, or ~9–10 generations between WT MAT and ΔREIIIHSS (Figure 4F). Therefore, REIII- is required for efficient recovery to the fully repressed state after heat perturbation. These data suggest that a central role of REIII is to ensure that epigenetic stability at MAT is maintained in the face of environmental perturbations in the wild.

Stability of heterochromatin in the absence of cenH and REIII trans-acting factors

To address dependence of the epigenetic maintenance of spreading on nucleation following heterochromatin establishment, we examined the behavior of cells following the removal of trans-acting factors required for the initial recruitment of nucleation factors such as Clr4, Swi6/HP1 and HDACs. This experiment is similar to the induced removal of the cis-acting sites in S. cerevisiae (Cheng and Gartenberg, 2000). ∆REIIIHSS and ∆KHSS-OFF isolate cells (see above, derived from nonselective plating of ∆KHSS) with established heterochromatin were crossed to mutants disrupting recruitment of nucleation factors at each element (Figure 5A). To impair REIII, we crossed the ∆KHSS-OFF reporter strain to ∆pcr1 (Noma et al., 2004). To impair ncRNA nucleation, we crossed the ∆REIIIHSS reporter strain to seb1-1, a mutant allele of the Seb1 RNA binding protein. Seb1 functions redundantly with the RNAi pathway in ncRNA nucleation, including binding cenH transcripts, and the mutant allele seb1-1 has defects in triggering nucleation at dh and dg pericentromeric elements (Marina et al., 2013). We focus on Seb1, as RNAi pathway mutants have little discernable effect on MAT when introduced after establishment (our unpublished data and [Hall et al., 2002]), indicating a stronger role for Seb1. Identifiable ∆REIIIHSSseb1-1 and ∆KHSS-OFF∆pcr1 colonies were grown for flow cytometry analysis immediately following mating and selection. The control cross mutant strains ∆REIIIHSS∆pcr1 and ∆KHSS-OFFseb1-1∆dcr1 (loss of all ncRNA-nucleation [Marina et al., 2013]) allowed us to assess any effects the trans-factor may have even in the absence of its cognate site of action.

Figure 5 with 1 supplement see all
Differential inheritance of ncRNA-dependent and independent spreading in the absence of nucleation factors.

(A) Scheme for removal of Pcr1 (REIII binding factor) in the ΔKHSS strain OFF isolate (ΔKHSS-OFF). Progeny of the cross was selected for ΔKHSS-OFFΔpcr1 genotype and identifiable colonies immediately grown for cytometry, and passaged for 456 hr. (B) Stable inheritance of repression in ΔKHSS-OFFΔpcr1. ΔKHSS-OFFΔpcr1 or ΔKHSS-OFF cells (dark blue lines) where analyzed by flow cytometry over consecutive days, the break indicating passaging without analysis. Δpcr1 had no significant effect on ΔREIIIHSS (light blue lines). (C) LEFT: scatter plots with partial point transparency of ΔKHSS-OFF or ΔKHSS-OFFΔpcr1 early and late in the time course. RIGHT: In the middle of the time course (asterisk in (B)), ΔKHSS-OFFΔpcr1 were struck for single colonies. The scatter plots for one of the isolates is shown. (D) Scheme for removal of functional Seb1 in ΔREIIIHSS strain. Selection and growth as in A., total passaging time 96 hr. (E) Weak inheritance of repression in ΔREIIIHSSseb1-1 (light blue lines). Analysis as above, total time course 96 hr. Removal of both Seb1 and RNAi pathways (ΔKHSS-OFFseb1-1Δpcr1) does not affect maintenance of silencing (dark blue lines). (F) Scatter plots of ΔREIIIHSS at 24 and 96 hr and through the entire time course for ΔREIIIHSSseb1-1. In these scatter plots, X and Y values of each cell are represented by purple dashes along the corresponding axis.

https://doi.org/10.7554/eLife.32948.018

Strikingly, most ∆KHSS-OFFpcr1 cells remains robustly repressed over 456 hr, around 200 generations (Figure 5B). However, removal of Pcr1 does have a small discernable effect, as the ∆KHSS-OFFpcr1 strain showed a small population of cells not completely in the OFF state compared to the ∆KHSS-OFF parent (Figure 5C LEFT). Further, by ~400 hr we detected a small fully ON population absent in the parent. This behavior is broadly consistent with the reported stability of intact ∆KOFF (switch rate of ~10−4 generation, Grewal and Klar, 1997; Thon and Friis, 1997), even though our assay appears to show even smaller ON populations. Very small ON populations are more apparent in a growth selection based assay as only the targeted population survives, as opposed to our assay, which captures all cells. We note a formal possilibty remains that selection of OFF colonies yields higher apparent stability. To get a closer view of the behavior of individual isolates from the population, after 288 hr of continuous passage, we streaked for single colonies and measured the resulting populations. While 5/6 isolates behaved like the broader population, we found 1/6 isolates that experienced more severe breakdown in its heterochromatic state (Figure 5C RIGHT). In this isolate heterochromatin collapsed in a manner not ordered with respect to REIII proximity and exhibited ‘green”ON/‘orange”OFF cells. In contrast, ∆REIIIHSSseb1-1 lost most spreading at the first measurement point (24 hr, Figure 5E) with progressively increasing de-repression of ‘orange’, but also some loss of ‘green’, over the next 72 hr (Figure 5F). This suggests that the epigenetic inheritance cenH-spreading requires continuous nucleation, at least via the Seb1 pathway, consistent with the behavior at synthetic nucleators (Audergon et al., 2015; Ragunathan et al., 2015).

REIII-, but not cenH-dependent heterochromatin suppresses histone turnover

It is known that REIII recruits the HDAC Clr3 (Yamada et al., 2005), which was later shown to repress the turnover of histones (Aygün et al., 2013). This suggested the intriguing possilibty that unstable epigenetic inheritance in the absence of REIII is linked to elevated histone turnover. To test this idea, we adopted the Recombination Induced Tag Exchange (RITE) system (Verzijlbergen et al., 2010) to assay replication-independent turnover of H3 in ∆REIIIHSS and ∆KHSS strains (Figure 6A). Tag switching (T7 for HA tag) in log phase growth was induced by administering β-estradiol concurrently with stalling replication with 15 mM hydroxyurea (HU) for 4 hr, during which time cells remain in early S phase (Figure 6—figure supplement 1). We compared the incorporation of T7 at 4 vs. 0 hr between ∆REIIIHSS, ∆KHSS-OFF and ∆KHSS-ON strains. First, we examined two euchromatic genes, pyk1 on chromosome 1, and mtd1, which is just outside the MAT locus. H3 turnover at these regions does not differ between the strains (Figure 6B) and is highest in in the strongly expressed pyk1 gene. We next examined sites in the MAT locus that are shared in sequence and genomic position between ∆REIIIHSS and ∆KHSS (probes indicated in diagram, Figure 6B). We note this includes also REIII, since this locus only differs between the strains by the 14bp containing the two Atf1/Pcr1 binding sites. In contrast to euchromatic loci, we observed that ∆KHSS-OFF experiences very low or no histone turnover at MAT targets by 4 hr HU compared to ∆KHSS-ON and ∆REIIIHSS, which experienced levels of H3 turnover more consistent with our euchromatic controls. This in not unexpected for ∆KHSS-ON, as is it effectively not heterochromatic (Figure 2E), and is consistent with previous results (Aygün et al., 2013). However, the observation that ∆REIIIHSS displays H3 exchange at levels similar to ∆KHSS-ON and euchromatin suggests that it is memory, rather than heterochromatin formation itself, that requires repressed histone turnover.

Figure 6 with 1 supplement see all
Histone turnover correlates with epigenetic stability in ncRNA-dependent and REIII-dependent heterochromatin.

(A) LEFT: Overview of the RITE system for histone 3.2. Cre recombinase allows tag exchange from HA to T7. RIGHT: experimental scheme for detecting replication-independent H3 turnover. Cells were grown to log phase and then grown for 4 hr in the presence of β-estradiol and 15 mM hydroxyurea. (B) Enrichment for H3-T7 at indicated loci in ΔKHSS-ON, ΔKHSS-OFFor ΔREIIIHSS strain. TOP: Location of amplicons for T7-ChIP indicated by bars. Dashed boxes in MAT indicated regions of genomic difference between ΔKHSS and ΔREIIIHSS. WT and m for REIII indicate presence or deletion of Atf1/Pcr1-binding sites, respectively. BOTTOM: Enrichment of T7 tag by ChIP at 4 hr in HU over 0 hr for indicated strains. one indicates no enrichment over 0 hr. Error bars indicate standard deviation of technical replicates. (C) Model for collaboration of cenH and REIII in establishing and maintaining the high fidelity MAT locus. (LEFT) During initial establishment, cenH heterochromatin raises the nucleation frequency at REIII (green arrow). A box right of REIII represents a putative additional nucleation element. (RIGHT) Labile cenH-nucleated spreading is disrupted, in part by de-stabilized nucleosomes, in a environmental perturbation or a stochastic event. REIII promotes reestablishment of the initial state by repressing histone turnover, limiting nucleosome loss (orange) and thus aiding spreading from cenH (light blue arrows, (1)), or promoting heterochromatin spreading from surrounding elements (dark blue arrows, (2)).

https://doi.org/10.7554/eLife.32948.020

Discussion

The patterning of the genome into regions of activity and inactivity underlies the formation of cellular identity. In many systems heterochromatin spreading is the dominant contributor to the pattern (Schultz, 1939; Schwartz et al., 2006; Wen et al., 2009). Maintaining identity requires the capacity to ‘remember’ the positional extent of heterochromatic spreading. Yet, how precise epigenetic memory is linked to the intrinsic properties of the spreading reaction itself has remained opaque. In this work, we were able to directly measure the heterochromatin spreading reaction in single S. pombe cells, separate from DNA-directed events at nucleation elements, and probe its behaviors and memory characteristics. The central principle that emerges form this work is that heterochromatin spreading in fission yeast, driven predominantly by ncRNA elements, is epigenetically unstable and requires stabilization by accessory elements for high fidelity epigenetic inheritance. At the MAT locus, which carries cell identity information, a separate type of heterochromatin, independent from nc-RNA elements and dependent on the REIII element, safeguards epigenetic propagation by repressing histone turnover.

ncRNA-triggered spreading is epigenetically unstable and labile in the face of perturbations

The dominant form of heterochromatin in S. pombe, triggered by ncRNA nucleators, leads to stochastic spreading of both silencing and H3K9 methylation that only occurs in some cells, and forms intermediate states (Figure 1 and Figure 1—figure supplement 1E, Figures 2C and 3D). This is consistent with position effect variegation in genetically disrupted systems (Elgin and Reuter, 2013; Nimmo et al., 1994). Additionally, the linear distance-dependent behavior we observe (Figure 1C) is reminiscent of the continuous spreading model in S. cerevisiae telomeres (Renauld et al., 1993; Talbert and Henikoff, 2006). This behavior of ncRNA spreading is not due to weak nucleation, as repressive histone marks accumulate to the same high extent at cenH in wild-type and ΔREIII and at Atf1/Pcr1 proximal region in ΔK cells. (Figures 1E and 2D and E). In a key result, we find ncRNA-triggered spreading to be epigenetically unstable. This is evidenced by highly dynamic behaviors over time and across generations, little discernable memory, and low resistance to chemical or environmental perturbations (Figures 35). Those behaviors are not necessarily predicted by the stochastic induction spreading, given that PEV in flies results in clonally inherited patches (Elgin and Reuter, 2013). This result opens the question how high fidelity can be achieved with ncRNA nucleators at loci that carry critical cell type specification information. The most likely cause for this instability is elevated and near-euchromatic levels of histone turnover (Figure 6B). This implies that while elevated histone turnover is compatible with heterochromatin formation per se, it is incompatible with epigenetic memory.

In contrast to the behavior of ΔREIII, ΔK cells, dependent on REIII for heterochromatin formation (Figure 2—figure supplement 2A&B), do not display stochasticity in spreading (Figures 2E and 3E), and instead repress MAT uniformly across nucleated cells in the population (Figure 2E). Under environmental perturbation, ΔK heterochromatin is extraordinarily resistant (Figure 4F and Figure 4—figure supplement 2A,F) and capable of high memory retention, even in the absence of the REIII-targeted Pcr1 protein, which attracts HDACs and Clr4/Swi6 (Jia et al., 2004; Kim et al., 2004) (Figure 5B,C). This is consistent with previously documented bistable behaviors ascribed to the overall locus (Dodd et al., 2007; Grewal and Klar, 1996). It, however, remains to be resolved whether heterochromatin in ΔK cells is nucleated by REIII and spreads outwards, or is nucleated at multiple sites, yielding apparent uniform heterochromatin formation. A REIII nucleated spreading model is favored by results presented here and by others (Jia et al., 2004; Wang and Moazed, 2017) that demonstrate that Atf1 and Pcr1 proteins or their binding sites are absolutely required for established of heterochromatin in ΔK cells, yet not for ΔREIII cells. However, unlike for cenH, where sufficiency has been clearly shown (Hall et al., 2002), we and others (Wang and Moazed, 2017) do not document significant heterochromatin formation by REIII when it is placed ectopically (Figure 2—figure supplement 1D). Thus, it cannot be differentiated whether the uniform heterochromatin formation in ΔK is the result of cooperation between different yet-to-be identified cis-acting elements, or a special property of REIII-driven spreading. Single site driven spreading of this ‘all or none’ type could be the result of looping, invoked for the polycomb system (Bantignies and Cavalli, 2011), predicted to improve spreading efficiency and memory in fission yeast (Erdel and Greene, 2016), or a unique molecular signature at REIII. For example, REIII recruits the HDAC Clr3 (Yamada et al., 2005), which promotes accumulation of the H3K9 trimethyl state, required for efficient spreading by Clr4 (Zhang et al., 2008; Al-Sady et al., 2013; Jih et al., 2017).

REIII stabilizes heterochromatin spreading by repressing histone turnover

Regulation of histone turnover has been linked to epigenetic memory in fission yeast (Taneja et al., 2017) and has been previously shown to be low at wild-type MAT (Aygün et al., 2013). Hence, the high histone turnover we observe in ΔREIII cells results from unaided cenH-spreading. REIII recruits the HDAC Clr3 (Yamada et al., 2005), which represses histone turnover (Aygün et al., 2013). Our finding that the ΔKOFF allele features very low histone turnover (Figure 6B), similar to the wild-type locus (Aygün et al., 2013), is thus consistent with REIII acting to repress histone turnover, when in a heterochromatic state. The extraordinary memory of repression we observe in ΔKOFF likely is explained by this repressed turnover, although we should note it is possible that isolation of ΔKOFF alleles, while consistent with the literature (for example [Grewal and Klar, 1996; Thon and Friis, 1997]), could bias the population to enhanced inheritance of repression. We speculate that reduced turnover increases retention of H3K9me3 nucleosomes, promoting methylation across nucleosomes by Clr4 via its H3K9me-dependent ‘read-write’ functionality (Al-Sady et al., 2013; Jih et al., 2017; Ragunathan et al., 2015; Zhang et al., 2008), thus facilitating re-establishment in the next generation. H3K9me3 is also directly promoted by Clr3, which is recruited to REIII (Yamada et al., 2005), further favoring reestablishment of methylation.

Collaboration of ncRNA-dependent and independent mechanisms in the maintenance of MAT heterochromatin

Repression of histone turnover and resulting epigenetic stability in ΔK strains requires cells to first nucleate and adopt a heterochromatic state (ΔKOFF, Figure 6B). However, since ΔK cells only nucleate infrequently (Figure 2E), how is REIII able to stabilize heterochromatin in most wild-type MAT cells (Figures 2B, 3B,C and 4B)? The independent action of cenH and REIII elements cannot account for this behavior, hence they must collaborate. We propose that in the context of wild-type MAT, cenH stimulates REIII nucleation (model, Figure 6C). Recent findings indicate that Atf1/Pcr1 are present at REIII even in non-silenced ΔK-type cells (Wang and Moazed, 2017). We speculate that since Atf1/Pcr1 recruits silencing factors such as Clr4 and HDACs (Jia et al., 2004; Kim et al., 2004; Yamada et al., 2005), heterochromatin originating from cenH might stabilize this recruitment. This hypothesis is supported by our observation for nucleation during TSA recovery. Although ΔKHSS cells very rarely renucleate (Figure 4D), REIII at the intact MAT locus must be active in most cells, as the heterochromatin reformed after erasure has much higher resistance to perturbation than that nucleated from cenH alone (red lines in Figure 4B and C).

Activated REIII in turn stabilizes the MAT locus most prominently when the heterochromatin state is perturbed. We infer this from the difference between the initial challenge and recovery from growth at high temperatures. When initially challenged, heterochromatin spreading at the wild-type MAT locus resembles that of ncRNA-nucleated heterochromatin (Figure 4—figure supplement 2A), suggesting that REIII or other nearby elements plays a minor role under normal circumstances at MAT. However, the heat recovery experiment suggests that changes in the REIII–dependent heterochromatin stabilization or assembly, not cenH nucleation (Figure 4—figure supplement 2B and C), takes on a major role in the accelerated recovery of the collapsed heterochromatin (Figure 4F). Thus, REIII is required under perturbation conditions to protect or quickly re-establish the heterochromatin state (Figure 4B,F and model Figure 6C). The relatively transient distal de-repression events experienced by wild-type MAT cells, which are much more pronounced in ΔREIII cells (Figure 3C and D), further points to REIII acting after stochastic loss of cenH spreading in steady state. It is possible that REIII does so by stabilizing existing heterochromatin via repression of histone turnover when the loss of these structures is sensed, or alternatively, that REIII-dependent structures expand or ‘fill-in’ the void left by collapse of cenH-spreading. In either case, we propose that REIII acts as a failsafe, ensuring the integrity, and ultimately epigenetic memory, of heterochromatin at MAT through perturbations.

In summary, we propose a model whereby the division of labor between cenH and REIII is uniquely suited for a cell type specification locus such as MAT, which requires silencing that is both robust and intergenerationally stable. ncRNA-nucleation is extremely robust but intrinsically too labile and stochastic to reliably control the cell type specification locus, thus requires support from an accessory element. The need for reliable control of cell type specification loci is shared in more complex systems. The nature of equivalent accessory elements to REIII and how they act in these cases remains to be determined.

Materials and methods

Strain construction

Plasmids and strain selection

Request a detailed protocol

Plasmids to generate constructs for genomic integration were generated by standard methods including Gibson assembly and in vivo recombination. S. pombe transformants were selected directly on dropout media for auxotrophic markers or onto rich media (YES) for 24 hr followed by selective media (YES + G418, YES + hygromycin or YES + nourseothricin). For all strains see Table 1.

Ura4 replacement method

Request a detailed protocol

To avoid interference of selection cassettes with heterochromatin function in our HSS, we produced ‘scarless’ genomic integrations, lacking selection markers. To do so, we marked the insertion site first with a ura4 cassette by genomic integration and then replaced this cassette either with a XFP cassette or altered genomic sequence for site mutations. ura4 replacements were isolated by 5-FOA counter-selection and confirmed by genomic PCR. This method was used to generate the atf/creb site deletions and sequence insertions. ura4 was targeted to the region between Mat3M and cenH, specifically including the two seven base atf/creb-binding sites (s1 and s2, and [Wang and Moazed, 2017]). The entire ura4 cassette was then replaced with a construct containing the two seven base pair deletions of s1 and s2 or a deletion of s1 with additional 700 bp of sequence from the sib1 open-reading frame. Desired point mutations and restoration of the pre-substitution locus was confirmed by PCR and sequencing.

Flow cytometry and FACS sorting

Request a detailed protocol

For standard flow cytometry experiments, cells were grown overnight in rich media (YES) and then diluted in the morning to OD = 0.1 in minimal media plus supplements (EMM complete) and grown 4–6 hr before analysis by flow cytometry. Flow cytometry was performed using Fortessa X20 Dual or LSRII instruments (Becton Dickinson, San Jose, CA). Samples sizes ranged from ~2000 to 100,000 cells depending on strain growth. Compensation was performed using cells expressing no XFPs and single-color controls expressing 1 XFP each. Compensated data was used for all downstream analyses. Fluorescence was detected for each color as described (Al-Sady et al., 2016).

For FACS sorting experiments, cells were grown overnight from OD = 0.025 in YES and then in the morning concentrated into a smaller volume to achieve a flow rate of ~5000 events/second on the cytometer. Sorting was performed using either Aria2 or Aria3u machines (Becton Dickinson). Prior to sorting cells were strained through a 35–40 μm mesh (Corning) to reduce clogs. Sorting criteria included a gate for size (forward (FSC) and side (SSC) scatter), removal of doublets, a gate for ‘green'OFF (‘green’ signal within the range of an unstained control) and then gated into Low, Intermediate, High ‘orange’ signal defined by the following: Low encompassed signal overlapping that of an unstained control and High encompassed signal overlapping that of the Δclr4 no heterochromatin control strain PAS355. Intermediate gate was set in between Low and High with about 100 fluorescence units of a gap (representing ~2% of the full range of captured fluorescence) to ensure reliable separation. The entire range of fluorescence detected was ~2.5 orders of magnitude. At least 8 × 106 cells were collected for each population for Chromatin Immunoprecipitation and 2 × 106 cells for RT-qPCR. Immediately after sorting, the final populations were subjected to the appropriate treatment for either Chromatin Immunoprecipitation or RT-qPCR. The R scripts for analysis is included as a text file, ‘Source Data 1’.

Sytox green staining and cell cycle analysis

Request a detailed protocol

Cell cycle analyses were performed essentially as described (Knutsen et al., 2011). Briefly, cells were fixed with 70% ethanol, washed with 20 mM EDTA pH 8.0, and treated with RNaseA for 3 hr at 37°C. Immediately before analysis by flow cytometry, 2 μM Sytox Green (Invitrogen) in 20 mM EDTA pH 8.0 was used to resuspend pelleted cells. Cells were excited with a 488 nm laser and Sytox Green signal was detected with a 505-nm longpass filter and a 530/30 bandpass filter. Cell cycle analysis was performed in the FlowJo Software (Tree Star Inc, Ashland, OR) The identification of cell populations and fraction of cells in each cell cycle phase (G2, S, and G1 + M) were determined as described (Knutsen et al., 2011).

Trichostatin A (TSA) gradient experiment

Request a detailed protocol

Cells were taken from fresh plates, and then grown overnight with shaking (Elmi) in 96-well plates containing 150 μL YES (Day −1). The next day (Day 0), cells were diluted into YES and measured by cytometry. At the end of Day 0, cells were passaged into YES + DMSO (0 μM TSA) or YES + 50 μM TSA overnight. The next day (Day 1), cells were diluted and grown briefly into the same pretreatment conditions and the 50 µM TSA pre-treated cells were checked for complete de-repression by flow cytometry. Complete de-repression was defined as a qualitative overlap of WT and Δclr4 profiles, with no evidence of repression. Both 0 and 50 μM TSA pretreated cells were then diluted into a gradient of TSA of 11 two-fold dilutions from 50 μM along with a twelfth 0 μM (DMSO) point. Cells were measured after ~6 hr and then passaged into the same TSA gradient conditions to continue growth.

The next day (Day 2), cells were diluted from overnight growth into the same gradient as above, measured ~6 hr later by flow cytometry and passaged into the same gradient again overnight. The same protocol was followed for Days 3 and 4. The full experiment was performed twice at different times (biological replicate). Given the lengthy continuous growth, contamination was occasionally observed in <1% of wells. The replicate shown was chosen based on lacking contamination.

Heat recovery experiment

Request a detailed protocol

Cells were taken from fresh plates, and then grown overnight with shaking (Elmi) at either 32°C or 38°C (Day-1) in 96-well plates containing 200 µL YES medium per well. In the morning, cells were diluted into 200 µL YES and grown ~6 hr at the same temperature before measurement by flow cytometry (Day 0). At the end of Day 0, all cells were all diluted again into YES and grown at 32°C. The next day (Day 1), cells were diluted from overnight growth into YES at 32°C, measured ~6 hr later by flow cytometry and passaged into the same temperature overnight. The same protocol was followed for Days 2, 3, and 4.

Nucleation factor removal experiment

Request a detailed protocol

HSS strains were crossed to parent strains lacking functional nucleation factors for REIII (Δpcr1) or cenH (Δdcr1 seb1-1). Cross progeny were identified via a random spore approach by growth on selective media 2 or 3 days after plating. Absence of pcr1 or dcr1 open-reading frames was confirmed by PCR. Presence of seb1-1 allele was confirmed by sequencing. Single colonies were grown in 96-well plates at 32°C containing 200 µL YES medium per well. In the morning, cells were diluted into 200 µL EMM and grown ~6 hr at the same temperature before measurement by flow cytometry. Cells were again diluted into 200 µL YES for overnight growth at 32°C and grown and measured similarly the subsequent days. For Δdcr1 and/or seb1-1 strains and their controls, this was continued for four days. For Δpcr1 strains and their controls this was continued for 5 days then resulting cells were plated onto selective media and allowed to grow 48 hr at 32°C. Patches were then passaged in bulk on selective plates every 36–48 hr for 7 additional days. On the 6th day, the passaged ΔKHSSΔpcr1 cells were additionally struck for singles. On the 8th day, patches of passaged cells and six single colonies of ΔKHSSΔpcr1 cells were grown in 96-well plates as above and measured by flow cytometry for 5 additional days.

Nucleosome turnover assay

Request a detailed protocol

Recombination Induced Tag Exchange (RITE) parent strain (HU2549) was crossed into HSS reporter strains. Resulting isolates were verified by growth on selective media. The cdc-25ts allele was crossed out. RITE was performed essentially as described (Audergon et al., 2015; Svensson et al., 2015) with the following exceptions. Given the labile nature of heterochromatin at elevated temperatures, replication stalling was performed with hydroxyurea as published (Aygün et al., 2013). Cells were grown to saturation overnight in YES supplemented with Hygromycin. In the morning cells were diluted to OD = 0.1 in 50 mL YES+Hygromycin and grown for 4 hr at 30°C, 225 rpm. After 4 hr of growth, 13 mL of cells were pelleted and processed for ChIP as the 0 hr time point. The remaining cells were washed twice in media devoid of Hygromycin and finally resuspended in YES supplemented with 15 mM Hydroxyurea (HU) and 1.5 μM β-Estradiol (ER) and incubated for 4 additional hours at 30°C, 225 rpm. After 4 hr incubation with HU and ER, 10 mL of cells were pelleted and processed for ChIP.

Chromatin immunoprecipitation (ChIP) and quantification

Request a detailed protocol

We found that sonication of a small number of cells such as can be collected by FACS leads to a marked increase in background signal from negative control regions that was absent when ChIP was performed with larger log phase cultures (>50 × 106 cells). To address this, ChIP in Figure 1E was performed on each of the FACS sorted populations with the addition of 42 × 106 formaldehyde fixed cells of S. cerevisiae W303 strain as a carrier. ChIP in Figure 2D was performed with 15 × 106 cells of each fission yeast strain and 50 × 106 additional W303. ChIPs for Figure 2E and Figure 2—figure supplement 2B were performed with 80 × 106 cells and no added W303. ChIPs for Figure 6B were performed with no added W303. ChIP was additionally performed on a sample of W303 alone, which only produced signal equivalent to background. S. pombe ChIP samples and W303 cells were fixed and pre-processed for ChIP separately, then mixed together immediately prior to lysis. Cells were cross-linked and lysates prepared for ChIP as described (Canzio et al., 2011) with the following exceptions: After lysis, the chromatin fraction was resuspended in 350 μL lysis buffer and sonication performed using a Diagenode Bioruptor Pico machine at 4°C, with 16–20 rounds of 30 s ON, 30 s rest. ChIP was essentially as described, with the total lysate split into 2–6 equal volumes (after ~8% set aside as input fraction) and ChIP performed in 600–800 μL per sample. Two or three technical replicates were performed across experiments. 1 μL of each of the following antibodies was added per ChIP replicate: anti-H3K9me2 (Abcam ab1220); anti-H3K4me3 (Active Motif 39159); anti-H3K9me3 (Millipore 07–442); anti-T7 (Novagen 69522–3). ChIP samples were agitated on a Nutator overnight at 4°C. Immune complexes were collected for 3 hr with 15–20 μL washed protein A Dynabead slurry (Invitrogen). Washing and downstream processing steps were essentially as described, except ‘wash buffer’ wash was performed once. Samples were purified using a Machery-Nagel PCR purification kit and NTB buffer for SDS containing samples. DNAs were quantified by RT-qPCR (see below). Enrichments were calculated as follows: For Figures 1E, 2D and E IP/input values for amplicons of interest were calculated and normalized to the IP/Input values for positive controls for each antibody, dh for H3K9me2 and H3K9me3 and the actin promoter for H3K4me3. For Figure 2—figure supplement 2B, ChIP signal was normalized to signal from a matched background Δclr4 strain. For Figure 6B IP/input values for the 4 hr time points were normalized to the IP/input values for the 0 hr time point.

RNA extraction and mRNA quantification

Request a detailed protocol

After sorting, samples were spun at 5000xg, supernatant decanted, and pellets flash frozen in liquid nitrogen and stored at −80°C. For the Δclr4 strain PAS335, cells were grown into log phase and then cell pellets were isolated in the same fashion. Total RNA was extracted in technical duplicates from the same cell pellets using the ‘MasterPure- Yeast RNA Purification Kit’ (Epicentre), including a 30 min DNAse treatment step post-RNA isolation. Reverse Transcription was performed with SuperScript III RT (Invitrogen), using the supplied protocol and 1.5–2 μg of RNA and an oligo dT primer. Following cDNA synthesis the reaction was treated with RNAse H (New England Biolabs). cDNA samples were quantified by RT-qPCR. For each sorted sample, mKO2 cDNA values were normalized to actin and then divided by the max value calculated similarly from PAS355 (Δclr4).

RT-qPCR

Request a detailed protocol

Real-time quantitative PCR was performed using a BioRad CFX-384 machine. 15 μL reactions were prepared, each containing 7.5 μL of Applied Biosystems SYBR Select Master Mix, 4.5 μL 3.3M betaine, 1.2 μL of 2.5 μM oligo mix, 0.8 μL water, and 1 μL template. The thermocycler protocol was: 2 min at 50°C then 2 min at 95°C followed by 40 cycles of 15 s at 95°C and then 1 min at 60°C followed by a plate read. Lastly a melt curve was generated. Standards were generated with five fold dilutions of genomic DNA containing templates for all PCR products.

Single-cell microscopy

Request a detailed protocol

Single cells of strains PAS 387, 389, 391 and 244 (see Table 1; E2Crimson under act1 promoter) were captured in microfluidic devices as described (Spivey et al., 2017). Multi-channel fission yeast lifespan microdissectors (multFYLM) contained six independent devices (channels), each of which is capable of capturing up to 392 cells (https://bio-protocol.org/e2783). In brief, the devices were cast in polydimethylsiloxane (PDMS, Sylgard 184, Dow Corning) using conventional soft lithography methods. Master structures were fabricated from P-doped silicon wafers (ID#452, University Wafers) and SU-8 photoresists 3005 and 2010 (Microchem, Westborough, MA). MultFYLMs were cleaned and adhered to glass coverslips (48 × 65 mm #1, Gold Seal), and then connected to syringes (60 mL, Becton-Dickson) containing YES 225 liquid media (Sunrise Science) via PFA tubing and microfluidic fittings (IDEX Health and Science). The multFYLM was maintained at 30˚C in a custom staged-mounted environmental chamber on an inverted microscope (Eclipse Ti, Nikon) equipped with NIS Elements software (Nikon), a 60X air objective (CFI Plan Apo λ, 0.95 NA, Nikon) fitted with an objective heater (Bioptechs), a motorized stage (Proscan III, Prior), and an active feedback-based focusing system (Perfect Focus System, Nikon). An LED lamp (Sola II, Lumencorp) and a scientific-grade CMOS camera (Zyla 5.5, Andor) were used for fluorescent imaging. Multi-color fluorescent imaging of sfGFP, mKO2 and E2Crimson fluorophores was carried out by alternating between three filter sets mounted in a computer-controlled filter ring (Chroma 49002, 49010 and 49015, respectively). To help with the semi-automated cell identification, each channel was imaged every ten minutes via brightfield imaging (100 ms exposure, both in focus and 4 μm below the focal plane). Fluorescent images of each of the three fluorophores were taken every 30 min (150 ms exposure). This illumination scheme was well below the phototoxicity limit, as described previously (Al-Sady et al., 2016). Raw images were saved as uncompressed 16 bit ND2 files and further analyzed using a custom-written image analysis pipeline (see below).

Cells were grown overnight (30°C with 225 rpm shaking) to saturation in YES media, then diluted in YES to an optical density at 600 nm (OD600) of 0.1 and allowed to grow for approximately 5 hr to reach an OD600 of 0.5. Cells (60 μL at OD 0.5 in YES + 2% Bovine Serum Albumin, BSA) were loaded at the entry port of the multFYLM. After cells entered individual channels, media lines were reattached and YES media was pumped through on a pulse cycle (14 min: 5 µLmin−1, 1 min: 55 µLmin−1) for the entire experiment. This flow regime was optimized to flush out occasional cell clumps that grew at the device inlets and other fluidic interfaces. Four genotypes were imaged simultaneously for 60 hr in each channel of a multFYLM device to ensure identical imaging and growth conditions. In all cases, we only analyze the innermost cell, which was the oldest cell pole (see below). Cells that were ejected or died within the first 12 hr after loading were not included in the downstream analysis.

Single-cell image analysis

Request a detailed protocol

Single-cell imaging data was processed using an updated version of the custom-written FYLM Critic analysis package (Spivey et al., 2017). The source-code is available via GitHub (https://github.com/finkelsteinlab/fylm; Rybarski et al., 2015; copy archived at https://github.com/elifesciences-publications/fylm). FYLM Critic performs the following automated processing on the raw images: (1) rotation; (2) jitter removal via a cross-correlation algorithm; and (3) generation of kymograph and individual cell images. The latter were used to create videos of individual cells in Fiji (Schindelin et al., 2012). The final outputs of FYLM critic are the position and contour of each dividing cell, as well as the time-dependent fluorescence intensities for each cell. These fluorescence intensities are obtained by averaging the intensity across all pixels that fall within the cell volume, as defined by the bright-field images. This normalization also ensures that the fluorescence intensity is corrected for the size of the rapidly dividing cells. Time-dependent fluorescent intensities were analyzed via custom-written MATLAB scripts (version 2017a Mathworks, available upon request). Background fluorescence from the PDMS device was subtracted using catch tubes that did not receive a cell. The maximum heterochromatin reporter (GFP, mKO2) fluorescence intensity was calculated using ∆clr4 cells in the same reporter construct background. To control for expression variation across the cell cycle, the fluorescence from heterochromatin reporters was also reported as a ratio of the control fluorophore, E2Crimson. Similarly, cells fluorescing in the clamp channel were removed from analysis for MAT-locus-derived strains (see Appendix 1-Supplemental Materials and methods).

Single-cell images generated by the FYLM Critic analysis were compiled into stacked movies using Fiji. Images in bright field and for each color channel were processed separately in batch and then later combined into a vertical stack. For each channel, 0.2% of pixels were allowed to become saturated and pixel values were normalized to the maximum range for the whole sequence in that channel. For bright field, every third image was included to match the imaging frequency of the fluorescent channels. Movies were edited for length to include contiguous imaging sequences without loss of focus and for size to remove non-cellular debris and cells from the opposite side of the channel that entered the field of view. After combining all color channels and bright field, the brightness and contrast were increased for cell 407 to match the red channel brightness of the other strains. Image sequences were saved as uncompressed. avi files with a rate of 15 frames/s.

Table 1
Yeast strains used in this study.
https://doi.org/10.7554/eLife.32948.022
StrainGenotype
PAS075Locus2::ade6p::3xE2C:hygMX at Locus2 (between SPBC1711.11 andSPBC1711.12)
PM03Wild-type strain: h(+); ura4-D18; leu1-32; ade6-M216; his7-366
PM1035ura4::natMX:dh fragment 1, clr4::KAN as in Marina et al. (2013)
PAS111ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 7 kb, ade6p:3xE2C: hygMX at Locus2
PAS112ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 7 kb, ade6p:3xE2C: hygMX at Locus2; clr4::kanMX
PAS133ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 1 kb, ade6p:3xE2C: hygMX at Locus2; clr4::kanMX
PAS134ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 1 kb, ade6p::3xE2C: hygMX at Locus2
PAS135ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, ade6p::3xE2C: hygMX at Locus2; clr4::kanMX
PAS136ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, ade6p::3xE2C: hygMX at Locus2
PAS141ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 5 kb, ade6p::3xE2C: hygMX at Locus2
PAS142ura4::natMX:dh:ade6p: SF-GFP, ade6p:mKO2 5 kb; ade6p::3xE2C: hygMX at Locus2; clr4::kanMX
PAS192ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2, h(-)
PAS193ΔK::ade6p:mKO2; ade6p:SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2; clr4::kanMX, h(-)
PAS214ΔK::ade6p:mKO2:ura4t; mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2; clr4::kanMX, h(-)
PAS215ΔK::ura4t:mKO2:ade6p; mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2; clr4::kanMX, h(-)
PAS216cenH::ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ade6p:3xE2C: hygMX at Locus2; clr4::kanMX, h90
PAS217cenH: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ade6p:3xE2C: hygMX at Locus2, h90
PAS218cenH::ade6p:mKO2 (Kint2); mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2; in clr4::kanMX, h90
PAS219cenH: ade6p:mKO2 (Kint2); mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2, h90
PAS231ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, leu1::ade6p:3xE2C: hygMX
PAS237ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, act1p::qxE2C: hygMX at Locus2; clr4::kanMX
PAS243ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, act1p::1xE2C: hygMX at Locus2; clr4::kanMX
PAS244ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, act1p::1xE2C: hygMX at Locus2
PAS264cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ade6p:3xE2C: hygMX at Locus2, pcr1::kanMX, h90
PAS268ΔK:: ade6p:mKO2; ade6p:SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2, REII::LEU2, h(-)
PAS269ΔK:: ade6p:mKO2; ade6p:SF-GFP between REIII and mat3M; ade6p:3xE2C:hygMX at Locus2; clr4::kanMX, REII::LEU2, h(-)
PAS331cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2;
ade6p:3xE2C:hygMX at Locus2; ΔREIII::REIII(Δs1, Δs2) in clr4::kanMX, h90
PAS332cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ade6p:3xE2C:hygMX at Locus2; ΔREIII::REIII(Δs1, Δs2), h90
PAS348ura4::hygMX:REIII:ade6p:SF-GFP; ade6p:mKO2 5 kb, ade6p:3xE2C:natMX at Locus2
PAS350ura4::hygMX:REIII:ade6p:SF-GFP; ade6p:mKO2 5 kb, ade6p:3xE2C:natMX at Locus2 dcr1::kanMX
PAS355ura4::natMX:dh:ade6p:SF-GFP, ade6p:mKO2 3 kb, leu1::ade6p:3xE2C:hygMX; clr4::kanMX
PAS385ΔK:: ade6p:mKO2; ade6p:SF-GFP between REIII and mat3M; act1p:1xE2C:hygMX at Locus2; clr4::kanMX, h(-)
PAS387ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; act1p:1xE2C: hygMX at Locus2, h(-)
PAS388cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ura4 at Locus2; leu1::act1p:1xE2C:hygMX, clr4::kanMX, h90
PAS389cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ura4 at Locus2; leu1::act1p:1xE2C:hygMX, h90
PAS390cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ura4 at Locus2;
ΔREIII::REIII(Δs1, Δs2), leu1::act1p:1xE2C:hygMX, in clr4::kanMX, h90
PAS391cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ura4 at Locus2;
ΔREIII::REIII(Δs1, Δs2), leu1::act1p:1xE2C:hygMX, h90
PAS398his1::natMX:dh:ade6p:mKO2; ade6p:SF-GFP 3 kb, ade6p::3xE2C:hygMX at Locus2, clr4::kanMX, ura4::phyB.
PAS399his1::natMX:dh:ade6p:mKO2; ade6p:SF-GFP 3 kb, ade6p::3xE2C:hygMX at Locus2, ura4::phyB.
PAS410ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M;
ade6p::3xE2C:hygMX at Locus2, natMX:clr4+, h(-); ‘OFF’ allele
PAS411ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M;
ade6p::3xE2C:hygMX at Locus2, natMX:clr4+, h(-); ‘ON’ allele
PAS464ΔK::ade6p:mKO2:ura4t; mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2; natMX:clr4+, h(-)
PAS465ΔK::ura4t:mKO2:ade6p; mat3m(EcoRV):: ade6p:SF-GFP; ade6p:3xE2C: hygMX at Locus2; natMX:clr4+, h(-)
PAS473ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; 700 bp sib1 ORF between REIII-s1 and mKO2;
ade6p:3xE2C: hygMX at Locus2, clr4::kanMX, h(-);
PAS474ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ΔREIII::REIII(Δs1, Δs2),
ade6p:3xE2C: hygMX at Locus2, clr4::kanMX, h(-);
PAS478ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; 700 bp sib1 ORF
between REIII-s1 and mKO2 ade6p:3xE2C: hygMX, natMX:clr4+, h(-);
PAS482ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2, h(-); ‘OFF’ allele
PAS483ΔK:: ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ΔREIII::REIII(Δs1, Δs2),
ade6p:3xE2C: hygMX at Locus2, natMX:clr4+, h(-);
PAS496cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2; ade6p:3xE2C:hygMX at Locus2; ΔREIII::REIII(Δs1, Δs2),
ars1::prad15:cre-EBD:LEU2; h3.2:lox:HA:hygMX:lox:T7; h90
PAS497ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2;
ars1::prad15:cre-EBD:LEU2; h3.2:lox:HA:hygMX:lox:T7; ‘OFF’ allele, h(-)
PAS498ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M; ade6p:3xE2C: hygMX at Locus2;
ars1::prad15:cre-EBD:LEU2; h3.2:lox:HA:hygMX:lox:T7; ‘ON’ allele; h(-)
PAS508ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M;
ade6p:3xE2C: hygMX at Locus2, ‘OFF’ allele; pcr1::kanMX
PAS510cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2;
ade6p:3xE2C:hygMX at Locus2; ΔREIII::REIII(Δs1, Δs2), pcr1::kanMX
PAS514ΔK::ade6p:mKO2; ade6p: SF-GFP between REIII and mat3M;
ade6p:3xE2C: hygMX at Locus2, ‘OFF’ allele; dcr1::kanMX; seb1-1:natMX
PAS515cenH:: ade6p:SF-GFP (Kint2); mat3m(EcoRV):: ade6p:mKO2;
ade6p:3xE2C:hygMX at Locus2; ΔREIII::REIII(Δs1, Δs2), seb1-1:natMX

Appendix 1

Supplemental discussion

Alternative formal model for cenH and REIII interaction

In the main discussion, we propose that cenH stimulates REIII nucleation to account for the high proportion of the spreading in wild-type MAT cells. Another formal possibility remains that non-nucleated, Atf1/Pcr1-bound, REIII raises cenH spreading efficiency. In the ΔREIII strain, Atf1/Pcr1 binding sites have been fully deleted, it is possible that when bound but not in a heterochromatin state, Atf1/Pcr1 acts to encourage more efficient spreading out of cenH, possibly by directing the locus to a more spreading competent location or via its recruitment of HDACs.

Supplemental materials and methods

Basic three color HSS analysis in R

Please note the relevant R scripts can be found in ‘Source Data 1’. Below a description of the method.

Reading in the data

Standard flow cytometry data files (.fcs) were read in with the R package flowCORE (Bioconductor, https://www.bioconductor.org/packages/release/bioc/html/flowCore.html).

Isolating successfully nucleated cells (Figures 1 and 2)

A strain closely matched in genetic background to HSS strains but containing no XFPs was analyzed under the same flow cytometry conditions in each experiment. This ‘unstained’ control was gated for cell size in the same manner as analysis strains and both the median fluorescence and standard deviation determined in green or orange channels (the signal from the ade6p:SF-GFP, ade6p:mKO2 or ade6p:E2C transcriptional units is referred to here as ‘green’, ‘orange’ or ‘red’). A nucleation cutoff was set for a value corresponding to the median of fluorescence units in the clamp channel plus two standard deviations from this unstained control. Only cells with signal less than this value were considered for post nucleation analysis.

Normalizing to max fluorescence values from Δclr4 strains (Figures 1 and 2)

Max values in Δclr4 strains were determined by calculating the median raw fluorescence in each color channel after gating for cell size. For each cell of each strain for analysis, the signal in each channel was divided by this max value for the corresponding Δclr4 strain.

Normalizing to constitutive red signal

For each cell of each strain, the ‘green’ and ‘orange’ values were divided by the ‘red’ value.

The red- and Δclr4-normalized values range from 0 to ~1.5, where one corresponds to the Δclr4 (max) value. As this value is derived from the mean of a cell distribution, with Δclr4 cells falling above and below the mean, we plotted out to 1.5 to capture cells with ratio values above 1.0.

Hexbin (2-D Histogram) Analysis

Normalized ‘green’ and ‘orange’ values (without any nucleation cutoff applied) were plotted as hexbin (or 2-D histogram) plot where density is color-coded in grayscale. Data points within x = 0–1.5 and y=0–1.5 were isolated and a hexbin plot was generated using n = 40–45 bins along each axis. Hexbin plots were generated using the R package hexbin (https://CRAN.R-project.org/package=hexbin)

Spreading Analysis with Nucleation Clamp

Cells within the FSC/SSC gate and with signal below the cutoff value for the nucleation color were plotted in a 1-D histogram with n = 50–200 bins where the points in the middle of the histogram bins were plotted connected by a line.

Modifications for Figure 4

For every TSA concentration and pre-condition (no TSA or 50 μM) each strain was normalized to the median fluorescence of a Δclr4 strain grown under the same treatment. The cytometer settings were adjusted so an unstained control had mean fluorescence in the green channel of 102. A ‘green’ cutoff for nucleation was assigned to be 400 fluorescence units. ‘orange’ cutoff values for each analysis strain were generated by determining mean and two standard deviations in PAS217 at 0 μM TSA from the no TSA precondition normalized to the appropriate Δclr4 strain for each analysis strain. Previous analysis demonstrated both colors in PAS217 to be fully repressed, as evident in the mean for each channel. For Figure 4B,C we calculated at each TSA concentration the fraction of cells with green signal below the ‘green’ cutoff that have an orange signal below the ‘orange’ cutoff. These values were normalized to the fraction calculated for cells of that strain at 0 μM TSA from the no TSA precondition. For Figure 4D we calculated at each concentration of TSA the fraction of all cells with orange signal below the ‘orange’ cutoff, because in the 50 μM TSA pre-condition, insufficient cells exist that are below the cutoff for ‘green’ to perform above analysis. These values were normalized to the fraction calculated for cells of that strain at 0 μM TSA from the no TSA precondition.

In Figure 4F, for each strain at each temperature, we calculated the fraction of cells that had ‘green’ signal less than the ‘green’ cutoff AND ‘orange’ signal less than the ‘orange’ cutoff. These values were normalized to the fraction calculated for cells of that strain at 32°C. The cutoff values were based on calculations of two standard deviations from the mean of a red only strain normalized to Δclr4 controls. For ‘green’ and ‘orange’ these values were approximately 0.35 and 0.4 respectively so these values were used for all strains to standardize the analysis.

Modifications for Figure 5

For Figure 5B,E for each strain we calculated at each time point the fraction of cells that had ‘green’ signal less than the ‘green’ cutoff AND ‘orange’ signal less than the ‘orange’ cutoff. The cutoff values for ‘green’ and ‘orange’ were each 0.4 as in Figure 4—figure supplement 2A. Samples in Figure 5B,C were gated for size because we were able to measure 105 cells as the strains grew well. Samples in Figure 5E,F were not gated for size because cell growth was poor in seb1-1 isolates. In Figure 5C 5000 cells were plotted in a scatter plot with point transparency. In Figure 5D ~1400 cells were plotted in a scatter plot with point transparency. For Figure 5F ~2500 cells were plotted in a scatter plot with point transparency. Data in Figure 5 Supplement A were plotted as in Figure 5C and in Figure 5 Supplement B as in Figure 5F.

Fitting t1/2 for 38°C spreading recovery

To derive a t1/2, which is the time required to recover to 50% of the full spreading observed at 32°C, we fit the data to a simple sigmoidal dose-response variable slope model:

fraction all cells with full spreading= Bottom+tn*(Top-Bottom)tn+t1/2 n

where Bottom is the starting fraction of cells with full spreading at t = 0, Top = 1, t is time in hrs. n represents a Hillslope.

In Figure 4—figure supplement 2A, we calculated the fraction of cells that had ‘green’ signal less than the ‘green’ cutoff AND ‘orange’ signal less than the ‘orange’ cutoff. The cutoff values were based on two standard deviations from the mean of a red only strain normalized to Δclr4 controls. For green and orange these values were approximately 0.4 and 0.4 respectively so these values were used for all strains to standardize the analysis. These cells were gated for FSC and SSC to isolate live cells because at the highest temperature many cells had died. In Figure 4—figure supplement 2B,C the 1-D histograms for all ‘green’ values were calculated as above with the modification that no size gate was called. In Figure 4—figure supplement 2C–E the 1-D histograms and hexbin plot insets were calculated as above with the modification that no size gate was called and the ‘green’ cutoff values were generated as in Figure 4F.

FYLM data analysis

Initial data calculations

Loss of focus was identified by red fluorescence measurements below a cutoff of one standard deviation from the mean of all collected values of red for all cells. This loss of focus data was removed from analysis. Background fluorescence from the PDMS device at each time point was then subtracted using catch tubes that did not receive a cell. For each MAT strain, its matched Δclr4 strain was also imaged and a mean and standard deviation were calculated. In each strain cells were normalized to this mean Δclr4 value (defined as ‘max’) and to their own red values as in the flow cytometry data analysis. An ON gate (used in Figure 3B) for cells that reached maximal de-repression was calculated for each strain from the Δclr4 strain mean less three standard deviations.

Calculating nucleation gates

As seen by flow cytometry and visual inspection of collected movies, the vast majority of MAT locus cells have a repressed nucleation reporter (‘green’), which allowed us to formulate a very strict nucleation cutoff from the collected FYLM data itself. This cutoff was the mean plus two standard deviations of all measured values of all cells. Only cells that maintained a green signal less than this cutoff for their entire measured lifespan were included for further analysis in Figure 3. We did not apply this nucleation gate to the ectopic strain, as few cells maintained ‘green’ tightly repressed throughout their measured lifespan. Instead, we show all the cells in the traces plot and highlight in grey example cells that remain nucleated or mostly nucleated throughout their measure lifespans.

Rescaling orange to fix negative values

Due to background subtraction (see above) a significant fraction of cells experienced negative values for their adjusted fluorescence in the orange channel. To account for this, the data for the MAT strains in Figure 3 was rescaled by determining the lowest value measured (minVal) and adding the difference between that value and 0 to every time point of every cell. Values for all timepoints were then divided by 1+minVal to rescale back to 1 = max. The ectopic strain was not rescaled.

Data smoothing

For trace plots and heatmaps data was smoothed using a moving average of the two-nearest neighbor data points before and after. This number was chosen as it represents the timeframe of one cell division and is on the timescale of the expression and maturation of the XFPs used in these strains.

Traces

Individual cell traces represent the red normalized and smoothed, green and orange fluorescence data plotted over time. Traces begin and end at whatever time a cell entered or exited the channel or died. Therefore, not all traces start at x = 0 or end at x = 60. Curated example cells were also plotted as overlays using gray lines. For these curated cells similar trace plots for orange divided by green was plotted in Figure 3—figure supplement 2C-BOTTOM, D-BOTTOM, and F.

Heatmaps

Points with red values greater than 50% of the mean were removed. For cells that remained nucleated throughout their measured lifespan, up to 36 hr of measurements of normalized green or orange was plotted as a heat map from blue (0) to yellow (1) for 30 of the longest imaged cells. White gaps indicate transient loss of focus of less than 2 hr (four time points). Heatmaps were no longer plotted for any cell that had a loss of focus event for more than four time points.

X-Y fluorescence line plots

For the cells in each behavior category (Figure 3—figure supplement 1A), an X-Y plot was generated that plots the unsmoothed, red-normalized ‘orange’ vs ‘green’ values for each cell as a line. Values were normalized to the mean in ‘green’ and ‘orange’ from a matched background Δclr4 strain to set the value of 1, while background fluorescence from empty channels were set to the value of 0. Line segments are colored from blue to pink generated based on the measurement time of that cell while the behavior was observed. The first measured point is represented in blue, the last in pink, and the color values in between are divided into increments by the total measured time for the cell.

X-Y fluorescence dot plots

For one or two selected cells per strain an X-Y fluorescence dot plot was generated that plots the smoothed ‘orange’ vs ‘green’ values for every third time point imaged over its entire measured lifespan. Points are colored in a grayscale that is generated based on the measured lifespan of that cell. The first measured point is represented in black, the last in white and the number of remaining points set by the total measured lifespan of the cell. In Figure 3—figure supplement 1B X-Y fluorescence dot plots were generated for four example category 1 cells. In this figure, time is colored from blue to pink as in the X-Y line plots and points represent 30 min intervals beginning when both colors are fully expressed and ending when both reach a local minimum.

Cell fate pie charts

The number of measured time points for each cell was determined and converted to hours (one fluorescence image every 30 min). Cells were then binned into lifetime groups of < 12 hr, 12–36 hr, or > 36 hr. Within these bins the cells were separated based on whether they died as annotated in FYLM Critic or if their traces were cut short due to late entry into or ejection from the catch channel.

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files. The code for analyzing live cell data is included in the submission. All reagents generated in this work are available upon request.

References

    1. Cheng TH
    2. Gartenberg MR
    (2000)
    Yeast heterochromatin is a dynamic structure that requires silencers continuously
    Genes & Development 14:452–463.
    1. Grewal SI
    2. Klar AJ
    (1997)
    A recombinationally repressed region between mat2 and mat3 loci shares homology to centromeric repeats and regulates directionality of mating-type switching in fission yeast
    Genetics 146:1221–1238.
    1. Guarente L
    (2000)
    Sir2 links chromatin silencing, metabolism, and aging
    Genes & Development 14:1021–1026.
    1. Nimmo ER
    2. Cranston G
    3. Allshire RC
    (1994)
    Telomere-associated chromosome breakage in fission yeast results in variegated expression of adjacent genes
    The EMBO Journal 13:3801–3811.
    1. Schultz J
    (1939)
    Proc Int Congr Genet
    The function of heterochromatin, Proc Int Congr Genet, Edinburgh, Cambridge University Press.
    1. Thon G
    2. Bjerling KP
    3. Nielsen IS
    (1999)
    Localization and properties of a silencing element near the mat3-M mating-type cassette of Schizosaccharomyces pombe
    Genetics 151:945–963.
    1. Thon G
    2. Bjerling P
    3. Bünner CM
    4. Verhein-Hansen J
    (2002)
    Expression-state boundaries in the mating-type region of fission yeast
    Genetics 161:611–622.
    1. Thon G
    2. Friis T
    (1997)
    Epigenetic inheritance of transcriptional silencing and switching competence in fission yeast
    Genetics 145:685–696.

Decision letter

  1. Edith Heard
    Reviewing Editor; Institut Curie, France
  2. Kevin Struhl
    Senior Editor; Harvard Medical School, United States

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

[Editors’ note: this article was originally rejected after discussions between the reviewers, but the authors were invited to resubmit after an appeal against the decision.]

Thank you for submitting your work entitled "A memory element imposes epigenetic behavior on intrinsically labile RNAi-induced heterochromatin spread" for consideration by eLife. Your article has been reviewed by a Senior Editor, a Reviewing Editor, and three reviewers.. The reviewers have opted to remain anonymous.

Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife. The overall conclusion is that this work is potentially of interest and could be worthy of publication in eLife in the future. However, in its current version it would require too much work and thus it cannot be deemed appropriate for revision within a reasonable time frame.

If you can incorporate the important critiques and suggestions made by the reviewers below, you are encouraged to resubmit your manuscript to eLife. In that case, the resubmitted manuscript would be considered a "new" submission and go through the same 2-stage process. However, unlike a true new submission, it would be typically seen by the editor and if needed the original reviewers.

Summary:

In this manuscript the authors set out to explore nature of heterochromatin spreading in fission yeast using a system that quantitatively measures heterochromatic silencing in living cells and follows the spreading reaction using fluorescent reporters. They report that spread of silencing from RNAi-nucleated elements is stochastic and multimodal and fluctuates over time. This is also supported by recent work by Obersriebnig et al. The authors also find another form of spreading, nucleated by the cis-acting element REIII, that they term deterministic, and that appears to show greater memory and fidelity in local spreading, being better buffered. This finding is novel. However, in the absence of a mechanism for REIII -mediated spreading, the work remains preliminary and there may be other interpretations of the data presented, as highlighted in the reviewers' comments. In summary, the authors do not really address the questions of heterochromatin spreading fidelity or memory as they set out to. In particular, evidence for actual "spreading" of heterochromatin from the so-called nucleation site and the genes assayed for silencing at some distance from it, is lacking. Based on the reviewers' comments, the overall conclusion is that this work is potentially of great interest and could be worthy of publication in eLife in the future. However, in its current version it would require too much work and thus it cannot be deemed appropriate for revision within a reasonable time frame.

Reviewer #1:

The fundamental problem conceptual with this paper is that the experiments are designed and interpreted in the context of a model in which spreading is assumed to be a process that emanates from a regulatory site and progressively reaches a detector site or sites. There is no evidence that "spreading" of heterochromatin of this type actually occurs in S. pombe or its relative S. cerevisiae. The crucial missing datum for all spreading models is some kind of time resolution between events at the cis-regulatory site, called the nucleation site in this study and others, and the genes at a distance that fall under the control of this regulatory site. Moreover, the way the experiments are conducted, for example by starting with cells and colonies that are preselected to be nucleated (green off) prevents the authors from the opportunity to observe the very data that could falsify the hypothesis (i.e. observing cells in which green is on and orange is off.)

The second major problem in the paper is the assertion of being able to distinguish between stochastic spreading and deterministic spreading. Let's for the moment imagine that spreading actually occurs. The distinction the authors make between stochastic and deterministic spreading is based upon the range of expression of the detector genes that are positioned to be under the influence of the nucleation site. If that range of expression is narrow, spreading is defined to be deterministic. If it is broad, spreading is defined to be stochastic. The competing model, never mentioned nor tested, is that heterochromatin formed in once case is more stable than heterochromatin formed in the other, but "spreading" is identical in both. Distinguishing between these two models would be easier than some of the experiments in the manuscript, but one needs to acknowledge the competing model before it can be tested.

The concept of fidelity of spreading is never defined and is less than obvious, yet experiments are interpreted in the context of this fidelity entity.

The "memory" experiments are interesting but again memory implies a time over which something is remembered. The experiments in Figure 4 may bear upon memory, but they are done in a rather complicated way show the sensitivity of those cells that have been pregrown for 25 generations to a subsequent challenge of various concentrations of TSA with respect to the stability of their expression state at that one time. A demonstration of memory requires a time component over which the thing of interest is remembered, but this experiments design precludes that. To be sure, there is some interesting biology behind the curves in this figure, but at present the experimental design would allow the authors to make interesting claims about the contribution of the various elements to the efficiency of re-establishing heterochromatin…. interesting in itself, but it isn't memory. It seems to me that the issue of memory would be resolved by the kind of clever experiments performed in S. cerevisiae by the Broach and Gartenberg labs in which the nucleation sites is removed and then the persistence of a state is measured over time in the absence of the element that "nucleated" it. Alternatively, the microfluidic technology behind Figure 3 could be used to address this issue.

Reviewer #2:

The assembly of heterochromatin domains involves nucleation and spreading of silencing factors. However, the degree of fidelity and heritable stability among heterochromatin domains throughout the genome can vary. Here, the authors dissect the origins of the intrinsic high-fidelity and memory capacity associated with stable epigenetic heterochromatin states. To do this, they utilize a system that simultaneously monitors nucleation and spreading driven by individual regulatory elements. They find that heterochromatin spreading driven by an RNAi nucleation site (cenH) is multimodal and occurs stochastically. This finding is consistent when measured at ectopic sites as well as at the endogenous MAT locus. By contrast, they argue that spreading driven solely from the endogenous REIII element at the MAT locus shows no intermediate levels and is complete ("deterministic"). They conclude that the REIII element is required to preserve the epigenetic memory of an ancestral state, as well as to ensure heterochromatin domain maintenance during environmental stress.

A major strength of this article is the use of a system that quantitatively measures heterochromatic silencing in living cells. The development of such tools has allowed them to tackle key questions in the field of heterochromatin assembly. The conclusion that the spreading of heterochromatin from an RNAi nucleated element is stochastic is also supported by recent work by Obersriebnig et al., 2016, which also used single cell assays. A key new finding here is that a cis-acting element, REIII, promotes deterministic spreading of heterochromatin that has high memory capacity. However, in the absence of a mechanism for the REIII -mediated spreading, this aspect of the work remains preliminary. Moreover, the evidence suggesting deterministic spreading from REIII needs further careful evaluation.

1) Regarding the mechanism for the role of REIII in epigenetic memory, previous studies have provided important insights that are highly relevant to this story and deserve to be addressed clearly as part of the main discussion. REIII (a binding site for ATF-CREB proteins Atf1-Pcr1) was initially shown to be critical for the stability of heterochromatin at MAT, particularly when heterochromatin nucleation by RNAi is compromised (Jia et al., 2004). A recent study confirmed this function of REIII (Wang and Moazed, 2017). However, more relevant to this paper are subsequent studies that revealed the mechanism as follows:

Yamada et al., 2005 found that REIII recruits an HDAC Clr3 (see Figure 2), the activity of which is critical for the stability and spread of heterochromatin. Further work showed that Clr3 is a component of the SHREC complex that localizes to REIII (see Figure 4 in Cam et al., 2008). These studies set the stage for a critical finding that explained the role of REIII /Clr3 in heterochromatin stability. Aygun et al., 2013 discovered Clr3 HDAC as a unique factor that promotes heterochromatin stability by suppressing turnover of histones. Indeed, loss of Clr3 is linked to loss of tri-H3K9me (see Figure 5A in Yamada et al., 2005), which provides the binding site for the chromodomain of Clr4 to reinforce a feedback loop in which tri-H3K9me recruits Clr4 HMTase to promote heterochromatin spreading/maintenance (Zhang et al. 2008; also, Al-Sady et al., and Jih et al., 2017).

In short: REIII → Atf1-Pcr1 → Clr3/SHREC → Suppresses histone turnover/stabilizes tri-H3K9me → Clr4 chromo binding to tri-H3K9me → Heterochromatin spreading/maintenance via feedback loop.

In light of these previous findings, an important possibility is that mutations in REIII impair the recruitment of Clr3 HDAC (e.g. see Figure 3 in Yamada et al., 2005), which in turn results in increased in turnover of H3K9 methylated histones that are required for the spreading and maintenance of heterochromatin. The authors shall test whether ∆REIII cells show an elevated level of histone turnover, hence affecting spreading and epigenetic memory.

2) An important assumption in this study is that cells carrying ORANGEON or ORANGEoff states contain comparable levels of heterochromatin at the nucleation site. However, it is possible that ORANGEON cells contain a relatively lower level of heterochromatin (but still sufficient to repress GREEN with a weak promoter) than ORANGEoff cells (e.g. in WT and ∆REIII cells). This should be tested experimentally by assaying H3K9me levels by ChIP at GREEN and ORANGE loci. Also, GREEN inserted in cenH is likely silenced post transcriptionally (transcripts of markers inserted near centromeric repeats are converted to siRNAs) and its silencing might not provide an accurate measure of silencing caused by heterochromatin modifications. If true, this would complicate the interpretation of the results.

3) This reviewer is concerned that the evidence that REIII causes deterministic spreading is an artifact of the system. Unlike WT or ∆REII cells that contain GREEN at a bona fide nucleation site, from which heterochromatin is known to spread outward, GREEN in ∆K cells is likely not placed at a nucleation site (there is no evidence that REIII alone can nucleate heterochromatin, and no evidence is presented to show that heterochromatin spreads from REIII toward ORANGE). In cells with GREEN located closer to the edge of the heterochromatin domain, the silencing of ORANGE located in the middle of the silenced domain might precede GREEN. In this situation, when silencing of ORANGE is assayed in GREENoff cells it would give the false impression of deterministic spreading. Consistent with this, markers inserted near mat3 are known to be silenced less efficiently. Similarly, if nucleation in ∆K cells occurred on the left side of ORANGE, all cells showing GREENoff will carry the ORANGEoff state. Can they confirm that heterochromatin in ∆K cells is indeed nucleated at REIII, or if REIII simply stabilizes heterochromatin (for example by recruiting Clr3 HDAC Yamada et al., 2005 and Aygun et al., 2013)? Moreover, it is also possible that heterochromatin in ∆K cells is formed in such as a way that it does not involve stepwise mechanism of nucleation and spreading and that the entire domain is silenced simultaneously.

4) The authors argue in the abstract that the epigenetic capacity contributing to heterochromatin spreading and maintenance is not encoded in the spreading reaction itself, but rather requires specialized memory elements. This statement is misleading. The factors involved in spreading and maintenance of heterochromatin, such as HDACs, are known to be loaded across silenced domains by HP1, which is also required for heterochromatin spreading. However, these same factors (HDACs) can be also loaded by so called specialized/memory elements, further enhancing the stability of heterochromatin. In this respect, the epigenetic capacity is at least partially encoded in the spreading reaction and can be further strengthened by elements that also load activities, such as HDACs, promoting heterochromatin stability.

Reviewer #3:

The fidelity and robustness are two essential characteristics to epigenetic inheritance, which is important both for chromatin structure and cell differentiation in many species. Its full establishment comprises of two stages: nucleation and spreading, which can be executed through RNA interference (RNAi)-dependent and RNAi-independent pathways.

In this study, Greenstein and colleagues set up an elegant and sensitive single-cell reporter system to dissect these two characteristics. They have used it for population studies and cell-lineage tracking for trans-generation memory. Consistent with previous findings, they showed spreading of heterochromatin nucleated by RNAi alone is robust but low on fidelity, as compared to that nucleated by wild-type locus. Importantly, they have attributed the fidelity in heterochromatin spreading to an RNAi-independent DNA element REIII. Further, they have tested the persistence of "memory" as nucleated by RNAi or REIII with and without perturbations. The biological importance and complicated nature should make it appealing to a broad audience.

In Figure 2B and 2C, where modified nucleators of RNAi-alone or REIII -alone are tested, the fidelity in spreading and robustness in nucleation segregated clearly in these population studies. However, one concern still remains. As spreading is greatly affected by the distance from the nucleator to the reporter, the actual distance for REIII -nucleated spreading should be measured between REIII and "orange" reporter, rather than between "green" and "orange" reporter as they are on different sides of the REIII nucleator, and to be consistent with Figure 2A and 2B. If the genome feature in Figure 2C is drawn to scale, it is around 1-2 kb in distance. At this distance, even RNAi-nucleated heterochromatin is demonstrating high fidelity in spreading as in the top panel of Figure 1C. At least one testing should be done with the reporter cassette moved further away from REIII nucleator by inserting "inert" DNA sequence (from a euchromatic region), to the similar 3-4 kb distance scale as Figure 2A and 2B. It is necessary to exclude the possibility that higher fidelity as in Figure 2C is due to a shorter loop for spreading. If the experiment is not possible, it shall be clearly explained in the text.

[Editors’ note: what now follows is the decision letter after the authors submitted for further consideration.]

Thank you for resubmitting your work entitled "A memory element imposes epigenetic behavior on intrinsically labile ncRNA-induced heterochromatin spread" for further consideration at eLife. Your revised article has been reviewed by Kevin Struhl (Senior Editor), a Reviewing Editor, and two reviewers.

The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below:

In the previous submission of this manuscript, the work presented was found to be potentially of great interest, particularly the discovery of a novel form of putative heterochromatin spreading, nucleated by the cis-acting element REIII that the authors claim has a deterministic role. However, direct evidence for REIII -mediated heterochromatin spreading was found to be lacking, and there was no mechanism. In their revised manuscript, the authors have provided several lines of new evidence to address some of the reviewers' comments. In particular, by deleting the REIII element, they show that this leads to increased turnover of histones which provides some insight into how REIII may be promoting epigenetic stability of heterochromatin. However, the conclusion that REIII promotes deterministic spreading of heterochromatin is still not really supported, as pointed out by reviewer 2. Furthermore, the argument that ∆K cells have better "memory" than WT cells is not supported by previous findings that ∆K cells contain metastable heterochromatin, with frequent conversion from the OFF to the ON state compared to WT. As both reviewers point out, it remains possible that bias may be introduced in the "memory" determinations, due to preselection for ∆KOFF colonies.

Despite these limitations, and given the new work included in this manuscript on this REIII -based heterochromatin process, we would be willing to consider a revised text, that focused on the single cell assay showing that loss of REIII affects the stability of heterochromatin. In this revised manuscript, specific references to « memory elements » and « deterministic » spreading should be removed, as suggested by reviewer 2.

Reviewer #2:

The authors performed several experiments that addressed the reviewers' comments and strengthened their conclusions. For example, the extensive chromatin immunoprecipitation analyses performed to assay H3K9me2 and H3K9me3 levels in different strains now correlates changes in gene silencing with heterochromatin levels. Also, this reviewer suggested that they investigate the effects of REIII deletion on histone turnover, which has implications for heterochromatin spreading and epigenetic memory. They do indeed find that deletion of REIII causes increased turnover of histones. Together with previously published findings showing that REIII recruits HDAC Clr3 (Yamada et al., 2005), which is involved in the suppression of histone turnover (Aygun et al., 2013), this finding provides a possible mechanistic insight into how REIII promotes epigenetic stability of heterochromatin. With the new additions, this part of the story has improved significantly.

However, their conclusion that REIII promotes deterministic spreading of heterochromatin remains ambiguous. In particular, unlike cenH that has been shown to nucleate heterochromatin, there is no clear evidence presented that heterochromatin spreads outward from the REIII element in ∆K cells. Rather, given the bistable nature of heterochromatin at mat in ∆K cells, the assembly of heterochromatin in ∆K cells cells might involve cooperativity between yet unknown elements, as opposed to reliance on the REIII element alone. In presenting their case for deterministic spreading from REIII, the authors mention previous work showing that loss of REIII affects heterochromatin in ∆K cells. However, it is important to note that previous studies (such as Jia, Noma and Grewal, 2004 and Wang and Moazed,.2017) and this article only show that REIII is required for the stability of heterochromatin (most likely via recruitment of HDAC Clr3 that suppresses histone turnover to prevent loss of H3K9me) but is not sufficient to nucleate heterochromatin. Similarly, their argument that ∆K cells have better "memory" than WT cells contrasts with a previous finding that ∆K cells contain metastable heterochromatin, as indicated by the frequent conversion from the OFF to the ON state at a rate many fold higher than in WT. Moreover, it remains possible that bias may be introduced in the "memory" determinations, due to preselection for ∆KOFF colonies (see subsection “Multi-generational single cell imaging reveals ncRNA-driven spreading to be unstable”), in which entire mat locus is known to be coated with heterochromatin (Nakayama, Klar and Grewal, 2000).

My recommendation is to prepare a revised manuscript for publication in eLife that focuses on the results from their elegant single cell assays showing that the loss of REIII affects the stability of heterochromatin and that this behavior of REIII is mechanistically linked to the suppression of histone turnover at mat. The discussion would be focused on succinctly describing this feature of REIII in the context of previous findings, including the results showing that REIII recruits the HDAC Clr3 (Yamada et al., 2005) that also suppresses histone turnover to promote epigenetic stability of heterochromatin (Aygun et al., 2013). Indeed, suppression of histone turnover by REIII /HDAC Clr3 may promote retention of methylated histones that are critical for the recruitment of Clr4; thus, promoting the efficient spreading of heterochromatin. Framing the results in this way would avoid the problematic issues raised by applying the terms "memory element" or "deterministic" spreading to REIII.

Reviewer #3:

The authors have addressed all the critiques raised by this reviewer. In particular, they have now provided more compelling evidence of "spreading" of heterochromatin from corresponding nucleation sites with a more complete testing cassette of different distances, directions, nucleation factors, etc. They have used these testing cassettes to demonstrate the stochasticity of ncRNA-mediated spreading and stability of REIII -mediated memory. In addition, they have now tested the anti-correlation of histone turnover and memory.

https://doi.org/10.7554/eLife.32948.029

Author response

[Editors’ note: the author responses to the first round of peer review follow.]

Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife. The overall conclusion is that this work is potentially of interest and could be worthy of publication in eLife in the future. However, in its current version it would require too much work and thus it cannot be deemed appropriate for revision within a reasonable time frame.

If you can incorporate the important critiques and suggestions made by the reviewers below, you are encouraged to resubmit your manuscript to eLife. In that case, the resubmitted manuscript would be considered a "new" submission and go through the same 2-stage process. However, unlike a true new submission, it would be typically seen by the editor and if needed the original reviewers. Summary: In this manuscript the authors set out to explore nature of heterochromatin spreading in fission yeast using a system that quantitatively measures heterochromatic silencing in living cells and follows the spreading reaction using fluorescent reporters. They report that spread of silencing from RNAi-nucleated elements is stochastic and multimodal and fluctuates over time. This is also supported by recent work by Obersriebnig et al. The authors also find another form of spreading, nucleated by the cis-acting element REIII, that they term deterministic, and that appears to show greater memory and fidelity in local spreading, being better buffered. This finding is novel. However, in the absence of a mechanism for REIII -mediated spreading, the work remains preliminary and there may be other interpretations of the data presented, as highlighted in the reviewers' comments. In summary, the authors do not really address the questions of heterochromatin spreading fidelity or memory as they set out to. In particular, evidence for actual "spreading" of heterochromatin from the so-called nucleation site and the genes assayed for silencing at some distance from it, is lacking. Based on the reviewers' comments, the overall conclusion is that this work is potentially of great interest and could be worthy of publication in eLife in the future. However, in its current version it would require too much work and thus it cannot be deemed appropriate for revision within a reasonable time frame.

We very much like to thank the editors and reviewers for their thoughtful feedback. We have made major additions to this work to address the reviewers’ comments. The experiments we have included here resulted in a more robust view of spreading in S. pombe and enable us to make mechanistic inferences about how the divergent spreading types we describeare enacted. Additionally, the data point to mechanisms underlying memory of repression that is required at a cell type specification locus like MAT. Together, they have made this a much stronger manuscript. The most important additions are as follows.

1) Establishing directionality of spreading: We have added an additional single cell tacking experiment in the ectopic strain, which together with our previous data, allows us to address the directionality of spreading using the FYLM device. We are happy to report that our data is consistent with a successive nucleation and spreading model. Additionally, we clarify our analysis method, which does not pre-filter cells of any nucleation proximal (“green”) or distal (“orange”) reporter values. We can now report that our cytometry without proximal or distal gating supports a nucleation seeded spreading model.

2) Tracking heterochromatin structure in our MAT locus strains: We have performed a suite of Chromatin Immuno-precipitation (ChIP) experiments to address the following issues: a. correlation of repression and heterochromatin structure, b. the stability of nucleation elements in strong and weak spreading strains and c. REIII -dependent heterochromatin formation. We can report that the new data support the conclusion that differences in spreading manifest as weaker heterochromatin but are not rooted in weak nucleation.

3) REIII mechanism and memory: We added an additional experiment, where we remove nucleation factors post-establishment, which more directly addresses the differences in memory over time in our strains. We provide a mechanistic explanation for increased or decreased memory by tracking histone 3 turnover at euchromatic and MAT locus heterochromatin sites. We can report that REIII confers a high degree of memory that is linked to very low histone turnover compared with strains nucleated by RNA-processes that show more euchromatic levels of histone turnover.

We have additionally adopted the reviewers’ suggestions to simplify and clarify the writing. Finally, both reviewer #1 and #2, compared our findings to those of Obersriebnig et al., a previous study of heterochromatin in single fission yeast cells. This was an important study and we would like to clarify the differences in experimental design and conclusions that distinguishes the two reports.

In the Obersriebnig work, heterochromatin formation was tracked either by following one color, at cenH, near REII or REIII, or with two colors, one at cenH and one near REIII in wild-type MAT. Since in all cases both REIII and cenH nucleators were intact, silencing of either color cannot be ascribed to spreading. In our work, we both inactivate one of the nucleators at a time and track nucleation behavior with one color and spreading with a second. The use of one color can yield silencing kinetics, but by design, the nucleation and spreading events cannot be untangled. The authors make claims about the stochasticity of heterochromatin establishment, which is novel, but not spreading, as is consistent with their experimental design. Conclusions they draw about establishment kinetics are noteworthy, yet those about spreading kinetics again are complicated by having both nucleators intact. This is important, because we believe the novelty of our work, beyond the novel behavior of REIII, is the ability to isolate any type of spreading event from the DNA/RNA sequence-driven nucleation component and define its characteristics, including in real time.

Reviewer #1:

The fundamental problem conceptual with this paper is that the experiments are designed and interpreted in the context of a model in which spreading is assumed to be a process that emanates from a regulatory site and progressively reaches a detector site or sites. There is no evidence that "spreading" of heterochromatin of this type actually occurs in S. pombe or its relative S. cerevisiae. The crucial missing datum for all spreading models is some kind of time resolution between events at the cis-regulatory site, called the nucleation site in this study and others, and the genes at a distance that fall under the control of this regulatory site.

We thank the reviewer for the comments which helped improve this work.

Spreading outward from nucleation sites, linearly or via loops, is the paradigm for propagating heterochromatin along the chromosome in S. pombe (Hall et al., 2002; Li et al., 2009; Noma et al., 2004; Zhang et al., 2008), S. cerevisiae (reviewed in (Kueng et al., 2013)), flies (reviewed in (Elgin and Reuter, 2013)), plants (Yang et al., 2017) and humans (Tchasovnikarova et al., 2015). However, we agree that experiments documenting kinetic separation of nucleation and spreading are rare. One example is vernalization of the FLC locus in Arabidopsis, where nucleation and spreading are separated in time (Yang et al., 2017). Our single cell tracking data already contained information that addresses this question directly, which we expanded in the revision:

Three lines of evidence establish that heterochromatin spreads outward from nucleation sites in our assays. First, we examined the kinetics of proximal and distal silencing in cells where the system switches from full derepression to repression. Our 3kb ectopic strain afforded us an opportunity to do so, (~5% of the cells in FYLM single cell tracking data of our 3kb ectopic strain, new Figure 3—figure supplement 1), while the MAT locus strains remain near uniformly repressed in the “green” nucleation proximal color. Full silencing of the “green” nucleation-proximal color almost exclusively precedes the distal “orange” spreading color in time (N=15/270 cells with repression events where “green” and “orange” start out fully de-repressed, see new Figure 3—figure supplement 1A,B). We find only one example (cat 4, 1/270) with “orange” repression preceding “green”, with the green closely following orange. This small kinetic gap therefore may suggest a fast spreading event.

Second, in the FYLM data, we did not observe any cells (N=0/270) where “green” becomes de-repressed while “orange” stays repressed (see new Figure 3—figure supplement 1A). In contrast we find many cells in which “orange” fluctuates while “green” is OFF for the entire measured lifespan (N=36/270, new Figure 3—figure supplement 1A,C) and more cells where “orange” is fully de-repressed while “green” is at least repressed partially for the entire measured lifespan (N=61/270, new Figure 3—figure supplement 1A,D). This reinforces the notion that silencing is directional with “green” preceding “orange”.

Third, using our steady state flow cytometry data for all RNA- nucleated strains, both ectopic and at MAT, gating for “green”OFF (nucleation site) yields a distribution of orange repression, OFF to ON, while gating in “orange”OFF (spreading site) only yields “green”OFF, supporting the above single-cell dynamics with a broader, population-based snapshot (see 2D-hexbins below as well as Figure 1C insets, Figure 2B,C,D). Therefore, we interpret this to indicate that spreading follows nucleation.

Moreover, the way the experiments are conducted, for example by starting with cells and colonies that are preselected to be nucleated (green off) prevents the authors from the opportunity to observe the very data that could falsify the hypothesis (i.e. observing cells in which green is on and orange is off.)

We apologize if our analysis method was not sufficiently clear, and we have now clarified the text. In fact, we do not preselect “green”OFF cells for analysis: All the 2D histogram and scatter plots in the paper show cells normalized only by “red”. No “green” or “orange” value filter is applied at this stage. Even for the ΔK strain, which can be isolated as ON or OFF alleles (Grewal and Klar, 1996), we describe the overall distribution of cells after introduction of clr4+ in Figure 2E. Since we established from our data (cytometry and single cell tracking) that nucleation indeed precedes spreading, we produce 1D “orange” histograms with an additional “green”OFF gate to isolate fully nucleated cells (Figure 1C, insets Figure 2B, C, E). Importantly, in the cytometry data, we do not observe cell populations that are “orange”OFF and “green”ON (Figure 1C inset plots are not pre-filtered and almost contain no cells in bottom right quadrant, highlighted for ectopic strains in Author response image 1). We believe our data is consistent with directional spreading from nucleation to distal sites and hope this clarification highlights that we include all cells with respect to “green” and “orange” expression for analysis.

Author response image 1
Exceedingly few cells are “green”ON and “orange”OFF.

2D hexbin plots of cells (max and min bin cell numbers indicated) with ectopic nucleation site (“green” nucleation proximal, “orange” distal at indicated distances, see Figure 1C). Fluorescence of all cells gated for cell size is plotted and normalized by “red” with no green or orange filtering. Note the very small number of cells detectable in the bottom right quadrant.

https://doi.org/10.7554/eLife.32948.027

The second major problem in the paper is the assertion of being able to distinguish between stochastic spreading and deterministic spreading. Let's for the moment imagine that spreading actually occurs. The distinction the authors make between stochastic and deterministic spreading is based upon the range of expression of the detector genes that are positioned to be under the influence of the nucleation site. If that range of expression is narrow, spreading is defined to be deterministic. If it is broad, spreading is defined to be stochastic. The competing model, never mentioned nor tested, is that heterochromatin formed in once case is more stable than heterochromatin formed in the other, but "spreading" is identical in both. Distinguishing between these two models would be easier than some of the experiments in the manuscript, but one needs to acknowledge the competing model before it can be tested.

The reviewer points out an alternative possible scenario for differences in spreading, which we now address in the text. It could be envisaged that the spreading probability depends primarily on the how stability of heterochromatin the nucleation site. In such a scenario the frequency of spreading would correlate precisely with the stability of the nucleation site. This model predicts that a narrow range of repression should therefore correlate with very robust heterochromatin, and conversely that spreading that has a wider range of states when the heterochromatic state is less stable or robust. Two lines of evidence argue against this possibility:

We examined the stability of the relevant nucleation sites to test this hypothesis, cenH in the wild-type (“orange” mostly repressed) or REIII -disrupted context (ΔREIII, stochastic spreading) and REIII in the cenH deleted context (ΔK, deterministic spreading). To do so, we performed Chromatin Immuno-precipitation.

(ChIP) experiments for the heterochromatin assembly mark H3K9me2 and the repressive/spreading mark H3K9me3 (see (Al-Sady et al., 2013; Jih et al., 2017)). Since factors assembling the heterochromatin structure are directly or indirectly attracted by these marks, we believe that they are the central features to examine. We observe that at all nucleation sites, the degree of both H3K9me2 and H3K9me3 methylation is extremely robust (around the value of the pericentromeric dh) and nearly invariant between the nucleation sites. Critically, there is no difference in H3K9me2 and me3 methylation at cenH between wildtype MAT and ΔREIII, strainswhich feature very different distributions of distal “orange” repression (Figure 2B, C). This new data is included in Figure 2D and E.

Author response image 2
“green” at cenH shows the same resistance to TSA irrespective of spreading behavior.

The same cells plotted in Figure 4B and C are plotted for the 0 TSA condition. Note that green lines represent fraction of total “green”OFF, orange lines represent fraction of “orange” repressed cells within “green”OFF.

https://doi.org/10.7554/eLife.32948.028

The TSA experiment in Figure 4 additionally afforded us an opportunity to compare directly the stability of the same cenH heterochromatin region in near-isogenic strains that spread robustly (wild-type MATHSS) or stochastically (∆REIIIHSS). To this end we plotted the fraction of “green”OFF cells along the [TSA] gradient (Author response image 2, relative to the 0 TSA point, =1) to determine if heterochromatin is differentially stable at the cenH locus between these strains. We found that the curves for “green” are essentially super-imposable, while the “orange” differs markedly, as shown in Figure 4B and C.

Together, we believe that both these new ChIP data and reanalysis of our existing TSA data shown here are not consistent with a model in which the differences in “orange” repression reflect stability of the heterochromatin domain.

Rather we believe that they more consistent with distinct heterochromatin spreading behaviors in these strains.

The concept of fidelity of spreading is never defined and is less than obvious, yet experiments are interpreted in the context of this fidelity entity.

We thank the reviewer for the helpful comments concerning terminology. We define spreading fidelity as the accurate inheritance of the positional and repressive extent of heterochromatin formed by spreading. High fidelity spreading should manifest as a tight distribution of states in the population and invariant inheritance when single cells are tracked. We have clarified and simplified terminology throughout the manuscript.

The "memory" experiments are interesting but again memory implies a time over which something is remembered. The experiments in Figure 4 may bear upon memory, but they are done in a rather complicated way show the sensitivity of those cells that have been pregrown for 25 generations to a subsequent challenge of various concentrations of TSA with respect to the stability of their expression state at that one time. A demonstration of memory requires a time component over which the thing of interest is remembered, but this experiments design precludes that. To be sure, there is some interesting biology behind the curves in this Figure, but at present the experimental design would allow the authors to make interesting claims about the contribution of the various elements to the efficiency of re-establishing heterochromatin…. interesting in itself, but it isn't memory. It seems to me that the issue of memory would be resolved by the kind of clever experiments performed in S. cerevisiae by the Broach and Gartenberg labs in which the nucleation sites is removed and then the persistence of a state is measured over time in the absence of the element that "nucleated" it. Alternatively, the microfluidic technology behind Figure 3 could be used to address this issue.

Measuring events over time is a central way to assess memory in a system, and we have provided additional measurements in this vein in the revised version. However, it is only one of several approaches to probe memory. The TSA experiments borrow conceptually from hysteresis measurements, which are a rigorous means to probe the memory capacity in the system, illustrated by Ferrell and colleagues (Bagowski and Ferrell, 2001). In this analysis, memory is manifested by the system behaving differently depending on its starting state. A stimulant-responsive system exhibits hysteresis (memory) if the behavior along the stimulant gradient is divergent between a fully stimulated or naïve starting condition. The similarity in molecular architectures, cooperativity and positive feedback, which drive both the biochemical systems traditionally subjected to hysteresis analysis and epigenetic systems (see (Angel et al., 2015; Dodd et al., 2007)), motivated us to perform hysteresis-type experiments. Here, we are probing whether the ancestral state biases the decision of system to undergo a spreading event. We measure the retention of this bias along a gradient of inhibitor, which is conceptually similar to measurements for a stimulant. We chose this course because stimulants are not readily available for heterochromatin. We believe this analysis is both rigorous and valid for probing memory capacity.

However, inspired by the reviewer’s invocation of the Gartenberg experiments, we probed the effect of removing trans-acting nucleation factors on the retention of the state originally formed by spreading. These experiments are conceptually similar to removal of cis-acting sites but avoid the caveat of low induced Cre-lox recombination frequency in S. pombe. We crossed our reporter strains to mutants required for either the cenH or REIII nucleator. Immediately following detection of cross progeny colonies, we proceeded with flow cytometry analysis. We find indeed that REIII nucleated heterochromatin can maintain itself after removal of a required trans-acting nucleation factor for ~ 200 generations, while cenH-nucleated heterochromatin cannot, and appears to require continual nucleation to maintain the state. These new data support our overall conclusions and are presented in the new Figure 5 and new Figure 5—figure supplement 1.

Finally, we also believe that our single cell tracking data in Figure 3 and Figure 3—figure supplement 2 lend independent support for the conclusion that REIII spreading is remembered, while cenH-spreading is not. Intergenerational tracking of spreading in the cenH, noncoding RNA (ncRNA)*-nucleatedstrains (Figure 3D) as well as ectopically ncRNA-nucleated strains (Figure 3—figure supplement 1 and Figure 3—figure supplement 2) show fluctuations of “orange” not compatible with strong memory, while REIII strains (Figure 3E) show no such fluctuations over similar generation times, strongly implying memory.

Together, our three independent approaches strongly support a differential in the ability to remember the state formed by spreading between cenH and REIII nucleated spreading.

We now use this more general descriptor for cenH and dh rather than RNAi.

Reviewer #2: The assembly of heterochromatin domains involves nucleation and spreading of silencing factors. However, the degree of fidelity and heritable stability among heterochromatin domains throughout the genome can vary. Here, the authors dissect the origins of the intrinsic high-fidelity and memory capacity associated with stable epigenetic heterochromatin states. To do this, they utilize a system that simultaneously monitors nucleation and spreading driven by individual regulatory elements. They find that heterochromatin spreading driven by an RNAi nucleation site (cenH) is multimodal and occurs stochastically. This finding is consistent when measured at ectopic sites as well as at the endogenous MAT locus. By contrast, they argue that spreading driven solely from the endogenous REIII element at the MAT locus shows no intermediate levels and is complete ("deterministic"). They conclude that the REIII element is required to preserve the epigenetic memory of an ancestral state, as well as to ensure heterochromatin domain maintenance during environmental stress. A major strength of this article is the use of a system that quantitatively measures heterochromatic silencing in living cells. The development of such tools has allowed them to tackle key questions in the field of heterochromatin assembly. The conclusion that the spreading of heterochromatin from an RNAi nucleated element is stochastic is also supported by recent work by Obersriebnig et al,. 2016, which also used single cell assays. A key new finding here is that a cis-acting element, REIII, promotes deterministic spreading of heterochromatin that has high memory capacity. However, in the absence of a mechanism for the REIII -mediated spreading, this aspect of the work remains preliminary. Moreover, the evidence suggesting deterministic spreading from REIII needs further careful evaluation.

We thank the reviewer for this assessment.

1) Regarding the mechanism for the role of REIII in epigenetic memory, previous studies have provided important insights that are highly relevant to this story and deserve to be addressed clearly as part of the main discussion. REIII (a binding site for ATF-CREB proteins Atf1-Pcr1) was initially shown to be critical for the stability of heterochromatin at MAT, particularly when heterochromatin nucleation by RNAi is compromised (Jia et al., 2004). A recent study confirmed this function of REIII (Wang and Moazed, 2017). However, more relevant to this paper are subsequent studies that revealed the mechanism as follows: Yamada et al., 2005 found that REIII recruits an HDAC Clr3 (see Figure 2), the activity of which is critical for the stability and spread of heterochromatin. Further work showed that Clr3 is a component of the SHREC complex that localizes to REIII (see Figure 4 in Cam et al., 2008). These studies set the stage for a critical finding that explained the role of REIII /Clr3 in heterochromatin stability. Aygun et al., 2013 discovered Clr3 HDAC as a unique factor that promotes heterochromatin stability by suppressing turnover of histones. Indeed, loss of Clr3 is linked to loss of tri-H3K9me (see Figure 5A in Yamada et al., 2005), which provides the binding site for the chromodomain of Clr4 to reinforce a feedback loop in which tri-H3K9me recruits Clr4 HMTase to promote heterochromatin spreading/maintenance (Zhang et al. 2008; also, Al-Sady et al., and Jih et al., 2017). In short: REIII → Atf1-Pcr1 → Clr3/SHREC → Suppresses histone turnover/stabilizes tri-H3K9me → Clr4 chromo binding to tri-H3K9me → Heterochromatin spreading/maintenance via feedback loop In light of these previous findings, an important possibility is that mutations in REIII impair the recruitment of Clr3 HDAC (e.g. see Figure 3 in Yamada et al..,2005), which in turn results in increased in turnover of H3K9 methylated histones that are required for the spreading and maintenance of heterochromatin. The authors shall test whether ∆REIII cells show an elevated level of histone turnover, hence affecting spreading and epigenetic memory.

We thank the reviewer for raising these points. We had included a description of possible mechanisms in our supplemental discussion in the original submission, including a number of the references cited by the reviewer. However, we recognize that the mechanistic arguments should be highlighted more prominently in the paper and have included them now in the main discussion.

We have now addressed the role of histone turnover as a possible mechanism for the differences in spreading and epigenetic memory by performing additional experiments. We adopted the published (Audergon et al., 2015; Svensson et al., 2015) Recombination Induced Tag Exchange system (RITE) to probe H3 turnover. Given the labile nature of heterochromatin at elevated temperatures (Figure 4E,F and Figure 4—figure supplement 2), replication stalling was performed with hydroxyurea as published (Aygun et al., 2013). We address turnover in the ∆KHSS-OFF, ∆KHSS-ON and ∆REIIIHSS strains to address the role of REIII, with ∆KHSS-ON serving as an isogenic non-heterochromatic control for ∆KHSS-OFF. We measured H3 turnover at two euchromatic control loci and find no difference between the strains. However, when we examined turnover at MAT using amplicons that are shared by all strains, we found significantly decreased histone turnover in ∆KHSS-OFF, while ∆KHSS-ON and ∆REIIIHSS consistently display elevated turnover, more consistent with euchromatic behavior. These new data are included in the new Figure 6 A,B. The data that histone turnover is elevated in the absence of nucleated REIII thus presents a mechanistic explanation for the strongly reduced memory we see in REIII and stochastic spreading, while the exceedingly low turnover in the presence of REIII accounts for the deterministic spreading and robust memory. The robust memory is underlined by a new experiment we include, where we remove by cross factors required for nucleation following establishment (new Figure 5).

2) An important assumption in this study is that cells carrying ORANGEON or ORANGEoff states contain comparable levels of heterochromatin at the nucleation site. However, it is possible that ORANGEON cells contain a relatively lower level of heterochromatin (but still sufficient to repress GREEN with a weak promoter) than ORANGEoff cells (e.g. in WT and ∆REIII cells). This should be tested experimentally by assaying H3K9me levels by ChIP at GREEN and ORANGE loci. Also, GREEN inserted in cenH is likely silenced post transcriptionally (transcripts of markers inserted near centromeric repeats are converted to siRNAs) and its silencing might not provide an accurate measure of silencing caused by heterochromatin modifications. If true, this would complicate the interpretation of the results.

We thank the reviewer for the comment. We have now performed extensive Chromatin Immuno-precipitation (ChIP) analysis for the heterochromatin assembly mark H3K9me2 and repressive/spreading mark H3K9me3 across the MAT locus.

Specifically, we find:

1) H3K9me2 and H3K9me3 marks accumulate to similar extents at the cenH nucleation site (“green” in WT and ∆REIIIHSS), even though this strain, unlike WT, features a high proportion of “orange”ON (new Figure 2D). We also find very similar proportion of H3K9me2 and me3 at “green” in ∆KHSS-OFF, see the new data in Figure 2E. H3K9me2 and me3 marks accumulate to very high extents at these sites in the three strains, which is not consistent with heterochromatin spreading behaviors being solely the result of different level of heterochromatin assembly at nucleation regions.

2) H3K9me2 and me3 marks decline rightwards of cenH only in the ∆REIIIHSS, but not wild type stain. This is evident even at the edge of cenH, and more strongly at the more distal “orange” (Figure 2D). However, H3K9me2 and me3 do not significantly decline at “orange” in ∆KHSS-OFF versus “green” (Figure 2E). We note that for H3K9me3 only there is a small, but not significant, difference between “green” and “orange”.

Together, these data are consistent with differences in heterochromatin spreading, as opposed to weakened nucleation, driving the differences in distal “orange” repression, and further that our single cell silencing data for “orange” accurately reflects differences in the spreading of heterochromatin outwards from nucleation sites.

In S. pombe, heterochromatin assembled via noncoding RNA (ncRNA)* processes triggers co-transcriptional silencing as the reviewer states. However, transcriptional silencing also plays important roles (Fischer et al., 2009) at these sites, and silencing broadly correlates well with H3K9me2 accumulation, a read-out of heterochromatin structure. Yet, it remains possible that placement within an RNAi generating element affects the behaviors we track. We believe the reviewer’s concern ought to be mitigated by a difference in the placement of “green” at ectopic loci versus the MAT locus. In ectopic strains, “green” (for ura4, “orange” for his1) is placed immediately adjacent but not within the dh element, reporting on silencing at the nucleation locus, while for MAT strains “green” is within cenH. As we find the behavior with respect to spreading to be similar (e.g. Figure 1C 3kb versus Figure 2C) we conclude that the nature of H3K9me- dependent silencing at the nucleation marker not to interfere with our interpretation.

*We now use this more general descriptor for cenH and dh rather than RNAi

3) This reviewer is concerned that the evidence that REIII causes deterministic spreading is an artifact of the system. Unlike WT or ∆REII cells that contain GREEN at a bona fide nucleation site, from which heterochromatin is known to spread outward, GREEN in ∆K cells is likely not placed at a nucleation site (there is no evidence that REIII alone can nucleate heterochromatin, and no evidence is presented to show that heterochromatin spreads from REIII toward ORANGE). In cells with GREEN located closer to the edge of the heterochromatin domain, the silencing of ORANGE located in the middle of the silenced domain might precede GREEN. In this situation, when silencing of ORANGE is assayed in GREENoff cells it would give the false impression of deterministic spreading. Consistent with this, markers inserted near mat3 are known to be silenced less efficiently. Similarly, if nucleation in ∆K cells occurred on the left side of ORANGE, all cells showing GREENoff will carry the ORANGEoff state. Can they confirm that heterochromatin in ∆K cells is indeed nucleated at REIII, or if REIII simply stabilizes heterochromatin (for example by recruiting Clr3 HDAC Yamada et al., 2005 and Aygun et al., 2013)? Moreover, it is also possible that heterochromatin in ∆K cells is formed in such as a way that it does not involve stepwise mechanism of nucleation and spreading and that the entire domain is silenced simultaneously.

We thank the reviewer for this comment.

Recent evidence suggests that in the absence of cenH, Atf1/Pcr1 binding sites at REIII are absolutely required for heterochromatin formation (Figure 3C in (Wang and Moazed, 2017)). Independently of these studies, three lines of evidence from our work suggest that ΔKHSS reports specifically on REIII -nucleated spreading, detected by “green”, towards “orange”:

1) Deletion of both Atf1/Pcr1 binding sites at REIII blocks any silencing in ΔKHSS (in 34/34 strains examined) and H3K9 methylation at “green” or “orange”. This new data is presented in Figure 2—figure supplement 2A,B. Because ΔREIII does not abolish heterochromatin at MAT, i.e. it is not required to stabilize heterochromatin formed elsewhere (Figure 2C), the above result points to REIII nucleating heterochromatin in ΔKHSS.

2) No repressive elements other than REII, cenH and REIII have been identified at MAT. Technically, it possible that in ΔKHSS, which lacks cenH, silencing is due to REII. We had addressed this in the original submission by probing ΔKHSS REII::LEU2, in which the REII element is deleted. This strain behaved similarly to ΔKHSS (Figure 2—figure supplement 1C), strongly suggesting that in ΔKHSS, REIII is the only functional nucleator.

3) We should note that in our 2D histogram and scatter plots (Figure 1, Figure 2 and Figure 5) cells are not gated for “green”OFF. Rather, all cells within the size gate are plotted normalized by “red”. In those 2D plots, no significant de-repression of “green” is evident in ΔKHSS in the repressed population. We only use the “green”OFF gate for 1D histograms, to plot spreading in cells that fully nucleate. To address the point concerning spreading from leftwards of the “orange” nucleator, we moved “green” 2.1 kb towards the right IR-R boundary in ΔKHSS. This did not lead to the appearance of significant de-repression of “green” in “orange”OFF cells, which might be expected if spreading emanated from a site leftwards of “orange” in ΔKHSS. This new data is presented Figure 2—figure supplement 2 C,D.

Together, we believe these lines of evidence and previously published data affirms that REIII is indeed the nucleating element in ΔKHSS and that it spreads heterochromatin outward.

As to the comment about a non-stepwise mechanism, this is entirely possible for ΔKHSS and we had discussed this possibility in the previous supplemental, and now moved it to the main discussion. Heterochromatin spreading can occur by a number of theoretical mechanisms (Talbert and Henikoff, 2006), such as linear “oozing”, skipping or looping. Looping mechanisms have been invoked for the polycomb system (Bantignies and Cavalli, 2011). Looping is one possible mechanism to account for deterministic behavior of ΔK, while linear “oozing” is more consistent with our observations for ncRNA-nucleated spreading (Figure 1, Figure 2C). In experiments outside the scope of this work, we look forward to testing this hypothesis.

4) The authors argue in the Abstract that the epigenetic capacity contributing to heterochromatin spreading and maintenance is not encoded in the spreading reaction itself, but rather requires specialized memory elements. This statement is misleading. The factors involved in spreading and maintenance of heterochromatin, such as HDACs, are known to be loaded across silenced domains by HP1, which is also required for heterochromatin spreading. However, these same factors (HDACs) can be also loaded by so called specialized/memory elements, further enhancing the stability of heterochromatin. In this respect, the epigenetic capacity is at least partially encoded in the spreading reaction and can be further strengthened by elements that also load activities, such as HDACs, promoting heterochromatin stability.

We thank the reviewer for raising this point. We have modified the language in the Abstract to reflect that the ncRNA-nucleated spreading reaction specifically does not confer the epigenetic fidelity required for robust and precise epigenetic inheritance, which additionally requires memory elements.

Reviewer #3: The fidelity and robustness are two essential characteristics to epigenetic inheritance, which is important both for chromatin structure and cell differentiation in many species. Its full establishment comprises of two stages: nucleation and spreading, which can be executed through RNA interference (RNAi)-dependent and RNAi-independent pathways. In this study, Greenstein and colleagues set up an elegant and sensitive single-cell reporter system to dissect these two characteristics. They have used it for population studies and cell-lineage tracking for trans-generation memory. Consistent with previous findings, they showed spreading of heterochromatin nucleated by RNAi alone is robust but low on fidelity, as compared to that nucleated by wild-type locus. Importantly, they have attributed the fidelity in heterochromatin spreading to an RNAi-independent DNA element REIII. Further, they have tested the persistence of "memory" as nucleated by RNAi or REIII with and without perturbations. The biological importance and complicated nature should make it appealing to a broad audience.

We thank the reviewer for the favorable assessment of this work.

In Figure 2B and 2C, where modified nucleators of RNAi-alone or REIII -alone are tested, the fidelity in spreading and robustness in nucleation segregated clearly in these population studies. However, one concern still remains. As spreading is greatly affected by the distance from the nucleator to the reporter, the actual distance for REIII -nucleated spreading should be measured between REIII and "orange" reporter, rather than between "green" and "orange" reporter as they are on different sides of the REIII nucleator, and to be consistent with Figure 2A and 2B. If the genome feature in Figure 2C is drawn to scale, it is around 1-2 kb in distance. At this distance, even RNAi-nucleated heterochromatin is demonstrating high fidelity in spreading as in the top panel of Figure 1C. At least one testing should be done with the reporter cassette moved further away from REIII nucleator by inserting "inert" DNA sequence (from a euchromatic region), to the similar 3-4 kb distance scale as Figure 2A and 2B. It is necessary to exclude the possibility that higher fidelity as in Figure 2C is due to a shorter loop for spreading. If the experiment is not possible, it shall be clearly explained in the text.

While REIII contains two potential nucleation sites (s1 and s2), only the distal (s2) sites has been implicated as the dominant nucleation element (Thon et al., 1999). The distance in ΔKHSS from the end of the spreading sensor ORF to s1 is ~1kb, and to s2 ~2.2kb. In comparison, the comparable distance in ΔREIII HSS is ~3.8kb. To mitigate the concern that the spreading distance in ΔK may be much shorter than in ΔREIII (i.e., from the ORF end to s1), we deleted s1 in the ΔK strain (ΔKΔs1) and inserted two different ~700bp sections of the sib1 gene ORF between the spreading sensor and the s2 site in the ΔKΔs1 background, increasing the spreading distance to ~3kb. Importantly, in both of these strains, we see the same deterministic behavior as in the ΔK parent. One of the ΔKΔs1 +700bp insertions, representative of two behaviorally indistinguishable insertion strains we constructed, is shown in the new Figure 2—figure supplement 2E. Larger insertions where not compatible with heterochromatin nucleation in ΔK. Given these data, and the observation that significant stochastic spreading and intermediate states are detected for the 3kb, and even 1kb spreading distance with ectopically placed ncRNA nucleators (Figure 1), we conclude that the deterministic behavior is ΔK in intrinsic to nucleation from the REIII element.

[Editors’ note: the author responses to the re-review follow.]

The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below: In the previous submission of this manuscript, the work presented was found to be potentially of great interest, particularly the discovery of a novel form of putative heterochromatin spreading, nucleated by the cis-acting element REIII that the authors claim has a deterministic role. However, direct evidence for REIII -mediated heterochromatin spreading was found to be lacking, and there was no mechanism. In their revised manuscript, the authors have provided several lines of new evidence to address some of the reviewers' comments. In particular, by deleting the REIII element, they show that this leads to increased turnover of histones which provides some insight into how REIII may be promoting epigenetic stability of heterochromatin. However, the conclusion that REIII promotes deterministic spreading of heterochromatin is still not really supported, as pointed out by reviewer 2. Furthermore, the argument that ∆K cells have better "memory" than WT cells is not supported by previous findings that ∆K cells contain metastable heterochromatin, with frequent conversion from the OFF to the ON state compared to WT. As both reviewers point out, it remains possible that bias may be introduced in the "memory" determinations, due to preselection for ∆KOFF colonies. Despite these limitations, and given the new work included in this manuscript on this REIII -based heterochromatin process, we would be willing to consider a revised text, that focused on the single cell assay showing that loss of REIII affects the stability of heterochromatin. In this revised manuscript, specific references to « memory elements » and « deterministic » spreading should be removed, as suggested by reviewer 2. Reviewer #2: The authors performed several experiments that addressed the reviewers' comments and strengthened their conclusions. For example, the extensive chromatin immunoprecipitation analyses performed to assay H3K9me2 and H3K9me3 levels in different strains now correlates changes in gene silencing with heterochromatin levels. Also, this reviewer suggested that they investigate the effects of REIII deletion on histone turnover, which has implications for heterochromatin spreading and epigenetic memory. They do indeed find that deletion of REIII causes increased turnover of histones. Together with previously published findings showing that REIII recruits HDAC Clr3 (Yamada et al., 2005), which is involved in the suppression of histone turnover (Aygun et al., 2013), this finding provides a possible mechanistic insight into how REIII promotes epigenetic stability of heterochromatin. With the new additions, this part of the story has improved significantly.

We thank the reviewer for this positive assessment.

However, their conclusion that REIII promotes deterministic spreading of heterochromatin remains ambiguous. In particular, unlike cenH that has been shown to nucleate heterochromatin, there is no clear evidence presented that heterochromatin spreads outward from the REIII element in ∆K cells. Rather, given the bistable nature of heterochromatin at mat in ∆K cells, the assembly of heterochromatin in ∆K cells cells might involve cooperativity between yet unknown elements, as opposed to reliance on the REIII element alone. In presenting their case for deterministic spreading from REIII, the authors mention previous work showing that loss of REIII affects heterochromatin in ∆K cells. However, it is important to note that previous studies (such as Jia, Noma and Grewal, 2004 and Wang and Moazed,.2017) and this article only show that REIII is required for the stability of heterochromatin (most likely via recruitment of HDAC Clr3 that suppresses histone turnover to prevent loss of H3K9me) but is not sufficient to nucleate heterochromatin.

The reviewer raises an important point. In ∆K cells, formation of H3K9me- marked heterochromatin is dependent on the action of REIII Atf1 and Pcr1 binding sites. Since there is no pre-existing H3K9me before introduction of clr4+ in our experiments, heterochromatin nucleation must occur in some fashion in ∆K cells. However, sufficiency of Atf1 and Pcr1 binding sites for nucleation remains to be demonstrated. Specifically, REIII, unlike cenH (e.g. Hall el al., 2002) cannot effect significant repression when placed ectopically in the genome (references cited and Figure 2—figure supplement 2). This result can be interpreted in at least two ways: REIII is not by itself sufficient to nucleate and requires other, yet to be defined elements, or, the local genomic context/or chromatin organization of REIII within MAT is required for its nucleation activity, hence it cannot function ectopically. Since it cannot be easily resolved which case is correct, we now soften the conclusion that REIII drives spreading in ∆K cells in the text and remove the term “deterministic spreading” in favor of REIII -dependent heterochromatin.

Similarly, their argument that ∆K cells have better "memory" than WT cells contrasts with a previous finding that ∆K cells contain metastable heterochromatin, as indicated by the frequent conversion from the OFF to the ON state at a rate many fold higher than in WT. Moreover, it remains possible that bias may be introduced in the "memory" determinations, due to preselection for ∆KOFF colonies (see subsection “Multi-generational single cell imaging reveals ncRNA-driven spreading to be unstable”), in which entire mat locus is known to be coated with heterochromatin (Nakayama, Klar and Grewal, 2000).

We thank the reviewer for this point. Here, we (1) Clarify our observations in the context of the literature, (2) Revise our interpretation of Figure 4D, and (3) Clarify our ∆K selection process.

Broadly, our results from ∆K cells are consistent with the published literature, namely a very low transition rate form OFF to ON per generation, about 10-4, (for example Thon and Friis, 1997; Grewal and Klar, 1998). The continuous growth in non-selective media of OFF cells over > 25 generations resulted in the appearance 1-2% ON colonies, consistent with these low transitions (Grewal and Klar, 1996). A key difference between our and some previous experiments is observing the entire population without selection (our assays) versus selection assays that permit survival of either the ON or OFF state. Thus, the appearance of very small ON populations may be less obvious in our experiments.

The memory experiments in Figure 4 reveal the key contribution of REIII to resistance and memory of MAT heterochromatin. The comparison of wild-type to ∆REIII cells here is central and valid in our view. However, we recognize that in Figure 4D, the results on its face point to stronger memory in ∆K cells than wild type. The “extreme” aspect of this memory is most likely due to the fact that ∆K cells, once ancestrally TSA treated, cannot efficiently re-nucleate. Given this caveat, we now refrain from interpreting Figure 4D as indicative of the degree of memory in ∆K. Overall, we did not intend to advance the notion that ∆K cells have more/stronger memory than wild-type and wherever this is ambiguous, have clarified the language. We note that unlike for ∆K cells, it is not possible to select for WT cells with active REIII, thus the population may contain some cells with REIII in an inactive state, which likely are more labile than the bulk of cells with active REIII. Hence, ∆KOFF populations can appear somewhat more stable as a result.

As to ∆KOFF selection, our approach is consistent with the literature (Grewal and Klar, 1996 and 1998; Jia et al., 2004; Thon and Friis, 1997 and others). We produce ∆KHSS strains in a ∆clr4 background and introduce clr4+ by cross. The population resulting from the germinated spore is often mixed, with both ON and OFF cells as shown in Figure 2E. Simply plating this population on nonselective media predominantly yields colonies that are ON or OFF. This reflects the metastable nature of ∆K cells, with some cells nucleating heterochromatin and some not, but then maintaining the state. In that sense, we do not isolate ∆KOFF cells using a directed selection process, promoting a given state, and do not feel that there should be any bias introduced by this method. Isolation of ON/OFF is required in some cases, for example, due to the insufficient number of OFF cells in the mixed population. However, we now qualify our results pertaining to memory clearly and state that bias may have been introduced by picking OFF colonies.

In general, as suggested, we re-focused our conclusions on ∆REIII. Together with our reinterpretation of the ∆K result in Figure 4, we believe we addressed remaining concerns.

My recommendation is to prepare a revised manuscript for publication in eLife that focuses on the results from their elegant single cell assays showing that the loss of REIII affects the stability of heterochromatin and that this behavior of REIII is mechanistically linked to the suppression of histone turnover at mat. The discussion would be focused on succinctly describing this feature of REIII in the context of previous findings, including the results showing that REIII recruits the HDAC Clr3 (Yamada et al., 2005) that also suppresses histone turnover to promote epigenetic stability of heterochromatin (Aygun et al., 2013). Indeed, suppression of histone turnover by REIII /HDAC Clr3 may promote retention of methylated histones that are critical for the recruitment of Clr4; thus, promoting the efficient spreading of heterochromatin. Framing the results in this way would avoid the problematic issues raised by applying the terms "memory element" or "deterministic" spreading to REIII.

We thank the reviewer for these suggestions. Indeed, we have re-written the manuscript, most prominently the Abstract and Discussion section, to focus on the overall result that ncRNA spreading is unstable and stabilized by REIII via suppression of histone turnover. We no longer use the terms “deterministic” or “memory element” with respect to REIII.

Reviewer #3: The authors have addressed all the critiques raised by this reviewer. In particular, they have now provided more compelling evidence of "spreading" of heterochromatin from corresponding nucleation sites with a more complete testing cassette of different distances, directions, nucleation factors, etc. They have used these testing cassettes to demonstrate the stochasticity of ncRNA-mediated spreading and stability of REIII -mediated memory. In addition, they have now tested the anti-correlation of histone turnover and memory.

We thank the reviewer for this positive assessment.

https://doi.org/10.7554/eLife.32948.030

Article and author information

Author details

  1. R A Greenstein

    1. Department of Microbiology & Immunology, George Williams Hooper Foundation, University of California San Francisco, San Francisco, United States
    2. TETRAD graduate program, University of California San Francisco, San Francisco, United States
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology, Writing—review and editing
    Competing interests
    No competing interests declared
  2. Stephen K Jones

    Department of Molecular Biosciences, The University of Texas at Austin, Austin, United States
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology, Writing—review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8387-1917
  3. Eric C Spivey

    Department of Molecular Biosciences, The University of Texas at Austin, Austin, United States
    Present address
    Department of Biomedical Engineering, Vanderbilt University, Nashville, United States
    Contribution
    Resources, Software, Formal analysis, Validation, Investigation, Methodology, Writing—review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4080-8616
  4. James R Rybarski

    Department of Molecular Biosciences, The University of Texas at Austin, Austin, United States
    Contribution
    Resources, Data curation, Software, Formal analysis, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9960-1161
  5. Ilya J Finkelstein

    1. Department of Molecular Biosciences, The University of Texas at Austin, Austin, United States
    2. Center for Systems and Synthetic Biology, The University of Texas at Austin, Austin, United States
    Contribution
    Conceptualization, Data curation, Supervision, Funding acquisition, Visualization, Methodology, Writing—review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9371-2431
  6. Bassem Al-Sady

    Department of Microbiology & Immunology, George Williams Hooper Foundation, University of California San Francisco, San Francisco, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing—original draft, Project administration
    For correspondence
    bassem.al-sady@ucsf.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8996-7941

Funding

National Institutes of Health (DP2GM123484)

  • Bassem Al-Sady

University of California, San Francisco (Program for Breakthrough Biomedical Research, New Frontier Research)

  • Bassem Al-Sady

American Federation for Aging Research (AFAR-020)

  • Ilya J Finkelstein

National Institutes of Health (F32 AG053051)

  • Stephen K Jones

Welch Foundation (F-l808)

  • Ilya J Finkelstein

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Shiv I Grewal, Karl Ekwall, and Hiten D Madhani for their generous gifts of fission yeast strains. We thank Graham A Anderson and Shengya Cao for stimulating discussions, especially on hysteresis, and Brandan La for the initial Matlab scripts for cytometry data analysis. In addition, we thank Carol A Gross for substantial help with writing the manuscript and Jonathan S Weissman and Sigurd Braun for critical comments. This work was supported by grants from the National Institutes of Health (DP2GM123484) and the UCSF Program for Breakthrough Biomedical Research (partially funded by the Sandler Foundation) to BA-S, American Federation of Aging Research (AFAR-020) and the Welch Foundation (F-l808) to IJF and the National Institute of Aging (F32 AG053051) to SKJ. Flow Cytometry data was generated in the UCSF Parnassus Flow Cytometry Core which is supported by the Diabetes Research Center (DRC) grant, NIH P30 DK063720.

Senior Editor

  1. Kevin Struhl, Harvard Medical School, United States

Reviewing Editor

  1. Edith Heard, Institut Curie, France

Version history

  1. Received: October 19, 2017
  2. Accepted: July 17, 2018
  3. Accepted Manuscript published: July 18, 2018 (version 1)
  4. Version of Record published: August 1, 2018 (version 2)

Copyright

© 2018, Greenstein et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,775
    Page views
  • 516
    Downloads
  • 20
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. R A Greenstein
  2. Stephen K Jones
  3. Eric C Spivey
  4. James R Rybarski
  5. Ilya J Finkelstein
  6. Bassem Al-Sady
(2018)
Noncoding RNA-nucleated heterochromatin spreading is intrinsically labile and requires accessory elements for epigenetic stability
eLife 7:e32948.
https://doi.org/10.7554/eLife.32948

Further reading

    1. Cell Biology
    2. Chromosomes and Gene Expression
    Maikel Castellano-Pozo, Georgios Sioutas ... Enrique Martinez-Perez
    Short Report Updated

    The cohesin complex plays essential roles in chromosome segregation, 3D genome organisation, and DNA damage repair through its ability to modify DNA topology. In higher eukaryotes, meiotic chromosome function, and therefore fertility, requires cohesin complexes containing meiosis-specific kleisin subunits: REC8 and RAD21L in mammals and REC-8 and COH-3/4 in Caenorhabditis elegans. How these complexes perform the multiple functions of cohesin during meiosis and whether this involves different modes of DNA binding or dynamic association with chromosomes is poorly understood. Combining time-resolved methods of protein removal with live imaging and exploiting the temporospatial organisation of the C. elegans germline, we show that REC-8 complexes provide sister chromatid cohesion (SCC) and DNA repair, while COH-3/4 complexes control higher-order chromosome structure. High-abundance COH-3/4 complexes associate dynamically with individual chromatids in a manner dependent on cohesin loading (SCC-2) and removal (WAPL-1) factors. In contrast, low-abundance REC-8 complexes associate stably with chromosomes, tethering sister chromatids from S-phase until the meiotic divisions. Our results reveal that kleisin identity determines the function of meiotic cohesin by controlling the mode and regulation of cohesin–DNA association, and are consistent with a model in which SCC and DNA looping are performed by variant cohesin complexes that coexist on chromosomes.

    1. Chromosomes and Gene Expression
    2. Developmental Biology
    Airat Ibragimov, Xin Yang Bing ... Paul Schedl
    Research Article Updated

    Though long non-coding RNAs (lncRNAs) represent a substantial fraction of the Pol II transcripts in multicellular animals, only a few have known functions. Here we report that the blocking activity of the Bithorax complex (BX-C) Fub-1 boundary is segmentally regulated by its own lncRNA. The Fub-1 boundary is located between the Ultrabithorax (Ubx) gene and the bxd/pbx regulatory domain, which is responsible for regulating Ubx expression in parasegment PS6/segment A1. Fub-1 consists of two hypersensitive sites, HS1 and HS2. HS1 is an insulator while HS2 functions primarily as an lncRNA promoter. To activate Ubx expression in PS6/A1, enhancers in the bxd/pbx domain must be able to bypass Fub-1 blocking activity. We show that the expression of the Fub-1 lncRNAs in PS6/A1 from the HS2 promoter inactivates Fub-1 insulating activity. Inactivation is due to read-through as the HS2 promoter must be directed toward HS1 to disrupt blocking.