HP1B is a euchromatic Drosophila HP1 homolog with links to metabolism

Heterochromatin Protein 1 (HP1) proteins are an important family of chromosomal proteins conserved among all major eukaryotic lineages. While HP1 proteins are best known for their role in heterochromatin, many HP1 proteins function in euchromatin as well. As a group, HP1 proteins carry out diverse functions, playing roles in the regulation of gene expression, genome stability, chromatin structure, and DNA repair. While the heterochromatic HP1 proteins are well studied, our knowledge of HP1 proteins with euchromatic distribution is lagging behind. We have created the first mutations in HP1B, a Drosophila HP1 protein with euchromatic function, and the Drosophila homolog most closely related to mammalian HP1α, HP1β, and HP1γ. We find that HP1B is a non-essential protein in Drosophila, with mutations affecting fertility and animal activity levels. In addition, animals lacking HP1B show altered food intake and higher body fat levels. Gene expression analysis of animals lacking HP1B demonstrates that genes with functions in various metabolic processes are affected primarily by HP1B loss. Our findings suggest that there is a link between the chromatin protein HP1B and the regulation of metabolism.


Introduction
Regulation of gene expression is an essential process. Errors in gene regulation lead to aberrant transcription or silencing, which in turn impact the fitness, health, and survival of the organism. In eukaryotes, chromatin, the composite of DNA and its associated proteins, mediates important aspects of gene regulation. Specific chromatin types are formed depending on the positioning of nucleosomes, the post-translational modification of the histones within the nucleosome, and the types of chromosomal proteins present. The simplest models distinguish two basic chromatin types, heterochromatin and euchromatin, but more elaborate models exist as well [1][2][3]. Generally speaking, euchromatin encompasses gene-rich, transcriptionally active sequences, while heterochromatin contains highly condensed, repeat-rich, transcriptionally inactive regions [4]. These chromatin types are biochemically distinct, and chromosomal proteins carry out a variety of functions in specifying the distinct chromatin structures. Some are histone modifying enzymes, nucleosome remodelers, or boundary proteins, while PLOS

Fly culture
All fly lines were cultivated on standard cornmeal-sucrose media [29] or Jazzmix Drosophila media (Fisher Scientific) at 25˚C and 70% humidity. Animals other than those used for the fertility and PEV assays, were raised under controlled 12h light/ 12h dark cycle. The starting stock to generate the new HP1b alleles carrying the P{PTT-GB}HP1b CB02849 insertion was obtained from the Carnegie Protein Trap Library [30]. The insertion is located in the 5' UTR of HP1b.

Generation of HP1b mutant alleles
The P element of HP1b CB02849 was mobilized by crossing homozygous females (HP1b CB02849-mw , y -/HP1b CB02849-mw , y -) to males containing the Δ2-3 P element transposase (w � ; ry 506 Sb 1 P{Δ2-3}99B/TM6B, Tb 1 ; Bloomington Drosophila Stock Center #1798). From the F1 population, males with the HP1b CB02849 allele and the Δ2-3 P element transposase were selected (by the Sb 1 phenotype) and crossed to balancer females (FM7j, B 1 , y -, w -/Df (1) ED6957 HP1b , y + , w mw ; derived from Bloomington Drosophila Stock Center #8033). The HP1b CB02849 transgene contains a copy of mini-white (mw), resulting in red eye pigmentation (mw). Thus, putative HP1b mutants with the FM7j balancer (HP1b � , y -/FM7j, B 1 , y -, w -) among the progeny from the cross were identified based on the lack of red eye color and the presence of a yphenotype, indicating the removal of the P{PTT-GB}HP1b CB02849 insertion. Putative mutants also were required to lack the Δ2-3 P element transposase by selecting against the Sb 1 phenotype. From 130 single male crosses set up in step 2, 89 female offspring lacking red eye color and with a yphenotype were identified. Stocks were established, and males (HP1b � /Y) were screened by PCR (5'-CTA GCG CAC GTG CAA CGC AAC-3'; 5'-ATA AAG GCG CTA GGT GGG CAG C -3') to identify deletions. In addition to many small insertions, we identified three deletion alleles: HP1b 16 (Δ539bp), HP1b 86 (Δ1078bp), and HP1b 49 (Δ1309bp). Prior to carrying out the analyses described below, the HP1b 16 (Δ539bp) and HP1b 86 (Δ1078bp) alleles that are characterized in detail here were backcrossed into to a yw -(yw) control strain for 6 generations to minimize genetic background effects.

Fertility assay-Female
Single virgin females (0-5 days old) from the yw, HP1b 16 and HP1b 86 strains were mated to two yw males in individual vials. After three days, these flies were moved to new vials, and again to a third set of vials on day six after the initial mating was set up. At day nine, the flies were discarded. From each of the three vials thus obtained from each female, offspring were counted four times, on days 10, 12, 14, and 16 after the vial was initially set. The offspring number per female was calculated by combining all offspring from the individual counts (12 total). Differences in offspring number between genotypes were analyzed by analysis of variance (ANOVA) and t-tests using R [31].

Fertility assay-Male
Single virgin males (0-5 days old) from the yw, HP1b 16 and HP1b 86 strains were mated to two yw females in individual vials. After three days, these flies were moved to new vials, and again to a third set of vials on day six after the initial mating was set up. At day nine, the flies were discarded. From each of the three vials thus obtained from each male, offspring were counted four times, on days 10, 12, 14, and 16 after the vial was initially set. The offspring number per male was calculated by combining all offspring from the individual counts (12 total). Differences in offspring number between genotypes were analyzed by ANOVA and t-tests using R [31].

Position effect variegation assay-White variegation
Eye pigment levels were measured by crossing males containing the P element reporter stock males to yw, HP1b 16 , or HP1b 86 homozygous virgin females. Offspring from the cross were aged to 3-5 days, and acid-extracted eye pigment was measured at OD 480 for male and female offspring separately [adapted from [32]]. For each assay, five flies were used, with 6-10 biological replicates per genotype. Due to the location of HP1b on the X chromosome, male offspring from this cross lack wildtype HP1b, while females have one dose. Differences in eye pigment levels between genotypes were analyzed by ANOVA and t-tests using R [31].

Position effect variegation assay-Stubble variegation
To assay the impact of HP1b mutations on Stubble variegation, the Sb V allele was used, similar to the assay described in Csink et al 1994 [33]. Offspring lacking the serrated wing phenotype associated with Ser 1 were examined, and the number of wildtype bristles was scored. In total, ten bristles per animal where assayed (four scutellar bristles, four dorsocentral bristles, and two posterior post-alar bristles), and the number of wildtype bristles was recorded for each animal. Due to the location of HP1b on the X chromosome, male offspring from this cross have lack wildtype HP1b, while females have one dose. Differences in bristle number between genotypes were analyzed by ANOVA, Tukey's HSD post hoc tests, and t-tests using R [31].

Animal activity-Third instar larvae
Sets of ten third instar larvae were placed in the center of a petri dish containing 1.5% agar. The petri dish was placed onto a diagram of concentric circles spaced 1cm apart to evaluate the movement of the larvae. The larvae were allowed to move freely, and pictures were taken with a digital camera every 15 seconds for 2 minutes. From the images, the number of animals in each concentric circle were recorded to provide an estimate of the distance traveled by the animals. For each genotype (yw, HP1b 16 , or HP1b 86 ), five trials with eight larvae each were performed. The data from the 60 second time point were recoded to divide all animals into two categories: those animals that remained in the center of the petri dish ("no movement") and those animals that left the center of the petri dish ("movement"). The fraction of animals moving versus not moving was compared between the HP1b mutant and yw control genotypes (chi-square test) using R [31].

Animal activity-Adult flies
Basal activity of the animals was monitored using a locomotion activity monitor (LAM25) from TriKinetics. 3-5 day old virgin flies from the yw, HP1b 16 , and HP1b 86 strains were placed into assay vials (n = 5 for each sex and genotype; ten flies per vial). Next, the vials were placed into the monitor, and the monitor was placed in the incubator (25˚C, 60% humidity, 12h light-dark cycle). The flies were allowed to recover for 1h prior to the start of the assay period. Animal activity was recorded every 5 minutes from 8AM to 10AM. We performed five independent replicates of this assay. Total activity during the 2hr period was used in the analyses. ANOVA (factors: sex, genotype, trial) and Tukey's HSD post hoc tests were performed in R to investigate the sources of the variation observed between the samples [31].

Stress assay-Starvation
3-5 day old virgin flies were moved onto starvation media (1.5ml of 1% agar to provide hydration) and placed into an incubator (25˚C, 60% humidity, 12 hr light-dark cycle). Deaths were recorded every 8 hours until all flies had died (n = 100 for each sex and genotype [yw, HP1b 16 , HP1b 86 ]; 20 flies per vial) [34]. Three independent replicates of this assay were performed. Survival analysis was carried out in R [31], using the "Survival" package, specifically its survival difference formula (Kaplan-Meier), to compare the survival curves between the three genotypes [31,35].

Stress assay-Heat stress
3-5 day old virgin flies were transferred to fresh Jazzmix medium (Fisher Scientific) and placed in a warm room maintained at 37˚C. Flies were monitored every hour until all flies had died (n = 100 for each sex and genotype; 20 flies per vial) [36]. We performed three independent replicates of this assay. Survival analysis was carried out in R [31], using the "Survival" package, specifically its survival difference formula (Kaplan-Meier) [31,35].

Stress assay-Oxidative stress
To assay oxidative stress, flies were exposed to the herbicide paraquat, an oxidizer. Filter paper disks saturated with paraquat solution (20mM paraquat, 5% sucrose) were placed into vials containing starvation media (1.5ml of 1% agar). 3-5 day old virgin flies were placed inside vials (n = 100 for each sex and genotype; 20 flies per vial), returned to the incubator (25˚C, 60% humidity, 12 hr light-dark cycle), and monitored every 4 hours until all flies had died [37]. We performed three independent replicates of this assay. Survival analysis was carried out in R [31], using the "Survival" package, specifically its survival difference formula (Kaplan-Meier) [31,35]. A second oxidative stress assay was carried out to determine if the increased survival of HP1b mutants was related to their lower food and thus paraquat intake. In this assay, the paraquat concentration was scaled based on food consumption; thus, the paraquat concentration for yw was decreased to 13.2mM paraquat. We performed this assay with 100 animals for each sex and genotype, again with 20 flies per vial. Data analysis was performed in R [31] using the R package "Survival," specifically, its survival difference formula (Kaplan-Meier) [31,35].

Lifespan assay
To determine lifespan, virgin flies were maintained on Jazzmix (Fisher Scientific) under standard conditions (25˚C, 60% humidity, 12 h light-dark cycle). Flies were moved to new food every three days without the use of CO 2 and monitored daily until all flies had died (starting population: n = 100 for each sex and genotype; 20 flies per vial) [36,38]. We performed three independent replicates of this assay. Survival analysis was carried out in R [31] using the "Survival" package, specifically its survival difference formula (Kaplan-Meier) [31,35].

Food consumption assay
Capillary Feeder (CAFÉ) assays were carried out as previously described [39]. Briefly, for the assay, animals were placed on starvation media (1.5ml of 1% agar to provide water) in standard fly vials that were cut in half, with four equally spaced holes around the perimeter. A nutrient solution (5% sucrose, 5% yeast solution) was provided through two capillaries inserted into the foam vial top. At 4PM of the day prior to the start of the assay, 3-5 day old virgin flies were placed into the assay vials (five flies per vial; n = 4-8 for each sex and genotype), the capillaries were filled with the nutrient solution, and the flies were returned to the incubator (25˚C, 60% humidity, 12 hr light-dark cycle) to acclimatize overnight to the new feeding method. At 8AM, the capillaries were refilled with the nutrient solution. Eight hours later, the drop in liquid levels within the capillaries was measured. We assayed the same flies twice over two eight-hour time periods (8AM to 4PM) on two consecutive days, and used average food consumption for further analysis. Evaporation was estimated using control vials lacking flies. We performed three independent replicates of this assay. Due to differences in the capillaries used, the three trials were analyzed separately. Data analysis was performed in R using ANOVA with a Tukey's honest significant difference (HSD) post hoc test and Kruskal-Wallis tests in case of non-normality [31].

Body composition analysis
Body fat content was measured using quantitative magnetic resonance (QMR) for flies at specific ages at the UAB Small Animal Phenotyping Core. Virgin flies were collected, sorted by sex and genotype, and aged for the desired time (3-5, 20, 30, 40 or 50 days old). At the appropriate age, flies were placed into a microcentrifuge tube (10 flies per tube, 4 tubes per sex and genotype for each time point), and stored on ice. QMR analysis was performed by the UAB Small Animal Phenotyping core with an EchoMRI 3-in-1 machine (Echo Medical Systems, Houston, TX). Flies were transferred into the biopsy tube and scanned using the biopsy setting with 9 primary accumulations. Data were analyzed in R using ANOVA with a Tukey's HSD post hoc test [31].

RNA-seq analysis
RNA was isolated using Trizol from third instar larvae of the HP1b 16 and HP1b 86 mutant stocks as well as of the yw control strain. Integrity of the RNA samples was confirmed by formaldehyde agarose gel electrophoresis. The RNA samples were prepared for RNA-seq analysis by the UAB Heflin Center for Genomic Science Genomics Core Lab. 30-40 million RNA-seq reads were collected per sample using the Illumina Sequencing Platform, and two samples per genotype were processed. The UAB Heflin Center for Genomic Science Genomics core processed the raw data using the Tuxedo suite software (Bowtie, Tophat, Cufflinks) and identified genes differentially expressed between the HP1b mutants and the yw control strain by Cuffdiff [40]. Only genes meeting an FDR (false discovery rate) cut-off of 0.05 were considered in the downstream analyses.

GO term enrichment analysis
The RNA-seq analysis revealed that while neither HP1b allele generated full-length mRNA or functional HP1B protein, an aberrant RNA was produced from the HP1b locus in the HP1b 16 stock. To minimize the potential impact of this aberrant mRNA on the interpretation of results, only genes misregulated in both HP1b mutant strains were included in the GO term analysis. Differentially expressed genes were classified as up-versus down-regulated, and genes misregulated in both HP1b mutant strains were identified using custom scripts. GO term enrichment analysis was carried out using PANTHER [41].

Krebs cycle metabolite analysis
Metabolite concentrations were measured by mass spectrometry for virgin flies aged 3-5 days. Flies were sorted by sex and genotype, aged 3-5 days, flash frozen in liquid nitrogen, homogenized in cold methanol, and delivered to the UAB Targeted Metabolomics and Proteomics Laboratory (TMPL) for analysis (n = 4 biological replicates per sex and genotype). At the TMPL, supernatants were stored at -80˚C until further processing. Standards were generated as a master mix of all compounds at 100 μg/mL in water and serial diluted to 10x of the final concentrations (0.05-10 μg/ml, 9 standards). Standards were further diluted to 1x in methanol to a total volume of 1 mL, and dried by a gentle stream of N 2 . For cell extracts, 1 mL of each were transferred to a glass tube and dried under a gentle stream of N 2 . Standards and samples were resuspended in 50μl of 5% acetic acid and vortexed for 15 seconds. Amplifex Keto Reagent (SCIEX, Concord, Ontario, Canada) (50 μL) was added to each sample and allowed to react for 1h at room temperature. Standards and samples were then dried under a gentle stream of N 2 and resuspended in 1ml of 0.1% formic acid.
Samples were analyzed by liquid-chromatography(LC)-multiple reaction ion monitoringmass spectrometry. Liquid chromatography was performed using a LC20AC HPLC system (Shimadzu, Columbia, MD) with a Synergi Hydro-RP 4 μm 80A 250 x 2 mm ID column (Phenomenex, Torrance, CA). Mobile phases were: A) 0.1% formic acid and B) methanol/0.1% formic acid. Compounds were eluted using a 5-40% linear gradient of B from 1 to 7 min, followed by a column wash of 40-100% B from 7 to 10 min, and re-equilibrated at 5% B from 10.5-15 min. Column eluant was passed into an electrospray ionization interface of an API 4000 triple-quadrupole mass spectrometer (SCIEX). The following mass transitions were monitored in the positive ion mode: m/z 261/118 for α-ketoglutarate, m/z 247/144 for oxaloacetate and m/z 204/144, 204/118 and 204/60 for pyruvate. In the negative mode, the following Peak areas of metabolites in the sample extracts were compared in MultiQuant software (SCIEX) to the those of the known standards to calculate metabolite concentrations. MetaboAnalyst 3.0 was used for the analysis [42], with additional data analysis performed in R [31].

Mitochondrial complex analysis
Mitochondrial complex 3 and citrate synthase activity were measured utilizing the resources at the UAB Bio-analytical Redox Biology (BARB) core. Thoraces were dissected from 3-5 day old male animals (five biological replicates per genotype) and placed in isolation buffer on ice [43]. The thoraces were homogenized, remaining solids removed by filtration, and protein concentration was measured with a Lowry assay [44]. To assay citrate synthase, a coupled reaction with acetyl-CoA, oxaloacetate, and 5,5-dithiobis-2-4 nitrobenzoic acid was used, and the results were measured using a spectrophotometer [45,46]. Complex III activity was measured by assaying the conversion of cytochrome C (its electron acceptor) from the oxidized to the reduced form using a spectrophotometer to detect a shift in absorbance [47]. Data analysis was performed in R [31].

HP1b is a non-essential gene
We generated loss-of-function mutations in the HP1 family gene HP1b by imprecise P element excision to study its function. Our screen generated three deletion alleles, HP1b 16 , HP1b 86 , and HP1b 49 (Fig 1A). The deletions span 539bp, 1,078bp, and 1,309bp respectively, with HP1b 16 and HP1b 86 deleting major portions of the HP1b coding sequence, and the HP1b 49 deletion extending into the intergenic region between HP1b and its downstream neighbor APC4 (Anaphase Promoting Complex 4). HP1b 16 and HP1b 86 are homozygous viable, while HP1b 49 is homozygous lethal early in development (no homozygous first instar larvae are observed). As HP1b 16 and HP1b 86 both lack full-length HP1b mRNA and protein ( Fig 1B and 1C), we attribute the lethality of HP1b 49 to an effect on APC4 due to the deletion of some of the regulatory sequences presumably present in the intergenic sequences between HP1b and APC4. While the full-length HP1b mRNA is absent in HP1b 16 and HP1b 86 mutant embryos, the two neighboring genes APC4 and CG7426 are expressed ( Fig 1B). Thus, our analysis here focuses on the HP1b 16 and HP1b 86 alleles. At 25˚C, homozygous stocks of HP1b 16 and HP1b 86 are viable and can be maintained without difficulties under standard conditions. Thus, we conclude that HP1b is a non-essential gene in contrast to the centromeric HP1 family member, HP1a, encoded by Su(var)205, which is homozygous lethal at the third instar stage [24,48,49].

Mutations in HP1b impact the number of offspring produced by females
Expression data generated by modENCODE available from Flybase (www.flybase.org, [50]) indicate that expression of HP1b is moderate throughout development with high expression during early embryonic development (2-4hr) and in the ovaries [51,52], suggesting that HP1b might impact the number of offspring produced per female. To investigate this hypothesis, we assayed the potential impact of HP1b mutations on offspring number. When the average offspring number of single females mated to two males of yw -(yw) is compared between HP1b mutants and control yw females, a reduction in offspring number is observed for the mutants (Fig 2A). While the yw females had on average 56 offspring, HP1b 16 and HP1b 86 had 31 and 36 offspring on average respectively, a drop that is statistically significant (p<0.01 for HP1b 16 and p<0.05 for HP1b 86 ; t-test). As a control, we also measured the average number of offspring produced by individual HP1b and yw males mated to two yw females. We found no significant differences in male offspring number between the HP1b mutant animals and the yw control strains (Fig 2B, p-value not significant). These data indicate that the lack of HP1B protein specifically impacts the offspring number of female animals.

HP1b is an enhancer of variegation
Position effect variegation (PEV) analysis in Drosophila is a powerful tool to investigate the impact of mutations on transgene expression in particular chromatin environments. Transgene expression from reporters within or close to heterochromatin tends to variegate due to silencing of the reporter in some cells but not others. If a mutation for example affects the chromatin structure, increased expression or suppression of variegation [Su(var) phenotype] is seen if the protein is required for heterochromatin formation or has repressive functions. If on the other hand decreased expression or enhancement of variegation is seen [E(var) phenotype], this result indicates that the protein is involved in counteracting heterochromatin formation or has activating functions. This assay has been used successfully to demonstrate that mutations in HP1a increase expression of a variety of reporter genes [Su(var) phenotype] [48], indicating it can act as a transcriptional repressor. In contrast, tethering HP1C to a reporter gene increased its expression, indicating that it functions as a transcriptional activator [53]. We used PEV analysis to investigate the role of HP1B in the regulation of gene expression. Males containing a P element reporter that fuses the hsp70 promoter to the white coding region were crossed to virgin females from the homozygous HP1b 16 and HP1b 86 stocks as well as to yw as a control. As PEV effects often depend on the specific reporter insertion used, five different reporter lines were tested: 39C-2 (chr. 2; pericentric), 39C-3 (chr. 2, pericentric), 118E-10 (chr. 4; pericentric), 39C-12 (chr. 4; arm), and 39C-27 (chr. 2R, telomeric) [54]. Assaying males, which lack HP1B entirely due to the gene's location on the X chromosome, we find that eye pigment levels are reduced in the presence of the HP1b mutations compared to the yw control ( Fig 3A). This E(var) phenotype is consistent for all five reporters assayed, including pericentric and telomeric reporter insertions (Fig 3B). In females, which are heterozygous for the HP1b mutations, the results differ depending on the reporter location and include E(var) responses (39C-2 and 39C-3) as well as Su(var) responses (118E-10 and 39C-12) (S1 Fig). Thus, HP1b mutations act as modifiers of PEV in both sexes for these hsp70-white reporters, and they have a strong E(var) effect in males in all genomic contexts assayed.
To ensure that the E(var) effect upon complete loss of HP1B is a general phenomenon and not restricted to the white reporter used in these assays, we investigated a second PEV system, variegation of the Stubble (Sb) allele Sb V . Sb V is a gain-of-function allele that causes a bristle defect. Animals heterozygous for Sb V have a mixture of wildtype and defective bristles. If an E Representative eye pictures are shown from males for two of the reporter insertions, 39C-3 located in the pericentric region of chromosome 2, and 118E-10 located in the pericentric region of chromosome 4. B. Quantitative eye pigment assays for males from the five reporter insertions assayed. All comparisons of mutant eye pigment values to the yw control are statistically significant (p<0.05, t-test). X-axis: reporter insertions. Y-axis: OD 480 measuring eye pigment. Results from the yw control are in red, from HP1b 16 in blue, and from HP1b 86 in green. Box plots: black bar-median; box-+25% and -25% quartiles; whiskers-maximum and minimum; circles-outliers; n = 6-10. C. Stubble variegation assays indicate that HP1b is an enhancer of variegation, as the presence of the HP1b mutant alleles significantly increases the number of wildtype bristles compared to the yw control (p<0.001, Tukey multiple comparisons of means; n = 19-39). Data from males. Results from the yw control are in red, from HP1b 16 in blue, and from HP1b 86 in green. Box plots: black bar-median; box-+25% and -25% quartiles; whiskers-maximum and minimum; circles-outliers. For data from females, see S1 (var) mutation is crossed to the Sb V reporter, the fraction of defective bristles will decrease, while in the presence of a Su(var) mutation, the fraction of defective bristles will increase. Males containing the Sb V reporter were crossed to virgin females from the homozygous HP1b 16 and HP1b 86 stocks as well as to yw as a control. In the offspring, ten bristles (four scutellar bristles, four dorsocentral bristles, and two posterior post-alar bristles) were examined, and the number of wildtype bristles was recorded. We find that the presence of HP1b null alleles significantly increases the number of wildtype bristles compared to the yw control (Fig 3C; Tukey's test, p<0.01). In males, for example, the wildtype bristle number increased approximately 2-fold from an average of 2.7 wildtype bristles per animal in yw to 5.3 and 6.1 in the HP1b 16 and HP1b 86 animals respectively. HP1b mutations also increase the wildtype bristle number in females (S1 Fig), and thus, HP1b mutations generally act as enhancers of variegation for the Sb V reporter. These findings confirm that the results obtained with the hsp70-white reporters and suggest that complete loss of HP1B results generally in enhancement of variegation as indicated by the effects in diverse genomic contexts with multiple reporters.

HP1b mutants have decreased activity levels
During the experiments investigating PEV and offspring number, we noticed that the HP1b mutants tended to be easier to handle than the yw controls as they appeared to move less. To determine if mutations in HP1b indeed affected the activity levels of the animals, experiments were carried out in third instar larvae and adult flies. In third instar larvae, we assayed crawling ability and determined the fraction of animals that moved or failed to move in a 60 second time period (Fig 4A). In this assay, HP1b mutant larvae tended to remain unmoving more often. For example, only 21% of the larvae from the HP1b 16 strain moved during the assay Drosophila HP1B has links to metabolism period (10/47), while 48% of the yw larvae moved (23/48). Altogether, significantly fewer HP1b mutant animals moved away from the center of the assay dish during the assay period (p< 0.001 for both mutant alleles, chi-square test, Fig 4A). These data indicate that lack of HP1B leads to lower activity levels in third instar larvae.
Next, activity levels of adult animals were assayed to determine if the lower activity level observed in larvae is a general feature of HP1b mutants consistent across different developmental stages. We utilized an activity monitor which measures fly activity levels by recording every instance a fly crosses the middle of the assay chamber. In this assay, HP1b mutants have significantly decreased levels of activity compared to the controls (Fig 4B). In a two-hour period, an average of 1142 crossings were recorded for HP1b 16 mutants, HP1b 86 mutants crossed the middle of the assay chamber on average 944 times, while the yw control flies logged 1441 crossings (Fig 4B). These differences are highly significant (Tukey's HSD; p<0.001 for both mutant alleles). While activity levels are significantly different between males and females (ANOVA, p<0.001), the decrease in activity associated with mutations in HP1b is observed in both males and females (S2 Fig). However, the decrease is more pronounced in females, where the decrease is significant for both HP1b 16 and HP1b 86 (Tukey's HSD, p<0.001), while in males, the activity of only HP1b 86 mutants is significantly decreased compared to yw (Tukey's HSD, p<0.001). Taken together with the results from the larval crawling assay, these results suggest that lowered animal activity levels are a general feature of HP1b mutants, irrespective of developmental stage.

Loss of HP1B alters stress resistance
Given that loss of HP1B impacts the number of offspring produced by the animals as well as their activity levels, we questioned whether other aspects of the animals' physiology might be impacted as well. First, we investigated stress resistance. In Drosophila, increased stress resistance is often observed in animals with increased lifespan [37,[55][56][57][58], and as stress assays require much less time, they are typically carried out prior to investigating lifespan. Three separate stress resistance assays were used to investigate resistance to starvation stress, oxidative stress, and heat stress. To assay resistance to starvation stress, flies were grown under standard conditions and transferred to vials containing 1% agar, providing hydration but no nutritional value. Deaths were recorded every eight hours, and survival curves were generated (Fig 5A). Survival curves of the HP1b mutant stocks are shifted to the right compared to the yw control, indicating increased stress resistance and longer survival. Living on average 51.65 and 49.37 hours under starvation conditions, HP1b mutants show significantly increased survival compared to yw, which lives an average of 45.60 hours (p = 1.022e-13 and p = 7.188e-06 for HP1b 16 and HP1b 86 , respectively, Kruskal-Wallis rank sum test). In addition, the survival curves are significantly different between HP1b mutant and yw control flies (Kaplan-Meier survival function: p = 1.11e-16 for HP1b 16 , and p = 9.83e-09 for HP1b 86 ). The increased survival under starvation conditions was observed also when the data from males and females were analyzed separately (S3 Fig). In females, both the survival curves and means are different between HP1b mutants and yw (Kruskal-Wallis rank sum test: p = 6.455e-09 and p = 4.672e-07 for HP1b 16 and HP1b 86 respectively; Kaplan-Meier survival function: p = 2.36e-10 for HP1b 16 , and p = 4.67e-08 for HP1b 86 ), while in males, the comparison of means is only significant for HP1b 16 (Kruskal-Wallis rank sum test: p = 2.054e-07), despite significantly shifted survival curves for both HP1b strains (Kaplan-Meier survival function: p = 1.86e-09 for HP1b 16 and p = 0.011 HP1b 86 ). Together, these data demonstrate that loss of HP1B allows flies to survive longer under starvation conditions. Next, we investigated resistance to oxidative stress using a paraquat feeding assay. Paraquat is an herbicide that induces oxidative stress when ingested, and it is widely used to measure oxidative stress resistance in flies [37]. Flies were fed a paraquat-laced sucrose solution and monitored to record time of death. Results from this assay revealed that the survival curve of the HP1b mutants were shifted to the right, and the control flies (red) showed a steeper dropoff in survival rate than the two HP1b mutant strains (blue and green, Fig 5B). In the combined analysis of males and females, HP1b mutants have significantly increased average lifespan compared to the yw control strain: HP1b 16 16 and p = 0 for HP1b 86 ; males: p = 6.96e-09 for HP1b 16 and p = 1.29e-14 for HP1b 86 ). Together, these data show that loss of HP1B increases oxidative stress resistance as measured by paraquat feeding assays.
Lastly, resistance to heat stress was investigated. Flies were raised under standard conditions and then moved to 37˚C. The flies' time of death was recorded hourly to generate survival curves. In contrast to the findings for starvation and oxidative stress, no increased resistance to heat stress was detected in the HP1b mutants, neither in males nor in females (Fig 5C and S3  Fig). In contrast, the HP1b strains actually showed decreased survival compared to the control, with HP1b 16 surviving on average 20.04 hours, HP1b 86 surviving 18.01 hours, and yw surviving 20.80 hours (p = 0.04158 for HP1b 16 and p = 1.874e-09 for HP1b 86 ; Kruskal-Wallis rank sum test). When males and females are analyzed separately, only the HP1b 86 strain shows significantly decreased survival, while HP1b 16 is not different from yw (females: p = 0.001803; males: p = 1.499e-08; Kruskal-Wallis rank sum test). Analysis of the survival curves shows similar results: only HP1b 86 has a survival curve that is significantly different from yw, be that in the combined analysis or in males/females (combined: p = 2.93e-09; females: p = 0.0002: males: p = 7.02e-7; Kaplan-Meier survival function). These results indicate that HP1b mutant flies do not have increased heat stress resistance. This result is unexpected, because most often, stress resistance in general is increased, irrespective of assay. Thus, HP1b mutants are unusual in showing resistance to starvation stress and oxidative stress but not showing resistance to heat stress.

Loss of HP1B impacts survival
In Drosophila, stress resistance phenotypes are often seen in animals that show increased longevity. Given the unusual stress resistance phenotypes we observed in the HP1b mutant strains, next, we investigated their longevity by generating survival curves. In three replicate experiments, we followed cohorts of 100 adult flies for each genotype/sex from hatching to death (Fig 6). Comparing the survival curves from all three strains, a statistically significant impact of "genotype" as a factor was detected (p = 0.0483, log rank test) when the data from both sexes is combined for the analysis. This difference is mainly driven by the results from females, which show a strong effect of genotype (p = 0.00273, log rank test), which is absent in males (S4 Fig). In the HP1b mutant strains, especially female animals experience less mortality early in life, until approximately day 40, when mortality rates increase. Towards the last quarter of the lifespan (~day 60), mortality rates in the yw control lines are lower than in the HP1b mutant strains, leading to a longer maximum lifespan (p = 1.11e-15 for HP1b 16 and p = 0 for HP1b 86 , log rank test on the last 10% of surviving animals). These results illustrate that loss of HP1b has a complex effect on survival, affecting the two sexes differently, and mainly impacting survival in early/mid-life.

Mutations in HP1b lower food intake
As dietary restrictions can result in altered survival curves [59], we investigated if a change in feeding behavior might contribute to the altered survival curves of HP1b mutants. We measured food intake using the CAFE assay [39], where the animals are fed a nutrient solution through glass capillaries, recording consumption over an 8hr time period. Representative results from one of three trials are shown in Fig 7A. Both genotype and sex affected the animals' food intake strongly in all trials (p<0.05, ANOVA). HP1b mutant animals consume less food than their yw counterparts (Fig 7A), a finding consistent among all three trials despite the fact that individual comparisons are not always significant (S1 Table). Across all trials, HP1b mutants consumed approximately 30% less food than the yw control flies (24-32% depending on trial), and this difference between the HP1b mutants and the yw control is highly significant in a combined analysis of all three trials (p = 0.0000305 for HP1b 16 and p = 0.0013900 for HP1b 86 ). The results from this assay indicate that loss of HP1b impacts food consumption, suggesting that the physiology and behavior of these animals might have been altered.

HP1b mutants have altered body fat levels
The lower activity levels observed in HP1b mutants together with the change in feeding behavior suggests that HP1b mutants might have altered body composition as well. We used quantitative magnetic resonance (QMR) to estimate fat and lean body mass. 3-5 day-old HP1b mutant flies showed significantly increased fat percentage compared to yw (15.5% vs 15.7% versus 10.0%; HP1b 16 and yw p = 0.003, HP1b 86 and yw p = 0.002; Tukey's multiple comparisons of means; Fig 7B). This increase in body fat content is observed in both male and female animals, but in males, only the HP1b 16 /yw comparison reaches statistical significance.
Next, we extended this assay to 50 days, measuring body fat content at 20, 30, 40, and 50 days of age. Specifically, we were interested in determining if the body fat levels remained elevated in the HP1b mutant animals, or if a loss of body fat might be responsible for the steep decline in survival seen for these animals around day 40 (Fig 6). Interestingly, the body fat content of yw animals had a trajectory different from animals lacking HP1B. In yw animals, body fat levels increased with age, from 10.0% at the earliest time point to approximately 18% at 40 and 50 days of age. In contrast, the body fat content of the HP1b mutant animals peaks at day 30 at 19% and drops to 14% at day 50. Thus, the level of body fat at day 50 is significantly lower

Oxidative stress resistance in HP1b mutants is due to altered feeding behavior
Given that our assays indicate that HP1b mutant animals have altered body fat content and altered behavior, both in terms of activity levels and in terms of food intake, we wondered if these changes might explain the increased stress resistance we had observed. Having a higher body fat content and lower energy expenditures due to decreased activity levels could explain the increased starvation resistance of HP1b mutants (Fig 5A), while the lower food intake could result in an overall lower dose of paraquat compared to yw, thus allowing the HP1b mutants to survive longer under oxidative stress (Fig 5B). We tested this hypothesis by carrying out a second oxidative stress assay, this time adjusting the paraquat concentration in the food to account for the reduced feeding rate in HP1b mutants. Compared to the initial assay (Fig 5A), the survival curve of yw shifts to the right, with yw animals now living slightly longer than the HP1b mutant animals (Fig 7C). Animals from the yw line die on average 83.74h after the initiation of oxidative stress, while HP1b 16 animals die on average after 75.08h and HP1b 86 after 74.52h (HP1b 16 : p = 9.693e-07; HP1b 86 : p = 6.548e-07; Kruskal-Wallis rank sum test; Fig 7C). This experiment suggests that the oxidative stress resistance observed for animals lacking HP1B in our initial experiment is due to the decreased food intake level of these animals. Loss of HP1B reduces food intake, alters body fat levels, and leads to oxidative stress resistance due to altered feeding behavior. A. Food intake measured using a CAFÉ assay (Y-axis, decrease in food levels over an 8hr time period as measured by drop in capillary meniscus in mm) reveals lower food intake in HP1b mutants (HP1b 16 : p = 0.0003155; HP1b 86 : p = 0.0095041; Tukey's multiple comparisons of means; n = 8 per sex and genotype). Data from one representative trial, both sexes combined are shown. For results from the additional trials and separated by sex, see S1 Table. B. Percent body fat (Y-axis) was measured by QMR in animals of different ages (X-axis). Body fat levels are increased in mutant animals early in life but decrease at a rapid rate after midlife (30-40 days). Red: yw; blue: HP1b 16 ; green: HP1b 86 . Error bars: SD. N = 4 per sex and genotype. C. After paraquat concentration is adjusted to take into account differences in feeding behavior between HP1b mutants (blue, green) and the yw control strain (red), HP1b mutants no longer exhibit increased resistance to oxidative stress. X-axis: Time to death in hours. Y-axis: Survival probability. For separate survival curves for males and females, see S5 Drosophila HP1B has links to metabolism

Expression analysis reveals misregulation of metabolic genes
To gain insights into the causes of the physiological and behavioral changes observed in the animals lacking HP1B, we performed an RNA-seq experiment to compare gene expression levels between the HP1b mutants and the yw control strain. Using a significance cut-off of p<0.05 (FDR, Benjamini-Hochberg), we find 2179 genes in HP1b 16 and 934 genes in HP1b 86 that differ significantly from yw in their expression levels (S6 Fig). Of these differentially expressed genes, 763 are significantly altered in both HP1b mutant strains (S6 Fig). Overall, if we consider all genes, gene expression changes observed between yw and HP1b mutants occur in roughly equal proportions in either direction (HP1b 16 : 46% up, 54% down; HP1b 86 : 56% up, 44% down; S6 Fig). In contrast, among the genes with significantly altered gene expression, upregulation is more common (S6 Fig). For example, in the gene set with significantly altered gene expression in both HP1b mutant strains, 85% of genes are upregulated, and only 15% are downregulated ( S6 Fig). Together, these data illustrate that loss of HP1B impacts a large number of genes and that while upregulation appears more common, gene expression changes can occur in both directions.
Next, we investigated the types of genes affected by loss of HP1B. By analyzing gene ontology (GO) terms associated with the differentially expressed genes using PANTHER [41], we identified a variety of GO terms enriched in the gene set affected by the loss of HP1B (Fig 8). While none of the GO SLIM terms were enriched significantly in the downregulated gene set, the upregulated genes were highly enriched for a variety of metabolism-related terms. These terms include the broad term "metabolic process", but also more specific terms such as "carbohydrate metabolic process", "lipid metabolic process", and "steroid metabolic process" (Fig 8,  S2 Table). This preponderance of metabolism-related terms is also evident in the full GO analysis (S3 and S4 Tables). Thus, the gene expression analysis suggests that loss of HP1B leads to a  Table. For the results from the basic GO term analysis see S3 (upregulated genes) and S4 (downregulated genes) Tables.
https://doi.org/10.1371/journal.pone.0205867.g008 profoundly altered metabolism which might be the root of the various other phenotypes of HP1b mutants noted above.

Levels of some Krebs cycle metabolites and mitochondrial complexes are lowered in HP1b mutants
Given that the expression analysis suggested that metabolic processes were strongly impacted by loss of HP1B, we carried out a targeted metabolite assay for Krebs cycle intermediates using mass spectrometry. Metabolites assayed included pyruvate, lactate, citrate, cis-aconitate, isocitrate, 2-hydroxyglutarate, α-keto-gluterate, succinate, fumarate, malate, oxalo-acetate, glutamate, glutamine, and aspartate. Of these compounds, cis-aconitate, isocitrate, α-keto-gluterate, and oxalo-acetate were below the limit of quantification. Among the remaining metabolites, citrate and malate levels are impacted by loss of HP1B (Fig 9). Citrate levels in extracts from HP1b mutant animals are lower than the citrate levels in the yw control samples (p<0.05, t-test, Fig 9A). Malate levels are also lower in the extracts from animals lacking HP1B, but the difference is not significant (p = 0.067159718, t-test, Fig 9B), most likely due to the relatively small sample size. Together, these results suggest that the levels of some Krebs cycle metabolites are indeed impacted by the loss of HP1B.
Because both our metabolite and RNA-seq analysis indicated that the mitochondria might be impacted by loss of HP1B, we next assayed two mitochondrial enzymes for their activity, citrate synthase (chosen based on the results from the metabolite analysis) and mitochondrial complex III, which has been linked to an inability to exercise in humans [60]. We found that, while citrate synthase activity levels were similar for HP1B mutant and yw control animals (ANOVA, p>0.05), HP1b mutants had decreased levels of complex III activity (Fig 9C,  ANOVA, HP1b 16 and yw p = 0.048, HP1b 86 and yw: p = 0.006). On average, the activity measured for complex III in yw was 5.4nmol/mg/min (+/-0.9), while the average for HP1b mutants was 3.3nmol/mg/min (+/-0.8, HP1b 16 ) and 2.5nmol/mg/min (+/-1.7; HP1b 86 ). Thus, complex III activity was reduced by approximately 48% in the mutant animals compared to the controls. These results demonstrate that loss of HP1B significantly impacts mitochondrial activity, which might explain the lower activity levels of the animals.

Discussion
In our study, we describe for the first time the phenotypes associated with null mutations in HP1b, a member of the HP1 family in Drosophila. Our results indicate that HP1B is not essential for survival, but nonetheless impacts important organism-level phenotypes. While we find that females lacking HP1B produce fewer adult offspring than control animals, homozygous HP1b mutants can be maintained as a healthy stock. These findings are consistent with the results of an RNAi knockdown study carried out in 2011, which found slightly decreased viability in flies with lower levels of HP1B [61]. With respect to viability, HP1B thus is similar to HP1C, another somatic HP1 homolog in Drosophila. While HP1c null alleles typically are maintained in a heterozygous stock, homozygous mutant flies survive to adulthood and are fertile [62]. The non-essential nature of HP1B and HP1C contrasts with the lethality seen for HP1a, which mainly functions in heterochromatin [24,25,48]. Thus, these observations indicate that the phenotypes of null mutations in the Drosophila HP1 proteins reflect their functional divergence.
Similar to mutations in many other chromatin proteins, mutations in HP1B affect PEV. However, as with the survival and fertility phenotypes, HP1B's impact on PEV is distinct from that of HP1a and HP1C. HP1a/Su(var)205 mutations are strong dominant suppressors of PEV for transgenes in pericentric heterochromatin and on chromosome 4 in both males and females [63][64][65][66][67]. Tethering HP1C to a transgene leads to increased transcription from the reporter [53], and analysis of an HP1C deficiency line suggests that it acts as dominant suppressor of variegation at the telomere (telomere position effect), but has no impact on PEV assayed by w m4 [68]. The results presented here for HP1B show that it acts as an enhancer of variegation in males in all sequences contexts examined (pericentric, telomeric, and chromosome 4). In females, HP1b mutations suppress variegation on chromosome 4, but enhance variegation for pericentric reporters. This discrepancy between the sexes is unusual, but most likely due to the fact that the HP1b gene is located on the X chromosome, and males thus lack HP1B completely, while females retain one dose, suggesting that impact of HP1B on PEV is dose-dependent. Additional experiments will be required to determine if females respond to a complete lack of HP1B the same way as males and how HP1B is able to impact sites at the telomere and centromere, as these domains are typically controlled by distinct sets of modifiers [66].
While not significantly impacting stock viability, phenotypic characterization of HP1b mutants revealed that loss of this chromatin protein has far-reaching consequences. The animals exhibit changed behaviors such as lowered food intake and lowered activity levels throughout development. These behavioral alterations lead to a shift in body composition, with the animals accumulating high body fat levels especially early in life. Together, these changes result in two stress resistance phenotypes, starvation stress resistance and oxidative stress resistance. Contrary to common belief, these stress resistance phenotypes are not due to an up-regulated stress response pathway; they are due to altered behavior and body composition: lower activity levels and higher initial body fat content allows the HP1b mutant animals to survive longer under starvation conditions, while the lower food intake allows them to better tolerate the oxidative stress conditions as they ingest less oxidizer.
Currently, it is unclear how often phenotypes like the ones described here for HP1b mutants occur in other chromatin mutants as few studies go beyond characterizing viability and fertility/offspring number. Fertility defects are seen in a number of chromatin mutants, including egg (the Drosophila SETDB1 homolog; [69,70]), su(Hw) (an insulator protein; [71,72]) and rhino (a female germline specific HP1 homolog; [73]). In mutants lacking the H3K9 histone methyltransferase SU(VAR)3-9 or the H3S10 kinase JIL-1, fewer than the expected offspring number are seen, suggesting that fertility and viability are affected [74,75]. Given the frequency in which fertility and survival are affected by mutations in chromatin components, it is likely that to date unreported phenotypes similar to those found in HP1b mutants might occur as well.
The RNA-seq analysis suggests that the various phenotypes of HP1b mutants originate in their altered metabolism. This hypothesis is supported by the large fraction of genes associated with metabolism GO terms that are found to be misregulated in the RNA-seq study, the results of the targeted metabolite study, and the mitochondria function assay. To date, no other metabolism studies have been reported for mutants lacking HP1 proteins. However, an RNAi knockdown study of HP1 proteins in Drosophila tissue culture cells found a significant enrichment of genes associated with metabolism (based on GO terms) among the misregulated gene set [76]. These results indicate that the link between HP1B and metabolism is conserved in the tissue culture system as well as animals.
Results from other model organisms further support a link between HP1 proteins and metabolism. Mice, as humans, have three HP1 homologs, HP1α, HP1β, and HP1γ (CBX5, CBX1, CBX3). A hypomorphic allele of the HP1 homolog Cbx3 (HP1γ) caused death in most mice prior to weaning and is associated with small pup size and altered energy homeostasis [77]. In C. elegans, loss of the HP1 homolog HP1L-2 leads to altered lipid metabolism in the animals [78]. These examples illustrate that the link identified in this study between an HP1 protein and metabolism might exist in other species as well.
The mechanism by which HP1B exerts its influence on metabolism genes is currently unclear. The finding that HP1B acts as E(var) in this study would suggest that it has a transcription activating function. However, all reporter genes tested in this study were located in heterochromatin, and when a lacI tagged version of HP1B was targeted to two euchromatic transgenes, silencing was reported [53]. Gene expression studies presented here as well as previous expression microarray [76] and qPCR studies [61] also show varying results in that loss of HP1B can lead to increased as well as decreased gene expression. Together, these findings suggest that the impact of HP1B loss on gene expression might be context-dependent.
Context-dependent impacts on gene expression have also been reported for another HP1 protein, HP1a in Drosophila. HP1a is mostly a transcriptional repressor [79], but it also acts as transcriptional activator for approximately 100 euchromatic loci [80] as well as several heterochromatic loci [49,[81][82][83]. Similar complex associations between HP1 proteins and transcription have been described in other species including mammals [reviewed in [84]], indicating that context-dependent roles in gene regulation are a common characteristic of HP1 proteins.
The complex impacts of HP1 proteins on gene expression might be explained in part by their diverse set of interacting proteins. Interactors include activators, repressors, and RNA pol II itself among many others, including other HP1 proteins [for data from Drosophila, see [85,86]; reviewed in [84]]. HP1B and HP1C specifically bind to many of the same genomic target sites [87][88][89], which raises the possibility that these two proteins could partially compensate for each other's function. However, in the RNA-seq analysis presented here, neither HP1a/Su (var)205 nor HP1c gene expressing are altered significantly. Currently, it is mostly unclear what happens to protein partners if individual HP1 proteins are depleted, and how this depletion affects the various protein complexes HP1 proteins participate in. Thus, knowledge of what additional protein interactions occur at individual HP1 protein binding sites might be necessary to predict the effects of HP1 proteins-and that of their removal-on gene expression at the target site. Further investigation of these interactors might also reveal how HP1B-and other HP proteins-are recruited to their target sites and how their regulation differs among HP1 paralogs. These studies will be essential to determine how the chromatin protein HP1B is linked to the control of metabolism as demonstrated in our studies.
Supporting information S1  Fig. HP1b mutations enhance Sb v variegation, but can act as Su(var) or E(var), depending on the reporter insertion site of a white reporter. A. Eye pigment assays for females from the five reporter insertions assayed. The response of an hsp70-white reporter to a single mutant HP1b allele in female flies depends on the reporter location. � : comparisons of mutant eye pigment levels to the ywcontrol that are statistically significant (p<0.05, t-test). Y-axis: OD 480 measuring eye pigment. X-axis: reporter insertions. Results from the ywcontrol are in red, from HP1b 16 in blue, and from HP1b 86 in green. Box plots: black bar-median; box-+25% and -25% quartiles; whiskers-maximum and minimum; circles-outliers; n = 6-10. B. Stubble variegation assays indicate that HP1b is an enhancer of variegation, as the presence of the HP1b mutant alleles significantly increases the number of wildtype bristles compared to the ywcontrol (p<0.01, Tukey multiple comparisons of means; n = 19-39). Data from females. Results from the ywcontrol are in red, from HP1b 16 in blue, and from HP1b 86 in green. Box plots: black bar-median; box-+25% and -25% quartiles; whiskers-maximum and minimum; circles-outliers. show significantly lower activity levels than animals of the yw control genotype (p<0.001 for both comparisons; Tukey's HSD). B. In males, the activity levels of animals lacking HP1B is reduced compared to yw, but the effect is less pronounced than in females and only significant for the HP1b 86 allele (p-value not significant for HP1b 16 , and p<0.001 for HP1b 86 , Tukey's HSD). Y-axis: activity level as measured by recorded beam crossings in an activity monitor for 10 flies within the 2hr assay period. Black diamond and bar: mean +/-SD. N = 25 per sex/genotype combination. (TIF)

S3 Fig. Loss of HP1B impacts survival during stress in both males and females. A. Female
HP1b mutant animals survive significantly longer during starvation condition than yw control animals (top; p = 6.455e-09 and p = 4.672e-07 for HP1b 16 and HP1b 86 respectively; Kruskal-Wallis rank sum test), while for males, only the comparison between HP1b 16 and yw is significant (bottom; p = 2.054e-07; Kruskal-Wallis rank sum test). B. Female HP1b mutant animals survive significantly longer after exposure to the oxidizer paraquat than yw control animals (top; p = 3.752e-10 and p < 2.2e-16 for HP1b 16 and HP1b 86 respectively; Kruskal-Wallis rank sum test), as do males (bottom; p = 1.252e-05 and p = 1.79e-14 for HP1b 16 and HP1b 86 respectively; Kruskal-Wallis rank sum test). C. HP1b mutant animals, both male and female, do not show improved survival during heat stress conditions (37˚C) compared to yw control animals. In contrast, HP1b 86 shows significantly lower survival (p = 0.001803 in females [top] and p = 1.499e-08in males [bottom]; Kruskal-Wallis rank sum test). For A-C: X-axis-time to death in hours. Y-axis-survival probability. Data shown are from three trials, each with 100 animals per genotype. (TIF) S4 Fig. Loss of HP1B impacts the survival curve of females and males differently. A. Survival curves differ significantly between the HP1b mutant strains and the yw control strain in females (HP1b 16 : p = 0.0259; HP1b 86 : p = 0.00563; log rank test). B. In males, the survival curves of HP1b mutant animals do not differ significantly from the survival curves of yw control animals (p not significant; log rank test). A+B: X-axis-time to death in days. Y-axissurvival probability. Data shown are combined from three trials, each with 100 animals per genotype/sex. (TIF) S5 Fig. Altered feeding behavior is the likely cause of increased oxidative stress resistance as measured by paraquat feeding in both males and females. The paraquat concentration was adjusted to take into account differences in feeding behavior between HP1b mutants (blue, green) and the yw control strain (red). A. Females. After adjusting the paraquat concentration, female HP1b mutant animals die earlier than their yw counterparts (HP1b 16 : p = 7.413e-11; HP1b 86 : p = 2.113e-14; Kruskal-Wallis rank sum test). B. Males. After adjusting the paraquat concentration, no difference in survival is detected between HP1b mutants and yw (HP1b 16 : p = 0.31; HP1b 86 : p = 0.9883; Kruskal-Wallis-rank sum test). Y-axis: survival probability; X-axis: Time to death in hours. N = 100 animals per genotype/sex.