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.

Resting state functional magnetic resonance (rs-fMRI) imaging offers insights into how different brain regions are connected into functional networks. It was recently shown that networks that are almost identical to the ones created from conventional correlation analysis can be obtained from a subset of high-amplitude data, suggesting that the functional networks may be driven by instantaneous co-activations of multiple brain regions rather than ongoing oscillatory processes. The rs-fMRI studies, however, rely on the blood oxygen level dependent (BOLD) signal, which is only indirectly sensitive to neural activity through neurovascular coupling. To provide more direct evidence that the neuronal co-activation events produce the time-varying network patterns seen in rs-fMRI studies, we examined the simultaneous rs-fMRI and local field potential (LFP) recordings in rats performed in our lab over the past several years. We developed complementary analysis methods that focus on either the temporal or spatial domain, and found evidence that the interaction between LFP and BOLD may be driven by instantaneous co-activation events as well. BOLD maps triggered on high-amplitude LFP events resemble co-activation patterns created from rs-fMRI data alone, though the co-activation time points are defined differently in the two cases. Moreover, only LFP events that fall into the highest or lowest thirds of the amplitude distribution result in a BOLD signal that can be distinguished from noise. These findings provide evidence of an electrophysiological basis for the time-varying co-activation patterns observed in previous studies.

A better understanding of how behavioral performance emerges from interacting brain systems may come from analysis of functional networks using functional magnetic resonance imaging. Recent studies comparing such networks with human behavior have begun to identify these relationships, but few have used a time scale small enough to relate their findings to variation within a single individual's behavior. In the present experiment we examined the relationship between a psychomotor vigilance task and the interacting default mode and task positive networks. Two time-localized comparative metrics were calculated: difference between the two networks' signals at various time points around each instance of the stimulus (peristimulus times) and correlation within a 12.3-s window centered at each peristimulus time. Correlation between networks was also calculated within entire resting-state functional imaging runs from the same individuals. These metrics were compared with response speed on both an intraindividual and an interindividual basis. In most cases, a greater difference or more anticorrelation between networks was significantly related to faster performance. While interindividual analysis showed this result generally, using intraindividual analysis it was isolated to peristimulus times 4 to 8 s before the detected target. Within that peristimulus time span, the effect was stronger for individuals who tended to have faster response times. These results suggest that the relationship between functional networks and behavior can be better understood by using shorter time windows and also by considering both intraindividual and interindividual variability.

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.

Functional magnetic resonance imaging (fMRI) is currently one of the most important neuroimaging methods in neuroscience. The image contrast in fMRI relies on the blood-oxygenation-level dependent (BOLD) signal, which indirectly reflects neural activity through neurovascular coupling. Because the mechanism that links the BOLD signal to neural activities involves multiple complicated processes, where neural activity, regional metabolism, hemodynamics, and the BOLD signal are all inter-connected, understanding the quantitative relationship between the BOLD signal and the underlying neural activities is crucial for interpreting fMRI data. Simultaneous local field potential (LFP) and fMRI recordings provide a method to study neurovascular coupling. There were a few studies that have shown non-linearities in stimulus related responses, but whether there is any non-linearity in LFP—BOLD relationship at rest has not been specifically quantified. In this study, we analyzed the simultaneous LFP and resting state-fMRI data acquired from rodents, and found that the relationship between LFP and BOLD is non-linear under isoflurane (ISO) anesthesia, but linear under dexmedetomidine (DMED) anesthesia. Subsequent analysis suggests that such non-linearity may come from the non-Gaussian distribution of LFP power and switching from LFP power to LFP amplitude can alleviate the problem to a degree. We also confirmed that, despite the non-linearity in the mean LFP—BOLD curve, the Pearson correlation between the two signals is relatively unaffected.

Resting state functional MRI (rs-fMRI) creates a rich four-dimensional data set that can be analyzed in a variety of ways. As more researchers come to view the brain as a complex dynamical system, tools are increasingly being drawn from other fields to characterize the complexity of the brain’s activity. However, given that the signal measured with rs-fMRI arises from the hemodynamic response to neural activity, the extent to which complexity metrics reflect neural complexity as compared to signal properties related to image quality remains unknown. To provide some insight into this question, correlation dimension, approximate entropy and Lyapunov exponent were calculated for different rs-fMRI scans from the same subject to examine their reliability. The metrics of complexity were then compared to several properties of the rs-fMRI signal from each brain area to determine if basic signal features could explain differences in the complexity metrics. Differences in complexity across brain areas were highly reliable and were closely linked to differences in the frequency profiles of the rs-fMRI signal. The spatial distributions of the complexity and frequency metrics suggest that they are both influenced by location-dependent signal properties that can obscure changes related to neural activity.

How do intrinsic brain dynamics interact with processing of external sensory stimuli? We sought new insights using functional magnetic resonance imaging to track spatiotemporal activity patterns at the whole brain level in lightly anesthetized mice, during both resting conditions and visual stimulation trials. Our results provide evidence that quasiperiodic patterns (QPPs) are the most prominent component of mouse resting brain dynamics. These QPPs captured the temporal alignment of anticorrelation between the default mode (DMN)- and task-positive (TPN)-like networks, with global brain fluctuations, and activity in neuromodulatory nuclei of the reticular formation. Specifically, the phase of QPPs prior to stimulation could significantly stratify subsequent visual response magnitude, suggesting QPPs relate to brain state fluctuations. This is the first observation in mice that dynamics of the DMN- and TPN-like networks, and particularly their anticorrelation, capture a brain state dynamic that affects sensory processing. Interestingly, QPPs also displayed transient onset response properties during visual stimulation, which covaried with deactivations in the reticular formation. We conclude that QPPs appear to capture a brain state fluctuation that may be orchestrated through neuromodulation. Our findings provide new frontiers to understand the neural processes that shape functional brain states and modulate sensory input processing.

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.