Epilepsy-related functional brain network alterations are already present at an early age in the GAERS rat model of genetic absence epilepsy.

Introduction: Genetic Absence Epilepsy Rats from Strasbourg (GAERS) represent a model of genetic generalized epilepsy. The present longitudinal study in GAERS and age-matched non-epileptic controls (NEC) aimed to characterize the epileptic brain network using two functional measures, resting state-functional magnetic resonance imaging (rs-fMRI) and manganese-enhanced MRI (MEMRI) combined with morphometry, and to investigate potential brain network alterations, following long-term seizure activity. Methods: Repeated rs-fMRI measurements at 9.4 T between 3 and 8 months of age were combined with MEMRI at the final time point of the study. We used graph theory analysis to infer community structure and global and local network parameters from rs-fMRI data and compared them to brain region-wise manganese accumulation patterns and deformation-based morphometry (DBM). Results: Functional connectivity (FC) was generally higher in GAERS when compared to NEC. Global network parameters and community structure were similar in NEC and GAERS, suggesting efficiently functioning networks in both strains. No progressive FC changes were observed in epileptic animals. Network-based statistics (NBS) revealed stronger FC within the cortical community, including regions of association and sensorimotor cortex, and with basal ganglia and limbic regions in GAERS, irrespective of age. Higher manganese accumulation in GAERS than in NEC was observed at 8 months of age, consistent with higher overall rs-FC, particularly in sensorimotor cortex and association cortex regions. Functional measures showed less similarity in subcortical regions. Whole brain volumes of 8 months-old GAERS were higher when compared to age-matched NEC, and DBM revealed increased volumes of several association and sensorimotor cortex regions and of the thalamus. Discussion: rs-fMRI, MEMRI, and volumetric data collectively suggest the significance of cortical networks in GAERS, which correlates with an increased fronto-central connectivity in childhood absence epilepsy (CAE). Our findings also verify involvement of basal ganglia and limbic regions. Epilepsy-related network alterations are already present in juvenile animals. Consequently, this early condition seems to play a greater role in dynamic brain function than chronic absence seizures.

Generation of a whole-brain hemodynamic response function and sex-specific differences in cerebral processing of mechano-sensation in mice detected by BOLD fMRI.

BOLD fMRI has become a prevalent method to study cerebral sensory processing in rodent disease models, including pain and mechanical hypersensitivity. fMRI data analysis is frequently combined with a general-linear-model (GLM) -based analysis, which uses the convolution of a hemodynamic response function (HRF) with the stimulus paradigm. However, several studies indicated that the HRF differs across species, sexes, brain structures, and experimental factors, including stimulation modalities or anesthesia, and hence might strongly affect the outcome of BOLD analyzes. While considerable work has been done in humans and rats to understand the HRF, much less is known in mice. As a prerequisite to investigate mechano-sensory processing and BOLD fMRI data in male and female mice, we (1) designed a rotating stimulator that allows application of two different mechanical modalities, including innocuous von Frey and noxious pinprick stimuli and (2) determined and statistically compared HRFs across 30 brain structures and experimental conditions, including sex and, stimulus modalities. We found that mechanical stimulation lead to brain-wide BOLD signal changes thereby allowing extraction of HRFs from multiple brain structures. However, we did not find differences in HRFs across all brain structures and experimental conditions. Hence, we computed a whole-brain mouse HRF, which is based on 88 functional scans from 30 mice. A comparison of this mouse-specific HRF with our previously reported rat-derived HRF showed significantly slower kinetics in mice. Finally, we detected pronounced differences in cerebral BOLD activation between male and female mice with mechanical stimulation, thereby exposing divergent processing of noxious and innocuous stimuli in both sexes.

The impact of vasomotion on analysis of rodent fMRI data.

Introduction: Small animal fMRI is an essential part of translational research in the cognitive neurosciences. Due to small dimensions and animal physiology preclinical fMRI is prone to artifacts that may lead to misinterpretation of the data. To reach unbiased translational conclusions, it is, therefore, crucial to identify potential sources of experimental noise and to develop correction methods for contributions that cannot be avoided such as physiological noise. Aim of this study was to assess origin and prevalence of hemodynamic oscillations (HDO) in preclinical fMRI in rat, as well as their impact on data analysis. Methods: Following the development of algorithms for HDO detection and suppression, HDO prevalence in fMRI measurements was investigated for different anesthetic regimens, comprising isoflurane and medetomidine, and for both gradient echo and spin echo fMRI sequences. In addition to assessing the effect of vasodilation on HDO, it was studied if HDO have a direct neuronal correlate using local field potential (LFP) recordings. Finally, the impact of HDO on analysis of fMRI data was assessed, studying both the impact on calculation of activation maps as well as the impact on brain network analysis. Overall, 303 fMRI measurements and 32 LFP recordings were performed in 71 rats. Results: In total, 62% of the fMRI measurements showed HDO with a frequency of (0.20 ± 0.02) Hz. This frequent occurrence indicated that HDO cannot be generally neglected in fMRI experiments. Using the developed algorithms, HDO were detected with a specificity of 95%, and removed efficiently from the signal time courses. HDO occurred brain-wide under vasoconstrictive conditions in both small and large blood vessels. Vasodilation immediately interrupted HDO, which, however, returned within 1 h under vasoconstrictive conditions. No direct neuronal correlate of HDO was observed in LFP recordings. HDO significantly impacted analysis of fMRI data, leading to altered cluster sizes and F-values for activated voxels, as well as altered brain networks, when comparing data with and without HDO. Discussion: We therefore conclude that HDO are caused by vasomotion under certain anesthetic conditions and should be corrected during fMRI data analysis to avoid bias.

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.

Fiber-based lactate recordings with fluorescence resonance energy transfer sensors by applying an magnetic resonance-informed correction of hemodynamic artifacts.

Significance: Fluorescence resonance energy transfer (FRET) sensors offer enormous benefits when studying neurophysiology through confocal microscopy. Yet, their use for fiber-based in vivo recordings is hampered by massive confounding effects and has therefore been scarcely reported. Aim: We aim to investigate whether in vivo fiber-based lactate recordings in the rodent brain are feasible with FRET sensors and implement a correction algorithm for the predominant hemodynamic artifact. Approach: We performed fiber-based FRET recordings of lactate (Laconic) and calcium (Twitch-2B) simultaneously with functional MRI and pharmacological MRI. MR-derived parameters were applied to correct hemodynamic artifacts. Results of FRET measurements were validated by local field potential, magnetic resonance spectroscopy, and blood analysis. Results: Hemodynamic artifacts dominated fiber-based in vivo FRET measurements with both Laconic and Twitch-2B. Our MR-based correction algorithm enabled to remove the artifacts and detect lactate and calcium changes during sensory stimulation or intravenous lactate injections. Conclusions: In vivo fiber-based lactate recordings are feasible using FRET-based sensors. However, signal corrections are required. MR-derived hemodynamic parameters can successfully be applied for artifact correction.

Retrosplenial Cortex Contributes to Network Changes during Seizures in the GAERS Absence Epilepsy Rat Model.

Resting state-fMRI was performed to explore brain networks in Genetic Absence Epilepsy Rats from Strasbourg and in nonepileptic controls (NEC) during monitoring of the brain state by simultaneous optical Ca2+-recordings. Graph theoretical analysis allowed for the identification of acute and chronic network changes and revealed preserved small world topology before and after seizure onset. The most prominent acute change in network organization during seizures was the segregation of cortical regions from the remaining brain. Stronger connections between thalamic with limbic regions compared with preseizure state indicated network regularization during seizures. When comparing between strains, intrathalamic connections were prominent in NEC, on local level represented by higher thalamic strengths and hub scores. Subtle differences were observed for retrosplenial cortex (RS), forming more connections beyond cortex in epileptic rats, and showing a tendency to lateralization during seizures. A potential role of RS as hub between subcortical and cortical regions in epilepsy was supported by increased numbers of parvalbumin-positive (PV+) interneurons together with enhanced inhibitory synaptic activity and neuronal excitability in pyramidal neurons. By combining multimodal fMRI data, graph theoretical methods, and electrophysiological recordings, we identified the RS as promising target for modulation of seizure activity and/or comorbidities.

A cortical rat hemodynamic response function for improved detection of BOLD activation under common experimental conditions.

For a reliable estimation of neuronal activation based on BOLD fMRI measurements an accurate model of the hemodynamic response is essential. Since a large part of basic neuroscience research is based on small animal data, it is necessary to characterize a hemodynamic response function (HRF) which is optimized for small animals. Therefore, we have determined and investigated the HRFs of rats obtained under a variety of experimental conditions in the primary somatosensory cortex. Measurements were performed on animals of different sex and strain, under different anesthetics, with and without ventilation and using different stimulation modalities. All modalities of stimulation used in this study induced neuronal activity in the primary somatosensory cortex or in subcortical regions. Since the HRFs of the BOLD responses in the primary somatosensory cortex showed a close concordance for the different conditions, we were able to determine a cortical rat HRF. This HRF is based on 143 BOLD measurements of 76 rats and can be used for statistical parametric mapping. It showed substantially faster progression than the human HRF, with a maximum after 2.8 ± 0.8 s, and a following undershoot after 6.1 ± 3.7 s. If the rat HRF was used statistical analysis of rat data showed a significantly improved detection performance in the somatosensory cortex in comparison to the commonly used HRF based on measurements in humans.

Functional MRI Readouts From BOLD and Diffusion Measurements Differentially Respond to Optogenetic Activation and Tissue Heating.

Functional blood-oxygenation-level-dependent (BOLD) MRI provides a brain-wide readout that depends on the hemodynamic response to neuronal activity. Diffusion fMRI has been proposed as an alternative to BOLD fMRI and has been postulated to directly rely on neuronal activity. These complementary functional readouts are versatile tools to be combined with optogenetic stimulation to investigate networks of the brain. The cell-specificity and temporal precision of optogenetic manipulations promise to enable further investigation of the origin of fMRI signals. The signal characteristics of the diffusion fMRI readout vice versa may better resolve network effects of optogenetic stimulation. However, the light application needed for optogenetic stimulation is accompanied by heat deposition within the tissue. As both diffusion and BOLD are sensitive to temperature changes, light application can lead to apparent activations confounding the interpretation of fMRI data. The degree of tissue heating, the appearance of apparent activation in different fMRI sequences and the origin of these phenomena are not well understood. Here, we disentangled apparent activations in BOLD and diffusion measurements in rats from physiological activation upon sensory or optogenetic stimulation. Both, BOLD and diffusion fMRI revealed similar signal shapes upon sensory stimulation that differed clearly from those upon heating. Apparent activations induced by high-intensity light application were dominated by T2 ∗-effects and resulted in mainly negative signal changes. We estimated that even low-intensity light application used for optogenetic stimulation reduces the BOLD response close to the fiber by up to 0.4%. The diffusion fMRI signal contained T2, T2 ∗ and diffusion components. The apparent diffusion coefficient, which reflects the isolated diffusion component, showed negative changes upon both optogenetic and electric forepaw stimulation. In contrast, positive changes were detected upon high-intensity light application and thus ruled out heating as a major contributor to the diffusion fMRI signal.