Functional connectivity in blood oxygenation level-dependent and cerebral blood volume-weighted resting state functional magnetic resonance imaging in the rat brain

Coupled pulsatile vascular and paravascular fluid dynamics in the human brain

BACKGROUND

Cardiac pulsation propels blood through the cerebrovascular network to maintain cerebral homeostasis. The cerebrovascular network is uniquely surrounded by paravascular cerebrospinal fluid (pCSF), which plays a crucial role in waste removal, and its flow is suspected to be driven by arterial pulsations. Despite its importance, the relationship between vascular and paravascular fluid dynamics throughout the cardiac cycle remains poorly understood in humans.

METHODS

In this study, we developed a non-invasive neuroimaging approach to investigate the coupling between pulsatile vascular and pCSF dynamics within the subarachnoid space of the human brain. Resting-state functional MRI (fMRI) and dynamic diffusion-weighted imaging (dynDWI) were retrospectively cardiac-aligned to represent cerebral hemodynamics and pCSF motion, respectively. We measured the time between peaks (∆TTP) ind d ϕ f M R I and dynDWI waveforms and measured their coupling by calculating the waveforms correlation after peak alignment (correlation at aligned peaks). We compared the ∆TTP and correlation at aligned peaks between younger [mean age: 27.9 (3.3) years, n = 9] and older adults [mean age: 70.5 (6.6) years, n = 20], and assessed their reproducibility within subjects and across different imaging protocols.

RESULTS

Hemodynamic changes consistently precede pCSF motion. ∆TTP was significantly shorter in younger adults compared to older adults (-0.015 vs. -0.069, p < 0.05). The correlation at aligned peaks were high and did not differ between younger and older adults (0.833 vs. 0.776, p = 0.153). The ∆TTP and correlation at aligned peaks were robust across fMRI protocols (∆TTP: -0.15 vs. -0.053, p = 0.239; correlation at aligned peaks: 0.813 vs. 0.812, p = 0.985) and demonstrated good to excellent within-subject reproducibility (∆TTP: intraclass correlation coefficient = 0.36; correlation at aligned peaks: intraclass correlation coefficient = 0.89).

CONCLUSION

This study proposes a non-invasive technique to evaluate vascular and paravascular fluid dynamics. Our findings reveal a consistent and robust cardiac pulsation-driven coupling between cerebral hemodynamics and pCSF dynamics in both younger and older adults.

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.

3D bioprinting of nanoparticle-laden hydrogel scaffolds with enhanced antibacterial and imaging properties

Biomaterial-associated microbial contaminations in biologically conducive three-dimensional (3D) tissue-engineered constructs have significantly limited the clinical applications of scaffold systems. To prevent such infections, antimicrobial biomaterials are rapidly evolving. Yet, the use of such materials in bioprinting-based approaches of scaffold fabrication has not been examined. This study introduces a new generation of bacteriostatic gelatin methacryloyl (GelMA)-based bioinks, incorporated with varying doses of antibacterial superparamagnetic iron oxide nanoparticles (SPIONs). The SPION-laden GelMA scaffolds showed significant resistance against the Staphylococcus aureus growth, while providing a contrast in magnetic resonance imaging. We simulated the bacterial contamination of cellular 3D GelMA scaffolds in vitro and demonstrated the significant effect of functionalized scaffolds in inhibiting bacterial growth, while maintaining cell viability and growth. Together, these results present a new promising class of functionalized bioinks to 3D bioprint tissue-engineered scaffold with markedly enhanced properties for the use in a variety of in vitro and clinical applications.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

Brain networks of rats under anesthesia using resting-state fMRI: comparison with dead rats, random noise and generative models of networks

OBJECTIVE

Connectivity networks are crucial to understand the brain resting-state activity using functional magnetic resonance imaging (rs-fMRI). Alterations of these brain networks may highlight important findings concerning the resilience of the brain to different disorders. The focus of this paper is to evaluate the robustness of brain network estimations, discriminate them under anesthesia and compare them to generative models.

APPROACH

The extraction of brain functional connectivity (FC) networks is difficult and biased due to the properties of the data: low signal to noise ratio, high dimension low sample size. We propose to use wavelet correlations to assess FC between brain areas under anesthesia using four anesthetics (isoflurane, etomidate, medetomidine, urethane). The networks are then deduced from the functional connectivity matrices by applying statistical thresholds computed using the number of samples at a given scale of wavelet decomposition. Graph measures are extracted and extensive comparisons with generative models of structured networks are conducted.

MAIN RESULTS

The sample size and filtering are critical to obtain significant correlations values and thereby detect connections between regions. This is necessary to construct networks different from random ones as shown using rs-fMRI brain networks of dead rats. Brain networks under anesthesia on rats have topological features that are mixing small-world, scale-free and random networks. Betweenness centrality indicates that hubs are present in brain networks obtained from anesthetized rats but locations of these hubs are altered by anesthesia.

SIGNIFICANCE

Understanding the effects of anesthesia on brain areas is of particular importance in the context of animal research since animal models are commonly used to explore functions, evaluate lesions or illnesses, and test new drugs. More generally, results indicate that the use of correlations in the context of fMRI signals is robust but must be treated with caution. Solutions are proposed in order to control spurious correlations by setting them to zero.

Sex differences in insular functional connectivity in response to noxious visceral stimulation in rats

Insular cortex (INS) plays a critical role in pain processing and shows sex differences in functional activation during noxious visceral stimulation. Less is known regarding functional interactions within the INS and between this structure and other parts of the brain. Cerebral blood flow mapping was performed using [¹⁴C]-iodoantipyrine perfusion autoradiography in male and female rats during colorectal distension (CRD) or no distension (controls). Forty regions of interest (ROIs) were defined anatomically to represent the granular, dysgranular, and agranular INS along the anterior-posterior (A-P) axis. Inter-ROI correlation matrices were calculated for each group to characterize intra-insular functional connectivity (FC). Results showed a clear FC segregation within the INS into an anterior (rostral to bregma +2.4 mm), a posterior (caudal to bregma -1.2 mm), and a mid INS subregion in between. Female controls showed higher FC density compared to males. During CRD, intra-insular FC density decreased greatly in females, but only modestly in males, with a loss of long-range connections between the anterior and mid INS noted in both sexes. New functional organization was characterized in both sexes by a cluster in the mid INS and primarily short-range FC along the A-P axis. Seed correlation analysis during CRD showed sex differences in FC of the anterior and mid agranular INS with the medial prefrontal cortex, thalamus, and brainstem areas (periaqueductal gray, parabrachial nucleus), suggesting sex differences in the modulatory aspect of visceral pain processing. Our findings suggest presence of substantial sex differences in visceral pain processing at the level of the insula.

Comparison of Resting-State Brain Activation Detected by BOLD, Blood Volume and Blood Flow

Resting-state brain activity has been widely investigated using blood oxygenation level dependent (BOLD) contrast techniques. However, BOLD signal changes reflect a combination of the effects of cerebral blood flow (CBF), cerebral blood volume (CBV), as well as the cerebral metabolic rate of oxygen (CMRO₂). In this study, resting-state brain activation was detected and compared using the following techniques: (a) BOLD, using a gradient-echo echo planar imaging (GE-EPI) sequence; (b) CBV-weighted signal, acquired using gradient and spin echo (GRASE) based vascular space occupancy (VASO); and (c) CBF, using pseudo-continuous arterial spin labeling (pCASL). Reliable brain networks were detected using VASO and ASL, including sensorimotor, auditory, primary visual, higher visual, default mode, salience and left/right executive control networks. Differences between the resting-state activation detected with ASL, VASO and BOLD could potentially be due to the different temporal signal-to-noise ratio (tSNR) and the short post-labeling delay (PLD) in ASL, along with differences in the spin-echo readout of VASO. It is also possible that the dynamics of spontaneous fluctuations in BOLD, CBV and CBF could differ due to biological reasons, according to their location within the brain.

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.

Dynamic resting state fMRI analysis in mice reveals a set of Quasi-Periodic Patterns and illustrates their relationship with the global signal

Time-resolved 'dynamic' over whole-period 'static' analysis of low frequency (LF) blood-oxygen level dependent (BOLD) fluctuations provides many additional insights into the macroscale organization and dynamics of neural activity. Although there has been considerable advancement in the development of mouse resting state fMRI (rsfMRI), very little remains known about its dynamic repertoire. Here, we report for the first time the detection of a set of recurring spatiotemporal Quasi-Periodic Patterns (QPPs) in mice, which show spatial similarity with known resting state networks. Furthermore, we establish a close relationship between several of these patterns and the global signal. We acquired high temporal rsfMRI scans under conditions of low (LA) and high (HA) medetomidine-isoflurane anesthesia. We then employed the algorithm developed by Majeed et al. (2011), previously applied in rats and humans, which detects and averages recurring spatiotemporal patterns in the LF BOLD signal. One type of observed patterns in mice was highly similar to those originally observed in rats, displaying propagation from lateral to medial cortical regions, which suggestively pertain to a mouse Task-Positive like network (TPN) and Default Mode like network (DMN). Other QPPs showed more widespread or striatal involvement and were no longer detected after global signal regression (GSR). This was further supported by diminished detection of subcortical dynamics after GSR, with cortical dynamics predominating. Observed QPPs were both qualitatively and quantitatively determined to be consistent across both anesthesia conditions, with GSR producing the same outcome. Under LA, QPPs were consistently detected at both group and single subject level. Under HA, consistency and pattern occurrence rate decreased, whilst cortical contribution to the patterns diminished. These findings confirm the robustness of QPPs across species and demonstrate a new approach to study mouse LF BOLD spatiotemporal dynamics and mechanisms underlying functional connectivity. The observed impact of GSR on QPPs might help better comprehend its controversial role in conventional resting state studies. Finally, consistent detection of QPPs at single subject level under LA promises a step forward towards more reliable mouse rsfMRI and further confirms the importance of selecting an optimal anesthesia regime.

Identifying Rodent Resting-State Brain Networks with Independent Component Analysis

Rodent models have opened the door to a better understanding of the neurobiology of brain disorders and increased our ability to evaluate novel treatments. Resting-state functional magnetic resonance imaging (rs-fMRI) allows for in vivo exploration of large-scale brain networks with high spatial resolution. Its application in rodents affords researchers a powerful translational tool to directly assess/explore the effects of various pharmacological, lesion, and/or disease states on known neural circuits within highly controlled settings. Integration of animal and human research at the molecular-, systems-, and behavioral-levels using diverse neuroimaging techniques empowers more robust interrogations of abnormal/ pathological processes, critical for evolving our understanding of neuroscience. We present a comprehensive protocol to evaluate resting-state brain networks using Independent Component Analysis (ICA) in rodent model. Specifically, we begin with a brief review of the physiological basis for rs-fMRI technique and overview of rs-fMRI studies in rodents to date, following which we provide a robust step-by-step approach for rs-fMRI investigation including data collection, computational preprocessing, and brain network analysis. Pipelines are interwoven with underlying theory behind each step and summarized methodological considerations, such as alternative methods available and current consensus in the literature for optimal results. The presented protocol is designed in such a way that investigators without previous knowledge in the field can implement the analysis and obtain viable results that reliably detect significant differences in functional connectivity between experimental groups. Our goal is to empower researchers to implement rs-fMRI in their respective fields by incorporating technical considerations to date into a workable methodological framework.

Quasi-periodic patterns of intrinsic brain activity in individuals and their relationship to global signal

Quasiperiodic patterns (QPPs) as reported by Majeed et al., 2011 are prominent features of the brain's intrinsic activity that involve important large-scale networks (default mode, DMN; task positive, TPN) and are likely to be major contributors to widely used measures of functional connectivity. We examined the variability of these patterns in 470 individuals from the Human Connectome Project resting state functional MRI dataset. The QPPs from individuals can be coarsely categorized into two types: one where strong anti-correlation between the DMN and TPN is present, and another where most areas are strongly correlated. QPP type could be predicted by an individual's global signal, with lower global signal corresponding to QPPs with strong anti-correlation. After regression of global signal, all QPPs showed strong anti-correlation between DMN and TPN. QPP occurrence and type was similar between a subgroup of individuals with extremely low motion and the rest of the sample, which shows that motion is not a major contributor to the QPPs. After regression of estimates of slow respiratory and cardiac induced signal fluctuations, more QPPs showed strong anti-correlation between DMN and TPN, an indication that while physiological noise influences the QPP type, it is not the primary source of the QPP itself. QPPs were more similar for the same subjects scanned on different days than for different subjects. These results provide the first assessment of the variability in individual QPPs and their relationship to physiological parameters.

Neural and metabolic basis of dynamic resting state fMRI

Resting state fMRI (rsfMRI) as a technique showed much initial promise for use in psychiatric and neurological diseases where diagnosis and treatment were difficult. To realize this promise, many groups have moved towards examining "dynamic rsfMRI," which relies on the assumption that rsfMRI measurements on short time scales remain relevant to the underlying neural and metabolic activity. Many dynamic rsfMRI studies have demonstrated differences between clinical or behavioral groups beyond what static rsfMRI measured, suggesting a neurometabolic basis. Correlative studies combining dynamic rsfMRI and other physiological measurements have supported this. However, they also indicate multiple mechanisms and, if using correlation alone, it is difficult to separate cause and effect. Hypothesis-driven studies are needed, a few of which have begun to illuminate the underlying neurometabolic mechanisms that shape observed differences in dynamic rsfMRI. While the number of potential noise sources, potential actual neurometabolic sources, and methodological considerations can seem overwhelming, dynamic rsfMRI provides a rich opportunity in systems neuroscience. Even an incrementally better understanding of the neurometabolic basis of dynamic rsfMRI would expand rsfMRI's research and clinical utility, and the studies described herein take the first steps on that path forward.

Functional circuit mapping of striatal output nuclei using simultaneous deep brain stimulation and fMRI

The substantia nigra pars reticulata (SNr) and external globus pallidus (GPe) constitute the two major output targets of the rodent striatum. Both the SNr and GPe converge upon thalamic relay nuclei (directly or indirectly, respectively), and are traditionally modeled as functionally antagonistic relay inputs. However, recent anatomical and functional studies have identified unanticipated circuit connectivity in both the SNr and GPe, demonstrating their potential as far more than relay nuclei. In the present study, we employed simultaneous deep brain stimulation and functional magnetic resonance imaging (DBS-fMRI) with cerebral blood volume (CBV) measurements to functionally and unbiasedly map the circuit- and network level connectivity of the SNr and GPe. Sprague-Dawley rats were implanted with a custom-made MR-compatible stimulating electrode in the right SNr (n=6) or GPe (n=7). SNr- and GPe-DBS, conducted across a wide range of stimulation frequencies, revealed a number of surprising evoked responses, including unexpected CBV decreases within the striatum during DBS at either target, as well as GPe-DBS-evoked positive modulation of frontal cortex. Functional connectivity MRI revealed global modulation of neural networks during DBS at either target, sensitive to stimulation frequency and readily reversed following cessation of stimulation. This work thus contributes to a growing literature demonstrating extensive and unanticipated functional connectivity among basal ganglia nuclei.

Noise and non-neuronal contributions to the BOLD signal: applications to and insights from animal studies

The BOLD signal reflects hemodynamic events within the brain, which in turn are driven by metabolic changes and neural activity. However, the link between BOLD changes and neural activity is indirect and can be influenced by a number of non-neuronal processes. Motion and physiological cycles have long been known to affect the BOLD signal and are present in both humans and animal models. Differences in physiological baseline can also contribute to intra- and inter-subject variability. The use of anesthesia, common in animal studies, alters neural activity, vascular tone, and neurovascular coupling. Most intriguing, perhaps, are the contributions from other processes that do not appear to be neural in origin but which may provide information about other aspects of neurophysiology. This review discusses different types of noise and non-neuronal contributors to the BOLD signal, sources of variability for animal studies, and insights to be gained from animal models.

Probing Intrinsic Resting-State Networks in the Infant Rat Brain

Resting-state functional magnetic resonance imaging (rs-fMRI) measures spontaneous fluctuations in blood oxygenation level-dependent (BOLD) signal in the absence of external stimuli. It has become a powerful tool for mapping large-scale brain networks in humans and animal models. Several rs-fMRI studies have been conducted in anesthetized and awake adult rats, reporting consistent patterns of brain activity at the systems level. However, the evolution to adult patterns of resting-state activity has not yet been evaluated and quantified in the developing rat brain. In this study, we hypothesized that large-scale intrinsic networks would be easily detectable but not fully established as specific patterns of activity in lightly anesthetized 2-week-old rats (N = 11). Independent component analysis (ICA) identified 8 networks in 2-week-old-rats. These included Default mode, Sensory (Exteroceptive), Salience (Interoceptive), Basal Ganglia-Thalamic-Hippocampal, Basal Ganglia, Autonomic, Cerebellar, as well as Thalamic-Brainstem networks. Many of these networks consisted of more than one component, possibly indicative of immature, underdeveloped networks at this early time point. Except for the Autonomic network, infant rat networks showed reduced connectivity with subcortical structures in comparison to previously published adult networks. Reported slow fluctuations in the BOLD signal that correspond to functionally relevant resting-state networks in 2-week-old rats can serve as an important tool for future studies of brain development in the settings of different pharmacological applications or disease.

Cluster-based analysis for characterizing dynamic functional connectivity

Different regions in the resting brain exhibit non-stationary functional connectivity (FC) over time. In this paper, a simple and efficient framework of clustering the variability in FC of a rat's brain at rest is proposed. This clustering process reveals areas that are always connected with a chosen region, called seed voxel, along with the areas exhibiting variability in the FC. This addresses an issue common to most dynamic FC analysis techniques, which is the assumption that the spatial extent of a given network remains constant over time. We increase the voxel size and reduce the spatial resolution to analyze variable FC of the whole resting brain. We hypothesize that the adjacent voxels in resting state functional magnetic resonance imaging (rsfMRI), just as in task-based fMRI, exhibit similar intensities, so they can be averaged to obtain larger voxels without any significant loss of information. Sliding window correlation is used to compute variable patterns of the rat's whole brain FC with the seed voxel in the sensorimotor cortex. These patterns are grouped based on their spatial similarities using binary transformed feature vectors in k-means clustering, not only revealing the variable and nonvariable portions of FC in the resting brain but also detecting the extent of the variability of these patterns.

Wavelet-based clustering of resting state MRI data in the rat

While functional connectivity has typically been calculated over the entire length of the scan (5-10min), interest has been growing in dynamic analysis methods that can detect changes in connectivity on the order of cognitive processes (seconds). Previous work with sliding window correlation has shown that changes in functional connectivity can be observed on these time scales in the awake human and in anesthetized animals. This exciting advance creates a need for improved approaches to characterize dynamic functional networks in the brain. Previous studies were performed using sliding window analysis on regions of interest defined based on anatomy or obtained from traditional steady-state analysis methods. The parcellation of the brain may therefore be suboptimal, and the characteristics of the time-varying connectivity between regions are dependent upon the length of the sliding window chosen. This manuscript describes an algorithm based on wavelet decomposition that allows data-driven clustering of voxels into functional regions based on temporal and spectral properties. Previous work has shown that different networks have characteristic frequency fingerprints, and the use of wavelets ensures that both the frequency and the timing of the BOLD fluctuations are considered during the clustering process. The method was applied to resting state data acquired from anesthetized rats, and the resulting clusters agreed well with known anatomical areas. Clusters were highly reproducible across subjects. Wavelet cross-correlation values between clusters from a single scan were significantly higher than the values from randomly matched clusters that shared no temporal information, indicating that wavelet-based analysis is sensitive to the relationship between areas.

Large-scale functional connectivity networks in the rodent brain

Resting-state functional Magnetic Resonance Imaging (rsfMRI) of the human brain has revealed multiple large-scale neural networks within a hierarchical and complex structure of coordinated functional activity. These distributed neuroanatomical systems provide a sensitive window on brain function and its disruption in a variety of neuropathological conditions. The study of macroscale intrinsic connectivity networks in preclinical species, where genetic and environmental conditions can be controlled and manipulated with high specificity, offers the opportunity to elucidate the biological determinants of these alterations. While rsfMRI methods are now widely used in human connectivity research, these approaches have only relatively recently been back-translated into laboratory animals. Here we review recent progress in the study of functional connectivity in rodent species, emphasising the ability of this approach to resolve large-scale brain networks that recapitulate neuroanatomical features of known functional systems in the human brain. These include, but are not limited to, a distributed set of regions identified in rats and mice that may represent a putative evolutionary precursor of the human default mode network (DMN). The impact and control of potential experimental and methodological confounds are also critically discussed. Finally, we highlight the enormous potential and some initial application of connectivity mapping in transgenic models as a tool to investigate the neuropathological underpinnings of the large-scale connectional alterations associated with human neuropsychiatric and neurological conditions. We conclude by discussing the translational potential of these methods in basic and applied neuroscience.

Considerations for resting state functional MRI and functional connectivity studies in rodents

Resting state functional MRI (rs-fMRI) and functional connectivity mapping have become widely used tools in the human neuroimaging community and their use is rapidly spreading into the realm of rodent research as well. One of the many attractive features of rs-fMRI is that it is readily translatable from humans to animals and back again. Changes in functional connectivity observed in human studies can be followed by more invasive animal experiments to determine the neurophysiological basis for the alterations, while exploratory work in animal models can identify possible biomarkers for further investigation in human studies. These types of interwoven human and animal experiments have a potentially large impact on neuroscience and clinical practice. However, impediments exist to the optimal application of rs-fMRI in small animals, some similar to those encountered in humans and some quite different. In this review we identify the most prominent of these barriers, discuss differences between rs-fMRI in rodents and in humans, highlight best practices for animal studies, and review selected applications of rs-fMRI in rodents. Our goal is to facilitate the integration of human and animal work to the benefit of both fields.

The neural basis of time-varying resting-state functional connectivity

Dynamic network analysis based on resting-state magnetic resonance imaging (rsMRI) is a fairly new and potentially powerful tool for neuroscience and clinical research. Dynamic analysis can be sensitive to changes that occur in psychiatric or neurologic disorders and can detect variations related to performance on individual trials in healthy subjects. However, the appearance of time-varying connectivity can also arise in signals that share no temporal information, complicating the interpretation of dynamic functional connectivity studies. Researchers have begun utilizing simultaneous imaging and electrophysiological recording to elucidate the neural basis of the networks and their variability in animals and in humans. In this article, we review findings that link changes in electrically recorded brain states to changes in the networks obtained with rsMRI and discuss some of the challenges inherent in interpretation of these studies. The literature suggests that multiple brain processes may contribute to the dynamics observed, and we speculate that it may be possible to separate particular aspects of the rsMRI signal to enhance sensitivity to certain types of neural activity, providing new tools for basic neuroscience and clinical research.

Exercise training reinstates cortico-cortical sensorimotor functional connectivity following striatal lesioning: development and application of a subregional-level analytic toolbox for perfusion autoradiographs of the rat brain

Current rodent connectome projects are revealing brain structural connectivity with unprecedented resolution and completeness. How subregional structural connectivity relates to subregional functional interactions is an emerging research topic. We describe a method for standardized, mesoscopic-level data sampling from autoradiographic coronal sections of the rat brain, and for correlation-based analysis and intuitive display of cortico-cortical functional connectivity (FC) on a flattened cortical map. A graphic user interface "Cx-2D" allows for the display of significant correlations of individual regions-of-interest, as well as graph theoretical metrics across the cortex. Cx-2D was tested on an autoradiographic data set of cerebral blood flow (CBF) of rats that had undergone bilateral striatal lesions, followed by 4 weeks of aerobic exercise training or no exercise. Effects of lesioning and exercise on cortico-cortical FC were examined during a locomotor challenge in this rat model of Parkinsonism. Subregional FC analysis revealed a rich functional reorganization of the brain in response to lesioning and exercise that was not apparent in a standard analysis focused on CBF of isolated brain regions. Lesioned rats showed diminished degree centrality of lateral primary motor cortex, as well as neighboring somatosensory cortex-changes that were substantially reversed in lesioned rats following exercise training. Seed analysis revealed that exercise increased positive correlations in motor and somatosensory cortex, with little effect in non-sensorimotor regions such as visual, auditory, and piriform cortex. The current analysis revealed that exercise partially reinstated sensorimotor FC lost following dopaminergic deafferentation. Cx-2D allows for standardized data sampling from images of brain slices, as well as analysis and display of cortico-cortical FC in the rat cerebral cortex with potential applications in a variety of autoradiographic and histologic studies.

Low-frequency calcium oscillations accompany deoxyhemoglobin oscillations in rat somatosensory cortex

Spontaneous low-frequency oscillations (LFOs) of blood-oxygen-level-dependent (BOLD) signals are used to map brain functional connectivity with functional MRI, but their source is not well understood. Here we used optical imaging to assess whether LFOs from vascular signals covary with oscillatory intracellular calcium (Ca(2+)i) and with local field potentials in the rat's somatosensory cortex. We observed that the frequency of Ca(2+)i oscillations in tissue (∼0.07 Hz) was similar to the LFOs of deoxyhemoglobin (HbR) and oxyhemoglobin (HbO2) in both large blood vessels and capillaries. The HbR and HbO2 fluctuations within tissue correlated with Ca(2+)i oscillations with a lag time of ∼5-6 s. The Ca(2+)i and hemoglobin oscillations were insensitive to hypercapnia. In contrast, cerebral-blood-flow velocity (CBFv) in arteries and veins fluctuated at a higher frequency (∼0.12 Hz) and was sensitive to hypercapnia. However, in parenchymal tissue, CBFv oscillated with peaks at both ∼0.06 Hz and ∼0.12 Hz. Although the higher-frequency CBFv oscillation (∼0.12 Hz) was decreased by hypercapnia, its lower-frequency component (∼0.06 Hz) was not. The sensitivity of the higher CBFV oscillations to hypercapnia, which triggers blood vessel vasodilation, suggests its dependence on vascular effects that are distinct from the LFOs detected in HbR, HbO2, Ca(2+)i, and the lower-frequency tissue CBFv, which were insensitive to hypercapnia. Hemodynamic LFOs correlated both with Ca(2+)i and neuronal firing (local field potentials), indicating that they directly reflect neuronal activity (perhaps also glial). These findings show that HbR fluctuations (basis of BOLD oscillations) are linked to oscillatory cellular activity and detectable throughout the vascular tree (arteries, capillaries, and veins).

Functional connectivity-based parcellation and connectome of cortical midline structures in the mouse: a perfusion autoradiography study

Rodent cortical midline structures (CMS) are involved in emotional, cognitive and attentional processes. Tract tracing has revealed complex patterns of structural connectivity demonstrating connectivity-based integration and segregation for the prelimbic, cingulate area 1, retrosplenial dysgranular cortices dorsally, and infralimbic, cingulate area 2, and retrosplenial granular cortices ventrally. Understanding of CMS functional connectivity (FC) remains more limited. Here we present the first subregion-level FC analysis of the mouse CMS, and assess whether fear results in state-dependent FC changes analogous to what has been reported in humans. Brain mapping using [¹⁴C]-iodoantipyrine was performed in mice during auditory-cued fear conditioned recall and in controls. Regional cerebral blood flow (CBF) was analyzed in 3-D images reconstructed from brain autoradiographs. Regions-of-interest were selected along the CMS anterior-posterior and dorsal-ventral axes. In controls, pairwise correlation and graph theoretical analyses showed strong FC within each CMS structure, strong FC along the dorsal-ventral axis, with segregation of anterior from posterior structures. Seed correlation showed FC of anterior regions to limbic/paralimbic areas, and FC of posterior regions to sensory areas–findings consistent with functional segregation noted in humans. Fear recall increased FC between the cingulate and retrosplenial cortices, but decreased FC between dorsal and ventral structures. In agreement with reports in humans, fear recall broadened FC of anterior structures to the amygdala and to somatosensory areas, suggesting integration and processing of both limbic and sensory information. Organizational principles learned from animal models at the mesoscopic level (brain regions and pathways) will not only critically inform future work at the microscopic (single neurons and synapses) level, but also have translational value to advance our understanding of human brain architecture.

Phase-amplitude coupling and infraslow (<1 Hz) frequencies in the rat brain: relationship to resting state fMRI

Resting state functional magnetic resonance imaging (fMRI) can identify network alterations that occur in complex psychiatric diseases and behaviors, but its interpretation is difficult because the neural basis of the infraslow BOLD fluctuations is poorly understood. Previous results link dynamic activity during the resting state to both infraslow frequencies in local field potentials (LFP) (<1 Hz) and band-limited power in higher frequency LFP (>1 Hz). To investigate the relationship between these frequencies, LFPs were recorded from rats under two anesthetics: isoflurane and dexmedetomidine. Signal phases were calculated from low-frequency LFP and compared to signal amplitudes from high-frequency LFP to determine if modulation existed between the two frequency bands (phase-amplitude coupling). Isoflurane showed significant, consistent phase-amplitude coupling at nearly all pairs of frequencies, likely due to the burst-suppression pattern of activity that it induces. However, no consistent phase-amplitude coupling was observed in rats that were anesthetized with dexmedetomidine. fMRI-LFP correlations under isoflurane using high frequency LFP were reduced when the low frequency LFP's influence was accounted for, but not vice-versa, or in any condition under dexmedetomidine. The lack of consistent phase-amplitude coupling under dexmedetomidine and lack of shared variance between high frequency and low frequency LFP as it relates to fMRI suggests that high and low frequency neural electrical signals may contribute differently, possibly even independently, to resting state fMRI. This finding suggests that researchers take care in interpreting the neural basis of resting state fMRI, as multiple dynamic factors in the underlying electrophysiology could be driving any particular observation.

Dynamic repertoire of intrinsic brain states is reduced in propofol-induced unconsciousness

The richness of conscious experience is thought to scale with the size of the repertoire of causal brain states, and it may be diminished in anesthesia. We estimated the state repertoire from dynamic analysis of intrinsic functional brain networks in conscious sedated and unconscious anesthetized rats. Functional resonance images were obtained from 30-min whole-brain resting-state blood oxygen level-dependent (BOLD) signals at propofol infusion rates of 20 and 40 mg/kg/h, intravenously. Dynamic brain networks were defined at the voxel level by sliding window analysis of regional homogeneity (ReHo) or coincident threshold crossings (CTC) of the BOLD signal acquired in nine sagittal slices. The state repertoire was characterized by the temporal variance of the number of voxels with significant ReHo or positive CTC. From low to high propofol dose, the temporal variances of ReHo and CTC were reduced by 78% ± 20% and 76%± 20%, respectively. Both baseline and propofol-induced reduction of CTC temporal variance increased from lateral to medial position. Group analysis showed a 20% reduction in the number of unique states at the higher propofol dose. Analysis of temporal variance in 12 anatomically defined regions of interest predicted that the largest changes occurred in visual cortex, parietal cortex, and caudate-putamen. The results suggest that the repertoire of large-scale brain states derived from the spatiotemporal dynamics of intrinsic networks is substantially reduced at an anesthetic dose associated with loss of consciousness.

Time-dependent effects of isoflurane and dexmedetomidine on functional connectivity, spectral characteristics, and spatial distribution of spontaneous BOLD fluctuations

Anesthesia is often necessary to perform fMRI experiments in the rodent model; however, commonly used anesthetic protocols may manifest changing brain conditions over the duration of the study. This possibility was explored in the current work. Eleven rats were anesthetized with 2% isoflurane anesthesia; four rats were anesthetized for a short period (30 min, simulating induction and fMRI setup) and seven rats were anesthetized for a long period (3 h, simulating surgical preparation). Following the initial anesthetic period, isoflurane was discontinued, and a dexmedetomidine bolus (0.025 mg/kg) and continuous subcutaneous infusion (0.05 mg/kg/h) were administered. Blood-oxygen-level dependent resting state imaging was performed every 30 min from 0.75 h post dexmedetomidine bolus until 5.75 h post-bolus. Evaluation of power spectra obtained from time courses in the primary somatosensory cortex revealed, in general, a monotonic increase in low-frequency power (0.05-0.3 Hz) in both groups over the duration of resting state imaging. Greater low-band spectral power (0.05-0.15 Hz) is present in the short isoflurane group for the first 2.75 h, but the spectra become highly uniform at 3.25 h. The emergence of a ~0.18 Hz peak, beginning at the 3.75 h time point, exists in both groups and evolves similarly, increasing in strength as the duration of dexmedetomidine sedation (and time since isoflurane cessation) extends. In the long isoflurane group only, bilateral functional connectivity strengthens with anesthetic duration, and correlation is linearly linked to low-band spectral power. Convergence of connectivity and spectral metrics between the short and long isoflurane groups occurs at ~3.25 h, suggesting the effects of isoflurane have subsided. Researchers using dexmedetomidine following isoflurane for functional studies should be aware of the duration specific effects of the pre-scan isoflurane durations as well as the continuing influences of long-term imaging under dexmedetomidine.

Remote brain network changes after unilateral cortical impact injury and their modulation by acetylcholinesterase inhibition

We explored whether cerebral cortical impact injury (CCI) effects extend beyond direct lesion sites to affect remote brain networks, and whether acetylcholinesterase (AChE) inhibition elicits discrete changes in functional activation of motor circuits following CCI. Adult male rats underwent unilateral motor-sensory CCI or sham injury. Physostigmine (AChE inhibitor) or saline were administered subcutaneously continuously via implanted minipumps (1.6 micromoles/kg/day) for 3 weeks, followed by cerebral perfusion mapping during treadmill walking using [(14)C]-iodoantipyrine. Quantitative autoradiographs were analyzed by statistical parametric mapping and functional connectivity (FC) analysis. CCI resulted in functional deficits in the ipsilesional basal ganglia, with increased activation contralesionally. Recruitment was also observed, especially contralesionally, of the red nucleus, superior colliculus, pedunculopontine tegmental nucleus, thalamus (ventrolateral n., central medial n.), cerebellum, and sensory cortex. FC decreased significantly within ipsi- and contralesional motor circuits and between hemispheres, but increased between midline cerebellum and select regions of the basal ganglia within each hemisphere. Physostigmine significantly increased functional brain activation in the cerebellar thalamocortical pathway (midline cerebellum→ventrolateral thalamus→motor cortex), subthalamic nucleus/zona incerta, and red nucleus and bilateral sensory cortex. In conclusion, CCI resulted in increased functional recruitment of contralesional motor cortex and bilateral subcortical motor regions, as well as recruitment of the cerebellar-thalamocortical circuit and contralesional sensory cortex. This phenomenon, augmented by physostigmine, may partially compensate motor deficits. FC decreased inter-hemispherically and in negative, but not positive, intra-hemispherical FC, and it was not affected by physostigmine. Circuit-based approaches into functional brain reorganization may inform future behavioral or molecular strategies to augment targeted neurorehabilitation.

Effects of severing the corpus callosum on electrical and BOLD functional connectivity and spontaneous dynamic activity in the rat brain

Functional networks, defined by synchronous spontaneous blood oxygenation level-dependent (BOLD) oscillations between spatially distinct brain regions, appear to be essential to brain function and have been implicated in disease states, cognitive capacity, and sensing and motor processes. While the topographical extent and behavioral function of these networks has been extensively investigated, the neural functions that create and maintain these synchronizations remain mysterious. In this work callosotomized rodents are examined, providing a unique platform for evaluating the influence of structural connectivity via the corpus callosum on bilateral resting state functional connectivity. Two experimental groups were assessed, a full callosotomy group, in which the corpus callosum was completely sectioned, and a sham callosotomy group, in which the gray matter was sectioned but the corpus callosum remained intact. Results indicated a significant reduction in interhemispheric connectivity in the full callosotomy group as compared with the sham group in primary somatosensory cortex and caudate-putamen regions. Similarly, electrophysiology revealed significantly reduced bilateral correlation in band limited power. Bilateral gamma Band-limited power connectivity was most strongly affected by the full callosotomy procedure. This work represents a robust finding indicating the corpus callosum's influence on maintaining integrity in bilateral functional networks; further, functional magnetic resonance imaging (fMRI) and electrophysiological connectivity share a similar decrease in connectivity as a result of the callosotomy, suggesting that fMRI-measured functional connectivity reflects underlying changes in large-scale coordinated electrical activity. Finally, spatiotemporal dynamic patterns were evaluated in both groups; the full callosotomy rodents displayed a striking loss of bilaterally synchronous propagating waves of cortical activity.

Neuroimaging evidence of major morpho-anatomical and functional abnormalities in the BTBR T+TF/J mouse model of autism

BTBR T+tf/J (BTBR) mice display prominent behavioural deficits analogous to the defining symptoms of autism, a feature that has prompted a widespread use of the model in preclinical autism research. Because neuro-behavioural traits are described with respect to reference populations, multiple investigators have examined and described the behaviour of BTBR mice against that exhibited by C57BL/6J (B6), a mouse line characterised by high sociability and low self-grooming. In an attempt to probe the translational relevance of this comparison for autism research, we used Magnetic Resonance Imaging (MRI) to map in both strain multiple morpho-anatomical and functional neuroimaging readouts that have been extensively used in patient populations. Diffusion tensor tractography confirmed previous reports of callosal agenesis and lack of hippocampal commissure in BTBR mice, and revealed a concomitant rostro-caudal reorganisation of major cortical white matter bundles. Intact inter-hemispheric tracts were found in the anterior commissure, ventro-medial thalamus, and in a strain-specific white matter formation located above the third ventricle. BTBR also exhibited decreased fronto-cortical, occipital and thalamic gray matter volume and widespread reductions in cortical thickness with respect to control B6 mice. Foci of increased gray matter volume and thickness were observed in the medial prefrontal and insular cortex. Mapping of resting-state brain activity using cerebral blood volume weighted fMRI revealed reduced cortico-thalamic function together with foci of increased activity in the hypothalamus and dorsal hippocampus of BTBR mice. Collectively, our results show pronounced functional and structural abnormalities in the brain of BTBR mice with respect to control B6 mice. The large and widespread white and gray matter abnormalities observed do not appear to be representative of the neuroanatomical alterations typically observed in autistic patients. The presence of reduced fronto-cortical metabolism is of potential translational relevance, as this feature recapitulates previously-reported clinical observations.

Dynamic functional connectivity: promise, issues, and interpretations

The brain must dynamically integrate, coordinate, and respond to internal and external stimuli across multiple time scales. Non-invasive measurements of brain activity with fMRI have greatly advanced our understanding of the large-scale functional organization supporting these fundamental features of brain function. Conclusions from previous resting-state fMRI investigations were based upon static descriptions of functional connectivity (FC), and only recently studies have begun to capitalize on the wealth of information contained within the temporal features of spontaneous BOLD FC. Emerging evidence suggests that dynamic FC metrics may index changes in macroscopic neural activity patterns underlying critical aspects of cognition and behavior, though limitations with regard to analysis and interpretation remain. Here, we review recent findings, methodological considerations, neural and behavioral correlates, and future directions in the emerging field of dynamic FC investigations.

Quasi-periodic patterns (QPP): large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity

Functional connectivity measurements from resting state blood-oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) are proving a powerful tool to probe both normal brain function and neuropsychiatric disorders. However, the neural mechanisms that coordinate these large networks are poorly understood, particularly in the context of the growing interest in network dynamics. Recent work in anesthetized rats has shown that the spontaneous BOLD fluctuations are tightly linked to infraslow local field potentials (LFPs) that are seldom recorded but comparable in frequency to the slow BOLD fluctuations. These findings support the hypothesis that long-range coordination involves low frequency neural oscillations and establishes infraslow LFPs as an excellent candidate for probing the neural underpinnings of the BOLD spatiotemporal patterns observed in both rats and humans. To further examine the link between large-scale network dynamics and infraslow LFPs, simultaneous fMRI and microelectrode recording were performed in anesthetized rats. Using an optimized filter to isolate shared components of the signals, we found that time-lagged correlation between infraslow LFPs and BOLD is comparable in spatial extent and timing to a quasi-periodic pattern (QPP) found from BOLD alone, suggesting that fMRI-measured QPPs and the infraslow LFPs share a common mechanism. As fMRI allows spatial resolution and whole brain coverage not available with electroencephalography, QPPs can be used to better understand the role of infraslow oscillations in normal brain function and neurological or psychiatric disorders.

Metabolic brain covariant networks as revealed by FDG-PET with reference to resting-state fMRI networks

The human brain is inherently organized as separate networks, as has been widely revealed by resting-state functional magnetic resonance imaging (fMRI). Although the large-scale functional connectivity can be partially explained by the underlying white-matter structural connectivity, the question of whether the underlying functional connectivity is related to brain metabolic factors is still largely unanswered. The present study investigated the presence of metabolic covariant networks across subjects using a set of fluorodeoxyglucose ((18)F, FDG) positron-emission tomography (PET) images. Spatial-independent component analysis was performed on the subject series of FDG-PET images. A number of networks that were mainly homotopic regions could be identified, including visual, auditory, motor, cerebellar, and subcortical networks. However, the anterior-posterior networks such as the default-mode and left frontoparietal networks could not be observed. Region-of-interest-based correlation analysis confirmed that the intersubject metabolic covariances within the default-mode and left frontoparietal networks were reduced as compared with corresponding time-series correlations using resting-state fMRI from an independent sample. In contrast, homotopic intersubject metabolic covariances observed using PET were comparable to the corresponding fMRI resting-state time-series correlations. The current study provides preliminary illustration, suggesting that the human brain metabolism pertains to organized covariance patterns that might partially reflect functional connectivity as revealed by resting-state blood oxygen level dependent (BOLD). The discrepancy between the PET covariance and BOLD functional connectivity might reflect the differences of energy consumption coupling and ongoing neural synchronization within these brain networks.

Neural correlates of time-varying functional connectivity in the rat

Functional connectivity between brain regions, measured with resting state functional magnetic resonance imaging, holds great potential for understanding the basis of behavior and neuropsychiatric diseases. Recently it has become clear that correlations between the blood oxygenation level dependent (BOLD) signals from different areas vary over the course of a typical scan (6-10 min in length), though the changes are obscured by standard methods of analysis that assume the relationships are stationary. Unfortunately, because similar variability is observed in signals that share no temporal information, it is unclear which dynamic changes are related to underlying neural events. To examine this question, BOLD data were recorded simultaneously with local field potentials (LFP) from interhemispheric primary somatosensory cortex (SI) in anesthetized rats. LFP signals were converted into band-limited power (BLP) signals including delta, theta, alpha, beta and gamma. Correlation between signals from interhemispheric SI was performed in sliding windows to produce signals of correlation over time for BOLD and each BLP band. Both BOLD and BLP signals showed large changes in correlation over time and the changes in BOLD were significantly correlated to the changes in BLP. The strongest relationship was seen when using the theta, beta and gamma bands. Interestingly, while steady-state BOLD and BLP correlate with the global fMRI signal, dynamic BOLD becomes more like dynamic BLP after the global signal is regressed. As BOLD sliding window connectivity is partially reflecting underlying LFP changes, the present study suggests it may be a valuable method of studying dynamic changes in brain states.

Infraslow LFP correlates to resting-state fMRI BOLD signals

The slow fluctuations of the blood-oxygenation-level dependent (BOLD) signal in resting-state fMRI are widely utilized as a surrogate marker of ongoing neural activity. Spontaneous neural activity includes a broad range of frequencies, from infraslow (<0.5 Hz) fluctuations to fast action potentials. Recent studies have demonstrated a correlative relationship between the BOLD fluctuations and power modulations of the local field potential (LFP), particularly in the gamma band. However, the relationship between the BOLD signal and the infraslow components of the LFP, which are directly comparable in frequency to the BOLD fluctuations, has not been directly investigated. Here we report a first examination of the temporal relation between the resting-state BOLD signal and infraslow LFPs using simultaneous fMRI and full-band LFP recording in rat. The spontaneous BOLD signal at the recording sites exhibited significant localized correlation with the infraslow LFP signals as well as with the slow power modulations of higher-frequency LFPs (1-100 Hz) at a delay comparable to the hemodynamic response time under anesthesia. Infraslow electrical activity has been postulated to play a role in attentional processes, and the findings reported here suggest that infraslow LFP coordination may share a mechanism with the large-scale BOLD-based networks previously implicated in task performance, providing new insight into the mechanisms contributing to the resting state fMRI signal.

Dynamic properties of functional connectivity in the rodent

Functional connectivity mapping with resting-state magnetic resonance imaging (MRI) has become an immensely powerful technique that provides insight into both normal cognitive function and disruptions linked to neurological disorders. Traditionally, connectivity is mapped using data from an entire scan (minutes), but it is well known that cognitive processes occur on much shorter time scales (seconds). Recent studies have demonstrated that the correlation between the blood oxygenation level-dependent (BOLD) MRI signal from different areas varies over time, motivating a further exploration of these fluctuations in apparent connectivity. However, it has also been shown that similar changes in correlation can arise when the timing relationships between voxels are randomized (Handwerker et al., 2012 ). In this work, we show that functional connectivity in the anesthetized rat exhibits dynamic properties that are similar to those previously observed in awake humans (Chang and Glover, 2010 ) and anesthetized monkeys (Hutchison et al., 2012 ). Sliding window correlation between BOLD time courses obtained from bilateral cortical and subcortical regions of interest results in periods of variable positive and negative correlation for most pairs of areas except homologous areas in opposite hemispheres, which exhibit a primarily positive correlation. A comparison with sliding window correlation of randomly matched time courses suggests that with the exception of homologous areas and sensorimotor connections, the dynamics cannot be distinguished from random fluctuations in correlation over time, supporting the idea that some of these dynamic patterns may be due to inherent properties of the signal rather than variations in neural coherence. Within the pairs of areas where the dynamics are most different from those of randomly matched time courses, ten common patterns of connectivity are identified, and their occurrence as a function of time is plotted for all animals. The observation of time-varying correlation in the rodent model will facilitate the future multimodal experiments needed to determine whether the changes in apparent connectivity are linked to underlying neural variability.

Track-weighted functional connectivity (TW-FC): a tool for characterizing the structural-functional connections in the brain

MRI provides a powerful tool for studying the functional and structural connections in the brain non-invasively. The technique of functional connectivity (FC) exploits the intrinsic temporal correlations of slow spontaneous signal fluctuations to characterise brain functional networks. In addition, diffusion MRI fibre-tracking can be used to study the white matter structural connections. In recent years, there has been considerable interest in combining these two techniques to provide an overall structural-functional description of the brain. In this work we applied the recently proposed super-resolution track-weighted imaging (TWI) methodology to demonstrate how whole-brain fibre-tracking data can be combined with FC data to generate a track-weighted (TW) FC map of FC networks. The method was applied to data from 8 healthy volunteers, and illustrated with (i) FC networks obtained using a seeded connectivity-based analysis (seeding in the precuneus/posterior cingulate cortex, PCC, known to be part of the default mode network), and (ii) with FC networks generated using independent component analysis (in particular, the default mode, attention, visual, and sensory-motor networks). TW-FC maps showed high intensity in white matter structures connecting the nodes of the FC networks. For example, the cingulum bundles show the strongest TW-FC values in the PCC seeded-based analysis, due to their major role in the connection between medial frontal cortex and precuneus/posterior cingulate cortex; similarly the superior longitudinal fasciculus was well represented in the attention network, the optic radiations in the visual network, and the corticospinal tract and corpus callosum in the sensory-motor network. The TW-FC maps highlight the white matter connections associated with a given FC network, and their intensity in a given voxel reflects the functional connectivity of the part of the nodes of the network linked by the structural connections traversing that voxel. They therefore contain a different (and novel) image contrast from that of the images used to generate them. The results shown in this study illustrate the potential of the TW-FC approach for the fusion of structural and functional data into a single quantitative image. This technique could therefore have important applications in neuroscience and neurology, such as for voxel-based comparison studies.

Thalamocortical dysconnectivity in schizophrenia

OBJECTIVE

The thalamus and cerebral cortex are connected via topographically organized, reciprocal connections. Previous studies have revealed thalamic abnormalities in schizophrenia; however, it is not known whether thalamocortical networks are differentially affected in the disorder. To explore this possibility, the authors examined functional connectivity in intrinsic low-frequency blood-oxygen-level-dependent (BOLD) signal fluctuations between major divisions of the cortex and thalamus using resting-state functional MRI (fMRI).

METHOD

Seventy-seven healthy subjects and 62 patients with schizophrenia underwent resting-state fMRI. To identify functional subdivisions of the thalamus, the authors parceled the cortex into six regions of interest: the prefrontal cortex, motor cortex/supplementary motor area, somatosensory cortex, temporal lobe, posterior parietal cortex, and occipital lobe. Mean BOLD time series were extracted for each region of interest and entered into a seed-based functional connectivity analysis.

RESULTS

Consistent with previous reports, activity in distinct cortical areas correlated with specific, largely nonoverlapping regions of the thalamus in both healthy comparison subjects and schizophrenia patients. Direct comparison between groups revealed reduced prefrontal-thalamic connectivity and increased motor/somatosensory-thalamic connectivity in schizophrenia. The changes in connectivity were unrelated to local gray matter content within the thalamus and to antipsychotic medication dosage. No differences were observed in temporal, posterior parietal, or occipital cortex connectivity with the thalamus.

CONCLUSIONS

These findings establish differential abnormalities of thalamocortical networks in schizophrenia. The etiology of schizophrenia may disrupt the development of prefrontal-thalamic connectivity and refinement of somatomotor connectivity with the thalamus that occurs during brain maturation.

Monkey in the middle: why non-human primates are needed to bridge the gap in resting-state investigations

Resting-state investigations based on the evaluation of intrinsic low-frequency fluctuations of the BOLD fMRI signal have been extensively utilized to map the structure and dynamics of large-scale functional network organization in humans. In addition to increasing our knowledge of normal brain connectivity, disruptions of the spontaneous hemodynamic fluctuations have been suggested as possible diagnostic indicators of neurological and psychiatric disease states. Though the non-invasive technique has been received with much acclamation, open questions remain regarding the origin, organization, phylogenesis, as well as the basis of disease-related alterations underlying the signal patterns. Experimental work utilizing animal models, including the use of neurophysiological recordings and pharmacological manipulations, therefore, represents a critical component in the understanding and successful application of resting-state analysis, as it affords a range of experimental manipulations not possible in human subjects. In this article, we review recent rodent and non-human primate studies and based on the examination of the homologous brain architecture propose the latter to be the best-suited model for exploring these unresolved resting-state concerns. Ongoing work examining the correspondence of functional and structural connectivity, state-dependency and the neuronal correlates of the hemodynamic oscillations are discussed. We then consider the potential experiments that will allow insight into different brain states and disease-related network disruptions that can extend the clinical applications of resting-state fMRI (RS-fMRI).

Structure-function relationship of cerebral networks in experimental neuroscience: contribution of magnetic resonance imaging

The analysis of neuronal networks, their interactions in resting condition as well as during brain activation have become of great interest for a better understanding of the signal processing of the brain during sensory stimulus or cognitive tasks. Parallel to the study of the functional networks and their dynamics, the underlying network structure is highly important as it provides the basis of the functional interaction. Moreover, under pathological conditions, some nodes in such a net may be impaired and the function of the whole network affected. Mechanisms such as functional deficit and improvement, and plastic reorganization are increasingly discussed in the context of existing structural and functional networks. While many of these aspects have been followed in human and clinical studies, the experimental range is limited for obvious reasons. Here, animal experimental studies are needed as they permit longer scan times and, moreover, comparison with invasive histology. Experimental non-invasive imaging modalities are now able to perform impressive contributions. In this review we try to highlight most recent new cutting-edge developments and applications in experimental neuroscience of functional and structural networks of the brain, relying on non-invasive imaging. We focus primarily on the potential of experimental Magnetic Resonance Imaging (MRI), but also touch upon micro positron emission tomography (μPET) and optical imaging developments where they are applicable to the topic of the present review.

Anxiolytic-like effect of pregabalin on unconditioned fear in the rat: an autoradiographic brain perfusion mapping and functional connectivity study

Clinical and preclinical evidence suggests anxiolytic-like efficacy of pregabalin (PGB, Lyrica). However, its mechanism of action remains under investigation. The current study applied [(14)C]-iodoantipyrine cerebral blood flow (CBF) mapping to examine the effect of PGB on neural substrates underlying unconditioned fear in a rat model of footshock-induced fear. Regional CBF (rCBF) was analyzed by statistical parametric mapping. Functional connectivity and graph theoretical analysis were used to investigate how footshock and PGB affect brain activation at the network level. Pregabalin significantly attenuated footshock-induced ultrasonic vocalization, but showed no significant effect on freezing behavior. Footshock compared to no-shock controls elicited significant increases in rCBF in limbic/paralimbic regions implicated in the processing of unconditioned fear and ultrasonic vocalization, including the amygdala, hypothalamus, lateral septum, dorsal periaqueductal gray, the anterior insular (aINS) and medial prefrontal cortex (mPFC). The activation pattern was similar in vehicle- and PGB-treated subjects, with PGB significantly attenuating activation in the amygdala, hypothalamus, and aINS. The vehicle/no-shock group showed strong, positive intra-structural correlations within the cortex, hypothalamus, amygdala, thalamus, and brainstem. The cortex was negatively correlated with the hypothalamus and brainstem. Footshock reduced the total number of significant correlations, but induced greater intra-cortical connectivity of the aINS and mPFC, and new positive correlations between the hypothalamus and amygdala. In no-shock controls, PGB significantly reduced the positive intra-structural correlations within the cortex and amygdala, as well as the negative cortico-subcortical correlations. Following footshocks, PGB disrupted both the network recruitment of aINS and mPFC, and the positive hypothalamic-amygdaloid correlations. Our findings suggest that PGB may exert anxiolytic effect by attenuating cortico-cortical and cortico-subcortical communication and inhibiting network recruitment of the aINS, mPFC, amygdala, and hypothalamus following a fear-inducing stimulus. Functional brain mapping in rodents may provide new endpoints for preclinical evaluation of anxiolytic drug candidates with potentially improved translational power compared to behavioral measurements alone.

An in vivo MRI Template Set for Morphometry, Tissue Segmentation, and fMRI Localization in Rats

Over the last decade, several papers have focused on the construction of highly detailed mouse high field magnetic resonance image (MRI) templates via non-linear registration to unbiased reference spaces, allowing for a variety of neuroimaging applications such as robust morphometric analyses. However, work in rats has only provided medium field MRI averages based on linear registration to biased spaces with the sole purpose of approximate functional MRI (fMRI) localization. This precludes any morphometric analysis in spite of the need of exploring in detail the neuroanatomical substrates of diseases in a recent advent of rat models. In this paper we present a new in vivo rat T2 MRI template set, comprising average images of both intensity and shape, obtained via non-linear registration. Also, unlike previous rat template sets, we include white and gray matter probabilistic segmentations, expanding its use to those applications demanding prior-based tissue segmentation, e.g., statistical parametric mapping (SPM) voxel-based morphometry. We also provide a preliminary digitalization of latest Paxinos and Watson atlas for anatomical and functional interpretations within the cerebral cortex. We confirmed that, like with previous templates, forepaw and hindpaw fMRI activations can be correctly localized in the expected atlas structure. To exemplify the use of our new MRI template set, were reported the volumes of brain tissues and cortical structures and probed their relationships with ontogenetic development. Other in vivo applications in the near future can be tensor-, deformation-, or voxel-based morphometry, morphological connectivity, and diffusion tensor-based anatomical connectivity. Our template set, freely available through the SPM extension website, could be an important tool for future longitudinal and/or functional extensive preclinical studies.

Translating glutamate: from pathophysiology to treatment

The neurotransmitter glutamate is the primary excitatory neurotransmitter in mammalian brain and is responsible for most corticocortical and corticofugal neurotransmission. Disturbances in glutamatergic function have been implicated in the pathophysiology of several neuropsychiatric disorders-including schizophrenia, drug abuse and addiction, autism, and depression-that were until recently poorly understood. Nevertheless, improvements in basic information regarding these disorders have yet to translate into Food and Drug Administration-approved treatments. Barriers to translation include the need not only for improved compounds but also for improved biomarkers sensitive to both structural and functional target engagement and for improved translational models. Overcoming these barriers will require unique collaborative arrangements between pharma, government, and academia. Here, we review a recent Institute of Medicine-sponsored meeting, highlighting advances in glutamatergic theories of neuropsychiatric illness as well as remaining barriers to treatment development.

Functional connectivity fMRI of the rodent brain: comparison of functional connectivity networks in rat and mouse

At present, resting state functional MRI (rsfMRI) is increasingly used in human neuropathological research. The present study aims at implementing rsfMRI in mice, a species that holds the widest variety of neurological disease models. Moreover, by acquiring rsfMRI data with a comparable protocol for anesthesia, scanning and analysis, in both rats and mice we were able to compare findings obtained in both species. The outcome of rsfMRI is different for rats and mice and depends strongly on the applied number of components in the Independent Component Analysis (ICA). The most important difference was the appearance of unilateral cortical components for the mouse resting state data compared to bilateral rat cortical networks. Furthermore, a higher number of components was needed for the ICA analysis to separate different cortical regions in mice as compared to rats.