Lhx1 maintains synchrony among circadian oscillator neurons of the SCN

Abstract

The robustness and limited plasticity of the master circadian clock in the suprachiasmatic nucleus (SCN) is attributed to strong intercellular communication among its constituent neurons. However, factors that specify this characteristic feature of the SCN are unknown. Here, we identified Lhx1 as a regulator of SCN coupling. A phase-shifting light pulse causes acute reduction in Lhx1 expression and of its target genes that participate in SCN coupling. Mice lacking Lhx1 in the SCN have intact circadian oscillators, but reduced levels of coupling factors. Consequently, the mice rapidly phase shift under a jet lag paradigm and their behavior rhythms gradually deteriorate under constant condition. Ex vivo recordings of the SCN from these mice showed rapid desynchronization of unit oscillators. Therefore, by regulating expression of genes mediating intercellular communication, Lhx1 imparts synchrony among SCN neurons and ensures consolidated rhythms of activity and rest that is resistant to photic noise.

https://doi.org/10.7554/eLife.03357.001

eLife digest

As anyone who has experienced jet lag can testify, our sleeping pattern is normally synchronized with the local day–night cycle. Nevertheless, if a person is made to live in constant darkness as part of an experiment, they still continue to experience daily changes in their alertness levels. In most individuals, this internal ‘circadian rhythm’ repeats with a period of just over 24 hr, and exposure to light brings it into line with the 24-hr clock.

The internal circadian rhythm is generated by a structure deep within the brain called the suprachiasmatic nucleus (SCN), which is essentially the ‘master clock’ of the brain. However, each cell within the SCN also contains its own clock, and can generate rhythmic activity independently of its neighbors. Cross-talk between these cells results in the production of a single circadian rhythm.

Now, Hatori et al. have identified the master regulator that controls this cross-talk. When mice living in 24-hr darkness were exposed to an hour of light in the early evening, they showed changes in the levels of proteins associated with many SCN genes. But one gene in particular, known as Lhx1, stood out because it was strongly suppressed by light.

Mice with a complete absence of Lhx1 die in the womb. However, mice that lose Lhx1 during embryonic development survive, although they struggle to maintain circadian rhythms when kept in complete darkness. This is not because their SCN cells fail to generate circadian rhythms. Instead, it is because the loss of Lhx1—a transcription factor that controls the expression of many other genes—means that the SCN cells do not produce the proteins they need to synchronize their outputs.

As well as identifying a key gene involved in the generation and maintenance of circadian rhythms, Hatori et al. have underlined the importance of cell-to-cell communication in these processes. These insights may ultimately have therapeutic relevance for individuals with sleep disturbances caused by jet lag, shift work or certain sleep disorders.

https://doi.org/10.7554/eLife.03357.002

Introduction

Circadian clocks generate ∼24 hr rhythms in behavior and physiology which allow an organism to anticipate and adjust to environmental changes accompanying the earth's day/night cycle. These rhythms are generated in a cell-autonomous manner by transcription–translation based feedback loops which are composed of clock proteins, such as PERIOD (PER1, PER2, and PER3), CRYPTOCHROME (CRY1 and CRY2), CLOCK, BMAL1, REV-ERB (REV-ERBα and REV-ERBβ), and ROR (RORα, RORβ, and RORγ) in mammals (Mohawk et al., 2012). These oscillatory loops reside in almost all tissue types and regulate their downstream effectors to generate oscillations in the steady-state mRNA levels of thousands of genes in a tissue-specific manner. Identification of tissue-specific circadian transcripts in peripheral organs has elucidated the mechanism by which circadian clocks dictate the temporal regulation of organ function. The tissue level clocks are organized in a hierarchical manner. The hypothalamic suprachiasmatic nucleus (SCN) composed of ∼20,000 densely packed neurons acts as the master clock by orchestrating molecular oscillations in peripheral tissues (Welsh et al., 2010). Tight intercellular communication among SCN neurons (coupling) drives synchronous oscillations. This, in turn, imparts overt rhythms in activity-rest and dependent rhythms in physiology and metabolism of the whole organism. While coupling between the SCN neurons buffers against the noise in oscillations of the constituent neurons, it is plastic enough to allow adaptive resetting of the phase of the SCN oscillator in response to changes in the environment.

Light is the principal cue for entraining the SCN circadian clock to environmental cycles. Light stimuli are perceived in the retina and transmitted to the SCN via melanopsin (OPN4)-expressing retinal ganglion cells (mRGCs) (Hatori and Panda, 2010). The time-of-the-day specific response (called ‘gating’) of the SCN to light pulses properly adjusts the phase of the circadian clock. In mice held under constant darkness, light administered at subjective daytime, subjective evening, or subjective late night causes no shift, phase delay or phase advance of the behavioral rhythm respectively. This phase-dependent light response is conserved across species (Pittendrigh, 1967; Zatz et al., 1988; Schwartz and Zimmerman, 1990). It is known that illumination at night triggers extensive chromatin remodeling in the mouse SCN (Crosio et al., 2000), which likely results in changes in the levels of a large number of transcripts. However, except for a few dozen transcripts including the clock genes Per1 and Per2 (Zhu et al., 2012; Jagannath et al., 2013), the extent of light-triggered transcriptional changes in the SCN is largely unknown. These acute transcriptional changes impinge on the molecular oscillator to adjust the phase of the mRNA rhythms in the SCN.

It is becoming increasingly clear that coupling among the SCN neurons buffers against phase shifts, and the transient weakening of such coupling facilitates large phase shifts. SCN neurons exhibit tight intercellular communication imposed by paracrine peptidergic signals such as VIP (Vasoactive intestinal polypeptide), AVP (Arginine vasopressin), and GRP (Gastrin-releasing peptide) (Welsh et al., 2010; Hogenesch and Herzog, 2011). As the deficiency of Vip or its receptor Vpac2r causes desynchronization among SCN neurons (Aton et al., 2005), this peptide-mediated coupling mechanism is a unique and necessary feature of the SCN in order to generate robust synchronous rhythms. Weaker coupling among the SCN neurons is suggested to facilitate rapid and large phase shifts of the overt rhythms (Herzog 2007; An et al., 2013). Therefore, transcription factors that regulate expression of the SCN coupling agents are central to the unique function of the SCN.

Here, using a combination of behavioral, genetic, and genomic tools, we vastly expand the understanding of the dynamic transcriptional landscape of the SCN. We conducted comprehensive analysis of the light-regulated, circadian, and tissue-enriched protein-coding transcriptome of the mouse SCN in order to understand the specificity of the master circadian clock. We found that the SCN-enriched transcription factor Lhx1 (LIM homeobox 1) is required for expression of a number of genes including Vip whose protein product participates in intercellular signaling. The SCN-specific loss of Lhx1 attenuates cell to cell coupling of cell-autonomous oscillators in the SCN, abolishing circadian behavioral activity consolidation in vivo.

Results

Light-regulated and circadian-controlled transcripts of the SCN

To comprehensively identify circadian-, light-regulated, and SCN-enriched protein-coding transcripts in the adult SCN, following 2 weeks of entrainment to 12 hr light:12 hr dark (LD) cycles male C57BL6/J mice were transferred to constant darkness (DD), and the SCN was collected every 2 hr over 48 hr. Light at subjective night, but not during subjective day is known to cause a behavioral phase shift. To assess the gene expression effect of a phase-shifting pulse of light, a subset of mice were exposed to a 1 hr white light pulse delivered at 30 hr, 40 hr, and 46 hr after the onset of DD (Figure 1A, Figure1—figure supplement 1) representing subjective daytime (CT6), early evening (CT16), and late night (CT22). We discovered 1412 genes (Figure 1—figure supplement 1; Supplementary file 1) that show circadian oscillation in transcript abundance (pMMCβ <0.05, pFGT <0.05, and median temporal expression >100). Gene expression measured at 1 hr, 2 hr, and 4 hr after the beginning of the 1 hr light pulse revealed 508 probesets whose levels changed (up- or down-regulated) in response to at least one of the three light pulses (Figure 1A,B, Figure 1—figure supplement 1; Supplementary file 2A). However, only 84 (17%) were also rhythmic (Figure 1—figure supplement 1; Supplementary file 2B).

Figure 1 with 3 supplements see all
Light-regulated transcripts of the SCN.

(A) Heatmap rendering of light-regulated SCN transcripts. For each time point, fold change between respective light treated and dark control was plotted. (B) Circadian gating of light-modulated transcripts. Cutoffs of two fold were set for up-regulation (blue) or suppression (red) after light pulse, and the number of probesets that satisfy each cutoff was plotted for each point. Quantitative RT-PCR (qRT-PCR) expression confirmation of genes detected as light-regulated by microarray. Examples of genes (C) induced or (D) repressed by light pulses at three different time points. (Mean +s.e.m., n = 4).

https://doi.org/10.7554/eLife.03357.003

We tested whether the number of light-modulated transcripts parallels the phase shifting effect of light. Light pulse delivered during subjective day (CT6) triggered changes in a small number of transcripts, while the same pulse at CT16 causes a large phase shift accompanied by a large number of transcripts changing >2 fold (Figure 1B, Figure 1—figure supplement 1; Supplementary file 2C). The magnitude of light modulation of transcripts also paralleled the known phase shifting effect of light in genetic models of light input perturbation. Compared to wild type (WT) C57BL6/J, the light response was unchanged in mice lacking rod/cone photoreceptors, attenuated in melanopsin (Opn4)-deficient mice and completely abolished in mice lacking all three photopigments or those with specific ablation of melanopsin expressing retinal ganglion cells (Opn4Cre/+;R26iDTR/+ + diphtheria toxin [DT]) (Figure 1—figure supplement 2). Overall, the light-induced transcriptional responses in the SCN closely correlate with the extent of light-induced behavioral phase shift (Hatori and Panda, 2010) and hence are informative of the mechanism and consequences of the phase shift.

The light-induced genes include the circadian clock components e.g., Per1, Per2, and Bhlhe40 (Dec1) and genes involved in the CREB and MAPK signaling pathways (Zhu et al., 2012; Jagannath et al., 2013), which are implicated in intracellular signaling leading to resetting the phase of cell autonomous oscillators (Figure 1C). The novel group of light repressed transcripts was enriched for those involved in inter-cellular communication including Vip, Avpr1a, Rasd1, and Npy6r (Figure 1D, Supplementary file 2C), suggesting that an effective phase shift of the SCN clock rests both on resetting the phase of cell autonomous clock and on light-induced relaxation of the intercellular coupling.

The data for circadian and light-dependent gene expression are also being made available in a user friendly searchable web interface at http://scn.salk.edu. The database can be queried using a gene symbol or a probeset identifier as the keyword.

Lhx1 is an SCN-enriched and light-modulated gene

Since the tight intercellular coupling of circadian oscillators is largely an SCN specific phenomenon (Herzog 2007), we reasoned that light likely suppresses the expression of an SCN-enriched factor that coordinates intercellular communication. We employed three-step enrichment criteria comparing the transcriptome of the SCN with that of 82 other mouse tissues including 14 different neural tissues (Su et al., 2004) (Figure 1—figure supplement 3). This analysis identified 213 SCN-enriched (not SCN-exclusive) genes (Figure 1—figure supplement 3, Supplementary file 3), including Rorα, Rorβ, Vip, Grp, Rgs16, and Prokr2 which are known to play important roles in SCN function (Welsh et al., 2010; Doi et al., 2011; Kasukawa et al., 2011; VanDunk et al., 2011).

Among the 13 SCN-enriched transcription factors discovered (Figure 1—figure supplement 3), Lhx1 mRNA was suppressed by light (Figure 1D), raising the possibility that it regulates the expression of SCN synchronizing agents and that the lack of Lhx1 might render the SCN prone to desynchrony. However, Lhx1 is a necessary factor for differentiation of several tissue types as Lhx1−/− embryos die ∼ E10 (Shawlot and Behringer, 1995). In the hypothalamus, Lhx1 expression begins at E11 and parallels that of Six6 (VanDunk et al., 2011), whose function is necessary for normal SCN development (Clark et al., 2013). Since Rorα expression is also SCN enriched (Figure 1—figure supplement 3) and its developmental expression follows Lhx1 expression in the SCN region, we generated RorαCre;Lhx1loxP mice for testing the role of Lhx1 in SCN function.

Loss of Lhx1 in the SCN affects circadian consolidation and light induced phase shifts of behavioral rhythms

The overlap between Rorα and Lhx1 expression in the adult brain is largely restricted to the SCN (Figure 1—figure supplement 3) and the dLGN (Chou et al., 2013). The Rorα/Lhx1 double positive dLGN neurons constitute the thalamocortical pathway for conveying visual information to the brain (Chou et al., 2013). Since this circuit is not implicated in circadian photoentrainment (Mure and Panda, 2012), we reasoned the RorαCre;Lhx1loxP mice will be appropriate for testing the role of Lhx1 in SCN.

Rorα is an essential component of the cell autonomous circadian oscillator (Sato et al., 2004). Its mRNA becomes detectable from E15 onward (VanDunk et al., 2011). Therefore, in the early developmental stage between E11 and E15, Lhx1 is expected to be expressed in the RorαCre;Lhx1loxP double mutant mice, which would allow normal Lhx1 function (if any) during early SCN differentiation, while uncovering the post-developmental role of Lhx1 in SCN function.

The RorαCre mouse has an IRES;Cre cassette knocked-in downstream of the Rorα locus (Chou et al., 2013) which permits normal expression of Rorα and co-expression of Cre. RorαCre;R26R mice (Figure 2A and Figure 2—figure supplement 1) or RorαCre;Z/AP mice (Figure 2B and Figure 2—figure supplement 1) showed robust Cre-dependent LacZ or ALPP (alkaline phosphatase) expression in the SCN region of the hypothalamus. ALPP staining of RorαCre;Z/AP mice revealed uniform staining of SCN neurons along dorso-ventral and rostro-caudal axes (Figure 2—figure supplement 1). As opposed to the developmental and circadian dysfunction in Rorα mutant (Sato et al., 2004), the RorαCre allele did not compromise Rorα function as the RorαCre mice exhibited normal development and circadian activity rhythm in light–dark cycles and in constant darkness that are indistinguishable from those of wild-type mice (see below).

Figure 2 with 2 supplements see all
Loss of Lhx1 expression in the SCN renders faster synchronization with change in LD regimes.

Enriched expression of a Rorα-driven marker in the SCN in RorαCre;R26R mice. (A) Ventral view of a whole brain (magnified view on the right) of adult RorαCre;R26R shows LacZ staining of the SCN. (B) Coronal section through the mid-SCN region (scale bar, 1 mm) and the magnified view of the SCN (scale bar, 100 µm) showing LacZ expression or (C) alkaline phosphatase expression in RorαCre;Z/AP mice. (D) qRT-PCR estimate of Lhx1 expression in the SCN (mean +s.e.m, n = 5). (E) Normal SCN innervation of the retinal ganglion cells in the WT mice as revealed by monocular injection of CTB-conjugated fluorescent marker is intact in (F) Lhx1SCN−KO mice. A 1 hr light pulse at CT16 causes (G) upregulation of light-induced genes (Per1, Per2, cFos, JunB), while (H) the light-suppressed transcripts (Lhx1, Vip, Avpr1a) in the WT SCN show reduced expression in the Lhx1SCN-KO mice. Mice were in DD for 2 days before the light pulse. Representative actograms of (I) RorαCre/Cre, (J) RorαCre/+;Lhx1loxP/loxP, and (K) RorαCre/Cre;Lhx1loxP/loxP mice subjected to 8 hr phase advance and 8 hr delay. (I) Average (+s.e.m., n = 5–8) activity onset and (K) average (+s.e.m.) number of days to re-entrain to advance or delay in light onset in three genotypes. Color codes in L and M correspond to the labels in IK.

https://doi.org/10.7554/eLife.03357.007

A single copy of Cre in RorαCre/+:Lhx1loxP/loxP mice reduced Lhx1 mRNA levels by nearly 40%, while in RorαCre/Cre:Lhx1loxP/loxP (Lhx1SCN-KO) mice Lhx1 expression in the SCN is severely reduced (Figure 2D). The gross morphology of the SCN in these mice remains intact (Figure 2—figure supplement 1) suggesting that the conditional loss of Lhx1 post early development in the SCN does not severely affect its differentiation unlike the loss of Six3, Six6, or Math5 (Wee et al., 2002; VanDunk et al., 2011; Clark et al., 2013). The activity pattern of these mice entrains normally to an imposed LD cycle (Figure 2—figure supplement 2) implying functional innervation of the SCN by the mRGCs, which is known to occur postnatally (McNeill et al., 2011). Accordingly, anterograde labeling using a Cholera Toxin B conjugated fluorescent marker indicated normal innervation of the SCN by RGCs (Figure 2E,F and Figure 2—figure supplement 1).

Light-induced upregulation of immediate early genes and clock components including c-Fos, JunB, Per1, and Per2 in the SCN of Lhx1SCN-KO mice was comparable to that in WT mice (Figure 2G). However, light suppressed transcripts involved in intercellular communication such as Vip and Avpr1a showed reduced basal expression under DD (Figure 2H). This acute induction of the phase-resetting branch of light input along with potentially weak intercellular coupling suggested that the Lhx1SCN-KO mice may be more susceptible to light-induced phase shifts. Accordingly, in response to an 8 hr phase advance or delay of the LD cycle, the activity onset of RorαCre/+;Lhx1loxP/loxP and RorαCre/Cre;Lhx1loxP/loxP mice readjusted much faster than the WT mice to the new LD regime irrespective of the direction of the shift (Figure 2I–L).

Next, we tested the consequence of the potentially weak intercellular communication in the absence of light. Under constant darkness (DD), the circadian locomotor activity rhythms of the Rorα+/+; Lhx1+/+, Rorα+/+;Lhx1loxP/loxP, RorαCre/+;Lhx1+/+, or RorαCre/Cre;Lhx1+/+ mice were comparable (Figure 3A–D, Table 1). The RorαCre/+;Lhx1loxP/loxP mice showed normal activity rhythm for up to 3 weeks under DD, after which the activity consolidation deteriorated with no apparent ∼24 hr rhythm (Figure 3E). The RorαCre/Cre;Lhx1loxP/loxP mice showed circadian activity rhythm for up to 4 days under constant darkness, after which they became arrhythmic (Figure 3F, Figure 3—figure supplement 1). The lack of circadian locomotor activity rhythm in Lhx1SCN-KO mice under DD does not result from the disruption of the cell autonomous circadian oscillator, since the median expression of core clock components Per1 and clock output gene Dbp largely remained equivalent in the SCN of Lhx1SCN-KO and wild-type cohorts (Figure 3G,H). Both transcripts showed a significantly dampened rhythm with reduced peak levels and increased expression at the trough, suggestive of oscillator desynchrony. On the other hand, transcripts participating in intercellular communication including Vip, Avpr1a, Rasd1, Pde7b, Creb3l1, and a cell matrix associated cell–cell interaction mediator Nov were significantly reduced in the Lhx1SCN-KO mice (Figure 3I).

Figure 3 with 2 supplements see all
Lhx1 sustains normal circadian activity rhythms by regulating expression of synchronizing factors.

(AF) Representative wheel running activity pattern over several days of LD followed by DD in wild type and mice lacking Lhx1 in the SCN. Double-plotted qRT-PCR quantification (average +s.e.m, n = 3–4 mice) of (G) Per1 and (H) Dbp in the SCN of DD adapted WT and Lhx1SCN−KO mice. (I) Average (+s.e.m. 8 time points every 3 hr over 24 hr) expression of several factors involved in intercellular communication or circadian clock in the SCN of dark adapted WT and Lhx1SCN−KO mice. (J) Average mRNA (+s.e.m., n = 3–6 mice, *p < 0.05) expression or (K) immunoreactivity of VIP is reduced in the SCN of Lhx1-deficient mice. (L) Transcriptional activation of mouse Vip promoter by mouse LHX1. pGL3-promoter vector was used as a control promoter vector. Values are mean +s.e.m, ANOVA **p<0.01, ***p<0.001 vs 0 ng (white bar).

https://doi.org/10.7554/eLife.03357.010
Table 1

Circadian running wheel activity period length of various mouse strains under constant darkness (n = 7–22)

https://doi.org/10.7554/eLife.03357.013
Average (h)SEM (h)
Rora+/+;Lhx1+/+23.780.23
Rora+/+;Lhx1loxP/loxP23.640.05
RoraCre/Cre;Lhx1+/+23.400.14
RoraCre/+;Lhx1+/+23.860.07
RoraCre/+;Lhx1loxP/loxP23.510.13
RoraCre/Cre;Lhx1loxP/loxPNANA
  1. Mice showing arrhythmic activity were excluded from the analysis. (NA = Not Applicable).

Lhx1 regulates expression of Vip

Among the downregulated genes, the mRNA levels of Vip were undetectable in the Lhx1SCN−KO mice (Figure 3J), and VIP protein level was also reduced in the SCN (Figure 3K). Hence, we tested whether LHX1 regulates Vip expression. Transcription from a 1 kb promoter region of mouse Vip driving a luciferase reporter was activated by wild-type mouse and human Lhx1 in a dose-dependent manner (Figure 3L and Figure 3—figure supplement 2). Such transcriptional regulation was dependent on the DNA binding function of LHX1 because the LHX1N230S mutant failed to activate transcription from Vip:luc reporter (Figure 3L). This mutant (Figure 3—figure supplement 2) carries a missense mutation in a highly conserved asparagine residue known to disrupt DNA binding (Thaler et al., 2002). Furthermore, wild-type LHX1 also failed to activate transcription from SV40 promoter, supporting the idea that Vip expression is specifically activated by LHX1 and that this requires the DNA-binding activity of LHX1. Collectively, these results demonstrate that LHX1 regulates expression of Vip, which along with sufficient expression of other synchronizing agents maintains synchrony among SCN neurons.

Lhx1 mediates synchrony among SCN oscillator neurons

Coupling among SCN neurons supports synchronous oscillations of individual oscillators that otherwise show variation in intrinsic period lengths (Aton et al., 2005). To directly test the role of Lhx1 in maintaining SCN synchrony, we recorded electrophysiological activity of SCN slices using a multielectrode array. SCNs of both WT and Lhx1SCN-KO mice showed no apparent difference in the firing frequency, supporting the notion that Lhx1SCN-KO has no discernible developmental defects in the SCN. The WT SCN from LD reared animals showed synchronous firing rhythm, as characterized by both robust oscillation of multiunit activity as well as synchrony of peak phase of activity among different channels. The Lhx1SCN-KO SCN, on the other hand, showed remarkable dampening of firing rhythm as well as phase dispersion on day 1 and became almost asynchronous within 3 days (Figure 4A–E). This parallels the timeline of emergence of behavioral arrhythmicity in Lhx1SCN-KO mice when they are released from LD to DD (Figure 3F). Similarly, in WT mice maintained in DD for 2 weeks, the SCN displayed circadian oscillation of the multi-unit activity with a similar peak phase among the different electrodes examined (Figure 4—figure supplement 1). In contrast, the slices from RorαCre/Cre;Lhx1loxP/loxP animals showed widely dispersed phases of the multi-unit activity peak, a finding coherent with the arrhythmic locomotor activity observed in these animals in DD (Figure 3A). Since VIP expression is severely reduced in the Lhx1SCN-KO mice, we tested if extrinsic supplementation of VIP can restore the synchrony of the SCN neurons. Daily application of VIP for 1 hr on SCN slices from DD adapted Lhx1SCN-KO mice restored the normal synchrony and waveforms of the SCN firing rhythm (Figure 4F,G). Finally, the gradual dampening of SCN multiunit activity and desynchrony in the LD-adapted SCN slice from Lhx1SCN-KO mouse can be reversed by the daily application of VIP (Figure 4H,I).

Figure 4 with 1 supplement see all
Lhx1 maintains synchrony among SCN neurons partly via VIP.

(A) Average (+s.e.m.) normalized multiunit activity (MUA) recorded from representative SCN slices of LD-adapted WT (n = 40, black) and Lhx1SCN−KO (n = 12, orange) mice. Data were binned every 60 min. Representative normalized MUAs and peak phase of activity from WT SCN (B and C, n = 40) and Lhx1SCN−KO mouse (D and E, n = 19). For C and E, left and right panels are respectively for days 1 and 3. (F) Average MUA of a DD adapted Lhx1SCN−KO SCN that received 1 hr perfusion of VIP daily for up to 7 days. (G) Peak phases of activity are gradually synchronized over 7 days. (H) Representative MUA from the SCN of an LD-adapted Lhx1SCN−KO mouse over several days. During the first 4 days, the activity dampened, which was rescued by daily application of VIP. Down arrows in F and H indicate the time of VIP application. (I) Peak phases of activity in WT and Lhx1SCN−KO SCN at the end of 7 days are shown.

https://doi.org/10.7554/eLife.03357.014

Discussion

Tissue-specific gene expression and temporal changes in transcript levels govern and facilitate tissue function. Although the SCN has long been recognized as the master circadian oscillator and the principal target for light modulation of circadian rhythms in mammals, the molecular basis for these at the protein-coding transcript level have not been comprehensively identified. In this study, we present a thorough analysis of transcriptional oscillations in the SCN, its response to light at different phases of the circadian oscillator and a glossary of SCN-enriched genes. This detailed description of the transcriptional landscape of the SCN has a variety of implications for understanding SCN function. In this study, we focused on the relevance of oscillator synchrony in circadian rhythms and light-induced phase shift in behavior.

Under natural conditions, time-of-the-day dependent interaction of the SCN with light input from the retina specifies the phase of overt activity-rest rhythms (Hatori and Panda, 2010). Light during the subjective night induces expression of several genes including immediate early genes c-fos, activates MAPK pathway, acutely upregulates the Per1 transcript and resets the phase of overt behavioral rhythm (Obrietan et al., 1998, 1999; Dziema et al., 2003; Cheng et al., 2004; Butcher et al., 2005). We found that global transcriptional changes in the SCN parallel the sensitivity of the activity rhythms to light. During the subjective day, when light is ineffective in resetting overt rhythms, a very small fraction of the transcriptome changed in response to light. Conversely, during the subjective night, extensive changes in transcripts correlated with the large phase shifting effect of light, suggesting that changes in multiple pathways accompany the shift in the phase of the SCN oscillator network. Among the transcripts that change both during daytime and nighttime, their magnitude of expression was attenuated during the daytime. The magnitude of transcript changes also reflected the effectiveness of retinal photoreceptors in entraining the clock. As seen in the Opn4−/− mice (Panda et al., 2002b), attenuated light-induced change in overt rhythms paralleled the reduced transcriptional responses to light. Mice lacking mRGCs or those lacking rod, cone and melanopsin photopigments exhibited no behavioral phase shift in response to light (Hatori and Panda, 2010). They also showed no significant mRNA expression changes in transcripts in the SCN (Figure 1—figure supplement 2). Collectively, these observations support the notion that a light signal perceived through retina photoreceptors causes transcriptional changes in the SCN that parallel the behavioral phase-shifting effect of light.

Included among the light-induced transcripts were several components of intracellular signaling cascades, kinases, phosphatases, and transcription factors (MAPK pathway genes, SIK1, PER1, PER2, DEC1, EGR1, and EGR2), which likely function at the different signaling steps needed to transduce light information received at the plasma membrane to generate an appropriate change in the phase of the core clock (Obrietan et al., 1998, 1999; Dziema et al., 2003; Cheng et al., 2004; Butcher et al., 2005; Jagannath et al., 2013). The light-suppressed transcripts on the other hand included intercellular signaling agents, cell surface receptors, and associated signaling components. This indicated that stimulus-induced downregulation of the tight coupling among SCN neurons is a likely mechanism to allow desynchronization and then resynchronization of cellular clocks to reset the SCN to a new phase. Since tight intercellular communication among the SCN neurons is a characteristic feature of the SCN, we further focused on the subset of transcripts that are light-suppressed as well as enriched in the SCN. The search for SCN-specific genes among 82 different mouse tissues did not yield any protein coding gene that is exclusively expressed in the SCN, suggesting that a combinatorial gene expression signature specifies SCN development and function. LHX1 is one of the SCN-enriched and light-suppressed transcription factors, and hence we reasoned that LHX1 might regulate transcription of at least a subset of the intercellular communication agents that specify SCN neuronal network features.

During the preparation of this manuscript, another study showed Six3-Cre dependent loss of Lhx1 expression severely affected terminal differentiation as well as peptidergic outputs of the SCN (Bedont et al., 2014). By using Rora-Cre in our study, the loss of Lhx1 expression likely occurred after the onset of Rorα expression at E14.5, thereby allowing any terminal differentiation function of Lhx1 in SCN development to occur normally. In both studies, the loss of Lhx1 in the SCN did not affect the overall expression level of any of the clock components, while the expression of several genes implicated in cellular synchrony including Vip and Avpr1a were downregulated. AVP, VIP, and their cognate receptors have been implicated in maintaining inter-cellular communication among SCN neurons (Herzog, 2007). Light pulse also reduced the expression of Avpr1a and Vip mRNA. Mice carrying loss-of-function alleles of Vip or vasopressin receptors exhibit a weakened SCN network, progressive desynchronization of SCN neurons and increased sensitivity to light-induced phase shifts (Herzog, 2007; Yamaguchi et al., 2013). Similarly, the loss of AVP receptor weakens the network and increases the sensitivity of the behavioral rhythm to changes in the light regime (Herzog, 2007). While these observations have indicated that AVP- and VIP-mediated intercellular communication constitutes the framework for the SCN network, the transcription factor(s) that determines this SCN-specific property was unknown. Direct transcriptional induction from the VIP promoter by LHX1, the severe loss of VIP mRNA and immunoreactivity in the Lhx1−/− SCN along with a desynchronous SCN now establishes Lhx1 as a critical regulator of VIP production in the SCN.

In summary, we have discovered that Lhx1 is a master regulator of multiple factors including VIP that maintain robust coupling among SCN neurons after their differentiation into oscillator neurons. While the loss of Vip (Colwell et al., 2003), Avpr1a (Wersinger et al., 2007; Li et al., 2009), or Rasd1 alone causes a mild alteration in the circadian organization or light responses (Cheng et al., 2004), the parallel perturbation of multiple intercellular signaling components in the Lhx1SCN-KO mice indicates a critical role for Lhx1 in determining the specific feature of SCN neurons that impart coupling among neurons. The intercellular coupling is thus as important as the cell-autonomous oscillations for maintaining the consolidated rhythm of activity-rest that can resist abrupt changes in the ambient light conditions.

Materials and methods

Animals

All animal experiments were carried out in accordance with the guidelines of the Institutional Animal Care and Use Committee of the Salk Institute. Mice were housed under 12 hr light: 12 hr dark (LD) cycles. Food and water were available ad libitum. C57BL/6J, C3H/HeJ strain (rd) carrying Pdebrd1 mutation, Cre-dependent lacZ reporter strain (R26R) (Soriano, 1999), and Cre-dependent human ALPP reporter strain (Z/AP) (Lobe et al., 1999) were obtained from the Jackson Laboratory. Opn4−/− mice (Panda et al., 2002b) were bred to rd/rd to generate rdrd;Opn4−/−. Both Opn4Cre and Opn4Cre;R26iDTR were described in Hatori et al. (2008). The floxed Lhx1 allele (Lhx1loxP) mice were originally generated in Kwan and Behringer (2002). RorαCre mouse was generated by knocking in an IRES;Cre cassette 3′ downstream of the Rorα locus (Chou et al., 2013). Both Lhx1loxP and RorαCre mice were back-crossed to C57BL/6J strain for at least eight generations.

DNA microarrays

Request a detailed protocol

168 male C57BL/6J mice of 6 weeks of age were maintained for 3 weeks on a 12 hr light:12 hr dark cycle. For circadian gene expression profiling, after being placed in DD for 2 days, four animals were sacrificed every 2 hr, beginning at hour 30 of DD, which corresponds to CT18, for two complete 24 hr cycles (Figure 1A). For light-regulated gene profiling, mice were maintained in DD then exposed to 1 hr light at CT30, 40, or 46, while control (no light pulse) mice were left in dark. After 1 hr, all mice were returned to DD, and four animals each were collected at 1, 2, and 4 hr from the beginning of light exposure from each CT. Mice were sacrificed by cervical dislocation, and the optic nerves were cut under dim red light. The SCN was quickly dissected and four individual SCNs were pooled to be rapidly frozen on dry ice. Total RNA was extracted by RNeasy mini column (Qiagen, CA, USA). For each time point, 100 ng of total RNA was used as starting material for Affymetrix MOE430 high density arrays (Affymetrix, CA, USA). For RorαCre;Lhx1loxP mice, the SCNs were collected every 4 hr. Total RNA was processed for qRT-PCR analyses following standard protocols.

DNA microarray data analysis

Identification of circadian transcriptome in the SCN

Request a detailed protocol

Data were analyzed by using COSOPT to identify the transcripts that cycle with a time period between 20 and 30 hr, as described previously (Panda et al., 2002a). To identify gene families that are overrepresented in the rhythmically expressing genes in the SCN, we used the DAVID pipeline (Dennis et al., 2003).

Identification of the light-regulated genes

Request a detailed protocol

Average expression levels of individual probesets were compared between respective light and dark samples. Those showing a fold change of 2 or higher at any of the 9 time-points were further examined.

Identification of SCN-enriched genes

Request a detailed protocol

We used mouse tissue expression database (Su et al., 2004) to identify SCN-enriched transcripts. We focused on probesets that had much higher signal in the SCN (in terms of the median temporal expression) than the remaining 82 tissues (duplicate values for each tissue), 14 neural tissues (duplicate), and the hypothalamus. Mathematically, this is interpreted as a high standard deviation for the given probeset expression level in the SCN relative to other tissues. We thus defined two Z values to measure how many standard deviations away the SCN expression is relative to the remainder of the data set, as follows:

Z (all) = ([Median temporal expression in the SCN] − [Average expression across 82 × 2 + 1 tissues])/(Standard deviation of expression across 82 × 2 + 1 tissues)

Z (neural) = ([Median temporal expression in the SCN] − [Average expression across 2 × 14 + 1 neural tissues])/(Standard deviation of expression across 2 × 14 + 1 neural tissues)

To provide a further layer of stringency, we performed hierarchical clustering of all tissues to identify the anatomical region of the brain that is closest to the SCN in terms of the gene expression signature. This region was the anatomically adjacent hypothalamus. We thus set our third criterion for SCN specificity of a gene as a high value of the SCN/Hypothalamus gene expression values.

We shortlisted probe sets that had Z (all tissues) ≥ 1.63, Z (neural tissues) ≥ 1.44, and SCN/hypothalamus ≥6. Further stringency was imposed by the requirement that median temporal SCN expression value for the probeset must exceed 200. 230 probe sets meet these criteria. Notably, many of these had previously been described as SCN-specific genes (Welsh et al., 2010; Kasukawa et al., 2011) and thus validated our approach.

Quantitative RT-PCR (qRT-PCR)

One hour light from CT40

Request a detailed protocol

Opn4Cre/+;R26iDTR/+ and control Opn4+/+;R26+/+ mice were individually housed in a wheel running cage and maintained in 12:12 LD cycle. Mice were intraperitoneally injected with diphtheria toxin to ablate melanopsin-expressing RGCs in Opn4Cre/+;R26iDTR/+ mice as described (Hatori et al., 2008). Ablation of melanopsin RGCs was confirmed by the loss of light entrainment of circadian activity rhythms. In addition to these mice, rd/rd;Opn4−/−, Opn4−/− and rd/rd mice were maintained in 12:12 LD, then transferred to DD. Mice were exposed to 1 hr light from CT40 or kept in dark, and sacrificed to collect the SCNs at 1 hr from the beginning of light under light or dark.

Wild type and Opn4−/− mice under 1 hr light from CT30, 40, or 46

Request a detailed protocol

C57BL/6J and Opn4−/− mice were maintained in LD, then transferred to DD. Mice were exposed to 1 hr light from CT30, 40, or 46 or kept in dark, and sacrificed to collect the SCNs at 0, 1, 2, or 4 hr from the beginning of light under light or dark.

Measurement of mRNA levels in RorαCre;Lhx1loxP mice

Request a detailed protocol

RorαCre/Cre;Lhx1loxP/loxP and the control mice were kept in LD cycle, then transferred to DD. Mice were collected at CT18, 22, 26, 28, 30, 32, 36, and 40. For circadian gene expression, total RNA from each time point was subject to qRT-PCR. For non-oscillating transcripts, data from all timepoints were pooled to evaluate levels of expression.

qRT-PCR

Request a detailed protocol

Total RNA from an individual mouse SCN was extracted following RNeasy mini column protocols (Qiagen, CA, USA). cDNA synthesis was carried out with Superscript III (Invitrogen, CA, USA) or qScript cDNA SuperMix (Quanta Biosciences, MD, USA). qRT-PCR was carried out in experimental (at least) triplicates using Power SYBR Green reagent in AB7900HT 384-well system (Applied Biosystems, CA, USA). Primer sequence information is available upon request. Abundance was calculated by normalization to beta Actin (Actb).

Histology

Request a detailed protocol

Procedures for alkaline phosphatase staining for RorαCre;Z/AP brain, X-gal staining for RoraCre;R26R brain, and anterograde tracing with fluorescent cholera toxin subunit B (CTB) have been described previously (Hatori et al., 2008; Chou et al., 2009; Brown et al., 2010).

Locomotor activity measurement

Request a detailed protocol

Daily locomotor activity of mice individually housed in wheel running cages was measured following standard methodology (Siepka and Takahashi, 2005). Typically, 6- to 10-week-old mice kept in cages were placed inside light tight boxes with independent illumination. During the light phase, the mice received ∼150 lux of white light from fluorescent light source. Wheel running activity in 5 min bins was continuously collected and later analyzed by ClockLab software (Actimetrics, Evanston, IL, USA). All routine animal husbandry care during the dark phase was performed under dim red light.

Multiunit activity recording of the SCN slices

Request a detailed protocol

Mice were sacrificed by cervical dislocation followed by rapid dissection. The brains were cooled down in an ice-cold ACSF solution (125 mM NaCl, 25 mM KCl, 1 mM MgCl2, 1.25 mM NaH2PO4, 2 mM CaCl2, 20 mM Glucose, 26 mM NaHCO3, Penicillin 5000 IU/ml, and Streptomycin 5000 µg/ml) saturated with 95% O2/5% CO2. Coronal slices (∼300 μm) were then prepared using a tissue chopper and trimmed to ∼10 mm2 slices containing both nuclei. Finally, the slices were transferred to the multi electrodes arrays (MEA). The MEA consists of a glass recording chamber, on the bottom of which are engraved 256 electrodes (10 μm in diameter, situated every 60 μm) and arranged in a 16 × 16 square grid (Multichannel Systems, Reutlingen, Germany). The chamber was continuously perfused with heated (35°C) and oxygenated ACSF-containing antibiotics. Extracellular electrical activity was continuously monitored (signal was acquired from all 256 channels, 10 kHz) and spikes crossing a threshold set at 3 times the standard deviation of the noise on each channel were recorded and stored for off-line analysis. Just prior to placing the MEA on the amplifier, a bright field picture of the slice position on the electrode was rapidly taken to assess SCN placement. The electrodes covering the SCN were then continuously recorded from 2 to 4 days. Channels displaying noise or monotonically decreasing activity were excluded from subsequent analysis. We recorded from animals housed for at least 2 weeks either in DD or LD (12/12) conditions. Animals were handled under dim red light until the optic nerve was cut. Data analysis and display were performed using NeuroExplorer (Plexon Denton, TX), Oriana (Kovach Computing Services, UK), and custom software written in MATLAB (MathWorks, Natick, MA). Peaks of firing were determined after smoothing of data (rloess, MATLAB) and fitting with a sinusoidal function.

Daily administration of VIP was realized by switching the perfusion from the tank containing the medium described earlier to a tank containing the same medium supplemented with VIP (Calbiochem, EMD Millipore, MA) at the concentration of 25 nM for 1 hr.

Plasmid construction

Request a detailed protocol

The coding regions corresponding to full-length mouse and human Lhx1 were amplified by PCR from pineal cDNA and subcloned into pcDNA3.1-TOPO (Invitrogen, CA, USA) to yield expression plasmids mouse Lhx1/pcDNA3.1 and human Lhx1/pcDNA3.1, respectively. Site-directed mutagenesis (Stratagene, CA, USA) was performed to construct mouse Lhx1/pcDNA3.1 to generate a point mutation changing asparagine (amino acid 230) to serine. A DNA fragment corresponding to ∼1 kb of mouse Vip promoter was amplified by PCR from mouse genomic DNA and cloned into pGL3 basic vector (Promega) to yield the Vip reporter vector.

Transcriptional assay

Request a detailed protocol

293T cells were cultured in DMEM supplemented with 10% FBS. The cells (40,000 cells) in 96-well plates were transfected by using TransIT-LT1 (Mirus Bio, WI, USA) with various amounts of expression plasmid (total amount was adjusted to 250 ng by adding empty vector pcDNA3.1), and 5 ng of firefly luciferase reporter plasmid. The cell lysates were prepared 46 hr after the transfection and subjected to dual-luciferase assay by luminometry (Promega, WI, USA).

References

    1. Schwartz WJ
    2. Zimmerman P
    (1990)
    Circadian timekeeping in BALB/c and C57BL/6 inbred mouse strains
    The Journal of Neuroscience 10:3685–3694.
    1. Wee R
    2. Castrucci AM
    3. Provencio I
    4. Gan L
    5. Van Gelder RN
    (2002)
    Loss of photic entrainment and altered free-running circadian rhythms in math5-/- mice
    The Journal of Neuroscience 22:10427–10433.

Article and author information

Author details

  1. Megumi Hatori

    Regulatory Biology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Present address
    School of Medicine, Keio University, Tokyo, Japan
    Contribution
    MH, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Shubhroz Gill
    For correspondence
    mhatori@a6.keio.jp
    Competing interests
    The authors declare that no competing interests exist.
  2. Shubhroz Gill

    Regulatory Biology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Contribution
    SG, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Megumi Hatori
    Competing interests
    The authors declare that no competing interests exist.
  3. Ludovic S Mure

    Regulatory Biology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Contribution
    LSM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  4. Martyn Goulding

    Molecular Neurobiology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Contribution
    MG, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  5. Dennis D M O'Leary

    Molecular Neurobiology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Contribution
    DDMO'L, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  6. Satchidananda Panda

    Regulatory Biology Laboratory, Salk Institute for Biological Studies, La Jolla, United States
    Contribution
    SP, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    panda@salk.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Eye Institute (EY016807)

  • Satchidananda Panda

National Institute of Diabetes and Digestive and Kidney Diseases (DK091618)

  • Satchidananda Panda

National Institute of Neurological Disorders and Stroke (NS31558)

  • Dennis D M O'Leary

National Institute of Mental Health (MH086147)

  • Dennis D M O'Leary

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Sheena Keding and Hiep Le for technical help. This work was partially supported by fellowships from the Japan Society for the Promotion of Science, and Messinger Healthy Living to MH and from Fyssen and Catharina Foundation to LSM. HA and Mary K Chapman Foundation, and The Leona M and Harry B Helmsley Charitable Trust grant #2012-PG-MED002 supported SG. NIH grants NS31558 and MH086147 to DDMO, NS80586 to MG and NIH grant EY016807 and Hearst Foundation to SP funded the study. Funding from Glenn Foundation, P30 CA014195 and P30 EY019005 facilitated circadian activity and gene expression analyses.

Author contributions: MH, SG, LM performed experiments and data analyses, MH, DDMO, MG, and SP conceptualized the experiments, MH and SP wrote the manuscript. Correspondence and requests for materials should be addressed to satchin@salk.edu or mhatori@a6.keio.jp.

Ethics

Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (#12-00026) of the Salk Institute for Biological Studies. The protocol was approved by the IACUC committee of the Salk Institute. All surgery was performed under IACUC approved anesthesia, and every effort was made to minimize suffering.

Copyright

© 2014, Hatori et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,638
    views
  • 538
    downloads
  • 68
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Megumi Hatori
  2. Shubhroz Gill
  3. Ludovic S Mure
  4. Martyn Goulding
  5. Dennis D M O'Leary
  6. Satchidananda Panda
(2014)
Lhx1 maintains synchrony among circadian oscillator neurons of the SCN
eLife 3:e03357.
https://doi.org/10.7554/eLife.03357

Share this article

https://doi.org/10.7554/eLife.03357

Further reading

    1. Computational and Systems Biology
    2. Genetics and Genomics
    Sanjarbek Hudaiberdiev, Ivan Ovcharenko
    Research Article

    Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

    1. Developmental Biology
    2. Genetics and Genomics
    Ignacy Czajewski, Bijayalaxmi Swain ... Daan MF van Aalten
    Research Article

    O-GlcNAcylation is an essential intracellular protein modification mediated by O-GlcNAc transferase (OGT) and O-GlcNAcase (OGA). Recently, missense mutations in OGT have been linked to intellectual disability, indicating that this modification is important for the development and functioning of the nervous system. However, the processes that are most sensitive to perturbations in O-GlcNAcylation remain to be identified. Here, we uncover quantifiable phenotypes in the fruit fly Drosophila melanogaster carrying a patient-derived OGT mutation in the catalytic domain. Hypo-O-GlcNAcylation leads to defects in synaptogenesis and reduced sleep stability. Both these phenotypes can be partially rescued by genetically or chemically targeting OGA, suggesting that a balance of OGT/OGA activity is required for normal neuronal development and function.