Abstract
In mammals, the suprachiasmatic nucleus (SCN) of the hypothalamus coordinates daily rhythms including sleep–wake, hormone release, and gene expression. The cells of the SCN must synchronize to each other to drive these circadian rhythms in the rest of the body. The ontogeny of circadian cycling and intercellular coupling in the SCN remains poorly understood. Recent in vitro studies have recorded circadian rhythms from the whole embryonic SCN. Here, we tracked the onset and precision of rhythms in PERIOD2 (PER2), a clock protein, within the SCN isolated from embryonic and postnatal mice of undetermined sex. We found that a few SCN cells developed circadian periodicity in PER2 by 14.5 d after mating (E14.5) with no evidence for daily cycling on E13.5. On E15.5, the fraction of competent oscillators increased dramatically corresponding with stabilization of their circadian periods. The cells of the SCN harvested at E15.5 expressed sustained, synchronous daily rhythms. By postnatal day 2 (P2), SCN oscillators displayed the daily, dorsal-ventral phase wave in clock gene expression typical of the adult SCN. Strikingly, vasoactive intestinal polypeptide (VIP), a neuropeptide critical for synchrony in the adult SCN, and its receptor, VPAC2R, reached detectable levels after birth and after the onset of circadian synchrony. Antagonists of GABA or VIP signaling or action potentials did not disrupt circadian synchrony in the E15.5 SCN. We conclude that endogenous daily rhythms in the fetal SCN begin with few noisy oscillators on E14.5, followed by widespread oscillations that rapidly synchronize on E15.5 by an unknown mechanism.
SIGNIFICANCE STATEMENT We recorded the onset of PER2 circadian oscillations during embryonic development in the mouse SCN. When isolated at E13.5, the anlagen of the SCN expresses high, arrhythmic PER2. In contrast, a few cells show noisy circadian rhythms in the isolated E14.5 SCN and most show reliable, self-sustained, synchronized rhythms in the E15.5 SCN. Strikingly, this synchrony at E15.5 appears before expression of VIP or its receptor and persists in the presence of blockers of VIP, GABA or neuronal firing. Finally, the dorsal-ventral phase wave of PER2 typical of the adult SCN appears ∼P2, indicating that multiple signals may mediate circadian synchrony during the ontogeny of the SCN.
Introduction
Daily patterns in behavior and physiology are coordinated by circadian pacemakers, which allow organisms to anticipate environmental cycles such as light/dark (Pittendrigh, 1960; Hut and Beersma, 2011). In mammals, the hypothalamic suprachiasmatic nucleus (SCN) is a master circadian pacemaker consisting of ∼20,000 neurons. Individual cells can generate circadian rhythms via a negative feedback loop involving core clock genes including Per2 (Hastings et al., 2014; Takahashi, 2017). How and when these rhythms arise during development is not known.
In mice, SCN neurogenesis occurs between embryonic days E10–E15 (where E1 is the day when a vaginal plug is detectable), peaking between E12–E14 and forming from ventrolateral to dorsomedial (Shimada and Nakamura, 1973; Kabrita and Davis, 2008; Shimogori et al., 2010). Transcription factors such as Lhx1, Shh, Six3, and Six6 participate in the embryonic specification of SCN lineages (Shimogori et al., 2010; VanDunk et al., 2011; Clark et al., 2013; Bedont et al., 2014). However, the markers of terminal differentiation and synaptogenesis peak postnatally (Moore and Bernstein, 1989; Shimogori et al., 2010). The rat SCN has fetal day/night differences in 2-deoxyglucose uptake in vivo ∼E19, in vivo daily oscillations in clock genes including Per2 by E20–E21 and in vitro cycling of electrical activity on the day of birth (E22; Reppert and Schwartz, 1983, 1984; Shibata and Moore, 1987; Sládek et al., 2004; Kováciková et al., 2006; Houdek and Sumová, 2014). The mouse SCN has detectable fetal daily rhythms in Per1 transcript levels at E17, and PERIOD1 (PER1) and PER2 proteins by E18 in vivo (Shimomura et al., 2001; Ansari et al., 2009). It is not clear whether these results reflect diurnal rhythms imposed by the mother in utero. Recent studies in vitro have found daily expression of the PER2::LUCIFERASE (Per2Luc) bioluminescent reporter in the whole SCN as early as E13 (Landgraf et al., 2015) or E15 (Wreschnig et al., 2014). It is unclear, however, if or when individual cells in the fetal SCN begin to oscillate and to synchronize their daily rhythms to each other.
To drive coherent rhythms in behavior, cells within the SCN must synchronize their genetic oscillations to an identical period and a consistent phase relationship (Mohawk and Takahashi, 2011; Herzog et al., 2015; Evans and Gorman, 2016). Vasoactive intestinal polypeptide (VIP) is required for normal synchrony in the SCN. Daily addition of VIP entrains SCN neurons, and the absence of VIP or its receptor, VPAC2R, reduces synchrony between neurons and consequently eliminates many daily rhythms (Harmar et al., 2002; Colwell et al., 2003; Aton et al., 2005; Maywood et al., 2006; Ciarleglio et al., 2009). It is not known which intercellular signals underlie circadian synchrony during development (Wang et al., 2014; Ono et al., 2016). In this study, we examined the embryonic onset of endogenous SCN circadian oscillation and intercellular synchrony.
Materials and Methods
Animals.
All mice were maintained on a C57BL/6JN background (WT) and housed in a 7:00 A.M. to 7:00 P.M. light/dark cycle in the Danforth Animal Facility at Washington University. The Per2Luc mouse line was generated by replacing the endogenous mouse Period2 gene locus with a Per2Luc reporter construct (Yoo et al., 2004). For immunochemistry, we compared male and female homozygous Per2Luc pups to Vip−/− and Vipr2−/− pups as negative controls. Vaginal plugs confirmed overnight mating of each female. We designated the morning after mating as E0.5 and the day of birth as postnatal age 0. All procedures were approved by the Animal Care and Use Committee of Washington University and followed National Institutes of Health guidelines.
Cultures and bioluminescence recording.
Pregnant mice were killed with CO2 and cervical dislocation and their embryos dissected into 4°C Hanks buffered saline solution (Sigma-Aldrich). All surgeries started at 1:00 P.M. (Zeitgeber time 06). We recorded bioluminescence from 300 μm coronal SCN slices from fetal (E13.5–E15.5 and E17.5) and postnatal day (P)2 homozygous Per2Luc mice (Landgraf et al., 2015). Briefly, brains were embedded in a block of 4% low-melting agarose and prepared with a vibraslicer. SCN were dissected with scalpels from sections of the ventral hypothalamus and placed on 0.4 mm membrane inserts (Millipore) in sealed 35 mm Petri dishes (BD Biosciences) with 1 ml of DMEM (Sigma-Aldrich, pH 7.2), supplemented with 25 U/ml penicillin, 25 μg/ml streptomycin (Invitrogen), 10 mm HEPES (Sigma-Aldrich), 2% B27 (Invitrogen), 0.35 g/L NaHCO3 (Sigma-Aldrich), and 0.15 mm beetle luciferin (Promega). SCN explants were transferred to a light-tight incubator at 36°C (Onyx, Stanford Photonics). We collected images (XR Mega-10AW camera, Stanford Photonics) every 6 min (6 s exposure) and then summed every 10 frames with ImageJ software (http://rsbweb.nih.gov/ij) to provide one image every hour. We applied adjacent frame minimization to each movie to filter out bright noise caused by the camera or cosmic rays.
Data analysis.
We analyzed the Per2Luc rhythms of the developing SCN at the level of the whole tissue and with 30 μm (pixel) resolution. The statistical analysis for each experiment is described in detail in its respective section. We assessed rhythms measured as the integrated intensity of the SCN tissue with cosinor analysis (software generous gift from B. Meier and A. Kramer, Charité University Berlin). We performed a pixel-based analysis of the same SCN movies. We excluded movies if the corners of the SCN slice moved (3 slices excluded at each of E15.5, E17.5, and P2). We measured periodicity of each pixel (1 px = 900 μm2) and defined pixels as circadian if they had a dominant period between 18 and 32 h (p < 0.05 by Lomb–Scargle periodogram analysis (Ruf, 1999), or a correlation coefficient >0.6 in cosinor analysis. Because the two methods yielded 90 ± 0.02% (mean ± SE; n = 30 brains) agreement in pixel classification, we reported results from the periodogram analysis only. The cycle to cycle variability in Per2Luc was evaluated by detrending rhythmic cells with a discrete wavelet transform (Leise, 2015) keeping only detail coefficients between 16 and 32 h, and measuring the interpeak distance (h) of each cycle. We evaluated the synchronization index using the Kuramato method (Kuramoto, 2003).
A radial distribution of mean peak times was constructed for each SCN to identify the presence of core-shell structure. The time of the initial peak for each pixel was calculated following detrending. A radial distribution of peak times was constructed by binning individual cell peaks based on radial distance from a reference point in the left core. A null distribution of radial peak times was calculated by bootstrapping: reassigning peak time values at random to pixels within each SCN and recalculating the radial distribution using 10,000 bootstrap runs. The SCN was determined to have core-shell structure if the radial peak time distribution showed significant (p < 0.05) peaks for the shell region (early peaking, middle distances) and at least one core region (late-peaking, near and far distances).
Software.
Data analysis was performed in the Python language, using packages scipy and numpy (Jones et al., 2014), PyWavelets, Statsmodels (Seabold and Perktold, 2010), and Matplotlib (Hunter, 2007). Parallelization of data processing was achieved via iPython (Pérez and Granger, 2007). Code for data processing is available on request.
Immunocytochemistry.
Brains of mice from ages E13.5–E17.5, P0, and P2 were removed and fixed in 4% paraformaldehyde in PBS overnight at 4°C. We transferred them to 30% sucrose in PBS for 24 h at 4°C and sectioned in 30-μm-thick coronal slices onto microscope slides. The sections were incubated for 1 h at 4°C in a blocking solution (10% nonfat dry milk, 10% bovine serum albumin and 0.3% Triton X-100 in PBS) and then for 48 h at 4°C in an anti-VIP (1:1000; ImmunoStar) or anti-VPAC2R (1:1000; Abcam) rabbit polyclonal antibody diluted in 3% BSA and 0.35% Triton X-100 in PBS. After washing in PBS, the slides next were incubated at room temperature in a biotinylated goat anti-rabbit IgG (1:200; Vector Laboratories) for 2 h and then ABC solution (1:200; Vectastain Elite ABC Kit, Vector Laboratories) for 2 h. Finally, we incubated each slide in 50 mm Tris-HCl using 3,3-diaminobenzidine kit (Sigma-Aldrich). After each step, we rinsed three times for 15 min in PBS. Mounted sections were dehydrated through a series of ethanol and xylene washes. Slides were imaged (NanoZoomer microscope, Leica) and the integrated optical density of the SCN was measured using ImageJ software (Schneider et al., 2012).
Drug treatments.
All drugs were diluted in deionized water and they remained in the recording medium throughout the entire recording without further medium changes. VIP antagonist (200 nm PG-99465) provided by Dr. P. Robberecht (Cutler et al., 2003), GABA antagonist (100 μm Gabazine, Tocris Bioscience), Meclofenamic acid (MEC; 10 μm; Sigma-Aldrich), or TTX (2.5 μm tetrodotoxin).
Results
SCN cells become reliable circadian oscillators after E14.5
Recent in vitro studies have demonstrated the presence of circadian rhythms in PER2 expression in whole SCN explants at E13 (Landgraf et al., 2015) or E15 (Wreschnig et al., 2014), around the times of peak or completed neurogenesis, respectively. To test whether cell-autonomous oscillations appear before tissue-level rhythms, we imaged bioluminescence in fetal Per2Luc/Per2Luc SCN slices. SCN explanted at E13.5 displayed no tissue-level circadian rhythms (n = 8 SCN from 2 litters; Fig. 1). At this stage, the ventricular zone of dorsal hypothalamus expressed high levels of PER2, sometimes in a circadian pattern (n = 2 of the 8 slices; Movie 1), but the bilateral SCN region had low to no PER2 expression and no circadian rhythms. In contrast, all tissues explanted on E14.5 or later displayed bilateral Per2Luc expression in the ventral hypothalamus, corresponding to the SCN. By E14.5, 45% of SCN were circadian (5 of 11 SCN from 4 litters; Fig. 1). Notably, some E14.5 explants (n = 4 of 11; Movie 2) also displayed no rhythmic Per2Luc expression outside the SCN in the ventricular zone of the dorsal hypothalamus. By E15.5, no hypothalamic explants showed high PER2 expression outside the SCN region and nearly all SCN were circadian (Fig. 1; Movie 3; n = 11 of 12 E15 SCN from 4 litters, 7 of 7 E17 SCN from 2 litters, and 5 of 6 P2 SCN from 1 litter). We conclude that PER2 expression is localized to the SCN and becomes circadian on E14.5.
Circadian rhythms in the SCN during development. Long-term, real-time recordings of PER2 expression from SCN harvested after E14.5 showed reliable circadian periods between 18 and 32 h. A, Representative records of Per2Luc bioluminescence from SCN starting on the day they were explanted. Note the weak circadian oscillations that were found in some E14.5 SCN. Insets, Two representative SCN from each age. Some SCN harvested at E13.5 expressed Per2Luc above the SCN region along the ventricular zone of dorsal hypothalamus (Extra SCN, blue arrowhead; n = 5/8). All others harvested at this age and older reliably expressed high levels of PER2 in the bilateral SCN (red arrow). Scale bar: 10 px, 300 μm. B, Nearly all SCN expressed significant circadian rhythms when harvested on or after E15.5 (top), with periods close to 24 h (mean ± SEM; middle) and durations of daily Per2Luc expression above the mean close to 12 h (alpha; mean ± SEM; bottom; p < 0.05, one-way ANOVA, Tukey's HSD). Note, eight slices were examined at E13.5, however, we excluded five because they did not express detectable PER2 in the SCN region.
Movie 1. Time-lapse movie of PER2 expression from SCN explanted at E13.5. Note that, at E13.5, the bilateral SCN region had little to no PER2 expression and no circadian rhythms. At this stage, the bilateral ventricular zone of the dorsal hypothalamus had high PER2 levels. PER2 expression increased in the E14.5 SCN, decreased in the dorsal hypothalamus, and showed circadian rhythms when explanted on or after E15.5. A daily dorsal-to-ventral wave of PER2 appeared around P2.
Time-lapse movie of PER2 expression from SCN explanted at E14.5. Note the increased expression in the SCN and decreased expression in the dorsal hypothalamus.
Time-lapse movie of PER2 expression from SCN explanted at E15.5. The bilateral SCN showed circadian rhythms when explanted on or after E15.5.
We next developed methods to culture and image the embryonic SCN for at least 5 d. From these movies, we quantified local circadian rhythms where each pixel included ∼10 or fewer cells (Fig. 2A,B). We found the fraction of circadian pixels (cells) increased approximately sevenfold from E14.5 (Fig. 2C,D; 0.10 ± 0.07, mean ± SEM; n = 11 SCN from 4 litters) to E15.5 [0.72 ± 0.11; n = 9 SCN from 4 litters; *p < 0.05, one-way Kruskal–Wallis H(3) = 10.62, with post hoc Tukey's Honestly Significant Difference (HSD)]. Nearly all regions were circadian in SCN explanted on E17.5 (0.96 ± 0.03; n = 4 SCN from 2 litters) and postnatal day 2 (P2, 0.97 ± 0.02; n = 3 SCN from 1 litter). The developmental increase in the number of circadian regions correlated with an increase in circadian period in the SCN (Fig. 3A,B). We conclude that a small number of SCN cells initiate endogenous circadian rhythmicity by E14.5, and by E15.5 SCN cells are rhythmic throughout the SCN.
Embryonic SCN cells express PER2 rhythmically beginning ∼E14.5 in vitro. A, Heat maps of two SCN explanted at each developmental stage show circadian pixels based on Lomb–Scargle periodogram (power > p < 0.05) and cosinor analysis (correlation coefficient >0.6). We therefore report results only from periodogram analyses for simplicity. B, Representative Per2Luc bioluminescence images of the same SCN explants (1 h integration, 2 × 2 binning). Scale bar: 10 px, 300 μm. C, Representative Per2Luc bioluminescence recordings from single rhythmic pixels from SCN cultured at different embryonic stages. Note the competence of the E15.5 SCN to generate and sustain circadian rhythms in PER2 protein expression. D, The fraction of circadian pixels in each SCN (mean ± SEM) was higher for explants at E15.5 and older [*p < 0.05, one-way ANOVA, Tukey's HSD; n = at least 3 SCN at each age where age is reported as days post-mating (E0.5) and post-birth (E20.5 = P0)].
Circadian period and precision of SCN cells increased with development. A, Heat maps from representative SCN show the mean period over 5 d of recording. B, The fraction of rhythmic pixels correlated with an increase in the mean SCN period (least-squares linear regression: p = 0.00019, Pearson's r = 0.74. SCN with <5 circadian pixels were considered arrhythmic and excluded from this analysis; E14.5: diamonds, n = 6 of 11; E15.5: circles, n = 9 of 9; E17.5: triangles, n = 4 of 4; P2: squares, n = 3 of 3). C, Representative mean cycle-to-cycle variability heat maps show the increase in circadian precision with age. Scale bar, 150 μm. D, The average variability of the daily period across all rhythmic pixels within each isolated SCN inversely correlated with the fraction of circadian cells (least-squares linear regression: p = 0.0006, Pearson's r = −0.70; symbols as in B), SCN with <5 circadian pixels were considered arrhythmic and excluded from this analysis.
To determine the precision of SCN rhythms during development, we measured the cycle-to-cycle difference in period for each pixel within the fetal SCN over 3 d of the recording. Interestingly the increase in the number of circadian oscillators correlated with an increase in circadian precision (Fig. 3C,D), suggesting that addition of circadian cells lengthens and stabilizes period in the SCN.
The phase wave typical of adult SCN appears ∼P2
In adults, PER2 expression progresses as a daily wave from dorsal to ventral across the SCN (Yan and Okamura, 2002; Quintero et al., 2003; Yamaguchi et al., 2003; Evans et al., 2011). We examined when this spatiotemporal patterning of circadian expression arises in the SCN. To do so, we identified the times of daily peak PER2 in each pixel of the SCN and tested for spatial organization of expression. We found reliable daily waves of PER2 from the dorsal to ventral SCN by P2, but not at earlier developmental stages (Movies 4, 5; Fig. 4). We conclude that phase relationships among SCN cells continue to mature after birth.
Time-lapse movie of PER2 expression from SCN explanted at E17.5. Note the high amplitude, in-phase circadian rhythms across the SCN.
Time-lapse movie of PER2 expression from SCN explanted at P2. Note the appearance of a daily dorsal-to-ventral wave of PER2 across the SCN.
The phase wave of PER2 expression in the SCN appeared only after birth. A, Heat maps of the time of peak PER2 expression across representative SCN. Note the dorsal-to-ventral distribution across the cells of the P2 SCN. Scale bar, 300 μm. B, Radial distribution of peak times revealed a core-shell structure in the P2, but not younger SCN explants. We compared the average differences in the times of peak PER2 expression as a function of distance from the ventral margin of each SCN with a null distribution of randomly shuffled peak times (shaded area = 95% CI for mean peak time calculated from 10,000 resamples). At P2, PER2 peaked significantly later in the ventral core than the dorsal shell, as seen in this representative example.
Onset of synchrony and intercellular communication within the SCN
Synchronous circadian oscillations in the SCN are required to coordinate daily rhythms in behavior (Schwartz et al., 1987; Aton et al., 2005; Ohta et al., 2005). To test how synchrony develops, we measured the cycle-to-cycle variability in period length and the synchrony index (SI), a quantity that ranges from 1.0 (all oscillators peak together) to 0.0 (all cells peak at uniformly different times of day) from 5 d of recording. We found that synchrony tends to increase between E14.5 (Fig. 5; SI = 0.80 ± 0.5, n = 6 SCN) and E15.5 (SI = 0.87 ± 0.05, n = 9 SCN) reaching a maximum ∼E17.5 (SI = 0.91 ± 0.02, n = 4) and P2 (SI = 0.94 ± 0.02, n = 3 SCN; 1 unilateral due to changes during recording between halves). Although there is no significant differences in the SI between ages, SI was positively correlated with the fraction of SCN cells displaying daily rhythms, indicating that intercellular synchrony is developed simultaneously with increased endogenously generated cellular oscillation (Fig. 5C).
The synchrony between SCN cells increased and became more stable as the fraction of circadian cells increased. A, Representative traces of the synchronization index (also called the Kuramato order parameter or Rayleigh statistic, r) among regions within the cultured SCN at different ages. B, The average synchronization index of SCN slices (mean ± SEM; 24 h excluded from either end of analysis) did not reliably increase with age, but (C) correlated with the fraction of circadian cells in the SCN (symbols as in Fig. 3; least-squares linear regression: p = 0.004, Pearson's r = 0.61, n = at least 3 SCN at each age).
We next examined the signaling pathways that could enable synchronization of SCN cells starting ∼E15.5. VIP signaling is necessary for coordinated circadian rhythms among SCN cells and in behavior of adult mice (Harmar et al., 2002; Colwell et al., 2003; Aton et al., 2005; Maywood et al., 2006; Ciarleglio et al., 2009). We measured the levels of VIP, and its cognate receptor, VPAC2R, by immunocytochemistry between days E14.5 and P2 (Fig. 6). We found VIP (12.5 ± 2.3, n = 7) and VPAC2R [15.0 ± 5.6 relative optical density (ROD), n = 6] were detectable postnatally above baseline (immunolabeling in the P2 SCN of Vip−/− = 2.1 ± 0.7 ROD, n = 5 brains; Vipr2−/− = 3.8 ± 2.1 ROD, n = 5, mean ± SEM), but not earlier stages. To further explore this surprising result, we found that a VIP antagonist did not affect synchrony of the E15.5 SCN (Fig. 7). Consistent with prior publications, we found that addition of the VPAC2R antagonist reduced the fraction of circadian cells in the isolated adult SCN to 30% from 84% (n = 1 adult SCN; Cutler et al., 2003; Aton et al., 2005). We conclude that VIP and its receptor VPAC2R appear in the SCN after the initiation of circadian rhythms and are not required for circadian synchrony among embryonic SCN cells.
VIP and VPAC2R expression matured in the postnatal SCN. A, Immunostaining intensity (mean relative optical density ± SEM; n = 9, 5, 7, 6, 7 and 7 brains at each age, respectively) of VIP in the SCN was significant by P2 compared with Vip knock-out controls (n = 5; *p < 0.05, one-way ANOVA with Tukey's HSD). B, Immunostaining intensity (mean ± SEM) of VPAC2R (n = 6, 6, 5, 5, 7, and 6 brains, respectively) was significant by P0 compared with Vipr2 knock-out controls (n = 5; *p < 0.05, one-way ANOVA with Tukey's HSD). C, Representative coronal SCN sections show the VIP (left) and VPAC2R (right) immunoreactivity at different ages. Scale bar, 200 μm.
Antagonists of GABA or VIP signaling, gap junctions or neuronal firing did not disrupt circadian synchrony in the developing SCN. Drugs applied during the first day of culture of E15.5 SCN explants did not change (A) the fraction of pixels with circadian Per2Luc (gray bars show the means across all SCN recorded at each developmental stage), or (B) the mean sync index over the 5 d of recording. C, The synchronization index within a representative cultured E15.5 SCN treated with Veh (Vehicle), Gabazine (Gz; 100 μm), MEC (10 μm), VIP antagonist (200 nm; PG-99465), or TTX (2.5 μm).
We then performed further pharmacological experiments to test whether GABA and/or gap junctions, other candidate coupling mechanisms, might drive synchrony in the SCN before VIP expression begins. Remarkably, antagonists against either or both GABA and gap junctions did not reduce the synchrony index of the E15.5 SCN (Fig. 7).
Discussion
Our results indicate that E15.5 is a critical time in the maturation of the SCN. Real-time imaging revealed that the clock gene, PER2, appears on E14.5 in the SCN and ∼10% of SCN cells are weakly circadian. By E15.5, nearly all cells within the SCN develop precise circadian rhythms. The cells of the E15.5 SCN are thus competent circadian oscillators. Previous work related this age with the end of neurogenesis (Shimada and Nakamura, 1973; Kabrita and Davis, 2008; Shimogori et al., 2010). At this age the expression of transcription factors as Lhx1 and RORα have a wider distribution across the SCN (VanDunk et al., 2011). Interestingly, RORα is a positive regulator of Bmal1 expression (Sato et al., 2004), Bmal1 is an essential transcription factor in circadian rhythmogenesis and activator of Per gene expression (Hastings et al., 2014; Takahashi, 2017). In adults, RORα knock-outs showed slight behavioral abnormalities in circadian rhythms. We therefore hypothesize that cell autonomous expression of RORα or another factor at E15.5 induces the onset of widespread, sustained circadian oscillations in SCN cells.
Studies using stem cells showed gradual development of circadian oscillations after inducing differentiation with retinoic acid (Yagita et al., 2010; Inada et al., 2014), suggesting that circadian rhythms can be initiated without direct maternal input early during differentiation. Consistent with a prior publication (Landgraf et al., 2015), we found little to no evidence that SCN cells become circadian spontaneously in vitro. For example, SCN explanted on 13.5 were arrhythmic and did not show an increase in the fraction of rhythmic pixels over the course of the experiment. SCN isolated at E14.5 or E15.5 similarly did not show changes in the amplitude, period, variability, or fraction of circadian cells during 7 d of recording. Although we cannot rule out that arrhythmic, fetal SCN cells cultured >7 d might develop daily cycling, our results indicate that the transition from arrhythmic to circadian is rapid and likely requires signals from outside the fetal SCN.
We found no circadian rhythms at E13.5, a small fraction of cells expressed noisy circadian rhythms at E14.5 and nearly all cells were circadian by E15.5. Notably, we found that PER2 was first expressed outside the SCN. At E14.5, extra-SCN PER2 expression was high and occasionally circadian. This extra-SCN expression diminishes by E15.5. It is not clear whether these cells move into the SCN or remain in the dorsal hypothalamus while reducing their PER2 expression. Inclusion of these extra-SCN regions in explants could explain the discrepancy between the two studies that reported the onset of PER2 daily oscillations in hypothalamic explants at E13 or E15 (Wreschnig et al., 2014; Landgraf et al., 2015). Importantly, the correlations in the increasing fraction of rhythmic cells, increasing circadian period and increasing intercellular and intracellular period precision during this time could reflect the consequences of cell–cell communication (Webb et al., 2009) and/or maturation of molecular components for circadian rhythm generation (VanDunk et al., 2011; Ono et al., 2016). Previous studies demonstrated that embryonic SCN have a strong phase-setting in response to perturbations such as surgical isolation (Wreschnig et al., 2014; Landgraf et al., 2015). Together, these data suggest that SCN cells become intrinsically circadian between E14.5 and E15.5. The presence of a functional circadian clock early in development may indicate the importance of precisely coordinating processes in the fetus such as cell division and differentiation, which are adversely affected in Clock mutant mice (Du Pré et al., 2014; Inada et al., 2014; Kobayashi et al., 2015).
The SCN in adult mice displays a dorsal to ventral phase wave in clock gene expression (Yan and Okamura, 2002; Quintero et al., 2003; Yamaguchi et al., 2003; Evans et al., 2011). Because this phase wave is not a fixed property (e.g., can adopt a different pattern after a temperature manipulation in vitro (Jeong et al., 2016), it is important to consider whether the results from the isolated SCN reflect its in vivo circadian properties. Critically, our results indicate that under constant conditions, the SCN phase wave first appears after birth. Because we isolated SCN at the same time of day, regardless of developmental age, we conclude that, although surgery can alter circadian phase (Wreschnig et al., 2014; Landgraf et al., 2015), it is unlikely to initiate rhythms in the fetal SCN or the spatial wave seen in the more mature SCN. This postnatal maturation of this SCN network likely involves synaptogenesis which begins ∼P2 when we found the onset of the phase wave (Bedont and Blackshaw, 2015). It will be intriguing to map the synaptic network topology that contributes to the phase relationships among SCN cells (Abel et al., 2016; Buijink et al., 2016).
Though the network continues to mature after birth, we found that intercellular synchrony is achieved much earlier. A prior report found that surgery can reset rhythms in the fetal SCN (Landgraf et al., 2015). Here, we found that synchrony of E14.5 and immature E15.5 SCN decreased over the recording, indicating weak or no synchronization between cells. Previous studies also relate an increase in the period variability with a loss of synchrony between cells (Honma et al., 2004; Webb et al., 2009). Our data showed an increase in synchrony by E15.5 (Fig. 5) corresponding to a decrease in cycle-to-cycle period variability, supporting this hypothesis. In adults, VIP and its receptor VPAC2R play a critical role in synchronizing oscillators in the SCN (Harmar et al., 2002; Colwell et al., 2003; Aton et al., 2005; Maywood et al., 2006; Ciarleglio et al., 2009). Although the expression of VIP mRNAs has been detected at E18 in rat (Ban et al., 1997; Houdek and Sumová, 2014) and VIP protein expression has been studied during postnatal ages (Herzog et al., 2000), our results show that the protein is not detectable by immunochemistry before birth. Interestingly, embryonic SCN explants from mice lacking VIP maintain circadian oscillations at the tissue level (Wreschnig et al., 2014). Moreover, our results demonstrate that synchrony occurs at E15.5, before the expression of VIP or receptor VPAC2R. That the exogenous application of VIP antagonist does not decrease synchrony suggests that other signals are likely involved in sustaining synchrony before the presence of VIP. There is evidence that dopamine, melatonin, or rhythmic feeding may be involved in establishing synchrony between the mother and fetal SCN to coordinate rhythmic cellular processes (Davis and Mannion, 1988; Weaver and Reppert, 1989; Viswanathan et al., 1994). The dopamine pathway in particular merits further investigation as a potential mechanism driving synchrony in the fetal SCN.
The increase in synchrony in the SCN at E15.5 does coincide with Lhx1 expression across the SCN (VanDunk et al., 2011). This transcription factor is important for the expression of diverse neuropeptide coupling signals aside from VIP, such as GRP or AVP (Bedont et al., 2014). Because these signals can mediate the synchrony in adults (Brown et al., 2005; Maywood et al., 2011) and in postnatal ages (Ono et al., 2016), they are candidates for mediating SCN synchrony in embryonic development. In rats the first synapses appear in perinatal ages and the sharp increase between P4 and P10 (Moore and Bernstein, 1989) suggests that other mechanisms could mediate synchrony before synapse formation. Gap junctions have recently been shown to not contribute to circadian synchrony in the adult SCN (Diemer et al., 2017). Our results, with a blocker of gap junctions, further indicate that they do not play a role in establishing synchrony during SCN development. We conclude that cellular circadian oscillations occur simultaneously with the onset of synchrony by unknown mechanisms.
Footnotes
This work was supported by postdoctoral fellowship 237875 Consejo Nacional de Ciencia y Tecnología (CONACYT), the March of Dimes Center for Prematurity, and NINDS Grant 095367. We thank members of the Herzog laboratory for valuable discussions.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Erik D. Herzog, Department of Biology, Box 1137, Washington University, St. Louis, MO 63130-4899. herzog{at}wustl.edu
