Author Response: H3K9me1/2 Methylation Limits the Lifespan of Daf-2 Mutants in C. Elegans

Meng Huang,Minjie Hong,Xinhao Hou,Chengming Zhu,Di Chen,Xiangyang Chen,Shouhong Guang,Xuezhu Feng

DOI: https://doi.org/10.7554/elife.74812.sa2

2022-01-01

Abstract:Article Figures and data Abstract Editor's evaluation Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract Histone methylation plays crucial roles in the development, gene regulation, and maintenance of stem cell pluripotency in mammals. Recent work shows that histone methylation is associated with aging, yet the underlying mechanism remains unclear. In this work, we identified a class of putative histone 3 lysine 9 mono/dimethyltransferase genes (met-2, set-6, set-19, set-20, set-21, set-32, and set-33), mutations in which induce synergistic lifespan extension in the long-lived DAF-2 (insulin growth factor 1 [IGF-1] receptor) mutant in Caenorhabditis elegans. These putative histone methyltransferase plus daf-2 double mutants not only exhibited an average lifespan nearly three times that of wild-type animals and a maximal lifespan of approximately 100 days, but also significantly increased resistance to oxidative and heat stress. Synergistic lifespan extension depends on the transcription factor DAF-16 (FOXO). mRNA-seq experiments revealed that the mRNA levels of DAF-16 Class I genes, which are activated by DAF-16, were further elevated in the daf-2;set double mutants. Among these genes, tts-1, F35E8.7, ins-35, nhr-62, sod-3, asm-2, and Y39G8B.7 are required for the lifespan extension of the daf-2;set-21 double mutant. In addition, treating daf-2 animals with the H3K9me1/2 methyltransferase G9a inhibitor also extends lifespan and increases stress resistance. Therefore, investigation of DAF-2 and H3K9me1/2 deficiency-mediated synergistic longevity will contribute to a better understanding of the molecular mechanisms of aging and therapeutic applications. Editor's evaluation This work presents important findings which pertain to how the length of an organism's lifespan can be affected by the presence of post-translational histone modifications, namely histone H3 lysine 9 mono- or di-methylation (H3K9me1/2). The authors find that the presence or absence of these marks plays a key role in modulating lifespan and stress resistance that synergizes with a long-lived insulin-pathway (daf-2) model. This work adds an important layer to aging regulation in organisms where a repressive chromatin state in specific genomic regions may limit the extent of lifespan. Hence, proper elucidation of this mechanism would be valuable to attain a comprehensive understanding of aging control and the reasons why these limits exist. https://doi.org/10.7554/eLife.74812.sa0 Decision letter Reviews on Sciety eLife's review process Introduction Lifespan is governed by complex interactions between genetic and environmental factors. The perturbation of insulin/insulin-like signaling, target of rapamycin (TOR) pathway, and mitochondrial functions have been shown to extensively modulate lifespan and health (Fontana et al., 2010; Kenyon, 2010; López-Otín et al., 2013). These genetic manipulations often lead to significant changes in gene expression at both the transcriptional and translational levels. Inhibition of DAF-2, the Caenorhabditis elegans ortholog of the insulin growth factor 1 (IGF-1) receptor, doubles adult lifespan by activating the DAF-16 (FOXO) transcription factor to regulate downstream genes involved in stress resistance, detoxification, and metabolism (Kenyon et al., 1993; Kimura et al., 1997; Lin et al., 1997; Ogg et al., 1997; Murphy et al., 2003; McElwee et al., 2004). In addition to genetic regulation, aging is also modulated by epigenetic processes. Epigenetic marks, including histone acetylation and methylation, as well as the associated chromatin states, are altered during aging. Aging-dependent loss of chromatin repression is correlated with alterations in gene expression patterns, which have been documented in species ranging from C. elegans to humans (Lee et al., 2000; Lund et al., 2002; Bennett-Baker et al., 2003; Lu et al., 2004). For example, the methylation of histone 3 lysine 4 methylation (H3K4me) is one of the posttranslational histone modifications that marks on regions of active transcription (Shilatifard, 2008). H3K4me is deposited by the MLL/COMPASS complex, which travels with elongating RNA Polymerase II during transcription (Wood et al., 2007). In C. elegans, animals with reductions in COMPASS complex subunits (wdr-5, ash-2, and set-2) live longer than wild-type individuals (Greer et al., 2010). However, little is known about the importance of repressive histone modification for aging marks. The H3K27me3 demethylase UTX-1 regulates lifespan independently of the presence of the germline but in a manner that depends on the insulin-FOXO signaling pathway (Maures et al., 2011). MET-2, a mammalian H3K9 methyltransferase SETDB1 homolog (Loyola et al., 2006), likely monomethylates and dimethylates H3K9 in C. elegans (Towbin et al., 2012). MET-2 is necessary both for a normal lifespan (Tian et al., 2016) and for the lifespan extension of wdr-5 mutants (Lee et al., 2019). SET-6 is a putative H3K9me2/3 methyltransferase but not H3K9me1. Although the depletion of SET-6 does not significantly increase the lifespan of C. elegans, it may increase healthy aging (Yuan et al., 2020). SET-25 is a putative trimethylase for H3K9 (Towbin et al., 2012; Padeken et al., 2021). However, the deletion of SET-25 does not significantly change the worm lifespan (Woodhouse et al., 2018). The enzymes responsible for histone lysine methylation are called histone methyltransferases (HMTs). HMTs typically contain a conserved catalytic domain called SET, which stems from Su(var)3–9, Enhancer of zeste, and Trithorax, the first HMTs known to carry this domain (Tschiersch et al., 1994). The C. elegans genome encodes 38 SET domain-containing proteins, of which 5 are essential for viability (Andersen and Horvitz, 2007; Ni et al., 2012). However, the biological roles of most SET proteins are largely unknown. To investigate the function of repressive histone methylation in lifespan regulation, we selected a number of putative H3K9 methyltransferases as well as methyltransferases that may indirectly regulate H3K9me and tested whether the loss-of-function of these SET proteins could change the lifespan of C. elegans. Interestingly, we found that the H3K9me1/2, but not H3K9me3, mutants exhibited a synergistic lifespan extension with daf-2 mutation. These animals showed an average lifespan of approximately 60 days, which is approximately 70% longer than that of daf-2 worms and is three times as long as that of wild-type N2 animals. The double mutants exhibited a maximal lifespan of approximately 100 days. The synergistic lifespan extension of DAF-2 and lifespan-limiting HMT mutants depend on the transcription factor DAF-16 (FOXO). mRNA-seq experiments showed that DAF-16 Class I genes (genes that are activated by DAF-16) are activated in long-lived worms. Therefore, we conclude that H3K9me1/2 may limit the lifespan of daf-2 animals by repressing the expression of DAF-16 Class I genes. Results The depletion of set-21 extends lifespan and enhances stress resistance in daf-2 mutant SET-21 is a homolog of human isoform 1 of histone-lysine N-methyltransferase SUV39H2, with functions unknown in C. elegans (Wenzel et al., 2011). We accidentally found that set-21(ok2320) significantly extended lifespan in daf-2 mutants, but did not change the lifespan of N2 animals (Figure 1A, Supplementary file 1 Statistical analyses of lifespan experiments). To confirm the function of SET-21 in lifespan regulation, we first generated a number of additional deletion alleles of set-21 by CRISPR/Cas9 technology (Figure 1—figure supplement 1A). The deletion of set-21 did not significantly extend the lifespan or change the brood size in N2 background (Figure 1A, Figure 1—figure supplement 1B). Strikingly, in daf-2(e1370) mutant worms, knocking out set-21 reduced the brood size and significantly extended lifespan (Figure 1A, Figure 1—figure supplement 1B). The average lifespan of daf-2(e1370);set-21(ust68) were 55% longer than that of daf-2(e1370) animals (Figure 1A). And the maximal lifespan of daf-2;set-21 animals achieved approximately 100 days, which is three times as long as that of N2 animals. Figure 1 with 1 supplement see all Download asset Open asset Synergistic lifespan extension and stress resistance in daf-2;set-21 mutants. (A) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of two independent experiments. The percentage of change was compared to the average lifespan of control animals. Asterisks indicate significant differences using log rank tests. ** 0.001< p < 0.01; ****p<0.0001. n.s., not significant, p>0.05. Lifespan data were summarized in Supplementary file 1. (B) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of two independent experiments. Asterisks indicate significant differences using log rank tests. ****p<0.0001. Lifespan data were summarized in Supplementary file 1. (C, D) Survival curves of the indicated animals. (C) Oxidative and (D) heat stress. Data in each panel are presented as the mean ± s.e.m. of five independent experiments. Asterisks indicate significant differences using two-tailed t tests. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; n.s., not significant, p>0.05. The same oxidative stress data of control N2 and daf-2 animals were used in (C) Figures 2 and 3 and Figure 3—figure supplement 3A. The same heat stress data of control N2 and daf-2 animals were used in (D) Figures 2 and 3 and Figure 3—figure supplement 3B. eat-2 mutant is a genetic model in dietary restriction (DR) research in C. elegans. The mutation of eat-2 renders a non-efficient pharynx in grinding bacteria and results in DR on regular media plates (Raizen et al., 1995). eat-2 animals exhibit phenotypes similar to those observed in other species subjected to DR, including an ~36% longer lifespan (Figure 1B; Lakowski and Hekimi, 1998). Similarly, knocking out set-21 further extended the lifespan of eat-2(ad465) mutant worms (Figure 1B). The average lifespan of eat-2(ad465);set-21(ust68) were 16% longer than that of eat-2(ad465) animals, suggesting a broader effect of these lifespan-limiting SET proteins. The extended lifespan of nematodes has been shown to correlate with the activation of stress response genes and increased stress resistance (Lithgow et al., 1995; Murphy et al., 2003; Samuelson et al., 2007). The deletion of set-21 did not significantly change the stress resistance in N2 background animals (Figure 1C–D). The daf-2 mutation led to a higher resistance to oxidative stress and heat-shock stress than N2 animals. Strikingly, daf-2;set-21 worms revealed a much higher resistance to oxidative stress via hydrogen peroxide treatment and heat-shock stress than daf-2 animals (Figure 1C–D). The depletion of met-2 extends lifespan in daf-2 mutant worms and enhances stress resistance The remarkable lifespan extension in daf-2;set-21 animals inspired us to re-investigate the role of histone H3 lysine 9 methylation in C. elegans. MET-2 is the mammalian SETDB1 homolog, which is involved in mono- and dimethylation of H3K9 (Figure 2A; Loyola et al., 2006; Towbin et al., 2012). The brood size of met-2 mutant is approximately half of that of wild-type N2 animals (Figure 2B). Interestingly, daf-2;met-2 double mutant exhibited synthetic defect and had a significantly reduced brood size. met-2 mutants exhibited a modest shorter lifespan than wild-type N2 animals (Figure 2C; Tian et al., 2016). Strikingly, daf-2;met-2 double mutants revealed an average lifespan of approximately 47 days, which is 30% longer than that of daf-2 mutation alone and is 2.3 times as long as that of wild-type N2 animals (Figure 2C). The depletion of met-2 enhanced the oxidative stress resistance and heat stress resistance in both N2 and daf-2 mutant worms (Figure 2D–E). These results suggest that met-2 and H3K9me1/2 may play different roles in lifespan regulation in N2 and daf-2 animals. Alternatively, the depletion of daf-2 may provide a sensitized genetic background to identify genes that are involved in aging and lifespan regulation. Figure 2 Download asset Open asset daf-2;met-2 mutants revealed extended lifespan and increased resistance to oxidative and heat stress compared to daf-2 animals. (A) Western blotting of L4 stage animals with the indicated antibodies. Numbers indicate the scanned density by ImageJ. (B) Brood size of indicated animals at 20°C. (C) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of three independent experiments. Asterisks indicate significant differences using log rank tests. ****p<0.0001. Lifespan data were summarized in Supplementary file 1. (D, E) Survival curves of the indicated animals. (D) Oxidative and (E) heat stress. Data in each panel are presented as the mean ± s.e.m. of five independent experiments. Asterisks indicate significant differences using two-tailed t tests. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; n.s., not significant, p>0.05. Figure 2—source data 1 Original files and figures of the blots in Figure 2A. https://cdn.elifesciences.org/articles/74812/elife-74812-fig2-data1-v2.zip Download elife-74812-fig2-data1-v2.zip SET-25 is a putative trimethylase for H3K9 (Figure 2A; Towbin et al., 2012). However, deletion of SET-25 did not significantly change the worm lifespan or stress resistance in either the wild-type N2 or daf-2 background animals (Figure 2C–E). The deletion of SET-25 only moderately enhanced the oxidative stress resistance in daf-2 mutant worms (Figure 2D–E). Identification of set genes required for lifespan limitation in daf-2 animals C. elegans’ genome encodes 38 proteins containing SET domains, most of which have unknown function (Andersen and Horvitz, 2007). To investigate the function of H3K9 methylation in lifespan regulation, we searched for other SET proteins that exhibit sequence similarity to SET-21 (Figure 3A, Figure 3—figure supplement 1). We selected seven additional SET proteins, including SET-6, SET-13, SET-15, SET-19, SET-20, SET-32, and SET-33, which are most closely related to the SET-21 in the dendrogram. All these proteins contain conserved catalytic SET domains and are annotated as putative H3K9 methyltransferases, as well as methyltransferases that may indirectly regulate H3K9me. Among them, SET-6 is required for H3K9me2/3 methylation (Yuan et al., 2020) and SET-32 is required for nuclear RNAi-directed H3K9me3 in C. elegans (Mao et al., 2015; Kalinava et al., 2017; Spracklin et al., 2017). Using N2 and the sensitized daf-2 background, we assayed the lifespan of these putative H3K9 methyltransferase mutants. Figure 3 with 3 supplements see all Download asset Open asset Synergistic lifespan extension of set-6, set-19, set-20, set-32, and set-33 with daf-2 animals. (A) Dendrogram comparing the protein sequences of SET proteins that are putative H3K9 methyltransferases or methyltransferases that may indirectly regulate H3K9me in Caenorhabditis elegans. The numbers indicate the length of each protein. (B) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of two independent experiments. The percentage of change was compared to the average lifespan of N2 animals. Asterisks indicate significant differences using log rank tests. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; n.s., not significant, p>0.05. Lifespan data were summarized in Supplementary file 1. (C) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of two independent experiments. The percentage of change was compared to the average lifespan of daf-2 animals. Asterisks indicate significant differences using log rank tests. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; n.s., not significant, p>0.05. (D, E) Survival curves of the indicated animals. (D) Oxidative and (E) heat stress. Data are presented as the mean ± s.e.m. of five independent experiments. Asterisks indicate significant differences using two-tailed t tests. *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; n.s., not significant, p>0.05. We acquired additional set mutants from CGC and also generated alleles by CRISPR/Cas9 technology (Figure 3—figure supplement 2A). Among them, the deletion of set-6, set-19, set-20, and set-32 modestly increased the lifespan compared to that of wild-type N2 animals (Figure 3B). However, in the daf-2 mutant background, mutations of set-6, set-19, set-20, set-32, and set-33 exhibited a striking synergistic lifespan extension. While the lifespan of daf-2;set-20 and daf-2;set-32 are approximately 60% longer than that of daf-2 worms, the daf-2;set-6 and daf-2;set-19 live 70% longer than that of daf-2 animals (Figure 3C). Especially, daf-2;set-19 exhibits a maximal lifespan of approximately 100 days, which is five times as long as the average lifespan of N2 animals. daf-2;set-6, daf-2;set-19, daf-2;set-20, daf-2;set-32, daf-2;set-33 had fewer brood sizes than the single set gene mutants (Figure 3—figure supplement 2B). The long-lived animals were also much more resistant than control animals to the oxidative stress induced by hydrogen peroxide (Figure 3D, Figure 3—figure supplement 3A) and heat stress (Figure 3E, Figure 3—figure supplement 3B) in daf-2 mutant worms but not in N2 background animals. set-32(ok1457) revealed a higher resistance to heat stress but not oxidative stress than N2 animals. To test whether SET-6, SET-19, SET-20, SET-21, SET-32, and SET-33 act in the same genetic pathway to regulate lifespan in daf-2 animals, we crossed set-21(ust68) to other set mutants. As expected, the triple mutants daf-2;set-21;set-6, daf-2;set-21;set-19, daf-2;set-21;set-20, daf-2;set-21;set-32, and daf-2;set-21;set-33 did not significantly further extend the lifespan than the double mutants (Figure 4A). Thus, these set genes probably act in the same genetic pathway to regulate lifespan. Figure 4 with 1 supplement see all Download asset Open asset Genetic pathway analysis of set genes in lifespan regulation in daf-2 mutants. (A–D) (Left) Survival curves of the indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Means ± s.e.m. of two independent experiments. The percentage change was compared to the lifespan of control animals. Asterisks indicate significant differences using log rank tests. ** 0.001< p < 0.01; ****p<0.0001; n.s., not significant, p>0.05. Lifespan data were summarized in Supplementary file 1. We speculated that there might be a correlation between the expression pattern and lifespan extension of the set genes Thus, we generated in situ single-copy 3xFLAG-GFP tagged SET-21, SET-25, and SET-32 transgenes by CRISPR/Cas9 technology and single-copy GFP-3xFlAG-tagged SET-13 transgenic strains using the Mos1-mediated single-copy insertion (MosSCI) technology. These animals were imaged in N2 and daf-2 background at embryonic and larva stages, respectively, to figure out the expression patterns and cellular localizations of SET proteins (Figure 4—figure supplement 1). However, we failed to find a clear correlation between their expression pattern with the lifespan extension phenotype. The depletion of SET-21 and SET-32 extended the lifespan of daf-2 animals, but GFP::SET-21 and GFP::SET-32 showed different expression patterns and subcellular localizations. GFP::SET-21 exclusively located in the nucleus but GFP::SET-32 located in both the nucleus and cytoplasm in embryos. At larval stage, GFP::SET-21 was expressed in soma but GFP::SET-32 was expressed in germline. Both set-13 and set-25 mutants did not show lifespan extension in N2 and daf-2 background worms. SET-13 and SET-25 were enriched in the nucleus. SET-13::GFP was expressed through eight-cell stage to late-embryonic stage, while SET-25::GFP was expressed in the whole embryonic stage and the germline in young adults. H3K9me3 is a hallmark of heterochromatin and is often enriched on repetitive sequences (Peters et al., 2003; Guenatri et al., 2004; Martens et al., 2005). The heterochromatin protein 1 (HP1) is a highly conserved small nonhistone protein, which was first identified in Drosophila, and recognizes H3K9me3 (James and Elgin, 1986). C. elegans has two HP1 paralogs: HP1 Like 1 and 2 (HPL-1 and HPL-2) (Couteau et al., 2002). The chromatin binding pattern of HPL-2 most closely resembles H3K9me1 and H3K9me2 and less closely resembles H3K9me3 (Garrigues et al., 2015). hpl-2(tm1489) mutation has been shown to extend lifespan (Meister et al., 2011). However, both daf-2;hpl-1 and daf-2;hpl-2 double mutants did not exhibit further lifespan extension compared to the daf-2 mutant (Figure 4B). SET-25 colocalized with both H3K9me3 and HPL-1, but did not colocalize with HPL-2 (Towbin et al., 2012). SET-25 is required for H3K9me3 methylation but is not involved in lifespan limitation (Figure 2A and C), further supporting that H3K9me3 may be dispensable for lifespan regulation in daf-2 animals. Previous work has shown that knocking out the H3K4me3 methyltransferase SET-2 extended worm lifespan of N2 animals (Greer et al., 2010). To test whether SET-2-dependent H3K4me3 and lifespan-limiting H3K9 methyltransferases act in the same genetic pathway to regulate lifespan, we crossed set-2 to daf-2;set-21 animals. The double mutants daf-2;set-2 and daf-2;set-21 and the triple mutants daf-2;set-2;set-21 exhibited similar lifespan extensions relative to that of daf-2 animals (Figure 4C), suggesting that H3K4me3 and H3K9me1/2 may function in the same genetic pathway or regulate the same cohorts of target genes for lifespan modulation. Inhibition of RSKS-1 (S6K), the TOR pathways, extends lifespan in C. elegans (Hansen et al., 2007; Pan et al., 2007), and the double mutant of daf-2;rsks-1 leads to synergistically prolonged longevity (Chen et al., 2013). However, the triple mutants daf-2(e1370);rsks-1(ok1255);set-21(ust68) revealed similar lifespan as long as those of daf-2(e1370);rsks-1(ok1255) and daf-2(e1370);set-21(ust68) double mutants (Figure 4D), suggesting that either rsks-1 and set-21 act in the same molecular pathway to regulate lifespan or there is likely an upceiling of maximal lifespan of approximately 100 days for C. elegans. Interesting, a report suggested that in a mutant of age-1, which encodes the class-I phosphatidylinositol 3-kinase catalytic subunit (PI3K(CS)), worms can survive to a median of 145–190 days at 20 degrees and a maximal lifespan of approximately 260 days, with nearly 10-fold extension of both median and maximum adult lifespan relative to control animals (Ayyadevara et al., 2008). Therefore, we concluded that a number of set genes limit the lifespan and stress resistance of daf-2 animals in C. elegans. Lifespan-limiting SET proteins regulate DAF-16 Class I genes in daf-2 animals The longevity phenotype of daf-2 animals depends on the downstream DAF-16 transcription factor (Kenyon et al., 1993). DAF-16 localizes in the cytoplasm in N2 worms but accumulates in the nucleus upon daf-2 mutation (Henderson and Johnson, 2001; Lin et al., 2001). To determine whether the expression pattern and subcellular localization of DAF-16 are altered upon the mutation of these lifespan-limiting HMTs, we crossed a widely used high-copied DAF-16::GFP transgene with daf-2, set-21, set-25, daf-2;set-21, and daf-2;set-25 animals. As reported, the mutation of daf-2 resulted in a pronounced nuclear accumulation of DAF-16::GFP (Figure 5—figure supplement 1A). In N2 background, the mutation of set-21 or set-25 did not detectably change the cytoplasmic localization of DAF-16::GFP (Figure 5—figure supplement 1A). In daf-2 background worms, the mutation of set-21 or set-25 did not further induce detectable changes in the expression pattern and subcellular localization of DAF-16::GFP. We generated a single-copied DAF-16::GFP transgene by inserting a GFP tag onto the C-terminal of DAF-16 via CRISPR/Cas9 technology (Figure 5A). Then, we crossed the single-copied DAF-16::GFP transgene with daf-2, set-21, set-25, daf-2;set-21, and daf-2;set-25 animals. DAF-16::GFP accumulated in the cytoplasm in N2 background animals. Unlike in the high-copied DAF-16::GFP transgenic animal, daf-2 mutation only modestly enriched DAF-16::GFP to the nucleus in the single-copied DAF-16::GFP transgenic strain (Figure 5B and Figure 5—figure supplement 1B). Strikingly, in daf-2;set-21 and daf-2;set-25 double mutant animals, DAF-16::GFP was strongly enriched in the nucleus in the single-copied transgenic animals. The reason for the subcellular localization difference between the high-copied and single-copied DAF-16::GFP transgenic animals is unclear. However, the single-copied DAF-16::GFP may be used to re-visit the mechanism of lifespan regulation pathways. In addition, the mechanism of enhanced nuclear accumulation of DAF-16::GFP in daf-2;set-21 and daf-2;set-25 animals is unknown. Given that daf-2;set-21 exhibited extended lifespan but daf-2;set-25 did not, and both of them revealed strong nuclear localization of DAF-16::GFP, it is quite possible that there is an additional level of regulation downstream of the nuclear localization of DAF-16. Figure 5 with 2 supplements see all Download asset Open asset The lifespan-limiting SET proteins regulate DAF-16 Class I genes in daf-2 mutants. (A) Construction of an in situ GFP-tagged DAF-16 transgene by CRISPR/Cas9 technology. The subcellular localization of DAF-16::GFP was scored at the pharynx in young adult animals. (B) (Left) Fluorescent images of the single-copied DAF-16::GFP transgene in indicated animals at 20°C. (Right) The subcellular localization of the single-copied DAF-16::GFP transgene was scored in indicated animals. (C, D) (Left) Survival curves of indicated animals. (Right) Histogram displaying the average lifespan of the indicated animals. Mean ± s.e.m. of three independent experiments. Asterisks indicate significant differences using log rank tests. ****p<0.0001. Lifespan data were summarized in Supplementary file 1. (E) Heatmap of the standardized fragments per kilobase of transcript per million mapped reads (FPKM) of reported top 50 DAF-16 Class I and Class II genes by mRNA-seq in the indicated animals. Statistical analysis was performed to obtain the expression levels of each gene in each strain and the average of the expression levels in the five strains. The gene expression levels in a single strain were compared with the average of the five strains, and the ratio obtained was processed by log2. The resulting value>0 means the expression level of the gene in the indicated strain is higher than the average expression levels of this gene in the five strains, as shown from yellow to red. The resulting value<0 means the expression level of the gene in the indicated strain is lower than the average expression levels of this gene in the five strains, as shown in blue. The expression levels are indicated by the color bar. The mRNA-seq was conducted in single replicate and validated by quantitative real-time PCR in Figure 7A. The 50 most significant positive (Class I) and negative (Class II) targets of DAF-16 were selected according to their responsiveness to DAF-16 (Tepper et al., 2013). To further test the role of DAF-16 in the daf-2;set-21-induced synergistic longevity, we constructed a daf-2(e1370);daf-16(mu86);set-21(ust68) triple mutant. The daf-16 mutation reverted the prolonged longevity phenotype of daf-2;set-21 to an average lifespan of 23 days, which is similar to that of set-21 or N2 alone (Figure 5C). In N2 background, daf-16 and daf-16;set-21 mutant worms lived shorter than wild-type N2 and set-21 animals (Figure 5D). The depletion of DAF-2 reduces the insulin signaling pathway and leads to both upregulation and downregulation of large sets of genes, referred to as DAF-16 Class I and II genes, respectively (Murphy et al., 2003). Class I genes are induced in daf-2 mutants but are repressed in daf-2;daf-16 double mutants. Class II genes are not induced in daf-2 mutants but are induced in daf-2; daf-16 double mutants; 1663 genes are classified as positive (Class I) DAF-16 targets and 1733 genes are classified as negative (Class II) DAF-16 targets of DAF-16 (Tepper et al., 2013). Class I genes are enriched for the Gene Ontology categories including oxidation, reduction, and energy metabolism, whereas Class II genes are enriched for genes involved in biosynthesis, growth, reproduction, and development. We performed mRNA-seq and analyzed Class I and II DAF-16 genes to identify the mis-regulated genes in daf-2;set mutants. Interestingly, the mRNA levels of DAF-16 Class I, but not Class II, genes are consistently activated in long-lived daf-2;set-19, daf-2;set-21, and daf-2;set-32 worms, than in control daf-2 and daf-2;set-25 animals (Figure 5E, Figure 5—figure supplement 2). Identification of the lifespan-limiting SET targeted genes in daf-2 animals To identify target genes in this group of lifespan extension daf-2;set mutants, we re-analyzed the mRNA expression profile in the double mutants (Figure 6, Figure 6—figure supplement 1). We identified 49 co-upregulated genes and 11 co-downregulated genes that are specifically enriched in the long-lived daf-2;set-19, daf-2,set-21, daf-2;set-32 double mutants, but not in short-lived control daf-2 and daf-2;daf-25 animals (Figure 6A–B and Figure 6—figure supplement 1A). Interestingly, among