BACKGROUND
spontaneous and sensory-evoked population activity in the mammalian neocortex is known to manifest in the form of traveling or propagating waves .
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  <cit> .
despite this proposed role of propagating waves in spatially coordinating neuronal computation  <cit> , the mechanisms governing wave propagation remain an issue of ongoing debate  <cit> .
in particular, little is known about how propagation dynamics are influenced by global network states  <cit> .
distinct global network states are a critical feature of neocortical dynamics, and characterize sleep and wakefulness  <cit> , as well as stages of arousal and anesthesia  <cit> .
these global network states can have profound effects on cortical information processing by biasing the dynamics of membrane potential synchrony  <cit>  and the interaction of sensory evoked responses with ongoing activity  <cit> .
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  <cit> .
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  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  spontaneously alternates between these two states , which have been shown to resemble, in some respects, forebrain states observed during sleep  <cit> .
we first used laminar current source density  analysis to investigate the cortical laminar organization of the spontaneous population activity patterns in each state.
using a temporospatial flow detection algorithm  <cit> , 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.
methods
all experiments were in compliance with the guidelines of the european community  and were approved by an ethics commission of the state of sachsen-anhalt.
twelve adult male wistar rats  were used for the experiments described here.
anesthesia
induction of anesthesia was performed by either intraperitoneal or intravenous infusion of an aqueous solution of urethane .
further anesthetic was given as necessary to maintain areflexia.
csd analysis
to analyze the layered activity distribution in both cortical states, three animals were implanted with a silicon microelectrode  into the primary visual cortex .
the back surface of the electrode was coated with the fluorescent carbocyanine dye diic <dig> to later localize the recording site in histological slices  <cit> .
laminar positioning of the electrode was controlled for by inserting it  <dig>  mm below the brain surface.
cortical layers were delineated based both on histology and on the typical laminar csd patterns following visual stimulation  <cit> .
electric field potentials were recorded with a biosignal recording system  and split into local field potentials  and multi-unit activity .
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  <cit> .
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  <cit> , the window was centered on each time point where the integrated mua  crossed a threshold of three standard deviations above baseline  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  <cit> .
we therefore filtered the lfp in a band of 3–6 hz, which typically shows a spectral peak in the desynchronized state .
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  and cortex was stained transdurally with the dye rh- <dig>  for  <dig>  h. optical signals were recorded and processed as described previously  <cit> .
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 .
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  <dig> –100 hz, prior to digitization.
the electrocardiogram was recorded as well, to allow triggered removal of heartbeat-related artifacts from the optical recording  <cit> .
classification of cortical state
in order to perform offline classification of cortical states, root mean square  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 .
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  and mathematica .
on average, episodes of each cortical state lasted between 2–10 min .
however, some episodes of the synchronized state lasted up to 50 min .
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 .
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  <cit> .
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  <cit> , 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 .
the performance characteristics and implementation details of this algorithm are discussed in detail elsewhere  <cit> .
in the present study, a correlation window of 200 ms  in combination with a maximum shift window of 25 ms  was used to calculate temporospatial flow; correlation windows of  <dig> and 250 ms showed similar results .
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 v <dig> as indicated in figure 2a.
the correlation window was moved in successive steps of  <dig> frames over the bandpass-filtered  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  <dig>  were discarded  <cit> .
flow trajectory distributions were tested for statistically significant differences using a circular kuiper two-sample test from the circstat toolbox  <cit> .
RESULTS
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  <cit> .
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  <cit> .
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  <cit> .
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 .
the variance of the total avrec which is explained by supragranular sink-source dipoles was  <dig> %  in the synchronized state and  <dig> %  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 , presumably reflecting the decreased spatial and temporal correlation of membrane potentials in more active forebrain states  <cit> .
local propagation speeds of spontaneous waves varied greatly across events and trials, apparently following poisson statistics , 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 , wilcoxon matched-pairs signed-ranks test, p <  <dig> ).
trajectories of spontaneous cortical waves are state-dependent
as for the directional properties of propagation, spontaneous waves measured at the center of v <dig> traveled in every direction , but with a highly anisotropic distribution in both states .
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 , differing by 180° for synchronized and desynchronized states.
notably, the respective directions of anisotropy were very stable across animals , 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.
discussion
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  <cit> , cell-type-specific firing patterns  <cit>  and membrane potential dynamics  <cit> .
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  <cit> .
in contrast, wave propagation in the urethane-desynchronized state most likely reflects the spread of more localized patterns of oscillatory synchrony  with a spectral peak around 4 hz .
despite these apparent spatiotemporal differences of the propagating waves, laminar csd profiles of synchronized and desynchronized states appear strikingly similar .
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  presumably reflect postsynaptic input to somata and dendrites of layer 2/ <dig> pyramidal neurons as well as to the distal dendrites of large layer  <dig> pyramidal neurons which ramify extensively into superficial layers  <cit> .
after being depolarized to firing threshold, pyramidal neurons in layers 2/ <dig> and  <dig> would then pass on the excitation to their efferent target neurons via intra- and interlaminar projections  <cit> , resulting in lateral spread of the depolarization as observed in supragranular and infragranular layers of cortical slices  <cit> .
of note, as there is accumulating evidence suggesting a leading role for layer  <dig> networks in cortical up state initiation  <cit>  and intracortical wave propagation  <cit> , 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/ <dig>  this scenario is supported by a recent study showing that layer 5a provides strong driving input to supragranular layers  <cit> .
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  <cit> .
however, the question whether lateral propagation of spontaneous waves of excitation is predominantly governed by intracortical mechanisms  <cit>  or whether there is an essential thalamic component involved  <cit>  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  <cit> .
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  <cit> .
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  <cit> , so that incoming excitatory postsynaptic potentials bring neurons faster to threshold as compared to more synchronized slow-wave states, where neurons are more hyperpolarized  <cit> .
this could lead to faster spread of activation across the network.
reported velocities for lateral propagation in the visual cortex differ greatly  <cit> , 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  <cit>  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 , and propagation parameters of evoked waves are not significantly affected by the global cortical state .
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 , the spread of these sensory-evoked waves is more strongly governed by subcortical-cortical interactions  <cit> .
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  <cit> .
to test this proposition requires experiments in a sensory area which can be stimulated with a point-source, such as the barrel cortex  <cit> .
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  <cit>  and animals  <cit> , 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  <cit>  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.
CONCLUSIONS
while lateral propagation may depend on local cortical circuits  <cit> , 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.
abbreviations
avrec: averaged rectified csd; csd: current source density; ecog: electrocorticogram; lfp: local field potential; mua: multi-unit activity; rms: root mean square; v1: primary visual cortex; v2: secondary visual cortex; v2m: secondary visual cortex, medial area; vsd: voltage-sensitive dye.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
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.
supplementary material
additional file 1
propagation of sensory-evoked activity does not depend on cortical state.
click here for file
acknowledgements
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- <dig>  dfg sfb-trr  <dig> , lin-founded special project 2007/ <dig> , cbbs-founded neuronetworks project ustr-3804r and alexander von humboldt stiftung .

 spontaneous and sensory-evoked population activity in the mammalian neocortex is known to manifest in the form of traveling or propagating waves. 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. despite this proposed role of propagating waves in spatially coordinating neuronal computation, 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. 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 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 spontaneously alternates between these two states, which have been shown to resemble, in some respects, forebrain states observed during sleep. we first used laminar current source density 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. methods all experiments were in compliance with the guidelines of the european community and were approved by an ethics commission of the state of sachsen-anhalt. twelve adult male wistar rats were used for the experiments described here. anesthesia induction of anesthesia was performed by either intraperitoneal or intravenous infusion of an aqueous solution of urethane. further anesthetic was given as necessary to maintain areflexia. csd analysis to analyze the layered activity distribution in both cortical states, three animals were implanted with a silicon microelectrode into the primary visual cortex. the back surface of the electrode was coated with the fluorescent carbocyanine dye diic to later localize the recording site in histological slices. laminar positioning of the electrode was controlled for by inserting it  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 and split into local field potentials and multi-unit activity. 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. 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 crossed a threshold of three standard deviations above baseline 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. 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 and cortex was stained transdurally with the dye rh- for  h. optical signals were recorded and processed as described previously. 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. 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 –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 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. 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 and mathematica. on average, episodes of each cortical state lasted between 2–10 min. however, some episodes of the synchronized state lasted up to 50 min. 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. 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. the performance characteristics and implementation details of this algorithm are discussed in detail elsewhere. in the present study, a correlation window of 200 ms in combination with a maximum shift window of 25 ms was used to calculate temporospatial flow; correlation windows of and 250 ms showed similar results. 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 v as indicated in figure 2a. the correlation window was moved in successive steps of frames over the bandpass-filtered 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 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. 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. the variance of the total avrec which is explained by supragranular sink-source dipoles was % in the synchronized state and % 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, presumably reflecting the decreased spatial and temporal correlation of membrane potentials in more active forebrain states. local propagation speeds of spontaneous waves varied greatly across events and trials, apparently following poisson statistics, 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, wilcoxon matched-pairs signed-ranks test, p <  ). trajectories of spontaneous cortical waves are state-dependent as for the directional properties of propagation, spontaneous waves measured at the center of v traveled in every direction, but with a highly anisotropic distribution in both states. 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, differing by 180° for synchronized and desynchronized states. notably, the respective directions of anisotropy were very stable across animals, 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. discussion 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 and membrane potential dynamics. 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. in contrast, wave propagation in the urethane-desynchronized state most likely reflects the spread of more localized patterns of oscillatory synchrony with a spectral peak around 4 hz. despite these apparent spatiotemporal differences of the propagating waves, laminar csd profiles of synchronized and desynchronized states appear strikingly similar. 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 presumably reflect postsynaptic input to somata and dendrites of layer 2/ pyramidal neurons as well as to the distal dendrites of large layer pyramidal neurons which ramify extensively into superficial layers. after being depolarized to firing threshold, pyramidal neurons in layers 2/ and 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. of note, as there is accumulating evidence suggesting a leading role for layer networks in cortical up state initiation and intracortical wave propagation, 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/ 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. however, the question whether lateral propagation of spontaneous waves of excitation is predominantly governed by intracortical mechanisms or whether there is an essential thalamic component involved 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, and propagation parameters of evoked waves are not significantly affected by the global cortical state. 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, 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 and animals, 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. abbreviations avrec: averaged rectified csd; csd: current source density; ecog: electrocorticogram; lfp: local field potential; mua: multi-unit activity; rms: root mean square; v1: primary visual cortex; v2: secondary visual cortex; v2m: secondary visual cortex, medial area; vsd: voltage-sensitive dye. competing interests the authors declare that they have no competing interests. authors’ contributions 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. supplementary material additional file 1 propagation of sensory-evoked activity does not depend on cortical state. click here for file acknowledgements 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- dfg sfb-trr, lin-founded special project 2007/, cbbs-founded neuronetworks project ustr-3804r and alexander von humboldt stiftung.
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. traveling wavewave propagationcortexstate-dependencevoltage-sensitive dye imagingurethanecurrent source densityrat visual cortex 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.
 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. traveling wavewave propagationcortexstate-dependencevoltage-sensitive dye imagingurethanecurrent source densityrat visual cortex 