Wave propagation of cortical population activity under urethane anesthesia is state dependent
© Wanger et al.; licensee BioMed Central Ltd. 2013
Received: 30 March 2013
Accepted: 3 July 2013
Published: 31 July 2013
Skip to main content
© Wanger et al.; licensee BioMed Central Ltd. 2013
Received: 30 March 2013
Accepted: 3 July 2013
Published: 31 July 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.
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 .
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).
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 .
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.
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.
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).
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
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).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.