Inhibition of histone H3-H4 chaperone pathways rescues C. elegans sterility by H2B loss

Oncohistone mutations are crucial drivers for tumorigenesis, but how a living organism governs the loss-of-function oncohistone remains unclear. We generated a histone H2B triple knockout (3KO) strain in Caenorhabditis elegans, which decreased the embryonic H2B, disrupted cell divisions, and caused animal sterility. By performing genetic suppressor screens, we uncovered that mutations defective in the histone H3-H4 chaperone UNC-85 restored H2B 3KO fertility by decreasing chromatin H3-H4 levels. RNA interference of other H3-H4 chaperones or H3 or H4 histones also rescued H2B 3KO sterility. We showed that blocking H3-H4 chaperones recovered cell division in C. elegans carrying the oncohistone H2BE74K mutation that distorts the H2B-H4 interface and induces nucleosome instability. Our results indicate that reducing chromatin H3-H4 rescues the dysfunctional H2B in vivo and suggest that inhibiting H3-H4 chaperones may provide an effective therapeutic strategy for treating cancers resulting from loss-of-function H2B oncohistone.


Introduction
Recent efforts to catalog histone mutations in cancer patients have identified many histone missense mutations as oncogenic drivers or contributors to tumor progression [1,2]. Histone globular domains form the core of the nucleosome, and the unstructured histone tails provide the target sites for various posttranslational modifications (PTM) [2][3][4]. The best-characterized oncohistone mutations distribute in the histone tail, including H3 K27M in diffuse intrinsic pontine gliomas, H3 G34V/R in pediatric glioblastomas, or H3 K36M in chondroblastomas [1,5,6]. In contrast to the well-known PTM-associated mutations, the most common oncohistone mutations localize in the globular histone fold domain, disrupting the structural integrity of the nucleosome [2,7]. Nucleosome assembly involves depositing a histone (H3-H4) 2 tetramer and rapidly following H2A-H2B heterodimer deposition [4]. The E76 residue of H2B interacts with the R92 residue of H4 through hydrogen bonding, and the H2B E76K oncohistone mutation profoundly destabilizes the interaction of H2B with H4 and inhibits histone octamer formation [2,7]. Despite the most frequent H2B oncohistone mutation, how a living organism controls H2B E76K remains elusive, and neither can we propose therapeutic strategies to treat H2B E76K tumors.
Functional redundancy is a significant technical hurdle for the genetic dissection of histone. A metazoan genome usually encodes 10-20 copies of the genes for each of the four core histone proteins [11,12]. We examined a complete inventory of green fluorescence protein (GFP)-tagged H2B knock-in C. elegans and uncovered four out of 17 H2B genes expressed in the germline and early embryos. By generating a series of H2B mutant alleles, we showed that deleting three of these four H2B genes (H2B 3KO) caused animal sterility with full penetration, which enabled us to perform forward genetic screens to isolate suppressors. Unexpectedly, we found that inhibiting H3-H4 histone chaperones or reducing H3-H4 histones restored the fertility of H2B 3KO animals. We provided evidence that the loss of H3-H4 histone chaperones rescued cell division failures in C. elegans bearing the H2B E76K mutation, suggesting that H3-H4 histone chaperones may be promising therapeutic targets for treating cancers resulting from dysfunctional H2B oncohistone mutations.

UNC-85 mutations restore the fertility of H2B triple knockout C. elegans
A metazoan genome usually encodes 10-20 copies of the genes for each histone protein [11,12]; the functional redundancy has been a significant technical hurdle for the genetic dissection of histone genes in animal cells. In C. elegans, 17 annotated H2B genes encode five distinct polypeptides (Figs 1A and S1) [12]. By examining a complete inventory of green Loss of H2B decreased brood size. N = 20-53, data are presented as mean ± SD (error bars). Statistical significance compared with fluorescence protein (GFP)-tagged H2B knock-in (KI) strains, we found that the C. elegans germline and early embryos only expressed four H2B genes encoding the HIS-58 subfamily protein (Figs 1B and S2 and S3). To study the four germline and embryonic H2B genes, we generated the single, double, or triple H2B knockout (H2B 3KO) strains (S4A and S4B Fig). Immunofluorescence detected an expected H2B reduction in H2B 3KO embryos (S4C and S4D Fig). The loss of H2B significantly decreased brood sizes in the single or double H2B KO strains, and the H2B 3KO animals induced the sterility phenotype with a 100% penetration ( Fig 1C).
We conducted forward genetic screens for suppressors that restore fertility in H2B 3KO (S5A Fig). Our genetic mapping, followed by whole-genome sequencing, identified the unc-85 (cas1100) mutant as a suppressor (Fig 1D-1F). The gene unc-85 encodes one of two C. elegans homologs of anti-silencing function 1 (Asf1) histone H3-H4 chaperone crucial for chromatin assembly [3,9,13]. The unc-85 (cas1100) allele alters a conserved glycine at residue 111 to arginine, mapped to the histone chaperone domain (Fig 1D and 1E). Using the structural information of the yeast ASF1, we showed that the UNC-85 G111 /Asf1 G110 localizes in β strand 7, a binding interface to H3-H4 histones (S5B and S5E Fig) [14]. Molecular dynamic simulation suggests that G111R mutation may perturb the interaction between the chaperone and H3-H4 (S5B- S5G Fig) and impair UNC-85 function. To understand the cellular behaviors of UNC-85 G111R , we compared the green fluorescence from the UNC-85::GFP and UNC-85 G111R ::GFP KI animals. We did not detect a noticeable change in the UNC-85::GFP G111R signal in the whole cell; however, the nucleus-to-cytoplasm ratios of UNC-85::GFP were consistently higher than that of UNC-85::GFP G111R in the early embryos or oocytes (Figs 2A and 2B and S6A and S6B and S1 Video), indicating that G111R mutation reduces the nuclear localization of UNC-85.

UNC-85 mutations rescue the meiosis defects of H2B 3KO animals
We followed spindle and chromosome dynamics to dissect the underlying cellular mechanism using GFP-tagged tubulin and mCherry-tagged histone H3 (HIS-72) in the C. elegans embryos the control with a matching color code is based on Student's t-test, n.s., not significant, �� p< 0.01, ��� p< 0.001. Three biological replicates calculated the fertile adult percentage. (D) Schematic of the C. elegans UNC-85 and ASFL-1 proteins. The amino acid changes (cas1067, cas1068, and cas1100) and the deletion allele of unc-85 (cas1571) or asfl-1 (cas1222) are indicated. The scale bar represents 20 amino acids. (E) Multi-sequence alignment of C. elegans (Ce) UNC-85 with ASFL-1 and its homologs in human (Hs), mouse (Mm), Drosophila (Dm), and Saccharomyces (Sc). The sequence alignment shows the protein sequence of UNC-85 histone chaperone domain flanking (± 20 aa) around the mutation isolated from the EMS screen, and the histone chaperone domain is highly acidic, which favors its interaction with histone. The red asterisk highlights the mutation isolated from the EMS screen. An asterisk ( � ) indicates positions 100% conserved in the alignment; a colon (:) indicates conservation between groups of strongly similar properties; a period (.) indicates conservation between groups of weakly similar properties. Hydrophobic residues are colored in blue, positive charge residues in yellow, negative charge residues in magenta, and polar residues in green. (F) Fertile adult percentage and brood size of the indicated strains. N = 10-70, values are presented as mean ± SD (error bars). Statistical significance compared with the control with a matching color code is based on Student's t-test, n.s., not significant, � p< 0.05, �� p< 0.01, ��� p< 0.001. Three biological replicates calculated the fertile adult percentage. https://doi.org/10.1371/journal.pgen.1010223.g001

PLOS GENETICS
Inhibition of Histone H3-H4 Chaperone Pathways Rescues C. elegans Sterility by H2B Loss (Fig 2C and 2D). The wild-type (WT) embryos undergo two consecutive chromosome segregation events (i.e., meiosis I and meiosis II) preceded by a single round of DNA replication. As the spindle microtubules (MTs) gradually decreased, inter-chromosomal MTs assembled between segregating homologous chromosomes, followed by polar body exclusion. The second meiotic spindles were similarly formed (Fig 2E and 2F and S2 Video) [15]. In H2B 3KO animals, when the oocyte underwent meiotic maturation, meiotic spindle assembly and anaphase I correctly occurred; however, the spindle midzone persisted, the first polar body did not form, and the meiosis II spindle failed to assemble (Fig 2E and 2F and S2 Video), leading to an ectopic replicative fate and embryonic lethality. In the H2B 3KO and unc-85 (cas1100) quadruple mutant animals, 83.6% of the examined oocytes completed meiosis and mitosis successfully ( Fig 2F and S3 Video). These results show that the loss of the H3-H4 chaperone UNC-85 rescued cell division failures in H2B 3KO animals.

Inhibition of other components in histone H3-H4 chaperone pathways suppresses the sterility in H2B 3KO
Histone chaperone pathways deliver nascent histone subunits from the cytoplasm to the nucleosome [8], and additional histone chaperones relay with ASF1 to transport H3-H4 and facilitate their deposition to chromatin ( Fig 3A) [9,16]. Next, we performed an RNAi screen to determine whether the loss of other histone chaperones rescues H2B 3KO sterility. After histone synthesis, the somatic nuclear autoantigenic sperm protein (NASP) presents newly folded histone H3-H4 dimers to acetyltransferase 1 HAT1 holoenzyme. Consistent with unc-85 and asfl-1 results, RNAi of nasp-1 or nasp-2, the orthologs of the human NASP, efficiently suppressed sterility of H2B 3KO animals by partially restoring the brood size ( Fig 3B). After ASF1 family proteins shuttled the acetylated H3-H4 dimers into the nucleus, the H3-H4 dimer in the Asf1-H3-H4 complex is transferred to the CAF-1 complex to form the (H3-H4) 2 tetramer, which deposits onto newly synthesized DNA. To deposit replication-independent histone variants, ASF1 interacts with the HIRA complex to hand over H3.3-H4 to assemble the nucleosome [9], and the chaperone DAXX-ATRX deposits H3.3-H4 at the pericentromeric and telomere regions [10]. We showed that RNAi of chaf-1 (the human CHAF1A) or hira-1 or dap-1 (the human HIRA and DAXX homologs) rescued brood size defects in H2B 3KO. Our data revealed that inhibition of the H3-H4 chaperones with distinct functions rescued the sterility of H2B 3KO. Like many other species, the C. elegans early embryos do not undergo active gene transcription until the eight-cell stage, excluding the possibility that these chaperones play distinct regulatory roles in gene expression. Instead, we suspect that the blockade of any chaperone will reduce the H3-H4 amount from chromatin, restoring the overall stoichiometry of histones within nucleosomes close to the normal level. By contrast, RNAi of H2A-H2B chaperone NAP1 or FACT orthologs nap-1 or spt-16 did not suppress the sterility phenotype ( Fig 3B). Thus, reducing components in the H3-H4 but not H2A-H2B chaperone pathways recovered the fertility of H2B 3KO animals.

PLOS GENETICS
Inhibition of Histone H3-H4 Chaperone Pathways Rescues C. elegans Sterility by H2B Loss

H3-H4 to H2A-H2B ratio is essential for cell divisions
How did the loss of H3-H4 chaperones restored fertility by the lack of H2B? Immunofluorescence of H2B 3KO and unc-85 quadruple mutant embryos did not reveal any increase in H2B compared to H2B 3KO embryos (S4C and S4D Fig), excluding the possibility of H2B recovery by inhibiting H3-H4 chaperones. Because the quadruple mutant embryos are viable but have a low H2B level similar to that in H2B 3KO embryos (S4C and S4D Fig), we argue that H2B reduction may not be the primary cause for H2B 3KO sterility. To determine whether the loss of H3-H4 chaperone reduced histone H3 proteins in the nucleus, we constructed the HIS-74 (H3)::GFP KI reporter. The whole-cell HIS-74::GFP protein fluorescence in the H2B 3KO and unc-85(cas1100) quadruple mutant strain was indistinguishable from WT in early embryos or oocytes (Figs 4A and 4B and S7A and S7B and S4 Video). However, the nucleus-to-cytoplasm ratio of H3 fluorescence intensity significantly decreased, indicating a reduction of H3 recruitment to the nuclei, consistent with the established role of UNC-85 as an H3-H4 chaperone. To explore the histone H3 content on chromatin, we generated genome-wide maps for H3 chromatin immunoprecipitation with the exogenous reference genome (ChIP-Rx), followed by massively parallel sequencing [17]. Compared to WT, we detected a marked reduction in the widespread distribution of the total H3 in the H2B 3KO and unc-85(cas1100) quadruple mutant animals (Figs 4C-4E and S7C). We performed immunostaining experiments to examine the histone H4 levels using an anti H4 antibody that recognizes all the 16 C. elegans histone H4 proteins with an identical amino acid sequence. We found that H2B and H4 protein levels were markedly reduced in the viable quadruple mutant embryos (S7D and S7E Fig), suggesting that the aberrant H2A-H2B to H3-H4 ratio might be responsible for the sterility in H2B 3KO. In support of this idea, RNAi of the H3 his-74 gene or H4 his-46 gene in H2B 3KO animals restored brood size to 16.3% or 9.8% of WT (Fig 3B), respectively. These results collectively show that reducing the chromatin H3 or H4 level rescued the H2B 3KO sterility.
To disrupt the H2A-H2B to H3-H4 balance alternatively, we used a heat-shock inducible system to express H2B or H3 in the embryos under the control of the Phsp-16.41 promoter (S8A Fig). The heat-shock induction led to a higher frequency of chromosome segregation defects, including anaphase bridges and lagging chromatin, and the embryos overexpressing H2B or H3 exhibited embryonic lethality with 38.2% or 52.4% penetrance, respectively (S8B- S8D Fig). Using the heat shock overexpression system, we observed that H2B or H3 overexpressed C. elegans embryos underwent abnormal chromatin segregation and defective cell division, consistent with the early results that overexpressed histone genes disrupted cell division in yeast and human cells [18][19][20]. Although we observed markedly higher penetrance of cell division failure in embryos overexpressing histone than the wild-type embryos treated by heat shock, we cannot exclude the effects of heat shock on embryonic development because heat shock of WT embryos caused 10~20% embryonic lethality [21,22]. These findings are consistent with the early yeast histone overexpression results [18] and support the notion that the proper H2A-H2B to H3-H4 ratio is essential for organizing chromatin and cell division.

Inhibition of H3-H4 chaperones rescues abnormalities caused by H2B oncohistone in C. elegans
We wondered whether inhibiting H3-H4 chaperones could rescue cellular abnormalities caused by the loss-of-function H2B E76K oncohistone mutation [1,2]. We generated a C. elegans strain that changed the germline H2B his-48(cas1172) to E74K, equivalent to the human E76K (Fig 5A and 5B). The his-58(cas1239) and his-66(cas1240) double knockouts (H2B 2KO) or HIS-48 E74K animals were fertile ( Fig 5C); however, introducing HIS-48 E74K into H2B 2KO decreased brood size and caused embryonic lethality in a temperature-dependent manner ( Fig  5C and 5D). Although the embryos completed oocyte meiosis I and assembled meiotic II spindle, their chromosomes appeared more dispersed with the spindle relative to WT during anaphase II. The spindle midzone persisted and failed to extrude the polar body at the end of meiosis II (Fig 5E and 5F and S5 Video), generating multiple maternal pronuclei that disrupted the subsequent mitotic divisions (S8E Fig). By genetically crossing unc-85 (cas1100) or asfl-1 (cas1222) into the strain carrying H2B 2KO and H2B E74K , we found an increased brood size and a reduction of embryonic lethality (Fig 5C and 5D and S6 Video). Live-cell imaging analyses showed that unc-85 (cas1100) successfully recovered meiosis II and subsequent mitoses in the embryos (Fig 5E and S6 Video), indicating the loss of H3-H4 chaperones rescues the defects caused by the dysfunctional oncohistone H2B E74K in C. elegans. We did not notice that the single HIS-48 E74K mutation causes apparent sterility (Fig 5C), likely due to the functional redundancy of histone genes. We speculate that the other three wild-type H2B genes generate normal H2B proteins that sustain embryo development in the his-48 E74K single mutant. However, if we deleted another two H2B genes, i.e., his-58 KO and his-66 KO , the only left WT H2B gene may not produce a sufficient amount of normal H2B protein to maintain the proper chromatin structure that supports early development in HIS-48 E74K , his-58 KO and his-66 KO triple mutant embryos. Besides, we deleted two germline H2B genes and changed the other two H2Bs into H2B E74K , and unc-85 (cas1100) did not restore the fecundity defects in the H2B quadruple mutants, indicating that H3-H4 chaperone-mediated restoration requires WT H2B. Because cancer cells are unlikely to mutate all the isoforms of their histone genes to oncohistone mutations, inhibiting histone chaperones holds the promise to treat oncohistone-induced cellular abnormalities.

Discussion
This study demonstrates the advantage of using a genetically tractable system to investigate dysfunctional oncohistone mutations in a living organism. By generating a series of H2B mutant alleles, we found that deleting three of four germline H2B genes caused animal sterility with full penetration, which enabled us to perform forward genetic screens to isolate suppressors. Unexpectedly, we uncovered that inhibiting H3-H4 histone chaperones or reducing H3-H4 histones restored the fertility of H2B 3KO animals. Inspired by the genetic results, we provided evidence that the loss of H3-H4 histone chaperones rescued cell division failures in C. elegans bearing the dysfunctional oncohistone H2B E74K mutation. Because the loss of histone H3-H4 chaperones caused a reduction of chromatin H3-H4 (Fig 4C-4E), the decreased H3-H4 amount might be comparable to the remaining functional H2A-H2B in chromatin, thereby restoring cell division and fertility. We argue that the loss of histone amount or activity can be tolerated in vivo as long as histone subunits are proportionally reduced.
Consistent with this argument, overexpression of histone subunits caused chromosome loss in diploid budding yeast by perturbing chromatin structure in Saccharomyces cerevisiae [18,19,23,24], but co-overexpression of all the core histone subunits reduced the frequency of chromosome loss [18]. Overexpression of histone resulted in embryo lethality and a high  (S8C and S8D Fig), demonstrating the correct histone subunit ratio's crucial role in animal cells. Biochemical results show that an equimolar ratio of H3-H4 to H2A-H2B is vital for chromatin formation [25]. Recognizing the importance of the balanced histone composition, we suggest that inhibition of H2A-H2B chaperone pathways could suppress the phenotypes caused by dysfunctional H3 and H4 oncohistone mutations. Oncohistone mutations are not ideal for drug development; but our study uncovers that the cellular and animal defects caused by oncohistones can be rescued by manipulating a previously unrecognized process independent of oncohistones. We anticipate that the development of small molecules that target the pathways other than the mutated histone genes holds the promise to devise effective strategies for treating oncohistone-induced tumorigenesis. In support of this idea, inhibiting the JNK pathway recovered the aberrant C. elegans oocytes resulting from H3K27M-oncohistone mutation and suppressed the abnormal cell replicative fate in the human H3.3 K27M tumor cells [26].

C. elegans strains, culture, and genetics
All the C. elegans strains were genetic derivatives of the wild-type (WT) strain Bristol N2, except for Hawaii strain CB4856 in SNP mapping. Strains were raised on NGM agar plates seeded with E. coli strain OP50 at 20˚C following standard procedures [27]. Some strains were obtained from the Caenorhabditis Genetics Center (CGC) funded by the NIH Office of Research Infrastructure Programs (P40 OD010440). Strains in this study are listed in S1 Table.

Genome editing and molecular biology
We perform genome editing according to established criteria. The sgRNA sequence was designed using the CRISPR design tool (https://zlab.bio/guide-design-resources) and was inserted into the plasmid pDD162 by linearizing the vector with 15 bp overlapped primers. To make knock-in and mutant strains, 1-to 1.5-kb upstream and downstream homology arms were amplified and inserted into pPD95.77. sgRNA construct and homology recombination (HR) template were co-injected into N2 with pRF4 [rol-6 (su1006)] and Podr-1::dsRed selection markers. Marker-positive F1 animals were singled and screened by PCR and Sanger sequencing [28]. S2 Table lists

Forward genetic screens
We used forward genetic screens to isolate H2B 3KO suppressors. The heterozygous H2B 3KO mutant animals (P0) were synchronized at the late L4 larval stage, collected in 4 mL M9 buffer, and incubated with 50 mM ethyl-methanesulfonate (EMS) for 4 hours at room temperature with constant rotation. Animals were then washed with M9 three times and cultured under standard conditions. After 20 hours, adult animals were bleached. Eggs (F1) were distributed indicated. Scale bar, 60 bp. (B) Multi-sequence alignment of C. elegans (Ce) H2B and its homologs in human (Hs), mouse (Mm), Drosophila (Dm), and Saccharomyces (Sc). A red asterisk highlights the oncohistone mutation. (C) Fertile adult percentage and brood size of WT, Mut I, Mut II, Mut II; unc-85(cas1100), Mut II; asfl-1(cas1222) and Mut III animals at 25˚C. Information of Mut I and Mut II is in (A). N = 9-33, Data are presented as mean ± SD (error bars). Statistical significance compared with the control with a matching color code is based on Student's t-test, ��� p < 0.001. Three biological replicates calculated the fertile adult percentage. (D) Quantification of embryonic lethality rates in Mut I, Mut II, Mut II; unc-85(cas1100), Mut II; asfl-1(cas1222) animals at 25˚C. n = 500-1612, from three biological replicates. Statistical significance compared with the control with a matching color code is based on Student's t-test, n.s., not significant, ��� p < 0.001. (E) Fluorescence time-lapse images of GFP-tagged tubulin with mCherry-tagged histone H3 and membrane (magenta) during oocyte meiosis in Mut I, Mut II, and Mut II; unc-85(cas1100) at 25˚C. Percentages of meiosis patterns are indicated. Scale bar, 5 μm. (F) Enlarged schematic timescales of oocyte meiosis in Mut II.

Germline dissection
Young adult hermaphrodite animals were cut with a syringe at the tail or head region in the gonad dissection buffer to expose their anterior or posterior gonads [29]. The dissected gonads were fixed in -20˚C pre-chilled methanol for 10 min and then transferred to poly-L-lysine coated slides. We rinsed slides with PBS, touched the edges of slides on Kimwipes to remove excess PBS, applied 30 μl mounting medium containing DAPI to the slides, and slowly tipped the coverslip onto the mounting medium to avoid air bubbles, dry at room temperature for 5-10 min and seal the edges of coverslip with nail polish.

RNA interference and candidates screening
RNAi by feeding was performed as described [30]. In brief, we performed a candidate RNAi screen using H2B 3KO worms. The RNAi clones were from the Ahringer C. elegans library and confirmed by sanger sequencing. Bacteria were cultured in LB plus carbenicillin at 50 μg ml -1 and tetracycline at 12.5 μg ml -1 overnight and then seeded on NGM plates containing 50 μg ml -1 carbenicillin, 12.5 μg ml -1 tetracycline, and 1 mM IPTG. Bacteria containing the empty vector L4440 were used as the control.

Live-cell imaging
C. elegans were anesthetized with 0.1 mmol/L levamisole in M9 buffer for live-cell imaging of meiosis and germline oocyte. We synchronized young adult worms to image early embryos, cut them with syringes to release embryos in M9 buffer, and then mounted them on 2% (wt/vol) agarose pads at 20˚C [31], time-lapse images were taken by μManager (www.micro-manager. org) with an exposure time of 300 ms every 30-60 s. Our imaging system includes an Axio Observer Z1 microscope (Carl Zeiss MicroImaging, Inc.) equipped with a 100X, 1.45 NA objective, an EM CCD camera (Andor iXon+ DU-897D-C00-#BV-500), and the 488-nm and 561-nm lines of a Sapphire CW CDRH USB Laser System to a spinning disk confocal scan head (Yokogawa CSU-X1 Spinning Disk Unit). ImageJ software (NIH) was used to process and quantify the images (http://rsbweb.nih.gov/ij/). The background fluorescence was subtracted for the intensity measurements of UNC-85 and HIS-74 distribution. The whole-cell and nucleus fluorescence intensity of UNC-85 and HIS-74 were measured, and the nucleus to cytoplasm fluorescence ratio (N/C ratio) was calculated to reflect the distribution of UNC-85 and HIS-74.

Self-brood size assay
Brood size counts were performed by picking individual late L4 stage hermaphrodites to separate plates and transferring daily to new plates until no more eggs were laid. Live progeny were counted at the late larva to adult stage [32].

Heat-shock treatment and viability
For C. elegans embryos viability analysis, 30 young adult worms were collected on a seeded NGM plate to lay approximately 100 eggs. The embryos were heat-shocked for 1 hour at 33˚C and cultured at 20˚C for three days. Viability is defined as the number of hatched worms divided by the original egg numbers.

Molecular dynamic simulation and structural analysis
Wild-type yeast Asf1-H3-H4 complex (PDB code: 2HUE) and free Asf1 (PDB code: 1ROC) structures were downloaded from PDB, and G110R mutation structures were created by USCF Chimera (RRID:SCR_004097) with energy minimization. Molecular dynamics simulations were performed using GROMACS (RRID:SCR_014565) version 4.6.7 using the OPL-S-AA force field parameter set. After dispersion correction, each system was solvated in 150 mM NaCl solvent with SPC/E water models in a cubic box. Subsequently, neutralizing counter ions were added, and the steepest descent energy minimization was performed, followed by a two-step equilibration consisting of 100 ps of isochoric-isothermal (NVT) equilibration and 100 ps of isothermal-isobaric (NPT) equilibration. All position restraints were removed, and simulations were performed for 5 ns. RMSF and RMSD analysis calculate the standard deviation (root mean square fluctuation/deviation) of the atomic positions of the specified amino acids compared to their starting position within the energy minimized and equilibrated structure. Simulations were performed in triplicate; values shown in figures are averages with a standard deviation of the mean highlighted as error bars. Structural analysis and images were performed and created using PyMol (RRID:SCR_000305). The electrostatic distribution surfaces were generated using the APBS tool (RRID:SCR_008387).

Chromatin immunoprecipitation
We performed chromatin immunoprecipitation as described previously with some modifications [33]. Synchronized young adults were harvested and washed three times with M9 to remove excess bacteria. Using a 200 μl pipette tip, we dripped the worm mix into a 100 ml beaker that contains liquid nitrogen to make frozen worm "popcorn". Then, we ground the worm "popcorn" into frozen powder with a mortar and crosslinked with 1.1% formaldehyde for 15 min at room temperature. Formaldehyde was quenched by the addition of glycine, and the fixed powder was washed with cold PBS 3 times. The fixed powder was resuspended and homogenized 2-3 strokes using a glass Dounce homogenizer that is assumed to disaggregate cells gently. Nuclei purification and chromatin digestion were conducted using SimpleChIP Enzymatic Chromatin IP Kit (CST, #9003). Briefly, we resuspend cells using buffer A and incubate them on ice for 10 min. Pellet nuclei were collected and resuspended with buffer B. We added the micrococcal nuclease and incubated the reaction for 4-12 min at 37˚C to digest DNA to the length of approximately 150-900 bp. We stopped the digestion by adding EDTA. Then samples were resuspended in ice-cold FA buffer with 0.1% SDS and sonicated using a Bioruptor sonication system at 70% power for 20 cycles of 30 s on and 30 s off. The sonicated samples were then centrifuged at 4˚C for 15 min at 16,000g, and we transferred the supernatant to a new tube. The crosslinked chromatin concentration was measured by BCA Protein Assay Kit (TIANGEN, #PA115), and 2% of the extract was saved as the input sample. 5 μg chromatin was precleaned with protein A/G agarose beads. We added 10 ng Drosophila line 2 (S2) cells spike-in chromatin (Active Motif, # 53083) to each sample and immunoprecipitated overnight using 2 μg H3 antibody (Abcam, #AB195277) and 0.5 μg spike-in antibody (Active Motif, #61686). We added 30 μl protein A/G agarose beads to each ChIP sample and rotated them at 4˚C for 2 h. We added 150 μl ChIP Elution Buffer and vortex gently at 1200 rpm for 30 min at 65˚C. Then we reversed the crosslink with proteinase K at 65˚C for 2h. DNA was purified using PCR Purification Kit (QIAGEN, #28106), and DNA library construction was conducted with DNA Library Prep Kit for Illumina V3 (Vazyme, #ND607). The library was sequenced on NovaSeq systems (Illumina) using 150 nucleotides kits and paired-end sequencing (Novogene).

ChIP-seq data analysis
Raw reads were filtered using cutadapt software (version 2.5) for clean reads and checked with fastqc software (version 0.11.5). Clean reads were aligned to a mixed reference worm genome using BWA (version 0.7.17). The reference genome was mixed in advance with Caenorhabditis elegans DNA toplevel assembly (WBcel235 from the ensemble) and Drosophila melanogaster DNA toplevel assembly (BDGP6.32 from the ensemble). The aligned reads were separated into two sequence-alignment map (SAM) files. The SAM file of Caenorhabditis elegans was then converted into binary alignment map (BAM) files with samtools software (version 1.9) for further analysis, and the SAM file of Drosophila melanogaster was used for calculating spike-in factors [34]. The sorted and indexed BAM file of Caenorhabditis elegans underwent filtration using sambamba software (version 0.6.6) with 'view -h -f bam -F "not unmapped and not duplicate" '. The filtered BAM file was then converted to bigwig(bw) file using deeptools software (version 3.3.1) with 'bamCoverage--binSize 10--scaleFactor α--effectiveGenomeSize 100286401--extendReads--ignoreDuplicates--bam 'where α stood for the Spike-in factors calculated from SAM file of Drosophila melanogaster [35]. The Spike-in factors were calculated as David A. described with this formula [17]: We generated heatmaps and profile maps of Histone H3 occupation using deeptools software (version 3.3.1). We used Integrative Genomics Viewer (IGV, version 2.7.2) to visualize signal tracks to view bigwig files generated from deeptools software [36]. The comparison heatmaps in IGV between samples were generated using 'bigwigCompare 'in deeptools software.