Insights in dose dependent effects of Isoflurane by analyzing static and dynamic functional connectivity in mice
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).

Figures

Figure 1: Loss of connectivity between homotopic regions in the two hemispheres as indicated by cross-subjects box plots for two ICs using FSL Nets. The statistical analysis shows the presence of connectivity for isoflurane 1.1% while it is absent for higher doses (1.2% and 1.5%). The results shown are after Bonferroni correction.

Figure 2: Change of correlation coefficients for two functional modules, the lateral cortical network (LCN) and default mode network (DMN) as function of isoflurane dose. The analysis revealed that correlations coefficients between homotopic brain areas decreased with the increasing dose of isoflurane. Values for the correlation coefficients have been normalized to 1.1% dose. The error bars show the standard deviation of correlation coefficients between different regions within these networks.

Figure 3: Dynamic functional networks for isoflurane doses of 1.1%, 1,2% and 1.5% (from left to right). The networks comprising of DMN-Limbic and Piriform-Limbic becomes more focused at 1.2% and looses all the information at 1.5%.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1749