ORIGINAL RESEARCH article

Front. Endocrinol., 26 April 2018

Sec. Neuroendocrine Science

Volume 9 - 2018 | https://doi.org/10.3389/fendo.2018.00199

Time-of-Day Effects on Metabolic and Clock-Related Adjustments to Cold

    FS

    Frederico Sander Mansur Machado 1,2,3

    ZZ

    Zhi Zhang 1,2,4

    YS

    Yan Su 2,4

    PD

    Paul de Goede 2,4

    RJ

    Remi Jansen 2,4

    EF

    Ewout Foppen 1,4

    CC

    Cândido Celso Coimbra 3

    AK

    Andries Kalsbeek 1,2,4*

  • 1. Department of Endocrinology and Metabolism, Academic Medical Center, University of Amsterdam, Amsterdam, Netherlands

  • 2. Hypothalamic Integration Mechanisms, Netherlands Institute for Neuroscience, Amsterdam, Netherlands

  • 3. Department of Physiology and Biophysics, Biological Sciences Institute, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil

  • 4. Laboratory of Endocrinology, Department of Clinical Chemistry, Academic Medical Centre, University of Amsterdam, Amsterdam, Netherlands

Article metrics

View details

17

Citations

5,4k

Views

1,8k

Downloads

Abstract

Background:

Daily cyclic changes in environmental conditions are key signals for anticipatory and adaptive adjustments of most living species, including mammals. Lower ambient temperature stimulates the thermogenic activity of brown adipose tissue (BAT) and skeletal muscle. Given that the molecular components of the endogenous biological clock interact with thermal and metabolic mechanisms directly involved in the defense of body temperature, the present study evaluated the differential homeostatic responses to a cold stimulus at distinct time-windows of the light/dark-cycle.

Methods:

Male Wistar rats were subjected to a single episode of 3 h cold ambient temperature (4°C) at one of 6 time-points starting at Zeitgeber Times 3, 7, 11, 15, 19, and 23. Metabolic rate, core body temperature, locomotor activity (LA), feeding, and drinking behaviors were recorded during control and cold conditions at each time-point. Immediately after the stimulus, rats were euthanized and both the soleus and BAT were collected for real-time PCR.

Results:

During the light phase (i.e., inactive phase), cold exposure resulted in a slight hyperthermia (p < 0.001). Light phase cold exposure also increased metabolic rate and LA (p < 0.001). In addition, the prevalence of fat oxidative metabolism was attenuated during the inactive phase (p < 0.001). These metabolic changes were accompanied by time-of-day and tissue-specific changes in core clock gene expression, such as DBP (p < 0.0001) and REV-ERBα (p < 0.01) in the BAT and CLOCK (p < 0.05), PER2 (p < 0.05), CRY1 (p < 0.05), CRY2 (p < 0.01), and REV-ERBα (p < 0.05) in the soleus skeletal muscle. Moreover, genes involved in substrate oxidation and thermogenesis were affected in a time-of-day and tissue-specific manner by cold exposure.

Conclusion:

The time-of-day modulation of substrate mobilization and oxidation during cold exposure provides a clear example of the circadian modulation of physiological and metabolic responses. Interestingly, after cold exposure, time-of-day mostly affected circadian clock gene expression in the soleus muscle, despite comparable changes in LA over the light–dark-cycle. The current findings add further evidence for tissue-specific actions of the internal clock in different peripheral organs such as skeletal muscle and BAT.

Introduction

Daily cyclic changes in environmental conditions are key signals for the adaptive and anticipatory activity of most living species, including mammals. The mammalian thermoregulatory system is fairly adapted to periodic changes in ambient temperatures that may reach high amplitudes depending on the geographic location (13). Intriguingly, in mammals, a role for daily body temperature cycles in the internal synchronization has been demonstrated, in vivo, ex vivo, as well as in vitro (1, 4, 5), indicating that the circadian timing system and the thermoregulatory system reciprocally influence each other.

In general, during the dark phase, ambient temperature decreases to its lowest daily levels. Lower ambient temperatures induce metabolic changes aimed to defend internal body temperature. Considering this, it is well known that a cold environment induces autonomic, cardiovascular, metabolic, and behavioral adjustments that depend on the synchronized activation of multiple independent pathways resulting in thermal adaptation/acclimation (6). These adjustments include physiological changes such as vasoconstriction (heat retention or storage) and thermogenesis from both the activation of brown adipose tissue (BAT) (non-shivering) and skeletal muscle involuntary contraction (shivering) (6). Synchronized activation of the autonomic innervation to BAT and white adipose tissue (WAT), liver, adrenal, and skeletal muscle is necessary to produce the necessary amount of energy and heat to keep body temperature within safe levels during cold exposure (7).

Contractile skeletal muscle activity acts as an important heat source during environmental cold exposure in rodents (6, 8) and humans (9, 10). Increased ADP/ATP ratio, 5′adenosine monophosphate-activated protein kinase (AMPK) activity, peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC1-α) activity, and consequently, free-fatty acid (FFA) uptake and oxidation, increases intracellular substrate availability for heat production (11). Notably, in BAT, the requirement of a functional molecular clock has been demonstrated, as a deficiency in the expression of the BMAL-1 and PER2 clock gene leads to activation of compensatory heat production mechanisms (12, 13). In humans, it has recently been demonstrated that BAT glucose uptake might be associated with the heat production rhythm (14). The latter data indicate that the central body clock has an important role in the tuning of the cold-evoked response, probably by setting the basal metabolic rate to a new reference level or balancing the autonomic tonus.

Metabolic rate increases in response to cold through shivering and non-shivering thermogenesis. Both processes are centrally regulated (6) and result in increased lipid and carbohydrate oxidation in involved tissues (mainly skeletal muscle and BAT) (8, 15). Therefore, lipid mobilization (from WAT lipolysis and liver de novo lipogenesis) increases to provide the main substrate used for the thermogenic activity. Interestingly, the molecular clock influences substrate oxidation (16, 17). Therefore, modulation of the central thermoregulatory pathways by the biological clock in the suprachiasmatic nuclei (SCN) might result in changes in cold exposure adjustments depending on the time-of-day. In line with this hypothesis, the master clock projects to the major brain areas involved in metabolic/thermal balance (18, 19).

Core components of the molecular circadian clock are expressed throughout the body, both centrally and peripherally and interact with intracellular pathways directly related to metabolism and heat production (1214, 16, 17, 20). In view of the above, the present study aimed to investigate how time-of-day modulates the peripheral adjustments induced by cold exposure at a physiological and a molecular level.

Materials and Methods

Animals

All experiments were performed in adult male Wistar rats (Charles River Breeding Laboratories, Sulzfeld, Germany). After arrival at the animal facility, animals were housed in individual cages (25 cm × 25 cm × 35 cm), with a 12/12-h light–dark (L/D) schedule [lights on at 0700 h, defined as Zeitgeber Time 0 (ZT0)]. Animals were allowed to adapt to the new environment for 1.5 weeks before the first experiments. All rats were kept under constant temperature (22 ± 2°C) and humidity (50 ± 5%) conditions. Food and water were available ad libitum. The animal care committee of the Royal Netherlands Academy of Arts and Sciences (DEC/KNAW) approved all experiments.

Experimental Procedures

To verify where would be a good site to monitor internal body temperature, a pilot study group was completed before the described experiments. In this preliminary group, each animal was implanted with two data loggers: one in the dorsal subcutaneous area, caudal to the BAT and the other one inside the peritoneal cavity, stitched to the abdominal wall. Since these loggers are not radio-telemetry based, interference between them was not observed, as it would be expected for other available models. After recovery from surgery, animals were exposed to episodes of lowered ambient temperatures during light and dark phases of the L/D cycle. With this initial study, we were able to identify that the intraperitoneal loggers produced the most stable and reliable results.

Experimental animals were anesthetized with isoflurane, which guaranteed a rapid recovery from the small surgery necessary to insert the temperature sensors. A ventral incision at the linea alba was made to introduce a data-logger probe (DST nano-T, StarOddi, Iceland) into the peritoneal cavity. Each logger was sutured to the inner musculature before the incision was closed. This procedure allowed continuous monitoring of core body temperature (Tcore) with a decreased risk of internal displacement of the sensor, which could cause misleading readings due to its position. These probes recorded internal Tcore with 5-min intervals.

For the main study (presented in Figure 1), on the fourth day after surgery, the basal 24-h locomotor activity (LA), energy expenditure (EE), and food and water intake were continuously recorded for each animal with an indirect calorimetry system (PhenoMaster/LabMaster, TSE Systems, Bad Homburg, Germany). LA was assessed as beam-breaks recorded during 15 min intervals. Oxygen consumption (VO2) and carbon dioxide production (VCO2) were assessed every 15 min for 100 s. Respiratory exchange ratio (RER) was calculated according to the formula: VCO2/VO2. EE was calculated with the formula: (CVO2 × VO2 + CVCO2 × VCO2)/1,000, considering CVO2 and CVCO2 preset reference values given by the manufacturer (CVO2 = 3.941 and CVCO2 = 1.106). Carbohydrate (CHO) and lipid oxidation were calculated with the adapted formulae (21, 22) previously used (23): (4.585 × VCO2 − 3.226 × VO2)/60,000 for CHO oxidation and (1.695 × VO2 − 1.701 × VCO2)/60,000 for lipid oxidation. After the basal recordings, the ambient temperature was lowered to 0–5°C at different times of the day: ZT23–2, ZT3–6, ZT7–10, ZT11–14, ZT15–18, and ZT19–22. The average rate of ambient cooling was −0.8°C/10 min during the period of 3 h of cooling activity. Peak rates of −2.3°C/10 min were observed during the first 45 min of the experimental protocol, while rates lower than −0.5°C/10 min were observed and maintained after 90 min of cooling. By the end of each episode of lowered ambient temperature, animals were rapidly anesthetized with 80% CO2 and immediately sacrificed by decapitation. For room temperature controls, animals were kept in the same experimental room, but outside the temperature controlled TSE chamber without any previous manipulation and were sacrificed at the same time points. After decapitation, the brain was removed, snap frozen on dry ice, and stored at −80°C. Soleus skeletal muscle and interscapular BAT were collected, frozen in liquid nitrogen, and then stored at −80°C.

Figure 1

Figure 1

Graphical representation of experimental paradigm.

Tissue Processing

RNA Extraction and cDNA Synthesis

Total RNA was extracted and purified with the TRIzol reagent protocol (Macherey-Nagel, Oensingen, Switzerland). The quality of RNA was examined by Agilent 2100 Bioanalyzer equipped with Nano chips (Agilent Technologies, Palo Alto, CA, USA) and concentrations were determined by Nanodrop spectrophotometer (ThermoScientific Technologies, Wilmington, DE, USA). A fixed amount of total RNA was reverse-transcribed with SensiFAST cDNA Synthesis Kit (Bioline, Taunton, MA, USA). For the control of genomic DNA contamination, we employed a minus reverse transcriptase sample (−RT).

Real-Time PCR (RT-PCR)

The expression of clock, metabolic, and thermogenesis-related genes was evaluated by RT-PCR (LightCycler ® 480, Roche) with the following reaction system: 2 µl of cDNA was incubated with 50 ng of both reverse and forward primer from gene of interest (see Table 1 for primer sequences) and SensiFAST no-ROX Mix (Bioline, Taunton, MA, USA) for a final volume of 10 µl. The relative amount of each gene was normalized against the geometric mean of three housekeeping genes: hypoxanthine-guanine phosphoribosyl transferase (HPRT), ribosomal protein S18 (S18), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) for the BAT; GAPDH, S18, and cyclophilin for the soleus muscle. The reference genes mentioned (HPRT, GAPDH, S18, and cyclophilin) were selected based on their constant expression under the various experimental conditions (i.e., time of the day and/or ambient temperature). The relative expression level of clock genes and other genes of interest in each sample was obtained by dividing the absolute amount of the target gene by the average of the reference genes values. Each gene/tissue RT-PCR was performed in a single plate.

Table 1

Genes Reference no. Reverse Forward
Clock genes
PER1 NM_001034125.1 TGGCCAGGATCTTGAACACTGCTA ATGCAGAAACAACAGCCACGGTTC
PER2 NM_031678.1 CAACGCCAAGGAGCTCAAGT CACCCTGAAAAGAAAGTGCGA
CRY1 NM_198750.2 TCATCATGGTCGTCGGACAGA AAGTCATCGTGCGCATTTCA
CRY2 NM_133405.2 TGTACAAGTCCCACAGGCGGTA TGGATAAGCACTTGGAACGGAA
ARNLT/BMAL-1 NM_024362.2 TGCAGTGTCCGAGGAAGATAGC CCGATGACGAACTGAAACACCT
CLOCK NM_021856.2 TTGCAGCTTGAGACATCGCT CGATCACAGCCCAACTCCTT
DBP NM_012543.3 TGCCTTCTTCATGATTGGCTG CCTTTGAACCTGATCCGGCT
REV-ERBα NM_001113422.1 CATGGGCATAGGTGAAGATTTCT ACAGCTGACACCACCCAGATC

Metabolic genes
CREB1 NM_031017.1 ACTCTGCTGGTTGTCTGCTC GCAGTGACTGAGGAGCTTGT
PGC1-α NM_031347.1 GGTCATTTGGTGACTCTGG TGCCATTGTTAAGACCGAG
PGC1-β NM_176075.2 AGGAGGGCTCATTGCGTTTT AAAAGGCCATCGGTGAAGGT
PPAR-α NM_013196.1 GGCCTTGACCTTGTTCATGT TCACACAATGCAATCCGTTT
PPAR-γ NM_013124.3 GGGGGTGATATGTTTGAACTTG CAGGAAAGACAACAGACAAATCA
HSP90 NM_001004082.3 ACCGAATCTTGTCCAGGGCATCA CGGGCCCACCCTGCTCTGTA
UCP1 NM_012682.2 GCTTTGTGCTTGCATTCTGA AATCAGCTTTGCTTCCCTCA
UCP2 NM_019354.3 GGGCACCTGTGGTGCTAC GACTCTGTAAAGCAGTTCTACACCAA
UCP3 NM_013167.2 ATAGTCAGGATGGTACCGAGCA GCACTGCAGCCTGTTTTGCTGA
CIRBP NM_031147.2 TAACCACCACCCCTCCAGAA GCGTTAGGAAGCTTGGGTGT
CPT1-α NM_031559.2 AAAGACTGGCGCTGCTCA ACAATGGGACATTCCAGGAG
CPT1-β NM_013200.1 TGCTTGACGGATGTGGTTCC GTGCTGGAGGTGGCTTTGGT-
AMPK NM_019142.2 TAGAGAATGACCCCGCTGCT TGTCACAGGCATATGGTGGTC
NAMPT NM_177928.3 TCGACACTATCAGGTGTCTCAG ACAGATACTGTGGCGGGAATTGCT
FAT/CD36 NM_001109218.1 CCTTGGCTAAATAACGAACTCTG ACAGTTTTGGATCTTTGACGTG
GLUT4 NM_012751.1 CAGCGAGGCAAGGCTAGA GGGCTGTGAGTGAGTGCTTTC
HSL NM_012859.1 CCACCCGTAAAGAGGGAACT TCACGCTACATAAAGGCTGCT
LPL NM_012598.2 AGCAATTCCCCGATGTCCA CAAAACAACCAGGCCTTCGA
ADR-β2 NM_012739.3 CGACCGCTATGAGCGTGTAG CGCTTCACGTTCGTGCTGGC
ADR-β3 NM_013108.2 CCTTGCTAGATCTCCATGG CTTCCCAGCTAGCCCTGTT
GR NM_012576 GGAGCAAAGCAGAGCAGGTTT ACCTGGATGACCAAATGACCC
FOXO-1 NM_001191846.2 GTAGGGACAGATTGTGGCGAA ACGAGTGGATGGTGAAGAGTG
ACC1 NM_022193.1 CAGGCTACCATGCCAATCTC GATGATCAAGGCCAGCTTGT
ACC2 NM_053922.1 GCTTCCGCTCCAGGGTAGAGT GCACGAGATTGCTTTCCTAG
mTOR NM_019906.1 CCCGAGGAATCATACAGGTG AGCAGCATGGGGTTTAGGT
CamK2a NM_012920.1 AAGGCTGTCATTCCAGGGTC TGGCGTGAAGGAATCCTCTG

Housekeeping genes
S18 NM_213557.1 TGGCCAGAACCTGGCTATACTTCC CTCTTCCACAGGAGGCCTACACG
HPRT1 NM_012583.2 AACAAAGTCTGGCCTGTATCCAA GCAGTACAGCCCCAAAATGG
GAPDH NM_017008.4 TCCACCACCCTGTTGCTGTA TGAACGGGAAGCTCACTGG
Cyclophilin NM_017101.1 GAAGGAATGGTTTGATGGGT ATGTGGTCTTTGGGAAGGTG

Information about gene primers.

Statistical Analysis

The effects of cold exposure and duration in Tcore, LA, VO2, RER, EE, CHO, and lipid oxidation, as well as food and water intake, were assessed with a repeated measures ANOVA two-way followed by an appropriate post hoc test. The net area under the curve (AUC) for each condition and parameter was calculated, i.e., the absolute change of the physiological variable during the analyzed experimental period compared to t = 0. “Net” AUC in this case means that the area of the negative peaks (i.e., decrease compared to t = 0) was subtracted from that of the positive peaks (i.e., increases compared to t = 0). To analyze the combined effect of cold exposure and time-of-day an ANOVA two-way was applied to the AUC data. The time-of-day effect was assessed with an independent ANOVA one-way for control and cold-exposed situations followed by post hoc analysis when appropriate. The effects of cold exposure and time-of-day on mRNA expression were assessed with ANOVA two-way followed by an appropriate post hoc test. Data are expressed as mean ± SEM. Differences were considered statistically significant at p < 0.05.

Results

Daily Oscillations of Body Temperature, LA, and Energy Metabolism

During the basal data collection, 24 h before the actual experiment, ambient temperature inside the calorimetric chamber was kept close to the typical animal facility temperature of 22.18 ± 0.21°C (Figure 2A), although a time-of-day effect was apparent (Table 2). Animals exhibited a clear day/night oscillation in all the studied variables, despite the high variability observed for both food and water intake (Table 2; Figure 2). Core body temperature (Figure 2B) and spontaneous LA (Figure 2C) peaked during the active phase. The increased activity was associated with increased food and water intake (Figures 2D,E) and resulted in augmented oxygen consumption (Figure 2F) and heat production (Figure 2G). During the active phase, substrate utilization was shifted toward carbohydrate oxidation, while during the rest phase, lipid metabolism was predominant as observed through RER (Figure 2H). The baseline levels of thermal, behavioral, and metabolic parameters at the onset of our cooling paradigm were similar to those observed during the previous day and both exhibited time of the day related changes (Table S1 in Supplementary Material).

Figure 2

Figure 2

Daily rhythms of thermal, behavioral, and metabolic parameters under regular animal facility conditions. Light phase begins at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) begins at ZT12. Regular daily oscillatory levels of the ambient temperature (A), body temperature (B), locomotor activity (C), food intake (D), water intake (E), oxygen consumption (F), heat production (G), and respiratory quotient (H) were monitored approximately 24 h before the cold episodes. Results are presented as mean ± SEM (n = 34–50).

Table 2

Mean24-h ± SD Meanlight phase ± SD Meandark phase ± SD n Time-of-day effect p Value
Ambient temperature (°C) 22.18 ± 0.21 21.99 ± 0.07 22.39 ± 0.06 50 F(96, 4.753) = 7.059 <0.0001
Core body temperature (°C) 37.20 ± 0.25 37.00 ± 0.14 37.40 ± 0.14 34 F(96, 3.201) = 2.695 <0.0001
Locomotor activity (Cnts) 337.1 ± 194.8 175.5 ± 81.1 502.1 ± 125.6 50 F(96, 4.753) = 12.54 <0.0001
Oxygen consumption (ml/h/kg) 1,500 ± 196.1 1,327 ± 86.65 1,677 ± 86.64 50 F(96, 4.753) = 14.44 <0.0001
Respiratory exchange ratio (VO2/CO2) 1.02 ± 0.02 1.03 ± 0.01 1.01 ± 0.01 50 F(96, 4.753) = 8.036 <0.0001
Heat Expenditure (kcal/h/kg) 7.56 ± 1.00 6.67 ± 0.44 8.47 ± 0.42 50 F(96, 4.753) = 15.51 <0.0001
Food intake (g/15 min) 0.24 ± 0.17 0.12 ± 0.10 0.36 ± 0.12 50 F(96, 4.749) = 4.241 <0.0001
Water intake (ml/15 min) 0.21 = 0.19 0.05 ± 0.05 0.37 ± 0.15 49 F(96, 4.749) = 5.263 <0.0001

Mean basal 24-h and light/dark levels of thermal, behavioral, and metabolic parameters under control conditions (24 h before environmental cooling).

Effect of Time-of-Day on Body Temperature and EE During Acute Environmental Cooling

Average ambient temperature after 2 h of cooling was 4.01 ± 0.07°C, regardless of time of the day. Body temperature was affected by cold exposure in a time-of-day-dependent fashion (Figures 3A,B). When experimental cooling started during the light phase (ZT3, 7, and 11), a slight and transient hyperthermia was observed (p < 0.001). Interestingly, during the ZT15–18 protocol, core body temperature stayed elevated during the final 120 min of cold exposure (p < 0.001), despite the normal decrease in this parameter observed in the control animals. Cold exposure during the ZT19–22 and ZT23–2 protocols did not affect core body temperature as compared to baseline values. The time-of-day-dependent body temperature responses were confirmed by the AUC analysis (Figure 3B; Table 3; Table S2 in Supplementary Material), showing a significant effect of both ZT and cold exposure.

Figure 3

Figure 3

Effect of time-of-day on body temperature (A,B), VO2(C,D), and energy expenditure (E,F) changes induced by environmental cooling. Results are presented as averages of 15 min bins during 2 h of environmental cooling (A,C,E) and as the area under the curve (AUC) (B,D,F) calculated as the net change from the basal levels before cold exposure. Light phase began at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) began on ZT12. Results are presented as mean ± SEM. * indicates differences between the first and the other time-points within the control or the cold-exposed group, p < 0.05. # indicates differences between control and cold conditions, p < 0.05. n = 5–11/group/ZT.

Table 3

Effect of ZT
Effect of cold
Interaction
F(DFn, DFd) p Value F(DFn, DFd) p Value F(DFn, DFd) p Value
Body temperature F(5, 56) = 3.734 0.0055 F(1, 56) = 4.487 0.0386 F(5, 56) = 2.183 0.0689
Locomotor activity F(5, 74) = 0.372 0.8665 F(1, 74) = 42.21 <0.0001 F(5, 74) = 0.839 0.5265
VO2 F(5, 88) = 3.743 0.0040 F(1, 88) = 201.2 <0.0001 F(5, 88) = 0.3105 0.9055
Respiratory exchange ratio F(5, 88) = 3.404 0.0074 F(1, 88) = 107.6 <0.0001 F(5, 88) = 5.558 0.0002
CHO oxidation F(5, 88) = 3.314 0.0087 F(1, 88) = 1.619 0.2067 F(5, 88) = 4.435 0.0012
Lipid oxidation F(5, 88) = 1.051 0.3931 F(1, 88) = 72.19 <0.0001 F(5, 88) = 2.017 0.0839
Energy expenditure F(5, 88) = 3.889 0.0031 F(1, 88) = 230.2 <0.0001 F(5, 88) = 2.405 0.0430
Food intake F(5, 86) = 0.074 0.9960 F(1, 86) = 2.434 0.1224 F(5, 86) = 1.018 0.4123
Water intake F(5, 88) = 1.065 0.3851 F(1, 88) = 0.037 0.8470 F(5, 88) = 1.433 0.2203

ANOVA two-way table (ZT × ambient temperature) demonstrating the effect of cold exposure and time of the day on net area under the curve (in arbitrary units for each parameter) for thermal, behavioral, and metabolic parameters.

Body temperature is the result of heat production and heat dissipation. Therefore, understanding heat production dynamics during cold exposure at different phases of the light/dark daily cycle would be of value for the present analysis. A good proxy for heat production is oxygen consumption. Oxygen consumption showed a strong increase during cold exposure (p < 0.0001; Figures 3C,D), however, despite the observed effect of ZT, there was no interaction between the two factors (Table 3). As soon as the ambient temperature dropped, VO2 rose until reaching a steady state, usually about 60 min after the start of the protocol (p < 0.0001). Together with the rise in metabolic rate (VO2), EE significantly increased with cold exposure (Figures 3E,F), independent of the time-of-day (Figure 3F; Table 3).

Effect of Time-of-Day on LA and Food/Water Intake During Acute Environmental Cooling

In view of the observed effects of time-of-day on the changes in body temperature and heat production in response to a cold environment, it was necessary to investigate whether cold-induced changes in LA or food and/or water intake could play a role. Despite their intrinsic relationship, the changes in Tcore were not strictly accompanied by changes in LA in a time-of-day-dependent manner (Figures 4A,B; Table 3). Cold exposure induced an increase in LA (Figures 4A,B) that was independent of Zeitgeber time (Table 3; Table S2 in Supplementary Material). Regarding food and water intake, regardless of time-of-day cold exposure did not affect these parameters (Figures 4C–F; Table 3).

Figure 4

Figure 4

Effect of time-of-day on locomotor activity (A,B), food intake (C,D), and water intake (E,F) changes induced by environmental cooling. Results are presented as averages of 15 min bins during 2 h of environmental cooling (A,C,E) and area under the curve (AUC) (B,D,F) calculated as the net change from the basal levels before cold exposure. Light phase began at ZT and dark phase (shaded areas) began on ZT12. Results are presented as mean ± SEM. * indicates differences between the first and the other time-points within the control or the cold-exposed group, p < 0.05. # indicates differences between control and cold conditions, p < 0.05. n = 7–11/group/ZT.

Effect of Time-of-Day on Respiratory Quotient and Substrate Oxidation During Acute Environmental Cooling

The respiratory quotient, an index of substrate oxidation prevalence, decreased during every cold exposure, an effect that was modulated by time-of-day (Figures 5A,B; Table 3). In line with a decrease in RER, cold exposure evoked a consistent increase in lipid oxidation (Figures 5E,F). Interestingly, the higher lipid oxidation’s AUC was apparent only at ZT7–9, 11–13, 15–17, and 19–21, indicating that this rise was accentuated during the dark phase (Figure 5F). Although carbohydrate (CHO) oxidation was not significantly affected by cold exposure, it showed a time-of-day association (Figures 5C,D; Table 3). Increased CHO oxidation was observed only during the ZT3–6 protocol and the first hour of the ZT11–14 protocol (p < 0.001). On the other hand, CHO oxidation decreased when animals were exposed to cold during the dark phase, especially at ZT15–18 and 19–22 (p < 0.05), when also the largest decreases in RER were observed.

Figure 5

Figure 5

Effect of time-of-day on respiratory quotient (A,B), carbohydrate (C,D), and lipid oxidation (E,F) changes induced by environmental cooling. Results are presented as averages of 15 min bins during 2 h of environmental cooling (A,C,E) and area under the curve (AUC) (B,D,F) calculated as the net change from the basal levels before cold exposure. Light phase began at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) began on ZT12. Results are presented as mean ± SEM * indicates differences between the first and the other time-points within the control or the cold-exposed group, p < 0.05. # indicates differences between control and cold conditions, p < 0.05. n = 7–11/group/ZT.

Effect of Time-of-Day on Clock Genes mRNA Expression in Skeletal Muscle and BAT After Acute Environmental Cooling

The mRNA expression of eight clock genes was investigated in both soleus muscle and BAT (Tables 4 and 5; Figure 5). Time-of-day affected the expression of five of eight and six of eight clock genes in soleus muscle and BAT, respectively. Expression of CRY2 and CLOCK genes was not significantly affected by time-of-day in either tissue, although CLOCK gene expression showed a tendency in the BAT (p = 0.054). Cold exposure caused a significant change in most (seven of eight) of the clock genes studied in BAT and half of the clock genes studied in the soleus muscle. In the BAT, time-of-day interacted with cold exposure only for DBP and REV-ERBα expression. In soleus muscle, time-of-day interacted with cold exposure for clock gene expression of five of eight of the genes studied (CLOCK, PER2, CRY1, CRY2, and REV-ERBα). BMAL-1 and DBP were significantly affected by cold exposure in the BAT (p < 0.001 and p < 0.0001), but not in soleus muscle (p = 0.479 and p = 0.313). While BMAL-1 was increased by cold exposure especially during the light phase, DBP decreased after the ZT7–10 and ZT11–14 protocols. The negative repressor loop of the core clock (PER1, PER2, CRY1, and CRY2) was upregulated in both BAT and soleus muscle. Interestingly, time-of-day interacted with the increased mRNA expression in three of four of the genes in soleus muscle, but with none of the genes in BAT. Finally, REV-ERBα mRNA expression showed a strong interaction effect of time-of-day and cold exposure in both tissues. In fact, except for ZT3–6, REV-ERBα expression was significantly decreased by cold exposure, with a steep reduction after ZT7–10 in both tissues.

Table 4

n Time-of-day Cold exposure Interaction
BMAL-1 80 <0.0001**** 0.0002*** 0.4584
CLOCK 79 0.0542 0.0018** 0.4676
PER1 82 0.0657 <0.0001**** 0.4003
PER2 80 <0.0001**** <0.0001**** 0.1492
CRY1 82 0.0003*** <0.0001**** 0.0580
CRY2 83 0.2550 <0.0001**** 0.3824
DBP 82 <0.0001*** <0.0001**** <0.0001****
REV-ERBα 8C <0.0001*** 0.3983 0.0018**
CREB 83 0.3688 0.0010** 0.2376
PGCl-α 83 0.0022** <0.0001**** 0.0024**
PGC1-β 81 0.0673 0.8981 0.0568
PPAR-α 82 0.0196* 0.3069 0.0592
PPAR-γ 83 0.1241 0.1606 0.1145
HSP90 83 0.5366 <0.0001**** 0.0321*
UCP1 82 0.0043** <0.0001**** 0.3767
CIRP 81 <0.0001**** <0.0001**** <0.0001****
CPT1-β 81 0.2195 0.0011** 0.2755
AMPK 82 0.0402* 0.0053** 0.2437
FAT/CD36 81 0.2195 0.0038** 0.3579
GLUT4 83 0.4685 0.0310* 0.4725
HSL 81 0.0203* 0.0390* 0.2161
LPL 80 0.0756 0.0001**** 0.4909
ADR-β3 79 0.0572 0.0412* 0.2272
GR 82 0.2009 0.0506 0.2960
ACC1 73 0.4955 0.0272* 0.5446
ACC2 79 0.4957 0.3617 0.0221*

Significance levels observed for the two-way ANOVA analysis of gene expression in the brown adipose tissue.

*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

Table 5

n Time-of-day Cold exposure Interaction
BMAL-1 80 <0.0001**** 0.4785 0.6391
CLOCK 83 0.1868 0.0817 0.0488*
PER1 80 0.0010** <0.0001**** 0.2543
PER2 78 <0.0001**** <0.0001**** 0.0183*
CRY1 81 <0.0001**** <0.0001**** 0.0199*
CRY2 82 0.3535 0.0002*** 0.0016**
DBP 81 <0.0001**** 0.3126 0.1582
REV-ERBα 76 <0.0001**** 0.7068 0.0308*
CREB 83 0.7026 0.0229* 0.0521
PGC1-α 80 <0.0001**** <0.0001**** 0.0006***
PGC1-β 79 0.3707 0.0297* 0.0984
PPAR-α 82 0.0340* 0.0169* 0.0136*
PPAR-γ 80 0.7800 0.0023** 0.3825
HSP90 78 0.0401* 0.1641 0.0015**
UCP2 65 0.5029 0.0189* 0.5634
UCP3 77 <0.0001**** 0.9447 0.0009***
CIRP 83 0.7462 0.4239 0.2013
CPT1-α 79 0.0084** 0.0009*** 0.0285*
CPT1-β 82 0.1321 0.3225 0.0004***
AMPK 83 0.4686 0.0008*** 0.0219*
NAMPT-1 82 0.0655 0.0071** 0.0234*
FAS/CD36 81 0.5285 0.0830 0.0067**
GLUT4 81 0.0071** <0.0001**** 0.0062**
HSL 80 0.3865 <0.0001**** 0.1707
LPL 83 0.4184 0.0200* 0.0125*
ADR-β2 82 0.0170* <0.0001**** 0.2198
GR 83 0.3263 <0.0001**** 0.0642
FOXO-1 80 0.3576 0.0005*** 0.3377
ACC2 13 0.3244 0.0300* 0.0919
mTOR 74 0.5338 0.3857 0.2628
CamK2a 79 0.7143 0.0104* 0.3953

Significance levels observed for the two-way ANOVA analysis of the gene expression in the soleus muscle.

*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

Effect of Time-of-Day on Metabolic Genes mRNA Expression in BAT After Acute Environmental Cooling

We studied 18 genes directly or indirectly related to metabolism or thermogenesis in the BAT (Table 4; Figure 6). Time-of-day affected the mRNA expression of 6 of 18 of the genes studied (PGC1-α, PPAR-α, UCP1, CIRBP, AMPK, and HSL) and just missed significance in 3/18 (PGC1-β, LPL, and ADR-β3). Cold exposure affected most of the BAT genes studied (13 of 18 genes) and just missed significance for GR mRNA expression (p = 0.051). Only 4 of 18 of the investigated genes showed a significant interaction between time-of-day and cold exposure (PGC1-α, HSP90, CIRBP, and ACC2), while PGC1-β and PPAR-α almost reached significance. The genes that were most clearly affected by cold exposure were PGC1-α, HSP90, UCP1, CIRBP, and LPL (p < 0.0001).

Figure 6

Figure 6

Effect of time-of-day on clock gene expression in both soleus muscle and brown adipose tissue (BAT) after exposure to environmental cooling. Light phase began at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) began at ZT12. Results are presented as mean ± SEM. p < 0.05. # indicates differences between control and cold conditions, p < 0.05. n = 5–7/group/ZT.

Effect of Time-of-Day on Metabolic Genes mRNA Expression in Skeletal Muscle After Acute Environmental Cooling

We studied 23 genes directly or indirectly related to metabolism or thermogenesis in the soleus muscle (Table 5; Figure 7). Time-of-day affected the mRNA expression of 7 of 23 genes (PGC1-α, PPAR-α, HSP90, UCP3, CPT1-α, GLUT4, and ADR-β2) and just missed significance for NAMPT-1 expression (p = 0.066). Cold exposure affected mRNA expression of 17 of 23 of the studied genes and almost reached significance for FAT/CD36 levels (p = 0.083). The interaction between time-of-day and cold exposure was significant for 11 of 23 of the studied genes and just missed significance for 4 of 23 genes (CREB, PGC1-β, GR, and ACC2). The genes that were most clearly affected by cold exposure were PGC1-α, GLUT4, HSL, ADRβ2, and GR (p < 0.0001).

Figure 7

Figure 7

Effect of time-of-day on metabolic gene expression in brown adipose tissue (BAT) after exposure to environmental cooling. Light phase began at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) began at ZT12. Results are presented as mean ± SEM. # indicates differences between control and cold conditions, p < 0.05. n = 5–7/group/ZT. Data for CREB, PGC1-β, PPAR-α, PPAR-γ, ACC1, and ACC2 are not shown (refer to Table 4 for results).

Discussion

The main finding of the present study is that the cold-induced metabolic response and changes in gene expression in BAT and muscle differ depending on the time-of-day of the cold exposure. For instance, the cold-induced increase in lipid oxidation was mainly observed during the dark phase (Figures 5E,F), indicating a daily modulation of the cold-induced metabolic adaptations. In the soleus skeletal muscle, particularly active during shivering thermogenesis, cold exposure increased the expression of clock genes in both the negative and positive loop of the core clock mechanism, PER/CRY and BMAL-1/CLOCK, respectively (Figure 6). In addition, besides upregulating clock gene expression in both regulatory loops of the BAT molecular clock as well, cold exposure inhibited the expression of DBP and REV-ERBα in BAT (Figure 6). The current findings add further evidence for a tissue-specific action of the internal clock in peripheral tissues such as the skeletal muscle and BAT. Whether this differential modulation in response to environmental stimuli relies on the activity of the central clock, extra-SCN sites within the central nervous system or peripheral inputs remains to be further elucidated.

Effects of Cold Exposure During Different Times-of-Day on Thermogenesis and Substrate Oxidation

The thermal and metabolic adjustments induced by an acute exposure to a cold environment are well known for mice (2429), rats (6, 8, 3032), and humans (9, 3335). Our experiments confirmed the effects of cold on the induction of increased heat production mostly through the concerted modulation of behavioral (LA, food, and water intake) and metabolic (Tcore, VO2, RER, EE, CHO, and lipid oxidation) processes. In accordance with our hypothesis, we found that time-of-day modulated the thermal, behavioral, and metabolic responses.

Significant decrements of ambient temperature pose a primitive threat to body temperature regulation. In order to keep body temperature within physiological levels, the central nervous system modulates a number of physiological processes involved in heat conservation and production. In the present experiment, we observed that rat Tcore is resilient to a temporally limited exposure to a reduced ambient temperature. In fact, after 120 min of environmental cooling, normal Tcore was preserved for each time-of-day exposure (Figure 3A). Curiously, we even detected a small rise in Tcore when the animals were challenged with a cold environment presented during the light phase, an effect that is not observed in smaller rodents, such as mice (2429). This mismatch of heat production and heat dissipation (leading to body heating) might be due to a heat defensive state potentiated by the sudden decrease in ambient temperature during the sleep period.

Indeed, a lower Tcore during the light phase is maintained through tail vasoconstriction and decreased EE (3639). In such a condition, the Tcore regulatory system likely presents a higher sensitivity to changes in locomotor and metabolic activity (40), resulting in increased Tcore during the first hour of cold exposure in the light phase, when heat production was increased (Figure 3A). Conversely, during the dark phase, basal Tcore is slightly raised as a function of the increased LA, EE, and circadian rhythm (39), making the immediate impact of the increased metabolic rate on Tcore less perceptible.

Indeed, the well-described circadian-dependent decrease in Tcore during the second half of dark phase seems to be counterbalanced by heat production induced by cold exposure (Figure 3A, ZT15–17). This cold-induced hyperthermia was also reported by others (6, 31, 4143). Initially, heat dissipation is reduced to minimize heat loss to the colder environment (6). As cold exposure is maintained, heat production is increased and heat loss and production reach a steady state in which Tcore can be successfully preserved. Therefore, as previously observed with other stressors (40, 4447), it seems that the time-of-day-dependent effects of cold exposure on body temperature reflect the transitory disturbance in heat loss and heat production mechanisms elicited by mixed signals from the internal circadian time and the thermoregulatory pathways that participate in the homeostatic responses.

In line with this perspective, the increased metabolic rate (Figures 3C,D) and LA (Figures 4A,B) in response to cold exposure might also reflect environmental temperature gradient (intensity of the stimulus) rather than a time-of-day-dependent modulation. Interestingly, LA spiked during the first hour of the cooling protocol but decreased thereafter (Figure 4A). Similar results were described by others (43, 48, 49), suggesting that this probably reflects coupled mechanisms of heat conservation (cold avoidance behavior reflecting increased LA) and heat production (shivering thermogenesis reflecting the decreased LA). In fact, as the environment gets colder, the animals lessen their LA without dampening metabolic rate and EE (Figures 3C,E), which points toward increased shivering activity of thermogenic pathways. Moreover, both oxygen consumption and calculated EE reached steady states of increased activity regardless of time-of-day. This is in accordance with previous experiments in mice showing that the relative changes in metabolic rate were similar when the cooling protocol started at the beginning of the light or dark phase (27, 29). Interestingly, Tokizawa and colleagues (29) observed that the threshold for increased heat production was elevated during the dark phase, while behavioral curling was elicited earlier during the light phase. In rats, a similar effect of time-of-day on thermoregulatory thresholds was also observed during exercise (47). To our knowledge, this is the first time that a time-of-day effect on the cold defensive response is reported for the rat model. With the present results, it seems that the circadian system regulates the basal settings of body temperature, thereby indirectly establishing the level of the homeostatic response required to minimize or even neutralize the physical challenge on thermal homeostasis posed by the reduced ambient temperature.

Even though thermal balance was successfully defended during cold exposure regardless of time-of-day, lipid and CHO utilization were affected by both environmental time and temperature (Figure 5). It has been previously shown that environmental cooling provides a key signal to substrate utilization for shivering (involuntary muscle contractions) and non-shivering (generated by proton gradient within mitochondrial oxidative activity and controlled by uncoupling proteins) thermogenesis from skeletal muscle and BAT, respectively (6, 8, 30, 5052). Environmental conditions and the intensity of cold exposure determine the metabolic rate and the prevalence of the thermogenic pathway activated to generate heat (6, 8, 51). Vaillancourt and colleagues (8) have previously shown that although affecting total EE; below a certain threshold (~15°C), the higher rates of CHO, lipid, and protein oxidations are not intensified by even lower environmental temperatures (10 and 5°C). To our knowledge, we are the first to report a clear time-of-day effect on substrate utilization during cold exposure.

Interestingly, we observed that during the active period, cold exposure caused a major switch by increasing lipid oxidation (Figures 5E,F) and lessening CHO oxidation (Figures 5C,D), whereas during the light phase both lipid and CHO oxidation rates were increased, although only transitory for CHO. This might reflect the first stage of cold defensive mechanisms in which in skeletal muscle shivering first consumes carbohydrates stocks in the glycolytic pathway to promptly generate ATP for muscular contraction and heat production (51, 52). As cold exposure persists, the participation of shivering in the thermal balance decreases and a primary role is appropriated by BAT non-shivering thermogenesis, with increased lipid uptake and oxidation (35, 52). Therefore, time-of-day influences the balance between CHO and lipids utilization during cold exposure, providing further evidence that the daily rhythms expressed through an internal clock in different body tissues might influence metabolic adjustments in heat-generating tissues.

Differences in the Metabolic and Thermal Adjustments to Cold Were Related to Distinct Changes in mRNA Expression of Clock and Clock-Controlled Genes

To increase metabolic heat production during exposure to a lower environmental temperature, activity in skeletal muscle and BAT is intensified through different pathways within the central nervous system (6). In the present study, we observed that after 3 h of cold exposure 13 of 18 and 17 of 23 of the chosen genes thought to be involved in cold-induced metabolic heat production were affected by our cold protocol in the BAT (Figure 7) and skeletal muscle (Figure 8), respectively.

Figure 8

Figure 8

Effect of time-of-day on metabolic gene expression in soleus muscle after exposure to environmental cooling. Light phase began at Zeitgeber Time 0 (ZT0) and dark phase (shaded areas) began at ZT12. Results are presented as mean ± SEM. # indicates differences between control and cold conditions, p < 0.05. n = 5–7/group/ZT. Data for CREB, PGC1-β, PPAR-γ, UCP2, CIRBP, CPT-β, NAMPT-1, FOXO-1, ACC2, CamK2a, and mTOR are not shown (refer to Table 5 for results).

In addition to the metabolic effects, cold exposure also affected clock and clock-controlled gene expression in BAT and soleus muscle, although to a different extent (Figure 6). In the BAT, seven of eight core clock genes studied were altered by the lower environmental temperature. Specifically, five of those genes were upregulated after cold exposure, without an interaction with time-of-day on this response. Others have already implicated PER1, PER2, CRY1, and CRY2 in BAT in cold exposure (12, 13). In fact, a functional clock is necessary to mediate molecular adjustments in the mitochondrial activity program (13, 17, 53, 54). However, the present study is the first to show this upregulation of BAT activity during cold exposure at different moments of the day, reinforcing the hypothesis that the thermogenic program, although it is intermingled with the core clock machinery (55), relies on the intensity of the stimulus rather than on time-of-day.

Interestingly, in BAT two of the clock genes, DBP and REV-ERBα, were downregulated by cold exposure with a time-of-day interaction. REV-ERBα is also downregulated after increased energetic demands and is acknowledged as an important repressor of mitochondrial activity and seems to be regulated by the positive loop of the core clock (16, 17, 53). In fact, we observed that the reduced REV-ERBα expression is accompanied by an augmented expression of BMAL-1 and PGC1-α, especially during the light phase (Figures 6 and 7), which points toward a pronounced shift in BAT mitochondrial activity in response to the light-to-dark transition (54) and change in ambient temperature, probably through a temperature compensation mechanism.

Regarding the soleus skeletal muscle, we observed that only the genes from the negative loop of the canonical core clock machinery were affected by cold exposure (Figure 6; Table 5), i.e., PER1, PER2, CRY1, and CRY2, that is half of the studied clock genes (four of eight). In line with this, resetting of the peripheral clock by cold exposure (13) and exercise (56) seems to be dependent of Per1/2 (57), as is resetting by glucocorticoids fluctuations (58). These findings reinforce the idea that cold exposure can reset the peripheral clock by changes in metabolic activity, since genes encoding for mitochondrial activity, such as CPT1, PGC1α, PPAR, and FOXO1 were affected. In addition, genes that regulate substrate uptake (LPL and GLUT4), cellular energetic state (AMPK), and mitochondrial activity (CPT1, NAMPT, PGC1α and PPARs) showed an interaction with time-of-day, suggesting that key components of the cellular energetic balance are dependent on temporal input to establish the necessary shifts in molecular programs required to survive in a colder environment. On the other hand, cold exposure also changed GR and ADRβ receptor expression in both soleus muscle and BAT, indicating that also hormonal and autonomic, i.e., non-metabolic, factors might be involved in this shifting process.

In the skeletal muscle (Figure 8), the increased expression of LPL, HSL, and FAT/CD36 points toward boosted lipid uptake, which is supported by previous physiological and molecular expression data (59). Glucose uptake was also elevated by cold exposure as suggested by increased GLUT4 expression in the present experiment and those performed by others (10, 11). Taken together, increased lipid and glucose uptake suggests an augmented substrate oxidative state within the soleus muscle. Reinforcing this assumption, the expression of key metabolic regulators such as AMPK, CamK2a, and NAMPT-1 was also increased after cold exposure (Table 5). It is thought that AMPK activation due to increased energetic demand regulates the key mitochondrial transporter for β-oxidation, CPT1-α/β, as well as key transcription factors for mitochondrial activity, such as PPARs, PGC1-α, and FOXO1 (6065). In fact, all these genes associated with mitochondrial activity were upregulated in the skeletal muscle. The present gene expression results thus support a role for the skeletal muscle in metabolic heat production. This idea is also supported by increased HSP90 mRNA expression, which might be related to higher oxidative activity demanding more chaperone protein content for protein stabilization (66).

Interestingly, despite the higher energetic demand, we did not observe a clear effect of cold exposure on UCP3 activity (Figure 8). It is thought that UCP3 could be directly involved in muscular heat production independent of shivering (35, 67). In line with our results, others (68) have shown that muscular UCP3 activity might not be involved in cold-induced metabolic heat production. Instead, uncoupling seems to constitute an intramuscular FFA buffering system during cold exposure, at least in murine models (12, 69). In addition, we also observed a possible influence of neuronal and hormonal stimulation on these molecular adjustments since expression of ADRβ2 and GR was upregulated after cold exposure, pointing to a direct modulation of the cold-induced adjustments by the sympatho-cortico-adrenal system (12, 7072).

In the BAT, cold exposure also upregulated the substrate uptake and mobilization program, as can be concluded from the increased mRNA expression of LPL, FAT/CD36, and GLUT4 (Figure 7; Table 4). Transcription of key enzymes for β-oxidation, such as CPT1-β and ACC1 was also upregulated by cold exposure. Aligned with this finding, mitochondrial activity, and biogenesis, expressed by the function of the key transcription factor PGC1α was increased by cold exposure. This increased mitochondrial activity state is reinforced by improved UCP1, increased HSP90 (66) and decreased CIRPB (73) expression after cold exposure. AMPK activity, which is thought to play a key role in mediating cellular metabolic flux, was further stimulated by cold exposure. Increased BAT activity is well supported by previous studies and accumulating evidence implicates both glucose and lipid oxidation in BAT non-shivering thermogenesis (14, 7477).

It is clear that the initial metabolic state, as defined by the circadian system, is crucial for the physiological response to environmental stimuli and therefore should be taken into account to understand the intrinsic capacity to adapt to changes in ambient temperature. In the current experiment, the cold-induced metabolic shift toward lipid oxidation provided a clear example of such a time-dependency and further proof that the internal clock plays an important role in shaping such physiological responses. In the soleus skeletal muscle, cold exposure upregulated specifically the negative loop of the canonical clock (PER/CRY), whereas in the BAT both the negative and positive loop (BMAL-1/CLOCK) were upregulated. Thereby, the present study for the first time showed a tissue-specific effect of cold exposure on clock and clock-controlled genes and thus demonstrated the existence of an interaction between time-of-day and homeostatic adjustments elicited by acute cold exposure.

Statements

Ethics statement

This study was carried out in accordance with the recommendations of the animal care committee of the Royal Netherlands Academy of Arts and Sciences (DEC/KNAW). The protocol was approved by the animal care committee of the Royal Netherlands Academy of Arts and Sciences (DEC/KNAW).

Author contributions

FM, CC, and AK designed the experiments. FM, ZZ, YS, PG, RJ, and EF acquired the data. FM, ZZ, EF, CC, and AK contributed to the analysis and interpretation of data. FM and AK participated in the elaboration of the manuscript and gave final approval for submission and publication, being accountable for all aspects of the present work.

Acknowledgments

We acknowledge Unga A. Unmehopa for her assistance on the quality control of RNA isolation and RT-qPCR. FM was supported by a CAPES/PDSE grant (BEX 9894/14-4). PG was supported by a ZonMW TOP grant (#91214047). ZZ was supported by the China Exchange Program of the Royal Netherlands Academy of Sciences, a Chinese Academy of Sciences Grant (11CDP001) and a Chinese Scholarship Council Grant (201206340004).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fendo.2018.00199/full#supplementary-material.

References

  • 1

    Brown SA Zumbrunn G Fleury-Olela F Preitner N Schibler U . Rhythms of mammalian body temperature can sustain peripheral circadian clocks. Curr Biol (2002) 12(18):157483.10.1016/S0960-9822(02)01145-4

  • 2

    Rensing L . Ruoff: temperature effect on entrainment, phase shifting, and amplitude of circadian clocks and its molecular bases. Chronobiol Int (2002) 19(5):80764.10.1081/CBI-120014569

  • 3

    El Allali K Achaâban MR Bothorel B Piro M Bouâouda H El Allouchi M et al Entrainment of the circadian clock by daily ambient temperature cycles in the camel (Camelus dromedarius). Am J Physiol Regul Integr Comp Physiol (2013) 304(11):R104452.10.1152/ajpregu.00466.2012

  • 4

    Buhr ED Yoo SH Takahashi JS . Temperature as a universal resetting cue for mammalian circadian oscillators. Science (2010) 330(6002):37985.10.1126/science.1195262

  • 5

    Saini C Morf J Stratmann M Gos P Schibler U . Simulated body temperature rhythms reveal the phase-shifting behavior and plasticity of mammalian circadian oscillators. Genes Dev (2012) 26(6):56780.10.1101/gad.183251.111

  • 6

    Nakamura K Morrison SF . Central efferent pathways for cold-defensive and febrile shivering. J Physiol (2011) 589(Pt 14):364158.10.1113/jphysiol.2011.210047

  • 7

    Morrison SF . Central neural control of thermoregulation and brown adipose tissue. Auton Neurosci (2016) 196:1424.10.1016/j.autneu.2016.02.010

  • 8

    Vaillancourt E Haman F Weber JM . Fuel selection in Wistar rats exposed to cold: shivering thermogenesis diverts fatty acids from re-esterification to oxidation. J Physiol (2009) 587(Pt 17):434959.10.1113/jphysiol.2009.175331

  • 9

    Weber JM Haman F . Fuel selection in shivering humans. Acta Physiol Scand (2005) 184(4):31929.10.1111/j.1365-201X.2005.01465.x

  • 10

    Hanssen MJ Hoeks J Brans B van der Lans AA Schaart G van den Driessche JJ et al Short-term cold acclimation improves insulin sensitivity in patients with type 2 diabetes mellitus. Nat Med (2015) 21(8):8635.10.1038/nm.3891

  • 11

    Oliveira RL Ueno M de Souza CT Pereira-da-Silva M Gasparetti AL Bezzera RM et al Cold-induced PGC-1alpha expression modulates muscle glucose uptake through an insulin receptor/Akt-independent, AMPK-dependent pathway. Am J Physiol Endocrinol Metab (2004) 287(4):E68695.10.1152/ajpendo.00103.2004

  • 12

    Li S Yu Q Wang GX Lin JD . The biological clock is regulated by adrenergic signaling in brown fat but is dispensable for cold-induced thermogenesis. PLoS One (2013) 8(8):e70109.10.1371/journal.pone.0070109

  • 13

    Chappuis S Ripperger JA Schnell A Rando G Jud C Wahli W et al Role of the circadian clock gene Per2 in adaptation to cold temperature. Mol Metab (2013) 2(3):18493.10.1016/j.molmet.2013.05.002

  • 14

    Lee P Bova R Schofield L Bryant W Dieckmann W Slattery A et al Brown adipose tissue exhibits a glucose-responsive thermogenic biorhythm in humans. Cell Metab (2016) 23(4):6029.10.1016/j.cmet.2016.02.007

  • 15

    Lowell BB Spiegelman BM . Towards a molecular understanding of adaptive thermogenesis. Nature (2000) 404(6778):65260.10.1038/35007527

  • 16

    Delezie J Dumont S Dardente H Oudart H Gréchez-Cassiau A Klosen P et al The nuclear receptor REV-ERB alpha is required for the daily balance of carbohydrate and lipid metabolism. FASEB J (2012) 26(8):332135.10.1096/fj.12-208751

  • 17

    Woldt E Sebti Y Solt LA Duhem C Lancel S Eeckhoute J et al Rev-erb-alpha modulates skeletal muscle oxidative capacity by regulating mitochondrial biogenesis and autophagy. Nat Med (2013) 19(8):103946.10.1038/nm.3213

  • 18

    Deurveilher S Semba K . Indirect projections from the suprachiasmatic nucleus to the median preoptic nucleus in rat. Brain Res (2003) 987(1):1006.10.1016/S0006-8993(03)03295-5

  • 19

    Saper CB . The central circadian timing system. Curr Opin Neurobiol (2013) 23(5):74751.10.1016/j.conb.2013.04.004

  • 20

    Reppert SM Weaver DR . Coordination of circadian timing in mammals. Nature (2002) 418(6901):93541.10.1038/nature00965

  • 21

    Frayn KN . Calculation of substrate oxidation rates in vivo from gaseous exchange. J Appl Physiol Respir Environ Exerc Physiol (1983) 55(2):62834.

  • 22

    Péronnet F Massicotte D . Table of nonprotein respiratory quotient: an update. Can J Sport Sci (1991) 16(1):239.

  • 23

    Su Y Foppen E Fliers E Kalsbeek A . Effects of intracerebroventricular administration of neuropeptide y on metabolic gene expression and energy metabolism in male rats. Endocrinology (2016) 157(8):307085.10.1210/en.2016-1083

  • 24

    Tokizawa K Uchida Y Nagashima K . Thermoregulation in the cold changes depending on the time-of-day and feeding condition: physiological and anatomical analyses of involved circadian mechanisms. Neuroscience (2009) 164(3):137786.10.1016/j.neuroscience.2009.08.040

  • 25

    Virtue S Even P Vidal-Puig A . Below thermoneutrality, changes in activity do not drive changes in total daily energy expenditure between groups of mice. Cell Metab (2012) 16(5):66571.10.1016/j.cmet.2012.10.008

  • 26

    Uchida Y Tokizawa K Nagashima K . Characteristics of activated neurons in the suprachiasmatic nucleus when mice become hypothermic during fasting and cold exposure. Neurosci Lett (2014) 579:17782.10.1016/j.neulet.2014.07.039

  • 27

    Abreu-Vieira G Xiao C Gavrilova O Reitman ML . Integration of body temperature into the analysis of energy expenditure in the mouse. Mol Metab (2015) 4(6):46170.10.1016/j.molmet.2015.03.001

  • 28

    Sato N Marui S Ozaki M Nagashima K . Cold exposure and/or fasting modulate the relationship between sleep and body temperature rhythms in mice. Physiol Behav (2015) 149:6975.10.1016/j.physbeh.2015.05.024

  • 29

    Tokizawa K Yoda T Uchida Y Kanosue K Nagashima K . Estimation of the core temperature control during ambient temperature changes and the influence of circadian rhythm and metabolic conditions in mice. J Therm Biol (2015) 51:4754.10.1016/j.jtherbio.2015.03.005

  • 30

    Davis TR Johnston DR Bell FC Cremer BJ . Regulation of shivering and non-shivering heat production during acclimation of rats. Am J Physiol (1960) 198:4715.

  • 31

    Morimoto A Murakami N Nakamori T Watanabe T . Suppression of non-shivering thermogenesis in the rat by heat-seeking behaviour during cold exposure. J Physiol (1986) 380:5419.10.1113/jphysiol.1986.sp016301

  • 32

    Gautier H . Body temperature regulation in the rat. J Therm Biol (2000) 25(4):2739.10.1016/S0306-4565(99)00097-2

  • 33

    Lennon JA Brech WJ Gordon ES . Effect of a short period of cold exposure on plasma FFA level in lean and obese humans. Metabolism (1967) 16(6):5036.10.1016/0026-0495(67)90078-9

  • 34

    Haman F . Shivering in the cold: from mechanisms of fuel selection to survival. J Appl Physiol (1985) (2006) 100(5):17028.10.1152/japplphysiol.01088.2005

  • 35

    Wijers SL Schrauwen P Saris WH van Marken Lichtenbelt WD . Human skeletal muscle mitochondrial uncoupling is associated with cold induced adaptive thermogenesis. PLoS One (2008) 3(3):e1777.10.1371/journal.pone.0001777

  • 36

    Fuller CA Sulzman FM Moore-Ede MC . Role of heat loss and heat production in generation of the circadian temperature rhythm of the squirrel monkey. Physiol Behav (1985) 34(4):5436.10.1016/0031-9384(85)90046-0

  • 37

    Gordon CJ Puckett E Padnos B . Rat tail skin temperature monitored noninvasively by radiotelemetry: characterization by examination of vasomotor responses to thermomodulatory agents. J Pharmacol Toxicol Methods (2002) 47(2):10714.10.1016/S1056-8719(02)00219-8

  • 38

    Refinetti R . Metabolic heat production, heat loss and the circadian rhythm of body temperature in the rat. Exp Physiol (2003) 88(3):4239.10.1113/eph8802521

  • 39

    Refinetti R . The circadian rhythm of body temperature. Front Biosci (Landmark Ed) (2010) 15:56494.10.2741/3634

  • 40

    Machado FS Rodovalho GV Coimbra CC . The time-of-day differently influences fatigue and locomotor activity: is body temperature a key factor?Physiol Behav (2015) 140:814.10.1016/j.physbeh.2014.11.069

  • 41

    Ishiwata T Hasegawa H Yazawa T Otokawa M Aihara Y . Functional role of the preoptic area and anterior hypothalamus in thermoregulation in freely moving rats. Neurosci Lett (2002) 325(3):16770.10.1016/S0304-3940(02)00266-5

  • 42

    Ishiwata T Saito T Hasegawa H Yazawa T Otokawa M Aihara Y . Changes of body temperature and extracellular serotonin level in the preoptic area and anterior hypothalamus after thermal or serotonergic pharmacological stimulation of freely moving rats. Life Sci (2004) 75(22):266575.10.1016/j.lfs.2004.04.040

  • 43

    Mohammed M Yanagisawa M Blessing W Ootsuka Y . Attenuated cold defense responses in orexin neuron-ablated rats. Temperature (Austin) (2016) 3(3):46575.10.1080/23328940.2016.1184366

  • 44

    Wachulec M Li H Tanaka H Peloso E Satinoff E . Suprachiasmatic nuclei lesions do not eliminate homeostatic thermoregulatory responses in rats. J Biol Rhythms (1997) 12(3):22634.10.1177/074873049701200304

  • 45

    Kalsbeek A Liu J Lei J Timmermans L Foppen E Cailotto C et al Differential involvement of the suprachiasmatic nucleus in lipopolysaccharide-induced plasma glucose and corticosterone responses. Chronobiol Int (2012) 29(7):83549.10.3109/07420528.2012.699123

  • 46

    Kalsbeek A Ruiter M La Fleur SE Van Heijningen C Buijs RM . The diurnal modulation of hormonal responses in the rat varies with different stimuli. J Neuroendocrinol (2003) 15(12):114455.10.1111/j.1365-2826.2003.01112.x

  • 47

    Machado FS Fóscolo DR Poletini MO Coimbra CC . Influence of time-of-day on maximal exercise capacity is related to daily thermal balance but not to induced neuronal activity in rats. Front Physiol (2016) 7:464.10.3389/fphys.2016.00464

  • 48

    Finger FW . Relation of general activity in rats to environmental temperature. Percept Mot Skills (1976) 43(3 Pt 1):87590.10.2466/pms.1976.43.3.875

  • 49

    Ishiwata T Saito T Hasegawa H Yazawa T Kotani Y Otokawa M et al Changes of body temperature and thermoregulatory responses of freely moving rats during GABAergic pharmacological stimulation to the preoptic area and anterior hypothalamus in several ambient temperatures. Brain Res (2005) 1048(1–2):3240.10.1016/j.brainres.2005.04.027

  • 50

    Banet M Hensel H Liebermann H . The central control of shivering and non-shivering thermogenesis in the rat. J Physiol (1978) 283:56984.10.1113/jphysiol.1978.sp012520

  • 51

    Gordon CJ . Thermal biology of the laboratory rat. Physiol Behav (1990) 47(5):96391.10.1016/0031-9384(90)90025-Y

  • 52

    Adán C Ardévol A Remesar X Alemany M Fernández-López JA . Hindleg muscle energy and substrate balances in cold-exposed rats. J Exp Biol (1995) 198(Pt 6):124351.

  • 53

    Gerhart-Hines Z Feng D Emmett MJ Everett LJ Loro E Briggs ER . The nuclear receptor Rev-Erb alpha controls circadian thermogenic plasticity. Nature (2013) 503(7476):4103.10.1038/nature12642

  • 54

    Orozco-Solis R Aguilar-Arnal L Murakami M Peruquetti R Ramadori G Coppari R et al The circadian clock in the ventromedial hypothalamus controls cyclic energy expenditure. Cell Metab (2016) 23(3):4678.10.1016/j.cmet.2016.02.003

  • 55

    Gerhart-Hines Z Lazar MA . Circadian metabolism in the light of evolution. Endocr Rev (2015) 36(3):289304.10.1210/er.2015-1007

  • 56

    Wolff G Esser KA . Scheduled exercise phase shifts the circadian clock in skeletal muscle. Med Sci Sports Exerc (2012) 44(9):166370.10.1249/MSS.0b013e318255cf4c

  • 57

    Tahara Y Aoyama S Shibata S . The mammalian circadian clock and its entrainment by stress and exercise. J Physiol Sci (2017) 67(1):110.10.1007/s12576-016-0450-7

  • 58

    Balsalobre A Brown SA Marcacci L Tronche F Kellendonk C Reichardt HM . Resetting of circadian time in peripheral tissues by glucocorticoid signaling. Science (2000) 89(5488):23447.10.1126/science.289.5488.2344

  • 59

    Koonen DP Benton CR Arumugam Y Tandon NN Calles-Escandon J Glatz JF et al Different mechanisms can alter fatty acid transport when muscle contractile activity is chronically altered. Am J Physiol Endocrinol Metab (2004) 286(6):E10429.10.1152/ajpendo.00531.2003

  • 60

    Wu Z Puigserver P Andersson U Zhang C Adelmant G Mootha V et al Mechanisms controlling mitochondrial biogenesis and respiration through the thermogenic coactivator PGC-1. Cell (1999) 98(1):11524.10.1016/S0092-8674(00)80611-X

  • 61

    Freyssenet D . Energy sensing and regulation of gene expression in skeletal muscle. J Appl Physiol (1985) (2007) 102(2):52940.10.1152/japplphysiol.01126.2005

  • 62

    Nickerson JG Momken I Benton CR Lally J Holloway GP Han XX . Protein-mediated fatty acid uptake: regulation by contraction, AMP-activated protein kinase, and endocrine signals. Appl Physiol Nutr Metab (2007) 32(5):86573.10.1139/H07-084

  • 63

    Rohas LM St-Pierre J Uldry M Jäger S Handschin C Spiegelman BM . A fundamental system of cellular energy homeostasis regulated by PGC-1alpha. Proc Natl Acad Sci U S A (2007) 104(19):79338.10.1073/pnas.0702683104

  • 64

    Schiaffino S Sandri M Murgia M . Activity-dependent signaling pathways controlling muscle diversity and plasticity. Physiology (Bethesda) (2007) 22:26978.10.1152/physiol.00009.2007

  • 65

    Lee Y Kim EK . AMP-activated protein kinase as a key molecular link between metabolism and clockwork. Exp Mol Med (2013) 45:e33.10.1038/emm.2013.65

  • 66

    Matz JM Blake MJ Tatelman HM Lavoi KP Holbrook NJ . Characterization and regulation of cold-induced heat shock protein expression in mouse brown adipose tissue. Am J Physiol (1995) 269(1 Pt 2):R3847.

  • 67

    Schrauwen P Hesselink M . Uncoupling protein 3 and physical activity: the role of uncoupling protein 3 in energy metabolism revisited. Proc Nutr Soc (2003) 62(3):63543.10.1079/PNS2003277

  • 68

    Boss O Samec S Kühne F Bijlenga P Assimacopoulos-Jeannet F Seydoux J et al Uncoupling protein-3 expression in rodent skeletal muscle is modulated by food intake but not by changes in environmental temperature. J Biol Chem (1998) 273(1):58.10.1074/jbc.273.1.5

  • 69

    Lin B Coughlin S Pilch F . Bidirectional regulation of uncoupling protein-3 and GLUT-4 mRNA in skeletal muscle by cold. Am J Physiol (1998) 275(3 Pt 1):E38691.

  • 70

    Navegantes LC Resano NM Migliorini RH Kettelhut IC . Catecholamines inhibit Ca(2+)-dependent proteolysis in rat skeletal muscle through beta(2)-adrenoceptors and cAMP. Am J Physiol Endocrinol Metab (2001) 281(3):E44954.10.1152/ajpendo.2001.281.3.E449

  • 71

    Braun TP Marks DL . Hypothalamic regulation of muscle metabolism. Curr Opin Clin Nutr Metab Care (2011) 14(3):23742.10.1097/MCO.0b013e328345bbcd

  • 72

    Khan MM Lustrino D Silveira WA Wild F Straka T Issop Y et al Sympathetic innervation controls homeostasis of neuromuscular junctions in health and disease. Proc Natl Acad Sci U S A (2016) 113(3):74650.10.1073/pnas.1524272113

  • 73

    Kornmann B Schaad O Bujard H Takahashi JS Schibler U . System-driven and oscillator-dependent circadian transcription in mice with a conditionally active liver clock. PLoS Biol (2007) 5(2):e34.10.1371/journal.pbio.0050034

  • 74

    Puigserver P Wu Z Park CW Graves R Wright M Spiegelman BM . A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell (1998) 92(6):82939.10.1016/S0092-8674(00)81410-5

  • 75

    Puigserver P Spiegelman BM . Peroxisome proliferator-activated receptor gamma coactivator 1 alpha (PGC-1 alpha): transcriptional coactivator and metabolic regulator. Endocr Rev (2003) 24(1):7890.10.1210/er.2002-0012

  • 76

    Bézaire V Seifert EL Harper ME . Uncoupling protein-3: clues in an ongoing mitochondrial mystery. FASEB J (2007) 21(2):31224.10.1096/fj.06-6966rev

  • 77

    Azzu V Jastroch M Divakaruni AS Brand MD . The regulation and turnover of mitochondrial uncoupling proteins. Biochim Biophys Acta (2010) 1797(6–7):78591.10.1016/j.bbabio.2010.02.035

Summary

Keywords

circadian, locomotor activity, thermoregulation, thermogenesis, gene expression, skeletal muscle, brown adipose tissue

Citation

Machado FSM, Zhang Z, Su Y, de Goede P, Jansen R, Foppen E, Coimbra CC and Kalsbeek A (2018) Time-of-Day Effects on Metabolic and Clock-Related Adjustments to Cold. Front. Endocrinol. 9:199. doi: 10.3389/fendo.2018.00199

Received

22 December 2017

Accepted

10 April 2018

Published

26 April 2018

Volume

9 - 2018

Edited by

Arturo Ortega, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV-IPN), Mexico

Reviewed by

Denis Richard, Laval University, Canada; Ricardo Orozco Solis, National Institute of Genomic Medicine, Mexico

Updates

Copyright

*Correspondence: Andries Kalsbeek,

Specialty section: This article was submitted to Neuroendocrine Science, a section of the journal Frontiers in Endocrinology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics