Sensory evoked fMRI paradigms in awake mice.

Mouse fMRI has become increasingly popular in the small animal imaging field. However, compared to the more commonly used rat fMRI, it is challenging for mouse fMRI to obtain robust and specific functional imaging results. In the meantime, in other neuroscience modalities such as optical imaging, functional recording in the awake mice is common and becoming standard. Therefore, in the current study we developed comprehensive setups and analysis pipeline for multi-sensory fMRI paradigms in the awake mice. Customized setups of somatosensory (whisker), auditory and olfactory stimulation were developed for use in the cryogenic coil in the awake mouse fMRI setting. After carefully evaluating head motion and motion artefacts, the nuisance regression approach was optimized for reducing the confounding effect of motion. The high temporal resolution data (TR = 0.35 s) revealed fast temporal dynamics (time-to-peak ~2 s) of evoked BOLD responses in most brain regions. Using the derived awake mouse specific hemodynamic response functions, high spatial resolution data revealed robust, specific and consistent cortical and subcortical activations in response to somatosensory, auditory and olfactory stimulations, respectively. Overall, we present comprehensive methods for acquiring and analyzing sensory evoked awake mouse fMRI data. The establishment of multi-sensory paradigms in awake mouse fMRI provides valuable tools for examining spatiotemporal characteristics and neural mechanisms of BOLD signals in the future.

Differential coupling between subcortical calcium and BOLD signals during evoked and resting state through simultaneous calcium fiber photometry and fMRI.

Task based and resting state fMRI has been widely utilized to study brain functions. As the foundation of fMRI, the underlying neural basis of the BOLD signal has been extensively studied, but the detailed mechanism remains elusive, particularly during the resting state. To examine the neurovascular coupling, it is important to simultaneously record neural and vascular signals. Here we developed a novel setup of camera based, scalable simultaneous calcium fiber photometry and fMRI in rats. Using this setup, we recorded calcium signals of superior colliculus (SC) and lateral geniculate nucleus (LGN) and fMRI simultaneously during visual stimulation and the resting state. Our results revealed robust, region-specific coupling between calcium and BOLD signals in the task state and weaker, whole brain correlation in the resting state. Interestingly, the spatial specificity of such correlation in the resting state was improved upon regression of white matter, ventricle signals and global signals in fMRI data. Overall, our results suggest differential coupling of calcium and BOLD signals for subcortical regions between evoked and resting states, and the coupling relationship in the resting state was related with resting state BOLD preprocessing strategies.

Awake and behaving mouse fMRI during Go/No-Go task.

Functional magnetic imaging (fMRI) has been widely used to examine the functional neural networks in both the evoked and resting states. However, most fMRI studies in rodents are performed under anesthesia, which greatly limits the scope of their application, and behavioral relevance. Efforts have been made to image rodents in the awake condition, either in the resting state or in response to sensory or optogenetic stimulation. However, fMRI in awake behaving rodents has not yet been achieved. In the current study, a novel fMRI paradigm for awake and behaving mice was developed, allowing functional imaging of the mouse brain in an olfaction-based go/no-go task. High resolution functional imaging with limited motion and image distortion were achieved at 9.4T with a cryogenic coil in awake and behaving mice. Distributed whole-brain spatiotemporal patterns were revealed, with drastically different activity profiles for go versus no-go trials. Therefore, we have demonstrated the feasibility of functional imaging of an olfactory behavior in awake mice. This fMRI paradigm in awake behaving mice could lead to novel insights into neural mechanisms underlying behaviors at a whole-brain level.

Aberrant resting-state co-activation network dynamics in major depressive disorder.

Major depressive disorder (MDD) is a globally prevalent and highly disabling disease characterized by dysfunction of large-scale brain networks. Previous studies have found that static functional connectivity is not sufficient to reflect the complicated and time-varying properties of the brain. The underlying dynamic interactions between brain functional networks of MDD remain largely unknown, and it is also unclear whether neuroimaging-based dynamic properties are sufficiently robust to discriminate individuals with MDD from healthy controls since the diagnosis of MDD mainly depends on symptom-based criteria evaluated by clinical observation. Resting-state functional magnetic resonance imaging (fMRI) data of 221 MDD patients and 215 healthy controls were shared by REST-meta-MDD consortium. We investigated the spatial-temporal dynamics of MDD using co-activation pattern analysis and made individual diagnoses using support vector machine (SVM). We found that MDD patients exhibited aberrant dynamic properties (such as dwell time, occurrence rate, transition probability, and entropy of Markov trajectories) in some transient networks including subcortical network (SCN), activated default mode network (DMN), de-activated SCN-cerebellum network, a joint network, activated attention network (ATN), and de-activated DMN-ATN, where some dynamic properties were indicative of depressive symptoms. The trajectories of other networks to deactivated DMN-ATN were more accessible in MDD patients. Subgroup analyses also showed subtle dynamic changes in first-episode drug-naïve (FEDN) MDD patients. Finally, SVM achieved preferable accuracies of 84.69%, 76.77%, and 88.10% in discriminating patients with MDD, FEDN MDD, and recurrent MDD from healthy controls with their dynamic metrics. Our findings reveal that MDD is characterized by aberrant dynamic fluctuations of brain network and the feasibility of discriminating MDD patients using dynamic properties, which provide novel insights into the neural mechanism of MDD.

Cell-type-specific optogenetic fMRI on basal forebrain reveals functional network basis of behavioral preference.

The basal forebrain (BF) is a complex structure that plays key roles in regulating various brain functions. However, it remains unclear how cholinergic and non-cholinergic BF neurons modulate large-scale functional networks and their relevance in intrinsic and extrinsic behaviors. With an optimized awake mouse optogenetic fMRI approach, we revealed that optogenetic stimulation of four BF neuron types evoked distinct cell-type-specific whole-brain BOLD activations, which could be attributed to BF-originated low-dimensional structural networks. Additionally, optogenetic activation of VGLUT2, ChAT, and PV neurons in the BF modulated the preference for locomotion, exploration, and grooming, respectively. Furthermore, we uncovered the functional network basis of the above BF-modulated behavioral preference through a decoding model linking the BF-modulated BOLD activation, low-dimensional structural networks, and behavioral preference. To summarize, we decoded the functional network basis of differential behavioral preferences with cell-type-specific optogenetic fMRI on the BF and provided an avenue for investigating mouse behaviors from a whole-brain view.

Instantaneous antidepressant effect of lateral habenula deep brain stimulation in rats studied with functional MRI.

The available treatments for depression have substantial limitations, including low response rates and substantial lag time before a response is achieved. We applied deep brain stimulation (DBS) to the lateral habenula (LHb) of two rat models of depression (Wistar Kyoto rats and lipopolysaccharide-treated rats) and observed an immediate (within seconds to minutes) alleviation of depressive-like symptoms with a high-response rate. Simultaneous functional MRI (fMRI) conducted on the same sets of depressive rats used in behavioral tests revealed DBS-induced activation of multiple regions in afferent and efferent circuitry of the LHb. The activation levels of brain regions connected to the medial LHb (M-LHb) were correlated with the extent of behavioral improvements. Rats with more medial stimulation sites in the LHb exhibited greater antidepressant effects than those with more lateral stimulation sites. These results indicated that the antidromic activation of the limbic system and orthodromic activation of the monoaminergic systems connected to the M-LHb played a critical role in the rapid antidepressant effects of LHb-DBS. This study indicates that M-LHb-DBS might act as a valuable, rapid-acting antidepressant therapeutic strategy for treatment-resistant depression and demonstrates the potential of using fMRI activation of specific brain regions as biomarkers to predict and evaluate antidepressant efficacy.

Sleep fMRI with simultaneous electrophysiology at 9.4 T in male mice.

Sleep is ubiquitous and essential, but its mechanisms remain unclear. Studies in animals and humans have provided insights of sleep at vastly different spatiotemporal scales. However, challenges remain to integrate local and global information of sleep. Therefore, we developed sleep fMRI based on simultaneous electrophysiology at 9.4 T in male mice. Optimized un-anesthetized mouse fMRI setup allowed manifestation of NREM and REM sleep, and a large sleep fMRI dataset was collected and openly accessible. State dependent global patterns were revealed, and state transitions were found to be global, asymmetrical and sequential, which can be predicted up to 17.8 s using LSTM models. Importantly, sleep fMRI with hippocampal recording revealed potentiated sharp-wave ripple triggered global patterns during NREM than awake state, potentially attributable to co-occurrence of spindle events. To conclude, we established mouse sleep fMRI with simultaneous electrophysiology, and demonstrated its capability by revealing global dynamics of state transitions and neural events.

Multimodal analysis demonstrating the shaping of functional gradients in the marmoset brain.

The discovery of functional gradients introduce a new perspective in understanding the cortical spectrum of intrinsic dynamics, as it captures major axes of functional connectivity in low-dimensional space. However, how functional gradients arise and dynamically vary remains poorly understood. In this study, we investigated the biological basis of functional gradients using awake resting-state fMRI, retrograde tracing and gene expression datasets in marmosets. We found functional gradients in marmosets showed a sensorimotor-to-visual principal gradient followed by a unimodal-to-multimodal gradient, resembling functional gradients in human children. Although strongly constrained by structural wirings, functional gradients were dynamically modulated by arousal levels. Utilizing a reduced model, we uncovered opposing effects on gradient dynamics by structural connectivity (inverted U-shape) and neuromodulatory input (U-shape) with arousal fluctuations, and dissected the contribution of individual neuromodulatory receptors. This study provides insights into biological basis of functional gradients by revealing the interaction between structural connectivity and ascending neuromodulatory system.

Automated Skull Stripping in Mouse Functional Magnetic Resonance Imaging Analysis Using 3D U-Net.

Skull stripping is an initial and critical step in the pipeline of mouse fMRI analysis. Manual labeling of the brain usually suffers from intra- and inter-rater variability and is highly time-consuming. Hence, an automatic and efficient skull-stripping method is in high demand for mouse fMRI studies. In this study, we investigated a 3D U-Net based method for automatic brain extraction in mouse fMRI studies. Two U-Net models were separately trained on T2-weighted anatomical images and T2*-weighted functional images. The trained models were tested on both interior and exterior datasets. The 3D U-Net models yielded a higher accuracy in brain extraction from both T2-weighted images (Dice > 0.984, Jaccard index > 0.968 and Hausdorff distance < 7.7) and T2*-weighted images (Dice > 0.964, Jaccard index > 0.931 and Hausdorff distance < 3.3), compared with the two widely used mouse skull-stripping methods (RATS and SHERM). The resting-state fMRI results using automatic segmentation with the 3D U-Net models are highly consistent with those obtained by manual segmentation for both the seed-based and group independent component analysis. These results demonstrate that the 3D U-Net based method can replace manual brain extraction in mouse fMRI analysis.

An integrated resource for functional and structural connectivity of the marmoset brain.

Comprehensive integration of structural and functional connectivity data is required to model brain functions accurately. While resources for studying the structural connectivity of non-human primate brains already exist, their integration with functional connectivity data has remained unavailable. Here we present a comprehensive resource that integrates the most extensive awake marmoset resting-state fMRI data available to date (39 marmoset monkeys, 710 runs, 12117 mins) with previously published cellular-level neuronal tracing data (52 marmoset monkeys, 143 injections) and multi-resolution diffusion MRI datasets. The combination of these data allowed us to (1) map the fine-detailed functional brain networks and cortical parcellations, (2) develop a deep-learning-based parcellation generator that preserves the topographical organization of functional connectivity and reflects individual variabilities, and (3) investigate the structural basis underlying functional connectivity by computational modeling. This resource will enable modeling structure-function relationships and facilitate future comparative and translational studies of primate brains.

Full activation pattern mapping by simultaneous deep brain stimulation and fMRI with graphene fiber electrodes.

Simultaneous deep brain stimulation (DBS) and functional magnetic resonance imaging (fMRI) constitutes a powerful tool for elucidating brain functional connectivity, and exploring neuromodulatory mechanisms of DBS therapies. Previous DBS-fMRI studies could not provide full activation pattern maps due to poor MRI compatibility of the DBS electrodes, which caused obstruction of large brain areas on MRI scans. Here, we fabricate graphene fiber (GF) electrodes with high charge-injection-capacity and little-to-no MRI artifact at 9.4T. DBS-fMRI with GF electrodes at the subthalamic nucleus (STN) in Parkinsonian rats reveal robust blood-oxygenation-level-dependent responses along the basal ganglia-thalamocortical network in a frequency-dependent manner, with responses from some regions not previously detectable. This full map indicates that STN-DBS modulates both motor and non-motor pathways, possibly through orthodromic and antidromic signal propagation. With the capability for full, unbiased activation pattern mapping, DBS-fMRI using GF electrodes can provide important insights into DBS therapeutic mechanisms in various neurological disorders.