What this is
- This research investigates the circadian clock in reindeer fibroblasts, exploring its rhythmicity and .
- Reindeer fibroblasts exhibit a shorter free-running period compared to mouse fibroblasts, with unique responses to temperature cycles.
- The study suggests weak coupling between core and secondary molecular clock loops in reindeer, affecting circadian regulation.
Essence
- Reindeer fibroblasts display rhythmic circadian activity and , yet their response to temperature cycles is atypical compared to mouse fibroblasts, indicating weak coupling in their molecular clock.
Key takeaways
- Reindeer fibroblasts have a free-running period of 21 h 19 min, shorter than the 25 h 20 min in mouse fibroblasts. This difference suggests a unique in reindeer adapted to Arctic conditions.
- Temperature cycles significantly influence clock gene expression in reindeer fibroblasts, which show a 9 to 12 h phase shift when returning to constant conditions. This behavior contrasts with mouse fibroblasts, which maintain their phase relationship.
- Weak coupling between core and secondary circadian feedback loops in reindeer fibroblasts may lead to direct temperature effects dominating over circadian transcription, affecting their overall circadian regulation.
Caveats
- The study's findings are based on in vitro experiments, which may not fully replicate the in vivo conditions of reindeer physiology.
- Differences in experimental protocols between studies may contribute to discrepancies in observed circadian behavior, warranting further investigation.
Definitions
- circadian rhythm: A biological process that displays an endogenous, entrainable oscillation of about 24 hours, influencing various physiological functions.
- temperature compensation: The ability of a circadian clock to maintain consistent period lengths across varying temperatures.
Simplified
Materials and Methods
Animals
All animals were kept in accordance with the EU directive 201/63/EU under licenses provided by the Norwegian food safety authority (Mattilsynet, FOTS 28929). Reindeer-derived biological material stemmed exclusively from a semi-domesticated herd held at the University of Tromsø, Norway (69°N).
Fibroblast Cell Cultures
Fibroblasts were cultivated from an ear skin sample of euthanized reindeer (Rangifer tarandus tarandus, female, 9 months old, collected on 27/02/2019 when natural photoperiod from start to end of civil twilight was 11 h 15 min) and from an ear puncture of a laboratory mouse (Mus musculus, C57Bl6/J, males, 10 weeks old, provided by Jaione Simón-Santamaría and Karolina Szafranska, UiT). Sampling and cultivation of fibroblasts broadly followed the protocols outlined in Du and Brown (2021).
In brief, tissue samples were collected in sampling medium containing Dulbecco’s Modified Eagle Medium (DMEM; Merck D5796), 50% fetal bovine serum (FBS; Biowest S181B), 1% penicillin-streptomycin (P/S; Merck P4333), and 1x amphotericin B on ice, then incubated in 2 mL digestion medium containing DMEM, 10% FBS, 1% P/S, 1x amphotericin B, and 3.5U liberase for 4 h at 37 °C and 5% CO2. Digested fragments were transferred into a 50 mL falcon tube filled with warm phosphate-buffered saline (PBS; Biowest L0615) and centrifuged for 5 min at 200 × g. The supernatant was removed and the formed pellet was resuspended in 0.2 mL culture medium (DMEM + 10% FBS + 1% P/S). The mix was placed in the center of a fresh 6-well plate, overlayed with a Millicell cell culture insert (Merck; PICM0RG50) and provided with 2 mL culture medium. The sample was incubated for ca. 2 weeks at 37 °C and 5% CO2 with a medium change every 3 to 4 days. Ready-to-harvest cells were trypsinized (Merck T4049) and replated. Cells were passaged until an appropriate volume was achieved. Active cultures were kept in DMEM, 10% FBS and 1% P/S and housed in cell incubators at 37 °C, 5% CO2 and high humidity. The primary fibroblasts from mouse and reindeer were not immortalized. Cultures not in active use were stored in liquid nitrogen in culture medium + 10% dimethyl sulfoxide (DMSO; Merck C6164).
When required, frozen cells were thawed quickly by incubating them in a water bath at 37 °C. Thawed cells were mixed with 10 mL culture medium, mixed and centrifuged for 5 min at 300 × g. The supernatant was discarded and the pellet was resuspended in 7 mL culture medium. Resuspended cells were transferred into a cell culture flask and incubated at 37 °C and 5% CO2.
Generation of theReporter Per2:luc
Circadian clock gene activity was recorded with a Bmal1:luc reporter (Addgene 68833, a gift from Steven Brown) and a Per2:luc reporter (Addgene 212035) contained in lentiviral transfer vectors (referred to as pLV6-Bmal-luc vector and pLV6-Per2-luc vector, respectively). The Bmal1:luc reporter contains a murine Bmal1 promoter with an adjacent luciferase sequence. The Per2:luc reporter was constructed in-house using a restriction enzyme strategy, which is efficient in constructing large plasmids (Zhang and Tandon, 2012). The Per2:luc reporter vector (pLV6-Per2-luc) consists of the backbone of the Bmal1:luc reporter vector (pLV6-Bmal-luc) and a murine Per2 promoter with adjacent luciferase sequence which is contained in the pGL3 basic E2 vector (Addgene 48747, a gift from Joseph Takahashi). The QuikChange Lightning Site-Directed Mutagenesis (SDM) kit (Agilent 210518) following the manufacturer’s manual was used to ligate the Per2:luc reporter with pLV6 backbone as follows.
The Bmal1:luc reporter vector contains one BamHI cutting site 3ʹ of the Bmal1:luc insert. The SDM kit and custom-designed oligonucleotides (Suppl. Table S1) were used to create an additional BamHI cutting site 5ʹ of the insert by a single-nucleotide mutation. After the SDM kit procedure, transformation of the mutated plasmid was done with stable competent Escherichia coli (NEB C3040H) according to the manufacturer’s manual. Cells were grown on LB agar plates (40 g/L LB agar, 100 µg/mL Ampicillin; Merck A9393) and single clones were cultivated in LB broth overnight (25 g LB, 100 µg/mL Ampicillin). DNA was extracted using the Qiagen miniprep (Qiagen 12123) and the mutation was checked by whole plasmid sequencing (Plasmidsaurus). After confirmation, the Bmal1:luc reporter vector was digested with BamHI (NEB R3136S), supplemented with calf intestinal alkaline phosphatase (Quick CIP, NEB M0525S) to avoid self-ligation. The digest was run on an 1% agarose gel (containing ethidium bromide) and the pLV6 backbone was extracted with the QIAquick Gel extraction kit (Qiagen 28704) according to the manufacturer’s manual.
Conveniently, the murine Per2 promoter and luciferase sequence in the pGL3-Per2-luc vector are by default flanked by BamHI cutting sites. Consequently, the pGL3-Per2-luc vector was directly digested with BamHI and the Per2:luc sequence was extracted from an agarose gel as explained above.
After successful digestions and extractions, the pLV6 backbone of the Bmal1:luc reporter vector and the Per2:luc fragment were ligated with T4 DNA ligase (NEB M0202S) according to the manufacturer’s manual. Subsequently, extracted vectors from transfected single colonies were sequenced by Plasmidsaurus and checked for correct ligation and correct sequence.
Generation of theWith Mutated RORE Sites Bmal1:luc
The two RORE elements (designated RORE1 and RORE2) of the Bmal1:luc reporter were mutated in series using the SDM kit following the manufacturer’s instructions. RORE-specific primers (Suppl. Table S1) were designed to mutate RORE sequences as previously reported (Guillaumond et al., 2005). Amplification of plasmids was achieved as described above. Confirmation of RORE mutation was confirmed by whole plasmid sequencing (Plasmidsaurus). The resulting reporter is subsequently referred to as Bmal1-mutated:luc and the whole vector is referred to as pLV6-Bmal1-mut-luc.
Lentiviral Vector Production and Transduction
The lentiviral vectors containing the respective reporters (Bmal1:luc, Per2:luc, and Bmal1-mutated:luc) were packed following broadly Du and Brown (2021) and Salmon and Trono (2007).
The day before transfection, HEK293T/17 cells (ATCC CRL-11268) were seeded at a density of 2 × 106/10 cm dish with the aim of a confluency of 40% to 50%. Packaging plasmids pMD2G and psPAX2 (Addgene 12259 and 12260, gifts from Didier Trono) and the respective transfer vector (pLV6-Bmal-luc, pLV6-Per2-luc, or pLV6-Bmal-mut-luc) were dissolved in TE buffer, pH 8.0 (Qiagen 12162) to a total volume of 250 µL. Next, 500 µL HEPES-buffered saline solution (Alfa Aesar AAJ62623AK) was added; 250 µL of 0.5M CaCl2 was added in an empty sterile 15 mL tube. The DNA/TE/HeB mix was added dropwise into the CaCl2-containing tube under constant vortexing. The mixture was then incubated for 20 to 30 min at room temperature during which a fine white translucent precipitate formed. One milliliter of the precipitate was added dropwise to the fibroblast culture and gently mixed. The next morning, the old medium was removed, washed twice with PBS, and 10 mL of fresh culture medium was added. For the next 2 days, 10 mL of culture medium was collected and replenished afterward. The pooled collection medium from two consecutive days (in total 20 mL) was centrifuged for 5 min at 500 × g and 4 °C. Thereafter, the supernatant was filtered through 0.45-μm syringe filters and then stored as 1 mL aliquots at −80 °C until further use.
Cultivated mouse and reindeer fibroblast were transduced with lentivirus containing either of the transfer vectors. The day before infection, confluent fibroblasts were split to achieve 50% confluency on the following day. The medium was replaced with warmed lentivirus-containing medium supplemented with 8 µg/mL protamine sulfate. After 48 h, the lentivirus-containing medium was renewed, and after 3 days, the medium was replaced with a medium containing blasticidin selection antibiotic (Merck 15205). The cells were then grown in T75 culture flasks before seeding into 35-mm dishes and were prepared for recording.
Two to three days before the recording, the cultures were transferred into a recording medium containing 5% FBS and 0.1 mM luciferin (Promega E1601) (Feeney et al., 2016; Yamazaki and Takahashi, 2005). Synchronization of the cultures was achieved by dexamethasone (DEX) (Merck D4902). Recording media in the culture dishes was exchanged for synchronization medium containing DEX (Recording media + 100 nM DEX) and incubated for ca. 30 min in the culture incubator (37 °C, 5% CO2). After incubation, the cultures were washed with PBS twice and a final volume of 2 to 3 mL recording media was filled in the 35-mm culture dishes. The dishes were sealed with parafilm and placed into a photomultiplier tube (PMT) (Hamamatsu Photonics LM-2400) placed inside a cell incubator. The PMT was connected to a PC with data acquisition software.
Cell Culture Experiments (Including Data Analysis)
Experiment 1: Different Temperatures, Luminescence Recording
Mouse and reindeer fibroblast cultures were transduced with either the clock gene promoter reporter Bmal1:luc (pLV6-Bmal1-luc) or Per2:luc (pLV6-Per2-luc). Fibroblast cultures (at passage number 6 for both species) were first measured under different temperatures in the following order 40 °C, 34 °C, 31 °C, and 37 °C. Luminescence was measured for 4 to 6 days under each temperature. Between each temperature recording, all cultures were re-synchronized as outlined above. Six replicates were used for each species-reporter combination, accumulating to a total of 24 cultures. However, two replicates of the reindeer-Per2 culture died early in the experiment resulting in four replicates for this group.
For each replicate, data for ca. 4 h following the DEX synchronization were omitted from the analysis and moving averages were calculated for 24 h around each respective data point, that is, 12 h before and after. Hence, a total of 16 h of data is not displayed in the graphs (4 h omitted + 12 h to calculate the first running mean). Phases were assessed by calculating the center of gravity of the bioluminescent peak on the first day 24 h after DEX synchronization. The center of gravity was generated by the CircWave program (version 1.4, programmed by Dr Roelof A. Hut). Phases were then normalized against the free-running period of the respective species and plotted in a circular graph in R studio (package: ggplot2). The displayed values represent means from all replicates recorded at 37 °C.
Periods, half-life times, and Q10 values were calculated based on damped sine wave analysis fitted with GraphPad prism (version 9.4.0) to Bmal1:luc recordings of mouse and reindeer cultures. Periods and half-life times were extracted from the damped sine wave analysis for each temperature. Half-life time is defined as the amount of time in hours for the amplitude of the rhythm to decay by half. Periods under different temperatures were used to calculate the Q10 value in R studio (R package: respirometry). The package calculates the Q10 based on the slope value of a fitted exponential model. Figures were produced in GraphPad prism. Statistical tests were performed in GraphPad prism.
Experiment 2: Temperature Entrainment, Luminescence Recording
In the second experiment, the ability of temperature entrainment for each species and each gene was measured. Both species were at passage number 8 at the time of recording. For each group, three replicates were used and transferred into the PMT under constant 37.2 °C. Approximately 12 h after the first synchronization, another three replicates of each group were synchronized. This resulted in half of each group being in antiphase to the other half. The cultures were recorded for 2 days under constant 37.2 °C. Thereafter, the cultures were transferred into a temperature cycle for 4 days with 12 h under 40.2 °C and 12 h under 37.2 °C. After those 4 days, the cultures were again transferred to 37.2 °C and measured for 2 days. Temperature was simultaneously recorded by a temperature logger (iButton, Maxim Integrated DS2922L) placed in the PMT and temperature was extracted with the program OneWire Viewer (Maxim Integrated Version 0.3.19.47).
Moving averages were calculated as described above. Centers of gravity were determined for each day of the recording in mouse and reindeer fibroblast and plotted as zeitgeber time (ZT = 0 refers to the first uprise in temperature). The values of the last day of the recording were compared by an unpaired t test. All graphs were produced in GraphPad prism.
We note that the PMT is sensitive to temperature; however, the amplitude of the changes in baseline within our temperature range is minimal and therefore does not confound our data (Suppl. Fig. S1).
Experiment 3: Serial Collection of Mouse and Reindeer Skin Fibroblasts for Quantitative Polymerase Chain Reaction Analysis
Mouse and reindeer skin fibroblasts were seeded into six-well plates and grown to confluency at 37.2 °C in DMEM + 10% FBS + P/S. Cultures (at passage number 9 for both species) were then synchronized using 100 nM DEX for 30 min and returned to the incubator for 48 h. After 50 h, we increased the temperature of the incubator to 40.2 °C for 12 h before returning the temperature to 37 °C until the end of the experiment. Plates were sampled starting at 48 h after synchronization and then every 4 h for the next 24 h. During sampling, the plates were rapidly retrieved from the incubator, the media was then removed before the dishes were washed 3× in ice-cold PBS. Finally, all excess PBS was removed before the plates were snap-frozen on dry ice and transferred to −80 °C freezer for storage.
RNA was extracted from each plate using RNeasy mini kit (QIAgen) and QIAshredder following the manufacturer’s instructions. We next converted 2 μg of RNA to cDNA using a high-capacity RNA to cDNA kit (Thermo Fisher) and performed quantitative polymerase chain reaction (qPCR) using GoTaq qPCR reagent. Gene and species-specific qPCR primers (Suppl. Table S1) were designed from genomic sequences (mouse, GRCm39; reindeer, GCA_949782905.1). We confirmed primer specificity by cloning and sequencing of PCR product and ensured primer efficiency (90%-105%) by standard dilutions. The housekeeping gene Ppib was used for data normalization.
Experiment 4: Temperature Entrainment of RORE-Mutated Bmal1 Reporter
In this experiment, U2OS cells, mouse fibroblasts and reindeer fibroblasts were transduced with either the normal Bmal1:luc (pLV6-Bmal1-luc) or the Bmal1-mutated:luc (pLV6-Bmal1-mut-luc). First, those cultures were recorded under constant conditions for 2 days. Thereafter, a three-stage protocol of temperature entrainment similar to experiment 2 was applied, with recordings under constant ambient temperature, followed by 4 days of a temperature cycle, followed by a release back into constant temperature. Moving averages were calculated as described above. Graphs were plotted in GraphPad.
Comparison of Reindeer and MousePromoter Sequences Bmal1
Genomic comparison was conducted by aligning the Bmal1 promoter sequence of the Bmal1:luc reporter (pLV6-Bmal1-luc) to the Bmal1 promoter sequence in reindeer (GCA_004026565.1) with Clustal Omega (https://www.ebi.ac.uk/↗). Furthermore, amino acid sequences of RORs and REV-ERBs between mouse and reindeer were compared by aligning them with Clustal Omega.
Results
Clock Gene Promoter Reporters Are Rhythmic and Temperature-compensated in Reindeer Skin Fibroblasts
We established primary cultures of mouse and reindeer skin fibroblasts and transduced them with either Bmal1:luc or Per2:luc clock gene promoter reporters. We then synchronized the cells with DEX and monitored bioluminescence of the reporter constructs for several days. In contrast to earlier work (Lu et al., 2010), both mouse and reindeer skin fibroblasts showed strong evidence of circadian oscillations, as shown by the reporter activity of Bmal1:luc and Per2:luc (Figure 1a). The free-running period (τ) of the Bmal1:luc was 25 h 20 min in mouse and 21 h 19 min in reindeer. Both species, however, had near-identical phase angles (Ψ) between the Per2 and Bmal1 reporters (Figure 1a, right panels).
Besides sustained rhythmicity under constant conditions, circadian oscillators are defined by their stable period length at different temperatures, a quality known as temperature compensation (Sweeney and Hastings, 1960). To test temperature compensation, and hence the robustness, of mouse and reindeer skin fibroblasts, we measured the oscillatory periods of cells transduced with the Bmal1:luc promoter reporter in cells incubated at 31 °C, 34 °C, 37 °C, or 40 °C (Suppl. Fig. S2). Although incubation temperature changed the dampening half-lives of the luciferase reporter (one-way analysis of variance [ANOVA], Fmouse = 97.73, Freindeer = 126.9, p for both <0.0001) (Figure 1b), both mouse and reindeer skin fibroblast cells remained rhythmic at all temperatures and maintained stable period lengths (Figure 1c). Q10 values were similar in mouse (0.916) and reindeer cells (0.907) (t test, p = 0.108, t = 1.766, df = 10), and well within the definitive 0.8 to 1.4 range of temperature compensation (Figure 1d) (Sweeney and Hastings, 1960).
Endogenous and temperature-compensated clock gene expression in mouse and reindeer fibroblast culture. (a) Representative luminescence data of mouse and reindeer fibroblasts under 37 °C ambient temperature. Expression of thereporter (black) is in approximate antiphase to the expressionreporter (gray) in both mouse and reindeer fibroblast cultures. Bioluminescence was measured in counts per minute (cpm). Phase angle difference (Ψ) betweenandfor the first day after DEX synchronization (from hour 24 to hour 48) are plotted as circular plots. Here, the:luc phase is calculated relative to the:luc phase which was set as 0° and expressed as a proportion of τ. Data are normalized for the respective free-running period (τ). (b) Half-life times ofexpression under different ambient temperatures. Half-life times were calculated by damped sine wave analysis and are defined as the amount of time for the amplitude to dampen by half. (c) Periods ofexpression under different ambient temperatures. Periods have been determined by damped sine wave analysis. (d) Q10 values ofin mouse and reindeer fibroblast cultures. Q10 values are based on periods as shown in panel c. Both Q10 values are within the acceptable range for temperature compensation (0.8-1.4) and do not differ significantly from each other (ns). Bmal1:luc of Per2:luc Bmal1:luc Per2:luc Bmal1 Per2 Bmal1:luc Bmal1:luc Bmal1:luc
Shows an Atypical Phase Relationship to Ambient Temperature Cycles in Reindeer Cells Bmal1
Daily fluctuations in body temperature are thought to be a universal entrainment cue in mammals (Buhr et al., 2010). To test temperature synchronization in mouse and reindeer skin fibroblasts, we ran a three-stage experiment. In stage 1, mouse and reindeer skin fibroblasts, transduced with either Bmal1:luc or Per2:luc reporters, were DEX-synchronized at two time points 12 h apart. Bioluminescence of those cultures was recorded under a constant temperature until stage 2. In stage 2, we exposed the cultures to a fluctuating temperature cycle with an amplitude of 3 °C for four 24 h cycles. In stage 3, the cultures were released back to constant temperature. The prediction for a temperature-entrainable circadian clock is that in stage 1, the cultures show two distinct phases 12 h apart, stage 2 will entrain cultures into the same phase, and stage 3 will reveal that this new phase persists in constant conditions.
In mouse fibroblasts, the Per2:luc reporter (Figure 2a) responded immediately to the temperature cycles (Figure 2c and Suppl. Fig. S3b) and the phase of the 12 h-apart synchronized cultures remained in phase to each other once back on constant temperature (t test, p = 0.0856, t = 2.271, df = 4). For the Bmal1:luc reporter (Figure 2b), the temperature cycle revealed a gradual entrainment (Figure 2c and Suppl. Fig. S3b) which emerged over several days. Release of the cultures from a temperature cycle to constant temperature showed near-synchrony of cultures whose acrophases (i.e., center of gravity) were 44 min ± 13 min apart from each other (Figure 2c) (t test, p = 0.0288, t = 3.341, df = 4). These data are consistent with previous temperature entrainment experiments in mouse cells (Saini et al., 2012).
In reindeer fibroblasts, the Per2:luc reporter activity under a temperature cycle (Figure 2d) was very similar to that seen in mouse skin fibroblasts, but with a remarkable 9 to 12 h phase jump when cells were released back into constant conditions. This surprising behavior was paralleled by unusual patterns of Bmal1:luc reporter activity (Figure 2e). Here, unlike the case in mouse fibroblasts, we observed an immediate alignment to the temperature cycle (Figure 2f and Suppl. Fig. S3b), with suppression by rising temperatures similar to the Per2:luc construct. When these cultures were released back into constant temperature conditions, Bmal1:luc and Per2:luc returned to an approximately antiphasic relationship with one another within the first cycle of free-running. Unlike the case for mouse fibroblasts, cultures that start the experiment in antiphase to one another still had quite distinctive phases (Bmal1:luc Ψ = 3 h 38 min ± 20 min; t test, p = 0.0004, t = 11.18, df = 4; Per2:luc Ψ = 2 h 24 min ± 36 min; t test, p = 0.0168, t = 3.951, df = 4).
We wondered if the observed acute sensitivity of Bmal1:luc to temperature was an artifact of the reporter construct. To test this, we cultured untransduced mouse and reindeer skin fibroblasts in six-well plates before synchronizing them with DEX and making serial collections over the first temperature cycle at a 4 h sampling resolution. We then measured the abundance of mouse and reindeer Per2 and Bmal1 mRNA during the time course using qPCR. In mouse skin fibroblasts, mouse Per2 and Bmal1 mRNA levels peaked in distinct phases during the temperature cycle (Figure 2g), whereas reindeer Per2 and Bmal1 mRNA levels were overlapping in phase (Figure 2h). These data are consistent with the reporter activity and provide independent evidence that, in contrast with mouse skin fibroblasts, temperature has a strong and immediate impact on both Bmal1 and Per2 gene expression in reindeer skin fibroblasts.
Temperature entrainment of clock genes expression in mouse and reindeer fibroblast cultures. Expression ofandin mouse (a-c) and reindeer fibroblasts (d-f) under constant and cycling ambient temperature (Ta). Two sets of cultures were DEX-synchronized 12 h apart from each other. Bioluminescence was measured in counts per minute (cpm). For reindeer, bioluminescence data after the temperature cycle were plotted in zoomed-in plots besides the corresponding graphs. Times of peak expressions (center of gravity) were determined with the CircWave program for each rhythm to track phase relationships and responses before, during, and after the temperature cycle (c, f). Phase relationships between the initially 12 h-apart synchronized cell cultures of the same gene were calculated after the temperature cycle and tested with an unpairedtest (asterisk indicates< 0.05, ns denotes> 0.05). (g, h) Relative gene expressions ofandin mouse and reindeer fibroblasts under one temperature cycle. Gene expression data were measured by qPCR with species-specific primers. Per2:luc Bmal1:luc t p p Bmal1 Per2
Mutation of RORE Elements in aPromoter Reporter Produces Reindeer-Like Reporter Rhythms Under Ambient Temperature Cycles in a Mouse Cell Context Bmal1
The REV-ERB and ROR families of nuclear hormone receptors are principal regulators of circadian Bmal1 expression. These factors activate (RORs) or repress (REV-ERBs) transcription through two conserved RORE sites located in the Bmal1 proximal promoter (Abe et al., 2022; Akashi and Takumi, 2005; Guillaumond et al., 2005; Preitner et al., 2002; Ueda et al., 2002), and, in the mouse, this stabilizes circadian rhythm expression (Abe et al., 2022).
We showed that Bmal1:luc in reindeer fibroblasts and reindeer Bmal1 mRNA abundance is highly temperature-sensitive (Figure 2). The proximal promoter of reindeer Bmal1 contains both ROREs deemed necessary for circadian transcription (Suppl. Fig. S5), and reindeer RORs and REV-ERBs show strong sequence homology across the functional domains required for their transcriptional activator/repressor activity (Suppl. Fig. S6). These findings suggest that the effect of temperature on Bmal1 promoter activity may be independent of the effect mediated via the circadian clock and its known early temperature-responding elements such as Per2 (Buhr et al., 2010; Saini et al., 2012; Tamaru et al., 2011).
To test this hypothesis, we mutated both RORE elements in the Bmal1:luc reporter (creating a new reporter: Bmal1-mutated:luc) (Figure 3a). To confirm that the circadian clock machinery was decoupled from our new Bmal1-mutated:luc reporter, we compared the reporter activity of the Bmal1:luc and Bmal1-mutated:luc reporters following DEX synchronization at a constant ambient temperature. As expected, transduction of U2OS cells (a human-derived stable cell line and common in vitro molecular clock model), mouse skin fibroblasts and reindeer skin fibroblasts with the Bmal1-mutated:luc reporter showed a profound loss of rhythmic reporter activity compared with cells transduced with the wild-type Bmal1 reporter (Suppl. Fig. S4a-c).
We next ran another temperature entrainment experiment but this time comparing the Bmal1:luc reporter with the Bmal1-mutated:luc reporter in mouse, reindeer, and U2OS cells. The responses of the wild-type Bmal1:luc reporter in mouse and reindeer skin fibroblasts (Figure 3b) were consistent with our previous experiment (Figure 2). Interestingly, whereas the Bmal1-mutated:luc reporter showed no cyclicity in mouse, reindeer, or U2OS cells following DEX synchronization, it showed a strong thermal response with rising temperatures causing a dramatic suppression of bioluminescence and falling temperatures increasing bioluminescence (Figure 3b and Suppl. Fig. S4d). Moreover, the relative phases of the bioluminescence profiles relative to the temperature cycle for the Bmal1-mutated:luc reporter were strikingly similar to those seen for the wild-type Bmal1:luc reporter in reindeer fibroblasts (Figure 3c). This result indicates that the proximal Bmal1 promoter contains elements that confer a high degree of direct sensitivity to temperature, independent of circadian coupled ROREs.
RORE mutation of thereporter and responses under temperature cycles. (a) RORE sites in thepromoter facilitate the response of REV-VERBs and RORs onexpression. Both RORE sites in the promoter reporter sequence were rendered non-functional by mutation (WT: un-mutated). (b) The mutated and un-mutated reporters were transduced into mouse and reindeer fibroblasts and were measured under an ambient temperature (Ta) cycle for 4 days. (c) Peak expression was determined by the CircWave program and plotted for both species and both promoter reporters. Bmal1:luc Bmal1 Bmal1
Discussion
Taken together, our data show that reindeer skin fibroblasts have a circadian rhythm that can be synchronized with DEX, are temperature-compensated, and undergo partial entrainment following 24 h temperature cycles. We show that reindeer fibroblasts have a comparatively short free-running period compared with mouse cells and that Bmal1 has heightened direct sensitivity to temperature in reindeer compared with mouse fibroblasts. This alteration of the relationship between temperature and circadian regulation of Bmal1 might be a consequence of weaker coupling between the core and secondary loops in the reindeer molecular clock, leading direct temperature effects on Bmal1 transcription to dominate over circadian transcription.
Previous work by Lu et al. (2010) reported that the reindeer molecular clock was non-functional. The striking differences between our results and those presented by Lu et al. (2010) were unexpected, especially as we ensured minimal differences between our experimental approaches. The reindeer skin biopsies were both collected from animals that originated from a herd held at the University of Tromsø, Norway (69°N) and the culture conditions and synchronization protocols between the studies were identical. The Per2:luc and Bmal1:luc reporters used in our study are of different origins, but ultimately, both sets of promoter reporters contain the key regulatory elements that allow mouse skin fibroblasts to drive rhythmic, phase-specific transcription. In contrast to the original study, we included an additional blasticidin selection step to ensure that all fibroblasts expressed the reporter constructs. However, it is unlikely that the previously reported arrhythmicity of reindeer fibroblasts is caused by a lower lentiviral transduction efficiency. Under such a scenario, weak rhythmicity would also be expected in the positive control, that is, mouse fibroblasts. We also ruled out if variable cell density could account for the differences between the two studies. Previous work has shown that low-density skin fibroblast cultures have low rhythmicity (Noguchi et al., 2013). In our own work, we ensured that the cells were confluent before synchronization. In the article by Lu et al. (2010), the confluency of the cells was not stated; however, the baseline of their raw data did not increase over time as would be expected in a low-density culture that it continues to proliferate over the course of the experiment. It is, therefore, unlikely that cell density accounts for the differences between our studies.
Although reindeer skin fibroblasts satisfy key criteria of circadian rhythms, we identify several distinctions between our mouse and reindeer model systems. Most notable were the data collected during our daily temperature cycle experiments. Here, the expressions of both Bmal1 and Per2 were temperature-sensitive in reindeer skin fibroblasts, leading to an unusual in-phase expression profile between the genes. We also observed a low amplitude and comparatively weak circadian alignment of the reindeer cultures once they were released back to constant temperature. This likely reflects the re-establishment of the cell-autonomous circadian oscillators within the reindeer skin fibroblast population. Here, the interaction of the positive (Bmal1) and negative (Per2) regulators of the core molecular clock resolves from their atypical in-phase state to their typical antiphase relationship evident as a large phase shift in Per2. The dampened amplitude might reflect phase variation between individual cells in the newly established circadian phase measured at the level of the whole cell population in a dish (Welsh et al., 2004). It would be useful to develop these observations in future studies through serial collection and measurement of Bmal1 and Per2 transcripts under constant conditions both before and following temperature cycles.
Intrigued by the unusual regulation of Bmal1 by daily temperature cycles in reindeer skin fibroblasts, we turned our attention to the molecular processes that underpinned the striking differences in daily temperature cycles between mouse and reindeer skin fibroblasts. To investigate the role that RORE-mediated transcription plays in temperature resetting, we experimentally uncoupled the circadian clock machinery from the Bmal1:luc reporter by mutating the two canonical RORE elements. Our experiments using this construct revealed a strikingly similar response to temperature cycles in mouse and U2OS cells as in reindeer with the wild-type reporter, namely an immediate phase shift with a trough at temperature rise (Figure 3).
These data show that RORE mutation, and therefore removal of the ROR/REV-ERB interaction at the Bmal1 promoter, leads to clock-independent temperature sensitivity of the reporter. It is clear, however, that endogenous circadian ROR/REV-ERB expression does regulate the Bmal1:luc reporter in reindeer skin fibroblasts as there is a robust circadian rhythm under constant conditions in reindeer skin fibroblasts which is disrupted by RORE mutation.
This leads us to suggest a model in which a species difference in the relative weighting between direct RORE-independent effects of temperature and indirect RORE-mediated circadian effects of temperature accounts for the different responses seen in temperature entrainment experiments (Figure 4). We believe that the most parsimonious interpretation of our data involves weakened coupling between the stabilizing (RORE-ROR/REV-ERB) and core (E-Box-BMAL1/PER) circadian clock loops in reindeer compared with mouse fibroblasts. Under this model, in mouse fibroblasts, the ROR/REV-ERB transcription factors strongly couple Bmal1 expression to the circadian clock, and it is this, probably via direct effects of temperature on Per2 expression (Buhr et al., 2010; Tamaru et al., 2011), which mediates the entrainment response to temperature cycles. In contrast, in reindeer fibroblasts, circadian clock–mediated effects of temperature on the ROR/REV-ERB transcription factors are seen as exerting a relatively weak regulatory effect on Bmal1 transcription, while direct, RORE-independent effects of temperature dominate. The consequence of this is that under a temperature cycling paradigm Bmal1 phasing immediately tracks temperature and the antiphasic relationship between Per2 and Bmal1, which is a hallmark for circadian coupling, breaks down. Interestingly, mathematical modeling of the circadian clock predicts that a weakened influence of the secondary loop on the core loop leads to shortened free-running periods (Yan et al., 2014), a prediction that is consistent with the comparatively short free-running period of the reindeer compared with mouse skin fibroblasts.
Our experiments leave open the question of how coupling between the reindeer core and secondary loops has become weaker than is the case in mice. We could find no obvious differences either in proximal Bmal1 promoter organization or in the major DNA binding, ligand binding, and protein-protein interaction domains in the REV-ERBs and RORs. It is possible that a recently described reindeer-specific mutation in the PER2 protein, which affects its interaction with CRY1 (Lin et al., 2019), is involved in the phenomena we describe, and genetic manipulation of this site is a desirable future objective.
The whole-organism implications of a reindeer-type compared with a mouse-type circadian oscillator are difficult to predict. However, given that thermoenergetic demands impact the circadian organization of 20 g nocturnal rodents radically differently compared with 100 kg diurnal ungulates of the Arctic (Hazlerigg and Tyler, 2019; van der Vinne et al., 2014), we speculate that our results reflect unappreciated diversity in the coupling of heat-sensitive signal transduction to circadian entrainment and metabolism.
In mouse fibroblasts, thereporter entrains gradually to an ambient temperature cycle eventually resuming its approximate antiphase relationship withexpression, suggesting that the-specific temperature effect is dominantly mediated via the circadian clock. In reindeer fibroblasts,is acutely responding to temperature cycles and moves into phase withexpression. We suggest that a weak coupling between circadian loops reduces the effect of temperature via the RORE elements on thepromoter; instead, an unknown RORE-independent pathway drivesexpression during temperature cycles in reindeer fibroblasts. Supporting this model, when we experimentally decoupled the circadian clock from thereporter by mutations of the ROREs, we rendered the promoter reporter acutely temperature-sensitive which led to a reindeer-like phase response under temperature cycles. Bmal1:luc Per2 Bmal1 Bmal1 Per2 Bmal1 Bmal1 Bmal1:luc
Supplemental Material
Acknowledgments
The authors would like to thank our animal technicians Hans Lian, Hans-Arne Solvang, and Renate Thorvaldsen. Only their dedication and hard work allows us studies on our fascinating Arctic model organisms. We also thank Jaione Simón-Santamaría and Karolina Szafranska for the mouse skin biopsies and Barbara Tomotani for useful discussions on the manuscript. We also thank Fredrik Markussen for his circular-plotting skills. The work was supported by grants from the Tromsø forskningsstiftelse (TFS) starter grant (TFS2016SW), the TFS infrastructure grant (IS3_17_SW) awarded to S.H.W. and Nansenfondet (ReinRytmer) awarded to A.C.W. The Arctic seasonal timekeeping initiative (ASTI) grant and UiT strategic funds support D.G.H., S.H.W., A.C.W., and D.A.