Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Articles, Systems/Circuits

Global Intracellular Slow-Wave Dynamics of the Thalamocortical System

Maxim Sheroziya and Igor Timofeev
Journal of Neuroscience 25 June 2014, 34 (26) 8875-8893; https://doi.org/10.1523/JNEUROSCI.4460-13.2014
Maxim Sheroziya
1Department of Psychiatry and Neuroscience, Université Laval, Québec G1J 2G3, Canada, and
2Centre de recherche de l'Institut universitaire en santé mentale de Québec, Université Laval, Québec G1J 2G3, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Igor Timofeev
1Department of Psychiatry and Neuroscience, Université Laval, Québec G1J 2G3, Canada, and
2Centre de recherche de l'Institut universitaire en santé mentale de Québec, Université Laval, Québec G1J 2G3, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

It is widely accepted that corticothalamic neurons recruit the thalamus in slow oscillation, but global slow-wave thalamocortical dynamics have never been experimentally shown. We analyzed intracellular activities of neurons either from different cortical areas or from a variety of specific and nonspecific thalamic nuclei in relation to the phase of global EEG signal in ketamine-xylazine anesthetized mice. We found that, on average, slow-wave active states started off within frontal cortical areas as well as higher-order and intralaminar thalamus (posterior and parafascicular nuclei) simultaneously. Then, the leading edge of active states propagated in the anteroposterior/lateral direction over the cortex at ∼40 mm/s. The latest structure we recorded within the slow-wave cycle was the anterior thalamus, which followed active states of the retrosplenial cortex. Active states from different cortical areas tended to terminate simultaneously. Sensory thalamic ventral posterior medial and lateral geniculate nuclei followed cortical active states with major inhibitory and weak tonic-like “modulator” EPSPs. In these nuclei, sharp-rising, large-amplitude EPSPs (“drivers”) were not modulated by cortical slow waves, suggesting their origin in ascending pathways. The thalamic active states in other investigated nuclei were composed of depolarization: some revealing “driver”- and “modulator”-like EPSPs, others showing “modulator”-like EPSPs only. We conclude that sensory thalamic nuclei follow the propagating cortical waves, whereas neurons from higher-order thalamic nuclei display “hub dynamics” and thus may contribute to the generation of cortical slow waves.

  • cortex
  • dynamics
  • intracellular
  • slow oscillation
  • thalamus

Introduction

Slow waves are the hallmark of slow-wave sleep. They usually start from frontal and/or midline cortical areas and travel to lateral and posterior cortical areas (Massimini et al., 2004; Mohajerani et al., 2010; Nir et al., 2011; Stroh et al., 2013). Within the thalamocortical (TC) system, the slow-wave activity can be recorded from both cortex and thalamus (Contreras and Steriade, 1995). Slow waves originate in the neocortex as they can be recorded from isolated cortical structures (Sanchez-Vives and McCormick, 2000; Timofeev et al., 2000; Cossart et al., 2003), but not from the thalamus of decorticated animals (Timofeev and Steriade, 1996). During slow-wave activity, cortical neurons are active, depolarized, and fire spikes during local field potential (LFP) depth-negative waves, and they are silent and hyperpolarized during depth-positive LFP waves (Contreras and Steriade, 1995; Steriade et al., 2001; Timofeev et al., 2001; Chauvette et al., 2011). Layer V pyramidal neurons play a leading role in the onset of active states (Sanchez-Vives and McCormick, 2000; Chauvette et al., 2010; Stroh et al., 2013). Intracellular recordings from ventral lateral and ventral posterior lateral thalamic nuclei demonstrated a passive hyperpolarization during cortical silent states and active inhibitory processes, originating from the reticular thalamic nucleus (RTN), during cortical active states (Contreras and Steriade, 1995; Contreras et al., 1996; Timofeev et al., 1996). Extracellular unit recordings from the ventrobasal complex and the posterior (PO) nucleus showed that these neurons fired at different phases of the slow oscillation (Slézia et al., 2011). Several studies demonstrate that the thalamus may actively contribute to the generation of slow-wave activity (Hughes et al., 2002; Crunelli and Hughes, 2010; David et al., 2013; Lemieux et al., 2014).

First-order thalamic nuclei receive inputs mainly from cortical layer VI and project into cortical layer IV and V/VI, whereas second-order thalamic nuclei receive inputs from cortical layer V and project into layer IV, superficial cortical layers, and into subcortical centers (Jones, 2009; Constantinople and Bruno, 2013). Axons from cortical layer VI (but not from layer V) send collaterals into the RTN (Deschênes et al., 1994; Bourassa et al., 1995). The inhibitory input into first-order thalamic nuclei is provided by the RTN, whereas higher-order thalamus receives inhibitory inputs mainly from the zona incerta (ZI) (Barthó et al., 2002; Trageser and Keller, 2004).

Large excitatory synapses (drivers) formed on proximal dendrites of first-order thalamic nuclei are formed by ascending axons, whereas higher-order nuclei receive drivers originating mainly from layer V pyramidal neurons (Sherman and Guillery, 1996, 1998). A majority of modulatory inputs on TC neurons (small synapses formed on distal dendrites) are formed by corticothalamic axons.

Given the distinct patterns of connectivity, particularly in the origin of “driver” inputs and inhibitory inputs to different thalamic nuclei, we hypothesized that TC neurons from variable thalamic regions differently contribute to the cortical slow oscillation.

Here, using phase analysis for EEG and simultaneously recorded intracellular activities, we describe the global slow-wave dynamics of the TC system in ketamine-xylazine anesthetized mice that is the accepted model of slow-wave sleep (Sharma et al., 2010; Chauvette et al., 2011).

Materials and Methods

Experiments on anesthetized mice.

Experiments were performed in accordance with the guideline of the Canadian Council on Animal Care and approved by the Université Laval Committee on Ethics and Animal Research. Adult CD1 mice of both sexes were used in experiments. First, mice were anesthetized with ketamine and xylazine (100 and 10 mg/kg, i.p.) and put in a stereotaxic frame. The scalp was infiltrated with lidocaine (0.5%) before the skin incision. Two screws (stainless steel, 1 mm diameter) were inserted in the bone over the right somatosensory cortex (S1, lateral: 2.5 mm; anteroposterior: −1.5 mm from bregma; according to Paxinos and Franklin, 2001) and cerebellum (reference electrode) to monitor the EEG slow oscillation. Additional anesthesia was administrated every 30–40 min when EEG slow waves started weakening. The temperature of the platform in the stereotaxic frame where the mice were placed on was maintained at 37°C. The gentle opening of the skull was done over the area of interest, and the dura was removed. Both cortical and thalamic slow-wave activity was very sensitive to cortical damage. For example, even light damage of the somatosensory cortex disrupted slow waves in ventral posterior medial (VPM) neurons. We found that it was practically impossible to remove the dura without cortical damage in mice. To avoid any contact of sharp tools with the cortex during removing the dura, we put a small drop of 3% peroxide on the dura. Peroxide itself did not penetrate under the dura but reacted with it. In 10–15 s, a bubble appeared under the dura. Then we washed out peroxide with saline and removed saline with mineral oil (Sigma). The space created by the bubble was used to remove the dura without a touch of the cortex. If detectable damage or bleeding from the cortex (not from the bone) were observed, obtained data were not used in phase analysis (see below). To avoid cortical surface drying, a mineral oil (Sigma) was applied over the exposed area. To reduce the brain pressure and pulsation, we opened the fourth ventricle. The stability of intracellular and LFP recordings was also reinforced by covering the opened cortical surface with 4% agar in saline. Since electrical shunting, 4% agar in saline was not used in experiments specifically investigating synchrony of LFP recorded from cortical surface. At the end of the recordings, the mice were killed by a lethal dose of sodium penthobarbital (120 mg/kg, i.p.) and perfused with 4% PFA.

Recordings and labeling.

AM 3000 amplifiers (A-M Systems) with custom modifications were used for EEG and LFP recordings (sampling rate 10 kHz). LFPs were recorded using monopolar tungsten electrodes (10–12 mΩ; FHC). In the majority of cases, the tips of LFP electrodes were located 0.8–0.9 mm from the cortical surface. The EEG recordings were obtained from a screw inserted in the bone. Intracellular and LFP recordings were made in the frontal association cortex (FrA, lateral: 1–1.5 mm, anteroposterior: 2.5–2.6), orbital (lateral: 1–1.5 mm, anteroposterior: 2.5–2.6, 1.2 mm < depth < 3 mm), motor (M1, lateral: 1.5–2 mm, anteroposterior: 0.7–1.2 mm), somatosensory (S1, lateral: 1–3 mm, anteroposterior: −1 to −2 mm), visual (V1, lateral: 2.5–3 mm, anteroposterior: −2.5 to −3 mm), visual and retrosplenial (V-RS, lateral: 1–1.5 mm, anteroposterior: −3.5 to −4 mm), and the most posterior visual and retrosplenial (V-RS posterior, lateral: 3 mm, anteroposterior: −4.5 to −5 mm) cortices. Intracellular recordings from the same cortical areas and from the thalamus were performed using glass micropipettes pulled from borosilicate glass capillaries (WPI; P-97, Sutter Instrument) and filled with a solution of 2 m potassium acetate (50–70 mΩ). For recordings in the thalamus, the pipettes were also filled with 1% biocytin (Sigma). A high-impedance amplifier (Neurodata IR-283 amplifiers; Cygnus Technology) with active bridge circuitry was used to record the membrane potential and to inject current into the neurons. To stain thalamic cells, positive current pulse (200 ms on-off duty cycle, 0.5–1.5 nA, 10 min) was applied through the recording pipettes. After 1–2 h of survival, deeply anesthetized mice were perfused through the heart with saline followed by 4% PFA in 0.1 m phosphate buffer (PB). The brains were then removed and stored overnight in 4% PFA in PB at 4°C. Then 50-μm-thick frontal sections were cut.

During intracellular recordings from sensory thalamus, we also stimulated contralateral whiskers and hair on a mouse muzzle by air-puff or the contralateral eye by light flashes. To stimulate whiskers or muzzle, we used a small tube (diameter 1 mm) connected to a compressed air balloon and equipped with a trigger. To stimulate an eye, we used white light-emitting diode directed into the contralateral eye also equipped with trigger.

Retrograde tracing.

Mice were deeply anesthetized with isoflurane, placed in a stereotaxic frame as above, and a small opening was done over targeted thalamic nuclei. Glass micropipettes with a tip diameter 30–40 μm were filled with 0.2% cholera toxin subunit B (CtB, List Biological Laboratories); and using microinjector (Drummond Scientific), 9.2 nl of CtB was delivered into the thalamus. CtB was injected either in VPM (lateral: 1.6 mm, anteroposterior: −1.6, depth: 3.2–3.3 mm) or PO (lateral: 1.2–1.3 mm, anteroposterior: −2.0 to −2.2 mm, depth: 3 mm) thalamic nucleus. After the injection, the micropipettes were left inside for 15 min to reduce the leakage along the injection tract. The openings were then covered with dental cement, and the mice were placed into their home cages. After 7 d of survival, the mice were deeply anesthetized with ketamine and xylazine and perfused through the heart with saline followed by ice-cold 4% PFA in PB. The brains were then removed and stored overnight in 4% PFA in PB at 4°C. Then, 50-μm-thick frontal sections were cut, and every other was collected.

Immunohistochemistry.

To visualize biocytin-labeled cells, free-floating sections were incubated in ABC kit (1:200; Vector Laboratories) in TBS for 2 h. After rinses, sections were incubated in TBS containing 0.05% DAB (Sigma), 0.005% CoCl2, and 0.00125% H2O2 for 5–10 min.

CtB-containing sections were first incubated in 2% normal rabbit serum, 0.5% BSA, and 0.3 Triton X-100 in PBS overnight at 4°C. The incubation in goat anti-CtB (1:4000; List Biological Laboratories), 2% normal rabbit serum, 0.5% BSA, and 0.3 Triton in PBS for 4 d at 4°C followed. After rinses, sections were incubated in biotinylated rabbit anti-goat IgG (1:200; List Biological Laboratories) for 2 h, and then ABC and DAB reaction followed like above. In the end, all sections were mounted onto gelatin-coated glass slides, dehydrated, and coverslipped. Microphotographs were taken, and Image Pro software (Media Cybernetics) was used for automatic quantification of the labeled cells (see below).

Electrophysiological data analysis.

To detect phases of the slow oscillation, we used the Hilbert transformation (Fig. 1A) for which the EEG was filtered (bandpass frequency 0.2–4 Hz). The extracted phase changes in a range of −1 to 1 of fraction of Pi. The instantaneous phase ±1 corresponds to the middle of the negative peak of the cortical EEG, whereas the instantaneous phase 0 corresponds to the middle of the positive peak of the cortical EEG (Fig. 1A). For phase analysis, we extracted spikes and postsynaptic potentials/noise from simultaneously recorded intracellular activities. Only stable intracellular recordings (at least 1 min recordings, but usually 2–5 min) were used for phase analysis (stable membrane potential, <−60 mV during silent states, spikes amplitude >45 mV). Spikes were detected using a threshold. IPSPs (negative peaks) were detected manually. To extract EPSPs, the intracellular activity was filtered between 0 and 1000 Hz, then differentiated, and the threshold was set on. To extract all EPSPs and/or synaptic noise, we used a threshold of 0.03–0.05 V/s, whereas to extract large-amplitude EPSPs in the thalamus, we used a threshold of 0.1 V/s (Fig. 1A). IPSPs were detected manually based on sharp negative deflection of >5 mV and their characteristic shape (Fig. 1A). As expected, in first-order thalamic nuclei, the time histogram of EPSPs with reference to the phase −0.5 was not modulated by the slow oscillation, but the IPSPs were clearly coupled to cortical active states (Fig. 1B). After instantaneous phases were determined, we plotted phase event histograms (Fig. 1C, top). The EEG signal is not completely sinusoidal, which leads to incorrect phase quantification when plotting phase histograms (Siapas et al., 2005; Mallet et al., 2008). Nonsinusoidal EEG waveform comes out from the time asymmetry of the waveform when the negative and/or positive phases of the EEG wave are longer in time than the transitions between them. To correct the phase analysis, we calculated how much time each phase bin occupied (Fig. 1C, middle, time-phase histogram). Then we divided the phase events histogram by the time-phase histogram (all the histograms were composed of 50 bins; i.e., normalized phase by time; Fig. 1C, bottom). All phase histograms were normalized by the total number of detected events divided by the number of bins (n = 50 in this study). A modification of the Rayleigh test for vector data (Moore–Rayleigh test) was used to select phase-locked neurons (significance level, p < 0.05) (Moore, 1980). For all phase-locked neurons, we calculated mean vectors. The Watson–Williams test was used to compare mean phases (angles) of groups of neurons (significance level, p < 0.05).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

The methodological approaches for the phase and time analysis of field potential and intracellular signals. A, Drawing represents a position of two electrodes used for EEG recordings. In all experiments, the EEG was recorded from the contralateral side compared with intracellular recordings. The EEG signal was filtered (bandpass 0.2–4 Hz), and then a Hilbert transform was performed. To extract EPSPs, we differentiated the intracellular potential over time and then we set up the threshold. Vertical red lines indicate identified EPSPs. IPSPs (three large-amplitude negative peaks are evidently seen here) were detected manually. B, Top, Averaged EEG. The zero time corresponds to the instantaneous phase −0.5. Bottom, Time histograms of EPSPs and IPSPs normalized by mean. C, Top, EPSP phase histogram for the same neuron before normalization by time. Because time spent in each phase is not identical, an artificial correlation between phase and EPSP occurrence is seen. Bottom, A time against phase histogram, which shows the average time duration of each phase bin (%). Because the amount of counts within each phase bin is proportional to its time duration, we divided the EPSP phase histogram by the time against phase histogram (i.e., we normalized the phase by time, bottom histogram). As expected, the EPSPs from ascending pathways were not phase-locked to the EEG (p > 0.05, Moore–Rayleigh test). D, Averaged EEG wave and intracellular potential plotted against phase using linear interpolation. EEG normalized by maximum values. Note averaged IPSPs on the intracellular wave.

We also plotted mean EEG, LFP, absolute value of γ LFP signal (bandpass frequency 50–300 Hz), and intracellular activities against phase by using linear interpolation (see Figs. 1D and 3). Considering each point of the wave plotted against phase as a vector (the amplitude of the wave gives the length of the vectors), we calculated the mean phase of the intracellular wave itself. We also calculated the first derivatives of the intracellular waves and EPSP histograms by phase and their mean vectors (see below). Before differentiation, waves and histograms plotted against phase were smoothed with a binomial filter.

Perievent time histograms were also plotted for extracted postsynaptic potentials. To quantify rising and descending slopes of the histograms, we measured the times (phases) and amplitudes of the histograms at their half amplitude and additionally at one-tenth amplitude for decay measured from minimum (see Figs. 9 and 11).

To reveal the delays between simultaneous LFP and intracellular recordings, the crossings at half of amplitude within each slow-wave cycle were detected (Fig. 2E). To analyze delays between LFP γ activities, we calculated their absolute values (γ2)1/2, and then the crossings with the threshold (7 SD of silent periods) corresponding to the transitions between active/silent states were detected (Fig. 2D).

Wave analysis and circular statistics were performed with Igor Pro version 6.2 software (WaveMetrics). Basic statistical and other data analysis was performed using Sigma Plot (Systat Software). Data are expressed as mean ± SD. or SEM when indicated. If the normality and equal variance tests passed (p > 0.05), then a two-tailed (parametric) t test was used; otherwise, a Mann–Whitney sum rank test (nonparametric) was used (significance level, p < 0.05).

Results

Slow-wave propagation over the cortex

We recorded slow-wave activity with monopolar LFP electrodes either simultaneously from the left M1, S1, V1, and V-RS cortices (Fig. 2A) or in pairs. In all experiments, EEG recorded from the contralateral somatosensory cortex was used for Hilbert transformation. Time and phase shifts between slow waves from M1, S1, V1, and V-RS cortices were apparent when LFP activities were recorded in deep layers (depth 0.7–0.9 mm, n = 6 animals; Fig. 2A,C). When LFP activities were recorded from the wet cortical surface, no apparent shifts between slow waves from M1, S1, and V1 cortical areas were found (Fig. 2B; n = 3 animals). The same result was obtained with EEG screws inserted in the skull on these cortical areas (n = 3 animals; data not shown). However, if LFP activities were recorded from the dry cortical surface, positive large-amplitude (4.5 ± 1 mV) propagating signals were detected (n = 3 animals; Fig. 2E). The active states clearly tended to terminate simultaneously within different cortical areas. Application of agar (based on saline) on the cortical surface diminished the amplitude of LFP signal (1.1 ± 0.2 mV) as well as eliminated the delays in the onset of active states (Fig. 2E). Because γ activity is a proven indicator of active states (Mukovski et al., 2007), we calculated the absolute values of LFP γ (50–300 Hz) activities recorded either from the surface or deep layers (Fig. 2D). Then we detected delays between the onsets of active or silent states (up delays or down delays correspondingly) as shown in Figure 2D, E between slow waves from different cortical areas and plotted delay as histograms (Fig. 3). We also plotted averaged LFP and absolute value of γ LFP activities against phase (Fig. 3). EEG and LFP of wet cortical surface (Fig. 3A,B) did not show propagation of slow-wave activity over the investigated cortical hemisphere. Normalized and plotted against phase, the LFP waves from wet cortical surface also coincided with the EEG wave from the contralateral hemisphere (Fig. 3B1). LFP slow wave from dry surface and γ LFP activities from dry surface or from deep layers revealed propagation of the leading edge of the surface positive (active state) wave in the anteroposterior direction and simultaneous termination on average of the wave (Fig. 3C–E). The γ LFP signal recorded from depth had higher amplitude compared with surface γ signal (Fig. 3D1,E1). Generally, slow waves started from the M1 area, propagated to S1, and then to V1 and V-RS areas (n = 6 animals; Fig. 3E1). The distributions of up delays (Fig. 3E2) revealed fluctuations of anteroposterior propagation (probability 30 ± 18% from S1 into M1, 8 ± 7% from V1 to M1, 3 ± 3% from V-RS to M1, n = 6 animals; Fig. 3E2). Although the probabilities of posterior–anterior propagation obtained from large-amplitude LFP signal of the dry surface was lower (probability 10 ± 8% from S1 into M1, 2%–3% from M1 to V1 and V-RS, n = 3 animals; Fig. 3C2). The mean up M1-S1, M1-V1, and M1 (V-RS) delays obtained from γ LFP signal recorded in deep layers were 48 ± 45, 98 ± 46, and 147 ± 42 ms, respectively (up M1-S1 vs up M1-V1, p = 0.054; up M1-S1 vs up M1-(V-RS), p = 0.003). The mean down M1-S1, M1-V1, and M1-(V-RS) delays were 14 ± 12, 24 ± 19, and 28 ± 14 ms, respectively (down M1-S1 vs down M1-(V-RS), p = 0.1). Taking into account a large scatter of down delays (almost ±200 ms centered around 0 ms; Fig. 3E3) compared with the relatively small (∼20 ms) and insignificant differences in mean down delays, we conclude that, in the slow-wave time scale, the active states tended to terminate simultaneously within different cortical areas.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

LFP/EEG recordings in the mouse cortex. A, Simultaneously recorded LFP and EEG activities in the motor, somatosensory, and visual cortex in deep layers (0.7–0.9 mm) as estimated from manipulator reading. The location of electrodes is schematically shown by colors on the mouse skull drawing. Left, Superposition of the beginning and the end of active states at a higher temporal resolution. B, LFP activities recorded from the wet cortical surface at the same points as in A. C, Averaged LFP waves from deep layers plotted against phase. In each animal (n = 6 mice), 2 min segments of LFPs were recorded simultaneously by pairs from M1 and successively in S1, and then V1, and RS (dotted line, the cortical area is shown by dotted point on the mouse skull drawing). Before averaging in a group, waves were normalized by maximum values. The positive part of SEM is shown with gray (negative part is not shown). D, An example of paired LFPs from deep layers in M1 and V1 and lower the corresponding absolute values of γ (50–300 Hz) activities are shown. The delays between γ LFP signals were detected using the crossings with the thresholds (7 SD of silent state activities). E, An example of paired LFPs from dry cortical surface in M1 and S1 before and after (bottom traces) of agar application on the cortex. The delays between LFP signals were detected using the crossings with half of wave amplitude. On average, the delays were insensitive to the threshold amplitude within significant range (shown with gray). F, The circuit that was used to measure a passive impedance of the cortex. The pipettes were filled with 3 m NaCl in 4% agar to avoid leakage, and chlorinated silver wires were used inside the pipette. Sinusoidal potential was applied for measurements (2 Hz, amplitude 0.4 V).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Propagation and termination of cortical LFP/EEG activities. Activities were recorded simultaneously by pairs in M1 and successively in S1, and then V1, and RS (color code as in Fig. 2). A, EEG activities recorded with screws were coherent across M1, S1, and V1. A1, Averaged EEG signals from M1, S1, and V1 are plotted against phase (n = 3 animals). Before grand averaging, waves were normalized by maximum values. The histograms for up (A2) and down (A3) M1–S1 and M1–V1 delays of the same EEG recordings are shown. The delays were detected as shown in Figure 2E. When S1 preceded M1 (V1 preceded M1), the time intervals calculated as negative. B, The similar plots for LFPs of wet cortical surface (n = 3 animals). B1, Black line indicates contralateral EEG. Other conventions as in A2 and A3. C, The similar plots for LFPs of dry cortical surface are shown (n = 3 animals). Other conventions as in A1 and A3. D, Plotted against phase, the absolute values of γ LFP signals (D1) recorded from the dry cortical surface, and up (D2) and down (D3) delay histograms are shown (n = 3 animals). Other conventions as in A1 and A3. The delays were detected as in Figure 2D. E, The similar plots for the absolute value of γ LFP signals recorded from the deep layers (n = 6 animals). In all panels, the positive part of SEM is shown with gray (negative part is not shown).

We did not detect neuronal correlations, which could explain phase and time shifts between wet surface and deep LFP slow-wave signals in mice (Fig. 4H). The simplest explanation for low variability of surface EEG (Fig. 3A) and LFP recorded from the cortical surface covered with liquid or agar (Fig. 3B) is a shunting effect, which leads to a spreading of the extracellular potential along the cortical surface compared with deep layers. To test this possibility, we measured the impedance between two pipettes, one of which was located on the wet cortical surface and the other was either on the surface or in deep layers (Fig. 2F). The impedance of the pipettes themselves was 190 ± 50 kΩ. When the second pipette was on the cortical surface and another one was inserted into the depth, the impedance of cortical tissue was 43 ± 25 kΩ (n = 5 pipette pairs). We failed to get a satisfactory estimation for the impedance of tissue on the surface and any distance dependence apparently because of pulsation of the cortex (measured impedance was 1 ± 7 kΩ). For comparison, the impedance between the same pipettes connected with saline-moistened napkin 5 × 10 mm was 6.5 ± 1.3 kΩ. These results indicate that the surface LFP signal depends on experimental conditions, whereas EEG signals recorded from screws do not represent the traveling wave but are reflecting the large-scale brain activity and therefore provide a reliable reference for the slow-wave phase detection.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Paired intracellular recordings from different cortical areas. A, The examples of simultaneously recorded intracellular activities of two neurons from the motor and somatosensory, and from the somatosensory and visual cortices. The cortical locations are indicated by colors. Note the time delays between active state onsets in the cell pairs. B, Superposition of averaged membrane potentials of n = 5 M1–S1 and 3 S1–V1 neuronal pairs are plotted against phase (left) and time (right). In the time plots, zero time corresponds to the instantaneous phase −0.5. C, The methodological approach for the delay analysis of paired intracellular recordings. The crossings at half of amplitude within each slow-wave cycle were detected. The delays were almost insensitive to the threshold amplitude within significant range (shown with gray). The time intervals between crossings were calculated and represent up and down delays as shown. When S1 preceded M1 (V1 preceded S1), the time intervals calculated as negative. Duration of active state calculated as the time interval between down and up crossings. D, Superposition of all M1–S1 up and down delays of 5 neuronal pairs (n delay plotted against n + 1 delay). E, Distribution of all M1–S1 (shown with black) and S1–V1 (shown with pink) up (left) and down (right) delays divided by the corresponding durations of M1 and S1 active states, respectively. F, All M1–S1 and S1–V1 up delays plotted against the duration of M1 and S1 active states, respectively. Black lines indicate linear fit for each neuronal pair. G, The mean up delays between neuronal pairs plotted against depth difference of neurons. The depth differences were calculated as the difference between M1 and S1 neurons (between S1 and V1, shown with gray). H, Left, Superposition of EPSP phase histograms normalized by mean representing the intracellular activity of 11 neurons recorded in the frontal cortex at different depths. Middle, Right, Mean vector phases and lengths calculated from the same EPSP histograms plotted against depth from which the neurons were recorded.

To average intracellular activities from different cortical and subcortical structures obtained in independent experiments, the most integrative and stereotypic process should be taken as a reference wave. This process should not be influenced by local fluctuations, and it should reflect as much as possible an averaged large-scale brain activity. In previous studies, to eliminate the influence of time variability of slow-wave duration, filtered EEG was taken for phase analysis with Hilbert transform (Mallet et al., 2008; Slézia et al., 2011). Because the delays of propagation of EEG were minimal (Fig. 3A), in the remaining part of the study, we used this signal as reference to align intracellular recordings obtained in different cortical or thalamic locations.

Paired intracellular recordings from different cortical areas

To investigate the cellular basis of the results obtained from LFP activities, we performed dual intracellular recordings of neuronal activities from M1 and S1 (n = 5 pairs of neurons) and from S1 and V1 (n = 3 pairs of neurons) cortices. In all recorded pairs of neurons, a majority of M1 active states started off before S1 active states, whereas active states in S1 also started before V1. Examples are shown in Figure 4A. Using Hilbert transform for simultaneously recorded EEG activities, we extracted the phase for each recording and plotted membrane potentials from pairs of neurons against the phase (Fig. 4B, left) and time (Fig. 4B, right). In time plots, zero time corresponded to the instantaneous phase −0.5. Averaged membrane potentials on phase and time plots revealed phase and time delays between neurons from different cortical areas. The mean time delays (measured at half rising times; Fig. 4B, right) in M1-S1 and S1-V1 neuronal pairs were 46 ± 11 ms and 55 ± 25 ms, respectively. Considering the mean distance from M1 to S1 and from S1 to V1 as ∼2 mm, the slow-wave propagation velocity was 40 ± 13 mm/s, which was similar to the velocity previously described in other studies (Mohajerani et al., 2010; Ruiz-Mejias et al., 2011).

To show the variability of the slow-wave propagation, we also calculated up and down delays between individual active states of simultaneously recorded neurons (Fig. 4C). Up and down delays fluctuated around the mean values and did not reveal any temporal/serial dependencies (Fig. 4D). The probabilities of posterior–anterior propagation were 14.5% (from 6% to 25%) and 9% (from 3% to 13%) in M1-S1 and S1-V1 pairs of neurons correspondingly. As in LFP experiments, the intracellular active states tended to terminate simultaneously (except one V1 neuron; Fig. 4B). Similar to a previous study on cats (Volgushev et al., 2006), individual down delays revealed large fluctuations centered around 0 ms (Fig. 4D,E). The mean values of SEM for up and down delays in M1-S1 pairs were 5.74 and 6.02 ms, respectively (the extent of fluctuations was similar). We did not detect any apparent correlation between the depth of neurons, the time delays between them, and the preferential phase of involvement of individual neurons in the slow oscillation (Fig. 4G,H). In all eight pairs of neurons recorded in this study, we found a positive correlation between the up delays and the durations of active states (Fig. 4F).

Intracellular slow-wave patterns in the mouse cortex and thalamus

To investigate how cortical and thalamocortical cells of different locations are involved in slow oscillation, we performed intracellular recordings within six cortical areas along the anterocaudal direction (see Materials and Methods) and from specific and nonspecific thalamic nuclei. Together, we recorded 18, 8, 11, 5, 7, and 7 neurons from FrA-O, M1, S1, V1, V-RS, and V-RS posterior cortical areas and 159 cells in the thalamus. Thalamocortical cells were identified either by their ability to respond with a characteristic low-threshold spike (LTS) in response to a hyperpolarizing current pulse or by morphological location. Within VPM nucleus and lateral geniculate nucleus (LGN), neurons were also usually identified by their response to sensory stimulation (whiskers and hair: air puff and eye–light flash, respectively).

As expected, all cortical neurons demonstrated an alternation of depolarized active (UP) and hyperpolarized silent (DOWN) states, whereas thalamic neurons revealed two distinct slow-wave patterns. Eighty-four thalamic neurons displayed unitary and/or spindle-like IPSPs, which followed cortical slow waves. The other 75 thalamic neurons displayed depolarized active states, which accompanied cortical slow waves, although occasionally we could also identify IPSPs. The examples of EEG and the two characteristic intracellular slow-wave patterns of TC neurons recorded simultaneously with cortical S1 neurons are shown in Figure 5.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Two intracellular slow-wave patterns recorded from the mouse thalamus. A, Simultaneously recorded intracellular activities of thalamic and cortical S1 neurons and contralateral EEG during the slow oscillation in anesthetized mice. Left, The thalamic neuron reveals a characteristic LTS in response to a hyperpolarizing (100 ms, 0.5 nA) current pulse. During cortical depolarizing active states, the TC neuron displayed unitary or repetitive IPSPs. Thalamic IPSPs follow the onset of cortical active states. B, A segment of EEG and intracellular activities recorded from another pair of thalamic (AD) and cortical neurons. A morphologically identified thalamic neuron (left) did not generate LTS in response to a hyperpolarizing current pulse (inset) of different amplitude (data not shown). The thalamic neuron generated active states composed of multiple depolarizing events that followed cortical active states. There is a time delay between the onset of intracellular cortical and thalamic active states (right).

Within particular thalamic nuclei, we identified the location of 51 neurons (Table 1; Figs. 5 and 6). All the neurons from the higher-order nuclei we identified histologically. Neurons from sensory nuclei were identified either by their response to sensory stimulation or both histologically and by their response to sensory stimuli. The neurons within VPM, LGN, PO, and AV nuclei generated LTS in response to a hyperpolarizing current pulse, whereas the neurons from PF, central lateral (CL), and AD did not generate LTS. Usually, we did not observe spontaneous LTS in the mouse's thalamus (Figs. 5 and 6), and only some PO neurons displayed spontaneous LTS when steady hyperpolarizing current was intracellularly injected (data not shown). Intralaminar PF neurons displayed the lowest input resistance (Table 1), which was significantly different from the input resistance of PO neurons (16 ± 5 mΩ vs 25 ± 6 mΩ, respectively, p = 0.005), of AV neurons (24 ± 8 mΩ, p < 0.05), and VPM and LGN neurons together (23 ± 7 mΩ, p < 0.05).

View this table:
  • View inline
  • View popup
Table 1.

Histologically identified neurons within the thalamusa

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Two thalamic slow-wave patterns recorded from the first- and higher-order nuclei. A, Left, Light microscopic images of the intracellularly stained VPM and LGN neurons. There is a specific first-order “bushy” morphology of the neurons. The neurons revealed LTS in response to hyperpolarizing (100 ms, 1 nA) current pulses. Middle, Simultaneously recorded cortical contralateral EEG and intracellular activities from VPM and LGN nuclei. Both neurons displayed unitary IPSPs during cortical active states and “fast-rising” EPSPs, which were not modulated by cortical states. Right, Selected IPSPs. B, Contralateral sensory stimulation (whiskers and hair on the muzzle-air puff, eye-flash) elicited synaptic responses consisting of IPSPs and fast-rising EPSPs in VPM and LGN neurons. Top and bottom traces represent the responses of the neurons to sensory stimulation at depolarized and hyperpolarized membrane potential respectively with injection of steady inward and outward current of 0.5 nA. C, Left, Light microscopic images of the intracellularly stained PO and PF neurons. In contrast to PO neurons, PF neurons did not reveal LTS in response to hyperpolarizing (100 ms) current pulses. Right, Simultaneously recorded cortical EEG and intracellular activities from PO and PF nuclei. Both neurons displayed the slow-wave pattern with depolarization during active states.

Using interpolation, we plotted membrane potentials against phase for all stable intracellular cortical and thalamic recordings (Fig. 7). At least 1 min segments of recordings with the robust slow oscillation of well-anesthetized mice were used. On average, the neuronal membrane potentials from PO (n = 7 neurons), parafascicular (PF, n = 6 neurons), CL (n = 1 neuron, data not shown), anterior dorsal (AD, n = 1 neuron; data not shown), and anterior ventral (AV, n = 7 neurons) nuclei displayed depolarized waves (averaged active states), which accompanied cortical depolarized waves (Fig. 7). In 3 of 7 analyzed PO cells, we also recorded occasional IPSPs that produced a local minimum during active states (Fig. 7, red dotted line). The averaged intracellular slow waves within VPM (n = 14 neurons) and LGN (n = 6 neurons) nuclei were dominated by inhibitory activities that produced averaged IPSPs with negative peaks <−65 mV. The membrane potential depolarization of all cortical cells as well as from all the neurons within thalamic PO, PF, CL, AD, and AV nuclei were phase-locked (p < 0.05, Moore–Rayleigh test).

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Membrane potentials of thalamic and cortical cells are modulated by slow oscillation. Averaged membrane potentials of neurons recorded in the cortex (FrA-O, M1, S1, V1, V-RS, and V-RS posterior: n = 18, 8, 11, 5, 7, and 7 neurons, respectively) and thalamus (PF, PO, AV, VPM, and LGN: n = 6, 7, 7, 14, and 6 neurons, respectively) plotted against phase. In all panels, the positive part of SEM is shown with gray (negative part is not shown). The cortical areas are schematically shown in the left.

Synaptic activities during slow oscillations in the mouse thalamus and cortex

The main excitatory inputs formed on thalamocortical neurons are usually divided to drivers and modulators based on the EPSP amplitude and rising time (Li et al., 2003; Reichova and Sherman, 2004). To our knowledge, no driver-like EPSPs were reported in the cortex. Typical examples of slow-wave cycles from different thalamic nuclei and from the somatosensory cortex are shown in Figure 8A. Large-amplitude, well-defined unitary EPSPs were a common feature of spontaneous activity of all the neurons recorded within PO, AV, VPM, and LGN nuclei (Fig. 8A, red arrows). Sensory responses within projecting VPM and LGN neurons were also composed of large-amplitude, well-defined unitary EPSPs as well IPSPs (Fig. 6B). Small-amplitude EPSPs occurring in groups were present in neurons from all recorded nuclei (Fig. 8A, blue arrows), but within projecting VPM and LGN nuclei they were weak, even if the membrane potential was hyperpolarized with injected current. We consider EPSPs called here “large-” and “small-” amplitude corresponding to thalamic drivers and modulators, respectively. We did not detect large-amplitude EPSPs in cortical neurons (Fig. 8A), so all cortical EPSPs displayed modulator-like behavior.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Drivers and modulators during slow oscillations in the mouse thalamus. A, Typical examples of neuronal activities recorded during slow-wave cycles from different thalamic nuclei and from the somatosensory cortex. The traces were aligned according to contralateral EEG (top). Red arrows and blue arrows indicate examples of unitary large-amplitude and merged together small EPSPs, respectively. Spontaneous large-amplitude EPSPs were a characteristic feature of all intracellular activities recorded within VPM, LGN, PO, and AV nuclei. AV neurons generated bursts of EPSPs. B, Averaged large-amplitude EPSPs recorded from four thalamic nuclei (VPM, LGN, PO, and AV: n = 9, 7, 11, and 10 neurons, respectively). The plots represent the mean amplitudes (right) and kinetics (bottom) of the EPSPs. *p < 0.05 (t test or Mann–Whitney sum rank test). **p < 0.01 (t test or Mann–Whitney sum rank test). ***p < 0.001 (t test or Mann–Whitney sum rank test). Large-amplitude EPSPs within VPM and LGN neurons, which were not modulated by cortical states and responded to sensory stimulation displayed the shortest rising time (fast-rising EPSPs).

First, we compared mean amplitudes and kinetics (half rise and half decay time counting from the EPSP peak; Fig. 8B) of the large-amplitude EPSPs within VPM, LGN, AV, and PO (n = 9, 7, 10, and 11 neurons, respectively). All measurements were performed without injection of intracellular steady current and with a membrane potential of ∼−60 to −75 mV. We found that EPSP amplitude within PO neurons was significantly higher than within VPM, LGN, and AV (7.37 ± 0.45 mV vs 4.15 ± 0.25, 4.32 ± 0.25, and 4.72 ± 0.47 mV, respectively, mean ± SEM, p < 0.001; Fig. 8B). Both rising and decaying phases of EPSPs within VPM, LGN, and AV nuclei were significantly (p < 0.05) faster than they were within second-order PO neurons (Fig. 8B). “All or none” behavior of driver EPSPs was reported previously (Li et al., 2003; Reichova and Sherman, 2004; Groh et al., 2008); therefore, we also calculated the ratio (SD)/(mean amplitude) for the large-amplitude EPSPs extracted in each neuron (an “all or none” behavior supposing an SD/mean ratio ≪ 1). We found that the ratio SD/mean amplitude was significantly smaller in VPM neurons compared with LGN, AV, and PO neurons (0.153 ± 0.038 vs 0.322 ± 0.126, 0.233 ± 0.032, and 0.397 ± 0.13, respectively, p ≤ 0.001). The ratios for AV and PO neurons were significantly different, too (p = 0.007). There were no evidently pronounced large-amplitude EPSPs within the other investigated nuclei (PF, CL, and AD), and we failed to detect unitary EPSPs in neurons from these nuclei either automatically or manually (Fig. 8A).

To reveal slow-wave modulation of the thalamic and cortical synaptic activities, we plotted phase and time PSP (and/or synaptic noise) histograms (Fig. 9). We also analyzed spiking activities. All detected cortical and thalamic PSPs were modulated by slow waves, with the exception of large-amplitude EPSPs (drivers) recorded within VPM and LGN nuclei (Fig. 9A). PSPs and spikes of all cortical cells as well as from all the neurons within thalamic PO (n = 7 neurons), PF (n = 6), AV (n = 7), CL (n = 1, data not shown), and AD (n = 1, data not shown) nuclei were phase-locked (p < 0.05, Moore–Rayleigh test). Extracted IPSPs from VPM (n = 14 neurons) and LGN (n = 6) neurons were also phase-locked (p < 0.05, Moore–Rayleigh test). In 4 of 14 analyzed VPM neurons, we extracted phase-locked small-amplitude EPSPs/noise corresponding to modulators (Figs. 8A and Fig. 9A, bottom, dotted black line, p < 0.05, Moore–Rayleigh test). However, when all large-amplitude EPSPs (corresponding to thalamic drivers) or small- and large-amplitude EPSPs were extracted together, they were not phase-locked in the VPM (Fig. 9A, bottom, solid black line). We failed to extract small-amplitude EPSP/noise modulated by cortical slow waves in LGN neurons, which were shown in Figure 8A. All EPSPs from that nucleus were not phase-locked to the slow oscillation (Fig. 9A, bottom, solid red line). Spike activities of only three VPM and no LGN neurons were phase-locked to the slow oscillation (p < 0.05, Moore–Rayleigh test; data not shown).

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Synaptic activities during slow oscillations in the mouse thalamus and cortex. A, Averaged PSP histograms of the neurons recorded in the cortex (FrA-O, M1, S1, V1, V-RS, and V-RS posterior: n = 18, 8, 11, 5, 7, and 7 neurons, respectively) and thalamus (PF, PO, AV, VPM, and LGN: n = 6, 7, 7, 14, and 6 neurons, respectively) plotted against phase. Here and below before grand averaging, the histogram for each neuron was normalized by the mean value (each bin of the histogram was divided by the value of mean bin of the histogram). Bottom, The histograms for small-amplitude EPSPs extracted in 4 VPM neurons are shown with a dotted line. The histograms for fast-rising EPSPs extracted in VPM and LGN neurons are shown with a solid line. B, Half rising amplitude of PSP histograms of each neuron is plotted against phase (averaged mean phases of the same points are shown in Fig. 11B). Bottom, Half rising amplitude of IPSP histograms of the VPM and LGN neurons is plotted against phase. C, Time PSP histograms. Zero time corresponds to the instantaneous phase −0.5. In all panels, the positive part of SEM is shown with gray (negative part is not shown).

Global slow oscillation dynamics in the mouse thalamocortical system

The averaged membrane potentials (Fig. 7) and PSP histograms (Fig. 9) revealed intacellular slow-wave dynamics in the thalamocortical system. The earliest structures involved in the generation of slow waves were the FrA-O and M1 cortices as well as PF and PO nuclei. Then, the intracellular active states propagated in the anteroposterior direction that was also shown with LFP recordings (Fig. 3E). The latest structures involved into the slow oscillation were V-RS posterior cortex and AV nucleus with an overall delay of approximately half of a cycle from the earliest structures. The mean phases (and significant differences between them) of membrane potentials, EPSPs, and spike histograms are shown in Tables 2 and 3. Both PF and PO half rising times were significantly different compared with rising times from S1, V1, V-RS, and V-RS posterior cortical neurons (Table 3). The mean PF phases obtained from EPSP histogram and from membrane potentials significantly differed from all the corresponding cortical mean phases (p < 0.05, Watson–Williams test, Table 2). The mean phases depend on the waveform and its duration. A more precise analysis representing the order of activation was obtained from the velocity of slow-wave activation (depolarization) plots (Fig. 10A,B; same data as in Figs. 7 and 9). The phase shifts were clearly revealed between thalamic and cortical neuron activation (Fig. 10A,B). The cortical mean phases of membrane potential derivatives (>0) and EPSP histogram derivatives (>0) were significantly different compared with each other respectively, except FrA-O versus M1 (p < 0.05, Watson–Williams test; mean phases of individual neurons are shown in Fig. 10A2,B2). The thalamic PF mean phases of membrane potential derivatives (>0) were significantly different from all the corresponding cortical mean phases (p < 0.05, Watson–Williams test). The thalamic PO mean phases of membrane potential derivatives (>0) were not significantly different from FrA-O and M1 phases. The thalamic PF and PO mean phases of EPSP histogram derivatives (>0) were significantly different from the corresponding mean S1, V1, V-RS, and V-RS posterior phases (p < 0.05, Watson–Williams test) but were not significantly different from FrA-O and M1 phases. Thus, the phases obtained from the membrane potential derivatives revealed that PF and PO led the slow-wave cycle, whereas the phases from EPSP histograms showed a synchrony between thalamic PF/PO and the frontal cortical areas. In Figure 10E, the membrane potential (left) and EPSP histogram (right) derivatives (>0) of PF and orbital cortex are plotted together, and their peaks are matched. The apparent statistical difference between PF and the frontal cortical areas obtained from the membrane potential derivatives is rather the result of the waveform and its duration since the peaks of activation occurred simultaneously (Fig. 10E). Taking into account the precision of the method (e.g., filtering, smoothing, Hilbert transformation), we conclude here that, on average, the thalamic PF and PO neurons started activation/depolarization simultaneously with the frontal cortical areas and before the other investigated parts of the cortex. Thalamic AV neurons were activated significantly after all other cortical and thalamic neurons, except V-RS posterior phase (Fig. 10A,B; Tables 2 and 3).

Table 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Table 2.

The number of neurons are indicated. Above the diagonal: significant differences between the mean phases from EPSP histograms (Fig. 9A) are shown by gray. Below the diagonal: significant differences between the mean phases from the averaged membrane potentials (Fig. 7) are shown by gray. p < 0.05 (Watson–Williams test).

Table 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Table 3.

The numbers of neurons are indicated. Above the diagonal: significant differences between the mean ½ rising times obtained from time histograms (represented as mean ± SEM, Figs. 9B and 11A). Below the diagonal: significant differences between the mean phases obtained from spike histograms (not shown). Differences are shown with different intensities of gray (for the phases, Watson–Williams test; for the times, two-tailed t test or Mann–Whitney test).

Figure 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 10.

Global slow-wave dynamics within the thalamocortical system of a mouse. A1, Differentiated by phase membrane potentials are plotted against phase of EEG slow oscillation. The averaged positive component of the derivative is shown, which represents the velocity of depolarization in cortical and thalamic neurons. Here and below before grand averaging, the plot for each neuron was normalized by the mean value. A2, The mean vectors (length) of the same components as in A1 are shown for individual neurons, which represent the mean phase of depolarization of the neurons during slow wave. B1, Differentiated by phase, phase EPSP histograms are plotted. The averaged positive component of the derivative represents the increments of EPSP number during transitions to the active states. B2, The mean vectors (length) of the same components as in B1 are shown for individual neurons, which represent the mean phase of the increments of EPSP number during transitions to the active states. C1, The averaged negative component of the derivative represents the velocity of hyperpolarization in cortical and thalamic neurons. C2, The mean vectors (length) of the same components as in C1 are shown for individual neurons, which represent the mean phase of hyperpolarization during the transition from active to silent state. D1, The velocity of termination of active states estimated from EPSPs histograms for cortical neurons (top) and thalamic neurons (bottom). D2, The mean vectors (length) of the same components as in D1 are shown for individual neurons, which represent the mean phase of EPSPs active state termination. E, The positive components of derivatives from the thalamic PF and orbital cortical neurons are plotted together, and their amplitudes are normalized by maximum of the first peak (left, obtained from membrane potentials; right, from EPSP histograms). In all panels, the positive part of SEM is shown with gray (negative part is not shown). F, The sequence of detection of slow waves in the cortex and thalamus is schematically shown.

We also found a significant phase shift between IPSPs within VPM and LGN nuclei (VPM vs LGN mean phases: −0.326 vs 0.145, p < 0.05, Watson–Williams test; Fig. 9A). The IPSPs from the VPM nucleus occurred significantly later than M1 (but not S1) cortex activation (mean phases, p < 0.05, Watson–Williams test) and IPSPs from the LGN nucleus occurred after the activation of both M1 and S1 (but not V1) cortices (mean phases, p < 0.05, Watson–Williams test). The inhibition within the sensory thalamus occurred statistically simultaneously with the activation of the corresponding primary cortical areas.

Despite an overt anteroposterior propagation of active state onsets, the termination of cortical active states occurred almost simultaneously in all investigated areas (Figs. 3E, 7, and 9). To support further this observation, we compared the velocities of active state termination in different cortical areas (negative components of the first derivatives; Fig. 10C,D). The averaged curves of cortical membrane potential and EPSP histogram derivatives (<0) did not reveal clear shifts between their peaks and in right parts indicating tendency to terminate simultaneously (Fig. 10C1,D1). Active state rising time or phase measured at half amplitude of histograms shown in Figure 9 revealed remarkable propagation, whereas decay time measured at half amplitude did not show significant propagation (Fig. 11A). Termination phase measured at half amplitude showed significant differences between frontal and posterior parts of the cortex, but the significance decreased when termination phase was measured at one-tenth amplitude that implies that measured differences depended on waveform (Fig. 11B). Similar to Figure 3E3 and the related analysis, comparison of segregate points of histograms revealed differences in termination (Fig. 11B), but taking into account the whole process (Fig. 10C1,D1), we conclude that the active states from different cortical areas tended to terminate simultaneously. Thus, during the propagation of the leading edge in the anteroposterior direction, active state duration became progressively shorter (Fig. 11A). The active state duration within the V-RS posterior versus the duration within frontal areas was almost twice shorter (p ≤ 0.001; Fig. 11A).

Figure 11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 11.

Duration (A) and proportion of phase (B) of active states within different cortical areas. A, Calculated parameters are shown schematically on top. Zero time was set at −0.5 of phase (Fig. 1). Active state duration was calculated as the time interval between half rising and half decay within EPSP time histograms. The slow-wave cycle duration was calculated as time between two consecutive crossings for decay. Active states become significantly shorter while propagating in the anteroposterior direction independently on the duration of slow-wave cycle. B, Calculated parameters are shown schematically on top. One-tenth decay phase was measured similar to half decay. Proportion of phase occupied by active states was calculated from the interval between the half decay and half rising phases within EPSP phase histograms using basic statistics as in A. *p < 0.05 (t test or Mann–Whitney sum rank test). **p < 0.01 (t test or Mann–Whitney sum rank test). ***p ≤ 0.001 (t test or Mann–Whitney sum rank test). Half, one-tenth decay, and half rising phases were compared using circular statistics. *p < 0.05 (Watson–Williams test). **p < 0.01 (Watson–Williams test). ***p ≤ 0.001 (Watson–Williams test). Data are mean ± SEM (phase duration and time) and mean ± confidence interval stated at the 95% confidence level (half decay and half rising phases).

To verify the results of phase analysis, we performed paired intracellular recordings within the thalamus and the cortex (Fig. 12). The neuron from the anterior thalamus started active states hundreds of milliseconds after the cortical S1 neuron (Fig. 12A), whereas the posterior thalamic neuron (presumably PO) was activated before the cortical S1 neuron recorded simultaneously (Fig. 12B). The results of time analysis obtained with dual intracellular recordings support the findings obtained with phase analysis of neuronal activation.

Figure 12.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 12.

The two examples of paired intracellular recordings revealed different time delays between cortical S1 cells and neurons from the anterior and posterior thalamus. A, The thalamic neuron was recorded within the dorsomedial part of anteroventral nucleus (AVDM). The thalamic active states started almost 200 ms after the cortical active states. B, The thalamic active states started ∼100 ms before the cortical active states. The thalamic neuron was not identified morphologically, but we targeted the posterior thalamus.

Comparison of the inhibitory inputs into projecting VPM and second-order PO nuclei

One of unexpected findings of the study was absence of well-defined inhibitory activities in most second-order thalamic nuclei. It is well documented that the main inhibitory input into first-order thalamic nuclei is provided by the RTN, whereas higher-order thalamus is also innervated with other subcortical nuclei (Barthó et al., 2002; Bokor et al., 2005). Particularly, the PO nucleus receives inhibitory input from the ZI (Trageser and Keller, 2004). As expected, in sensory thalamic nuclei, we found IPSPs phase-locked to the slow oscillation, but we did not find profound IPSPs within higher-order thalamic nuclei during slow oscillations. Therefore, using retrograde labeling by CtB, we investigated differences in the inputs from RTN and ZI into VPM and PO nuclei, respectively. CtB injected in the VPM nucleus (n = 7 mice) retrogradely labeled neurons within the cortical layer VI of the somatosensory cortex as well as neurons in the RTN (Fig. 13A), whereas, when injected into PO (n = 6 mice), CtB labeled neurons mainly within the cortical layer V and in the ZI (Fig. 13B). The labeled cells within the RTN tended to be grouped in clusters, whereas labeled neurons within the ZI tended to stain as individual neurons (Fig. 13). We used an unbiased automatic quantification of all labeled objects (individual neurons or clusters of neurons, area of >40 μm2) from the sections containing the RTN or the ZI (Fig. 13C). To quantify differences between local densities of stained neurons within VPM and ZI nuclei, we calculated the ratio between the integrated areas of clusters (object length >25 μm) to the total area of all labeled objects for each animal. The ratios for RTN and ZI nuclei were 0.74 ± 0.15 versus 0.20 ± 0.09, respectively (p ≤ 0.001), reflecting the staining in clusters of neurons within the RTN versus the staining of individual neurons in the ZI.

Figure 13.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 13.

Comparison of the main inhibitory inputs projecting into first-order VPM and second-order PO nuclei. A, CtB injected in the VPM retrogradely labeled neurons within the cortical layer VI as well as cells in the RTN. B, When CtB was injected into the PO nucleus, it labeled neurons within the cortical layer V and in ZI. Unlike ZI neurons, RTN-labeled cells tend to group together. C, Quantification of labeled object size within ZI and RTN nuclei. Labeled RTN neurons tend to form clusters.

Discussion

In the present study, we show that, in the mouse thalamus, during active phases of the slow oscillation, the inhibition dominates in sensory nuclei, whereas the excitation dominates in higher-order nuclei. Applying phase and time analysis to intracellular and LFP activities within different cortical and thalamic areas, this study is the first to provide information on the global structure of the slow oscillation within the TC system of anesthetized mice. (1) Similar to previous studies (Massimini et al., 2004; Mohajerani et al., 2010; Nir et al., 2011; Stroh et al., 2013), we confirm that most of the slow oscillation cycles started from the frontal cortical areas and the leading edge of the slow waves propagated over the cortex in the anteroposterior direction, whereas the end of active states occurred nearly simultaneously in different cortical sites (Volgushev et al., 2006) (Figs. 3, 7, 9, and 10); (2) the first-order (sensory) thalamic nuclei followed the corresponding cortical areas during the slow oscillation; (3) the higher-order PO and intralaminar PF thalamic neurons started off the slow-wave active states approximately simultaneously; and (4) the latest structure to become involved in the slow-wave cycle was the anterior thalamus, which approximately followed the retrosplenial cortex.

Neurons within the sensory thalamus revealed IPSPs modulated by cortical slow waves as well as large-amplitude, fast-rising EPSPs originating likely from ascending structures and therefore not modulated by the slow oscillation. When “driver” inputs arrived to TC neurons during IPSP, we never observed spikes triggered by the driver, which is consistent with the gating role of the slow oscillation (Timofeev et al., 1996). The inhibitory potentials in VPM and LGN neurons during the slow oscillation originated likely because of burst firing of RTN neurons (Contreras and Steriade, 1995; Timofeev and Steriade, 1996). Instead, the nonsensory thalamus revealed depolarizing active states and therefore was actively involved in the cortical slow oscillation by providing a phase-locked positive feedback and likely contributing to the synchronization of cortical active states.

The two intracellular slow-wave patterns we detected in the sensory versus higher-order thalamus of anesthetized mice could be explained by the patterns of their connectivity. (1) The drivers within sensory nuclei originate in ascending sensory pathways, whereas the drivers in higher-order nuclei originate from layer V pyramidal cells (Groh et al., 2008), which by themselves are the neurons playing a leading role in the generation of slow oscillation (Chauvette et al., 2010; Beltramo et al., 2013). (2) Consistent with the known ratio of large versus small excitatory synapses in different thalamic nuclei (Van Horn and Sherman, 2007), the activities of modulators were much more pronounced within higher-order nuclei than within sensory VPM and LGN nuclei. (3) It appears that the RTN provides a stronger inhibition into sensory thalamic nuclei during cortical active states than the ZI onto higher-order nuclei because a majority of RTN neurons fire bursts of action potentials (Avanzini et al., 1989; Contreras et al., 1992; Fuentealba et al., 2005), but ZI neurons fire tonically (Trageser et al., 2006). Therefore, the postsynaptic action of RTN neurons should be stronger compared with ZI neurons. Also, the RTN neurons form clusters (Fig. 13). It is known that RTN neurons are interconnected by chemical and electrical synapses (Landisman et al., 2002; Shu and McCormick, 2002). Clusters of electrically coupled RTN neurons generate synchronous activities (Long et al., 2004), suggesting that tightly interconnected bursting RTN neurons within clusters would induce stronger IPSPs in their targets than dispersed tonically firing neurons from ZI.

Consistent with morphological observations in primates (Rovó et al., 2012), we did not observe “driver” EPSPs within intralaminar PF neurons, but we observed it in the PO nucleus. The thalamic PO drivers provided a strong relay between the primary and the secondary somatosensory cortices (Theyel et al., 2010), contributing to the long-range cortical synchronization. The “driver” EPSPs within thalamic AV neurons probably originated from the hippocampal system and/or retrosplenial cortex (Shibata, 1998; Petrof and Sherman, 2009).

The neurons within PF, CL, and AD nucleus did not display LTS in response to a hyperpolarizing current pulse, corroborating previous extracellular unit and whole-cell in vitro studies (Lacey et al., 2007; Beatty et al., 2009). Interestingly, the same thalamic neurons in our study did not display driver EPSPs, suggesting their very particular role in the thalamic functions.

One of the intriguing findings of the study was that, on average, the onset of active states in PO and in PF neurons occurred simultaneously with the onset of earliest cortical active states. The onset of depolarizing states appeared as a buildup of synaptic “modulator” events in the PF nucleus or as large-amplitude EPSPs in the PO nucleus (Fig. 8), suggesting its origin in presynaptic structures. Cortical inputs to PF nucleus of rat arrive from primary and secondary motor cortex, primary somatosensory cortex, and parietal cortex (Cornwall and Phillipson, 1988; Van der Werf et al., 2002). The frontal cortical areas as well as the thalamus, forming long-distance and multiple connections, are considered to be hubs within the functional brain network (Honey et al., 2007; van den Heuvel and Sporns, 2011). If one assumes a probabilistic mechanism of slow-wave initiation at any cortical site (Timofeev et al., 2000; Chauvette et al., 2010), and the PF nucleus represents a hub connecting these cortical areas, which integrates multiple inputs from long-distance sites, the PF nucleus would be at least the second to activate during each slow-wave cycle and thus, in relation to all cycles together (average), it would be detected as the first structure. The PF nucleus also receives inputs from laterodorsal and pedunculopontine tegmentum, locus ceruleus, and raphe nuclei, most densely the dorsal raphe (Cornwall and Phillipson, 1988; Van der Werf et al., 2002). During slow-wave sleep, locus ceruleus neurons tend to fire before the onset of active states in the prefrontal cortex (Eschenko et al., 2012); therefore, they can actively contribute to the early onset of active states in PF neurons. Phase-locked firing of other neuromodulatory structures may also contribute to shaping of the slow oscillation.

Despite the fact that the PF neurons demonstrated the early onset of slow-wave active states, their role in the global slow-wave generation is unclear. The majority of projections from the PF nucleus reach striatum (Deschênes et al., 1996; Van der Werf et al., 2002). The cortical projections of these neurons form small terminal puffs in layer VI of the frontal motor cortex (Deschênes et al., 1996). These data indicate that PF neurons can boost the slow oscillation mainly in frontal areas.

Neurons of the PO nucleus receive ascending inputs from somatosensory, auditory, visual, and vestibular systems (Groenewegen and Witter, 2004), which should not be synchronized with the slow oscillation. The cortical (driver) inputs arrive mainly from somatosensory (S1 and S2) areas, but to a lesser extent from medial frontal cortices (Veinante et al., 2000; Guandalini, 2001); and frontal cortex is the first to display the onset of slow waves, which potentially could drive these neurons and show the early onset of activity in the PO nucleus. However, the density of corticothalamic terminals originating from the frontal cortex to PO and AV nuclei is similar (Guandalini, 2001); the PO nucleus was among the first to enter into active phases of the slow oscillation and the AV nucleus was the last one (Fig. 10).

This and other studies (Massimini et al., 2004; Volgushev et al., 2006; Ruiz-Mejias et al., 2011) show that cortical active state onsets propagate in the anteroposterior direction as well as demonstrate simultaneous termination of active states across neocortex. Although the exact mechanisms of this simultaneous termination remain unknown, several investigations point to a possible role of inhibitory processes in such termination (Steriade et al., 1993; Puig et al., 2008; Chen et al., 2012). Active contribution of higher-order thalamocortical neurons in active state maintenance suggests that some of them can provide a small, but widespread, influence over distant cortical areas and therefore control termination of active states.

In this study using stereotypic global slow-wave pattern in ketamine-xylazine anesthetized mice, we show that the active phases of slow waves occur first in frontal cortical areas and in second-order (posterior) thalamic nuclei (PF and PO). The slow waves propagate in the rostrocaudal direction in the cortex (from frontal to retrosplenial cortex). The sequence of detection of slow waves in the thalamus was in the caudorostral direction (from posterior PF and PO nuclei to the anterior AV nucleus; Fig. 10F). The AV nucleus forms reciprocal connections with the limbic system and, in particular, with the retrosplenial cortex (Groenewegen and Witter, 2004) as well as the medial prefrontal cortex (Guandalini, 2001), potentially making a link between cortical areas with the earliest and the latest involvement in the slow oscillation. Our results demonstrate that only the first-order nuclei might contribute to the sleep-related gating of peripheral inputs, whereas higher-order thalamic nuclei may facilitate an intracortical dialog during sleep, which is possibly a key component of sleep-dependent memory consolidation (Rasch and Born, 2013).

Footnotes

  • This work was supported by National Science and Engineering Research Council of Canada Grant 298475, Canadian Institutes of Health Research Grants MOP-37862 and MOP-67175, and National Institute of Neurological Disorders and Stroke Grants 1R01-NS060870 and 1R01-NS059740. We thank Serge Ftomov for technical assistance, Josée Seigneur for management assistance, and Sylvain Chauvette for comments on the manuscript.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Dr. Igor Timofeev, Centre de recherche de l'Institut universitaire en santé mentale de Québec, 2601 de la Canardière Québec, Quebec G1J 2G3, Canada. Igor.Timofeev{at}phs.ulaval.ca

References

  1. ↵
    1. Avanzini G,
    2. de Curtis M,
    3. Panzica F,
    4. Spreafico R
    (1989) Intrinsic properties of nucleus reticularis thalami neurones of the rat studied in vitro. J Physiol 416:111–122, pmid:2558172.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Barthó P,
    2. Freund TF,
    3. Acsády L
    (2002) Selective gabaergic innervation of thalamic nuclei from zona incerta. Eur J Neurosci 16:999–1014, doi:10.1046/j.1460-9568.2002.02157.x, pmid:12383229.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Beatty JA,
    2. Sylwestrak EL,
    3. Cox CL
    (2009) Two distinct populations of projection neurons in the rat lateral parafascicular thalamic nucleus and their cholinergic responsiveness. Neuroscience 162:155–173, doi:10.1016/j.neuroscience.2009.04.043, pmid:19393292.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Beltramo R,
    2. D'Urso G,
    3. Dal Maschio M,
    4. Farisello P,
    5. Bovetti S,
    6. Clovis Y,
    7. Lassi G,
    8. Tucci V,
    9. De Pietri Tonelli D,
    10. Fellin T
    (2013) Layer-specific excitatory circuits differentially control recurrent network dynamics in the neocortex. Nat Neurosci 16:227–234, doi:10.1038/nn.3306, pmid:23313909.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bokor H,
    2. Frère SG,
    3. Eyre MD,
    4. Slézia A,
    5. Ulbert I,
    6. Lüthi A,
    7. Acsády L
    (2005) Selective GABAergic control of higher-order thalamic relays. Neuron 45:929–940, doi:10.1016/j.neuron.2005.01.048, pmid:15797553.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Bourassa J,
    2. Pinault D,
    3. Deschênes M
    (1995) Corticothalamic projections from the cortical barrel field to the somatosensory thalamus in rats: a single-fibre study using biocytin as an anterograde tracer. Eur J Neurosci 7:19–30, pmid:7711933.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Chauvette S,
    2. Volgushev M,
    3. Timofeev I
    (2010) Origin of active states in local neocortical networks during slow sleep oscillation. Cereb Cortex 20:2660–2674, doi:10.1093/cercor/bhq009, pmid:20200108.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Chauvette S,
    2. Crochet S,
    3. Volgushev M,
    4. Timofeev I
    (2011) Properties of slow oscillation during slow-wave sleep and anesthesia in cats. J Neurosci 31:14998–15008, doi:10.1523/JNEUROSCI.2339-11.2011, pmid:22016533.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Chen JY,
    2. Chauvette S,
    3. Skorheim S,
    4. Timofeev I,
    5. Bazhenov M
    (2012) Interneuron-mediated inhibition synchronizes neuronal activity during slow oscillation. J Physiol 590:3987–4010, doi:10.1113/jphysiol.2012.227462, pmid:22641778.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Constantinople CM,
    2. Bruno RM
    (2013) Deep cortical layers are activated directly by thalamus. Science 340:1591–1594, doi:10.1126/science.1236425, pmid:23812718.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Contreras D,
    2. Steriade M
    (1995) Cellular basis of EEG slow rhythms: a study of dynamic corticothalamic relationships. J Neurosci 15:604–622, pmid:7823167.
    OpenUrlAbstract
  12. ↵
    1. Contreras D,
    2. Curró Dossi R,
    3. Steriade M
    (1992) Bursting and tonic discharges in two classes of reticular thalamic neurons. J Neurophysiol 68:973–977, pmid:1432063.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Contreras D,
    2. Timofeev I,
    3. Steriade M
    (1996) Mechanisms of long-lasting hyperpolarizations underlying slow sleep oscillations in cat corticothalamic networks. J Physiol 494:251–264, pmid:8814619.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Cornwall J,
    2. Phillipson OT
    (1988) Afferent projections to the parafascicular thalamic nucleus of the rat, as shown by the retrograde transport of wheat germ agglutinin. Brain Res Bull 20:139–150, doi:10.1016/0361-9230(88)90171-2, pmid:2836036.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Cossart R,
    2. Aronov D,
    3. Yuste R
    (2003) Attractor dynamics of network up states in the neocortex. Nature 423:283–288, doi:10.1038/nature01614, pmid:12748641.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Crunelli V,
    2. Hughes SW
    (2010) The slow (<1 Hz) rhythm of non-REM sleep: a dialogue between three cardinal oscillators. Nat Neurosci 13:9–17, doi:10.1038/nn.2445, pmid:19966841.
    OpenUrlCrossRefPubMed
  17. ↵
    1. David F,
    2. Schmiedt JT,
    3. Taylor HL,
    4. Orban G,
    5. Di Giovanni G,
    6. Uebele VN,
    7. Renger JJ,
    8. Lambert RC,
    9. Leresche N,
    10. Crunelli V
    (2013) Essential thalamic contribution to slow waves of natural sleep. J Neurosci 33:19599–19610, doi:10.1523/JNEUROSCI.3169-13.2013, pmid:24336724.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Deschênes M,
    2. Bourassa J,
    3. Pinault D
    (1994) Corticothalamic projections from layer V cells in rat are collaterals of long-range corticofugal axons. Brain Res 664:215–219, pmid:7895031.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Deschênes M,
    2. Bourassa J,
    3. Doan VD,
    4. Parent A
    (1996) A single-cell study of the axonal projections arising from the posterior intralaminar thalamic nuclei in the rat. Eur J Neurosci 8:329–343, doi:10.1111/j.1460-9568.1996.tb01217.x, pmid:8714704.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Eschenko O,
    2. Magri C,
    3. Panzeri S,
    4. Sara SJ
    (2012) Noradrenergic neurons of the locus coeruleus are phase locked to cortical up-down states during sleep. Cereb Cortex 22:426–435, doi:10.1093/cercor/bhr121, pmid:21670101.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Fuentealba P,
    2. Timofeev I,
    3. Bazhenov M,
    4. Sejnowski TJ,
    5. Steriade M
    (2005) Membrane bistability in thalamic reticular neurons during spindle oscillations. J Neurophysiol 93:294–304, pmid:15331618.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Groenewegen HJ,
    2. Witter MP
    (2004) in The rat nervous system, Thalamus, ed Paxinos G (Elsevier Academic, New York), Ed 3, pp 407–453.
  23. ↵
    1. Groh A,
    2. de Kock CP,
    3. Wimmer VC,
    4. Sakmann B,
    5. Kuner T
    (2008) Driver or coincidence detector: modal switch of a corticothalamic giant synapse controlled by spontaneous activity and short-term depression. J Neurosci 28:9652–9663, pmid:18815251.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Guandalini P
    (2001) The efferent connections to the thalamus and brainstem of the physiologically defined eye field in the rat medial frontal cortex. Brain Res Bull 54:175–186, doi:10.1016/S0361-9230(00)00444-5, pmid:11275407.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Honey CJ,
    2. Kötter R,
    3. Breakspear M,
    4. Sporns O
    (2007) Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc Natl Acad Sci U S A 104:10240–10245, doi:10.1073/pnas.0701519104, pmid:17548818.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Hughes SW,
    2. Cope DW,
    3. Blethyn KL,
    4. Crunelli V
    (2002) Cellular mechanisms of the slow (<1 Hz) oscillation in thalamocortical neurons in vitro. Neuron 33:947–958, doi:10.1016/S0896-6273(02)00623-2, pmid:11906700.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Jones EG
    (2009) Synchrony in the interconnected circuitry of the thalamus and cerebral cortex. Ann N Y Acad Sci 1157:10–23, doi:10.1111/j.1749-6632.2009.04534.x, pmid:19351352.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Lacey CJ,
    2. Bolam JP,
    3. Magill PJ
    (2007) Novel and distinct operational principles of intralaminar thalamic neurons and their striatal projections. J Neurosci 27:4374–4384, doi:10.1523/JNEUROSCI.5519-06.2007, pmid:17442822.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Landisman CE,
    2. Long MA,
    3. Beierlein M,
    4. Deans MR,
    5. Paul DL,
    6. Connors BW
    (2002) Electrical synapses in the thalamic reticular nucleus. J Neurosci 22:1002–1009, pmid:11826128.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Lemieux M,
    2. Chen JY,
    3. Lonjers P,
    4. Bazhenov M,
    5. Timofeev I
    (2014) The impact of cortical deafferentation on the neocortical slow oscillation. J Neurosci 34:5689–5703, doi:10.1523/JNEUROSCI.1156-13.2014, pmid:24741059.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Li J,
    2. Guido W,
    3. Bickford ME
    (2003) Two distinct types of corticothalamic EPSPs and their contribution to short-term synaptic plasticity. J Neurophysiol 90:3429–3440, doi:10.1152/jn.00456.2003, pmid:12890796.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Long MA,
    2. Landisman CE,
    3. Connors BW
    (2004) Small clusters of electrically coupled neurons generate synchronous rhythms in the thalamic reticular nucleus. J Neurosci 24:341–349, doi:10.1523/JNEUROSCI.3358-03.2004, pmid:14724232.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Mallet N,
    2. Pogosyan A,
    3. Márton LF,
    4. Bolam JP,
    5. Brown P,
    6. Magill PJ
    (2008) Parkinsonian β oscillations in the external globus pallidus and their relationship with subthalamic nucleus activity. J Neurosci 28:14245–14258, doi:10.1523/JNEUROSCI.4199-08.2008, pmid:19109506.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Massimini M,
    2. Huber R,
    3. Ferrarelli F,
    4. Hill S,
    5. Tononi G
    (2004) The sleep slow oscillation as a traveling wave. J Neurosci 24:6862–6870, doi:10.1523/JNEUROSCI.1318-04.2004, pmid:15295020.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Mohajerani MH,
    2. McVea DA,
    3. Fingas M,
    4. Murphy TH
    (2010) Mirrored bilateral slow-wave cortical activity within local circuits revealed by fast bihemispheric voltage-sensitive dye imaging in anesthetized and awake mice. J Neurosci 30:3745–3751, doi:10.1523/JNEUROSCI.6437-09.2010, pmid:20220008.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Moore BR
    (1980) A modification of the Rayleigh Test for vector data. Biometrika 67:175–180, doi:10.1093/biomet/67.1.175.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Mukovski M,
    2. Chauvette S,
    3. Timofeev I,
    4. Volgushev M
    (2007) Detection of active and silent states in neocortical neurons from the field potential signal during slow-wave sleep. Cereb Cortex 17:400–414, doi:10.1093/cercor/bhj157, pmid:16547348.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Nir Y,
    2. Staba RJ,
    3. Andrillon T,
    4. Vyazovskiy VV,
    5. Cirelli C,
    6. Fried I,
    7. Tononi G
    (2011) Regional slow waves and spindles in human sleep. Neuron 70:153–169, doi:10.1016/j.neuron.2011.02.043, pmid:21482364.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Paxinos G,
    2. Franklin K
    (2001) The mouse brain in stereotaxic coordinates (Academic, New York), Ed 2.
  40. ↵
    1. Petrof I,
    2. Sherman SM
    (2009) Synaptic properties of the mammillary and cortical afferents to the anterodorsal thalamic nucleus in the mouse. J Neurosci 29:7815–7819, pmid:19535593.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Puig MV,
    2. Ushimaru M,
    3. Kawaguchi Y
    (2008) Two distinct activity patterns of fast-spiking interneurons during neocortical up states. Proc Natl Acad Sci U S A 105:8428–8433, doi:10.1073/pnas.0712219105, pmid:18550841.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Rasch B,
    2. Born J
    (2013) About sleep's role in memory. Physiol Rev 93:681–766, pmid:23589831.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Reichova I,
    2. Sherman SM
    (2004) Somatosensory corticothalamic projections: distinguishing drivers from modulators. J Neurophysiol 92:2185–2197, doi:10.1152/jn.00322.2004, pmid:15140908.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Rovó Z,
    2. Ulbert I,
    3. Acsády L
    (2012) Drivers of the primate thalamus. J Neurosci 32:17894–17908, doi:10.1523/JNEUROSCI.2815-12.2012, pmid:23223308.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Ruiz-Mejias M,
    2. Ciria-Suarez L,
    3. Mattia M,
    4. Sanchez-Vives MV
    (2011) Slow and fast rhythms generated in the cerebral cortex of the anesthetized mouse. J Neurophysiol 106:2910–2921, doi:10.1152/jn.00440.2011, pmid:21880935.
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Sanchez-Vives MV,
    2. McCormick DA
    (2000) Cellular and network mechanisms of rhythmic recurrent activity in neocortex. Nat Neurosci 3:1027–1034, doi:10.1038/79848, pmid:11017176.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Sharma AV,
    2. Wolansky T,
    3. Dickson CT
    (2010) A comparison of sleeplike slow oscillations in the hippocampus under ketamine and urethane anesthesia. J Neurophysiol 104:932–939, doi:10.1152/jn.01065.2009, pmid:20538775.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Sherman SM,
    2. Guillery RW
    (1996) Functional organization of thalamocortical relays. J Neurophysiol 76:1367–1395, pmid:8890259.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Sherman SM,
    2. Guillery RW
    (1998) On the actions that one nerve cell can have on another: distinguishing “drivers” from “modulators.” Proc Natl Acad Sci U S A 95:7121–7126.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Shibata H
    (1998) Organization of projections of rat retrosplenial cortex to the anterior thalamic nuclei. Eur J Neurosci 10:3210–3219, doi:10.1046/j.1460-9568.1998.00328.x, pmid:9786214.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Shu Y,
    2. McCormick DA
    (2002) Inhibitory interactions between ferret thalamic reticular neurons. J Neurophysiol 87:2571–2576, pmid:11976393.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Siapas AG,
    2. Lubenov EV,
    3. Wilson MA
    (2005) Prefrontal phase locking to hippocampal θ oscillations. Neuron 46:141–151, doi:10.1016/j.neuron.2005.02.028, pmid:15820700.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Slézia A,
    2. Hangya B,
    3. Ulbert I,
    4. Acsády L
    (2011) Phase advancement and nucleus-specific timing of thalamocortical activity during slow cortical oscillation. J Neurosci 31:607–617, doi:10.1523/JNEUROSCI.3375-10.2011, pmid:21228169.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Steriade M,
    2. Nuñez A,
    3. Amzica F
    (1993) A novel slow (<1 Hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components. J Neurosci 13:3252–3265, pmid:8340806.
    OpenUrlAbstract
  55. ↵
    1. Steriade M,
    2. Timofeev I,
    3. Grenier F
    (2001) Natural waking and sleep states: a view from inside neocortical neurons. J Neurophysiol 85:1969–1985, pmid:11353014.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Stroh A,
    2. Adelsberger H,
    3. Groh A,
    4. Rühlmann C,
    5. Fischer S,
    6. Schierloh A,
    7. Deisseroth K,
    8. Konnerth A
    (2013) Making waves: initiation and propagation of corticothalamic Ca2+ waves in vivo. Neuron 77:1136–1150, doi:10.1016/j.neuron.2013.01.031, pmid:23522048.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Theyel BB,
    2. Llano DA,
    3. Sherman SM
    (2010) The corticothalamocortical circuit drives higher-order cortex in the mouse. Nat Neurosci 13:84–88, doi:10.1038/nn.2449, pmid:19966840.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Timofeev I,
    2. Steriade M
    (1996) Low-frequency rhythms in the thalamus of intact-cortex and decorticated cats. J Neurophysiol 76:4152–4168, pmid:8985908.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Timofeev I,
    2. Contreras D,
    3. Steriade M
    (1996) Synaptic responsiveness of cortical and thalamic neurones during various phases of slow sleep oscillation in cat. J Physiol 494:265–278, pmid:8814620.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Timofeev I,
    2. Grenier F,
    3. Bazhenov M,
    4. Sejnowski TJ,
    5. Steriade M
    (2000) Origin of slow cortical oscillations in deafferented cortical slabs. Cereb Cortex 10:1185–1199, doi:10.1093/cercor/10.12.1185, pmid:11073868.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Timofeev I,
    2. Grenier F,
    3. Steriade M
    (2001) Disfacilitation and active inhibition in the neocortex during the natural sleep-wake cycle: an intracellular study. Proc Natl Acad Sci U S A 98:1924–1929, doi:10.1073/pnas.98.4.1924, pmid:11172052.
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Trageser JC,
    2. Keller A
    (2004) Reducing the uncertainty: gating of peripheral inputs by zona incerta. J Neurosci 24:8911–8915, doi:10.1523/JNEUROSCI.3218-04.2004, pmid:15470158.
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Trageser JC,
    2. Burke KA,
    3. Masri R,
    4. Li Y,
    5. Sellers L,
    6. Keller A
    (2006) State-dependent gating of sensory inputs by zona incerta. J Neurophysiol 96:1456–1463, doi:10.1152/jn.00423.2006, pmid:16775205.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. van den Heuvel MP,
    2. Sporns O
    (2011) Rich-club organization of the human connectome. J Neurosci 31:15775–15786, doi:10.1523/JNEUROSCI.3539-11.2011, pmid:22049421.
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. Van der Werf YD,
    2. Witter MP,
    3. Groenewegen HJ
    (2002) The intralaminar and midline nuclei of the thalamus: anatomical and functional evidence for participation in processes of arousal and awareness. Brain Res Rev 39:107–140, doi:10.1016/S0165-0173(02)00181-9, pmid:12423763.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Van Horn SC,
    2. Sherman SM
    (2007) Fewer driver synapses in higher order than in first order thalamic relays. Neuroscience 146:463–470, doi:10.1016/j.neuroscience.2007.01.026, pmid:17320295.
    OpenUrlCrossRefPubMed
  67. ↵
    1. Veinante P,
    2. Lavallée P,
    3. Deschênes M
    (2000) Corticothalamic projections from layer 5 of the vibrissal barrel cortex in the rat. J Comp Neurol 424:197–204, doi:10.1002/1096-9861(20000821)424:2<197::AID-CNE1>3.0.CO%3B2-6, pmid:10906697.
    OpenUrlCrossRefPubMed
  68. ↵
    1. Volgushev M,
    2. Chauvette S,
    3. Mukovski M,
    4. Timofeev I
    (2006) Precise long-range synchronization of activity and silence in neocortical neurons during slow-wave sleep. J Neurosci 26:5665–5672, doi:10.1523/JNEUROSCI.0279-06.2006, pmid:16723523.
    OpenUrlAbstract/FREE Full Text
View Abstract
Back to top

In this issue

The Journal of Neuroscience: 34 (26)
Journal of Neuroscience
Vol. 34, Issue 26
25 Jun 2014
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Global Intracellular Slow-Wave Dynamics of the Thalamocortical System
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Global Intracellular Slow-Wave Dynamics of the Thalamocortical System
Maxim Sheroziya, Igor Timofeev
Journal of Neuroscience 25 June 2014, 34 (26) 8875-8893; DOI: 10.1523/JNEUROSCI.4460-13.2014

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Global Intracellular Slow-Wave Dynamics of the Thalamocortical System
Maxim Sheroziya, Igor Timofeev
Journal of Neuroscience 25 June 2014, 34 (26) 8875-8893; DOI: 10.1523/JNEUROSCI.4460-13.2014
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • cortex
  • dynamics
  • intracellular
  • slow oscillation
  • thalamus

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Memory Retrieval Has a Dynamic Influence on the Maintenance Mechanisms That Are Sensitive to ζ-Inhibitory Peptide (ZIP)
  • Neurophysiological Evidence for a Cortical Contribution to the Wakefulness-Related Drive to Breathe Explaining Hypocapnia-Resistant Ventilation in Humans
  • Monomeric Alpha-Synuclein Exerts a Physiological Role on Brain ATP Synthase
Show more Articles

Systems/Circuits

  • Oddball Evoked Deviant Responses Reflect Complex Context-Dependent Expectations in Mouse V1
  • V2b Neurons Act via Multiple Targets to Produce in Phase Inhibition during Locomotion
  • Exploring Neural Dynamics in the Auditory Telencephalon of Crows Using Functional Ultrasound Imaging
Show more Systems/Circuits
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.