Qasim Bukhari1, Aileen Schröter1, and Markus Rudin1,2
1Department of Information Technology and Electrical Engineering, Institute of Biomedical Engineering, ETH and University of Zürich, Zürich, Switzerland, 2Institute of Pharmacology and Taxicology, University of Zürich, Zürich, Switzerland
Synopsis
The neurophysiological
effects of anesthetics on brain functional networks are not completely
understood. In this work we investigated the resting state functional brain networks
under different doses of isoflurane in mice. We used static and dynamic
functional connectivity (dFC) analysis to get insights in dose dependent
effects of isoflurane. The results from dFC analysis show that spatial
segregation across brain functional networks is lost with the increasing dose
of anesthesia thus it may be indicative of a deep anesthetic state. Static
network analysis using dual regression revealed loss of functional connectivity
between the bilateral regions, that is also supported with further results
showing decrease in functional correlations with increased dose of isoflurane.Purpose
Anesthesia
affects brain functional connectivity (Stuart et. al. 2006, Desai et. al. 2009). Yet, the link between neurophysiological effects of
anesthetics such as alteration of local/global neural activity and functional
networks and its clinical effects, i.e. loss of consciousness is not understood.
EEG as readout of neural activity has been widely used for studying underlying mechanisms
(Patrick et. al. 2013). However, EEG provides poor spatial resolution and thus
does not yield information on anesthesia-induced changes in cerebral functional
connectivity. In contrast, resting-state fMRI enables studying functional
networks and their modulation by pharmacological interventions, though at
inferior temporal resolution. In this study, we evaluated the effects of a
commonly used anesthetic, isoflurane, on the brain functional networks in mice using
resting-state fMRI. Both static and dynamic functional connectivity analysis have
been applied to capture the effect of the varying the anesthetic dose on
stationary and transient network patterns.
Method
Data
acquisition: Resting-state fMRI data sets of anesthetized mice were collected
on a Bruker BioSpec 94/30 system operating at 9.4T using a gradient-echo
echo-planar imaging (GE-EPI) sequence with 1s temporal resolution. Anesthesia
regimen included isoflurane with doses of 1.1%, 1.2% and 1.5% in an air/oxygen
mixture (4:1). Animals were mechanically ventilated throughout the experiment.
Data
analysis: After preprocessing and realignment, concat-ICA was applied using the
MELODIC toolbox of FSL followed by dual regression analysis and ‘randomize’. These
results were used for Network Modeling (Smith et. al. 2011) of the brain
networks using FSLNets. We also calculated the dynamic functional connectivity
networks using the dictionary learning approach (Leonardi et. al. 2013,
Leonardi et. al. 2014). Time courses were extracted from the ICA maps and fed
into the dFC algorithm. Simple building blocks (referred as ‘atoms’) of whole
brain connectivity were estimated using dictionary learning algorithm (Leonardi
et. al. 2014) with 30 folds each made of 200 iterations. The patterns thus
extracted explained more than 50% of the variance in the sliding window
correlation matrix. These building blocks were estimated for each anesthetic
dose separately.
Results
Dual
regression analysis revealed static functional connectivity between homotopic
regions in the two hemispheres. Network analysis was performed for selected ICs
using FSL Nets that revealed significant loss of functional connectivity
between the homotopic regions in the two hemispheres with the increased dose of
isoflurane (Fig. 1). Increasing
the isoflurane dose led to a general decrease in static correlation values all
over the brain, though the susceptibility varied across brain regions. This is
illustrated for two brain networks, the lateral cortical network (LCN) and the
default mode network (DMN), for which an average correlation values is given
(averaged over the various ROIs contained in the network, Fig.2).
The
dynamic functional connectivity analysis revealed significant interactions
among functional networks that were not apparent from the conventional analysis
(Fig. 3). Louvain
algorithm was used for estimating the modular functional connectivity matrix. dFC
analysis using dictionary-learning method (Leonardi et. al. 2013) displayed
significant interaction between the various modules identified (LCN, DMN,
thalamic and limbic networks). Increasing the isoflurane dose led to decreased interactions
among these modules (comparison 1.1 and 1.2% isoflurane, Fig. 3). At 1.5%
isoflurane spatial segregation is largely lost, i.e. the individual atoms (dynamic functional
connectivity states) displayed no structure, i.e. they
were rather ‘all off’ or ‘all on’. We might speculate that this loss of spatial
segregation might be associated to loss of consciousness as suggested in
earlier studies, which observed ‘global’ synchronization across the brain
during deep anesthesia (Grandjean et al. 2014).
Conclusion
Brain
functional networks are dependent on the anesthetic doses. While static FC
analysis with increasing dose of anesthetic revealed a steady loss of
functional correlation among homotopic brain regions, dynamic FC analysis
revealed a loss of spatial segregation across brain functional networks,
indicative of a deep anesthetic state. Dynamic
functional network analysis allowed the analysis of network changes that were
not apparent in the conventional static analysis.
Acknowledgements
No acknowledgement found.References
1. Stuart R. Hameroff
(2006) The Entwined Mysteries of Anesthesia and Consciousness: Is There
a Common Underlying Mechanism? Anesthesiology 8 2006, Vol.105, 400-412
2. Desai M, Kahn I, Ishizawa Y, Brown EN, Graybiel AM, Moore
CI. The impact of sevofluorane anesthesia on tactile activation of cortex and
striatum: high-resolution 9.4T fMRI studies in squirrel monkeys.Society for
Neuroscience, Chicago, IL, October 17-21, 2009
3. Patrick L. Purdon et al. (2013) Electroencephalogram signatures of loss and recovery
of consciousness from propofol
PNAS
110(12)
4. Smith SM, et al. (2011) Network modelling methods for
FMRI. Neuroimage. 54:875-91
5. Grandjean J,
Schroeter A, et al. (2014) Optimization of anesthesia
protocol for resting-state fMRI studies in mice based on their differential
effects on the functional connectivity pattern. NeuroImage 102: 838-47
6. Leonardi N, et al.
(2013) Principal components of functional connectivity: a new approach to study
dynamic brain connectivity during rest. (Translated from eng) Neuroimage
83:937-950 (in eng).
7. Leonardi N, Shirer
WR, Greicius MD, & Van De Ville D (2014) Disentangling dynamic networks:
Separated and joint expressions of functional connectivity patterns in time.
(Translated from eng) Hum Brain Mapp 35(12):5984-5995 (in eng).