Despite its sequence identity with canonical H4, Drosophila H4r product is enriched at specific chromatin regions

Histone variants are different from their canonical counterparts in structure and are encoded by solitary genes with unique regulation to fulfill tissue or differentiation specific functions. A single H4 variant gene (His4r or H4r) that is located outside of the histone cluster and gives rise to a polyA tailed messenger RNA via replication-independent expression is preserved in Drosophila strains despite that its protein product is identical with canonical H4. In order to reveal information on the possible role of this alternative H4 we epitope tagged endogenous H4r and studied its spatial and temporal expression, and revealed its genome-wide localization to chromatin at the nucleosomal level. RNA and immunohistochemistry analysis of H4r expressed under its cognate regulation indicate expression of the gene throughout zygotic and larval development and presence of the protein product is evident already in the pronuclei of fertilized eggs. In the developing nervous system a slight disequibrium in H4r distribution is observable, cholinergic neurons are the most abundant among H4r-expressing cells. ChIP-seq experiments revealed H4r association with regulatory regions of genes involved in cellular stress response. The data presented here indicate that H4r has a variant histone function.

The Histone 4 replacement gene (H4r) encodes a protein identical in amino acid sequence with its canonical histone H4 counterpart. H4r has been identified in 14 out of 22 sequenced Drosophila species so far 1 . Unlike canonical H4 genes, which are found in multiple copies within the histone cluster on the second chromosome 2 , H4r, similarly to the other histone variants, is located outside of the canonical histone cluster in a single copy. The third chromosomal H4r gene contains an intron, expressed independently of replication and produces polyadenylated mRNA product 2 . H4r mRNA is only weakly expressed in the germline and shows much higher expression in terminally differentiated cells. This observation led to the hypothesis on H4r replaces canonical H4 in postmitotic cells due to its replication independent expression 2 . Another hypothesis on the role of the alternative H4 is that it may play a role in environmental stress response, as mRNA expression from H4r changes upon ethanol treatment 3 . Since H4r encodes the same amino acid sequence as the canonical H4 in 14 Drosophila species with a codon usage different from that of H4, a further theory for the function of H4r is that it is incorporated to the chromatin with distinct co-or post-translational modifications in distinct cell types or in different environmental conditions 1 . Deletion of H4r does not cause a visible phenotypic change 4 , however, the loss of H4r causes reduced viability; female H4r mutants showing lower viability than males. Loss of H4r is also coupled to increased heat-stress resistance, supposedly due to the less condensed chromatin that allows a quicker and stronger response to heat stress 5 .
Although little data are available on the expression of H4r, on NCBI and on FlyBase it is indicated that the amount of H4r may differ significantly in distinct cell types of a tissue. For instance, experiments performed on larval and adult brains at various stages show that the amount of repressive marker PcG increases on H4r by the progression of differentiation whereas the amount of RNA polymerase II decreases, suggesting that H4r is expressed in undifferentiated cells, but not or at a lower level in mature neurons 6 . There are currently no data Scientific Reports | (2022) 12:5007 | https://doi.org/10.1038/s41598-022-09026-x www.nature.com/scientificreports/ available in the literature about the genomic distribution and interaction partners of H4r. Intrigued by the fact that a replacement histone gene with a gene product of identical structure but different expression pattern as the corresponding canonical histone is preserved in Drosophila species, we performed a detailed analysis of H4r expression and H4r localization in order to gain information on the function of this so far enigmatic gene and its product.

Results
Tagging H4r with 3xFlag-tag. To overcome the problem that the amino acid sequences of the products of H4r and H4 are identical making the two proteins indistinguishable, we fused a 3xFlag epitope tag encoding DNA sequence to the genomic H4r gene retaining its original expression pattern (Supplementary Figure S1). As a result, with the use of a Flag-specific antibody the product of H4r can be detected among the structurally identical H4 proteins. We created fly strains in which the 3xFlag-tag is located at the N-terminus of H4r and analysed it from a developmental, cellular and molecular perspective. To facilitate identification of labelled H4r animals, we included a dsRed marker gene flanked by two loxP recombination sites downstream of H4r. Because expression of the red fluorescent marker protein dsRed in brain interferes with analysis of images of immunostained brains, it was necessary to remove this marker gene before sample preparation. We achieved this by Cre-mediated recombination.
Analysis of H4r expression. Taking advantage of the Flag epitope tagged H4r expressed under its canonical control, we performed immunohistochemical stainings to detect H4r presence in different stages of development in different tissues and cell types. We found that H4r was detectable in both male and female pronuclei of embryos, and the expression of H4r remained ubiquitous throughout embryonic development ( Fig. 1a-d). In the brain of wandering larvae and adults different levels of H4r expression was observed ( Fig. 1e-g). Anti-Flag and DAPI double staining did not show clearly observable differences in the chromatin compaction of the Flagpositive cells. In order to identify cell types expressing H4r in larval brain we generated transgenic lines in which cell type-specific Lamin B-GFP expression could be achieved using the Gal4-UAS system. For the Drosophila stocks used for this experiment see Supplementary Table S3. By using these lines we sought to identify cell types in which colocalization of H4r with Lamin B-GFP could be detected. Analysing expression by immunohistochemistry in Drosophila lines expressing 3xFlag-H4r and Lamin B-GFP under the control of OK371, Gad1, ChAT, elav and Insc permitting the identification of glutaminergic, GABAergic, cholinergic neurons, mature neurons in general and neuroblasts, respectively, we found partial colocalization of lamin and Flag positive nuclei in every Drosophila line. This observation indicates that the expression of H4r is not restricted strictly and characteristic uniquely for any of the investigated neuron types. Nonetheless, we found that although the colocalization was not perfect, the majority (53.38% of cells in the whole brain, 69.2% in the eye discs) of H4r accumulating cells was cholinergic neurons (Fig. 2). Weak level of colocalization were found with mature neurons in general, with glutaminergic and GABAergic neurons and with neuroblasts (Table 1 and Supplementary Figure.  The expression of H4r does not change significantly upon heat stress. H4r was suggested to play a role in the chromatin formation at the loci of inducible genes such as heat shock genes 5 . We wondered whether there are changes in H4r expression upon heat shock and following recovery, which would further clarify the role of H4r in the expression of inducible genes. Therefore, we measured the changes in H4r expression at mRNA and protein levels upon heat stress and recovery. In addition we also determined if change in the ratio of soluble and chromatin associated H4r was detectable under the above conditions. We found no significant change in H4r expression at either mRNA (normalised to tubulin; n = 2; p ≥ 0.1357; one-way ANOVA with Sidak's multiple comparison test) or protein level, nor did we find alteration in the distribution of H4r protein between free (n = 2; p ≥ 0.8675; two-way ANOVA with Tukey's multiple comparison test) and chromatin associated forms (n = 2; p ≥ 0.7529; two-way ANOVA with Tukey's multiple comparison test) (Fig. 3).
H4r preferentially binds to specific chromosomal loci. One way of obtaining hints on H4r specific function could be the determination of its association with loci of particular chromosomal regions or genomewide. Staining of polytene chromosomes gave no information in this respect since H4r seemed to be associated with regions of giant chromosomes throughout roughly in inverse proportion of RNA polymerase II. In order to examine genomic distribution of H4r in diploid cells with a better resolution we performed chromatin-immunoprecipitation followed by sequencing (ChIP-seq) assays for epitope tagged H4r with Flag antibody. In parallel we also detected H3 distribution. Our aim was to detect the genome-wide nucleosome occupancy and analyze H4r distribution. As H3 is present in nucleosomes with H4 in equimolar quantity the two can provide identical information about nucleosome occupancy. Therefore instead of H4 antibody, we used an H3-specific antibody in these experiments as the latter one gave more consistent results in ChIP assays. In addition, we also detected the localization of the H3 variant H3.3. To achieve this, for chromatin preparation we used a Drosophila line which expresses transgenic H3.3 fused with a 3xFlag tag under the control of elav-Gal4. Thus, as the H3-antibody recognises both the canonical H3 and H3.3 but the Flag-antibody only H3.3, we could determine -using the same samples -H3 localization reflecting chromatin compaction and also H3.3 distribution. We found extensive similarities between H4r and H3.3 localizations and significant differences in the localization of these variants and that of canonical H3 (Fig. 4a).
Results of repeated ChIP-seq experiments indicated a reproducible non-random and non-uniform localization of H4r throughout the genome (Fig. 4b). Overall, H4r and H3.3 shared extensive similarity in their genome-wide distribution (Spearman's rank correlation coefficient ρ = 0.78, p = 0.0041, Fig. 4c)   www.nature.com/scientificreports/ overlap between the localization of the two histone variants was far from being perfect, suggesting functional differences for the two alternative histone forms. Notably, H4r distribution did not show strong correlation with that of H3, which can be assumed to reflect canonical H4 distribution (Spearman's rank correlation coefficient ρ = 0.35, p = 0.0187, see Fig. 4c), but was similar to H3.3 (as it was expected, the correlation between H3.3 and canonical H3 was higher, Spearman's rank correlation coefficient ρ = 0.47, p = 0.1695, see Fig. 4c). Identification of the regions at which H4r was specifically found in higher abundance revealed that this variant histone was most frequently bound to promoter regions (61.3% of total H4r were found at promoter regions, 14.0% at gene bodies and 24.7% at distal intergenic regions; 61.4% of total H3.3 were found at promoter regions, 15.8% at gene bodies and 22.8% at distal intergenic regions). On the contrary, a lower enrichment of H4r was detected in intergenic regions compared to the canonical H3 (22.0% of total H3 were found at promoter regions, 27.9% at gene bodies and 50.1% at distal intergenic regions) (Fig. 5). Altogether, these results indicated H4r to be more similar functionally to a histone variant, such as H3.3 than to a canonical histone.
Next, we analyzed which functional groups of genes show H4r enrichment. According to a PANTHER GO-Slim Biological Process analysis H4r, similarly to H3.3, shows enrichment on genes coupled with differentiation or necessary for normal functions -including genes inducible by various stimuli -and did not show significant enrichment on genes associated with cell cycle and cell division (Supplementary Table S1).
Interestingly, out of 2297 genes enriched in H4r and from the 2294 enriched with H3.3, 1479 genes were identical. According to expressional data of FlyBase RNA-seq expression profile, most genes showing H4r and H3.3 localization are highly expressed in the larval central nervous system and in adult brain. We found that the amount of H4r and H3.3 relative to that of the canonical H3 is significantly higher on genes that are inducible or highly expressed in the adult brain than on those genes that are weakly or not expressed (Fig. 6a). H4r/H3 ratio was 1.833 at inducible genes, 1.907 at highly expressed genes and 0.82 at weakly or not expressed genes; H3.3/H3 ratio was 1.667 at inducible genes, 1.859 at highly expressed genes and 0.937 at weakly or not expressed genes. The difference was not significant in case of the H4r/H3 ratio between the inducible and highly expressed genes (p > 0.9999) but it was significant between inducible or highly expressed genes and weakly/not expressed genes (p < 0.0001). The differences in the H3.3/H3 ratios were significant between the inducible and highly expressed genes (p = 0.0147), and between highly or inducible expressed genes and weakly/not expressed genes (p < 0.0001).
On genes showing enrichment for H4r and H3.3 the amount of these variants changes differently upon heat stress and recovery. Increased resistance to heat stress of H4r knock-out flies have been described recently 5 , and we found H4r enrichment on genes involved in environmental stress responses. These observations prompted us to examine how H4r localization changes upon activation of the heat stress response and after a period of recovery time following that. We performed ChIP-seq experiments using chromatin samples obtained from heads of flies exposed to heat stress at 37 °C for 20 min and from heads of animals, which were allowed to a 3 h recovery at 25 °C after an identical heat stress. The changes in the amount of H4r and H3.3 on the Hsp genes upon heat shock and following recovery are shown on Fig. 6b, c. The level of H4r  www.nature.com/scientificreports/ and H3.3 decreased, however, not significantly upon heat shock (p = 0.1266 and p = 0.1662, respectively), and increased significantly upon recovery (p < 0.0001 and p = 0.0006, respectively). The ratio of H4r to canonical H3 visibly increases upon heat shock (p = 0.0069) and shows a mild, but not significant increase after recovery as well (p > 0.9999). In contrast to H4r, although the ratio of H3.3 to canonical H3 increases upon heat shock (p < 0.0001), it mildly decreases after recovery (p = 0.1266). The results show a similar distribution for the two examined alternative histones upon transcription activation of Hsp genes, but indicate a difference in their distribution upon transcription silencing following the activation. According this result, nucleosomes built in during and following recovery contain more canonical H3 than H3.3, but the increasing H4r/H3 ratio suggests that these nucleosomes contain mostly H4r.

Discussion
In this study we created an experimental system in which the alternative histone H4r and canonical H4 can be distinguished despite the same amino acid sequences they have. We used this experimental tool to examine the expression and genomic localization of H4r and to compare it to another alternative histone, to the more extensively studied H3. 3. In previous studies based on RNA expression of H4r and deletion mutant phenotype analysis, two hypothesis arose about the function of H4r. One of these states that H4r replaces canonical H4 in the post-mitotic cells where canonical H4 is not expressed, and the other hypothesis states that H4r may play role in environmental stress responses.
Here we demonstrated that H4r is present already in the pronuclei of the fertilized embryos, which means that H4r is transferred to the eggs as a maternal gene product. The presence of H4r remains ubiquitous during the process of embryogenesis, and becomes to some extent cell type-specific upon neuronal differentiation. H4r differs from the H3.3 histone variant in its presence in nuclei during early embryogenesis as the latter has been  Figure S2). Results of ChIP-seq experiments performed with Flag-antibody specific for the 3xFlag-tagged H4r showed that the distribution of H4r throughout the genome is more similar to the other examined variant histone H3.3 than to the canonical histone H3. Both histone variants were detected preferentially bound to promoter regions of genes whereas canonical histone H3 was more abundant in intergenic regions. These findings support the assumption that H4r has a specific histone variant function. Moreover, compared to canonical H3 both histone variants were significantly more abundant on genes that are highly expressed or inducible than on weakly expressed and constitutively inactive genes. www.nature.com/scientificreports/ On induced Hsp gene promoters, similarly to H3.3, the amount of H4r increases upon heat shock relative to canonical H3. This notion supports the assumption that these histone variants play role in/during transcription activation. H3.3 has already been described to be involved in transcription activation 8 . H4r seems to behave similarly upon heat induction. While the amount of promoter localized H4r decreases upon heat shock, the ratio of nucleosomes that contain H4r increases, as it is seen in the case of H3.3 as well. We noticed that following activation, during recovery when transcription on Hsp genes is silenced, some H3.3 remains incorporated in nucleosomes at Hsp gene promoters but most of the newly integrated nucleosomes contain canonical H3 (Fig. 6c). H4r remains incorporated in newly inbuilt nucleosomes with higher frequency. It has been shown that H3.3 expression increased upon heat shock and most of H3.3 got incorporated into the chromatin, but during recovery H3.3 association to chromatin decreased 9 . Interestingly, H4r gene is not induced during heat stress and the level present in the cells is sufficient to fulfill H4 replacement connected to transcription activation of Hsp genes. Taken together these findings indicate that H4r and H3.3 may be involved in the transcription activation of inducible genes and these alternative histones might play role in the establishment of transcriptional memory. Nucleosomes containing alternative histones around promoter regions are a characteristic of inducible genes 10 . It is believed that these nucleosomes provide higher flexibility to the chromatin structure around promoters of inducible genes allowing a quicker and stronger expression upon stimulus. Our findings that H4r is abundant at promoters of genes that are inducible or associated with developmental processes and gets incorporated to promoters of Hsp genes after recovery from heat stress give rise to the assumption that H4r might play a role in establishing a transcriptional memory. Incorporation of H4r to the promoters of inducible genes might also be part of a priming for easier transcription.
H4r loss of function mutants show reduced viability whereas the fertility and longevity of mutants is not affected. The lack of H4r leads to minor changes in the transcriptional pattern, mostly upregulating genes 5 . Our results of H4r localization on genes associated with developmental processes and genes implicated in response to various stimuli, might be interpreted as sign of a possible role of the H4r gene to provide a H4 pool that can be involved in dynamic alteration of chromatin structure after transcriptional changes that emerges with differentiation and environmental stimuli. In differentiating cells, upon entry into G0 phase, genes associated with the differentiated state are transcribed, whereas genes that drive the cell cycle remain silenced, packed by canonical histones. In contrast, genes associated with differentiation are packed in nucleosomes by histone variants, which have replication-independent expression and provide easier access for further transcription. In the absence of H4r, maintenance of the normal transcription pattern in differentiated cells could be impaired, leading to developmental defects and consequently reduced viability.

Materials and methods
Generation of transgenic fly lines. We used CRISPR/Cas9 system for generation of flies expressing H4r fused with a 3xFlag tag. We followed the protocol described by Port et al., 2014 11 andHenn et al., 2020 12 . For the guide-RNA sequences cloned into pCFD4 plasmid see Supplementary Table S2.
For modifying the H4r locus we amplified an extended H4r genomic region (primers used: H4r extended genomic region Fw and Rev, see Supplementary Table S2), and ligated the NheI-digested amplicon to pBlueScript II KS ( +) plasmid, creating donor plasmids for the CRIPSR/Cas9 system mediated homologous recombination. For tagging H4r, we created a subclone by amplifying the gene using H4r subcloning Fw and Rev primers, and cloning it to a modified pBlueScript II KS ( +) plasmid via BglII and KpnI. We fused 3xFlag-tag to the H4r gene via PCR using N-terminal Flag-tag Fw and Rev primers, generating a daughter plasmid. We made the generated daughter plasmid circularised by digesting both end with ClaI, and ligated the digested ends by T4 ligase. We cloned the tagged H4r gene to the donor plasmid using HpaI and BglII restriction endonuclases, generating a daughter donor plasmids. Then a new subclone plasmids was generated by cloning the modified H4r gene and www.nature.com/scientificreports/ the adjacent 265 bp upstream and 430 bp downstream sequences using SacII and PstI enzymes. On this subclone, the PAM sequences were mutated by Sequence and Ligation Independent Cloning (SLIC) 13 : the amplification of H4r gene surrounded by mutated PAM sequences was made by the PAM mutation insert Fw and Rev primers, the plasmid sequences bounded by the mutated PAM sequences were amplified with the PAM mutation plasmid Fw and Rev. The loxP-dsRed-loxP (derived from pHD-dsRed plasmid) sequence was cloned to the NdeI recognition sequence of the plasmids after blunting. The modified H4r gene and dsRed marker gene surrounded by loxP sequences were cloned to the donor plasmid by using SacII and PstI enzymes. For the sequences of the primers used see Supplementary Table S2. An outline of the cloning steps for creating donor plasmids is shown on Supplementary Figure S1. For generation of a line expressing 3xFlag-tagged H4r and dsRed marker, we injected the donor plasmids as it is described by Henn et al., 2020 12 .
For immunohistochemical staining we created a derivative of the above described line expressing 3xFlag-H4r and dsRed. For knocking out dsRed in order to avoid high background on immunohistochemistry samples, we used Cre-mediated recombination of loxP sites.
Transgenic line carrying UAS-H3.3-3xFlag gene (w; UAS-H3.3-3xFlag) 14 was crossed with elav-Gal4 (BL458) for ChIP-seq experiments. Chromatin samples were prepared from the heads of adult offspring of this crossing. www.nature.com/scientificreports/ Immunohistochemistry and western blotting. Immunostaining of embryos was performed as described by Henn et al., 12 . Larval and adult brains were dissected in Ringer's solution and then fixed in PBS containing 4% formaldehyde by rotating for 20 min at room temperature. After removing the fixative solution, samples were washed three times for 5 min at room temperature in PBT (0.1% Triton-X-100 in PBS). They were blocked in PBT-N (0.1% Triton-X-100, 1% BSA and 5% FBS in PBS) at room temperature for 1 h. Samples were incubated overnight in PBT-N containing primary antibodies, then they were washed three times at room temperature in PBT for 10 min, then incubated in PBT-N containing secondary antibodies and DAPI for 1 h at room temperature. After another three washes with 10 min of PBT, samples were placed on microscope slide and mounted in Fluoromount-G (Invitrogen). For positive and negative controls of the immunohistochemical experiments see Supplementary Figure S3.
To prepare complete protein extracts for western blot 15 adult heads were homogenized in RIPA buffer (150.4 mM NaCl; 0.1% SDS; 0.5% Na-DOC; 0.01% Triton-X-100; 1 × PIC; 0.025 M Tris-HCl pH 7.5) (4 μl/ adult head) and then incubated on ice for 30 min. At the end of incubation time, samples were centrifuged at 13,000 rpm for 10 min, and the supernatants were transferred to clean tubes. Pellets were resuspended in an amount of RIPA buffer equal to the amount of supernatants, and the supernatants containing the soluble proteins and the resuspended pellets containing the chromatin-bound proteins were boiled with SDS Loading buffer for 10 min.
Proteins were separated via Tricin-SDS-PAGE 15 . Prior to blotting, membranes were washed in methanol for 15 s, then in water for 2 min, then in blotting buffer (0.02 M Tris-HCl pH 8.0; 0.15 M glycine; 20% methanol). Following blotting, membranes were blocked in TBST (10 mM Tris-HCl pH 8.0; 150 mM NaCl; 0.05% Tween-20) containing 5% non-fat dry milk for 1 h at room temperature, then membranes were incubated in TBST containing 0.02% BSA and primary antibodies at 4 °C overnight. After removing TBST with BSA and primary antibodies, membranes were washed in TBST four times for 10 min and then they were incubated in TBST containing 0.02% BSA and secondary antibodies for 1 h at room temperature. After washing them four times for 10 min in TBST, membranes were incubated in 10 × diluted ECL reagent (Millipore) at room temperature for 5 min then signals were detected via Li-Cor C-DiGit Scanner and measured with ImageJ software. For statistical analysis one-way ANOVA with Tukey's multiple comparison test was performed via GraphPad Prism 8.0.1.
RNA extraction and qPCR. RNA was extracted from adult heads homogenized in NE buffer (15 mM HEPES pH 7.6, 10 mM KCl, 5 mM MgCl 2 , 0.1 mM EDTA, 0.5 mM EGTA, 350 mM sucrose, 0.1% Tween 20, 1 mM DTT, 1 × PIC (Proteinase Inhibitor Cocktail (Roche))). RNA extraction was performed by using TRIzol Reagent (Thermo) according the recommendations of the manufacturer. DNA contamination in the extracted RNA samples were removed with DNaseI (Thermo), and reverse transcription was performed via TaqMan Reverse Transcription Reagents (Invitrogen). qPCR was performed by using Promega GoTaq qPCR Master Mix (Thermo). The primers used for qPCR are shown in Supplementary Table S2. For statistical analysis one-way ANOVA with Sidak's multiple comparison test was performed via GraphPad Prism 8.0.1.

Chromatin-immunoprecipitation and sequencing.
Chromatin-immunoprecipitation was performed as described by Schauer et al., 2013 16 , with the following modifications: 0.954-2.575 µg chromatin was used to each ChIP assay. Dynabeads protein G (Invitrogen) matrix was equilibrated in RIPA buffer without carrier DNA and BSA. After immunopurification, beads were washed five times with RIPA buffer, once with LiCl buffer and once with 10 mM Tris-HCl pH 8.0 (without EDTA), 5 min ea. Beads were resuspended in 10 mM Tris-HCl pH 8.0. Samples were incubated with Proteinase K (Serva) for 3 h at 50 °C after removing RNA contamination and reverse crosslinking. Immunoprecipitated DNA was purified by phenol-chloroform-isoamilalcohol extraction followed by precipitation with ethanol. DNA samples were resuspended in 10 mM Tris-HCl pH 8.0.
Library preparation was performed by using NEBNext Ultra II DNA Library Prep Kit for Illumina according the recommendations of the manufacturer, without size selection. MiSeq Reagent Kit v3 was used for sequencing. For each sample, two biological replicates were sequenced two times.
Since total input controls were not sequenced, peaks and reads of variant histones were compared with normalization to H3 reads of the given samples. Genomic annotations were performed on Drosophila melanogaster assembly release BDGP6.28.101. Read numbers were determined using BedCov Galaxy Version 2.0.2 24 , and read numbers were normalized following this formula: number of reads on a given region in a given sample / total number of reads in a given sample. For determining read counts on entire genes, annotation file was downloaded from FlyMine (https:// www. flymi ne. org/ flymi ne/ regio ns). For gene ontology annotations PANTHER GO-Slim (http:// panth erdb. org/) with Fisher's Exact test type and FDR correction was used. Expressional data for counting reads on highly, inducible or weakly/not expressed genes derive from FlyBase expression profile (http:// flyba se. org/ rnaseq/ profi le_ search): highly expressed genes are determined as genes showing peak expression levels not less than 'very high' in adults in the age of 1-5 days and peak expression levels not less than 'moderately high' in adult heads. Weakly or not expressed genes are determined as genes showing peak expression levels not higher than 'low' in adults in the age of 1-5 days and peak expression levels not higher than 'no/extremely low' in adult heads. Inducible genes are determined as genes showing peak expression levels not less than 'very high' in adults in the age of 1-5 days and peak expression level not less than 'moderately high' in adult heads in case of any treatment with available expressional data. Statistical analysis of the comparisons for the genome-wide distribution of H4r, H3.3 and H3 was performed using two-way ANOVA with Tukey's multiple comparison test was performed via GraphPad Prism 8.0.1. For the statistical analysis of H4r and H3.3 abundance relative to H3 on inducible, highly and weakly/not expressed genes Kruskal-Wallis test followed with Dunn's multiple comparison test was used via GraphPad Prism 8.0.1. For statistical analysis of the changes in the absolute and relative abundance of H4r and H3.3 on the Hsp genes in the distinct conditions Friedman-test followed with Dunn's multiple comparison test was used via GraphPad Prism 8.0.1.

Data availability
Datasets generated and used in this study are available on the National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) under accession GSE197256.