Infraslow LFP correlates to resting-state fMRI BOLD signals

Short timescale abnormalities in the states of spontaneous synchrony in the functional neural networks in Alzheimer's disease

Alzheimer's disease (AD) is a prevalent neurodegenerative condition that can lead to severe cognitive and functional deterioration. Functional magnetic resonance imaging (fMRI) revealed abnormalities in AD in intrinsic synchronization between spatially separate regions in the so-called default mode network (DMN) of the brain. To understand the relationship between this disruption in large-scale synchrony and the cognitive impairment in AD, it is critical to determine whether and how the deficit in the low frequency hemodynamic fluctuations recorded by fMRI translates to much faster timescales of memory and other cognitive processes. The present study employed magnetoencephalography (MEG) and a Hidden Markov Model (HMM) approach to estimate spontaneous synchrony variations in the functional neural networks with high temporal resolution. In a group of cognitively healthy (CH) older adults, we found transient (mean duration of 150-250 ms) network activity states, which were visited in a rapid succession, and were characterized by spatially coordinated changes in the amplitude of source-localized electrophysiological oscillations. The inferred states were similar to those previously observed in younger participants using MEG, and the estimated cortical source distributions of the state-specific activity resembled the classic functional neural networks, such as the DMN. In patients with AD, inferred network states were different from those of the CH group in short-scale timing and oscillatory features. The state of increased oscillatory amplitudes in the regions overlapping the DMN was visited less often in AD and for shorter periods of time, suggesting that spontaneous synchronization in this network was less likely and less stable in the patients. During the visits to this state, in some DMN nodes, the amplitude change in the higher-frequency (8-30 Hz) oscillations was less robust in the AD than CH group. These findings highlight relevance of studying short-scale temporal evolution of spontaneous activity in functional neural networks to understanding the AD pathophysiology. Capacity of flexible intrinsic synchronization in the DMN may be crucial for memory and other higher cognitive functions. Our analysis yielded metrics that quantify distinct features of the neural synchrony disorder in AD and may offer sensitive indicators of the neural network health for future investigations.

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.

Effects of Tissue-Specific Functional Magnetic Resonance Imaging Signal Regression on Resting-State Functional Connectivity

Neuroimaging studies typically consider white matter as unchanging in different neural and metabolic states. However, a recent study demonstrated that white matter signal regression (WMSR) produced a similar loss of neurometabolic information to global (whole-brain) signal regression (GSR) in resting-state functional magnetic resonance imaging (R-fMRI) data. This was unexpected as the loss of information would normally be attributed to neural activity within gray matter correlating with the global R-fMRI signal. Indeed, WMSR has been suggested as an alternative to avoid such pitfalls in GSR. To address these concerns about tissue-specific regression in R-fMRI data analysis, we performed GSR, WMSR, and gray matter signal regression (GMSR) on R-fMRI data from the 1000 Functional Connectomes Project. We describe several regional and motion-related differences between different types of regressions. However, the overall effects of concern, particularly network-specific alteration of correlation coefficients, are present for all regressions. This suggests that tissue-specific regression is not an adequate strategy to counter pitfalls of GSR. Conversely, if GSR is desired, but the studied disease state excludes either gray matter or white matter from analysis (e.g., due to tissue atrophy), our results indicate that WMSR or GMSR may reproduce the gross effects of GSR.

Fluctuations of the EEG-fMRI correlation reflect intrinsic strength of functional connectivity in default mode network

Both functional magnetic resonance imaging (fMRI) and electrophysiological recordings have revealed that resting-state functional connectivity is temporally variable in human brain. Combined full-band electroencephalography-fMRI (fbEEG-fMRI) studies have shown that infraslow (<.1 Hz) fluctuations in EEG scalp potential are correlated with the blood-oxygen-level-dependent (BOLD) fMRI signals and that also this correlation appears variable over time. Here, we used simultaneous fbEEG-fMRI to test the hypothesis that correlation dynamics between BOLD and fbEEG signals could be explained by fluctuations in the activation properties of resting-state networks (RSNs) such as the extent or strength of their activation. We used ultrafast magnetic resonance encephalography (MREG) fMRI to enable temporally accurate and statistically robust short-time-window comparisons of infra-slow fbEEG and BOLD signals. We found that the temporal fluctuations in the fbEEG-BOLD correlation were dependent on RSN connectivity strength, but not on the mean signal level or magnitude of RSN activation or motion during scanning. Moreover, the EEG-fMRI correlations were strongest when the intrinsic RSN connectivity was strong and close to the pial surface. Conversely, weak fbEEG-BOLD correlations were attributable to periods of less coherent or spatially more scattered intrinsic RSN connectivity, or RSN activation in deeper cerebral structures. The results thus show that the on-average low correlations between infra-slow EEG and BOLD signals are, in fact, governed by the momentary coherence and depth of the underlying RSN activation, and may reach systematically high values with appropriate source activities. These findings further consolidate the notion of slow scalp potentials being directly coupled to hemodynamic fluctuations.

Spontaneous Infra-slow Brain Activity Has Unique Spatiotemporal Dynamics and Laminar Structure

Systems-level organization in spontaneous infra-slow (<0.1Hz) brain activity, measured using blood oxygen signals in fMRI and optical imaging, has become a major theme in the study of neural function in both humans and animal models. Yet the neurophysiological basis of infra-slow activity (ISA) remains unresolved. In particular, is ISA a distinct physiological process, or is it a low-frequency analog of faster neural activity? Here, using whole-cortex calcium/hemoglobin imaging in mice, we show that ISA in each of these modalities travels through the cortex along stereotypical spatiotemporal trajectories that are state dependent (wake versus anesthesia) and distinct from trajectories in delta (1-4 Hz) activity. Moreover, mouse laminar electrophysiology reveals that ISA travels through specific cortical layers and is organized into unique cross-laminar temporal dynamics that are different from higher frequency local field potential activity. These findings suggest that ISA is a distinct neurophysiological process that is reflected in fMRI blood oxygen signals.

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.

Spatial-temporal-spectral EEG patterns of BOLD functional network connectivity dynamics

OBJECTIVE

Growing interest in the examination of large-scale brain network functional connectivity dynamics is accompanied by an effort to find the electrophysiological correlates. The commonly used constraints applied to spatial and spectral domains during electroencephalogram (EEG) data analysis may leave part of the neural activity unrecognized. We propose an approach that blindly reveals multimodal EEG spectral patterns that are related to the dynamics of the BOLD functional network connectivity.

APPROACH

The blind decomposition of EEG spectrogram by parallel factor analysis has been shown to be a useful technique for uncovering patterns of neural activity. The simultaneously acquired BOLD fMRI data were decomposed by independent component analysis. Dynamic functional connectivity was computed on the component's time series using a sliding window correlation, and between-network connectivity states were then defined based on the values of the correlation coefficients. ANOVA tests were performed to assess the relationships between the dynamics of between-network connectivity states and the fluctuations of EEG spectral patterns.

MAIN RESULTS

We found three patterns related to the dynamics of between-network connectivity states. The first pattern has dominant peaks in the alpha, beta, and gamma bands and is related to the dynamics between the auditory, sensorimotor, and attentional networks. The second pattern, with dominant peaks in the theta and low alpha bands, is related to the visual and default mode network. The third pattern, also with peaks in the theta and low alpha bands, is related to the auditory and frontal network.

SIGNIFICANCE

Our previous findings revealed a relationship between EEG spectral pattern fluctuations and the hemodynamics of large-scale brain networks. In this study, we suggest that the relationship also exists at the level of functional connectivity dynamics among large-scale brain networks when no standard spatial and spectral constraints are applied on the EEG data.

Ultra-Slow Single-Vessel BOLD and CBV-Based fMRI Spatiotemporal Dynamics and Their Correlation with Neuronal Intracellular Calcium Signals

Functional MRI has been used to map brain activity and functional connectivity based on the strength and temporal coherence of neurovascular-coupled hemodynamic signals. Here, single-vessel fMRI reveals vessel-specific correlation patterns in both rodents and humans. In anesthetized rats, fluctuations in the vessel-specific fMRI signal are correlated with the intracellular calcium signal measured in neighboring neurons. Further, the blood-oxygen-level-dependent (BOLD) signal from individual venules and the cerebral-blood-volume signal from individual arterioles show correlations at ultra-slow (<0.1 Hz), anesthetic-modulated rhythms. These data support a model that links neuronal activity to intrinsic oscillations in the cerebral vasculature, with a spatial correlation length of ∼2 mm for arterioles. In complementary data from awake human subjects, the BOLD signal is spatially correlated among sulcus veins and specified intracortical veins of the visual cortex at similar ultra-slow rhythms. These data support the use of fMRI to resolve functional connectivity at the level of single vessels.

Spontaneous activity forms a foundation for odor-evoked activation maps in the rat olfactory bulb

Fluctuations in spontaneous activity have been observed by many neuroimaging techniques, but because these resting-state changes are not evoked by stimuli, it is difficult to determine how they relate to task-evoked activations. We conducted multi-modal neuroimaging scans of the rat olfactory bulb, both with and without odor, to examine interaction between spontaneous and evoked activities. Independent component analysis of spontaneous fluctuations revealed resting-state networks, and odor-evoked changes revealed activation maps. We constructed simulated activation maps using resting-state networks that were highly correlated to evoked activation maps. Simulated activation maps derived by intrinsic optical signal (IOS), which covers the dorsal portion of the glomerular sheet, significantly differentiated one odor's evoked activation map from the other two. To test the hypothesis that spontaneous activity of the entire glomerular sheet is relevant for representing odor-evoked activations, we used functional magnetic resonance imaging (fMRI) to map the entire glomerular sheet. In contrast to the IOS results, the fMRI-derived simulated activation maps significantly differentiated all three odors' evoked activation maps. Importantly, no evoked activation maps could be significantly differentiated using simulated activation maps produced using phase-randomized resting-state networks. Given that some highly organized resting-state networks did not correlate with any odors' evoked activation maps, we posit that these resting-state networks may characterize evoked activation maps associated with odors not studied. These results emphasize that fluctuations in spontaneous activity form a foundation for active processing, signifying the relevance of resting-state mapping to functional neuroimaging.

Macroscale variation in resting-state neuronal activity and connectivity assessed by simultaneous calcium imaging, hemodynamic imaging and electrophysiology

Functional imaging of spontaneous activity continues to play an important role in the field of connectomics. The most common imaging signal used for these experiments is the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal, but how this signal relates to spontaneous neuronal activity remains incompletely understood. Genetically encoded calcium indicators represent a promising tool to study this problem, as they can provide brain-wide measurements of neuronal activity compared to point measurements afforded by electrophysiological recordings. However, the relationship between the calcium signal and neurophysiological parameters at the mesoscopic scale requires further systematic characterization. Therefore, we collected simultaneous resting-state measurements of electrophysiology, along with calcium and hemodynamic imaging, in lightly anesthetized mice to investigate two aims. First, we examined the relationship between each imaging signal and the simultaneously recorded electrophysiological signal in a single brain region, finding that both signals are better correlated with multi-unit activity compared to local field potentials, with the calcium signal possessing greater signal-to-noise ratio and regional specificity. Second, we used the resting-state imaging data to model the relationship between the calcium and hemodynamic signals across the brain. We found that this relationship varied across brain regions in a way that is consistent across animals, with delays increasing by600 ms towards posterior cortical regions. Furthermore, while overall functional connectivity (FC) measured by the hemodynamic signal is significantly correlated with FC measured by calcium, the two estimates were found to be significantly different. We hypothesize that these differences arise at least in part from the observed regional variation in the hemodynamic response. In total, this work highlights some of the caveats needed in interpreting hemodynamic-based measurements of FC, as well as the need for improved modeling methods to reduce this potential source of bias.

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.

How networks communicate: propagation patterns in spontaneous brain activity

Initially regarded as 'noise', spontaneous (intrinsic) activity accounts for a large portion of the brain's metabolic cost. Moreover, it is now widely known that infra-slow (less than 0.1 Hz) spontaneous activity, measured using resting state functional magnetic resonance imaging of the blood oxygen level-dependent (BOLD) signal, is correlated within functionally defined resting state networks (RSNs). However, despite these advances, the temporal organization of spontaneous BOLD fluctuations has remained elusive. By studying temporal lags in the resting state BOLD signal, we have recently shown that spontaneous BOLD fluctuations consist of remarkably reproducible patterns of whole brain propagation. Embedded in these propagation patterns are unidirectional 'motifs' which, in turn, give rise to RSNs. Additionally, propagation patterns are markedly altered as a function of state, whether physiological or pathological. Understanding such propagation patterns will likely yield deeper insights into the role of spontaneous activity in brain function in health and disease.This article is part of the themed issue 'Interpreting blood oxygen level-dependent: a dialogue between cognitive and cellular neuroscience'.

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.

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.

Graph Theoretic Analysis of Resting State Functional MR Imaging

Graph theoretic analyses applied to examine the brain at rest have played a critical role in clarifying the foundations of the brain's intrinsic and task-related activity. There are many opportunities for clinical scientists to describe and predict dysfunction using a network perspective. This primer describes the theoretic basis and practical application of graph theoretic analysis to resting state functional MR imaging data. Major practices, concepts, and findings are concisely reviewed. The theoretic and practical frontiers of resting state functional MR imaging are highlighted with observations about major avenues for conceptual advances and clinical translation.

Brain functional connectivity alterations in a rat model of excessive alcohol drinking: A resting-state network analysis

Alcohol use disorders (AUD) are a major public health concern. Understanding the brain network alterations is of the utmost importance to diagnose and develop treatment strategies. Employing resting-state functional magnetic resonance imaging, we have performed a longitudinal study in a rat model of chronic excessive alcohol consumption, to identify functional alterations in brain networks triggered by alcohol drinking. Two time points were considered: 1) before alcohol consumption (control condition) and 2) after 30 days of alcohol drinking (alcohol condition). We first identified nine resting-state networks with group independent component analysis. Afterwards, dual regression was applied to obtain subject-specific time courses and spatial maps. L2-regularized partial correlation analysis between pairs of networks showed that functional connectivity (FC) between the retrosplenial-visual and striatal networks decreases due to alcohol consumption, whereas FC between the prefrontal-cingulate and striatal networks increases. Analysis of subject-specific spatial maps revealed FC decreases within networks after alcohol drinking, including the striatal, motor-parietal, prefrontal-cingulate, retrosplenial-visual and left motor-parietal networks. Overall, our results unveil a generalized decrease in brain FC induced by alcohol drinking in genetically predisposed animals, even after a relatively short period of exposure (1 month). The only exception to this hypo-connectivity state is the functional association between the striatal and prefrontal-cingulate networks, which increases after drinking, supporting evidence in human alcoholics.

Concurrent tACS-fMRI Reveals Causal Influence of Power Synchronized Neural Activity on Resting State fMRI Connectivity

Resting state fMRI (rs-fMRI) is commonly used to study the brain's intrinsic neural coupling, which reveals specific spatiotemporal patterns in the form of resting state networks (RSNs). It has been hypothesized that slow rs-fMRI oscillations (<0.1 Hz) are driven by underlying electrophysiological rhythms that typically occur at much faster timescales (>5 Hz); however, causal evidence for this relationship is currently lacking. Here we measured rs-fMRI in humans while applying transcranial alternating current stimulation (tACS) to entrain brain rhythms in left and right sensorimotor cortices. The two driving tACS signals were tailored to the individual's α rhythm (8-12 Hz) and fluctuated in amplitude according to a 1 Hz power envelope. We entrained the left versus right hemisphere in accordance to two different coupling modes where either α oscillations were synchronized between hemispheres (phase-synchronized tACS) or the slower oscillating power envelopes (power-synchronized tACS). Power-synchronized tACS significantly increased rs-fMRI connectivity within the stimulated RSN compared with phase-synchronized or no tACS. This effect outlasted the stimulation period and tended to be more effective in individuals who exhibited a naturally weak interhemispheric coupling. Using this novel approach, our data provide causal evidence that synchronized power fluctuations contribute to the formation of fMRI-based RSNs. Moreover, our findings demonstrate that the brain's intrinsic coupling at rest can be selectively modulated by choosing appropriate tACS signals, which could lead to new interventions for patients with altered rs-fMRI connectivity.SIGNIFICANCE STATEMENT Resting state fMRI (rs-fMRI) has become an important tool to estimate brain connectivity. However, relatively little is known about how slow hemodynamic oscillations measured with fMRI relate to electrophysiological processes. It was suggested that slowly fluctuating power envelopes of electrophysiological signals synchronize across brain areas and that the topography of this activity is spatially correlated to resting state networks derived from rs-fMRI. Here we take a novel approach to address this problem and establish a causal link between the power fluctuations of electrophysiological signals and rs-fMRI via a new neuromodulation paradigm, which exploits these power synchronization mechanisms. These novel mechanistic insights bridge different scientific domains and are of broad interest to researchers in the fields of Medical Imaging, Neuroscience, Physiology, and Psychology.

Altered dynamics of brain connectivity in major depressive disorder at-rest and during task performance

Major depressive disorder (MDD) has been associated with alterations in several functional brain networks. Previous studies investigating brain networks in MDD during the performance of a task have yielded inconsistent results with the function of the brain at rest. In this study, we used functional magnetic resonance imaging at rest and during a goal-directed task to investigate dynamics of functional connectivity in 19 unmedicated patients with MDD and 19 healthy controls across both experimental paradigms. Patients had spatial differences in the default mode network (DMN), in the executive network (EN), and in the dorsal attention network (DAN) compared to controls at rest and during task performance. In patients the amplitude of the low frequency (LFO) oscillations was reduced in the motor and in the DAN networks during both paradigms. There was a diagnosis by paradigm interaction on the LFOs amplitude of the salience network, with increased amplitude change between task and rest in patients relative to controls. Our findings suggest that the function of several networks could be intrinsically affected in MDD and this could be viable phenotype for the investigation on the neurobiological mechanisms of this disorder and its treatment.

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.

A new neuroinformatics approach to personalized medicine in neurology: The Virtual Brain

PURPOSE OF REVIEW

An exciting advance in the field of neuroimaging is the acquisition and processing of very large data sets (so called 'big data'), permitting large-scale inferences that foster a greater understanding of brain function in health and disease. Yet what we are clearly lacking are quantitative integrative tools to translate this understanding to the individual level to lay the basis for personalized medicine.

RECENT FINDINGS

Here we address this challenge through a review on how the relatively new field of neuroinformatics modeling has the capacity to track brain network function at different levels of inquiry, from microscopic to macroscopic and from the localized to the distributed. In this context, we introduce a new and unique multiscale approach, The Virtual Brain (TVB), that effectively models individualized brain activity, linking large-scale (macroscopic) brain dynamics with biophysical parameters at the microscopic level. We also show how TVB modeling provides unique biological interpretable data in epilepsy and stroke.

SUMMARY

These results establish the basis for a deliberate integration of computational biology and neuroscience into clinical approaches for elucidating cellular mechanisms of disease. In the future, this can provide the means to create a collection of disease-specific models that can be applied on the individual level to personalize therapeutic interventions. VIDEO ABSTRACT.

A Hitchhiker's Guide to Functional Magnetic Resonance Imaging

Functional Magnetic Resonance Imaging (fMRI) studies have become increasingly popular both with clinicians and researchers as they are capable of providing unique insights into brain functions. However, multiple technical considerations (ranging from specifics of paradigm design to imaging artifacts, complex protocol definition, and multitude of processing and methods of analysis, as well as intrinsic methodological limitations) must be considered and addressed in order to optimize fMRI analysis and to arrive at the most accurate and grounded interpretation of the data. In practice, the researcher/clinician must choose, from many available options, the most suitable software tool for each stage of the fMRI analysis pipeline. Herein we provide a straightforward guide designed to address, for each of the major stages, the techniques, and tools involved in the process. We have developed this guide both to help those new to the technique to overcome the most critical difficulties in its use, as well as to serve as a resource for the neuroimaging community.

Scale-Free Neural and Physiological Dynamics in Naturalistic Stimuli Processing

Neural activity recorded at multiple spatiotemporal scales is dominated by arrhythmic fluctuations without a characteristic temporal periodicity. Such activity often exhibits a 1/f-type power spectrum, in which power falls off with increasing frequency following a power-law function: [Formula: see text], which is indicative of scale-free dynamics. Two extensively studied forms of scale-free neural dynamics in the human brain are slow cortical potentials (SCPs)-the low-frequency (<5 Hz) component of brain field potentials-and the amplitude fluctuations of α oscillations, both of which have been shown to carry important functional roles. In addition, scale-free dynamics characterize normal human physiology such as heartbeat dynamics. However, the exact relationships among these scale-free neural and physiological dynamics remain unclear. We recorded simultaneous magnetoencephalography and electrocardiography in healthy subjects in the resting state and while performing a discrimination task on scale-free dynamical auditory stimuli that followed different scale-free statistics. We observed that long-range temporal correlation (captured by the power-law exponent β) in SCPs positively correlated with that of heartbeat dynamics across time within an individual and negatively correlated with that of α-amplitude fluctuations across individuals. In addition, across individuals, long-range temporal correlation of both SCP and α-oscillation amplitude predicted subjects' discrimination performance in the auditory task, albeit through antagonistic relationships. These findings reveal interrelations among different scale-free neural and physiological dynamics and initial evidence for the involvement of scale-free neural dynamics in the processing of natural stimuli, which often exhibit scale-free dynamics.

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.

Physiological characterization of a robust survival rodent fMRI method

Anesthetics are commonly used in preclinical functional MRI studies. It is well-appreciated that proper choice of anesthetics is of critical importance for maintaining a physiologically normal range of autonomic functioning. A recent study, using a low dose of dexmedetomidine (active isomer of medetomidine) in combination with a low dose of isoflurane, suggested stable measurements across repeated fMRI experiments in individual animals with each session lasting up to several hours. The rat default mode network has been successfully identified using this preparation, indicating that this protocol minimally disturbs brain network functions. However, medetomidine is known to cause peripheral vasoconstriction, respiratory suppression, and bradycardia, each of which could independently confound the BOLD signal. The goal of this study was to systematically characterize physiological conditions for fMRI experiments under this anesthetic regimen. To this end, we acquired somatosensory stimulation "task-evoked" and resting-state fMRI to evaluate the integrity of neurovascular coupling and brain network function during three time windows (0-30min, 30-90min, and 90-150min) following dexmedetomidine initiation. Results demonstrate that both evoked BOLD response and resting-state fMRI signal remained stable during the 90-150min time window, while autonomic physiological parameters maintained near-normal conditions during this period. Our data suggest that using a spontaneously-inhaled, low dose of isoflurane in combination with a continuous low dose of dexmedetomidine is a viable option for longitudinal imaging studies in rats.

Correlated inter-regional variations in low frequency local field potentials and resting state BOLD signals within S1 cortex of monkeys

The hypothesis that specific frequency components of the spontaneous local field potentials (LFPs) underlie low frequency fluctuations of resting state fMRI (rsfMRI) signals was tested. The previous analyses of rsfMRI signals revealed differential inter-regional correlations among areas 3a, 3b, and 1 of primary somatosensory cortex (S1) in anesthetized monkeys (Wang et al. [2013]: Neuron 78:1116-1126). Here LFP band(s) which correlated between S1 regions, and how these inter-regional correlation differences covaried with rsfMRI signals were examined. LFP signals were filtered into seven bands (delta, theta, alpha, beta, gamma low, gamma high, and gamma very high), and then a Hilbert transformation was applied to obtain measures of instantaneous amplitudes and temporal lags between regions of interest (ROI) digit-digit pairs (areas 3b-area 1, area 3a-area 1, area 3a-area 3b) and digit-face pairs (area 3b-face, area 1-face, and area 3a-face). It was found that variations in the inter-regional correlation strengths between digit-digit and digit-face pairs in the delta (1-4 Hz), alpha (9-14 Hz), beta (15-30 Hz), and gamma (31-50 Hz) bands parallel those of rsfMRI signals to varying degrees. Temporal lags between digit-digit area pairs varied across LFP bands, with area 3a mostly leading areas 1/2 and 3b. In summary, the data demonstrates that the low and middle frequency range (1-50 Hz) of spontaneous LFP signals similarly covary with the low frequency fluctuations of rsfMRI signals within local circuits of S1, supporting a neuronal electrophysiological basis of rsfMRI signals. Inter-areal LFP temporal lag differences provided novel insights into the directionality of information flow among S1 areas at rest. Hum Brain Mapp 37:2755-2766, 2016. © 2016 Wiley Periodicals, Inc.

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.

Transient neuronal coactivations embedded in globally propagating waves underlie resting-state functional connectivity

Resting-state functional connectivity (FC), which measures the correlation of spontaneous hemodynamic signals (HemoS) between brain areas, is widely used to study brain networks noninvasively. It is commonly assumed that spatial patterns of HemoS-based FC (Hemo-FC) reflect large-scale dynamics of underlying neuronal activity. To date, studies of spontaneous neuronal activity cataloged heterogeneous types of events ranging from waves of activity spanning the entire neocortex to flash-like activations of a set of anatomically connected cortical areas. However, it remains unclear how these various types of large-scale dynamics are interrelated. More importantly, whether each type of large-scale dynamics contributes to Hemo-FC has not been explored. Here, we addressed these questions by simultaneously monitoring neuronal calcium signals (CaS) and HemoS in the entire neocortex of mice at high spatiotemporal resolution. We found a significant relationship between two seemingly different types of large-scale spontaneous neuronal activity-namely, global waves propagating across the neocortex and transient coactivations among cortical areas sharing high FC. Different sets of cortical areas, sharing high FC within each set, were coactivated at different timings of the propagating global waves, suggesting that spatial information of cortical network characterized by FC was embedded in the phase of the global waves. Furthermore, we confirmed that such transient coactivations in CaS were indeed converted into spatially similar coactivations in HemoS and were necessary to sustain the spatial structure of Hemo-FC. These results explain how global waves of spontaneous neuronal activity propagating across large-scale cortical network contribute to Hemo-FC in the resting state.

Hemodynamic signals in fNIRS

Near-infrared spectroscopy (NIRS) was originally designed for clinical monitoring of tissue oxygenation, and it has also been developed into a useful tool in neuroimaging studies, with the so-called functional NIRS (fNIRS). With NIRS, cerebral activation is detected by measuring the cerebral hemoglobin (Hb), where however, the precise correlation between NIRS signal and neural activity remains to be fully understood. This can in part be attributed to the situation that NIRS signals are inherently subject to contamination by signals arising from extracerebral tissue. In recent years, several approaches have been investigated to distinguish between NIRS signals originating in cerebral tissue and signals originating in extracerebral tissue. Selective measurements of cerebral Hb will enable a further evolution of fNIRS. This chapter is divided into six sections: first a summary of the basic theory of NIRS, NIRS signals arising in the activated areas, correlations between NIRS signals and fMRI signals, correlations between NIRS signals and neural activities, and the influence of a variety of extracerebral tissue on NIRS signals and approaches to this issue are reviewed. Finally, future prospects of fNIRS are described.

N-methyl-D-aspartate receptor encephalitis mediates loss of intrinsic activity measured by functional MRI

Spontaneous brain activity is required for the development and maintenance of normal brain function. Many disease processes disrupt the organization of intrinsic brain activity, but few pervasively reduce the amplitude of resting state blood oxygen level dependent (BOLD) fMRI fluctuations. We report the case of a female with anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis, longitudinally studied during the course of her illness to determine the contribution of NMDAR signaling to spontaneous brain activity. Resting state BOLD fMRI was measured at the height of her illness and 18 weeks following discharge from hospital. Conventional resting state networks were defined using established methods. Correlation and covariance matrices were calculated by extracting the BOLD time series from regions of interest and calculating either the correlation or covariance quantity. The intrinsic activity was compared between visits, and to expected activity from 45 similarly aged healthy individuals. Near the height of the illness, the patient exhibited profound loss of consciousness, high-amplitude slowing of the electroencephalogram, and a severe reduction in the amplitude of spontaneous BOLD fMRI fluctuations. The patient's neurological status and measures of intrinsic activity improved following treatment. We conclude that NMDAR-mediated signaling plays a critical role in the mechanisms that give rise to organized spontaneous brain activity. Loss of intrinsic activity is associated with profound disruptions of consciousness and cognition.

Evaluation of sliding window correlation performance for characterizing dynamic functional connectivity and brain states

A promising recent development in the study of brain function is the dynamic analysis of resting-state functional MRI scans, which can enhance understanding of normal cognition and alterations that result from brain disorders. One widely used method of capturing the dynamics of functional connectivity is sliding window correlation (SWC). However, in the absence of a "gold standard" for comparison, evaluating the performance of the SWC in typical resting-state data is challenging. This study uses simulated networks (SNs) with known transitions to examine the effects of parameters such as window length, window offset, window type, noise, filtering, and sampling rate on the SWC performance. The SWC time course was calculated for all node pairs of each SN and then clustered using the k-means algorithm to determine how resulting brain states match known configurations and transitions in the SNs. The outcomes show that the detection of state transitions and durations in the SWC is most strongly influenced by the window length and offset, followed by noise and filtering parameters. The effect of the image sampling rate was relatively insignificant. Tapered windows provide less sensitivity to state transitions than rectangular windows, which could be the result of the sharp transitions in the SNs. Overall, the SWC gave poor estimates of correlation for each brain state. Clustering based on the SWC time course did not reliably reflect the underlying state transitions unless the window length was comparable to the state duration, highlighting the need for new adaptive window analysis techniques.

Disconnection and hyper-connectivity underlie reorganization after TBI: A rodent functional connectomic analysis

While past neuroimaging methods have contributed greatly to our understanding of brain function after traumatic brain injury (TBI), resting state functional MRI (rsfMRI) connectivity methods have more recently provided a far more unbiased approach with which to monitor brain circuitry compared to task-based approaches. However, current knowledge on the physiologic underpinnings of the correlated blood oxygen level dependent signal, and how changes in functional connectivity relate to reorganizational processes that occur following injury is limited. The degree and extent of this relationship remain to be determined in order that rsfMRI methods can be fully adapted for determining the optimal timing and type of rehabilitative interventions that can be used post-TBI to achieve the best outcome. Very few rsfMRI studies exist after experimental TBI and therefore we chose to acquire rsfMRI data before and at 7, 14 and 28 days after experimental TBI using a well-known, clinically-relevant, unilateral controlled cortical impact injury (CCI) adult rat model of TBI. This model was chosen since it has widespread axonal injury, a well-defined time-course of reorganization including spine, dendrite, axonal and cortical map changes, as well as spontaneous recovery of sensorimotor function by 28 d post-injury from which to interpret alterations in functional connectivity. Data were co-registered to a parcellated rat template to generate adjacency matrices for network analysis by graph theory. Making no assumptions about direction of change, we used two-tailed statistical analysis over multiple brain regions in a data-driven approach to access global and regional changes in network topology in order to assess brain connectivity in an unbiased way. Our main hypothesis was that deficits in functional connectivity would become apparent in regions known to be structurally altered or deficient in axonal connectivity in this model. The data show the loss of functional connectivity predicted by the structural deficits, not only within the primary sensorimotor injury site and pericontused regions, but the normally connected homotopic cortex, as well as subcortical regions, all of which persisted chronically. Especially novel in this study is the unanticipated finding of widespread increases in connection strength that dwarf both the degree and extent of the functional disconnections, and which persist chronically in some sensorimotor and subcortically connected regions. Exploratory global network analysis showed changes in network parameters indicative of possible acutely increased random connectivity and temporary reductions in modularity that were matched by local increases in connectedness and increased efficiency among more weakly connected regions. The global network parameters: shortest path-length, clustering coefficient and modularity that were most affected by trauma also scaled with the severity of injury, so that the corresponding regional measures were correlated to the injury severity most notably at 7 and 14 days and especially within, but not limited to, the contralateral cortex. These changes in functional network parameters are discussed in relation to the known time-course of physiologic and anatomic data that underlie structural and functional reorganization in this experiment model of TBI.

Spontaneous Neural Dynamics and Multi-scale Network Organization

Spontaneous neural activity has historically been viewed as task-irrelevant noise that should be controlled for via experimental design, and removed through data analysis. However, electrophysiology and functional MRI studies of spontaneous activity patterns, which have greatly increased in number over the past decade, have revealed a close correspondence between these intrinsic patterns and the structural network architecture of functional brain circuits. In particular, by analyzing the large-scale covariation of spontaneous hemodynamics, researchers are able to reliably identify functional networks in the human brain. Subsequent work has sought to identify the corresponding neural signatures via electrophysiological measurements, as this would elucidate the neural origin of spontaneous hemodynamics and would reveal the temporal dynamics of these processes across slower and faster timescales. Here we survey common approaches to quantifying spontaneous neural activity, reviewing their empirical success, and their correspondence with the findings of neuroimaging. We emphasize invasive electrophysiological measurements, which are amenable to amplitude- and phase-based analyses, and which can report variations in connectivity with high spatiotemporal precision. After summarizing key findings from the human brain, we survey work in animal models that display similar multi-scale properties. We highlight that, across many spatiotemporal scales, the covariance structure of spontaneous neural activity reflects structural properties of neural networks and dynamically tracks their functional repertoire.

Low- but Not High-Frequency LFP Correlates with Spontaneous BOLD Fluctuations in Rat Whisker Barrel Cortex

Resting-state magnetic resonance imaging (rsMRI) is thought to reflect ongoing spontaneous brain activity. However, the precise neurophysiological basis of rsMRI signal remains elusive. Converging evidence supports the notion that local field potential (LFP) signal in the high-frequency range correlates with fMRI response evoked by a task (e.g., visual stimulation). It remains uncertain whether this relationship extends to rsMRI. In this study, we systematically modulated LFP signal in the whisker barrel cortex (WBC) by unilateral deflection of rat whiskers. Results show that functional connectivity between bilateral WBC was significantly modulated at the 2 Hz, but not at the 4 or 6 Hz, stimulus condition. Electrophysiologically, only in the low-frequency range (<5 Hz) was the LFP power synchrony in bilateral WBC significantly modulated at 2 Hz, but not at 4- or 6-Hz whisker stimulation, thus distinguishing these 2 experimental conditions, and paralleling the findings in rsMRI. LFP power synchrony in other frequency ranges was modulated in a way that was neither unique to the specific stimulus conditions nor parallel to the fMRI results. Our results support the hypothesis that emphasizes the role of low-frequency LFP signal underlying rsMRI.

Frequency-dependent changes in amplitude of low-frequency oscillations in depression: A resting-state fMRI study

OBJECTIVE

We conducted this fMRI study to examine whether the alterations in amplitudes of low-frequency oscillation (LFO) of major depressive disorder (MDD) patients were frequency dependent.

MATERIALS AND METHODS

The LFO amplitudes (as indexed by amplitude of low-frequency fluctuation [ALFF] and fractional ALFF [fALFF]) within 4 narrowly-defined frequency bands (slow-5: 0.01-0.027Hz, slow-4: 0.027-0.073Hz, slow-3: 0.073-0.198Hz, and slow-2: 0.198-0.25Hz) were computed using resting-state fMRI data of 35 MDD patients and 32 healthy subjects. Repeated-measures analysis of variance (ANOVA) was performed on ALFF and fALFF both within the low frequency bands of slow-4 and slow-5 and within all of the four bands.

RESULTS

We observed significant main effects of group and frequency on ALFF and fALFF in widely distributed brain regions. Importantly, significant group and frequency interaction effects were observed in the ventromedial prefrontal cortex, inferior frontal gyrus, precentral gyrus, in a left-sided fashion, the bilateral posterior cingulate and precuneus, during ANOVA both within slow-4 and slow-5 bands and within all the frequency bands.

CONCLUSIONS

The results suggest that the alterations of LFO amplitudes in specific brain regions in MDD patients could be more sensitively detected in the slow-5 rather than the slow-4 bands. The findings may provide guidance for the frequency choice of future resting-state fMRI studies of MDD.

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

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

Laminar Distribution of Phase-Amplitude Coupling of Spontaneous Current Sources and Sinks

Although resting-state functional connectivity is a commonly used neuroimaging paradigm, the underlying mechanisms remain unknown. Thalamo-cortical and cortico-cortical circuits generate oscillations at different frequencies during spontaneous activity. However, it remains unclear how the various rhythms interact and whether their interactions are lamina-specific. Here we investigated intra- and inter-laminar spontaneous phase-amplitude coupling (PAC). We recorded local-field potentials using laminar probes inserted in the forelimb representation of rat area S1. We then computed time-series of frequency-band- and lamina-specific current source density (CSD), and PACs of CSD for all possible pairs of the classical frequency bands in the range of 1–150 Hz. We observed both intra- and inter-laminar spontaneous PAC. Of 18 possible combinations, 12 showed PAC, with the highest measures of interaction obtained for the pairs of the theta/gamma and delta/gamma bands. Intra- and inter-laminar PACs involving layers 2/3–5a were higher than those involving layer 6. Current sinks (sources) in the delta band were associated with increased (decreased) amplitudes of high-frequency signals in the beta to fast gamma bands throughout layers 2/3–6. Spontaneous sinks (sources) of the theta and alpha bands in layers 2/3–4 were on average linked to dipoles completed by sources (sinks) in layer 6, associated with high (low) amplitudes of the beta to fast-gamma bands in the entire cortical column. Our findings show that during spontaneous activity, delta, theta, and alpha oscillations are associated with periodic excitability, which for the theta and alpha bands is lamina-dependent. They further emphasize the differences between the function of layer 6 and that of the superficial layers, and the role of layer 6 in controlling activity in those layers. Our study links theories on the involvement of PAC in resting-state functional connectivity with previous work that revealed lamina-specific anatomical thalamo-cortico-cortical connections.

On the origins of the 'global signal' determined by functional magnetic resonance imaging in the resting state

OBJECTIVE

Functional magnetic resonance imaging blood oxygen level dependent (BOLD) determinations of correlations between 'resting-state' neuronal activity in different regions of cortex have generated much interest. Determination of these correlations requires regressing out signals that are correlated in all parts of the cortex and are taken to be artefactual, such as those due to movement, respiration and cardiovascular activity. However when these are removed there still remains a 'global signal' (GS), which is taken to be of unknown physiological origin, and is regressed out by some researchers but not by others.

APPROACH

We have investigated the origin of this GS using cortical models consisting of coupled networks of modules representing regions of interest.

MAIN RESULTS

We show that the GS has an amplitude that is linearly related to the average correlation between the modules/voxels in the network over a large range of such correlations. The GS arises as a consequence of feedback between the modules/voxels leading to correlations in their BOLD signals. Given the relationship between the GS and the average correlations it might be anticipated that regressing out the GS during preprocessing will significantly modify the correlations subsequently determined. This is shown to be the case when comparing the connections of individual modules with that predicted by the correlations.

SIGNIFICANCE

The present model shows that such correlations can arise as a consequence of the intermodular feedback connectivity without recourse to imposing a GS independent of the connectivity. Our model indicates that the GS reflects the extent of feedback pathways provided by the intermodular/inter-regional connections and hence the average correlation between modules or regions of cortex. However the model has not been used to elucidate the possible contributions of a GS independent of the connectivity, which might indeed contribute to the GS of the cortex.

Cortical Network Models of Firing Rates in the Resting and Active States Predict BOLD Responses

Measurements of blood oxygenation level dependent (BOLD) signals have produced some surprising observations. One is that their amplitude is proportional to the entire activity in a region of interest and not just the fluctuations in this activity. Another is that during sleep and anesthesia the average BOLD correlations between regions of interest decline as the activity declines. Mechanistic explanations of these phenomena are described here using a cortical network model consisting of modules with excitatory and inhibitory neurons, taken as regions of cortical interest, each receiving excitatory inputs from outside the network, taken as subcortical driving inputs in addition to extrinsic (intermodular) connections, such as provided by associational fibers. The model shows that the standard deviation of the firing rate is proportional to the mean frequency of the firing when the extrinsic connections are decreased, so that the mean BOLD signal is proportional to both as is observed experimentally. The model also shows that if these extrinsic connections are decreased or the frequency of firing reaching the network from the subcortical driving inputs is decreased, or both decline, there is a decrease in the mean firing rate in the modules accompanied by decreases in the mean BOLD correlations between the modules, consistent with the observed changes during NREM sleep and under anesthesia. Finally, the model explains why a transient increase in the BOLD signal in a cortical area, due to a transient subcortical input, gives rises to responses throughout the cortex as observed, with these responses mediated by the extrinsic (intermodular) connections.

Steady-state BOLD Response to Higher-order Cognition Modulates Low-Frequency Neural Oscillations

Steady-state responses (SSRs) reflect the synchronous neural oscillations evoked by noninvasive and consistently repeated stimuli at the fundamental or harmonic frequencies. The steady-state evoked potentials (SSEPs; the representative form of the SSRs) have been widely used in the cognitive and clinical neurosciences and brain–computer interface research. However, the steady-state evoked potentials have limitations in examining high-frequency neural oscillations and basic cognition. In addition, synchronous neural oscillations in the low frequency range (&lt;1 Hz) and in higher-order cognition have received a little attention. Therefore, we examined the SSRs in the low frequency range using a new index, the steady-state BOLD responses (SSBRs) evoked by semantic stimuli. Our results revealed that the significant SSBRs were induced at the fundamental frequency of stimuli and the first harmonic in task-related regions, suggesting the enhanced variability of neural oscillations entrained by exogenous stimuli. The SSBRs were independent of neurovascular coupling and characterized by sensorimotor bias, an indication of regional-dependent neuroplasticity. Furthermore, the amplitude of SSBRs may predict behavioral performance and show the psychophysiological relevance. Our findings provide valuable insights into the understanding of the SSRs evoked by higher-order cognition and how the SSRs modulate low-frequency neural oscillations.

Propagated infra-slow intrinsic brain activity reorganizes across wake and slow wave sleep

Propagation of slow intrinsic brain activity has been widely observed in electrophysiogical studies of slow wave sleep (SWS). However, in human resting state fMRI (rs-fMRI), intrinsic activity has been understood predominantly in terms of zero-lag temporal synchrony (functional connectivity) within systems known as resting state networks (RSNs). Prior rs-fMRI studies have found that RSNs are generally preserved across wake and sleep. Here, we use a recently developed analysis technique to study propagation of infra-slow intrinsic blood oxygen level dependent (BOLD) signals in normal adults during wake and SWS. This analysis reveals marked changes in propagation patterns in SWS vs. wake. Broadly, ordered propagation is preserved within traditionally defined RSNs but lost between RSNs. Additionally, propagation between cerebral cortex and subcortical structures reverses directions, and intra-cortical propagation becomes reorganized, especially in visual and sensorimotor cortices. These findings show that propagated rs-fMRI activity informs theoretical accounts of the neural functions of sleep.

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.

The power of using functional fMRI on small rodents to study brain pharmacology and disease

Functional magnetic resonance imaging (fMRI) is an excellent tool to study the effect of pharmacological modulations on brain function in a non-invasive and longitudinal manner. We introduce several blood oxygenation level dependent (BOLD) fMRI techniques, including resting state (rsfMRI), stimulus-evoked (st-fMRI), and pharmacological MRI (phMRI). Respectively, these techniques permit the assessment of functional connectivity during rest as well as brain activation triggered by sensory stimulation and/or a pharmacological challenge. The first part of this review describes the physiological basis of BOLD fMRI and the hemodynamic response on which the MRI contrast is based. Specific emphasis goes to possible effects of anesthesia and the animal’s physiological conditions on neural activity and the hemodynamic response. The second part of this review describes applications of the aforementioned techniques in pharmacologically induced, as well as in traumatic and transgenic disease models and illustrates how multiple fMRI methods can be applied successfully to evaluate different aspects of a specific disorder. For example, fMRI techniques can be used to pinpoint the neural substrate of a disease beyond previously defined hypothesis-driven regions-of-interest. In addition, fMRI techniques allow one to dissect how specific modifications (e.g., treatment, lesion etc.) modulate the functioning of specific brain areas (st-fMRI, phMRI) and how functional connectivity (rsfMRI) between several brain regions is affected, both in acute and extended time frames. Furthermore, fMRI techniques can be used to assess/explore the efficacy of novel treatments in depth, both in fundamental research as well as in preclinical settings. In conclusion, by describing several exemplary studies, we aim to highlight the advantages of functional MRI in exploring the acute and long-term effects of pharmacological substances and/or pathology on brain functioning along with several methodological considerations.

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.

Mesoscale infraslow spontaneous membrane potential fluctuations recapitulate high-frequency activity cortical motifs

Neuroimaging of spontaneous, resting-state infraslow (<0.1 Hz) brain activity has been used to reveal the regional functional organization of the brain and may lead to the identification of novel biomarkers of neurological disease. However, these imaging studies generally rely on indirect measures of neuronal activity and the nature of the neuronal activity correlate remains unclear. Here we show, using wide-field, voltage-sensitive dye imaging, the mesoscale spatiotemporal structure and pharmacological dependence of spontaneous, infraslow cortical activity in anaesthetized and awake mice. Spontaneous infraslow activity is regionally distinct, correlates with electroencephalography and local field potential recordings, and shows bilateral symmetry between cortical hemispheres. Infraslow activity is attenuated and its functional structure abolished after treatment with voltage-gated sodium channel and glutamate receptor antagonists. Correlation analysis reveals patterns of infraslow regional connectivity that are analogous to cortical motifs observed from higher-frequency spontaneous activity and reflect the underlying framework of intracortical axonal projections.