1717

Reorganization of Rodent Resting-state Functional Networks after mild Traumatic Brain Injury
Chia-Feng Lu1,2,3, Yu-Chieh Jill Kao1,2, Huai-Lu Chen1,4, Fei-Ting Hsu1,5, Ping-Huei Tsai1,2,5, Hua-Shan Liu1,6, Li-Chun Hsieh1,5, Gilbert Aaron Lee1,4, and Cheng-Yu Chen1,2,5

1Translational Imaging Research Center, College of Medicine, Taipei Medical University, Taipei, Taiwan, 2Department of Radiology, School of Medicine, College of Medicine, Taipei Medical University, Taipei, Taiwan, 3Department of Biomedical Imaging and Radiological Sciences, National Yang-Ming University, Taipei, Taiwan, 4Department of Medical Research, Taipei Medical University Hospital, Taipei, Taiwan, 5Department of Medical Imaging, Taipei Medical University Hospital, Taipei, Taiwan, 6School of Biomedical Engineering, College of Biomedical Engineering, Taipei Medical University, Taipei, Taiwan

Synopsis

Complex network analysis unraveled the mTBI-related functional reorganization with the disruption of long-distance connections of thalamocortical interactions and enhanced local connectivity segregation in the regions of neocortex and thalamic nuclei.

Background and Purpose

Mild traumatic brain injury (mTBI) is a complex disease with temporary or persistent cognitive and behavioral disturbance while the gross brain structures appear to be normal in the early phase. Recent studies revealed the abnormalities of functional connectivity after mTBI based on the resting-state fMRI (rsfMRI)1-3. However, it is still unclear how the intricate functional networks reorganize in response to the brain injury, and how the injury-induced functional reorganization can be correlated to the extent of mTBI. In this study, an experimental mTBI rodent model was employed to investigate the functional reorganization after injury. We anticipated the occurrence of altered functional connectivity related to the impact site, and aimed to investigate the corresponding network reorganization and its relations to the extent of injury. Exploratory network analysis based on the graph theory was employed to globally (whole-brain network) and locally (region-wise) assess the characteristics of functional network measured by the blood oxygenation level dependent (BOLD) rsfMRI.

Materials and Methods

The animal study was approved by the Institutional Animal Care and Use Committee with the agreement of the 3Rs principle. Male Sprague–Dawley rats (weighing 250–400g) were anesthetized and placed in a stereotactic frame. An mTBI model modified from weight-drop impact acceleration injury4 with a weight of 600 gw dropping through a 1m-height tube into the metal disc cemented on top of the skull of the left primary somatosensory cortex was used. Two experimental conditions, a single impact (6 rats) and repetitive impacts with an interval of 1 hour (6 rats), was employed to create different mTBI extents. Animal MRI data were acquired using a Brucker 7T PharmaScan scanner with a circular surface coil for signal receiving. A RARE sequence (TR/TE=2650/40 ms; voxel size= 0.13x0.13x1.00 mm3) was performed to acquire anatomical images. BOLD rsfMRI was acquired using a single-shot GE-EPI (TR/TE=1000/15 ms; voxel size= 0.55x0.55x1.00 mm3; 300 volumes) under a combined anesthetic regime of low-dose isoflurane and dexmedetomidine infusion5. Physiological parameters were continuously monitored and maintained within normal ranges during the experiment. In total, data at pre-injury (n=12) and 24 h after injury (each n=6 after single and repetitive injury) were used in this study.

fMRI preprocessing steps included the correction of slice timing, realignment, co-registration between subjects, and spatial smoothing using SPM8 (http://www.fil.ion.ucl.ac.uk/spm/). The gray matter was partitioned into 3 categories, including neocortex (NE, 10 regions), subcortical regions (SC, 12 regions), and thalamic nuclei (TN, 10 regions) for each hemisphere based on Waxholm Space atlas (Table 1)6,7. The Pearson cross-correlation coefficients were computed between regional BOLD signals (between 0.01-0.3 Hz)8 with regression of motion parameters followed by Fisher’s r-to-z transform. Topological parameters based on the graph theory were then applied to characterize functional networks9. One-sample t-test was employed to determine the significance of functional connectivity within a condition (p<0.05, with FDR correction), and two-sample t-test was used to determine the differences between single-/repetitive-impact and pre-injury condition (p<0.05).

Results and Discussion

The applied animal model successfully mimicked the radiological traits of mTBI, i.e., no skull fracture and intact brain structure on RARE T2-weighted images. Persistent declines in neurological and general locomotor performance were observed in animals after injury. Our results of network analysis also revealed the alterations of functional organization. Substantial functional connections, with mostly positive connectivity for adjacent regions and negative values for long-distance connections between NC, SC, and TN, were presented on the connectivity matrix before injury (Fig.1A). When injury extent increased, the functional organization tended to lose the long-distance connections, especially between NC and TN (thalamocortical circuits) and elevated the connectivity strength within categories (Fig.1B and 1C). We further investigated the transitions of network hubs (pivotal regions with 15% top rank of connectivity strength) between conditions. Before injury, the hubs were located at bilateral somatosensory, motor, striatum, and basal forebrain regions (Fig.2A) and transferred to the TN regions, including sensorimotor-associated nuclei (VPM, VPL), mediodorsal nuclei (MD) and posterior nuclear group (Po) (Fig.2C and 2E). The global topological parameters of functional networks confirmed an impact-related increment of connectivity strength and clustering coefficient (local segregation) while preserving the network efficiency (Fig.3). Finally, the local topological parameters further revealed a significantly decreased degree (number of functionally connected neighbors) and efficiency at the impacted left somatosensory region (Fig.4). Several TN regions showed contrary increases of degree, larger increments of strength and clustering coefficient (Fig.4). Our results demonstrated a compensatory mechanism in mTBI involving the regulation of thalamocortical circuits presented by a functional reorganization after injury.

Acknowledgements

This study was funded in part by the Taipei Medical University (TMU103-AE1-B20) and the Ministry of Science and Technology (MOST 105-2314-B-038-014, MOST 104-2923-B-038-003-MY3, MOST 105-2628-B-038-002-MY2), Taipei, Taiwan.

References

1. Mayer AR, Mannell MV, Ling J, Gasparovic C, Yeo RA. Functional connectivity in mild traumatic brain injury. Human brain mapping. 2011;32(11):1825-35.

2. Shumskaya E, Andriessen TM, Norris DG, Vos PE. Abnormal whole-brain functional networks in homogeneous acute mild traumatic brain injury. Neurology. 2012;79(2):175-82.

3. Zhou Y, Milham MP, Lui YW, Miles L, Reaume J, Sodickson DK, Grossman RI, Ge Y. Default-mode network disruption in mild traumatic brain injury. Radiology. 2012;265(3):882-92.

4. Marmarou A, Foda MA, Brink WV, Campbell J, Kita H, Demetriadou K. A new model of diffuse brain injury in rats: Part I: Pathophysiology and biomechanics. Journal of neurosurgery. 1994;80(2):291-300.

5. Lu H, Zou Q, Gu H, Raichle ME, Stein EA, Yang Y. Rat brains also have a default mode network. Proceedings of the National Academy of Sciences. 2012;109(10):3979-84.

6. Papp EA, Leergaard TB, Calabrese E, Johnson GA, Bjaalie JG. Waxholm Space atlas of the Sprague Dawley rat brain. NeuroImage. 2014;97:374-86.

7. Paxinos G, Watson C. The Rat Brain in Stereotaxic Coordinates, 6th ed. Academic Press, 2007.

8. Grandjean J, Schroeter A, Batata I, Rudin M. Optimization of anesthesia protocol for resting-state fMRI in mice based on differential effects of anesthetics on functional connectivity patterns. Neuroimage. 2014;102:838-47.

9. Lu CF, Soong BW, Wu HM, Teng S, Wang PS, Wu YT. Disrupted cerebellar connectivity reduces whole-brain network efficiency in multiple system atrophy. Movement Disorders. 2013;28(3):362-9.

Figures

Table 1 The employed parcellation of gray matter into 32 regions for each hemisphere.

Figure 1 The mean 64 x 64 connectivity matrices (constructed by the correlations between 64 brain regions) for the pre-injury (A), single-impact (B), and repetitive-impact (C) conditions. The connectivity strength is presented by Fisher’s z values, and only the significant connectivity (one-sample t-test, p<0.05 with FDR correction) for each condition was displayed in the figures. NC: neocortical regions; SC: subcortical regions; TN: thalamic nuclei.

Figure 2 The topographic maps of functional connectivity in the conditions of pre-injury (A and B), single-impact (C and D), and repetitive-impact (E and F) status. The impact sites, i.e., on top of the left primary sensorimotor cortex, are labeled by the gray columns. The network hubs (pivotal nodes) are labeled with red. A, B, and C subplots show that the positive connectivity is associated with the bilateral connections and within NC, SC, or TN connections. Whereas B, D, and F subplots present the distribution of negative connectivity is related to the long-distance connections cross three categories.

Figure 3 The global topological parameters in characterizing the functional organization. An asterisk (*) represents significant difference (p<0.05) between pre-injury and repetitive-impact conditions, and a pound sign (#) represents significant difference (p<0.05) between pre-injury and single-impact conditions

Figure 4 The local topological parameters, including degree (number of functionally connected neighbors), strength (average absolute z values), clustering coefficient (local segregation), and network efficiency. Only the properties of neocortical regions and thalamic nuclei are displayed. Green and blue asterisks represent significant difference (p<0.05) between pre-injury and repetitive-impact conditions and between pre-injury and single-impact conditions, respectively.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
1717