- Research article
- Open Access
Wave propagation of cortical population activity under urethane anesthesia is state dependent
BMC Neurosciencevolume 14, Article number: 78 (2013)
Propagating waves of excitation have been observed extensively in the neocortex, during both spontaneous and sensory-evoked activity, and they play a critical role in spatially organizing information processing. However, the state-dependence of these spatiotemporal propagation patterns is largely unexplored. In this report, we use voltage-sensitive dye imaging in the rat visual cortex to study the propagation of spontaneous population activity in two discrete cortical states induced by urethane anesthesia.
While laminar current source density patterns of spontaneous population events in these two states indicate a considerable degree of similarity in laminar networks, lateral propagation in the more active desynchronized state is approximately 20% faster than in the slower synchronized state. Furthermore, trajectories of wave propagation exhibit a strong anisotropy, but the preferred direction is different depending on cortical state.
Our results show that horizontal wave propagation of spontaneous neural activity is largely dependent on the global activity states of local cortical circuits.
Spontaneous and sensory-evoked population activity in the mammalian neocortex is known to manifest in the form of traveling or propagating waves (e.g. [1–8]). Among other functions, these waves are postulated to provide nonspecific background depolarization to neuronal assemblies, thus modulating the likelihood of action potential initiation within spatially segregated populations of cells [9, 10]. Despite this proposed role of propagating waves in spatially coordinating neuronal computation [11–15], the mechanisms governing wave propagation remain an issue of ongoing debate . In particular, little is known about how propagation dynamics are influenced by global network states .
Distinct global network states are a critical feature of neocortical dynamics, and characterize sleep and wakefulness , as well as stages of arousal and anesthesia . These global network states can have profound effects on cortical information processing by biasing the dynamics of membrane potential synchrony  and the interaction of sensory evoked responses with ongoing activity [21, 22]. Thus, changes in cortical state should lead to pronounced changes in the operation of cortical networks underlying wave propagation, but investigations to date have mainly focused on wave propagation in single states . Therefore, in order to comprehensively account for the mechanisms of cortical wave propagation and to understand differences in observations from different experimental and clinical preparations, insight into the state-dependence of propagation patterns is crucial.
Here, we address the state-dependence of lateral intracortical networks by analyzing spatiotemporal patterns of propagating population activity as recorded with voltage-sensitive dye (VSD) imaging in rat visual cortex, under a high-amplitude ECoG-synchronized state and a low-amplitude ECoG-desynchronized state induced by urethane anesthesia. Under moderate planes of urethane anesthesia, the ECoG (electrocorticogram) spontaneously alternates between these two states (Figure 1A), which have been shown to resemble, in some respects, forebrain states observed during sleep [23, 24]. We first used laminar current source density (CSD) analysis to investigate the cortical laminar organization of the spontaneous population activity patterns in each state. Using a temporospatial flow detection algorithm , we then quantified local propagation patterns of spontaneous activity in single trials. Our results show that, despite a large degree of similarity in laminar networks, propagation of spontaneous excitable waves is very susceptible to cortical state.
All experiments were in compliance with the guidelines of the European Community (EUVD 86/609/EEC) and were approved by an ethics commission of the state of Sachsen-Anhalt. Twelve adult male Wistar rats (250–400 g) were used for the experiments described here.
Induction of anesthesia was performed by either intraperitoneal or intravenous infusion of an aqueous solution of urethane (1.25–1.5 g/kg). Further anesthetic was given as necessary to maintain areflexia.
To analyze the layered activity distribution in both cortical states, three animals were implanted with a silicon microelectrode (32 recording sites spaced 50 μm apart, each 400 μm2, NeuroNexus Technologies, Inc.) into the primary visual cortex (V1). The back surface of the electrode was coated with the fluorescent carbocyanine dye DiIC18 to later localize the recording site in histological slices . Laminar positioning of the electrode was controlled for by inserting it 1.6 mm below the brain surface. Cortical layers were delineated based both on histology and on the typical laminar CSD patterns following visual stimulation . Electric field potentials were recorded with a biosignal recording system (MAP System, Plexon Inc.) and split into local field potentials (LFP, 1–150 Hz) and multi-unit activity (MUA, 0.9–8.8 kHz). LFPs were digitized continuously, and triggered spike data were saved as waveforms and timestamps for offline analysis.
CSDs were obtained from LFP signals as described previously [26, 27]. The calculation of the characteristic CSD profile for both states was done within a window of 1000 ms. In the synchronized state, which continuously cycles between phases of generalized network activity and network silence , the window was centered on each time point where the integrated MUA (binned at 5 ms) crossed a threshold of three standard deviations above baseline (i.e. the preceding bin window) and remained higher for at least 500 ms. We found this triggering to be considerably more accurate than triggering off of the low-pass filtered LFP. For the desynchronized state, thresholding based on MUA cannot be used, since neurons are constantly discharging . We therefore filtered the LFP in a band of 3–6 Hz, which typically shows a spectral peak in the desynchronized state (Figure 1C). After band-pass filtering, we calculated the Hilbert-phase of the signal and triggered the average calculation at zero degree phase angle. This procedure results in both cases in an average CSD profile of thousands of events, locked to either down-up state transitions or phase crossings. Averaged rectified CSDs (AVRECs,  were calculated from 200 ms windows around the trigger.
For VSD imaging, methods were adapted from our previous work [3, 29]. Briefly, craniotomy of approximately 5 mm diameter was performed over the visual cortex of the left hemisphere in nine animals (center: bregma: 7 mm, lateral: 4 mm) and cortex was stained transdurally with the dye RH-1691 (1 mg/ml in normal saline; Optical Imaging) for 1.5 h. Optical signals were recorded and processed as described previously [3, 29].
ECoGs were recorded simultaneously with a 1 mm Ag/AgCl ball electrode placed on the edge of the exposed cortical surface in the lower left corner of the imaging field (plan view, rostrocaudal axis). A second Ag/AgCl electrode placed over the cerebellum served as reference. Signals were amplified with a custom-made 1000-gain amplifier and filtered in a band of 0.5–100 Hz, prior to digitization. The electrocardiogram was recorded as well, to allow triggered removal of heartbeat-related artifacts from the optical recording .
Classification of cortical state
In order to perform offline classification of cortical states, root mean square (RMS) values of 2 s segments were determined from filtered LFP or ECoG recordings. The resulting bimodal distribution was then divided at the saddle point to classify the segment into distinct brain states (Figure 1B). Segments with an RMS value larger than criterion were classified as synchronized state, whereas segments with values smaller than criterion were classified as desynchronized state. Power-spectral densities were estimated from 2 s long, Hamming-windowed segments, and normalized by the integral of power over frequency within the specified window. Data were processed in Matlab (Mathworks) and Mathematica (Wolfram Research).
On average, episodes of each cortical state lasted between 2–10 min (Figure 1D). However, some episodes of the synchronized state lasted up to 50 min (Figure 1D). Such persistent episodes of the urethane synchronized state were found to be particularly prevalent during the first few hours of each recording session, when anesthesia levels are presumably at their deepest. At later stages in the experiment, alternations between synchronized and desynchronized episodes became more frequent and the average duration of each discrete state shortened (data not shown).
Quantification of propagation patterns
In order to characterize the lateral propagation velocity and trajectories of spontaneous activity in VSD recordings, we applied a temporospatial flow detection algorithm . Briefly, this algorithm compares signal traces from neighboring detectors and extracts the time shift necessary to obtain best correlation coefficients between signal traces. The algorithm then aggregates the spatial pattern of these time shifts from neighboring detector pairs, to obtain the propagation pattern within the local area. Although this algorithm can detect the three known patterns of spatial propagation , translation waves, source waves and spiral waves, in our experiments, translation waves were predominant and therefore the other patterns were not analyzed in further detail. Software to perform these calculations is freely available from (http://www.sourceforge.net/projects/nounou). The performance characteristics and implementation details of this algorithm are discussed in detail elsewhere .
In the present study, a correlation window of 200 ms (320 frames) in combination with a maximum shift window of 25 ms (40 frames) was used to calculate temporospatial flow; correlation windows of 150 and 250 ms showed similar results (see Additional file 1). These windows are sufficiently large to encompass the main features of the activity pattern, while not so large as to increase the risk of aliasing across cycles of oscillation. Of note, given the subsequent spatial pattern matching in our flow algorithm, such aliasing would be rejected from the dataset and would therefore lead to lower sensitivity of the algorithm, but would not bias the results. For detection of spontaneous flow events, the algorithm was applied to signal traces from a hexagonal ring of detectors centered in V1 as indicated in Figure 2A. The correlation window was moved in successive steps of 4 frames over the bandpass-filtered (0.5–35 Hz) VSD signal, which consisted of 2 s long traces of spontaneous activity from each detector. Prior to flow detection, VSD traces were classified as synchronized or desynchronized based on the criterion described above. The obtained flow vectors were than matched to templates and detected events with flow reliability lower than 0.85 were discarded . Flow trajectory distributions were tested for statistically significant differences using a circular Kuiper two-sample test from the CircStat toolbox .
We investigated the state-dependence of spontaneous laterally propagating waves in the urethane-anesthetized rat visual cortex. We first examined the laminar organization of these events using CSD analysis of silicon microelectrode recordings to investigate the network population activity associated with each state. We then used VSD imaging to record spontaneous activity in both states, and compared the characteristics of spatiotemporal propagation, i.e. propagation velocities and trajectories.
Laminar networks of spontaneous population activity are similar for synchronized and desynchronized states
Desynchronized forebrain states are characterized by higher spontaneous firing rates of thalamic neurons than synchronized slow-wave states . Therefore, it could be hypothesized that wave propagation in the desynchronized state may potentially involve a stronger thalamic component, whereas propagation patterns in synchronized states are thought to be mainly governed by intracortical inputs [10, 16]. Such a network difference would obscure comparison of propagation patterns across states. As one means to address this concern and to investigate which laminar cortical layers are most involved in these spontaneous activity patterns, we used laminar CSD analysis .
For the synchronized state, CSD profiles were obtained by averaging spontaneous events centered on the transition from silent to active network states. For the desynchronized state, averaging was triggered off of the instantaneous phase of its main oscillatory band. In all three animals examined, the resulting CSD profiles were qualitatively similar across states, in terms of the laminar arrangement and the relative amplitudes of respective current sinks and sources. In particular, both states featured characteristic supragranular sink-source dipoles, the amplitude of which was considerably larger compared to respective granular or infragranular components (Figure 1E). The variance of the total AVREC which is explained by supragranular sink-source dipoles was 68.4% (± 4.2%, SEM) in the synchronized state and 69.9% (± 2.2%, SEM) in the desynchronized state, respectively.
Since laminar CSD profiles of synchronized and desynchronized states show no obvious discrepancies, we interpret this as supportive of the proposition that comparison of propagation patterns in these two states is meaningful with respect to the laminar networks involved.
Horizontal propagation of spontaneous waves is faster in the desynchronized state
Spontaneous propagating waves of depolarization were observed in both cortical states. Waves in the desynchronized state, however, appeared more fragmented spatially (Figure 2B), presumably reflecting the decreased spatial and temporal correlation of membrane potentials in more active forebrain states [20, 33].
Local propagation speeds of spontaneous waves varied greatly across events and trials, apparently following Poisson statistics (Figure 2C), with desynchronized probability distributions being significantly biased towards faster speeds. A comparison of the medians of the respective speed distributions from nine animals showed that, on average, propagation speeds were about 20% faster in the desynchronized state as compared to the synchronized state (Figure 2D, meanSYNC = 98 mm/s ± 7 mm/s and meanDESYNC = 118 mm/s ± 9 mm/s (± SEM), Wilcoxon matched-pairs signed-ranks test, p < 0.02).
Trajectories of spontaneous cortical waves are state-dependent
As for the directional properties of propagation, spontaneous waves measured at the center of V1 traveled in every direction (Figure 2E), but with a highly anisotropic distribution in both states (Figure 2E). While the dominant axis of propagation was always the lateral-medial axis regardless of cortical state, preferred propagation directions within this biased axis were clearly state dependent (Kuiper’s test, p < 0.001 for every animal), differing by 180° for synchronized and desynchronized states. Notably, the respective directions of anisotropy were very stable across animals (Figure 2E), suggesting commonalities in the network organization of these spontaneous waves.
Thus, both the distribution of propagation speeds and of propagation directions showed a clear dependence on the global cortical state.
Based on laminar CSD analyses of the predominant ongoing activity patterns in each state, we first demonstrate that laminar networks in both states demonstrate a considerable degree of similarity. Despite these similarities in network patterns, spontaneous population activity in synchronized and desynchronized states differs, by definition, with respect to its spectral amplitude, as well as in MUA activity patterns , cell-type-specific firing patterns [34, 35] and membrane potential dynamics [20, 36]. In this report, we further demonstrate that this activity also differs in its temporal and spatial characteristics: namely, spontaneous waves of depolarization as recorded by VSD imaging in rat visual cortex tend to propagate significantly faster and with different anisotropies in an activated globally desynchronized state as compared to a globally synchronized slow-wave state.
Waves in synchronized slow-wave states are known to reflect the globalized spread of cellular up states of depolarization across the cortical network [37–39]. In contrast, wave propagation in the urethane-desynchronized state most likely reflects the spread of more localized patterns of oscillatory synchrony (cf. ) with a spectral peak around 4 Hz (Figure 1C). Despite these apparent spatiotemporal differences of the propagating waves, laminar CSD profiles of synchronized and desynchronized states appear strikingly similar (Figure 1E).
In particular, the dipolar arrangement of sink-source pairs in supragranular layers, which was observed across states, indicates that the bulk of synaptic transmission is mediated via projections to somata and/or neural processes located in supragranular layers. Mechanistically, the observed supragranular dipoles (cf. ) presumably reflect postsynaptic input to somata and dendrites of layer 2/3 pyramidal neurons as well as to the distal dendrites of large layer 5 pyramidal neurons which ramify extensively into superficial layers . After being depolarized to firing threshold, pyramidal neurons in layers 2/3 and 5 would then pass on the excitation to their efferent target neurons via intra- and interlaminar projections , resulting in lateral spread of the depolarization as observed in supragranular and infragranular layers of cortical slices [44, 45].
Of note, as there is accumulating evidence suggesting a leading role for layer 5 networks in cortical up state initiation [46–50] and intracortical wave propagation [39, 45, 50], the left panel in Figure 1E should not be misinterpreted as evidence for a supragranular origin of propagating cortical up states. Rather, and within the mechanistic framework described above, we reckon that the prominent supragranular dipoles in Figure 1E are evoked by postsynaptic input coming from infragranular neurons which project onto somata and neural processes in layer 2/3. This scenario is supported by a recent study showing that layer 5A provides strong driving input to supragranular layers .
Regardless of exact circuit mechanisms, our laminar CSD patterns do argue against a fundamental functional network difference in the two states, as far as the relative degree of thalamic contribution to the respective current sinks and sources is concerned. For instance, if the contribution of thalamic components to spontaneous population activity would be significantly stronger in the desynchronized state as compared to the synchronized state, one would expect the right panel in Figure 1E to more closely resemble the classical CSD pattern of a prominent granular sink followed by activation of other cortical layers [27, 32, 49]. However, the question whether lateral propagation of spontaneous waves of excitation is predominantly governed by intracortical mechanisms [15, 50–52] or whether there is an essential thalamic component involved [53–55] remains to be answered.
Lateral or horizontal spread of mass depolarization is thought of as propagation along a network of interconnected neurons or as a consequence of the phase delay between coupled local oscillators . In either case, the main factor determining the velocity of propagation is the strength of excitatory-excitatory connections in the network and the proximity of the membrane potential of excitatory neurons to threshold . The difference in lateral propagation velocity of spontaneous cortical waves between the two states can thus be interpreted as an indicator of stronger population synaptic strength, i.e. lateral coupling, in cortical networks under desynchronized conditions. Such changes in lateral coupling may result from the fact that during desynchronized states, the membrane potential of neurons is more depolarized , so that incoming excitatory postsynaptic potentials bring neurons faster to threshold as compared to more synchronized slow-wave states, where neurons are more hyperpolarized . This could lead to faster spread of activation across the network. Reported velocities for lateral propagation in the visual cortex differ greatly , and our results demonstrate that differences in brain state can contribute considerably to this discrepancy. With respect to cortical information processing, a faster spread of depolarizing waves would facilitate communication between spatially segregated neural ensembles by depolarizing more neurons in a shorter time window, thus enhancing the gain for incoming stimuli in a larger patch of cortical tissue. This, in turn, might be a factor contributing to efficient integration of information in desynchronized forebrain states  which accompany attentive modes of processing like dreaming or wakefulness.
Of note, sensory-evoked cortical waves propagate about twice as fast as spontaneous waves in the same cortical state (see Additional file 1), and propagation parameters of evoked waves are not significantly affected by the global cortical state (Additional file 1). These pronounced differences in propagation velocity between spontaneous and sensory-evoked waves are readily explained by the qualitatively different nature of spontaneous and evoked waves in our study. Since our visual stimulus is not a point-source (which is impractical in rodent visual cortex), the spread of these sensory-evoked waves is more strongly governed by subcortical-cortical interactions . Mechanistically, in this case, a patterned feedforward input from the thalamus might excite neighboring cortical regions with a time delay, the result of which would be an apparent wave motion pattern, corresponding to a subcortical oscillator model where a phase delay in excitatory drive results in fictive wave motion . To test this proposition requires experiments in a sensory area which can be stimulated with a point-source, such as the barrel cortex .
Since wave propagation is thought to have functional relevance in spatial integration, the state-dependence of preferred propagation trajectories is of particular note. Whereas anisotropic propagation of spontaneous cortical waves has been observed in both humans [2, 5] and animals [12, 15, 59], a brain state-dependent change in the directions of spontaneous traveling waves has, to our knowledge, not been demonstrated to date. Given that preferred trajectories associated with each state were fairly stable across animals, our data suggests a strong systematic modulation of effective network pathways by cortical state. A reversal in the bias of the direction of information transfer between higher- and lower-order visual areas  could potentially associate with different modes of "top-down" or "bottom-up" processing, and an important future direction will be to test the behavioral effects of selective inhibition of directional propagation.
While lateral propagation may depend on local cortical circuits , our results demonstrate that the global activity state of the cortical network also has a significant impact on the velocity and preferred direction of propagating waves. Our findings shed light on how propagation parameters, and thus spatiotemporal integration, changes as a function of cortical state.
Averaged rectified CSD
Current source density
Local field potential
Root mean square
Primary visual cortex
Secondary visual cortex
Secondary visual cortex, medial area
Slovin H, Arieli A, Hildesheim R, Grinvald A: Long-term voltage-sensitive dye imaging reveals cortical dynamics in behaving monkeys. J Neurophysiol. 2002, 88: 3421-3438. 10.1152/jn.00194.2002.
Massimini M, Huber R, Ferrarelli F, Hill S, Tononi G: The sleep slow oscillation as a traveling wave. J Neurosci. 2004, 24: 6862-6870. 10.1523/JNEUROSCI.1318-04.2004.
Lippert MT, Takagaki K, Xu W, Huang X, Wu J-Y: Methods for voltage-sensitive dye imaging of rat cortical activity with high signal-to-noise ratio. J Neurophysiol. 2007, 98: 502-512. 10.1152/jn.01169.2006.
Xu W, Huang X, Takagaki K, Wu J: Compression and reflection of visually evoked cortical waves. Neuron. 2007, 55: 119-129. 10.1016/j.neuron.2007.06.016.
Takahashi K, Saleh M, Penn RD, Hatsopoulos NG: Propagating waves in human motor cortex. Front Hum Neurosci. 2011, 5: 40.
Reimer A, Hubka P, Engel AK, Kral A: Fast propagating waves within the rodent auditory cortex. Cereb Cortex. 2011, 21: 166-177. 10.1093/cercor/bhq073.
Nauhaus I, Busse L, Ringach DL, Carandini M: Robustness of traveling waves in ongoing activity of visual cortex. J Neurosci. 2012, 32: 3088-3094. 10.1523/JNEUROSCI.5827-11.2012.
Freeman W, Quiroga RQ: Imaging Brain Function With EEG: Advanced Temporal and Spatial Analysis of Electroencephalographic Signals. 2012, New York: Springer, 2013
Ermentrout GB, Kleinfeld D: Traveling electrical waves in cortex: insights from phase dynamics and speculation on a computational role. Neuron. 2001, 29: 33-44. 10.1016/S0896-6273(01)00178-7.
Wu J-Y, Huang X, Zhang C: Propagating waves of activity in the neocortex: what they are, what they do. Neuroscientist. 2008, 14: 487-502.
Roland PE, Hanazawa A, Undeman C, Eriksson D, Tompa T, Nakamura H, Valentiniene S, Ahmed B: Cortical feedback depolarization waves: a mechanism of top-down influence on early visual areas. Proc Natl Acad Sci U S A. 2006, 103: 12586-12591. 10.1073/pnas.0604925103.
Rubino D, Robbins KA, Hatsopoulos NG: Propagating waves mediate information transfer in the motor cortex. Nat Neurosci. 2006, 9: 1549-1557. 10.1038/nn1802.
Benucci A, Frazor RA, Carandini M: Standing waves and traveling waves distinguish two circuits in visual cortex. Neuron. 2007, 55: 103-117. 10.1016/j.neuron.2007.06.017.
Han F, Caporale N, Dan Y: Reverberation of recent visual experience in spontaneous cortical waves. Neuron. 2008, 60: 321-327. 10.1016/j.neuron.2008.08.026.
Takagaki K, Zhang C, Wu J-Y, Lippert MT: Crossmodal propagation of sensory-evoked and spontaneous activity in the rat neocortex. Neurosci Lett. 2008, 431: 191-196. 10.1016/j.neulet.2007.11.069.
Sato TK, Nauhaus I, Carandini M: Traveling waves in visual cortex. Neuron. 2012, 75: 218-229. 10.1016/j.neuron.2012.06.029.
Muller L, Destexhe A: Propagating waves in thalamus, cortex and the thalamocortical system: Experiments and models. J Physiol Paris. 2012, 106: 222-238. 10.1016/j.jphysparis.2012.06.005.
Steriade MM, McCarley RW: Brain Control of Wakefulness and Sleep. 2005, New York: Springer, 2
Harris KD, Thiele A: Cortical state and attention. Nat Rev Neurosci. 2011, 12: 509-523. 10.1038/nrn3084.
Poulet JFA, Petersen CCH: Internal brain state regulates membrane potential synchrony in barrel cortex of behaving mice. Nature. 2008, 454: 881-885. 10.1038/nature07150.
Castro-Alamancos MA: Absence of rapid sensory adaptation in neocortex during information processing states. Neuron. 2004, 41: 455-464. 10.1016/S0896-6273(03)00853-5.
Curto C, Sakata S, Marguet S, Itskov V, Harris KD: A simple model of cortical dynamics explains variability and state dependence of sensory responses in urethane-anesthetized auditory cortex. J Neurosci. 2009, 29: 10600-10612. 10.1523/JNEUROSCI.2053-09.2009.
Lincoln DW: Correlation of unit activity in the hypothalamus with EEG patterns associated with the sleep cycle. Exp Neurol. 1969, 24: 1-18. 10.1016/0014-4886(69)90002-8.
Clement EA, Richard A, Thwaites M, Ailon J, Peters S, Dickson CT: Cyclic and sleep-like spontaneous alternations of brain state under urethane anaesthesia. PLoS One. 2008, 3: e2004-10.1371/journal.pone.0002004.
Takagaki K, Zhang C, Wu J-Y, Ohl FW: Flow detection of propagating waves with temporospatial correlation of activity. J Neurosci Methods. 2011, 200: 207-218. 10.1016/j.jneumeth.2011.05.023.
Lippert MT, Takagaki K, Kayser C, Ohl FW: Asymmetric Multisensory Interactions of Visual and Somatosensory Responses in a Region of the Rat Parietal Cortex. PLoS One. 2013, 8: e63631-10.1371/journal.pone.0063631.
Happel MFK, Jeschke M, Ohl FW: Spectral integration in primary auditory cortex attributable to temporally precise convergence of thalamocortical and intracortical input. J Neurosci. 2010, 30: 11114-11127. 10.1523/JNEUROSCI.0689-10.2010.
Steriade M, Timofeev I, Grenier F: Natural waking and sleep states: a view from inside neocortical neurons. J Neurophysiol. 2001, 85: 1969-1985.
Takagaki K, Lippert MT, Dann B, Wanger T, Ohl FW: Normalization of voltage-sensitive dye signal with functional activity measures. PLoS One. 2008, 3: e4041-10.1371/journal.pone.0004041.
Berens P: CircStat: A MATLAB Toolbox for Circular Statistics. J Stat Softw. 2009, 31: 1-21.
Castro-Alamancos MA, Oldford E: Cortical sensory suppression during arousal is due to the activity-dependent depression of thalamocortical synapses. J Physiol (Lond). 2002, 541: 319-331. 10.1113/jphysiol.2002.016857.
Mitzdorf U: Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. Physiol Rev. 1985, 65: 37-100.
Destexhe A, Contreras D, Steriade M: Spatiotemporal analysis of local field potentials and unit discharges in cat cerebral cortex during natural wake and sleep states. J Neurosci. 1999, 19: 4595-4608.
Sakata S, Harris KD: Laminar-dependent effects of cortical state on auditory cortical spontaneous activity. Front Neural Circuits. 2012, 6: 109.
Niell CM, Stryker MP: Modulation of visual responses by behavioral state in mouse visual cortex. Neuron. 2010, 65: 472-479. 10.1016/j.neuron.2010.01.033.
Gentet LJ, Avermann M, Matyas F, Staiger JF, Petersen CCH: Membrane potential dynamics of GABAergic neurons in the barrel cortex of behaving mice. Neuron. 2010, 65: 422-435. 10.1016/j.neuron.2010.01.006.
Petersen CCH, Hahn TTG, Mehta M, Grinvald A, Sakmann B: Interaction of sensory responses with spontaneous depolarization in layer 2/3 barrel cortex. Proc Natl Acad Sci U S A. 2003, 100: 13638-13643. 10.1073/pnas.2235811100.
Ferezou I, Bolea S, Petersen CCH: Visualizing the cortical representation of whisker touch: voltage-sensitive dye imaging in freely moving mice. Neuron. 2006, 50: 617-629. 10.1016/j.neuron.2006.03.043.
Luczak A, Barthó P, Marguet SL, Buzsáki G, Harris KD: Sequential structure of neocortical spontaneous activity in vivo. Proc Natl Acad Sci U S A. 2007, 104: 347-352. 10.1073/pnas.0605643104.
Mohajerani MH, McVea DA, Fingas M, Murphy TH: Mirrored bilateral slow-wave cortical activity within local circuits revealed by fast bihemispheric voltage-sensitive dye imaging in anesthetized and awake mice. J Neurosci. 2010, 30: 3745-3751. 10.1523/JNEUROSCI.6437-09.2010.
Csercsa R, Dombovári B, Fabó D, Wittner L, Eross L, Entz L, Sólyom A, Rásonyi G, Szucs A, Kelemen A, Jakus R, Juhos V, Grand L, Magony A, Halász P, Freund TF, Maglóczky Z, Cash SS, Papp L, Karmos G, Halgren E, Ulbert I: Laminar analysis of slow wave activity in humans. Brain. 2010, 133: 2814-2829. 10.1093/brain/awq169.
Markram H: A network of tufted layer 5 pyramidal neurons. Cereb Cortex. 1997, 7: 523-533. 10.1093/cercor/7.6.523.
Bannister AP: Inter- and intra-laminar connections of pyramidal cells in the neocortex. Neurosci Res. 2005, 53: 95-103. 10.1016/j.neures.2005.06.019.
Tanifuji M, Sugiyama T, Murase K: Horizontal propagation of excitation in rat visual cortical slices revealed by optical imaging. Science. 1994, 266: 1057-1059. 10.1126/science.7973662.
Wester JC, Contreras D: Columnar interactions determine horizontal propagation of recurrent network activity in neocortex. J Neurosci. 2012, 32: 5454-5471. 10.1523/JNEUROSCI.5006-11.2012.
Sanchez-Vives MV, McCormick DA: Cellular and network mechanisms of rhythmic recurrent activity in neocortex. Nat Neurosci. 2000, 3: 1027-1034. 10.1038/79848.
Beltramo R, D’Urso G, Dal Maschio M, Farisello P, Bovetti S, Clovis Y, Lassi G, Tucci V, De Pietri TD, Fellin T: Layer-specific excitatory circuits differentially control recurrent network dynamics in the neocortex. Nat Neurosci. 2013, 16: 227-234. 10.1038/nn.3306.
Chauvette S, Volgushev M, Timofeev I: Origin of active states in local neocortical networks during slow sleep oscillation. Cereb Cortex. 2010, 20: 2660-2674. 10.1093/cercor/bhq009.
Sakata S, Harris KD: Laminar structure of spontaneous and sensory-evoked population activity in auditory cortex. Neuron. 2009, 64: 404-418. 10.1016/j.neuron.2009.09.020.
Stroh A, Adelsberger H, Groh A, Rühlmann C, Fischer S, Schierloh A, Deisseroth K, Konnerth A: Making waves: initiation and propagation of corticothalamic Ca2+ waves in vivo. Neuron. 2013, 77: 1136-1150. 10.1016/j.neuron.2013.01.031.
Wright JJ, Liley DT: Simulation of electrocortical waves. Biol Cybern. 1995, 72: 347-356. 10.1007/BF00202790.
Timofeev I, Steriade M: Low-frequency rhythms in the thalamus of intact-cortex and decorticated cats. J Neurophysiol. 1996, 76: 4152-4168.
Farley BJ, Noreña AJ: Spatiotemporal coordination of slow-wave ongoing activity across auditory cortical areas. J Neurosci. 2013, 33: 3299-3310. 10.1523/JNEUROSCI.5079-12.2013.
Crunelli V, Hughes SW: The slow (<1 Hz) rhythm of non-REM sleep: a dialogue between three cardinal oscillators. Nat Neurosci. 2010, 13: 9-17. 10.1038/nn.2445.
Ushimaru M, Ueta Y, Kawaguchi Y: Differentiated participation of thalamocortical subnetworks in slow/spindle waves and desynchronization. J Neurosci. 2012, 32: 1730-1746. 10.1523/JNEUROSCI.4883-11.2012.
Goulet J, Ermentrout GB: The mechanisms for compression and reflection of cortical waves. Biol Cybern. 2011, 105: 253-268. 10.1007/s00422-011-0465-3.
Alkire MT, Hudetz AG, Tononi G: Consciousness and anesthesia. Science. 2008, 322: 876-880. 10.1126/science.1149213.
Devonshire IM, Grandy TH, Dommett EJ, Greenfield SA: Effects of urethane anaesthesia on sensory processing in the rat barrel cortex revealed by combined optical imaging and electrophysiology. Eur J Neurosci. 2010, 32: 786-797. 10.1111/j.1460-9568.2010.07322.x.
Fucke T, Suchanek D, Nawrot MP, Seamari Y, Heck DH, Aertsen A, Boucsein C: Stereotypical spatiotemporal activity patterns during slow-wave activity in the neocortex. J Neurophysiol. 2011, 106: 3035-3044. 10.1152/jn.00811.2010.
We thank Lydia Löw and Kathrin Ohl for expert technical assistance. We also thank Prof. Jian-Young Wu and Prof. Christoph Kayser for scientific discussion. This study was supported by Deutsche Forschungsgemeinschaft DFG SFB-779, DFG SFB-TRR 31 (FWO), LIN-founded Special Project 2007/08 (TW, JG and FWO), CBBS-founded NeuroNetworks Project USTR-3804R and Alexander von Humboldt Stiftung (KT).
The authors declare that they have no competing interests.
TW carried out the acquisition of VSD and ECoG data, the analysis of VSD and ECoG data and wrote the manuscript. KT carried out the acquisition of VSD and ECoG data, and participated in the analysis of VSD and ECoG data, in the preparation of the manuscript and in the design of the study. MTL carried out the acquisition of CSD data, the analysis of CSD data, and participated in the acquisition of VSD and ECoG data, and in the preparation of the manuscript. JG participated in the design of the study. FWO participated in the analysis of CSD data, in the preparation of the manuscript and in the design of the study. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
- Traveling wave
- Wave propagation
- Voltage-sensitive dye imaging
- Current source density
- Rat visual cortex