- 1International Institute for Integrative Sleep Medicine (WPI-IIIS), University of Tsukuba, Tsukuba, Japan
- 2Department of Psychiatry & Neuroscience, Peter O’Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, United States
Non-rapid eye movement (NREM) sleep is tightly homeostatically regulated and essential for survival. In the electroencephalogram (EEG), oscillations in the delta (0.5–4 Hz) range are prominent during NREM sleep. These delta oscillations are, to date, the best indicator for homeostatic sleep regulation; they are increased after prolonged waking and fade during NREM sleep. The precise mechanisms underlying sleep homeostasis and the generation of EEG delta oscillations are still being investigated. Activity-dependent neuronal calcium influx has been hypothesized to play an important role in generating delta oscillations and might be involved in downstream signaling that mediates sleep function. Dihydropyridine blockers of L-type voltage-gated calcium channels (VGCCs) are in wide clinical use to treat hypertension and other cardiovascular disorders and are readily blood-brain-barrier penetrant. We therefore, wanted to investigate their potential effects on EEG delta oscillation and homeostatic NREM sleep regulation in freely behaving mice. In vivo two-photon imaging of cortical neurons showed larger spontaneous calcium transients in NREM sleep compared to waking. Application of the dihydropyridine calcium blocker nicardipine significantly reduced cortical calcium transients without affecting the generation of delta oscillations. Nicardipine also did not affect EEG delta oscillations over 24 h following application. The time spent in NREM sleep and NREM episode duration was also not affected. Thus, acute block of calcium entry through L-type VGCCs does not interfere with EEG delta oscillations or their homeostatic regulation, despite prior evidence from calcium channel knockout mice.
Introduction
In this study, we aimed to test two hypotheses. First, whether an activity-dependent influx of calcium through L-type voltage-gated calcium channels (VGCCs) in cortical neurons is necessary for the generation of the characteristic non-rapid eye movement (NREM) sleep electroencephalogram (EEG) delta oscillations (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018). Second, whether L-type VGCC activity is necessary for the gradual reduction of delta oscillations that occur during normal NREM sleep. Sleep is necessary for survival in all animals with a central nervous system and ensures its proper long-term function (Roth, 2009; Zielinski et al., 2016). Rapid eye movement (REM) and NREM sleep can be distinguished in birds and mammals. NREM sleep is characterized by high activity in the low-frequency EEG delta (0.5–4 Hz) band. These delta oscillations are caused by synchronous transitions of cortical neurons between a hyperpolarized, silent OFF state and a depolarized, active ON state, with each state lasting from a few tens to a few hundred milliseconds (Contreras and Steriade, 1995; Amzica and Steriade, 1998; Destexhe et al., 1999). Recent evidence points to a central role of intracellular calcium in cortical neurons in the generation of ON-OFF transitions (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018); according to this model, ON state action potential firing in cortical neurons triggers calcium influx, the rise in intracellular calcium concentration causes the opening of calcium-dependent potassium channels, inhibiting cortical neurons and triggering ON-OFF transitions; calcium extrusion during OFF state silence closes calcium-dependent potassium channels in cortical neurons and triggers OFF-ON transitions. The power of delta oscillations is, to date, the best marker for the homeostatic regulation of NREM sleep (Aeschbach and Borbely, 1993; Achermann et al., 1993; Borbély and Achermann, 1999; Dijk, 2009). High sleep need (e.g., after prolonged waking) is accompanied by high EEG delta power, and as sleep need is reduced during NREM sleep, delta oscillations fade (Aeschbach et al., 1997; Bjorness et al., 2016; Borbély et al., 2016). The mechanisms responsible for the generation of delta oscillations are still poorly understood. Cortical ON/OFF transitions can be observed in the absence of thalamocortical input (Steriade et al., 1993); however, in physiological NREM sleep, thalamic afferents may trigger and synchronize OFF-ON transitions (Gent et al., 2019). Moreover, it is still unknown whether delta oscillations reduce sleep need and exert the beneficial effects of NREM sleep on the brain or whether they are simply an epiphenomenon of the process.
In an effort to find more measures for successful sleep (i.e., sleep that reduces sleep pressure), several recent studies have investigated the effects of waking and sleep on the neural transcriptome, proteome, and phospho-proteome and have found clear patterns that correlate with sleep need (Wang et al., 2018; Brüning et al., 2019; Noya et al., 2019; Suzuki et al., 2020). However, no direct link has been established between neural activity, which generates delta oscillations, and the signaling that leads to the changes in gene expression, protein, and phospho-protein abundance observed after NREM sleep. An influx of calcium ions into neurons is one of the most established mechanisms linking activity to downstream cellular signals, and it has been hypothesized that this plays a role in slow wave generation and sleep (Llinás and Jahnsen, 1982; Deschênes et al., 1984; Steriade et al., 1985; Domich et al., 1986; Bal et al., 1995; Hughes et al., 2004; Lee et al., 2004; Anderson et al., 2005). A block of L-type VGCCs might therefore interfere with the signaling pathways linking neural activity to sleep need resolution. L-type VGCCs are widely expressed in cortical neurons (Hell et al., 1993) and are responsible for a postsynaptic influx of calcium, which can lead to synaptic plasticity and activity-induced regulation of gene expression (Catterall, 2000; Dolmetsch et al., 2001; Schafe and LeDoux, 2008). Dihydropyridines block L-type VGCCs and cause a significant reduction in activity-dependent calcium influx in neurons (Borst and Sakmann, 1996). They are widely used to treat cardiovascular disorders, mainly hypertension, and typically penetrate the blood-brain barrier (Ritz et al., 2010; He et al., 2018). Sleep disturbances are not frequently reported under dihydropyridine treatment and may be limited to patients suffering from obstructive sleep apnea (Nerbass et al., 2011). However, milder symptoms might have been overlooked. Therefore, we were motivated to investigate the effects of the dihydropyridine nicardipine on NREM sleep. We monitored in vivo calcium transients in cortical neurons and observed the duration and intensity of NREM sleep following intracerebroventricular (ICV) treatment in mice. Despite robust suppression of cortical calcium influx, dihydropyridine treatment did not significantly affect NREM sleep.
Materials and methods
Animals and surgeries
Animals
Experimental procedures were carried out in accordance with local and national regulations and after approval by the animal care and use committee of the University of Tsukuba. C57BL/6J male mice, 8–25 weeks old (Jackson Laboratory, Japan; N = 7), were used for sleep recordings. For in vivo calcium imaging, we crossed Ai148 (TITL2-GC6F-ICL-tTA2) mice (Allen Institute, WA, USA) with mice expressing Cre under the CaMKII promoter (Tsien et al., 1996), which led to the expression of GCaMP6f in forebrain excitatory neurons (N = 3). Mice were singly housed after surgery in a 12-h light and 12-h dark cycle with food and water ad libitum. Zeitgeber time 0 (ZT0) is light onset.
Electrode implantation
Mice were anesthetized with isoflurane (3–4% for induction and 2–2.5% for maintenance) and placed in a stereotaxic frame (David Kopf Instruments). The body temperature was maintained at 34–35°C with a heating pad. The scalp was removed, and the skull surface was cleaned with a scalpel blade. Two EEG screw electrodes were inserted into the skull for sleep recording. Two wire-electrodes (Cooner Wire, CA, USA) were implanted into the neck muscles bilaterally for electromyogram (EMG) recordings. Electrodes were connected to a standard 0.1-inch four-pin connector. This connector was fixed to the skull with cyanoacrylate and dental cement.
Microdrive implantation
Surgery followed with a previously published protocol (Matsumoto et al., 2020), briefly: tetrodes, EEG screws, and EMG wires were implanted using a custom-built microdrive. Mice were anesthetized with isoflurane; the surface of the brain was exposed by scalp incision followed by craniotomy. The dura was removed, and a microdrive system with six tetrodes (KANTHAL Precision Technology, nichrome, 14 μm in diameter) was implanted together with the EEG screws and EMG wires (Cooner Wire, CA, USA). Three tetrodes were placed in layer V of the primary motor cortex (M1) of the right and left hemispheres, and one tetrode into the ventral hippocampal commissure.
Intracerebroventricular cannula implantation
A single guide cannula of 27 gauge (RWD Life Science, CA, USA) was inserted with a tip distance of 1.9 mm from the surface of the brain at a 10° angle, targeting the lateral ventricle. During the ICV injection, injector cannulas of 33 gauge (RWD Life Science, CA, USA) were inserted through the guide cannula.
Optical window insertion
Optical windows were inserted into the skulls of three animals after EEG and EMG electrode implantation and ICV cannula insertion. A craniotomy was performed, covering the M1 and secondary somatosensory (S2) regions of the cortex using a dental drill. A hooked needle was used to remove a circular (2 mm) flap of the skull, leaving the dura intact. A Lipidure-CM5206 (NOF corporation) treated circular glass coverslip was held against the brain surface and glued to the skull with ultraviolet light hardening glue (BONDIC Starter Kit Complete; BONDIC, NY, USA). A custom-made metal head plate with a sufficient central opening was then secured to the skull with dental cement (SuperBond C&B set; Sun Medical, Shiga, Japan).
Pharmacology
Mice were injected ICV with a 10% sorbitol vehicle solution (Sigma-Aldrich, USA) dissolved in ultrapure H2O (MilliporeSigma, Sigma-Aldrich, Inc., USA) or with 100 nM of nicardipine hydrochloride (N7510, Sigma-Aldrich, USA) dissolved in vehicle solution. During the ICV injection, an injector cannula of 33G was inserted through the guide cannula to within 0.1 mm of the tip of the guide cannula. The solution was carried via polyethylene tubing (RWD Life Science, CA, USA) from a 1 ml syringe (Terumo, Japan).
Sleep recording and analysis
Recording of electrical signals
EEG and EMG electrodes were connected to a low-noise amplifier and digitizer board (Intan Technologies, CA, USA). Signals were recorded at a sampling rate of 1 kHz with the provided open-source software (Intan Technologies, USA). EEG and EMG were extracted, and band-pass filtered (EEG 0.1–500 Hz, EMG 10–500 Hz) using MATLAB (Matlab 2018a, The MathWorks, Inc., MA, USA) and provided import functions (Intan Technologies, CA, USA). After the electrode implantation, mice were given time to recover and acclimatize to recording chambers for 7 days, followed by 3 days of baseline recordings, 1 day of vehicle injection, followed by 1 day of nicardipine injection, with one rest day between injections.
For EEG and EMG recording during imaging, the signals were band-pass filtered and amplified using an analog amplifier (MEG-5200, NIHON KOHDEN, Japan), digitized at 2000 Hz (Digidata 1440A, Molecular Devices, USA), and acquired using Clampex 10.3 (Molecular Devices, USA). Each scan was digitized and recorded to match the imaging data.
The procedure for habituating the animals after tetrode implantation and obtaining neuronal spiking data, EEG, and EMG recordings followed published methods (Matsumoto et al., 2020).
Sleep deprivation
Mice were deprived of sleep by gentle handling and cage changes every hour for 4 h (Suzuki et al., 2013). Sleep deprivation started at ZT0 (onset of the light period). Mice were allowed to fall asleep undisturbed after sleep deprivation.
Sleep scoring
For sleep scoring, EEG and EMG signals were divided into regular-spaced epochs. Epoch length was set to 4 s for EEG data obtained during in vivo imaging because of the needed high resolution. Epochs of 10-s length were used for all other sleep analyses. Each epoch was scored as either wake, NREM, or REM sleep. EEG signals underwent fast Fourier transformation by custom MATLAB scripts (MATLAB 2018a, The MathWorks, Inc., MA, USA). Scoring was based on the power in the EEG delta band (0.5–4 Hz), the EEG theta band (6–10 Hz), and the ratio of delta to theta band power, as well as the integral of the EMG signal. Epochs with a high integral of EMG signal and low delta power were scored as wake. Epochs with a low integral of EMG signal and high delta power were scored as NREM. Epochs with a very low integral of EMG signal and a high theta to delta power ratio were scored as REM. The scoring was conducted using custom MATLAB programs (MATLAB 2018a, The MathWorks, Inc., MA, USA) (Bjorness et al., 2016).
Spike sorting
Tetrode signals were recorded with a Neuralynx Digital Lynx system using Neuralynx Cheetah software (Neuralynx, MT, USA) (Matsumoto et al., 2020). Spikes were detected online, and the multi-unit activity was sorted into single-unit activity offline with Spike Sort 3D software (Neuralynx, MT, USA). Cluster classification was performed using previous methods (Matsumoto et al., 2020). Briefly, clusters were considered to represent single units if the isolation distance was ≥ 20 and L-ratio ≤ 0.3. Only units with a peak-to-peak width greater than 250 μs were classified as putative excitatory neurons and were used for this study.
In vivo two-photon imaging and analysis
After the optical window surgery, mice were given 7 days to recover/habituate to the head plate. For imaging, they were placed on an air-suspended spherical treadmill. The mice could undergo spontaneous wake and NREM transitions on the spherical treadmill during imaging. Due to the short recording duration of 2 h or less, we did not reliably observe REM episodes. The calcium activity was observed with a two-photon microscope (Axio Examiner Z1, Zeiss) (Miyazaki et al., 2020). Imaging of GCaMP6f expressing neurons was conducted in cortical layer II/III, with excitation at 910 nm with a Ti: Sa laser (Mai Tai DeepSee, Spectra-Physics). For each session, 1,400 frames were acquired at a rate of 440 ms per scan, which took about 11 min. In total, there were 13 sessions per animal: three sessions for baseline, three sessions after vehicle injection, and seven sessions after nicardipine injection.
The scanned images were loaded into a custom-made MATLAB program (MATLAB 2018a, The MathWorks, Inc., MA, USA), where they were analyzed semi-automatically. The extracted raw data were motion corrected (Pnevmatikakis and Giovannucci, 2017). Regions of interest (ROI) were manually drawn to match the somata of the neurons. Their calcium activity was detected as a change in fluorescence over baseline fluorescence traces (dF/F), following this function: dF/F = (Ft–F0)/F0, where F0 is the baseline mean of the four frames previous to the Ft frame. The integral of dF/F over the 0.5 thresholds was calculated to compare the calcium activity among different conditions. To determine the threshold, we calculated the average standard deviation of the dF/F signal under nicardipine for all neurons when the relative contribution of calcium transients to noise was the lowest. The value of 0.5 represents two standard deviations above the mean (rounded from a value of 0.502). Thus, the thresholding removes the background noise and shows only the difference in calcium activity during baseline, vehicle, and nicardipine sessions. The temporal resolution of our system (440 ms per scan) and the decay time constant of the indicator are insufficient (Chen et al., 2013; Umpierre et al., 2020; Huang et al., 2021; Zhang et al., 2021) for a proper analysis of underlying spiking patterns. The average of the above 0.5 dF/F integrals of the three baseline scans was set to 100%. The integral obtained during scans following vehicle injection and nicardipine injection was normalized to that value. Scans were obtained 40 min after vehicle injection and 100 min after nicardipine injection.
The statistical analyses were performed in MATLAB (MATLAB 2018a, The MathWorks, Inc., MA, USA) and IBM SPSS Statistics (SPSS for Windows, Version 25.0. Armonk, NY: IBM Corp.) for two-way repeated measures ANOVA. Statistical significance was defined as p < 0.05. Results are presented as mean ± standard error of the mean (SEM) unless otherwise noted.
Results
In this study, we could assess in vivo calcium dynamics together with EEG/EMG recordings in mice freely transitioning between vigilance states. These measures provide optimal quantitative measures for the potential role of L-type VGCC-mediated neuronal calcium influx in sleep function.
Calcium influx into excitatory neurons is larger in non-rapid eye movement sleep compared to waking
Burst firing is an effective way to produce large intracellular calcium transients. In shorter recordings (Ohyama et al., 2020), we have shown that the propensity for burst firing is larger in NREM sleep compared to waking. In 24 h recordings in which the animals naturally transitioned between vigilance states, we again found that burst frequency in NREM was significantly larger than in waking (NREM = 0.362 ± 0.14 min–1; waking = 0.295 ± 0.12 min–1, ten neurons, t = -3.502, p < 0.007, paired, two-tailed T-test). This is compatible with larger calcium transients in NREM sleep; we observed intracellular calcium transients directly to further explore this. Using in vivo two-photon imaging of the genetically encoded calcium sensor GCaMP6f, we observed calcium signals in naturally occurring wake and NREM sleep in head-fixed animals that could freely walk on air-suspended Styrofoam balls. Calcium influx into excitatory neurons was larger during NREM sleep compared to waking (Figure 1).
Figure 1. Cortical calcium transients during waking and NREM sleep. Somatic calcium transients are larger in NREM sleep compared to waking in simultaneous EEG and EMG recordings and in vivo two-photon GCaMP6f calcium imaging. (A) Sample of one recording session. Top trace: hypnogram derived from the simultaneously recorded EEG and EMG. Middle graph: Pseudo-colored spectrogram of the EEG recording. The green line below indicates the time during which the imaging data was obtained (SCAN). The graph below shows power in the EEG delta band (0.5–4 Hz) over time. (B) Enlarged time window from panel (A) (indicated by the dotted lines), representing one imaging episode of 1,400 frames lasting 11 min. Top trace: Hypnogram, times in NREM sleep are marked with blue bars, showing two long wake episodes. Blue trace: EEG delta power in 4 s epochs. Black traces: sample of dF/F traces from five different regions of interest (ROI), drawn over somata of layer II/III excitatory neurons. Notice the decrease in activity during the wake episodes. (C) Average data from a total of 100 neurons in 3 mice are shown (mean ± SEM). The dF/F traces were thresholded at 0.5 for all ROIs (see section “Materials and methods”), and the integral over the threshold was calculated. Given the calcium fluorescence variability between ROIs, we calculated the activity in its percentage during NREM and its percentage during wakefulness, corrected by the relative durations of NREM and wakefulness. P values less than 0.001 are indicated with three asterisks.
The higher calcium transients in NREM sleep are compatible with previously published results (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018), indicating that intracellular calcium transients contribute to the characteristic ON/OFF oscillation in NREM sleep. We, therefore, wanted to test the effects of acutely perturbing calcium influx into cortical neurons and measure the effect on NREM sleep.
Intracerebroventricular nicardipine reduces calcium influx into excitatory neurons
Therefore, we next investigated the inhibition of calcium activity under the influence of nicardipine, a dihydropyridine L-type VGCC blocker, using in vivo two-photon calcium imaging of GCaMP6f labeled neurons. The calcium activity of individual neurons decreased significantly 100 min after 100 nM ICV nicardipine injection (Figure 2). We calculated the integral over time of the dF/F signals above a 0.5 dF/F threshold across three animals. As shown in Figure 2B, this integral was significantly reduced under nicardipine, indicating that it successfully blocked calcium entry through L-type VGCCs (Wu et al., 1999). As expected (Wu et al., 1999), nicardipine did not eliminate calcium transients (Figure 2C), indicating that not all calcium influx was blocked.
Figure 2. Intracerebroventricular application of Nicardipine reduces cortical somatic calcium transients. Layer II/III neuron somatic calcium transients in in vivo two-photon imaging of GCaMP6f under baseline conditions, after ICV vehicle injection, and after ICV nicardipine injection. (A) Representative grayscale image of a field of view from one experiment. (B) Comparison of calcium signals under baseline conditions (black), after 40 min of ICV vehicle injection (blue), and after 100 min of ICV nicardipine injection. A threshold of 0.5 was applied to the dF/F signals (see section “Materials and methods”), which were then separately integrated over time for the three conditions. (N = 3 animals; 63 neurons, 24 neurons, and 13 neurons), error bars are SEM. P values less than 0.001 are indicated with three asterisks. (C) Representative dF/F traces over time of two neurons in the three conditions (black, baseline; blue after vehicle injection; red after nicardipine injection).
We have thus shown that activity-dependent calcium influx into cortical neurons is significantly reduced by ICV nicardipine application. We next investigated the effects of such reduced calcium influx on NREM sleep, including the generation of EEG delta waves, and on the efficiency with which the need for sleep was resolved, as indicated by the reduction of EEG delta power during NREM sleep.
Intracerebroventricular nicardipine has no significant effect on non-rapid eye movement sleep
To observe an impairment in NREM sleep, we administered vehicle or nicardipine via ICV injection under three different sleep conditions and at different zeitgeber times (ZT). First, at ZT0, when mice are at their normal peak sleep need. Second, at ZT12, when sleep need is typically lowest, and third, after four h of sleep deprivation starting at ZT0, resulting in very high sleep need at ZT4.
There was no difference in NREM EEG delta (0.5–4 Hz) energy between the two treatment conditions, nicardipine, compared to control injections, both immediately after the injection and a few hours later, when the nicardipine effect had dissipated (Figure 3A).
Figure 3. Application of Nicardipine ICV does not Affect EEG Delta Oscillations. The power of EEG oscillations in the delta (0.5–4 Hz) band during NREM sleep did not differ when the animals received ICV vehicle injections or ICV nicardipine injections. (A) Power in the EEG delta band during NREM sleep is integrated over time to obtain delta energy. Delta energy production was measured for 3 h in three different conditions [at ZT0, at ZT12, and at ZT4 after 4 h of sleep deprivation (SlDep)]. Each animal was measured after vehicle injection (blue) and after nicardipine injection (red). Delta energy is plotted over time in 0.5-h bins. Nicardipine did not affect EEG delta energy in any of the three conditions at any time point. (B) The effectiveness of NREM sleep episodes can be measured by the reduction in EEG delta power over the duration of the episode. We analyzed NREM sleep episodes longer than 5 min in the recordings from ZT0, ZT12, and ZT4 after sleep deprivation. Peak EEG delta power within the first 3 min is plotted on the left-hand side, and the last epoch’s delta power is plotted on the right-hand side. After both vehicle treatment (blue) and nicardipine treatment (red), EEG delta power was lower at the end of the episode compared to the beginning. (C) Box plots (orange line, mean; blue box 25th percentile; black whiskers 90th percentile) of comparison between vehicle treatment and nicardipine treatment on the average reduction in delta power during NREM sleep episodes longer than 5 min. In all three conditions, there was no significant difference between vehicle and nicardipine. P-values were obtained with paired two-tailed Student’s t-test (N = 7).
Therefore, there was no correlation between the generation of NREM sleep EEG delta oscillations (Supplementary material Appendix 1 and Supplementary Figure 1) and the drastic reduction of activity-dependent calcium influx in cortical neurons. Moreover, sleep-need reduction during NREM sleep episodes (Bjorness et al., 2016) was not affected by the reduced calcium influx (Figures 3B,C and Table 1). The reduction in the power of EEG delta oscillations during NREM episodes was not different between treatment conditions (Figure 3C).
We also investigated the sleep patterns of the mice under different treatment conditions by measuring the time spent in wake, NREM sleep, and REM sleep. We did not observe any significant difference between nicardipine and vehicle-treated conditions in the overall durations of wake, NREM sleep, and REM sleep (Table 2 and Supplementary Figure 2).
In summary, ICV injection of nicardipine drastically reduced activity-dependent calcium influx into cortical neurons but did not affect the generation of NREM EEG delta oscillations or the dissipation of sleep need during NREM sleep.
Discussion
In this study, we demonstrated a strong reduction of calcium influx into cortical neurons after applying the L-type VGCC blocker nicardipine in vivo. This reduced calcium influx did not interfere with the generation of EEG delta waves and did not block the gradual reduction of EEG delta power during NREM sleep. NREM EEG delta power has been studied as a strong indicator of sleep need (Borbély and Achermann, 1999; Dijk, 2009). However, the link between sleep need and the cortical ON/OFF oscillation that underlies the generation of delta oscillations and cortical activity is not well elucidated. One compelling argument posits that calcium activity is involved in generating delta oscillations (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018). Thus, we sought to investigate the role of activity-dependent calcium influx through L-type VGCCs using pharmacological block by nicardipine. L-type VGCCs are widely expressed in cortical neurons, and 50% of action potential-driven postsynaptic calcium influx is mediated by L-type VGCCs (Wu et al., 1999). A dihydropyridine blocker was selected since they are a widely used class of drugs that readily penetrate the blood-brain barrier and might have side effects that have been overlooked so far.
In agreement with previous results, we found higher incidences of burst firing (Hubel, 1959; Hobson and McCarley, 1971; Ohyama et al., 2020) and larger calcium transients (Ohyama et al., 2020) in cortical neurons in NREM sleep compared to waking. Tetrode recordings lasting 24 h show that increased bursting is a stable characteristic of NREM sleep. Some groups have reported lower calcium activity in cortical neurons in NREM sleep compared to waking (Niethard et al., 2016, 2018). However, a complex picture emerges with different cell types and compartments showing different changes in calcium activity between waking and NREM sleep (Sigl-Glöckner and Seibt, 2019). At this stage, we cannot fully resolve the differences between our findings and some prior results, but previous studies have found increased calcium activity in subtypes of excitatory cortical neurons (Niethard et al., 2016, 2021).
It was recently proposed that increased burst firing and the related calcium influx were fundamental for generating ON/OFF state oscillations and NREM sleep (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018). Our results do not support an absolute need for activity-driven calcium influx in the generation of EEG delta waves because, as aforementioned, a large percentage of the calcium influx is through L-type VGCCs. Alternative sources of intracellular calcium might still play a role. Other VGCCs such as N, P/Q or T-type, NMDA receptors, or phospholipase C signaling are other possible sources for activity-dependent calcium entry. Calcium influx through T-type channels was shown to characterize high-frequency bursts of action potentials during NREM sleep and to initiate UP states (Crunelli et al., 2014). Because of their presynaptic physiology—coupling neuronal excitation to the secretion of neurotransmitters (Fletcher et al., 2001; Wormuth et al., 2016), investigation of other VGCCs remains challenging. Several studies have shown that systematic application of the NMDA receptor blocker ketamine results in particularly intense EEG delta oscillations (Fontanini et al., 2003; Chauvette et al., 2011; Ahnaou et al., 2017), arguing against disruption of ON/OFF oscillations and delta waves by blocking neuronal calcium influx. Thus, an obligatory link between EEG delta oscillations and calcium influx is unlikely. Further investigations into the calcium activities of different receptors and signaling processes during the ON/OFF phases are necessary.
A more fundamental and less addressed issue concerns the potential down streams signals generated by delta wave activity mediating sleep need resolution. Does neural activity underlying delta waves contribute to sleep need reduction? If so, how? In our study, L-type VGCC blockade produced a substantial reduction of calcium transients; nevertheless, the decay of delta wave activity was unchanged from the control. Thus, calcium influx does not link neuronal activity and sleep need resolution. However, we would like to emphasize the need to investigate other potential links between EEG delta power and sleep need resolution.
This might be achieved by the recent advances in the proteomic and transcriptomic investigation of sleep, which provide us with good markers for sleep need buildup and resolution (Bjorness et al., 2009, 2016; Wang et al., 2018; Suzuki et al., 2020). For example, wake-expressed c-fos levels rapidly drop during subsequent NREM sleep (Basheer et al., 1999), and reduced levels of immediate early genes may constitute a signaling pathway linking NREM delta power to changes in gene expression patterns. In the future, these markers can be used to determine whether NREM delta power is causally related to sleep need resolution or is an epiphenomenon of the process.
Our results show that neuronal firing and calcium activity increase during NREM sleep compared to wakefulness. Acute inhibition of calcium influx through L-type VGCC did not affect EEG delta power or sleep need resolution. This result agrees with the clinical literature on dihydropyridines, which readily cross the blood-brain barrier. Clinical studies showed no patients complaining of side effects in terms of mood and sleep (Leonetti and Salvetti, 1994; He et al., 2018). A recent clinical study showed that hypertensive patients receiving the dihydropyridine amlodipine less frequently suffered from insomnia than when they received other medications (Isaveya et al., 2020). In light of the proposed importance of calcium transients in NREM sleep and NREM sleep delta power (Tatsuki et al., 2016, 2017; Shi and Ueda, 2018), the lack of significant effects on NREM sleep of even high doses of dihydropyridines is reassuring for the large number of patients who take these substances as a daily medication.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by the animal care and use committee of the University of Tsukuba.
Author contributions
RG and KV contributed to the conception and design of the study. JD performed the statistical analysis. SM contributed to the database. GH performed the experiments and wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This research was supported by AMED under Grant Number JP21zf0127005.
Acknowledgments
We would like to thank the Japan Science and Technology Agency, Core Research for Evolutional Science and Technology (A3A28043), for supporting KV.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.969712/full#supplementary-material
Supplementary Figure 1. Sample EEG Delta Power and Hypnogram Plots over 24 h for Vehicle and Nicardipine Treatment. Power in the EEG delta (0.5–4 Hz) band (top) and hypnogram (bottom) from one mouse plotted against time for two non-consecutive 24 h recordings. Light (inactive) phase from ZT0 to ZT12 and dark (active) phase from ZT12 to ZT24. The graph on the left is the vehicle condition (indicated by the blue vertical bar), and the graph on the right is the nicardipine condition (indicated by the yellow vertical bar). (A) Injection of vehicle and nicardipine at ZT0. (B) Injection of vehicle and nicardipine at ZT 12. (C) Injection of vehicle and nicardipine at ZT4 after 4 h of sleep deprivation from ZT0 to ZT4.
Supplementary Figure 2. Vigilance State Plots at 1-h Intervals for Vehicle and Nicardipine Treatment. Plots of vigilance state duration in minutes per 1 h bin (top, wake; middle, NREM sleep; bottom, REM sleep) over 24 h. Data is mean ± SEM for N = 7 animals. ZT0–ZT12 is the light (inactive) phase, and ZT12–ZT24 is the dark (active) phase. Data from ICV vehicle treatment in orange, from ICV nicardipine treatment in blue. Treatment is indicated by the gray vertical bar. (A) Treatment time ZT0. (B) Treatment time ZT12. (C) Treatment time ZT4 after sleep deprivation from ZT0 to ZT4. P-values on top of graphs are from two-way repeated measures ANOVA for treatment effect. There was no significant difference between treatment conditions for any vigilance state.
References
Achermann, P., Dick, D. J., Brunner, D. P., and Borbely, A. (1993). A model of human sleep homeostasis based on EEG slow-wave activity: quantitative comparison of data and simulations. Brain Res. Bull. 31, 97–113. doi: 10.1016/0361-9230(93)90016-5
Aeschbach, D., and Borbely, A. A. (1993). All-night dynamics of the human sleep EEG. J. Sleep Res. 2, 70–81.
Aeschbach, D., Dijk, D.-J., and Borbély, A. A. (1997). Dynamics of EEG spindle frequency activity during extended sleep in humans: relationship to slow-wave activity and time of day. Brain Res. 748, 131–136. doi: 10.1016/s0006-8993(96)01275-9
Ahnaou, A., Huysman, H., Biermans, R., Manyakov, N. V., and Drinkenburg, W. H. I. M. (2017). Ketamine: differential neurophysiological dynamics in functional networks in the rat brain. Transl. Psychiatry 7:e1237.
Amzica, F., and Steriade, M. (1998). Electrophysiological correlates of sleep delta waves. Electroencephalogr. Clin. Neurophysiol. 107, 69–83. doi: 10.1016/S0013-4694(98)00051-0
Anderson, M. P., Mochizuki, T., Xie, J., Fischler, W., Manger, J. P., Talley, E. M., et al. (2005). Thalamic Cav3.1 T-type Ca2+ channel plays a crucial role in stabilizing sleep. Proc. Natl. Acad. Sci. U.S.A. 102, 1743–1748. doi: 10.1073/pnas.0409644102
Bal, T., von Krosigk, M., and McCormick, D. A. (1995). Role of the ferret perigeniculate nucleus in the generation of synchronized oscillations in vitro. J. Physiol. 483(Pt 3), 665–685. doi: 10.1113/jphysiol.1995.sp020613
Basheer, R., Sherin, J., Saper, C., Morgan, J. I., McCarley, R. W., and Shiromani, P. J. (1999). Effects of sleep on wake-induced c-fos expression. J. Neurosci. 17, 9746–9750. doi: 10.1523/JNEUROSCI.17-24-09746.1997
Bjorness, T., Kelly, C., Gao, T., Poffenberger, V., and Greene, R. W. (2009). Control and function of the homeostatic sleep response by adenosine A1 receptors. J. Neurosci. 29, 1267–1276. doi: 10.1523/JNEUROSCI.2942-08.2009
Bjorness, T. E., Dale, N., Mettlach, G., Sonneborn, A., Sahin, B., Fienberg, A. A., et al. (2016). An adenosine-mediated glial-neuronal circuit for homeostatic sleep. J. Neurosci. 36, 3709–3721. doi: 10.1523/JNEUROSCI.3906-15.2016
Borbély, A. A., and Achermann, P. (1999). Sleep homeostasis and models of sleep regulation. J. Biol. Rhythms 14, 557–568. doi: 10.1177/074873099129000894
Borbély, A. A., Daan, S., Wirz-Justice, A., and Deboer, T. (2016). The two-process model of sleep regulation: a reappraisal. J. Sleep Res. 25, 131–143. doi: 10.1111/jsr.12371
Borst, J., and Sakmann, B. (1996). Calcium influx and transmitter release in a fast CNS synapse. Nature 383, 431–434. doi: 10.1038/383431a0
Brüning, F., Noya, S. B., Bange, T. J., Koutsouli, S., Rudolph, J. D., Tyagarajan, S. K., et al. (2019). Sleep-wake cycles drive daily dynamics of synaptic phosphorylation. Science 366:eaav3617.
Catterall, W. A. (2000). Structure and regulation of voltage-gated Ca2+ channels. Ann. Rev. Cell Dev. Biol. 16, 521–555. doi: 10.1146/annurev.cellbio.16.1.521
Chauvette, S., Crochet, S., Volgushev, M., and Timofeev, I. (2011). Properties of slow oscillation during slow-wave sleep and anesthesia in cats. J. Neurosci. 19, 14998–15008.
Chen, T. W., Wardill, T. J., Sun, Y., Pulver, S. R., Renninger, S. L., Baohan, A., et al. (2013). Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499, 295–300. doi: 10.1038/nature12354
Contreras, D., and Steriade, M. (1995). Cellular basis of EEG slow rhythms: a study of dynamic corticothalamic relationships. J. Neurosci. 15, 604–622. doi: 10.1523/jneurosci.15-01-00604.1995
Crunelli, V., David, F., Leresche, N., and Lambert, R. C. (2014). Role for T-type Ca2+ channels in sleep waves. Pflugers Arch. Eur. J. Physiol. 466, 735–745. doi: 10.1007/s00424-014-1477-3
Deschênes, M., Paradis, M., Roy, J. P., and Steriade, M. (1984). Electrophysiology of neurons of lateral thalamic nuclei in cat: resting properties and burst discharges. J. Neurophysiol. 51, 1196–1219. doi: 10.1152/jn.1984.51.6.1196
Destexhe, A., Contreras, D., and Steriade, M. (1999). Spatiotemporal analysis of local field potentials and unit discharges in cat cerebral cortex during natural wake and sleep states. J. Neurosci. 19, 4595–4608. doi: 10.1523/jneurosci.19-11-04595.1999
Dijk, D. J. (2009). Regulation and functional correlates of slow wave sleep. J. Clin. Sleep Med. 5(2 Suppl.), S6–S15.
Dolmetsch, R. E., Pavani, U., Fife, K., Spotts, J. M., and Greenberg, M. E. (2001). Signaling to the nucleus by an L-type calcium channel-calmodulin complex through the MAP kinase pathway. Science 294, 333–339.
Domich, L., Oakson, G., and Steriade, M. (1986). Thalamic burst patterns in the naturally sleeping cat: a comparison between cortically projecting and reticularis neurones. J. Physiol. 379, 429–449. doi: 10.1113/jphysiol.1986.sp016262
Fletcher, C. F., Tottene, A., Lennon, V. A., Wilson, S. M., Dubel, S. J., Paylor, R., et al. (2001). Dystonia and cerebellar atrophy in Cacna1a null mice lacking P/Q calcium channel activity. FASEB J. 15, 1288–1290. doi: 10.1096/fj.00-0562fje
Fontanini, A., Spano, P. F., and Bower, J. (2003). Ketamine-xylazine-induced slow (<1.5Hz) oscillations in the rat piriform (olfactory) cortex are functionally correlated with respiration. J. Neurosci. 23, 7993–8001.
Gent, T. C., Bandarabadi, M., Herrera, C. G., and Adamantidis, A. R. (2019). Thalamic dual control of sleep and wakefulness. Nat. Neurosci. 21, 974–984. doi: 10.1038/s41593-018-0164-7
He, Q., Liu, J., Liang, J., Liu, X., Li, W., Liu, Z., et al. (2018). Towards improvements for penetrating the blood-brain barrier—recent progress fro a material and pharmaceutical perspective. Cells 7:24. doi: 10.3390/cells7040024
Hell, J. W., Westenbroek, R. E., Warner, C., Ahlijanian, M. K., Prystay, W., Gilbert, M. M., et al. (1993). Identification and differential subcellular localization of the neuronal class C and class D L-type calcium channel α1 subunits. J. Cell Biol. 123, 949–962. doi: 10.1083/jcb.123.4.949
Hobson, J. A., and McCarley, R. W. (1971). Cortical unit activity in sleep and waking. Electrophysiol. Clin. Neurophys. 30, 97–112.
Huang, L., Ledochowitsch, P., Knoblich, U., Lecoq, J., Murphy, G. J., Reid, R. C., et al. (2021). Relationship between simultaneously recorded spiking activity and fluorescence signal in GCaMP6 transgenic mice. eLife 10:e51675.
Hubel, D. H. (1959). Single unit activity in striate cortex of unrestrained cats. J. Physiol. 147, 26–38.
Hughes, S. W., Lörincz, M., Cope, D. W., Blethyn, K. L., Kékesi, K. A., Parri, H. R., et al. (2004). Synchronized oscillations at alpha and theta frequencies in the lateral geniculate nucleus. Neuron 42, 253–268. doi: 10.1016/s0896-6273(04)00191-6
Isaveya, G., Buriakovska, O., and Shalimova, A. (2020). Effect of amlodipine on the manifestations of chronic insomnia in hypertensive patients with type 2 diabetes mellitus. J. Polish Soc. Hypertens. 24, 192–199.
Lee, J., Kim, D., and Shin, H. S. (2004). Lack of delta waves and sleep disturbances during non-rapid eye movement sleep in mice lacking alpha1G-subunit of T-type calcium channels. Proc. Natl. Acad. Sci. U.S.A. 101, 18195–18199. doi: 10.1073/pnas.0408089101
Leonetti, G., and Salvetti, A. (1994). Effects of cliazapril and nitrendipine on blood pressure, mood, sleep, and cognitive function in elderly hypertensive patients: an Italian multicenter study. J. Cardiovasc. Pharmacol. 24, (Suppl. 3), S73–S77.
Llinás, R., and Jahnsen, H. (1982). Electrophysiology of mammalian thalamic neurones in vitro. Nature 297, 406–408. doi: 10.1038/297406a0
Matsumoto, S., Ohyama, K., Díaz, J., Yanagisawa, M., Greene, R. W., and Vogt, K. E. (2020). Enhanced cortical responsiveness during natural sleep in freely behaving mice. Sci. Rep. 10:2278. doi: 10.1038/s41598-020-59151-8
Miyazaki, T., Kanda, T., Tsujino, N., Ishii, R., Nakatsuka, D., Kizuka, M., et al. (2020). Dynamics of cortical local connectivity during sleep–wake states and the homeostatic process. Cerebr. Cortex 30, 3977–3990. doi: 10.1093/cercor/bhaa012
Nerbass, F. B., Pedrosa, R. P., Genta, P. R., Drager, L. F., and Lorenzi-Filho, G. (2011). Calcium channel blockers are independently associated with short sleep duration in hypertensive patients with obstructive sleep apnea. J. Hypertens. 29, 1236–1241. doi: 10.1097/HJH.0b013e3283462e8b
Niethard, N., Brodt, S., and Born, J. (2021). Cell-type-specific dynamics of calcium activity in cortical circuits over the course of slow-wave sleep and rapid eye movement sleep. J. Neurosci. 41, 4212–4222. doi: 10.1523/JNEUROSCI.1957-20.2021
Niethard, N., Hasegawa, M., Itokazu, T., Oyanedel, C. N., Born, J., and Sato, T. R. (2016). Sleep-stage-specific regulation of cortical excitation and inhibition. Curr. Biol. 26, 2739–2749. doi: 10.1016/j.cub.2016.08.035
Niethard, N., Ngo, H. V., Ehrlich, I., and Born, J. (2018). Cortical circuit activity underlying sleep slow oscillations and spindles. Proc. Natl. Acad. Sci. U.S.A. 115, E9220–E9229. doi: 10.1073/pnas.1805517115
Noya, S. B., Colameo, D., Brüning, F., Spinnler, A., Mircsof, D., Opitz, L., et al. (2019). The forebrain synaptic transcriptome is organized by clocks but its proteome is driven by sleep. Science 366:eaav2642. doi: 10.1126/science.aav2642
Ohyama, K., Kanda, T., Miyazaki, T., Tsujino, N., Ishii, R., Ishikawa, Y., et al. (2020). Structure of cortical network activity across natural wake and sleep states in mice. PLoS One 15:e0233561. doi: 10.1371/journal.pone.0233561
Pnevmatikakis, E. A., and Giovannucci, A. (2017). NoRMCorre: an online algorithm for piecewise rigid motion correction of calcium imaging data. J. Neurosci. 291, 83–94.
Ritz, B., Rhodes, S., Qian, L., Schernhammer, E., Olsen, J. H., and Friis, S. (2010). L-type calcium channel blocker snad Parkinson’s disease in Denmark. Ann. Neurol. 67, 600–606. doi: 10.1002/ana.21937
Schafe, G. E., and LeDoux, J. E. (2008). Learning and memory: a comprehensive reference; 4.11 Neural and molecular mechanisms of fear memory. Acad. Press 4, 157–192. doi: 10.1016/B978-012370509-9.00045-0
Shi, S., and Ueda, H. R. (2018). Ca2+ -dependent hyperpolarization pathways in sleep homeostasis and mental disorders. Bioessays 40:1700105. doi: 10.1002/bies.201700105
Sigl-Glöckner, J., and Seibt, J. (2019). Peeking into the sleeping brain: using in vivo imaging in rodents to understand the relationship between sleep and cognition. J. Neurosci. Methods 316, 71–82. doi: 10.1016/j.jneumeth.2018.09.011
Steriade, M., Deschênes, M., Domich, L., and Mulle, C. (1985). Abolition of spindle oscillations in thalamic neurons disconnected from nucleus reticularis thalami. J. Neurophysiol. 54, 1473–1497. doi: 10.1152/jn.1985.54.6.1473
Steriade, M., Nunez, A., and Amzica, F. (1993). Intracellular analysis of relations between the slow (< 1 Hz) neocortical oscillation and other sleep rhythms of the electroencephalogram. J. Neurosci. 13, 3266–3283. doi: 10.1523/JNEUROSCI.13-08-03266.1993
Suzuki, A., Sinton, C. M., Greene, R. W., and Yanagisawa, M. (2013). Behavioral and biochemical dissociation of arousal and homeostatic sleep need influenced by prior wakeful experience in mice. PNAS 110, 10288–10293.
Suzuki, A., Yanagisawa, M., and Greene, R. (2020). Loss of Arc attenuates the behavioral and molecular responses for sleep homeostasis in mice. PNAS 117, 10547–10553. doi: 10.1073/pnas.1906840117
Tatsuki, F., Ode, K., and Ueda, H. (2017). Ca2+-dependent hyperpolarization hypothesis for mammalian sleep. Neurosci. Res. 118, 48–55. doi: 10.1016/j.neures.2017.03.012
Tatsuki, F., Sunagawa, G. A., Shi, S., Susaki, E. A., Yukinaga, H., Perrin, D., et al. (2016). Involvement of Ca(2+)-dependent hyperpolarization in sleep duration in mammals. Neuron 90, 70–85. doi: 10.1016/j.neuron.2016.02.032
Tsien, J. Z., Chen, D. F., Gerber, D., Tom, C., Mercer, E. H., Anderson, D. J., et al. (1996). Subregion- and cell type–restricted gene knockout in mouse brain. Cell 87, 1317–1326.
Umpierre, A. D., Bystrom, L. L., Ying, Y., Liu, Y. U., Worrell, G., and Wu, L. J. (2020). Microglial calcium signaling is attuned to neuronal activity in awake mice. eLife 9:e56502.
Wang, Z., Ma, J., Miyoshi, C., Li, Y., Sato, M., Ogawa, Y., et al. (2018). Quantitative phophoproteomic analysis of the molecular substrates of sleep need. Nature 558, 435–439. doi: 10.1038/s41586-018-0218-8
Wormuth, C., Lundt, A., Henseler, C., Müller, R., Broich, K., Papazoglou, A., et al. (2016). Review: Cav2.3 R-type voltage-gated Ca2+ channels - functional implications in convulsive and non-convulsive seizure activity. Open Neurol. J. 10, 99–126. doi: 10.2174/1874205X01610010099
Wu, L. G., Westenbroek, R. E., Borst, J. G., Catterall, W. A., and Sakmann, B. (1999). Calcium channel types with distinct presynaptic localization couple differentially to transmitter release in single calyx-type synapses. J. Neurosci. 19, 726–736.
Zhang, Y., Rózsa, M., Liang, Y., Bushey, D., Wei, Z., Zheng, J., et al. (2021). Fast and sensitive GCaMP calcium indicators for imaging neural populations. bioRxiv [Preprint]. doi: 10.1101/2021.11.08.467793v2
Keywords: non-rapid eye movement wave sleep, delta activity, NREM, voltage gated calcium channel, dihydropyridine, calcium imaging
Citation: Han G, Matsumoto S, Diaz J, Greene RW and Vogt KE (2022) Dihydropyridine calcium blockers do not interfere with non-rapid eye movement sleep. Front. Neurosci. 16:969712. doi: 10.3389/fnins.2022.969712
Received: 16 June 2022; Accepted: 21 September 2022;
Published: 19 October 2022.
Edited by:
Ritchie Edward Brown, United States Department of Veterans Affairs, United StatesReviewed by:
Laura M. J. Fernandez, Université de Lausanne, SwitzerlandChristine M. Muheim, Washington State University, United States
Copyright © 2022 Han, Matsumoto, Diaz, Greene and Vogt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kaspar E. Vogt, vogt.kaspar.fu@u.tsukuba.ac.jp