Abstract
Circadian rhythms of mammalian physiology and behavior are coordinated by the suprachiasmatic nucleus (SCN) in the hypothalamus. Within SCN neurons, various aspects of cell physiology exhibit circadian oscillations, including circadian clock gene expression, levels of intracellular Ca2+ ([Ca2+]i), and neuronal firing rate. [Ca2+]i oscillates in SCN neurons even in the absence of neuronal firing. To determine the causal relationship between circadian clock gene expression and [Ca2+]i rhythms in the SCN, as well as the SCN neuronal network dependence of [Ca2+]i rhythms, we introduced GCaMP3, a genetically encoded fluorescent Ca2+ indicator, into SCN neurons from PER2::LUC knock-in reporter mice. Then, PER2 and [Ca2+]i were imaged in SCN dispersed and organotypic slice cultures. In dispersed cells, PER2 and [Ca2+]i both exhibited cell autonomous circadian rhythms, but [Ca2+]i rhythms were typically weaker than PER2 rhythms. This result matches the predictions of a detailed mathematical model in which clock gene rhythms drive [Ca2+]i rhythms. As predicted by the model, PER2 and [Ca2+]i rhythms were both stronger in SCN slices than in dispersed cells and were weakened by blocking neuronal firing in slices but not in dispersed cells. The phase relationship between [Ca2+]i and PER2 rhythms was more variable in cells within slices than in dispersed cells. Both PER2 and [Ca2+]i rhythms were abolished in SCN cells deficient in the essential clock gene Bmal1. These results suggest that the circadian rhythm of [Ca2+]i in SCN neurons is cell autonomous and dependent on clock gene rhythms, but reinforced and modulated by a synchronized SCN neuronal network.
Significance Statement
Intracellular calcium ([Ca2+]i) oscillates with a circadian (∼24-h) rhythm within the SCN, the brain’s master circadian pacemaker. However, what drives these [Ca2+]i rhythms is not well understood. To test how rhythmic circadian clock gene expression and the SCN neuronal network influence [Ca2+]i rhythms, we monitored [Ca2+]i and PER2 (a rhythmic clock gene) in individual SCN neurons in dispersed and slice cultures by fluorescence and bioluminescence imaging, respectively. We also tested the effects of blocking neuronal firing or knocking out Bmal1, an essential clock gene. We found that the [Ca2+]i rhythm requires intact clock gene rhythms, and that it does not require the SCN neuronal network (so is cell autonomous) but is reinforced and modulated by the neuronal network.
Introduction
In mammals, circadian rhythmicity of behavior and physiology is controlled by a master circadian pacemaker in the hypothalamus: the suprachiasmatic nucleus (SCN). Circadian rhythm generation within SCN neurons and other cells depends on a core transcription-translation negative feedback loop (TTFL). This core loop involves activation of Period (Per1, Per2, Per3) and Cryptochrome (Cry1 and Cry2) gene transcription by a BMAL1/CLOCK heterodimer, and delayed inhibition of this process by complexes containing PER and CRY proteins (Mohawk and Takahashi, 2011). Null mutation of Bmal1 alone is sufficient to abolish circadian rhythms of behavior (Bunger et al., 2000) or single SCN neurons (Ko et al., 2010).
In SCN neurons, various cellular processes exhibit circadian rhythms, including clock gene expression, Ca2+, neuronal firing rate, and neuropeptide release (Welsh et al., 2010). SCN neurons communicate through synapses (Yamaguchi et al., 2003), diffusible messengers (Silver et al., 1996; Maywood et al., 2011), and possibly gap junctions (Colwell, 2000b) to create coherent rhythms. Although individual SCN neurons can function as independent circadian oscillators (Welsh et al., 1995), the SCN network contributes to the strength of cellular rhythmicity (Webb et al., 2009).
Ca2+ plays important roles in both generation of circadian rhythms in SCN neurons and their synchronization by retinal input (Colwell, 2011). Previous studies found adequate extracellular and intracellular Ca2+ levels to be essential for generation and regulation of cellular circadian rhythms (Lundkvist et al., 2005; Harrisingh et al., 2007). In retinorecipient SCN neurons, [Ca2+]i is also an important mediator in the signaling pathway from postsynaptic glutamate receptors to clock gene expression (Colwell, 2011).
Circadian oscillation of [Ca2+]i has been demonstrated in plants (Johnson et al., 1995), Drosophila (Liang et al., 2016), and rodent SCN (Colwell, 2000a; Ikeda et al., 2003; Enoki et al., 2012; Brancaccio et al., 2013). One early study using the Ca2+-sensitive dye fura2 (Colwell, 2000a) found a circadian [Ca2+]i rhythm that is dependent on both neuronal firing and voltage-gated Ca2+ channel activity. In contrast, studies using genetically encoded Ca2+ reporters have revealed [Ca2+]i rhythms that persist in the absence of neuronal firing (Enoki et al., 2012; Brancaccio et al., 2013) but are suppressed by blockers of ryanodine-sensitive intracellular [Ca2+] stores (Ikeda et al., 2003). Although the relationship between clock gene and [Ca2+]i rhythms has not been studied at a single-cell level in dispersed culture, [Ca2+]i peak times precede PER2 peak times in SCN slices (Brancaccio et al., 2013; Enoki et al., 2017a) and show a spatiotemporal wave pattern similar to that previously reported for Per1 and PER2 (Yamaguchi et al., 2003; Foley et al., 2011; Enoki et al., 2012). In Cry1/Cry2 double knockout mice, SCN [Ca2+]i and neuronal firing rhythms are abolished, especially in dorsal SCN (Enoki et al., 2017a). Phase relationships between [Ca2+]i and neuronal firing rhythms are different in dorsal and ventral mouse SCN (Enoki et al., 2017b). In Drosophila, the [Ca2+]i rhythm is attenuated in per01 null mutant flies, which lack molecular and behavioral rhythmicity, and its phase is modulated by neuropeptide signaling (Liang et al., 2016). These studies raise intriguing questions about how the SCN network regulates [Ca2+]i–PER2 phase relationships spatiotemporally, and how [Ca2+]i rhythms are related to neuronal firing, neuronal networks, and other important clock genes.
In this study, we evaluated three possible mechanisms that might contribute to SCN [Ca2+]i rhythms: (1) the TTFL molecular clock, (2) spontaneous action potentials within single SCN neurons, and (3) the synchronized neuronal network coupled by neurotransmitters. To monitor [Ca2+]i and clock gene expression rhythms simultaneously in individual SCN neurons, the genetically encoded Ca2+ reporter GCaMP3 was introduced by adeno-associated virus (AAV) to SCN slices or dispersed cells cultured from PER2::LUC knock-in reporter mice, which express a fusion of PER2 and firefly luciferase in all cells (Welsh et al., 2004). Furthermore, we used computational models and image analysis to reveal neuronal network contributions to the generation and spatiotemporal patterning of clock gene and [Ca2+]i rhythms.
Materials and Methods
Experimental design and statistical analysis
Animals
Generation of mPer2Luciferase (PER2::LUC) knock-in mice was described previously (Yoo et al., 2004). For this study, we used an alternative PER2::LUC mouse line incorporating an SV40 polyadenylation site to enhance expression levels (Welsh et al., 2004). The mice were developed at Northwestern University using the same methodology as the original strain of knock-in mice. PER2::LUC mice were backcrossed for >10 generations with wild-type (WT) mice (C57BL/6J) purchased from the Jackson Laboratory. Bmal1 –/– mice were donated by Dr. Pamela Mellon (University of California, San Diego). This line was created by crossing Bmal1-floxed mice (The Jackson Laboratory, stock number 7668; Storch et al., 2007) and zona pellucida 3 (Zp3)-cre mice (Lewandoski et al., 1997). PER2::LUC;Bmal1 –/– mice were created by crossing PER2::LUC mice and Bmal1 –/– mice. PER2::LUC mice were bred as homozygotes, and PER2::LUC;Bmal1 –/– mice were bred as heterozygotes at the Bmal1 allele (Bmal1 +/–). Genotyping was performed by multiplex PCR as described in a previous report (Storch et al., 2007).
Mice were maintained in LD 12:12-h cycles (12 h light, 12 h dark) throughout gestation and from birth until used for experiments. Mouse studies were conducted in accordance with regulations of the Institutional Animal Care and Use Committee at University of California, San Diego.
SCN slice and dispersed culture
Both male and female mice were used in this study. For SCN slice cultures, 4- to 7-d-old PER2::LUC knock-in mice or PER2::LUC;Bmal1 –/– mice were anesthetized by ketamine (Ketaset, Fort Dodge Animal Health) and decapitated. The hypothalamus block was dissected from the whole brain, and coronal slices were cut by a tissue chopper (Stoelting) to a thickness of 250 µm. Each SCN slice was cultured on a Millicell-CM membrane insert (EMD Millipore PICMORG50; Merck KGaA) in explant medium with low sodium bicarbonate (EML) [DMEM (Thermo Fisher Scientific, 12100046), supplemented with 350 mg/L sodium bicarbonate, 10 mM HEPES, 25 U/mL penicillin, 25 µg/mL streptomycin, 2% B-27 (Thermo Fisher Scientific, 17504-044), and 200 µM d-luciferin potassium salt (BioSynth L-8220; no phenol red, serum-free; pH 7.4].
Dispersed SCN neurons were obtained from neonatal (0- to 4-d-old) PER2::LUC knock-in mice. 1 × 105 cells were plated in a well of a glass bottom dish (MatTek P35G-1.5-10-C) and cultured as described previously (Welsh et al., 1995; Liu et al., 2007). Tetrodotoxin (TTX) and ryanodine were purchased from Sigma-Aldrich (T5651).
Viral vector transduction
AAV-2/1-hSynapsin-GCaMP3 was purchased from the vector core laboratory of the University of Pennsylvania. For dispersed cell culture, 3 µl virus solution (5 × 1011 genomic copies/0.1 mL) was mixed with dispersed cells (∼200 µl) during the first 2 h of plating. For slice culture, 0.5 µl virus solution was dropped onto an SCN slice immediately after the start of culture. To allow GCaMP3 expression levels to saturate, SCN dispersed cells and slices were cultured for at least 2 wks before imaging. Strong GCaMP3 expression was observed in most neurons.
Fluorescence and bioluminescence imaging
Simultaneous fluorescence (FL) and bioluminescence (BL) imaging was performed by adding FL imaging to previously described methods for BL imaging (Welsh et al., 2005; Welsh and Noguchi, 2012). The culture was placed on the stage of a microscope (Olympus IX-71). Light from the sample was collected with a 10× objective lens (10× UPlanSApo; Olympus) and transmitted to a cooled CCD camera (iKon-M; Andor) mounted on the bottom port of the microscope. The image was focused using GCaMP3 fluorescence, lights were turned off, and stray light was eliminated by covering the dish with a small black lucite box and draping the microscope with black plastic sheeting (Thorlabs BK5). GCaMP3 was excited by a 500-nm LED (pE-2; coolLED), and emitted light was filtered by a yellow fluorescent protein (YFP) filter cube (YFP-2427A; Semrock). LED light intensity was reduced to either 10% or 50% of maximal level, depending on GCaMP3 expression levels. FL images were collected at intervals of 30 min or 1 h with 1-s exposure duration, no binning, pre-amplifier gain ×4, and readout speed 5 MHz. To avoid autofluorescence of luciferin, excitation and filter settings were optimized for YFP rather than for green fluorescent protein (GFP), which has absorbance and emission spectra more similar to GCaMP3 (Tian et al., 2009). For BL image collection, the LED light was turned off, and the filter cube turret was moved to an open (no filter cube) position. BL images were collected at intervals of 30 min or 1 h, alternating with FL images, with 29-min exposure duration, binning 1 × 1 or 4 × 4, pre-amplifier gain ×4, and readout speed 50 kHz. To align with FL images, BL images were shifted 1 pixel in the x direction and –9 pixels in the y direction. MetaMorph (Molecular Devices) was used for control of camera and LED shutter and for image analysis. Time series images were typically collected for 5–7 d.
Image processing for manual cell tracking
In MetaMorph, stack files of time series of either FL or BL images were created. To remove cosmic ray artifacts, pairs of consecutive images were compared pixel-wise, and the minimum values of all pixels were used to construct a new image from every pair of consecutive images (running minimum algorithm; Welsh et al., 2004). The same algorithm was run on a FL stack file to make it comparable to a BL stack file. Then, a combined stack file was constructed by interleaving consecutive FL and BL image planes. In dispersed cell cultures, neurons were morphologically identified by their round cell bodies and long processes in FL images, and matching of neurons in alternating FL and BL images ensured identification of single cells. Bright PER2::LUC signals also confirmed that neurons were derived from the SCN. Both FL and BL intensity were measured within a region of interest (ROI) defined manually for each cell. For dispersed cultures, smaller (3 × 3-pixel) ROIs were used for FL images, selecting the brightest part of a cell, and larger (∼10 × 10-pixel) ROIs about the same size as a cell body were used for BL images. In slice cultures, identification of single cells was difficult, especially in BL images with lower resolution, because of high cell density and multiple layers of cells in a slice. For slice cultures, small (3 × 3-pixel) ROIs were used for both FL and BL, selecting the brightest part of a cell in FL images, and measuring the intensity of both FL and BL in the same position in a combined stack file. The position of the ROI was adjusted to accommodate movements of cells during the experiment, but its size was kept constant across the time series and across cells. Data were logged to Microsoft Excel files for plotting and further analysis.
BL intensity values were corrected for background level in each experiment by subtracting a single minimum intensity of a background region devoid of cells. Background-subtracted values were normalized to an average BL intensity (PER2 level) of neurons in the same dish. After subtracting 1 from each normalized value, the resulting value was plotted in relative units (i.e., 0 = average PER2 level of neurons in the dish; Figs. 1B, C, 2D, 3A, B, and 4A). See figures in the Results section.
Figure 1-1
Figure 1-2
Figure 3-1
FL intensity of background regions devoid of cells varied with position and time, reflecting variations in the LED light source, CCD camera operation, and a slow upward trend in GCaMP3 expression over time. Therefore, FL intensity of a cell was normalized using a set of multiple background regions devoid of cells, and then detrended by subtracting a 24-h running average. In each dissociated cell experiment, six different ROIs in areas devoid of cells were selected as backgrounds. For slice cultures, nonrhythmic tissue areas surrounding SCN were included as background ROIs. Time series were plotted of the intensities of these six backgrounds (BG1 to BG6, named from low to high intensity). At each time point, we assumed that intensities of these background regions were linearly related: where BGn are background intensity values, An are coefficients, and Bn are constants. An and Bn (n = 1, 2, 3, 4, 5) could then be obtained using BGn values at two time points (minimum and maximum values in the time series). Both An and Bn (n = 1, 2, 3, 4, 5) were linearly related to the average intensity of BGn. Using these linear relationships, Acell and Bcell could then be obtained for any cell-containing region with a particular average intensity, and the expected background intensity for each cell (BGcell) estimated by the following equation:
Thus, FL intensity of each cell at each time point was normalized to an expected background intensity by dividing by BGcell, and long-term trend was removed by subtracting a 24-h running average. For data in the first or last 12 h, the first or last 24-h average was subtracted, respectively. Processed FL intensity ([Ca2+]i) values are shown on the right axes of data plots (Figs. 1B, C, 2D, 3A, B, and 4A).
Rhythm data analysis for manually tracked data
Single-cell circadian rhythmicity was assessed in two ways. Autocorrelation was applied as a simple and robust metric of rhythmicity, using xcorr in Matlab (MathWorks). This method has the advantage of being independent of waveform, so it is equally applicable to PER2 rhythms (which resemble a sine wave) and [Ca2+]i rhythms (which do not). BL and FL data were processed as described above, and then the autocorrelation was calculated for 3- or 4-d segments of each time series after linear detrending. The rhythmicity index (RI) was defined as the value of autocorrelation at 24-h lag, which is applicable even for neurons without significant rhythmicity, such as Bmal1 –/– neurons. RI varies from 0 to 1, where 1 indicates a perfectly periodic time series. To assess for the presence of statistically significant rhythmicity, we used a test developed previously (Leise et al., 2012) based on regression of the log-log plot of the power spectral density and motivated by the analysis of arrhythmic Bmal1–/– SCN cells reported by Ko et al. (2010). Here, we call this test the PSD-regression test. Cells with p < 0.05 (PSD-regression test) were categorized as rhythmic. Circular statistics were run in Matlab using the CircStat package (Berens, 2009).
In Figs. 1, 2, and 4, 4D segments of time series were used to assess rhythmicity with both RI and the PSD-regression test. In Fig. 3, 3D segments before, during, and after TTX application were used. In control conditions, RI was affected by duration of data (p < 0.01, t test; Figs. 2E and 3C, E), so for TTX experiments, we ensured that all segments were of similar duration. For cross-correlations in TTX experiments (Fig. 3E, H), average PER2 intensity over the 2.0 d before TTX application (before), from 0.5 to 2.5 d after TTX application (TTX), and from 0.5 to 2.5 d after washout (after) were used.
For cells categorized as rhythmic, coherence between PER2 and [Ca2+]i rhythms was calculated as the peak of cross-correlation for linear-detrended time series using xcorr in Matlab. The time lag associated with this peak was used to estimate the phase difference between the PER2 and [Ca2+]i rhythms of each cell. To test for statistical differences of RI and coherence between SCN slice cultures and dispersed cell cultures, Fisher z-transformation followed by hierarchical (mixed effect) model (Bates et al., 2015) was applied using R (Bates et al., 2015; Brockhoff et al., 2016).
Spatiotemporal imaging data analysis
After removing cosmic rays by a running minimum algorithm in MetaMorph, FL and BL stack files were analyzed using ImageJ and Matlab scripts using methods described previously (Evans et al., 2011a). Briefly, ROIs of cell-like regions were identified as small bright disks in the 2D wavelet filtered sum of days 0.5–2.5 of the recording. Cell-like regions were judged to be significantly rhythmic if they passed the Siegel test at significance level 0.01. Peak times of PER2 and [Ca2+]i for each cell-like region were calculated after detrending using a discrete wavelet transform.
Mathematical models
The effect of TTX on the SCN was studied using previously developed mathematical models of the SCN (Vasalou and Henson, 2010, 2011; Kingsbury et al., 2016). All simulations were performed in Matlab. The model had four basic components: core molecular clock, electrophysiology, network topology, and neurotransmitter signaling. The previous model contained 21 ordinary differential equations (ODEs): the 16 ODEs of the Leloup–Goldbeter core molecular clock model (Leloup and Goldbeter, 2003), two ODEs governing both cytoplasmic ([Ca2+]i) and ryanodine storage calcium concentrations (Vasalou and Henson, 2010), one ODE for the phosphorylated cAMP response element-binding protein (CREB) concentration (To et al., 2007), and two ODEs for the neurotransmitters vasoactive intestinal peptide (VIP) and γ-aminobutyric acid (GABA), which controlled all cell-to-cell signaling. The relationship between [Ca2+]i and the kinetics of CREB phosphorylation, as well as activated CREB’s feedback on the molecular clock, have been modeled previously (To et al., 2007) based on principles from experimental observation (Obrietan et al., 1999; Tischkau et al., 2003).
To generate heterogeneity within the network population, we randomly distributed the parameters determining the baseline rates of Per mRNA transcription (vsP0; 0.94 ± 0.06), Bmal1 mRNA transcription (vsB; 1.00 ± 0.01), and Bmal1 mRNA degradation (vmB; 0.80 ± 0.01), as described previously (Vasalou et al., 2009), as well as a parameter scaling the entire ODE governing [Ca2+]i (distributed evenly between 0 and 1). All other parameter values were preserved from previous modeling studies (Vasalou and Henson, 2010; Vasalou et al., 2011; Kingsbury et al., 2016).
For simulations in which cells were coupled together, the GABA network topology was designed (Kingsbury et al., 2016) to conform to experimental results in which GABA connections were mapped using high-density multi-electrode arrays (Freeman et al., 2013). As described in Kingsbury et al. (2016), the VIP network was taken as a subset of the GABA network by assigning a random variable to each cell that removed outgoing VIP connections from 80% of the network. For dispersed cells, neurons did not receive any VIP or GABA signals, whereas in the coupled network, all neurons were forced to have at least one incoming VIP signal. No VIP or GABA signals were received during simulated TTX application; this was achieved by setting firing rates to zero and increasing the values of parameters governing GABA and VIP degradation (ndGABA and ndVIP, respectively) from 0.2 to 1.0.
Results
In dispersed SCN neurons, both PER2 and [Ca2+]i rhythms are cell autonomous, but the strength of [Ca2+]i rhythmicity is limited by PER2 rhythmicity
To study whether [Ca2+]i rhythms are cell autonomous or require a rhythmic SCN neuronal network, we measured [Ca2+]i in SCN neurons dissociated from PER2::LUC mice and maintained in dispersed cultures. Both PER2 and [Ca2+]i dynamics were monitored by imaging in single neurons (Fig. 1A). Cells were manually tracked in imaging data over several days (4.7–7.8 d). Rhythmicity of both PER2 expression and [Ca2+]i were analyzed for 87 cells in 3 cultures over a 4-d window (Fig. 1). Strength of rhythmicity (RI) was assessed by autocorrelation, and the presence of statistically significant rhythmicity was assessed by the PSD-regression test for both PER2 and [Ca2+]i, as described in Materials and Methods. There was no statistical correlation between GCaMP3 expression level and [Ca2+]i RI. Even cells with the lowest GCaMP3 expression still exhibited robust circadian rhythm, suggesting that weak expression of CCaMP3 was not responsible for lack of rhythmicity in other cells.
Most cells (80%) had robust circadian rhythms of PER2 expression, as reported previously (Liu et al., 2007; Webb et al., 2009; Ko et al., 2010; Fig. 1C, D). [Ca2+]i rhythms, like PER2 rhythms, were independently phased across cells, and therefore cell autonomous. However, only 40% of dispersed cells had significant [Ca2+]i rhythms; even in cells with rhythmic PER2, only 46% had significant [Ca2+]i rhythms (Fig. 1B, D). In contrast to the nearly sinusoidal PER2 waveforms, [Ca2+]i rhythm waveforms were more variable, often showing spike-like waveforms with narrow peaks (Fig. 1C). All cells with significant [Ca2+]i rhythms also had significant PER2 rhythms, except for three cells with marginal rhythmicity values (Fig. 1F). Peaks of [Ca2+]i preceded peaks of PER2 by an average of 7.3 ± 2.2 h (circular mean ± SD, n = 32; Fig. 1E).
In dispersed SCN neurons, the RI value for [Ca2+]i rarely exceeded that for PER2 (Fig. 1F). To better understand the relationship of PER2 and [Ca2+]i rhythmicity, we ran simulations of dispersed neurons in a mathematical model that included both TTFL and electrophysiology. Three parameters influencing the period and amplitude of TTFL and one parameter influencing [Ca2+]i concentration were independently and randomly modulated in 400 dispersed neurons (Vasalou and Henson, 2010), as described in Materials and Methods. Similar to experimental data (Fig. 1F), the amplitude of [Ca2+]i was limited by the amplitude of PER (Fig. 1G). Two mechanisms encoded in the model enabled it to achieve behavior similar to that observed experimentally (Fig. 1G): (1) [Ca2+]i rhythms were dependent on TTFL rhythms, and (2) TTFL rhythms persisted at baseline levels even in the absence of [Ca2+]i oscillations (To et al., 2007).
A rhythmic SCN neuronal network reinforces both PER2 and [Ca2+]i rhythmicity
To examine the effects of the SCN network on PER2 and [Ca2+]i rhythmicity, we compared the rhythmicity of single SCN neurons oscillating independently in dispersed cell cultures to those oscillating synchronously in slice cultures. Because we sliced 250-μm coronal sections, the SCN was usually divided into two slices (Noguchi and Watanabe, 2008). The slice containing a larger part of the SCN was used for imaging. Cells located in various parts of SCN slices were evenly selected and manually tracked in imaging data over several days (4.7–6.8 d). Both PER2 and [Ca2+]i rhythmicity were analyzed for 45 cells in 3 SCN slices for 4 d (Fig. 2A–C). Single-cell PER2 and [Ca2+]i rhythms appeared more robust in slice cultures than in dispersed cell cultures (Fig. 2D). In SCN slices, as in dispersed cells, [Ca2+]i rhythm waveforms were more variable than PER2 rhythm waveforms. In SCN slices, PER2 and [Ca2+]i patterns were rhythmic in 100% and 71% of cells, respectively (PSD-regression test, p < 0.05). Using RI as a measure of rhythm strength, single-cell PER2 and [Ca2+]i rhythms were both significantly stronger in slices than in dispersed cultures (p < 0.05, mixed effect model; Fig. 2E). To measure coherence between PER2 and [Ca2+]i waveforms, cross-correlations were performed for both dispersed cells and cells in slices. The cross-correlation values were significantly higher for cells in slices (Fig. 2F, p < 0.05, mixed effect model), suggesting higher similarity between PER2 and [Ca2+]i waveforms or improved rhythmicity. Peaks of [Ca2+]i preceded peaks of PER2 by an average of 4.8 ± 2.6 h (circular mean ± SD, n = 32; Fig. 2G), which was significantly different from dispersed cells (7.3 ± 2.2 h; circular mean ± SD, n = 32; p < 0.01, Watson–Williams test; Fig. 1E). Compared with dispersed cells, cells in slices also exhibited a significantly larger variability of this difference between PER2 and [Ca2+]i peak time (p < 0.05, Wallraff’s nonparametric test of circular homoscedasticity).
Spatiotemporal relationships between PER2 and [Ca2+]i rhythms
To analyze spatiotemporal patterns of PER2 and [Ca2+]i activation in the SCN and how the neuronal network may influence those patterns, peak times of PER2 and [Ca2+]i across the SCN were computationally detected in five SCN slices. Cell-like regions exhibiting significant rhythmicity in either PER2 or [Ca2+]i are depicted as small circles. Spatiotemporal patterns of PER2 and [Ca2+]i peaks and phase difference between PER2 and [Ca2+]i peaks in all five slices are shown in Fig. 2H–K. Generally, peak times of PER2 occurred earlier in the dorsomedial region and later in the ventrolateral region, as reported in previous studies (Evans et al., 2011b; Foley et al., 2011; Fig. 2H). [Ca2+]i and PER2 phases generally showed a similar spatial pattern (Fig. 2H–J). Phase differences between PER2 and [Ca2+]i were relatively uniform across SCN, with [Ca2+]i peaks preceding PER2 peaks by ∼4–6 h (Fig. 2K). Phase relationships between PER2 and [Ca2+]i peaks in each slice were stable during the entire 3- to 7-d durations of imaging experiments.
A minority of cells in the dorsomedial region expressed very different phase relationships, including [Ca2+]i peaks antiphase to or lagging PER2 peaks. These distinct phase differences were caused by changes in [Ca2+]i phasing rather than PER2 phasing (Fig. 2H–K). We hypothesized that this unusual phase relationship might arise from a difference in rostrocaudal position of the slice within the SCN. Additional SCN slices containing either a rostral or a caudal tip of the SCN (Noguchi and Watanabe, 2008), were imaged. However, unusual phasing was not observed in these extreme rostral (n = 1) or caudal (n = 3) slices.
Synchronous SCN network electrical activity reinforces rhythmic PER2 expression
To determine whether reinforcement of PER2 or [Ca2+]i rhythms requires a rhythmic, synchronized SCN network, we compared effects of TTX (blocking neuronal firing) on rhythmicity in slice and dispersed cultures. Both PER2 and [Ca2+]i rhythmicity were analyzed for 69 SCN neurons in 5 slices and 45 SCN neurons in 3 dispersed cell cultures. Cultures were imaged in control medium for 3–4 d, then in medium containing 1 µM TTX for 3.0 d, and finally in fresh control medium for 3.8 d.
For SCN neurons in slice cultures, on average, TTX significantly decreased PER2 RI to 70% of the value before TTX, and refreshing medium resulted in partial recovery to 84% of the pre-TTX control value (Fig. 3A, C). TTX had more variable effects on [Ca2+]i rhythms: e.g., no effect or loss of rhythmicity. On average, TTX decreased [Ca2+]i RI to 60% of the value before TTX, and refreshing medium resulted in recovery to 95% of the pre-TTX control value (Fig. 3A, C). TTX also significantly reduced average PER2 expression to 42% of the value before TTX, and refreshing medium restored it to a level similar to that before TTX, consistent with previous reports (Yamaguchi et al., 2003; Fig. 3A, D). Similarly, the coherence between PER2 and [Ca2+]i waveforms significantly decreased during TTX application and recovered after refreshing medium (Fig. 3E).
For SCN neurons in dispersed cultures, unlike in slice cultures, TTX did not significantly affect PER2 or [Ca2+]i RI (Fig. 3B, F); nor did TTX significantly affect average PER2 expression or coherence (Fig. 3G, H). For most individual cell traces, TTX had no discernible effect on PER2 rhythmicity (Fig. 3B). Just as in slice cultures, TTX had more variable effects on [Ca2+]i rhythms: e.g., no effect (Fig. 3B), loss of rhythmicity, or gain of rhythmicity. However, such changes in [Ca2+]i rhythmicity were also sometimes observed spontaneously in control conditions without TTX, and on average TTX had no significant effect on [Ca2+]i RI (Fig. 3F).
Ryanodine receptors have been suggested to regulate circadian [Ca2+]i release (Ikeda et al., 2003). We applied 100 µM ryanodine to dispersed SCN cultures (n = 20 cells, in 2 cultures), but no significant difference was observed in either PER2 or [Ca2+]i RI (Fig. 1–3).
Mathematical models predict that PER and [Ca2+]i rhythm amplitude are reduced by TTX in a multicellular SCN network but not in a single SCN neuron
To explore possible reasons that TTX suppressed PER2 and [Ca2+]i rhythmicity in SCN cells in slices but not in dispersed cultures, effects of TTX were simulated in both single-cell and multicellular SCN mathematical models (Vasalou and Henson, 2010, 2011; Vasalou et al., 2011). Neuronal firing rate was set to 0 Hz during TTX application, abolishing VIP and GABA signaling. In the multicellular model, desynchronization of individual cells and the absence of VIP receptor activation during TTX application reduced the amplitude of PER and [Ca2+]i oscillations (Fig. 3I). In contrast, in the single-cell model, there is no VIP or GABA intercellular signaling even at baseline, in the absence of TTX; consequently, applying TTX had no effect on amplitude of PER and [Ca2+]i oscillations (Fig. 3J).
[Ca2+]i rhythmicity requires an intact transcriptional clock
To test whether clock gene rhythms are necessary to generate [Ca2+]i rhythms in SCN neurons, we imaged [Ca2+]i rhythms in SCN slice cultures from mice deficient in the essential clock gene Bmal1 (Bunger et al., 2000; Ko et al., 2010). Cells were manually tracked in imaging data (4.9–6.7 d). Both PER2 and [Ca2+]i rhythms were analyzed for 35 cells in 4 slice cultures for 4 d. Both PER2 and [Ca2+]i showed stochastic fluctuations rather than constant levels or circadian rhythms (Fig. 4A). PER2 RI, [Ca2+]i RI, and coherence between PER2 and [Ca2+]i traces were significantly and dramatically decreased compared with WT controls (45 cells in 3 SCN slices; Fig. 4B, C).
Discussion
To study the causal relationships among SCN clock gene rhythms, [Ca2+]i rhythms, neuronal firing rhythms, and the SCN neuronal network, we simultaneously imaged single-cell PER2 and [Ca2+]i rhythmicity in both SCN slice and dispersed cell cultures. Our results demonstrate that [Ca2+]i rhythms in SCN neurons (1) are cell autonomous and do not require a synchronized network (Fig. 1), yet (2) are reinforced by the rhythmic neuronal network in SCN slices (Fig. 2), (3) do not require circadian rhythms of neuronal firing (Fig. 3), (4) have flexible phasing modulated by the neuronal network (Fig. 2), and (5) require an intact TTFL (Fig. 4).
In dispersed SCN neurons, circadian rhythms of spontaneous firing or PER2 expression free-run independently with their own intrinsic periods (Welsh et al., 1995; Webb et al., 2009), establishing that firing and PER2 rhythms are cell autonomous. The presence of [Ca2+]i rhythms in independently oscillating dispersed SCN neurons, as shown here (Fig. 1), establishes that [Ca2+]i rhythms are also cell-autonomous.
Many dispersed SCN neurons did not show [Ca2+]i rhythms yet exhibited robust PER2 rhythms. The strength of [Ca2+]i rhythmicity was always limited by PER2 rhythmicity, i.e., never clearly stronger than PER2 rhythmicity in a given cell. These experimental results are consistent with mathematical simulations, in which [Ca2+]i rhythms were not necessary to generate clock gene rhythms in single cells. We propose that this analysis method, plotting a rhythmicity strength of two parameters, can be used to test the causality of two rhythmic parameters. Cytoplasmic Ca2+ is, of course, necessary for cell survival and various physiologic functions including circadian rhythm generation (Lundkvist et al., 2005; Colwell, 2011). However, it appears that only tonic levels of [Ca2+]i (and not circadian rhythmicity of [Ca2+]i) are necessary for these processes. Our data suggest that [Ca2+]i oscillation works as a clock-controlled messenger to regulate circadian cellular functions rather than as a clock component to generate circadian rhythms.
Both PER2 and [Ca2+]i rhythms, but especially the latter, were greatly strengthened in SCN slice cultures relative to dispersed cultures, suggesting that the SCN neuronal network amplifies [Ca2+]i rhythms. Coherence between PER2 and [Ca2+]i rhythms was also improved in slice cultures, suggesting that control of PER2 over [Ca2+]i is stronger in slices than in dispersed cells. We cannot, however, rule out the possibility that improved rhythmicity in SCN slices was at least partly artifactual, owing to scattered light from neighboring cells smoothing the rhythm wave form and reducing noise. This concern could be addressed by future experiments restricting [Ca2+]i reporter expression to only a few cells per SCN slice, e.g., by gene gun (Ikeda et al., 2003).
Several lines of studies consistently suggested that TTX application reduces Per1 and Per2 mRNA and protein levels in cultured SCN slices (Yamaguchi et al., 2003; Webb et al., 2009). However, previous studies have found variable effects of TTX on [Ca2+]i rhythms in SCN slices (Colwell, 2000a; Ikeda et al., 2003; Enoki et al., 2012; Hong et al., 2012; Brancaccio et al., 2013). In our SCN slice cultures, TTX decreased both PER2 and [Ca2+]i rhythmicity, as well as PER2 expression levels. However, interestingly, in dispersed cell cultures, which already lack coupling sufficient to synchronize cellular rhythms, TTX did not significantly affect either PER2 or [Ca2+]i rhythms. Mathematical models of single-cell and multicellular SCN correctly simulated these distinct TTX effects, indicating that acute reduction of Pers by TTX within a day are caused by loss of network reinforcement. They also explain why TTX induces nonphotic-type phase shifts in neonatal rat SCN slices (Noguchi and Watanabe, 2005) but not in dispersed cells (Welsh et al., 1995).
On average, [Ca2+]i peaked ∼7.3 h before PER2 for single SCN neurons in dispersed cell cultures, and ∼4.8 h before PER2 for single SCN neurons in slice cultures. These results are similar to the 7.2-h phase difference observed previously in SCN slices (Brancaccio et al., 2013). In our study, the average phase difference between PER2 and [Ca2+]i rhythms in single cells was significantly different for SCN slices versus dispersed cell cultures (Figs. 1E and 2G), with a wider distribution of phase differences among cells in slices than among dispersed cells. These results suggest that the SCN neuronal network modulates the phase relationship between PER2 and [Ca2+]i rhythms.
More detailed computational spatiotemporal peak time analysis showed that [Ca2+]i peaked ∼4–6 h before PER2 in most cells in SCN slices (Fig. 2H–S), as suggested in previous work (Enoki et al., 2012; Brancaccio et al., 2013). However, in some slices, a small dorsal SCN region stably maintained [Ca2+]i peak time nearly antiphasic to the PER2 peak (Fig. 2J). These divergent groups of cells were not seen in rostral or caudal slices. Possibly, the divergent phase relationship may occur selectively in slices containing VIP-positive cells, which are more abundant in the central part of the SCN (Abrahamson and Moore, 2001; Noguchi and Watanabe, 2008).
Compared with the stable sinusoidal waveform of PER2 rhythms, [Ca2+]i rhythm waveforms were more variable. For example, [Ca2+]i rhythms showed narrow spike type waveforms (Figs. 1C and 2D ), sinusoidal waveforms, two daily [Ca2+]i peaks, or arrhythmicity (Fig. 1B). The diverse phases and waveforms of [Ca2+]i rhythms in dispersed SCN neurons suggest that [Ca2+]i rhythms can be modulated not only by the intercellular neuronal network, but also by intracellular modulators linking the TTFL to [Ca2+]i. Such intracellular mechanisms could involve store-operated Ca2+ entry (SOCE; Ikeda et al., 2003; Shen et al., 2011). Further investigations are necessary to determine the specific mechanisms and possible significance of variable [Ca2+]i waveforms and variable phase relationships between PER2 and [Ca2+]i rhythms.
It has been reported that circadian rhythms of [Ca2+]i and electrical activity can persist despite treatment with Per1 + Per2 antisense oligonucleotides, which significantly reduced PER immunoreactivity (Sugiyama et al., 2004). Thus, we considered the possibility that [Ca2+]i rhythms could be driven by TTFL-independent mechanisms. In this study, we showed that [Ca2+]i rhythms were definitively abolished in Bmal1–/– SCN slices, in which PER2 is also nonrhythmic. This is consistent with the hypothesis that an intact TTFL drives [Ca2+]i rhythms in SCN neurons (Enoki et al., 2017a). In the previous experiment using Per1 + Per2 antisense oligonucleotides (Sugiyama et al., 2004), [Ca2+]i rhythms may have been rescued by a few cells in which the TTFL was not completely suppressed by antisense oligonucleotides (Liu et al., 2007). In some cases, Bmal1–/– SCN slices (but not dispersed cells) have been observed to generate stochastic PER2 oscillations with variable periods overlapping the circadian range (Ko et al., 2010). We did not observe such unstable oscillations in this study, possibly because of toxicity of viral vectors or [Ca2+]i buffering by GCaMP3 (Looger and Griesbeck, 2012).
Studies in SCN slices have found that Ca2+ rhythms depend on Gq signaling (Brancaccio et al., 2013) and ryanodine-sensitive intracellular Ca2+ stores (Ikeda et al., 2003), but not on gap junctions (Enoki et al., 2012). Ryanodine receptors have been shown to modulate temporal [Ca2+]i release in the mouse SCN and to be necessary for rat behavioral rhythm output (Mercado et al., 2009; Aguilar-Roblero et al., 2016). In this study, ryanodine had little effect on circadian rhythms of [Ca2+]i or PER2. Differences in experimental conditions, such as developmental stage and reporters, might have led to different [Ca2+]i responses, but, our report is consistent with previous work in finding that ryanodine does not affect PER2 rhythms (Aguilar-Roblero et al., 2016). Possible mechanisms by which clock genes control spontaneous firing of SCN neurons include regulation of ion channels at the transcriptional level (e.g., BK calcium-activated potassium channels; Meredith et al., 2006), via cAMP [e.g., hyperpolarization-activated, cyclic nucleotide-gated (HCN) ion channels; O’Neill et al., 2008], or via redox status (Wang et al., 2012). However, the mechanism by which the TTFL exerts rhythmic control over [Ca2+]i release from intracellular stores is completely unknown. SCN microarray studies reveal circadian oscillation of several interesting Ca2+ related genes (e.g., calreticulin and calumenin, which are Ca2+-binding proteins located in ER, and NCX1, which is a sodium/calcium exchanger expressed in plasma membrane and mitochondria; Panda et al., 2002; Jung and Kim, 2004; Barbuti and DiFrancesco, 2008). Further investigation may reveal interesting connections between the TTFL and Ca2+ regulatory proteins.
In summary, in this study, we demonstrate autonomy of circadian [Ca2+]i rhythms in dispersed SCN neurons, SCN network reinforcement of both PER2 and [Ca2+]i rhythms and modulation of the phase relationship between the two rhythms, and dependence of SCN [Ca2+]i rhythms on the transcriptional circadian clock but not on neuronal firing within individual cells. How the transcriptional circadian clock drives [Ca2+]i rhythms, and how SCN neuronal networks modulate [Ca2+]i rhythm phase relative to PER2, are intriguing questions for future study. Much work remains to be done to elucidate the complete mechanism of the SCN circadian clock, including cytosolic signaling pathways and intercellular networks.
Acknowledgments
Acknowledgments: We thank Dr. Loren Looger (Janelia Farm, HHMI) for advice on using GCaMP3 and Dr. Pamela Mellon (UCSD) for donating Bmal1–/– mice. We thank Kevin Lo for technical assistance.
Footnotes
Authors report no conflict of interest.
Funding sources: Supported by NIH (R01 MH082945 to DKW) and a V.A. Career Development Award (DKW). TLL gratefully acknowledges support from the Amherst College Faculty Research Award Program funded by the H. Axel Schupf ‘57 Fund for Intellectual Life.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.