Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Differential modulation of thalamo-parietal interactions by varying depths of isoflurane anesthesia

  • Dongrae Cho ,

    Contributed equally to this work with: Dongrae Cho, Teo Jeon Shin

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Teo Jeon Shin ,

    Contributed equally to this work with: Dongrae Cho, Teo Jeon Shin

    Affiliation Department of Pediatric Dentistry and Dental Research Institute, School of Dentistry, Seoul National University, Seoul, Korea

  • Jinsil Ham,

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Dong-Hyuk Choi,

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Seonghyun Kim,

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Seongwook Jeong,

    Affiliation Department of Anesthesiology and Pain Medicine, Chonnam National University Medical School, Gwangju, South Korea

  • Hyoung-Ihl Kim,

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Jae Gwan Kim ,

    leebr@gist.ac.kr (BL); jaekim@gist.ac.kr (JGK)

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

  • Boreom Lee

    leebr@gist.ac.kr (BL); jaekim@gist.ac.kr (JGK)

    Affiliation Department of Biomedical Science and Engineering, Institute of Integrated Technology, Gwangju Institute of Science and Technology, Gwangju, South Korea

Abstract

The thalamus is thought to relay peripheral sensory information to the somatosensory cortex in the parietal lobe. Long-range thalamo-parietal interactions play an important role in inducing the effect of anesthetic. However, whether these interaction changes vary with different levels of anesthesia is not known. In the present study, we investigated the influence of different levels of isoflurane-induced anesthesia on the functional connectivity between the thalamus and the parietal region. Microelectrodes were implanted in rats to record local field potentials (LFPs). The rats underwent different levels of isoflurane anesthesia [deep anesthesia: isoflurane (ISO) 2.5 vol%, light anesthesia (ISO 1 vol%), awake, and recovery state] and LFPs were recorded from four different brain areas (left parietal, right parietal, left thalamus, and right thalamus). Partial directed coherence (PDC) was calculated for these areas. With increasing depth of anesthesia, the PDC in the thalamus-to-parietal direction was significantly increased mainly in the high frequency ranges; however, in the parietal-to-thalamus direction, the increase was mainly in the low frequency band. For both directions, the PDC changes were prominent in the alpha frequency band. Functional interactions between the thalamus and parietal area are augmented proportionally to the anesthesia level. This relationship may pave the way for better understanding of the neural processing of sensory inputs from the periphery under different levels of anesthesia.

Introduction

General Anesthesia is defined as the state in which reversible loss of consciousness, sensation, motor, and reflex function are induced. In order to be able to perform surgery associated with highly noxious stimuli, general anesthetics should have hypnotic, amnesic, and analgesic effect. Although the mechanism of unconsciousness during anesthesia has been extensively investigated in previous studies [13], the neurophysiological mechanisms of the analgesic effect induced by anesthesia have attracted little interest.

Recently, the functional role of the thalamus as a consciousness switch has gained attention as a putative process underlying the hypnotic and amnesic effect of anesthetics [4]. However, a large body of evidence has shown that the thalamus processes sensory information from the periphery and relays it to the sensory cortex [5]. The integration of information between two brain areas is necessary to information processing [6,7]. Therefore, we infer that many anesthetics can modulate the interaction between the thalamus and the sensory cortex. In agreement with this hypothesis, several lines of evidence suggest that anesthetics alter thalamo-cortical functional connectivity [3,8,9]. High concentrations of anesthetics are usually administered before highly invasive surgical procedures because it is believed that a deep level of anesthesia is needed to deal with the elevated surgical stress. This procedure suggests that the level of anesthetic depth may affect sensory processing, therefore affecting the responses to noxious stimuli. However, to our knowledge, the effects of varying depths of anesthesia on the interaction between the thalamus and the parietal area have not been fully investigated.

Therefore, in the present study, we investigate how different levels of anesthesia induced by a volatile anesthetic (e.g. isoflurane) affect the interaction between the thalamus and the parietal region.

Materials and methods

Microelectrode implantation

All animal experiments complied with the institutional guidelines of the Gwangju Institute of Science and Technology (GIST). All procedures were approved by the Institutional Animal Care and Use Committee of GIST. Ten male Long-Evans rats that weighted 504 ± 106 g (mean ± standard deviation) were used for this experiment. The rats were single-caged in a controlled animal facility at 21 ± 1°C with water ad libitum. All animal care units were preserved on a 12-hour light-dark cycle with light on at 7:00 h. The microelectrodes were implanted in the parietal lobe and thalamic region in order to measure local field potentials (LFP). The rat was anesthetized using a mixture of ketamine hydrochloride (100mg/kg) and xylazine (7mg/kg). And then, its head was fixed inside a stereotaxic instrument. The incisor bar was moved to bring the bregma and lambda at the same height. We performed a durotomy (i.e. piercing a tiny hole in the dura) to make a descending path for the microelectrodes. The following four brain areas were selected as regions of interest for acquiring LFPs: left parietal lobe (LP), right parietal lobe (RP), left thalamic (LT) and right thalamic (RT) regions. Stereotaxic coordinates of the target areas (in mm) mentioned in the order (anterior or posterior to bregma, lateral to midline, ventral to dura mater) were as follows: (LP) (−4, +5.5, 0), RP (−4, −5.5, 0), LT (−3, +3, −6.5), and RT (−3, −3, −6.5). Fig 1B shows the brain histology indicating that the deep electrodes reached both the thalami. The implanted microelectrodes (2 skull microelectrodes for parietal lobe and 2 deep microelectrodes for thalamic regions) were tightly fixed to the skull using bone cement. The wound on the scalp was sutured after removing the rat from the stereotaxic frame. After the suture, the rats were placed in a chamber for recovering from the anesthesia and ketoprofen (2.5mg/kg, Uni Biotech, Korea) was injected intramuscularly to reduce pain. After the surgery, antibiotics (Ceftazol 20 mg/kg, Guju Pharma Co, Korea) was also injected intramuscularly for 3 days. All rats had a one-week postoperative recovery period. In addition, we frequently observed the recovery process in this period.

thumbnail
Fig 1. LFP recording under the four different stages and the locations of deep microelectrode in the thalamus.

(A): LFPs were recorded from the right parietal (RP), left parietal (LP), right thalamus (RT) and left thalamus (LT) in the awake, deep anesthesia, light anesthesia and recovery stages. (B): Deep electrodes (red arrows) reached bilateral thalami (black dash line).

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

LFP recording and data acquisition

A week after the implantation of the microelectrodes, we performed LFP recording on rats under isoflurane inhalational anesthesia. The overall protocol of this experiment has been previously described [9]. Briefly, the microelectrodes on the skull were connected to a customized LFP recording device and the acquired LFP was digitized with a Ni-DAQ card (NI PXIe-6368, National Instruments Corp.; 300 Hz sampling frequency). Initially, LFPs from the target areas were measured for 30 minutes in freely moving rats. The acquired LFPs were used as individual baseline references. LFPs were subsequently measured at different levels of isoflurane anesthesia. In order to induce a rapid loss of consciousness (LOC), we anesthetized the rat by delivering isoflurane (3 vol%) to an induction chamber for 3 minutes. We removed the unconscious rat from the induction chamber and placed the mouth of the rat in the nasal cone for delivering the inhalation anesthetic. Following this procedure, we administrated isoflurane at a concentration of 2.5% for 7 minutes. In particular, the first 2 minutes at the 2.5 vol% isoflurane stage were sufficient for reaching a pseudo state-equilibrium between blood plasma and the effective isoflurane concentration. The last 5 minutes at the 2.5 vol% isoflurane stage maintained a stable concentration of the anesthetic. After sustaining 2.5-vol% isoflurane for 5 minutes, the concentration of inhalational anesthetic was decreased to 1 vol% and this concentration was further administrated for 5 minutes. After being exposed to 1 vol% isoflurane for 5 minutes, the rats were allowed to recover from anesthesia for 30 minutes by removing nasal cone from the mouth of the rat. The LFPs in the target areas (i.e. parietal lobe and thalamus) were continuously recorded during the administration of isoflurane with a custom-made software developed in the LabVIEW (LabVIEW 2014, National Instruments Cop.) environment. Before data preprocessing, we chose the LFPs according to the level of anesthetic concentration. In order to analyze the changes of connectivity direction and strength between the regions of interest induced by different anesthesia levels, we selected four different 40 s LFPs, respectively. Since the selected LFP segments were artifact-free, each LEP segment have an information about stable anesthetic concentration. The selected segments were divided into the following four stages: awake (pre-anesthetic period), deep anesthesia (2.5 vol% isoflurane administration), light anesthesia (1 vol% isoflurane administration), recovery (post-anesthetic period) in Fig 1A. Therefore, the total number of segments was 40 (i.e. 10[subjects]×4[stages]).

Data preprocessing

The selected LFPs from the regions of interest were preprocessed before the application of the connectivity algorithm in order to remove the external noise induced by several external factors. Therefore, we conducted notch filtering to reduce the effect of the 60 Hz line noise generated by commercial electricity. Furthermore, we performed downsampling, i.e. the process of reducing the sampling rate of a signal in order to decrease the effect of added noise originating from the external environment. Since our brain wave range of interest was from 0.5 Hz to 60 Hz, we reduced the sampling rate of LFPs from 300 Hz to 150 Hz. In addition, to enhance the signal-to-noise ratio, we used averaged LFPs from the parietal lobe (P) and thalamus (T). These LFPs were calculated by averaging the signals from left (LP, LT) and right (RP, LT) sides of each area, respectively. Lastly, we divided the 40 s averaged LFPs (P, T) into 20 epochs (20[epochs]×2[seconds]) to apply the connectivity algorithm. The detailed description of epoch length and its effect are described in the next section. Therefore, the total number of epochs is 800 (10[subjects]×4[stages]×20[epochs]).

Partial Directed Coherence (PDC).

Partial Directed Coherence (PDC) derived by the multivariate autoregressive model (MVAR) is a popular connectivity method and was defined by Baccala and sameshiam [10]. Unlike Granger causality (GC), which is used to represent the causal relationship between two different signals in the time domain, PDC can describe the direction of information flow between two or more signals in the frequency domain. In other words, PDC extends the concept of directed influence from the time domain to the frequency domain and from two dimensional signals to multidimensional signals. Therefore, PDC is particularly well-suited to analyze brain connectivity using LFPs. Indeed, LFPs consist of multichannel signals and the characteristics of each wave (e.g. gamma, beta, alpha, theta and delta) are highly variable.

In the present study, we assume that the set of Wide Sense Stationary (WSS) signal x is estimated from multichannel LFPs and is described as (1) where N and t are the number of the LFP channel and the time index, respectively. The superscript T stands for the transpose operation. Let us assume that the MVAR process mentioned below is a proper description of the signal set x(t) (2)

The p donated model order of MVAR, residual noise e(t) = [e1(t), e2(t),…, eN(t)]T is the vector of size N and A(l) is the N×N model coefficients matrix. Since the performance of MVAR depends on the model order p, the determination of the appropriate mod’el order p can be obtained using an information theory approach [11]. Among the information theory methods available for selecting the proper criteria, we used the Akaike Information Criteria (AIC) for MVAR process [12]. The matrix of model coefficients A(l) are calculated by minimizing the residual e(t). In brief, Eq (2) is transformed to matrix form in the frequency domain and it is described as (3) where ajk(f), ej(f) and xj(f) are calculated using a Fourier transform. Among them, an element ajk(f) of the model coefficient A(f) represents the influence that signal xj(f) exerts on another signal xj(f) in the frequency f. The matrix of model coefficient A(f) is obtained as (4) where I is a N×N identity matrix. Eventually, estimated PDC is calculated as (5) where θkj(f) represents the directional flow of information from k-th signal to j-th signal in the frequency f. Therefore, the result of PDC is a useful index that can be used to analyze the causality between two different signals in the frequency domain. In particular, the PDC process is well-suited to analyze LFPs because these signals have different characteristics depending on the frequency range. However, since the classical PDC based on MVAR is a stationary process, it is not appropriate to apply the classical PDC to a nonstationary signal such as the LFPs because it may produce unexpected results. In order to overcome this issue, the PDC process is applied to LFPs using a narrow window where an optimal window size has been selected [11]. Previous neurophysiological experiments and signal analyses have revealed that the selection of a proper window size for applying MVAR processes can estimate the appropriate stationary stochastic process as well. In the present study, we used a 2 s window size.

Statistical analysis.

The results of PDC were divided in the following frequency bands: delta (1–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), and gamma (30–60 Hz) and was averaged for each respective frequency band. Therefore, the total number of PDC feature is 4000 (10 [subjects]×4 [stages]×20 [epoch]×5 [bands]). The overall PDC process structure is represented in Fig 2. The statistical differences in PDC for each frequency band at different anesthetic states were calculated by a Friedman test followed by Dunn’s post-hoc test using a free statistical software R (R 2.13.1) [13]. A P-value less than 0.05 was considered statistically significant.

thumbnail
Fig 2. Block diagram for calculating the Partial Directed Coherence (PDC).

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

Results

To confirm the direction and strength of interactions between thalamus and parietal area affected by different levels of anesthesia, PDC was calculated and averaged for all subjects. As shown in Figs 3 and 4, PDC changed at different anesthesia levels for both the thalamus-to-parietal and parietal-to-thalamus directions. The PDC for both directions (thalamus → parietal, parietal → thalamus) showed similar trends as anesthesia depth increased. The PDC in the thalamus-to-parietal direction significantly increased with increasing anesthesia mainly in the high frequency range. PDC increase associated with anesthesia depth changes was most prominent in the alpha frequency band (F = 18.72, P<0.001), suggesting that interactions between the thalamus and parietal areas mainly occur in the alpha frequency range. However, anesthesia-related changes in the parietal-to-thalamus PDC were only significant in in the low frequency range. PDC increase associated with anesthesia depth changes was most prominent in the alpha frequency band (F = 12.77, P = 0.0052), similar to what we observed in the reverse direction. This result suggests that the parietal-to-thalamus interactions mainly occur in the alpha frequency range. The detailed description of change in PDC between parietal and thalamus associated with anesthesia depth changes is present in S1 Table (see supporting information).

thumbnail
Fig 3. Averaged PDC computed in the thalamus-to-parietal direction at each frequency band according to the anesthesia depth.

Error bars indicate standard error. GA: general anesthesia, P: parietal area, T: thalamus.

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

thumbnail
Fig 4. Averaged PDC computed in the parietal-to-thalamus direction at each frequency band according to the anesthesia depth.

Error bars indicate standard error. GA: general anesthesia, P: parietal area, T: thalamus.

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

Taken together, these results suggest that the pattern of connectivity between parietal and thalamic regions is influenced by the depth of anesthesia (Fig 5).

thumbnail
Fig 5. Estimated connectivity of the parietal area and thalamus using PDC.

The solid line indicates an increase in connectivity between the two areas. The frequency shown on the line indicates the frequency bands between the two areas that are statistically significant (P < 0.05). T: thalamus, P: parietal area.

https://doi.org/10.1371/journal.pone.0175191.g005

Discussion

We have shown that PDC preferentially changes in the thalamus-to-cortex direction depending on the anesthesia depth. To the best of our knowledge, this is the first study to investigate large-scale thalamo-parietal interactions across the brain at different levels of anesthesia.

Several lines of evidence have shown that the functional connections between the thalamus and cortical areas are critical in maintaining consciousness. Information flow between the thalamus and the cortex is abruptly changed during the transitional phase of consciousness [2,3]. However, until now, the thalamo-cortical interaction changes induced by anesthesia depth have not been thoroughly investigated. Recently, we have found that different levels of anesthesia have a different impact on the interaction between the thalamus and the frontal cortex, a brain region critical in regulating cognitive function [9]. Our previous findings lead to the hypothesis that anesthesia depth per se may affect the functional connectivity between brain areas.

In agreement with our hypothesis and previously reported findings, different levels of anesthesia significantly affect the thalamo-parietal interaction. Considering that the analgesic effect usually increases proportionally to the depth of anesthesia, we predicted that bottom-up sensory processing may be impaired during anesthesia. However, contrary to our prediction, thalamo-parietal synchrony was proportionally increased, rather than being decreased, during anesthesia. Sensory perception can be understood as a hierarchical predictive coding between bottom-up and top down information processing [1416]. Anesthesia also affects the function of the frontal area, a region that plays a critical role in executive processing [1718]. If fragmented sensory inputs are not integrated within the framework of sensory information processing, the senses cannot be perceived. Taken together, our results and previous findings lead to the speculation that feedforward connections may be impaired during anesthesia while bottom-up sensory processing appears to be improved.

In this context, we speculate that thalamus-to-parietal bottom-up connectivity may not be reduced because this connectivity is not directed to a high order brain region. However, it is unclear why bottom-up processing is aroused. One conclusion would be that as anesthesia deepens, isoflurane-induced uneven modulation of cortico-thalamic loop might make these robust connections more prominent. Under anesthesia, spontaneous thalamic activity is regulated by feedback connections driven by cortical neuron activity [19]. If feedback connection attenuates the input to the thalamus associated with different anesthetic levels, the inherent and robust connections between thalamus and parietal area may be more prominent.

We have shown that parietal-thalamus interactions are rather increased contrary to frontal thalamus interactions in our previous study. It is not clear why thalamocortical interactions differ from one brain area to another with different anesthetic depths. However, these different results may be explained by considering that sensory stimulation is also an important factor controlling the level of consciousness. Indeed, sensory processing is also affected by modulation of thalamocortical interaction depending on attentional states [20]. Also, it has been known that the detection of complex sound sequences is not significantly impaired in the absence of consciousness [21], suggesting that sensory processing may be more preserved than consciousness processing. Increased thalamo-parietal interaction may be thought as an increase in sensory processing.

However, from the perspective of information integration theory, a synchronized network does not respond efficiently to incoming inputs. A study investigating the mechanisms underlying unconsciousness in generalized epilepsy reported an increase in the synchrony between cortical and subcortical structures [22]. The authors suggested that, on a global neuronal scale, enhanced synchrony paradoxically might block the integration of information relayed from other areas. It has also been reported that connections between thalamic nuclei involved in transferring sensory input from the periphery are relatively maintained while connections between nonspecific nuclei mediating cortical integration are depressed [23]. The parietal lobe receives major sensory inputs from the periphery through the thalamus. Although it is not clear whether the increase in interaction is due to a complementary mechanism that increases arousal or by robust interactions of itself, increased interactions may enhance the paradoxical inherent inter-brain interactions, rendering them less flexible to respond to various types of external stimuli.

The thalamic reticular nucleus is important in regulating the interaction between the thalamus and the cortex [24]. Recent studies have shown that TRN controls sensory processing associated with attention states [20,25]. Isoflurane, a volatile anesthetic, modulates the gamma-aminobutyric acid (GABA) receptors of TRN differently under different anesthetic conditions [26]. The TRN constitue of thin GABAergic neurons [27]. Considering that TRN is important for thalamocortical interaction and regulation of thalamocortical function in relation to sensory processing as well as the function control of thalamic nucleus for sensory processing, the alteration of GABAergic TRN neuron associated with different levels of anesthesia may affect reciprocal thalamo-parietal interaction although a further study warrants to clarify the mechanism of increased thalamocortical interactions from the perspective of neurotransmission.

Interestingly, modulations of thalamo-cortical and cortico-thalamic interactions were statistically significant in the alpha frequency band. Several studies have suggested that brain oscillation in the alpha band may be a core feature of communication between brain areas. The thalamus is known to be involved in generating alpha oscillation [2830]. Thalamo-cortical circuits are mainly coordinated by an alpha rhythm generated from the thalamus [31,32]. Therefore, it is surprising that reciprocal thalamo-parietal connections increase proportionally with the anesthesia depth in this study. Rather, thalamic activity is augmented during propofol induced ROC unlike depressed cortical activity [33,34]. Taken together, these findings suggest that the robust thalamo-cortical connections in the alpha frequency band are aroused as anesthesia deepens.

Given the large variety of anesthetic agents, it is likely that other types of anesthetics may exert different effects on thalamo-parietal connectivity at different levels of anesthesia. This variability should be taken into account when extending our findings to other types of anesthetics. While our study provides a plausible explanation of what happens to thalamo-parietal connections under different levels of anesthesia, it is unclear whether the increased thalamo-parietal interactions observed under anesthesia is the consequence of a direct action of the anesthetic on the thalamus and cortex or rather this effect is mediated by other brain areas not analyzed in the present study.

In conclusion, we have found that the functional interaction between the thalamus and parietal area is augmented proportionally to the anesthesia levels. The present findings provide a new insight into the processing of sensory inputs from the periphery under different levels of anesthesia. Indeed, the mechanisms underlying this anesthetic-dependent modulation deserve further investigations from the perspective of a global brain network.

Supporting information

S1 Table. Changes in thalamo-parietal connectivity at each frequency band and the time domain depending on anesthetic depth.

T: thalamus, P: frontal area, GA: general anesthesia, + indicates statistical significance (P<0.05), ns: statistically insignificant.

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

(DOCX)

Author Contributions

  1. Conceptualization: DRC TS JK BL.
  2. Data curation: DRC TS.
  3. Formal analysis: DRC TS.
  4. Funding acquisition: DRC TS JH DHC SK SJ HK JK BL.
  5. Investigation: DRC TS JH DHC SK.
  6. Methodology: DRC TS BL.
  7. Project administration: TS JK BL.
  8. Resources: DRC JH DHC SK.
  9. Software: DRC.
  10. Supervision: JK BL.
  11. Visualization: DRC TS.
  12. Writing – review & editing: TS BL.

References

  1. 1. Xu X, Wang G, Tian X. Effects of Anesthesia on Effective Connectivity in the Brain. World Journal of Neuroscience. 2015;5(2):99–107.
  2. 2. White NS, Alkire MT. Impaired thalamocortical connectivity in humans during general-anesthetic-induced unconsciousness. Neuroimage. 2003;19(2):402–411.
  3. 3. Mashour George A., Alkire Michael T.. Consciousness, Anesthesia, and the Thalamocortical System. Anesthesiology. 2013;118(1):13–15. pmid:23208518
  4. 4. Alkire M. T., Haier R. J., Fallon J. H.. Toward a Unified Theory of Narcosis: Brain Imaging Evidence for a Thalamocortical Switch as the Neurophysiologic Basis of Anesthetic-Induced Unconsciousness. Consciousness and Cognition. 2000;9(3):370–386. pmid:10993665
  5. 5. Sherman S. Murray. Thalamic relays and cortical functioning. Progress in Brain Research. 2005;149:107–126. pmid:16226580
  6. 6. Tononi G. Consciousness, information integration, and the brain. Progress in Brain Research. 2005;150:109–126. pmid:16186019
  7. 7. Tononi G. Information integration: its relevance to brain function and consciousness. Archives Italiennes de Biologie. 2010;148(3):299–322. pmid:21175016
  8. 8. Crone JS, Soddu A, Holler Y, Vanhaudenhuyse A, Schurz M, Bergmann J, et al. Altered network properties of the fronto-parietal network and the thalamus in impaired consciousness. NeuroImage: Clinical. 2014;4:240–248
  9. 9. Shin TJ,Cho D,Ham J,Choi DH,Kim S,Jeong S et al. Changes in thalamo-frontal interaction under different levels of anesthesia in rats. Neuroscience Letters. 2016;627:18–23. pmid:27230989
  10. 10. Baccala LA, Sameshima K. Partial directed coherence: a new concept in neural structure determination. Biological Cybernetics. 2001;84(6):463–474. pmid:11417058
  11. 11. Ding M, Bressler SL, Yang W, Liang H. Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment. Biological Cybernetics. 2000;83(1):35–45. pmid:10933236
  12. 12. Sun Y,Zhang H,Feng T,Qiu Y,Zhu Y,Tong S. Early Cortical Connective Network Relating to Audiovisual Stimulation by Partial Directed Coherence Analysis. IEEE Transactions on Biomedical Engineering. 2009;56(11):2721–2724.
  13. 13. R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2015.
  14. 14. Clark A. Whatever Next? Predictive Brains, Situated Agents, and the Future of Cognitive Science. Behavioral and Brain Sciences. 2013;36(3):181–204. pmid:23663408
  15. 15. Rao RP, Ballard DH. Predictive coding in the visual cortex: a functional interpretation of some extra-classical receptive-field effects. Nature Neuroscience. 1999;2(1):79–87. pmid:10195184
  16. 16. Gilbert CD, Sigman M. Brain states: Top-down influences in sensory processing. Neuron. 2007;54(5):677–696. pmid:17553419
  17. 17. Ku SW, Lee U, Noh GJ, Jun IG, Mashour GA. Preferential inhibition of frontal-to-parietal feedback connectivity is a neurophysiologic correlate of general anesthesia in surgical patients. PLoS ONE. 2011;6(10):e25155. pmid:21998638
  18. 18. Lee U, Kim S, Noh GJ, Choi BM, Hwang E, Mashour GA. The directionality and functional organization of frontoparietal connectivity during consciousness and anesthesia in humans. Consciousness and Cognition. 2009;18(4):1069–1078. pmid:19443244
  19. 19. Vahle-Hinz C, Detsch O, Siemers M, Kochs E. Contributions of GABAergic and glutamatergic mechanisms to isoflurane-induced suppression of thalamic somatosensory information transfer. Experimental Brain Research. 2007;176(1):159–172. pmid:16847609
  20. 20. Halassa MM, Chen Z, Wimmer RD, Brunetti PM, Zhao S, Zikopoulos B, et al. State-dependent architecture of thalamic reticular subnetworks. Cell. 2014;158(4):808–821. pmid:25126786
  21. 21. Tzovara A, Simonin A, Oddo M, Rossetti AO, Lucia MD. Neural detection of complex sound sequences in the absence of consciousness. Brain. 2015;138(5):1160–1166.
  22. 22. Arthuis M,Valton L,Régis J,Chauvel P,Wendling F, Naccache L et al. Impaired consciousness during temporal lobe seizures is related to increased long-distance cortical-subcortical synchronization. Brain. 2009;132(Pt 8):2091–2101. pmid:19416952
  23. 23. Liu X, Lauer KK, Ward BD, Li SJ, Hudetz AG. Differential effects of deep sedation with propofol on the specific and nonspecific thalamocortical systems: a functional magnetic resonance imaging study. Anesthesiology. 2013;118(1):59–69. pmid:23221862
  24. 24. Pratt JA and Morris BJ. The thalamic reticular nucleus: a functional hub for thalamocortical network dysfunction in schizophrenia and a target for drug discovery. Journal of Psychopharmacology. 2015;29(2):127–137. pmid:25586397
  25. 25. McAlonan K, Cavanaugh J, Wurtz RH. Attentional modulation of thalamic reticular neurons. Journal of Neuroscience. 2006;26(16):4444–4450. pmid:16624964
  26. 26. Gottschalk A and Miotke SA. Volatile Anesthetic Action in a Computational Model of the Thalamic Reticular Nucleus. Anesthesiology. 2009;110:996–1010. pmid:19387174
  27. 27. Houser CR, Vaughn JE, Barber RP, Roberts E. GABA neurons are the major cell type of the nucleus reticularis thalami. Brain Research. 1980;200(2):341–354. pmid:7417821
  28. 28. Hughes SW, Crunelli V. Thalamic mechanisms of EEG alpha rhythms and their pathological implications. Neuroscientist. 2005;11(4):357–372. pmid:16061522
  29. 29. Bollimunta A, Mo J, Schroeder CE, Ding M. Neuronal mechanisms and attentional modulation of corticothalamic α oscillations. Journal of Neuroscience. 2011;31(13):4935–4943. pmid:21451032
  30. 30. Jahnsen H, Llinás R. Ionic basis for the electro-responsiveness and oscillatory properties of guinea-pig thalamic neurones in vitro. The Journal of Physiology. 1984;349:227–247. pmid:6737293
  31. 31. Lorincz ML, Kekesi KA, Juhasz G, Crunelli V, Hughes SW. Temporal framing of thalamic relay-mode firing by phasic inhibition during the alpha rhythm. Neuron. 2009;63(5):683–696. pmid:19755110
  32. 32. Vijayan S, Kopell N.J. Thalamic model of awake alpha oscillations and implications for stimulus processing. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(45):18553–18558. pmid:23054840
  33. 33. Boly M, Moran R, Murphy M, Boveroux P, Bruno MA, Noirhomme Q et al. Connectivity changes underlying spectral EEG changes during propofol-induced loss of consciousness. The Journal of Neuroscience. 2013;32(20):7082–7090.
  34. 34. Velly LJ, Rey MF, Bruder NJ, Gouvitsos FA, Witjas T, Regis JM et al. Differential dynamic of action on cortical and subcortical structures of anesthetic agents during induction of anesthesia. Anesthesiology. 2007;107(2):202–212. pmid:17667563