Previous studies have demonstrated inconsistent neurophysiologic effects of ketamine, although discrepant findings might relate to differences in doses studied, brain regions analyzed, coadministration of other anesthetic medications, and resolution of the electroencephalograph. The objective of this study was to characterize the dose-dependent effects of ketamine on cortical oscillations and functional connectivity.
Ten healthy human volunteers were recruited for study participation. The data were recorded using a 128-channel electroencephalograph during baseline consciousness, subanesthetic dosing (0.5 mg/kg over 40 min), anesthetic dosing (1.5 mg/kg bolus), and recovery. No other sedative or anesthetic medications were administered. Spectrograms, topomaps, and functional connectivity (weighted and directed phase lag index) were computed and analyzed.
Frontal theta bandwidth power increased most dramatically during ketamine anesthesia (mean power ± SD, 4.25 ± 1.90 dB) compared to the baseline (0.64 ± 0.28 dB), subanesthetic (0.60 ± 0.30 dB), and recovery (0.68 ± 0.41 dB) states; P < 0.001. Gamma power also increased during ketamine anesthesia. Weighted phase lag index demonstrated theta phase locking within anterior regions (0.2349 ± 0.1170, P < 0.001) and between anterior and posterior regions (0.2159 ± 0.1538, P < 0.01) during ketamine anesthesia. Alpha power gradually decreased with subanesthetic ketamine, and anterior-to-posterior directed connectivity was maximally reduced (0.0282 ± 0.0772) during ketamine anesthesia compared to all other states (P < 0.05).
Ketamine anesthesia correlates most clearly with distinct changes in the theta bandwidth, including increased power and functional connectivity. Anterior-to-posterior connectivity in the alpha bandwidth becomes maximally depressed with anesthetic ketamine administration, suggesting a dose-dependent effect.
Ketamine is a unique anesthetic with neural effects that are distinct from more commonly-used γ-aminobutyric acid agonists
Although various electroencephalography studies of ketamine have been performed, none have assessed spectral power and connectivity at subanesthetic and anesthetic doses
Ketamine had dose-dependent effects on spectral power, functional connectivity, and directed connectivity
Anesthetic doses of ketamine resulted in markedly increased theta power across the cortex as well as increased gamma and delta power
Increased anterior-posterior connectivity in the theta bandwidth and decreased connectivity in the alpha bandwidth were specific for ketamine anesthesia
KETAMINE is a general anesthetic with unique features at the molecular, neural, and behavioral levels. Unlike the intravenous and inhaled anesthetics in common clinical use, ketamine is thought to work primarily by antagonizing N-methyl-d-aspartate (NMDA) receptors1–4 and hyperpolarization-activated, cyclic-nucleotide-gated channel 1,5,6 without strong agonist effects on synaptic γ-aminobutyric acid (GABA) receptors. Furthermore, electroencephalographic characteristics of ketamine anesthesia are distinct from those associated with anesthetics that act primarily via GABA-receptor agonism. For example, coherent frontal alpha oscillations are observed with propofol-induced unconsciousness7,8 that are believed to result from thalamocortical hypersynchrony9 and that possibly inhibit corticocortical communication.8 These same patterns are noted with ether-based volatile anesthetics during surgical levels of unconsciousness.10 Ketamine anesthesia, however, is not associated with these electroencephalographic characteristics.11,12 During anesthetic induction with ketamine, an increase in gamma power has been consistently reported,10,12,13 and when ketamine is administered in the presence of propofol or ether-based volatile anesthetics, there is a power shift to the beta bandwidth.14,15 These effects may be due to anti-NMDA-mediated disinhibition of pyramidal neurons16 and hyperpolarization-activated, cyclic-nucleotide-gated channel 1 inhibition.15
Previous studies11–13 have examined the electroencephalographic effects of ketamine, but the precise neurophysiologic correlates remain unclear for a number of reasons. For example, ketamine has often been studied in the presence of other concomitantly administered sedative and anesthetic medications,13,14,17 potentially altering or obfuscating neurophysiologic signatures specific to ketamine. Additionally, electroencephalogram studies to date have been constrained by either regionally limited analyses or low-resolution acquisition.11–13 Studies have also focused on either subanesthetic18,19 or anesthetic dosing,11–13 without being able to directly compare dose-dependent electroencephalographic characteristics. Last, there is a dearth of information regarding recovery from ketamine-induced anesthesia, as studies have tended to focus on the beginning of surgical cases, with GABAergic anesthetics often coadministered or administered directly after ketamine induction. Improving the understanding of ketamine-specific effects on the electroencephalogram may provide insights into the neuroscientific correlates of consciousness, particularly given the unique molecular, neural, and phenomenologic features of ketamine that differ from other anesthetics. Thus, in this study, we used spectral and connectivity analyses of high-density electroencephalographic recordings to characterize neurophysiologic changes associated with ketamine as a single agent during subanesthetic administration, anesthetic dosing, and a recovery period. We hypothesized that ketamine would induce distinct, dose-dependent effects on spectral and functional connectivity patterns that would distinguish between these arousal states.
Materials and Methods
This study was approved by the University of Michigan Medical School Institutional Review Board, Ann Arbor, Michigan (HUM00061087), and written informed consent was obtained from all participants before the study. All study procedures were conducted at the University of Michigan Medical School, Ann Arbor, Michigan. Ten volunteers were recruited from September 2015 through January 2016 using recruitment flyers posted throughout the medical school and University Hospital of the University of Michigan, Ann Arbor, Michigan. Phone screening was conducted by a member of the research team to review inclusion and exclusion criteria before enrollment.
Participants were considered eligible if they were American Society of Anesthesiologists physical status class I, between the ages of 20 and 40 yr, had a body mass index less than 30, and had no predictors of a difficult airway. Candidates were excluded from participation if they had cardiovascular disease, cardiac conduction abnormalities, hypertension, obstructive sleep apnea, asthma, ongoing respiratory illness, gastroesophageal reflux, history of drug use (or positive drug screen before the experiment), family history of problems with anesthesia, neurologic disorders, psychiatric disorders, or current pregnancy. The number of volunteers recruited was based on previous studies that investigated neurophysiologic correlates of anesthetic-induced unconsciousness.7,20
Experiments were conducted between the hours of 8:00 am and 1:00 pm for all participants except for one (participant No. 2, studied between 4:00 pm and 7:00 pm). Before initiation of the experiment, a full medical and anesthetic history was obtained, and a physical examination was performed. All participants fasted from food and drink for 8 h before the experiment. Peripheral intravenous catheters were placed, and American Society of Anesthesiologists standard monitors were applied before drug administration. At least two anesthesiologists were present for the full duration of the experiment. All participants underwent the following stepwise protocol, with electroencephalogram data collection throughout each period:
a 5-min eye-closed resting period before ketamine administration (baseline condition)
subanesthetic ketamine infusion (0.5 mg/kg) over 40 min (subanesthetic condition) with eyes closed, followed by 8 mg ondansetron for nausea and vomiting prophylaxis
break for completion of questionnaire (data not reported here)
anesthetic (1.5 mg/kg) bolus dose (anesthetic condition) with eyes closed
recovery period (recovery condition) with eyes closed
Participants were instructed to keep their eyes closed throughout each of these recording periods. Rather than using target-controlled infusions, we chose dosing strategies that were directly relevant to clinical care for either depression (0.5 mg/kg over 40 min) or anesthetic induction (1.5 mg/kg bolus) to enhance the translational relevance of this study. Participants were monitored through both loss of consciousness (LOC) and return of consciousness (ROC) using a previously described protocol20 ; we acknowledge that our definition of LOC relates to consciousness of the environment and does not exclude the possibility of endogenous experiences such as dreams or hallucinations. In brief, volunteers held an object in each hand that makes a sound when squeezed. An audio loop then started playing that instructed participants, in a random manner, to squeeze their right or left hand. This audio loop would deliver a command every 30 s, and the loop played before, during, and after anesthetic ketamine bolus administration. LOC was marked when subjects ceased to respond to two audio commands in a row. A scopolamine patch (1.5 mg) was also placed after ROC for further nausea and vomiting prophylaxis. Participants were then evaluated and monitored after the conclusion of the experiment. Additional antiemetics were administered if needed. Participants were discharged once they were awake, alert, and responsive; had no significant nausea or vomiting; were able to ambulate with minimal assistance; and had a responsible adult to accompany them home. Follow-up calls were made to patients both the evening after the experiment and the next day.
Electroencephalogram data were acquired with 128-channel HydroCel nets, Net Amps 400 amplifiers, and Net Station 4.5 software (Electrical Geodesics, Inc., USA). The electroencephalogram was digitized continuously at 500 Hz with a vertex reference. Per manufacturer recommendations, channel impedances were kept at less than 50 kΩ, and the net was wrapped with gauze to optimize contact between the electrodes and scalp.
Spectral Processing and Analysis.
Data processing was performed with Chronux (http://chronux.org/)21,22 and custom MATLAB (MathWorks, USA) scripts and toolboxes. All visual inspection was performed by one of the investigators (T.B.-B.) trained in dense electroencephalography review. After each recording session, the data were bandpass-filtered at 0.5 to 55 Hz (eegiltnew, firfilt plugin, zero-phase, Hamming-windowed finite impulse response filter, 3,301 points (6.6 s), 0.5 Hz transition band, 0.5 and 55 pass-band edges, −6 dB cutoff frequencies: 0.25 and 55.25) to decrease the influence of low-frequency drift and high-frequency artifacts. Electrodes on the lowest parts of the face and head were removed, leaving 98 remaining channels. Bad channels were then detected and removed using visual inspection and the rejchan and clean rawdata functions, with 86 to 91 channels retained. The channels were then average-referenced to the mean of the voltage across all remaining channels. Conditions of interest—the baseline eyes-closed period (5 min), subanesthetic infusion (0 to 38 min), anesthetic bolus dose (from LOC to 5 min after LOC, 3.8 min for one volunteer), and recovery period (5 min, collected 8 to 15 min after ROC)—were extracted from the data for each participant. These specific time epochs were chosen to include data relatively preserved from artifact and to represent stable neurophysiologic periods for each of the four experimental conditions. Large artifacts were removed after visual inspection and the rejcont function. Remaining periods of continuous data were segmented into 3-s epochs. Data epochs with remaining large artifacts were removed via a combination of visual inspection and the rejkurt, rejtrend, and eegthresh functions. For each participant, the remaining epochs from the four periods of interest were submitted to independent component analysis (runica function, infomax, extended). Independent components representing eye blink, lateral eye, muscle, facial electromyography, focal channel noise, and focal trial noise were removed from all epochs using visual inspection and the following component classification plugins: SASICA (semi-automated selection of independent components of the electroencephalogram for artifact correction),23 IC-MARC (independent components of electroencephalogram into multiple artifact classes),24 and ADJUST (automatic electroencephalogram artifact detection based on the joint use of spatial and temporal features).25 All remaining epochs for each period and each participant were then visually inspected again to remove any remaining epochs with excessive artifacts. For the four conditions (baseline, subanesthetic, anesthetic, and recovery), the mean remaining trial counts were 54 (SD = 17), 574 (SD = 48), 66 (SD = 8), and 73 (SD = 11), respectively. The previously removed channels in each data set were interpolated to the 98-channel subset. Spectral power was computed with multitaper spectral analyses (mtspecgramc function; time window: 3 s, overlap: 0.5 s, number of tapers: 3, time-bandwidth product: 5, spectral resolution: 0.25 Hz). The median absolute power (10*log10[μV2/Hz]) was calculated for each of the four experimental periods at each of five frequency bands (delta: 1 to 4 Hz, theta: 4 to 8 Hz, alpha: 8 to 13 Hz, beta: 13 to 30 Hz, gamma: 30 to 48 Hz) for all 98 channels and then for eight frontal channels centered on the Fz site (figure in Supplemental Digital Content 1, http://links.lww.com/ALN/B458). The final data presented in the spectrograms are the sequential, 3-s epochs (grand-averaged across all participants) that remained after the data-cleaning steps described above were implemented. The data were temporally sequential but not necessarily contiguous given that some epochs were removed after cleaning and artifact removal. Images presented are log-transformed (10*log10 transform of the grand average of the single-subject data). The recovery period data for one participant were unavailable, and therefore recovery period calculations were completed for only nine participants.
The topographic maps of spectral power for each experimental condition and electroencephalogram bandwidth were constructed using the topoplot function in the MATLAB (MathWorks) toolbox EEGLAB.26
Connectivity Processing and Analysis.
The selected continuous data epochs from the spectral analyses were used for connectivity analysis without channel interpolation or independent component analysis; this was done to preserve the phase information of the original signals.27,28 Additionally, eight occipital channels were removed due to excessive artifact presence, which was deemed necessary for connectivity analyses, leaving 90 remaining channels. The functional connectivity between electroencephalogram channels was examined across all experimental conditions. The undirected and directed connectivity were measured with debiased weighted phase lag index (wPLI)29 and directed phase lag index (dPLI),30 respectively. The wPLI is a measure of how the instantaneous phases of two electroencephalogram signals are phase-locked to each other. If the instantaneous phase of one signal is consistently ahead or behind of the other signal, the phases are considered locked, and wPLI = 1. However, if the signals randomly alternate between a phase lead and a phase lag relationship, there is no phase locking, and wPLI = 0. The directionality of functional connectivity between two electroencephalogram signals was determined with the asymmetry of the phase lead/lag relationship. If the instantaneous phase of one signal consistently leads that of the other signal, dPLI = 1, and with the inverse case, dPLI = −1. If there is no a bias in the phase lead/lag relationship, then dPLI = 0. In practice, both measures are robust to the volume conduction problem of scalp electroencephalography. Electroencephalogram signals were divided into 1-min-long epochs with 50% overlap, which was further divided into 2-s nonoverlapping subepochs. For each subepoch, the cross-spectral density was estimated using the multitaper method, with time-bandwidth product = 2 and the number of tapers = 3,21,22 and from these repetitions the averaged wPLI and dPLI values at variable frequencies were estimated using a custom-written function adapted from the Fieldtrip toolbox.31 To remove the bias of the measures for a given electroencephalogram data set, the shuffled-data method was used. A series (N = 20) of shuffled signal pairs were generated for each pair of electroencephalogram signals, and the wPLI and dPLI measures were calculated with the shuffled data, the mean of which was subtracted from the raw wPLI and dPLI values as the final estimation of undirected and directed connectivity.32 Final baseline and recovery time epochs were chosen as 2-min segments selected in the middle of each respective recording. Subanesthetic epochs were selected as two sequential 2-min epochs (noted in the figures as Subanes-1 and Subanes-2) that were representative of the observed spectral power during the subanesthetic infusion period. Anesthetic epochs (noted in the figures as Anes-1 and Anes-2) were chosen and displayed in the same manner.
We tested for differences between the experimental periods (baseline, subanesthetic, anesthetic, and recovery) by entering the single-subject electroencephalogram-derived values into a linear-mixed model analysis in SPSS 22 (IBM, USA). Linear-mixed model analysis allowed for statistical comparisons in the event of partially missing data from a given participant (as was the case with missing recovery period data from one volunteer, noted above). Each experimental period served as the fixed factor with participant data serving as a correlated random effect. We used restricted maximum likelihood estimation and diagonal covariance structure with heterogenous variances and zero correlation between elements. Direct pairwise comparisons were adjusted using Bonferroni’s method. We did not model repeated covariance effects, and a P value less than 0.05 was considered statistically significant. The average of median connectivity values was reported for connectivity data, as has been previously described methodologically.7 All variability estimates are in SD.
All study participants successfully completed subanesthetic and anesthetic dosing protocols; 6 of 10 (60%) experienced nausea and vomiting after ketamine anesthesia, and 5 of 10 (50%) required additional antiemetic treatment. One volunteer (1 of 10, 10%) briefly required a chin lift and jaw thrust to maintain airway patency shortly after LOC. No adverse clinical events were otherwise noted. All 10 participants experienced LOC after the 1.5 mg/kg ketamine bolus; the mean time between bolus administration and LOC was 1.23 (±0.35) min, and LOC lasted an average of 10.68 (±3.51) min.
Dose-dependent Global Effects of Ketamine—Spectral and Topographic Analysis
Prominent changes were apparent across theta, alpha, and gamma bandwidths. There was a marked increase in theta power during ketamine anesthesia (figs. 1 and 2, and Supplemental Digital Content 2, http://links.lww.com/ALN/B459, for individual participants), and this power localized to both frontal and posterior channels (fig. 2). There was a reduction in posterior alpha power from baseline through anesthetic dosing, and no increase or anteriorization of alpha power was noted with either ketamine dose (fig. 2). Increased gamma power was also apparent during the anesthetic epoch compared to other states (figs. 1 and 2).
Frontal Channel Cluster—Ketamine-induced Changes across Brain States
Frontal channels are of particular interest given their accessibility and use in the operating room. During the subanesthetic infusion, a gradual dissipation of alpha power was noted compared to baseline (mean power ± SD, 0.46 ± 0.34 and 0.88 ± 0.80 dB, respectively), and alpha power remained decreased during the periods of ketamine anesthesia (0.49 ± 0.22 dB) and recovery (0.34 ± 0.38 dB); F(3,26) = 3.679, P < 0.05 (fig. 3, A and B). With anesthetic dosing, there was a significant increase in theta bandwidth power (4.25 ± 1.90 dB) compared to baseline (0.64 ± 0.28 dB), subanesthetic (0.60 ± 0.30 dB), and recovery (0.68 ± 0.41 dB) periods; F(3,26) = 41.01, P < 0.001 (fig. 3, A and B; table 1). Theta power was not otherwise significantly modulated during nonanesthetic periods, and relative increases in theta power during ketamine anesthesia were higher than that of any other bandwidth (table 1).
Functional Connectivity Analyses
Given the significant increases in theta power and decreases in alpha power with progressive ketamine dosing, we undertook connectivity analyses in these bandwidths across each condition. Reduced functional connectivity in the alpha bandwidth has also been previously demonstrated with anesthetic ketamine dosing,11 further motivating our choice of bandwidths for analysis. Although we noted marked increases in gamma bandwidth power, electromyographic artifact might contribute to gamma activity despite our data cleaning.33,34 Thus, we did not examine gamma connectivity in this study. Furthermore, subanesthetic ketamine has been shown to increase gamma power in other paradigms,35,36 whereas increased theta power seems to correlate specifically to the anesthetized state.19,37,38
Weighted Phase Lag Index—Theta Phase Locking during Ketamine Anesthesia
Figure 4 presents the global and regional wPLI changes across each experimental condition. Figure 4A shows the average median wPLI between anterior and posterior regions in the time and frequency domains (see Supplemental Digital Content 3, http://links.lww.com/ALN/B460, for individual participants). Alpha wPLI was dominant at baseline and increased during the initial phase of the subanesthetic infusion (fig. 4A). During the anesthetic period, however, wPLI connectivity strength shifted to the theta bandwidth (fig. 4A). Figure 4B presents the topographic changes of global wPLI networks for the theta and alpha bandwidths across experimental conditions. Within the anterior region, theta wPLI was highest during the anesthetic period (mean wPLI ± SD; 0.2349 ± 0.1170) compared to all other states; F(5,44) = 7.996, P < 0.001 (fig. 3B). Theta wPLI was also highest during ketamine anesthesia between anterior and posterior regions (0.2159 ± 0.1538); F(5,44) = 4.192, P < 0.01 (fig. 3B). Alpha wPLI within the anterior region remained relatively unchanged throughout each study condition, F(5,44) = 0.697, P = 0.628, although significant changes were noted between anterior and posterior regions, with mean wPLI lowest (0.0833 ± 0.1083) during the anesthetic period; F(5,44) = 3.147, P < 0.05 (fig. 4B).
Directed Phase Lag Index—Directional Connectivity Changes during Ketamine Anesthesia
Figure 5 demonstrates how the directionality of functional connectivity changed across the experimental conditions. Figure 5A shows the average median dPLI between anterior and posterior electrodes in the time and frequency domains (see Supplemental Digital Content 3, http://links.lww.com/ALN/B460, for individual participants). An anterior-to-posterior directionality (the positive dPLIs denoted with red color) was observed in the alpha bandwidth during the baseline and subanesthetic periods. Theta shows a consistently opposite directionality from posterior-to-anterior regions (the negative dPLI denoted with blue color). Figure 5B shows the frequency-specific and state-specific topographies of electroencephalographic directionality. Alpha dPLI was maximally reduced between anterior and posterior channels during ketamine anesthesia (mean dPLI ± SD, 0.0282 ± 0.0772); F(5,44) = 2.484; P < 0.05 (fig. 4B). Alpha dPLI was also examined between prefrontal and frontal channels and was again maximally reduced during ketamine anesthesia (−0.0224 ± 0.0400); F(5,44) = 4.137; P < 0.01 (fig. 4B). Compared to baseline (0.0427 ± 0.0717), direct pairwise comparisons (with Bonferroni correction) revealed a statistically significant reduction in alpha dPLI only during each of the anesthetic epochs (Anes-1, −0.0170 ± 0.0481, P < 0.05; Anes-2, −0.0224, ± 0.0400, P < 0.05). Thus, an overall maximally reduced anterior-to-posterior connectivity pattern was demonstrated with alpha dPLI during ketamine anesthesia among these channel regions. Although a posterior-to-anterior pattern was noted with theta dPLI during ketamine anesthesia (fig. 5B), these associations did not reach statistically significant differences between anterior and posterior channels, F(5,44) = 0.177, P = 0.970, or for prefrontal and frontal channels, F(5, 44) = 0.768, P = 0.578.
This study characterized ketamine-induced electroencephalographic changes during subanesthetic, anesthetic, and recovery periods (relative to baseline) in healthy volunteers using high-density electroencephalography. During subanesthetic ketamine administration, spectral power gradually (and modestly) shifts from the alpha bandwidth to the theta bandwidth, with anterior-to-posterior connectivity (as measured by alpha dPLI) maintained. During ketamine anesthesia, however, the power shifts dramatically to the theta bandwidth, theta wPLI increases in anterior and posterior regions, and anterior-to-posterior alpha connectivity (as measured by alpha dPLI) is significantly reduced. Upon recovery, these connectivity patterns return to near-baseline levels in each bandwidth. Ketamine anesthesia was also associated with increased gamma and delta power, as previously reported,10,12,13 and frontal gamma power remained slightly elevated compared to baseline.
Our results demonstrate increased theta power correlating with ketamine anesthesia. Of note, increased theta oscillatory activity during ketamine anesthesia was first documented more than five decades ago. In the seminal human study that first examined the pharmacologic effects of ketamine in 1965, Domino et al.37 remarked, “Characteristic thetalike waves were recorded at dosage levels which produced coma; these were distributed over most of the cerebral hemisphere.” Increased theta activity has since been documented during anesthetic ketamine administration.13,15,38 A strictly receptor- and channel-based explanation for this theta resonance may be difficult to describe, because ketamine has a considerably diverse array of molecular targets.1,2,5,6,39–42 A network-level consideration of the phenomenon might provide insight regarding the consequences for information processing. In the awake state, theta and alpha frequencies interact in a circular, reciprocating pattern of “information flow” between anterior and posterior brain regions.43,44 Although oscillatory amplitude and phase are theoretically independent, amplitude (and thus power) may associate with phase to facilitate interactions with other neuronal populations within a network.45 As alpha power and anterior-to-posterior connectivity become maximally depressed during ketamine anesthesia, increased theta power and phase locking between the anterior and posterior regions may represent large-scale network disequilibrium. This is, however, a speculative interpretation, and additional studies are needed to further explore the neurobiologic relevance of such connectivity patterns. From a clinical perspective, the distinct appearance of theta power during anesthesia may be informative for clinicians who routinely administer ketamine. The appearance of increasing theta power could alert clinicians to anesthetic ketamine dosing, which may be particularly helpful in settings where subanesthetic dosing is desired (as in the treatment of depression).46,47 As the global and frontal spectrograms in this study both illustrated marked increases in theta power during ketamine anesthesia, the use of frontal electroencephalogram channels alone may be sufficient to monitor for this theta signature.
It is notable that ketamine anesthesia is associated with a maximal depression of anterior-to-posterior connectivity in the alpha bandwidth. Inhibition of frontal-to-parietal connectivity (with certain measures) has correlated with general anesthesia in humans across a diverse range of anesthetics,12,32,48,49 and our group has previously demonstrated reduced alpha directional connectivity during ketamine anesthesia with low-resolution electroencephalography.11 Recently, disruptions in NMDA– and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid–mediated frontoparietal connectivity patterns have been demonstrated during subanesthetic ketamine administration in healthy volunteers,35 who remained conscious during these experiments. This raises the possibility that disrupted frontal-to-parietal connectivity may not strictly correlate with anesthetic-induced unconsciousness. An alternative explanation, however, is that such connectivity was not maximally depressed during these subanesthetic infusion periods. Findings from our study demonstrate that alpha directional connectivity remained relatively unchanged during subanesthetic dosing, decreased dramatically during anesthetic ketamine administration, and returned to baseline levels at recovery. The connectograms (figs. 4A and 5A) demonstrate significant functional connectivity in the alpha bandwidth during baseline and subanesthetic states—despite a gradual reduction of power—and this connectivity diminishes significantly during ketamine anesthesia. Thus, there may be a dose-dependent relationship between connectivity and loss of consciousness. In fact, Untergehrer et al.50 illustrated dynamic fluctuations in functional connectivity patterns, whereby maximum functional connectivity occurs during varying time intervals on a continuous scale, with surrogates of information transfer times being fastest during unconsciousness. Accordingly, connectivity is a dynamic process and may need to cross certain real-time, continuous thresholds to correlate with varying levels of consciousness.
Multiple and unique strengths of this study are worth highlighting. No other sedative, opioid, or hypnotic medications were administered; thus, electroencephalographic effects of ketamine were not confounded by coadministration of other psychoactive medications. Additionally, an advanced-data cleaning process was used for spectral analysis. In addition to visual inspection, data from each time epoch underwent independent component analysis blind-source separation, and independent components representing eye blink, eye muscle, facial muscle, channel noise, and single-trial artifact were removed. Analytic comparisons were also made across multiple states—we were able to compare effects of ketamine before administration, at both subanesthetic and anesthetic doses, and during recovery. High-density electroencephalography was utilized, allowing for high-resolution and multiregional analyses. Last, functional connectivity based on phase relationships was also assessed with the use of wPLI and dPLI. To our knowledge, this is the first study to combine all of these methodologic strengths for the electroencephalographic analysis of ketamine in humans (table 2).
There are significant limitations to the study as well. Administering the subanesthetic ketamine infusion immediately before anesthetic dosing may influence neurophysiologic patterns observed during the anesthetic period. Our spectral data demonstrate a gradual reduction in alpha bandwidth power and increased theta power during the subanesthetic infusion. Whether this preceding subanesthetic exposure modulates subsequent neurophysiology and network connectivity during anesthetic dosing is unclear. Additionally, although a 128-channel electroencephalogram system was used, 86 to 91 channels were ultimately retained for spectral analyses after removal of artifact and bad channels. These artifact-removal strategies may also yield spectral patterns that differ from a real-world, clinical setting. For example, ocular artifact removal may have attenuated changes in delta power that might otherwise be present.51 As we wanted to optimize electroencephalogram data for connectivity analyses, we did not employ the same artifact-removal strategies (e.g., independent component analysis) as the spectral analyses, because they may interfere with phase-synchronization data.27,28 Thus, spectral and connectivity data results emerged from different data-processing steps. The administration of both ondansetron and scopolamine could potentially have modulated central neurophysiology, but we regarded this as unlikely, and both of these medications were ultimately warranted for participant comfort. Participants were instructed to keep their eyes closed during each experimental period, but sporadic eye opening was noted throughout the recording sessions. This could have impacted certain epochs of the neurophysiologic data acquisition and analysis.
There are also significant limitations to consider regarding connectivity analyses. Although connectivity measures (in this case, wPLI and dPLI) are used as surrogates for functional interactions across the cortex, this is merely an assumption. Our recent data in nonhuman primates confirm a breakdown of corticocortical information transfer during ketamine anesthesia,52 but different connectivity measures or experimental conditions can yield disparate results.43 The phase lag index can also be associated with significant intersubject variability, particularly given the complex neural processes associated with these phase-based methods; as a result, phase lag index is often not amenable to statistical averaging or traditional artifact reduction strategies.53 In terms of the general interpretation of directional connectivity results, it must be kept in mind that these measures reflect a large spatial scale and a long temporal scale, with an unclear relationship to the underlying neural spiking networks. As an intuitive example, the Dow Jones Industrial Average samples the performance of only a select number of publically traded companies. As a single value that rises or falls, it does not capture the rich dynamics and widespread stock exchange on the trading-room floor. Similarly, directed connectivity measures sample features that might reflect large-scale cortical interactions, but they do not capture the neuronal dynamics and widespread information exchange in cortical and subcortical systems. With that being said, changes in directed connectivity measures and surrogates of information transfer observed in the brain during anesthetic-induced unconsciousness have been shown to reflect fundamental network properties.45
In summary, this study used high-resolution electroencephalographic data to characterize ketamine-induced changes across multiple doses and associated levels of consciousness. Notably, increased theta power and phase locking occur between anterior and posterior regions during ketamine anesthesia, returning to baseline upon recovery. Anterior-to-posterior connectivity in the alpha bandwidth was maximally suppressed during ketamine anesthesia, and this connectivity is also restored to baseline levels during recovery. This constellation of power and connectivity results (increased theta power, as seen during rapid-eye-movement sleep, coupled with loss of anterior-to-posterior alpha connectivity, as seen during anesthesia) might contribute to the unique qualities of ketamine anesthesia.
The authors thank Yumeng Li, M.S., Center for Statistical Consultation and Research, University of Michigan, Ann Arbor, Michigan, for statistical consultation.
Supported by grant Nos. T32GM103730 (to Drs. Vlisides and Mashour) and R01GM111293 (to Dr. Mashour) from the National Institutes of Health, Bethesda, Maryland, and by the Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan.
The authors declare no competing interests.