Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Articles, Behavioral/Systems/Cognitive

Spatiotemporal Reconfiguration of Large-Scale Brain Functional Networks during Propofol-Induced Loss of Consciousness

Manuel S. Schröter, Victor I. Spoormaker, Anna Schorer, Afra Wohlschläger, Michael Czisch, Eberhard F. Kochs, Claus Zimmer, Bernhard Hemmer, Gerhard Schneider, Denis Jordan and Rüdiger Ilg
Journal of Neuroscience 12 September 2012, 32 (37) 12832-12840; https://doi.org/10.1523/JNEUROSCI.6046-11.2012
Manuel S. Schröter
1Max Planck Institute of Psychiatry, 80804 Munich, Germany,
2Department of Neurology,
3Department of Anesthesiology, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Victor I. Spoormaker
1Max Planck Institute of Psychiatry, 80804 Munich, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anna Schorer
2Department of Neurology,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Afra Wohlschläger
4Department of Neuroradiology, Klinikum rechts der Isar, Technische Universität München, 81675 Munich, Germany, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Czisch
1Max Planck Institute of Psychiatry, 80804 Munich, Germany,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eberhard F. Kochs
3Department of Anesthesiology, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Claus Zimmer
4Department of Neuroradiology, Klinikum rechts der Isar, Technische Universität München, 81675 Munich, Germany, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bernhard Hemmer
2Department of Neurology,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gerhard Schneider
3Department of Anesthesiology, and
5Department of Anesthesiology, Witten/Herdecke University, HELIOS Clinic Wuppertal, 42283 Wuppertal, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Denis Jordan
3Department of Anesthesiology, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rüdiger Ilg
2Department of Neurology,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Applying graph theoretical analysis of spontaneous BOLD fluctuations in functional magnetic resonance imaging (fMRI), we investigated whole-brain functional connectivity of 11 healthy volunteers during wakefulness and propofol-induced loss of consciousness (PI-LOC). After extraction of regional fMRI time series from 110 cortical and subcortical regions, we applied a maximum overlap discrete wavelet transformation and investigated changes in the brain's intrinsic spatiotemporal organization. During PI-LOC, we observed a breakdown of subcortico-cortical and corticocortical connectivity. Decrease of connectivity was pronounced in thalamocortical connections, whereas no changes were found for connectivity within primary sensory cortices. Graph theoretical analyses revealed significant changes in the degree distribution and local organization metrics of brain functional networks during PI-LOC: compared with a random network, normalized clustering was significantly increased, as was small-worldness. Furthermore we observed a profound decline in long-range connections and a reduction in whole-brain spatiotemporal integration, supporting a topological reconfiguration during PI-LOC. Our findings shed light on the functional significance of intrinsic brain activity as measured by spontaneous BOLD signal fluctuations and help to understand propofol-induced loss of consciousness.

Introduction

The majority of the brain's energy consumption is associated with metabolism inherent to intrinsic neuronal activity (Ames, 2000). Intrinsic brain activity, as reflected in spontaneous slow fluctuations in BOLD signals in functional magnetic resonance imaging (fMRI), has been proposed to provide important clues on brain organization (Fox and Raichle, 2007) and mapping of its spatiotemporal correlation structure seems a promising approach to assess basic properties of brain function (Biswal et al., 2010). Patterns of coherent signal fluctuations in fMRI have been associated with large-scale functional brain networks (Biswal et al., 1995; Cordes et al., 2001; Greicius et al., 2003; Fox et al., 2006; Seeley et al., 2007) and suggested to contribute to variability in human behavioral task-performance and perception (Hampson et al., 2006; Fox et al., 2007).

Recent research has probed the functional significance of intrinsic brain activity to explain altered states of consciousness. Positron-emission tomography (PET) studies reported decreases in global brain metabolism throughout non-rapid eye movement (NREM) sleep (Braun et al., 1997) and in pharmacologically modulated (Alkire et al., 1995; Fiset et al., 1999) and pathologically altered states of consciousness (Laureys et al., 2004). For hypnotic anesthetics, such as propofol, hypometabolism during loss of consciousness was reported for parts of the ascending reticular activation system, in particular thalamus, the association cortices, such as posterior cingulate cortex (PCC), and frontoparietal cortices (Kaisti et al., 2003). Since findings overlap largely with metabolic decreases during NREM sleep, anesthesia and sleep have been suggested to share some common functional circuits to induce loss of consciousness (Lydic and Baghdoyan, 2005).

But how does the organization of large-scale functional networks change during propofol-induced loss of consciousness (PI-LOC)? While early fMRI studies of BOLD signal fluctuations across levels of consciousness reported persistent intrinsic activity and relatively preserved spatial organization during light sleep and sedation (Greicius et al., 2008; Horovitz et al., 2008; Larson-Prior et al., 2009), recent studies point to a state-dependent modulation: during propofol-anesthesia, a linear association between functional connectivity of frontoparietal networks and levels of consciousness (Boveroux et al., 2010), a reorganization of key nodes of the default mode network (Stamatakis et al., 2010), and decreased subcortical-cortical connectivity have been reported (Mhuircheartaigh et al., 2010). State-dependent changes in functional connectivity have also been observed during deep NREM sleep (Horovitz et al., 2009; Sämann et al., 2011). How these findings relate to recent models of spontaneous brain activity and conscious processing (Dehaene and Changeux, 2005) or information integration (Tononi, 2004) has yet to be tested systematically.

Assuming potential similarities between some physiologic characteristics of deep NREM sleep and anesthesia-induced loss of consciousness (Lydic and Baghdoyan, 2005) and motivated by previous observations of state-dependent connectivity changes during NREM sleep (Spoormaker et al., 2010), we applied graph theory to investigate changes in topological properties of whole-brain functional networks. During PI-LOC, we hypothesized to find reduced functional connectivity strength, particularly of subcortical-cortical connections, and a reduced probability of high-degree hubs. Furthermore, we expected differential changes in brain topology, such as a decrease in spatiotemporal integration and a more insulated pattern of activity, possibly related to a reduction in long-range connections.

Materials and Methods

Subjects.

Fifteen male subjects participated in the study. The study protocol was in accordance with the Declaration of Helsinki and approved by the local ethics committee of the medical department of the Technische Universität München, Germany. All subjects were right-handed, drug-free, and had no history of neurological or psychiatric disorders. Subjects had to give written informed consent and were reimbursed for participation. Individuals with excessive head movements during fMRI scanning (>2 mm between two subsequent EPI images) were excluded, leaving a sample of 11 subjects (mean age ± SD: 25.8 ± 3.0 years). In the present analyses, we focus on changes from wakefulness to PI-LOC.

Sedation protocol.

Propofol [2,6 di-isopropylphenol; (James and Glen, 1980)] is an anesthetic agent that is thought not to interfere with cerebral blood flow response and flow-metabolism coupling (Johnston et al., 2003; Veselis et al., 2005). Propofol acts mainly on GABAA receptors, potentiating GABAergic interneurons in the cortex and enhancing neuronal inhibition (Antkowiak, 1999; Franks, 2008; Brown et al., 2011). The effect depends mainly on the dosage and speed of administration. During simultaneous electroencephalography (EEG) and fMRI, propofol was administered in the scanner using a target controlled infusion (TCI) pump (Open TCI, Space Infusion System; Braun Medical) to obtain constant effect-site concentrations, as estimated by the pharmacokinetic model of Marsh et al. (1991). TCI concentrations were increased in 0.4 μg/ml steps beginning at 1.2 μg/ml until volunteers stopped responding to the verbal command “squeeze my hand” accompanied by a firm handshake [PI-LOC corresponding to a Ramsay sedation scale of 5–6 (Ramsay et al., 1974), i.e., sluggish or no response to tactile or auditory stimulation]. TCI concentration at this point was maintained and kept stable for another 10 min to ensure equilibration of the estimated effect-site concentration. Having reached this equilibrium, simultaneous EEG and fMRI measurements were performed subsequently for ∼10 min. After that, a phase of sedation at 50% of the PI-LOC concentration was recorded; however, movement artifacts prevented a reliable analysis of these fMRI data. Throughout the experiment, standard monitoring parameters, such as breathing frequency, blood pressure, cardiac rhythm, and pulse oxymentry, were measured. Oxygen administration (via nose tube) and expiratory CO2 concentrations were continuously monitored.

EEG analysis.

Nonlinear permutation entropy (PeEn) of EEG electrodes has been established to reliably separate normal wakefulness from PI-LOC (Jordan et al., 2008). PeEn is a measure of the irregularity of signals and is based on a comparison of the neighboring order of signal values (Bandt and Pompe, 2002). It has been shown that PeEn is largely unaffected by signal disturbances and can be used to analyze time series generated by high-dimensional dynamic systems with low stationarity. In the present investigation, we used PeEn as an additional measure to indicate levels of vigilance and to verify the effect of state on brain functional networks. EEG recordings were performed using a 64-channel electrode cap with nonmagnetic Ag/AgCl electrodes (Easycap) and two 32-channel nonmagnetic battery-operated EEG amplifiers (BrainAmp MR; Brain Products) designed for MRI in-scanner measurement. One EEG channel was dedicated to register the electrocardiogram (ECG) and placed over the subject's chest. The resolution of the EEG acquisition system was 100 and 500 nV for ECG and signals were recorded with a sampling rate of 5 kHz (BrainVision Recorder; Brain Products). EEG signal preprocessing was performed with BrainVision Analyzer 2 (Brain Products). Gradient artifacts produced by fMRI were averaged over a sliding window and subtracted from the EEG signals. Cardioballistic artifacts caused by cardiomechanic electrode induction were removed using a template-detection method. Templates were based on detected local maxima (R-peak) of the ECG and subtracted from the EEG using sliding windows of 21 epochs. Basic artifact rejection (EEG with amplitudes exceeding 250 μV), average reference, and independent component analysis for blind source separation of noncortical signal components were computed. Finally, PeEn of 18 frontal EEG electrodes (Fp1–FC6, average reference, total bandwidth 0.5–30 Hz) was calculated.

fMRI acquisition and preprocessing.

Resting-state fMRI recordings were acquired at a 3T MR system (Achieva Quasar Dual; Philips) using an eight-channel phased-array head coil and a gradient echo EPI sequence (TR, 1838 ms; TE, 30 ms; flip angle, 75°; 72 × 72 matrix; FOV, 220 × 220 mm; 32 slices oriented along AC-PC; slice thickness, 3 mm; 1 mm interslice gap; scan duration 10 min). The following fMRI preprocessing steps were performed using SPM (version 5; www.fil.ion.ucl.ac.uk/spm): slice time correction to account for interleaved slice acquisition, realignment to the first image to correct for small head movements, and spatial normalization of all images to a standard EPI template in Montreal Neurological Institute (MNI) standard space. Only scans with head movements <2 mm between two subsequent volumes were accepted, resulting in 200 scans per subject for the final connectivity analysis. Finally, images were resampled to a voxel resolution of 2 × 2 × 2 mm3 using a fifth-degree spline interpolation. To remove unspecific signal variation from the time series, the following regressors were defined in a multiple regression model: six head movement parameters derived from the realignment step (1–6) and global signals within white matter (WM) and ventricular CSF (7–8). Time series for WM and CSF were calculated using 3 mm radius spheres in deep WM and CSF (Chang and Glover, 2009). No global signal removal was applied, since it was observed to suppress meaningful neural activity (Zhang et al., 2008; Schölvinck et al., 2010) and to artificially increase anticorrelations (Chang and Glover, 2009; Weissenbacher et al., 2009). Resultant residual images only contained signals not explained by regressors 1–8 and were used for extraction of regional time courses. To ensure that changes in functional connectivity were not biased by differences in head motion, we compared mean movements of experimental conditions using the root-mean-square of translational parameters (Van Dijk et al., 2012).

Brain parcellation and time course extraction.

Multiple brain atlases were used to extract time series for our connectivity analyses. The main analysis was performed on the Harvard–Oxford Probabilistic Atlas (threshold at 25%), as implemented in FSL [FMRIB Software Library, www.fmrib.ox.ac.uk/fsl (Smith et al., 2004)], which subdivides the brain into 110 cortical and subcortical regions. Regional time series were extracted using the average signal of all voxels inside each region using the Marsbar toolbox (Brett et al., 2002). Time series then underwent a maximum overlap discrete wavelet transform (MODWT) by applying the WMTSA wavelet toolkit (www.atmos.washington.edu/∼wmtsa/) with a Daubechies wavelet (Percival and Walden, 2000). This resulted in four frequency scales [0.13–0.27 Hz (scale 1), 0.07–0.13 Hz (scale 2), 0.03–0.07 Hz (scale 3), and 0.01–0.03 Hz (scale 4)] with frequency-band-specific correlation matrices. To increase comparability with previous studies (Achard et al., 2006; Spoormaker et al., 2010) and to reduce confounds through nonneuronal physiological variations (Birn et al., 2006; Shmueli et al., 2007), we focus our main analyses on frequency band 0.03–0.07 Hz (scale 3). To control for possible effects of brain parcellation (Zalesky et al., 2010), we additionally performed graph analyses on the widely used automatic anatomical labeling (AAL) template (Tzourio-Mazoyer et al., 2002) and a recently proposed functional brain parcellation by Dosenbach et al. (2010), in which all 160 cerebral spheres have the exact same size and do not overlap. To allow comparisons of our analysis with classical resting-state fMRI approaches, we further analyzed the frequency band 0.01–0.10 Hz using a third-order Butterworth filter.

Graph construction.

Graph theoretical analyses were run on binary adjacency matrices A. A is gained by thresholding the weighted matrix M at an individual correlation threshold R, ensuring statistics on graphs with the same number of edges and nodes. Statistical analysis on binary graphs was restricted to a range of thresholds including connection probabilities of 36–50% (1% increments, 15 values in total). This procedure follows the rationale that some networks may become disconnected at lower thresholds (for the present graphs, this was the case for thresholds <36% connection probability) and that network configurations become increasingly unspecific and similar to those of random graphs at very high thresholds [e.g., a connection probability >50% (Humphries and Gurney, 2008)]. Comparable approaches have been applied in previous studies (Achard and Bullmore, 2007; Bassett et al., 2008; Liu et al., 2008; van den Heuvel et al., 2009; Lynall et al., 2010). Network characterizations based on binary networks gained by thresholding approaches are likely to neglect at least some information, such as the role of negative correlations. However, since the biological relevance of anticorrelations in fMRI networks is still under debate (Chang and Glover, 2009; Fox et al., 2009) and relies to some degree on the choice of fMRI preprocessing strategy (Murphy et al., 2009; Weissenbacher et al., 2009; Schwarz and McGonigle, 2011), we ground our main functional network analysis on graphs that are fully connected and maximally sparse (i.e., in the described range of connection probabilities).

Graph theoretical properties.

The following functional connectivity properties were estimated on the weighted matrices M of wavelet correlations r between the time courses of nodes i and j: regional connectivity strength, diversity, and integration. The regional connectivity strength S(i) is the average of all correlations between the ith node and all other nodes in the brain network. Nodal connection diversity D(i) is defined as the column variance over the connectivity strength S(i) of each node i (Campbell et al., 1986; Bassett et al., 2011). Interareal integration I of the functional brain network was estimated by a principal component analysis using the ratio of the first eigenvalue to the sum of all eigenvalues (Lynall et al., 2010). Functional network metrics degree k, clustering coefficient C, path length L, clustering ratio γ, path length ratio λ, and small-worldness σ were estimated for the range of connection probabilities 36–50% on adjacency matrices A: the nodal degree k of a node is the number of edges that link a node to the rest of the network; it is calculated on the adjacency matrix A. The regional clustering coefficient C is a measure of functional segregation, which is calculated by the ratio of the number of edges that exist between direct neighbors of a node and the maximum number of possible connections between these neighbors (Watts and Strogatz, 1998); only nodes with degree ≥2 were considered. Compared with complex networks, random networks are expected to have a lower average clustering. Characteristic path length L of a network is calculated as the harmonic mean of all geodesic distances di,j between node i and j; it measures the average minimum number of steps that must be traversed to get from one node to another (Watts and Strogatz, 1998; Hayasaka and Laurienti, 2010). Path length has been associated to a network's global efficiency (Latora and Marchiori, 2001; Achard and Bullmore, 2007). Finally, small-worldness σ is given by the clustering ratio γ, defined as C/Crandom, and the path length ratio λ, defined as L/Lrandom. A network is assumed to possess small-world organization if γ > 1 and λ ∼ 1, resulting in σ values >1. For λ and γ, average Crandom and Lrandom were computed by applying a Markov chain algorithm to randomize all edges of the original networks while preserving their degree distribution. For each individual brain graph and all thresholds in the chosen range of connection probabilities, 10 randomly rewired matrices were generated (Maslov and Sneppen, 2002; Stam et al., 2007). Graph metrics were calculated with the Brain Connectivity Toolbox [www.brain-connectivity-toolbox.net (Rubinov and Sporns, 2010)]. Small-world organization seems an intuitively attractive model to describe brain networks since it enables segregated information processing in specialized modules as well as globally coordinated integration across the brain (Bassett and Bullmore, 2006). Small-worldness was found on different scales and for different modes of neuronal networks (Bullmore and Sporns, 2009). Dense local clustering of functionally related cortices and relatively sparse global interconnectedness, i.e., a short average path length, is present in functional and structural brain networks (Salvador et al., 2005; Achard et al., 2006; Hagmann et al., 2008). Modeling work supported the advantages of this network topology for efficient information processing (Latora and Marchiori, 2001), learning (Simard et al., 2005), and dynamic adaptivity and emergence of complexity (Sporns and Zwi, 2004).

Subsystem connectivity analysis.

We compared subsystem-specific connectivity changes across wakefulness and PI-LOC by allocating all 110 regions of the Harvard–Oxford Atlas into five major subsystems as proposed by Mesulam (2000): association, limbic, paralimbic, primary sensory, and subcortical regions. Differences in intra-subsystem and inter-subsystem connectivity were computed on the regional MODWT correlation coefficients and grouped into 15 categories containing all possible subsystem interactions (e.g., association–association, association–limbic, association–paralimbic, association–primary, association–subcortical, etc.).

Topography analysis.

We computed differential connection distance during PI-LOC by splitting all MODWT correlation coefficients into two groups according to their Euclidean distance in the anatomical brain graph (MNI space). The length of a connection was estimated by subtracting the respective x, y, and z coordinates of two nodes i and j from each other and by square root transforming the sum of the squared differential scores. The average connection distance of all brain regions in the graph, as measured from their respective center of mass, was 78 mm. Long distant connections were defined as correlations with >78 mm length; local connections were defined as having ≤78 mm length.

Statistical analysis.

Statistical analysis focused on differences in functional connectivity and topology between wakefulness and PI-LOC. Analyses on graph measures across the range of connection probabilities (36–50%), weighted graph measures, and subsystem connectivity were performed using paired t tests; a false discovery rate (FDR) was applied to correct for multiple comparisons (q set at 0.05). Paired t tests were also calculated on the differential scores of topographical distance measures. Changes in the degree distribution were calculated on mean graphs of the range of connection probabilities using χ2 tests; distributions were split into high-degree nodes (≥45 connections) and low-degree nodes (≤45 connections). MATLAB (version 2009a; MathWorks) and SPSS software version 18.0 were used for all statistical and computational operations. Effects of propofol on PeEn were tested using 95% percentile bootstrap confidence intervals (threshold p < 0.05) based on the area under the receiver operating characteristic curve (AUC) by comparing results in clustered frontal electrodes according to the 10-20 scheme (Jordan et al., 2010).

Visualization.

Changes in connectivity strength were rendered on cortical surface with Caret v5.61 software [http://brainvis.wustl.edu/wiki (Van Essen et al., 2001)]. Cortical coloring is derived from the respective group values for each region of the Harvard–Oxford Atlas and mapped to the Caret PALS-B12 average template. Additionally, BrainNet Viewer v.1.0 (www.nitrc.org/projects/bnv) was used to visualize changes in connection density.

Results

Physiological measures

During PI-LOC, subjects showed a significant drop in systolic blood pressure (Table 1). We found no significant effect on heart rate or respiration by propofol. There was a decrease in oxyhemoglobin saturation during PI-LOC, but no statistically significant change for CO2 saturation compared with baseline. In the selected 11 subjects, no significant difference in individual's mean head movement was observed between baseline and PI-LOC. Due to problems in EEG signal recording, one subject had to be excluded from the PeEn analysis. PeEn of frontal EEG electrodes of remaining participants robustly separated wakefulness from PI-LOC: an AUC of 0.94 (0.76–1.00) of PeEn (mean values of 18 EEG channels) confirmed significant changes in EEG from wakefulness to PI-LOC, where maximum of sensitivity and specificity was found at a threshold of 3.45.

View this table:
  • View inline
  • View popup
Table 1.

Physiological variables

Effects of PI-LOC on brain functional connectivity

Paired t tests yielded significant decreases in average connectivity strength S between wakefulness and PI-LOC across all scales of the wavelet decomposition, including MODWT scale 3 (Table 2). Significant reduction of regional connectivity strength was found in several regions of occipital, temporal, and parietal lobes (bilateral inferior, middle and superior temporal gyri, occipital, lingual and fusiform gyri, and superior parietal lobule) and in bilateral putamen and thalamus (Fig. 1). We further observed a significant reduction in whole-brain integration I (t(10) = −2.36, p = 0.039), i.e., the percentage of variance accounted for by the first principal component decreased from 36% in the awake resting condition to 27% in PI-LOC (Table 2).

View this table:
  • View inline
  • View popup
Table 2.

Changes in functional connectivity and topological metrics

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Effects of PI-LOC on functional connectivity strength. A, Average regional connectivity strength S during wakefulness and after loss of consciousness (PI-LOC); yellow to red colors indicate increasing connectivity strength. B, Lower figure shows significant p values of between-condition decreases in connectivity strength S of left and right hemisphere (LH and RH); light to dark blue colors indicate increasing significance. The dotted line depicts the FDR threshold (q = 0.05; see Materials and Methods).

Effects of PI-LOC on brain topological properties

In both conditions, the brain graphs showed a classical small-world topology (σ > 1). The following differences in the network's global/local organization were observed: Using one-sided paired t tests, we found a significant decrease of clustering coefficient C during PI-LOC (t(10) = −2.19, p = 0.026; Fig. 2; Table 2). Since clustering coefficient Crandom of comparable random graphs decreased relatively more strongly (t(10) = −2.48, p = 0.016), we found a significant increase in the clustering ratio γ (t(10) = 1.86, p = 0.046). There was no significant effect on path length L and path length ratio λ for the applied range of connection probabilities. Accordingly, the increase in small-worldness σ (t(10) = 1.90, p = 0.043) was driven by differences in the clustering ratio γ. Furthermore, we found significant differences in mean degree (Fig. 3A) and degree distribution of average graphs across the range of connection probabilities (χ2(9) = 44.10, p < 0.001), as well as a significant decrease in the proportion of high-degree nodes during PI-LOC [χ2(1) = 17.50, p < 0.001; degree distribution split into high-degree nodes (≥45 connections) and low-degree nodes (≤45 connections)]. Findings were robust across different brain parcellations (AAL template, regions of interest derived from Dosenbach et al., 2010), frequency bands of the MODWT, and an analysis with 0.01–0.1 Hz bandpass-filtered time series (data not shown).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Effects of PI-LOC on brain topological properties. Figure shows topological changes of brain functional networks across wakefulness (black) and PI-LOC (gray). A–C, Summary statistics for clustering (C, Crandom, γ = C/Crandom) (A), path length (L, Lrandom, λ = L/Lrandom) (B), and small-worldness (σ = γ/λ) (C) were performed on MODWT scale 3 (0.03–0.07 Hz) throughout a range of 36–50% connection probability (Table 2). Error bars indicate 95% CI.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Changes in mean degree and short/long-range connections. A, Group average networks for wakefulness and loss of consciousness. During PI-LOC, brain functional networks decreased in mean degree; depicted at correlation threshold r = 0.3. B, Ratio of long- to short-range connections at PI-LOC compared with wakefulness (mean and SEM values of all long/short-range correlations). Both long- and short-range correlations decreased significantly during PI-LOC; * indicates significance of wakefulness versus PI-LOC for one-sample t test (p < 0.05, one-sided). The reduction in long-range connections was significantly stronger; * indicates significance of wakefulness versus PI-LOC for paired t test (p < 0.05).

Effects of PI-LOC on subsystem connectivity

Connectivity within and between functionally defined subsystems (Mesulam, 2000) largely decreased during PI-LOC (Table 3). We found significant alterations of functional connectivity between association and subcortical subsystems (t(10) = −4.19, pFDR < 0.05). Moreover, nominally significant differences were noted for paralimbic and subcortical subsystems (t(10) = −2.49, p = 0.033, uncorrected) and primary sensory and subcortical subsystems (t(10) = −2.87, p = 0.016, uncorrected). We also observed decreases in connectivity strength within association subsystems (t(10) = −2.89, p = 0.015, uncorrected), between association and primary subsystems (t(10) = −2.62, p = 0.025, uncorrected), and limbic and primary sensory subsystems (t(10) = −2.33, p = 0.041, uncorrected). There was no significant change of connectivity within the primary sensory areas or between regions of the limbic subsystem. A repeated-measures ANOVA revealed a significant effect of subsystem on connectivity strength (F(14,140) = 24.58, p < 0.001) and a trend for a significant effect of state on connectivity strength (F(1,10) = 4.68, p = 0.056). The interaction between subsystem and state was significant (F(14,140) = 3.72, p < 0.001). Post hoc paired t tests on the delta scores (PI-LOC minus wakefulness) of each subsystem revealed that connectivity between the subcortical and association subsystem decreased more than connectivity within the association subsystem (t(10) = −1.98, one-sided p = 0.038) and than connectivity within the subcortical subsystem (t(10) = −2.21, one-sided p = 0.025).

View this table:
  • View inline
  • View popup
Table 3.

Differences in subsystem connectivity

Effects of PI-LOC on long- and short-range connectivity

There was a significant effect of PI-LOC on the ratio of long- to short-range connectivity (Fig. 3B). We found a significant decrease in long-range connections (>78 mm; delta scores: mean(wakefulness) − mean(PI-LOC); one-sample t(10) = 2.86, p = 0.016) and in short-range connections (<78 mm; one-sample t(10) = 2.22, p = 0.050). A paired t test on the mean delta scores revealed a stronger effect on long-range than on short-range connections (t(10) = 2.83, p = 0.016).

Discussion

In the present study, we applied graph theoretical analysis of spontaneous BOLD fluctuations in fMRI during PI-LOC, highlighting significant quantitative changes in the brain's spatiotemporal organization. First, we observed a decrease in general connectivity strength, particularly between higher-association cortices and subcortical regions. Second, changes in local and global connectivity metrics during PI-LOC revealed a decrease of whole-brain integration I, an increase in clustering ratio γ, and a reduced probability of high-degree nodes. Third, we found a pronounced reduction of long-range connections during PI-LOC.

Effects of PI-LOC on subcortical and cortical functional connectivity

Our results indicate that even though propofol acts rather globally on the brain's intrinsic connectivity, effects are not uniform and can be associated with structures that have been linked to functional circuits of arousal and awareness (Franks, 2008). This is in line with PET studies showing a global, yet locally graded, reduction in brain metabolism during propofol anesthesia (Fiset et al., 1999; Kaisti et al., 2003).

Consistent with recent fMRI studies on the effect of propofol, we observed systemic decreases in subcortico-cortical and corticocortical connectivity strength (Boveroux et al., 2010). Regionally, we found significantly reduced connectivity strength of the thalamus, putamen, and several multimodal associative and sensory cortices, primarily located in posterior regions of the brain. The exact role of thalamus and thalamic functional connectivity in anesthesia-induced loss of consciousness is still debated (Alkire et al., 2008). Although thalamus is very likely involved in PI-LOC, observed changes in connectivity may be attributed to direct effects of propofol on thalamic nuclei or indirect effects via thalamocortical feedback loops. A study in Parkinson's patients reported that propofol-induction lead to slowing in cortical before subcortical electrophysiologic patterns, which may indicate a potentially secondary role of thalamus in initiating loss of consciousness (Velly et al., 2007). In contrast, for the onset of NREM sleep, deactivation of thalamus was found to precede that of cortex (Magnin et al., 2010). Different effects in the sensitivity of cortical and subcortical structures to propofol may also be related to the disparity in the distribution of GABAA receptors (Alkire and Haier, 2001). Ex vivo and PET studies revealed a lower GABAA receptor density in subcortical structures, such as the thalamus, than in structures of the neocortex (Braestrup et al., 1977; Gründer et al., 2001). To estimate the significance of changes in thalamic functional connectivity between states of vigilance, a better understanding of the nature of infra-slow oscillations (ISOs) seems required (Hughes et al., 2011). In particular, to what degree ISOs modulate behaviorally relevant EEG oscillations and whether their spatiotemporal organization promotes functional integration between distant brain regions requires further attention (Leopold et al., 2003; He and Raichle, 2009, He et al., 2010).

In addition to subcortico-cortical connectivity changes, significant connectivity decreases of parieto-occipital brain structures support the view that this complex may be important for conscious processing (Alkire et al., 2008). Coordinated interplay between posterior and frontal association cortices has been related to modulation of internal and external awareness (Boly et al., 2008). Evidence for reduced coupling between posterior, parietal, and frontal cortices was reported in recent studies during deep NREM sleep (Horovitz et al., 2009; Sämann et al., 2011) and anesthesia (Boveroux et al., 2010; Schrouff et al., 2011). Reorganized connectivity of posterior nodes, such as PCC, has been reported already for moderate levels of propofol (Stamatakis et al., 2010).

In accordance with previous studies, the connectivity within primary sensory cortices was not affected (Boveroux et al., 2010). The decrease in connectivity strength between low-level primary sensory and high-level association cortices seems to support the notion that the brain's ability to process sensory information up the cortical hierarchy is likely to be impaired during PI-LOC (Heinke and Schwarzbauer, 2001). It is not clear whether this observation also relates to a decrease in cortico-thalamocortical coupling, as disconnection of the thalamus from the functional brain graph may not only prohibit flow of ascending information to the cortex, but also alter integration of information between cortical areas (Theyel et al., 2010).

Effects of PI-LOC on topology of brain functional graphs

We observed that brain graphs demonstrate small-world topology during both wakefulness and PI-LOC. A relative maintenance of global principles may reflect an adaptive reconfiguration of brain organization (Bassett et al., 2006), a process that has also been suggested for propofol-induced transitions of consciousness in a recent EEG study (Lee et al., 2010). A significant increase in values of the clustering ratio γ was found during PI-LOC compared with wakefulness. Of note, this change resulted from a proportionally stronger reduction of Crandom than C. The larger decrease in Crandom may be related to the random rewiring of a more locally connected graph and to differences in the degree distributions between the two states. The increase in the clustering ratio γ could favor more local and less globally integrated information processing. This global-to-local switch is further supported by the significant decrease in long-range connections during PI-LOC.

Our data also suggests that anesthetics particularly affect functional connectivity hubs, as found during normal wakefulness (Tomasi and Volkow, 2011), showing a decreased incidence of high-degree nodes and a significant reduction in connectivity strength for structures such as thalamus and multimodal association cortices during PI-LOC. These findings were also accompanied by altered whole-brain integration, that is, a higher spatiotemporal decomposition of the functional network. Similar results have been reported by high-density EEG studies during deep NREM sleep (Massimini et al., 2005) and midazolam-induced sedation (Ferrarelli et al., 2010).

Taking the view of major theories of consciousness, an increase of local clustering and a decrease in global coherence could mirror a decoupling of functionally related cortical processing modules, possibly relevant for large-scale information integration in the brain (Mashour, 2006). Whether changes in levels of consciousness may be represented in ISOs of neuronal activity and how modifications in the organization and temporal coherence of ISOs affect brain function should be further investigated (He and Raichle, 2009).

Methodological considerations

Some limitations regarding our data have to be carefully considered. First, spontaneous BOLD signal fluctuations are an indirect measure of neural activity. Although there is mounting evidence on the neuronal origin of the BOLD signal and its correlation with local field potentials (LFPs), the exact nature and behavioral relevance of its fluctuations is not yet fully understood (Goense et al., 2012). Results in favor of a neuronal origin of BOLD signal fluctuations have been reported by He et al. (2008), showing a similar spatial correlation structure of BOLD signal fluctuations and LFPs for the frequency range of slow cortical potentials. Conversely, research has shown that BOLD signal fluctuations may contain diverse nonneuronal contributors, such as fluctuations from peripheral physiological sources, movement, or artifacts during data collection (Birn et al., 2006; Chang and Glover, 2009). A better understanding of the neurovascular component of BOLD signal fluctuations is still needed. Even though we controlled for nuisance in our data and restricted our analysis to a frequency band <0.10 Hz, we cannot fully exclude that wavelet correlations used in our study have been influenced by other nonneuronal variables. One also has to consider secondary effects of propofol, which may mediate regional specific BOLD responses and complicate inferences from functional connectivity. To directly test whether differential effects on connectivity are due to drug exposure or loss of consciousness, future studies should investigate changes in network metrics at different dosages of propofol. However, our results overlap largely with summaries on the molecular effect of propofol (Franks, 2008) and studies that used PET to investigate its effect on brain metabolism (Kaisti et al., 2003). It has to be stressed, however, that other dosages or types of anesthetics may show quite different effects on functional connectivity (Stamatakis et al., 2010).

Second, one has to note that present graph theoretical analysis provides only a macroscopic view on brain functional organization and previous research has shown that how brain graphs are constructed, compared, and normalized to reference models may change results considerably (Zalesky et al., 2010; Wang et al., 2011). To control for a methodological bias, we compared metrics across different brain parcellations and performed analyses on degree-preserved graph partitions. Since correlation strength of functional connectivity was significantly reduced during PI-LOC, the usage of strength as a single criterion to gain binary graphs would have lead to significant differences in the number of edges, and hence, different results with regard to measures such as the network's characteristic path length. The exact functional relationship between connection strength and topology requires further investigation (Rubinov and Sporns, 2011).

Conclusion

Our findings emphasize the significance of brain functional networks as measured by spontaneous BOLD fluctuations during resting fMRI. They support the notion that PI-LOC may be associated with a disintegration of the spatiotemporal architecture of brain activity. Results indicate that PI-LOC is associated with decreased subcortico-cortical connectivity and a breakdown of connectivity within higher-order association cortices and between higher-order association and primary sensory cortices. Mirroring observations during deep NREM sleep, topological results indicate that intrinsic activity during PI-LOC is more locally segregated and less globally integrated.

Footnotes

  • We thank the Brain Mechanisms of Consciousness and Anesthesia Research Group: Adem Omerovic, Stefan Hartl, Eva Pichlmaier, Sabine Paprotny, Christine Preibisch, and Susanne Neufang.

  • Correspondence should be addressed to Dr. Rüdiger Ilg, Klinikum rechts der Isar, Department of Neurology, Technische Universität München, Ismaningerstrasse 22, 81675 Munich, Germany. ilg{at}lrz.tu-muenchen.de

References

  1. ↵
    1. Achard S,
    2. Bullmore E
    (2007) Efficiency and cost of economical brain functional networks. PLoS Comput Biol 3:e17.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Achard S,
    2. Salvador R,
    3. Whitcher B,
    4. Suckling J,
    5. Bullmore E
    (2006) A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J Neurosci 26:63–72.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Alkire MT,
    2. Haier RJ
    (2001) Correlating in vivo anaesthetic effects with ex vivo receptor density data supports a GABAergic mechanism of action for propofol, but not for isoflurane. Br J Anaesth 86:618–626.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Alkire MT,
    2. Haier RJ,
    3. Barker SJ,
    4. Shah NK,
    5. Wu JC,
    6. Kao YJ
    (1995) Cerebral metabolism during propofol anesthesia in humans studied with positron emission tomography. Anesthesiology 82:393–403, discussion 27A.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Alkire MT,
    2. Hudetz AG,
    3. Tononi G
    (2008) Consciousness and anesthesia. Science 322:876–880.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Ames A 3rd.
    (2000) CNS energy metabolism as related to function. Brain Res Brain Res Rev 34:42–68.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Antkowiak B
    (1999) Different actions of general anesthetics on the firing patterns of neocortical neurons mediated by the GABA(A) receptor. Anesthesiology 91:500–511.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bandt C,
    2. Pompe B
    (2002) Permutation entropy: a natural complexity measure for time series. Phys Rev Lett 88:174102.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bassett DS,
    2. Bullmore E
    (2006) Small-world brain networks. Neuroscientist 12:512–523.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Bassett DS,
    2. Meyer-Lindenberg A,
    3. Achard S,
    4. Duke T,
    5. Bullmore E
    (2006) Adaptive reconfiguration of fractal small-world human brain functional networks. Proc Natl Acad Sci U S A 103:19518–19523.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Bassett DS,
    2. Bullmore E,
    3. Verchinski BA,
    4. Mattay VS,
    5. Weinberger DR,
    6. Meyer-Lindenberg A
    (2008) Hierarchical organization of human cortical networks in health and schizophrenia. J Neurosci 28:9239–9248.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Bassett DS,
    2. Brown JA,
    3. Deshpande V,
    4. Carlson JM,
    5. Grafton ST
    (2011) Conserved and variable architecture of human white matter connectivity. Neuroimage 54:1262–1279.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Birn RM,
    2. Diamond JB,
    3. Smith MA,
    4. Bandettini PA
    (2006) Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31:1536–1548.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Biswal B,
    2. Yetkin FZ,
    3. Haughton VM,
    4. Hyde JS
    (1995) Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med 34:537–541.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Biswal BB,
    2. Mennes M,
    3. Zuo XN,
    4. Gohel S,
    5. Kelly C,
    6. Smith SM,
    7. Beckmann CF,
    8. Adelstein JS,
    9. Buckner RL,
    10. Colcombe S,
    11. Dogonowski AM,
    12. Ernst M,
    13. Fair D,
    14. Hampson M,
    15. Hoptman MJ,
    16. Hyde JS,
    17. Kiviniemi VJ,
    18. Kötter R,
    19. Li SJ,
    20. Lin CP,
    21. et al.
    (2010) Toward discovery science of human brain function. Proc Natl Acad Sci U S A 107:4734–4739.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Boly M,
    2. Phillips C,
    3. Tshibanda L,
    4. Vanhaudenhuyse A,
    5. Schabus M,
    6. Dang-Vu TT,
    7. Moonen G,
    8. Hustinx R,
    9. Maquet P,
    10. Laureys S
    (2008) Intrinsic brain activity in altered states of consciousness: how conscious is the default mode of brain function? Ann N Y Acad Sci 1129:119–129.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Boveroux P,
    2. Vanhaudenhuyse A,
    3. Bruno MA,
    4. Noirhomme Q,
    5. Lauwick S,
    6. Luxen A,
    7. Degueldre C,
    8. Plenevaux A,
    9. Schnakers C,
    10. Phillips C,
    11. Brichant JF,
    12. Bonhomme V,
    13. Maquet P,
    14. Greicius MD,
    15. Laureys S,
    16. Boly M
    (2010) Breakdown of within- and between-network resting state functional magnetic resonance imaging connectivity during propofol-induced loss of consciousness. Anesthesiology 113:1038–1053.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Braestrup C,
    2. Albrechtsen R,
    3. Squires RF
    (1977) High densities of benzodiazepine receptors in human cortical areas. Nature 269:702–704.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Braun AR,
    2. Balkin TJ,
    3. Wesenten NJ,
    4. Carson RE,
    5. Varga M,
    6. Baldwin P,
    7. Selbie S,
    8. Belenky G,
    9. Herscovitch P
    (1997) Regional cerebral blood flow throughout the sleep-wake cycle: an H2(15)O PET study. Brain 120:1173–1197.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Brett M,
    2. Anton JL,
    3. Valabregue R,
    4. Poline JP
    (2002) Presented at the 8th International Conference on Functional Mapping of the Human Brain (Sendai, Japan), Region of interest analysis using a SPM toolbox.
  21. ↵
    1. Brown EN,
    2. Purdon PL,
    3. Van Dort CJ
    (2011) General anesthesia and altered states of arousal: a systems neuroscience analysis. Annu Rev Neurosci 34:601–628.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Bullmore E,
    2. Sporns O
    (2009) Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci 10:186–198.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Campbell K,
    2. Marsden P,
    3. Hurlbert J
    (1986) Social resources and socio-economic status. Soc Networks 8:97–117.
    OpenUrlCrossRef
  24. ↵
    1. Chang C,
    2. Glover GH
    (2009) Effects of model-based physiological noise correction on default mode network anti-correlations and correlations. Neuroimage 47:1448–1459.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Cordes D,
    2. Haughton VM,
    3. Arfanakis K,
    4. Carew JD,
    5. Turski PA,
    6. Moritz CH,
    7. Quigley MA,
    8. Meyerand ME
    (2001) Frequencies contributing to functional connectivity in the cerebral cortex in “resting-state” data. AJNR Am J Neuroradiol 22:1326–1333.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Dehaene S,
    2. Changeux JP
    (2005) Ongoing spontaneous activity controls access to consciousness: a neuronal model for inattentional blindness. PLoS Biol 3:e141.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Dosenbach NU,
    2. Nardos B,
    3. Cohen AL,
    4. Fair DA,
    5. Power JD,
    6. Church JA,
    7. Nelson SM,
    8. Wig GS,
    9. Vogel AC,
    10. Lessov-Schlaggar CN,
    11. Barnes KA,
    12. Dubis JW,
    13. Feczko E,
    14. Coalson RS,
    15. Pruett JR Jr.,
    16. Barch DM,
    17. Petersen SE,
    18. Schlaggar BL
    (2010) Prediction of individual brain maturity using fMRI. Science 329:1358–1361.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Ferrarelli F,
    2. Massimini M,
    3. Sarasso S,
    4. Casali A,
    5. Riedner BA,
    6. Angelini G,
    7. Tononi G,
    8. Pearce RA
    (2010) Breakdown in cortical effective connectivity during midazolam-induced loss of consciousness. Proc Natl Acad Sci U S A 107:2681–2686.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Fiset P,
    2. Paus T,
    3. Daloze T,
    4. Plourde G,
    5. Meuret P,
    6. Bonhomme V,
    7. Hajj-Ali N,
    8. Backman SB,
    9. Evans AC
    (1999) Brain mechanisms of propofol-induced loss of consciousness in humans: a positron emission tomographic study. J Neurosci 19:5506–5513.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Fox MD,
    2. Raichle ME
    (2007) Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci 8:700–711.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Fox MD,
    2. Corbetta M,
    3. Snyder AZ,
    4. Vincent JL,
    5. Raichle ME
    (2006) Spontaneous neuronal activity distinguishes human dorsal and ventral attention systems. Proc Natl Acad Sci U S A 103:10046–10051.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Fox MD,
    2. Snyder AZ,
    3. Vincent JL,
    4. Raichle ME
    (2007) Intrinsic fluctuations within cortical systems account for intertrial variability in human behavior. Neuron 56:171–184.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Fox MD,
    2. Zhang D,
    3. Snyder AZ,
    4. Raichle ME
    (2009) The global signal and observed anticorrelated resting state brain networks. J Neurophysiol 101:3270–3283.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Franks NP
    (2008) General anaesthesia: from molecular targets to neuronal pathways of sleep and arousal. Nat Rev Neurosci 9:370–386.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Goense J,
    2. Whittingstall K,
    3. Logothetis NK
    (2012) Neural and BOLD responses across the brain. WIREs Cogn Sci 3:75–86.
    OpenUrlCrossRef
  36. ↵
    1. Greicius MD,
    2. Krasnow B,
    3. Reiss AL,
    4. Menon V
    (2003) Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A 100:253–258.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Greicius MD,
    2. Kiviniemi V,
    3. Tervonen O,
    4. Vainionpää V,
    5. Alahuhta S,
    6. Reiss AL,
    7. Menon V
    (2008) Persistent default-mode network connectivity during light sedation. Hum Brain Mapp 29:839–847.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Gründer G,
    2. Siessmeier T,
    3. Lange-Asschenfeldt C,
    4. Vernaleken I,
    5. Buchholz HG,
    6. Stoeter P,
    7. Drzezga A,
    8. Lüddens H,
    9. Rösch F,
    10. Bartenstein P
    (2001) [18F]fluoroethylflumazenil: a novel tracer for PET imaging of human benzodiazepine receptors. Eur J Nucl Med 28:1463–1470.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Hagmann P,
    2. Cammoun L,
    3. Gigandet X,
    4. Meuli R,
    5. Honey CJ,
    6. Wedeen VJ,
    7. Sporns O
    (2008) Mapping the structural core of human cerebral cortex. PLoS Biol 6:e159.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Hampson M,
    2. Driesen NR,
    3. Skudlarski P,
    4. Gore JC,
    5. Constable RT
    (2006) Brain connectivity related to working memory performance. J Neurosci 26:13338–13343.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Hayasaka S,
    2. Laurienti PJ
    (2010) Comparison of characteristics between region-and voxel-based network analyses in resting-state fMRI data. Neuroimage 50:499–508.
    OpenUrlCrossRefPubMed
  42. ↵
    1. He BJ,
    2. Raichle ME
    (2009) The fMRI signal, slow cortical potential and consciousness. Trends Cogn Sci 13:302–309.
    OpenUrlCrossRefPubMed
  43. ↵
    1. He BJ,
    2. Snyder AZ,
    3. Zempel JM,
    4. Smyth MD,
    5. Raichle ME
    (2008) Electrophysiological correlates of the brain's intrinsic large-scale functional architecture. Proc Natl Acad Sci U S A 105:16039–16044.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. He BJ,
    2. Zempel JM,
    3. Snyder AZ,
    4. Raichle ME
    (2010) The temporal structures and functional significance of scale-free brain activity. Neuron 66:353–369.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Heinke W,
    2. Schwarzbauer C
    (2001) Subanesthetic isoflurane affects task-induced brain activation in a highly specific manner: a functional magnetic resonance imaging study. Anesthesiology 94:973–981.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Horovitz SG,
    2. Fukunaga M,
    3. de Zwart JA,
    4. van Gelderen P,
    5. Fulton SC,
    6. Balkin TJ,
    7. Duyn JH
    (2008) Low frequency BOLD fluctuations during resting wakefulness and light sleep: a simultaneous EEG-fMRI study. Hum Brain Mapp 29:671–682.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Horovitz SG,
    2. Braun AR,
    3. Carr WS,
    4. Picchioni D,
    5. Balkin TJ,
    6. Fukunaga M,
    7. Duyn JH
    (2009) Decoupling of the brain's default mode network during deep sleep. Proc Natl Acad Sci U S A 106:11376–11381.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Hughes SW,
    2. Lorincz ML,
    3. Parri HR,
    4. Crunelli V
    (2011) Infraslow (<0.1 Hz) oscillations in thalamic relay nuclei basic mechanisms and significance to health and disease states. Prog Brain Res 193:145–162.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Humphries MD,
    2. Gurney K
    (2008) Network ‘small-world-ness’: a quantitative method for determining canonical network equivalence. PLoS One 3:e0002051.
    OpenUrlCrossRefPubMed
  50. ↵
    1. James R,
    2. Glen JB
    (1980) Synthesis, biological evaluation, and preliminary structure-activity considerations of a series of alkylphenols as intravenous anesthetic agents. J Med Chem 23:1350–1357.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Johnston AJ,
    2. Steiner LA,
    3. Chatfield DA,
    4. Coleman MR,
    5. Coles JP,
    6. Al-Rawi PG,
    7. Menon DK,
    8. Gupta AK
    (2003) Effects of propofol on cerebral oxygenation and metabolism after head injury. Br J Anaesth 91:781–786.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Jordan D,
    2. Stockmanns G,
    3. Kochs EF,
    4. Pilge S,
    5. Schneider G
    (2008) Electroencephalographic order pattern analysis for the separation of consciousness and unconsciousness: an analysis of approximate entropy, permutation entropy, recurrence rate, and phase coupling of order recurrence plots. Anesthesiology 109:1014–1022.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Jordan D,
    2. Steiner M,
    3. Kochs EF,
    4. Schneider G
    (2010) A program for computing the prediction probability and the related receiver operating characteristic graph. Anesth Analg 111:1416–1421.
    OpenUrlCrossRefPubMed
  54. ↵
    1. Kaisti KK,
    2. Långsjö JW,
    3. Aalto S,
    4. Oikonen V,
    5. Sipilä H,
    6. Teräs M,
    7. Hinkka S,
    8. Metsähonkala L,
    9. Scheinin H
    (2003) Effects of sevoflurane, propofol, and adjunct nitrous oxide on regional cerebral blood flow, oxygen consumption, and blood volume in humans. Anesthesiology 99:603–613.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Larson-Prior LJ,
    2. Zempel JM,
    3. Nolan TS,
    4. Prior FW,
    5. Snyder AZ,
    6. Raichle ME
    (2009) Cortical network functional connectivity in the descent to sleep. Proc Natl Acad Sci U S A 106:4489–4494.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Latora V,
    2. Marchiori M
    (2001) Efficient behavior of small-world networks. Phys Rev Lett 87:198701.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Laureys S,
    2. Perrin F,
    3. Faymonville ME,
    4. Schnakers C,
    5. Boly M,
    6. Bartsch V,
    7. Majerus S,
    8. Moonen G,
    9. Maquet P
    (2004) Cerebral processing in the minimally conscious state. Neurology 63:916–918.
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Lee U,
    2. Oh G,
    3. Kim S,
    4. Noh G,
    5. Choi B,
    6. Mashour GA
    (2010) Brain networks maintain a scale-free organization across consciousness, anesthesia, and recovery: evidence for adaptive reconfiguration. Anesthesiology 113:1081–1091.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Leopold DA,
    2. Murayama Y,
    3. Logothetis NK
    (2003) Very slow activity fluctuations in monkey visual cortex: implications for functional brain imaging. Cereb Cortex 13:422–433.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Liu Y,
    2. Liang M,
    3. Zhou Y,
    4. He Y,
    5. Hao Y,
    6. Song M,
    7. Yu C,
    8. Liu H,
    9. Liu Z,
    10. Jiang T
    (2008) Disrupted small-world networks in schizophrenia. Brain 131:945–961.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Lydic R,
    2. Baghdoyan HA
    (2005) Sleep, anesthesiology, and the neurobiology of arousal state control. Anesthesiology 103:1268–1295.
    OpenUrlCrossRefPubMed
  62. ↵
    1. Lynall ME,
    2. Bassett DS,
    3. Kerwin R,
    4. McKenna PJ,
    5. Kitzbichler M,
    6. Muller U,
    7. Bullmore E
    (2010) Functional connectivity and brain networks in schizophrenia. J Neurosci 30:9477–9487.
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Magnin M,
    2. Rey M,
    3. Bastuji H,
    4. Guillemant P,
    5. Mauguière F,
    6. Garcia-Larrea L
    (2010) Thalamic deactivation at sleep onset precedes that of the cerebral cortex in humans. Proc Natl Acad Sci U S A 107:3829–3833.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Marsh B,
    2. White M,
    3. Morton N,
    4. Kenny GN
    (1991) Pharmacokinetic model driven infusion of propofol in children. Br J Anaesth 67:41–48.
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. Mashour GA
    (2006) Integrating the science of consciousness and anesthesia. Anesth Analg 103:975–982.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Maslov S,
    2. Sneppen K
    (2002) Specificity and stability in topology of protein networks. Science 296:910–913.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Massimini M,
    2. Ferrarelli F,
    3. Huber R,
    4. Esser SK,
    5. Singh H,
    6. Tononi G
    (2005) Breakdown of cortical effective connectivity during sleep. Science 309:2228–2232.
    OpenUrlAbstract/FREE Full Text
  68. ↵
    1. Mesulam MM
    (2000) Principles of behavioral and cognitive neurology (Oxford UP, Oxford), 2nd edition.
  69. ↵
    1. Mhuircheartaigh RN,
    2. Rosenorn-Lanng D,
    3. Wise R,
    4. Jbabdi S,
    5. Rogers R,
    6. Tracey I
    (2010) Cortical and subcortical connectivity changes during decreasing levels of consciousness in humans: a functional magnetic resonance imaging study using propofol. J Neurosci 30:9095–9102.
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Murphy K,
    2. Birn RM,
    3. Handwerker DA,
    4. Jones TB,
    5. Bandettini PA
    (2009) The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage 44:893–905.
    OpenUrlCrossRefPubMed
  71. ↵
    1. Percival DB,
    2. Walden AT
    (2000) Wavelet methods for time series analysis (Cambridge UP, Cambridge).
  72. ↵
    1. Ramsay MA,
    2. Savege TM,
    3. Simpson BR,
    4. Goodwin R
    (1974) Controlled sedation with alphaxalone-alphadolone. Br Med J 2:656–659.
    OpenUrlAbstract/FREE Full Text
  73. ↵
    1. Rubinov M,
    2. Sporns O
    (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52:1059–1069.
    OpenUrlCrossRefPubMed
  74. ↵
    1. Rubinov M,
    2. Sporns O
    (2011) Weight-conserving characterization of complex functional brain networks. Neuroimage 56:2068–2079.
    OpenUrlCrossRefPubMed
  75. ↵
    1. Sämann PG,
    2. Wehrle R,
    3. Hoehn D,
    4. Spoormaker VI,
    5. Peters H,
    6. Tully C,
    7. Holsboer F,
    8. Czisch M
    (2011) Development of the brain's default mode network from wakefulness to slow wave sleep. Cereb Cortex 21:2082–2093.
    OpenUrlAbstract/FREE Full Text
  76. ↵
    1. Salvador R,
    2. Suckling J,
    3. Schwarzbauer C,
    4. Bullmore E
    (2005) Undirected graphs of frequency-dependent functional connectivity in whole brain networks. Philos Trans R Soc Lond B Biol Sci 360:937–946.
    OpenUrlAbstract/FREE Full Text
  77. ↵
    1. Schölvinck ML,
    2. Maier A,
    3. Ye FQ,
    4. Duyn JH,
    5. Leopold DA
    (2010) Neural basis of global resting-state fMRI activity. Proc Natl Acad Sci U S A 107:10238–10243.
    OpenUrlAbstract/FREE Full Text
  78. ↵
    1. Schrouff J,
    2. Perlbarg V,
    3. Boly M,
    4. Marrelec G,
    5. Boveroux P,
    6. Vanhaudenhuyse A,
    7. Bruno MA,
    8. Laureys S,
    9. Phillips C,
    10. Pélégrini-Issac M,
    11. Maquet P,
    12. Benali H
    (2011) Brain functional integration decreases during propofol-induced loss of consciousness. Neuroimage 57:198–205.
    OpenUrlCrossRefPubMed
  79. ↵
    1. Schwarz AJ,
    2. McGonigle J
    (2011) Negative edges and soft thresholding in complex network analysis of resting state functional connectivity data. Neuroimage 55:1132–1146.
    OpenUrlCrossRefPubMed
  80. ↵
    1. Seeley WW,
    2. Menon V,
    3. Schatzberg AF,
    4. Keller J,
    5. Glover GH,
    6. Kenna H,
    7. Reiss AL,
    8. Greicius MD
    (2007) Dissociable intrinsic connectivity networks for salience processing and executive control. J Neurosci 27:2349–2356.
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Shmueli K,
    2. van Gelderen P,
    3. de Zwart JA,
    4. Horovitz SG,
    5. Fukunaga M,
    6. Jansma JM,
    7. Duyn JH
    (2007) Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. Neuroimage 38:306–320.
    OpenUrlCrossRefPubMed
  82. ↵
    1. Simard D,
    2. Nadeau L,
    3. Kroger H
    (2005) Fastest learning in small-world neuronal networks. Phys Lett A 336:8–15.
    OpenUrlCrossRef
  83. ↵
    1. Smith SM,
    2. Jenkinson M,
    3. Woolrich MW,
    4. Beckmann CF,
    5. Behrens TE,
    6. Johansen-Berg H,
    7. Bannister PR,
    8. De Luca M,
    9. Drobnjak I,
    10. Flitney DE,
    11. Niazy RK,
    12. Saunders J,
    13. Vickers J,
    14. Zhang Y,
    15. De Stefano N,
    16. Brady JM,
    17. Matthews PM
    (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23(Suppl 1):S208–S219.
    OpenUrlCrossRefPubMed
  84. ↵
    1. Spoormaker VI,
    2. Schröter MS,
    3. Gleiser PM,
    4. Andrade KC,
    5. Dresler M,
    6. Wehrle R,
    7. Samann PG,
    8. Czisch M
    (2010) Development of a large-scale functional brain network during human non-rapid eye movement sleep. J Neurosci 30:11379–11387.
    OpenUrlAbstract/FREE Full Text
  85. ↵
    1. Sporns O,
    2. Zwi JD
    (2004) The small world of the cerebral cortex. Neuroinformatics 2:145–162.
    OpenUrlCrossRefPubMed
  86. ↵
    1. Stam CJ,
    2. Jones BF,
    3. Nolte G,
    4. Breakspear M,
    5. Scheltens P
    (2007) Small-world networks and functional connectivity in Alzheimer's disease. Cereb Cortex 17:92–99.
    OpenUrlAbstract/FREE Full Text
  87. ↵
    1. Stamatakis EA,
    2. Adapa RM,
    3. Absalom AR,
    4. Menon DK
    (2010) Changes in resting neural connectivity during propofol sedation. PLoS One 5:e14224.
    OpenUrlCrossRefPubMed
  88. ↵
    1. Theyel BB,
    2. Llano DA,
    3. Sherman SM
    (2010) The corticothalamocortical circuit drives higher-order cortex in the mouse. Nat Neurosci 13:84–88.
    OpenUrlCrossRefPubMed
  89. ↵
    1. Tomasi D,
    2. Volkow ND
    (2011) Association between functional connectivity hubs and brain networks. Cereb Cortex 21:2003–2013.
    OpenUrlAbstract/FREE Full Text
  90. ↵
    1. Tononi G
    (2004) An information integration theory of consciousness. BMC Neurosci 5:42.
    OpenUrlCrossRefPubMed
  91. ↵
    1. Tzourio-Mazoyer N,
    2. Landeau B,
    3. Papathanassiou D,
    4. Crivello F,
    5. Etard O,
    6. Delcroix N,
    7. Mazoyer B,
    8. Joliot M
    (2002) Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15:273–289.
    OpenUrlCrossRefPubMed
  92. ↵
    1. van den Heuvel MP,
    2. Stam CJ,
    3. Kahn RS,
    4. Hulshoff Pol HE
    (2009) Efficiency of functional brain networks and intellectual performance. J Neurosci 29:7619–7624.
    OpenUrlAbstract/FREE Full Text
  93. ↵
    1. Van Dijk KR,
    2. Sabuncu MR,
    3. Buckner RL
    (2012) The influence of head motion on intrinsic functional connectivity MRI. Neuroimage 59:431–438.
    OpenUrlCrossRefPubMed
  94. ↵
    1. Van Essen DC,
    2. Drury HA,
    3. Dickson J,
    4. Harwell J,
    5. Hanlon D,
    6. Anderson CH
    (2001) An integrated software suite for surface-based analyses of cerebral cortex. J Am Med Inform Assoc 8:443–459.
    OpenUrlAbstract/FREE Full Text
  95. ↵
    1. Velly LJ,
    2. Rey MF,
    3. Bruder NJ,
    4. Gouvitsos FA,
    5. Witjas T,
    6. Regis JM,
    7. Peragut JC,
    8. Gouin FM
    (2007) Differential dynamic of action on cortical and subcortical structures of anesthetic agents during induction of anesthesia. Anesthesiology 107:202–212.
    OpenUrlCrossRefPubMed
  96. ↵
    1. Veselis RA,
    2. Feshchenko VA,
    3. Reinsel RA,
    4. Beattie B,
    5. Akhurst TJ
    (2005) Propofol and thiopental do not interfere with regional cerebral blood flow response at sedative concentrations. Anesthesiology 102:26–34.
    OpenUrlCrossRefPubMed
  97. ↵
    1. Wang JH,
    2. Zuo XN,
    3. Gohel S,
    4. Milham MP,
    5. Biswal BB,
    6. He Y
    (2011) Graph theoretical analysis of functional brain networks: test-retest evaluation on short- and long-term resting-state functional MRI data. PLoS One 6:e21976.
    OpenUrlCrossRefPubMed
  98. ↵
    1. Watts DJ,
    2. Strogatz SH
    (1998) Collective dynamics of ‘small-world’ networks. Nature 393:440–442.
    OpenUrlCrossRefPubMed
  99. ↵
    1. Weissenbacher A,
    2. Kasess C,
    3. Gerstl F,
    4. Lanzenberger R,
    5. Moser E,
    6. Windischberger C
    (2009) Correlations and anticorrelations in resting-state functional connectivity MRI: a quantitative comparison of preprocessing strategies. Neuroimage 47:1408–1416.
    OpenUrlCrossRefPubMed
  100. ↵
    1. Zalesky A,
    2. Fornito A,
    3. Harding IH,
    4. Cocchi L,
    5. Yücel M,
    6. Pantelis C,
    7. Bullmore ET
    (2010) Whole-brain anatomical networks: does the choice of nodes matter? Neuroimage 50:970–983.
    OpenUrlCrossRefPubMed
  101. ↵
    1. Zhang D,
    2. Snyder AZ,
    3. Fox MD,
    4. Sansbury MW,
    5. Shimony JS,
    6. Raichle ME
    (2008) Intrinsic functional relations between human cerebral cortex and thalamus. J Neurophysiol 100:1740–1748.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

The Journal of Neuroscience: 32 (37)
Journal of Neuroscience
Vol. 32, Issue 37
12 Sep 2012
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Spatiotemporal Reconfiguration of Large-Scale Brain Functional Networks during Propofol-Induced Loss of Consciousness
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Spatiotemporal Reconfiguration of Large-Scale Brain Functional Networks during Propofol-Induced Loss of Consciousness
Manuel S. Schröter, Victor I. Spoormaker, Anna Schorer, Afra Wohlschläger, Michael Czisch, Eberhard F. Kochs, Claus Zimmer, Bernhard Hemmer, Gerhard Schneider, Denis Jordan, Rüdiger Ilg
Journal of Neuroscience 12 September 2012, 32 (37) 12832-12840; DOI: 10.1523/JNEUROSCI.6046-11.2012

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Spatiotemporal Reconfiguration of Large-Scale Brain Functional Networks during Propofol-Induced Loss of Consciousness
Manuel S. Schröter, Victor I. Spoormaker, Anna Schorer, Afra Wohlschläger, Michael Czisch, Eberhard F. Kochs, Claus Zimmer, Bernhard Hemmer, Gerhard Schneider, Denis Jordan, Rüdiger Ilg
Journal of Neuroscience 12 September 2012, 32 (37) 12832-12840; DOI: 10.1523/JNEUROSCI.6046-11.2012
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Conclusion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Memory Retrieval Has a Dynamic Influence on the Maintenance Mechanisms That Are Sensitive to ζ-Inhibitory Peptide (ZIP)
  • Neurophysiological Evidence for a Cortical Contribution to the Wakefulness-Related Drive to Breathe Explaining Hypocapnia-Resistant Ventilation in Humans
  • Monomeric Alpha-Synuclein Exerts a Physiological Role on Brain ATP Synthase
Show more Articles

Behavioral/Systems/Cognitive

  • Gravin Orchestrates Protein Kinase A and β2-Adrenergic Receptor Signaling Critical for Synaptic Plasticity and Memory
  • Generation of Intensity Selectivity by Differential Synaptic Tuning: Fast-Saturating Excitation But Slow-Saturating Inhibition
  • Episodic Reinstatement in the Medial Temporal Lobe
Show more Behavioral/Systems/Cognitive
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.