1343

Diffusion functional MRI characterizes dynamical brain function in a neuropsychiatric disease model mouse
Yoshifumi Abe1
1Keio University School of Medicine, Tokyo, Japan

Synopsis

We propose an analytical framework to characterize dynamic brain function in neuropsychiatric conditions by taking advantage of the technical aspects of diffusion functional MRI (DfMRI). The pipeline consists of local activity analysis with apparent diffusion coefficient (ADC) data, functional connectivity (FC) analysis with diffusion-weighted data (Sb1800), and ignition-driven mean integration (IDMI) analysis combining both. We illustrated its utility by analyzing model mice with an obsessive–compulsive disorder (OCD)-related behavior. The framework was successful in detecting hyperactivation and biased connectivity across the cortico-striato-thalamic circuitry. The IDMI analysis found unseen local activity-initiated propagation to the global network.

Introduction

Functional MRI (fMRI) based on a blood oxygenation-level dependent (BOLD) mechanism is a powerful tool for investigating whole brain connectivity as well as brain activity. However, the BOLD signal sometimes fails when neurovascular coupling is disrupted (e.g. cerebrovascular disease, astrocytic dysfunction, intake of drugs, such as anesthetized condition)1-3. Thus, BOLD-fMRI under these conditions may exhibited false-positive or false-negative findings of brain network and activity.
An essential tool to resolve the drawback of BOLD-fMRI is diffusion fMRI (DfMRI) which has proposed a functional imaging technique as an alternative to BOLD-fMRI4. DfMRI has been shown to reflect neural activity more directly, compared with BOLD-fMRI, as DfMRI signals are not of vascular origin3-5. It suggests that DfMRI has an advantage to evaluate brain function regardless of the status of hemodynamic responses.
Here, by employing the technical advantages of DfMRI, we proposed the analytical framework to characterize dynamical brain function in neuropsychiatric disease model mice with an obsessive–compulsive disorder (OCD)-related behavior. The framework consists of three analyses; local activity analysis with water apparent diffusion coefficients (ADC) data, functional connectivity (FC) analysis with diffusion-weighted signal obtained with b=1800 s/mm² (Sb1800) data, and ignition-driven mean integration (IDMI) analysis with both data (Fig. 1).

Methods

Animal: We recruited astrocytic glutamate transporter (GLT1)-knockdown (KD) mice (GLASTCreERT2/+;GLT1 flox/flox, n=14)6. GLAST+/+;GLT1flox/flox mice were used as the control (ctl) mice (n=14). All animals were intraperitoneally injected with 100 mg/kg tamoxifen (TAM) at postnatal days 17-21 for consecutive five days to induce the deficit of GLT1.
MRI acquisition: MRI were conducted under awake condition at 11.7T equipped with cryoprobe and a gradient system allowing a maximum gradient strength of 1000 mT/m. (BioSpec, Bruker). A detail method of awake MRI was previously described7. DfMRI images were acquired with the following parameters: diffusion-sensitized double spine echo (SE)-echo planner imaging (EPI) sequence; TR/TE=2000/37 ms, Resolution=0.2×0.2×0.8 mm3, Slice number=10, Repetition number=150, Scan time=10 min, and b-values=1000 and 1800 mm2/s along 1 directions; [X=1,Y=1,Z=1].
Analysis: Pre-processing (realignment, slice timing, normalization, and smoothing) and statistics were performed using SPM12. ADCs were calculated as ADC=ln(Sb1000/Sb1800)/800 and compared a magnitude of ADC for local activity3. FC was calculated using Sb1800 by CONN. ROIs of 59 bilateral brain loci were defined from the Allen Mouse Brain Atlas. IDMI analysis were referred as the previous literatures8,9. Briefly, intrinsic ignition events for each given brain region were determined using time points of up-regulation triggers of neural activation (ADC decrease). Then we calculated the IDMI value from the event-related phase synchronized connectivity of Sb1800 data (See Fig. 5A and B). The IDMI value reflects amounts of information propagation from local activity defined by the driving events to global network.

Results and Discussion

ADC values were used as a quantitative marker of neuronal activity without hemodynamic confound3,10. We compared ADC magnitude as local activity between the ctl and KD mice. Voxel based and ROI-based analysis found that ADCs were significantly decreased in motor cortex (MO), somatosensory cortex (SS), prelimbic area (PL), infralimbic area (ILA), orbitofrontal area (ORB), angular insular cortex (AI), Caudate-Putamen (CP), nuclear accumbens (ACB), ventral globus pallidus (PALv), medial group of the dorsal thalamus (MED) in the KD mice (Fig. 2). According to previous evidence that the ADC decrease reflected an increase of neural activity3, this result indicated that the KD mice exhibited hyperactivation of the cortico-striato-thalamic (CST) circuit. The previous electrophysiological study using the same KD mice confirmed hyperactivation at the striatum6.
Next, we calculated FC as a pairwise Peason’s correlation of time-series Sb1800 between brain nodes, and then generated FC maps and matrix with the seeded bilateral ROIs (Fig. 3A and B). We found FC alterations in the CST circuit (Fig. 3C). From these results, local activity and FC analysis confirmed the hyperactivation and abnormal connectivity of CST circuit in the KD mice, as suggested from previous findings in OCD patients11-13.
Finally, we performed IDMI analysis to elucidate the relationship between hyperactivation and abnormal FC patterns. As a result, we found the IDMI values were significantly increased at MO, SS, ORB, and RSP in the KD mice while they were decreased at VIS, ACA, PL, ILA, CP, ACB, PALd, PALv, MED (Fig. 5C). The IDMI increase reflected an increase in neural propagation from local activity of the given brain region to global network. In addition, we investigated relationship between the IDMI value and ADC or FC (Fig. 5). In ORB and SS, the ADC decrease was related with the IDMI increase and the enhancement of their connectivity, On the other hand, in PL and dorsal globus pallidus (PALd), the ADC decrease was related with the IDMI decrease and the weakening of their connectivity. These results indicated that the hyperactivation of the CST circuit was associated with the abnormal connectivity and neural propagation.

Conclusion

We proposed DfMRI analytical framework and characterized dynamical brain function in the neuropsychiatric conditions. The hemodynamic-independent measurement by DfMRI is a key to succeed in detecting these brain function. We anticipate that our DfMRI approach is suitable to measure brain function in conditions of neurovascular coupling impairment. Applying our approach to clinical studies with geriatric patients who commonly comorbid dysfunction of vascular response will be a worthwhile challenge.

Acknowledgements

This work was supported by a Grant-in-Aid for Research Fellowships of the Japan Society for the Promotion of Science (JSPS Research Fellow) under grant number 18J00922 and Early-Career Scientists of JSPS KAKENHI under grant number 19K16298.

References

1. Shih, Y. Y. et al. A new scenario for negative functional magnetic resonance imaging signals: endogenous neurotransmission. J Neurosci 29, 3036-3044 (2009).

2. Veldsman, M., Cumming, T. & Brodtmann, A. Beyond BOLD: optimizing functional imaging in stroke populations. Hum Brain Mapp 36, 1620-1636 (2015).

3. Abe, Y., Tsurugizawa, T. & Le Bihan, D. Water diffusion closely reveals neural activity status in rat brain loci affected by anesthesia. PLoS Biol 15, e2001494, (2017).

4. Le Bihan, D., Urayama, S., Aso, T., Hanakawa, T. & Fukuyama, H. Direct and fast detection of neuronal activation in the human brain with diffusion MRI. Proc Natl Acad Sci U S A. 103, 8263-8268 (2006).

5. Tsurugizawa, T., Ciobanu, L. & Le Bihan, D. Water diffusion in brain cortex closely tracks underlying neuronal activity. Proc Natl Acad Sci U S A 110, 11636-11641 (2013).

6. Aida, T. et al. Astroglial glutamate transporter deficiency increases synaptic excitability and leads to pathological repetitive behaviors in mice. Neuropsychopharmacology 40, 1569-1579 (2015).

7. Yoshida, K. et al. Physiological effects of a habituation procedure for functional MRI in awake mice using a cryogenic radiofrequency probe. J Neurosci Methods 274, 38-48 (2016).

8. Deco, G., Tagliazucchi, E., Laufs, H., Sanjuan, A. & Kringelbach, M. L. Novel Intrinsic Ignition Method Measuring Local-Global Integration Characterizes Wakefulness and Deep Sleep. eNeuro 4 (2017).

9. Deco, G. & Kringelbach, M. L. Hierarchy of Information Processing in the Brain: A Novel 'Intrinsic Ignition' Framework. Neuron 94, 961-968 (2017).

10. Abe, Y., Van Nguyen, K., Tsurugizawa, T., Ciobanu, L. & Le Bihan, D. Modulation of water diffusion by activation-induced neural cell swelling in Aplysia Californica. Sci Rep 7 (2017).

11. Hou, J. M. et al. Resting-state functional connectivity abnormalities in patients with obsessive-compulsive disorder and their healthy first-degree relatives. J Psychiatry Neurosci 39, 304-311 (2014).

12. Jung, W. H. et al. Altered functional network architecture in orbitofronto-striato-thalamic circuit of unmedicated patients with obsessive-compulsive disorder. Hum Brain Mapp 38 (2017).

13. Yang, X. et al. Intrinsic brain abnormalities in drug-naive patients with obsessive-compulsive disorder: A resting-state functional MRI study. J Affect Disord 245, 861-868 (2019).

Figures

A schematic diagram of the DfMRI analytical framework. Two type of time-series signals of DfMRI with b=1000 and 1800 mm 2/s (Sb1000 and Sb1800, respectively) were obtained from the GLT1-KD mice and control mice in resting state and awake condition. Time-series ADCs were computed with Sb1000 and Sb1800. ADC was used for investigating local activity and defining the driving events ignited by neural excitation in IDMI. Sb1800 was used for calculating functional connectivity (FC) and phase connectivity in IDMI.

ADC-DfMRI detected hyperexcitation of the CST circuit in KD mice. (A) A voxel-based comparison of ADC between the ctl and KD mice. Blue color represents a significant ADC decreases coupled with an increase in neural activity, while hot color represents the opposite (FWE corrected p<0.05). The numbers below the maps represent distances from the bregma. (B) ADC differences were compared by ROI-based analysis at each area of the CST circuit. * p<0.05, ** p<0.01 (Bonferroni-corrected p-value of Student’s t-test).

FC analysis found biased connectivity across the CST circuit. (A1-9) Representative averaged FC maps with 9 seeded ROIs (green regions) were shown. All connectivity was calculated with only ctl mice. (B) FC matrix of each connectivity over the CST circuit were shown at below (ctl) and upper (KD). (C) Comparisons of z-scored FC values of specific pathway of the CST circuit between the ctl and KD mice. * p<0.05, ** p<0.01 (Bonferroni-corrected p-value of Student’s t-test).

IDMI revealed alterations of local activity-initiated propagation across the CST circuit. (A) A scheme of IDMI analysis as an example of a given brain region to estimating the driving event. (B) Phase lock connectivity was computed from time window (gray). From the binarized matrix, IDMI was calculated by finding the largest subcomponent. (C) IDMI of each given region were compared between ctl and KD groups. Green color IDMI was calculated from shuffled Sb1800 data to validate the empirical data. * p<0.05 (FDR p-value, ctl vs KD). # p<0.001 (two-way ANOVA, empirical vs shuffle).

Local activity-initiated propagation had relationships with local activity and functional connectivity. (A-D) Correlations between IDMI and ADC (left panels), IDMI and FC of the given pathway (middle panels), and ADC and FC (right panels) at ORB (A), SS (B), CP (C), and PALd (D). ORB and SS showed that the IDMI increases were correlated with the increases in local activity (ADC decrease) and the enhanced connectivity. On the other hand, CP and PALd showed that the IDMI decreases were correlated with the increases in local activity and the weaken connectivity.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
1343