Abstract
Circadian clocks may mediate lifespan extension by caloric or dietary restriction (DR). We find that the core clock transcription factor Clock is crucial for a robust longevity and fecundity response to DR in Drosophila. To identify clock-controlled mediators, we performed RNA-sequencing from abdominal fat bodies across the 24 h day after just 5 days under control or DR diets. In contrast to more chronic DR regimens, we did not detect significant changes in the rhythmic expression of core clock genes. Yet we discovered that DR induced de novo rhythmicity or increased expression of rhythmic clock output genes. Network analysis revealed that DR increased network connectivity in one module comprised of genes encoding proteasome subunits. Adult, fat body specific RNAi knockdown demonstrated that proteasome subunits contribute to DR-mediated lifespan extension. Thus, clock control of output links DR-mediated changes in rhythmic transcription to lifespan extension.
Introduction
Circadian (~24 h) clocks regulate a wide range of rhythmic metabolic, physiological and behavioral parameters to acclimate to environmental changes in light, temperature, and food availability (Patke et al., 2020). Circadian clock disruption has been implicated in advanced aging and the longevity response to caloric or dietary restriction (CR or DR) (Froy, 2018; Galikova and Flatt, 2010; Manoogian and Panda, 2017; Nakahata and Fukada, 2022; Zhu et al., 2022). DR, reduction in food intake without causing malnutrition, robustly extends longevity in various animal models including yeast, worms, flies, and monkeys (Green et al., 2022; Mc Auley, 2022). Yet, the molecular mechanisms by which DR delays aging are not fully understood. Understanding how the clock impacts aging and DR sensitivity may provide novel avenues to understanding aging.
The circadian clock consists of a widely conserved transcriptional feedback loop that drives 24 h molecular oscillations. In flies, the heterodimer transcription factor CLK/CYC forms the positive arm of the loop and activates their repressors, PER and TIM. The PER-TIM complex functions as the negative arm of the loop and inhibits CLK-CYC activity (Allada and Chung, 2010). This feedback loop drives core clock gene rhythms and controls rhythmic physiological, metabolic, and behavioral parameters via clock control of output genes (Patke et al., 2020). Genetically hybrid mice with a deviation of the circadian period from 24 h by over seven minutes showed a higher mortality rate than the mice with less deviated periods (Libert et al., 2012). However, whether the altered circadian period is correlated with or causes the increased mortality is not clear. Genetic inactivation of CYC ortholog Bmal1 as well as other circadian clock mutants also significantly reduced lifespan in mice (Dubrovsky et al., 2010; Fu et al., 2002; Kondratov et al., 2006; Lee et al., 2010). Yet when Bmal1 knockout was restricted to adulthood, lifespan was normal (Yang et al., 2016). While a lifelong DR did not significantly extend lifespan of Bmal1 knockout mice (Patel et al., 2016a), chronic (~2 mo) DR exposure increased core clock amplitude in mice (Patel et al., 2016b; Sato et al., 2017). This suggests that the circadian clock may be among the molecular mechanisms of DR. However, mice under DR restrict their feeding behavior to a narrow temporal window (Acosta-Rodriguez et al., 2017). Thus, DR induced changes in core clocks may instead be due to the well-known effects of time-restricted feeding (Hatori et al., 2012). Indeed, core clock genes are important for age-dependent cardiac function and lifespan extending effects of time-restricted feeding (Gill et al., 2015; Ulgherait et al., 2021). A recent study revealed that a basal level lifespan extension by DR is further increased when DR is temporally aligned with mice’s natural meal timing (i.e., during the night) {Acosta-Rodriguez, 2022 #117}. Thus, it remains unclear whether disruption of the circadian clock itself or other factors, such as a defect during development, results in lifespan reduction and is responsible for the lack of DR response.
Circadian clocks have also been implicated in aging and the DR longevity response in flies as well. DR mortality effects are rapid, fully evident within just 2 ~ 4 days of a diet shift in flies (Mair et al., 2003; McCracken et al., 2020), making them an attractive model organism for DR studies. Loss-of-function mutants in the activator and repressor complexes that “fix” the clock at different points in the cycle have tested the functional significance of the clock in aging and DR. Inhibition of neuronal Clk appears to reduce the lifespan extending DR effects, where flies were tested for DR effects with two (ad libitum and DR) diets (Hodge et al., 2022). However, this observation is inconclusive to the role of Clk for DR effects as inhibition of neuronal Clk decreases food intake {Xu, 2008 #103}. Reduction of food intake can decrease lifespan under DR while increasing lifespan under ad libitum, masking the true DR response {Flatt, 2014 #74}. per01 and tim01 mutant flies exhibit inconsistent DR longevity responses perhaps due to differences in microbial content (Katewa et al., 2016; Ulgherait et al., 2020; Ulgherait et al., 2016; Ulgherait et al., 2021). Thus, the role of core clock genes in mediating DR effects could be further clarified. Notably, the rhythmic amplitude of core clock genes of flies is enhanced after chronic (>10 days) DR (Katewa et al., 2016). Knockdown of modestly CLK- and DR-regulated genes in the eye modulate lifespan without apparent effects on DR-dependent longevity (Hodge et al., 2022). tim overexpression increased the amplitude of core clock oscillations and extended lifespan under control but not DR diets (Katewa et al., 2016). While clock oscillation amplitudes between tim overexpression on a control diet and wild-type flies under DR are comparable, their lifespan remain quite different, suggesting that core clock effects may not be required for lifespan extension (Katewa et al., 2016). Thus, it remains unclear if lifespan extension functions via the circadian clock or instead through clock output genes and what the role of specific clock output genes is in DR-dependent longevity. Using multiple diets, we demonstrate that Clk mutants suppress DR longevity and fecundity responses, providing more definitive demonstration of the role of the core clock. Nonetheless, using a shorter-term DR strategy, we reveal that primary DR effects on the circadian transcriptome spare core clock genes, suggesting a primary effect on circadian clock output. Network analysis suggests that a diet-dependent effect on a gene module containing proteasome subunit genes. Moreover, suppression of proteasome subunit expression, predominantly in the abdominal fat body, limits lifespan extension by DR. These results provide crucial genetic evidence that circadian clock output pathways, specifically those involving the proteasome, link DR-mediated changes in rhythmic transcription to lifespan extension. These studies raise the possibility of using chronotherapy to combat aging and age-related diseases.
Results
The Effects of Dietary Restriction on Lifespan and Fecundity Is Dramatically Suppressed in Mutants of the Core Clock Transcription Factor Clk
To tease apart the role of the circadian clock in the DR longevity response, we first evaluated the roles of the positive and negative arm of the feedback loop by testing per01 as in prior reports (Katewa et al., 2016; Ulgherait et al., 2020; Ulgherait et al., 2016), and ClkJrk, a dominant negative allele of Clk (Allada et al., 1998), which had not been previously examined for DR studies. We applied a common DR regimen where the concentration of both yeast and sucrose are diluted (whole food dilution: Control: 15% [w/v] Sucrose and Yeast, 15SY; DR: 5% [w/v] Sucrose and Yeast, 5SY) in 12hr light: 12hr dark (LD) cycles (Bass et al., 2007; Kabil et al., 2011). We used female flies where DR responses are more robust (Magwere et al., 2004). We observed robust lifespan extension by DR in wild-type iso31 (w1118, 29%) and comparable extension in per01 flies (25%) (Fig. S1; diet*genotype interaction p > 0.05 between iso31 and per01). These results are consistent with one report which showed that per01 flies display normal DR extension effects (Ulgherait et al., 2016). However, only ~ 10% lifespan extension by DR was observed in ClkJrk mutants (Fig. S1; diet*genotype interaction p < 0.0001 between iso31 and ClkJrk). Interestingly, the DR responses in per01 and ClkJrk flies were significantly different from each other (Fig. S1; diet*genotype interaction p < 0.0001 between per01 and ClkJrk). We hypothesize that ClkJrk and per01 arrest the clock at opposite points in the cycle, only one of which impacts the DR response. Although the DR response is typically tested by comparing lifespan with just two diets (Solovev et al., 2019), it can potentially mis-assign the effects of DR or diet (Flatt, 2014; Tatar, 2007). For example, chico mutants show little effect of DR looking at just two diet concentrations but showed an almost identical response if their lifespan was measured across seven diet concentrations (Clancy et al., 2002). To distinguish between this possibility and a “true” DR response, we performed a reaction norm analysis by comparing the mean lifespan of ClkJrk flies to that of wild-type over serially diluted diets (1, 5, 10, 15, and 20SY) (Bass et al., 2007; Flatt, 2014; Tatar, 2007). As expected in wild-type flies, we observed a reduction of lifespan as food concentration increased from 5SY to 20SY. A reduction of lifespan was also observed when going from 5SY to 1SY, presumably due to malnutrition (Fig. 1A & C). However, mean lifespan of ClkJrk mutant flies showed a more flattened reaction curve to diets (Fig. 1C, diet*genotype interaction p < 0.0001 between wild-type iso31 and ClkJrk), suggesting that Clk is a “true” DR gene. For example, ClkJrk flies only show a 14% increase in lifespan between 5SY and 15SY while iso31 flies show a 34% increase. In addition, while ClkJrk flies are short-lived relative to wild-type at 5SY, they significantly outlived wild-type flies in 1SY (Fig. 1 & S2, p < 0.0005 by log-rank test). Similar results were obtained in an independent trial (Fig. S2). It has been suggested that a whole food dilution may cause dehydration effects, especially on high concentrations of yeast (Ja et al., 2009) and sucrose (van Dam et al., 2020), which may lead to a false conclusion on lifespan response to diet. To eliminate this possibility, we additionally tested the DR response of ClkJrk mutant flies using the yeast-restriction strategy, where varying yeast concentration with a fixed sucrose concentration(Bass et al., 2007; Ja et al., 2009; McCracken et al., 2020). We first tested the DR response of ClkJrk mutant flies using the same experimental diets used for Fig.1 and Fig. S3 except with a fixed sucrose concentration at 5%. Although the reaction pattern of ClkJrk flies in this yeast-restriction DR regimen was qualitatively somewhat different from that of whole food dilution, we confirmed that the DR response of the mutants was strongly suppressed as in whole food dilution protocols (Fig. 1 and Fig. S2). We then further confirmed this observation using a different yeast-restriction protocol where purified yeast extract is used (Katewa et al., 2016; Ulgherait et al., 2016) instead of whole-cell lysates of yeast (McCracken et al., 2020; Min et al., 2007). Although overall mean lifespan and lifespan response to the yeast extract diets (Fig. S4, Table S1) was qualitatively different from those of whole food restriction (Fig. 1 and S2) and whole cell lysates yeast restriction (Fig. S3), DR response of ClkJrk mutant flies in yeast extract diets was significantly impaired compared to wild-type iso31 control flies. Thus, across several diet conditions, these data indicate that ClkJrk robustly suppresses the DR longevity response. While increasing diet concentration has a negative impact on lifespan, it has a strong positive correlation with female fecundity (e.g., egg laying) primarily due to increased protein sources in yeast (Bass et al., 2007; Skorupa et al., 2008). To determine if ClkJrk also impacted this diet response, we measured egg laying in both iso31 and ClkJrk mutant flies over 7 days. As expected, iso31 flies increased egg production with increased diet concentrations (Fig. 2). However, egg production was strongly decreased in ClkJrk mutant flies (Fig. 2, p < 0.0001 by regression analysis). More importantly, diet-dependent increases in egg laying were significantly suppressed in ClkJrk mutants (Fig. 2, p < 0.01 for pair-wise comparison in each diet by t-test). A similar trend was observed in an independent trial (Fig. S5). Taken together, these data indicate that ClkJrk strongly disrupts how flies respond to DR at both levels of longevity and fecundity.
Shorter Term Dietary Restriction Selectively Reprograms Circadian Output Genes in the Abdominal Fat Body
Given the role of ClkJrk in mediating responses to DR, we then asked how the circadian transcriptome responds to DR. A recent RNA-Seq analysis using whole-fly lysates of female flies after seven days of DR or ad libitum diets (a yeast extract restriction protocol) suggested that DR alters the circadian transcriptome in whole flies (Hodge et al., 2022). We focused our studies on the fat body, an analog of the mammalian liver and adipose tissue, given its critical role in mediating the effect of DR (Bai et al., 2012; Banerjee et al., 2012; Dobson et al., 2018; Katewa et al., 2016). Importantly, the fat body also has its own core clock system, including Clk, regulating energy metabolism, feeding, and egg-laying (Xu et al., 2011; Xu et al., 2008). To assess DR-dependent circadian rhythms, we performed a fine scale (every 2 h over 24 h) RNA-Seq analysis from dissected abdominal fat body tissues (DiAngelo and Birnbaum, 2009; Xu et al., 2011; Xu et al., 2008) of young (8 days old) iso31 females. To prepare dissected fat body samples for RNA-Seq analysis, we entrained flies for ~ 5 days under 12h:12h light-dark (12LD) cycles on either control or DR diets. Previous studies showed that ~ 3 days in LD cycles are sufficient to entrain the fat body clock in Drosophila (Erion et al., 2016; Xu et al., 2011; Xu et al., 2008). Moreover, DR reshapes mortality of flies within 2 ~ 4 days of diet shift in flies (Mair et al., 2003; McCracken et al., 2020). We also observed that 5 days on DR diet was sufficient to affect metabolism and physiology evidenced by changes in fecundity (Fig. 2 & S3). Thus, our environmental settings in light schedule and diet are sufficient to capture significant diet- and circadian-dependent transcriptional changes important for lifespan extension by DR. We then analyzed the samples from control and DR diets separately to examine if and how DR changes the circadian gene expression pattern in the fat body. For this diet-dependent analysis, we used BooteJTK (Hutchison et al., 2018) to identify rhythmic genes, using a fold-change threshold ≥ 1.5 at an FDR < 0.25. This cutoff corresponds to p values of 0.011 and 0.014 for control and DR, respectively, making the p value cutoff for this study similar to or more stringent than several recent studies (Abruzzi et al., 2017; Eckel-Mahan et al., 2013; Kuintzle et al., 2017; Sato et al., 2017). We found a significant reorganization of the circadian transcriptome by DR (Fig. 3B-C, Fig. 3E-G). Collectively, we identified 623 oscillating in either one or both conditions, of which 136 cycle in both conditions, 188 cycle only in the control diet, and 299 cycle in the DR diet. Thus, there is a net increase of 50% in cycling genes under DR (Fig. 3B-C). Remarkably, core clock genes were not significantly impacted by DR in both phase and peak expression (using thresholds FDR < 0.1 and log2 (fold change) > 0.5 from differential analysis, see methods) (Fig. 3D). This indicates that rhythmicity of core clock genes is largely resistant to short term (~5 days) diet changes. Although peak phases of common rhythmic genes remained largely unaffected (Fig. 3D), DR significantly, albeit modestly, increased overall expression (10-70% in TPM) of many of these common cyclers, although notably none of the core clock genes were increased (Fig. 3D, F & G). From an averaged expression comparison across all time points between control and DR in the 136 common rhythmic genes, 54 genes were significantly increased (t-test, FDR < 0.05) while only one was downregulated (Fig. 3E-G). It shows that short-term DR (5 days) increases expression of robustly rhythmic genes without affecting expression of the core clock genes, arguing that DR impacts rhythmic output while sparing core clocks.
Weighted Gene Coexpression Network Analysis Identifies a DR-Specific Cycling Proteasome Module
In order to identify novel genes and pathways that are associated with or even causal to the lifespan extension under DR and also to understand circadian transcriptomic organization in the fat body under DR, we took a network approach (Zhang et al., 2013). First, we performed Weighted Gene Coexpression Network Analysis (WGCNA) to reconstruct gene coexpression networks from our time-series RNA-Seq data collected under DR diet (See Methods). We identified 41 network modules of co-regulated genes (Table S2). Notably, 12 modules (29%) were “cycling modules”, i.e., those enriched with rhythmic genes (FDR < 0.05). Genes in each of these “cycling modules” exhibited highly similar phases and waveforms, revealing coordinated circadian gene expression. In order to further examine how the transcriptomic organization is altered by DR, we computed the modular differential connectivity (MDC) (Zhang et al., 2013), which was expressed as a ratio to reflect the difference in gene co-expression strength (i.e., network connectivity) of a module between DR and control diets (see Methods and Supplementary Information). At FDR < 0.05, we identified 14 network modules among the 41 that gained connectivity (MDC > 1) and one network module that lost connectivity (MDC < 1) under DR diet compared to control diet (Table S2). Importantly, three of the differentially connected modules were also cycling modules identified from WGCNA (Fig. S6 A-B), exhibiting higher network connectivity under DR. This suggests that DR may increase the circadian coordination of gene expression in these modules. One of these modules, which we term the “proteasome module”, was of special interest, as the protein products from many of the genes in this module are the core and auxiliary components of the proteasome complex (Fig. 4A-C), showing a dramatically higher pathway enrichment scores than the other two modules (Fig. S6C). This observation also agrees with the primary analysis (Fig. 3B, 3C) that, out of the 33 subunits of the proteasome complex in Drosophila (Belote and Zhong, 2009), one and 12 proteasome subunits were defined as cycling (fold-change in TPM ≥ 1.5 at an FDR < 0.25 BooteJTK analysis) on control and DR diets, respectively. The proteasome complex functions as one of the major proteolytic degradation machines (von Mikecz et al., 2008). Remarkably, 25 of the 33 subunit genes that comprise the proteasome complex (Belote and Zhong, 2009) were found in this cycling module which gained network connectivity under DR, suggesting a DR-specific circadian coordination in the gene expression. In addition to differential connectivity, DR mildly (~20%) but significantly elevated the expression of 21 genes (FDR < 0.05) in the proteasome module (Fig. 4D). Thus, this observation suggests that circadian clocks modulate daily proteasome subunit gene expression in a nutrition-dependent manner.
Knockdown of proteasome subunit genes in fat body limited normal lifespan and the effects of dietary restriction
Our finding of DR-induced changes in proteasome gene expression and cycling suggests a potential mechanism by which the clock could mediate DR effects. Loss of protein homeostasis (proteostasis) is a hallmark of aging and can determine lifespan in both humans and model organisms, including flies (Kaushik and Cuervo, 2015; Koyuncu et al., 2021; Meller and Shalgi, 2021; Santra et al., 2019; Yang et al., 2019; Yu and Hyun, 2021). The proteasome plays a central role in proteostasis by clearing, recycling, and breaking down up to ~ 90% of cellular proteins (Jang, 2018; von Mikecz et al., 2008). In multiple animal models, activation of the proteasome system extends lifespan and healthspan (Anderson et al., 2022; Chondrogianni et al., 2015; Kruegel et al., 2011; Munkacsy et al., 2019; Tonoki et al., 2009; Vilchez et al., 2012). For example, in Drosophila, global overexpression of rpn11, a regulatory subunit of the proteasome complex, extends lifespan while knock-down of rpn11 decreases lifespan (Tonoki et al., 2009). Similarly, adult-specific ubiquitous overexpression of prosϐ5, a catalytic subunit of the proteasome, extends lifespan in flies (Nguyen et al., 2019). Moreover, pharmacological inhibition of the proteasome reduces lifespan in a dosage-dependent manner (Tsakiri et al., 2013). Intriguingly, a recent study in worms shows that DR extends lifespan through promoting proteostasis, suggesting a link between DR and the proteasome (Matai et al., 2019). However, in flies, whether fat body proteasome function is linked to aging and/or the DR response has not been reported. To gain insight into the role of the fat body proteasome, we used the mifepristone (RU486)-inducible GAL4/UAS Gene-Switch (GS) system (Osterwalder et al., 2001; Roman et al., 2001) to knockdown the expression of 11 proteasome subunits predominantly in the adult abdominal fat body using the S106-GS-Gal4 driver (Bai et al., 2012; Jin et al., 2020; Poirier et al., 2008; Roman et al., 2001; Taylor et al., 2022) and assess its effects on lifespan (Fig. S7). Among those genes that displayed the most robust (>25% change) effects on lifespan, we focused on two genes, prosβ3 and rpn7. While we observed some variability among independent trials (Fig. S8), which has previously been observed in longevity assays in Drosophila (Bai et al., 2015; Katewa et al., 2016), we found significant reductions in lifespan as well as suppressions of DR-mediated lifespan extension for knockdown of both of these subunits (p < 0.0001 and p = 0.0222 for prosβ3 and rpn7, respectively for the gene*diet interaction effect from pooled data)(Fig. 5 & S8, Table S1). In the case of prosβ3, three out of the four trials demonstrated statistically significant gene*diet interactions (Table S1). Although the S106-GS-Gal4 is the most widely used inducible fat body driver, it also mis-expresses in the digestive system (Poirier et al., 2008). To test whether the reduced DR response by prosβ3 and rpn7 knockdown with S106-GS-Gal4 is solely from the abdominal fat body or in the intestine or both, we performed two rounds (trial 1 and 2) of adult- and tissue-specific knockdown experiments using the gut-specific inducible TIGS-2 (TIGS-Gal4) (Poirier et al., 2008; Ulgherait et al., 2020). Although the DR effect was weaker in TIGS-Gal4 control flies (without RU486) compared to other wild-type control flies tested in this study (Fig.1, S2, S3, S7, S8), presumably due to genetic background effects (Jin et al., 2020; Liao et al., 2013; Wilson et al., 2020), knockdown of prosβ3 and rpn7 in the gut consistently displayed a stronger lifespan extension by DR (Fig. S9 and Table S1) despite some variation in the lifespan pattern between the two trials. Notably, compared to S106-GS-Gal4 experiment (Fig. 5 and S8), flies with either of the subunits knocked down in the gut suffered much stronger reductions in lifespan regardless of diet types (Fig. S9), implying that the S106-GS-Gal4 results are not significantly affected by leaky misexpression in the gut. Thus, our data indicate that proteasome function in the adult abdominal fat body is important for DR-mediated lifespan extension.
Discussion
While circadian rhythms dampen during aging, environmental and genetic perturbations leading to high-amplitude circadian rhythms correlate with many health benefits (Froy, 2018; He et al., 2016). Here we demonstrate that the master circadian clock transcription factor Clk is important for DR effects on lifespan and fecundity. We also discovered that DR acutely (~ 5 days) alters the circadian transcriptome in the fat body while sparing the core clock suggesting that the primary effects of DR are on circadian output genes. Using adult- and tissue-specific RNAi, we show that diet sensitive clock-controlled proteasome subunit genes in the abdominal fat body are important for the lifespan extending effects of DR. Our data suggest a role of the clock in DR effects on lifespan and reveal a molecular pathway, the proteasome, through which the clock may exert some of these effects.
Our data demonstrate a profound role for the transcription factor Clk in mediating the effects of DR on two independent phenotypes: lifespan and fecundity. The effects of ClkJrk on DR-mediated lifespan expansion were robust, replicable, and, importantly, exhibited across a wide range of diets. Testing DR effects by using just two diets can mask a shift in the diet-dependent lifespan curve (Clancy et al., 2002; Flatt, 2014), for example, due to changes in feeding. The fact that we observe a suppressed response to DR across a wide range of diets suggests a strong DR phenotype. In fact, Clk is one of just a handful of fly genes (Banerjee et al., 2012; Katewa et al., 2016; Ulgherait et al., 2016; Wang et al., 2009; Zid et al., 2009) for which this more rigorous standard has been achieved, suggesting a unique and central role of Clk in DR. ClkJrk mutants are also not simply adversely affected by restricted nutrient intake, as may be the case for Bmal1 mutant mice (Kondratov et al., 2006), as ClkJrk mutants are much more long-lived than iso31 flies under the 1SY and 1Y malnutrition diet (Fig. 1, S2, S3). Female fecundity in Drosophila is strongly correlated with diet concentration (Bass et al., 2007). ClkJrk flies also exhibited a suppressed diet-dependent fecundity response especially at higher diet concentrations (Fig. 2, S5), resulting in a reduced fecundity response to a range of diets. These data suggest a more general role for Clk in dietary sensitivity beyond lifespan. Reduction of fecundity in female ClkJrk mutants is consistent with previous observations under standard diets (Beaver et al., 2002). While we observed reduced DR effects in ClkJrk mutants, we observed a relatively robust DR lifespan response in arrhythmic per01 mutants. Notably, our results reproduced one of the prior reports (Ulgherait et al., 2016) with per01. Loss of the major activator (Clk) and repressor (per) stop the clock at different stages of the cycle (Claridge-Chang et al., 2001; Emery et al., 1998; Glossop et al., 1999). As a result,per and Clk mutants can often yield distinct, even opposing, phenotypes (Keene et al., 2010). While we cannot exclude the possibility that ClkJrk may not act via control of oscillatory gene expression (see (McDonald and Rosbash, 2001)), we favor the idea that the clock drives daily oscillations between DR-sensitive (low PER) and DR-insensitive (low CLK) states which may be adapted to daily feeding (Xu et al., 2008) and/or the sleep/wake rhythm.
In addition to demonstrating a critical role for Clk in mediating the effects of DR on lifespan, we also find that DR reprograms the circadian transcriptome, not by changing core clocks, but rather by inducing or altering rhythmicity of a key set of clock-controlled output genes. In accordance with the recent guidelines (Hughes et al., 2017), we collected samples every 2 h from two independent cycles of 24 h (every 2h for 48 h), increasing the statistical power for rhythm detection, and computed the false discovery rate. To determine how the circadian clock may mediate DR effects we assessed the circadian transcriptome under DR conditions in the fat body, a tissue important for metabolism, longevity, and DR. Importantly we assessed the transcriptome after just 5 days of DR sufficient time for DR to induce changes in mortality rate (Mair et al., 2003; McCracken et al., 2020). We found that this short-term DR was sufficient to produce ~ 35% more rhythmic genes in total compared to control diet. DR also increases overall expression of the rhythmic genes in the control diet that remain rhythmic under DR (common rhythmic genes) (Fig. 3B-F). This is consistent with the observation in mouse liver that DR increases the number of rhythmic genes as well as their amplitude (Sato et al., 2017), indicating this feature of circadian DR sensitivity is widely conserved. Strikingly, although we discovered that DR for 5 days is sufficient to initiate reprogramming of the circadian transcriptome, core clocks remained virtually unaffected (Fig. 3D), a time at which DR mortality effects are observed (Whitaker et al., 2014). In contrast, increased amplitudes of core clocks in mice liver is seen after 2~ 6 months of DR (Patel et al., 2016b; Sato et al., 2017) and in flies after 10 days (Katewa et al., 2016).
Genetic induction of rhythmic amplitude of the core clock gene timeless altered lifespan in a DR-sensitive manner (Katewa et al., 2016). Yet, this tim induction was not accompanied by downstream changes in other components of the feedback loop, suggesting the core clock per se was not involved. Thus, we hypothesize that on the time scale of DR-induced changes in mortality rate, core clocks are not affected and that later changes in core clocks reflect an indirect and delayed effect of DR. It will be of interest to determine if short-term DR in mammals also spares core clock genes. Together, we postulate that core clocks are resistant to amplitude changes by short-term DR regimens while a longer term, depending on species, gradually increases their amplitude, which can further reprogram diet-dependent CCGs. This also implies that the length of DR shapes the pattern of circadian transcriptome reprogramming. Our data suggest a model by which the circadian clock gates the response to dietary restriction to control the complement and amplitude of clock regulated gene expression (Fig. 6). First, the circadian clock rhythmically controls the daily expression of multiple components of the proteasome. Second, DR increases the expression levels and enhances coordinated rhythmic expression of proteasome subunit genes. WGCNA followed by MDC analysis discovered that the transcriptomic organization in the proteasome module was altered by DR. This change is due to enhanced and coordinated rhythmic expression of proteasome subunits by DR (Fig. 4A-C & S6). In addition, DR mildly but significantly increased average expression (t-test with the time-averaged expression analysis, FDR < 0.05) of many of the proteasome subunits (Fig. 4D). In line with the observation by Katewa et al that amplitude changes by DR in core clocks is gradual and requires a minimum of 6 ~ 10 days in a yeast restriction DR (Katewa et al., 2016), we speculate that acute DR-responsive genes/pathways can enhance CLK driven rhythmic processes. Importantly, we provide in vivo evidence that the diet and clock sensitive subunits of the proteasome in the fat body are important for diet sensitive effects on lifespan, providing a pathway of linking clocks, DR, and aging. Although it is well-established in multiple species that proteasome activity decreases during aging and is generally positively correlated with lifespan (Anderson et al., 2022; Chondrogianni et al., 2015; Huang et al., 2019; Kruegel et al., 2011; Nguyen et al., 2019; Pickering et al., 2015; Tonoki et al., 2009; Vilchez et al., 2012), little is known about its tissue-specific contribution to aging and the link to DR. We found that the knock-down of several cycling proteasome subunits predominantly in the adult fat body, using the S106-GS gal4 (Poirier et al., 2008), significantly reduces lifespan, consistent with whole organism manipulations (Fig. S7). Moreover, further testing for potential DR effects with two selected subunits (prosβ3 and rpn7) revealed that knock-down of these subunits reduced DR effects, providing evidence that suppression of proteasome function in the fat body limits DR-mediated lifespan extension (Fig. 5). While some variability was observed, perhaps due to inconsistent delivery of RU486 to flies (Yamada et al., 2017), in the case of prosβ3, significant effects on DR response were observed in three out of four trials. Moreover, combined data from independent trials confirmed knockdown of these subunits reduces the DR effect. We favor the idea that impaired proteasome function contributes to amino acids imbalance, which is known to be critical for DR-mediated lifespan extension in flies (Grandison et al., 2009), leading to reduced DR response (Fig. 6). As DR is suggested to be the most promising intervention to delay aging, the work presented here has important implications for integrating timing into anti-aging therapies. In conjunction with our observations, a recent study also demonstrated that the beneficial effects of time-restricted feeding on longevity is strongly abolished in core clock mutant flies including ClkJrk flies (Ulgherait et al., 2021), emphasizing the roles of circadian clocks in dietary interventions for health and longevity. We propose that time-of-day activation of proteasome may, at least partially, mediate the beneficial effect of DR. Thus, the daily timing of anti-aging therapies may be crucial for lifespan and healthspan extension.
Materials and Methods
Fly Rearing and Media
All the flies used for experiments were raised on a standard yeast-cornmeal-molasses based diet under a light-dark (LD) 12:12h cycle. The following flies were used in this study. ClkJrk and per01 flies were backcrossed to the wild type (w1118) iso31 line (Bloomington stock number: 5905) 6 times. S106-GeneSwitch (S106-GS) (Bloomington stock number: 8151) and RNAi lines were obtained from Bloomington stock center (See Table S1).
Lifespan and Fecundity Assay
For the lifespan assay, young adult female flies (~ 48 hours cohorts) were separated under light CO2 anesthesia after mating with males for ~ 2 days in food bottles. Separated female flies were kept in groups of 20~25 flies in plastic vials on the Sucrose-Yeast (SY) diet and transferred to fresh food vials every 2~3 days. For the lifespan experiment with Mifepristone (RU486)-inducible GeneSwitch system, flies were kept in the vials containing either vehicle (1% EtOH) or RU486 (200 uM). Dead flies were recorded at each transfer. For the fecundity assay, single females (~ 3 days old) were placed with two males of similar age in vials containing different SY diets. Vials were changed daily at ZT0 (8 AM) for 7 days and stored at 4°C until the eggs were counted. Vials with dead females or sterile females (no eggs laid) during the assay were removed from the analysis. For both lifespan and fecundity assays, flies were kept at 25°C, 12hr light : 12 dark (12L:12D) and 60% relative humidity.
Survival Statistics
Survival analysis, log-rank test to evaluate statistical differences between survival curves, and Cox proportional hazards analysis to evaluate ability of tested genes to modify lifespan in the specified range of diets were performed with JMP® statistical package (version 14, SAS Institute Inc.) with data from replicate vials combined. P values from the Cox proportional hazards analysis represent the probability of main (Gene (G) as a nominal variable and Diet (D) as a continuous variable) and interaction effects (G x D) by the likelihood ratio chi-square test.
Fat Body Dissection
Young mated female flies (~ 3 day old) were entrained under either DR diet (5SY) or control diet (15SY) for 5 days in 12L:12D cycles at 25°C. At every 2 hours, flies were directly dissected without dry ice to harvest fat tissues in the abdomen. Pinned flies were cut to remove organs in the abdomen (intestine, ovaries, malpighian tubules, etc). Fat tissues attached to the epidermis were collected (DiAngelo and Birnbaum, 2009; Katewa et al., 2016; Xu et al., 2011; Xu et al., 2008). Fat body from ~ 10 flies were harvested within 10 minutes for each time point of RNA-Seq analysis.
RNA-Seq
Dissected fat body tissues from ~8 days old mated female flies were homogenized in pH 7.4 PBS for 2 min using a Kontes motor and pestle, and were incubated TRIzol LS Reagent (Thermo Fisher, Waltham, MA) for 15 min. RNA was extracted according to the manufacturer’s instructions and residual DNA in the homogenized samples was removed by RNAse free DNase I (Thermo Fisher Scientific). Quality of RNA samples were checked with Agilent 2100 Bioanalyzer. cDNA library was constructed with poly(A) selected mRNA using Truseq RNA library preparation kit and then sequenced at the Genomics Core Facility at the University of Chicago on Illumina HiSeq 2000 System.
Quantification of Transcript, Normalization, and Batch Correction
RNA-seq data were quantified at transcript level using Kallisto (Bray et al., 2016), reference transcriptome used FlyBase_r6.14 (Gramates et al., 2017). Quantified transcripts were summed up to the gene level using tximport library (Soneson et al., 2015). The resulting gene set was initially filtered based on the TPM level; genes with less than 1 TPM across more than 40% of time points per condition (including replicates), were removed from the further analysis. Quantified and filtered samples were normalized within condition with RUVSeq library (Risso et al., 2014), under the RUVg protocol using a set of “in-silico empirical” negative controls (the least significantly DE genes based on a first-pass DE analysis performed prior to RUVg normalization). Technical batch correction between conditions was performed with EDASeq protocol, using upper-quartile (UQ) normalization (Risso et al., 2011).
Rhythmicity Detection
Rhythm detection was performed using BooteJTK (Hutchison et al., 2018) on filtered, TPM level data, with parameters set to: period detection for 24 hr, sampling interval 2hr, model function cosine, phases 00-22hr by 2 hr and asymmetries 02-22hr by 2 hr. Genes with FDR corrected p-value < 0.25 and 1.5 > fold change in raw TPM were assumed to be cycling. Rhythm detection was performed on the commonly expressed genes between analyzed conditions (conditions pooled, 7772 IDs) for detection of the most robust cycling in fat body, additionally condition specific detection was performed (2 replicates per conditions, 7937 in control diet, 7865 in DR diet). Differential gene expression analysis was performed using DESeq2 (Love et al., 2014).
WGCNA and MDC
We reconstructed gene coexpression networks under the DR condition using the WGCNA/r package (Langfelder and Horvath, 2008; Zhang and Horvath, 2005). Briefly, we first computed the network adjacency matrix as kij = [0.5 * (1 – rij)] ^ β. kij is the network connectivity and rij is the Pearson correlation coefficient between a pair of genes i and j. Soft power threshold β was chosen so that the topology of the network was scale-free. A topological overlap matrix (TOM) was then computed to evaluate the neighborhood similarities between genes and to classify network genes into modules using hierarchical clustering and dynamic tree cut. We used DAVID (v6.8) to functionally annotate the identified network modules. To identify network modules that were organized by circadian rhythmicity, we tested the enrichment of rhythmically expressed genes (FDR < 0.1 identified using concatenated data from both DR and control diets) in network modules using fisher’s exact test. Since we defined the network as a “signed” network (i.e., using kij = [0.5 * (1 – rij)] ^ β, instead of the default choice of kij = |rij| ^ β), genes in the cycling network modules shared highly similar circadian phases (as opposed to also including genes with the exact opposite phases) as well as cycling waveforms. To evaluate changes in network connectivity between DR and control diets, we implemented in R a method to compute modular differential connectivity (MDC) described by Zhang et al., 2013 (Zhang et al., 2013). MDC is defined as the ratio of the summed pairwise connectivity among genes in a module under DR conditions and that of the same genes under control conditions (i.e., MDC = ∑i ∑j kijDR / ∑i ∑j kijcontrol). Statistical significance was determined using a permutation-based FDR approach. Two types of FDR estimates were computed, one based on randomly permuted samples to generate networks with nonrandom nodes but random connections, and the other based on randomly permuted gene labels to generate networks with random nodes but nonrandom connections. The final FDR was determined as the larger of the two estimates. 1000 permutations were computed for each type of FDR estimates, and FDRs for modules that gained connectivity (MDC > 1) or lost connectivity (MDC < 1) under DR conditions were estimated separately as described by Zhang et.al., 2013 (Zhang et al., 2013).
Differential Expression Analysis
Differential gene expression analysis was performed using DESeq2 library (Love et al., 2014). We used est. count data generated by kallisto, and a union of pre-filtered gene lists for both conditions (gene lists were in agreement with the TPM level pre-filtering, as described above) resulting in 8030 gene IDs. The full model included time, in a form of 3rd degree polynomial, and diet type as factors. We assumed genes with FDR < 0.2 and absolute log2 fold change > 0.5 as differentially expressed between the two dietary conditions, DR (5SY) and control (15SY) diets. Gene functional classification of differentially expressed genes was performed using DAVID (v6.8) (Huang da et al., 2009).
Acknowledgements
We thank Besim Becoja, Eric Ho Cheung, Alejandra Diaz, Gwang Min Han, Alexandra Raymond, Benjamin Green, and Tova Beeber for technical assistance for lifespan assays and the Bloomington Stock Center for RNAi lines. This work was supported by Defense Advanced Research Projects Agency (DARPA) (D12AP00023 to RA) and the Data Science Initiative (DSI) Research Support Program (to RA) at Northwestern University. This effort was in part sponsored by DARPA; the content of the information does not necessarily reflect the position or the policy of the government, and no official endorsement should be inferred. This work was supported by the NSF-Simons Center for Quantitative Biology at Northwestern University. This work was supported by a grant from the Simons Foundation (597491-RWC) and the National Science Foundation (1764421). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Science Foundation and Simons Foundation. DSH was supported by the National Institutes of Health T32 Institutional Training Grant (Northwestern Univ: grant NIH T32HL007909) and the National Institute of General Medical Sciences of the National Institutes of Health under Award Number P20GM103436 (Univ. of Louisville). ALH was supported by the National Institutes of Health Medical Scientist Training program at the University of Chicago (grant NIGMS T32GM07281).
Declaration of Interests
The authors declare no competing interests.
References
- RNA-seq analysis of Drosophila clock and non-clock neurons reveals neuron-specific cycling and novel candidate neuropeptidesPLoS Genet 13
- Mice under Caloric Restriction Self-Impose a Temporal Restriction of Food Intake as Revealed by an Automated Feeder SystemCell Metab 26:267–277
- Circadian organization of behavior and physiology in DrosophilaAnnu Rev Physiol 72:605–624
- A mutant Drosophila homolog of mammalian Clock disrupts circadian rhythms and transcription of period and timelessCell 93:791–804
- Hyperactivation of the proteasome in Caenorhabditis elegans protects against proteotoxic stress and extends lifespanJ Biol Chem 298
- Drosophila insulin-like peptide-6 (dilp6) expression from fat body extends lifespan and represses secretion of Drosophila insulin-like peptide-2 from the brainAging cell 11:978–985
- Drosophila Longevity Assurance Conferred by Reduced Insulin Receptor Substrate Chico Partially Requires d4eBPPloS one 10
- dSir2 in the adult fat body, but not in muscles, regulates life span in a diet-dependent mannerCell Rep 2:1485–1491
- Optimization of dietary restriction protocols in DrosophilaJ Gerontol A Biol Sci Med Sci 62:1071–1081
- Loss of circadian clock function decreases reproductive fitness in males of Drosophila melanogasterProc Natl Acad Sci U S A 99:2134–2139
- Duplicated proteasome subunit genes in Drosophila and their roles in spermatogenesisHeredity (Edinb) 103:23–31
- Near-optimal probabilistic RNA-seq quantificationNat Biotechnol 34:525–527
- 20S proteasome activation promotes life span extension and resistance to proteotoxicity in Caenorhabditis elegansFASEB J 29:611–622
- Dietary restriction in long-lived dwarf fliesScience 296
- Circadian regulation of gene expression systems in the Drosophila headNeuron 32:657–671
- Regulation of fat cell mass by insulin in Drosophila melanogasterMol Cell Biol 29:6341–6352
- Tissue-specific transcriptome profiling of Drosophila reveals roles for GATA transcription factors in longevity by dietary restrictionNPJ Aging Mech Dis 4
- Deficiency of circadian protein CLOCK reduces lifespan and increases age-related cataract development in miceAging (Albany NY) 2:936–944
- Reprogramming of the circadian clock by nutritional challengeCell 155:1464–1478
- CRY, a Drosophila clock and light-regulated cryptochrome, is a major contributor to circadian rhythm resetting and photosensitivityCell 95:669–679
- Neural clocks and Neuropeptide F/Y regulate circadian gene expression in a peripheral metabolic tissueElife 5
- Plasticity of lifespan: a reaction norm perspectiveProc Nutr Soc 73:532–542
- Circadian rhythms, nutrition and implications for longevity in urban environmentsProc Nutr Soc 77:216–222
- The circadian gene Period2 plays an important role in tumor suppression and DNA damage response in vivoCell 111:41–50
- Dietary restriction and other lifespan extending pathways converge at the activation of the downstream effector takeoutAging (Albany NY) 2:387–389
- Time-restricted feeding attenuates age-related cardiac decline in DrosophilaScience (New York, N.Y.) 347:1265–1269
- Interlocked feedback loops within the Drosophila circadian oscillatorScience 286:766–768
- FlyBase at 25: looking to the futureNucleic Acids Res 45:D663–D671
- Amino-acid imbalance explains extension of lifespan by dietary restriction in DrosophilaNature 462:1061–1064
- Molecular mechanisms of dietary restriction promoting health and longevityNat Rev Mol Cell Biol 23:56–73
- Time-restricted feeding without reducing caloric intake prevents metabolic diseases in mice fed a high-fat dietCell Metab 15:848–860
- The Small Molecule Nobiletin Targets the Molecular Oscillator to Enhance Circadian Rhythms and Protect against Metabolic SyndromeCell Metab 23:610–621
- Dietary restriction and the transcription factor clock delay eye aging to extend lifespan in Drosophila MelanogasterNat Commun 13
- Systematic and integrative analysis of large gene lists using DAVID bioinformatics resourcesNat Protoc 4:44–57
- RiboTag translatomic profiling of Drosophila oenocytes under aging and induced oxidative stressBMC Genomics 20
- Guidelines for Genome-Scale Analysis of Biological RhythmsJ Biol Rhythms 32:380–393
- Bootstrapping and Empirical Bayes Methods Improve Rhythm Detection in Sparsely Sampled DataJ Biol Rhythms 33:339–349
- Water-and nutrient-dependent effects of dietary restriction on Drosophila lifespanProc Natl Acad Sci U S A 106:18633–18637
- Regulation of Protein Degradation by Proteasomes in CancerJ Cancer Prev 23:153–161
- Genetic and metabolomic architecture of variation in diet restriction-mediated lifespan extension in DrosophilaPLoS Genet 16
- Increased transsulfuration mediates longevity and dietary restriction in DrosophilaProc Natl Acad Sci U S A 108:16831–16836
- Peripheral Circadian Clocks Mediate Dietary Restriction-Dependent Changes in Lifespan and Fat Metabolism in DrosophilaCell Metab 23:143–154
- Proteostasis and agingNat Med 21:1406–1415
- Clock and cycle limit starvation-induced sleep loss in DrosophilaCurr Biol 20:1209–1215
- Early aging and age-related pathologies in mice deficient in BMAL1, the core componentof the circadian clockGenes & development 20:1868–1873
- Rewiring of the ubiquitinated proteome determines ageing in Celegans. Nature 596:285–290
- Elevated proteasome capacity extends replicative lifespan in Saccharomyces cerevisiaePLoS Genet 7
- Circadian deep sequencing reveals stress-response genes that adopt robust rhythmic expression during agingNat Commun 8
- WGCNA: an R package for weighted correlation network analysisBMC Bioinformatics 9
- Disrupting circadian homeostasis of sympathetic signaling promotes tumor development in micePloS one 5
- Genetic variation in responses to dietary restriction--an unbiased tool for hypothesis testingExp Gerontol 48:1025–1029
- Deviation of innate circadian period from 24 h reduces longevity in miceAging cell 11:794–800
- Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2Genome Biol 15
- Sex differences in the effect of dietary restriction on life span and mortality rates in female and male Drosophila melanogasterJ Gerontol A Biol Sci Med Sci 59:3–9
- Demography of dietary restriction and death in DrosophilaScience 301:1731–1733
- Circadian rhythms, time-restricted feeding, and healthy agingAgeing Res Rev 39:59–67
- Dietary restriction improves proteostasis and increases life span through endoplasmic reticulum hormesisProc Natl Acad Sci U S A 116:17383–17392
- Dietary restriction and ageing: Recent evolutionary perspectivesMech Ageing Dev 208
- The hidden costs of dietary restriction: Implications for its evolutionary and mechanistic originsSci Adv 6
- Microarray analysis and organization of circadian gene expression in DrosophilaCell 107:567–578
- The aging proteostasis decline: From nematode to humanExp Cell Res 399
- Counting calories in Drosophila diet restrictionExp Gerontol 42:247–251
- Neuronal-specific proteasome augmentation via Prosbeta5 overexpression extends lifespan and reduces age-related cognitive declineAging cell 18
- Molecular connections between circadian clock and health/ageingJ Biochem 171:473–476
- Proteasome beta5 subunit overexpression improves proteostasis during aging and extends lifespan in Drosophila melanogasterSci Rep 9
- A conditional tissue-specific transgene expression system using inducible GAL4Proc Natl Acad Sci U S A 98:12596–12601
- Circadian clocks govern calorie restriction-mediated life span extension through BMAL1-and IGF-1-dependent mechanismsFASEB J 30:1634–1642
- Calorie restriction regulates circadian clock gene expression through BMAL1 dependent and independent mechanismsSci Rep 6
- Molecular mechanisms and physiological importance of circadian rhythmsNature Reviews Molecular Cell Biology 21:67–84
- Lifespan of mice and primates correlates with immunoproteasome expressionJ Clin Invest 125:2059–2068
- Characterization of the Drosophila gene-switch system in aging studies: a cautionary taleAging cell 7:758–770
- Normalization of RNA-seq data using factor analysis of control genes or samplesNat Biotechnol 32:896–902
- GC-content normalization for RNA-Seq dataBMC Bioinformatics 12
- P[Switch], a system for spatial and temporal control of gene expression in Drosophila melanogasterProc Natl Acad Sci U S A 98:12602–12607
- Proteostasis collapse is a driver of cell aging and deathProc Natl Acad Sci U S A 116:22173–22178
- Circadian Reprogramming in the Liver Identifies Metabolic Pathways of AgingCell 170:664–677
- Dietary composition specifies consumption, obesity, and lifespan in Drosophila melanogasterAging cell 7:478–490
- Circadian clock genes’ overexpression in Drosophila alters diet impact on lifespanBiogerontology 20:159–170
- Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferencesF1000Res 4
- Diet restriction in Drosophila melanogasterDesign and analysis. Interdisciplinary topics in gerontology 35:115–136
- Sirt6 regulates lifespan in Drosophila melanogasterProc Natl Acad Sci U S A 119
- Genetic evidence linking age-dependent attenuation of the 26S proteasome with the aging processMol Cell Biol 29:1095–1106
- Proteasome dysfunction in Drosophila signals to an Nrf2-dependent regulatory circuit aiming to restore proteostasis and prevent premature agingAging cell 12:802–813
- Circadian regulation of mitochondrial uncoupling and lifespanNat Commun 11
- Dietary Restriction Extends the Lifespan of Circadian Mutants tim and perCell Metab 24:763–764
- Circadian autophagy drives iTRF-mediated longevityNature 598:353–358
- Sugar-Induced Obesity and Insulin Resistance Are Uncoupled from Shortened Survival in DrosophilaCell Metab 31:710–725
- RPN-6 determines Celegans longevity under proteotoxic stress conditions. Nature 489:263–268
- The nuclear ubiquitin-proteasome system: visualization of proteasomes, protein aggregates, and proteolysis in the cell nucleusMethods Mol Biol 463:191–202
- Long-lived Indy and calorie restriction interact to extend life spanProc Natl Acad Sci U S A 106:9262–9267
- Dietary switch reveals fast coordinated gene expression changes in Drosophila melanogasterAging (Albany NY) 6:355–368
- GWAS for Lifespan and Decline in Climbing Ability in Flies upon Dietary Restriction Reveal decima as a Mediator of Insulin-like Peptide ProductionCurr Biol 30:2749–2760
- The circadian clock interacts with metabolic physiology to influence reproductive fitnessCell Metab 13:639–654
- Regulation of feeding and metabolism by neuronal and peripheral clocks in DrosophilaCell Metab 8:289–300
- Mifepristone Reduces Food Palatability and Affects Drosophila Feeding and LifespanJ Gerontol A Biol Sci Med Sci 72:173–180
- Timing of expression of the core clock gene Bmal1 influences its effects on aging and survivalSci Transl Med 8
- Multidimensional Proteomics Identifies Declines in Protein Homeostasis and Mitochondria as Early Signals for Normal Aging and Age-associated Disease in DrosophilaMol Cell Proteomics 18:2078–2088
- Proteostasis-associated aging: lessons from a Drosophila modelGenes Genomics 43:1–9
- Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s diseaseCell 153:707–720
- A general framework for weighted gene co-expression network analysisStat Appl Genet Mol Biol 4
- Deciphering clock genes as emerging targets against agingAgeing Res Rev 81
- 4E-BP extends lifespan upon dietary restriction by enhancing mitochondrial activity in DrosophilaCell 139:149–160
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Copyright
© 2023, Hwangbo et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 612
- downloads
- 18
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.