Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Attenuation of High-Frequency (50-200 Hz) Thalamocortical EEG Rhythms by Propofol in Rats Is More Pronounced for the Thalamus than for the Cortex

Correction

5 Oct 2015: Reed SJ, Plourde G (2015) Correction: Attenuation of High-Frequency (50-200 Hz) Thalamocortical EEG Rhythms by Propofol in Rats Is More Pronounced for the Thalamus than for the Cortex. PLOS ONE 10(10): e0140087. https://doi.org/10.1371/journal.pone.0140087 View correction

Abstract

Background

Thalamocortical EEG rhythms in gamma (30-80 Hz) and high-gamma (80-200 Hz) ranges have been linked to arousal and conscious processes. To test the hypothesis that general anesthetics attenuate these rhythms, we characterized the concentration-effect relationship of propofol on the spectral power of these rhythms. In view of the ongoing debate about cortex versus thalamus as the primary site of anesthetic action for unconsciousness, we also compared the relative sensitivity of cortex and thalamus to this effect propofol.

Methods

Adult male Long-Evans rats were chronically implanted with electrodes in somatosensory (barrel) cortex and ventroposteromedial thalamus. Propofol was delivered by a computer-controlled infusion using real-time pharmacokinetic modeling to obtain the desired plasma concentration. Spectral power was assessed during baseline, at four stable propofol plasma-concentrations (0, 3,6,9,12 μg/ml) and during recovery over four frequency ranges (30-50, 51-75, 76-125, 126-200 Hz). Unconsciousness was defined as complete loss of righting reflex. Multiple regression was used to model the change of power (after logarithmic transformation) as a function of propofol concentration and recording site.

Results

Unconsciousness occurred at the 9 μg/ml concentration in all animals. Propofol caused a robust linear concentration-dependent attenuation of cortical power in the 76-200 Hz range and of thalamic power in the 51-200 Hz range. In all instances the concentration-effect slope for the thalamus was markedly steeper than for the cortex. Furthermore the lowest concentration causing unconsciousness significantly reduced cortical power in the 126-200 Hz range and thalamic power in the 51-200 Hz range.

Conclusions

Propofol causes a concentration-dependent attenuation of the power of thalamocortical rhythms in the 51-200 Hz range and this effect is far more pronounced for the thalamus, where the attenuation provides a robust correlate of the hypnotic action of propofol.

Introduction

The mechanisms by which general anesthetics suppress consciousness remain elusive [1, 2]. The hypothesis that general anesthetics impair consciousness through, in part, attenuation of high-frequency thalamocortical EEG rhythms in the gamma (30–80 Hz) range has been proposed [36]. Gamma rhythms are an important element of on-going activity of the waking brain, reflecting the depolarization of thalamic and cortical neurons [7]. They occur in all sensory areas and have been associated with a number of high-level neurocognitive processes including general arousal, sensory, and attention [8, 9]. Early observations of high-frequency rhythms were confined to the gamma range but recent developments have drawn attention to the broadband high-gamma (80–200 Hz) range [10, 11].

There is emerging evidence that anesthetics such as isoflurane and propofol attenuate high-frequency cortical rhythms above 50 Hz in humans and animal models [1214] as well as thalamocortical network dynamics [15, 16]. We recently reported [17] that propofol anesthesia in rats is associated with a decrease of 80–200 Hz power in the cortex and thalamus and that physostigmine given during continuous infusion of propofol causes behavioral arousal that correlates with an increase in 80–200 Hz power in the thalamus. Furthermore, there was a suggestion in this study that the attenuation of power in the 50–200 Hz range during propofol was more pronounced for the thalamus than for the cortex but the evidence was not sufficiently strong to warrant an explicit statement. Thus further evaluation of the effect of propofol on gamma and high-gamma rhythms on cortex and thalamus appears warranted. Demonstration of a concentration-effect relationship is the standard method to assess the effect of a drug and to provide evidence of causality. A formal comparison of the cortex and thalamus in regard to their sensitivity to propofol could inform the issue of whether anesthetic-induced unconsciousness primarily results for impairment of the cortex, thalamus or both [1].

Our primary aim was to characterize the concentration-effect relationship of propofol for the attenuation of high frequency (30–200 Hz) thalamocortical EEG rhythms and to compare the cortex and thalamus in regard to their sensitivity to propofol. The 30–200 Hz range was divided in the following bands chosen on the basis of prior studies [12, 14]: 30–50 Hz; 51–75 Hz; 76–125 Hz and 126–200 Hz. For each frequency band, we evaluated to concentration effect relationship of propofol on the power of the EEG recorded from the somatosensory (barrel) cortex and ventroposteromedial thalamic nucleus (VPM) of the rats with chronically implanted electrodes and we compared the steepness of the cortical and thalamic regression lines. We chose the ventroposteromedial thalamic nucleus (VPM) because of its documented sensitivity to general anesthetics [18] and because it is suitable for recording spontaneous fast oscillations [19]. The barrel cortex was selected because it receives inputs from the VPM and also displays fast spontaneous oscillation [20, 21].

Our secondary aim was to determine whether the attenuation of high-frequency power in the cortex or thalamus could reflect neural events that play a role in causing unconsciousness. We used the following criterion as a necessary (but not sufficient) condition: if the attenuation of power with a frequency band is linked in any way to unconsciousness, then it should be possible to demonstrate that the highest level of power for this frequency band observed during unconsciousness (i.e. the 9 and 12 μg/ml concentrations) is significantly less than the lowest level of power observed during any period when the animal is conscious. Any attenuation of power that fulfills this stringent criterion can be considered a potential correlate of the hypnotic effect of propofol and would deserve further investigation.

Materials and Methods

Ethics Statement

This study was carried out in strict accordance with the guidelines of the Canadian Council on Animal Care. The protocol was approved by the Montreal Neurological Institute Animal Care Committee. All surgery was performed under general anesthesia with ketamine and xylazine. All efforts were made to minimize suffering. Male Long-Evans rats (300–320g, n = 9) were acquired from Charles River Laboratories (Senneville, Quebec), and arrived with a right jugular vein catheter. Animals were housed individually in an enriched environment with food and water available, ad libitum. A normal light-cycle (light on from 07:00 to 19:00 Hrs) was maintained throughout the duration of the experiment. After completion of the experiments, the animals were euthanized under deep urethane anesthesia by perfusion via the ascending aorta.

Surgery

Rats were anesthetized with ketamine (50 mg/kg, ip) / xylazine (5 mg/kg, ip) and secured to a stereotaxic frame which was equipped with a heating pad and rectal temperature monitor. A midline incision was made along the scalp and the periosteum was retracted with hemostats. Bipolar electrodes made of twisted, teflon-coated stainless-steel wires (125μm diameter; vertical separation between the wire tips: 0.5–1.0 mm; AM Systems, WA) were then positioned in the VPM (Paxinos atlas [22] coordinates: A/P: -3.5, L/M: 2.7, V/D: -6.6, relative to bregma) and sensory (barrel) cortex (A/P: -2.3, L/M: 5.0, V/D: -3.6) for local field potentials (LFPs) recordings. LFPs consist of EEG activity recorded with micro-electrodes inserted in brain tissue [23]. The cortical electrode was adjusted vertically to optimize the amplitude of the field excitatory post-synaptic potentials evoked by VPM stimulation. Thicker wires (0.011μm exposed tips, AM Systems, WA) soldered to stainless steel screws were placed in the contralateral parietal bone and the ipsilateral frontal bone to serve as the reference and ground electrodes, respectively. Electrode leads were then fastened to gold-plated Amphenol pins and inserted into a nine-pin ABS connector (Ginder Scientific, Nepean ON, Canada). This assembly was then secured to the skull using acrylic dental cement. All animals were administered ketoprofen (5 mg/kg S.C. at end of surgery and repeated once a day for 2 days) and buprenorphine (0.5 mg/kg S.C. at end of surgery and repeated every 12 hours for up to 3 days if needed). Animals were routinely monitored to minimize any suffering or post-operative complications. One week was allotted for recovery and observation, whilst the catheter was flushed every 3–5 days with a gentamycin / heparin solution. No animals exhibited clinical endpoints warranting removal from the study, including significant weight loss, behavior changes, or other adverse reactions to anesthesia and surgery.

Design

The experimental sessions were preceded by a forty-minute acclimation session, which occurred daily for three days, allowing the animals to explore the testing chamber. During this time, the animal’s head assembly was attached to the commutator, but with no recording. This reduced any anxiety-provoked excitement that would confound the results.

Each rat underwent a single testing session with six conditions (Baseline—target plasma propofol concentrations of 3 μg /ml, 6 μg/ml, 9 μg/ml, 12.0 μg/ml, and recovery after return of spontaneous ambulation). The choice of these concentrations was based on previous experiments [17] and designed to cover both the sub-hypnotic (3 and 6 μg/ml) and hypnotic (9 and 12 μg/ml) ranges. Our previous work [17] had shown that unconsciousness (defined as loss of righting reflex and absence of any attempts at righting) occurred at concentrations between 6.5 and 9 μg/ml. The initial baseline condition consisted of recording LFPs and observing behavior (ambulation, grooming…) while the rat explored the environment. This sequence was repeated for the 3, 6, 9, and 12 μg/ml concentrations. The animals were closely monitored during the step-wise administration of propofol to detect any unintended effects. Finally, the propofol administration was terminated and the animals were given time to regain their righting reflex and spontaneous ambulation before the last recording (recovery).

The anesthetic was administered in the right jugular vein catheter with a Harvard-22 syringe pump controlled by the Stanpump software developed by Steven L. Shafer and colleagues (Department of Anesthesiology, Standford University, CA) using pharmacokinetic parameters derived from Knibbe et al. [24]. After a 15 minutes drug equilibration period, LFPs were recorded and spontaneous behavior was observed. If the animal was not ambulating, the righting reflex was tested by placing the animal on its side and by observing attempts to right. The animal was considered unconscious if no righting attempt was made.

Electrophysiological Recording

Local field potentials from the cortical and thalamic contacts were recorded with a common reference on the contralateral parietal bone (referential montage), amplified (0.1 Hz to 475 Hz band pass), digitized at 3030 Hz, and stored for offline analysis. For each period, two minutes of high quality data devoid of artifacts was obtained. The recordings were subsequently reformatted offline to obtain a bipolar recording for each recording site (S1Bf and VPM) consisting of the difference between the recordings of the two contacts (spaced 1 mm from each other) of each electrode. All analyses were conducted with the bipolar recordings because they are less subject to contamination from volume-conducted myogenic artifacts.

Spectral power was computed with the Welch’s method [25] using 2-second long non-overlapping segments and a Hamming window (Matlab Signal Processing toolbox, version 6, MathWorks Inc. Sherborn, MA). Segments containing more than 5% outliers (defined as values outside mean ± three standard deviations of the entire recording of about 2 minutes) were excluded from analysis. The frequency bands were defined as: 30–50, 51–75, 76–125, 126–200 Hz. To minimize the impact of interference from external electrical sources, power at 59–61 Hz, 119–121 Hz, and 179–181 Hz were excluded from analysis. Since baseline power differs between animals, normalization of the spectra is required to avoid an undue influence of animals with high baseline power. The power spectra were normalized for each animal and recording site by dividing the power values by the average power from 0.1–200 Hz during baseline. Thus, for each animal and recording site, the average power during baseline was equal to unity and the spectra for the other periods were scaled with the same parameter.

Perfusion and Histology

The animals were euthanized under deep urethane anesthesia (2000 mg/kg, I.P) by perfusion via the ascending aorta of 200 ml of heparinized saline followed by 500 ml of 10% neutral buffered formalin. The isolated brains were post-fixed in 4% paraformaldehyde followed by a paraffin embedding. Histological preparations for the first five rats yielded tissue unsuitable for microscopy. The other four brains were processed by the Goodman Cancer Research Centre’s Histology Facility at McGill University. Slices were made in 50 μm step-sections at a thickness of 6 μm and stained with haematoxylin and eosin.

Statistical Analysis

All power values were transformed with log base 10. Multiple regression was used to model the changes of log power as a linear function of concentration with a common slope for all animals. The recovery period was excluded from regression analysis. A separate analysis was conducted for each frequency band (30–50, 51–75, 76–125, 126–200 Hz). For each analysis, there were 5 repeated measurements (corresponding to the five 5 propofol concentrations; xj, j = 1, 2, 3, 4, 5 corresponding to 0, 3, 6, 9, 12 μg/ml) for each recording site (k = 1for cortex and k = 2 for thalamus).

We modeled each outcome yijk as: < yijk = αik + β0k + β1kxj + εijk > where yijk denotes the spectral power of rat i (i = 1,2,…,9) at concentration j for site k; αik is the difference between the random intercept of rat i and β0k; β0k is the regression constant for site k; β1k is the regression coefficient for the linear term (i.e. slope) for site k; εijk is the error term, which accounts for unexplained errors (e.g., measurement error).

The model contains both fixed effects (the influence of concentration and recording site) and a random effect (the baseline values of each rat). This mixed model was estimated with PROC MIXED in SAS Statistical Software (version 9.2, SAS Institute Inc., Cary, NC) using an autoregressive covariance structure. The choice of a mixed, more complex regression model is required because the experimental design involves serial (repeated) measures in then same animals (see reference [26] for a detailed account).

Inspection of the residuals revealed no anomalies. The p values for the linear term were corrected for multiple comparisons (8 tests) with Hommel's procedure [27]. The goodness of fit was assessed with the concordance correlation coefficient (CCC), which was developed for mixed models as an alternative to the traditional R2 in linear regression [28, 29]. This analysis was done for each site separately. The value of CCC ranges from -1 to 1 with a value of -1 reflecting a perfect discordance (similar to a correlation coefficient of -1), a value of 1, perfect concordance and a value of 0, no relationship at all. In the present case we are only interested in the 0 to 1 range.

We compared the cortical and thalamic regression coefficients for the linear term (β1) for each frequency range to assess the significance of differences in the slope to concentration-effect plots by modeling signals from two sites simultaneously based on approaches for multivariate longitudinal data [30]. We tested the null hypothesis that (β1Cortex—β1Thalamus) = 0 assuming that (β1Cortexβ1Thalamus) / SE (β1Cortexβ1Thalamus) follows a normal distribution and where SE (β1Cortexβ1Thalamus) is the standard error of the difference of the coefficients. (We use the underline (instead of the caret (^)) to denote the estimate of a parameter e.g. β1Cortex is the estimate, derived from the data, of β1Cortex.) To obtain the standard error, we used the variance of (β1Cortexβ1Thalamus) computed with the standard formula [31] su-v2 = su2 + sv2–2 • r • su • sv where su-v2 is the variance of the difference between the variables u and v; su and sv, their standard deviation and r the correlation coefficient between u and v. These tests were adjusted for multiple comparisons (4 tests) with Hommel’s procedure [27] implemented in SAS Statistical Software (version 9.2, SAS Institute Inc., Cary, NC).

To identify power changes that could potentially be linked to loss of consciousness, we used paired t-tests to compare the concentration with highest mean spectral power (estimated from the entire LFP epoch of 2 minutes and averaged across animals) observed during unconsciousness (i.e. the 9 and 12 μg/ml concentrations) with the period (baseline, 3 or 6 μg/ml concentrations or recovery) with the lowest mean level of power observed when the animal is conscious. This was done for each frequency band and recording site separately, with adjustment for multiple comparisons (8 tests). Paired t-tests were also used to assess the difference in power between the 12 μg/ml concentration and recovery (transition from anesthetized to awake) and were similarly corrected for multiple comparisons (8 tests). Criterion for significance was p≤0.05 with the two-tailed hypothesis.

Results

Behavioral Observations

Propofol had a consistent effect on behavior. At the 3 μg/ml concentration, the animals were lethargic with reduced ambulation compared with baseline. The 6 μg/ml concentration was associated with increased ambulation compared with baseline combined with frequent rearings and occasional jumps. The 9 μg/ml abolished righting attempts in all animals, which remained anesthetized until the termination of the propofol infusion. The animals regained consciousness with spontaneous ambulation 8–24 minutes (mean 14.8; std 6.1) after the termination of the propofol infusion.

Concentration-effect relationships

Fig 1 shows that propofol caused a significant, linear decrease of logarithmic power as a function of propofol concentration for the cortex in the 76–125 Hz and 125–200 Hz frequency ranges as well as for the thalamus in all frequency ranges. More importantly, the slope of the regression line (β1 parameter in Table 1) for the thalamus is significantly steeper than for the cortex in all frequency ranges (Table 1). This was a robust effect with p values of 0.0018 for the 126–200 Hz ranges at p<0.0001 for the other frequency ranges. It should be noted that comparisons between slopes for power in the 30–50 Hz and 51–75 Hz ranges are valid even if cortical power in these ranges did not change as a function of propofol concentration (i.e. the slope of the regression line for cortical power did not differ from zero).

thumbnail
Fig 1. Power as a function of propofol concentration for each frequency band for the cortex and thalamus.

Each animal represented by a thin line. The thick black line shows the mean. The dash red line shows the regression fit. Recovery period was not included in the regression. The significance of the regression model is shown in the ellipses: NS: non-significant; *: p<0.05; **: p<0.01; ***: p<0.001; see Table 1 for exact p values and regression parameters estimates.

https://doi.org/10.1371/journal.pone.0123287.g001

Table 1 presents in detail the results of the regression analysis including the regression coefficients for the linear term (β1, i.e. the slope of the regression line), along with their standard error and the probability that they differ from zero, as well as the concordance correlation coefficient (CCC), which provides a measure of how well the regression parameters fit the data. As expected when the regression analysis yields negative results, the CCCs are low (<0.2) for cortical power in the 30–50 Hz and 51–75 Hz ranges. In all other instances, the CCCs are all above 0.5, an indication that the regression analysis provided a satisfactory fit of the data.

Close inspection of Fig 1 suggests that there may have been a peak increase in power at the 6 μg/ml. This was most evident for the 30–50 Hz and 51–75 Hz bands for both cortex and thalamus as well as in the 76–125 Hz band for the cortex. To determine if this increase was significant, power at the 6 μg/ml concentration was compared with the average of baseline and 12 μg/ml power. There were significant outcomes before correction multiple comparisons for the 30–50 Hz bands for both recording sites as well as for cortical power in the 51–75 Hz band (Table 2). However, no difference remained significant after comparison for multiple comparisons (Table 2).

thumbnail
Table 2. Test hypothesis of a peak of log power at 6 μg/ml: 6 μg/ml versus mean of baseline and 12 μg/ml (mean ± std).

https://doi.org/10.1371/journal.pone.0123287.t002

The concentration-dependent effect of propofol on the power of fast rhythms is readily noticeable on the average spectra with all animals combined (Fig 2 top panels) as well as on the power spectra from a single animal (Fig 2 bottom panels). The distance between the traces in these spectra is greater for the thalamus than for the cortex, an observation that is in line with the steeper negative slope of the concentration-effect plots for the thalamus.

thumbnail
Fig 2. Power spectra (log-log scales).

Top panel. Average spectra for all 9 animals for the cortex and thalamus. Bars show standard error. The blanks represent power excluded from analysis to minimize the impact of interference from external electrical sources, at 60, 120 and 180 Hz. Recovery data not shown to avoid distracting overlap. Bottom panel. Power spectra from one animal for the cortex and thalamus.

https://doi.org/10.1371/journal.pone.0123287.g002

Segments of the original recordings from one animal are shown in Fig 3. The cortical traces show an increase in fast activity at the 3 μg/ml concentration. The other concentrations caused a progressive decrease of fast activity and a progressive increase in slow activity The thalamic traces show minimal changes between baseline and the 3 μg/ml concentration. The 6, 9, and 12 μg/ml concentrations were associated with a progressive increase in slow activity and decrease of fast activity.

thumbnail
Fig 3. Segments of the original recordings from the same animal as Fig 2.

The cortical traces show an increase in fast activity at the 3 μg/ml concentration. The 6, 9, and 12 μg/ml concentrations were associated with a progressive increase in slow activity and decrease of fast activity. The thalamic traces show minimal changes between baseline and the 3 μg/ml concentration. The 6, 9, and 12 μg/ml concentrations were associated with a progressive increase in slow activity and decrease of fast activity.

https://doi.org/10.1371/journal.pone.0123287.g003

Relationship with level of consciousness

Our criterion for considering that the attenuation of power within a frequency band is related to the hypnotic effect is that the highest level of power within this frequency band observed during the periods when the animal is unconscious (i.e. the 9 and 12 μg/ml concentrations) should be less than the lowest level of power observed during any period when the animal is conscious. For cortical power, only the 126–200 Hz range satisfied this condition. Its highest value during unconsciousness (9 μg/ml concentration) was significantly less (p = 0.0476) than the lowest power during consciousness (recovery) (Table 3). For the 30–50 Hz range, highest cortical power during unconsciousness was unexpectedly significantly higher (p = 0.0160) than lowest power during consciousness (baseline). In the thalamus, highest power during unconsciousness (invariably at the 9 μg/ml concentration) was significantly reduced compared with the lowest power during consciousness (6 μg/ml concentration) for all frequency bands except the 30–50 Hz band (0.0112≤p≤0.0240) (Table 3).

Recovery

The transition from the last concentration studied (12 μg/ml) to recovery was associated with a significant (0.0021≤p≤0.0162) increase in power in all instances (Table 4). Recovery recordings were obtained to confirm a return towards baseline power values to rule out other possible causes for the observed power changes (duration of experiment, equipment malfunction).

thumbnail
Table 4. To test the hypothesis that power during recovery is higher that during the 12 μg/ml concentration.

https://doi.org/10.1371/journal.pone.0123287.t004

Histology

Accurate electrode placement was confirmed by microscopy in the four brains available (Fig 4).

thumbnail
Fig 4. Histology.

The filled circles denote the end of the electrode tracks on an overlay from the Paxinos atlas [22] and confirm correct electrode placement in the 4 available specimens. Distances are relative to bregma.

https://doi.org/10.1371/journal.pone.0123287.g004

Discussion

The main finding of this study is that propofol causes a more extensive and pronounced attenuation of high-frequency rhythms in the sensory thalamus than in the sensory cortex. The attenuation of power of the high-frequency rhythms occurred over a wider frequency range in the thalamus (where all four frequency bands were affected) than in the cortex (where only the 76–125 Hz and 126–200 Hz bands were affected). More importantly, the slope of the concentration-effect line was steeper for the thalamus than for the cortex for all frequency bands and the differences were highly significant (p < 0.0018 for the 126–200 Hz range; p < 0.0001 for the other bands.

Another important finding of this study is that the attenuation of thalamic power appears to provide a more consistent marker of unconsciousness than cortical power. Thalamic power the 51–200 Hz range passed a fairly stringent test requiring that there be a significant difference between the highest level of power observed during the periods when the animal is unconscious (i.e. the 9 and 12 μg/ml concentrations) and the lowest level of power observed during the periods when the animal is conscious (baseline, propofol 3μg/ml or 6 μg/ml, recovery). For cortical power, only the 126–200 Hz range met these criteria. Further work on these rhythms appears warranted to assess their usefulness as potential markers of anesthetic-induced unconsciousness.

In regard to the issue of whether anesthetic-induced unconsciousness primarily results for impairment of the cortex, thalamus or both, [1] the present findings support the view that both structures are affected. This conclusion is in line with our recent report showing, on the basis of EEG recordings from the motor cortex and sensory thalamus in three patients, that induction of anesthesia with propofol is associated with concurrent alterations of cortical and thalamic activity [32].

Hudetz et al [12] had showed in rats that isoflurane attenuates spectral power in the 70–140 Hz range in neocortical and hippocampal recordings in a concentration-dependent manner. The present report confirms and extends these observations by showing that propofol exerts a similar effect on neocortical recordings. Breashears et al [14] had reported that propofol anesthesia is associated with a decrease of power in the 75–205 Hz range in human electrocorticograms. The present report shows that the same can be observed in an animal model. The importance of providing such confirmatory observations should not be underestimated. There is an increasing awareness of the need for replication in biomedical science in view of the alarmingly high proportion of findings that cannot be reproduced even though they were published in good faith by honest and competent investigators [33].

It is easy to forget that the results described above were obtained on data after logarithmic transform (base 10). Thus a significant linear term should not be interpreted as evidence that “power decreases linearly as a function of propofol concentration” but rather that “the logarithm of power decreases linearly as a function of propofol concentration”.

Our observation of a greater effect of propofol on the thalamus is in line the report of Ying et al [34]. These authors measured the concentration of propofol required to decrease synaptic responsiveness of thalamocortical relay neurons in the ventrobasal complex (VB) mouse brain slices and noted that this concentration was much less than for potentiation of GABAA receptor-elicited responses in the hippocampus or olfactory cortex, suggesting that the ventro-basal thalamic neurons are more sensitive to propofol.

The attenuation of the power of high-frequency rhythms by propofol likely reflects a decrease of neuronal firing rates. Propofol decreases spontaneous neuronal firing rates in the both cortex [35, 36] and thalamus [35] by 30% or more and high-gamma power in LFPs is strongly correlated with the average firing rate[11, 37]. Concurrent recordings of both LFPs and of action potentials with the same electrode would allow verification of this interpretation.

Propofol exerted a puzzling effect on power in the 30–50 Hz range. It caused a concentration-dependent attenuation of thalamic, but not cortical, power in this range (Table 1). During unconsciousness, however, cortical power in the 30–50 Hz range was significantly increased compared with the period of consciousness with lowest power (Table 3). This is a clear indication that power in this band does not provide a good reflection of anesthetic effect, as noted by others [12]. The inconsistent effects of general anesthetics of the 30–50 Hz band may perhaps be justified to some extent by the fact that this frequency band occupies the border zone between the low frequency range where anesthetics increase power and the high-gamma range where anesthetics decrease power [38].

The inconsistent effects of propofol (and isoflurane [12]) on 30–50 Hz power are surprising given numerous observations that anesthesia in human subjects is associated with attenuation of “40 Hz” oscillations present in sensory evoked responses [3941]. Furthermore, this effect is also concentration-dependent [3, 42, 43]. A key element to consider is that these studies did not look at spontaneous “40 Hz” activity but relied on the 40 Hz component of sensory evoked potentials [44]. The effects of isoflurane on these responses have been documented in the rat and are similar to those seen in humans [45].This confirms that there are important differences between spontaneous and evoked gamma oscillations [46].

Limitations

We acknowledge two main limitations. First, we did not assay the propofol concentrations in blood to confirm the accuracy of the target controlled infusion protocol. We believe that collecting sequential arterial blood samples for each propofol concentration was not compatible with our testing protocol. Instead, we relied on predicted concentrations based on the Stanpump software and are thus dependent on the accuracy of the kinetic model [24]. In the worst-case scenario, we can conclude that we have demonstrated a dose-dependent (as opposed to concentration-dependent) effect since the dose of propofol and the rate of administration were similar across animals. While we are aware of prior use of target controlled infusion of propofol in rodents [47], we found no such studies where blood concentrations were measured. Second, histological confirmation of proper electrode placement is only available in four animals. Problems with tissue processing made the other specimens unusable. We nevertheless feel confident that the electrodes were properly positioned in all animals because we used electrophysiological guidance for fine adjustment of the antero-posterior position of the cortical electrode by measuring the cortical field excitatory post-synaptic potentials in response to electrical stimulation of the VPM.

Conclusions

We conclude that propofol causes a concentration-dependent attenuation of the spectral power of thalamocortical rhythms in the 51–200 Hz range and that this attenuation is of sufficient magnitude to cause statistically significant reductions of power during unconsciousness. This effect is more consistent and more pronounced for the thalamus than for the cortex, suggesting greater thalamic sensitivity to propofol. Further work is warranted to determine to what extent and how this attenuation could contribute to the hypnotic action of propofol.

Supporting Information

S1 Histology. A histology slide showing the electrode tracks for cortex (matched to middle cut on Fig 4) and VPM (corresponding to the most posterior cut on Fig 4).

https://doi.org/10.1371/journal.pone.0123287.s001

(TIF)

Acknowledgments

We thank Steven Shafer, MD (Department of Anesthesiology, Columbia University, New York, NY) for help with the Stanpump software; René Deroussent for technical assistance; Xianming Tan (Biostatistics Core Facilities, McGill University Health Center Research Institute) for advice on regression analysis and SAS programming; Jean Gotman for his comments on the manuscript and his suggestion to explain inconsistent effects of general anesthetics on 30–50 Hz range; Jean-François Hardy for his comments on the manuscript.

Author Contributions

Conceived and designed the experiments: SJR GP. Performed the experiments: SJR GP. Analyzed the data: SJR GP. Wrote the paper: SJR GP. Surgery and testing: SJR. Assisted surgery and testing: GP. Wrote first draft: SJR. Contributed to final version: SJR GP.

References

  1. 1. Alkire MT, Hudetz AG, Tononi G. Consciousness and anesthesia. Science 2008;322: 876–880. pmid:18988836
  2. 2. Franks NP. General anaesthesia: from molecular targets to neuronal pathways of sleep and arousal. Nat Rev Neurosci 2008;9: 370–386. pmid:18425091
  3. 3. Madler C, Keller I, Schwender D, Pöppel E. Sensory information processing during general anaesthesia: Effect of isoflurane on auditory evoked neuronal oscillations. Br J Anaesth 1991;66: 81–87. pmid:1997064
  4. 4. Plourde G. Depth of anesthesia. Clinical use of the 40-Hz auditory steady state response. Int Anesthesiol Clin 1993;31: 107–120. pmid:8270322
  5. 5. John ER, Prichep LS. The anesthetic cascade: a theory of how anesthesia suppresses consciousness. Anesthesiology 2005;102: 447–471. pmid:15681963
  6. 6. Mashour GA. Consciousness unbound: Toward a paradigm of general anesthesia. Anesthesiology 2004;100: 428–433. pmid:14739820
  7. 7. Steriade M. Central core modulation of spontaneous oscillations and sensory transmission in thalamocortical systems. Curr Opin Neurobiol 1993;3: 619–625. pmid:8219730
  8. 8. Bartos M, Vida I, Jonas P. Synaptic mechanisms of synchronized gamma oscillations in inhibitory interneuron networks. Nat Rev Neurosci 2007;8: 45–56. pmid:17180162
  9. 9. Fries P. Neuronal gamma-band synchronization as a fundamental process in cortical computation. Annu Rev Neurosci 2009;32: 209–224. pmid:19400723
  10. 10. Crone NE, Korzeniewska A, Franaszczuk PJ. Cortical gamma responses: searching high and low. Int J Psychophysiol 2011;79: 9–15. pmid:21081143
  11. 11. Ray S, Maunsell JH. Different origins of gamma rhythm and high-gamma activity in macaque visual cortex. PLoS Biol 2011;9: e1000610. pmid:21532743
  12. 12. Hudetz AG, Vizuete JA, Pillay S. Differential effects of isoflurane on high-frequency and low-frequency gamma oscillations in the cerebral cortex and hippocampus in freely moving rats. Anesthesiology 2011;114: 588–595. pmid:21293253
  13. 13. Kortelainen J, Jia X, Seppanen T, Thakor N. Increased electroencephalographic gamma activity reveals awakening from isoflurane anaesthesia in rats. Br J Anaesth 2012;109: 782–789. pmid:22907339
  14. 14. Breshears JD, Roland JL, Sharma M, Gaona CM, Freudenburg ZV, Tempelhoff R, et al. Stable and dynamic cortical electrophysiology of induction and emergence with propofol anesthesia. Proc Natl Acad Sci U S A 2010;107: 21170–21175. pmid:21078987
  15. 15. Ni Mhuircheartaigh R, Warnaby C, Rogers R, Jbabdi S, Tracey I. Slow-wave activity saturation and thalamocortical isolation during propofol anesthesia in humans. Sci Transl Med 2013;5: 208ra148.
  16. 16. Ching S, Cimenser A, Purdon PL, Brown EN, Kopell NJ. Thalamocortical model for a propofol-induced alpha-rhythm associated with loss of consciousness. Proc Natl Acad Sci U S A 2010;107: 22665–22670. pmid:21149695
  17. 17. Reed SJ, Plourde G, Tobin S, Chapman CA. Partial antagonism of propofol anaesthesia by physostigmine in rats is associated with potentiation of fast (80–200 Hz) oscillations in the thalamus. Br J Anaesth 2013;110: 646–653. pmid:23213035
  18. 18. 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. pmid:10350532
  19. 19. Steriade M, Contreras D, Amzica F, Timofeev I. Synchronization of fast (30–40 Hz) spontaneous oscillations in intrathalamic and thalamocortical networks. J Neurosci 1996;16: 2788–2808. pmid:8786454
  20. 20. MacDonald KD, Barth DS. High frequency (gamma-band) oscillating potentials in rat somatosensory and auditory cortex. Brain Res 1995;694: 1–12. pmid:8974632
  21. 21. Jones MS, Barth DS. Spatiotemporal organization of fast (>200 Hz) electrical oscillations in rat Vibrissa/Barrel cortex. J Neurophysiol 1999;82: 1599–1609. pmid:10482773
  22. 22. Paxinos G, Watson C. The Rat Brain in Stereotaxic Coordinates. San Diego: Academic Press; 2004.
  23. 23. Buzsaki G, Anastassiou CA, Koch C. The origin of extracellular fields and currents—EEG, ECoG, LFP and spikes. Nat Rev Neurosci 2012;13: 407–420. pmid:22595786
  24. 24. Knibbe CA, Zuideveld KP, Aarts LP, Kuks PF, Danhof M. Allometric relationships between the pharmacokinetics of propofol in rats, children and adults. Br J Clin Pharmacol 2005;59: 705–711. pmid:15948935
  25. 25. Welch PD. The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging Over Short, Modified Periodograms. IEEE Trans Audio Electroacoust 1967;15: 70–73.
  26. 26. Glantz SA, Slinker BK. Primer of applied regression & analysis of variance. Toronto: McGraw-Hill Inc.; 2001.
  27. 27. Hommel G. A Stagewise rejective multiple test procedure based on a modified Bonferroni test. Biometrika 1988;75: 383–386.
  28. 28. Lin LI. A concordance correlation coefficient to evaluate reproducibility. Biometrics 1989;45: 255–268. pmid:2720055
  29. 29. Vonesh EF, Chinchilli VM, Pu K. Goodness-of-fit in generalized nonlinear mixed-effects models. Biometrics 1996;52: 572–587. pmid:10766504
  30. 30. Gao F, Thompson P, Xiong C, Miller JP. Analyzing multivariate longitudinal data using SAS. In: K. J. LeBouton, editor editors. Thirty-first Annual SAS Users Group International (SUGI) Conference, 2006, San Francisco.
  31. 31. Sokal RR, Rohlf FJ. Biometry. San Francisco: Freeman; 1969.
  32. 32. Verdonck O, Reed SJ, Hall J, Gotman J, Plourde G. The sensory thalamus and cerebral motor cortex are affected concurrently during induction of anesthesia with propofol: a case series with intracranial electroencephalogram recordings. Can J Anaesth 2014;61: 254–262. pmid:24449402
  33. 33. Russell JF. If a job is worth doing, it is worth doing twice. Nature 2013;496: 7. pmid:23552910
  34. 34. Ying SW, Goldstein PA. Propofol suppresses synaptic responsiveness of somatosensory relay neurons to excitatory input by potentiating GABA(A) receptor chloride channels. Mol Pain 2005;1: 2. pmid:15813991
  35. 35. Andrada J, Livingston P, Lee BJ, Antognini J. Propofol and etomidate depress cortical, thalamic, and reticular formation neurons during anesthetic-induced unconsciousness. Anesth Analg 2012;114: 661–669. pmid:22190559
  36. 36. Lewis LD, Weiner VS, Mukamel EA, Donoghue JA, Eskandar EN, Madsen JR, et al. Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness. Proc Natl Acad Sci U S A 2012;109: E3377–3386. pmid:23129622
  37. 37. Ray S, Crone NE, Niebur E, Franaszczuk PJ, Hsiao SS. Neural correlates of high-gamma oscillations (60–200 Hz) in macaque local field potentials and their potential implications in electrocorticography. J Neurosci 2008;28: 11526–11536. pmid:18987189
  38. 38. Pichlmayr I, Lips U, Kunkel H. The electroencephalogram in anesthesia. Fundamentals, practical applications, examples (translation by Bonatz E. and Masyk-Iversen T.). Berlin: Springer-Verlag; 1984.
  39. 39. Madler C, Pöppel E. Auditory evoked potentials indicate the loss of neuronal oscillations during general anesthesia. Naturwissenschaften 1987;74: 42–43. pmid:3561520
  40. 40. Plourde G, Picton TW. Human auditory steady-state response during general anesthesia. Anesth Analg 1990;71: 460–468. pmid:2221405
  41. 41. Kulli J, Koch C. Does anesthesia cause loss of consciousness? Trends Neurosci 1991;14: 6–10. pmid:1709532
  42. 42. Heneghan CPH, Thornton C, Navaratnarajah M, Jones JG. Effect of isoflurane on the auditory evoked response in man. Br J Anaesth 1987;59: 277–282. pmid:3828175
  43. 43. Plourde G, Villemure C, Fiset P, Bonhomme V, Backman SB. Effect of isoflurane on the auditory steady-state response and on consciousness in human volunteers. Anesthesiology 1998;89: 844–851. pmid:9778001
  44. 44. Plourde G. Auditory evoked potentials. Best Pract Res Clin Anaesthesiol 2006;20: 129–139. pmid:16634420
  45. 45. Santarelli R, Carraro L, Conti G, Capello M, Plourde G, Arslan E. Effects of isoflurane on auditory middle latency (MLRs) and steady-state (SSRs) responses recorded from the temporal cortex of the rat. Brain Res 2003;973: 240–251. pmid:12738068
  46. 46. Bullock TH. Introduction to induced rhythms: A widespread, heterogeneous class of oscillations. In: Basar E. and Bullock T. H., editors. Induced rhythms in the brain (Brain dynamics series). Boston: Birkhauser. 1992 pp. 1–26.
  47. 47. Chemali JJ, Van Dort CJ, Brown EN, Solt K. Active emergence from propofol general anesthesia is induced by methylphenidate. Anesthesiology 2012;116: 998–1005. pmid:22446983