Abstract
In this functional magnetic resonance imaging study, we examined the effect of mild propofol sedation and propofol-induced unconsciousness on resting state brain connectivity, using graph analysis based on independent component analysis and a classical seed-based analysis. Contrary to previous propofol research, which mainly emphasized the importance of connectivity in the default mode network (DMN) and external control network (ECN), we focused on the salience network, thalamus, and brainstem. The importance of these brain regions in brain arousal and organization merits a more detailed examination of their connectivity response to propofol. We found that the salience network disintegrated during propofol-induced unconsciousness. The thalamus decreased connectivity with the DMN, ECN, and salience network, while increasing connectivity with sensorimotor and auditory/insular cortices. Brainstem regions disconnected from the DMN with unconsciousness, while the pontine tegmental area increased connectivity with the insulae during mild sedation. These findings illustrate that loss of consciousness is associated with a wide variety of decreases and increases of both cortical and subcortical connectivity. It furthermore stresses the necessity of also examining resting state connectivity in networks representing arousal, not only those associated with awareness.
Introduction
Consciousness is a phenomenon that has so far evaded detailed description. With the advent of modern brain imaging techniques, it has become possible to better locate changes in brain activity associated with altered consciousness. In the healthy awake brain, structured brain activity occurs even in the absence of a specific task. This spontaneous activity is organized into resting state networks (RSNs) (Damoiseaux et al., 2006). Each RSN consists of brain regions that have a comparable temporal pattern of spontaneous neuronal activity, which are said to be functionally connected. Resting state functional magnetic resonance imaging (resting state fMRI) is an established method used to examine changes in the brain's connectivity when consciousness is altered, assuming that low frequency changes in blood oxygenation reflect neuronal activity. It has been demonstrated that especially higher-order RSNs show decreases in connectivity when consciousness is reduced or lost. The default mode network (DMN), external control network (ECN), and salience network are considered to represent the major portion of higher-order RSNs. The DMN is known to show changes in the spatial pattern during altered states of consciousness like deep sleep, vegetative state, and general anesthesia (Boveroux et al., 2010; Guldenmund et al., 2012; Horovitz et al., 2009; Schrouff et al., 2011; Vanhaudenhuyse et al., 2009). More recent studies also show substantial disintegration of the ECN, an RSN with spontaneous neuronal activity that is anticorrelated to that of the DMN (Boveroux et al., 2010; Vanhaudenhuyse et al., 2011), with decreasing consciousness. However, although the DMN and ECN are associated with the mediation of internal and external awareness, respectively, and are thus logical networks of interest for anesthesia studies, these are not the only brain networks showing changes during loss of consciousness. For instance, little is known about connectivity of the salience network during mild sedation and unconsciousness. The salience network is an RSN that consists of the bilateral frontal insulae, anterior cingulate cortex, and presupplementary motor area (Seeley et al., 2007). Activity in the salience RSN is associated with the detection of novel, salient stimuli, and is thought to play a role in coordinating an adequate response by recruiting appropriate brain networks. Therefore, some authors suggest it might play a role in coordinating between DMN and ECN activity (Bonnelle et al., 2012; Menon and Uddin, 2010).
Furthermore, although some resting state fMRI literature discusses changing connectivity between higher-order, highly integrative networks, and the thalamus, much remains to be discovered about changing connectivity between thalamus, brainstem and lower-order networks (like the auditory and sensorimotor RSNs). The propofol study by Boveroux et al. (2010) showed an involvement of the thalamus and brainstem regions in propofol-induced unconsciousness. However, no separate seed-based analysis with these loci as seed regions was performed to document their connectivity changes in detail. The brainstem pontine tegmental area (PTA), mesencephalon, and thalamus are brain structures closely associated with brain arousal via the ascending reticular arousal system (Angeles Fernandez-Gil et al., 2010; Giacino et al., 2012). In fact, microinjection of pentobarbital in the mesopontine tegmental region in rats has been shown to result in a condition resembling general anesthesia (Devor and Zalkind, 2001). Furthermore, the thalamus is regarded as one of the key loci for information integration and brain RSN coordination (Brown et al., 2010; Pinault, 2011; Schiff, 2008; Schiff and Fins, 2007; Tang et al., 2011). In general, propofol-induced unconsciousness is thought to result from increased GABA-ergic transmission in the cortex and at inhibitory projections from the preoptic area of the hypothalamus to arousal regions in the brainstem (Brown et al., 2011). Physiological effects of propofol, such as atonia and apnea, can also be traced back to the drug's actions on the brainstem (Brown et al., 2011).
These findings stress the potential importance of brainstem, thalamus and salience RSN activity in propofol-induced mild sedation and loss of consciousness. Therefore, we further examined their functional connectivity with resting state fMRI. Potential thalamic and brainstem involvement in the DMN and auditory RSN was first tested using a recently developed analysis method based on independent component analysis (ICA). Next, we used a classical seed-based technique to further examine connectivity changes with regions in the salience network, thalamus, and brainstem. A seed was also placed in the hippocampus, to examine expected changes in connectivity with the DMN other than the previously described frontoparietal disconnection (Boveroux et al., 2010). Given its role in the detection of salient stimuli and thus direction of attention, we expected decreases in salience network connectivity with unconsciousness, while the PTA, mesencephalon and thalamus, being involved in brain arousal and information integration, are also likely to show a propofol-induced change in connectivity.
Materials and Methods
Subjects
We used previously published fMRI data from 20 healthy right-handed volunteers (Boveroux et al., 2010; Schrouff et al., 2011). One subject was excluded from the analysis due to the occurrence of hyperventilation, while data from two other subjects were left out as they were acquired during pilot sessions. Data from the remaining 17 volunteers (13 women and 4 men; mean age: 21.9 years; standard deviation: 1.9 years; none was under medication or pregnant) were used for our analysis. The study was approved by the Ethics Committee of the Faculty of Medicine of the University of Liège (University Hospital, Liège, Belgium) and subjects gave written informed consent.
Sedation protocol
Subjects fasted for at least 6 h for solids and 2 h for liquids before the sedation. They wore headphones and earplugs in the scanner. Propofol infusion, using a target controlled infusion device (Diprifusor©-algorithm, Pharmacokinetics and Pharmacodynamics Software Server, Department of Anesthesia, Stanford University, USA) to obtain constant effect-site concentrations, occurred via an intravenous catheter placed into a vein of the right forearm or hand. During all four levels of consciousness, the blood pressure, electrocardiogram, breathing frequency, and pulse oxymetry (Sp02) were continuously monitored. For the whole duration of the experiment, subjects were breathing spontaneously, while additional oxygen was delivered at 5 L/min via a loosely fitting plastic facemask. The level of consciousness was assessed using the Ramsay scale (Olson et al., 2007). The subject was asked twice per consciousness level assessment to strongly squeeze the hand of the investigator. The awake states before sedation and after recovery of consciousness were Ramsay 2 (strong squeezing of the hand), mild sedation was Ramsay 3 (clear but slow squeezing), and propofol-induced unconsciousness was Ramsay 5–6 (no response). In addition, a reaction time task was also given to the subject before and after each session to help define the level of consciousness. This reaction task consisted of a block of 20 beeps delivered via the headphones, and the subjects were asked to press a keypad as fast as they could. After reaching the desired effect-site concentration, a 5 min equilibration period was established. Mean propofol plasma concentrations for wakefulness, mild sedation, unconsciousness, and recovery were 0 μg/mL (standard deviation: 0), 1.71 μg/mL (standard deviation: 0.72), 3.02 μg/mL (standard deviation: 1.03), and 0.59 μg/mL (standard deviation: 0.28), respectively. These propofol measurements were based on arterial blood samples taken directly before and after each scan. Two certified anesthesiologists and complete resuscitation equipment were present throughout the experiment [for supplementary protocol information, see Boveroux et al. (2010)].
fMRI data acquisition
Functional images were acquired on a 3T Siemens Allegra scanner (Siemens AG, Munich, Germany; Echo Planar Imaging sequence using 32 slices, repetition time=2460 ms, echo time=40 ms, field of view=220 mm, voxel size=3.45×3.45×3 mm, matrix size=64×64×32). Ten minute acquisitions were made (mean: 253 scans, standard deviation: 74 scans) during four different states of consciousness: normal wakefulness, mild propofol sedation, propofol-induced unconsciousness, and recovery. The temporal order of mild sedation and propofol-induced unconsciousness sessions was randomized. A high-resolution T1 image was also made for each subject for coregistration purposes. Total scan time per condition was around 30 min.
ICA-based analysis
The first analysis method we used was an automated ICA-based process combined with graph theory. This approach was used to analyze propofol-induced connectivity changes of whole networks: the DMN and the auditory RSN. These RSNs were chosen as they, and regions anticorrelating with the DMN and auditory RSN, are well-described and robust. This is important for the selection of regions of interest (ROIs) representing the RSNs, a vital part of our selection procedure as will be described below. Furthermore, these two RSNs are thought to be good representatives of a “higher-order,” highly integrative RSN and a sensory “lower-order” RSN, respectively. Structural and functional magnetic resonance images were realigned, normalized, smoothed (8 mm kernel), and corrected for 3D motion using BrainVoyager QX software (R. Goebel, Brain Innovation, Maastricht, The Netherlands). ICA with 30 components yielded statistically independent RSNs and possible artifacts (Hyvarinen et al., 2001; McKeown et al., 1998). Independent components corresponding to the DMN were selected using a previously published automated selection technique (Maudoux et al., 2012; Soddu et al., 2012). For this technique, using an average DMN template calculated on the basis of resting state fMRI scans from 12 healthy subjects [4 women, 8 men; mean age 21 years; standard deviation: 3 years (Demertzi et al., 2011)], who were scanned for a previous study on a 3T Siemens Allegra scanner (10 min resting state with eyes closed), we defined 14 ROIs (10×10×10 mm) that were considered to be most representative of the DMN based on previous research (Fair et al., 2008; Fox et al., 2005) (Supplementary Material S1; Supplementary Data are available online at

Default mode network (DMN) integrity during wakefulness, mild sedation, unconsciousness, and recovery. Scalar maps (left in each panel), fingerprints (bottom left in each panel) and connectivity graphs (right in each panel) in four conditions: wakefulness (the right inferior parietal region of the DMN is just below the z-axis of the image)

Auditory resting state network (RSN) integrity during wakefulness, mild sedation, unconsciousness, and recovery. Scalar maps (left in each panel), fingerprints (bottom left in each panel), and connectivity graphs (right in each panel) in four conditions: wakefulness
Seed-based analysis
Although our automatic ICA method is a robust tool for analyzing the well-described DMN and auditory RSN, seed analysis was chosen for further analysis of the connective repertoire of ROIs arising from previous research in the context of consciousness. Regions selected for the seed-based analysis were the right anterior insular cortex (33, 22, 6), thalamus (−7, −16, 6 and 7, −16, 6 combined), mesencephalon (−2, −24, −8), PTA (−3, −18, −27), and hippocampus (25, −16, −15). Literature confirms strong relationships of these regions with cognition and consciousness: the PTA, mesencephalon, and thalamus are of interest because of their roles in cortex activation and coordination (Abulafia et al., 2009; Saper et al., 2005; Schiff, 2008), while the hippocampus has a pivotal role in memory processing and has a close association with the DMN, providing additional information about DMN integrity (Wang and Orser, 2011). The right anterior insular cortex has been associated with the making of go/no-go decisions (Aron et al., 2004; Duann et al., 2009) and is part of the salience network, which is thought to enable switching between internally-oriented thought (DMN-driven) and externally-oriented activity (ECN-driven) (Seeley et al., 2007). For the same 17 subjects that we used for our ICA-based analysis, and for all four consciousness conditions, we used SPM8 software (statistical parametric mapping, Wellcome Trust Centre for Neuroimaging,
Although whole brain family-wise error correction, and to lesser extent whole brain false discovery rate correction, are often employed in studies where relatively great changes in brain activity can be expected (such as sensory and motor processes in classical fMRI studies), more subtle changes will go unnoticed. In fMRI research, there is a trend towards dramatically trying to minimize the chance of finding false positives (type I errors), resulting in an increase in type II errors (missing true effects). However, a great number of highly reproducible studies have used p<0.001 uncorrected thresholds (Lieberman and Cunningham, 2009). This threshold is generally considered to be a good balance between the chance of finding type I and II errors (Lieberman and Cunningham, 2009). Nevertheless, carefulness is advised with the interpretation of clusters consisting of less than 10 voxels. The employment of knowledge of the region's behavior in comparable experimental setups outlined in previous studies, used for applying small volume correction at a family-wise error of p<0.05 in this study, is a commonly employed procedure in the field of fMRI and further aids in minimizing potential type I errors.
Results
ICA of the DMN
During wakefulness before administration of propofol, DMN connectivity could be observed in all subjects (Fig. 2A). The DMN consisted of medial prefrontal cortex, posterior cingulate cortex/precuneus, and bilateral inferior parietal cortices (Beckmann et al., 2005; Damoiseaux et al., 2006). With mild propofol sedation, no significant changes were observed in DMN connectivity strength as compared to wakefulness (Fig. 2B, see Supplementary Material S2, DMN, for contrasts). The scalar contrast map showed that anticorrelation between the bilateral anterior insular cortices and the DMN had decreased from wakefulness to mild sedation (Supplementary Material S2, DMN).
During propofol-induced unconsciousness (Fig. 2C), the only difference found with the graph analysis between wakefulness and unconsciousness that survived bonferroni multiple comparisons correction was the loss of frontoparietal connectivity (Fig. 2C & Supplementary Material S2, DMN). Disconnection of the frontal part of the DMN was seen on the scalar map. As with mild sedation, the bilateral anterior insular cortices were less anticorrelated to the DMN as compared to wakefulness (Supplementary Material S2, DMN).
Recovery of consciousness was accompanied by the re-establishment of frontoparietal connectivity (Fig. 2D). No significant differences were observed between wakefulness before sedation and recovery, apart from weaker anticorrelation of the right anterior insular cortex with the DMN during recovery, as seen on the contrast scalar map (Supplementary Material S2, DMN). On the scalar map, two brainstem areas (the PTA and a region of the mesencephalon) showed up as being functionally connected to the DMN, whereby the mesencephalon demonstrated a significant increase in connectivity to the DMN during recovery from unconsciousness (Supplementary Material S2, DMN). We also found the thalamus and the cerebellum to be connected to the DMN during recovery (Fig. 2D).
The fingerprints from the chosen components of all four conditions were similar to the average fingerprint obtained from healthy, nonsedated controls (De Martino et al., 2007; Demertzi et al., 2011), and were thus likely to be neuronal. Movement parameters showed that there were no significant differences in patient movement between patients, or between scanning sessions.
ICA of the auditory RSN
The identified auditory RSN during wakefulness consisted of bilateral insular regions, the (supplementary) motor area and thalamus (Fig. 2A). An increase in connectivity between thalamus and insula was observed during mild sedation (Fig. 2B). No other significant differences were observed during either of the four consciousness conditions (Fig. 2C, D; Supplementary Material S2, auditory RSN). During wakefulness, we observed on the scalar maps an anticorrelation between the auditory RSN and two bilateral homotopic parts of the secondary visual cortex, while during mild sedation; connectivity was seen between the auditory RSN and a medial part of the secondary visual cortex.
Seed-based analysis
1. The first seed was placed in the anterior insular cortex (Supplementary Material S3). Decreased connectivity between the anterior insular cortex and bilateral inferior parietal cortices, precuneus, secondary sensorimotor area, and primary auditory cortex was found during mild sedation. Proceeding from mild sedation to unconsciousness, the anterior insular cortex also disconnected from the the salience network, thalamus, and superior temporal gyrus.
2. The second seed region of our analysis was the thalamus (Fig. 3 & Supplementary Material S3). With mild sedation, an increase in connectivity between the thalamus and the auditory cortex, insular cortex, primary somatosensory cortex, primary motor cortex, and supplementary motor area was observed. Decreased connectivity was found between thalamus and cerebellum. Compared to mild sedation, unconsciousness was accompanied by a decrease in thalamic connectivity with insular regions. Connectivity between thalamus and both the salience network regions and bilateral inferior parietal lobules was found to decrease from mild sedation to unconsciousness. Thalamic connectivity with regions overlapping with those of the DMN decreased during unconsciousness as compared to wakefulness. Connectivity between thalamus and supplementary motor and primary motor areas remained stronger than during full wakefulness.
3. In the mesencephalon seed region (Supplementary Material S3), mild sedation was associated with decreased connectivity with the posterior cingulate cortex/precuneus and the mesencephalon itself. Proceeding from mild sedation to unconsciousness, the seed region decreased connectivity with the posterior cingulate cortex/precuneus and left inferior parietal cortex. Compared to wakefulness, unconsciousness was associated with decreased connectivity with the posterior cingulate cortex/precuneus, inferior parietal cortices, and anterior cingulate cortex.
4. The fourth seed region was the PTA (Fig. 4 & Supplementary Material S3). During mild sedation, the PTA became more connected to the auditory and other insular regions. From mild sedation to unconsciousness, a decrease in connectivity compared to mild sedation was observed between PTA and auditory/insular areas, until preanesthesia levels of connectivity. Compared to wakefulness, during unconsciousness, a decrease in PTA connectivity with the posterior cingulate cortex/precuneus and right superior frontal cortex was found, while connectivity between PTA and putamen increased. During both mild sedation and unconsciousness, a decreased connectivity between PTA and thalamus was found.
5. The fifth seed region was the hippocampus (Supplementary Material S3). With mild sedation, we found a connectivity decrease between hippocampus and both posterior cingulate cortex/precuneus and primary/secondary visual cortex. From mild sedation to unconsciousness, a decrease in connectivity between hippocampus and the inferior parietal region was observed.

Changes in connectivity with the thalamus.

Changes in connectivity with the PTA.
Discussion
We examined the changing connectivity of the thalamus, brainstem and salience network during mild propofol sedation and propofol-induced unconsciousness. The results of this analysis complement the research performed by Boveroux et al. (2010), who focused on DMN and ECN connectivity changes associated with a decreasing level of consciousness. This previous research found propofol-induced decreased frontoparietal connectivity in the DMN, decreased anticorrelation between DMN and ECN, and decreased connectivity between thalamus and both DMN and ECN. Furthermore, they reported decreased connectivity between visual and auditory RSNs, while connectivity within the auditory and visual RSNs did not change significantly. We here present results from our analyses focusing on propofol-induced salience network, thalamus, and brainstem connectivity changes.
Salience network
Disconnection within the salience network and between salience network regions and thalamus was found with the thalamus and anterior insular cortex seeds. This disconnection has not been previously described in anesthesia and could play a major role in loss of consciousness. The influence of correct salience network functioning on brain connectivity was shown in a study with patients with cognitive impairment following mild traumatic brain injury (Bonnelle et al., 2012). The patients were performing a stop-signal task. Rapid deactivation of DMN regions was associated with efficient inhibitory control. Failure of this deactivation correlated with white matter damage in tracts connecting the salience network nodes. Menon and Uddin (2010) therefore propose that the salience network organizes DMN and ECN activity to account for an appropriate response to salient stimuli (Menon and Uddin, 2010). Disrupted salience network functioning has also been associated with a multitude of other cognitive disorders, among which are schizophrenia (Ellison-Wright and Bullmore, 2010), psychosis (Palaniyappan and Liddle, 2012), bipolar disorder (Ellison-Wright and Bullmore, 2010), autism (Menon and Uddin, 2010), and frontotemporal dementia (Seeley, 2008, 2010). Further research should elucidate how the disruptions in DMN and ECN functioning during propofol anesthesia relate to disrupted salience network functioning.
Thalamus
Seed analysis from the thalamus showed decreases in connectivity between thalamus and regions overlapping with the DMN and ECN, similar to those found in the study of Boveroux et al. (2010), as well as the salience RSN. Therefore, unconsciousness seems to correlate with disintegration of all three major higher-order networks, as well as their connectivity with the thalamus. Using positron emission tomography, thalamic metabolism has been shown to decrease significantly during anesthesia-induced unconsciousness (Fiset et al., 1999). Furthermore, a model has been suggested in which the thalamus orchestrates the commonly observed increased and coherent alpha frequency activity in the frontal cortex during propofol-induced unconsciousness (Ching et al., 2010; Cimenser et al., 2011). The authors suggest that this steady thalamic alpha rhythm could impede conduction and thus responsiveness to external stimuli. The fact that such a relatively active thalamus during unconsciousness has not been described for positron emission tomography and fMRI experiments could be because the signals obtained with these imaging methods are thought to be mostly associated with high-frequency activity, rather than alpha activity (Lachaux et al., 2007). It is not unlikely that links might exist between increased frontal cortex alpha activity and the disconnection of frontal cortex components of the DMN and salience network described in this paper. Indeed, new evidence supports the view that alpha activity in a certain brain area might be associated with inhibition in that area (Jensen and Mazaheri, 2010). In the awake brain, alpha activity is thought to serve maintenance and change of attentional focus by inhibition of distractors (Ben-Simon et al., 2013; Foxe and Snyder, 2011; Hanslmayr et al., 2011; Mathewson et al., 2011). In the unconscious brain, alpha activity might reflect inhibition of perception (Ching et al., 2010).
We observed increased connectivity between thalamus and both sensorimotor and auditory/insular cortices during mild sedation (Fig. 5A). Although increased thalamic connectivity with these regions has previously been observed during mild sevoflurane sedation (Martuzzi et al., 2010), to our knowledge, no previous resting state fMRI reports exist on the persistence of such increased connectivity during unconsciousness (Fig. 5B). One possible partial explanation for this connectivity increase might be that brain disturbance by propofol favors functional connectivity between brain regions that make up an above average contribution to the functional and possibly structural connective repertoire of the central thalamus. This is the case for the sensorimotor cortices (Fair et al., 2010). However, further research should elucidate the origin of this increased connectivity.

Summary of propofol-induced connectivity changes.
Brainstem
The brainstem contains a major part of the ascending reticular arousal system (de Lecea et al., 2012; Saper et al., 2005). Using our ICA of the DMN, we detected involvement of the mesencephalon within the DMN with recovery from unconsciousness. With seed analysis, we found propofol-induced decreased connectivity between DMN regions, especially the posterior cingulate cortex/precuneus, and both PTA and mesencephalon. The posterior cingulate cortex/precuneus is one of the most active regions in the human brain during wakefulness, and diminished activity here has been associated with loss of consciousness due to anesthesia and disorders of consciousness (Laureys, 2008; Xie et al., 2011). The importance of the PTA in the ascending reticular arousal system is well illustrated in the case of brainstem coma, a state in which there is a total absence of brain arousal and awareness (Laureys, 2008), where bilateral damage in the PTA alone can induce coma (Parvizi and Damasio, 2003). Recovery of brainstem functioning during recovery from coma shows similarities with recovery from anesthesia (Brown et al., 2010).
The mesencephalon has been shown to have solid connections with the posterior cingulate cortex/precuneus and anterior cingulate cortex in a multitude of studies and is considered to be a key region in pain saliency processing (Linnman et al., 2012). Moreover, two nuclei in the mesencephalon, the periaquaductal gray matter and the parabrachial nucleus, are thought to have a prominent role in modulating global brain activity. Projections have been found from these nuclei to the brain arousal areas of the reticular nuclei and basal forebrain, as well as intralaminar thalamic nuclei (Parvizi and Damasio, 2001).
Considering the brain arousal modulating capacities of the pontine and mesencephalic nuclei, their decreased connectivity with DMN regions during unconsciousness and the increasing connectivity between DMN and mesencephalon during recovery suggest brainstem-driven changes in brain arousal to partly underlie propofol-induced DMN connectivity decreases. These DMN decreases include disconnection of the frontal cortex and the hippocampus. A decreased connectivity between thalamus and PTA was also found; two pivotal entities in brain arousal that are known to work closely together (Saper et al., 2005). Decreased connectivity between them could potentially also represent a change in coordination of brain arousal.
We observed increased connectivity between PTA and insular cortices with mild sedation. These PTA connectivity changes have not been previously described and further research should indicate whether they could potentially underlie changes in brain processing associated with the insulae, such as auditory processing, salience detection, and self-perception (Menon and Uddin, 2010; Modinos et al., 2009; Saper et al., 2005). No increased PTA connectivity with the insular cortex was found during unconsciousness as compared to wakefulness. This suggests a nonlinear behavior during decreasing levels of consciousness, with increased connectivity during mild sedation but not unconsciousness. The reason for the absence of this increased connectivity between PTA and insulae during unconsciousness might be related to the state of general depression in the brain during unconsciousness, impeding salience detection and execution of behavior.
Methodological considerations
In our study, we used both ICA- and seed-based methods to try to examine connectivity decreases and increases during mild sedation and unconsciousness. For this, we used resting state fMRI, which assumes that low frequency changes in blood oxygenation reflect neuronal activity. ICA is an established data-driven method that is an adequate tool for exploring data and discriminating neuronal activity from artifact (Hyvarinen et al., 2001). It relies on the method of maximum discrimination of brain activity frequencies within a dataset given a certain number of components. However, potential bias might occur during the component selection step, as well as from the assumption that in every consciousness condition, remnants of the RSN of interest (e.g., DMN) are still present. The seed-based technique does not have these kinds of biases. However, given that proper preprocessing has been conducted, seed-based techniques might still suffer from the bias of ROI selection. Seed-based analysis shows the total connective repertoire of a chosen ROI. It looks at the degree of similarity of brain activity with all frequency patterns found in the ROI. Therefore, the techniques might be considered to be complementary, although comparing findings obtained using ICA with those gathered with seed-based methods requires some caution, especially when using different software packages for each method. However, we consider our results as being robust, as both the ICA- and seed-based methods found increased connectivity between thalamus and insular/auditory areas during mild sedation. Inspection of component fingerprints (De Martino et al., 2007) showed that it is possible to reliably find and study RSNs during propofol-induced unconsciousness.
Anesthesia could potentially affect the number of RSNs present in the brain. If this number changes, and we still only divide into 30 independent components with ICA, possibly RSNs are grouped together that would normally not be incorporated in the same independent component, or RSN regions could stop being associated with a certain RSN and become associated with another RSN in another independent component. A seeming disconnection or connectivity increase may thus be the result of a changing number of discriminable RSNs available in the brain. However, previously published research using seed-based methodology also showed propofol-induced disconnection of the frontal part of the DMN (Boveroux et al., 2010). Increased connectivity in the auditory RSN seen during mild sedation with ICA was also found with our seed-based method. Therefore, although it is imaginable that propofol might confound by changing discriminability of RSNs, we do not expect that this will affect our results greatly.
A possible explanation for the increased connectivity between thalamus and primary sensorimotor areas might be increased body movement during propofol sedation. Although we did not observe significant head motion, we did not perform electromyography to exclude significant body movement. However, a recent high-density electroencephalography study with propofol (Murphy et al., 2011) examined subjects at Ramsey 3 (mild sedation) and 5 (unconsciousness) levels of sedation; comparable to the sedation levels used in our study. Using electromyography, they found no significantly increased body movement during sedation stages. However, they report an increase in spindle (12–15 Hz) and beta (15–25 Hz) power during mild sedation, while increased gamma (25–40 Hz) power was observed in both mild sedation and unconsciousness. The authors excluded the possibility that these increases originated from ocular or muscular sources. Using electroencephalography source reconstruction, another report showed increased thalamic excitability during mild sedation and unconsciousness (Boly et al., 2012).
A further point of caution is the discrimination between connectivity and correlated activity. Although the central thalamic nuclei are often regarded as part of a vital brain arousal pathway with a broad connective range (Brown et al., 2010; Schiff, 2008), our findings cannot single out the possibility that increased connectivity with sensorimotor/insular regions might rather be a form of coincidental correlation. An in vivo cell recording study in rats with lesions in the central lateral intralaminar nucleus, destroying thalamic afferent connections, showed that brain arousal of local cortical networks might not be associated with thalamic afferent input (Constantinople and Bruno, 2011). Synchronous activity of thalamus and sensorimotor/insular cortices could therefore also be a result of both regions receiving similar input from another brain region, or being connected via that other region. However, other studies have shown significant decreases in brain arousal associated with central thalamic injury and significant deafferentation of its neurons due to diffuse brain insults in humans (Schiff and Plum, 2000; Shah and Schiff, 2010). Furthermore, sevoflurane-induced unconsciousness has been reversed by microinjection of nicotine, a cholinergic agonist, into the rat central thalamus (Alkire et al., 2007). These results emphasize the close direct relationship of the thalamus with the cortex, suggesting propofol influences direct connectivity between thalamus and cortical regions. It is possible that complementary arousal pathways provide an indirect connection between thalamus and cortical regions. An example would be the part of the ascending arousal pathway that does not run through the thalamus, but is influenced by it (Saper et al., 2005). Further research should elucidate the exact role of thalamic nuclei in cortical excitation during states of altered consciousness.
Another possible reason for carefulness when interpreting fMRI data might be the potential influence of pCO2 levels on the BOLD signal. However, it has been shown that pCO2 levels do not seem to change the BOLD response to neuronal activity (Birn et al., 2006; Wise et al., 2004). Furthermore, for our analysis, we are interested in correlations rather than specific regional effects and are therefore confident pCO2 levels do not significantly influence our results (Corfield et al., 2001).
Conclusion
Our study shows that propofol-induced mild sedation and unconsciousness is associated with decreased connectivity within the salience network, in addition to previously reported DMN and ECN breakdown. Thalamic connectivity with these three higher-order RSNs was also found to be reduced. Furthermore, brainstem connectivity with DMN-related regions decreased. In contrast, connectivity increases were observed between the thalamus and sensorimotor/insular/auditory cortices, as well as between the PTA and insular/auditory cortices. These findings give insight into how connectivity changes in DMN and ECN internal and external awareness networks might be related to salience network integrity and arousal components located in the thalamus and brainstem. It emphasizes the necessity of examining resting state connectivity of brain regions associated with both brain arousal and awareness together to make inferences upon mechanisms responsible for loss of consciousness.
Funding and Acknowledgments
This research was funded by the Belgian National Funds for Scientific Research (Brussels, Belgium), the European Commission (Brussels, Belgium), the James McDonnell Foundation (Saint Louis, Missouri, USA), the Mind Science Foundation (San Antonio, Texas, USA), the French Speaking Community Concerted Research Action (ARC, 06/11-340, Brussels, Belgium), the Fondation Médicale Reine Elisabeth (Brussels, Belgium), the University of Liège (Liège, Belgium) and the University Hospital of Liège (Liège, Belgium).
Footnotes
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
