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

The field of brain connectomics develops our understanding of the brain's intrinsic organization by characterizing trends in spontaneous brain activity. Linear correlations in spontaneous blood-oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) fluctuations are often used as measures of functional connectivity (FC), that is, as a quantity describing how similarly two brain regions behave over time. Given the natural spectral scaling of BOLD-fMRI signals, it may be useful to represent BOLD-fMRI as multiple processes occurring over multiple scales. The wavelet domain presents a transform space well suited to the examination of multiscale systems as the wavelet basis set is constructed from a self-similar rescaling of a time and frequency delimited kernel. In the present study, we utilize wavelet transforms to examine fluctuations in whole-brain BOLD-fMRI connectivity as a function of wavelet spectral scale in a sample (N = 31) of resting healthy human volunteers. Information theoretic criteria measure relatedness between spectrally-delimited FC graphs. Voxelwise comparisons of between-spectra graph structures illustrate the development of preferential functional networks across spectral bands.

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.

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.

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.

Measures of whole-brain activity, from techniques such as functional Magnetic Resonance Imaging, provide a means to observe the brain's dynamical operations. However, interpretation of whole-brain dynamics has been stymied by the inherently high-dimensional structure of brain activity. The present research addresses this challenge through a series of scale transformations in the spectral, spatial, and relational domains. Instantaneous multispectral dynamics are first developed from input data via a wavelet filter bank. Voxel-level signals are then projected onto a representative set of spatially independent components. The correlation distance over the instantaneous wavelet-ICA state vectors is a graph that may be embedded onto a lower-dimensional space to assist the interpretation of state-space dynamics. Applying this procedure to a large sample of resting-state and task-active data (acquired through the Human Connectome Project), we segment the empirical state space into a continuum of stimulus-dependent brain states. Upon observing the local neighborhood of brain-states adopted subsequent to each stimulus, we may conclude that resting brain activity includes brain states that are, at times, similar to those adopted during tasks, but that are at other times distinct from task-active brain states. As task-active brain states often populate a local neighborhood, back-projection of segments of the dynamical state space onto the brain's surface reveals the patterns of brain activity that support many experimentally-defined states.

While functional connectivity has typically been calculated over the entire length of the scan (5-10. min), 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.

To examine the neural basis of the blood oxygenation level dependent (BOLD) magnetic resonance imaging (MRI) signal, we have developed a rodent model in which functional MRI data and in vivo intracortical recording can be performed simultaneously. The combination of MRI and electrical recording is technically challenging because the electrodes used for recording distort the MRI images and the MRI acquisition induces noise in the electrical recording. To minimize the mutual interference of the two modalities, glass microelectrodes were used rather than metal and a noise removal algorithm was implemented for the electrophysiology data. In our studies, two microelectrodes were separately implanted in bilateral primary somatosensory cortices (SI) of the rat and fixed in place. One coronal slice covering the electrode tips was selected for functional MRI. Electrode shafts and fixation positions were not included in the image slice to avoid imaging artifacts. The removed scalp was replaced with toothpaste to reduce susceptibility mismatch and prevent Gibbs ringing artifacts in the images. The artifact structure induced in the electrical recordings by the rapidly-switching magnetic fields during image acquisition was characterized by averaging all cycles of scans for each run. The noise structure during imaging was then subtracted from original recordings. The denoised time courses were then used for further analysis in combination with the fMRI data. As an example, the simultaneous acquisition was used to determine the relationship between spontaneous fMRI BOLD signals and band-limited intracortical electrical activity. Simultaneous fMRI and electrophysiological recording in the rodent will provide a platform for many exciting applications in neuroscience in addition to elucidating the relationship between the fMRI BOLD signal and neuronal activity.

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.