Electrophysiological low-frequency coherence and cross-frequency coupling contribute to BOLD connectivity

Intracranial Electrophysiology of the Human Default Network

The human default network (DN) plays a critical role in internally directed cognition, behavior, and neuropsychiatric disease. Despite much progress with functional neuroimaging, persistent questions still linger concerning the electrophysiological underpinnings, fast temporal dynamics, and causal importance of the DN. Here, we review how direct intracranial recording and stimulation of the DN provides a unique combination of high spatiotemporal resolution and causal information that speaks directly to many of these outstanding questions. Our synthesis highlights the electrophysiological basis of activation, suppression, and connectivity of the DN, each key areas of debate in the literature. Integrating these unique electrophysiological data with extant neuroimaging findings will help lay the foundation for a mechanistic account of DN function in human behavior and cognition.

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.

Maturation trajectories of cortical resting-state networks depend on the mediating frequency band

The functional significance of resting state networks and their abnormal manifestations in psychiatric disorders are firmly established, as is the importance of the cortical rhythms in mediating these networks. Resting state networks are known to undergo substantial reorganization from childhood to adulthood, but whether distinct cortical rhythms, which are generated by separable neural mechanisms and are often manifested abnormally in psychiatric conditions, mediate maturation differentially, remains unknown. Using magnetoencephalography (MEG) to map frequency band specific maturation of resting state networks from age 7 to 29 in 162 participants (31 independent), we found significant changes with age in networks mediated by the beta (13-30 Hz) and gamma (31-80 Hz) bands. More specifically, gamma band mediated networks followed an expected asymptotic trajectory, but beta band mediated networks followed a linear trajectory. Network integration increased with age in gamma band mediated networks, while local segregation increased with age in beta band mediated networks. Spatially, the hubs that changed in importance with age in the beta band mediated networks had relatively little overlap with those that showed the greatest changes in the gamma band mediated networks. These findings are relevant for our understanding of the neural mechanisms of cortical maturation, in both typical and atypical development.

Brain Resting-State Functional Connectivity Is Preserved Under Sevoflurane Anesthesia in Patients with Pervasive Developmental Disorders: A Pilot Study

Functional connectivity studies play a huge role in understanding the relationship between the network connections and the behavioral phenotype of patients with pervasive developmental disorders (PDD). Some patients with PDD may not be able to tolerate the imaging procedure while they are awake, and, hence, they often need general anesthesia. General anesthesia is a confounding factor in functional imaging studies due to its effect on the functional connectivity. The objective of this study is to look at the resting-state functional connectivity (RS-FC) under sevoflurane anesthesia in patients with PDDs. Thirteen adults with PDD scheduled for magnetic resonance imaging (MRI) of the brain under general anesthesia were recruited for the study. Resting-state functional MRI (fMRI) scans were acquired at 1 minimum alveolar concentration (MAC) of sevoflurane. Spontaneous blood oxygenation level-dependent fluctuations were measured, and a seed-voxel analysis was done to identify the resting-state networks. Subjects' data were compared with data from 16 nonanesthetized healthy controls. Six networks (default mode network [DMN], executive control network [ECN], salience network [SN], auditory, visual, and sensorimotor) were investigated. At 1 MAC sevoflurane anesthesia, RS-FC was preserved in all the networks. Secondary analysis of connectivity showed a decrease in connectivity within the thalamus and an increase in DMN-ECN and DMN-SN cross-network connectivity in the anesthetized patient group compared to healthy controls. Previous reports suggested that even mild levels of anesthesia could reduce overall fluctuation levels in the major brain. However, our results provide strong evidence that most networks can sustain detectable levels of activity in patients with PDDs even under deep levels of anesthesia.

Machine Learning Based Classification of Deep Brain Stimulation Outcomes in a Rat Model of Binge Eating Using Ventral Striatal Oscillations

Neuromodulation-based interventions continue to be evaluated across an array of appetitive disorders but broader implementation of these approaches remains limited due to variable treatment outcomes. We hypothesize that individual variation in treatment outcomes may be linked to differences in the networks underlying these disorders. Here, Sprague-Dawley rats received deep brain stimulation separately within each nucleus accumbens (NAc) sub-region (core and shell) using a within-animal crossover design in a rat model of binge eating. Significant reductions in binge size were observed with stimulation of either target but with significant variation in effectiveness across individuals. When features of local field potentials (LFPs) recorded from the NAc were used to classify the pre-defined stimulation outcomes (response or non-response) from each rat using a machine-learning approach (lasso), stimulation outcomes could be classified with greater accuracy than expected by chance (effect sizes: core = 1.13, shell = 1.05). Further, these LFP features could be used to identify the best stimulation target for each animal (core vs. shell) with an effect size = 0.96. These data suggest that individual differences in underlying network activity may relate to the variable outcomes of circuit based interventions, and measures of network activity could have the potential to individually guide the selection of an optimal stimulation target to improve overall treatment response rates.

Increased thalamic resting-state connectivity as a core driver of LSD-induced hallucinations

OBJECTIVE

It has been proposed that the thalamocortical system is an important site of action of hallucinogenic drugs and an essential component of the neural correlates of consciousness. Hallucinogenic drugs such as LSD can be used to induce profoundly altered states of consciousness, and it is thus of interest to test the effects of these drugs on this system.

METHOD

100 μg LSD was administrated orally to 20 healthy participants prior to fMRI assessment. Whole brain thalamic functional connectivity was measured using ROI-to-ROI and ROI-to-voxel approaches. Correlation analyses were used to explore relationships between thalamic connectivity to regions involved in auditory and visual hallucinations and subjective ratings on auditory and visual drug effects.

RESULTS

LSD caused significant alterations in all dimensions of the 5D-ASC scale and significantly increased thalamic functional connectivity to various cortical regions. Furthermore, LSD-induced functional connectivity measures between the thalamus and the right fusiform gyrus and insula correlated significantly with subjective auditory and visual drug effects.

CONCLUSION

Hallucinogenic drug effects might be provoked by facilitations of cortical excitability via thalamocortical interactions. Our findings have implications for the understanding of the mechanism of action of hallucinogenic drugs and provide further insight into the role of the 5-HT_2A -receptor in altered states of consciousness.

A Meta-Analysis Suggests Different Neural Correlates for Implicit and Explicit Learning

A meta-analysis of non-human primates performing three different tasks (Object-Match, Category-Match, and Category-Saccade associations) revealed signatures of explicit and implicit learning. Performance improved equally following correct and error trials in the Match (explicit) tasks, but it improved more after correct trials in the Saccade (implicit) task, a signature of explicit versus implicit learning. Likewise, error-related negativity, a marker for error processing, was greater in the Match (explicit) tasks. All tasks showed an increase in alpha/beta (10-30 Hz) synchrony after correct choices. However, only the implicit task showed an increase in theta (3-7 Hz) synchrony after correct choices that decreased with learning. In contrast, in the explicit tasks, alpha/beta synchrony increased with learning and decreased thereafter. Our results suggest that explicit versus implicit learning engages different neural mechanisms that rely on different patterns of oscillatory synchrony.

Modulation of Long-Range Connectivity Patterns via Frequency-Specific Stimulation of Human Cortex

There is increasing interest in how the phase of local oscillatory activity within a brain area determines the long-range functional connectivity of that area. For example, increasing convergent evidence from a range of methodologies suggests that beta (20 Hz) oscillations may play a vital role in the function of the motor system [1, 2, 3, 4, 5]. The “communication through coherence” hypothesis posits that the precise phase of coherent oscillations in network nodes is a determinant of successful communication between them [6, 7]. Here we set out to determine whether oscillatory activity in the beta band serves to support this theory within the cortical motor network in vivo. We combined non-invasive transcranial alternating-current stimulation (tACS) [8, 9, 10, 11, 12] with resting-state functional MRI (fMRI) [13] to follow both changes in local activity and long-range connectivity, determined by inter-areal blood-oxygen-level-dependent (BOLD) signal correlation, as a proxy for communication in the human cortex. Twelve healthy subjects participated in three fMRI scans with 20 Hz, 5 Hz, or sham tACS applied separately on each scan. Transcranial magnetic stimulation (TMS) at beta frequency has previously been shown to increase local activity in the beta band [14] and to modulate long-range connectivity within the default mode network [15]. We demonstrated that beta-frequency tACS significantly changed the connectivity pattern of the stimulated primary motor cortex (M1), without changing overall local activity or network connectivity. This finding is supported by a simple phase-precession model, which demonstrates the plausibility of the results and provides emergent predictions that are consistent with our empirical findings. These findings therefore inform our understanding of how local oscillatory activity may underpin network connectivity.

•

tACS does not alter overall functional connectivity between major network nodes

•

However, tACS modulates the connectivity pattern of the stimulated motor cortex

•

These data directly support the “communication through coherence” hypothesis

•

We provide evidence for how disordered connectivity arises from oscillatory changes

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.

Complementary contributions of concurrent EEG and fMRI connectivity for predicting structural connectivity

While averaged dynamics of brain function are known to estimate the underlying structure, the exact relationship between large-scale function and structure remains an unsolved issue in network neuroscience. These complex functional dynamics, measured by EEG and fMRI, are thought to arise from a shared underlying structural architecture, which can be measured by diffusion MRI (dMRI). While simulation and data transformation (e.g. graph theory measures) have been proposed to refine the understanding of the underlying function-structure relationship, the potential complementary and/or independent contribution of EEG and fMRI to this relationship is still poorly understood. As such, we explored this relationship by analyzing the function-structure correlation in fourteen healthy subjects with simultaneous resting-state EEG-fMRI and dMRI acquisitions. We show that the combination of EEG and fMRI connectivity better explains dMRI connectivity and that this represents a genuine model improvement over fMRI-only models for both group-averaged connectivity matrices and at the individual level. Furthermore, this model improves the prediction within each resting-state network. The best model fit to underlying structure is mediated by fMRI and EEG-δ connectivity in combination with Euclidean distance and interhemispheric connectivity with more local contributions of EEG-γ at the scale of resting-state networks. This highlights that the factors mediating the relationship between functional and structural metrics of connectivity are context and scale dependent, influenced by topological, geometric and architectural features. It also suggests that fMRI studies employing simultaneous EEG measures may characterize additional and essential parts of the underlying neuronal activity of the resting-state, which might be of special interest for both clinical studies and the investigation of resting-state dynamics.

Low-frequency hippocampal-cortical activity drives brain-wide resting-state functional MRI connectivity

The hippocampus, including the dorsal dentate gyrus (dDG), and cortex engage in bidirectional communication. We propose that low-frequency activity in hippocampal-cortical pathways contributes to brain-wide resting-state connectivity to integrate sensory information. Using optogenetic stimulation and brain-wide fMRI and resting-state fMRI (rsfMRI), we determined the large-scale effects of spatiotemporal-specific downstream propagation of hippocampal activity. Low-frequency (1 Hz), but not high-frequency (40 Hz), stimulation of dDG excitatory neurons evoked robust cortical and subcortical brain-wide fMRI responses. More importantly, it enhanced interhemispheric rsfMRI connectivity in various cortices and hippocampus. Subsequent local field potential recordings revealed an increase in slow oscillations in dorsal hippocampus and visual cortex, interhemispheric visual cortical connectivity, and hippocampal-cortical connectivity. Meanwhile, pharmacological inactivation of dDG neurons decreased interhemispheric rsfMRI connectivity. Functionally, visually evoked fMRI responses in visual regions also increased during and after low-frequency dDG stimulation. Together, our results indicate that low-frequency activity robustly propagates in the dorsal hippocampal-cortical pathway, drives interhemispheric cortical rsfMRI connectivity, and mediates visual processing.

Impact of Visual Corticostriatal Loop Disruption on Neural Processing within the Parahippocampal Place Area

The caudate nucleus is a part of the visual corticostriatal loop (VCSL), receiving input from different visual areas and projecting back to the same cortical areas via globus pallidus, substantia nigra, and thalamus. Despite perceptual and navigation impairments in patients with VCSL disruption due to caudate atrophy (e.g., Huntington's disease, HD), the relevance of the caudate nucleus and VCSL on cortical visual processing is not fully understood. In a series of fMRI experiments, we found that the caudate showed a stronger functional connection to parahippocampal place area (PPA) compared with adjacent regions (e.g., fusiform face area, FFA) within the temporal visual cortex. Consistent with this functional link, the caudate showed a higher response to scenes compared with faces, similar to the PPA. Testing the impact of VCSL disruption on neural processes within PPA, HD patients showed reduced scene-selective activity within PPA compared with healthy matched controls. In contrast, the level of selective activity in adjacent cortical and subcortical face-selective areas (i.e., FFA and amygdala) remained intact. These results show some of the first evidence for the direct impact and potential clinical significance of VCSL on the generation of "selective" activity within PPA.

SIGNIFICANCE STATEMENT

Visual perception is often considered the product of a multistage feedforward neural processing between visual cortical areas, ignoring the likely impact of corticosubcortical loops on this process. Here, we provide evidence for the contribution of visual corticostriatal loop and the caudate nucleus on generating selective response within parahippocampal place area (PPA). Our results show that disruption of this loop in Huntington's disease patients reduces the level of selective activity within PPA, which may lead to related perceptual impairments in these patients.

Organization and hierarchy of the human functional brain network lead to a chain-like core

The brain is a paradigmatic example of a complex system: its functionality emerges as a global property of local mesoscopic and microscopic interactions. Complex network theory allows to elicit the functional architecture of the brain in terms of links (correlations) between nodes (grey matter regions) and to extract information out of the noise. Here we present the analysis of functional magnetic resonance imaging data from forty healthy humans at rest for the investigation of the basal scaffold of the functional brain network organization. We show how brain regions tend to coordinate by forming a highly hierarchical chain-like structure of homogeneously clustered anatomical areas. A maximum spanning tree approach revealed the centrality of the occipital cortex and the peculiar aggregation of cerebellar regions to form a closed core. We also report the hierarchy of network segregation and the level of clusters integration as a function of the connectivity strength between brain regions.

The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation.

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.

Evaluation of the effectiveness of transcranial direct current stimulation (tDCS) and psychosensory stimulation through DOCS scale in a minimally conscious subject

The aim of our study was to assess the effectiveness of transcranial direct current stimulation (tDCS) on alertness improvement in a patient in a minimally conscious state (MCS) by means of disorders of consciousness scale combined with psycho-sensory stimulation. The effects of tDCS on muscle hypertonia through the Ashworth scale were also examined. tDCS was performed through a two-channel intra-cephalic stimulator. After stimulation, the patient followed a psychosensory stimulation training. Results pointed out an increase in DOCunit score, as well as an increase in alertness maintenance and an improvement in muscle hypertonia, although a MCS state persisted.

Acetylcholine Release in Prefrontal Cortex Promotes Gamma Oscillations and Theta-Gamma Coupling during Cue Detection

The capacity for using external cues to guide behavior ("cue detection") constitutes an essential aspect of attention and goal-directed behavior. The cortical cholinergic input system, via phasic increases in prefrontal acetylcholine release, plays an essential role in attention by mediating such cue detection. However, the relationship between cholinergic signaling during cue detection and neural activity dynamics in prefrontal networks remains unclear. Here we combined subsecond measures of cholinergic signaling, neurophysiological recordings, and cholinergic receptor blockade to delineate the cholinergic contributions to prefrontal oscillations during cue detection in rats. We first confirmed that detected cues evoke phasic acetylcholine release. These cholinergic signals were coincident with increased neuronal synchrony across several frequency bands and the emergence of theta-gamma coupling. Muscarinic and nicotinic cholinergic receptors both contributed specifically to gamma synchrony evoked by detected cues, but the effects of blocking the two receptor subtypes were dissociable. Blocking nicotinic receptors primarily attenuated high-gamma oscillations occurring during the earliest phases of the cue detection process, while muscarinic (M1) receptor activity was preferentially involved in the transition from high to low gamma power that followed and corresponded to the mobilization of networks involved in cue-guided decision making. Detected cues also promoted coupling between gamma and theta oscillations, and both nicotinic and muscarinic receptor activity contributed to this process. These results indicate that acetylcholine release coordinates neural oscillations during the process of cue detection.SIGNIFICANCE STATEMENT The capacity of learned cues to direct attention and guide responding ("cue detection") is a key component of goal-directed behavior. Rhythmic neural activity and increases in acetylcholine release in the prefrontal cortex contribute to this process; however, the relationship between these neuronal mechanisms is not well understood. Using a combination of in vivo neurochemistry, neurophysiology, and pharmacological methods, we demonstrate that cue-evoked acetylcholine release, through distinct actions at both nicotinic and muscarinic receptors, triggers a procession of neural oscillations that map onto the multiple stages of cue detection. Our data offer new insights into cholinergic function by revealing the temporally orchestrated changes in prefrontal network synchrony modulated by acetylcholine release during cue detection.

Frequency-specific electrophysiologic correlates of resting state fMRI networks

Resting state functional MRI (R-fMRI) studies have shown that slow (< 0.1 Hz), intrinsic fluctuations of the blood oxygen level dependent (BOLD) signal are temporally correlated within hierarchically organized functional systems known as resting state networks (RSNs) (Doucet et al., 2011). Most broadly, this hierarchy exhibits a dichotomy between two opposed systems (Fox et al., 2005). One system engages with the environment and includes the visual, auditory, and sensorimotor (SMN) networks as well as the dorsal attention network (DAN), which controls spatial attention. The other system includes the default mode network (DMN) and the fronto-parietal control system (FPC), RSNs that instantiate episodic memory and executive control, respectively. Here, we test the hypothesis, based on the spectral specificity of electrophysiologic responses to perceptual vs. memory tasks (Klimesch, 1999; Pfurtscheller and Lopes da Silva, 1999), that these two large-scale neural systems also manifest frequency specificity in the resting state. We measured the spatial correspondence between electrocorticographic (ECoG) band-limited power (BLP) and R-fMRI correlation patterns in awake, resting, human subjects. Our results show that, while gamma BLP correspondence was common throughout the brain, theta (4–8 Hz) BLP correspondence was stronger in the DMN and FPC, whereas alpha (8–12 Hz) correspondence was stronger in the SMN and DAN. Thus, the human brain, at rest, exhibits frequency specific electrophysiology, respecting both the spectral structure of task responses and the hierarchical organization of RSNs.

Studying Brain Circuit Function with Dynamic Causal Modeling for Optogenetic fMRI

Defining the large-scale behavior of brain circuits with cell type specificity is a major goal of neuroscience. However, neuronal circuit diagrams typically draw upon anatomical and electrophysiological measurements acquired in isolation. Consequently, a dynamic and cell-type-specific connectivity map has never been constructed from simultaneous measurements across the brain. Here, we introduce dynamic causal modeling (DCM) for optogenetic fMRI experiments-which uniquely allow cell-type-specific, brain-wide functional measurements-to parameterize the causal relationships among regions of a distributed brain network with cell type specificity. Strikingly, when applied to the brain-wide basal ganglia-thalamocortical network, DCM accurately reproduced the empirically observed time series, and the strongest connections were key connections of optogenetically stimulated pathways. We predict that quantitative and cell-type-specific descriptions of dynamic connectivity, as illustrated here, will empower novel systems-level understanding of neuronal circuit dynamics and facilitate the design of more effective neuromodulation therapies.

Long-range projections coordinate distributed brain-wide neural activity with a specific spatiotemporal profile

One challenge in contemporary neuroscience is to achieve an integrated understanding of the large-scale brain-wide interactions, particularly the spatiotemporal patterns of neural activity that give rise to functions and behavior. At present, little is known about the spatiotemporal properties of long-range neuronal networks. We examined brain-wide neural activity patterns elicited by stimulating ventral posteromedial (VPM) thalamo-cortical excitatory neurons through combined optogenetic stimulation and functional MRI (fMRI). We detected robust optogenetically evoked fMRI activation bilaterally in primary visual, somatosensory, and auditory cortices at low (1 Hz) but not high frequencies (5-40 Hz). Subsequent electrophysiological recordings indicated interactions over long temporal windows across thalamo-cortical, cortico-cortical, and interhemispheric callosal projections at low frequencies. We further observed enhanced visually evoked fMRI activation during and after VPM stimulation in the superior colliculus, indicating that visual processing was subcortically modulated by low-frequency activity originating from VPM. Stimulating posteromedial complex thalamo-cortical excitatory neurons also evoked brain-wide blood-oxygenation-level-dependent activation, although with a distinct spatiotemporal profile. Our results directly demonstrate that low-frequency activity governs large-scale, brain-wide connectivity and interactions through long-range excitatory projections to coordinate the functional integration of remote brain regions. This low-frequency phenomenon contributes to the neural basis of long-range functional connectivity as measured by resting-state fMRI.

Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach

Phase-amplitude coupling (PAC) plays an important role in neural communication and computation. Interestingly, recent studies have indicated the presence of ubiquitous PAC phenomenon even during the resting state. Despite the importance of PAC phenomenon, estimation of significant physiological PAC is challenging because of the lack of appropriate surrogate measures to control false positives caused by non-physiological PAC. Therefore, in the present study, we evaluated PAC phenomenon during resting-state magnetoencephalography (MEG) signal and considered various surrogate measures and computational approaches widely used in the literature in addition to proposing new ones. We evaluated PAC phenomenon over the entire length of the MEG signal and for multiple shorter time segments. The results indicate that the extent of PAC phenomenon mainly depends on the surrogate measures and PAC computational methods used, as well as the evaluation approach. After a careful and critical evaluation, we found that resting-state MEG signals failed to exhibit ubiquitous PAC phenomenon, contrary to what has been suggested previously.

Topology, Cross-Frequency, and Same-Frequency Band Interactions Shape the Generation of Phase-Amplitude Coupling in a Neural Mass Model of a Cortical Column

Phase-amplitude coupling (PAC), a type of cross-frequency coupling (CFC) where the phase of a low-frequency rhythm modulates the amplitude of a higher frequency, is becoming an important indicator of information transmission in the brain. However, the neurobiological mechanisms underlying its generation remain undetermined. A realistic, yet tractable computational model of the phenomenon is thus needed. Here we analyze a neural mass model of a cortical column, comprising fourteen neuronal populations distributed across four layers (L2/3, L4, L5 and L6). A control analysis showed that the conditional transfer entropy (cTE) measure is able to correctly estimate the flow of information between neuronal populations. Then, we computed cTE from the phases to the amplitudes of the oscillations generated in the cortical column. This approach provides information regarding directionality by distinguishing PAC from APC (amplitude-phase coupling), i.e. the information transfer from amplitudes to phases, and can be used to estimate other types of CFC such as amplitude-amplitude coupling (AAC) and phase-phase coupling (PPC). While experiments often only focus on one or two PAC combinations (e.g., theta-gamma or alpha-gamma), we found that a cortical column can simultaneously generate almost all possible PAC combinations, depending on connectivity parameters, time constants, and external inputs. PAC interactions with and without an anatomical equivalent (direct and indirect interactions, respectively) were analyzed. We found that the strength of PAC between two populations was strongly correlated with the strength of the effective connections between the populations and, on average, did not depend on whether the PAC connection was direct or indirect. When considering a cortical column circuit as a complex network, we found that neuronal populations making indirect PAC connections had, on average, higher local clustering coefficient, efficiency, and betweenness centrality than populations making direct connections and populations not involved in PAC connections. This suggests that their interactions were more effective when transmitting information. Since approximately 60% of the obtained interactions represented indirect connections, our results highlight the importance of the topology of cortical circuits for the generation of the PAC phenomenon. Finally, our results demonstrated that indirect PAC interactions can be explained by a cascade of direct CFC and same-frequency band interactions, suggesting that PAC analysis of experimental data should be accompanied by the estimation of other types of frequency interactions for an integrative understanding of the phenomenon.

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.

Activation of Direct and Indirect Pathway Medium Spiny Neurons Drives Distinct Brain-wide Responses

A central theory of basal ganglia function is that striatal neurons expressing the D1 and D2 dopamine receptors exert opposing brain-wide influences. However, the causal influence of each population has never been measured at the whole-brain scale. Here, we selectively stimulated D1 or D2 receptor-expressing neurons while visualizing whole-brain activity with fMRI. Excitation of either inhibitory population evoked robust positive BOLD signals within striatum, while downstream regions exhibited significantly different and generally opposing responses consistent with-though not easily predicted from-contemporary models of basal ganglia function. Importantly, positive and negative signals within the striatum, thalamus, GPi, and STN were all associated with increases and decreases in single-unit activity, respectively. These findings provide direct evidence for the opposing influence of D1 and D2 receptor-expressing striatal neurons on brain-wide circuitry and extend the interpretability of fMRI studies by defining cell-type-specific contributions to the BOLD signal.

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.

Aberrant Spontaneous and Task-Dependent Functional Connections in the Anxious Brain

A number of brain regions have been implicated in the anxiety disorders, yet none of these regions in isolation has been distinguished as the sole or discrete site responsible for anxiety disorder pathology. Therefore, the identification of dysfunctional neural networks as represented by alterations in the temporal correlation of blood-oxygen level dependent (BOLD) signal across several brain regions in anxiety disorders has been increasingly pursued in the past decade. Here, we review task-independent (e.g., resting state) and task-induced functional connectivity magnetic resonance imaging (fcMRI) studies in the adult anxiety disorders (including trauma- and stressor-related and obsessive compulsive disorders). The results of this review suggest that anxiety disorder pathophysiology involves aberrant connectivity between amygdala-frontal and frontal-striatal regions, as well as within and between canonical "intrinsic" brain networks - the default mode and salience networks, and that evidence of these aberrations may help inform findings of regional activation abnormalities observed in the anxiety disorders. Nonetheless, significant challenges remain, including the need to better understand mixed findings observed using different methods (e.g., resting state and task-based approaches); the need for more developmental work; the need to delineate disorder-specific and transdiagnostic fcMRI aberrations in the anxiety disorders; and the need to better understand the clinical significance of fcMRI abnormalities. In meeting these challenges, future work has the potential to elucidate aberrant neural networks as intermediate, brain-based phenotypes to predict disease onset and progression, refine diagnostic nosology, and ascertain treatment mechanisms and predictors of treatment response across anxiety, trauma-related and obsessive compulsive disorders.

Underconnected, But Not Broken? Dynamic Functional Connectivity MRI Shows Underconnectivity in Autism Is Linked to Increased Intra-Individual Variability Across Time

Autism spectrum disorder (ASD) is characterized by core sociocommunicative impairments. Atypical intrinsic functional connectivity (iFC) has been reported in numerous studies of ASD. A majority of findings has indicated long-distance underconnectivity. However, fMRI studies have thus far exclusively examined static iFC across several minutes of scanning. We examined temporal variability of iFC, using sliding window analyses in selected high-quality (low-motion) consortium datasets from 76 ASD and 76 matched typically developing (TD) participants (Study 1) and in-house data from 32 ASD and 32 TD participants. Mean iFC and standard deviation of the sliding window correlation (SD-iFC) were computed for regions of interest (ROIs) from default mode and salience networks, as well as amygdala and thalamus. In both studies, ROI pairings with significant underconnectivity (ASD<TD) were identified. Mediation analyses showed that decreased mean iFC in the ASD groups was significantly affected by increased SD-iFC. Our study is the first to identify temporal variability across time as a significant contributing factor to the common finding of static underconnectivity in ASD. Since peak connectivity across time was not significantly reduced in ASD, static underconnectivity findings may have to be reinterpreted, suggesting that connections are not actually "broken" in ASD, but subject to greater intra-individual variability across time. Our findings indicate the need for dynamic approaches to iFC in clinical functional connectivity MRI (fcMRI) investigations.

Oxytocin reduces amygdala activity, increases social interactions, and reduces anxiety-like behavior irrespective of NMDAR antagonism

Standard dopamine therapies for schizophrenia are not efficacious for negative symptoms of the disease, including asociality. This reduced social behavior may be due to glutamatergic dysfunction within the amygdala, leading to increased fear and social anxiety. Several studies have demonstrated the prosocial effects of oxytocin in schizophrenia patients. Therefore, this study evaluates the effect of subchronic oxytocin on EEG activity in amygdala of mice during performance of the three-chamber social choice and open field tests following acute ketamine as a model of glutamatergic dysfunction. Oxytocin did not restore social deficits introduced by ketamine but did significantly increase sociality in comparison to the control group. Ketamine had no effect on time spent in the center during the open field trials, whereas oxytocin increased overall center time across all groups, suggesting a reduction in anxiety. Amygdala activity was consistent across all drug groups during social and nonsocial behavioral trials. However, oxytocin reduced overall amygdala EEG power during the two behavioral tasks. Alternatively, ketamine did not significantly affect EEG power throughout the tasks. Decreased EEG power in the amygdala, as caused by oxytocin, may be related to both reduced anxiety and increased social behaviors. Data suggest that separate prosocial and social anxiety pathways may mediate social preference.

Dynamic regional phase synchrony (DRePS): An Instantaneous Measure of Local fMRI Connectivity Within Spatially Clustered Brain Areas

Dynamic functional brain connectivity analysis is a fast expanding field in computational neuroscience research with the promise of elucidating brain network interactions. Sliding temporal window based approaches are commonly used in order to explore dynamic behavior of brain networks in task-free functional magnetic resonance imaging (fMRI) data. However, the low effective temporal resolution of sliding window methods fail to capture the full dynamics of brain activity at each time point. These also require subjective decisions regarding window size and window overlap. In this study, we introduce dynamic regional phase synchrony (DRePS), a novel analysis approach that measures mean local instantaneous phase coherence within adjacent fMRI voxels. We evaluate the DRePS framework on simulated data showing that the proposed measure is able to estimate synchrony at higher temporal resolution than sliding windows of local connectivity. We applied DRePS analysis to task-free fMRI data of 20 control subjects, revealing ultra-slow dynamics of local connectivity in different brain areas. Spatial clustering based on the DRePS feature time series reveals biologically congruent local phase synchrony networks (LPSNs). Taken together, our results demonstrate three main findings. Firstly, DRePS has increased temporal sensitivity compared to sliding window correlation analysis in capturing locally synchronous events. Secondly, DRePS of task-free fMRI reveals ultra-slow fluctuations of ∼0.002-0.02 Hz. Lastly, LPSNs provide plausible spatial information about time-varying brain local phase synchrony. With the DRePS method, we introduce a framework for interrogating brain local connectivity, which can potentially provide biomarkers of human brain function in health and disease. Hum Brain Mapp 37:1970-1985, 2016. © 2016 Wiley Periodicals, Inc.

Recovery of slow-5 oscillations in a longitudinal study of ischemic stroke patients

Functional networks in resting-state fMRI are identified by characteristics of their intrinsic low-frequency oscillations, more specifically in terms of their synchronicity. With advanced aging and in clinical populations, this synchronicity among functionally linked regions is known to decrease and become disrupted, which may be associated with observed cognitive and behavioral changes. Previous work from our group has revealed that oscillations within the slow-5 frequency range (0.01–0.027 Hz) are particularly susceptible to disruptions in aging and following a stroke. In this study, we characterized longitudinally the changes in the slow-5 oscillations in stroke patients across two different time-points. We followed a group of ischemic stroke patients (n = 20) and another group of healthy older adults (n = 14) over two visits separated by a minimum of three months (average of 9 months). For the stroke patients, one visit occurred in their subacute window (10 days to 6 months after stroke onset), the other took place in their chronic window (> 6 months after stroke). Using a mid-order group ICA method on 10-minutes eyes-closed resting-state fMRI data, we assessed the frequency distributions of a component's representative time-courses for differences in regards to slow-5 spectral power. First, our stroke patients, in their subacute stage, exhibited lower amplitude slow-5 oscillations in comparison to their healthy counterparts. Second, over time in their chronic stage, those same patients showed a recovery of those oscillations, reaching near equivalence to the healthy older adult group. Our results indicate the possibility of an eventual recovery of those initially disrupted network oscillations to a near-normal level, providing potentially a biomarker for stroke recovery of the cortical system. This finding opens new avenues in infra-slow oscillation research and could serve as a useful biomarker in future treatments aimed at recovery.

•

Slow-5 oscillation amplitudes are reduced in stroke patients at the subacute stage.

•

Slow-5 oscillation amplitudes correlate with cognitive performance.

•

Slow-5 oscillations recover in the same patients at the chronic stage.

•

Findings support the high implication of slow-5 oscillations in network disruption.

•

Slow-5 oscillations may serve as a bio-marker of functional network health.

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.

Pain Related Cortical Oscillations: Methodological Advances and Potential Applications

Alongside the time-locked event-related potentials (ERPs), nociceptive somatosensory inputs can induce modulations of ongoing oscillations, appeared as event-related synchronization or desynchronization (ERS/ERD) in different frequency bands. These ERD/ERS activities are suggested to reflect various aspects of pain perception, including the representation, encoding, assessment, and integration of the nociceptive sensory inputs, as well as behavioral responses to pain, even the precise details of their roles remain unclear. Previous studies investigating the functional relevance of ERD/ERS activities in pain perception were normally done by assessing their latencies, frequencies, magnitudes, and scalp distributions, which would be then correlated with subjective pain perception or stimulus intensity. Nevertheless, these temporal, spectral, and spatial profiles of stimulus induced ERD/ERS could only partly reveal the dynamics of brain oscillatory activities. Indeed, additional parameters, including but not limited to, phase, neural generator, and cross frequency couplings, should be paid attention to comprehensively and systemically evaluate the dynamics of oscillatory activities associated with pain perception and behavior. This would be crucial in exploring the psychophysiological mechanisms of neural oscillation, and in understanding the neural functions of cortical oscillations involved in pain perception and behavior. Notably, some chronic pain (e.g., neurogenic pain and complex regional pain syndrome) patients are often associated with the occurrence of abnormal synchronized oscillatory brain activities, and selectively modulating cortical oscillatory activities has been showed to be a potential therapy strategy to relieve pain with the application of neurostimulation techniques, e.g., repeated transcranial magnetic stimulation (rTMS) and transcranial alternating current stimulation (tACS). Thus, the investigation of the oscillatory activities proceeding from phenomenology to function, opens new perspectives to address questions in human pain psychophysiology and pathophysiology, thereby promoting the establishment of rational therapeutic strategy.

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.

Intrinsic Functional Connectivity Patterns Predict Consciousness Level and Recovery Outcome in Acquired Brain Injury

For accurate diagnosis and prognostic prediction of acquired brain injury (ABI), it is crucial to understand the neurobiological mechanisms underlying loss of consciousness. However, there is no consensus on which regions and networks act as biomarkers for consciousness level and recovery outcome in ABI. Using resting-state fMRI, we assessed intrinsic functional connectivity strength (FCS) of whole-brain networks in a large sample of 99 ABI patients with varying degrees of consciousness loss (including fully preserved consciousness state, minimally conscious state, unresponsive wakefulness syndrome/vegetative state, and coma) and 34 healthy control subjects. Consciousness level was evaluated using the Glasgow Coma Scale and Coma Recovery Scale-Revised on the day of fMRI scanning; recovery outcome was assessed using the Glasgow Outcome Scale 3 months after the fMRI scanning. One-way ANOVA of FCS, Spearman correlation analyses between FCS and the consciousness level and recovery outcome, and FCS-based multivariate pattern analysis were performed. We found decreased FCS with loss of consciousness primarily distributed in the posterior cingulate cortex/precuneus (PCC/PCU), medial prefrontal cortex, and lateral parietal cortex. The FCS values of these regions were significantly correlated with consciousness level and recovery outcome. Multivariate support vector machine discrimination analysis revealed that the FCS patterns predicted whether patients with unresponsive wakefulness syndrome/vegetative state and coma would regain consciousness with an accuracy of 81.25%, and the most discriminative region was the PCC/PCU. These findings suggest that intrinsic functional connectivity patterns of the human posteromedial cortex could serve as a potential indicator for consciousness level and recovery outcome in individuals with ABI.

SIGNIFICANCE STATEMENT

Varying degrees of consciousness loss and recovery are commonly observed in acquired brain injury patients, yet the underlying neurobiological mechanisms remain elusive. Using a large sample of patients with varying degrees of consciousness loss, we demonstrate that intrinsic functional connectivity strength in many brain regions, especially in the posterior cingulate cortex and precuneus, significantly correlated with consciousness level and recovery outcome. We further demonstrate that the functional connectivity pattern of these regions can predict patients with unresponsive wakefulness syndrome/vegetative state and coma would regain consciousness with an accuracy of 81.25%. Our study thus provides potentially important biomarkers of acquired brain injury in clinical diagnosis, prediction of recovery outcome, and decision making for treatment strategies for patients with severe loss of consciousness.

Structural and functional connectivity between the lateral posterior-pulvinar complex and primary visual cortex in the ferret

The role of higher-order thalamic structures in sensory processing remains poorly understood. Here, we used the ferret (Mustela putorius furo) as a novel model species for the study of the lateral posterior (LP)-pulvinar complex and its structural and functional connectivity with area 17 [primary visual cortex (V1)]. We found reciprocal anatomical connections between the lateral part of the LP nucleus of the LP-pulvinar complex (LPl) and V1. In order to investigate the role of this feedback loop between LPl and V1 in shaping network activity, we determined the functional interactions between LPl and the supragranular, granular and infragranular layers of V1 by recording multiunit activity and local field potentials. Coherence was strongest between LPl and the supragranular V1, with the most distinct peaks in the delta and alpha frequency bands. Inter-area interaction measured by spike-phase coupling identified the delta frequency band being dominated by the infragranular V1 and multiple frequency bands that were most pronounced in the supragranular V1. This inter-area coupling was differentially modulated by full-field synthetic and naturalistic visual stimulation. We also found that visual responses in LPl were distinct from those in V1 in terms of their reliability. Together, our data support a model of multiple communication channels between LPl and the layers of V1 that are enabled by oscillations in different frequency bands. This demonstration of anatomical and functional connectivity between LPl and V1 in ferrets provides a roadmap for studying the interaction dynamics during behaviour, and a template for identifying the activity dynamics of other thalamo-cortical feedback loops.

Attention Decreases Phase-Amplitude Coupling, Enhancing Stimulus Discriminability in Cortical Area MT

Local field potentials (LFPs) in cortex reflect synchronous fluctuations in the synaptic activity of local populations of neurons. The power of high frequency (>30 Hz) oscillations in LFPs is locked to the phase of low frequency (<30 Hz) oscillations, an effect known as phase-amplitude coupling (PAC). While PAC has been observed in a variety of cortical regions and animal models, its functional role particularly in primate visual cortex is largely unknown. Here, we document PAC for LFPs recorded from extra-striate area MT of macaque monkeys, an area specialized for the processing of visual motion. We further show that directing spatial attention into the receptive field of MT neurons decreases the coupling between the low frequency phase and high frequency power of LFPs. This attentional suppression of PAC increases neuronal discriminability for attended visual stimuli. Therefore, we hypothesize that visual cortex uses PAC to regulate inter-neuronal correlations and thereby enhances the coding of relevant stimuli.

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.

The modular and integrative functional architecture of the human brain

Network-based analyses of brain imaging data consistently reveal distinct modules and connector nodes with diverse global connectivity across the modules. How discrete the functions of modules are, how dependent the computational load of each module is to the other modules' processing, and what the precise role of connector nodes is for between-module communication remains underspecified. Here, we use a network model of the brain derived from resting-state functional MRI (rs-fMRI) data and investigate the modular functional architecture of the human brain by analyzing activity at different types of nodes in the network across 9,208 experiments of 77 cognitive tasks in the BrainMap database. Using an author-topic model of cognitive functions, we find a strong spatial correspondence between the cognitive functions and the network's modules, suggesting that each module performs a discrete cognitive function. Crucially, activity at local nodes within the modules does not increase in tasks that require more cognitive functions, demonstrating the autonomy of modules' functions. However, connector nodes do exhibit increased activity when more cognitive functions are engaged in a task. Moreover, connector nodes are located where brain activity is associated with many different cognitive functions. Connector nodes potentially play a role in between-module communication that maintains the modular function of the brain. Together, these findings provide a network account of the brain's modular yet integrated implementation of cognitive functions.

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.

Modeling the Generation of Phase-Amplitude Coupling in Cortical Circuits: From Detailed Networks to Neural Mass Models

Phase-amplitude coupling (PAC), the phenomenon where the amplitude of a high frequency oscillation is modulated by the phase of a lower frequency oscillation, is attracting an increasing interest in the neuroscience community due to its potential relevance for understanding healthy and pathological information processing in the brain. PAC is a diverse phenomenon, having been experimentally detected in at least ten combinations of rhythms: delta-theta, delta-alpha, delta-beta, delta-gamma, theta-alpha, theta-beta, theta-gamma, alpha-beta, alpha-gamma, and beta-gamma. However, a complete understanding of the biophysical mechanisms generating this diversity is lacking. Here we review computational models of PAC generation that range from detailed models of neuronal networks, where each cell is described by Hodgkin-Huxley-type equations, to neural mass models (NMMs) where only the average activities of neuronal populations are considered. We argue that NMMs are an appropriate mathematical framework (due to the small number of parameters and variables involved and the richness of the dynamics they can generate) to study the PAC phenomenon.

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.

Task-sensitive reconfiguration of corticocortical 6-20 Hz oscillatory coherence in naturalistic human performance

Electrophysiological oscillatory coherence between brain regions has been proposed to facilitate functional long-range connectivity within neurocognitive networks. This notion is supported by intracortical recordings of coherence in singled-out corticocortical connections in the primate cortex. However, the manner in which this operational principle manifests in the task-sensitive connectivity that supports human naturalistic performance remains undercharacterized. Here, we demonstrate task-sensitive reconfiguration of global patterns of coherent connectivity in association with a set of easier and more demanding naturalistic tasks, ranging from picture comparison to speech comprehension and object manipulation. Based on whole-cortex neuromagnetic recording in healthy behaving individuals, the task-sensitive component of long-range corticocortical coherence was mapped at spectrally narrow-band oscillatory frequencies between 6 and 20 Hz (theta to alpha and low-beta bands). This data-driven cortical mapping unveiled markedly distinct and topologically task-relevant spatiospectral connectivity patterns for the different tasks. The results demonstrate semistable oscillatory states relevant for neurocognitive processing. The present findings decisively link human behavior to corticocortical coherence at oscillatory frequencies that are widely thought to convey long-range, feedback-type neural interaction in cortical functional networks.

Top-down control of the phase of alpha-band oscillations as a mechanism for temporal prediction

The physiological state of the brain before an incoming stimulus has substantial consequences for subsequent behavior and neural processing. For example, the phase of ongoing posterior alpha-band oscillations (8-14 Hz) immediately before visual stimulation has been shown to predict perceptual outcomes and downstream neural activity. Although this phenomenon suggests that these oscillations may phasically route information through functional networks, many accounts treat these periodic effects as a consequence of ongoing activity that is independent of behavioral strategy. Here, we investigated whether alpha-band phase can be guided by top-down control in a temporal cueing task. When participants were provided with cues predictive of the moment of visual target onset, discrimination accuracy improved and targets were more frequently reported as consciously seen, relative to unpredictive cues. This effect was accompanied by a significant shift in the phase of alpha-band oscillations, before target onset, toward each participant's optimal phase for stimulus discrimination. These findings provide direct evidence that forming predictions about when a stimulus will appear can bias the phase of ongoing alpha-band oscillations toward an optimal phase for visual processing, and may thus serve as a mechanism for the top-down control of visual processing guided by temporal predictions.

Dynamic network communication as a unifying neural basis for cognition, development, aging, and disease

Perception, cognition, and social interaction depend upon coordinated neural activity. This coordination operates within noisy, overlapping, and distributed neural networks operating at multiple timescales. These networks are built upon a structural scaffolding with intrinsic neuroplasticity that changes with development, aging, disease, and personal experience. In this article, we begin from the perspective that successful interregional communication relies upon the transient synchronization between distinct low-frequency (<80 Hz) oscillations, allowing for brief windows of communication via phase-coordinated local neuronal spiking. From this, we construct a theoretical framework for dynamic network communication, arguing that these networks reflect a balance between oscillatory coupling and local population spiking activity and that these two levels of activity interact. We theorize that when oscillatory coupling is too strong, spike timing within the local neuronal population becomes too synchronous; when oscillatory coupling is too weak, spike timing is too disorganized. Each results in specific disruptions to neural communication. These alterations in communication dynamics may underlie cognitive changes associated with healthy development and aging, in addition to neurological and psychiatric disorders. A number of neurological and psychiatric disorders-including Parkinson's disease, autism, depression, schizophrenia, and anxiety-are associated with abnormalities in oscillatory activity. Although aging, psychiatric and neurological disease, and experience differ in the biological changes to structural gray or white matter, neurotransmission, and gene expression, our framework suggests that any resultant cognitive and behavioral changes in normal or disordered states or their treatment are a product of how these physical processes affect dynamic network communication.

Hippocampal-dorsolateral prefrontal coupling as a species-conserved cognitive mechanism: a human translational imaging study

Hippocampal-prefrontal cortex (HC-PFC) interactions are implicated in working memory (WM) and altered in psychiatric conditions with cognitive impairment such as schizophrenia. While coupling between both structures is crucial for WM performance in rodents, evidence from human studies is conflicting and translation of findings is complicated by the use of differing paradigms across species. We therefore used functional magnetic resonance imaging together with a spatial WM paradigm adapted from rodent research to examine HC-PFC coupling in humans. A PFC-parietal network was functionally connected to hippocampus (HC) during task stages requiring high levels of executive control but not during a matched control condition. The magnitude of coupling in a network comprising HC, bilateral dorsolateral PFC (DLPFC), and right supramarginal gyrus explained one-fourth of the variability in an independent spatial WM task but was unrelated to visual WM performance. HC-DLPFC coupling may thus represent a systems-level mechanism specific to spatial WM that is conserved across species, suggesting its utility for modeling cognitive dysfunction in translational neuroscience.