- Research article
- Open Access
Chronic circadian misalignment results in reduced longevity and large-scale changes in gene expression in Drosophila
© The Author(s). 2019
- Received: 19 September 2018
- Accepted: 20 December 2018
- Published: 7 January 2019
Circadian clocks are found in nearly all organisms, from bacteria to mammals, and ensure that behavioral and physiological processes occur at optimal times of day and in the correct temporal order. It is becoming increasingly clear that chronic circadian misalignment (CCM), such as occurs in shift workers or as a result of aberrant sleeping and eating schedules common to modern society, has profound metabolic and cognitive consequences, but the proximate mechanisms connecting CCM with reduced organismal health are unknown. Furthermore, it has been difficult to disentangle whether the health effects are directly induced by misalignment or are secondary to the alterations in sleep and activity levels that commonly occur with CCM. Here, we investigated the consequences of CCM in the powerful model system of the fruit fly, Drosophila melanogaster. We subjected flies to daily 4-h phase delays in the light-dark schedule and used the Drosophila Activity Monitoring (DAM) system to continuously track locomotor activity and sleep while simultaneously monitoring fly lifespan.
Consistent with previous results, we find that exposing flies to CCM leads to a ~ 15% reduction in median lifespan in both male and female flies. Importantly, we demonstrate that the reduced longevity occurs independent of changes in overall sleep or activity. To uncover potential molecular mechanisms of CCM-induced reduction in lifespan, we conducted whole body RNA-sequencing to assess differences in gene transcription between control and misaligned flies. CCM caused progressive, large-scale changes in gene expression characterized by upregulation of genes involved in response to toxic substances, aging and oxidative stress, and downregulation of genes involved in regulation of development and differentiation, gene expression and biosynthesis.
Many of these gene expression changes mimic those that occur during natural aging, consistent with the idea that CCM results in premature organismal decline, however, we found that genes involved in lipid metabolism are overrepresented among those that are differentially regulated by CCM and aging. This category of genes is also among the earliest to exhibit CCM-induced changes in expression, thus highlighting altered lipid metabolism as a potentially important mediator of the negative health consequences of CCM.
- Circadian misalignment
The daily rotation of the Earth on its axis produces predictable cycles of light and temperature. This has led to the evolution of endogenous circadian clocks, which allow organisms to anticipate and adapt to environmental cycles to ensure that behavioral and physiological processes occur at optimal times of day. Perhaps equally importantly, circadian clocks orchestrate the temporal sequencing of internal events such that they occur at appropriate times with respect to one another. Together, these extrinsic and intrinsic advantages of circadian clocks are thought to promote organismal health . For example, the circadian system generates rhythms in foraging and feeding behaviors such that they occur during times of maximal food availability, and concurrently upregulates metabolic pathways in anticipation of increased food intake.
Research undertaken over the past ~ 50 years has identified many of the genes and molecules that confer upon cells the ability to keep time. In most organisms, the molecular circadian clock consists of a transcriptional-translational feedback loop that produces rhythmic oscillations in gene expression, and this molecular clock is present both in master circadian pacemaker cells in the brain as well as in cells in peripheral tissues [2, 3]. A key feature of circadian clocks is their ability to maintain rhythmic oscillations under constant conditions; however, because they oscillate with an intrinsic period that is not exactly 24-h, they must be reset on a daily basis through the process of entrainment in order to remain synchronized to environmental cues . More extensive clock resetting is necessary following abrupt shifts in lighting schedules, such as occur with trans-meridian travel, shiftwork, or simply as a result of the aberrant eating and sleeping schedules that are common in modern society. The latter in particular, which has been termed ‘social jetlag’, has led to a situation in which as much as 2/3 of the population live under a state of chronic circadian misalignment (CCM) typified by a behavioral schedule that runs outside the dictates of their endogenous circadian rhythms [5, 6]. These irregular schedules likely result in a misalignment between circadian clocks and external environmental signals as well as internal desynchrony between different clock-containing tissues in the body.
Evidence suggests that circadian disruption has severe consequences. Epidemiological studies in humans have demonstrated that various forms of CCM are correlated with higher incidences of a number of disorders including obesity, diabetes, cardiovascular disease, cancer, and mood disorders, along with increased overall mortality rates. [4, 7, 8]. Controlled studies using animal models have corroborated these findings and demonstrated a causal link between CCM and decreased health and longevity, both in invertebrate and vertebrate systems [7–17]. However, despite a growing knowledge of the importance of circadian function for organismal health, the cellular and molecular mechanisms through which CCM negatively impacts health and longevity have yet to be established. Furthermore, it is unclear whether the negative health effects are directly caused by misalignment or whether they are secondary to behavioral changes, such as reduced sleep, that often accompany CCM .
Here we employed a model of CCM in the fruit fly, Drosophila melanogaster, in which we exposed flies to daily 4-h phase delays in lighting conditions. We used the Drosophila Activity Monitoring (DAM) system to continuously monitor fly locomotor activity and sleep while simultaneously conducting longevity analysis. This allowed us to assess the effect of long-term CCM on aging-associated changes in locomotor activity and sleep levels, and to correlate these measures with fly lifespan. We demonstrate that CCM leads to an ~ 15% reduction in median lifespan of both male and female flies, in agreement with previous studies [10, 12]. While misaligned flies exhibited aberrant patterns of locomotor activity, as evidenced by reduced rest:activity rhythm strength under light-dark (LD) conditions, overall sleep and activity levels were largely unchanged and did not account for the reduction in longevity produced by CCM. We further conducted RNA sequencing analysis and identified several categories of genes that are up- or downregulated by CCM. Interestingly, many of the gene expression changes induced by CCM mirror those observed with natural aging, suggesting that CCM results in cellular stress and premature organismal decline. However, we find that genes involved in lipid metabolism, which are among the earliest to undergo CCM-induced changes in expression, are also differentially regulated by CCM and aging, highlighting this category of genes as a potentially important mediator of the negative consequences of circadian disruption.
CCM reduces lifespan
It should be noted that all flies were maintained on minimal sucrose-agar food during longevity analysis, resulting in a relatively short lifespan of even control flies in these experiments compared with those performed with standard food. This food source is typical for use with the DAM system and thus allowed us to more directly compare activity and sleep levels in our experiments to those observed in other studies using this system. Nevertheless, the ~ 15% reduction in median lifespan we observed agrees well with results obtained with flies kept on standard food , which demonstrates that the negative health impacts of CCM are similar regardless of diet.
CCM causes aberrant locomotor activity patterns
Despite the fact that rest:activity rhythms remained time-locked to phase-delaying schedules, we did observe activity patterns in CCM-exposed flies indicative of mistimed locomotor activity. In flies, dawn and dusk lighting transitions are normally preceded by anticipatory increases in locomotor activity that depend on a functional circadian clock . In young control flies, we found that morning and evening anticipatory activity began ~ 8–9 h following the previous lighting transition, with more prominent activity increases occurring in anticipation of lights turning off (Fig. 2a and c). CCM flies also began to ramp up their activity ~ 8–9 h following the previous lighting transition; however, because of the delayed lighting schedule, this led to a substantial mistiming of anticipatory activity which resulted in prolonged activity peaks that started to resolve prior to the actual lighting transitions (Fig. 2b and d). This mistimed anticipation of lighting conditions persisted throughout the duration of the experiment. Thus, in weeks 2 and 3, although all flies displayed lower amplitude activity rhythms, CCM-exposed flies continued to show prolonged activity peaks that preceded lighting transitions by several hours more than controls (Fig. 2a-d). In addition to this prolonged anticipatory activity preceding lighting changes, CCM flies exhibited a modest mid-afternoon activity peak during the normal siesta time. This was most prominently observed in males during the second week of CCM exposure (Fig. 2d). This increased daytime activity may account in part for the increase in week 2 and lifetime activity levels in male CCM flies compared to controls. In contrast, overall activity was unchanged in CCM females (Fig. 2e-f). Together, these results demonstrate that circadian clocks in CCM-exposed flies have to be constantly adjusted to changing environmental cycles, leading to chronically mistimed activity.
To determine whether CCM affects rhythm amplitude in addition to timing of activity, we conducted chi2 periodogram analysis, which demonstrated that CCM-exposed flies displayed reduced locomotor activity rhythm strength compared to controls (Fig. 3a and b). Both experimental and control flies showed expected reductions in rest:activity rhythm strength as they aged  (Fig. 3a-b), typified by less robust activity peaks and reduced anticipatory activity, including an apparent lack of lights-on anticipation in control flies by week 3 (Fig. 2a and c). However, CCM flies had lower rhythm power than controls at all time points, perhaps as a result of mistimed morning and evening anticipation relative to lights-on and lights-off, leading to prolonged activity bouts.
This reduced rhythmicity could be due either to the misalignment between internal and external rhythms or to CCM-induced damage to core clock neurons or molecular cycling. To test for the latter, we assessed locomotor rhythmicity of flies in constant darkness (DD) following exposure to varying amounts of time in either CCM or control conditions (Fig. 3e-h). Importantly, we saw no differences in free-running rhythm power or period between control and experimental flies. Furthermore, CCM flies immediately adopted normal free-running periods upon transition to DD conditions, showing no residual effect of their prolonged 28-h day exposure (Fig. 3i). The fact that CCM flies have normal free-running rhythms indicates that the central clock and associated output pathways maintain proper function following CCM exposure and further suggests that the reduced rhythm strength observed under LD conditions is a result of the discrepancy between the endogenous period of the fly and environmental cycles.
CCM decreases longevity independent of changes in locomotor activity or sleep
Effect of CCM and Sleep on Survival
Effect of CCM on Survival in Univariate Cox Proportional Hazard Model
p < 0.001 ***
Effect of Sleep Time on Survival in Univariate Cox Proportional Hazard Model
Sleep (Per Hour)
p < 0.001 ***
Effect of CCM and Sleep Time on Survival in Multivariate Cox Proportional Hazard Model
p < 0.01 **
Sleep (Per Hour)
p < 0.001 ***
Effect of CCM on Survival in Sleep-Matched Univariate Cox Proportional Hazard Model
p < 0.05 *
Effect of CCM on Survival in Univariate Cox Proportional Hazard Model
p < 0.001 ***
Effect of Sleep Time on Survival in Univariate Cox Proportional Hazard Model
Sleep (Per Hour)
p = 0.32
Because CCM shortens daily sleep time, we wanted to know if the effect of CCM on survival was due to the shortened sleep of these flies or if CCM was more detrimental than would be predicted by its effects on sleep time alone. We took two approaches. In one approach, we used a multivariate Cox Proportional Hazard Model with both treatment and week 1 sleep time in the model. The hazard ratio for CCM in the multivariate model was 1.49 (Table 1), indicating that misaligned flies were roughly 49% more likely to die than controls flies at a given point in time, even with sleep time held constant. Put another way, the risk of death for a CCM fly was roughly equivalent to the risk of death for a control fly with an hour and a half less sleep. In the other approach, we matched each control fly with a CCM fly with a similar sleep time and calculated a hazard ratio for CCM in this smaller, sleep-matched data set, which importantly also normalized overall activity levels in CCM flies. This produced a similar result, with a hazard ratio of 1.37 for CCM (Table 1). In female flies, the estimated hazard ratio for CCM in a univariate model was 1.65 (Table 1), although we note that the model assumption of proportional hazards is not upheld for female flies, as CCM has a stronger effect on survival at earlier time points. Sleep time did not significantly affect survival in female flies, so was not a potential confounding factor. It is therefore clear that the effect of CCM on survival in female flies is not due to its effects on sleep time. Together, these analyses demonstrate that the reduced longevity resulting from CCM is independent of the minor changes in activity and sleep that are associated with the daily 4-h phase delay schedule. This idea is further supported by the fact that CCM equivalently reduced lifespan in male and female flies (Fig. 1a and b) despite the fact that males and females have differences in baseline sleep and that CCM more strongly reduced sleep amount in males than females (Fig. 4e and f).
CCM induces large-scale changes in gene expression
Having determined that CCM reduces lifespan independent of overt behavioral changes, we sought to identify changes at the molecular level that could underlie the physiological consequences of CCM. We therefore used whole-fly RNA sequencing to assess changes in gene transcription associated with prolonged circadian misalignment. For these experiments, we pooled male and female flies, which we reasoned would allow us to identify central, common mechanisms through which circadian misalignment reduces lifespan, even though this may obscure sex-specific gene expression changes.
Although expression changes in week 3 downregulated genes as a whole were not correlated with week 2 changes, there were several individual genes that were similarly affected following different CCM exposure times. Notably, 3/6 genes found to be downregulated following 2 weeks of CCM exposure were also identified as significantly differentially expressed in week 3, providing independent confirmation of these changes. The other 3 week 2 downregulated genes showed a trend towards decreased expression in week 3 samples; however, these genes also showed substantial age-related decreases in expression in control flies that somewhat mitigated the expression differences between control and CCM flies at the week 3 time point (Fig. 6c and d).
CCM-induced gene expression changes indicate cellular stress and organismal decline
Overrepresented GO Terms Associated with Week 2 Differentially Expressed Genes
Number of Genes
fatty acid metabolic process
organic acid metabolic process
lipid metabolic process
organic acid biosynthetic process
Overrepresented GO Terms Associated with Week 3 Differentially Expressed Genes
Number of Genes
--Terms Associated with Downregulated Genes--
cellular developmental process
multicellular organism development
animal organ development
nervous system development
regulation of gene expression
cellular component morphogenesis
regulation of biosynthetic process
pattern specification process
regulation of cellular component organization
--Terms Associated with Upregulated Genes--
response to toxic substance
cofactor metabolic process
sulfur compound metabolic process
chaperone-mediated protein folding
determination of adult lifespan
chaperone cofactor-dependent protein refolding
response to DDT
response to heat
response to oxidative stress
response to temperature stimulus
cellular response to unfolded protein
Overrepresented Nervous System Related GO Terms Associated with Week 3 Downregulated Genes
Number of Genes
nervous system development
learning or memory
central nervous system development
chemical synaptic transmission
GO term analysis of week 3 upregulated genes suggested that CCM-exposed flies were experiencing high levels of cellular stress. We identified 17 overrepresented GO terms associated with week 3 upregulated genes, and these included functions such as response to toxic substances, unfolded protein response, response to oxidative stress, and determination of lifespan (Table 3 and Additional file 6). These findings demonstrating increased expression of stress-related genes, together with the downregulation of genes involved in developmental and biosynthetic processes, are consistent with the idea that CCM accelerates organismal decline.
In line with this, we noted that the majority of differentially expressed genes following 3 weeks of CCM exposure were also found to be similarly regulated by aging in control samples. Thus, 145/245 CCM downregulated genes and 50/104 CCM upregulated genes were similarly affected by aging (using a Benjamini-Hochberg FDR cutoff of 0.1 to identify genes differentially expressed with aging). Furthermore, an additional 50 CCM downregulated genes and 22 CCM upregulated showed a trend towards aging-related changes in the same direction (genes in which aging-related log2FC was in the same direction as the CCM-induced log2FC and had a Benjamini-Hochberg adjusted p-value between 0.1 and 0.3) (Additional file 2). This indicates that CCM may accelerate or exacerbate normal age-related changes in gene expression.
Genes Oppositely Regulated by Circadian Misalignment and Aging
week 3 log2FC
week 3 p-adj
--Downregulated with CCM; Upregulated with Aging--
dipeptidyl-peptidase activity; serine-type peptidase activity
fatty acid elongase activity
carbohydrate:proton symporter activity; glucose transmembrane transporter activity
urate oxidase activity
carbohydrate:proton symporter activity; glucose transmembrane transporter activity
transmembrane transporter activity
fatty acid elongase activity
--Upregulated with CCM; Downregulated with Aging--
asparagine synthase (glutamine-hydrolyzing) activity; protein homodimerization activity
We have chronically exposed flies to aberrant lighting schedules in order to better characterize the behavioral, physiological and molecular consequences of circadian misalignment. Our findings of reduced longevity demonstrate the profound consequences of CCM, and are in line with a number of other studies conducted in flies and mammals. In one of the earliest studies to test the functional significance of circadian rhythms in animals, Pittendrigh and Minis demonstrated that exposure of fruit flies to long (27-h) or short (21-h) days (the equivalent of daily 3-h phase delays or advances) resulted in reduced lifespan , and these findings were later corroborated by similar studies in blowflies . In mammalian models, as in invertebrate systems, studies have consistently demonstrated negative health consequences of CCM, particularly for aged or disease-prone animals. Chronic misalignment has been found to increase mortality in aged mice , immune-challenged mice  and mice and hamsters prone to cardiovascular disease  and cancer . In addition to these studies demonstrating that CCM can exacerbate or hasten disease progression, misalignment has also been found to decrease longevity in otherwise healthy mammals. For example, reduced lifespan was observed in mice exposed to an extremely short 8-h light-dark cycle, an effect that was absent in mutant mice that lacked a functional molecular clock . More recently, it was demonstrated that wildtype mice exposed to weekly phase shifts exhibited signs of neurodegeneration, kidney disease, ulcerative dermatitis, and cancer, and these were associated with reduced lifespan . It should be noted that CCM paradigms in mammals have typically involved less frequent phase shifts than paradigms used in flies. It is likely that the severity of the consequences of CCM will vary depending on the extent and duration of misalignment.
The effects of exposing organisms to non-24-h light cycles provided the basis for Pittendrigh’s theory of circadian resonance, which stated that organisms function best when the difference between endogenous circadian period and environmental cycles is small, thus minimizing the magnitude of daily circadian clock adjustments . Other studies have directly tested the resonance theory by showing that the negative consequences of exposing animals to non-standard length days were at least partially ameliorated in mutant animals whose endogenous period matched that of the aberrant lighting schedule [12, 27–29]. In general, these studies demonstrate that large daily adjustments to the circadian clock are detrimental to organismal health.
CCM could negatively influence organismal health for several reasons, which are not mutually exclusive. First, misalignment could result in a loss of the extrinsic advantage of circadian rhythms to properly time behavioral and physiological processes with respect to external environmental conditions. This possibility is supported by studies of wild-caught chipmunks, in which surgical destruction of the central circadian clock led to premature death, likely due to mistimed activity leading to higher rates of predation . However, in controlled laboratory experiments such as those we have conducted, which include ad libitum food availability and a lack of predation, it is unlikely that extrinsic factors play a prominent role in producing adverse health effects such as decreased longevity. One exception, however, is that some physiological processes may be optimized to occur in either the light or dark. If this were the case, then we would expect our CCM flies to be negatively affected since they show evidence of mistimed behavior, particularly in the hours preceding lighting transitions, which suggests that other physiological processes may be similarly mistimed.
A second possibility is that CCM causes internal desynchrony such that physiological and behavioral processes become disordered with respect to one another. This is thought to be a major cause of jetlag symptoms in mammals, as studies have shown that clocks in different tissues resynchronize to phase shifts at different rates . Furthermore, studies in humans have shown that gene expression rhythms in peripheral tissues quickly dampen under conditions of imposed circadian misalignment, likely due to a loss of synchronization between behavioral schedules and the central circadian clock [32, 33]. It is unclear to what extent this type of internal desynchrony occurs in flies. Our locomotor activity tracking demonstrated apparent entrainment of rest:activity rhythms to the 28-h day produced by daily 4-h phase delays; however, further studies are necessary to assay clock protein cycling in central and peripheral tissues under these conditions to determine whether molecular clocks are indeed entrained to long or short environmental cycles and to test whether central and peripheral oscillators remain synchronized to one another. Even if entrainment is occurring, the substantial discrepancy between the 28-h environmental cycles and the ~ 24-h endogenous period of the fly may be challenging the efficiency of the endogenous clock in CCM flies, as substantial daily adjustment to the phase of the clock would be required. Furthermore, CCM conditions could still result in transient desynchrony occurring on a daily basis, especially towards the end of the extended dark period because clock cycling cannot be adjusted by entrainment mechanisms until subsequent light exposure.
Circadian misalignment could also indirectly affect health secondary to changes in overall activity or sleep levels, which are known to accompany CCM [4, 7, 8]. In general, we found only very modest effects of CCM on sleep duration in flies. Importantly, there was no correlation between sleep amount and longevity in female flies, which indicates that the slightly reduced sleep exhibited by CCM females in the first week of the experiment could not account for their reduced longevity. Furthermore, although sleep amount and longevity were correlated in males, multivariate models demonstrated that CCM decreased longevity to a further extent than what would be expected based on the effect of CCM on sleep. Thus, while sleep disturbances undoubtedly contribute to the negative consequences associated with CCM in humans, our results clearly demonstrate that CCM can reduce lifespan independent of its effects on sleep.
To better understand the proximal causes of the adverse health effects of CCM, we conducted RNA sequencing of flies following 2 or 3 weeks of circadian misalignment. Notably, the 2-week time point occurs prior to most CCM-induced death, which we reasoned would allow us to identify early molecular changes directly associated with CCM as opposed to nonspecific changes occurring as a consequence of CCM-induced organismal decline. Consistent with the relatively healthy state of CCM-exposed flies at 2 weeks, we found few genes that exhibited significant up- or downregulation at this time. Nevertheless, analysis of differentially expressed genes demonstrated an overrepresentation of those associated with lipid metabolism. The relatively early expression changes in these genes suggests that alterations in lipid metabolism could underlie the initial consequences of CCM.
In addition to investigating gene expression changes associated with CCM, we analyzed control samples to identify genes regulated by aging independent of circadian misalignment. Notably, we found that many more genes were differentially expressed in response to aging compared to CCM. The aging-related genes we observed included many known to be upregulated with aging. For example, we found that the heat-shock genes hsp22 and hsp70 were upregulated with aging and were both among the top 5 differentially expressed genes in our analysis. These genes have been found previously to be robustly upregulated with aging in fly heads and whole flies [25, 34].
In order to identify gene expression changes specific to CCM exposure, we isolated genes in which the consequences of CCM were opposite those of aging. Interestingly, genes involved in lipid binding and metabolism were also overrepresented among this dataset, further pointing towards an important contribution of altered lipid homeostasis in mediating the negative effects of CCM. Given the wide range of functions associated with lipid metabolism, more focused experiments will be necessary to determine the specific contribution of the genes implicated in our studies to the negative consequences of CCM. However, it is of note that several candidate genes are involved in fatty acid elongation, which is essential for the creation of very long-chain fatty acids and sphingolipids that make up integral parts of cell membranes and have been implicated in diverse processes including cell growth and cell signaling . This includes CG18609, which is a fatty acid elongase that we found to be significantly downregulated following both 2 and 3 weeks of CCM. CG18609 normally increases in expression as flies age, which indicates that the gene expression changes are specific to CCM and are not simply occurring alongside CCM-induced organismal decline. Interestingly, lipid metabolism has been linked to aging processes . Furthermore, genes involved in lipid metabolism were also recently found to be acutely dysregulated in human skeletal muscle following a 12-h phase shift  as well as in livers of mice exposed to weekly 8-h phase shifts , which suggests that lipid metabolism may be affected by CCM across multiple species.
In order to control for circadian variation in gene expression, we extracted our RNA samples between ZT0–3 in both CCM and control flies. At this time, most flies in both groups are awake and active. However, we cannot rule out that some gene expression changes observed in CCM flies reflect the increased activity in the preceding few hours due to the prolonged morning anticipation in these flies, or to potential differences in feeding behavior associated with this prolonged activity. An assessment of gene expression in CCM flies across multiple circadian time points and behavioral states would help to differentiate between these possibilities.
Many more genes exhibited differential expression following 3 weeks of CCM compared to only 2 weeks. This difference may be due to a progressive dysregulation of expression such that changes are underway by week 2 that are too subtle to detect with whole fly RNA sequencing. The fact that there was an overall correlation between gene expression changes in weeks 2 and 3 supports this idea. On the other hand, it is possible that many of the gene expression changes observed in week 3 reflect general organismal decline, which appears to be accelerated in CCM-exposed flies. In fact, we saw substantial overlap between genes altered by CCM and those that exhibited expression changes due to aging. These included upregulation of stress-response genes and downregulation of genes involved in development and regulation of gene expression. In addition, we observed dysregulation of a number of genes involved in neuronal function, including those responsible for generation of circadian rhythms. This may indicate that molecular clock function is compromised, which could in turn exacerbate the negative health consequences of CCM, especially as the circadian clock protects against oxidative stress in aged flies [25, 38]. However, given that we found no effect of CCM on free-running behavioral rhythms, which are controlled by the central clock, our results are more consistent with reduced circadian function of peripheral clocks. Tissue-specific RNA sequencing in the face of circadian misalignment would help to determine whether this is the case. If peripheral clock function is indeed altered, then the reduced longevity produced by CCM could be a consequence both of the acute physiological effects of misalignment and a reduced ability to effectively mitigate those effects.
Our results confirm that CCM is associated with severe, adverse consequences, including reduced lifespan. We furthermore show that the negative health effects of CCM can occur independent of changes in sleep amount, which indicates that frequent circadian phase shifts are inherently unhealthy, likely due to behavioral and physiological processes occurring at sub-optimal times. Finally, our RNA sequencing analysis demonstrates that CCM produces large-scale changes in gene expression, and in particular identifies genes involved in lipid metabolism as being potentially important mediators of the negative health consequences of CCM. These studies identify candidate molecules and pathways that can be targeted in future studies to further delineate the etiology of CCM-induced organismal decline. Our findings have important implications for the large number of people who routinely follow behavioral schedules that run outside the dictates of their internal circadian clocks.
Prior to behavioral experiments, fly stocks were maintained in 6 oz. square-bottom polypropylene bottles (Fisher Scientific) with standard cornmeal-molasses food on a 12:12 light:dark (LD) cycle at 25 °C. For longevity experiments, male and female iso31 flies  were collected under light CO2 anesthesia within 1 day of eclosion. Individual flies were loaded into glass tubes containing a 5% sucrose and 2% agar food source for locomotor activity analysis with the Drosophila Activity Monitoring (DAM) system (Trikinetics, Waltham, MA). Humidity and temperature-controlled incubators (Percival, Perry IA) were used to expose flies to either a 24-h schedule (LD 12:12; control group), or a 28-h schedule (LD 14:14; CCM group), which results in daily 4-h phase delays. A 28-h schedule was chosen because previous studies indicated a more severe consequence of long versus short days on Drosophila lifespan  and because a daily 4-h phase shift results in CCM-exposed flies completing 6 full LD cycles in a week such that total weekly light and dark exposure was equivalent between CCM and control flies. Incubator temperature was held constant at 25 °C and humidity levels were kept between 70 and 80%. Temperature, humidity and light intensity were checked regularly to ensure uniformity between experimental groups. Flies were transferred to new tubes each week to supply fresh food. Occasionally, flies escaped during transfer to new tubes, and these flies were excluded from all analyses. Longevity analysis was conducted in two independent experiments, each of which included 64 male and female flies for each condition. Because consistent results were observed between these two experiments, we combined data for analysis.
The DAM system records each time a fly breaks an infrared beam that bisects the center of each glass tube. Longevity was determined by identifying the last minute that a beam break occurred. Occasionally we observed “ghost” readings, where single beam breaks were detected after extended periods of inactivity. Thus, we ignored isolated activity bins that occurred > 12 h after previous activity. Data were collected until all flies in the experiment were dead. Fly lifespans were analyzed with the Log-rank test (GraphPad Prism 7), with p < 0.05 considered significant.
Locomotor activity and sleep analysis
Analysis of locomotor activity rhythmicity was performed with ClockLab software (Actimetrics, Wilmette, IL). Locomotor activity rhythm period and strength (power) were determined with chi2 periodogram analysis. Rhythm power was defined as the amplitude of the periodogram line at the dominant period minus the chi2 significance line at a significance of p < 0.01. Weekly LD rhythmicity was concurrently determined for the same group of flies that were used for longevity experiments. In a separate group of flies, the period and power of free-running rhythms under conditions of constant darkness (DD) were calculated to assess the health of the endogenous circadian clock. These flies were exposed to control or CCM conditions for 1–3 weeks and then transferred to DD for an additional week for rhythmicity analysis. Because fly activity is affected by transfer between LD and DD incubators, the first day of DD was not included in rhythmicity analysis. For all rhythmicity analyses the bin length was 1 min. 2-Way ANOVA followed by Sidak’s posthoc test (GraphPad Prism 7) were used to compare rhythm strength between control and CCM flies at the different time points, with p < 0.05 considered significant.
The average (normalized) function in ClockLab was used to create mean normalized activity plots in Fig. 3i. Flies were kept in LD 14:14 (CCM) or 12:12 (control) conditions for 2 weeks prior to transfer to DD, and graphs show the final 3 days in LD followed by the first 3 days in DD.
Sleep was defined as 5 consecutive minutes of inactivity  as determined by a custom-developed Excel formula. 2-Way ANOVA followed by Sidak’s posthoc test (GraphPad Prism 7) were used to separately compare activity (in mean beam breaks/min) and sleep amount (in mean min/hr) between control and CCM flies for week 1, week 2, week 3, and full-life, with p < 0.05 considered significant. For full-life measurements, we removed the last three days of life from analysis because flies exhibited reduced activity during this time, making it difficult to separate sleep from an aging-induced decreases in locomotor activity.
For weekly rest:activity rhythm period and power as well as activity amount and sleep duration analysis, only flies that survived through the entire week were included and thus the number of flies analyzed varied between weeks. For female control flies, 121 flies survived through week 1, 118 survived through week 2, 92 through week 3, and 26 through week 4. For female CCM flies, 123 survived through week 1, 119 through week 2, 52 through week 3, and 13 through week 4. For control male flies, 125 survived through week 1, 124 through week 2, 84 through week 3, and 27 through week 4. For CCM male flies, 115 survived through week 1, 107 through week 2, 44 through week 3, and 6 through week 4. Thus, these are the ns for experiments depicted in Figs. 3 and 4. For Fig. 2, ns are slightly smaller since analysis was conducted starting from ZT18 on day 1 to ZT18 of day 8 for each week instead of starting at ZT0, and a handful of flies died between ZT0 and ZT18 at the end of the week.
Spearman’s rank test (GraphPad Prism 7) was used to calculate the correlation between sleep and longevity, with p < 0.05 considered significant. To further assess the effect of CCM and sleep time on longevity the Cox Proportional Hazard model was used. Hazard ratios for CCM and sleep time were first calculated in univariate models. For male flies, a multivariate Cox Proportional Hazard Model was also used with both treatment and week 1 sleep time in the model to determine if CCM reduced lifespan independent of its effects on sleep duration. In a separate approach, each control fly was matched with a CCM-exposed fly with a similar sleep time (within 0.5 min/hr) and a hazard ratio for CCM was calculated in this smaller data set.
Male and female iso31 flies were exposed to 2 or 3 weeks of control or CCM conditions in DAM monitors as described above. 10 male and 10 female flies were collected for each sample (3 CCM and 3 control samples at each time point, 12 total) and total RNA was extracted at ZT0–3 with TRIreagent (Molecular Research Center, Inc., Cincinnati, OH) according to the manufacturer’s recommendations. Following phase separation, samples were additionally subjected to DNAse treatment and purification with the Zymo Direct-zol RNA kit (Zymo Research, Irvine, CA). All samples passed quality control with RIN > 7.2 and a 260/280 ratio > 2.0. Library preparation with the NEB Next Ultra II kit and subsequent 150 base pair, paired-end RNA sequencing on the Hi-Seq platform were conducted by Novogene (Davis, CA).
Approximately 21–28 million raw reads were obtained per sample. Novogene filtered raw reads to remove reads containing adapters, reads containing N > 10% (N represents the bases that cannot be determined), and reads containing Phred quality score ≥ 5 over more than 50% of the read. Following filtering, > 94.78% of reads remained for each sample.
Differential gene expression analysis
Reads were mapped to Dmel Reference Genome (dm6) with RNA Star  in Galaxy (using the public Galaxy server at usegalaxy.org). Drosophila_melanogaster.BDGP6.87.gtf was used for the annotation file and filtered to remove all non-coding genes because we were interested in identifying differentially expressed genes that encode for proteins. The aligned BAM files were processed with MMquant  to quantify the number of reads that mapped to each gene. MMquant allows for quantification of duplicated genes by creating merged genes out of different genes for which the sequence is multiply mapped.
DEseq2  in Galaxy was used to identify genes exhibiting differential expression between CCM and control samples at the 2- and 3-week time points. We used a Benjamini-Hochberg FDR cutoff of 0.1 to determine differentially expressed genes, but used no cutoff for fold-change (FC), which was generally in the moderate range (log2FC ranged from − 1.26 to 0.88). DEseq2 uses shrinkage estimate to generate a conservative estimate of FC.
We also used DEseq2 to assess aging-related changes by comparing gene expression in control samples between week 2 and 3. We then identified genes among the 349 that were found to be up- or downregulated following 3 weeks of CCM that were oppositely regulated by aging, which included those with a log2FC in the opposite direction as that produced by CCM and that were significant with a FDR cutoff of 0.3.
We used the Pearson correlation coefficient (GraphPad Prism 7) to assess the correlation between the CCM-induced log2FC in gene expression between week 2 and week 3 for all genes as well as the subsets of genes found to be up- or downregulated following 3 weeks of CCM. p < 0.05 was considered significant.
Gene ontology (GO) term analysis
GO term analysis for Biological Function was performed with Princeton University’s GOTermFinder (http://go.princeton.edu/cgi-bin/GOTermFinder). Analysis was done separately for up- and down-regulated genes from each week. Because some genes have unique as well as overlapping sequences and therefore showed up more than once in the MMquant output, we manually curated the differentially expressed gene list to include only one incidence of each gene for gene ontology analysis. We provided the background gene list which consisted of all genes tested for differential expression. The Bonferroni-corrected p-value cutoff was set to 0.05 to determine over- or underrepresented GO terms. The resulting lists were then passed through REVIGO (http://revigo.irb.hr)  with an allowed similarity of 0.7 to remove redundant terms. Because there were > 350 GO terms associated with week 3 downregulated genes, we further filtered that list to remove highly generic terms (those with an average distance to root of < 4). Table 3 lists only the top 15 most significant of these filtered GO terms.
Quantitative PCR (qPCR) analysis
Male and female iso31 flies were exposed to 2 or 3 weeks of control or CCM conditions in DAM monitors and RNA extraction was performed as described above. Flies used for qPCR analysis were from a separate cohort from those used for RNA sequencing. cDNA synthesis was conducted with the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems) from 1 μg total RNA and qPCR was conducted with the PowerUp SYBR Green Master Mix (Applied Biosystems) using a StepOne Plus Real-Time PCR System (Applied Biosystems). We included 3 samples per group and 3 technical replicates per sample. Primers were as follows: CG18609 forward: 5′- CTCTTCATCGTGGGCATTTGC – 3′; CG18609 reverse: 5′- TCGAGGACCTTGTTCAGCAG -3′; CG3739 forward: 5′- AGGATAGTGACTTGCCCTGG – 3′; CG3739 reverse: 5′- GTGGATCCAAAGGGTTGGCT -3′. Quantification was done with the relative standard curve method, and gene expression was normalized to actin5c levels for each sample. RNA sequencing showed no expression changes in actin5c with CCM or aging. ANOVA with Sidak’s posthoc test was used to compare normalized gene expression levels. p < 0.05 was considered significant.
We thank Hallie Chen for help with experiments and members of the Cavanaugh Lab for helpful discussions.
This work was supported by a NARSAD Young Investigator Award from the Brain and Behavior Research Institute to DJC. The funding agency had no role in the design of the study, collection, analysis, or interpretation of data, or in writing the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. RNAseq data have been uploaded to Gene Expression Omnibus. Accession # GSE124065.
ACB contributed to conception and design of experiments, data acquisition, analysis and interpretation of data, and was involved in drafting and revising the manuscript. GDS contributed to data acquisition and analysis and interpretation of data. ACS contributed to analysis and interpretation of data. SDH contributed to data acquisition and analysis and interpretation of data. PP contributed to data acquisition and analysis and interpretation of data. HEW contributed to analysis and interpretation of data and was involved in revising the manuscript. CMD contributed to analysis and interpretation of data and was involved in revising the manuscript. DJC contributed to conception and design of experiments, data acquisition, analysis and interpretation of data, and was involved in drafting and revising the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Vaze KM, Sharma VK. On the adaptive significance of circadian clocks for their owners. Chronobiol Int. 2013;30:413–33.View ArticleGoogle Scholar
- Mohawk JA, Green CB, Takahashi JS. Central and peripheral circadian clocks in mammals. Annu Rev Neurosci. 2012;35:445–62. https://doi.org/10.1146/annurev-neuro-060909-153128.View ArticlePubMedPubMed CentralGoogle Scholar
- Dubowy C, Sehgal A. Circadian rhythms and sleep in Drosophila melanogaster. Genetics. 2017;205:1373–97. https://doi.org/10.1534/genetics.115.185157.View ArticlePubMedPubMed CentralGoogle Scholar
- Roenneberg T, Merrow M. The circadian clock and human health. Curr Biol. 2016;26:R432–43.View ArticleGoogle Scholar
- Wittmann M, Dinich J, Merrow M, Roenneberg T. Social jetlag: misalignment of biological and social time. Chronobiol Int. 2006;23:497–509.View ArticleGoogle Scholar
- Foster RG, Peirson SN, Wulff K, Winnebeck E, Vetter C, Roenneberg T. Sleep and circadian rhythm disruption in social jetlag and mental illness. Prog Mol Biol Transl Sci. 2013;119:325–46.View ArticleGoogle Scholar
- Golombek DA, Casiraghi LP, Agostino PV, Paladino N, Duhart JM, Plano SA, et al. The times they’re a-changing: effects of circadian desynchronization on physiology and disease. J Physiol Paris. 2013;107:310–22. https://doi.org/10.1016/j.jphysparis.2013.03.007.View ArticlePubMedGoogle Scholar
- Evans JA, Davidson AJ. Health consequences of circadian disruption in humans and animal models. Prog Mol Biol Transl Sci. 2013;119:283–323.View ArticleGoogle Scholar
- Kettner NM, Voicu H, Finegold MJ, Coarfa C, Sreekumar A, Putluri N, et al. Circadian homeostasis of liver metabolism suppresses Hepatocarcinogenesis. Cancer Cell. 2016;30:909–24.View ArticleGoogle Scholar
- Pittendrigh CS, Minis DH. Circadian systems: longevity as a function of circadian resonance in Drosophila melanogaster. Proc Natl Acad Sci U S A. 1972;69:1537–9.View ArticleGoogle Scholar
- von Saint Paul U, Aschoff J. Longevity among blowfliesPhormia terraenovae R.D. kept in non-24-hour light-dark cycles. J Comp Physiol. 1978;127:191–5. https://doi.org/10.1007/BF01350109.View ArticleGoogle Scholar
- Vaccaro A, Birman S, Klarsfeld A. Chronic jet lag impairs startle-induced locomotion in Drosophila. Exp Gerontol. 2016;85:24–7.View ArticleGoogle Scholar
- Castanon-Cervantes O, Wu M, Ehlen JC, Paul K, Gamble KL, Johnson RL, et al. Dysregulation of inflammatory responses by chronic circadian disruption. J Immunol. 2010;185:5796–805.View ArticleGoogle Scholar
- Davidson AJ, Sellix MT, Daniel J, Yamazaki S, Menaker M, Block GD. Chronic jet-lag increases mortality in aged mice. Current biology. 2006;16:R914–6.View ArticleGoogle Scholar
- Penev PD, Kolker DE, Zee PC, Turek FW. Chronic circadian desynchronization decreases the survival of animals with cardiomyopathic heart disease. Am J Phys. 1998;275(6 Pt 2):H2334–7.Google Scholar
- Park N, Cheon S, Son GH, Cho S, Kim K. Chronic circadian disturbance by a shortened light-dark cycle increases mortality. Neurobiol Aging. 2012;33:1122.e11–22.View ArticleGoogle Scholar
- Filipski E, Delaunay F, King VM, Wu M-W, Claustrat B, Grechez-Cassiau A, et al. Effects of chronic jet lag on tumor progression in mice. Cancer Res. 2004;64:7879–85.View ArticleGoogle Scholar
- Aschoff J, Pohl H. Phase relations between a circadian rhythm and its zeitgeber within the range of entrainment. Naturwissenschaften. 1978;65:80–4.View ArticleGoogle Scholar
- Koh K, Evans JM, Hendricks JC, Sehgal A. A Drosophila model for age-associated changes in sleep:wake cycles. Proc Natl Acad Sci U S A. 2006;103:13843–7.View ArticleGoogle Scholar
- Shaw PJ, Cirelli C, Greenspan RJ, Tononi G. Correlates of sleep and waking in Drosophila melanogaster. Science. 2000;287:1834–7.View ArticleGoogle Scholar
- Kayser MS, Yue Z, Sehgal A. A critical period of sleep for development of courtship circuitry and behavior in Drosophila. Science. 2014;344:269–74.View ArticleGoogle Scholar
- Huber R, Hill SL, Holladay C, Biesiadecki M, Tononi G, Cirelli C. Sleep homeostasis in Drosophila melanogaster. Sleep. 2004;27:628–39.View ArticleGoogle Scholar
- Bushey D, Hughes KA, Tononi G, Cirelli C. Sleep, aging, and lifespan in Drosophila. BMC Neurosci. 2010;11:56.View ArticleGoogle Scholar
- Koudounas S, Green EW, Clancy D. Reliability and variability of sleep and activity as biomarkers of ageing in Drosophila. Biogerontology. 2012;13:489–99.View ArticleGoogle Scholar
- Kuintzle RC, Chow ES, Westby TN, Gvakharia BO, Giebultowicz JM, Hendrix DA. Circadian deep sequencing reveals stress-response genes that adopt robust rhythmic expression during aging. Nat Commun. 2017;8:14529.View ArticleGoogle Scholar
- Weiss R, Bartok O, Mezan S, Malka Y, Kadener S. Synergistic interactions between the molecular and neuronal circadian networks drive robust behavioral circadian rhythms in Drosophila melanogaster. PLoS Genet. 2014;10:e1004252.View ArticleGoogle Scholar
- Klarsfeld A, Rouyer F. Effects of circadian mutations and LD periodicity on the life span of Drosophila melanogaster. J Biol Rhythm. 1998;13:471–8.View ArticleGoogle Scholar
- Martino TA, Oudit GY, Herzenberg AM, Tata N, Koletar MM, Kabir GM, et al. Circadian rhythm disorganization produces profound cardiovascular and renal disease in hamsters. Am J Physiol Regul Integr Comp Physiol. 2008;294:R1675–83.View ArticleGoogle Scholar
- Takasu NN, Nakamura TJ, Tokuda IT, Todo T, Block GD, Nakamura W. Recovery from age-related infertility under environmental light-dark cycles adjusted to the intrinsic circadian period. Cell Rep. 2015;12:1407–13.View ArticleGoogle Scholar
- DeCoursey PJ, Walker JK, Smith SA. A circadian pacemaker in free-living chipmunks: essential for survival? J Comp Physiol A. 2000;186:169–80.View ArticleGoogle Scholar
- Yamazaki S, Numano R, Abe M, Hida A, Takahashi R, Ueda M, et al. Resetting central and peripheral circadian oscillators in transgenic rats. Science. 2000;288:682–5.View ArticleGoogle Scholar
- Archer SN, Laing EE, Moller-Levet CS, van der Veen DR, Bucca G, Lazar AS, et al. Mistimed sleep disrupts circadian regulation of the human transcriptome. Proc Natl Acad Sci U S A. 2014;111:E682–91.View ArticleGoogle Scholar
- Kervezee L, Cuesta M, Cermakian N, Boivin DB. Simulated night shift work induces circadian misalignment of the human peripheral blood mononuclear cell transcriptome. Proc Natl Acad Sci U S A. 2018;115:5540–5.View ArticleGoogle Scholar
- Landis G, Shen J, Tower J. Gene expression changes in response to aging compared to heat stress, oxidative stress and ionizing radiation in Drosophila melanogaster. Aging (Albany NY). 2012;4:768–89.View ArticleGoogle Scholar
- Jung A, Hollmann M, Schafer MA. The fatty acid elongase NOA is necessary for viability and has a somatic role in Drosophila sperm development. J Cell Sci. 2007;120(Pt 16):2924–34.View ArticleGoogle Scholar
- Bustos V, Partridge L. Good Ol’ fat: links between lipid signaling and longevity. Trends Biochem Sci. 2017;42:812–23.View ArticleGoogle Scholar
- Wefers J, van Moorsel D, Hansen J, Connell NJ, Havekes B, Hoeks J, et al. Circadian misalignment induces fatty acid metabolism gene profiles and compromises insulin sensitivity in human skeletal muscle. Proc Natl Acad Sci U S A. 2018;115:7789–94.View ArticleGoogle Scholar
- Giebultowicz JM. Circadian regulation of metabolism and healthspan in Drosophila. Free Radic Biol Med. 2018;119:62–8.View ArticleGoogle Scholar
- Ryder E, Blows F, Ashburner M, Bautista-Llacer R, Coulson D, Drummond J, et al. The DrosDel collection: a set of P-element insertions for generating custom chromosomal aberrations in Drosophila melanogaster. Genetics. 2004;167:797–813.View ArticleGoogle Scholar
- Hendricks JC, Finn SM, Panckeri KA, Chavkin J, Williams JA, Sehgal A, et al. Rest in Drosophila is a sleep-like state. Neuron. 2000;25:129–38.View ArticleGoogle Scholar
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.View ArticleGoogle Scholar
- Zytnicki M. Mmquant: how to count multi-mapping reads? BMC Bioinformatics. 2017;18:411.View ArticleGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.View ArticleGoogle Scholar
- Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.View ArticleGoogle Scholar