Multiplex Networks of Cortical and Hippocampal Neurons Revealed at Different Timescales

Recent studies have emphasized the importance of multiplex networks – interdependent networks with shared nodes and different types of connections – in systems primarily outside of neuroscience. Though the multiplex properties of networks are frequently not considered, most networks are actually multiplex networks and the multiplex specific features of networks can greatly affect network behavior (e.g. fault tolerance). Thus, the study of networks of neurons could potentially be greatly enhanced using a multiplex perspective. Given the wide range of temporally dependent rhythms and phenomena present in neural systems, we chose to examine multiplex networks of individual neurons with time scale dependent connections. To study these networks, we used transfer entropy – an information theoretic quantity that can be used to measure linear and nonlinear interactions – to systematically measure the connectivity between individual neurons at different time scales in cortical and hippocampal slice cultures. We recorded the spiking activity of almost 12,000 neurons across 60 tissue samples using a 512-electrode array with 60 micrometer inter-electrode spacing and 50 microsecond temporal resolution. To the best of our knowledge, this preparation and recording method represents a superior combination of number of recorded neurons and temporal and spatial recording resolutions to any currently available in vivo system. We found that highly connected neurons (“hubs”) were localized to certain time scales, which, we hypothesize, increases the fault tolerance of the network. Conversely, a large proportion of non-hub neurons were not localized to certain time scales. In addition, we found that long and short time scale connectivity was uncorrelated. Finally, we found that long time scale networks were significantly less modular and more disassortative than short time scale networks in both tissue types. As far as we are aware, this analysis represents the first systematic study of temporally dependent multiplex networks among individual neurons.     

Fast two-photon in vivo imaging with three-dimensional random-access scanning in large tissue volumes

The understanding of brain computations requires methods that read out neural activity on different spatial and temporal scales. Following signal propagation and integration across a neuron and recording the concerted activity of hundreds of neurons pose distinct challenges, and the design of imaging systems has been mostly focused on tackling one of the two operations. We developed a high-resolution, acousto-optic two-photon microscope with continuous three-dimensional (3D) trajectory and random-access scanning modes that reaches near-cubic-millimeter scan range and can be adapted to imaging different spatial scales. We performed 3D calcium imaging of action potential backpropagation and dendritic spike forward propagation at sub-millisecond temporal resolution in mouse brain slices. We also performed volumetric random-access scanning calcium imaging of spontaneous and visual stimulation-evoked activity in hundreds of neurons of the mouse visual cortex in vivo. These experiments demonstrate the subcellular and network-scale imaging capabilities of our system.     

Intrinsic connectivity reveals functionally distinct cortico-hippocampal networks in the human brain

Episodic memory depends on interactions between the hippocampus and interconnected neocortical regions. Here, using data-driven analyses of resting-state functional magnetic resonance imaging (fMRI) data, we identified the networks that interact with the hippocampus—the default mode network (DMN) and a “medial temporal network” (MTN) that included regions in the medial temporal lobe (MTL) and precuneus. We observed that the MTN plays a critical role in connecting the visual network to the DMN and hippocampus. The DMN could be further divided into 3 subnetworks: a “posterior medial” (PM) subnetwork comprised of posterior cingulate and lateral parietal cortices; an “anterior temporal” (AT) subnetwork comprised of regions in the temporopolar and dorsomedial prefrontal cortex; and a “medial prefrontal” (MP) subnetwork comprised of regions primarily in the medial prefrontal cortex (mPFC). These networks vary in their functional connectivity (FC) along the hippocampal long axis and represent different kinds of information during memory-guided decision-making. Finally, a Neurosynth meta-analysis of fMRI studies suggests new hypotheses regarding the functions of the MTN and DMN subnetworks, providing a framework to guide future research on the neural architecture of episodic memory.

Episodic memory depends on interactions between the hippocampus and interconnected neocortical regions. This study uses network analyses of intrinsic brain networks at rest to identify and characterize brain networks that interact with the hippocampus and have distinct functions during memory-guided decision making.     

Noise during rest enables the exploration of the brain's dynamic repertoire

Traditionally brain function is studied through measuring physiological responses in controlled sensory, motor, and cognitive paradigms. However, even at rest, in the absence of overt goal-directed behavior, collections of cortical regions consistently show temporally coherent activity. In humans, these resting state networks have been shown to greatly overlap with functional architectures present during consciously directed activity, which motivates the interpretation of rest activity as day dreaming, free association, stream of consciousness, and inner rehearsal. In monkeys, it has been shown though that similar coherent fluctuations are present during deep anesthesia when there is no consciousness. Here, we show that comparable resting state networks emerge from a stability analysis of the network dynamics using biologically realistic primate brain connectivity, although anatomical information alone does not identify the network. We specifically demonstrate that noise and time delays via propagation along connecting fibres are essential for the emergence of the coherent fluctuations of the default network. The spatiotemporal network dynamics evolves on multiple temporal scales and displays the intermittent neuroelectric oscillations in the fast frequency regimes, 1–100 Hz, commonly observed in electroencephalographic and magnetoencephalographic recordings, as well as the hemodynamic oscillations in the ultraslow regimes, <0.1 Hz, observed in functional magnetic resonance imaging. The combination of anatomical structure and time delays creates a space–time structure in which the neural noise enables the brain to explore various functional configurations representing its dynamic repertoire.     

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic–haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

What is the relationship between functional network architectures inferred from electromagnetic and haemodynamic data? This study shows that superposition of electromagnetic networks at canonical frequency bands manifests as highly structured patterns of haemodynamic functional connectivity in the human brain, reflecting systematic variation in cytoarchitecture.     

Altered Resting State Brain Dynamics in Temporal Lobe Epilepsy Can Be Observed in Spectral Power,Functional Connectivity and Graph Theory Metrics

Despite a wealth of EEG epilepsy data that accumulated for over half a century, our ability to understand brain dynamics associated with epilepsy remains limited. Using EEG data from 15 controls and 9 left temporal lobe epilepsy (LTLE) patients, in this study we characterize how the dynamics of the healthy brain differ from the “dynamically balanced” state of the brain of epilepsy patients treated with anti-epileptic drugs in the context of resting state. We show that such differences can be observed in band power, synchronization and network measures, as well as deviations from the small world network (SWN) architecture of the healthy brain. The θ (4–7 Hz) and high α (10–13 Hz) bands showed the biggest deviations from healthy controls across various measures. In particular, patients demonstrated significantly higher power and synchronization than controls in the θ band, but lower synchronization and power in the high α band. Furthermore, differences between controls and patients in graph theory metrics revealed deviations from a SWN architecture. In the θ band epilepsy patients showed deviations toward an orderly network, while in the high α band they deviated toward a random network. These findings show that, despite the focal nature of LTLE, the epileptic brain differs in its global network characteristics from the healthy brain. To our knowledge, this is the only study to encompass power, connectivity and graph theory metrics to investigate the reorganization of resting state functional networks in LTLE patients.     

Motor Network Plasticity and Low-Frequency Oscillations Abnormalities in Patients with Brain Gliomas: A Functional MRI Study

Brain plasticity is often associated with the process of slow-growing tumor formation, which remodels neural organization and optimizes brain network function. In this study, we aimed to investigate whether motor function plasticity would display deficits in patients with slow-growing brain tumors located in or near motor areas, but who were without motor neurological deficits. We used resting-state functional magnetic resonance imaging to probe motor networks in 15 patients with histopathologically confirmed brain gliomas and 15 age-matched healthy controls. All subjects performed a motor task to help identify individual motor activity in the bilateral primary motor cortex (PMC) and supplementary motor area (SMA). Frequency-based analysis at three different frequencies was then used to investigate possible alterations in the power spectral density (PSD) of low-frequency oscillations. For each group, the average PSD was determined for each brain region and a nonparametric test was performed to determine the difference in power between the two groups. Significantly reduced inter-hemispheric functional connectivity between the left and right PMC was observed in patients compared with controls (P<0.05). We also found significantly decreased PSD in patients compared to that in controls, in all three frequency bands (low: 0.01–0.02 Hz; middle: 0.02–0.06 Hz; and high: 0.06–0.1 Hz), at three key motor regions. These findings suggest that in asymptomatic patients with brain tumors located in eloquent regions, inter-hemispheric connection may be more vulnerable. A comparison of the two approaches indicated that power spectral analysis is more sensitive than functional connectivity analysis for identifying the neurological abnormalities underlying motor function plasticity induced by slow-growing tumors.     

Altered Functional Connectivity and Small-World in Mesial Temporal Lobe Epilepsy

Background

The functional architecture of the human brain has been extensively described in terms of functional connectivity networks, detected from the low–frequency coherent neuronal fluctuations that can be observed in a resting state condition. Little is known, so far, about the changes in functional connectivity and in the topological properties of functional networks, associated with different brain diseases.

Methodology/Principal Findings

In this study, we investigated alterations related to mesial temporal lobe epilepsy (mTLE), using resting state functional magnetic resonance imaging on 18 mTLE patients and 27 healthy controls. Functional connectivity among 90 cortical and subcortical regions was measured by temporal correlation. The related values were analyzed to construct a set of undirected graphs. Compared to controls, mTLE patients showed significantly increased connectivity within the medial temporal lobes, but also significantly decreased connectivity within the frontal and parietal lobes, and between frontal and parietal lobes. Our findings demonstrated that a large number of areas in the default-mode network of mTLE patients showed a significantly decreased number of connections to other regions. Furthermore, we observed altered small-world properties in patients, along with smaller degree of connectivity, increased n-to-1 connectivity, smaller absolute clustering coefficients and shorter absolute path length.

Conclusions/Significance

We suggest that the mTLE alterations observed in functional connectivity and topological properties may be used to define tentative disease markers.     

Frequency Dependent Alterations in Regional Homogeneity of Baseline Brain Activity in Schizophrenia

Low frequency oscillations are essential in cognitive function impairment in schizophrenia. While functional connectivity can reveal the synchronization between distant brain regions, the regional abnormalities in task-independent baseline brain activity are less clear, especially in specific frequency bands. Here, we used a regional homogeneity (ReHo) method combined with resting-state functional magnetic resonance imaging to investigate low frequency spontaneous neural activity in the three different frequency bands (slow-5∶0.01–0.027 Hz; slow-4∶0.027–0.08 Hz; and typical band: 0.01–0.08 Hz) in 69 patients with schizophrenia and 62 healthy controls. Compared with controls, schizophrenia patients exhibited decreased ReHo in the precentral gyrus, middle occipital gyrus, and posterior insula, whereas increased ReHo in the medial prefrontal cortex and anterior insula. Significant differences in ReHo between the two bands were found in fusiform gyrus and superior frontal gyrus (slow-4> slow-5), and in basal ganglia, parahippocampus, and dorsal middle prefrontal gyrus (slow-5> slow-4). Importantly, we identified significant interaction between frequency bands and groups in the inferior occipital gyrus and caudate body. This study demonstrates that ReHo changes in schizophrenia are widespread and frequency dependent.     

Effective Connectivity of Depth-Structure–Selective Patches in the Lateral Bank of the Macaque Intraparietal Sulcus

Extrastriate cortical areas are frequently composed of subpopulations of neurons encoding specific features or stimuli, such as color, disparity, or faces, and patches of neurons encoding similar stimulus properties are typically embedded in interconnected networks, such as the attention or face-processing network. The goal of the current study was to examine the effective connectivity of subsectors of neurons in the same cortical area with highly similar neuronal response properties. We first recorded single- and multi-unit activity to identify two neuronal patches in the anterior part of the macaque intraparietal sulcus (IPS) showing the same depth structure selectivity and then employed electrical microstimulation during functional magnetic resonance imaging in these patches to determine the effective connectivity of these patches. The two IPS subsectors we identified—with the same neuronal response properties and in some cases separated by only 3 mm—were effectively connected to remarkably distinct cortical networks in both dorsal and ventral stream in three macaques. Conversely, the differences in effective connectivity could account for the known visual-to-motor gradient within the anterior IPS. These results clarify the role of the anterior IPS as a pivotal brain region where dorsal and ventral visual stream interact during object analysis. Thus, in addition to the anatomical connectivity of cortical areas and the properties of individual neurons in these areas, the effective connectivity provides novel key insights into the widespread functional networks that support behavior.     

Altered Resting State Brain Networks in Parkinson’s Disease

Parkinson’s disease (PD) is a neurodegenerative disorder affecting dopaminergic neurons in the substantia nigra leading to dysfunctional cortico-striato-thalamic-cortical loops. In addition to the characteristic motor symptoms, PD patients often show cognitive impairments, affective changes and other non-motor symptoms, suggesting system-wide effects on brain function. Here, we used functional magnetic resonance imaging and graph-theory based analysis methods to investigate altered whole-brain intrinsic functional connectivity in PD patients (n = 37) compared to healthy controls (n = 20). Global network properties indicated less efficient processing in PD. Analysis of brain network modules pointed to increased connectivity within the sensorimotor network, but decreased interaction of the visual network with other brain modules. We found lower connectivity mainly between the cuneus and the ventral caudate, medial orbitofrontal cortex and the temporal lobe. To identify regions of altered connectivity, we mapped the degree of intrinsic functional connectivity both on ROI- and on voxel-level across the brain. Compared to healthy controls, PD patients showed lower connectedness in the medial and middle orbitofrontal cortex. The degree of connectivity was also decreased in the occipital lobe (cuneus and calcarine), but increased in the superior parietal cortex, posterior cingulate gyrus, supramarginal gyrus and supplementary motor area. Our results on global network and module properties indicated that PD manifests as a disconnection syndrome. This was most apparent in the visual network module. The higher connectedness within the sensorimotor module in PD patients may be related to compensation mechanism in order to overcome the functional deficit of the striato-cortical motor loops or to loss of mutual inhibition between brain networks. Abnormal connectivity in the visual network may be related to adaptation and compensation processes as a consequence of altered motor function. Our analysis approach proved sensitive for detecting disease-related localized effects as well as changes in network functions on intermediate and global scale.     

Physiological and pathological oscillatory networks in the human motor system.

Human brain functions are heavily contingent on neural interactions both at the single neuron and the neural population or system level. Accumulating evidence from neurophysiological studies strongly suggests that coupling of oscillatory neural activity provides an important mechanism to establish neural interactions. With the availability of whole-head magnetoencephalography (MEG) macroscopic oscillatory activity can be measured non-invasively from the human brain with high temporal and spatial resolution. To localise, quantify and map oscillatory activity and interactions onto individual brain anatomy we have developed the 'dynamic imaging of coherent sources' (DICS) method which allows to identify and analyse cerebral oscillatory networks from MEG recordings. Using this approach we have characterized physiological and pathological oscillatory networks in the human sensorimotor system. Coherent 8 Hz oscillations emerge from a cerebello-thalamo-premotor-motor cortical network and exert an 8 Hz oscillatory drive on the spinal motor neurons which can be observed as a physiological tremulousness of the movement termed movement discontinuities. This network represents the neurophysiological substrate of a discrete mode of motor control. In parkinsonian resting tremor we have identified an extensive cerebral network consisting of primary motor and lateral premotor cortex, supplementary motor cortex, thalamus/basal ganglia, posterior parietal cortex and secondary somatosensory cortex, which are entrained in the tremor or twice the tremor rhythm. This low frequency entrapment of motor areas likely plays an important role in the pathophysiology of parkinsonian motor symptoms. Finally, studies on patients with postural tremor in hepatic encephalopathy revealed that this type of tremor results from a pathologically slow thalamocortical and cortico-muscular coupling during isometric hold tasks. In conclusion, the analysis of oscillatory cerebral networks provides new insights into physiological mechanisms of motor control and pathophysiological mechanisms of tremor disorders.     

Development of Large-Scale Functional Brain Networks in Children

The ontogeny of large-scale functional organization of the human brain is not well understood. Here we use network analysis of intrinsic functional connectivity to characterize the organization of brain networks in 23 children (ages 7–9 y) and 22 young-adults (ages 19–22 y). Comparison of network properties, including path-length, clustering-coefficient, hierarchy, and regional connectivity, revealed that although children and young-adults' brains have similar “small-world” organization at the global level, they differ significantly in hierarchical organization and interregional connectivity. We found that subcortical areas were more strongly connected with primary sensory, association, and paralimbic areas in children, whereas young-adults showed stronger cortico-cortical connectivity between paralimbic, limbic, and association areas. Further, combined analysis of functional connectivity with wiring distance measures derived from white-matter fiber tracking revealed that the development of large-scale brain networks is characterized by weakening of short-range functional connectivity and strengthening of long-range functional connectivity. Importantly, our findings show that the dynamic process of over-connectivity followed by pruning, which rewires connectivity at the neuronal level, also operates at the systems level, helping to reconfigure and rebalance subcortical and paralimbic connectivity in the developing brain. Our study demonstrates the usefulness of network analysis of brain connectivity to elucidate key principles underlying functional brain maturation, paving the way for novel studies of disrupted brain connectivity in neurodevelopmental disorders such as autism.     

Time-Dependent Increase in Network Response to Stimulation

In vitro neuronal cultures have become a popular method with which to probe network-level neuronal dynamics and phenomena in controlled laboratory settings. One of the key dynamics of interest in these in vitro studies has been the extent to which cultured networks display properties indicative of learning. Here we demonstrate the effects of a high frequency electrical stimulation signal in training cultured networks of cortical neurons. Networks receiving this training signal displayed a time-dependent increase in the response to a low frequency probing stimulation, particularly in the time window of 20–50 ms after stimulation. This increase was found to be statistically significant as compared to control networks that did not receive training. The timing of this increase suggests potentiation of synaptic mechanisms. To further investigate this possibility, we leveraged the powerful Cox statistical connectivity method as previously investigated by our group. This method was used to identify and track changes in network connectivity strength.     

Frequency-Dependent Modulation of Regional Synchrony in the Human Brain by Eyes Open and Eyes Closed Resting-States

The eyes-open (EO) and eyes-closed (EC) states have differential effects on BOLD-fMRI signal dynamics, affecting both the BOLD oscillation frequency of a single voxel and the regional homogeneity (ReHo) of several neighboring voxels. To explore how the two resting-states modulate the local synchrony through different frequency bands, we decomposed the time series of each voxel into several components that fell into distinct frequency bands. The ReHo in each of the bands was calculated and compared between the EO and EC conditions. The cross-voxel correlations between the mean frequency and the overall ReHo of each voxel’s original BOLD series in different brain areas were also calculated and compared between the two states. Compared with the EC state, ReHo decreased with EO in a wide frequency band of 0.01–0.25 Hz in the bilateral thalamus, sensorimotor network, and superior temporal gyrus, while ReHo increased significantly in the band of 0–0.01 Hz in the primary visual cortex, and in a higher frequency band of 0.02–0.1 Hz in the higher order visual areas. The cross-voxel correlations between the frequency and overall ReHo were negative in all the brain areas but varied from region to region. These correlations were stronger with EO in the visual network and the default mode network. Our results suggested that different frequency bands of ReHo showed different sensitivity to the modulation of EO-EC states. The better spatial consistency between the frequency and overall ReHo maps indicated that the brain might adopt a stricter frequency-dependent configuration with EO than with EC.     

A Statistical Method of Identifying Interactions in Neuron–Glia Systems Based on Functional Multicell Ca2+ Imaging

Crosstalk between neurons and glia may constitute a significant part of information processing in the brain. We present a novel method of statistically identifying interactions in a neuron–glia network. We attempted to identify neuron–glia interactions from neuronal and glial activities via maximum-a-posteriori (MAP)-based parameter estimation by developing a generalized linear model (GLM) of a neuron–glia network. The interactions in our interest included functional connectivity and response functions. We evaluated the cross-validated likelihood of GLMs that resulted from the addition or removal of connections to confirm the existence of specific neuron-to-glia or glia-to-neuron connections. We only accepted addition or removal when the modification improved the cross-validated likelihood. We applied the method to a high-throughput, multicellular in vitro Ca2+ imaging dataset obtained from the CA3 region of a rat hippocampus, and then evaluated the reliability of connectivity estimates using a statistical test based on a surrogate method. Our findings based on the estimated connectivity were in good agreement with currently available physiological knowledge, suggesting our method can elucidate undiscovered functions of neuron–glia systems.     

Visual working memory load-related changes in neural activity and functional connectivity

Background

Visual working memory (VWM) helps us store visual information to prepare for subsequent behavior. The neuronal mechanisms for sustaining coherent visual information and the mechanisms for limited VWM capacity have remained uncharacterized. Although numerous studies have utilized behavioral accuracy, neural activity, and connectivity to explore the mechanism of VWM retention, little is known about the load-related changes in functional connectivity for hemi-field VWM retention.

Methodology/Principal Findings

In this study, we recorded electroencephalography (EEG) from 14 normal young adults while they performed a bilateral visual field memory task. Subjects had more rapid and accurate responses to the left visual field (LVF) memory condition. The difference in mean amplitude between the ipsilateral and contralateral event-related potential (ERP) at parietal-occipital electrodes in retention interval period was obtained with six different memory loads. Functional connectivity between 128 scalp regions was measured by EEG phase synchronization in the theta- (4–8 Hz), alpha- (8–12 Hz), beta- (12–32 Hz), and gamma- (32–40 Hz) frequency bands. The resulting matrices were converted to graphs, and mean degree, clustering coefficient and shortest path length was computed as a function of memory load. The results showed that brain networks of theta-, alpha-, beta-, and gamma- frequency bands were load-dependent and visual-field dependent. The networks of theta- and alpha- bands phase synchrony were most predominant in retention period for right visual field (RVF) WM than for LVF WM. Furthermore, only for RVF memory condition, brain network density of theta-band during the retention interval were linked to the delay of behavior reaction time, and the topological property of alpha-band network was negative correlation with behavior accuracy.

Conclusions/Significance

We suggest that the differences in theta- and alpha- bands between LVF and RVF conditions in functional connectivity and topological properties during retention period may result in the decline of behavioral performance in RVF task.     

Bursting Reverberation as a Multiscale Neuronal Network Process Driven by Synaptic Depression-Facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5–30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2–5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.     

Determination of Vascular Dementia Brain in Distinct Frequency Bands with Whole Brain Functional Connectivity Patterns

Recent studies have shown that multivariate pattern analysis (MVPA) can be useful for distinguishing brain disorders into categories. Such analyses can substantially enrich and facilitate clinical diagnoses. Using MPVA methods, whole brain functional networks, especially those derived using different frequency windows, can be applied to detect brain states. We constructed whole brain functional networks for groups of vascular dementia (VaD) patients and controls using resting state BOLD-fMRI (rsfMRI) data from three frequency bands - slow-5 (0.01∼0.027 Hz), slow-4 (0.027∼0.073 Hz), and whole-band (0.01∼0.073 Hz). Then we used the support vector machine (SVM), a type of MVPA classifier, to determine the patterns of functional connectivity. Our results showed that the brain functional networks derived from rsfMRI data (19 VaD patients and 20 controls) in these three frequency bands appear to reflect neurobiological changes in VaD patients. Such differences could be used to differentiate the brain states of VaD patients from those of healthy individuals. We also found that the functional connectivity patterns of the human brain in the three frequency bands differed, as did their ability to differentiate brain states. Specifically, the ability of the functional connectivity pattern to differentiate VaD brains from healthy ones was more efficient in the slow-5 (0.01∼0.027 Hz) band than in the other two frequency bands. Our findings suggest that the MVPA approach could be used to detect abnormalities in the functional connectivity of VaD patients in distinct frequency bands. Identifying such abnormalities may contribute to our understanding of the pathogenesis of VaD.