PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcneusBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Neuroscience
 
BMC Neurosci. 2010; 11: 122.
Published online Sep 23, 2010. doi:  10.1186/1471-2202-11-122
PMCID: PMC2955690
Cross-approximate entropy of cortical local field potentials quantifies effects of anesthesia - a pilot study in rats
Matthias Kreuzer,#1 Harald Hentschke,corresponding author#2,3 Bernd Antkowiak,2 Cornelius Schwarz,3,4 Eberhard F Kochs,1 and Gerhard Schneider1
1Department of Anesthesiology, Klinikum rechts der Isar, Technische Universität München, München, Germany
2Section of Experimental Anaesthesiology, University Hospital Tübingen, Tübingen, Germany
3Department of Cognitive Neurology, Hertie-Institute for Clinical Brain Research, University Tübingen, Tübingen, Germany
4Werner Reichardt Centre for Integrative Neuroscience, University Tübingen, Tübingen, Germany
corresponding authorCorresponding author.
#Contributed equally.
Matthias Kreuzer: m.kreuzer/at/lrz.tu-muenchen.de; Harald Hentschke: harald.hentschke/at/uni-tuebingen.de; Bernd Antkowiak: bernd.antkowiak/at/uni-tuebingen.de; Cornelius Schwarz: cornelius.schwarz/at/uni-tuebingen.de; Eberhard F Kochs: e.f.kochs/at/lrz.tu-muenchen.de; Gerhard Schneider: g.schneider/at/lrz.tu-muenchen.de
Received February 26, 2010; Accepted September 23, 2010.
Background
Anesthetics dose-dependently shift electroencephalographic (EEG) activity towards high-amplitude, slow rhythms, indicative of a synchronization of neuronal activity in thalamocortical networks. Additionally, they uncouple brain areas in higher (gamma) frequency ranges possibly underlying conscious perception. It is currently thought that both effects may impair brain function by impeding proper information exchange between cortical areas. But what happens at the local network level? Local networks with strong excitatory interconnections may be more resilient towards global changes in brain rhythms, but depend heavily on locally projecting, inhibitory interneurons. As anesthetics bias cortical networks towards inhibition, we hypothesized that they may cause excessive synchrony and compromise information processing already on a small spatial scale. Using a recently introduced measure of signal independence, cross-approximate entropy (XApEn), we investigated to what degree anesthetics synchronized local cortical network activity. We recorded local field potentials (LFP) from the somatosensory cortex of three rats chronically implanted with multielectrode arrays and compared activity patterns under control (awake state) with those at increasing concentrations of isoflurane, enflurane and halothane.
Results
Cortical LFP signals were more synchronous, as expressed by XApEn, in the presence of anesthetics. Specifically, XApEn was a monotonously declining function of anesthetic concentration. Isoflurane and enflurane were indistinguishable; at a concentration of 1 MAC (the minimum alveolar concentration required to suppress movement in response to noxious stimuli in 50% of subjects) both volatile agents reduced XApEn by about 70%, whereas halothane was less potent (50% reduction).
Conclusions
The results suggest that anesthetics strongly diminish the independence of operation of local cortical neuronal populations, and that the quantification of these effects in terms of XApEn has a similar discriminatory power as changes of spontaneous action potential rates. Thus, XApEn of field potentials recorded from local cortical networks provides valuable information on the anesthetic state of the brain.
It has long been known that anesthetics alter neuronal activity of the brain, yet it is still unclear which of the alterations of brain activity are instrumental in causing sedation, amnesia and unconsciousness. Anesthetics likely disrupt communication between cortical areas [1], and thus impair large-scale integration of information hypothesized to be a prerequisite to proper brain function, particularly conscious perception [2-5]. A number of experimental observations are in agreement with this concept. In rats stimulated with light flashes, volatile anesthetics disrupted anterior-posterior phase synchronization of field responses [6] and depressed long-latency spike responses in visual cortex, thought to arise from cortico-cortical interactions [7]. In humans, during the transition from waking to loss of consciousness, various general anesthetics decoupled gamma rhythms between anterior and posterior cortical areas as well as between homologous areas in different hemispheres [8]. Employing transcranial magnetic stimulation, Ferrarelli et al. recently demonstrated a breakdown of cortical effective connectivity upon loss of consciousness induced by midazolam [9].
A complementary insight into network effects of anesthetics was recently provided by Hudetz and colleagues, who demonstrated that in rats volatile anesthetics diminished the independence of spontaneous electroencephalographic (EEG) signals recorded from different hemispheres [10]: isoflurane, and halothane above 0.4%, reduced cross-approximate entropy (XApEn), a nonlinear information statistical parameter which quantifies the independence (or dissimilarity) of signals. Fittingly, cholinergic stimulation reversed the isoflurane-induced decrease of XApEn [11]. These findings suggest that anesthetics induce changes of thalamocortical networks and slow predominant frequencies. This impairs brain function by promoting uniformity of signals and impeding an exchange of independent information between cortical areas, as a contrast to the transient synchronization of fast oscillations which has been suggested as the mechanism underlying conscious perception [2].
It is less clear at present in which manner anesthetics impair information processing on a much smaller spatial scale. In local cortical networks, subnetworks defined by strong excitatory connections exist which may operate quite independently of each other [12,13]. With anesthetics, their independence of operation may be compromised by a general decrease of cortical neuronal excitability [14-18]. Furthermore, inhibitory interneurons project locally and have a great potential to pace their postsynaptic targets [19,20]. Under conditions of pharmacologically enhanced GABAergic transmission, they may coerce independent subnetworks into more synchronous, uniform activity patterns [21,22].
In the present pilot study, we investigated to what degree anesthetics alter signal independence in the somatosensory ('barrel') cortex of the rat. We recorded spontaneous local field potential (LFP) activity from multiple, closely spaced electrodes, which sample from a much smaller subset of neurons than EEG electrodes [23,24]. The effects of the volatile anesthetics isoflurane, enflurane and halothane were evaluated in three animals. Our results show that the independence of activity patterns across recording sites as quantified by XApEn is a monotonously declining function of anesthetic concentration, suggesting that volatile anesthetics strongly promote uniform activity already on the level of local neocortical circuits.
Figure Figure11 illustrates the location of the electrode arrays within neocortex and shows exemplary recordings from an experiment in which isoflurane was administered. Under control conditions, spectral power of the LFP peaked at 6-7 Hz and gradually declined towards higher frequencies (Figure (Figure1C,1C, top). The signals recorded from different electrodes seemed similar to varying degrees. In particular, similarity between signals depended on the spatial arrangement of the electrodes: within each row, adjacent electrodes recorded signals which were in general more similar to each other than electrodes placed further apart. This is evident in the pattern of cross correlation values in Figure Figure1D1D (top). Moreover, the signals from electrode row 1 had a higher amplitude than those from row 2 and qualitatively appeared somewhat different. The qualitative difference in signal appearance between rows was reflected in cross correlation values which were on average lower than within-row correlations (Figure (Figure1D,1D, top; compare upper right rectangular area with the two triangular areas, separated by white broken lines). We attribute the differences in signal amplitude and appearance between electrode rows to the different layers in which they were located (see Methods; [25,26]) as well as to their separation by more than twice the interelectrode distance (500 versus 180 um).
Figure 1
Figure 1
Signal characteristics of neocortical local field potentials. A, top view (schematic) of a rat cortical hemisphere and approximate location of the recording electrode array in one animal (rat1). The electrodes were arranged in two rows and placed in the (more ...)
At all concentrations applied, isoflurane greatly enhanced the amplitude of low-frequency signal components (up to ~20 Hz) and reduced high frequency components (above 100 Hz, Figure 1B, C). At the highest concentration (1.45%), high-amplitude LFP 'spikes' occurred, and in one animal burst suppression patterns appeared. Enflurane and halothane led to signal changes which were similar at low concentrations, but at concentrations close to 1 MAC the anesthetics clearly differed in their tendency to produce high-amplitude LFP spikes and burst suppression patterns (enflurane > isoflurane > halothane). Both the shifts in the spectral composition of neural activity common to all three anesthetics and the agent-specific differences are in agreement with electroencephalographic observations in rats [10,11,18,27-29].
XApEn quantifies the regularity of patterns in a pair of time series (Figure (Figure2).2). Before embarking on the time-consuming computation of the entire data set, we were interested in finding out whether it related, qualitatively, to cross correlation, a classical measure of the similarity of signals (see Figure Figure1D).1D). Figure Figure3A3A shows LFP excerpts recorded from three electrodes under control conditions. In contrast to all other data analyzed and presented in this study, this excerpt was chosen to include an episode of 'high voltage rhythmic spikes' (HVRS, see Methods). XApEn values of all three channel combinations (Figure (Figure3B)3B) were very similar throughout the length of the segment although their cross correlations (Figure (Figure3C)3C) differed substantially, in agreement with physical electrode distance. Moreover, XApEn reacted very sensitively to HVRS, declining during this episode. Cross correlation values also changed during HVRS, but the relative change was much smaller than that of XApEn. Thus, this analysis illustrates that XApEn captures characteristics of cortical LFP activity beyond similarity as quantified by cross correlation, and emphasizes that XApEn is well suited to detect changes in cortical network patterns.
Figure 2
Figure 2
Calculation of XapEn. Short excerpt of raw data from three channels shown to illustrate the computation of XApEn. Parameters were set to r = 23 μV (tolerance) and m = 1 (sequence length). Thus, the 'sequence' length is one data point. For illustration (more ...)
Figure 3
Figure 3
XApEn compared to cross correlation. A, raw data excerpts of 30 s length from channels 1,2 and 7 of the same recording as shown in Figure 1 (control condition). Note the sequence of 'high voltage rhythmic spikes' (HVRS) starting at ~12 s. B, XApEn values (more ...)
XApEn values were subsequently computed for all channel combinations and depicted in color-coded matrices. Figure Figure44 compares two different experimental sessions with one animal, separated by one day, with isoflurane and halothane administration, respectively. As seen under control conditions, the matrices formed fingerprint-like patterns which did not seem to depend on the spatial relations between electrodes. Although the baseline values of XApEn differed, the patterns were quite similar in the two control measurements shown, indicating that no major changes of the cortical network had occurred in the interim. In general, XApEn patterns were remarkably stable over the period of several days in which the experiments were conducted: cross correlations of the linearized XApEn values of the control conditions ranged from 0.86 to 0.99 (median values 0.94, 0.99 and 0.99 for the three animals, each of which was subjected to three recording sessions on separate days).
Figure 4
Figure 4
Anesthetic-induced changes of XapEn. Each color-coded matrix contains values obtained from non-redundant pairwise combinations of recording sites (cf. Figure 1D). The animal was exposed, in separate experimental sessions, to increasing concentrations (more ...)
Volatile anesthetics strongly and reversibly decreased XApEn. In the experiments shown in Figure Figure4,4, the lowest doses of isoflurane (0.3%) and halothane (0.65%) reduced XApEn values to approximately 50%. Concentration increases to 1 MAC produced only minor additional changes. This was confirmed in dose-response relationships shown in Figure Figure5.5. The data points, normalized to control values, were fitted with the monoexponential function y = exp(-x/tau). The free parameter tau, describing the decrease of XApEn with anesthetic concentration expressed in MAC, was as follows [95% confidence intervals]: 0.72 [0.58 0.86] for isoflurane, 0.72 [0.58 0.86] for enflurane and 1.70 [1.49 1.90] for halothane. Thus, halothane affected the cortical network less strongly than isoflurane and halothane. However, all three anesthetics were about equally potent at concentrations corresponding to human MACawake, which render 50% of a human test population unconscious [30-34] and which should be at the lower end of concentrations affecting cortical function of rats (blue rectangles in Figure Figure5).5). Both observations are consistent with the action potential-depressing properties of the agents [16]. XApEn computed from strongly lowpass-filtered versions of the data segments (cutoff frequency of 30 Hz as opposed to 200 Hz) yielded qualitatively comparable, but numerically inferior results (Additional file 1, Figure S1A), suggesting that signal changes in both low and high frequency ranges contribute to the anesthetic-induced decrease of XApEn.
Figure 5
Figure 5
Dose-response curves of XapEn. Dependence of cross-approximate entropy on anesthetic concentration expressed in MAC. The data points represent values from all non-redundant pairwise combinations of recording sites, normalized to control and averaged (± (more ...)
XApEn calculated from the LFP of closely spaced intracortical sites showed significant changes with anesthetic concentrations. This finding demonstrates that volatile anesthetics coerce small cortical sub-networks, here represented by rat barrel cortex, into uniform, synchronized activity patterns.
Barrel cortex forms a large part of rodent somatosensory cortex, characterized by a one-to-one correspondence between the sensory organs (follicles at the base of the large facial whiskers) and cytoarchitectonically segregated structures in layer 4 termed 'barrels' [35-37]. This columnar, somatotopic organization results in an orderly bottom-up spread of sensory-evoked activity which, in the initial stage of processing, is spatially restricted to the discrete termination zones of the major thalamic afferents [38-40]. Yet, barrel cortex is also characterized by a large degree of synaptic divergence and interconnectivity, characteristic of neocortex in general [41]. Axons of pyramidal neurons may span several barrels, especially in layers 2/3 [40,42], and other cortical areas including contralateral barrel cortex form reciprocal, spatially dispersed connections [37,43]. Probably owing to this high degree of cortical interconnectivity, the permanent, ongoing activity therein, unrelated to sensory input, has a spatiotemporal profile which is largely independent of the barrel architecture [44,45]. In the present study, we recorded and analyzed this activity. Signal characteristics showed a dramatic and fundamental change during administration of volatile anesthetics.
The decrease of XApEn with all three volatile anesthetics suggests that these agents transform the diversity of synaptic inputs impinging on closely spaced pyramidal cells into a more uniform, synchronous pattern. These patterns may arise from an enhancement of GABAergic currents and a weakening of glutamatergic currents [46], which bias synaptic communication towards inhibition. Specifically, GABAergic interneurons likely gain in influence on cortical activity patterns by entraining local networks to common rhythms. Both the relative insensitivity of some interneuron classes to GABAergic inhibition [47] and the finding that inhibitory inputs in neighboring pyramidal neurons are more synchronous than excitatory inputs [20] are consistent with this idea. Furthermore, input via long-range connections from other cortical areas is impaired under anesthesia [9,48,49]. Experimental findings on figure-ground separation in monkey visual cortex could also be interpreted along the lines of a functional disconnection of cortical areas with anesthetics [50]. Therefore, it seems likely that the long-range excitatory synaptic inputs emanating from various cortical areas have an overall desynchronizing influence on local network activity in barrel cortex, and, by extension, that the impairment of this input by volatile anesthetics contributes to more local synchrony and thus a decrease of XApEn as observed.
Given that intracortical connections outnumber subcortical afferents [41] and that volatile anesthetics alter network activity in isolated cortical networks in vitro [16,51] we argue that the decline of XApEn was to a substantial part due to intracortical effects of the anesthetics. In addition, decreases in XApEn probably also reflect the anesthetic-induced transformation of activity in subcortical areas projecting to cortex. In particular thalamus, with its intricate reciprocal connections to cortex [36,37], must be considered [52]. Shown to be sensitive to volatile anesthetics and prone to bursting behavior, it may imprint its activity patterns on cortex [53-55]. Other likely candidates include the basal forebrain, which modulates cortical activity via cholinergic, GABAergic and glutamatergic afferents [56-58] as well as hypothalamic sleep pathways [59].
We found that the three volatile anesthetics differed in their potency to alter cortical activity patterns as quantified by XApEn. While isoflurane and enflurane were indistinguishable, halothane had significantly weaker effects. This finding fits the profile of this anesthetic, which has previously been found to exert weaker effects than isoflurane on XApEn computed from interhemispheric EEG [10] and on spontaneous action potential activity in neocortex in vitro and in vivo [16,51]. In rat visual cortex, Imas et al. found an enhancement of event-related gamma oscillations at intermediate concentrations of halothane [60], a finding which underlines the particular characteristic of this anesthetic. A potential limitation of our results is the fact that body temperature of the animals was not controlled, possibly leading to hypothermia during anesthesia. To minimize the influence of this phenomenon, the animals were placed into a plastic housing which provided some thermal insulation, and therefore at least it seems unlikely that they experienced severe hypothermia.
It is surprising that isoflurane- and halothane-induced changes of XApEn computed from interhemispheric EEG and LFP match qualitatively so well (Hudetz et al. used a different set of filter frequencies and computational parameters, so that a quantitative comparison of values may not lead to valid results [10,11]). First, although the basal mechanisms underlying EEG and LFP signals are identical - synaptic population currents in pyramidal cells giving rise to extracellular potential gradients - the populations of cells sampled from are not. EEG electrodes record potentials from a much larger population than intracortical electrodes. Furthermore, EEG signals are dominated by the largest dipole-generating contributors, pyramidal cells of layer 5, which extend their dendrites to layer 1. Intracortical signals, by contrast, are sampled from a restricted spatial volume [24] and are thus lamina-dependent [25,26]. Second, cortical sites separated by less than 1-2 mm as in our experiments receive a great deal of common synaptic input [20,61], and consequently exhibit more synchronous activity than cortical networks in different hemispheres. Thus, the qualitatively similar depression of XApEn values reported by Hudetz et al. and here show that anesthetics unfold synchronizing effects on widely different spatial scales.
Anesthetics affect intracortical connections as well as corticocortical connections and inputs from subcortical areas. A disconnection or suppression of communication between different cortical areas is hypothesized to be a key player in the process of unconsciousness [1,62]. Buzsáki [4] suggested that local computations register in large parts of the cortex through long range-connections, and that this local-global wiring is necessary for subjective experiences. Our results suggest that local intracortical communication suffers from the uniformity of signals induced by volatile anesthetics, in parallel to a suppression or disconnection of long-range interareal connections as evident in EEG studies. Thus, our results underline the usefulness of multisite electrophysiological recordings - be it LFP in experimental animals or EEG in humans - in combination with XApEn and related parameters [29,63] for the quantification of anesthetic effects. By the same token, such approaches may serve to probe specific hypotheses on the neurophysiological correlate of consciousness, particularly those postulating a precisely timed convergence of synaptic inputs in neocortex [64], and possibly also to detect pathological connectivity between brain areas.
In vivo surgery & recording
All procedures described were in accordance with the policy on the use of animals in neuroscience research of the Society for Neuroscience and approved by the Ethical Committee on Animal Care and Use of the Government of Baden Württemberg, Germany. The experimental procedures were identical to those described in Hentschke et al [16] and the field potential data analyzed here were won from three of the four animals in the same study. Briefly, male or female Sprague-Dawley rats aged 12-16 weeks were anesthetized with ketamine/xylazine (100 mg/kg and 15 mg/kg, respectively). A craniectomy over the right hemisphere was performed, the dura was removed and a steel ring (O.D. 5 mm) placed on the pial surface. 4-16 custom-made etched tungsten microelectrodes (impedance 1-5 MΩ, tip separation ~180 μm) arranged as single or double linear arrays (one animal, 1 × 4; two animals, 2 × 8) were implanted through the steel ring into the neocortical somatosensory representation of the mystacial vibrissae ('barrel cortex'). Stereotactical target coordinates were Bregma -3.2 mm and lateral 4 mm. The array was oriented at an angle of 30-40° relative to the cortical surface such that the medial electrode row touched the cortical surface first and was thus located in a deeper cortical lamina than the lateral electrode row; the latter was aimed at the border of layers 4 and 5 (depth of penetration for this row from the point of contact with the pia was ~1100 um). Proper location of the electrodes within barrel cortex was verified by mapping of neuronal responses elicited by manual deflection of individual vibrissae. The electrode array and a head post were then fixed with dental cement and via 7-9 stainless steel screws driven into the skull. Recordings commenced after a recovery period of 2-4 d. The animals were head-restrained, sitting in a plastic housing, and placed in a sealable plexiglass box into which vaporizers (Drägerwerke, Lübeck, Germany) driven by air pumps delivered the anesthetics. For any given concentration, animals were exposed for periods of 27-39 min to the anesthetic. After discontinuation of the anesthetic, the animals were exposed to air for 40-60 min. A maximum of one recording session was performed per day.
Voltage traces recorded from the electrodes were referenced to the steel ring (rat1 and rat2) or to one of the electrodes (rat3), amplified and passed through a bandpass filter (-3 dB passband 0.5-200 Hz), digitized at 20 kHz and stored with a multichannel extracellular recording system (MultiChannelSystems, Reutlingen, Germany). Deteriorated, low-impedance electrodes with a small signal amplitude were excluded from analysis. The number of useable electrodes were 3 (rat3), 6 (rat2) and 12 (rat1). Anesthetic concentration is expressed in volume-% or as MAC (minimum alveolar concentration) as given in [65]: isoflurane 1.46%, enflurane 2.21%, halothane 1.03%.
Selection of data segments
Under control conditions, the field potential data contained episodes of oscillatory activity with a dominant peak at 8-9 Hz (theta component) and an additional (non-harmonic) peak at 13-16 Hz (beta component; an example is given in Figure Figure3).3). The nature of these oscillations, also called 'high voltage rhythmic spikes' (HVRS), is a matter of debate. They appear in resting animals, often in conjunction with low-amplitude whisker movements [66] (whisker 'twitching') and may reflect a specific kind of idling, responsive state of the whisker sensory system [67-69]. Others consider HVRS manifestations of absence epilepsies [70-72]. We observed that at subanesthetic to hypnotic concentrations (isoflurane, 0.3-0.75%; enflurane, 0.5-1.2%; halothane, 0.65-1.05%) characteristics of the oscillations changed in several ways. Most notably, they were less coherent across channels and lasted for shorter periods. Furthermore, the dominant frequency was between 13 and 16 Hz. At the highest concentrations, which were equal to or above MAC (isoflurane, 1.45%; enflurane, 2.4%; halothane, 1.6%) the oscillations subsided and low frequency components and/or burst suppression patterns dominated.
We decided to exclude data with HVRS from analysis due to the unresolved nature of this activity and because the computational load of our analysis restricted the amount of data that could be analyzed in reasonable time to about 10 seconds per recording. To this end, the data acquired under control conditions were divided into segments of 2048 points overlapping by 730 points. From these, the spectral power of the signals in the range 7-9 Hz was computed for all channels. Of the resulting segment-wise power values, the 75th and 90th percentile were determined for each channel. We rejected segments with a power greater than the 90th percentile on any of the recorded channels and/or with a power larger than the 75th percentile on 50% or more of all recorded channels. The procedure was repeated for recordings with anesthetics, but with power determined in the range 10-16 Hz. Finally, within each recording, groups of four consecutive (overlapping) intervals which satisfied the criteria above were combined to yield segments of 6002 points length, corresponding to ~3 s. Three of such segments per recording, picked randomly from the beginning, middle and end of each recording were subjected to the subsequent analysis.
Data analysis: Pearson's correlation and cross-approximate entropy
Cross-approximate entropy (XApEn) was introduced by Pincus for the analysis of hormone secretion processes [73]. It is an extension of approximate entropy which has been widely used in the analysis of biosignals such as EEG, electrocardiogram and, recently, local field potentials [29,74-76]. XApEn quantifies the predictability of the appearance of common patterns in two time series, not heeding the order in which they appear. Higher values of XApEn indicate higher degrees of independence or dissimilarity of the signals. As pointed out by Hudetz et al [10], XApEn is sensitive to changes in signal independence over a wide frequency spectrum and does not require stationarity of signals; thus it is very well suited to the analysis of local field potentials. In a preprocessing step, the selected data sequences were low pass filtered by a MATLAB 6.5 (The MathWorks, Inc., Natick, MA, USA) double reverse butterworth filtering routine (-3dB frequency of 200 Hz), which allows filtering with zero phase shift, and subsequently downsampled to 500 Hz. Next, data from each channel were normalized by subtracting the mean, followed by a division by the standard deviation. XApEn was calculated according to Pincus et al. [73]:
equation M7
Φm(r)(x || y)is the average of equation M8. equation M9 is the number of times which a sequence of defined length m in signal x starting at data point i has a similar counterpart anywhere in signal y, divided by N-m+1 (the number of comparisons possible). Two sequences are defined as similar if none of their scalar component differences (xi - yi) exceeds tolerance r (see Figure Figure22 for a detailed illustration of the computations). For analysis of the presented data, length m was set to 1 and tolerance r was 20% of the standard deviation of the signal in the channel combination with the lower channel number, i.e. if the channel combination was [1,4], SD was calculated for the signal recorded from channel 1. These settings are in accordance with the settings recommended by Pincus [73]. XApEn calculation was performed with MATLAB 6.5.
In addition to XApEn, Pearson cross correlation coefficient defined as
equation M10
where equation M11 and equation M12 are the mean of the intervals x = [x1, ..., xn] and y = [y1, ..., yn], was calculated using LabView 6i (National Instruments, Austin, Tx, USA). R [77]http://www.r-project.org/ was used for generating Figure Figure22.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
BA and GS conceived of the study and participated in its coordination. HH and CS performed animal surgery and the experiments. CS wrote data conversion code. MK performed cross-approximate entropy and cross correlation analysis and together with HH prepared the figures and performed statistical analyses. GS, BA, EK and CS participated in the preparation of the manuscript and interpretation of data. HH and MK wrote the paper. All authors read and approved the final manuscript.
Additional file 1
Figure S1A depicts the dose response relationships of XApEn for the same data as used in Figure Figure5,5, but filtered with an upper cutoff frequency of 30 Hz and resampled at 100 Hz. Same conventions as for Figure Figure55 apply, including monoexponential fits represented by solid lines. Qualitatively, the same observations as with the original data set, containing frequencies up to 200 Hz, can be made: decline of XApEn with anesthetic concentration, significant difference between halothane on the one hand and enflurane and isoflurane on the other hand. However, the data are more variable among animals and show less decline with anesthetic frequency. The free parameter tau was [95% confidence intervals]: 2.37 [1.34 3.41] for isoflurane, 2.30 [1.50 3.10] for enflurane and 6.31 [4.67 7.94] for halothane. R2 values were 0.52, 0.68 and 0.51 (same order of agents). Figure S1B shows dose response relations of cross correlation values (normalized to control) computed from the same data as used for Figure Figure5.5. While on average the cross correlation increased with anesthetic concentration, there was little dose dependence, and intra- and interindividual variability was extremely high (cf. Figure Figure33).
Acknowledgements
This work was supported by the German Ministry for Education and Research (BMBF 0311858). Further support was provided by the Hertie Foundation and the Hermann and Lilly Schilling Foundation. The authors thank Kuno Kirschfeld for generous support of this study, and Susanne Kramer for excellent technical support.
  • Alkire MT, Hudetz AG, Tononi G. Consciousness and anesthesia. Science. 2008;322:876–880. doi: 10.1126/science.1149213. [PMC free article] [PubMed] [Cross Ref]
  • Varela F, Lachaux J-P, Rodriguez E, Martinerie J. The brainweb: phase synchronization and large-scale integration. Nat Rev Neurosci. 2001;2:229–239. doi: 10.1038/35067550. [PubMed] [Cross Ref]
  • Lamme VAF. Towards a true neural stance of consciousness. Trends Cogn Sci. 2006;10:494–501. doi: 10.1016/j.tics.2006.09.001. [PubMed] [Cross Ref]
  • Buzsáki G. The structure of consciousness. Nature. 2007;446:267. doi: 10.1038/446267a. [PubMed] [Cross Ref]
  • Uhlhaas PJ, Pipa G, Lima B, Melloni L, Neuenschwander S, Nikolic D, Singer W. Neural synchrony in cortical networks: history, concept and current status. Frontiers in Integrative Neuroscience. 2009;3:17. doi: 10.3389/neuro.07.017.2009. [PMC free article] [PubMed] [Cross Ref]
  • Imas OA, Ropella KM, Wood JD, Hudetz AG. Isoflurane disrupts anterio-posterior phase synchronization of flash-induced field potentials in the rat. Neurosci Lett. 2006;402:216–221. doi: 10.1016/j.neulet.2006.04.003. [PubMed] [Cross Ref]
  • Hudetz AG, Vizuete JA, Imas OA. Desflurane selectively suppresses long-latency cortical neuronal response to flash in the rat. Anesthesiology. 2009;111:231–239. doi: 10.1097/ALN.0b013e3181ab671e. [PMC free article] [PubMed] [Cross Ref]
  • John ER, Prichep LS, Kox W, Valdés-Sosa P, Bosch-Bayard J, Aubert E, Tom E, diMichele F, Gugino LD. Invariant reversible QEEG effects of anesthetics. Conscious Cogn. 2001;10:165–183. doi: 10.1006/ccog.2001.0507. [PubMed] [Cross Ref]
  • Ferrarelli F, Massimini M, Sarasso S, Casali A, Riedner BA, Angelini G, Tononi G, Pearce RA. Breakdown in cortical effective connectivity during midazolam-induced loss of consciousness. Proc Natl Acad Sci USA. 2010;107:2681–2686. doi: 10.1073/pnas.0913008107. [PubMed] [Cross Ref]
  • Hudetz AG. Effect of volatile anesthetics on interhemispheric EEG cross-approximate entropy. Brain Res. 2002;954:123–131. doi: 10.1016/S0006-8993(02)03358-9. [PubMed] [Cross Ref]
  • Hudetz AG, Wood JD, Kampine JP. Cholinergic reversal of isoflurane anesthesia in rats as measured by cross-approximate entropy of the electroencephalogram. Anesthesiology. 2003;99:1125–1131. doi: 10.1097/00000542-200311000-00019. [PubMed] [Cross Ref]
  • Song S, Sjöström PJ, Reigl M, Nelson S, Chklovskii DB. Highly nonrandom features of synaptic connectivity in local cortical circuits. PLoS Biol. 2005;3:e68–e80. doi: 10.1371/journal.pbio.0030068. [PMC free article] [PubMed] [Cross Ref]
  • Kampa BM, Letzkus JJ, Stuart GJ. Cortical feed-forward networks for binding different streams of sensory information. Nat Neurosci. 2006;9:1472–1473. doi: 10.1038/nn1798. [PubMed] [Cross Ref]
  • Simons DJ, Carvell GE, Hershey AE, Bryant P. Responses of barrel cortex neurons in awake rats and effects of urethane anesthesia. Exp Brain Res. 1992;91:259–272. doi: 10.1007/BF00231659. [PubMed] [Cross Ref]
  • Erchova IA, Lebedev MA, Diamond ME. Somatosensory cortical neuronal population activity across states of anaesthesia. Eur J Neurosci. 2002;15:744–752. doi: 10.1046/j.0953-816x.2002.01898.x. [PubMed] [Cross Ref]
  • Hentschke H, Schwarz C, Antkowiak B. Neocortex is the major target of sedative concentrations of volatile anesthetics: strong depression of firing rates and increase of GABAA receptor-mediated inhibition. Eur J Neurosci. 2005;21:93–102. doi: 10.1111/j.1460-9568.2004.03843.x. [PubMed] [Cross Ref]
  • Berger T, Borgdorff A, Crochet S, Neubauer FB, Lefort S, Fauvet B, Ferezou I, Carleton A, Lüscher H-R, Petersen CCH. Combined voltage and calcium epiflourescence imaging in vitro and in vivo reveals subthreshold and suprathreshold dynamics of mouse barrel cortex. J Neurophysiol. 2007;97:3751–3762. doi: 10.1152/jn.01178.2006. [PubMed] [Cross Ref]
  • Sleigh JW, Vizuete JA, Voss L, Steyn-Ross A, Steyn-Ross M, Marcuccilli CJ, Hudetz AG. The electrocortical effects of enflurane: experiment and theory. Anesth Analg. 2009;109:1253–1262. doi: 10.1213/ANE.0b013e3181add06b. [PMC free article] [PubMed] [Cross Ref]
  • Blatow M, Rozov A, Katona I, Hormuzdi SG, Meyer AH, Whittington MA, Caputi A, Monyer H. A novel network of multipolar bursting interneurons generates theta frequency oscillations in neocortex. Neuron. 2003;38:805–817. doi: 10.1016/S0896-6273(03)00300-3. [PubMed] [Cross Ref]
  • Hasenstaub A, Shu Y, Haider B, Kraushaar U, Duque A, McCormick DA. Inhibitory postsynaptic potentials carry synchronized frequency information in active cortical networks. Neuron. 2005;47:423–435. doi: 10.1016/j.neuron.2005.06.016. [PubMed] [Cross Ref]
  • Merriam EB, Netoff TI, Banks MI. Bistable network behavior of layer I interneurons in auditory cortex. J Neurosci. 2005;25:6175–6186. doi: 10.1523/JNEUROSCI.0512-05.2005. [PMC free article] [PubMed] [Cross Ref]
  • Drexler B, Roether CL, Jurd R, Rudolph U, Antkowiak B. Opposing actions of etomidate on cortical theta oscillations are mediated by different γ-aminobutyric acid type A receptor subtypes. Anesthesiology. 2005;102:346–352. doi: 10.1097/00000542-200502000-00017. [PubMed] [Cross Ref]
  • Creutzfeldt OD, Fuster JM, Lux HD, Nacimiento A. Experimenteller Nachweis von Beziehungen zwischen EEG-Wellen und Aktivitaet corticaler Nervenzellen. Naturwissenschaften. 1964;51:166–167. doi: 10.1007/BF00622299. [Cross Ref]
  • Katzner S, Nauhaus I, Benucci A, Bonini V, Ringach DL, Carandini M. Local origin of field potentials in visual cortex. Neuron. 2009;61:35–41. doi: 10.1016/j.neuron.2008.11.016. [PMC free article] [PubMed] [Cross Ref]
  • Kandel A, Buzsáki G. Cellular-synaptic generation of sleep spindles, spike-and-wave discharges, and evoked thalamocortical responses in the neocortex of the rat. J Neurosci. 1997;17:6783–6797. [PubMed]
  • Haslinger R, Ulbert I, Moore CI, Brown EN, Devor A. Analysis of LFP phase predicts sensory response of barrel cortex. J Neurophysiol. 2006;96:1658–1663. doi: 10.1152/jn.01288.2005. [PubMed] [Cross Ref]
  • Antunes LM, Golledge HD, Roughan JV, Flecknell PA. Comparison of electroencephalogram activity and auditory evoked responses during isoflurane and halothane anaesthesia in the rat. Veterinary Anaesthesia and Analgesia. 2003;30:15–23. doi: 10.1046/j.1467-2995.2003.00085.x. [PubMed] [Cross Ref]
  • Orth M, Bravo E, Barter L, Carstens E, Antognini JF. The differential effects of halothane and isoflurane on electroencephalographic responses to electrical microstimulation of the reticular formation. Anesthesia & Analgesia. 2006;102:1709–1714. [PubMed]
  • Silva A, Cardoso-Cruz H, Silva F, Galhardo V, Antunes LM. Comparison of anesthetic depth indexes based on thalamocortical local field potentials in rats. Anesthesiology. 2010;112:355–363. doi: 10.1097/ALN.0b013e3181ca3196. [PubMed] [Cross Ref]
  • Stoelting RK, Longnecker DE, Eger EI. Minimum alveolar concentrations in man on awakening from methoxyflurane, halothane, ether and fluroxene anesthesia: MAC awake. Anesthesiology. 1970;33:5–9. doi: 10.1097/00000542-197007000-00004. [PubMed] [Cross Ref]
  • Katoh T, Suguro Y, Kimura T, Ikeda K. Cerebral awakening concentration of sevoflurane and isoflurane predicted during slow and fast alveolar washout. Anesth Analg. 1993;77:1012–1017. doi: 10.1213/00000539-199311000-00024. [PubMed] [Cross Ref]
  • Gaumann DM, Mustaki JP, Tassonyi E. MAC-awake of isoflurane, enflurane and halothane evaluated by slow and fast alveolar washout. Br J Anaesth. 1992;68:81–84. doi: 10.1093/bja/68.1.81. [PubMed] [Cross Ref]
  • Alkire MT, Haier RJ, Shah NK, Anderson CT. Positron emission tomography study of regional cerebral metabolism in humans during isoflurane anesthesia. Anesthesiology. 1997;86:549–557. doi: 10.1097/00000542-199703000-00006. [PubMed] [Cross Ref]
  • Alkire MT, Pomfrett CJ, Haier RJ, Gianzero MV, Chan CM, Jacobsen BP, Fallon JH. Functional brain imaging during anesthesia in humans: effects of halothane on global and regional cerebral glucose metabolism. Anesthesiology. 1999;90:701–709. doi: 10.1097/00000542-199903000-00011. [PubMed] [Cross Ref]
  • Woolsey TA, van der Loos H. The structural organization of layer IV in the somatosensory region of mouse cerebral cortex. The description of a cortical field composed of discrete cytoarchitectonic units. Brain Res. 1970;17:205–242. doi: 10.1016/0006-8993(70)90079-X. [PubMed] [Cross Ref]
  • Chmielowska J, Carvell GE, Simons DJ. Spatial organization of thalamocortical and corticothalamic projection systems in rat SmI barrel cortex. J Comp Neurol. 1989;285:325–338. doi: 10.1002/cne.902850304. [PubMed] [Cross Ref]
  • Welker E, Hoogland PV, van der Loos H. Organization of feedback and feedforward projections of the barrel cortex: a PHA-L study in the mouse. Exp Brain Res. 1988;73:411–435. doi: 10.1007/BF00248234. [PubMed] [Cross Ref]
  • Petersen RS, Diamond ME. Spatial-temporal distribution of whisker-evoked activity in rat somatosensory cortex and the coding of stimulus location. J Neurosci. 2000;20:6135–6143. [PubMed]
  • Laaris N, Keller A. Thalamic-evoked synaptic interactions in barrel cortex revealed by optical imaging. J Neurosci. 2000;20:1529–1537. [PubMed]
  • Petersen CCH, Grinvald A, Sakmann B. Spatiotemporal dynamics of sensory responses in layer 2/3 of rat barrel cortex measured in vivo by voltage-sensitive dye imaging combined with whole-cell voltage recordings and neuron reconstructions. J Neurosci. 2003;23:1298–1309. [PubMed]
  • Braitenberg V, Schüz A. Cortex: Statistics and Geometry of Neuronal Connectivity. Berlin, Heidelberg: Springer-Verlag; 1998.
  • Hoeflinger BF, Bennett-Clarke CA, Chiaia NL, Killackey HP, Rhoades RW. Patterning of local intracortical projections within the vibrissae representation of rat primary somatosensory cortex. J Comp Neurol. 1995;354:551–563. doi: 10.1002/cne.903540406. [PubMed] [Cross Ref]
  • Cauller LJ, Clancy B, Connors BW. Backward cortical projections to primary somatosensory cortex in rats extend long horizontal axons in layer I. J Comp Neurol. 1998;390:297–310. doi: 10.1002/(SICI)1096-9861(19980112)390:2<297::AID-CNE11>3.0.CO;2-V. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/j.neuron.2006.03.043. [PubMed] [Cross Ref]
  • Kerr JN, de Kock CPJ, Greenberg DS, Bruno RM, Sakmann B, Helmchen F. Spatial organization of neuronal population responses in layer 2/3 of rat barrel cortex. J Neurosci. 2007;27:13316–13328. doi: 10.1523/JNEUROSCI.2210-07.2007. [PubMed] [Cross Ref]
  • MacIver MB, Mikulec AA, Amagasu SM, Monroe FS. Volatile anesthetics depress glutamate transmission via presynaptic actions. Anesthesiology. 1996;85:823–834. doi: 10.1097/00000542-199610000-00018. [PubMed] [Cross Ref]
  • Vardya I, Drasbek KR, Dósa Z, Jensen K. Cell type-specific GABAA receptor-mediated tonic inhibition in mouse neocortex. J Neurophysiol. 2008;100:526–532. doi: 10.1152/jn.01224.2007. [PubMed] [Cross Ref]
  • Chen LM, Friedmann RM, Roe AW. Optical imaging of SI topography in anesthetized and awake squirrel monkeys. J Neurosci. 2005;25:7648–7659. doi: 10.1523/JNEUROSCI.1990-05.2005. [PubMed] [Cross Ref]
  • Kuhn B, Denk W, Bruno RM. In vivo two-photon voltage-sensitive dye imaging reveals top-down control of cortical layers 1 and 2 during wakefulness. Proc Natl Acad Sci USA. 2008;105:7588–7593. doi: 10.1073/pnas.0802462105. [PubMed] [Cross Ref]
  • Lamme VAF, Zipser K, Spekreijse H. Figure-ground activity in primary visual cortex is suppressed by anesthesia. Proc Natl Acad Sci USA. 1998;95:3263–3268. doi: 10.1073/pnas.95.6.3263. [PubMed] [Cross Ref]
  • Antkowiak B, Helfrich-Förster C. Effects of small concentrations of volatile anesthetics on action potential firing of neocortical neurons in vitro. Anesthesiology. 1998;88:1592–1605. doi: 10.1097/00000542-199806000-00024. [PubMed] [Cross Ref]
  • Alkire MT, Haier RJ, Fallon JH. Toward a unified theory of narcosis: brain imaging evidence for a thalamocortical switch as the neurophysiologic basis of anesthetic-induced unconsciousness. Conscious Cogn. 2000;9:370–386. doi: 10.1006/ccog.1999.0423. [PubMed] [Cross Ref]
  • Steriade M, McCormick DA, Sejnowski TJ. Thalamocortical oscillations in the sleeping and aroused brain. Science. 1993;262:679–685. doi: 10.1126/science.8235588. [PubMed] [Cross Ref]
  • Detsch O, Vahle-Hinz C, Kochs E, Siemers M, Bromm B. Isoflurane induces dose-dependent changes of thalamic somatosensory information transfer. Brain Res. 1999;829:77–89. doi: 10.1016/S0006-8993(99)01341-4. [PubMed] [Cross Ref]
  • Gottschalk A, Miotke SA. Volatile anesthetic action in a computational model of the thalamic reticular nucleus. Anesthesiology. 2009;110:996–1010. doi: 10.1097/ALN.0b013e31819db923. [PubMed] [Cross Ref]
  • Détári L, Rasmusson DR, Semba K. Phasic relationship between the activity of basal forebrain neurons and cortical EEG in urethane-anesthetized rat. Brain Res. 1997;759:112–121. doi: 10.1016/S0006-8993(97)00252-7. [PubMed] [Cross Ref]
  • Manns ID, Alonso A, Jones B. Rhythmically discharging basal forebrain units comprise cholinergic, GABAergic, and putative glutamatergic cells. J Neurophysiol. 2003;89:1057–1066. doi: 10.1152/jn.00938.2002. [PubMed] [Cross Ref]
  • Goard M, Dan Y. Basal forebrain activation enhances cortical coding of natural scenes. Nat Neurosci. 2009;12:1440–1445. doi: 10.1038/nn.2402. [PMC free article] [PubMed] [Cross Ref]
  • Zecharia AY, Nelson LE, Gent TC, Schumacher M, Jurd R, Rudolph U, Brickley SG, Maze M, Franks NP. The involvement of hypothalamic sleep pathways in general anesthesia: testing the hypothesis using the GABAA receptor β3N265 M knock-in mouse. J Neurosci. 2009;29:2177–2187. doi: 10.1523/JNEUROSCI.4997-08.2009. [PubMed] [Cross Ref]
  • Imas OA, Ropella KM, Wood JD, Hudetz AG. Halothane augments event-related γ oscillations in rat visual cortex. Neuroscience. 2004;123:269–278. doi: 10.1016/j.neuroscience.2003.09.014. [PubMed] [Cross Ref]
  • Lampl I, Reichova I, Ferster D. Synchronous membrane potential fluctuations in neurons of the cat visual cortex. Neuron. 1999;22:361–374. doi: 10.1016/S0896-6273(00)81096-X. [PubMed] [Cross Ref]
  • Tononi G, Edelman GM. Consciousness and complexity. Science. 1998;282:1846–1851. doi: 10.1126/science.282.5395.1846. [PubMed] [Cross Ref]
  • Jordan D, Stockmanns G, Kochs EF, Pilge S, Schneider G. Electroencephalographic order pattern analysis for the separation of consciousness and unconsciousness: an analysis of approximate entropy, permutation entropy, recurrence rate, and phase coupling of order recurrence plots. Anesthesiology. 2008;109:1014–1022. doi: 10.1097/ALN.0b013e31818d6c55. [PubMed] [Cross Ref]
  • John ER, Prichep LS. The anesthetic cascade. Anesthesiology. 2005;102:447–471. doi: 10.1097/00000542-200502000-00030. [PubMed] [Cross Ref]
  • Mazze RI, Rice SA, Baden JM. Halothane, isoflurane, and enflurane MAC in pregnant and nonpregnant female and male mice and rats. Anesthesiology. 1985;62:339–341. doi: 10.1097/00000542-198503000-00021. [PubMed] [Cross Ref]
  • Semba K, Komisaruk BR. Neural substrates of two different rhythmical vibrissal movements in the rat. Neuroscience. 1984;12:761–774. doi: 10.1016/0306-4522(84)90168-4. [PubMed] [Cross Ref]
  • Kaplan BJ. The epileptic nature of rodent electrocortical polyspiking is still unproven. Exp Neurol. 1985;88:425–436. doi: 10.1016/0014-4886(85)90204-3. [PubMed] [Cross Ref]
  • Nicolelis MAL, Fanselow EE. Thalamocortical optimization of tactile processing according to behavioral state. Nat Neurosci. 2002;5:517–523. doi: 10.1038/nn0602-517. [PubMed] [Cross Ref]
  • Wiest MC, Nicolelis MAL. Behavioral detection of tactile stimuli during 7-12 Hz cortical oscillations in awake rats. Nat Neurosci. 2003;6:913–914. doi: 10.1038/nn1107. [PubMed] [Cross Ref]
  • Meeren HKM, Pijn JPM, Luijtelaar ELJM, Coenen AML, Lopes da Silva F. Cortical focus drives widespread corticothalamic networks during spontaneous absence seizures in rats. J Neurosci. 2002;22:1480–1495. [PubMed]
  • Shaw F-Z. Is spontaneous high-voltage rhythmic spike discharge in Long Evans rats an absence-like seizure activity? J Neurophysiol. 2004;92:63–77. [PubMed]
  • Polack P-O, Charpier S. Intracellular activity of cortical and thalamic neurones during high-voltage rhythmic spike discharge in Long-Evans rats in vivo. J Physiol. 2006;571:461–476. doi: 10.1113/jphysiol.2005.100925. [PubMed] [Cross Ref]
  • Pincus SM, Mulligan T, Iranmanesh A, Gheorghiu S, Godschalk M, Veldhuis JD. Older males secrete luteinizing hormone and testosterone more irregularly, and jointly more asynchronously, than younger males. PNAS. 1996;93:14100–14105. doi: 10.1073/pnas.93.24.14100. [PubMed] [Cross Ref]
  • Bruhn J, Röpcke H, Hoeft A. Approximate entropy as an electroencephalographic measure of anesthetic drug effect during desflurane anesthesia. Anesthesiology. 2000;92:715–726. doi: 10.1097/00000542-200003000-00016. [PubMed] [Cross Ref]
  • Pincus SM. Approximate entropy as a measure of system complexity. Proc Natl Acad Sci USA. 1991;88:2297–2301. doi: 10.1073/pnas.88.6.2297. [PubMed] [Cross Ref]
  • Beckers F, Ramaekers D, Aubert AE. Approximate Entropy of Heart Rate Variability: Validation of Methods and Application in Heart Failure. Cardiovascular Engineering. 2001;1:177–182. doi: 10.1023/A:1015212328405. [Cross Ref]
  • R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2005.
Articles from BMC Neuroscience are provided here courtesy of
BioMed Central