Broadband slow-wave modulation in posterior and anterior cortex tracks distinct states of propofol-induced unconsciousness

A controversy has developed in recent years over the roles of frontal and posterior cortices in mediating consciousness and unconsciousness. Disruption of posterior cortex during sleep appears to suppress the contents of dreaming, yet activation of frontal cortex appears necessary for perception and can reverse unconsciousness under anesthesia. We used anesthesia to study how regional cortical disruption, mediated by slow wave modulation of broadband activity, changes during unconsciousness in humans. We found that broadband slow-wave modulation enveloped posterior cortex when subjects initially became unconscious, but later encompassed both frontal and posterior cortex when subjects were more deeply anesthetized and likely unarousable. Our results suggest that unconsciousness under anesthesia comprises several distinct shifts in brain state that disrupt the contents of consciousness distinct from arousal and awareness of those contents.


Results
We quantified cross-frequency coupling using a correlation between the band-pass filtered slow voltage (0.1-4 Hz) and the instantaneous amplitude of the band-pass filtered high-frequency signal (Fig. 1a).This correlation provides a one-dimensional summary of cross-frequency coupling: if the correlation is positive, it means that the high frequency amplitude is higher when the slow voltage is positive (peakmax coupling); if the correlation is negative, it means that the high frequency amplitude is higher when the slow voltage is negative (troughmax coupling).We can then vary the frequency band used for the high frequency amplitude, to see how different frequencies relate to the slow-wave.Using this approach, we can look at how peakmax and troughmax dynamics vary across electrodes and frequencies during the transition to unconsciousness.
Figure 1b shows how this plays out for a representative subject: the subject was administered propofol in increasing doses every 14 min while performing an auditory task with a button-click response.When the propofol dose was sufficiently high, the subject stopped responding to the stimulus (LOC).After 5 increasing levels of propofol, the dose was reduced and eventually the subject started responding to the stimulus again (ROC).
The correlation-based cross-frequency coupling over frontal electrodes reaffirms the basic result from previous work on the alpha rhythm: at 10 Hz, coupling to the slow-wave begins as troughmax dynamics, switches to peakmax at higher doses, and transitions back through troughmax during recovery.By looking beyond 10 Hz, however, we can see that the troughmax dynamics begin well before LOC at higher frequencies up to about 20 Hz, frequencies that show higher power during propofol sedation in which subjects remain conscious 11 .This is also consistent with the observation that patients may be conscious during troughmax dynamics 7,11,12 .
Perhaps the most striking result in the cross-frequency coupling is that the peakmax pattern extends across all frequencies (5-50 Hz).Since there are no narrow-band rhythms in the signal other than alpha and slow, this result suggests that non-rhythmic broadband activity is coupled to the slow-wave.This is consistent with an interpretation in which peakmax dynamics reflect underlying cortical up-and down-states, wherein cortical activity, whether rhythmic or not, is shut down during the trough of the slow-wave [14][15][16] .
Comparing the cross-frequency coupling between the frontal and posterior electrodes, we can see that a broadband peakmax pattern develops in posterior channels shortly after loss of consciousness, at the same time as frontal troughmax.That is, even before frontal channels transition into peakmax dynamics, posterior channels already show broadband coupling to the slow-wave.In all subjects, the onset of posterior peakmax dynamics occurred after LOC and before the onset of frontal peakmax dynamics (See Supplementary Figures S2-S10).This result is further shown in the topographic maps of the 8-16 Hz cross-frequency coupling by level: shortly after loss of consciousness, a collection of frontal electrodes still have troughmax coupling while electrodes over posterior electrodes show peakmax coupling.As the propofol level is increased, the frontal electrodes begin to participate in the peakmax coupling.This suggests that broadband silencing of cortical activity by the slow-wave begins first over posterior areas before extending frontally at higher doses of propofol.
It has been noted that many traditional phase-amplitude coupling metrics could produce spurious coupling in the presence of non-sinusoidal waveforms 17 .The correlation metric used here avoids estimating the phase of the slow wave, making our metric less prone to these errors (see Methods).We also visually inspected the time-domain waveforms that generated the effects we report to ensure that waveform features did not generate spurious coupling (for example, see Supplementary Figure S13).The raw voltage traces show that time periods with broadband peakmax coupling (Low Dose for posterior electrodes and High Dose for both frontal and posterior electrodes) do not have the kind of sharp transients that would result in spurious broadband coupling 17 .
For the rest of the analyses, we identified four levels of interest: (1) Baseline, before the administration of propofol; (2) Sedation, the level before the subject stopped responding; (3) Unconscious Low Dose, the first level after the subject stopped responding; and (4) Unconscious High Dose, the highest level of propofol that did not result in burst suppression, a coma-like state beyond what is required for unconsciousness.This allowed us to combine information across subjects, given that the timecourse of LOC and ROC were different for each subject.
Since the patterns of cross-frequency coupling across frequency varied by electrode, propofol level, and subject, we wanted to identify the coupling patterns that best summarized this activity.To do this, we used a non-centered principal component analysis (PCA) to decompose the cross-frequency coupling patterns into principal frequency modes (Fig. 2).The principal modes are orthogonal patterns across frequency that capture successively smaller proportions of the total energy in the matrix representing cross frequency coupling by frequency, electrode, propofol level, and subject (Supplementary Figure S11).
The first principal mode is positive for all frequencies, so it represents broadband effects in the cross-frequency coupling, and it captures 78% of the total energy.The second and third principal modes, which have narrow-band peaks in the alpha and beta frequency ranges, respectively, capture much less of the energy.Hence the entrainment of broadband cortical activity to the slow-wave is by far the biggest contributor to the cross-frequency coupling during propofol anesthesia.
This dominant first frequency mode is of particular interest, because it can capture the broadband peakmax pattern observed in the individual subject summaries (Fig. 1 and Supplementary Figures S2-S10).In particular, we would like to use the first principal mode to describe the spatial distribution of the broadband peakmax phenomenon among the cortical generators of the EEG during progressively higher doses of propofol.In order to do so, we performed the cross-frequency coupling analysis on the source-localized EEG signals, and projected the resulting cross-frequency coupling patterns across cortical location onto the first principal mode from the sensor-space analysis described above.Figure 3 shows the resulting projections, morphed to the Freesurferaverage surface 18 and averaged across subjects.On these surfaces, a positive projection (red) represents broadband peakmax coupling to the slow-wave., "x"): "troughmax" (left) and "peakmax" (right).For each, the broadband EEG voltage (black) is bandpass filtered into low-frequency activity (LFA, orange) and high frequency activity ( α , light green).The instantaneous amplitude of the high frequency signal is computed using the analytic signal ( α AMP , dark green).The cross-frequency coupling is quantified by the correlation between the low frequency activity and the instantaneous amplitude of the high frequency signal.On the left the correlation is negative (troughmax, blue), and on the right the correlation is positive (peakmax, red).(b) Summary of cross-frequency coupling for one subject (28 years old).Top Panel: the target effect site concentration of propofol administered over the course of the session.Vertical lines indicate Loss of Consciousness (LOC) and Return of Consciousness (ROC).Second Panel: the probability of response curves for the two types of auditory stimuli (Verbal and Clicks).LOC was defined when the probability of response fell below 5% and remained for at least 5 min, and ROC was defined when the probability of response went above 5% and remained there for at least 5 min.Third and Fourth Panels: the cross-frequency coupling between the LFA and a range of amplitude frequencies (y-axis), for every 30 s interval in the session, for a set of five frontal electrodes (third panel) and a set of six posterior electrodes (fourth panel) with scalp positions indicated on the right.Fifth Panel: Scalp distribution of the cross-frequency coupling between the LFA and the 8-16 Hz amplitude.
Consistent with the individual subject sensor-space results, broadband peakmax coupling to the slow-wave is present over posterior cortical areas during the Unconscious Low Dose condition.It strengthens and spreads to frontal cortical areas during the Unconscious High Dose condition.In Fig. 4, we can see that during the Unconscious Low Dose condition the parietal, temporal, and occipital lobes all have broadband peakmax coupling, but the frontal lobe has virtually no broadband coupling to the slow-wave.During the Unconscious High Dose condition, in contrast, the frontal lobe joins the other lobes in broadband peakmax coupling.
Together these results suggest that broadband peakmax coupling to the slow-wave is a major contributor to cross-frequency coupling dynamics during propofol anesthesia (Fig. 2).Furthermore, posterior cortical areas are the first to exhibit broadband peakmax coupling after loss of consciousness, followed by frontal areas at high doses of propofol (Figs. 3, 4).

Discussion
We propose that the broadband coupling to the slow-wave that we see in the EEG signal is a macro-scale indicator of regional cortical up-and down-states.Population spiking activity has a broadband effect on the neural power spectrum 15,16 , so the entrainment of population spiking activity to the slow-wave during cortical up-and down-states should cause broadband power to be strongest at a particular phase of the slow-wave.This idea is supported by micro-scale evidence in electrocorticography and local field potential data showing that up-states are associated with simultaneous increases in both population spiking activity and broadband power (< 50 Hz) 14 .While it is common to assume that increased slow power in the EEG always corresponds to up-and down-states in cortical firing, our data suggests that the slow-wave must reach a critical amplitude, which may vary by location, before it entrains the population spiking activity enough to produce broadband modulation.Posterior areas cross this threshold at lower doses of propofol than frontal areas.The idea of a critical threshold for the slow-wave is consistent with the slow-wave activity saturation hypothesis 19 , in which saturation of the slow-wave power accompanies an attenuation of cortical and thalamic BOLD fMRI response to sensory stimulation.
We hypothesize that the cortical up-and down-states themselves provide a mechanism for the disruption of cortical sensory processing.Classical studies have shown that a midpontine trigeminal brainstem transection produces desynchronized "wakeful" EEG patterns in which animals can exhibit visual tracking and ocular accommodation.However, in these same animals, these visual behaviors cease when slow oscillations develop, as they periodically did during these studies [20][21][22] .More recently, many studies have demonstrated how slow wave dynamics are associated with disrupted cortical sensory function 19,[23][24][25] .Under this interpretation, broadband peakmax dynamics signify that local cortical activity is being disrupted by up-and down-states on a scale large enough to be detected in the EEG.
The ability to monitor a patient's state of unconsciousness during general anesthesia has obvious practical benefits, but identifying principled methods for doing so remains an open question.We can make a distinction between unconscious states in which patients can be aroused to consciousness by external stimuli (arousable), versus those in which patients cannot be aroused (unarousable).Some recent data suggest indirectly that patients can be arousable when frontal alpha rhythms are strongest at the trough of the slow-wave (troughmax), but are unarousable when the frontal alpha rhythms are strongest at the peak of the slow-wave (peakmax) 7,11,12 .Here we show that when frontal alpha rhythms are coupled to the peak of the slow-wave (peakmax), broadband spectral power over both frontal and posterior cortices is also coupled to the slow-wave.Hence, unarousability appears Here we show the results of a non-centered PCA analysis on the cross-frequency coupling patterns (i.e., to identify frequency-dependent functions describing the modulation of the high frequency signal by the LFA).The analysis was performed over all sensors, subjects, and propofol levels of interest (baseline, sedation, unconscious low dose, and unconscious high dose).Left: The first three principal modes, representing the patterns across frequency that capture the greatest percentage of total energy in the coupling patterns.Right: the percentage of total energy captured by each principal mode, in order from greatest to least.The first principal mode captures 77% of the total energy of the coupling patterns.
to be associated with broadband coupling to the slow-wave over both frontal and posterior cortex.At lower anesthetic doses in which patients are unconscious but potentially arousable, we find that the posterior cortex engages in broadband peakmax dynamics while the prefrontal cortex participates in (alpha-band) troughmax dynamics.Thus, these different spatial patterns of cross-frequency coupling appear to indicate different states of unconsciousness.Knowledge of these states could be used by anesthesiologists to establish different arousable or unarousable states of unconsciousness appropriate to the clinical circumstances.
The "posterior hot zone" [1][2][3] and fronto-posterior [4][5][6] hypotheses propose starkly different roles for posterior and frontal circuits in mediating consciousness.Anesthesia introduces another dimension to this problem in that unconscious patients may or may not be arousable by external stimuli depending on the drug dose.Along this dimension, the prefrontal cortex appears to play a central role, in which stimulation of prefrontal cortex www.nature.com/scientificreports/can jump-start consciousness 10 .Our results would seem to harmonize these viewpoints, in that the disruption of frontal versus posterior cortices coincides with different states of unconsciousness.First, both the posterior hot zone and fronto-posterior hypotheses predict that disruption of posterior areas should degrade the contents of consciousness.In our results, posterior areas are the first to show broadband peakmax at low doses of propofol, when subjects are unconscious but may be arousable.This is consistent with the interpretation that posterior cortical regions mediate the contents of consciousness 2 , which are disrupted by alternating up-and down-states during unconsciousness 3,14 .Second, frontal areas develop broadband peakmax coupling at higher doses of propofol, in a state associated with unarousability.This is consistent with the interpretation that prefrontal cortex is required for behavioral arousal to consciousness 10 and that alternating up-and down-states in prefrontal cortex prevent this arousal.Our results suggest that these two consciousness-supporting subsystems are disrupted at different levels of anesthesia-induced unconsciousness, leading to two empirically distinguishable states of unconsciousness.A limitation of this work is that unconsciousness was determined using only an auditory response task; arousability of the subjects by stronger and possibly noxious stimuli was not assessed.Further work will be required to more firmly establish the connections between the brain states observed in this study and behavioral arousability.This future work could enhance the clinical utility of these results by identifying the precise neural signatures that can distinguish arousable from unarousable states.Based on the present results and previous work 10,12 , we predict that patients would not be arousable to consciousness in the global (frontal plus posterior) peakmax state, but could be arousable when posterior peakmax is present in the absence of frontal peakmax.
The main result of this study is that brain-wide cross-frequency coupling analysis captures reliable but spatially-varying changes in brain dynamics during the transition to unconsciousness and unarousability under propofol anesthesia.We interpret broadband coupling to the peak of the slow-wave to be a macro-scale indicator of micro-scale cortical up-and down-states that serve to disrupt cortical function.This functional disruption encompasses both frontal and posterior cortical areas at different states of unconsciousness.Our results suggest that unconsciousness under anesthesia is not a singular phenomenon but rather involves several distinct shifts in brain state that disrupt posterior cortical processing of the "contents" of consciousness distinct from prefrontal arousal and awareness of those contents.We believe this distinction can help answer the fundamental clinical question of whether or not a patient can be aroused to consciousness during surgery.Moreover, this work identifies EEG signatures that could be used practically to identify and monitor those states of unconsciousness.From a neuroscientific perspective, our results bridge two historically separate perspectives in consciousness research, one focused on the neural correlates of consciousness and the other studying levels of consciousness and unconsciousness, by identifying distinct cortical substrates that could mediate these phenomena.In addition, our work suggests that other forms of anesthesia-or drug-induced altered consciousness could be studied under a similar paradigm to better understand the interplay of the level and contents of consciousness.

Experimental design
This study analyzes data that has been previously described in 11 and 13 .All procedures were approved by the Human Research Committee (IRB) at the Massachusetts General Hospital.All research was performed in accordance with relevant guidelines and regulations, and informed consent was obtained from all participants.
After at least 14 min of eyes-closed rest, subjects were administered propofol anesthesia, increasing propofol dosage every 14 min to target 1, 2, 3, 4, or 5 µg mL −1 effect-site concentration.During emergence, propofol dose was decreased in three 14-min steps targeting effect-site concentrations of 0.5 µg mL −1 , 1.0 µg mL −1 , and 1.5 µg mL −1 less than the dose at which the subject stopped responding to the auditory cues.Electroencepha- logram (EEG) recordings were collected with a high density (64 channel) BrainVision MRI Plus system (Brain Products) with a sampling rate of 5,000 samples per second and the electrode positions were digitized (Polhemus FASTRACK 3D).
Throughout the study, subjects performed a task in which they responded to an auditory click or word with a finger-press of a button.Loss of consciousness (LOC) was defined when a state-space model estimated a less than 5% probability of a correct response and remained so for at least 5 min 11 .Similarly, return of consciousness (ROC) was defined as the first time during emergence when the state-space model reached a 5% probability of a correct response and remained so for at least 5 min.

EEG preprocessing
Channels with severe, persistent artifacts were identified by visual inspection and removed from the analysis.We also removed electrodes that showed signs of being electrically bridged according to the eBridge software package 26 .
The cross-frequency coupling analysis described below uses versions of the signal that have been bandpass filtered in the slow band (0.1-4 Hz) and a set of higher bands between 4 and 50 Hz, each 2 Hz in bandwidth with a 1 Hz transition band.We filtered the data for the entire session into these bands using a zero-phase FIR filter and downsampled to 200 samples per second.We chose the upper limit of the higher bands, 50 Hz, in order to avoid line noise and capture the frequencies with the largest EEG power 27 .
For sensor-level analyses, we applied a Laplacian reference based on first nearest neighbors after filtering 28 .
Because of the stepped structure of the propofol dosage, subjects had approximately 12 min at each dosage with constant predicted effect-site concentration.We call these periods "levels", and defined four levels of interest: Baseline, Sedation, Unconscious Low Dose, and Unconscious High Dose.Baseline was defined as the period before the onset of propofol, Sedation was defined as the level prior to LOC, Unconscious Low Dose was defined as the first full level after LOC, and Unconscious High Dose was defined as the highest level of propofol or the level before burst suppression (for subjects who entered burst suppression).For analyses performed on these levels, we selected 10 artifact-free epochs 30 s long from each level for each subject.Two of the subjects stopped responding during the first level of propofol, so they did not have a sedation period.As a result, figures summarizing across subjects represent 10 subjects for Baseline, Low Dose, and High Dose, and 8 subjects for Sedation.

Cross-frequency coupling analysis
Under propofol, high frequency amplitude is much more likely to couple to the peak or the trough of the slowwave than to the rising or falling phases 11 (Supplementary Figure S1).Hence we summarize cross-frequency coupling using a metric that is either positive (peakmax: high frequency amplitude is higher during the peak of the slow-wave) or negative (troughmax: high frequency amplitude is higher during the trough of the slow-wave).In particular, we use the Pearson correlation coefficient between the instantaneous amplitude of the band-passed high frequency signal and the band-passed slow voltage (0.1-4 Hz) 29 : where V is a column vector representing a timeseries of the slow voltage and A is a vector representing the corresponding timeseries of the instantaneous amplitude of the high frequency band.The instantaneous amplitude was computed using the magnitude of the analytic signal of the bandpass filtered signal.Before computing the correlation, the instantaneous amplitude for each 30 s epoch was centered by subtracting the mean.V was not explicitly centered because its expected value is zero.
In order to find the average correlation across space (electrodes, source locations) and/or time (multiple 30 s epochs), the dimensions to be averaged are stacked vertically in the V and A column vectors in Equation 1. (1) Since the correlation represents the covariance between the signals (in the numerator, V T A ) normalized by the standard deviations of each signal (in the denominator, √ V T V and √ A T A ), stacking the signals vertically has the effect of estimating the covariance and standard deviations separately before performing the normalization.
Note that this analysis does not involve estimating the phase of the slow-wave, a typical step in phase-amplitude coupling analysis 29 .We chose not to estimate phase for three reasons.First, phase-amplitude coupling in propofol clusters around zero (slow-wave peak, when the signal is positive) and π (slow-wave trough, where the signal is negative) 11 (see also Supplementary Figure S1).Second, the slow-wave is thought to reflect up states and down states, regardless of whether it is well-described by a sinusoid.Early in propofol-induced unconsciousness, the slow-wave is often quite irregular in waveform shape and frequency.Since phase estimation can be sensitive to waveform shape 17 , it was important to choose a metric that would be robust to these effects.Hence, we did not estimate the phase of the slow wave, but rather filtered the signal in a relatively wide band, 0.1-4 Hz, in order to capture more of the non-sinusoidal features of the waveform.Under propofol anesthesia, there is no narrow-band delta frequency rhythmic activity (1-4 Hz), so the contribution of this frequency range to the band-passed signal was non-rhythmic, affecting the waveform shape.By using correlation to describe the relationship of the high frequency amplitude to the non-sinusoidal low frequency waveform, we are therefore better able to capture the relationship between high frequency activity and the up and down fluctuations that are apparent in the broadband signal (Fig. 1, Panel A).Finally, reducing the problem to one dimension (positive and negative coupling) simplified the analysis of a single location and frequency, allowing us to expand the analysis to other frequencies and locations.
Panel A in Fig. 1 shows the computation for a frontal electrode, computing the correlation between the slow voltage (low-frequency activity, LFA) and the alpha band ( α , 8-16 Hz) during two time windows.For the time window on the left, the alpha amplitude is higher when the slow voltage is negative, leading to a negative correlation (troughmax).For the time window on the right, the alpha amplitude is higher when the slow voltage is positive, leading to a positive correlation (peakmax).We can perform this analysis for all 30 s epochs in the session and for a range of amplitude frequencies (2 Hz bands between 4 and 50 Hz): this results in a modulogram showing the cross-frequency coupling between each frequency and the slow-wave across time in the session for the selected electrode.Similarly, we can compute the correlations for each electrode within propofol levels, and show the spatial distribution of cross-frequency coupling between the alpha band and the slow-wave as a function of propofol level.
We performed the cross-frequency coupling analysis across amplitude bands, time, electrodes, and subjects.We characterize time in two ways.The first was a session-based analysis, in which the coupling was assessed for every 30 s epoch in the session separately, resulting in modulograms such as the frontal and posterior summaries in Panel B of Fig. 1: these modulograms represent the average correlations across several electrodes, indicated in the scalp maps to the right.The second way to characterize time was based on propofol level, in which the coupling was assessed over 10 artifact-free 30 s epochs at each propofol level.Hence, each coupling value represents the correlation over 300 s of data, and we have one such value for each electrode, level, and subject.This level-based analysis was used for the topoplots in Panel B of Fig. 1 and for the Principal Component Analysis described below.

Source localization
Source localization was performed by minimum norm estimation 30 using the MNE-python toolbox 31 , www.martinos.org/mne/].The noise covariance was assumed to be diagonal, estimated from the EEG using the spectral power above 65 Hz.The source space was an ico-3 decimation of the Freesurfer reconstruction of the cortical surface and the sources were constrained to be perpendicular to the cortical surface.The forward model was computed using a 3-layer boundary element model, eliminating any sources within 5mm of the inner skull surface.
To perform the cross-frequency coupling analysis in source space, the band-passed signals in the slow band and each 2 Hz band between 4 and 50 Hz were source localized using the methods described above.For the cortical surface representations (Fig. 3), the cross-frequency coupling analysis described above was applied to each source location for each subject, for the four levels of interest.The average cross-frequency coupling within lobes (Fig. 4) used the same technique, but the signals for all of the sources in a lobe were combined: the 10 epochs (for each level) for all of the sources in the lobe were stacked vertically in the V and A column vectors in Equation 1.

Non-centered principal component analysis and projection into source space
To identify patterns of coupling across frequency that explained most of the observed features, we performed a non-centered principal component analysis 32 on the sensor-space coupling patterns across propofol levels and subjects.The level-based coupling results for the four levels of interest (baseline, sedation, unconscious low dose, and unconscious high dose) were combined into an aggregate matrix A(f, i) where f indexes the amplitude band and i indexes the propofol level, electrode, and subject (i.e. the rows contain all of the levels, electrodes, and subjects stacked).Supplementary Figure S11 shows the cross-frequency coupling results for two electrodes during the four levels of interest, for all subjects: the goal of the analysis was to decompose these patterns into components that capture their major features across frequencies.Principal component analysis involves a singular value decomposition: In the decomposition, U is a matrix whose columns are the principal modes.The first column is the first principal mode, meaning the pattern across frequencies that captures the most energy in the A matrix.The second column is the second principal mode, which captures the most energy in the data after the first mode has been

Figure 1 .
Figure 1.Peakmax coupling is broadband, and it begins after loss of consciousness on posterior electrodes before frontal electrodes.(a) Quantification of cross-frequency coupling.Two example traces from a frontal electrode (middle, "x"): "troughmax" (left) and "peakmax" (right).For each, the broadband EEG voltage (black) is bandpass filtered into low-frequency activity (LFA, orange) and high frequency activity ( α , light green).The instantaneous amplitude of the high frequency signal is computed using the analytic signal ( α AMP , dark green).The cross-frequency coupling is quantified by the correlation between the low frequency activity and the instantaneous amplitude of the high frequency signal.On the left the correlation is negative (troughmax, blue), and on the right the correlation is positive (peakmax, red).(b) Summary of cross-frequency coupling for one subject (28 years old).Top Panel: the target effect site concentration of propofol administered over the course of the session.Vertical lines indicate Loss of Consciousness (LOC) and Return of Consciousness (ROC).Second Panel: the probability of response curves for the two types of auditory stimuli (Verbal and Clicks).LOC was defined when the probability of response fell below 5% and remained for at least 5 min, and ROC was defined when the probability of response went above 5% and remained there for at least 5 min.Third and Fourth Panels: the cross-frequency coupling between the LFA and a range of amplitude frequencies (y-axis), for every 30 s interval in the session, for a set of five frontal electrodes (third panel) and a set of six posterior electrodes (fourth panel) with scalp positions indicated on the right.Fifth Panel: Scalp distribution of the cross-frequency coupling between the LFA and the 8-16 Hz amplitude.

Figure 2 .
Figure 2. Broadband effects dominate the cross-frequency coupling patterns across subjects, electrodes, and anesthetic drug levels.Here we show the results of a non-centered PCA analysis on the cross-frequency coupling patterns (i.e., to identify frequency-dependent functions describing the modulation of the high frequency signal by the LFA).The analysis was performed over all sensors, subjects, and propofol levels of interest (baseline, sedation, unconscious low dose, and unconscious high dose).Left: The first three principal modes, representing the patterns across frequency that capture the greatest percentage of total energy in the coupling patterns.Right: the percentage of total energy captured by each principal mode, in order from greatest to least.The first principal mode captures 77% of the total energy of the coupling patterns.

Figure 3 .
Figure 3. Broadband peakmax occurs first over posterior cortical regions during unconsciousness at lower doses of propofol, and extends to encompass to frontal cortical regions at higher propofol doses.These figures show the projection of the source-space cross-frequency coupling patterns onto the first principal mode for the four levels of interest.Positive values reflect broadband cross-frequency coupling to the peak of the slow-wave, while negative values represent broadband coupling to the trough of the slow-wave.We present the average coupling across all subjects.

Figure 4 .
Figure 4.During unconsciousness at low propofol doses, the broadband peakmax effect occurs over the occipital, parietal, and temporal lobes, but not the frontal lobe.Here we show cross-frequency coupling patterns estimated across each cortical lobe using source localization, projected onto the first principal mode and averaged across subjects.Error bars represent 95% confidence intervals, estimated using a bootstrap across subjects.