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.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo.

PURPOSE: Recent studies indicate that T1 in white matter (WM) is influenced by fiber orientation in B0 . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T1 relaxation time in humans in vivo as well as in rat brain ex vivo. METHODS: Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T1 plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T1 and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T1 within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B0 and T1 , and diffusion MRI images acquired at 9.4 T. T1 angular plots were determined at several rotation angles in B0 . RESULTS: Angular T1 plots from global WM provided reference for estimated fiber orientation-linked T1 changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T1 , matching that estimated from WM T1 data. In CC, where large and giant axons are numerous, the measured T1 change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T1 plots at 9.4 T, matching those observed at 7 T in vivo. CONCLUSION: These data causally link axon fiber orientation in B0 to the T1 relaxation anisotropy in WM.

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.

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.

MRI-Guided Electrode Implantation for Chronic Intracerebral Recordings in a Rat Model of Post-Traumatic Epilepsy-Challenges and Gains.

Brain atrophy induced by traumatic brain injury (TBI) progresses in parallel with epileptogenesis over time, and thus accurate placement of intracerebral electrodes to monitor seizure initiation and spread at the chronic postinjury phase is challenging. We evaluated in adult male Sprague Dawley rats whether adjusting atlas-based electrode coordinates on the basis of magnetic resonance imaging (MRI) increases electrode placement accuracy and the effect of chronic electrode implantations on TBI-induced brain atrophy. One group of rats (EEG cohort) was implanted with two intracortical (anterior and posterior) and a hippocampal electrode right after TBI to target coordinates calculated using a rat brain atlas. Another group (MRI cohort) was implanted with the same electrodes, but using T2-weighted MRI to adjust the planned atlas-based 3D coordinates of each electrode. Histological analysis revealed that the anterior cortical electrode was in the cortex in 83% (25% in targeted layer V) of the EEG cohort and 76% (31%) of the MRI cohort. The posterior cortical electrode was in the cortex in 40% of the EEG cohort and 60% of the MRI cohort. Without MRI-guided adjustment of electrode tip coordinates, 58% of the posterior cortical electrodes in the MRI cohort will be in the lesion cavity, as revealed by simulated electrode placement on histological images. The hippocampal electrode was accurately placed in 82% of the EEG cohort and 86% of the MRI cohort. Misplacement of intracortical electrodes related to their rostral shift due to TBI-induced cortical and hippocampal atrophy and caudal retraction of the brain, and was more severe ipsilaterally than contralaterally (p < 0.001). Total lesion area in cortical subfields targeted by the electrodes (primary somatosensory cortex, visual cortex) was similar between cohorts (p > 0.05). MRI-guided adjustment of coordinates for electrodes improved the success rate of intracortical electrode tip placement nearly to that at the acute postinjury phase (68% vs. 62%), particularly in the posterior brain, which exhibited the most severe postinjury atrophy. Overall, MRI-guided electrode implantation improved the quality and interpretation of the origin of EEG-recorded signals.

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 > 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.

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.

Relaxation anisotropy of quantitative MRI parameters in biological tissues.

Quantitative MR relaxation parameters vary in the sensitivity to the orientation of the tissue in the magnetic field. In this study, the orientation dependence of multiple relaxation parameters was assessed in various tissues. Ex vivo samples of each tissue type were prepared either from bovine knee (tendon, cartilage) or mouse (brain, spinal cord, heart, kidney), and imaged at 9.4 T MRI with T1, T2, continuous wave (CW-) T1ρ, adiabatic T1ρ and T2ρ, and Relaxation along fictitious field (RAFF2-4) sequences at five different orientations with respect to the main magnetic field. Relaxation anisotropy of the measured parameters was quantified and compared. The highly ordered collagenous tissues, i.e. cartilage and tendon, presented the highest relaxation anisotropy for T2, CW-T1ρ with spin-lock power < 1 kHz, Ad-T2ρ and RAFF2-4. Maximally anisotropy was 75% in cartilage and 30% in tendon. T1 and adiabatic T1ρ did not exhibit observable anisotropy. In the other measured tissue types, anisotropy was overall less than 10% for all the parameters. The results confirm that highly ordered collagenous tissues have properties that induce very clearly observable relaxation anisotropy, whereas in other tissues the effect is not as prominent. Quantitative comparison of anisotropy of different relaxation parameters highlights the importance of sequence choice and design in MR imaging.

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.

Hippocampal position and orientation as prognostic biomarkers for posttraumatic epileptogenesis: An experimental study in a 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, T 2 ∗ -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.

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.

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.

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.

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.