Histone variant H2A.Z mutant suppresses the senescence-associated secretory phenotype

Histone variant H2A.Z mutant suppresses the senescence-associated secretory phenotype | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Histone variant H2A.Z mutant suppresses the senescence-associated secretory phenotype Zong Ming Chua, Hiroshi Tanaka, Adrianna Abele, Adarsh Rajesh, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8818498/v1 This work is licensed under a CC BY 4.0 License Status: Under Revision Version 1 posted 13 You are reading this latest preprint version Abstract Cellular senescence features a durable cell-cycle arrest and a pro-inflammatory senescence-associated secretory phenotype (SASP) driven in part by chromatin remodeling. The histone variant H2A.Z plays an essential role in regulating gene expression through modulating nucleosome dynamics and is known to regulate the expression of cell cycle genes during the early stages of cellular senescence. However, how the intrinsic stability of H2A.Z-containing nucleosomes plays a role in the establishment of the senescent phenotype remains unexplored. To investigate this, we employed H2A.Z R80C, a H2A.Z mutant that destabilizes nucleosomes by disrupting histone-DNA interactions. We observed that expression of H2A.Z R80C causes suppression of SASP in senescent primary human fibroblasts, without affecting expression of cell cycle genes. H2A.Z knockdown did not suppress the SASP, demonstrating that SASP suppression is likely due to altered stability of H2A.Z-containing nucleosomes rather than loss of H2A.Z function. Mechanistically, SASP suppression is linked to decreased H3K27ac at SASP gene loci. These findings offer a novel avenue for understanding and manipulating the SASP during aging and other senescence-related pathologies. Biological sciences/Cell biology Biological sciences/Genetics Biological sciences/Molecular biology Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 INTRODUCTION Cellular senescence, a state of irreversible cell cycle arrest triggered by both extrinsic and intrinsic stimuli, is a complex biological process marked by significant alterations in chromatin structure, epigenetic modifications, metabolism and gene expression 1 – 6 . Senescent cells remain viable after stable cell cycle arrest, and their burden is known to increase with age in many vertebrates 7 . Senescent cells produce a specific pro-inflammatory senescence-associated secretory phenotype (SASP) which consists of proinflammatory cytokines such as IL-6 and IL-1, chemokines such as IL-8, and extracellular proteases such as MMP-1. This SASP can improve tissue repair and regeneration, anti-tumor immunity, and cellular plasticity 8 – 13 . However, sustained and chronic expression of the SASP, especially in older organisms, can result in the spread of senescence and impaired tissue function through paracrine signaling 14 . Chronic SASP signaling can also lead to a pro-inflammatory and immunosuppressive microenvironment that drives tumorigenesis 15 . Accordingly, removal of senescent cells has been shown to improve tissue function and promote healthy aging 16 , 17 . Induction of cellular senescence and the SASP takes place concurrently with changes in nucleosome composition and localization of histone variants, which in turn regulate key age-related and senescence phenotypes 18 – 23 . Among the histone variants that play a role in regulating the senescent phenotype is H2A.Z, a variant of core histone H2A, which has been shown to accumulate with age and is linked to age-related cognitive decline in mice 24 . H2A.Z has a roughly 60% degree of sequence similarity with canonical H2A, with a higher proportion of residues in the N- and C-terminal tail domain showing differences compared to the more conserved core regions 25 . In some contexts, H2A.Z has a pro-proliferative, anti-senescence role. The p53-dependent localization of the histone variant H2A.Z by the p400 complex represses basal p21 transcription in proliferating cells, and upon DNA damage, H2A.Z is evicted to permit Tip60 recruitment and p21-mediated proliferation arrest 26 . Additionally, H2A.Z supports the expression of proliferation-associated genes such as Myc and Ki-67, further reinforcing its role in cell cycle progression and senescence avoidance 27 . In pancreatic ductal adenocarcinoma (PDAC), overexpression of H2A.Z enhances tumor growth and confers chemoresistance by suppressing senescence, again highlighting a pro-proliferative function 28 . On the other hand, H2A.Z is expressed throughout the cell cycle and in non-proliferating cells, unlike canonical H2A which is predominantly expressed during S phase 29 . Accordingly, it has been shown that the growth-arrested senescent phenotype is associated with the site-specific deposition of H2A.Z, suggesting a role for H2A.Z in non-proliferating senescent cells. The histone variant H2A.Z directly impacts the structure and the stability of assembled nucleosomes, altering the activity of chromatin remodelers and the accessibility of the underlying DNA, thereby affecting diverse DNA-related processes such as repair, recombination, and transcription 30 – 33 . H2A.Z is localized at the transcription start site (TSS) of actively transcribed genes, and its presence plays a role in regulating the activity of poised genes which are marked by both activating and repressive histone marks 34 , 35 . Homotypic nucleosomes containing pairs of H2A.Z histones show decreased nucleosomal stability, while H2A.Z containing nucleosomes show altered dynamics, mobility, and increased accessibility of the local chromatin to transcriptional machinery 36 – 39 . Given both the importance of cellular senescence and the accompanying SASP in aging, and the profound changes in transcriptome that take place during the onset of cellular senescence and the proposed role of H2A.Z in control of gene expression through nucleosome stability, we sought to investigate how modulating the stability of H2A.Z containing nucleosomes would impact the SASP, potentially providing new avenues to ameliorate age-associated disease phenotypes. To understand how the stability of H2A.Z-containing nucleosomes plays a role in the induction of cellular senescence, we examined the effects of a specific histone mutant H2A.Z R80C in senescent cells. H2A.Z R80C was previously identified as a histone mutant recurrently found in some cancers (a so-called oncohistone mutation) and has been shown to decrease the stability of the core nucleosomal complex 40 , 41 . Modeling of the H2A.Z R80C mutant shows that removal of the arginine 80 residue reduces the affinity between the H2A.Z histone and the adjacent DNA strand, as the arginine interacts favorably with the negatively charged DNA backbone while also extending into the minor groove of the DNA strand to stabilize the nucleosome 41 . We used H2A.Z R80C as a tool to probe the relationship between nucleosome stability and phenotypes of cell senescence. RESULTS H2A.Z R80C alters nucleosome stability To investigate the relationship between H2A.Z stability and phenotypes of senescence, we employed the H2A.Z R80C mutant that subtly destabilizes nucleosomes 40 . The H2A.Z R80C substitution mutation replaces an arginine residue with cysteine, which is modeled to destabilize the nucleosome due to the loss of the favorable arginineDNA interaction in an assembled nucleosome (Fig. 1 A). We carried out biochemical characterization of lentivirus-mediated ectopically expressed H2A.Z wild type (WT) or mutant H2A.Z R80C, both with an N-terminal HA-tag or GFP-tag, in IMR-90 primary human fibroblasts. We confirmed that ectopic H2A.Z WT and R80C were expressed at comparable levels, based on western blotting and quantification of the immunofluorescent signals (Fig. 1 B, 1 C). More than 99% of cells showed clear ectopic expression of H2A.Z WT or R80C. Several lines of evidence indicated that both H2A.Z WT and R80C are incorporated into nucleosomal chromatin. First, the fluorescent signal of the HA-tag co-localized with DAPI and H3K27me3 in both interphase and mitotic chromatin (Fig. 1 B, 1 D). Co-localization of the WT and mutant histones on mitotic chromosomes strongly suggests incorporation of both variants into chromatin. A DNA digestion assay targeted by anti-HA antibodies using Tn5 transposase in IMR90 cells expressing either HA-tagged H2A.Z WT or R80C histones revealed clear nucleosomal ladder patterns, indicating that the ectopically expressed histones are both incorporated into chromatin (Fig. 1 E). To assess the protein mobility of H2A.Z WT and R80C in the nucleus, we carried out fluorescence recovery after photobleaching (FRAP) to measure the kinetics of nuclear diffusion of GFP-tagged H2A.Z WT and R80C. The H2A.Z R80C mutant showed a more rapid fluorescence recovery rate compared to H2A.Z WT, indicating that H2A.Z R80C is more mobile within the nucleus, consistent with decreased nucleosome stability 40 (Fig. 1 F). In sum, H2A.Z WT and R80C are comparably incorporated into nucleosomes and chromatin at steady state, but R80C is more mobile consistent with decreased nucleosome stability. H2A.Z R80C mutant suppresses induction of a subset of the SASP Having characterized H2A.Z WT and R80C mutant in proliferating IMR-90 cells, we next investigated their impact on senescence-associated phenotypes. Upon irradiation, H2A.Z WT and H2A.Z R80C exhibited a comparable large and flattened morphology characteristic of senescence within 14 days (Fig. 2 A). Consistent with the induction of senescence, for both H2A.Z WT and R80C cells we observed a clear downregulation of genes associated with cellular proliferation and cell cycle progression following irradiation (Fig. 2 B). A principal component (PC) analysis of mRNA-seq at six time points, 0, 5, 8, 10, 12, and 14 days post-irradiation, revealed that roughly a quarter of transcriptomic variance (25.9%) was explained by PC1 with a clear separation of samples by time post-irradiation (Supplementary Fig. 1). H2A.Z R80C tended to separate from empty and H2A.Z WT along PC2. An unbiased clustering analysis of differentially expressed genes at day 14 post-irradiation identified four distinct gene clusters when comparing empty and H2A.Z WT controls and R80C-expressing cells (Fig. 2 C). Among these clusters, cluster 1 represented genes with expression levels that were unaffected during senescence in both empty and H2A.Z WT cells, but were upregulated in H2A.Z R80C cells. Cluster 2 included genes that were upregulated during senescence in empty, H2A.Z WT and H2A.Z R80C cells, albeit more so and less so in R80C in some cases. Cluster 3 contained genes downregulated during senescence in empty and H2A.Z WT cells, with further decrease in H2A.Z R80C. Finally, cluster 4 represented genes markedly upregulated during senescence in empty and H2A.Z WT cells, but less so in H2A.Z R80C cells. Functional enrichment analysis of each gene cluster revealed distinct biological associations: cluster 1 was enriched with genes related to glycerolipid and cholesterol metabolism; cluster 2 with genes involved in cytokine and inflammatory responses; cluster 3 with genes related to extracellular matrix organization; and cluster 4 with genes encoding cell surface proteins (Supplementary Fig. 2). Further analysis of these differentially expressed genes using the STRING database supports the results of this functional enrichment analysis in identifying both a group of pro-inflammatory genes and a cluster of collagen and matrix metalloproteinase genes that are downregulated in H2A.Z R80C expressing cells (Supplementary Fig. 3). We validated these functional enrichment analyses with a gene set enrichment analysis (GSEA) on the same set of differentially expressed genes and observed that expression of genes in pathways related to inflammation was repressed in H2A.Z R80C samples at day 14 post irradiation compared to H2A.Z WT (Supplementary Fig. 4). Repression of genes related to inflammation was also apparent in parts of cluster 2 (IL1A, IL1B, IL6, CXCL3, CXCL5) and cluster 4 (CXCL1, CXCL8) (Supplementary Fig. 2). These functional enrichment results indicate that there are two dominant repressive gene expression changes associated with expression of the H2A.Z R80C mutant during senescence, with one centered on the repression of a pro-inflammatory gene signature and the other on the repression of genes associated with the extracellular matrix. Given that the SASP is enriched in both inflammatory cytokines/chemokines and extracellular matrix components, we examined whether the gene signatures from our unbiased analysis reflect the known SASP signature 13 . We examined the expression of several key SASP genes and found that their expression in senescence was lower in H2A.Z R80C expressing cells compared to H2A.Z WT (Fig. 2 D). Expression of these key SASP genes was significantly downregulated in senescent H2A.Z R80C expressing cells (Fig. 2 E). Importantly, the expression of H2A.Z R80C did not affect the expression of canonical interferon-stimulated genes (ISGs) that were upregulated during senescence (Fig. 2 F), confirming that H2A.Z R80C specifically affects a subset of the SASP. We verified suppression of key SASP markers IL1A , IL6 , and IL8 / CXCL8 in H2A.Z R80C expressing cells by qPCR, while the expression of cell cycle genes CDKN2A and CDKN1A was not significantly impaired (Fig. 2 G). In summary, these results indicate that expression of H2A.Z R80C specifically suppresses genes associated with the SASP, including chemokines and cytokines involved in the proinflammatory phenotype of senescence. SASP suppression by H2A.Z R80C is specific to senescence and independent of loss of H2A.Z function Next, we set out to investigate the mechanism by which H2A.Z R80C suppresses the SASP. First, we asked whether suppression of SASP is specific to the R80C substitution or can also be achieved by other amino-acid substitutions. To do this, we tested various H2A.Z R80 substitutions for their effect on senescence, including bulky residues (R80W), and those with altered (R80D, R80E) or similar (R80K) charge properties to the original arginine residue. Expression of all mutants was confirmed by RT-qPCR (Fig. 3 A). We observed similar reductions in SASP gene expression during senescence across all mutations (Fig. 3 B), indicating that SASP suppression is agnostic to the specific identity of the residue in the substitution. This is consistent with structural studies showing that R80 makes specific contact with DNA, and suggests that this interaction is disrupted by any mutation 40 . Second, we asked whether H2A.Z R80C blocks expression of the canonical SASP genes only when they were induced as part of the senescence program. To investigate this, we added recombinant IL1A, a pro-inflammatory cytokine, into the media of senescent IMR90 cells and observed expression of SASP components IL1A , IL6 and IL8 by RT-qPCR. Under these conditions, we observed the reduction in IL1A , IL6 and IL8 expression by H2A.Z R80C was largely absent (Fig. 3 C). This indicates that H2A.Z R80C does not affect the capacity of cells to express IL6 and IL8 in response to other non-senescence triggers. Its suppressive effect on these SASP components is specific to the senescence program. Finally, we wanted to know whether the suppression of the SASP by H2A.Z R80C is due to a loss of H2A.Z function. To test this, we asked whether suppression of the SASP was also achieved by H2A.Z knockdown. We performed siRNA knockdown of H2A.Z in senescent cells using a pool of three different siRNAs targeting H2A.Z.1. H2A.Z.2 (H2AFV) was not included in this analysis due to its substantially lower expression levels relative to H2A.Z.1 in our system, making it less likely to contribute significantly to the observed effects (Supplementary Fig. 5). Using RT-qPCR, we verified knockdown efficiency of ~ 98% based on transcript level of both total H2A.Z as well as ectopically expressed HA-tagged H2A.Z (Fig. 3 D). As expected, knockdown of ectopic H2A.Z R80C abrogated the reduction of the SASP by H2A.Z R80C, confirming that H2A.Z R80C in senescent cells represses the SASP (Fig. 3 E). However, knockdown of H2A.Z did not suppress the SASP gene expression in both Empty and H2A.Z WT cells. Indeed, knockdown of H2A.Z increased expression of SASP genes in several cases, suggesting a transcription repressive role in some cases. Regardless, these results indicate that the suppression of the SASP during senescence by H2A.Z R80C is not due to a loss of H2A.Z function. H2A.Z R80C shows altered epigenetic signatures in image-based deep learning classification Given the role of histone modifications in control of gene expression, we hypothesized that the altered transcriptomic program of H2A.Z R80C-expressing cells is linked to changes in epigenetic signatures. We employed deep learning techniques to examine differences in various epigenetic signatures in proliferating and senescent IMR-90 cells using immunofluorescent images of key regulatory histone modifications: H3K27ac, H3K27me3, H3K4me1, and H3K56ac 42 , 43 . Comparisons among different neural networks revealed that the ResNet50 pre-trained convolutional neural networks performed best to differentiate between proliferating cells based on immunofluorescence staining of H3K27ac, H3K27me3, H3K4me1, H3K56ac and DAPI staining (Fig. 4 A). The trained ResNet50 model robustly distinguished various nuclear epigenetic marks, while highlighting a visually similar signature between the enhancer marks H3K27ac and H3K4me1 compared to other modifications, demonstrating the power of this machine learning approach to detect biological differences in epigenetic signatures between individual cells (Fig. 4 B). We observed that, based solely on DAPI staining, a separate trained ResNet model was able to robustly distinguish between proliferating and senescent cells (Fig. 4 C). While DAPI alone could not clearly distinguish R80C cells from others, H3K27ac, H3K27me3, and H3K4me1 clearly separated it from both Empty and H2A.Z WT cells (Fig. 4 D). These results suggest that specific chromatin marks such as H3K27ac and H3K27me3 reflect distinct epigenetic states between H2A.Z WT and R80C-expressing cells. Epigenetic signatures at SASP genes are altered by H2A.Z R80C expression Based on this finding, we wondered whether the repression of the SASP genes by H2A.Z R80C might be associated with altered histone modifications at those genes. To test this, we performed CUT&Tag for histone modifications H3K27ac, H3K27me3 and H3K4me1 as well as for HA-tagged H2A.Z WT and R80C in proliferating and senescent cells. Clustering of samples based on genome-wide distribution of HA-tagged H2A.Z wild type and R80C mutant showed greater correlation with distribution of H3K27ac than with H3K27me3, consistent with preferential enrichment of H2A.Z at active gene loci (Fig. 5 A). The enrichment level of H2A.Z R80C did not significantly differ from that of H2A.Z WT across proliferating and senescent cells in H2A.Z WT and R80C expressing cells (Supplementary Fig. 6). By contrast, H3K27ac, an activating histone mark, was significantly enriched at differentially expressed SASP gene loci during senescence, but this was not the case in H2A.Z R80C expressing cells (Fig. 5 B). These findings suggest that H2A.Z R80C may suppress SASP gene expression by attenuating recruitment of activating H3K27ac. Together, these results demonstrate that H3K27ac enrichment at SASP gene loci depends on the stability of H2A.Z. These findings collectively suggest that H2A.Z R80C modulates SASP gene expression by interfering with H3K27ac-associated chromatin activation. DISCUSSION Our results reveal a novel role of the stability of H2A.Z-containing nucleosomes in regulating expression of the SASP during cellular senescence. Expression of the H2A.Z R80C mutant decreased the stability of the nucleosome core particle and specifically downregulated the SASP during senescence, without affecting upregulation of other genes in senescence, such as ISGs, CDKN1A and CDKN2A . The effect of H2A.Z R80C mutation was recapitulated by other substitutions of R80, was specific to senescence over other inducers of cytokine/chemokine expression, and was not due to loss of H2A.Z function. Rather, the effect of H2A.Z R80C appears to be due to impaired marking of SASP genes by H3K27ac. H2A.Z R80C was first identified as a histone mutant that is enriched in cancer. While many histone mutants are seemingly random in nature (likely ‘passenger’ mutations), a small subset of them are recurrent and over-represented within cancer genomes, and are therefore considered potential ‘drivers’ of disease 40 , 44 . Such ‘oncohistone’ mutations are predicted to play a crucial role in affecting nucleosome structure and diverse cellular processes such as chromatin accessibility, DNA damage repair, replication, and histone post-translational modifications - ultimately predisposing cell fate towards transformation. Given the tumor suppressive role of the SASP in senescent cells through its ability to recruit the innate immune system, it is tempting to speculate that the oncogenic effect of H2A.Z R80C is, at least in part, due to suppression of the SASP 45 . Our biochemical analysis of cells ectopically expressing either wild type or mutant H2A.Z proteins suggests that these ectopically expressed histones are incorporated into chromatin and nucleosomes in a manner that is similar to endogenous histones based on imaging and CUT&Tag data. We further observed through bulk RNA-seq that a subset of key SASP genes was repressed when cells ectopically express H2A.Z R80C. This included genes associated with the pro-inflammatory SASP, including IL1A , IL1B , IL6 , CXCL1 , CXCL3 , CXCL5 , and CXCL8 46 . However expression of ISGs which represent a separate arm of the senescence program that is distinct from the SASP as well as cell cycle genes CDKN1A and CDKN2A was not affected 47 , 48 . Notably numerous genes associated with the extracellular matrix and matrix metalloproteinases were downregulated during senescence and were further repressed in H2A.Z R80C-expressing cells, including COL3A1 , COL6A1 , COL6A2 , and PCOLCE . This suggests that H2A.Z R80C selectively modulates senescence-associated genes that influence extracellular environment, rather than globally altering the senescence transcriptional program. We have observed a similar phenotype in cells expressing different substitutions at the R80 site, specifically cysteine (C), aspartic acid (D), glutamic acid (E), lysine (K) and tryptophan (W). These results indicate that the arginine residue at H2A.Z R80 is essential for establishing the SASP phenotype. The residue R80 in histone H2A.Z plays an important role in influencing the structure and stability of the nucleosome due to its close interaction with the adjacent DNA strand that wraps around the nucleosome. Crystal structure modeling of the nucleosome including the H2A.Z histone variant indicates that the positively charged H2A.Z R80 arginine residue closely fits the minor groove of the adjacent negatively charged DNA strand 25 . This suggests that the R80 residue is important in modulating histone-DNA interactions within the assembled nucleosome, and that mutating this residue would potentially impact the kinetics of the formation/disassembly of H2A.Z-containing nucleosomes 41 , 49 . Results reported here indicate that substitution of R80 with any amino-acid is sufficient to disrupt its function. Interestingly, substitution of cysteine for arginine can occur through a simple C to T transition at the first base of the DNA codon, potentially through deamination of methylated C me GT (arginine) to TGT (cysteine) in H2A.Z.1. Together, these results and observations suggest that the frequency of R80C over other substitutions is in part a function of the genetic code, not a specific advantage for cancer of cysteine over substitutions. We further tested if the suppression of the SASP by H2A.Z R80C was specific to the context of cellular senescence or if this mutant would also suppress the same cytokine/chemokines induced by a different trigger. We observed that, upon stimulation of pro-inflammatory genes with IL1A supplementation in media, the activation of genes was comparable between H2A.Z WT and R80C expressing cells, indicating the effect of R80C is specific to the context of senescence 50 . Notably, knockdown of H2A.Z did not recapitulate the suppression of the SASP observed with R80C expression, suggesting that the effect is independent of H2A.Z abundance. This rules out a simple loss-of-function model where the R80C mutation merely leads to functional inactivation of H2A.Z. Instead, our data supports a model where the presence and incorporation of the destabilized H2A.Z R80C protein is required to actively suppress the SASP. Potentially, the mutant protein functions in a dominant-negative or neomorphic fashion. The induction of the SASP is accompanied by redistribution of histone modifications to establish a pro-inflammatory gene expression program. We observed distinct visual epigenetic signatures in senescent cells expressing the R80C mutant using an unbiased, deep-learning approach analyzing immunofluorescent images. Models trained on specific histone modifications—particularly the active enhancer and promoter mark H3K27ac—could successfully discriminate between the two cell states. This initial finding pointed to a defect in the establishment of active chromatin in senescent cells. We confirmed this directly at target gene loci using CUT&Tag. During senescence, H2A.Z WT cells exhibited a significant enrichment of H3K27ac at SASP gene loci, consistent with their robust transcriptional activation. In stark contrast, this gain of H3K27ac was notably reduced in H2A.Z R80C-expressing cells. In conclusion, we have observed a previously unremarked relationship between nucleosome stability and senescence phenotype. Decreasing the stability of H2A.Z containing nucleosomes results in a significant and specific reduction in the expression of pro-inflammatory SASP genes. This finding advances our understanding of cellular senescence, and how it can be modulated. More remains to be explored, especially in terms of the mechanism underlying this phenomenon and its relevance to aging biology. Modulating nucleosome stability may offer a strategy to intervene in senescence and pave the way for new therapeutics in the future that ameliorate the detrimental effects of cellular senescence. METHODS Cell culture IMR90 fibroblasts and HEK293T cells were cultured in growth media DMEM (1x) (10% Fetal Bovine Serum, 2mM L-glutamine, 1% pen-strep) (Gibco cat # 10313021) and incubated at 37 o C, 5% CO2, 3.5% O2 for the former and 37 o C, 5% CO2 for the latter. IMR90s and HEK293Ts were split every 2 to 3 days or when confluent by incubating with 0.25% Trypsin-EDTA (Gibco cat # 25200056) for 5 mins within the cell incubator to detach adherent cells. A minimum of 3x volume of growth media is added to detached cells before cells are spun down at 300g for 3 minutes. Liquid above the cell pellet is then aspirated before the cell pellet is resuspended with growth media with gentle pipetting. Resuspended cells are replated on appropriate surface-treated polystyrene. Frozen cell stocks are made from freezing down growing cells by pelleting cells as described above. Cell pellets are then resuspended in Fetal Bovine Serum with 10% DMSO. Resuspended cells are placed in cryovials in a styrofoam box with the lid closed in -80 o C freezer for at least 6 hours. Cells are transferred to a liquid nitrogen tank afterwards for long term storage. Senescence induction IR induced senescence was carried out by irradiating cells at a dose of 20Gy before an immediate refresh of the appropriate cell media. Irradiation was carried out within an X Ray Irradiator machine (Rad Source Technologies RS 2000 Small Animal Irradiator) that is regularly calibrated. Cells were split 1:3 three days post irradiation, and media of cells are refreshed every 2 to 3 days after becoming senescent and have stopped proliferating. Plasmid preparation DNA sequences for the sequences encoding H2A.Z WT, H2A.Z R80C, H2A.Z R80D, H2A.Z R80E, H2A.K R80K, and H2A.Z R80W were ordered from Integrated DNA Technologies (IDT). pLenti-PGK plasmids used for subsequent lentivirus preparation were obtained from a collaborator and the plasmid sequence checked using Sanger sequencing. Cloning of DNA sequences into plasmid was carried out using the restriction endonucleases BamHI (R3136S) and SalI (R3138S) for digestion of both insert and plasmid sequences while T4 ligase (M0202L) was used for the subsequent ligation reaction. All enzymes used were obtained from New England Biolabs and utilized based on recommended incubation time and temperature. Cloned plasmids were transformed into stable competent C3040 E. coli using heat shock for one minute at 42 o C and then immediately chilling the sample on ice before streaking on sterile Agar plates containing Carbenicillin for selection of transformed E. coli. After incubation of plate for 24 hours at 37 o C single colonies are picked using a sterile pipette tip and grown in Luria-Bertani (LB) broth for 12 hours. LB broth was prepared with LB medium capsules (3002011) from MP Biomedicals and ultrapure water from a Barnstead™ GenPure™ xCAD Plus Ultrapure Water Purification System (Thermo Fisher Scientific, cat # 50136151). Plasmid DNA was subsequently isolated using a ZymoPURE II Plasmid Midiprep Kit (Zymo Research, cat # D4201) and sequence verified using Sanger sequencing. Lentivirus generation and use HEK293T cells were grown to 50–60% confluence and transfected using Lipofectamine 2000 (Thermo Fisher Scientific, cat # 11668030) with 2.4ug/ml of lentiviral expression vector, 1ug/ml psPAX2 packaging plasmid and 0.6ug/ml VSVG plasmid over 6 hrs. Viral laden media was collected after 24 hrs post transfection and aliquoted and stored at -80 C for long term storage. Cells were infected by lentivirus by thawing frozen viral laden media and adding it to fresh growth media of cells after titrating viral MOI. After 24 hours have passed media is refreshed with fresh growth media containing 1ug/ml of puromycin for selection. After 72 hours have passed, cells are split and replated with fresh growth media. Antibodies For western blots and CUT&TAG against GFP we used Abcam anti-GFP antibody (ab6556), while anti-H3K27ac (ab177178) from Abcam and anti-H3K27me3 (9733S) from Cell Signaling Technology was used for CUT&TAG. We further used antibodies against alpha tubulin from Santa Cruz (sc-8035), actin from Sigma-Aldrich (A1978) and H2A.Z from Cell Signaling Technology (2718S) in western blots. For CUT&TAG and immunofluorescence staining against specific epigenetic markers for deep learning based epigenetic analysis, we used the antibodies anti-H3K56ac from Active Motif (39281), anti-H3K4me1 from Abcam (ab8895), anti-H3K27ac from Abcam (ab4729) and from Active Motif (39133), anti-H3K27me3 from Cell Signaling (9733S), anti-H3K4me2 from Sigma-Aldrich (07–030). We used as an IgG control, antibodies from Vector Laboratories (I-1000, I-2000). We used the antibodies against FLAG tag from Cell Signaling Technology (14793S) and against HA tag from Santa Cruz (sc-7392) for western blotting, CUT&TAG, and immunofluorescence staining. Western Blot Proteins are extracted from cells using sample buffer (0.5M TrisHCl ph6.8, Glycerol, SDS, BPB) and ran in denaturing gels. Proteins are blotted on to Immobilon-P PVDF Membrane (MiliporeSigma cat # IPVH00010) using BioRad transblot Turbo (BioRad cat # 1704150) before staining with Ponceau S (Sigma-Aldrich cat # P7170S-1L). Membrane is blocked at room temperature for 30 mins in TBST (Tris-buffered saline, 0.1% Tween 20) with 5% milk and left on a shaker in the dark overnight at 4C with antibody targeting protein of interest. After 3x 10 min washes in TBST, secondary antibodies are diluted to 1:10000 in 5% milk with TBST and incubated with membrane for 1hr at room temperature. Protein bands are imaged using SuperSignal™ West Femto Maximum Sensitivity Substrate (Thermo Fisher Scientific, cat # 34094). An alternative Western Blot protocol using Bullet Blocking buffer (Nacalai USA cat # 13779-01) instead of 5% milk was also employed to achieve similar results. RNA-seq RNA-seq was carried out on the Illumina platform and all analysis was performed in R. Alignment was carried out using pseudo-alignment with Kallisto against hg38 cDNA. R packages used include tidyverse, edgeR, and matrixStats for statistical calculations, as well as plotly, ggrepel, gt, ggplot, gplots, RColorBrewer, and enrichplot, for tables, heatmaps, and plots. GSEA and GO term analysis was carried out in R using the packages GSEABase, msigdbr, and gprofiler2 Immunofluorescence Cells were washed with DPBS (Gibco, cat # 14190144) before fixing with 10% Neutral Buffered Formalin (Thermo Fisher Scientific, cat # 9400-1) for 10 mins at room temperature. Cells were incubated in blocking solution (4% bovine serum albumin, 1% goat serum in DPBS) before antibody staining in blocking solution according to the manufacturer's recommended concentration of antibodies. Secondary antibodies used are Alexa Fluor Goat anti-Rabbit IgG 488 (A11008) and Goat anti-Mouse (A11029) as well as Invitrogen Goat anti-Rabbit IgG 555 (A32732) and Goat anti-Mouse 555 (A32727). DAPI staining was carried out using DAPI solution (100ng/ml DAPI in DPBS). RNA quantification After a DPBS wash step RNA is extracted directly from cells through addition of RNA lysis buffer on to plate and proceeding with steps outlined in Quick-RNA MiniPrep kit (Zymo Research, cat #11–328). RNA is reverse transcribed using ReverAid RT (Thermo Fisher Scientific, cat # K1691) to obtain cDNA which is diluted 1:10 in water and quantified in 384 well format with SYBR™ Green PCR Master Mix (Thermo Fisher Scientific, cat # 4309155) with appropriate primers (full list in Supplementary Table 1). siRNA knockdown siRNA knockdown was performed using a 20nM pool of three different siRNAs targeting H2A.Z using the Trifect RNAi hs.Ri.H2AFZ.13 by IDT. Cells were treated with 3ul RNAiMAX (Thermo Fisher Scientific, cat #13778150) per 6 well volume along with pooled siRNA for 24hr before media was refreshed with fresh growth media. Knockdown is measured at least 48 hrs after treatment of cells with siRNA. CUT&TAG CUT&TAG was carried out according to the original protocol outlined online by Steven Henikoff ( https://www.protocols.io/view/cut-amp-tag-home-bd26i8he ), using native nuclei isolated from proliferating IMR90 cells. 100,000 nuclei were used for each sample. Single amendment was made to use Halt™ Protease Inhibitor Cocktail (100X) (cat # 78430) instead of Roche Complete Protease Inhibitor EDTA-Free tablets. CUT&TAG pAG-Tn5 protein was from Epicypher (SKU: 15-1017). Analysis was carried out using bowtie2 for alignment and SEACR for peak calling. Peak annotation was carried out using annotatePeaks.pl from HOMER. Alignment and annotations were carried out using the latest version of hg38. Heatmaps were generated using deepTools. Epigenetic characterization with deep learning All algorithms and model training were developed in Jupyter with Python 3.10. Pre-trained ResNet models were downloaded using Pytorch and weights were frozen before unfreezing the last layer during training. All code used to carry out cell segmentation, training of ResNet models, and classification is publicly available on GitHub ( https://github.com/zongmingchua/Epigenetic_image_ML ) Declarations DATA AVAILABILITY The sequencing datasets generated and analyzed during the current study are not yet publicly available due to ongoing processing in the Gene Expression Omnibus repository but will be available from the corresponding author on reasonable request or upon publication of this manuscript. AUTHORS CONTRIBUTIONS P.D. A conceived the concept of the study, procured the funding, and was in charge of the overall direction. Z.M and H.T produced the results shown in Fig 1. Z.M produced the results shown in Figs 2-5. A.A assisted in producing plasmids encoding different substitution mutants at the H2A.Z R80 site. A.R and X.L assisted with advice in conducting bioinformatics and machine learning analysis. M.G.T, K.M, A.D, C.C.M, L.H, and N.D assisted with refining experimental protocols and general technical analysis. All authors including were involved in the preparation of the manuscript. FUNDING / FINANCIAL SUPPORT This research was supported by the National Institute on Aging of the National Institutes of Health under Award Number R01AG071464. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. CONFLICT OF INTEREST The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. References Hernandez-Segura, A., Nehme, J. & Demaria, M. Hallmarks of Cellular Senescence. Trends Cell Biol. 28, 436–453 (2018). Rai, T. S. et al. HIRA orchestrates a dynamic chromatin landscape in senescence and is required for suppression of neoplasia. Genes Dev. 28, 2712–2725 (2014). Ivanov, A. et al. Lysosome-mediated processing of chromatin in senescence. J. Cell Biol. 202, 129–143 (2013). Campisi, J. & d’Adda Di Fagagna, F. Cellular senescence: when bad things happen to good cells. Nat. Rev. Mol. Cell Biol. 8, 729–740 (2007). Hayflick, L. THE LIMITED IN VITRO LIFETIME OF HUMAN DIPLOID CELL STRAINS. Exp. Cell Res. 37, 614–636 (1965). Hayflick, L. & Moorhead, P. S. The serial cultivation of human diploid cell strains. Exp. Cell Res. 25, 585–621 (1961). Wang, C. et al. DNA damage response and cellular senescence in tissues of aging mice. Aging Cell 8, 311–323 (2009). Chen, H.-A. et al. Senescence rewires microenvironment sensing to facilitate anti-tumor immunity. Cancer Discov. CD-22-0528 (2022) doi: 10.1158/2159-8290.CD-22-0528 . Reyes, N. S. et al. Sentinel p16INK4a+ cells in the basement membrane form a reparative niche in the lung. Science 378, 192–201 (2022). Ritschka, B. et al. The senescence-associated secretory phenotype induces cellular plasticity and tissue regeneration. Genes Dev. 31, 172–183 (2017). Michaloglou, C. et al. BRAFE600-associated senescence-like cell cycle arrest of human naevi. Nature 436, 720–724 (2005). Xue, W. et al. Senescence and tumour clearance is triggered by p53 restoration in murine liver carcinomas. Nature 445, 656–660 (2007). Coppé, J.-P., Desprez, P.-Y., Krtolica, A. & Campisi, J. The Senescence-Associated Secretory Phenotype: The Dark Side of Tumor Suppression. Annu. Rev. Pathol. Mech. Dis. 5, 99–118 (2010). Hernandez-Gonzalez, F. et al. Human senescent fibroblasts trigger progressive lung fibrosis in mice. Aging 15, 6641–6657 (2023). Ruhland, M. K. et al. Stromal senescence establishes an immunosuppressive microenvironment that drives tumorigenesis. Nat. Commun. 7, 11762 (2016). Baker, D. J. et al. Naturally occurring p16(Ink4a)-positive cells shorten healthy lifespan. Nature 530, 184–189 (2016). Tchkonia, T., Zhu, Y., van Deursen, J., Campisi, J. & Kirkland, J. L. Cellular senescence and the senescent secretory phenotype: therapeutic opportunities. J. Clin. Invest. 123, 966–972 (2013). Yi, S.-J. & Kim, K. New Insights into the Role of Histone Changes in Aging. Int. J. Mol. Sci. 21, 8241 (2020). Paluvai, H., Di Giorgio, E. & Brancolini, C. The Histone Code of Senescence. Cells 9, 466 (2020). Rogakou, E. P. & Sekeri–Pataryas, K. E. Histone variants of H2A and H3 families are regulated during in vitro aging in the same manner as during differentiation. Exp. Gerontol. 34, 741–754 (1999). Tvardovskiy, A., Schwämmle, V., Kempf, S. J., Rogowska-Wrzesinska, A. & Jensen, O. N. Accumulation of histone variant H3.3 with age is associated with profound changes in the histone methylation landscape. Nucleic Acids Res. 45, 9272–9289 (2017). Contrepois, K. et al. Histone variant H2A.J accumulates in senescent cells and promotes inflammatory gene expression. Nat. Commun. 8, 14995 (2017). Duarte, L. F. et al. Histone H3.3 and its proteolytically processed form drive a cellular senescence programme. Nat. Commun. 5, 5210 (2014). Stefanelli, G. et al. Learning and Age-Related Changes in Genome-wide H2A.Z Binding in the Mouse Hippocampus. Cell Rep. 22, 1124–1131 (2018). Suto, R. K., Clarkson, M. J., Tremethick, D. J. & Luger, K. Crystal structure of a nucleosome core particle containing the variant histone H2A.Z. Nat. Struct. Biol. 7, 1121–1124 (2000). Gevry, N., Chan, H. M., Laflamme, L., Livingston, D. M. & Gaudreau, L. p21 transcription is regulated by differential localization of histone H2A.Z. Genes Amp Dev. 21, 1869–1881 (2007). Sales-Gil, R. et al. Non-redundant functions of H2A.Z.1 and H2A.Z.2 in chromosome segregation and cell cycle progression. EMBO Rep. 22, e52061 (2021). Ávila-López, P. A. et al. H2A.Z overexpression suppresses senescence and chemosensitivity in pancreatic ductal adenocarcinoma. Oncogene 40, 2065–2080 (2021). Belotti, E. et al. H2A.Z is dispensable for both basal and activated transcription in post-mitotic mouse muscles. Nucleic Acids Res. 48, 4601–4613 (2020). Buschbeck, M. & Hake, S. B. Variants of core histones and their roles in cell fate decisions, development and cancer. Nat. Rev. Mol. Cell Biol. 18, 299–314 (2017). Luger, K., Dechassa, M. L. & Tremethick, D. J. New insights into nucleosome and chromatin structure: an ordered state or a disordered affair? Nat. Rev. Mol. Cell Biol. 13, 436–447 (2012). Andrews, A. J., Chen, X., Zevin, A., Stargell, L. A. & Luger, K. The Histone Chaperone Nap1 Promotes Nucleosome Assembly by Eliminating Nonnucleosomal Histone DNA Interactions. Mol. Cell 37, 834–842 (2010). Angelov, D. et al. The Histone Variant MacroH2A Interferes with Transcription Factor Binding and SWI/SNF Nucleosome Remodeling. Mol. Cell 11, 1033–1041 (2003). Cui, K. et al. Chromatin signatures in multipotent human hematopoietic stem cells indicate the fate of bivalent genes during differentiation. Cell Stem Cell 4, 80–93 (2009). Valdés-Mora, F. et al. Acetylation of H2A.Z is a key epigenetic modification associated with gene deregulation and epigenetic remodeling in cancer. Genome Res. 22, 307–321 (2012). Hu, G. et al. H2A.Z facilitates access of active and repressive complexes to chromatin in embryonic stem cell self-renewal and differentiation. Cell Stem Cell 12, 180–192 (2013). Horikoshi, N., Arimura, Y., Taguchi, H. & Kurumizaka, H. Crystal structures of heterotypic nucleosomes containing histones H2A.Z and H2A. Open Biol. 6, 160127. Li, S., Wei, T. & Panchenko, A. R. Histone variant H2A.Z modulates nucleosome dynamics to promote DNA accessibility. Nat. Commun. 14, 769 (2023). Jin, C. & Felsenfeld, G. Nucleosome stability mediated by histone variants H3.3 and H2A.Z. Genes Dev. 21, 1519–1529 (2007). Nacev, B. A. et al. The expanding landscape of ‘oncohistone’ mutations in human cancers. Nature 567, 473–478 (2019). Arimura, Y. et al. Cancer-associated mutations of histones H2B, H3.1 and H2A.Z.1 affect the structure and stability of the nucleosome. Nucleic Acids Res. 46, 10007–10018 (2018). Berryman, S., Matthews, K., Lee, J. H., Duffy, S. P. & Ma, H. Image-Based Cell Phenotyping Using Deep Learning . http://biorxiv.org/lookup/doi/10.1101/817544 (2019) doi: 10.1101/817544 . Heckenbach, I. et al. Nuclear morphology is a deep learning biomarker of cellular senescence. Nat. Aging 2, 742–755 (2022). Schwartzentruber, J. et al. Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma. Nature 482, 226–231 (2012). Kang, T.-W. et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature 479, 547–551 (2011). Coppé, J.-P. et al. Senescence-Associated Secretory Phenotypes Reveal Cell-Nonautonomous Functions of Oncogenic RAS and the p53 Tumor Suppressor. PLOS Biol. 6, e301 (2008). Cecco, M. D. et al. L1 drives IFN in senescent cells and promotes age-associated inflammation. Nature 566, 73–78 (2019). Serrano, M., Lin, A. W., McCurrach, M. E., Beach, D. & Lowe, S. W. Oncogenic ras Provokes Premature Cell Senescence Associated with Accumulation of p53 and p16INK4a. Cell 88, 593–602 (1997). Bagert, J. D. et al. Oncohistone mutations enhance chromatin remodeling and alter cell fates. Nat. Chem. Biol. 17, 403–411 (2021). Orjalo, A. V., Bhaumik, D., Gengler, B. K., Scott, G. K. & Campisi, J. Cell surface-bound IL-1α is an upstream regulator of the senescence-associated IL-6/IL-8 cytokine network. Proc. Natl. Acad. Sci. 106, 17031–17036 (2009). Additional Declarations No competing interests reported. Supplementary Files Feb2026R80CManuscriptCZMHTPDASupplementary.docx Feb2026R80CUncroppedgelimageCZMHTPDA.docx Cite Share Download PDF Status: Under Revision Version 1 posted Editorial decision: Revision requested 17 Mar, 2026 Reviews received at journal 16 Mar, 2026 Reviews received at journal 08 Mar, 2026 Reviews received at journal 05 Mar, 2026 Reviewers agreed at journal 03 Mar, 2026 Reviewers agreed at journal 03 Mar, 2026 Reviewers agreed at journal 24 Feb, 2026 Reviewers agreed at journal 23 Feb, 2026 Reviewers agreed at journal 23 Feb, 2026 Reviewers invited by journal 23 Feb, 2026 Editor assigned by journal 23 Feb, 2026 Submission checks completed at journal 19 Feb, 2026 First submitted to journal 07 Feb, 2026 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8818498","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":597191756,"identity":"aba63d5d-65bf-4376-a733-143fa10608bb","order_by":0,"name":"Zong Ming Chua","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Zong","middleName":"Ming","lastName":"Chua","suffix":""},{"id":597191757,"identity":"fadf37a9-5d7e-4fa7-8112-99cd86467e67","order_by":1,"name":"Hiroshi Tanaka","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Hiroshi","middleName":"","lastName":"Tanaka","suffix":""},{"id":597191758,"identity":"76e4f2bc-6206-4f73-a0f6-931875d53810","order_by":2,"name":"Adrianna Abele","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Adrianna","middleName":"","lastName":"Abele","suffix":""},{"id":597191759,"identity":"ab493290-6d12-4686-96c4-6c893727edbc","order_by":3,"name":"Adarsh Rajesh","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Adarsh","middleName":"","lastName":"Rajesh","suffix":""},{"id":597191760,"identity":"768b4f59-49a8-46b0-a420-da27c08fa02d","order_by":4,"name":"Marcos Garcia Teneche","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Marcos","middleName":"Garcia","lastName":"Teneche","suffix":""},{"id":597191761,"identity":"65085fc3-e29d-4f93-b810-f772e17025b4","order_by":5,"name":"Xue Lei","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Xue","middleName":"","lastName":"Lei","suffix":""},{"id":597191762,"identity":"d2aa12d6-da3b-4c2b-994e-0788db613959","order_by":6,"name":"Karl Miller","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Karl","middleName":"","lastName":"Miller","suffix":""},{"id":597191763,"identity":"f35a5183-6b8d-4390-8b99-a1949528ebe7","order_by":7,"name":"Andrew Davis","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Andrew","middleName":"","lastName":"Davis","suffix":""},{"id":597191764,"identity":"35391ab5-53a1-4960-89bd-c0ceb1c8b0b5","order_by":8,"name":"Carolina Cano Macip","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Carolina","middleName":"Cano","lastName":"Macip","suffix":""},{"id":597191765,"identity":"f635613c-f519-432b-9ed4-87cdd095983a","order_by":9,"name":"Laurence Haddadin","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Laurence","middleName":"","lastName":"Haddadin","suffix":""},{"id":597191766,"identity":"f805ecfc-3d98-4cf8-93ee-385ced8a2749","order_by":10,"name":"Nirmalya Dasgupta","email":"","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":false,"prefix":"","firstName":"Nirmalya","middleName":"","lastName":"Dasgupta","suffix":""},{"id":597191767,"identity":"0fd1ebfe-e334-4dc4-9240-009a970d9ad4","order_by":11,"name":"Peter D. Adams","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA70lEQVRIiWNgGAWjYHACAyC2AZEGB0jRkka6lsNgLcSpN2c/vPHTjT/n7c2lmzce/FFjk9gvfcaA4eOeWpxaLHvSiqVz224n7pxzrOAwz7G0xJl9OQaMM54dx+2qAzkG0rkNtxMMbuQYHGZsOGxscIbHgJnnwDHcWs6/Mf6d8+ecPUjLwZ9ALfYEtdzIMZPOYTvAuAGo5QBvw2E5Ax6wlhrcfpnxrMw6ty0Z7hc5iTNsBQdnHDiAU4s5f/Lm2zl/7EAhtvkjMMR4+HuYNz74cKAOt8PgLAk4iwMUp4dJ0sL+AEjgtmUUjIJRMApGHAAAptddNdjjouQAAAAASUVORK5CYII=","orcid":"","institution":"Sanford Burnham Prebys Medical Discovery Institute","correspondingAuthor":true,"prefix":"","firstName":"Peter","middleName":"D.","lastName":"Adams","suffix":""}],"badges":[],"createdAt":"2026-02-08 01:08:16","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8818498/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8818498/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":103473805,"identity":"4b9e4b01-570a-488b-aa97-433860bd1f1d","added_by":"auto","created_at":"2026-02-26 06:27:43","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1097255,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eBiochemical characterization of wild type and mutants of histone variant H2A.Z in cells. \u003c/strong\u003e(A) Location of H2A.Z R80 within the nucleosome core particle (NCP) (B) Representative images of cells co-stained with DAPI, anti-H3K27me3, and anti-HA antibody including quantification of HA signal intensity in cell nuclei (n=4) (C) Western blotting of cells ectopically expressing HA-H2A.Z R80 in proliferating IMR90s using an anti-HA antibody. (D) Representative images of dividing cells in anaphase showing the co-localization of DAPI, HA, and H3K27me3 signal from immunofluorescence. (E) High sensitivity DNA tapestation shows distinct nucleosomal ladder of ectopically expressed HA-H2A.Z histones in IMR90 cells (F) Fluorescence Recovery after Photobleaching (FRAP) was carried out in order to assess kinetics of histone mutants as measured by fluorescence recovery over 90 mins after photobleaching.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/d73d222b4d567b988b47556f.png"},{"id":103508097,"identity":"ec412bc2-6402-4cc5-8a8f-c5c319570905","added_by":"auto","created_at":"2026-02-26 13:47:11","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1399697,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eR80C histone variant H2A.Z specifically suppresses the SASP of senescent cells \u003c/strong\u003e(A) Images of senescent cells expressing H2A.Z WT and R80C mutant histone proteins (B) Heatmap showing the expression dynamics of genes associated with proliferation and cell cycle progression over a 14-day time course following irradiation (C) Heatmap of overlapping differentially expressed genes between H2A.Z WT vs H2A.Z R80C and Empty vs H2A.Z R80C expressing cells 14 days post IR. Gene signature clusters 1-4 are marked. (D) Expression and fold change of selected components of the SASP for senescent cells 14 days post irradiation (E) Volcano plot of H2A.Z WT vs H2A.Z R80C expressing cells at day 14 post irradiation. (F) A panel of interferon stimulated genes show broad upregulation during senescence but no significant changes in expression between H2A.Z WT and R80C mutant expressing cells (G) RT-qPCR of key SASP markers \u003cem\u003eIL1A\u003c/em\u003e, \u003cem\u003eIL6\u003c/em\u003e, and \u003cem\u003eIL8 \u003c/em\u003eshow downregulation in H2A.Z R80C compared to H2A.Z WT expressing senescent cells without affecting expression of \u003cem\u003eCDKN1A \u003c/em\u003e(p21) and \u003cem\u003eCDKN2A \u003c/em\u003e(p16). Statistical significance was assessed using unpaired t-tests; *\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.05, **\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.01, ***\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/159838ce66a13510203e7087.png"},{"id":103473804,"identity":"537007ba-8295-4660-a102-05722caad392","added_by":"auto","created_at":"2026-02-26 06:27:43","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":661622,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMechanistic analysis of SASP suppression by R80C histone variant H2A.Z \u003c/strong\u003e(A) Changing R80 mutant does not notably change reduction of SASP in R80 mutant expressing senescent cells. (B) SASP gene activation at varying concentrations of IL1A compared to untreated cells (NT) in media does not differ significantly between H2A.Z WT and H2A.Z R80C expressing proliferating IMR90 cells. (C) SASP gene activation is not significantly different between H2A.Z WT and H2A.Z R80C expressing senescent IMR90 cells with exposure to 100pg/ml of IL1A for 6 hours in media in comparison to untreated cells (NT). (D) Validation of total H2A.Z and ectopically expressed HA-H2A.Z using RT-qPCR showed a highly efficient knockdown of H2A.Z histones. (E) Knockdown of H2A.Z abrogates reduction of SASP in H2A.Z R80C expressing compared to H2A.Z WT expressing senescent IMR90 cells. Statistical significance was assessed using unpaired t-tests for all panels; *\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.05, **\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.01, ***\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/8a2c69ab36780a5fc660b962.png"},{"id":103473800,"identity":"f294be55-add8-4e57-9184-ca2381fc3288","added_by":"auto","created_at":"2026-02-26 06:27:42","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":587352,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eImage classification using deep learning models to investigate consequences of wild type and mutant histone variant H2A.Z \u003c/strong\u003e(A) Binary Cross Entropy Loss and Accuracy on test data during training of different model architectures after training for 80 to 100 epochs (B) Cells stained with different epigenetic markers can be robustly classified using a trained deep learning model (C) A trained ResNet101 was able to robustly classify proliferating and senescent cells based solely on images of DAPI staining of cell nuclei (D) Cells ectopically expressing H2A.Z R80C show distinct epigenetic patterns that can distinguished by deep learning models.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/b27b5598a6ad85fcff2378ea.png"},{"id":103507907,"identity":"4e88b2e7-eb17-4a27-a054-373f2d7e9d6a","added_by":"auto","created_at":"2026-02-26 13:46:24","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1916841,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eEpigenome distribution of wild type and mutant histone variant H2A.Z \u003c/strong\u003e(A) Heatmap of signal correlation between HA tag and H3K27ac and HA tag and H3K27me3 signal between proliferating and senescent conditions of IMR90 cells expressing both H2A.Z wild type and mutant R80C histones. (B) Enrichment at key SASP genes of HA tag and H3K27ac signal as assessed using CUT\u0026amp;TAG of senescent IMR90 cells expressing wild type and mutant H2A.Z histones at 11 days post irradiation. Statistical significance was assessed using unpaired t-tests; *\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.05, **\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.01, ***\u003cem\u003eP \u003c/em\u003e\u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/f161000f163fc18b7ebdae97.png"},{"id":104397701,"identity":"c7bf2f81-da93-4e8e-bda2-ac6eaab40c11","added_by":"auto","created_at":"2026-03-11 11:54:44","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":7834666,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/10f8be6a-c79f-43c9-a630-f4b922220593.pdf"},{"id":103473806,"identity":"42226d48-c282-42a7-b09d-c7402fe6b522","added_by":"auto","created_at":"2026-02-26 06:27:43","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":5581211,"visible":true,"origin":"","legend":"","description":"","filename":"Feb2026R80CManuscriptCZMHTPDASupplementary.docx","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/e524fb1589c71761f68b0d5a.docx"},{"id":103508005,"identity":"b7274ac4-a84d-46f9-a818-a2f4ec9a71de","added_by":"auto","created_at":"2026-02-26 13:46:50","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":462133,"visible":true,"origin":"","legend":"","description":"","filename":"Feb2026R80CUncroppedgelimageCZMHTPDA.docx","url":"https://assets-eu.researchsquare.com/files/rs-8818498/v1/84e4863929254f287046546c.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Histone variant H2A.Z mutant suppresses the senescence-associated secretory phenotype","fulltext":[{"header":"INTRODUCTION","content":"\u003cp\u003eCellular senescence, a state of irreversible cell cycle arrest triggered by both extrinsic and intrinsic stimuli, is a complex biological process marked by significant alterations in chromatin structure, epigenetic modifications, metabolism and gene expression\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e1\u003c/span\u003e–\u003cspan class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e. Senescent cells remain viable after stable cell cycle arrest, and their burden is known to increase with age in many vertebrates\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e. Senescent cells produce a specific pro-inflammatory senescence-associated secretory phenotype (SASP) which consists of proinflammatory cytokines such as IL-6 and IL-1, chemokines such as IL-8, and extracellular proteases such as MMP-1. This SASP can improve tissue repair and regeneration, anti-tumor immunity, and cellular plasticity\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e8\u003c/span\u003e–\u003cspan class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. However, sustained and chronic expression of the SASP, especially in older organisms, can result in the spread of senescence and impaired tissue function through paracrine signaling\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e. Chronic SASP signaling can also lead to a pro-inflammatory and immunosuppressive microenvironment that drives tumorigenesis\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e. Accordingly, removal of senescent cells has been shown to improve tissue function and promote healthy aging\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e16\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eInduction of cellular senescence and the SASP takes place concurrently with changes in nucleosome composition and localization of histone variants, which in turn regulate key age-related and senescence phenotypes\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e18\u003c/span\u003e–\u003cspan class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. Among the histone variants that play a role in regulating the senescent phenotype is H2A.Z, a variant of core histone H2A, which has been shown to accumulate with age and is linked to age-related cognitive decline in mice\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. H2A.Z has a roughly 60% degree of sequence similarity with canonical H2A, with a higher proportion of residues in the N- and C-terminal tail domain showing differences compared to the more conserved core regions\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eIn some contexts, H2A.Z has a pro-proliferative, anti-senescence role. The p53-dependent localization of the histone variant H2A.Z by the p400 complex represses basal p21 transcription in proliferating cells, and upon DNA damage, H2A.Z is evicted to permit Tip60 recruitment and p21-mediated proliferation arrest\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. Additionally, H2A.Z supports the expression of proliferation-associated genes such as Myc and Ki-67, further reinforcing its role in cell cycle progression and senescence avoidance\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. In pancreatic ductal adenocarcinoma (PDAC), overexpression of H2A.Z enhances tumor growth and confers chemoresistance by suppressing senescence, again highlighting a pro-proliferative function\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. On the other hand, H2A.Z is expressed throughout the cell cycle and in non-proliferating cells, unlike canonical H2A which is predominantly expressed during S phase\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Accordingly, it has been shown that the growth-arrested senescent phenotype is associated with the site-specific deposition of H2A.Z, suggesting a role for H2A.Z in non-proliferating senescent cells.\u003c/p\u003e \u003cp\u003eThe histone variant H2A.Z directly impacts the structure and the stability of assembled nucleosomes, altering the activity of chromatin remodelers and the accessibility of the underlying DNA, thereby affecting diverse DNA-related processes such as repair, recombination, and transcription\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e30\u003c/span\u003e–\u003cspan class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. H2A.Z is localized at the transcription start site (TSS) of actively transcribed genes, and its presence plays a role in regulating the activity of poised genes which are marked by both activating and repressive histone marks\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e34\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. Homotypic nucleosomes containing pairs of H2A.Z histones show decreased nucleosomal stability, while H2A.Z containing nucleosomes show altered dynamics, mobility, and increased accessibility of the local chromatin to transcriptional machinery\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e36\u003c/span\u003e–\u003cspan class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eGiven both the importance of cellular senescence and the accompanying SASP in aging, and the profound changes in transcriptome that take place during the onset of cellular senescence and the proposed role of H2A.Z in control of gene expression through nucleosome stability, we sought to investigate how modulating the stability of H2A.Z containing nucleosomes would impact the SASP, potentially providing new avenues to ameliorate age-associated disease phenotypes.\u003c/p\u003e \u003cp\u003eTo understand how the stability of H2A.Z-containing nucleosomes plays a role in the induction of cellular senescence, we examined the effects of a specific histone mutant H2A.Z R80C in senescent cells. H2A.Z R80C was previously identified as a histone mutant recurrently found in some cancers (a so-called oncohistone mutation) and has been shown to decrease the stability of the core nucleosomal complex\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e40\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. Modeling of the H2A.Z R80C mutant shows that removal of the arginine 80 residue reduces the affinity between the H2A.Z histone and the adjacent DNA strand, as the arginine interacts favorably with the negatively charged DNA backbone while also extending into the minor groove of the DNA strand to stabilize the nucleosome\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. We used H2A.Z R80C as a tool to probe the relationship between nucleosome stability and phenotypes of cell senescence.\u003c/p\u003e"},{"header":"RESULTS","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eH2A.Z R80C alters nucleosome stability\u003c/h2\u003e \u003cp\u003eTo investigate the relationship between H2A.Z stability and phenotypes of senescence, we employed the H2A.Z R80C mutant that subtly destabilizes nucleosomes\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e. The H2A.Z R80C substitution mutation replaces an arginine residue with cysteine, which is modeled to destabilize the nucleosome due to the loss of the favorable arginineDNA interaction in an assembled nucleosome (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eA). We carried out biochemical characterization of lentivirus-mediated ectopically expressed H2A.Z wild type (WT) or mutant H2A.Z R80C, both with an N-terminal HA-tag or GFP-tag, in IMR-90 primary human fibroblasts. We confirmed that ectopic H2A.Z WT and R80C were expressed at comparable levels, based on western blotting and quantification of the immunofluorescent signals (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eB, \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). More than 99% of cells showed clear ectopic expression of H2A.Z WT or R80C.\u003c/p\u003e \u003cp\u003eSeveral lines of evidence indicated that both H2A.Z WT and R80C are incorporated into nucleosomal chromatin. First, the fluorescent signal of the HA-tag co-localized with DAPI and H3K27me3 in both interphase and mitotic chromatin (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eB, \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eD). Co-localization of the WT and mutant histones on mitotic chromosomes strongly suggests incorporation of both variants into chromatin. A DNA digestion assay targeted by anti-HA antibodies using Tn5 transposase in IMR90 cells expressing either HA-tagged H2A.Z WT or R80C histones revealed clear nucleosomal ladder patterns, indicating that the ectopically expressed histones are both incorporated into chromatin (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eE).\u003c/p\u003e \u003cp\u003eTo assess the protein mobility of H2A.Z WT and R80C in the nucleus, we carried out fluorescence recovery after photobleaching (FRAP) to measure the kinetics of nuclear diffusion of GFP-tagged H2A.Z WT and R80C. The H2A.Z R80C mutant showed a more rapid fluorescence recovery rate compared to H2A.Z WT, indicating that H2A.Z R80C is more mobile within the nucleus, consistent with decreased nucleosome stability\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e1\u003c/span\u003eF). In sum, H2A.Z WT and R80C are comparably incorporated into nucleosomes and chromatin at steady state, but R80C is more mobile consistent with decreased nucleosome stability.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eH2A.Z R80C mutant suppresses induction of a subset of the SASP\u003c/h3\u003e\n\u003cp\u003eHaving characterized H2A.Z WT and R80C mutant in proliferating IMR-90 cells, we next investigated their impact on senescence-associated phenotypes. Upon irradiation, H2A.Z WT and H2A.Z R80C exhibited a comparable large and flattened morphology characteristic of senescence within 14 days (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eA). Consistent with the induction of senescence, for both H2A.Z WT and R80C cells we observed a clear downregulation of genes associated with cellular proliferation and cell cycle progression following irradiation (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eB). A principal component (PC) analysis of mRNA-seq at six time points, 0, 5, 8, 10, 12, and 14 days post-irradiation, revealed that roughly a quarter of transcriptomic variance (25.9%) was explained by PC1 with a clear separation of samples by time post-irradiation (Supplementary Fig.\u0026nbsp;1). H2A.Z R80C tended to separate from empty and H2A.Z WT along PC2.\u003c/p\u003e \u003cp\u003eAn unbiased clustering analysis of differentially expressed genes at day 14 post-irradiation identified four distinct gene clusters when comparing empty and H2A.Z WT controls and R80C-expressing cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eC). Among these clusters, cluster 1 represented genes with expression levels that were unaffected during senescence in both empty and H2A.Z WT cells, but were upregulated in H2A.Z R80C cells. Cluster 2 included genes that were upregulated during senescence in empty, H2A.Z WT and H2A.Z R80C cells, albeit more so and less so in R80C in some cases. Cluster 3 contained genes downregulated during senescence in empty and H2A.Z WT cells, with further decrease in H2A.Z R80C. Finally, cluster 4 represented genes markedly upregulated during senescence in empty and H2A.Z WT cells, but less so in H2A.Z R80C cells.\u003c/p\u003e \u003cp\u003eFunctional enrichment analysis of each gene cluster revealed distinct biological associations: cluster 1 was enriched with genes related to glycerolipid and cholesterol metabolism; cluster 2 with genes involved in cytokine and inflammatory responses; cluster 3 with genes related to extracellular matrix organization; and cluster 4 with genes encoding cell surface proteins (Supplementary Fig.\u0026nbsp;2). Further analysis of these differentially expressed genes using the STRING database supports the results of this functional enrichment analysis in identifying both a group of pro-inflammatory genes and a cluster of collagen and matrix metalloproteinase genes that are downregulated in H2A.Z R80C expressing cells (Supplementary Fig.\u0026nbsp;3).\u003c/p\u003e \u003cp\u003eWe validated these functional enrichment analyses with a gene set enrichment analysis (GSEA) on the same set of differentially expressed genes and observed that expression of genes in pathways related to inflammation was repressed in H2A.Z R80C samples at day 14 post irradiation compared to H2A.Z WT (Supplementary Fig.\u0026nbsp;4). Repression of genes related to inflammation was also apparent in parts of cluster 2 (IL1A, IL1B, IL6, CXCL3, CXCL5) and cluster 4 (CXCL1, CXCL8) (Supplementary Fig.\u0026nbsp;2). These functional enrichment results indicate that there are two dominant repressive gene expression changes associated with expression of the H2A.Z R80C mutant during senescence, with one centered on the repression of a pro-inflammatory gene signature and the other on the repression of genes associated with the extracellular matrix.\u003c/p\u003e \u003cp\u003eGiven that the SASP is enriched in both inflammatory cytokines/chemokines and extracellular matrix components, we examined whether the gene signatures from our unbiased analysis reflect the known SASP signature\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. We examined the expression of several key SASP genes and found that their expression in senescence was lower in H2A.Z R80C expressing cells compared to H2A.Z WT (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eD). Expression of these key SASP genes was significantly downregulated in senescent H2A.Z R80C expressing cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eE). Importantly, the expression of H2A.Z R80C did not affect the expression of canonical interferon-stimulated genes (ISGs) that were upregulated during senescence (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eF), confirming that H2A.Z R80C specifically affects a subset of the SASP. We verified suppression of key SASP markers \u003cem\u003eIL1A\u003c/em\u003e, \u003cem\u003eIL6\u003c/em\u003e, and \u003cem\u003eIL8\u003c/em\u003e/\u003cem\u003eCXCL8\u003c/em\u003e in H2A.Z R80C expressing cells by qPCR, while the expression of cell cycle genes \u003cem\u003eCDKN2A\u003c/em\u003e and \u003cem\u003eCDKN1A\u003c/em\u003e was not significantly impaired (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e2\u003c/span\u003eG). In summary, these results indicate that expression of H2A.Z R80C specifically suppresses genes associated with the SASP, including chemokines and cytokines involved in the proinflammatory phenotype of senescence.\u003c/p\u003e \u003cp\u003e \u003cem\u003eSASP suppression by H2A.Z R80C is specific to senescence and independent of loss of H2A.Z function\u003c/em\u003e \u003c/p\u003e \u003cp\u003eNext, we set out to investigate the mechanism by which H2A.Z R80C suppresses the SASP. First, we asked whether suppression of SASP is specific to the R80C substitution or can also be achieved by other amino-acid substitutions. To do this, we tested various H2A.Z R80 substitutions for their effect on senescence, including bulky residues (R80W), and those with altered (R80D, R80E) or similar (R80K) charge properties to the original arginine residue. Expression of all mutants was confirmed by RT-qPCR (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e3\u003c/span\u003eA). We observed similar reductions in SASP gene expression during senescence across all mutations (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e3\u003c/span\u003eB), indicating that SASP suppression is agnostic to the specific identity of the residue in the substitution. This is consistent with structural studies showing that R80 makes specific contact with DNA, and suggests that this interaction is disrupted by any mutation\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eSecond, we asked whether H2A.Z R80C blocks expression of the canonical SASP genes only when they were induced as part of the senescence program. To investigate this, we added recombinant IL1A, a pro-inflammatory cytokine, into the media of senescent IMR90 cells and observed expression of SASP components \u003cem\u003eIL1A\u003c/em\u003e, \u003cem\u003eIL6\u003c/em\u003e and \u003cem\u003eIL8\u003c/em\u003e by RT-qPCR. Under these conditions, we observed the reduction in \u003cem\u003eIL1A\u003c/em\u003e, \u003cem\u003eIL6\u003c/em\u003e and \u003cem\u003eIL8\u003c/em\u003e expression by H2A.Z R80C was largely absent (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e3\u003c/span\u003eC). This indicates that H2A.Z R80C does not affect the capacity of cells to express \u003cem\u003eIL6\u003c/em\u003e and \u003cem\u003eIL8\u003c/em\u003e in response to other non-senescence triggers. Its suppressive effect on these SASP components is specific to the senescence program.\u003c/p\u003e \u003cp\u003eFinally, we wanted to know whether the suppression of the SASP by H2A.Z R80C is due to a loss of H2A.Z function. To test this, we asked whether suppression of the SASP was also achieved by H2A.Z knockdown. We performed siRNA knockdown of H2A.Z in senescent cells using a pool of three different siRNAs targeting H2A.Z.1. H2A.Z.2 (H2AFV) was not included in this analysis due to its substantially lower expression levels relative to H2A.Z.1 in our system, making it less likely to contribute significantly to the observed effects (Supplementary Fig.\u0026nbsp;5). Using RT-qPCR, we verified knockdown efficiency of ~\u0026thinsp;98% based on transcript level of both total H2A.Z as well as ectopically expressed HA-tagged H2A.Z (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e3\u003c/span\u003eD). As expected, knockdown of ectopic H2A.Z R80C abrogated the reduction of the SASP by H2A.Z R80C, confirming that H2A.Z R80C in senescent cells represses the SASP (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e3\u003c/span\u003eE). However, knockdown of H2A.Z did not suppress the SASP gene expression in both Empty and H2A.Z WT cells. Indeed, knockdown of H2A.Z increased expression of SASP genes in several cases, suggesting a transcription repressive role in some cases. Regardless, these results indicate that the suppression of the SASP during senescence by H2A.Z R80C is not due to a loss of H2A.Z function.\u003c/p\u003e\n\u003ch3\u003eH2A.Z R80C shows altered epigenetic signatures in image-based deep learning classification\u003c/h3\u003e\n\u003cp\u003eGiven the role of histone modifications in control of gene expression, we hypothesized that the altered transcriptomic program of H2A.Z R80C-expressing cells is linked to changes in epigenetic signatures. We employed deep learning techniques to examine differences in various epigenetic signatures in proliferating and senescent IMR-90 cells using immunofluorescent images of key regulatory histone modifications: H3K27ac, H3K27me3, H3K4me1, and H3K56ac\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e,\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e. Comparisons among different neural networks revealed that the ResNet50 pre-trained convolutional neural networks performed best to differentiate between proliferating cells based on immunofluorescence staining of H3K27ac, H3K27me3, H3K4me1, H3K56ac and DAPI staining (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). The trained ResNet50 model robustly distinguished various nuclear epigenetic marks, while highlighting a visually similar signature between the enhancer marks H3K27ac and H3K4me1 compared to other modifications, demonstrating the power of this machine learning approach to detect biological differences in epigenetic signatures between individual cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). We observed that, based solely on DAPI staining, a separate trained ResNet model was able to robustly distinguish between proliferating and senescent cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e4\u003c/span\u003eC). While DAPI alone could not clearly distinguish R80C cells from others, H3K27ac, H3K27me3, and H3K4me1 clearly separated it from both Empty and H2A.Z WT cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e4\u003c/span\u003eD). These results suggest that specific chromatin marks such as H3K27ac and H3K27me3 reflect distinct epigenetic states between H2A.Z WT and R80C-expressing cells.\u003c/p\u003e\n\u003ch3\u003eEpigenetic signatures at SASP genes are altered by H2A.Z R80C expression\u003c/h3\u003e\n\u003cp\u003eBased on this finding, we wondered whether the repression of the SASP genes by H2A.Z R80C might be associated with altered histone modifications at those genes. To test this, we performed CUT\u0026amp;Tag for histone modifications H3K27ac, H3K27me3 and H3K4me1 as well as for HA-tagged H2A.Z WT and R80C in proliferating and senescent cells. Clustering of samples based on genome-wide distribution of HA-tagged H2A.Z wild type and R80C mutant showed greater correlation with distribution of H3K27ac than with H3K27me3, consistent with preferential enrichment of H2A.Z at active gene loci (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e5\u003c/span\u003eA). The enrichment level of H2A.Z R80C did not significantly differ from that of H2A.Z WT across proliferating and senescent cells in H2A.Z WT and R80C expressing cells (Supplementary Fig.\u0026nbsp;6). By contrast, H3K27ac, an activating histone mark, was significantly enriched at differentially expressed SASP gene loci during senescence, but this was not the case in H2A.Z R80C expressing cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e5\u003c/span\u003eB). These findings suggest that H2A.Z R80C may suppress SASP gene expression by attenuating recruitment of activating H3K27ac. Together, these results demonstrate that H3K27ac enrichment at SASP gene loci depends on the stability of H2A.Z. These findings collectively suggest that H2A.Z R80C modulates SASP gene expression by interfering with H3K27ac-associated chromatin activation.\u003c/p\u003e"},{"header":"DISCUSSION","content":"\u003cp\u003eOur results reveal a novel role of the stability of H2A.Z-containing nucleosomes in regulating expression of the SASP during cellular senescence. Expression of the H2A.Z R80C mutant decreased the stability of the nucleosome core particle and specifically downregulated the SASP during senescence, without affecting upregulation of other genes in senescence, such as ISGs, \u003cem\u003eCDKN1A\u003c/em\u003e and \u003cem\u003eCDKN2A\u003c/em\u003e. The effect of H2A.Z R80C mutation was recapitulated by other substitutions of R80, was specific to senescence over other inducers of cytokine/chemokine expression, and was not due to loss of H2A.Z function. Rather, the effect of H2A.Z R80C appears to be due to impaired marking of SASP genes by H3K27ac.\u003c/p\u003e \u003cp\u003eH2A.Z R80C was first identified as a histone mutant that is enriched in cancer. While many histone mutants are seemingly random in nature (likely ‘passenger’ mutations), a small subset of them are recurrent and over-represented within cancer genomes, and are therefore considered potential ‘drivers’ of disease\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e40\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e. Such ‘oncohistone’ mutations are predicted to play a crucial role in affecting nucleosome structure and diverse cellular processes such as chromatin accessibility, DNA damage repair, replication, and histone post-translational modifications - ultimately predisposing cell fate towards transformation. Given the tumor suppressive role of the SASP in senescent cells through its ability to recruit the innate immune system, it is tempting to speculate that the oncogenic effect of H2A.Z R80C is, at least in part, due to suppression of the SASP\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eOur biochemical analysis of cells ectopically expressing either wild type or mutant H2A.Z proteins suggests that these ectopically expressed histones are incorporated into chromatin and nucleosomes in a manner that is similar to endogenous histones based on imaging and CUT\u0026amp;Tag data. We further observed through bulk RNA-seq that a subset of key SASP genes was repressed when cells ectopically express H2A.Z R80C. This included genes associated with the pro-inflammatory SASP, including \u003cem\u003eIL1A\u003c/em\u003e, \u003cem\u003eIL1B\u003c/em\u003e, \u003cem\u003eIL6\u003c/em\u003e, \u003cem\u003eCXCL1\u003c/em\u003e, \u003cem\u003eCXCL3\u003c/em\u003e, \u003cem\u003eCXCL5\u003c/em\u003e, and \u003cem\u003eCXCL8\u003c/em\u003e\u003csup\u003e46\u003c/sup\u003e. However expression of ISGs which represent a separate arm of the senescence program that is distinct from the SASP as well as cell cycle genes \u003cem\u003eCDKN1A\u003c/em\u003e and \u003cem\u003eCDKN2A\u003c/em\u003e was not affected\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e47\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e. Notably numerous genes associated with the extracellular matrix and matrix metalloproteinases were downregulated during senescence and were further repressed in H2A.Z R80C-expressing cells, including \u003cem\u003eCOL3A1\u003c/em\u003e, \u003cem\u003eCOL6A1\u003c/em\u003e, \u003cem\u003eCOL6A2\u003c/em\u003e, and \u003cem\u003ePCOLCE\u003c/em\u003e. This suggests that H2A.Z R80C selectively modulates senescence-associated genes that influence extracellular environment, rather than globally altering the senescence transcriptional program.\u003c/p\u003e \u003cp\u003eWe have observed a similar phenotype in cells expressing different substitutions at the R80 site, specifically cysteine (C), aspartic acid (D), glutamic acid (E), lysine (K) and tryptophan (W). These results indicate that the arginine residue at H2A.Z R80 is essential for establishing the SASP phenotype. The residue R80 in histone H2A.Z plays an important role in influencing the structure and stability of the nucleosome due to its close interaction with the adjacent DNA strand that wraps around the nucleosome. Crystal structure modeling of the nucleosome including the H2A.Z histone variant indicates that the positively charged H2A.Z R80 arginine residue closely fits the minor groove of the adjacent negatively charged DNA strand\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e. This suggests that the R80 residue is important in modulating histone-DNA interactions within the assembled nucleosome, and that mutating this residue would potentially impact the kinetics of the formation/disassembly of H2A.Z-containing nucleosomes\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e41\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. Results reported here indicate that substitution of R80 with any amino-acid is sufficient to disrupt its function. Interestingly, substitution of cysteine for arginine can occur through a simple C to T transition at the first base of the DNA codon, potentially through deamination of methylated C\u003csup\u003eme\u003c/sup\u003eGT (arginine) to TGT (cysteine) in H2A.Z.1. Together, these results and observations suggest that the frequency of R80C over other substitutions is in part a function of the genetic code, not a specific advantage for cancer of cysteine over substitutions.\u003c/p\u003e \u003cp\u003eWe further tested if the suppression of the SASP by H2A.Z R80C was specific to the context of cellular senescence or if this mutant would also suppress the same cytokine/chemokines induced by a different trigger. We observed that, upon stimulation of pro-inflammatory genes with IL1A supplementation in media, the activation of genes was comparable between H2A.Z WT and R80C expressing cells, indicating the effect of R80C is specific to the context of senescence\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eNotably, knockdown of H2A.Z did not recapitulate the suppression of the SASP observed with R80C expression, suggesting that the effect is independent of H2A.Z abundance. This rules out a simple loss-of-function model where the R80C mutation merely leads to functional inactivation of H2A.Z. Instead, our data supports a model where the presence and incorporation of the destabilized H2A.Z R80C protein is required to actively suppress the SASP. Potentially, the mutant protein functions in a dominant-negative or neomorphic fashion.\u003c/p\u003e \u003cp\u003eThe induction of the SASP is accompanied by redistribution of histone modifications to establish a pro-inflammatory gene expression program. We observed distinct visual epigenetic signatures in senescent cells expressing the R80C mutant using an unbiased, deep-learning approach analyzing immunofluorescent images. Models trained on specific histone modifications—particularly the active enhancer and promoter mark H3K27ac—could successfully discriminate between the two cell states. This initial finding pointed to a defect in the establishment of active chromatin in senescent cells. We confirmed this directly at target gene loci using CUT\u0026amp;Tag. During senescence, H2A.Z WT cells exhibited a significant enrichment of H3K27ac at SASP gene loci, consistent with their robust transcriptional activation. In stark contrast, this gain of H3K27ac was notably reduced in H2A.Z R80C-expressing cells.\u003c/p\u003e \u003cp\u003eIn conclusion, we have observed a previously unremarked relationship between nucleosome stability and senescence phenotype. Decreasing the stability of H2A.Z containing nucleosomes results in a significant and specific reduction in the expression of pro-inflammatory SASP genes. This finding advances our understanding of cellular senescence, and how it can be modulated. More remains to be explored, especially in terms of the mechanism underlying this phenomenon and its relevance to aging biology. Modulating nucleosome stability may offer a strategy to intervene in senescence and pave the way for new therapeutics in the future that ameliorate the detrimental effects of cellular senescence.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003c/div\u003e \u003c/div\u003e\n\n "},{"header":"METHODS","content":"\u003ch2\u003eCell culture\u003c/h2\u003e\u003cp\u003eIMR90 fibroblasts and HEK293T cells were cultured in growth media DMEM (1x) (10% Fetal Bovine Serum, 2mM L-glutamine, 1% pen-strep) (Gibco cat # 10313021) and incubated at 37\u003csup\u003eo\u003c/sup\u003eC, 5% CO2, 3.5% O2 for the former and 37\u003csup\u003eo\u003c/sup\u003eC, 5% CO2 for the latter. IMR90s and HEK293Ts were split every 2 to 3 days or when confluent by incubating with 0.25% Trypsin-EDTA (Gibco cat # 25200056) for 5 mins within the cell incubator to detach adherent cells. A minimum of 3x volume of growth media is added to detached cells before cells are spun down at 300g for 3 minutes. Liquid above the cell pellet is then aspirated before the cell pellet is resuspended with growth media with gentle pipetting. Resuspended cells are replated on appropriate surface-treated polystyrene.\u003c/p\u003e\u003cp\u003eFrozen cell stocks are made from freezing down growing cells by pelleting cells as described above. Cell pellets are then resuspended in Fetal Bovine Serum with 10% DMSO. Resuspended cells are placed in cryovials in a styrofoam box with the lid closed in -80\u003csup\u003eo\u003c/sup\u003eC freezer for at least 6 hours. Cells are transferred to a liquid nitrogen tank afterwards for long term storage.\u003c/p\u003e\u003ch3\u003eSenescence induction\u003c/h3\u003e\u003cp\u003eIR induced senescence was carried out by irradiating cells at a dose of 20Gy before an immediate refresh of the appropriate cell media. Irradiation was carried out within an X Ray Irradiator machine (Rad Source Technologies RS 2000 Small Animal Irradiator) that is regularly calibrated. Cells were split 1:3 three days post irradiation, and media of cells are refreshed every 2 to 3 days after becoming senescent and have stopped proliferating.\u003c/p\u003e\u003ch2\u003ePlasmid preparation\u003c/h2\u003e\u003cp\u003eDNA sequences for the sequences encoding H2A.Z WT, H2A.Z R80C, H2A.Z R80D, H2A.Z R80E, H2A.K R80K, and H2A.Z R80W were ordered from Integrated DNA Technologies (IDT). pLenti-PGK plasmids used for subsequent lentivirus preparation were obtained from a collaborator and the plasmid sequence checked using Sanger sequencing.\u003c/p\u003e\u003cp\u003eCloning of DNA sequences into plasmid was carried out using the restriction endonucleases BamHI (R3136S) and SalI (R3138S) for digestion of both insert and plasmid sequences while T4 ligase (M0202L) was used for the subsequent ligation reaction. All enzymes used were obtained from New England Biolabs and utilized based on recommended incubation time and temperature.\u003c/p\u003e\u003cp\u003eCloned plasmids were transformed into stable competent C3040 E. coli using heat shock for one minute at 42\u003csup\u003eo\u003c/sup\u003eC and then immediately chilling the sample on ice before streaking on sterile Agar plates containing Carbenicillin for selection of transformed E. coli. After incubation of plate for 24 hours at 37\u003csup\u003eo\u003c/sup\u003eC single colonies are picked using a sterile pipette tip and grown in Luria-Bertani (LB) broth for 12 hours. LB broth was prepared with LB medium capsules (3002011) from MP Biomedicals and ultrapure water from a Barnstead™ GenPure™ xCAD Plus Ultrapure Water Purification System (Thermo Fisher Scientific, cat # 50136151). Plasmid DNA was subsequently isolated using a ZymoPURE II Plasmid Midiprep Kit (Zymo Research, cat # D4201) and sequence verified using Sanger sequencing.\u003c/p\u003e\u003ch2\u003eLentivirus generation and use\u003c/h2\u003e\u003cp\u003eHEK293T cells were grown to 50–60% confluence and transfected using Lipofectamine 2000 (Thermo Fisher Scientific, cat # 11668030) with 2.4ug/ml of lentiviral expression vector, 1ug/ml psPAX2 packaging plasmid and 0.6ug/ml VSVG plasmid over 6 hrs. Viral laden media was collected after 24 hrs post transfection and aliquoted and stored at -80 C for long term storage.\u003c/p\u003e\u003cp\u003eCells were infected by lentivirus by thawing frozen viral laden media and adding it to fresh growth media of cells after titrating viral MOI. After 24 hours have passed media is refreshed with fresh growth media containing 1ug/ml of puromycin for selection. After 72 hours have passed, cells are split and replated with fresh growth media.\u003c/p\u003e\u003ch2\u003eAntibodies\u003c/h2\u003e\u003cp\u003eFor western blots and CUT\u0026amp;TAG against GFP we used Abcam anti-GFP antibody (ab6556), while anti-H3K27ac (ab177178) from Abcam and anti-H3K27me3 (9733S) from Cell Signaling Technology was used for CUT\u0026amp;TAG. We further used antibodies against alpha tubulin from Santa Cruz (sc-8035), actin from Sigma-Aldrich (A1978) and H2A.Z from Cell Signaling Technology (2718S) in western blots.\u003c/p\u003e\u003cp\u003eFor CUT\u0026amp;TAG and immunofluorescence staining against specific epigenetic markers for deep learning based epigenetic analysis, we used the antibodies anti-H3K56ac from Active Motif (39281), anti-H3K4me1 from Abcam (ab8895), anti-H3K27ac from Abcam (ab4729) and from Active Motif (39133), anti-H3K27me3 from Cell Signaling (9733S), anti-H3K4me2 from Sigma-Aldrich (07–030). We used as an IgG control, antibodies from Vector Laboratories (I-1000, I-2000).\u003c/p\u003e\u003cp\u003eWe used the antibodies against FLAG tag from Cell Signaling Technology (14793S) and against HA tag from Santa Cruz (sc-7392) for western blotting, CUT\u0026amp;TAG, and immunofluorescence staining.\u003c/p\u003e\u003ch2\u003eWestern Blot\u003c/h2\u003e\u003cp\u003eProteins are extracted from cells using sample buffer (0.5M TrisHCl ph6.8, Glycerol, SDS, BPB) and ran in denaturing gels. Proteins are blotted on to Immobilon-P PVDF Membrane (MiliporeSigma cat # IPVH00010) using BioRad transblot Turbo (BioRad cat # 1704150) before staining with Ponceau S (Sigma-Aldrich cat # P7170S-1L). Membrane is blocked at room temperature for 30 mins in TBST (Tris-buffered saline, 0.1% Tween 20) with 5% milk and left on a shaker in the dark overnight at 4C with antibody targeting protein of interest. After 3x 10 min washes in TBST, secondary antibodies are diluted to 1:10000 in 5% milk with TBST and incubated with membrane for 1hr at room temperature. Protein bands are imaged using SuperSignal™ West Femto Maximum Sensitivity Substrate (Thermo Fisher Scientific, cat # 34094).\u003c/p\u003e\u003cp\u003eAn alternative Western Blot protocol using Bullet Blocking buffer (Nacalai USA cat # 13779-01) instead of 5% milk was also employed to achieve similar results.\u003c/p\u003e\u003ch2\u003eRNA-seq\u003c/h2\u003e\u003cp\u003eRNA-seq was carried out on the Illumina platform and all analysis was performed in R. Alignment was carried out using pseudo-alignment with Kallisto against hg38 cDNA. R packages used include tidyverse, edgeR, and matrixStats for statistical calculations, as well as plotly, ggrepel, gt, ggplot, gplots, RColorBrewer, and enrichplot, for tables, heatmaps, and plots. GSEA and GO term analysis was carried out in R using the packages GSEABase, msigdbr, and gprofiler2\u003c/p\u003e\u003ch2\u003eImmunofluorescence\u003c/h2\u003e\u003cp\u003eCells were washed with DPBS (Gibco, cat # 14190144) before fixing with 10% Neutral Buffered Formalin (Thermo Fisher Scientific, cat # 9400-1) for 10 mins at room temperature. Cells were incubated in blocking solution (4% bovine serum albumin, 1% goat serum in DPBS) before antibody staining in blocking solution according to the manufacturer's recommended concentration of antibodies. Secondary antibodies used are Alexa Fluor Goat anti-Rabbit IgG 488 (A11008) and Goat anti-Mouse (A11029) as well as Invitrogen Goat anti-Rabbit IgG 555 (A32732) and Goat anti-Mouse 555 (A32727). DAPI staining was carried out using DAPI solution (100ng/ml DAPI in DPBS).\u003c/p\u003e\u003ch2\u003eRNA quantification\u003c/h2\u003e\u003cp\u003eAfter a DPBS wash step RNA is extracted directly from cells through addition of RNA lysis buffer on to plate and proceeding with steps outlined in Quick-RNA MiniPrep kit (Zymo Research, cat #11–328). RNA is reverse transcribed using ReverAid RT (Thermo Fisher Scientific, cat # K1691) to obtain cDNA which is diluted 1:10 in water and quantified in 384 well format with SYBR™ Green PCR Master Mix (Thermo Fisher Scientific, cat #\u003c/p\u003e\u003cp\u003e4309155) with appropriate primers (full list in Supplementary Table\u0026nbsp;1).\u003c/p\u003e\u003ch2\u003esiRNA knockdown\u003c/h2\u003e\u003cp\u003esiRNA knockdown was performed using a 20nM pool of three different siRNAs targeting H2A.Z using the Trifect RNAi hs.Ri.H2AFZ.13 by IDT. Cells were treated with 3ul RNAiMAX (Thermo Fisher Scientific, cat #13778150) per 6 well volume along with pooled siRNA for 24hr before media was refreshed with fresh growth media. Knockdown is measured at least 48 hrs after treatment of cells with siRNA.\u003c/p\u003e\u003ch2\u003eCUT\u0026amp;TAG\u003c/h2\u003e\u003cp\u003eCUT\u0026amp;TAG was carried out according to the original protocol outlined online by Steven Henikoff (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.protocols.io/view/cut-amp-tag-home-bd26i8he\u003c/span\u003e\u003cspan class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), using native nuclei isolated from proliferating IMR90 cells. 100,000 nuclei were used for each sample. Single amendment was made to use Halt™ Protease Inhibitor Cocktail (100X) (cat # 78430) instead of Roche Complete Protease Inhibitor EDTA-Free tablets. CUT\u0026amp;TAG pAG-Tn5 protein was from Epicypher (SKU: 15-1017).\u003c/p\u003e\u003cp\u003eAnalysis was carried out using bowtie2 for alignment and SEACR for peak calling. Peak annotation was carried out using annotatePeaks.pl from HOMER. Alignment and annotations were carried out using the latest version of hg38. Heatmaps were generated using deepTools.\u003c/p\u003e\u003ch2\u003eEpigenetic characterization with deep learning\u003c/h2\u003e\u003cp\u003eAll algorithms and model training were developed in Jupyter with Python 3.10. Pre-trained ResNet models were downloaded using Pytorch and weights were frozen before unfreezing the last layer during training.\u003c/p\u003e\u003cp\u003eAll code used to carry out cell segmentation, training of ResNet models, and classification is publicly available on GitHub (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/zongmingchua/Epigenetic_image_ML\u003c/span\u003e\u003cspan class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eDATA AVAILABILITY\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe sequencing datasets generated and analyzed during the current study are not yet publicly available due to ongoing processing in the Gene Expression Omnibus repository but will be available from the corresponding author on reasonable request or upon publication of this manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAUTHORS CONTRIBUTIONS\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eP.D. A conceived the concept of the study, procured the funding, and was in charge of the overall direction. Z.M and H.T produced the results shown in Fig 1. Z.M produced the results shown in Figs 2-5. A.A assisted in producing plasmids encoding different substitution mutants at the H2A.Z R80 site. A.R and X.L assisted with advice in conducting bioinformatics and machine learning analysis. M.G.T, K.M, A.D, C.C.M, L.H, and N.D assisted with refining experimental protocols and general technical analysis. All authors including were involved in the preparation of the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFUNDING / FINANCIAL SUPPORT\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research was supported by the National Institute on Aging of the National Institutes of Health under Award Number R01AG071464. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCONFLICT OF INTEREST\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eHernandez-Segura, A., Nehme, J. \u0026amp; Demaria, M. Hallmarks of Cellular Senescence. \u003cem\u003eTrends Cell Biol.\u003c/em\u003e 28, 436\u0026ndash;453 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRai, T. S. \u003cem\u003eet al.\u003c/em\u003e HIRA orchestrates a dynamic chromatin landscape in senescence and is required for suppression of neoplasia. \u003cem\u003eGenes Dev.\u003c/em\u003e 28, 2712\u0026ndash;2725 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIvanov, A. \u003cem\u003eet al.\u003c/em\u003e Lysosome-mediated processing of chromatin in senescence. \u003cem\u003eJ. Cell Biol.\u003c/em\u003e 202, 129\u0026ndash;143 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCampisi, J. \u0026amp; d\u0026rsquo;Adda Di Fagagna, F. Cellular senescence: when bad things happen to good cells. \u003cem\u003eNat. Rev. Mol. Cell Biol.\u003c/em\u003e 8, 729\u0026ndash;740 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHayflick, L. THE LIMITED IN VITRO LIFETIME OF HUMAN DIPLOID CELL STRAINS. \u003cem\u003eExp. Cell Res.\u003c/em\u003e 37, 614\u0026ndash;636 (1965).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHayflick, L. \u0026amp; Moorhead, P. S. The serial cultivation of human diploid cell strains. \u003cem\u003eExp. Cell Res.\u003c/em\u003e 25, 585\u0026ndash;621 (1961).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, C. \u003cem\u003eet al.\u003c/em\u003e DNA damage response and cellular senescence in tissues of aging mice. \u003cem\u003eAging Cell\u003c/em\u003e 8, 311\u0026ndash;323 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen, H.-A. \u003cem\u003eet al.\u003c/em\u003e Senescence rewires microenvironment sensing to facilitate anti-tumor immunity. \u003cem\u003eCancer Discov.\u003c/em\u003e CD-22-0528 (2022) doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1158/2159-8290.CD-22-0528\u003c/span\u003e\u003cspan address=\"10.1158/2159-8290.CD-22-0528\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eReyes, N. S. \u003cem\u003eet al.\u003c/em\u003e Sentinel p16INK4a+ cells in the basement membrane form a reparative niche in the lung. \u003cem\u003eScience\u003c/em\u003e 378, 192\u0026ndash;201 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRitschka, B. \u003cem\u003eet al.\u003c/em\u003e The senescence-associated secretory phenotype induces cellular plasticity and tissue regeneration. \u003cem\u003eGenes Dev.\u003c/em\u003e 31, 172\u0026ndash;183 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMichaloglou, C. \u003cem\u003eet al.\u003c/em\u003e BRAFE600-associated senescence-like cell cycle arrest of human naevi. \u003cem\u003eNature\u003c/em\u003e 436, 720\u0026ndash;724 (2005).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXue, W. \u003cem\u003eet al.\u003c/em\u003e Senescence and tumour clearance is triggered by p53 restoration in murine liver carcinomas. \u003cem\u003eNature\u003c/em\u003e 445, 656\u0026ndash;660 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCopp\u0026eacute;, J.-P., Desprez, P.-Y., Krtolica, A. \u0026amp; Campisi, J. The Senescence-Associated Secretory Phenotype: The Dark Side of Tumor Suppression. \u003cem\u003eAnnu. Rev. Pathol. Mech. Dis.\u003c/em\u003e 5, 99\u0026ndash;118 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHernandez-Gonzalez, F. \u003cem\u003eet al.\u003c/em\u003e Human senescent fibroblasts trigger progressive lung fibrosis in mice. \u003cem\u003eAging\u003c/em\u003e 15, 6641\u0026ndash;6657 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRuhland, M. K. \u003cem\u003eet al.\u003c/em\u003e Stromal senescence establishes an immunosuppressive microenvironment that drives tumorigenesis. \u003cem\u003eNat. Commun.\u003c/em\u003e 7, 11762 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaker, D. J. \u003cem\u003eet al.\u003c/em\u003e Naturally occurring p16(Ink4a)-positive cells shorten healthy lifespan. \u003cem\u003eNature\u003c/em\u003e 530, 184\u0026ndash;189 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTchkonia, T., Zhu, Y., van Deursen, J., Campisi, J. \u0026amp; Kirkland, J. L. Cellular senescence and the senescent secretory phenotype: therapeutic opportunities. \u003cem\u003eJ. Clin. Invest.\u003c/em\u003e 123, 966\u0026ndash;972 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYi, S.-J. \u0026amp; Kim, K. New Insights into the Role of Histone Changes in Aging. \u003cem\u003eInt. J. Mol. Sci.\u003c/em\u003e 21, 8241 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePaluvai, H., Di Giorgio, E. \u0026amp; Brancolini, C. The Histone Code of Senescence. \u003cem\u003eCells\u003c/em\u003e 9, 466 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRogakou, E. P. \u0026amp; Sekeri\u0026ndash;Pataryas, K. E. Histone variants of H2A and H3 families are regulated during in vitro aging in the same manner as during differentiation. \u003cem\u003eExp. Gerontol.\u003c/em\u003e 34, 741\u0026ndash;754 (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTvardovskiy, A., Schw\u0026auml;mmle, V., Kempf, S. J., Rogowska-Wrzesinska, A. \u0026amp; Jensen, O. N. Accumulation of histone variant H3.3 with age is associated with profound changes in the histone methylation landscape. \u003cem\u003eNucleic Acids Res.\u003c/em\u003e 45, 9272\u0026ndash;9289 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eContrepois, K. \u003cem\u003eet al.\u003c/em\u003e Histone variant H2A.J accumulates in senescent cells and promotes inflammatory gene expression. \u003cem\u003eNat. Commun.\u003c/em\u003e 8, 14995 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDuarte, L. F. \u003cem\u003eet al.\u003c/em\u003e Histone H3.3 and its proteolytically processed form drive a cellular senescence programme. \u003cem\u003eNat. Commun.\u003c/em\u003e 5, 5210 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStefanelli, G. \u003cem\u003eet al.\u003c/em\u003e Learning and Age-Related Changes in Genome-wide H2A.Z Binding in the Mouse Hippocampus. \u003cem\u003eCell Rep.\u003c/em\u003e 22, 1124\u0026ndash;1131 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSuto, R. K., Clarkson, M. J., Tremethick, D. J. \u0026amp; Luger, K. Crystal structure of a nucleosome core particle containing the variant histone H2A.Z. \u003cem\u003eNat. Struct. Biol.\u003c/em\u003e 7, 1121\u0026ndash;1124 (2000).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGevry, N., Chan, H. M., Laflamme, L., Livingston, D. M. \u0026amp; Gaudreau, L. p21 transcription is regulated by differential localization of histone H2A.Z. \u003cem\u003eGenes Amp Dev.\u003c/em\u003e 21, 1869\u0026ndash;1881 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSales-Gil, R. \u003cem\u003eet al.\u003c/em\u003e Non-redundant functions of H2A.Z.1 and H2A.Z.2 in chromosome segregation and cell cycle progression. \u003cem\u003eEMBO Rep.\u003c/em\u003e 22, e52061 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e\u0026Aacute;vila-L\u0026oacute;pez, P. A. \u003cem\u003eet al.\u003c/em\u003e H2A.Z overexpression suppresses senescence and chemosensitivity in pancreatic ductal adenocarcinoma. \u003cem\u003eOncogene\u003c/em\u003e 40, 2065\u0026ndash;2080 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBelotti, E. \u003cem\u003eet al.\u003c/em\u003e H2A.Z is dispensable for both basal and activated transcription in post-mitotic mouse muscles. \u003cem\u003eNucleic Acids Res.\u003c/em\u003e 48, 4601\u0026ndash;4613 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBuschbeck, M. \u0026amp; Hake, S. B. Variants of core histones and their roles in cell fate decisions, development and cancer. \u003cem\u003eNat. Rev. Mol. Cell Biol.\u003c/em\u003e 18, 299\u0026ndash;314 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLuger, K., Dechassa, M. L. \u0026amp; Tremethick, D. J. New insights into nucleosome and chromatin structure: an ordered state or a disordered affair? \u003cem\u003eNat. Rev. Mol. Cell Biol.\u003c/em\u003e 13, 436\u0026ndash;447 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAndrews, A. J., Chen, X., Zevin, A., Stargell, L. A. \u0026amp; Luger, K. The Histone Chaperone Nap1 Promotes Nucleosome Assembly by Eliminating Nonnucleosomal Histone DNA Interactions. \u003cem\u003eMol. Cell\u003c/em\u003e 37, 834\u0026ndash;842 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAngelov, D. \u003cem\u003eet al.\u003c/em\u003e The Histone Variant MacroH2A Interferes with Transcription Factor Binding and SWI/SNF Nucleosome Remodeling. \u003cem\u003eMol. Cell\u003c/em\u003e 11, 1033\u0026ndash;1041 (2003).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCui, K. \u003cem\u003eet al.\u003c/em\u003e Chromatin signatures in multipotent human hematopoietic stem cells indicate the fate of bivalent genes during differentiation. \u003cem\u003eCell Stem Cell\u003c/em\u003e 4, 80\u0026ndash;93 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVald\u0026eacute;s-Mora, F. \u003cem\u003eet al.\u003c/em\u003e Acetylation of H2A.Z is a key epigenetic modification associated with gene deregulation and epigenetic remodeling in cancer. \u003cem\u003eGenome Res.\u003c/em\u003e 22, 307\u0026ndash;321 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHu, G. \u003cem\u003eet al.\u003c/em\u003e H2A.Z facilitates access of active and repressive complexes to chromatin in embryonic stem cell self-renewal and differentiation. \u003cem\u003eCell Stem Cell\u003c/em\u003e 12, 180\u0026ndash;192 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHorikoshi, N., Arimura, Y., Taguchi, H. \u0026amp; Kurumizaka, H. Crystal structures of heterotypic nucleosomes containing histones H2A.Z and H2A. \u003cem\u003eOpen Biol.\u003c/em\u003e 6, 160127.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi, S., Wei, T. \u0026amp; Panchenko, A. R. Histone variant H2A.Z modulates nucleosome dynamics to promote DNA accessibility. \u003cem\u003eNat. Commun.\u003c/em\u003e 14, 769 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJin, C. \u0026amp; Felsenfeld, G. Nucleosome stability mediated by histone variants H3.3 and H2A.Z. \u003cem\u003eGenes Dev.\u003c/em\u003e 21, 1519\u0026ndash;1529 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNacev, B. A. \u003cem\u003eet al.\u003c/em\u003e The expanding landscape of \u0026lsquo;oncohistone\u0026rsquo; mutations in human cancers. \u003cem\u003eNature\u003c/em\u003e 567, 473\u0026ndash;478 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eArimura, Y. \u003cem\u003eet al.\u003c/em\u003e Cancer-associated mutations of histones H2B, H3.1 and H2A.Z.1 affect the structure and stability of the nucleosome. \u003cem\u003eNucleic Acids Res.\u003c/em\u003e 46, 10007\u0026ndash;10018 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBerryman, S., Matthews, K., Lee, J. H., Duffy, S. P. \u0026amp; Ma, H. \u003cem\u003eImage-Based Cell Phenotyping Using Deep Learning\u003c/em\u003e. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://biorxiv.org/lookup/doi/10.1101/817544\u003c/span\u003e\u003cspan address=\"http://biorxiv.org/lookup/doi/10.1101/817544\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2019) doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1101/817544\u003c/span\u003e\u003cspan address=\"10.1101/817544\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeckenbach, I. \u003cem\u003eet al.\u003c/em\u003e Nuclear morphology is a deep learning biomarker of cellular senescence. \u003cem\u003eNat. Aging\u003c/em\u003e 2, 742\u0026ndash;755 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchwartzentruber, J. \u003cem\u003eet al.\u003c/em\u003e Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma. \u003cem\u003eNature\u003c/em\u003e 482, 226\u0026ndash;231 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKang, T.-W. \u003cem\u003eet al.\u003c/em\u003e Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. \u003cem\u003eNature\u003c/em\u003e 479, 547\u0026ndash;551 (2011).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCopp\u0026eacute;, J.-P. \u003cem\u003eet al.\u003c/em\u003e Senescence-Associated Secretory Phenotypes Reveal Cell-Nonautonomous Functions of Oncogenic RAS and the p53 Tumor Suppressor. \u003cem\u003ePLOS Biol.\u003c/em\u003e 6, e301 (2008).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCecco, M. D. \u003cem\u003eet al.\u003c/em\u003e L1 drives IFN in senescent cells and promotes age-associated inflammation. \u003cem\u003eNature\u003c/em\u003e 566, 73\u0026ndash;78 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSerrano, M., Lin, A. W., McCurrach, M. E., Beach, D. \u0026amp; Lowe, S. W. Oncogenic ras Provokes Premature Cell Senescence Associated with Accumulation of p53 and p16INK4a. \u003cem\u003eCell\u003c/em\u003e 88, 593\u0026ndash;602 (1997).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBagert, J. D. \u003cem\u003eet al.\u003c/em\u003e Oncohistone mutations enhance chromatin remodeling and alter cell fates. \u003cem\u003eNat. Chem. Biol.\u003c/em\u003e 17, 403\u0026ndash;411 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOrjalo, A. V., Bhaumik, D., Gengler, B. K., Scott, G. K. \u0026amp; Campisi, J. Cell surface-bound IL-1α is an upstream regulator of the senescence-associated IL-6/IL-8 cytokine network. \u003cem\u003eProc. Natl. Acad. Sci.\u003c/em\u003e 106, 17031\u0026ndash;17036 (2009).\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"npj-aging","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"Learn more about [npj Aging](https://www.nature.com/npjamd/)","snPcode":"41514","submissionUrl":"https://submission.springernature.com/new-submission/41514/3","title":"npj Aging","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-8818498/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8818498/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eCellular senescence features a durable cell-cycle arrest and a pro-inflammatory senescence-associated secretory phenotype (SASP) driven in part by chromatin remodeling. The histone variant H2A.Z plays an essential role in regulating gene expression through modulating nucleosome dynamics and is known to regulate the expression of cell cycle genes during the early stages of cellular senescence. However, how the intrinsic stability of H2A.Z-containing nucleosomes plays a role in the establishment of the senescent phenotype remains unexplored. To investigate this, we employed H2A.Z R80C, a H2A.Z mutant that destabilizes nucleosomes by disrupting histone-DNA interactions. We observed that expression of H2A.Z R80C causes suppression of SASP in senescent primary human fibroblasts, without affecting expression of cell cycle genes. H2A.Z knockdown did not suppress the SASP, demonstrating that SASP suppression is likely due to altered stability of H2A.Z-containing nucleosomes rather than loss of H2A.Z function. Mechanistically, SASP suppression is linked to decreased H3K27ac at SASP gene loci. These findings offer a novel avenue for understanding and manipulating the SASP during aging and other senescence-related pathologies.\u003c/p\u003e","manuscriptTitle":"Histone variant H2A.Z mutant suppresses the senescence-associated secretory phenotype","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-02-26 06:27:35","doi":"10.21203/rs.3.rs-8818498/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2026-03-17T13:29:54+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2026-03-16T15:41:35+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2026-03-08T13:17:35+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2026-03-05T16:27:12+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"263791656954893923157925577749278285975","date":"2026-03-03T21:11:06+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"327059443416181975869975463189396673521","date":"2026-03-03T16:14:02+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"288025275816362856217951916967021335839","date":"2026-02-24T16:14:53+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"271390456693019115975622440405192130668","date":"2026-02-24T00:04:14+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"111527800128300082497810089146851775503","date":"2026-02-23T23:28:15+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2026-02-23T22:55:58+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2026-02-23T18:59:47+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2026-02-20T04:40:43+00:00","index":"","fulltext":""},{"type":"submitted","content":"npj Aging","date":"2026-02-08T01:00:24+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"npj-aging","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"Learn more about [npj Aging](https://www.nature.com/npjamd/)","snPcode":"41514","submissionUrl":"https://submission.springernature.com/new-submission/41514/3","title":"npj Aging","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"49baaeea-8c7e-46e7-a7bb-c9534453464c","owner":[],"postedDate":"February 26th, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"in-revision","subjectAreas":[{"id":63548591,"name":"Biological sciences/Cell biology"},{"id":63548592,"name":"Biological sciences/Genetics"},{"id":63548593,"name":"Biological sciences/Molecular biology"}],"tags":[],"updatedAt":"2026-03-17T13:41:01+00:00","versionOfRecord":[],"versionCreatedAt":"2026-02-26 06:27:35","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-8818498","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8818498","identity":"rs-8818498","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}