ComBating inter-site differences in field strength: harmonizing preclinical traumatic brain injury MRI data

Integrating datasets from multiple sites and scanners can increase statistical power for neuroimaging studies but can also introduce significant inter-site confounds. We evaluated the effectiveness of ComBat, an empirical Bayes approach, to combine longitudinal preclinical MRI data acquired at 4.7 or 9.4 T at two different sites in Australia. Male Sprague Dawley rats underwent MRI on Days 2, 9, 28, and 150 following moderate/severe traumatic brain injury (TBI) or sham injury as part of Project 1 of the NIH/NINDS-funded Centre Without Walls EpiBioS4Rx project. Diffusion-weighted and multiple-gradient-echo images were acquired, and outcomes included QSM, FA, and ADC. Acute injury measures including apnea and self-righting reflex were consistent between sites. Mixed-effect analysis of ipsilateral and contralateral corpus callosum (CC) summary values revealed a significant effect of site on FA and ADC values, which was removed following ComBat harmonization. Bland-Altman plots for each metric showed reduced variability across sites following ComBat harmonization, including for QSM, despite appearing to be largely unaffected by inter-site differences and no effect of site observed. Following harmonization, the combined inter-site data revealed significant differences in the imaging metrics consistent with previously reported outcomes. TBI resulted in significantly reduced FA and increased susceptibility in the ipsilateral CC, and significantly reduced FA in the contralateral CC compared with sham-injured rats. Additionally, TBI rats also exhibited a reversal in ipsilateral CC ADC values over time with significantly reduced ADC at Day 9, followed by increased ADC 150 days after injury. Our findings demonstrate the need for harmonizing multi-site preclinical MRI data and show that this can be successfully achieved using ComBat while preserving phenotypical changes due to TBI.

Epilepsy phenotype and its reproducibility after lateral fluid percussion-induced traumatic brain injury in rats: Multicenter EpiBioS4Rx study project 1

OBJECTIVE

This study was undertaken to assess reproducibility of the epilepsy outcome and phenotype in a lateral fluid percussion model of posttraumatic epilepsy (PTE) across three study sites.

METHODS

A total of 525 adult male Sprague Dawley rats were randomized to lateral fluid percussion-induced brain injury (FPI) or sham operation. Of these, 264 were assigned to magnetic resonance imaging (MRI cohort, 43 sham, 221 traumatic brain injury [TBI]) and 261 to electrophysiological follow-up (EEG cohort, 41 sham, 220 TBI). A major effort was made to harmonize the rats, materials, equipment, procedures, and monitoring systems. On the 7th post-TBI month, rats were video-EEG monitored for epilepsy diagnosis.

RESULTS

A total of 245 rats were video-EEG phenotyped for epilepsy on the 7th postinjury month (121 in MRI cohort, 124 in EEG cohort). In the whole cohort (n = 245), the prevalence of PTE in rats with TBI was 22%, being 27% in the MRI and 18% in the EEG cohort (p > .05). Prevalence of PTE did not differ between the three study sites (p > .05). The average seizure frequency was .317 ± .725 seizures/day at University of Eastern Finland (UEF; Finland), .085 ± .067 at Monash University (Monash; Australia), and .299 ± .266 at University of California, Los Angeles (UCLA; USA; p < .01 as compared to Monash). The average seizure duration did not differ between UEF (104 ± 48 s), Monash (90 ± 33 s), and UCLA (105 ± 473 s; p > .05). Of the 219 seizures, 53% occurred as part of a seizure cluster (≥3 seizures/24 h; p >.05 between the study sites). Of the 209 seizures, 56% occurred during lights-on period and 44% during lights-off period (p > .05 between the study sites).

SIGNIFICANCE

The PTE phenotype induced by lateral FPI is reproducible in a multicenter design. Our study supports the feasibility of performing preclinical multicenter trials in PTE to increase statistical power and experimental rigor to produce clinically translatable data to combat epileptogenesis after TBI.

Successful harmonization in EpiBioS4Rx biomarker study on post-traumatic epilepsy paves the way towards powered preclinical multicenter studies

OBJECTIVE

Project 1 of the Preclinical Multicenter Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx) consortium aims to identify preclinical biomarkers for antiepileptogenic therapies following traumatic brain injury (TBI). The international participating centers in Finland, Australia, and the United States have made a concerted effort to ensure protocol harmonization. Here, we evaluate the success of harmonization process by assessing the timing, coverage, and performance between the study sites.

METHOD

We collected data on animal housing conditions, lateral fluid-percussion injury model production, postoperative care, mortality, post-TBI physiological monitoring, timing of blood sampling and quality, MR imaging timing and protocols, and duration of video-electroencephalography (EEG) follow-up using common data elements. Learning effect in harmonization was assessed by comparing procedural accuracy between the early and late stages of the project.

RESULTS

The animal housing conditions were comparable between the study sites but the postoperative care procedures varied. Impact pressure, duration of apnea, righting reflex, and acute mortality differed between the study sites (p < 0.001). The severity of TBI on D2 post TBI assessed using the composite neuroscore test was similar between the sites, but recovery of acute somato-motor deficits varied (p < 0.001). A total of 99% of rats included in the final cohort in UEF, 100% in Monash, and 79% in UCLA had blood samples taken at all time points. The timing of sampling differed on day (D)2 (p < 0.05) but not D9 (p > 0.05). Plasma quality was poor in 4% of the samples in UEF, 1% in Monash and 14% in UCLA. More than 97% of the final cohort were MR imaged at all timepoints in all study sites. The timing of imaging did not differ on D2 and D9 (p > 0.05), but varied at D30, 5 months, and ex vivo timepoints (p < 0.001). The percentage of rats that completed the monthly high-density video-EEG follow-up and the duration of video-EEG recording on the 7th post-injury month used for seizure detection for diagnosis of post-traumatic epilepsy differed between the sites (p < 0.001), yet the prevalence of PTE (UEF 21%, Monash 22%, UCLA 23%) was comparable between the sites (p > 0.05). A decrease in acute mortality and increase in plasma quality across time reflected a learning effect in the TBI production and blood sampling protocols.

SIGNIFICANCE

Our study is the first demonstration of the feasibility of protocol harmonization for performing powered preclinical multi-center trials for biomarker and therapy discovery of post-traumatic epilepsy.

Image data harmonization tools for the analysis of post-traumatic epilepsy development in preclinical multisite MRI studies

Preclinical MRI studies have been utilized for the discovery of biomarkers that predict post-traumatic epilepsy (PTE). However, these single site studies often lack statistical power due to limited and homogeneous datasets. Therefore, multisite studies, such as the Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx), are developed to create large, heterogeneous datasets that can lead to more statistically significant results. EpiBioS4Rx collects preclinical data internationally across sites, including the United States, Finland, and Australia. However, in doing so, there are robust normalization and harmonization processes that are required to obtain statistically significant and generalizable results. This work describes the tools and procedures used to harmonize multisite, multimodal preclinical imaging data acquired by EpiBioS4Rx. There were four main harmonization processes that were utilized, including file format harmonization, naming convention harmonization, image coordinate system harmonization, and diffusion tensor imaging (DTI) metrics harmonization. By using Python tools and bash scripts, the file formats, file names, and image coordinate systems are harmonized across all the sites. To harmonize DTI metrics, values are estimated for each voxel in an image to generate a histogram representing the whole image. Then, the Quantitative Imaging Toolkit (QIT) modules are utilized to scale the mode to a value of one and depict the subsequent harmonized histogram. The standardization of file formats, naming conventions, coordinate systems, and DTI metrics are qualitatively assessed. The histograms of the DTI metrics were generated for all the individual rodents per site. For inter-site analysis, an average of the individual scans was calculated to create a histogram that represents each site. In order to ensure the analysis can be run at the level of individual animals, the sham and TBI cohort were analyzed separately, which depicted the same harmonization factor. The results demonstrate that these processes qualitatively standardize the file formats, naming conventions, coordinate systems, and DTI metrics of the data. This assists in the ability to share data across the study, as well as disseminate tools that can help other researchers to strengthen the statistical power of their studies and analyze data more cohesively.

Orientation selective stimulation with tetrahedral electrodes of the rat infralimbic cortex to indirectly target the amygdala

Introduction

Deep brain stimulation (DBS) is a rapidly developing therapeutic intervention with constantly expanding neurological and psychiatric indications. A major challenge for the approach is the precise targeting and limitation of the effect on the desired neural pathways. We have introduced a new approach, orientation selective stimulation (OSS) that allows free rotation of the induced electric field on a plane when using a probe with three parallel electrodes forming an equilateral triangle at the tip. Here, we expand the technique by introducing a tetrahedral stimulation probe that enables adjustment of the primary electric field direction freely at any angle in a 3D space around the stimulating probe. OSS in 3D will enable better targeting of the electric field according to the local brain anatomy. We tested its utility in a rat model of DBS for treatment-resistant depression. The stimulation directed to the subgenual anterior cingulate cortex (sgACC) has yielded dramatic improvement in individual patients suffering from therapy resistant depression, but no consistent benefit in larger series. This failure has been ascribed to the challenging anatomy of sgACC with several crossing neural tracts and individual differences in the local anatomy.

Methods

We stimulated infralimbic cortex (IL), the rat analog of sgACC, and recorded local electrical responses in amygdala (AMG) that is monosynaptically connected to IL and plays a central role in emotional states. We further traced AMG-IL connections using a viral vector and tractography using diffusion magnetic resonance imaging (MRI). Finally, we mimicked the clinical situation by delivering sustained 130 Hz stimulation at IL at the most effective field orientation and followed changes in resting-state functional connectivity with IL using functional MRI. To help interpretation of responses in functional connectivity, we stimulated only the left IL, which we did not expect to evoke measurable changes in the rat behavior.

Results

The AMG evoked responses depended systematically on the IL stimulation field orientation and yielded the maximum response in near vertical field orientation in accordance with tractography. Sustained 130 Hz stimulation at a field orientation yielding the strongest AMG evoked responses increased functional connectivity between IL and AMG on the stimulation side.

Conclusion

These findings suggest that OSS in 3D provides a new approach to optimize the DBS for every individual patient with a single stimulation probe implantation.

A consensus protocol for functional connectivity analysis in the rat brain

Task-free functional connectivity in animal models provides an experimental framework to examine connectivity phenomena under controlled conditions and allows for comparisons with data modalities collected under invasive or terminal procedures. Currently, animal acquisitions are performed with varying protocols and analyses that hamper result comparison and integration. Here we introduce StandardRat, a consensus rat functional magnetic resonance imaging acquisition protocol tested across 20 centers. To develop this protocol with optimized acquisition and processing parameters, we initially aggregated 65 functional imaging datasets acquired from rats across 46 centers. We developed a reproducible pipeline for analyzing rat data acquired with diverse protocols and determined experimental and processing parameters associated with the robust detection of functional connectivity across centers. We show that the standardized protocol enhances biologically plausible functional connectivity patterns relative to previous acquisitions. The protocol and processing pipeline described here is openly shared with the neuroimaging community to promote interoperability and cooperation toward tackling the most important challenges in neuroscience.

Diffusion MRI approaches for investigating microstructural complexity in a rat model of traumatic brain injury

Our study explores the potential of conventional and advanced diffusion MRI techniques including diffusion tensor imaging (DTI), and single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD) to investigate complex microstructural changes following severe traumatic brain injury in rats at a chronic phase. Rat brains after sham-operation or lateral fluid percussion (LFP) injury were scanned ex vivo in a 9.4 T scanner. Our region-of-interest-based approach of tensor-, and SS3T-CSD derived fixel-, 3-tissue signal fraction maps were sensitive to changes in both white matter (WM) and grey matter (GM) areas. Tensor-based measures, such as fractional anisotropy (FA) and radial diffusivity (RD), detected more changes in WM and GM areas as compared to fixel-based measures including apparent fiber density (AFD), peak FOD amplitude and primary fiber bundle density, while 3-tissue signal fraction maps revealed distinct changes in WM, GM, and phosphate-buffered saline (PBS) fractions highlighting the complex tissue microstructural alterations post-trauma. Track-weighted imaging demonstrated changes in track morphology including reduced curvature and average pathlength distal from the primary lesion in severe TBI rats. In histological analysis, changes in the diffusion MRI measures could be associated to decreased myelin density, loss of myelinated axons, and increased cellularity, revealing progressive microstructural alterations in these brain areas five months after injury. Overall, this study highlights the use of combined conventional and advanced diffusion MRI measures to obtain more precise insights into the complex tissue microstructural alterations in chronic phase of severe brain injury.

Plasma Neurofilament Light Chain (NF-L) Is a Prognostic Biomarker for Cortical Damage Evolution but Not for Cognitive Impairment or Epileptogenesis Following Experimental TBI

Plasma neurofilament light chain (NF-L) levels were assessed as a diagnostic biomarker for traumatic brain injury (TBI) and as a prognostic biomarker for somatomotor recovery, cognitive decline, and epileptogenesis. Rats with severe TBI induced by lateral fluid-percussion injury (n = 26, 13 with and 13 without epilepsy) or sham-operation (n = 8) were studied. During a 6-month follow-up, rats underwent magnetic resonance imaging (MRI) (day (D) 2, D7, and D21), composite neuroscore (D2, D6, and D14), Morris-water maze (D35−D39), and a 1-month-long video-electroencephalogram to detect unprovoked seizures during the 6th month. Plasma NF-L levels were assessed using a single-molecule assay at baseline (i.e., naïve animals) and on D2, D9, and D178 after TBI or a sham operation. Plasma NF-L levels were 483-fold higher on D2 (5072.0 ± 2007.0 pg/mL), 89-fold higher on D9 (930.3 ± 306.4 pg/mL), and 3-fold higher on D176 32.2 ± 8.9 pg/mL after TBI compared with baseline (10.5 ± 2.6 pg/mL; all p < 0.001). Plasma NF-L levels distinguished TBI rats from naïve animals at all time-points examined (area under the curve [AUC] 1.0, p < 0.001), and from sham-operated controls on D2 (AUC 1.0, p < 0.001). Plasma NF-L increases on D2 were associated with somatomotor impairment severity (ρ = −0.480, p < 0.05) and the cortical lesion extent in MRI (ρ = 0.401, p < 0.05). Plasma NF-L increases on D2 or D9 were associated with the cortical lesion extent in histologic sections at 6 months post-injury (ρ = 0.437 for D2; ρ = 0.393 for D9, p < 0.05). Plasma NF-L levels, however, did not predict somatomotor recovery, cognitive decline, or epileptogenesis (p > 0.05). Plasma NF-L levels represent a promising noninvasive translational diagnostic biomarker for acute TBI and a prognostic biomarker for post-injury somatomotor impairment and long-term structural brain damage.

Acute Hippocampal Damage as a Prognostic Biomarker for Cognitive Decline but Not for Epileptogenesis after Experimental Traumatic Brain Injury

It is necessary to develop reliable biomarkers for epileptogenesis and cognitive impairment after traumatic brain injury when searching for novel antiepileptogenic and cognition-enhancing treatments. We hypothesized that a multiparametric magnetic resonance imaging (MRI) analysis along the septotemporal hippocampal axis could predict the development of post-traumatic epilepsy and cognitive impairment. We performed quantitative T2 and T2* MRIs at 2, 7 and 21 days, and diffusion tensor imaging at 7 and 21 days after lateral fluid-percussion injury in male rats. Morris water maze tests conducted between 35-39 days post-injury were used to diagnose cognitive impairment. One-month-long continuous video-electroencephalography monitoring during the 6th post-injury month was used to diagnose epilepsy. Single-parameter and regularized multiple linear regression models were able to differentiate between sham-operated and brain-injured rats. In the ipsilateral hippocampus, differentiation between the groups was achieved at most septotemporal locations (cross-validated area under the receiver operating characteristic curve (AUC) 1.0, 95% confidence interval 1.0-1.0). In the contralateral hippocampus, the highest differentiation was evident in the septal pole (AUC 0.92, 95% confidence interval 0.82-0.97). Logistic regression analysis of parameters imaged at 3.4 mm from the contralateral hippocampus's temporal end differentiated between the cognitively impaired rats and normal rats (AUC 0.72, 95% confidence interval 0.55-0.84). Neither single nor multiparametric approaches could identify the rats that would develop post-traumatic epilepsy. Multiparametric MRI analysis of the hippocampus can be used to identify cognitive impairment after an experimental traumatic brain injury. This information can be used to select subjects for preclinical trials of cognition-improving interventions.

A companion to the preclinical common data elements and case report forms for in vivo rodent neuroimaging: A report of the TASK3-WG3 Neuroimaging Working Group of the ILAE/AES Joint Translational Task Force

The International League Against Epilepsy/American Epilepsy Society (ILAE/AES) Joint Translational Task Force established the TASK3 working groups to create common data elements (CDEs) for various aspects of preclinical epilepsy research studies, which could help improve the standardization of experimental designs. In this article, we discuss CDEs for neuroimaging data that are collected in rodent models of epilepsy, with a focus on adult rats and mice. We provide detailed CDE tables and case report forms (CRFs), and with this companion manuscript, we discuss the methodologies for several imaging modalities and the parameters that can be collected.

Histopathological modeling of status epilepticus-induced brain damage based on in vivo diffusion tensor imaging in rats

Non-invasive magnetic resonance imaging (MRI) methods have proved useful in the diagnosis and prognosis of neurodegenerative diseases. However, the interpretation of imaging outcomes in terms of tissue pathology is still challenging. This study goes beyond the current interpretation of in vivo diffusion tensor imaging (DTI) by constructing multivariate models of quantitative tissue microstructure in status epilepticus (SE)-induced brain damage. We performed in vivo DTI and histology in rats at 79 days after SE and control animals. The analyses focused on the corpus callosum, hippocampal subfield CA3b, and layers V and VI of the parietal cortex. Comparison between control and SE rats indicated that a combination of microstructural tissue changes occurring after SE, such as cellularity, organization of myelinated axons, and/or morphology of astrocytes, affect DTI parameters. Subsequently, we constructed a multivariate regression model for explaining and predicting histological parameters based on DTI. The model revealed that DTI predicted well the organization of myelinated axons (cross-validated R = 0.876) and astrocyte processes (cross-validated R = 0.909) and possessed a predictive value for cell density (CD) (cross-validated R = 0.489). However, the morphology of astrocytes (cross-validated R &gt; 0.05) was not well predicted. The inclusion of parameters from CA3b was necessary for modeling histopathology. Moreover, the multivariate DTI model explained better histological parameters than any univariate model. In conclusion, we demonstrate that combining several analytical and statistical tools can help interpret imaging outcomes to microstructural tissue changes, opening new avenues to improve the non-invasive diagnosis and prognosis of brain tissue damage.

Light sedation with short habituation time for large-scale functional magnetic resonance imaging studies in rats

Traditionally, preclinical resting state functional magnetic resonance imaging (fMRI) studies have been performed in anesthetized animals. Nevertheless, as anesthesia affects the functional connectivity (FC) in the brain, there has been a growing interest in imaging in the awake state. Obviously, awake imaging requires resource- and time-consuming habituation prior to data acquisition to reduce the stress and motion of the animals. Light sedation has been a less widely exploited alternative for awake imaging, requiring shorter habituation times, while still reducing the effect of anesthesia. Here, we imaged 102 rats under light sedation and 10 awake animals to conduct an FC analysis. We established an automated data-processing pipeline suitable for both groups. Additionally, the same pipeline was used on data obtained from an openly available awake rat database (289 measurements in 90 rats). The FC pattern in the light sedation measurements closely resembled the corresponding patterns in both onsite and offsite awake datasets. However, fewer datasets had to be excluded due to movement in rats with light sedation. The temporal analysis of FC in the lightly sedated group indicated a lingering effect of anesthesia that stabilized after the first 5 min. In summary, our results indicate that the light sedation protocol is a valid alternative for large-scale studies where awake protocols may become prohibitively resource-demanding, as it provides similar results to awake imaging, preserves more scans, and requires shorter habituation times. The large amount of fMRI data obtained in this work are openly available for further analyses.

Event-recurring multiband SWIFT functional MRI with 200-ms temporal resolution during deep brain stimulation and isoflurane-induced burst suppression in rat

PURPOSE

To develop a high temporal resolution functional MRI method for tracking repeating events in the brain.

METHODS

We developed a novel functional MRI method using multiband sweep imaging with Fourier transformation (SWIFT), termed event-recurring SWIFT (EVER-SWIFT). The method is able to image similar repeating events with subsecond temporal resolution. Here, we demonstrate the use of EVER-SWIFT for detecting functional MRI responses during deep brain stimulation of the medial septal nucleus and during spontaneous isoflurane-induced burst suppression in the rat brain at 9.4 T with 200-ms temporal resolution.

RESULTS

The EVER-SWIFT approach showed that the shapes and time-to-peak values of the response curves to deep brain stimulation significantly differed between downstream brain regions connected to the medial septal nucleus, resembling findings obtained with traditional 2-second temporal resolution. In contrast, EVER-SWIFT allowed for detailed temporal measurement of a spontaneous isoflurane-induced bursting activity pattern, which was not achieved with traditional temporal resolution.

CONCLUSION

The EVER-SWIFT technique enables subsecond 3D imaging of both stimulated and spontaneously recurring brain activities, and thus holds great potential for studying the mechanisms of neuromodulation and spontaneous brain activity.

Spatial signatures of anesthesia-induced burst-suppression differ between primates and rodents

During deep anesthesia, the electroencephalographic (EEG) signal of the brain alternates between bursts of activity and periods of relative silence (suppressions). The origin of burst-suppression and its distribution across the brain remain matters of debate. In this work, we used functional magnetic resonance imaging (fMRI) to map the brain areas involved in anesthesia-induced burst-suppression across four mammalian species: humans, long-tailed macaques, common marmosets, and rats. At first, we determined the fMRI signatures of burst-suppression in human EEG-fMRI data. Applying this method to animal fMRI datasets, we found distinct burst-suppression signatures in all species. The burst-suppression maps revealed a marked inter-species difference: in rats, the entire neocortex engaged in burst-suppression, while in primates most sensory areas were excluded—predominantly the primary visual cortex. We anticipate that the identified species-specific fMRI signatures and whole-brain maps will guide future targeted studies investigating the cellular and molecular mechanisms of burst-suppression in unconscious states.

The development of anesthesia was a significant advance in medicine. It allows individuals to undergo surgery without feeling pain or remembering the experience. But scientists still do not know how anesthesia works. During anesthesia, scientists have measured brain activity using electroencephalograms (EEG) and found that the brain appears to turn on and off. Comatose patients also have similar switches between bursts of electrical activity and periods of silence. This burst-suppression pattern may be related to unconsciousness.

But scientists still have many questions about how anesthesia causes burst-suppression. One challenge is that while an EEG can tell scientists when the brain turns on and off, it does not show exactly where this occurs. Another imaging method called functional Magnetic Resonance Imaging (fMRI) may fill this gap by allowing scientists to map where the brain activity occurs.

Sirmpilatze et al. have created detailed maps of burst-suppression in humans, primates, and rats under anesthesia by analyzing brain scans using fMRI. In rats, the entire outer layer or cortex of the brain underwent a synchronized pattern of burst-suppression. In humans and primates, areas of the brain like those responsible for eyesight did not follow the rest of the cortex in switching on and off.

The experiments reveal crucial differences in how rats and humans and other primates respond to anesthesia. The fMRI mapping technique Sirmpilatze et al. created may help scientists learn more about these differences and why some parts of human brains do not undergo burst-suppression. This may help scientists learn more about unconsciousness and help improve anesthesia or the care of comatose patients.

Orientation selective DBS of entorhinal cortex and medial septal nucleus modulates activity of rat brain areas involved in memory and cognition

The recently introduced orientation selective deep brain stimulation (OS-DBS) technique freely controls the direction of the electric field's spatial gradient by using multiple contacts with independent current sources within a multielectrode array. The goal of OS-DBS is to align the electrical field along the axonal track of interest passing through the stimulation site. Here we utilized OS-DBS with a planar 3-channel electrode for stimulating the rat entorhinal cortex (EC) and medial septal nucleus (MSN), two promising areas for DBS treatment of Alzheimer's disease. The brain responses to OS-DBS were monitored by whole brain functional magnetic resonance imaging (fMRI) at 9.4 T with Multi-Band Sweep Imaging with Fourier Transformation (MB-SWIFT). Varying the in-plane OS-DBS stimulation angle in the EC resulted in activity modulation of multiple downstream brain areas involved in memory and cognition. Contrary to that, no angle dependence of brain activations was observed when stimulating the MSN, consistent with predictions based on the electrode configuration and on the main axonal directions of the targets derived from diffusion MRI tractography and histology. We conclude that tuning the OS-DBS stimulation angle modulates the activation of brain areas relevant to Alzheimer's disease, thus holding great promise in the DBS treatment of the disease.

Chronic hypometabolism in striatum and hippocampal network after traumatic brain injury and their relation with memory impairment - [18F]-FDG-PET and MRI 4 months after fluid percussion injury in rat

Hippocampal and thalamo-cortico-striatal networks are critical for memory function as well as execution of a variety of learning strategies. In subjects with memory impairment as a sequel of traumatic brain injury (TBI), the contribution of late metabolic depression across these networks to memory deficit is poorly understood. We used [18F]-FDG-PET to measure chronic post-TBI glucose uptake in the striatum and connected brain areas (septal and temporal hippocampus, thalamus, entorhinal cortex, frontoparietal cortex and amygdala) in rats with lateral fluid-percussion injury (LFPI). Then we assessed a link between network hypometabolism and memory impairment. At 4 months post TBI, glucose uptake was decreased in ipsilateral striatum (10%, p = 0.027), frontoparietal cortex (17%, p = 0.00009), and hippocampus (22%, p = 0.027) as compared to sham operated controls. Thalamic uptake was 6% lower ipsilaterally than contralaterally, p = 0.00004). At 5 months, Morris water maze (MWM) showed memory impairment in 83% of the rats with TBI. The lower the hippocampal or striatal [18F]-FDG uptake, the poorer the MWM performance (hippocampus: r = -0.471, p < 0.05; striatum: r = -0.696, p < 0.001). Striatal [18F]-FDG-PET identified the injured animals with memory impairment with 100% specificity and sensitivity (AUC = 1.000, p = 0.009). Interestingly, the low striatal glucose uptake was a better diagnostic biomarker for memory impairment than the reduced hippocampal (AUC = 0.806, p = 0.112) or entorhinal (AUC = 0.528, p = 0.885) glucose uptake. The volumetric atrophy assessed in T2 weighted MRI or the gliotic area in Nissl staining did not correlate with glucose uptake. Arterial spin labeling did not indicate any reduction in the striatal blood flow. Our study suggests that TBI-induced chronic hypometabolism in striatum contributes to the cognitive deficits.

Hippocampal position and orientation as prognostic biomarkers for post-traumatic epileptogenesis - an experimental study in rat lateral fluid-percussion model

Objective

This study was undertaken to identify prognostic biomarkers for posttraumatic epileptogenesis derived from parameters related to the hippocampal position and orientation.

Methods

Data were derived from two preclinical magnetic resonance imaging (MRI) follow‐up studies: EPITARGET (156 rats) and Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx; University of Eastern Finland cohort, 43 rats). Epileptogenesis was induced with lateral fluid percussion‐induced traumatic brain injury (TBI) in adult male Sprague Dawley rats. In the EPITARGET cohort, T2∗‐weighted MRI was performed at 2, 7, and 21 days and in the EpiBioS4Rx cohort at 2, 9, and 30 days and 5 months post‐TBI. Both hippocampi were segmented using convolutional neural networks. The extracted segmentation mask was used for a geometric construction, extracting 39 parameters that described the position and orientation of the left and right hippocampus. In each cohort, we assessed the parameters as prognostic biomarkers for posttraumatic epilepsy (PTE) both individually, using repeated measures analysis of variance, and in combination, using random forest classifiers.

Results

The extracted parameters were highly effective in discriminating between sham‐operated and TBI rats in both the EPITARGET and EpiBioS4Rx cohorts at all timepoints (t; balanced accuracy > .9). The most discriminating parameter was the inclination of the hippocampus ipsilateral to the lesion at t = 2 days and the volumes at t ≥ 7 days after TBI. Furthermore, in the EpiBioS4Rx cohort, we could effectively discriminate epileptogenic from nonepileptogenic animals with a longer MRI follow‐up, at t = 150 days (area under the curve = .78, balanced accuracy = .80, p = .0050), based on the orientation of both hippocampi. We found that the ipsilateral hippocampus rotated outward on the horizontal plane, whereas the contralateral hippocampus rotated away from the vertical direction.

Significance

We demonstrate that assessment of TBI‐induced hippocampal deformation by clinically translatable MRI methodologies detects subjects with prior TBI as well as those at high risk of PTE, paving the way toward subject stratification for antiepileptogenesis studies.

Convolutional Neural Networks Enable Robust Automatic Segmentation of the Rat Hippocampus in MRI After Traumatic Brain Injury

Registration-based methods are commonly used in the automatic segmentation of magnetic resonance (MR) brain images. However, these methods are not robust to the presence of gross pathologies that can alter the brain anatomy and affect the alignment of the atlas image with the target image. In this work, we develop a robust algorithm, MU-Net-R, for automatic segmentation of the normal and injured rat hippocampus based on an ensemble of U-net-like Convolutional Neural Networks (CNNs). MU-Net-R was trained on manually segmented MR images of sham-operated rats and rats with traumatic brain injury (TBI) by lateral fluid percussion. The performance of MU-Net-R was quantitatively compared with methods based on single and multi-atlas registration using MR images from two large preclinical cohorts. Automatic segmentations using MU-Net-R and multi-atlas registration were of excellent quality, achieving cross-validated Dice scores above 0.90 despite the presence of brain lesions, atrophy, and ventricular enlargement. In contrast, the performance of single-atlas segmentation was unsatisfactory (cross-validated Dice scores below 0.85). Interestingly, the registration-based methods were better at segmenting the contralateral than the ipsilateral hippocampus, whereas MU-Net-R segmented the contralateral and ipsilateral hippocampus equally well. We assessed the progression of hippocampal damage after TBI by using our automatic segmentation tool. Our data show that the presence of TBI, time after TBI, and whether the hippocampus was ipsilateral or contralateral to the injury were the parameters that explained hippocampal volume.

Metabolism of hyperpolarised [1-13 C]pyruvate in awake and anaesthetised rat brains

The use of hyperpolarised ¹³ C pyruvate for nononcological neurological applications has not been widespread so far, possibly due to delivery issues limiting the visibility of metabolites. First proof-of-concept results have indicated that metabolism can be detected in human brain, and this may supersede the results obtained in preclinical settings. One major difference between the experimental setups is that preclinical MRI/MRS routinely uses anaesthesia, which alters both haemodynamics and metabolism. Here, we used hyperpolarised [1-¹³ C]pyruvate to compare brain metabolism in awake rats and under isoflurane, urethane or medetomidine anaesthesia. Spectroscopic [1-¹³ C]pyruvate time courses measured sequentially showed that pyruvate-to-bicarbonate and pyruvate-to-lactate labelling rates were lower in isoflurane animals than awake animals. An increased bicarbonate-to-lactate ratio was observed in the medetomidine group compared with other groups. The study shows that hyperpolarised [1-¹³ C]pyruvate experiments can be performed in awake rats, thus avoiding anaesthesia-related issues. The results suggest that haemodynamics probably dominate the observed pyruvate-to-metabolite labelling rates and area-under-time course ratios of referenced to pyruvate. On the other hand, the results obtained with medetomidine suggest that the ratios are also modulated by the underlying cerebral metabolism. However, the ratios between intracellular metabolites were unchanged in awake compared with isoflurane-anaesthetised rats.

Whole-brain studies of spontaneous behavior in head-fixed rats enabled by zero echo time MB-SWIFT fMRI

Understanding the link between the brain activity and behavior is a key challenge in modern neuroscience. Behavioral neuroscience, however, lacks tools to record whole-brain activity in complex behavioral settings. Here we demonstrate that a novel Multi-Band SWeep Imaging with Fourier Transformation (MB-SWIFT) functional magnetic resonance imaging (fMRI) approach enables whole-brain studies in spontaneously behaving head-fixed rats. First, we show anatomically relevant functional parcellation. Second, we show sensory, motor, exploration, and stress-related brain activity in relevant networks during corresponding spontaneous behavior. Third, we show odor-induced activation of olfactory system with high correlation between the fMRI and behavioral responses. We conclude that the applied methodology enables novel behavioral study designs in rodents focusing on tasks, cognition, emotions, physical exercise, and social interaction. Importantly, novel zero echo time and large bandwidth approaches, such as MB-SWIFT, can be applied for human behavioral studies, allowing more freedom as body movement is dramatically less restricting factor.

Infantile status epilepticus disrupts myelin development

Temporal lobe epilepsy (TLE) is the most prevalent type of epilepsy in adults; it often starts in infancy or early childhood. Although TLE is primarily considered to be a grey matter pathology, a growing body of evidence links this disease with white matter abnormalities. In this study, we explore the impact of TLE onset and progression in the immature brain on white matter integrity and development utilising the rat model of Li-pilocarpine-induced TLE at the 12th postnatal day (P). Diffusion tensor imaging (DTI) and Black-Gold II histology uncovered disruptions in major white matter tracks (corpus callosum, internal and external capsules, and deep cerebral white matter) spreading through the whole brain at P28. These abnormalities were mostly not present any longer at three months after TLE induction, with only limited abnormalities detectable in the external capsule and deep cerebral white matter. Relaxation Along a Fictitious Field in the rotating frame of rank 4 indicated that white matter changes observed at both timepoints, P28 and P72, are consistent with decreased myelin content. The animals affected by TLE-induced white matter abnormalities exhibited increased functional connectivity between the thalamus and medial prefrontal and somatosensory cortex in adulthood. Furthermore, histological analyses of additional animal groups at P15 and P18 showed only mild changes in white matter integrity, suggesting a gradual age-dependent impact of TLE progression. Taken together, TLE progression in the immature brain distorts white matter development with a peak around postnatal day 28, followed by substantial recovery in adulthood. This developmental delay might give rise to cognitive and behavioural comorbidities typical for early-onset TLE.

Microstructural Tissue Changes in a Rat Model of Mild Traumatic Brain Injury

Our study investigates the potential of diffusion MRI (dMRI), including diffusion tensor imaging (DTI), fixel-based analysis (FBA) and neurite orientation dispersion and density imaging (NODDI), to detect microstructural tissue abnormalities in rats after mild traumatic brain injury (mTBI). The brains of sham-operated and mTBI rats 35 days after lateral fluid percussion injury were imaged ex vivo in a 11.7-T scanner. Voxel-based analyses of DTI-, fixel- and NODDI-based metrics detected extensive tissue changes in directly affected brain areas close to the primary injury, and more importantly, also in distal areas connected to primary injury and indirectly affected by the secondary injury mechanisms. Histology revealed ongoing axonal abnormalities and inflammation, 35 days after the injury, in the brain areas highlighted in the group analyses. Fractional anisotropy (FA), fiber density (FD) and fiber density and fiber bundle cross-section (FDC) showed similar pattern of significant areas throughout the brain; however, FA showed more significant voxels in gray matter areas, while FD and FDC in white matter areas, and orientation dispersion index (ODI) in areas most damage based on histology. Region-of-interest (ROI)-based analyses on dMRI maps and histology in selected brain regions revealed that the changes in MRI parameters could be attributed to both alterations in myelinated fiber bundles and increased cellularity. This study demonstrates that the combination of dMRI methods can provide a more complete insight into the microstructural alterations in white and gray matter after mTBI, which may aid diagnosis and prognosis following a mild brain injury.

Acute thalamic damage as a prognostic biomarker for post-traumatic epileptogenesis

OBJECTIVE

To identify magnetic resonance imaging (MRI) biomarkers for post-traumatic epilepsy.

METHODS

The EPITARGET (Targets and biomarkers for antiepileptogenesis, epitarget.eu) animal cohort completing T₂ relaxation and diffusion tensor MRI follow-up and 1-month-long video-electroencephalography monitoring included 98 male Sprague-Dawley rats with traumatic brain injury and 18 controls. T₂ imaging was performed on day (D) 2, D7, and D21 and diffusion tensor imaging (DTI) on D7 and D21 using a 7-Tesla Bruker PharmaScan MRI scanner. The mean and standard deviation (SD) of the T₂ relaxation rate, multiple diffusivity measures, and diffusion anisotropy at each time-point within the ventroposterolateral and ventroposteromedial thalamus were used as predictor variables in multi-variable logistic regression models to distinguish rats with and without epilepsy.

RESULTS

Twenty-nine percent (28/98) of the rats with traumatic brain injury (TBI) developed epilepsy. The best-performing logistic regression model utilized the D2 and D7 T₂ relaxation time as well as the D7 diffusion tensor data. The model distinguished rats with and without epilepsy (Bonferroni-corrected p-value < .001) with a cross-validated concordance statistic of 0.74 (95% confidence interval [CI] 0.60-0.84). In a cross-validated classification test, the model exhibited 54% sensitivity and 91% specificity, enriching the epilepsy rate within the study population from the expected 29% to 71%. A model using the D2 T₂ data only resulted in a 73% enriched epilepsy rate (regression p-value .007, cross-validated concordance 0.70, 95% CI 0.56-0.80, sensitivity 29%, specificity 96%).

SIGNIFICANCE

An MRI parameter set reporting on acute and subacute neuropathologic changes common to experimental and human TBI presents a diagnostic biomarker for post-traumatic epileptogenesis. Significant enrichment of the study population was achieved even when using a single time-point measurement, producing an expected epilepsy rate of 73%.

Spinal cord fMRI with MB-SWIFT for assessing epidural spinal cord stimulation in rats

Purpose

Electrical epidural spinal cord stimulation (SCS) is used as a treatment for chronic pain as well as to partially restore motor function after a spinal cord injury. Monitoring the spinal cord activity during SCS with fMRI could provide important and objective measures of integrative responses to treatment. Unfortunately, spinal cord fMRI is severely challenged by motion and susceptibility artifacts induced by the implanted electrode and bones. This pilot study introduces multi‐band sweep imaging with Fourier transformation (MB‐SWIFT) technique for spinal cord fMRI during SCS in rats. Given the close to zero acquisition delay and high bandwidth in 3 dimensions, MB‐SWIFT is demonstrated to be highly tolerant to motion and susceptibility‐induced artifacts and thus holds promise for fMRI during SCS.

Methods

MB‐SWIFT with 0.78 × 0.78 × 1.50 mm³ spatial resolution and 3‐s temporal resolution was used at 9.4 Tesla in rats undergoing epidural SCS at different frequencies. Its performance was compared with spin echo EPI. The origin of the functional contrast was also explored using suppression bands.

Results

MB‐SWIFT was tolerant to electrode‐induced artifacts and respiratory motion, leading to substantially higher fMRI sensitivity than spin echo fMRI. Clear stimulation frequency‐dependent responses to SCS were detected in the rat spinal cord close to the stimulation site. The origin of MB‐SWIFT fMRI signals was consistent with dominant inflow effects.

Conclusion

fMRI of the rat spinal cord during SCS can be consistently achieved with MB‐SWIFT, thus providing a valuable experimental framework for assessing the effects of SCS on the central nervous system.

Isoflurane affects brain functional connectivity in rats 1 month after exposure

Isoflurane, the most commonly used preclinical anesthetic, induces brain plasticity and long-term cellular and molecular changes leading to behavioral and/or cognitive consequences. These changes are most likely associated with network-level changes in brain function. To elucidate the mechanisms underlying long-term effects of isoflurane, we investigated the influence of a single isoflurane exposure on functional connectivity, brain electrical activity, and gene expression. Male Wistar rats (n = 22) were exposed to 1.8% isoflurane for 3 h. Control rats (n = 22) spent 3 h in the same room without exposure to anesthesia. After 1 month, functional connectivity was evaluated with resting-state functional magnetic resonance imaging (fMRI; n = 6 + 6) and local field potential measurements (n = 6 + 6) in anesthetized animals. A whole genome expression analysis (n = 10+10) was also conducted with mRNA-sequencing from cortical and hippocampal tissue samples. Isoflurane treatment strengthened thalamo-cortical and hippocampal-cortical functional connectivity. Cortical low-frequency fMRI power was also significantly increased in response to the isoflurane treatment. The local field potential results indicating strengthened hippocampal-cortical alpha and beta coherence were in good agreement with the fMRI findings. Furthermore, altered expression was found in 20 cortical genes, several of which are involved in neuronal signal transmission, but no gene expression changes were noted in the hippocampus. Isoflurane induced prolonged changes in thalamo-cortical and hippocampal-cortical function and expression of genes contributing to signal transmission in the cortex. Further studies are required to investigate whether these changes are associated with the postoperative behavioral and cognitive symptoms commonly observed in patients and animals.

Brain fMRI during orientation selective epidural spinal cord stimulation

Epidural spinal cord stimulation (ESCS) is widely used for chronic pain treatment, and is also a promising tool for restoring motor function after spinal cord injury. Despite significant positive impact of ESCS, currently available protocols provide limited specificity and efficiency partially due to the limited number of contacts of the leads and to the limited flexibility to vary the spatial distribution of the stimulation field in respect to the spinal cord. Recently, we introduced Orientation Selective (OS) stimulation strategies for deep brain stimulation, and demonstrated their selectivity in rats using functional MRI (fMRI). The method achieves orientation selectivity by controlling the main direction of the electric field gradients using individually driven channels. Here, we introduced a similar OS approach for ESCS, and demonstrated orientation dependent brain activations as detected by brain fMRI. The fMRI activation patterns during spinal cord stimulation demonstrated the complexity of brain networks stimulated by OS-ESCS paradigms, involving brain areas responsible for the transmission of the motor and sensory information. The OS approach may allow targeting ESCS to spinal fibers of different orientations, ultimately making stimulation less dependent on the precision of the electrode implantation.

RAFF-4, Magnetization Transfer and Diffusion Tensor MRI of Lysophosphatidylcholine Induced Demyelination and Remyelination in Rats

Remyelination is a naturally occurring response to demyelination and has a central role in the pathophysiology of multiple sclerosis and traumatic brain injury. Recently we demonstrated that a novel MRI technique entitled Relaxation Along a Fictitious Field (RAFF) in the rotating frame of rank n (RAFFn) achieved exceptional sensitivity in detecting the demyelination processes induced by lysophosphatidylcholine (LPC) in rat brain. In the present work, our aim was to test whether RAFF4, along with magnetization transfer (MT) and diffusion tensor imaging (DTI), would be capable of detecting the changes in the myelin content and microstructure caused by modifications of myelin sheets around axons or by gliosis during the remyelination phase after LPC-induced demyelination in the corpus callosum of rats. We collected MRI data with RAFF4, MT and DTI at 3 days after injection (demyelination stage) and at 38 days after injection (remyelination stage) of LPC (n = 12) or vehicle (n = 9). Cell density and myelin content were assessed by histology. All MRI metrics detected differences between LPC-injected and control groups of animals in the demyelination stage, on day 3. In the remyelination phase (day 38), RAFF4, MT parameters, fractional anisotropy, and axial diffusivity detected signs of a partial recovery consistent with the remyelination evident in histology. Radial diffusivity had undergone a further increase from day 3 to 38 and mean diffusivity revealed a complete recovery correlating with the histological assessment of cell density attributed to gliosis. The combination of RAFF4, MT and DTI has the potential to differentiate between normal, demyelinated and remyelinated axons and gliosis and thus it may be able to provide a more detailed assessment of white matter pathologies in several neurological diseases.

Data-Driven Regularization Parameter Selection in Dynamic MRI

In dynamic MRI, sufficient temporal resolution can often only be obtained using imaging protocols which produce undersampled data for each image in the time series. This has led to the popularity of compressed sensing (CS) based reconstructions. One problem in CS approaches is determining the regularization parameters, which control the balance between data fidelity and regularization. We propose a data-driven approach for the total variation regularization parameter selection, where reconstructions yield expected sparsity levels in the regularization domains. The expected sparsity levels are obtained from the measurement data for temporal regularization and from a reference image for spatial regularization. Two formulations are proposed. Simultaneous search for a parameter pair yielding expected sparsity in both domains (S-surface), and a sequential parameter selection using the S-curve method (Sequential S-curve). The approaches are evaluated using simulated and experimental DCE-MRI. In the simulated test case, both methods produce a parameter pair and reconstruction that is close to the root mean square error (RMSE) optimal pair and reconstruction. In the experimental test case, the methods produce almost equal parameter selection, and the reconstructions are of high perceived quality. Both methods lead to a highly feasible selection of the regularization parameters in both test cases while the sequential method is computationally more efficient.

Quantitative susceptibility mapping of the rat brain after traumatic brain injury

The primary lesion arising from the initial insult after traumatic brain injury (TBI) triggers a cascade of secondary tissue damage, which may also progress to connected brain areas in the chronic phase. The aim of this study was, therefore, to investigate variations in the susceptibility distribution related to these secondary tissue changes in a rat model after severe lateral fluid percussion injury. We compared quantitative susceptibility mapping (QSM) and R₂ * measurements with histological analyses in white and grey matter areas outside the primary lesion but connected to the lesion site. We demonstrate that susceptibility variations in white and grey matter areas could be attributed to reduction in myelin, accumulation of iron and calcium, and gliosis. QSM showed quantitative changes attributed to secondary damage in areas located rostral to the lesion site that appeared normal in R₂ * maps. However, combination of QSM and R₂ * was informative in disentangling the underlying tissue changes such as iron accumulation, demyelination, or calcifications. Therefore, combining QSM with R₂ * measurement can provide a more detailed assessment of tissue changes and may pave the way for improved diagnosis of TBI, and several other complex neurodegenerative diseases.

Automated joint skull-stripping and segmentation with Multi-Task U-Net in large mouse brain MRI databases

Skull-stripping and region segmentation are fundamental steps in preclinical magnetic resonance imaging (MRI) studies, and these common procedures are usually performed manually. We present Multi-task U-Net (MU-Net), a convolutional neural network designed to accomplish both tasks simultaneously. MU-Net achieved higher segmentation accuracy than state-of-the-art multi-atlas segmentation methods with an inference time of 0.35 s and no pre-processing requirements. We trained and validated MU-Net on 128 T2-weighted mouse MRI volumes as well as on the publicly available MRM NeAT dataset of 10 MRI volumes. We tested MU-Net with an unusually large dataset combining several independent studies consisting of 1782 mouse brain MRI volumes of both healthy and Huntington animals, and measured average Dice scores of 0.906 (striati), 0.937 (cortex), and 0.978 (brain mask). Further, we explored the effectiveness of our network in the presence of different architectural features, including skip connections and recently proposed framing connections, and the effects of the age range of the training set animals. These high evaluation scores demonstrate that MU-Net is a powerful tool for segmentation and skull-stripping, decreasing inter and intra-rater variability of manual segmentation. The MU-Net code and the trained model are publicly available at https://github.com/Hierakonpolis/MU-Net.

RatLesNetv2: A Fully Convolutional Network for Rodent Brain Lesion Segmentation

We present a fully convolutional neural network (ConvNet), named RatLesNetv2, for segmenting lesions in rodent magnetic resonance (MR) brain images. RatLesNetv2 architecture resembles an autoencoder and it incorporates residual blocks that facilitate its optimization. RatLesNetv2 is trained end to end on three-dimensional images and it requires no preprocessing. We evaluated RatLesNetv2 on an exceptionally large dataset composed of 916 T2-weighted rat brain MRI scans of 671 rats at nine different lesion stages that were used to study focal cerebral ischemia for drug development. In addition, we compared its performance with three other ConvNets specifically designed for medical image segmentation. RatLesNetv2 obtained similar to higher Dice coefficient values than the other ConvNets and it produced much more realistic and compact segmentations with notably fewer holes and lower Hausdorff distance. The Dice scores of RatLesNetv2 segmentations also exceeded inter-rater agreement of manual segmentations. In conclusion, RatLesNetv2 could be used for automated lesion segmentation, reducing human workload and improving reproducibility. RatLesNetv2 is publicly available at https://github.com/jmlipman/RatLesNetv2.

Assessment of the structural complexity of diffusion MRI voxels using 3D electron microscopy in the rat brain

Validation and interpretation of diffusion magnetic resonance imaging (dMRI) requires detailed understanding of the actual microstructure restricting the diffusion of water molecules. In this study, we used serial block-face scanning electron microscopy (SBEM), a three-dimensional electron microscopy (3D-EM) technique, to image seven white and grey matter volumes in the rat brain. SBEM shows excellent contrast of cellular membranes, which are the major components restricting the diffusion of water in tissue. Additionally, we performed 3D structure tensor (3D-ST) analysis on the SBEM volumes and parameterised the resulting orientation distributions using Watson and angular central Gaussian (ACG) probability distributions as well as spherical harmonic (SH) decomposition. We analysed how these parameterisations described the underlying orientation distributions and compared their orientation and dispersion with corresponding parameters from two dMRI methods, neurite orientation dispersion and density imaging (NODDI) and constrained spherical deconvolution (CSD). Watson and ACG parameterisations and SH decomposition captured well the 3D-ST orientation distributions, but ACG and SH better represented the distributions due to its ability to model asymmetric dispersion. The dMRI parameters corresponded well with the 3D-ST parameters in the white matter volumes, but the correspondence was less evident in the more complex grey matter. SBEM imaging and 3D-ST analysis also revealed that the orientation distributions were often not axially symmetric, a property neatly captured by the ACG distribution. Overall, the ability of SBEM to image diffusion barriers in intricate detail, combined with 3D-ST analysis and parameterisation, provides a step forward toward interpreting and validating the dMRI signals in complex brain tissue microstructure.

Long-lasting blood-brain barrier dysfunction and neuroinflammation after traumatic brain injury

BACKGROUND

Traumatic brain injury (TBI) causes 10-20% of acquired epilepsy, which typically develops within 2 years post-injury with poorly understood mechanisms. We investigated the location, severity, evolution and persistence of blood-brain barrier (BBB) dysfunction and associated neuroinflammation after TBI, and their contribution to post-traumatic seizure susceptibility.

METHODS

TBI was induced with lateral fluid-percussion in adult male Sprague-Dawley rats (6 sham, 12 TBI). Permeability of the BBB was assessed using T1-weighted magnetic resonance imaging (MRI) with gadobutrol (Gd) contrast enhancement at 4 days, 2 weeks, 2 months, and 10 months post-injury and with intravenously administered fluorescein at 11 months post-TBI. Continuous (24/7) video-EEG monitoring was performed for 3 weeks at 11 months post-injury followed by the pentylenetetrazol (PTZ) seizure-susceptibility test. In the end, rats were perfused for histology to assess albumin extravasation, iron deposits, calcifications, reactive astrocytes, microglia and monocytes. To investigate the translational value of the data obtained, BBB dysfunction and neuroinflammation were investigated immunohistochemically in autopsy brain tissue from patients with TBI and PTE.

RESULTS

MRI indicated persistent Gd leakage in the impacted cortex and thalamus of variable severity in all rats with TBI which correlated with fluorescein extravasation. In the impacted cortex BBB dysfunction was evident from 4 days post-injury onwards to the end of the 10-months follow-up. In the ipsilateral thalamus, leakage was evident at 2 and 10 months post-injury. The greater the BBB leakage in the perilesional cortex at 10 months after the injury, the greater the expression of the endothelial cell antigen RECA-1 (r = 0.734, p < 0.01) and the activated macrophages/monocytes/microglia marker CD68 (r = 0.699, p < 0.05) at 11 months post-injury. Seven of the 12 rats with TBI showed increased seizure susceptibility in the PTZ-test. Unlike expected, we did not find any association between increased Gd-leakage or neuroinflammation with seizure susceptibility at 11 months post-TBI. Analysis of human autopsy tissue indicated that similar to the animal model, chronic BBB dysfunction was also evident in the perilesional cortex and thalamus of patients with PTE, characterized by presence of albumin, iron deposits and calcifications as well as markers of neuroinflammation, including reactive astrocytes, microglia and monocytes.

CONCLUSIONS

Rats and humans with TBI have long-lasting cortical BBB dysfunction and neuroinflammation. Focal Gd-enhancement matched with loci of neuroinflammation, particularly in the thalamus. Although BBB leakage did not associate with increased seizure susceptibility after TBI, our data suggest that for treatments aimed to mitigate BBB damage and its secondary pathologies like chronic neuroinflammation, there is a region-specific, long-lasting therapeutic time window.

Early Increase in Cortical T2 Relaxation Is a Prognostic Biomarker for the Evolution of Severe Cortical Damage, but Not for Epileptogenesis, after Experimental Traumatic Brain Injury

Prognostic biomarkers for post-injury outcome are necessary for the development of neuroprotective and antiepileptogenic treatments for traumatic brain injury (TBI). We hypothesized that T₂ relaxation magnetic resonance imaging (MRI) predicts the progression of perilesional cortical pathology and epileptogenesis. The EPITARGET animal cohort used for MRI analysis included 120 adult male Sprague-Dawley rats with TBI induced by lateral fluid-percussion injury and 24 sham-operated controls. T₂ MRI was performed at days 2, 7, and 21 post-TBI. The lesioned cortex was outlined, and the T₂ value of each imaging voxel within the lesion area was scored using a five-grade pathology classification. Analysis of 1-month video-electroencephalography recordings initiated 5 months post-TBI indicated that 27% (31 of 114) of the animals with TBI developed epilepsy. Multiple linear regression analysis indicated that T₂-based classification of lesion volume at day 2 and day 7 post-TBI explained the necrotic lesion volume with greatly increased T₂ (>102 ms) at day 21 post-TBI (F_(13,103) = 52.5; p < 0.001; R² =  0.87; adjusted R² = 0.85). The volume of moderately increased (78-102 ms) T₂ at day 7 post-TBI predicted the evolution of large (>12 mm³) cortical lesions (area under the curve, 0.92; p < 0.001; cutoff, 1.9 mm³; false positive rate, 0.10; true positive rate, 0.62). Logistic regression analysis, however, showed that the different severities of T₂ lesion volumes at days 2, 7, and 21 post-TBI did not explain the development of epilepsy (χ²_(18,95) = 18.4; p = 0.427). In addition, the location of the T₂ abnormality within the cortex did not correlate with epileptogenesis. A single measurement of T₂ relaxation MRI in the acute post-TBI phase is useful for identifying post-TBI subjects at highest risk of developing large cortical lesions, and thus, in the greatest need of neuroprotective therapies after TBI, but not the development of post-traumatic epilepsy.

Orientation selective deep brain stimulation of the subthalamic nucleus in rats

Deep brain stimulation (DBS) has become an important tool in the management of a wide spectrum of diseases in neurology and psychiatry. Target selection is a vital aspect of DBS so that only the desired areas are stimulated. Segmented leads and current steering have been shown to be promising additions to DBS technology enabling better control of the stimulating electric field. Recently introduced orientation selective DBS (OS-DBS) is a related development permitting sensitization of the stimulus to axonal pathways with different orientations by freely controlling the primary direction of the electric field using multiple contacts. Here, we used OS-DBS to stimulate the subthalamic nucleus (STN) in healthy rats while simultaneously monitoring the induced brain activity with fMRI. Maximal activation of the sensorimotor and basal ganglia-thalamocortical networks was observed when the electric field was aligned mediolaterally in the STN pointing in the lateral direction, while no cortical activation was observed with the electric field pointing medially to the opposite direction. Such findings are consistent with mediolateral main direction of the STN fibers, as seen with high resolution diffusion imaging and histology. The asymmetry of the OS-DBS dipolar field distribution using three contacts along with the potential stimulation of the internal capsule, are also discussed. We conclude that OS-DBS offers an additional degree of flexibility for optimization of DBS of the STN which may enable a better treatment response.

Multiscale Imaging Approach for Studying the Central Nervous System: Methodology and Perspective

Non-invasive imaging methods have become essential tools for understanding the central nervous system (CNS) in health and disease. In particular, magnetic resonance imaging (MRI) techniques provide information about the anatomy, microstructure, and function of the brain and spinal cord in vivo non-invasively. However, MRI is limited by its spatial resolution and signal specificity. In order to mitigate these shortcomings, it is crucial to validate MRI with an array of ancillary ex vivo imaging techniques. These techniques include histological methods, such as light and electron microscopy (EM), which can provide specific information on the tissue structure in healthy and diseased brain and spinal cord, at cellular and subcellular level. However, these conventional histological techniques are intrinsically two-dimensional (2D) and, as a result of sectioning, lack volumetric information of the tissue. This limitation can be overcome with genuine three-dimensional (3D) imaging approaches of the tissue. 3D highly resolved information of the CNS achievable by means of other imaging techniques can complement and improve the interpretation of MRI measurements. In this article, we provide an overview of different 3D imaging techniques that can be used to validate MRI. As an example, we introduce an approach of how to combine diffusion MRI and synchrotron X-ray phase contrast tomography (SXRPCT) data. Our approach paves the way for a new multiscale assessment of the CNS allowing to validate and to improve our understanding of in vivo imaging (such as MRI).

Multi-band SWIFT enables quiet and artefact-free EEG-fMRI and awake fMRI studies in rat

Functional magnetic resonance imaging (fMRI) studies in animal models provide invaluable information regarding normal and abnormal brain function, especially when combined with complementary stimulation and recording techniques. The echo planar imaging (EPI) pulse sequence is the most common choice for fMRI investigations, but it has several shortcomings. EPI is one of the loudest sequences and very prone to movement and susceptibility-induced artefacts, making it suboptimal for awake imaging. Additionally, the fast gradient-switching of EPI induces disrupting currents in simultaneous electrophysiological recordings. Therefore, we investigated whether the unique features of Multi-Band SWeep Imaging with Fourier Transformation (MB-SWIFT) overcome these issues at a high 9.4 T magnetic field, making it a potential alternative to EPI. MB-SWIFT had 32-dB and 20-dB lower peak and average sound pressure levels, respectively, than EPI with typical fMRI parameters. Body movements had little to no effect on MB-SWIFT images or functional connectivity analyses, whereas they severely affected EPI data. The minimal gradient steps of MB-SWIFT induced significantly lower currents in simultaneous electrophysiological recordings than EPI, and there were no electrode-induced distortions in MB-SWIFT images. An independent component analysis of the awake rat functional connectivity data obtained with MB-SWIFT resulted in near whole-brain level functional parcellation, and simultaneous electrophysiological and fMRI measurements in isoflurane-anesthetized rats indicated that MB-SWIFT signal is tightly linked to neuronal resting-state activity. Therefore, we conclude that the MB-SWIFT sequence is a robust preclinical brain mapping tool that can overcome many of the drawbacks of conventional EPI fMRI at high magnetic fields.

MRS Reveals Chronic Inflammation in T2w MRI-Negative Perilesional Cortex - A 6-Months Multimodal Imaging Follow-Up Study

Sustained inflammation in the injured cortex is a promising therapeutic target for disease-modification after traumatic brain injury (TBI). However, its extent and dynamics of expansion are incompletely understood which challenges the timing and placement of therapeutics to lesioned area. Our aim was to characterize the evolution of chronic inflammation during lesion expansion in lateral fluid-percussion injury (FPI) rat model with focus on the MRI-negative perilesional cortex. T2-weighted MR imaging (T2w MRI) and localized magnetic resonance spectroscopy (MRS) were performed at 1, 3, and 6 months post-injury. End-point histology, including Nissl for neuronal death, GFAP for astrogliosis, and Prussian Blue for iron were used to assess perilesional histopathology. An additional animal cohort was imaged with a positron emission tomography (PET) using translocator protein 18 kDa (TSPO) radiotracer [¹⁸F]-FEPPA. T2w MRI assessed lesion growth and detected chronic inflammation along the lesion border while rest of the ipsilateral cortex was MRI-negative (MRI-). Instead, myo-inositol that is an inflammatory MRS marker for gliosis, glutathione for oxidative stress, and choline for membrane turnover were elevated throughout the 6-months follow-up in the MRI- perilesional cortex (all p < 0.05). MRS markers revealed chronically sustained inflammation across the ipsilateral cortex but did not indicate the upcoming lesion expansion. Instead, the rostral expansion of the cortical lesion was systematically preceded by a hyperintense band in T2w images months earlier. Histologic analysis of the hyperintensity indicated scattered astrocytes, incomplete glial scar, and intracellularly packed and free iron. Yet, the band was negative in [¹⁸F]-FEPPA-PET. [¹⁸F]-FEPPA also showed no cortical TSPO expression within the MRS voxel in MRI- perilesional cortex or anywhere along glial scar when assessed at 2 months post-injury. However, [¹⁸F]-FEPPA showed a robust signal increase, indicating reactive microgliosis in the ipsilateral thalamus at 2 months post-TBI. We present evidence that MRS reveals chronic posttraumatic inflammation in MRI-negative perilesional cortex. The mismatch in MRS, MRI, and PET measures may allow non-invasive endophenotyping of beneficial and detrimental inflammatory processes to aid targeting and timing of anti-inflammatory therapeutics.

Impaired hippocampal-cortical coupling but preserved local synchrony during sleep in APP/PS1 mice modeling Alzheimer's disease

Sleep, in addition to its brain restorative processes, plays an important role in memory transfer from its temporary store in the hippocampus to the more permanent storage in the neocortex. Alzheimer’s disease (AD) affects memory and sleep. The aim of this study was to explore disturbances in global and local synchrony patterns between brain regions in the APP/PS1 mouse model of the AD during natural sleep. We used 8 male APPswe/PS1dE9 mice and 6 wild-type littermates, aged 5–6 months, with multiple electrode bundles implanted into cortical regions, thalamus and hippocampus. We measured video-EEG in freely moving animals and analyzed synchrony during NREM vs REM sleep. Global synchrony between medial frontal cortex and hippocampus measured with magnitude-squared coherence was slightly decreased in delta range during NREM stage of sleep in APP/PS1 mice. In contrast, local hippocampal synchrony measured with cross-frequency coupling remained intact. Ripple structure or frequency did not differ between the genotypes. However, the coupling of the spindle-band power peak in the medial prefrontal cortex to hippocampal ripples was significantly decreased compared to wild-type animals. The delicate timing of hippocampal ripples, frontal delta, and corticothalamic spindle oscillations may be the first sign of impaired memory in amyloid plaque-forming transgenic mice.

Mechanosensitive meningeal nociception via Piezo channels: Implications for pulsatile pain in migraine?

BACKGROUND

Recent discovery of mechanosensitive Piezo receptors in trigeminal ganglia suggested the novel molecular candidate for generation of migraine pain. However, the contribution of Piezo channels in migraine pathology was not tested yet. Therefore, in this study, we explored a potential involvement of Piezo channels in peripheral trigeminal nociception implicated in generation of migraine pain.

METHODS

We used immunohistochemistry, calcium imaging, calcitonin gene related peptide (CGRP) release assay and electrophysiology in mouse and rat isolated trigeminal neurons and rat hemiskulls to study action of various stimulants of Piezo receptors on migraine-related peripheral nociception.

RESULTS

We found that essential (35%) fraction of isolated rat trigeminal neurons responded to chemical Piezo1 agonist Yoda1 and about a half of Yoda1 positive neurons responded to hypo-osmotic solution (HOS) and a quarter to mechanical stimulation by focused ultrasound (US). In ex vivo hemiskull preparation, Yoda1 and HOS largely activated persistent nociceptive firing in meningeal branches of trigeminal nerve. By using our novel cluster analysis of pain spikes, we demonstrated that 42% of fibers responded to Piezo1 agonist and 20% of trigeminal fibers were activated by Yoda1 and by capsaicin, suggesting expression of Piezo receptors in TRPV1 positive peptidergic nociceptive nerve fibers. Consistent with this, Yoda1 promoted the release of the key migraine mediator CGRP from hemiskull preparation.

CONCLUSION

Taken together, our data suggest the involvement of mechanosensitive Piezo receptors, in particular, Piezo1 subtype in peripheral trigeminal nociception, which provides a new view on mechanotransduction in migraine pathology and suggests novel molecular targets for anti-migraine medicine.

Harmonization of pipeline for preclinical multicenter MRI biomarker discovery in a rat model of post-traumatic epileptogenesis

Preclinical imaging studies of posttraumatic epileptogenesis (PTE) have largely been proof-of-concept studies with limited animal numbers, and thus lack the statistical power for biomarker discovery. Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx) is a pioneering multicenter trial investigating preclinical imaging biomarkers of PTE. EpiBios4Rx faced the issue of harmonizing the magnetic resonance imaging (MRI) procedures and imaging data metrics prior to its execution. We present here the harmonization process between three preclinical MRI facilities at the University of Eastern Finland (UEF), the University of Melbourne (Melbourne), and the University of California, Los Angeles (UCLA), and evaluate the uniformity of the obtained MRI data. Adult, male rats underwent a lateral fluid percussion injury (FPI) and were followed by MRI 2 days, 9 days, 1 month, and 5 months post-injury. Ex vivo scans of fixed brains were conducted 7 months post-injury as an end point follow-up. Four MRI modalities were used: T2-weighted imaging, multi-gradient-echo imaging, diffusion-weighted imaging, and magnetization transfer imaging, and acquisition parameters for each modality were tailored to account for the different field strengths (4.7 T and 7 T) and different MR hardwares used at the three participating centers. Pilot data collection resulted in comparable image quality across sites. In interim analysis (of data obtained by April 30, 2018), the within-site variation of the quantified signal properties was low, while some differences between sites remained. In T2-weighted images the signal-to-noise ratios were high at each site, being 35 at UEF, 48 at Melbourne, and 32 at UCLA (p < 0.05). The contrast-to-noise ratios were similar between the sites (9, 10, and 8, respectively). Magnetization transfer ratio maps had identical white matter/ gray matter contrast between the sites, with white matter showing 15% higher MTR than gray matter despite different absolute MTR values (MTR both in white and gray matter was 3% lower in Melbourne than at UEF, p < 0.05). Diffusion-weighting yielded different degrees of signal attenuation across sites, being 83% at UEF, 76% in Melbourne, and 80% at UCLA (p < 0.05). Fractional anisotropy values differed as well, being 0.81 at UEF, 0.73 in Melbourne, and 0.84 at UCLA (p < 0.05). The obtained values in sham animals showed low variation within each site and no change over time, suggesting high repeatability of the measurements. Quality control scans with phantoms demonstrated stable hardware performance over time. Timing of post-TBI scans was designed to target specific phases of the dynamic pathology, and the execution at different centers was highly accurate. Besides a few outliers, the 2-day scans were done within an hour from the target time point. At day 9, most animals were scanned within an hour from the target time point, and all but 2 outliers within 24 h from the target. The 1-month post-TBI scans were done within 31 ± 3 days. MRI procedures and animal physiology during scans were similar between the sites. Taken together, the 10% inter-site difference in FA and 3% difference in MTR values should be included into analysis as a covariate or balanced out in post-processing in order to detect disease-related effects on brain structure at the same scale. However, for a MRI biomarker for post-traumatic epileptogenesis to have realistic chance of being successfully translated to validation in clinical trials, it would need to be a robust TBI-induced structural change which tolerates the inter-site methodological variability described here.

Sleep-State Dependent Alterations in Brain Functional Connectivity under Urethane Anesthesia in a Rat Model of Early-Stage Parkinson's Disease

Parkinson's disease (PD) is characterized by the gradual degeneration of dopaminergic neurons in the substantia nigra, leading to striatal dopamine depletion. A partial unilateral striatal 6-hydroxydopamine (6-OHDA) lesion causes 40-60% dopamine depletion in the lesioned rat striatum, modeling the early stage of PD. In this study, we explored the connectivity between the brain regions in partially 6-OHDA lesioned male Wistar rats under urethane anesthesia using functional magnetic resonance imaging (fMRI) at 5 weeks after the 6-OHDA infusion. Under urethane anesthesia, the brain fluctuates between the two states, resembling rapid eye movement (REM) and non-REM sleep states. We observed clear urethane-induced sleep-like states in 8/19 lesioned animals and 8/18 control animals. 6-OHDA lesioned animals exhibited significantly lower functional connectivity between the brain regions. However, we observed these differences only during the REM-like sleep state, suggesting the involvement of the nigrostriatal dopaminergic pathway in REM sleep regulation. Corticocortical and corticostriatal connections were decreased in both hemispheres, reflecting the global effect of the lesion. Overall, this study describes a promising model to study PD-related sleep disorders in rats using fMRI.

Tuning Neuromodulation Effects by Orientation Selective Deep Brain Stimulation in the Rat Medial Frontal Cortex

Previous studies that focused on treating major depressive disorder with conventional deep brain stimulation (DBS) paradigms produced inconsistent results. In this proof-of-concept preclinical study in rats (n = 8), we used novel paradigms of orientation selective DBS for stimulating the complex circuitry crossing the infralimbic cortex, an area considered analogous to human subgenual cingulate cortex. Using functional MRI at 9.4 T, we monitored whole brain responses to varying the electrical field orientation of DBS within the infralimbic cortex. Substantial alterations of functional MRI responses in the amygdala, a major node connected to the infralimbic cortex implicated in the pathophysiology of depression, were observed. As expected, the activation cluster near the electrode was insensitive to the changes of the stimulation orientation. Hence, our findings substantiate the ability of orientation selective stimulation (OSS) to recruit neuronal pathways of distinct orientations relative to the position of the electrode, even in complex circuits such as those involved in major depressive disorder. We conclude that OSS is a promising approach for stimulating brain areas that inherently require individualisation of the treatment approach.

Awake Rat Brain Functional Magnetic Resonance Imaging Using Standard Radio Frequency Coils and a 3D Printed Restraint Kit

Functional magnetic resonance imaging (fMRI) is a powerful noninvasive tool for studying spontaneous resting state functional connectivity (RSFC) in laboratory animals. Brain function can be significantly affected by generally used anesthetics, however, rendering the need for awake imaging. Only a few different awake animal habituation protocols have been presented, and there is a critical need for practical and improved low-stress techniques. Here we demonstrate a novel restraint approach for awake rat RSFC studies. Our custom-made 3D printed restraint kit is compatible with a standard Bruker Biospin MRI rat bed, rat brain receiver coil, and volume transmitter coil. We also implemented a progressive habituation protocol aiming to minimize the stress experienced by the rats, and compared RSFC between awake, lightly sedated, and isoflurane-anesthetized rats. Our results demonstrated that the 3D printed restraint kit was suitable for RSFC studies of awake rats. During the short 4-day habituation period, the plasma corticosterone concentration, movement, and heart rate, which were measured as stress indicators, decreased significantly, indicating adaptation to the restraint protocol. Additionally, 10 days after the awake MRI session, rats exhibited no signs of depression or anxiety based on open-field and sucrose preference behavioral tests. The RSFC data revealed significant changes in the thalamo-cortical and cortico-cortical networks between the awake, lightly sedated, and anesthetized groups, emphasizing the need for awake imaging. The present work demonstrates the feasibility of our custom-made 3D printed restraint kit. Using this kit, we found that isoflurane markedly affected brain connectivity compared with that in awake rats, and that the effect was less pronounced, but still significant, when light isoflurane sedation was used instead.

Clinical deep brain stimulation strategies for orientation-selective pathway activation

OBJECTIVE

This study investigated stimulation strategies to increase the selectivity of activating axonal pathways within the brain based on their orientations relative to clinical deep brain stimulation (DBS) lead implants.

APPROACH

Previous work has shown how varying electrode shape and controlling the primary electric field direction through preclinical electrode arrays can produce orientation-selective axonal stimulation. Here, we significantly extend those results using computational models to evaluate the degree to which clinical DBS leads can direct stimulus-induced electric fields and generate orientation-selective activation of fiber pathways in the brain. Orientation-selective pulse paradigms were evaluated in conceptual models and in patient-specific models of subthalamic nucleus (STN)-DBS for treating Parkinson's disease.

MAIN RESULTS

Single-contact monopolar or two-contact bipolar stimulation through clinical DBS leads with cylindrical electrodes primarily activated axons orientated parallel to the lead. Conversely, multi-contact monopolar stimulation with a cathode-leading pulse waveform selectively activated axons perpendicular to the DBS lead. Clinical DBS leads with segmented rows of electrodes and a single current source provided additional angular resolution for activating axons oriented 0°, ±22.5°, ±45°, ±67.5°, or 90° relative to the lead shaft. Employing multiple independent current sources to deliver unequal amounts of current through these leads further increased the angular resolution of activation relative to the lead shaft. The patient-specific models indicated that multi-contact cathode configurations, which are rarely used in clinical practice, could increase activation of the hyperdirect pathway collaterals projecting into STN (a putative therapeutic target), while minimizing direct activation of the corticospinal tract of internal capsule, which can elicit sensorimotor side-effects when stimulated.

SIGNIFICANCE

When combined with patient-specific tissue anisotropy and patient-specific anatomical morphologies of neural pathways responsible for therapy and side effects, orientation-selective DBS approaches show potential to significantly improve clinical outcomes of DBS therapy for a range of existing and investigational clinical indications.