Monitoring of switches in heterochromatin-induced silencing shows incomplete establishment and developmental instabilities

Farah Bughioa, Gary R. Huckellb,1, and Keith A. Maggerta,2

aDepartment of Cellular and Molecular Medicine, University of Arizona, Tucson, AZ 85724; and bPrivate address, San Diego, CA 92120

Edited by James A. Birchler, University of Missouri, Columbia, MO, and approved August 28, 2019 (received for review June 6, 2019) Position effect variegation (PEV) in Drosophila results from new While the natural factors that influence establishment have juxtapositions of euchromatic and heterochromatic chromosomal only recently begun to be explored (34, 35), multiple models for regions, and manifests as striking bimodal patterns of gene ex- maintenance have been proposed, including DNA cytosine pression. The semirandom patterns of PEV, reflecting clonal rela- and histone modifications stably inherited through tionships between cells, have been interpreted as gene-expression S-phase (36–39). The former does not exist in the common yeast states that are set in development and thereafter maintained model systems, in Caenorhabditis elegans or in Drosophila (40, 41; without change through subsequent cell divisions. The rate of in- cf. refs. 42 and 43), and so it is widely accepted that epigenetic stability of PEV is almost entirely unexplored beyond the final gene-regulatory (e.g., heterochromatin-induced silencing) infor- expression of the modified gene; thus the origin of the expressiv- mation must be encoded by the histone modifications in these ity and patterns of PEV remain unexplained. Many properties of model systems. Histone modifications (and DNA methylation) PEV are not predicted from currently accepted biochemical and are part of the mechanisms of activation or repression, although theoretical models. In this work we investigate the time at which they do not necessarily encode information instructive of gene expressivity of silencing is set, and find that it is determined be- expression. Histone modifications and DNA methylation pat- fore heterochromatin exists. We employ a mathematical simula- terns are easily changed by introduction of transcriptional re- GENETICS tion and a corroborating experimental approach to monitor pressors or activators (44–48). Intentional modification or switching (i.e., gains and losses of silencing) through develop- recruitment of heterochromatin affect ment. In contrast to current views, we find that gene silencing is incompletely set early in embryogenesis, but nevertheless is repeat- (49, 50), but it is not clear that modification of histones alone is necessary and sufficient to delimit self-sustaining (i.e., “epige- edly lost and gained in individual cells throughout development. ” Our data support an alternative to locus-specific “epigenetic” silenc- netic ) heterochromatin. Furthermore, many experimental ob- “ ing at variegating gene promoters that more fully accounts for the servations are at odds with the simple model of establish and ” final patterns of PEV. maintain epigenetics, suggesting a need to refine our under- standing of these phenomena. heterochromatin | position effect variegation | Drosophila | gene For example, in Drosophila (51) and Saccharomyces cerevisiae silencing | epigenetics (52), excision of a silenced reporter gene from a heterochromatic locus completely derepresses silencing, even in interphase, nalyses of position effect variegation (PEV) in Drosophila showing that the silencing information is not maintained on Amelanogaster and yeasts are the foundation of modern thegeneinquestion.Rather,silencing must be continually models of epigenetic inheritance. In Drosophila, the normally imposed by the juxtaposing heterochromatin. Furthermore, + euchromatic cell-autonomous white gene necessary for eye pig- mentation is subject to heterochromatin-induced gene silencing Significance through transposition or chromosome rearrangement (Fig. 1A) + (1–4). PEV of white manifests as clonal patches of contiguous The genetic properties of heterochromatin can create strikingly ommatidia that share an expression state (ON or OFF) (Fig. 1 B random patterns of inherited gene-expression states in other- and C), resulting in hallmark “clonal expression variation.” This wise identical cells. The origin of these patterns of clonal ex- developmental clonality of expression has led to concepts of es- pression variation are not understood: How and when in tablishment and maintenance. At some point in development an development it is determined how many cells will express, how asymmetry in otherwise identical sister cells was randomly estab- expression states are inherited through DNA replication, and lished, leading to expression in 1 lineage and gene silencing in the how often instabilities arise. We developed genetic, mathe- other (5). Thenceforth, expression states were maintained through matical, and molecular–genetic approaches to investigate these subsequent cell divisions, including through 1 or many S-phases questions. Our results indicate that these characteristics are when all chromosome-bound proteins are removed from chro- predetermined before chromatin structure is set, and the in- mosomes during replication. This behavior of PEV demanded that heritance displays remarkable instability through development. the biochemical processes of DNA replication be appended to account for locus-specific regulatory “epigenetic” information Author contributions: F.B. and K.A.M. designed research; F.B. and K.A.M. performed re- being maintained through S-phase in parallel to the replication of search; F.B., G.R.H., and K.A.M. contributed new reagents/analytic tools; F.B. and K.A.M. DNA sequence. These concepts of epigenetic establishment and analyzed data; and F.B. and K.A.M. wrote the paper. maintenance have been definitively demonstrated for some The authors declare no conflict of interest. phenomena—X inactivation (6–8), genomic imprinting (9), chromo- This article is a PNAS Direct Submission. some imprinting (10, 11), centromere identity (12–14)—but the Published under the PNAS license. concept now broadly pervades theories of gene expression in 1Retired. general (15–21), in ecology (15, 22–25), in evolution (26–28), and 2To whom correspondence may be addressed. Email: [email protected]. in disease development (29–31), which imagine many or all gene- This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. expression states to be induced and inherited in much the same 1073/pnas.1909724116/-/DCSupplemental. way as heterochromatin-induced silencing (cf. refs. 20, 32, and 33). First published September 16, 2019.

www.pnas.org/cgi/doi/10.1073/pnas.1909724116 PNAS | October 1, 2019 | vol. 116 | no. 40 | 20043–20053 Downloaded by guest on September 28, 2021 A w+ 2 - transposition w+ 1 - chromosome inversion

BCD

E

Fig. 1. PEV expresses a range of phenotypes. (A) Chromosome inversions or transpositions create new heterochromatin–euchromatin breakpoints, and genes + + near these new breakpoints undergo PEV with different clonal expression patterns. Red circle (“w ”) represents the white gene, thick stippled line rep- resents heterochromatin, thin line represents euchromatin, dumbbell shape is transposed DNA, cross-hatching indicates variegation. (B) Patched (or sectored) variegation of Y,10C.(C) Salt-and-pepper variegation of Y,B840.1 (63). (D) Low expressivity of heterochromatin-induced gene silencing resulting in high + expression of the variegating white gene in In(1)wm4. Left and right eyes of the same individual are shown. (E) High expressivity of silencing, and thus low + white expression, in an isogenic sibling of the fly in D. (Magnification, B and C,36×; D and E,12×.)

derepression is correlated with high-frequency natural changes showed both patterns expressed a remarkably similar appearance to ribosomal DNA (rDNA) copy number, which is known to throughout development (60). affect heterochromatin function (53, 54), suggesting unmapped We sought an alternative, heterochromatin-focused, explana- genome changes may underlie some of the PEV pattern. The tion that could account for the same clonal patterns of PEV. This manifest expression versus nonexpression in PEV may therefore alternative would need to explain 4 additional features of PEV not be due to established and maintained silencing information that current models do not. First, how can expressing and non- at gene promoters, but rather fluctuations in overall hetero- expressing patches be demarcated before the variegated gene chromatin “strength” in the nucleus. A straightforward way to promoter is activated (58, 61)? Second, why are the temperature- test this possibility would be to determine whether multiple sensitive periods of PEV limited to very early and very late in variegating genes fluctuate in concert when independently jux- development when heterochromatin first forms and later un- taposed to different heterochromatic regions. Indeed, when dergoes developmentally induced relaxation (58, 60)? Third, how tested at the chromatin of the silent mating-type loci in S. cer- is maintenance of silencing independent of the cell cycle (60), evisiae, 2 stochastically inactivated marker genes were perfectly despite S-phase being the salient critical time for maintenance? correlated in their expression (55). Fourth, how can different heterochromatic regions of the ge- In Drosophila, centric or Y-linked heterochromatin often in- nome induce different expressivities of PEV, strongly repressed duce large patches (or sectors) of expression (Fig. 1B), whereas to near wild-type (62, 63), and how can they induce different chromosome 4 or telomeric heterochromatin often result in clonal expression patterns (2, 60)? “salt-and-pepper” variegation (Fig. 1C), although exceptions In this work, we address the fundamental question of when occur in both directions. The reasons for these different mani- expressivity (i.e., the extent of silencing as a fraction of pig- festations are not understood. It was initially imagined that early mented and nonpigmented ommatidia) of silencing is deter- stochastic establishment of gene silencing by heterochromatin mined. We find that it is generally set before, at, or around the would create large patches, while late establishment would cre- time of fertilization (likely during oogenesis), which is inconsis- ate salt-and-pepper patterns (2, 56–58). The fact that variegating tent with the currently accepted model of silencing being set + white alleles that possess the natural enhancer/promoter of the early then lost from individual genes through development. This + white gene [e.g., In(1)wm4, In(1)wm4h] are not expressed until led us to investigate the time of derepression. We created a late pupae (59) would seem prima facie to invalidate these mathematical simulation of clonal variation to test the accepted models. However, work by Joel Eissenberg and colleagues (58, model and an alternative model: That the final patterns of PEV 60) refuted the hypothesis of differential establishment of si- can be better understood as a uniformly early event that affects lencing, instead supporting an alternative view where silencing whether a clone of descendent cells manifests relatively stable or was uniformly established early and stochastically lost at differ- unstable silencing (rare or frequent switches both ON-to-OFF ent times in development. These 2 models differ very little, as and OFF-to-ON). We recapitulated both patched and salt-and- they both envision clonal patterns as similar errors in establish- pepper patterns by simulating this alternative model. The sim- ment or maintenance that differ only in their developmental ulation also indicated a strongly discriminating test: The ac- timing. This timing framework still stands, despite observations cepted model does not predict gains of silencing through by Janice Spofford that the different patterns may differ in ge- development, while the alternative relies on such occurrences. netic properties (2), and by Eissenberg and colleagues, who Gains of silencing have not been directly observed in the past,

20044 | www.pnas.org/cgi/doi/10.1073/pnas.1909724116 Bughio et al. Downloaded by guest on September 28, 2021 mostly for methodological reasons. We used a variegating gal80 A BC reporter gene in conjunction with a lineage marking system (64) to monitor and observe switches directly in living precursor cells of the variegating eye. Taken together, our data suggest that different expressivities and clonal patterns of PEV result from differences in the early establishment of silencing, and in failures of heterochromatin to continually reestablish silencing continu- ally through development. Results and Discussion DEF The extent of silencing in the adult eye (i.e., the expressivity of PEV) is a combination of at least 2 things: Variation in the amount of silencing in cell lineages within an individual, and variation in the amount of silencing between isogenic siblings (e.g., the individuals in Fig. 1 C and D are isogenic siblings) (63). The time at which this latter contribution to expressivity is established is not known. G Bilaterally Correlated PEV Shows That Expressivity Is Set Around Fertilization. We reasoned that if the prevailing model for PEV is correct (65–70), then every S-phase represents an independent potential loss of silencing, andsolossesofsilencinginany lineage would be independent of losses in any other lineage. We tested this prediction by analyzing the variegation in both halves of individuals from recently isogenized strains of 3 variegating + − − white reporter genes: whitemottled 4, whitemottled 4h, and Y,10C.

The first 2 are chromosomes bearing inversions that juxtapose GENETICS the natural white enhancer/promoter with centric heterochro- Fig. 2. Clonal variation in left and right eyes is correlated. (A) Categories + matin of the X chromosome. Y,10C is a Y-linked P{white }, (≤25% of the eye pigmented, <50% [but more than 25%], <75%, <100%, + fully wild-type) of pigmentation of left (abscissa) and right (ordinate) eyes which possesses a white transgene controlled by the minimal m4 hs from individual In(1)w male flies. Areas of the red circles represent frac- heat-shock-70 gene promoter (w ) (71, 72) and a GFP gene tion of the total population (n = 180 eyes) falling into those categories. expressed by a promoter (73), both of which variegate Areas of the black circles represent expectation if left and right eyes were (54). We quantified expression categorically (as in ref. 63) and + independent, determined by binomial probabilities of population distribu- found a strong correlation in expressivity of white between left tion derived from population pigmentation data. The observed and expected 2 and right eyes of the same individuals from all 3 strains (Fig. 2 A–C; data were compared using χ analysis. Rejected H0 states that the eyes are see figure legend for statistical analyses). This correlation was independent in their pigment categories, χ2 = 79.21, df = 14. (B)AsinA but clear even when the flies were raised at high (29 °C) or low with In(1)wm4h, n = 168, χ2 = 51.25. (C)AsinA but with Y,10C, n = 194, χ2 = m4 2 (12 °C) temperatures, which suppress and enhance variegation, 48.84. (D)AsinA but with In(1)w raised at 29 °C, n = 152, χ = 80.21. (E)Asin m4 = χ2 = respectively, or were raised on medium containing a chemical A but with In(1)w raised at 12 °C, n 126, 226.95. (F)AsinA but with In(1)wm4 raised in the presence of 50 mM sodium butyrate, n = 102, χ2 = 30.81. suppressor of variegation (74) (Fig. 2 D–F). Rare extreme excep- + (G) Distributions of pigment categories from F2 offspring of flies with high tions to the preponderant range of white expression were addi- degree of silencing (pink, n = 152), medium silencing (red, n = 190), or low tionally informative. Eyes with almost no expression (strongly silencing (dark red, n = 210). P < 0.05 for graphs in A–F. repressed, in some cases merely single ommatidia expressing + white )werealways(n = 10 of 10) found on a fly head opposite a second such rare eye. We confirmed that these exceptions were or placodal determination, before cytological condensation or not the consequence of cryptic polymorphisms by collecting low- heterochromatin formation in nuclear division cycle 14 (78–82), + expressing exceptions, flies with medial ranges of expression, and before the activation of the hsp70 (58) or white promoters (59, high-expressing flies, and mating them pairwise. Even in the sec- 73), before the onset of any zygotic transcription at the mid- ond subsequent generation, these 3 subpopulations each re- blastula transition (83), and before heterochromatin-induced capitulated normal and statistically indistinguishable variation in silencing is first genetically detectable (58, 84). This suggests a expression (Fig. 2G)(H0 states that these 3 populations of pig- common factor that affects expressivity between individuals is set ment categories derive from 1 population, χ2 = 7.15, df = 8, P > down in eggs, sperm, or at fertilization (85), is influential prior to 0.5). We conclude that the variation in expression is not due to and at the time heterochromatin is first formed, and its influence + spontaneous second-site enhancers or suppressors, either domi- persists through development until as late as white activation nant or recessive. in pupae. This correlation between left and right indicates that a com- To determine if this correlation persists throughout develop- mon factor controls expressivity in both halves of the organism. ment, we analyzed third-instar larval males bearing Y,10C Syncytial nuclei, and tissues separated by morphogenetic move- chromosomes and scored for expression of variegating GFP in- ments after cellularization, generally do not cross the left–right dependently in the left and right eye imaginal discs. Silencing + division plane; thus, whatever factor correlates PEV expressivity of the white and GFP transgenes on a normally paternally in adults acts prior to the separation of the presumptive eye inherited Y,10C was near-complete, so we removed the paternal anlagen in embryos. Analysis of patterns of X chromosome loss imprint on the Y (63) by crossing chromosomes through females in gynandromorphs (75–77) indicates that the nuclear division of and analyzing the male progeny of a C(1)DX/Y,10C mother. the embryo that separates left and right imaginal cuticular an- When total fluorescence in dissected imaginal discs was quanti- lagen is usually the very first zygotic division. Thus, we interpret fied, the levels of expression of GFP varied between isogenic the correlations between the left and right eyes of adults to re- individuals but were not statistically different between left and veal a common silencing potential between the first 2 zygotic right halves of the organism (Fig. 3A) (regression coefficient 1.00, nuclei. This occurs well before cellularization, before germ-layer 95% confidence interval 0.71 to 1.30, Pearson’s R2 = 0.734;

Bughio et al. PNAS | October 1, 2019 | vol. 116 | no. 40 | 20045 Downloaded by guest on September 28, 2021 AB

CD

+ Fig. 3. Relationship of clonal variation of 2 linked variegating genes: Ubiquitin-GFP and white .(A) Correlation between GFP fluorescence in left- and right- wing imaginal discs dissected from third-instar males bearing Y,10C. Fluorescence units are arbitrary, but consistent across all individuals (n = 21) of the + experiment. (B) Correlation between different white pigment categories (abscissa) and GFP expression in the thorax and abdomen (ordinate), taken from adult males bearing Y,10C selected on the basis of the former. (C) Illustrative examples for B.(D) Illustrative examples of individual Y,10C-bearing male third- instar larvae, selected on the basis of low GFP expression (Upper) or high GFP expression (Lower) in nongerm tissue, then scored for white+ expression as adults (Right). (Magnification, eyes in C and D,25×; adults in C,12×; and larvae in D,15×.)

rejected H0 states that there is no correlation between fluorescence GFP expression, 17 (26%) were white or nearly so, while 48 + of the left and right discs, specifically that the aforementioned (74%) had high levels of white pigmentation as adults. We + coefficient is not unity, F = 52.44, P < 0.05). These results in- conclude that the adult expression of white is statistically pre- dicate that the influence of an early-acting factor can be detected dicted by larval GFP expression (rejected H0 states that a priori throughout development. selection based on GFP fluorescence would have no bearing on + white pigmentation as adults, χ2 = 48.0305, df = 1, P < 0.05). Covariegation of Linked Genes Confirms the Persistence of This effect is unlikely to be cross-talk between enhancers and Expressivity through Development. Our finding that the propensity promoters of the 2 reporter genes because it has been previously to derepress is set before repression exists challenges current shown that induction of 1 variegating gene does not affect var- concepts for heterochromatin-induced silencing that posit that iegation of a closely linked variegating gene (60, 86). epigenetic information is stored at a gene promoter, and requires These 3 correlations—GFP expression between different tis- + us to consider PEV as a property of heterochromatin strength sues of larvae (Fig. 3A), white expression to GFP expression in (i.e., its ability to act as a transcriptional repressor). This has a adults (Fig. 3 B and C), and larval GFP expression to adult + strong prediction not made by the current model, namely that white expression (Fig. 3D)—all argue that the variation seen in expressivity in 1 part of the body should correlate with expressivity PEV is not reflecting variation in epigenetic information at gene in other parts of the body. To test this prediction, we again ana- loci, but rather is reflecting organism-by-organism variation in lyzed males bearing matroclinous Y,10C chromosomes. We sorted + some common factor, set at the time of fertilization, and gen- adult males with strong and weak silencing of white and then erally maintained in the soma throughout development. quantified the level of GFP expression in thoraces and abdomens of whole animals. Expression of the former was strongly correlated Development of an Alternative Model for PEV. The currently ac- with the latter (Fig. 3 B and C) (Kendall’s τ = 0.695, rejected H0 cepted model for PEV posits that large patches of expression states that there is no correlation between pigmentation category result from early derepression events while smaller clones of salt- and GFP fluorescence, P < 0.05). and-pepper variegation result from late events of derepression. It Finally, we confirmed that expressivity was stable through is implicit in this model that, regardless of when derepression development by hand-sorting male third-instar larvae by cate- occurs, silencing once lost is permanently lost. Although many gorical (“high” and “low”) GFP fluorescence, and returning mathematical modelings of epigenetic silencing have hypotheti- them to food and allowing them to grow to adulthood (Fig. 3D). cal gains and losses with each mitosis (50, 87–90), we are not Of 69 larvae with low GFP expression, 59 (86%) had white or aware of any study that has explained any mechanism for gains of near-white eyes as adults, the remainder (14%) possessed eyes silencing. Thus, the assumption that silencing once lost is never + with many white -expressing ommatidia. Of 65 larvae with high gained has remained almost entirely unchallenged.

20046 | www.pnas.org/cgi/doi/10.1073/pnas.1909724116 Bughio et al. Downloaded by guest on September 28, 2021 We envisioned another possibility for the origin of these 2 indicates that Gal80 has never been absent because even a short distinct patterns of clonal variation that can accommodate gains period of time of Gal4 activity, in the absence of Gal80, is suf- in silencing. We reasoned that it would be possible to create ficient to cause the FLP-dependent irreversible activation of clonal patches or individual salt-and-pepper patterns by leaving GFP expression. fixed the time of “establishment” of heterochromatin-induced The presence of these nonfluorescent patches of cells indicates silencing (i.e., early embryo), and instead alter the stability (the that those cell clones expressing a variegating gene late in de- probability of accurate reestablishment) of silencing in sub- velopment are the result of progenitor cells in which silencing of + sequent divisions. This alternative hypothesis might explain the gal80 was never established. For example, white clones in adult 2 extreme patterns of PEV. Large patches could arise from early eyes may best be understood not as cells in which silencing was establishment/nonestablishment and relatively low frequencies of established and lost, but as cells in which establishment was gains and losses in silencing (decisions once made are faithfully defective. This observation is consistent with the presence of maintained through development). Salt-and-pepper variegation a very early “sensitive” period of heterochromatin-mediated si- could arise from equally early establishment/nonestablishment lencing and our finding that factors controlling expressivity may followed by relatively high frequencies in switching (many sub- act around fertilization. sequent losses and gains of silencing through development). To Green and green-yellow fluorescent cell clones indicate losses of silencing, establish the validity of this hypothesis, we first created a sim- and yellow clones indicate persistent silencing. Cells expressing only plified simulation to determine if stability rather than timing GFP have had Gal4 activity at some point in development, and could create patched and salt-and-pepper patterns of PEV. The then lost it due to derepression of gal80. Gal80 binds and re- −9 simulation and our analyses can be found in SI Appendix. A key presses preformed Gal4-DNA complex with a Kd < 10 , a 10- finding was that the alternative stability hypothesis stated above fold higher affinity than that of Gal4 for DNA (95, 96), so the requires the occurrence of gains of silencing throughout devel- gal4-dependent G-TRACE is very sensitive to gal80 expression opment, while the conventional timing hypothesis does not. regardless of the level of gal4 expression. In solely green fluo- Therefore, we sought to create a live reporter system sensitive to rescent clones (Fig. 5C), derepression must have occurred some such gain events. hours or days before, to provide Gal80 time to repress Gal4- mediated transcription, and Gal4-dependent RFP mRNA and Switch Monitoring to Characterize PEV through Development. In time to decay below detection limits. Some clones had

order to monitor switches of heterochromatin-induced silencing low (Fig. 5D) or moderate (Fig. 5E) levels of RFP at the time of GENETICS throughout development, we considered heterochromatin func- analysis, presumably from more-recent losses of silencing and tion as a developmental stage, and adopted the G-TRACE de- incomplete decay of RFP gene products (e.g., a cell shortly after velopmental lineage tracer (64) to monitor it. G-TRACE contains point “3” in Fig. 4C). For the most part, low or moderate levels a UAS-RFP that reports on current gal4 activity, and a paired of RFP occurred in patches of contiguous cells with comparable UAS-FLP and ubiquitin-FRT-STOP-FRT-GFP that reports on RFP fluorescence (Fig. 5A), which we interpret to be loss events historic gal4 activity. To create a system particularly sensitive to in the single progenitor cell to that clone. The sizes of such acquisitions of silencing, we made gal4 expression ubiquitous and clones modally varied from about 15 to 60 cells, indicating losses embedded the Gal4 repressor, gal80, in silencing heterochromatin of silencing 4 to 6 cell divisions prior to analysis. Cell clones that (Fig. 4A). had both high-RFP and high-GFP (Fig. 5F) likely indicate per- In this configuration a red fluorescent cell in a nonfluorescent sistent silencing or losses of silencing very recently. clone would indicate recent gains of silencing of gal80, the point We observed rare contiguous clusters of yellow fluorescent of discrimination between the current and the alternative hy- cells that contained individual green fluorescent cells in their potheses. Yellow fluorescence would indicate persistent silenc- midst (Fig. 5G). We believe these individuals to be cells that lost ing, green fluorescence the loss of silencing, and nonfluorescent silencing individually without passage through the cell cycle. We cells those in which silencing was not established (Fig. 4B). Be- acknowledge that those individual cells are sisters to cells that cause RFP expression precedes GFP expression (Fig. 4C, events might also have lost silencing, perhaps as a product of their ulti- 1 and 2) and because RFP decays once gal4 is repressed (Fig. 4C, mate shared cell division; however, we see no evidence of paired between events 3 and 4), whereas GFP will not, the RFP-to-GFP sister cells with similarly reduced RFP. The half-life of RFP is ratio and size of clones can further be used to infer the timing of reported to be 90 h (97), although it is not known if this is true of gains and losses of silencing. Cells in which red fluorescence is all cell types and in all developmental stages, so loss of silencing matched with submaximal green fluorescence are those in which may have happened much earlier or more recently in develop- silencing has been gained [and Gal80 decayed, with a half-life of ment. Regardless, we believe that green fluorescent cells within ∼2 h (91)], and the degree of green fluorescence will indicate yellow fluorescent patches represent bona fide loss of silencing how long ago that happened (Fig. 4D). Cells in which green events. fluorescence is maximal and red fluorescence is not are those in In terminally differentiated examples of clonal expression + which silencing was lost, and the degree of RFP decay will in- variation, such as white ommatidia in the Drosophila eye or + dicate how long ago derepression occurred. We used this “switch yellow cuticular spots, clones of expressing cells would corre- monitoring” (SwiM) system to detect losses and gains of silenc- spond to a sum of those cells that never established silencing ing through development of living eye-antennal imaginal discs. (i.e., nonfluorescent cells in the SwiM system) and those in which silencing was established and then lost (i.e., green fluorescence Analysis of Eye-Antennal Imaginal Discs Using the SwiM System. or those cells in which green fluorescence is maximal and red Because of the high degree of variation in PEV between individ- fluorescence is dwindling). Application of the SwiM system has uals (Figs. 1 D and E and 2) (92–94) it was impossible for us to allowed us to discriminate between the different histories of meaningfully determine frequencies of these events in populations heterochromatin-induced silencing that these cells experienced of flies. We therefore sought to test the stability hypothesis by through development. Based on our studies, we conclude that close analysis of individual imaginal discs (Fig. 5A). clonal variation is best viewed as originating from 2 different Nonfluorescent cell clones indicate incomplete establishment of silencing. sources: Failure to establish silencing and failure to maintain it. In larval eye-antennal and wing discs we observed clones of cells We view PEV as a “heterochromatin” phenomenon rather than that were nonfluorescent (Fig. 5B). These correspond to cells in a “promoter” phenomenon; therefore, it is not certain that these which the ubiquitous gal4 activity was inhibited by the presence are mechanistically different events (see below in Conclu- of Gal80. Notably, the lack of green fluorescence in these cells sions and Further Thoughts).Ashasbeenreportedearlier

Bughio et al. PNAS | October 1, 2019 | vol. 116 | no. 40 | 20047 Downloaded by guest on September 28, 2021 A

B

C

D

Fig. 4. Schematic diagram of the SwiM system. (A) Graphical representation of the components of the SwiM system. Cross-hatched GAL80 (purple) rep- resents variegating gal80 gene. Red uRFP (UAS-RFP), blue uFLP (UAS-FLP), aGFP (actin-GFP), and circle with triangle (excised FRT-Terminator) extrachromo- somal circle represent components of the G-TRACE lineage tracer system. (B) Fluorescence outcomes of different repression and derepression states for the SwiM system. Asterisks (*) indicate the time of either the loss-of-silencing or gain-of-silencing events. (C) Delays between RFP activation (1) and GFP activation (2), and between RFP inactivation and the beginning of detectable mRNA and protein decay (3), and the final state without detectable RFP (4), allow for comparison of GFP and RFP fluorescence in individual cells to determine whether a cell has undergone gal80 repression (1, 2) or derepression (3, 4). (D) Quick- reference table summarizing (B and C).

(51, 60, 98), we see no evidence that cell division represents a We detected red fluorescent cells appearing in 2 patterns. particularly sensitive point in maintaining heterochromatin- First, they appeared as individual cells within fields of non- induced silencing. fluorescent cells (Fig. 5H). These cells have lost gal80 expression Red cell clones indicate gains of silencing. The stability hypothesis for the first time in development. Such red fluorescent cells will predicted that loss events would be observed throughout devel- mature into yellow fluorescence as GFP accumulates, and close opment, but also that we would detect gains of silencing. The analysis of red and green fluorescence in these cells indicates SwiM system was designed to be sensitive to gains, which are that there is indeed a very low level of green fluorescence in expected to appear as red fluorescent cells with no or low levels of many of them at the time of analysis. GFP expression. The time frame in which RFP accumulates to The presence of these cells supports the hypothesis that the detectability prior to the build-up of FLP, the FLP-catalyzed ge- final expressivity and pattern of clonal variation includes gains as nome rearrangement, and the build-up and maturation of GFP to well as losses of heterochromatin-induced gene silencing. It is detectability is relatively narrow, less than 2 h. We therefore therefore preferable to think of PEV as fluctuations in hetero- expected red fluorescent cells to be rare, but critically informative. chromatin strength through development, perhaps similar to the

20048 | www.pnas.org/cgi/doi/10.1073/pnas.1909724116 Bughio et al. Downloaded by guest on September 28, 2021 A Green We also observed a few GFP-expressing and non-RFP cells that were multiple cell diameters away from any other fluorescent cell (Fig. 5I). These are likely to be cells in which silencing was first lost, then reacquired, all without passage through S-phase. We interpret these cells as direct indicators of cell-autonomous fluc- tuations of heterochromatin strength. The second pattern we detected was of a few red fluorescent cells occurring in clusters (Fig. 5J). We can imagine 2 possible explanations for this observation. First, they could represent artifacts of the gain of gal80 silencing at some point in the past Red (e.g., 2 cell divisions for clones of 4 cells), and the recent decay of Gal80 to a level incapable of repressing Gal4. Second, they could represent synchronized gains of silencing in sisters, cousins, and second-cousins. This scenario would indicate that cells are pre-

D F determined to gain (or lose) silencing at some time in the future. J

C I Supportive evidence that this is the correct interpretation comes

E G from the preponderance of small green fluorescent clusters in B nonfluorescing clones, suggesting that losses of silencing (which

H are not subject to issues of perdurance of Gal80)arealsopre- determined some divisions before manifestation. The nature of such a predeterminant is unknown. However, these predetermined B CD instabilities may be explained as a consequence of fluctuations in the efficacy of a heterochromatin component in cell clones. For example, if the concentration of some protein necessary for het- erochromatin function were to decrease in a single cell, one might expect that the entire clone descended from it may synchronously

EFG lose silencing. GENETICS A different lineage tracer more sensitive to changes. To confirm the dynamism of switching we infer from the G-TRACE components, we combined the gal4 and variegating gal80 transgene components of the SwiM system with I-TRACE (101). This lineage tracer uses a similar UAS-RFP, but also contains a ubiquitously expressed HI J GFP and a UAS-driven short-hairpin RNA directed at GFP. With I-TRACE in place of G-TRACE, red fluorescent cells would be those in which gal80 is currently repressed, green cells indicative of current gal80 expression and consequent repression of gal4,and yellow cells revealing new repression or derepression events Fig. 5. SwiM system analysis of a third-instar eye-antennal imaginal disk. as GFP and RFP fluorescence replace each other (Fig. 6A). (A) Entire disk with red- and green-channel separations. Separations are I-TRACE is therefore more sensitive to changes in expression of presented as inverted black-on-white for clarity, bordered with green and gal80, but unable to resolve the direction of the switch. In imaginal red. In all cases the green separation is on the top and the red on the bot- discs from flies bearing ubiquitously expressed gal4, variegating tom. (Inset) Location of subsequent regions (B–J). (B) Nonfluorescent cells gal80, and I-TRACE, we observe most cells express red or green indicate cells in which silencing was never established. Red and green fluorescence. However, many contiguous or individual cells exhibit channels are enhanced past linearity to establish the lack of either fluores- both red and green fluorescence (Fig. 6B). These are cells expe- cence in the cells. (C) Green fluorescence (within dotted border) indicates riencing new gains and losses (Fig. 6C), confirming the ongoing losses of silencing some time before analysis. (D) Green with a trace of red fluorescence indicates losses of silencing more recently than in C.(E)Red instability of heterochromatin-induced gene silencing. fluorescence in the 2 patches (arrows) are greater than in D and differ from WhilewehavepresentedadetailedanalysisoftheSwiM each other, indicating progressively more-recent losses of silencing. (F) Bright system in a single eye-antennal disk, the same types of cells are (maximal) red and green fluorescence indicate ongoing silencing. (G)Asingle found in all such imaginal discs we observed (SI Appendix, Figs. cell (arrow) expressing green fluorescence but not red fluorescence indicates a S2 and S3), and we believe the SwiM system is revealing common loss of silencing within a clone of cells that are still silencing. (H) Array of cells features of the fluctuating nature of heterochromatin-induced expressing red and not green fluorescence, indicating recent gains of silencing. gene silencing through development. Note that they are individual, but nonetheless cluster within the imaginal disk. (I) Individual green fluorescing cells (arrows) without red fluorescence in the Photoswitchable Fluorescent Proteins Confirm Gains and Losses Both middle of nonfluorescent clones indicate cells that were not silenced, then Exist and Are Postmitotic. Finally, we combined ubiquitously acquired silencing, then lost it again. (J) Contiguous cluster of cells expressing red fluorescence and not green fluorescence (within dotted border) suggest expressed gal4 and variegating gal80 with a gal4-controlled Kaede that all of the cells in this clone gained silencing at the same time, despite photoswitchable protein. Kaede is translated and matures as a being many cell divisions descendent from a common ancestor. (Magnifica- green fluorescing protein, but brief exposure to UV radiation tion, A,150×; Insets indicate magnification of B–J.) converts it to red fluorescence. The former state is relatively un- stable (half-life of 6 h), whereas the red form is very stable (half- life greater than 7 d) (102). We photoconverted intact animals, returned them to culture, and allowed them 24 h before dissection idea proposed by Ahmad and Henikoff (98), rather than in locus- and analysis. We once again observed the same types of cell events specific modulation. The existence of such gains of silencing is that we detected using G-TRACE. Yellow fluorescence using inarguable in purely undifferentiated mitotic systems, such as in Kaede occurs in cells in which gal80 silencing is complete, and S. cerevisiae and Schizosaccharomyces pombe, where individual corresponds to yellow fluorescence using G-TRACE. Non- cells with either expression state of reporter genes are capable of fluorescent cells using Kaede corresponds to no or green fluo- recreating colonies with both states (99, 100). rescence with G-TRACE. Red fluorescence indicates losses of

Bughio et al. PNAS | October 1, 2019 | vol. 116 | no. 40 | 20049 Downloaded by guest on September 28, 2021 AD

B E

CF

Fig. 6. Use of I-TRACE and Kaede with variegating gal80 to reveal recent switches in repression state. (A) Graphical representation of I-TRACE components in conjunction with variegating gal80 and ubiquitously expressed gal4. Temporal schematic shows that I-TRACE is sensitive to changes in the repression state of gal80, although it poorly discriminates between gains and losses of silencing. (B) Illustrative image of a third-instar eye-antennal imaginal disk, including region enlarged in C.(C) Highlight of region of active switching (dotted outline). (D) Graphical representation of photoswitchable Kaede in conjunction with variegating gal80 and ubiquitously expressed gal4. Temporal schematic shows that Kaede can detect and discriminate derepression of gal80 prior to the time of photoconversion (purple line) and new repression of gal80 after photoconversion. (E) Illustrative example including region enlarged in F.(F) Highlight of a clone of cells (dotted outline) or individual cells (arrows) exhibiting new gains of gal80 silencing. (Magnification, B, 100×; E, 120×; Insets indicate magnifi- cation of C and F.)

gal80 silencing at some point in development prior to UV pho- Conclusions and Further Thoughts. Our results indicate that the toconversion, just as green (or green-yellow) fluorescence does potential strength of heterochromatin is set at or before fertil- when using G-TRACE. With Kaede, strictly green fluorescence ization, incompletely established around gastrulation (60), then indicates gains of silencing at some point after UV photo- lost and gained stochastically and repeatedly during develop- conversion (Fig. 6D). The time between UV treatment and dis- ment. Such fluctuations of heterochromatin strength can occur section was longer than the time between RFP detection and GFP without consequence to the variegating genes, at least until the detection in G-TRACE, and so we expected to see more green genes are “queried” for silencing at the time of gene expression. fluorescent cells when using Kaede. We observed individual and In the case of white or yellow variegation this occurs in late pu- small clones of green fluorescent cells (Fig. 6 E and F), confirming pae, with the final PEV expressivity and patterns reflecting the the presence of gains of silencing in imaginal discs. Individual final outcome of gains and losses through development up until green fluorescent cells are unequivocally gains of silencing without that point. This is the logical way in which fields of various sizes mitosis (especially in those cells immediately behind the mor- (i.e., the amount of the eye to be pigmented, the patched vs. salt- phogenetic furrow and the final cell division in these tissues), just and-pepper patterns) can exist prior to when the gene expression as individual red fluorescent cells are individuals in which silencing that defines them has initiated. In our SwiM system, where the was lost without cell division. Small clones of green fluorescent ubiquitin promoter queries gene expression continuously, we cells no doubt share ancestry and are coordinately acting in their directly observe these fluctuations creating fields of cells, related gain of silencing. Again, this implies that these cells have had some by lineage, whose silencing potential has been recorded through change and are determined to gain silencing some hours (or cell development. The cause of the heterochromatin fluctuation is divisions) hence, suggesting the cause of gains and losses of si- yet unknown, but could conceivably be 1 or more of a number of lencing prefigures the actual gains and losses of silencing. factors. Fluctuations may be influenced by changes at the rDNA,

20050 | www.pnas.org/cgi/doi/10.1073/pnas.1909724116 Bughio et al. Downloaded by guest on September 28, 2021 which are somatically unstable (103), known to impact PEV Materials and Methods (54, 104, 105), and can cause the same phenomenon in yeast (53). Strains and Husbandry. The main SwiM strain contained variegating GAL80 + Loss of other repeat DNAs (106) that alter PEV expressivity in and a ubiquitous GAL4,wasy1 w67c23; Tp(P{w mC=tubP-GAL80ts}10)PEV.4; +mC 1 trans (107, 108) might also affect heterochromatin fluctuations. P{w =Act5C-gal4}17bF01/TM6B, Tb . The individual in Fig. 5 was y ts− PEV.4 + It is also consistent with our data that fluctuations in concen- w/Y; P{tubP-GAL80 }/ ; P{Act5C-gal4}/P{UAS-RFP}, P{UAS-FLP}, P{Ubi-(FRT.STOP.FRT)GFP}, the product of a cross between the SwiM strain + + trations or activities of soluble heterochromatin-establishing fac- and a G-TRACE strain (w*; P{w mC=UAS-RedStinger}6, P{w mC=UAS-FLP.Exel}3, + tors, such as Piwi-interacting RNAs, small-interfering RNAs, or P{w mC=Ubi-p63E(FRT.STOP)Stinger}15F2). In those cell clones in which gal80 is protein components of heterochromatin (e.g., concentrations of active, the phenotype is RFP− GFP− because of lack of activation of any gal4- dependent transcription. In those cell clones in which gal80 is repressed, the maternally deposited or zygotically expressed HP1) could underlie + + the patterns of PEV. Some or all of these factors might have been phenotype is RFP GFP , because of direct activation of RFP by Gal4 and FLP- mediated rearrangement and activation of GFP. identified previously as dose-limiting genetic modifiers of PEV, Full genotypes of other strains, and details on culture conditions, are although a simple mass action model is likely insufficient to fully available in SI Appendix. account for the behaviors of PEV that we and others have de- scribed (see ref. 109 for a detailed treatment). Dissection and Microscopy. Images of whole flies were taken with a Sony a7iii Viewing PEV as fluctuations in heterochromatin function attached to a Nikon SMZ-1500 microscope, illuminated with a Peak Plus rather than epigenetic information at gene promoters is a rela- Tactical LED Flashlight. Larvae were dissected in PBS and were visualized with a Zeiss AxioZoom.v16 and images captured with the Zeiss Axiocam 506-mono tively simple shift in our view, but is capable of explaining camera. Photography and quantification are described in SI Appendix. heretofore difficult-to-explain PEV phenomena. These include the following: 1) The instability of silencing in G1 (51) despite Dissection, Microscopy, and Fluorescence Detection. Adult male flies were the evident stability through G1 and S, 2) the correlations be- killed by exposure to ether. They were rinsed in 70% ethanol and placed in tween 2 variegating reporters in the same nucleus (55), 3) the PBS. Wings and legs were removed for better visualization on a a Zeiss correlation between PEV and the instability of repeat genes AxioZoom.v16 and images captured with the Zeiss Axiocam 506-mono camera. Quantitation for correlation experiments was performed using (110, 111), 4) why the temperature sensitive periods for PEV NIH image. Larvae were dissected in PBS, and were visualized with pho- suppression/enhancement do not correspond to the times of the tography and quantification are described in SI Appendix. Live third-instar presumed losses of silencing, and 5) our observation here that larvae were collected into a Petri dish containing PBS and visualized to select expressivity precedes activity of a variegating promoter. If this male larvae. These males were further sorted based on high and low GFP, GENETICS model is true, it renders moot the search for biochemical mech- excluding the gonads which were uniformly bright between siblings. Indi- vidual larvae were imaged as above and returned to individually labeled anisms to maintain locus-specific histone modification patterns vials. Adults were then scored for white+ eye pigment and correlated to through S-phase, because this type of maintenance may not exist their GFP fluorescence. for PEV. What appears as maintenance may instead be the suc- Nonfluorescence images of whole flies were taken with a Sony a7iii at- cesses or failures of heterochromatin to persistently reestablish tached to a Nikon SMZ-1500 microscope, illuminated with a Peak Plus Tactical silencing on closely neighboring genes. LED Flashlight. Notably, in a study of mouse PEV (50), ectopic recruitment of Statistical Analyses. In all cases, α was set to 0.05 prior to the study, and we HP1 induced silencing (49) that could be monitored for stability report whether the null hypotheses (H0) could be rejected based on that in the form of self-sustained mitotic maintenance. The empirical criterion. Specific test statistics are reported in the text and figure legends. rates derived from these experiments did not themselves explain Justifications for specific tests are described in SI Appendix. the stability of silencing in either the transgene or natural het- erochromatin, and some hypothetical factor was required to ACKNOWLEDGMENTS. This work is supported by National Institutes of Health Office of the Director’s Transformative Research Award R01 properly model self-sustaining heterochromatin-induced gene GM076092. Special support was provided by Dr. Patrick Lyons. Core support silencing (50, 87). We suggest that it is not possible to model was given by the University of Arizona Cancer Center, Grant P30 CA023074, stable heterochromatin-induced gene silencing without a bona and core services at the University of Arizona provided by the Office of the Vice President for Research. We thank K. G. Golic for suggesting the name fide block of heterochromatin that is constantly reestablishing “SwiM,” and S. L. Eckert, K. G. Golic, M. Golic, and C. Huckell for support silencing (51, 52). during the finalization of this work.

1. J. W. Gowen, E. H. Gay, Chromosome constitution and behavior in eversporting and 12. G. H. Karpen, R. C. Allshire, The case for epigenetic effects on centromere identity mottling in Drosophila melanogaster. Genetics 19, 189–208 (1934). and function. Trends Genet. 13, 489–496 (1997). 2. J. B. Spofford, “Position-effect variegation in Drosophila” in The Genetics and Bi- 13. B. C. Williams, T. D. Murphy, M. L. Goldberg, G. H. Karpen, Neocentromere activity ology of Drosophila, M. Ashburner, E. Novitski, Eds. (Academic Press, 1976), vol. 1c, of structurally acentric mini-chromosomes in Drosophila. Nat. Genet. 18,30–37 chap. 24, pp. 955–1019. (1998). 3. K. S. Weiler, B. T. Wakimoto, Heterochromatin and gene expression in Drosophila. 14. K. A. Maggert, G. H. Karpen, Acquisition and metastability of centromere identity Annu. Rev. Genet. 29, 577–605 (1995). and function: Sequence analysis of a human neocentromere. Genome Res. 10, 725– 4. S. C. Elgin, G. Reuter, Position-effect variegation, heterochromatin formation, 728 (2000). and gene silencing in Drosophila. Cold Spring Harb. Perspect. Biol. 5, a017780 15. M. D. Anway, A. S. Cupp, M. Uzumcu, M. K. Skinner, Epigenetic transgenerational (2013). actions of endocrine disruptors and male fertility. Science 308, 1466–1469 (2005). 5. J. Locke, M. A. Kotarski, K. D. Tartof, Dosage-dependent modifiers of position effect 16. C. Lang-Mladek et al., Transgenerational inheritance and resetting of stress-induced variegation in Drosophila and a mass action model that explains their effect. Ge- loss of epigenetic gene silencing in Arabidopsis. Mol. Plant 3, 594–602 (2010). netics 120, 181–198 (1988). 17. R. J. Schmitz et al., Transgenerational epigenetic instability is a source of novel 6. M. F. Lyon, A. G. Searle, C. E. Ford, S. Ohno, A mouse translocation suppressing sex- methylation variants. Science 334, 369–373 (2011). linked variegation. Cytogenetics 3, 306–323 (1964). 18. C. Becker, D. Weigel, Epigenetic variation: Origin and transgenerational inheritance. 7. S. Ohno, M. F. Lyon, Cytological study of searle’s X-autosome translocation in Mus Curr. Opin. Plant Biol. 15, 562–567 (2012). musculus. Chromosoma 16,90–100 (1965). 19. D. Crews et al., Epigenetic transgenerational inheritance of altered stress responses. 8. P. S. Harper, Mary Lyon and the hypothesis of random X chromosome inactivation. Proc. Natl. Acad. Sci. U.S.A. 109, 9143–9148 (2012). Hum. Genet. 130, 169–174 (2011). 20. E. Heard, R. A. Martienssen, Transgenerational epigenetic inheritance: Myths and 9. A. C. Ferguson-Smith, Genomic imprinting: The emergence of an epigenetic para- mechanisms. Cell 157,95–109 (2014). digm. Nat. Rev. Genet. 12, 565–575 (2011). 21. M. Trerotola, V. Relli, P. Simeone, S. Alberti, Epigenetic inheritance and the missing 10. U. Nur, Heterochromatization and euchromatization of whole genomes in scale heritability. Hum. Genomics 9, 17 (2015). insects (Coccoidea: Homoptera). Dev. Suppl. 1990,29–34 (1990). 22. T. M. Valtonen, K. Kangassalo, M. Pölkki, M. J. Rantala, Transgenerational effects of 11. G. Herrick, J. Seger, Imprinting and paternal genome elimination in insects. Results parental larval diet on offspring development time, adult body size and pathogen Probl. Cell Differ. 25,41–71 (1999). resistance in Drosophila melanogaster. PLoS One 7, e31611 (2012).

Bughio et al. PNAS | October 1, 2019 | vol. 116 | no. 40 | 20051 Downloaded by guest on September 28, 2021 23. F. C. Zucchi et al., Maternal stress induces epigenetic signatures of psychiatric and 61. W. K. Baker, A clonal system of differential gene activity in Drosophila. Dev. Biol. 16, neurological diseases in the offspring. PLoS One 8, e56967 (2013). 1–17 (1967). 24. A. Soubry, C. Hoyo, R. L. Jirtle, S. K. Murphy, A paternal environmental legacy: Evi- 62. C. P. Bishop, Evidence for intrinsic differences in the formation of chromatin do- dence for epigenetic inheritance through the male germ line. Bioessays 36, 359–371 mains in Drosophila melanogaster. Genetics 132, 1063–1069 (1992). (2014). 63. K. A. Maggert, K. G. Golic, The Y chromosome of Drosophila melanogaster exhibits 25. A. Soubry, Epigenetic inheritance and evolution: A paternal perspective on dietary chromosome-wide imprinting. Genetics 162, 1245–1258 (2002). influences. Prog. Biophys. Mol. Biol. 118,79–85 (2015). 64. C. J. Evans et al., G-TRACE: Rapid Gal4-based cell lineage analysis in Drosophila. Nat. 26. B. Lemos, A. T. Branco, D. L. Hartl, Epigenetic effects of polymorphic Y chromosomes Methods 6, 603–605 (2009). modulate chromatin components, immune response, and sexual conflict. Proc. Natl. 65. C. Alabert, A. Groth, Chromatin replication and epigenome maintenance. Nat. Rev. Acad. Sci. U.S.A. 107, 15826–15831 (2010). Mol. Cell Biol. 13, 153–167 (2012). 27. M. Pigliucci, G. Müller, Konrad Lorenz Institute for Evolution and Cognition Re- 66. T. Rohlf et al., Modeling the dynamic epigenome: From histone modifications to- search, Evolution, the Extended Synthesis (MIT Press, Cambridge, MA, 2010). wards self-organizing chromatin. Epigenomics 4, 205–219 (2012).   28. J. P. Lim, A. Brunet, Bridging the transgenerational gap with epigenetic memory. 67. A. Boskovic, O. J. Rando, Transgenerational epigenetic inheritance. Annu. Rev. Trends Genet. 29, 176–186 (2013). Genet. 52,21–41 (2018).  29. J. G. Falls, D. J. Pulford, A. A. Wylie, R. L. Jirtle, Genomic imprinting: Implications for 68. K. Skvortsova, N. Iovino, O. Bogdanovic, Functions and mechanisms of epigenetic human disease. Am. J. Pathol. 154, 635–647 (1999). inheritance in animals. Nat. Rev. Mol. Cell Biol. 19, 774–790 (2018). 30. D. C. Dolinoy, J. R. Weidman, R. L. Jirtle, Epigenetic gene regulation: Linking early 69. R. A. Greenstein, B. Al-Sady, Epigenetic fates of gene silencing established by het- – developmental environment to adult disease. Reprod. Toxicol. 23, 297–307 (2007). erochromatin spreading in cell identity and genome stability. Curr. Genet. 65, 423 31. J. Bohacek, I. M. Mansuy, Epigenetic inheritance of disease and disease risk. Neu- 428 (2019). ropsychopharmacology 38, 220–236 (2013). 70. M. F. Perez, B. Lehner, Intergenerational and transgenerational epigenetic in- – 32. K. A. Maggert, Genetics: Polymorphisms, epigenetics, and something in between. heritance in animals. Nat. Cell Biol. 21, 143 151 (2019). Genet. Res. Int. 2012, 867951 (2012). 71. K. G. Golic, S. Lindquist, The FLP recombinase of yeast catalyzes site-specific re- – 33. C. Deans, K. A. Maggert, What do you mean, “epigenetic”? Genetics 199, 887–896 combination in the Drosophila genome. Cell 59, 499 509 (1989). (2015). 72. K. G. Golic, M. M. Golic, Engineering the Drosophila genome: Chromosome re- – 34. K. H. Seong, T. Maekawa, S. Ishii, Inheritance and memory of stress-induced epige- arrangements by design. Genetics 144, 1693 1711 (1996). nome change: Roles played by the ATF-2 family of transcription factors. Genes Cells 73. R. Drysdale; FlyBase Consortium, FlyBase : A database for the Drosophila research – 17, 249–263 (2012). community. Methods Mol. Biol. 420,45 59 (2008). 35. J. C. Aldrich, K. A. Maggert, Transgenerational inheritance of diet-induced genome 74. G. Reuter, R. Dorn, H. J. Hoffmann, Butyrate sensitive suppressor of position-effect – rearrangements in Drosophila. PLoS Genet. 11, e1005148 (2015). variegation mutations in Drosophila melanogaster. Mol. Gen. Genet. 188, 480 485 36. V. Tran, C. Lim, J. Xie, X. Chen, Asymmetric division of Drosophila male germline (1982). stem cell shows asymmetric histone distribution. Science 338, 679–682 (2012). 75. C. W. Hinton, The behavior of an unstable ring chromosome of Drosophila mela- – 37. K. Ahmad, S. Henikoff, No strand left behind. Science 361, 1311–1312 (2018). nogaster. Genetics 40, 951 961 (1955). 38. N. Petryk et al., MCM2 promotes symmetric inheritance of modified histones during 76. W. Janning, Gynandromorph fate maps in Drosophila. Results Probl. Cell Differ. 9, 1–28 (1978). DNA replication. Science 361, 1389–1392 (2018). 77. E. M. Underwood, F. R. Turner, A. P. Mahowald, Analysis of cell movements and fate 39. C. Yu et al., A mechanism for preventing asymmetric histone segregation onto mapping during early embryogenesis in Drosophila melanogaster. Dev. Biol. 74, replicating DNA strands. Science 361, 1386–1389 (2018). 286–301 (1980). 40. S. Urieli-Shoval, Y. Gruenbaum, J. Sedat, A. Razin, The absence of detectable 78. A. P. Mahowald, P. A. Hardy, Genetics of Drosophila embryogenesis. Annu. Rev. methylated bases in Drosophila melanogaster DNA. FEBS Lett. 146, 148–152 (1982). Genet. 19, 149–177 (1985). 41. M. Schaefer, F. Lyko, Lack of evidence for DNA methylation of Invader4 retroele- 79. J. C. Eissenberg, Position effect variegation in Drosophila: Towards a genetics of ments in Drosophila and implications for Dnmt2-mediated epigenetic regulation. chromatin assembly. Bioessays 11,14–17 (1989). Nat. Genet. 42, 920–921, author reply 921 (2010). 80. M. Demerec, The Biology of Drosophila (Cold Spring Harbor Laboratory Press, 42. S. Phalke et al., Retrotransposon silencing and telomere integrity in somatic cells of Plainview, NY, ed. Facsim., 1994), p. x, 632 pp. Drosophila depends on the cytosine-5 methyltransferase DNMT2. Nat. Genet. 41, 81. K. Yuan, P. H. O’Farrell, TALE-light imaging reveals maternally guided, H3K9me2/3- 696–702 (2009). independent emergence of functional heterochromatin in Drosophila embryos. 43. V. Krauss, G. Reuter, DNA methylation in Drosophila—A critical evaluation. Prog. Genes Dev. 30, 579–593 (2016). Mol. Biol. Transl. Sci. 101, 177–191 (2011). 82. C. A. Seller, C. Y. Cho, P. H. O’Farrell, Rapid embryonic cell cycles defer the estab- 44. M. Hampsey, D. Reinberg, Tails of intrigue: Phosphorylation of RNA polymerase II lishment of heterochromatin by Eggless/SetDB1 in Drosophila. Genes Dev. 33, 403– mediates histone methylation. Cell 113, 429–432 (2003). 417 (2019). 45. U. Grossniklaus et al., Transgenerational epigenetic inheritance: How important is it? 83. B. A. Edgar, G. Schubiger, Parameters controlling transcriptional activation during Nat. Rev. Genet. 14, 228–235 (2013). early Drosophila development. Cell 44, 871–877 (1986). 46. M. Ptashne, Epigenetics: Core misconcept. Proc. Natl. Acad. Sci. U.S.A. 110, 7101– 84. J. C. Eissenberg, Y. W. Ge, T. Hartnett, Increased phosphorylation of HP1, a 7103 (2013). heterochromatin-associated protein of Drosophila, is correlated with heterochro- 47. S. Henikoff, S. Ramachandran, Pioneers invade the nucleosomal landscape. Mol. Cell matin assembly. J. Biol. Chem. 269, 21315–21321 (1994). – 71, 193 194 (2018). 85. U. Nur, Reversal of heterochromatization and the activity of the paternal chromo- 48. S. Ramachandran, S. Henikoff, MINCE-Seq: Mapping in vivo nascent chromatin with some set in the male mealy bug. Genetics 56, 375–389 (1967). – EdU and sequencing. Methods Mol. Biol. 1832, 159 168 (2018). 86. H. Steller, V. Pirrotta, Expression of the Drosophila white gene under the control of 49. Y. Li, J. R. Danzer, P. Alvarez, A. S. Belmont, L. L. Wallrath, Effects of tethering HP1 the hsp70 heat shock promoter. EMBO J. 4, 3765–3772 (1985). – to euchromatic regions of the Drosophila genome. Development 130, 1817 1824 87. C. Hodges, G. R. Crabtree, Dynamics of inherently bounded histone modification (2003). domains. Proc. Natl. Acad. Sci. U.S.A. 109, 13296–13301 (2012). 50. N. A. Hathaway et al., Dynamics and memory of heterochromatin in living cells. Cell 88. D. Richards, S. Berry, M. Howard, Illustrations of mathematical modeling in biology: – 149, 1447 1460 (2012). Epigenetics, meiosis, and an outlook. Cold Spring Harb. Symp. Quant. Biol. 77, 175– 51. K. Ahmad, K. G. Golic, Somatic reversion of chromosomal position effects in Dro- 181 (2012). – sophila melanogaster. Genetics 144, 657 670 (1996). 89. P. A. Steffen, J. P. Fonseca, L. Ringrose, Epigenetics meets mathematics: Towards a 52. T. H. Cheng, M. R. Gartenberg, Yeast heterochromatin is a dynamic structure that quantitative understanding of chromatin biology. Bioessays 34, 901–913 (2012). – requires silencers continuously. Genes Dev. 14, 452 463 (2000). 90. H. Zhang, X. J. Tian, A. Mukhopadhyay, K. S. Kim, J. Xing, Statistical mechanics model 53. A. H. Michel, B. Kornmann, K. Dubrana, D. Shore, Spontaneous rDNA copy number for the dynamics of collective epigenetic histone modification. Phys. Rev. Lett. 112, – variation modulates Sir2 levels and epigenetic gene silencing. Genes Dev. 19, 1199 068101 (2014). 1210 (2005). 91. K. Ang et al., Mediator acts upstream of the transcriptional activator Gal4. PLoS Biol. 54. S. Paredes, K. A. Maggert, Ribosomal DNA contributes to global chromatin regula- 10, e1001290 (2012). tion. Proc. Natl. Acad. Sci. U.S.A. 106, 17829–17834 (2009). 92. B. McClintock, The significance of responses of the genome to challenge. Science 55. Y. Mano, T. J. Kobayashi, J. Nakayama, H. Uchida, M. Oki, Single cell visualization of 226, 792–801 (1984). yeast gene expression shows correlation of epigenetic switching between multiple 93. L. Piacentini et al., Transposons, environmental changes, and heritable induced heterochromatic regions through multiple generations. PLoS Biol. 11, e1001601 phenotypic variability. Chromosoma 123, 345–354 (2014). (2013). 94. L. Fanti, L. Piacentini, U. Cappucci, A. M. Casale, S. Pimpinelli, Canalization by se- 56. W. K. Baker, Position-effect variegation. Adv. Genet. 14, 133–169 (1968). lection of de Novo induced mutations. Genetics 206, 1995–2006 (2017). 57. K. D. Tartof, M. Bremer, Mechanisms for the construction and developmental control 95. K. Melcher, H. E. Xu, Gal80-Gal80 interaction on adjacent Gal4p binding sites is re- of heterochromatin formation and imprinted chromosome domains. Dev. Suppl. quired for complete GAL gene repression. EMBO J. 20, 841–851 (2001). 1990,35–45 (1990). 96. O. Egriboz et al., Self-association of the Gal4 inhibitor protein Gal80 is impaired by 58. B. Y. Lu, C. P. Bishop, J. C. Eissenberg, Developmental timing and tissue specificity of Gal3: Evidence for a new mechanism in the GAL gene switch. Mol. Cell. Biol. 33, heterochromatin-mediated silencing. EMBO J. 15, 1323–1332 (1996). 3667–3674 (2013). 59. A. Fjose, L. C. Polito, U. Weber, W. J. Gehring, Developmental expression of the 97. V. V. Verkhusha et al., High stability of discosoma DsRed as compared to aequorea white locus of Drosophila melanogaster. EMBO J. 3, 2087–2094 (1984). EGFP. Biochemistry 42, 7879–7884 (2003). 60. B. Y. Lu, J. Ma, J. C. Eissenberg, Developmental regulation of heterochromatin- 98. K. Ahmad, S. Henikoff, Modulation of a transcription factor counteracts hetero- mediated gene silencing in Drosophila. Development 125, 2223–2234 (1998). chromatic gene silencing in Drosophila. Cell 104, 839–847 (2001).

20052 | www.pnas.org/cgi/doi/10.1073/pnas.1909724116 Bughio et al. Downloaded by guest on September 28, 2021 99. L. Pillus, J. Rine, Epigenetic inheritance of transcriptional states in S. cerevisiae. Cell 105. J. Zhou et al., Y chromosome mediates ribosomal DNA silencing and modulates the 59, 637–647 (1989). chromatin state in Drosophila. Proc. Natl. Acad. Sci. U.S.A. 109, 9941–9946 (2012). 100. L. L. Sandell, D. E. Gottschling, V. A. Zakian, Transcription of a yeast telomere alle- 106. K. H. Wei, J. K. Grenier, D. A. Barbash, A. G. Clark, Correlated variation and pop- viates telomere position effect without affecting chromosome stability. Proc. Natl. ulation differentiation in satellite DNA abundance among lines of Drosophila – Acad. Sci. U.S.A. 91, 12061–12065 (1994). melanogaster. Proc. Natl. Acad. Sci. U.S.A. 111,18793 18798 (2014). 101. J. A. Bosch, T. M. Sumabat, I. K. Hariharan, Persistence of RNAi-mediated knockdown 107. J. B. Spofford, R. DeSalle, Nucleolus organizer-suppressed position-effect variegation in Drosophila melanogaster. Genet. Res. 57, 245–255 (1991). in Drosophila complicates mosaic analysis yet enables highly sensitive lineage trac- 108. M. Gatti, S. Pimpinelli, Functional elements in Drosophila melanogaster hetero- ing. Genetics 203, 109–118 (2016). chromatin. Annu. Rev. Genet. 26, 239–275 (1992). 102. M. Tomura et al., Monitoring cellular movement in vivo with photoconvertible 109. S. Henikoff, Dosage-dependent modification of position-effect variegation in Dro- “ ” fluorescence protein Kaede transgenic mice. Proc. Natl. Acad. Sci. U.S.A. 105, sophila. Bioessays 18, 401–409 (1996). – 10871 10876 (2008). 110. J. C. Peng, G. H. Karpen, Heterochromatic genome stability requires regulators of 103. S. Cohen, K. Yacobi, D. Segal, Extrachromosomal circular DNA of tandemly repeated histone H3 K9 methylation. PLoS Genet. 5, e1000435 (2009). genomic sequences in Drosophila. Genome Res. 13, 1133–1145 (2003). 111. J. C. Aldrich, K. A. Maggert, Simple quantitative PCR approach to reveal naturally 104. K. Larson et al., Heterochromatin formation promotes longevity and represses ri- occurring and mutation-induced repetitive sequence variation on the Drosophila Y bosomal RNA synthesis. PLoS Genet. 8, e1002473 (2012). chromosome. PLoS One 9, e109906 (2014). GENETICS

Bughio et al. PNAS | October 1, 2019 | vol. 116 | no. 40 | 20053 Downloaded by guest on September 28, 2021