High Spatiotemporal Resolution Radial Encoding Single-Vessel fMRI.

High-field preclinical functional MRI (fMRI) is enabled the high spatial resolution mapping of vessel-specific hemodynamic responses, that is single-vessel fMRI. In contrast to investigating the neuronal sources of the fMRI signal, single-vessel fMRI focuses on elucidating its vascular origin, which can be readily implemented to identify vascular changes relevant to vascular dementia or cognitive impairment. However, the limited spatial and temporal resolution of fMRI is hindered hemodynamic mapping of intracortical microvessels. Here, the radial encoding MRI scheme is implemented to measure BOLD signals of individual vessels penetrating the rat somatosensory cortex. Radial encoding MRI is employed to map cortical activation with a focal field of view (FOV), allowing vessel-specific functional mapping with 50 × 50 µm2 in-plane resolution at a 1 to 2 Hz sampling rate. Besides detecting refined hemodynamic responses of intracortical micro-venules, the radial encoding-based single-vessel fMRI enables the distinction of fMRI signals from vessel and peri-vessel voxels due to the different contribution of intravascular and extravascular effects.

Characterizing brain stage-dependent pupil dynamics based on lateral hypothalamic activity.

Pupil dynamics presents varied correlation features with brain activity under different vigilant levels. The modulation of brain dynamic stages can arise from the lateral hypothalamus (LH), where diverse neuronal cell types contribute to arousal regulation in opposite directions via the anterior cingulate cortex (ACC). However, the relationship of the LH and pupil dynamics has seldom been investigated. Here, we performed local field potential (LFP) recordings at the LH and ACC, and whole-brain fMRI with simultaneous fiber photometry Ca2+ recording in the ACC, to evaluate their correlation with brain state-dependent pupil dynamics. Both LFP and functional magnetic resonance imaging (fMRI) data showed various correlations to pupil dynamics across trials that span negative, null, and positive correlation values, demonstrating brain state-dependent coupling features. Our results indicate that the correlation of pupil dynamics with ACC LFP and whole-brain fMRI signals depends on LH activity, suggesting a role of the latter in brain dynamic stage regulation.

Cortical depth-dependent human fMRI of resting-state networks using EPIK.

Introduction: Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear. Methods: Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers. Results: Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest). Discussion: The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

Mapping of whole-cerebrum resting-state networks using ultra-high resolution acquisition protocols.

Resting-state functional magnetic resonance imaging (fMRI) has been used in numerous studies to map networks in the brain that employ spatially disparate regions. However, attempts to map networks with high spatial resolution have been hampered by conflicting technical demands and associated problems. Results from recent fMRI studies have shown that spatial resolution remains around 0.7 × 0.7 × 0.7 mm3 , with only partial brain coverage. Therefore, this work aims to present a novel fMRI technique that was developed based on echo-planar-imaging with keyhole (EPIK) combined with repetition-time-external (TR-external) EPI phase correction. Each technique has been previously shown to be effective in enhancing the spatial resolution of fMRI, and in this work, the combination of the two techniques into TR-external EPIK provided a nominal spatial resolution of 0.51 × 0.51 × 1.00 mm3 (0.26 mm3 voxel) with whole-cerebrum coverage. Here, the feasibility of using half-millimetre in-plane TR-external EPIK for resting-state fMRI was validated using 13 healthy subjects and the corresponding reproducible mapping of resting-state networks was demonstrated. Furthermore, TR-external EPIK enabled the identification of various resting-state networks distributed throughout the brain from a single fMRI session, with mapping fidelity onto the grey matter at 7T. The high-resolution functional image further revealed mesoscale anatomical structures, such as small cerebral vessels and the internal granular layer of the cortex within the postcentral gyrus.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications.

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Contribution of animal models toward understanding resting state functional connectivity.

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

Decoding the brain state-dependent relationship between pupil dynamics and resting state fMRI signal fluctuation.

Pupil dynamics serve as a physiological indicator of cognitive processes and arousal states of the brain across a diverse range of behavioral experiments. Pupil diameter changes reflect brain state fluctuations driven by neuromodulatory systems. Resting-state fMRI (rs-fMRI) has been used to identify global patterns of neuronal correlation with pupil diameter changes; however, the linkage between distinct brain state-dependent activation patterns of neuromodulatory nuclei with pupil dynamics remains to be explored. Here, we identified four clusters of trials with unique activity patterns related to pupil diameter changes in anesthetized rat brains. Going beyond the typical rs-fMRI correlation analysis with pupil dynamics, we decomposed spatiotemporal patterns of rs-fMRI with principal component analysis (PCA) and characterized the cluster-specific pupil-fMRI relationships by optimizing the PCA component weighting via decoding methods. This work shows that pupil dynamics are tightly coupled with different neuromodulatory centers in different trials, presenting a novel PCA-based decoding method to study the brain state-dependent pupil-fMRI relationship.

Mapping the Brain-Wide Network Effects by Optogenetic Activation of the Corpus Callosum.

Optogenetically driven manipulation of circuit-specific activity enables causality studies, but its global brain-wide effect is rarely reported. Here, we applied simultaneous functional magnetic resonance imaging (fMRI) and calcium recording with optogenetic activation of the corpus callosum (CC) connecting barrel cortices (BC). Robust positive BOLD was detected in the ipsilateral BC due to antidromic activity, spreading to the ipsilateral motor cortex (MC), and posterior thalamus (PO). In the orthodromic target, positive BOLD was reliably evoked by 2 Hz light pulses, whereas 40 Hz light pulses led to reduced calcium, indicative of CC-mediated inhibition. This presumed optogenetic CC-mediated inhibition was further elucidated by pairing light pulses with whisker stimulation at varied interstimulus intervals. Whisker-induced positive BOLD and calcium signals were reduced at intervals of 50/100 ms. The calcium-amplitude-modulation-based correlation with whole-brain fMRI signal revealed that the inhibitory effects spread to contralateral BC, ipsilateral MC, and PO. This work raises the need for fMRI to elucidate the brain-wide network activation in response to optogenetic stimulation.

Indexing brain state-dependent pupil dynamics with simultaneous fMRI and optical fiber calcium recording.

Pupillometry, a noninvasive measure of arousal, complements human functional MRI (fMRI) to detect periods of variable cognitive processing and identify networks that relate to particular attentional states. Even under anesthesia, pupil dynamics correlate with brain-state fluctuations, and extended dilations mark the transition to more arousable states. However, cross-scale neuronal activation patterns are seldom linked to brain state-dependent pupil dynamics. Here, we complemented resting-state fMRI in rats with cortical calcium recording (GCaMP-mediated) and pupillometry to tackle the linkage between brain-state changes and neural dynamics across different scales. This multimodal platform allowed us to identify a global brain network that covaried with pupil size, which served to generate an index indicative of the brain-state fluctuation during anesthesia. Besides, a specific correlation pattern was detected in the brainstem, at a location consistent with noradrenergic cell group 5 (A5), which appeared to be dependent on the coupling between different frequencies of cortical activity, possibly further indicating particular brain-state dynamics. The multimodal fMRI combining concurrent calcium recordings and pupillometry enables tracking brain state-dependent pupil dynamics and identifying unique cross-scale neuronal dynamic patterns under anesthesia.

MRI-guided robotic arm drives optogenetic fMRI with concurrent Ca2+ recording.

Optical fiber-mediated optogenetic activation and neuronal Ca2+ recording in combination with fMRI provide a multi-modal fMRI platform. Here, we developed an MRI-guided robotic arm (MgRA) as a flexible positioning system with high precision to real-time assist optical fiber brain intervention for multi-modal animal fMRI. Besides the ex vivo precision evaluation, we present the highly reliable brain activity patterns in the projected basal forebrain regions upon MgRA-driven optogenetic stimulation in the lateral hypothalamus. Also, we show the step-wise optical fiber targeting thalamic nuclei and map the region-specific functional connectivity with whole-brain fMRI accompanied by simultaneous calcium recordings to specify its circuit-specificity. The MgRA also guides the real-time microinjection to specific deep brain nuclei, which is demonstrated by an Mn-enhanced MRI method. The MgRA represents a clear advantage over the standard stereotaxic-based fiber implantation and opens a broad avenue to investigate the circuit-specific functional brain mapping with the multi-modal fMRI platform.

Multimodal assessment of recovery from coma in a rat model of diffuse brainstem tegmentum injury.

Despite the association between brainstem lesions and coma, a mechanistic understanding of coma pathogenesis and recovery is lacking. We developed a coma model in the rat mimicking human brainstem coma, which allowed multimodal analysis of a brainstem tegmentum lesion's effects on behavior, cortical electrophysiology, and global brain functional connectivity. After coma induction, we observed a transient period (∼1h) of unresponsiveness accompanied by cortical burst-suppression. Comatose rats then gradually regained behavioral responsiveness concurrent with emergence of delta/theta-predominant cortical rhythms in primary somatosensory cortex. During the acute stage of coma recovery (∼1-8h), longitudinal resting-state functional MRI revealed an increase in functional connectivity between subcortical arousal nuclei in the thalamus, basal forebrain, and basal ganglia and cortical regions implicated in awareness. This rat coma model provides an experimental platform to systematically study network-based mechanisms of coma pathogenesis and recovery, as well as to test targeted therapies aimed at promoting recovery of consciousness after coma.

Identifying Respiration-Related Aliasing Artifacts in the Rodent Resting-State fMRI.

Resting-state functional magnetic resonance imaging (rs-fMRI) combined with optogenetics and electrophysiological/calcium recordings in animal models is becoming a popular platform to investigate brain dynamics under specific neurological states. Physiological noise originating from the cardiac and respiration signal is the dominant interference in human rs-fMRI and extensive efforts have been made to reduce these artifacts from the human data. In animal fMRI studies, physiological noise sources including the respiratory and cardiorespiratory artifacts to the rs-fMRI signal fluctuation have typically been less investigated. In this article, we demonstrate evidence of aliasing effects into the low-frequency rs-fMRI signal fluctuation mainly due to respiration-induced B0 offsets in anesthetized rats. This aliased signal was examined by systematically altering the fMRI sampling rate, i.e., the time of repetition (TR), in free-breathing conditions and by adjusting the rate of ventilation. Anesthetized rats under ventilation showed a significantly narrower frequency bandwidth of the aliasing effect than free-breathing animals. It was found that the aliasing effect could be further reduced in ventilated animals with a muscle relaxant. This work elucidates the respiration-related aliasing effects on the rs-fMRI signal fluctuation from anesthetized rats, indicating non-negligible physiological noise needed to be taken care of in both awake and anesthetized animal rs-fMRI studies.

High magnetic field induced otolith fusion in the zebrafish larvae.

Magnetoreception in animals illustrates the interaction of biological systems with the geomagnetic field (geoMF). However, there are few studies that identified the impact of high magnetic field (MF) exposure from Magnetic Resonance Imaging (MRI) scanners (>100,000 times of geoMF) on specific biological targets. Here, we investigated the effects of a 14 Tesla MRI scanner on zebrafish larvae. All zebrafish larvae aligned parallel to the B0 field, i.e. the static MF, in the MRI scanner. The two otoliths (ear stones) in the otic vesicles of zebrafish larvae older than 24 hours post fertilization (hpf) fused together after the high MF exposure as short as 2 hours, yielding a single-otolith phenotype with aberrant swimming behavior. The otolith fusion was blocked in zebrafish larvae under anesthesia or embedded in agarose. Hair cells may play an important role on the MF-induced otolith fusion. This work provided direct evidence to show that high MF interacts with the otic vesicle of zebrafish larvae and causes otolith fusion in an "all-or-none" manner. The MF-induced otolith fusion may facilitate the searching for MF sensors using genetically amenable vertebrate animal models, such as zebrafish.