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

Dynamic effective connectivity

Metastability is a key source of itinerant dynamics in the brain; namely, spontaneous spatiotemporal reorganization of neuronal activity. This itinerancy has been the focus of numerous dynamic functional connectivity (DFC) analyses - developed to characterize the formation and dissolution of distributed functional patterns over time, using resting state fMRI. However, aside from technical and practical controversies, these approaches cannot recover the neuronal mechanisms that underwrite itinerant (e.g., metastable) dynamics-due to their descriptive, model-free nature. We argue that effective connectivity (EC) analyses are more apt for investigating the neuronal basis of metastability. To this end, we appeal to biologically-grounded models (i.e., dynamic causal modelling, DCM) and dynamical systems theory (i.e., heteroclinic sequential dynamics) to create a probabilistic, generative model of haemodynamic fluctuations. This model generates trajectories in the parametric space of EC modes (i.e., states of connectivity) that characterize functional brain architectures. In brief, it extends an established spectral DCM, to generate functional connectivity data features that change over time. This foundational paper tries to establish the model's face validity by simulating non-stationary fMRI time series and recovering key model parameters (i.e., transition probabilities among connectivity states and the parametric nature of these states) using variational Bayes. These data are further characterized using Bayesian model comparison (within and between subjects). Finally, we consider practical issues that attend applications and extensions of this scheme. Importantly, the scheme operates within a generic Bayesian framework - that can be adapted to study metastability and itinerant dynamics in any non-stationary time series.

Arousal Contributions to Resting-State fMRI Connectivity and Dynamics

Resting-state functional magnetic resonance imaging (rsfMRI) is being widely used for charting brain connectivity and dynamics in healthy and diseased brains. However, the resting state paradigm allows an unconstrained fluctuation of brain arousal, which may have profound effects on resting-state fMRI signals and associated connectivity/dynamic metrics. Here, we review current understandings of the relationship between resting-state fMRI and brain arousal, in particular the effect of a recently discovered event of arousal modulation on resting-state fMRI. We further discuss potential implications of arousal-related fMRI modulation with a focus on its potential role in mediating spurious correlations between resting-state connectivity/dynamics with physiology and behavior. Multiple hypotheses are formulated based on existing evidence and remain to be tested by future studies.

Propagations of spontaneous brain activity in awake rats

Slow propagations of spontaneous brain activity have been reported in multiple species. However, systematical investigation of the organization of such brain activity is still lacking. In this study, we analyzed propagations of spontaneous brain activity using a reference library of characteristic resting-state functional connectivity (RSFC) patterns in awake rodents. We found that transitions through multiple distinct RSFC patterns were reproducible not only in transition sequences but also in transition time delays. In addition, the organization of these transitions and their spatiotemporal dynamic patterns were revealed using a graphical model. We further identified prominent brain regions involved in these transitions. These results provide a comprehensive framework of brainwide propagations of spontaneous activity in awake rats. This study also offers a new tool to study the spatiotemporal dynamics of activity in the resting brain.

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

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

DC Shifts-fMRI: A Supplement to Event-Related fMRI

Event-related fMRI have been widely used in locating brain regions which respond to specific tasks. However, activities of brain regions which modulate or indirectly participate in the response to a specific task are not event-related. Event-related fMRI can't locate these regulatory regions, detrimental to the integrity of the result that event-related fMRI revealed. Direct-current EEG shifts (DC shifts) have been found linked to the inner brain activity, a fusion DC shifts-fMRI method may have the ability to reveal a more complete response of the brain. In this study, we used DC shifts-fMRI to verify that even when responding to a very simple task, (1) The response of the brain is more complicated than event-related fMRI generally revealed and (2) DC shifts-fMRI have the ability of revealing brain regions whose responses are not in event-related way. We used a classical and simple paradigm which is often used in auditory cortex tonotopic mapping. Data were recorded from 50 subjects (25 male, 25 female) who were presented with randomly presented pure tone sequences with six different frequencies (200, 400, 800, 1,600, 3,200, 6,400 Hz). Our traditional fMRI results are consistent with previous findings that the activations are concentrated on the auditory cortex. Our DC shifts-fMRI results showed that the cingulate-caudate-thalamus network which underpins sustained attention is positively activated while the dorsal attention network and the right middle frontal gyrus which underpin attention orientation are negatively activated. The regional-specific correlations between DC shifts and brain networks indicate the complexity of the response of the brain even to a simple task and that the DC shifts can effectively reflect these non-event-related inner brain activities.

The Not-So-Global Blood Oxygen Level-Dependent Signal

Global signal regression is a controversial processing step for resting-state functional magnetic resonance imaging, partly because the source of the global blood oxygen level-dependent (BOLD) signal remains unclear. On the one hand, nuisance factors such as motion can readily introduce coherent BOLD changes across the whole brain. On the other hand, the global signal has been linked to neural activity and vigilance levels, suggesting that it contains important neurophysiological information and should not be discarded. Any widespread pattern of coordinated activity is likely to contribute appreciably to the global signal. Such patterns may include large-scale quasiperiodic spatiotemporal patterns, known also to be tied to performance on vigilance tasks. This uncertainty surrounding the separability of the global BOLD signal from concurrent neurological processes motivated an examination of the global BOLD signal's spatial distribution. The results clarify that although the global signal collects information from all tissue classes, a diverse subset of the BOLD signal's independent components contribute the most to the global signal. Further, the timing of each network's contribution to the global signal is not consistent across volunteers, confirming the independence of a constituent process that comprises the global signal.

Resting State fMRI in Rodents

Resting state functional MRI (fMRI) and functional connectivity are widely applied in humans to examine the role of brain networks in normal function and dysfunction. A similar approach can be taken in rodents, either to obtain translational measures in models of brain disorders or to more carefully examine the neurophysiological underpinnings of the networks. A protocol for resting state functional connectivity in the anesthetized rat, from animal setup to data acquisition to possible pipelines for data analysis, is described. © 2018 by John Wiley & Sons, Inc.

Sampling Rate Effects on Resting State fMRI Metrics

Low image sampling rates used in resting state functional magnetic resonance imaging (rs-fMRI) may cause aliasing of the cardiorespiratory pulsations over the very low frequency (VLF) BOLD signal fluctuations which reflects to functional connectivity (FC). In this study, we examine the effect of sampling rate on currently used rs-fMRI FC metrics. Ultra-fast fMRI magnetic resonance encephalography (MREG) data, sampled with TR 0.1 s, was downsampled to different subsampled repetition times (sTR, range 0.3–3 s) for comparisons. Echo planar k-space sampling (TR 2.15 s) and interleaved slice collection schemes were also compared against the 3D single shot trajectory at 2.2 s sTR. The quantified connectivity metrics included stationary spatial, time, and frequency domains, as well as dynamic analyses. Time domain methods included analyses of seed-based functional connectivity, regional homogeneity (ReHo), coefficient of variation, and spatial domain group level probabilistic independent component analysis (ICA). In frequency domain analyses, we examined fractional and amplitude of low frequency fluctuations. Aliasing effects were spatially and spectrally analyzed by comparing VLF (0.01–0.1 Hz), respiratory (0.12–0.35 Hz) and cardiac power (0.9–1.3 Hz) FFT maps at different sTRs. Quasi-periodic pattern (QPP) of VLF events were analyzed for effects on dynamic FC methods. The results in conventional time and spatial domain analyses remained virtually unchanged by the different sampling rates. In frequency domain, the aliasing occurred mainly in higher sTR (1–2 s) where cardiac power aliases over respiratory power. The VLF power maps suffered minimally from increasing sTRs. Interleaved data reconstruction induced lower ReHo compared to 3D sampling (p < 0.001). Gradient recalled echo-planar imaging (EPI BOLD) data produced both better and worse metrics. In QPP analyses, the repeatability of the VLF pulse detection becomes linearly reduced with increasing sTR. In conclusion, the conventional resting state metrics (e.g., FC, ICA) were not markedly affected by different TRs (0.1–3 s). However, cardiorespiratory signals showed strongest aliasing in central brain regions in sTR 1–2 s. Pulsatile QPP and other dynamic analyses benefit linearly from short TR scanning.

Quasi-periodic patterns contribute to functional connectivity in the brain

Functional connectivity is widely used to study the coordination of activity between brain regions over time. Functional connectivity in the default mode and task positive networks is particularly important for normal brain function. However, the processes that give rise to functional connectivity in the brain are not fully understood. It has been postulated that low-frequency neural activity plays a key role in establishing the functional architecture of the brain. Quasi-periodic patterns (QPPs) are a reliably observable form of low-frequency neural activity that involve the default mode and task positive networks. Here, QPPs from resting-state and working memory task-performing individuals were acquired. The spatiotemporal pattern, strength, and frequency of the QPPs between the two groups were compared and the contribution of QPPs to functional connectivity in the brain was measured. In task-performing individuals, the spatiotemporal pattern of the QPP changes, particularly in task-relevant regions, and the QPP tends to occur with greater strength and frequency. Differences in the QPPs between the two groups could partially account for the variance in functional connectivity between resting-state and task-performing individuals. The QPPs contribute strongly to connectivity in the default mode and task positive networks and to the strength of anti-correlation seen between the two networks. Many of the connections affected by QPPs are also disrupted during several neurological disorders. These findings contribute to understanding the dynamic neural processes that give rise to functional connectivity in the brain and how they may be disrupted during disease.

Dynamic properties of simulated brain network models and empirical resting-state data

Brain network models (BNMs) have become a promising theoretical framework for simulating signals that are representative of whole-brain activity such as resting-state fMRI. However, it has been difficult to compare the complex brain activity obtained from simulations to empirical data. Previous studies have used simple metrics to characterize coordination between regions such as functional connectivity. We extend this by applying various different dynamic analysis tools that are currently used to understand empirical resting-state fMRI (rs-fMRI) to the simulated data. We show that certain properties correspond to the structural connectivity input that is shared between the models, and certain dynamic properties relate more to the mathematical description of the brain network model. We conclude that the dynamic properties that explicitly examine patterns of signal as a function of time rather than spatial coordination between different brain regions in the rs-fMRI signal seem to provide the largest contrasts between different BNMs and the unknown empirical dynamical system. Our results will be useful in constraining and developing more realistic simulations of whole-brain activity.

Quasi-periodic patterns of brain activity in individuals with attention-deficit/hyperactivity disorder

Individuals with attention-deficit/hyperactivity disorder have disrupted functional connectivity in the default mode and task positive networks. Traditional fMRI analysis techniques that focus on 'static' changes in functional connectivity have been successful in identifying differences between healthy controls and individuals with ADHD. However, such analyses are unable to explain the mechanisms behind the functional connectivity differences observed. Here, we study dynamic changes in functional connectivity in individuals with ADHD through investigation of quasi-periodic patterns (QPPs). QPPs are reliably recurring low-frequency spatiotemporal patterns in the brain linked to infra-slow electrical activity. They have been shown to contribute to functional connectivity observed through static analysis techniques. We find that QPPs contribute to functional connectivity specifically in regions that are disrupted in individuals with ADHD. Individuals with ADHD also show differences in the spatiotemporal pattern observed within the QPPs. This difference results in a weaker contribution of QPPs to functional connectivity in the default mode and task positive networks. We conclude that quasi-periodic patterns provide insight into the mechanisms behind functional connectivity differences seen in individuals with ADHD. This allows for a better understanding of the etiology of the disorder and development of effective treatments.

Cognitive and Physiologic Impacts of the Infraslow Oscillation

Brain states are traditionally recognized via sleep-wake cycles, but modern neuroscience is beginning to identify many sub-states within these larger arousal types. Multiple lines of converging evidence now point to the infraslow oscillation (ISO) as a mediator of brain sub-states, with impacts ranging from the creation of resting state networks (RSNs) in awake subjects to interruptions in neural activity during sleep. This review will explore first the basic characteristics of the ISO in human subjects before reviewing findings in sleep and in animals. Networks of consistently correlated brain regions known as RSNs seen in human functional neuroimaging studies oscillate together at infraslow frequencies. The infraslow rhythm subdivides nonREM in a manner that may correlate with plasticity. The mechanism of this oscillation may be found in the thalamus and may ultimately come from glial cells. Finally, I review the functional impacts of ISOs on brain phenomena ranging from higher frequency oscillations, to brain networks, to information representation and cognitive performance. ISOs represent a relatively understudied phenomenon with wide effects on the brain and behavior.

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.

Spontaneous cognitive processes and the behavioral validation of time-varying brain connectivity

In cognitive neuroscience, focus is commonly placed on associating brain function with changes in objectively measured external stimuli or with actively generated cognitive processes. In everyday life, however, many forms of cognitive processes are initiated spontaneously, without an individual's active effort and without explicit manipulation of behavioral state. Recently, there has been increased emphasis, especially in functional neuroimaging research, on spontaneous correlated activity among spatially segregated brain regions (intrinsic functional connectivity) and, more specifically, on intraindividual fluctuations of such correlated activity on various time scales (time-varying functional connectivity). In this Perspective, we propose that certain subtypes of spontaneous cognitive processes are detectable in time-varying functional connectivity measurements. We define these subtypes of spontaneous cognitive processes and review evidence of their representations in time-varying functional connectivity from studies of attentional fluctuations, memory reactivation, and effects of baseline states on subsequent perception. Moreover, we describe how these studies are critical to validating the use of neuroimaging tools (e.g., fMRI) for assessing ongoing brain network dynamics. We conclude that continued investigation of the behavioral relevance of time-varying functional connectivity will be beneficial both in the development of comprehensive neural models of cognition, and in informing on best practices for studying brain network dynamics.

Definition and Classification of Periodic and Rhythmic Patterns

The growing use of continuous video-EEG recording in the inpatient setting, in particular in patients with varying degrees of encephalopathy, has yielded a window to the brain with an excellent temporal resolution. This increasingly available tool has become more than an instrument to detect nonconvulsive seizures (its primary use), and clinical indications span from ischemia detection in acute brain injuries, neuroprognostication of comatose patients, to monitoring the degree of encephalopathy. In this context, abnormal findings such as periodic discharges and rhythmic delta activity were increasingly recognized; however, significant subjectivity remained in the interpretation of these findings pertaining to key features regarding their spatial involvement, prevalence of occurrence, duration, associated morphologic features, and behavior. In 2005, the American Clinical Neurophysiology Society proposed standardized definitions and classification of electroencephalographic rhythmic and periodic patterns. This was subsequently revised in 2011 and in 2012 and is now being used by centers worldwide, with the final version published in early 2013 as an official guideline of the ACNS. The resulting uniform terminology has allowed for significant advances in the understanding of the pathophysiology, epileptogenic potential, and overall clinical implication of these patterns. Investigators across multiple institutions are now able to collaborate while exploring diagnostic and therapeutic algorithms to these patterns, an effort that may soon provide definitive evidence guiding treating clinicians on the management of these patients.

A randomised, double-blind, placebo-controlled parallel trial of closed-loop infraslow brain training in food addiction

The posterior cingulate cortex (PCC) is involved in food craving in obese food addicted individuals. This randomised, double-blind, placebo-controlled parallel study explored the potential therapeutic effects of infraslow neurofeedback (ISF-NF) on food craving targeting the PCC in obese women with symptoms of food addiction. Participants received six sessions of either ISF-NF (n = 11) or placebo (n = 10) over a three-week period. There were no reported adverse effects. Electrophysiologically, there were significant increases in infraslow activity (p = 0.0002) and infraslow/beta nesting (p < 0.001) in the PCC in the ISF-NF group (mean r = 0.004 ± 0.002) compared to placebo (mean r = 0.02 ± 0.002) two days after the last intervention. Also, there was a significant decrease in different dimensions of state food craving compared to baseline and to placebo. Findings suggest that source localized IFS-NF results in electrophysiological changes and may be associated with reduced food craving. This trial is registered at www.anzctr.org.au , identifier, ACTRN12617000601336. This study was funded by the Otago Medical Research Grant: CT375.

Quasi-Periodic Patterns of Neural Activity improve Classification of Alzheimer's Disease in Mice

Resting state (rs)fMRI allows measurement of brain functional connectivity and has identified default mode (DMN) and task positive (TPN) network disruptions as promising biomarkers for Alzheimer's disease (AD). Quasi-periodic patterns (QPPs) of neural activity describe recurring spatiotemporal patterns that display DMN with TPN anti-correlation. We reasoned that QPPs could provide new insights into AD network dysfunction and improve disease diagnosis. We therefore used rsfMRI to investigate QPPs in old TG2576 mice, a model of amyloidosis, and age-matched controls. Multiple QPPs were determined and compared across groups. Using linear regression, we removed their contribution from the functional scans and assessed how they reflected functional connectivity. Lastly, we used elastic net regression to determine if QPPs improved disease classification. We present three prominent findings: (1) Compared to controls, TG2576 mice were marked by opposing neural dynamics in which DMN areas were anti-correlated and displayed diminished anti-correlation with the TPN. (2) QPPs reflected lowered DMN functional connectivity in TG2576 mice and revealed significantly decreased DMN-TPN anti-correlations. (3) QPP-derived measures significantly improved classification compared to conventional functional connectivity measures. Altogether, our findings provide insight into the neural dynamics of aberrant network connectivity in AD and indicate that QPPs might serve as a translational diagnostic tool.

Time-Resolved Resting-State Functional Magnetic Resonance Imaging Analysis: Current Status, Challenges, and New Directions

Time-resolved analysis of resting-state functional magnetic resonance imaging (rs-fMRI) data allows researchers to extract more information about brain function than traditional functional connectivity analysis, yet a number of challenges in data analysis and interpretation remain. This article briefly summarizes common methods for time-resolved analysis and presents some of the pressing issues and opportunities in the field. From there, the discussion moves to interpretation of the network dynamics observed with rs-fMRI and the role that rs-fMRI can play in elucidating the large-scale organization of brain activity.

Infraslow Electroencephalographic and Dynamic Resting State Network Activity

A number of studies have linked the blood oxygenation level dependent (BOLD) signal to electroencephalographic (EEG) signals in traditional frequency bands (δ, θ, α, β, and γ), but the relationship between BOLD and its direct frequency correlates in the infraslow band (<1 Hz) has been little studied. Previously, work in rodents showed that infraslow local field potentials play a role in functional connectivity, particularly in the dynamic organization of large-scale networks. To examine the relationship between infraslow activity and network dynamics in humans, direct current (DC) EEG and resting state magnetic resonance imaging data were acquired simultaneously. The DC EEG signals were correlated with the BOLD signal in patterns that resembled resting state networks. Subsequent dynamic analysis showed that the correlation between DC EEG and the BOLD signal varied substantially over time, even within individual subjects. The variation in DC EEG appears to reflect the time-varying contribution of different resting state networks. Furthermore, some of the patterns of DC EEG and BOLD correlation are consistent with previous work demonstrating quasiperiodic spatiotemporal patterns of large-scale network activity in resting state. These findings demonstrate that infraslow electrical activity is linked to BOLD fluctuations in humans and that it may provide a basis for large-scale organization comparable to that observed in animal studies.

Co-activation patterns in resting-state fMRI signals

The brain is a complex system that integrates and processes information across multiple time scales by dynamically coordinating activities over brain regions and circuits. Correlations in resting-state functional magnetic resonance imaging (rsfMRI) signals have been widely used to infer functional connectivity of the brain, providing a metric of functional associations that reflects a temporal average over an entire scan (typically several minutes or longer). Not until recently was the study of dynamic brain interactions at much shorter time scales (seconds to minutes) considered for inference of functional connectivity. One method proposed for this objective seeks to identify and extract recurring co-activation patterns (CAPs) that represent instantaneous brain configurations at single time points. Here, we review the development and recent advancement of CAP methodology and other closely related approaches, as well as their applications and associated findings. We also discuss the potential neural origins and behavioral relevance of CAPs, along with methodological issues and future research directions in the analysis of fMRI co-activation patterns.

Incorporating spatial constraint in co-activation pattern analysis to explore the dynamics of resting-state networks: An application to Parkinson's disease

The dynamics of the brain's intrinsic networks have been recently studied using co-activation pattern (CAP) analysis. The CAP method relies on few model assumptions and CAP-based measurements provide quantitative information of network temporal dynamics. One limitation of existing CAP-related methods is that the computed CAPs share considerable spatial overlap that may or may not be functionally distinct relative to specific network dynamics. To more accurately describe network dynamics with spatially distinct CAPs, and to compare network dynamics between different populations, a novel data-driven CAP group analysis method is proposed in this study. In the proposed method, a dominant-CAP (d-CAP) set is synthesized across CAPs from multiple clustering runs for each group with the constraint of low spatial similarities among d-CAPs. Alternating d-CAPs with less overlapping spatial patterns can better capture overall network dynamics. The number of d-CAPs, the temporal fraction and spatial consistency of each d-CAP, and the subject-specific switching probability among all d-CAPs are then calculated for each group and used to compare network dynamics between groups. The spatial dissimilarities among d-CAPs computed with the proposed method were first demonstrated using simulated data. High consistency between simulated ground-truth and computed d-CAPs was achieved, and detailed comparisons between the proposed method and existing CAP-based methods were conducted using simulated data. In an effort to physiologically validate the proposed technique and investigate network dynamics in a relevant brain network disorder, the proposed method was then applied to data from the Parkinson's Progression Markers Initiative (PPMI) database to compare the network dynamics in Parkinson's disease (PD) and normal control (NC) groups. Fewer d-CAPs, skewed distribution of temporal fractions of d-CAPs, and reduced switching probabilities among final d-CAPs were found in most networks in the PD group, as compared to the NC group. Furthermore, an overall negative association between switching probability among d-CAPs and disease severity was observed in most networks in the PD group as well. These results expand upon previous findings from in vivo electrophysiological recording studies in PD. Importantly, this novel analysis also demonstrates that changes in network dynamics can be measured using resting-state fMRI data from subjects with early stage PD.

Stereotypical modulations in dynamic functional connectivity explained by changes in BOLD variance

Spontaneous activity measured in human subject under the absence of any task exhibits complex patterns of correlation that largely correspond to large-scale functional topographies obtained with a wide variety of cognitive and perceptual tasks. These "resting state networks" (RSNs) fluctuate over time, forming and dissolving on the scale of seconds to minutes. While these fluctuations, most prominently those of the default mode network, have been linked to cognitive function, it remains unclear whether they result from random noise or whether they index a nonstationary process which could be described as state switching. In this study, we use a sliding windows-approach to relate temporal dynamics of RSNs to global modulations in correlation and BOLD variance. We compare empirical data, phase-randomized surrogate data, and data simulated with a stationary model. We find that RSN time courses exhibit a large amount of coactivation in all three cases, and that the modulations in their activity are closely linked to global dynamics of the underlying BOLD signal. We find that many properties of the observed fluctuations in FC and BOLD, including their ranges and their correlations amongst each other, are explained by fluctuations around the average FC structure. However, we also report some interesting characteristics that clearly support nonstationary features in the data. In particular, we find that the brain spends more time in the troughs of modulations than can be expected from stationary dynamics.

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.

PEAR: PEriodic And fixed Rank separation for fast fMRI

PURPOSE

In functional MRI (fMRI), faster acquisition via undersampling of data can improve the spatial-temporal resolution trade-off and increase statistical robustness through increased degrees-of-freedom. High-quality reconstruction of fMRI data from undersampled measurements requires proper modeling of the data. We present an fMRI reconstruction approach based on modeling the fMRI signal as a sum of periodic and fixed rank components, for improved reconstruction from undersampled measurements.

METHODS

The proposed approach decomposes the fMRI signal into a component which has a fixed rank and a component consisting of a sum of periodic signals which is sparse in the temporal Fourier domain. Data reconstruction is performed by solving a constrained problem that enforces a fixed, moderate rank on one of the components, and a limited number of temporal frequencies on the other. Our approach is coined PEAR - PEriodic And fixed Rank separation for fast fMRI.

RESULTS

Experimental results include purely synthetic simulation, a simulation with real timecourses and retrospective undersampling of a real fMRI dataset. Evaluation was performed both quantitatively and visually versus ground truth, comparing PEAR to two additional recent methods for fMRI reconstruction from undersampled measurements. Results demonstrate PEAR's improvement in estimating the timecourses and activation maps versus the methods compared against at acceleration ratios of R = 8,10.66 (for simulated data) and R = 6.66,10 (for real data).

CONCLUSIONS

This paper presents PEAR, an undersampled fMRI reconstruction approach based on decomposing the fMRI signal to periodic and fixed rank components. PEAR results in reconstruction with higher fidelity than when using a fixed-rank based model or a conventional Low-rank + Sparse algorithm. We have shown that splitting the functional information between the components leads to better modeling of fMRI, over state-of-the-art methods.

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.

PEAR: PEriodic and ApeRiodic signal separation for fast FMRI

Undersampling of functional MRI (fMRI) data leads to increased temporal resolution, as it allows shorter acquisition time per frame. High quality reconstruction of fMRI data from undersampled measurements requires proper modeling of the fMRI data. Recent publications suggest that the fMRI signal is a superposition of periodic and aperiodic signals. In this paper we develop an fMRI reconstruction approach based on this modeling. The fMRI data is assumed to be composed of two components: a component that holds a sum of periodic signals which is sparse in the temporal Fourier domain and an component that holds the remaining imaging information (consisting of the background and aperiodic signals) which has low rank. Data reconstruction is done by solving a constrained problem that enforces a fixed, moderate rank on one of the components, and a limited number of temporal frequencies on the other. Our approach is coined PEAR - PEriodic and ApeRiodic signal separation for fast fMRI. Experimental results are based on fMRI reconstruction using realistic timecourses. Evaluation was performed both quantitatively and visually versus ground truth. Results demonstrate PEAR's improvement in estimating the realistic timecourses versus state-of-the-art approaches at acceleration ratio of R=16.6.

The behavioral and cognitive relevance of time-varying, dynamic changes in functional connectivity

Recent advances in neuroimaging methods and analysis have led to an expanding body of research that investigates how large-scale brain network organization dynamically adapts to changes in one's environment, including both internal state changes and external stimulation. It is now possible to detect changes in functional connectivity that occur on the order of seconds, both during an unconstrained resting state and during the performance of constrained cognitive tasks. It is thought that these dynamic, time-varying changes in functional connectivity, often referred to as dynamic functional connectivity (dFC), include features that are relevant to behavior and cognition. This review summarizes four aspects of the nascent literature directly testing that assumption: 1) how changes in functional network organization on the order of task blocks relate to differences in task demands and to cognitive ability; 2) how differences in dFC variability between different contexts relate to cognitive demands and behavioral performance; 3) how ongoing fluctuations in dFC impact perception and attention; and 4) how different patterns of dFC correspond to individual differences in cognition. The review ends by discussing promising directions for future research in this field. First, it comments on how dFC analyses can help to elucidate the mechanisms of healthy cognition. Next, it describes how dFC processes may be disrupted in disease, and how probing such dysfunction can increase understanding of neural etiology, as well as behavioral and cognitive impairments, observed in psychiatric and neurologic populations. Last, it considers the potential for computational models to uncover neuronal mechanisms of dFC, and how both healthy cognition and disease emerge from network dynamics.

Cortex-wide BOLD fMRI activity reflects locally-recorded slow oscillation-associated calcium waves

Spontaneous slow oscillation-associated slow wave activity represents an internally generated state which is characterized by alternations of network quiescence and stereotypical episodes of neuronal activity - slow wave events. However, it remains unclear which macroscopic signal is related to these active periods of the slow wave rhythm. We used optic fiber-based calcium recordings of local neural populations in cortex and thalamus to detect neurophysiologically defined slow calcium waves in isoflurane anesthetized rats. The individual slow wave events were used for an event-related analysis of simultaneously acquired whole-brain BOLD fMRI. We identified BOLD responses directly related to onsets of slow calcium waves, revealing a cortex-wide BOLD correlate: the entire cortex was engaged in this specific type of slow wave activity. These findings demonstrate a direct relation of defined neurophysiological events to a specific BOLD activity pattern and were confirmed for ongoing slow wave activity by independent component and seed-based analyses.

When a person is in a deep non-dreaming sleep, neurons in their brain alternate slowly between periods of silence and periods of activity. This gives rise to low-frequency brain rhythms called slow waves, which are thought to help stabilize memories. Slow wave activity can be detected on multiple scales, from the pattern of electrical impulses sent by an individual neuron to the collective activity of the brain’s entire outer layer, the cortex. But does slow wave activity in an individual group of neurons in the cortex affect the activity of the rest of the brain?

To find out, Schwalm, Schmid, Wachsmuth et al. took advantage of the fact that slow waves also occur under general anesthesia, and placed anesthetized rats inside miniature whole-brain scanners. A small region of cortex in each rat had been injected with a dye that fluoresces whenever the neurons in that region are active. An optical fiber was lowered into the rat’s brain to transmit the fluorescence signals to a computer. Monitoring these signals while the animals lay inside the scanner revealed that slow-wave activity in any one group of cortical neurons was accompanied by slow-wave activity across the cortex as a whole. This relationship was seen only for slow waves, and not for other brain rhythms.

Slow waves seem to occur in all species of animal with a backbone, and in both healthy and diseased brains. While it is not possible to inject fluorescent dyes into the human brain, it is possible to monitor neuronal activity using electrodes. Comparing local electrode recordings with measures of whole-brain activity from scanners could thus allow similar experiments to be performed in people. There is growing evidence – from animal models and from studies of patients – that slow waves may be altered in Alzheimer’s disease. Further work is required to determine whether detecting these changes could help diagnose disease at earlier stages, and whether reversing them may have therapeutic potential.

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.

Frequency-phase analysis of resting-state functional MRI

We describe an analysis method that characterizes the correlation between coupled time-series functions by their frequencies and phases. It provides a unified framework for simultaneous assessment of frequency and latency of a coupled time-series. The analysis is demonstrated on resting-state functional MRI data of 34 healthy subjects. Interactions between fMRI time-series are represented by cross-correlation (with time-lag) functions. A general linear model is used on the cross-correlation functions to obtain the frequencies and phase-differences of the original time-series. We define symmetric, antisymmetric and asymmetric cross-correlation functions that correspond respectively to in-phase, 90° out-of-phase and any phase difference between a pair of time-series, where the last two were never introduced before. Seed maps of the motor system were calculated to demonstrate the strength and capabilities of the analysis. Unique types of functional connections, their dominant frequencies and phase-differences have been identified. The relation between phase-differences and time-delays is shown. The phase-differences are speculated to inform transfer-time and/or to reflect a difference in the hemodynamic response between regions that are modulated by neurotransmitters concentration. The analysis can be used with any coupled functions in many disciplines including electrophysiology, EEG or MEG in neuroscience.

The global signal in fMRI: Nuisance or Information?

The global signal is widely used as a regressor or normalization factor for removing the effects of global variations in the analysis of functional magnetic resonance imaging (fMRI) studies. However, there is considerable controversy over its use because of the potential bias that can be introduced when it is applied to the analysis of both task-related and resting-state fMRI studies. In this paper we take a closer look at the global signal, examining in detail the various sources that can contribute to the signal. For the most part, the global signal has been treated as a nuisance term, but there is growing evidence that it may also contain valuable information. We also examine the various ways that the global signal has been used in the analysis of fMRI data, including global signal regression, global signal subtraction, and global signal normalization. Furthermore, we describe new ways for understanding the effects of global signal regression and its relation to the other approaches.

Combined spatiotemporal ICA (stICA) for continuous and dynamic lag structure analysis of MREG data

This study investigated lag structure in the resting-state fMRI by applying a novel independent component (ICA) method to magnetic resonance encephalography (MREG) data. Briefly, the spatial ICA (sICA) was used for defining the frontal and back nodes of the default mode network (DMN), and the temporal ICA (tICA), which is enabled by the high temporal resolution of MREG (TR=100ms), was used to separate both neuronal and physiological components of these two spatial map regions. Subsequently, lag structure was investigated between the frontal (DMNvmpf) and posterior (DMNpcc) DMN nodes using both conventional method with all-time points and a sliding-window approach. A rigorous noise exclusion criterion was applied for tICs to remove physiological pulsations, motion and system artefacts. All the de-noised tICs were used to calculate the null-distributions both for expected lag variability over time and over subjects. Lag analysis was done for the three highest correlating denoised tICA pairs. Mean time lag of 0.6s (± 0.5 std) and mean absolute correlation of 0.69 (± 0.08) between the highest correlating tICA pairs of DMN nodes was observed throughout the whole analyzed period. In dynamic 2min window analysis, there was large variability over subjects as ranging between 1-10sec. Directionality varied between these highly correlating sources an average 28.8% of the possible number of direction changes. The null models show highly consistent correlation and lag structure between DMN nodes both in continuous and dynamic analysis. The mean time lag of a null-model over time between all denoised DMN nodes was 0.0s and, thus the probability of having either DMNpcc or DMNvmpf as a preceding component is near equal. All the lag values of highest correlating tICA pairs over subjects lie within the standard deviation range of a null-model in whole time window analysis, supporting the earlier findings that there is a consistent temporal lag structure across groups of individuals. However, in dynamic analysis, there are lag values exceeding the threshold of significance of a null-model meaning that there might be biologically meaningful variation in this measure. Taken together the variability in lag and the presence of high activity peaks during strong connectivity indicate that individual avalanches may play an important role in defining dynamic independence in resting state connectivity within networks.

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.

Dynamic Connectivity Patterns in Conscious and Unconscious Brain

Brain functional connectivity undergoes dynamic changes from the awake to unconscious states. However, how the dynamics of functional connectivity patterns are linked to consciousness at the behavioral level remains elusive. In this study, we acquired resting-state functional magnetic resonance imaging data during wakefulness and graded levels of consciousness in rats. Data were analyzed using a dynamic approach combining the sliding window method and k-means clustering. Our results demonstrate that whole-brain networks contained several quasi-stable patterns that dynamically recurred from the awake state into anesthetized states. Remarkably, two brain connectivity states with distinct spatial similarity to the structure of anatomical connectivity were strongly biased toward high and low consciousness levels, respectively. These results provide compelling neuroimaging evidence linking the dynamics of whole-brain functional connectivity patterns and states of consciousness at the behavioral level.

Resting-state functional MRI as a tool for evaluating brain hemodynamic responsiveness to external stimuli in rats

PURPOSE

Anesthesia is a major confounding factor in functional MRI (fMRI) experiments attributed to its effects on brain function. Recent evidence suggests that parameters obtained with resting-state fMRI (rs-fMRI) are coupled with anesthetic depth. Therefore, we investigated whether parameters obtained with rs-fMRI, such as functional connectivity (FC), are also directly related to blood-oxygen-level-dependent (BOLD) responses.

METHODS

A simple rs-fMRI protocol was implemented in a pharmacological fMRI study to evaluate the coupling between hemodynamic responses and FC under five anesthetics (α-chloralose, isoflurane, medetomidine, thiobutabarbital, and urethane). Temporal change in the FC was evaluated at 1-hour interval. Supplementary forepaw stimulation experiments were also conducted.

RESULTS

Under thiobutabarbital anesthesia, FC was clearly coupled with nicotine-induced BOLD responses. Good correlation values were also obtained under isoflurane and medetomidine anesthesia. The observations in the thiobutabarbital group were supported by forepaw stimulation experiments. Additionally, the rs-fMRI protocol revealed significant temporal changes in the FC in the α-chloralose, thiobutabarbital, and urethane groups.

CONCLUSION

Our results suggest that FC can be used to estimate brain hemodynamic responsiveness to stimuli and evaluate the level and temporal changes of anesthesia. Therefore, analysis of the fMRI baseline signal may be highly valuable tool for controlling the outcome of preclinical fMRI experiments. Magn Reson Med 78:1136-1146, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

Ongoing Slow Fluctuations in V1 Impact on Visual Perception

The human brain’s ongoing activity is characterized by intrinsic networks of coherent fluctuations, measured for example with correlated functional magnetic resonance imaging signals. So far, however, the brain processes underlying this ongoing blood oxygenation level dependent (BOLD) signal orchestration and their direct relevance for human behavior are not sufficiently understood. In this study, we address the question of whether and how ongoing BOLD activity within intrinsic occipital networks impacts on conscious visual perception. To this end, backwardly masked targets were presented in participants’ left visual field only, leaving the ipsi-lateral occipital areas entirely free from direct effects of task throughout the experiment. Signal time courses of ipsi-lateral BOLD fluctuations in visual areas V1 and V2 were then used as proxies for the ongoing contra-lateral BOLD activity within the bilateral networks. Magnitude and phase of these fluctuations were compared in trials with and without conscious visual perception, operationalized by means of subjective confidence ratings. Our results show that ipsi-lateral BOLD magnitudes in V1 were significantly higher at times of peak response when the target was perceived consciously. A significant difference between conscious and non-conscious perception with regard to the pre-target phase of an intrinsic-frequency regime suggests that ongoing V1 fluctuations exert a decisive impact on the access to consciousness already before stimulation. Both effects were absent in V2. These results thus support the notion that ongoing slow BOLD activity within intrinsic networks covering V1 represents localized processes that modulate the degree of readiness for the emergence of visual consciousness.

Multiscale network activity in resting state fMRI

The brain is inherently multiscalar in both space and time. We argue that this multiscalar nature is reflected in the blood oxygenation level dependent (BOLD) fluctuations used to map functional connectivity. We present evidence that global fluctuations in activity, quasiperiodic spatiotemporal patterns, and aperiodic time-varying activity coexist within the BOLD signal. These processes can be separated using careful analysis and appear to reflect electrical activity on similar scales, suggesting that the BOLD signal fluctuations can provide novel insight into the functional architecture of the brain.

Ultra-fast magnetic resonance encephalography of physiological brain activity - Glymphatic pulsation mechanisms?

The theory on the glymphatic convection mechanism of cerebrospinal fluid holds that cardiac pulsations in part pump cerebrospinal fluid from the peri-arterial spaces through the extracellular tissue into the peri-venous spaces facilitated by aquaporin water channels. Since cardiac pulses cannot be the sole mechanism of glymphatic propulsion, we searched for additional cerebrospinal fluid pulsations in the human brain with ultra-fast magnetic resonance encephalography. We detected three types of physiological mechanisms affecting cerebral cerebrospinal fluid pulsations: cardiac, respiratory, and very low frequency pulsations. The cardiac pulsations induce a negative magnetic resonance encephalography signal change in peri-arterial regions that extends centrifugally and covers the brain in ≈1 Hz cycles. The respiratory ≈0.3 Hz pulsations are centripetal periodical pulses that occur dominantly in peri-venous areas. The third type of pulsation was very low frequency (VLF 0.001-0.023 Hz) and low frequency (LF 0.023-0.73 Hz) waves that both propagate with unique spatiotemporal patterns. Our findings using critically sampled magnetic resonance encephalography open a new view into cerebral fluid dynamics. Since glymphatic system failure may precede protein accumulations in diseases such as Alzheimer's dementia, this methodological advance offers a novel approach to image brain fluid dynamics that potentially can enable early detection and intervention in neurodegenerative diseases.

Low frequency steady-state brain responses modulate large scale functional networks in a frequency-specific means

Neural oscillations are essential for brain functions. Research has suggested that the frequency of neural oscillations is lower for more integrative and remote communications. In this vein, some resting-state studies have suggested that large scale networks function in the very low frequency range (<1 Hz). However, it is difficult to determine the frequency characteristics of brain networks because both resting-state studies and conventional frequency tagging approaches cannot simultaneously capture multiple large scale networks in controllable cognitive activities. In this preliminary study, we aimed to examine whether large scale networks can be modulated by task-induced low frequency steady-state brain responses (lfSSBRs) in a frequency-specific pattern. In a revised attention network test, the lfSSBRs were evoked in the triple network system and sensory-motor system, indicating that large scale networks can be modulated in a frequency tagging way. Furthermore, the inter- and intranetwork synchronizations as well as coherence were increased at the fundamental frequency and the first harmonic rather than at other frequency bands, indicating a frequency-specific modulation of information communication. However, there was no difference among attention conditions, indicating that lfSSBRs modulate the general attention state much stronger than distinguishing attention conditions. This study provides insights into the advantage and mechanism of lfSSBRs. More importantly, it paves a new way to investigate frequency-specific large scale brain activities.

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.

Different dynamic resting state fMRI patterns are linked to different frequencies of neural activity

Resting state functional magnetic resonance imaging (rsfMRI) results have indicated that network mapping can contribute to understanding behavior and disease, but it has been difficult to translate the maps created with rsfMRI to neuroelectrical states in the brain. Recently, dynamic analyses have revealed multiple patterns in the rsfMRI signal that are strongly associated with particular bands of neural activity. To further investigate these findings, simultaneously recorded invasive electrophysiology and rsfMRI from rats were used to examine two types of electrical activity (directly measured low-frequency/infraslow activity and band-limited power of higher frequencies) and two types of dynamic rsfMRI (quasi-periodic patterns or QPP, and sliding window correlation or SWC). The relationship between neural activity and dynamic rsfMRI was tested under three anesthetic states in rats: dexmedetomidine and high and low doses of isoflurane. Under dexmedetomidine, the lightest anesthetic, infraslow electrophysiology correlated with QPP but not SWC, whereas band-limited power in higher frequencies correlated with SWC but not QPP. Results were similar under isoflurane; however, the QPP was also correlated to band-limited power, possibly due to the burst-suppression state induced by the anesthetic agent. The results provide additional support for the hypothesis that the two types of dynamic rsfMRI are linked to different frequencies of neural activity, but isoflurane anesthesia may make this relationship more complicated. Understanding which neural frequency bands appear as particular dynamic patterns in rsfMRI may ultimately help isolate components of the rsfMRI signal that are of interest to disorders such as schizophrenia and attention deficit disorder.

Functional connectivity arises from a slow rhythmic mechanism

The mechanism underlying temporal correlations among blood oxygen level-dependent signals is unclear. We used oxygen polarography to better characterize oxygen fluctuations and their correlation and to gain insight into the driving mechanism. The power spectrum of local oxygen fluctuations is inversely proportional to frequency raised to a power (1/f) raised to the beta, with an additional positive band-limited component centered at 0.06 Hz. In contrast, the power of the correlated oxygen signal is band limited from ∼ 0.01 Hz to 0.4 Hz with a peak at 0.06 Hz. These results suggest that there is a band-limited mechanism (or mechanisms) driving interregional oxygen correlation that is distinct from the mechanism(s) driving local (1/f) oxygen fluctuations. Candidates for driving interregional oxygen correlation include rhythmic or pseudo-oscillatory mechanisms.

Neuroplasticity in Human Alcoholism: Studies of Extended Abstinence with Potential Treatment Implications

Alcoholism is characterized by a lack of control over excessive alcohol consumption despite significant negative consequences. This impulsive and compulsive behavior may be related to functional abnormalities within networks of brain regions responsible for how we make decisions. The abnormalities may result in strengthened networks related to appetitive drive-or the need to fulfill desires-and simultaneously weakened networks that exercise control over behaviors. Studies using functional magnetic resonance imaging (fMRI) in abstinent alcoholics suggest that abstinence is associated with changes in the tone of such networks, decreasing resting tone in appetitive drive networks, and increasing resting tone in inhibitory control networks to support continued abstinence. Identifying electroencephalographic (EEG) measures of resting tone in these networks initially identified using fMRI, and establishing in longitudinal studies that these abstinence-related changes in network tone are progressive would motivate treatment initiatives to facilitate these changes in network tone, thereby supporting successful ongoing abstinence.

Assessing dynamic brain graphs of time-varying connectivity in fMRI data: application to healthy controls and patients with schizophrenia

Graph theory-based analysis has been widely employed in brain imaging studies, and altered topological properties of brain connectivity have emerged as important features of mental diseases such as schizophrenia. However, most previous studies have focused on graph metrics of stationary brain graphs, ignoring that brain connectivity exhibits fluctuations over time. Here we develop a new framework for accessing dynamic graph properties of time-varying functional brain connectivity in resting-state fMRI data and apply it to healthy controls (HCs) and patients with schizophrenia (SZs). Specifically, nodes of brain graphs are defined by intrinsic connectivity networks (ICNs) identified by group independent component analysis (ICA). Dynamic graph metrics of the time-varying brain connectivity estimated by the correlation of sliding time-windowed ICA time courses of ICNs are calculated. First- and second-level connectivity states are detected based on the correlation of nodal connectivity strength between time-varying brain graphs. Our results indicate that SZs show decreased variance in the dynamic graph metrics. Consistent with prior stationary functional brain connectivity works, graph measures of identified first-level connectivity states show lower values in SZs. In addition, more first-level connectivity states are disassociated with the second-level connectivity state which resembles the stationary connectivity pattern computed by the entire scan. Collectively, the findings provide new evidence about altered dynamic brain graphs in schizophrenia, which may underscore the abnormal brain performance in this mental illness.

Dynamic resting state functional connectivity in awake and anesthetized rodents

Since its introduction, resting-state functional magnetic resonance imaging (rsfMRI) has been a powerful tool for investigating functional neural networks in both normal and pathological conditions. When measuring resting-state functional connectivity (RSFC), most rsfMRI approaches do not consider its temporal variations and thus only provide the averaged RSFC over the scan time. Recently, there has been a surge of interest to investigate the dynamic characteristics of RSFC in humans, and promising results have been yielded. However, our knowledge regarding the dynamic RSFC in animals remains sparse. In the present study we utilized the single-volume co-activation method to systematically study the dynamic properties of RSFC within the networks of infralimbic cortex (IL) and primary somatosensory cortex (S1) in both awake and anesthetized rats. Our data showed that both IL and S1 networks could be decomposed into several spatially reproducible but temporally changing co-activation patterns (CAPs), suggesting that dynamic RSFC was indeed a characteristic feature in rodents. In addition, we demonstrated that anesthesia profoundly impacted the dynamic RSFC of neural circuits subserving cognitive and emotional functions but had less effects on sensorimotor systems. Finally, we examined the temporal characteristics of each CAP, and found that individual CAPs exhibited consistent temporal evolution patterns. Together, these results suggest that dynamic RSFC might be a general phenomenon in vertebrate animals. In addition, this study has paved the way for further understanding the alterations of dynamic RSFC in animal models of brain disorders.