静息状态下脑功能连接的磁共振成像研究

静息状态下脑功能连接的磁共振成像研究近年来取得了迅猛发展. 通过对fMRI信号低频涨落成分的同步性分析,可以得到大脑静息态任意脑区的功能连接和多套网络系统,其中“默认网络”的发现可能为人脑固有网络的研究提供新的思路. 而静息态网络与解剖连接之间可能存在的对应,以及在神经精神疾病患者脑中性质和连接的异常改变,使其具有重要的研究和临床应用价值. 该文总结了静息状态功能磁共振成像的主要研究成果,对静息状态脑功能网络的发现和发展、研究方法、各网络及其特点以及在临床方面的应用进行简单的介绍和分析.     

An algorithm to estimate anatomical connectivity between brain regions using diffusion MRI

The study of anatomical connectivity is essential for interpreting functional MRI data and for establishing how brain areas are linked together into networks to support higher-order functions. Diffusion-weighted MR images (DWI) and tractography provide a unique noninvasive tool to explore the connectional architecture of the brain. The identification of anatomical circuits associated with a specific function can be better accomplished by the joint application of diffusion and functional MRI. In this article, we propose a simple algorithm to identify the set of pathways between two regions of interest. The method is based upon running deterministic tractography from all possible starting positions in the brain and selecting trajectories that intersect both regions. We compare results from single-fiber tractography using diffusion tensor imaging and from multi-fiber tractography using reduced-encoding persistent angular structure (PAS) MRI on standard DWI datasets from healthy human volunteers. Our results show that, in comparison with single-fiber tractography, the multi-fiber technique reveals additional putative routes of connection. We demonstrate highly consistent results of the proposed technique over a cohort of 16 healthy subjects.     

Use of random matrix theory in the discovery of resting state brain networks

Connectomics identifies brain networks in vivo in resting state functional MRI. However, the presence of noise produces spurious identification of brain networks, which have low test-retest reliability. A Network Based Statistics approach to network identification has been previously proposed that affords much better statistical power relative to Bonferroni method but nevertheless provides a sufficiently conservative, family-wise control for false positives. We propose the use of Random Matrix Theory (RMT) to discover brain networks and to associate those networks with demographic and clinical variables. We parcellated the brain into cortical and subcortical regions using either an anatomical or a functional brain atlas. We applied RMT to study functional connectivity across brain regions by first computing the correlation matrix for time courses in those brain regions and then identifying eigenvalues that deviate from the theoretical random distribution that RMT predicts, on the assumption that real brain networks would produce eigenvalues that differ significantly from the random distribution. We assessed the specificity and test-retest reliability of identified networks through application of this RMT-based approach to (1) synthetic data generated under the null-hypothesis, (2) resting state functional MRI data from 4 real-world cohorts of patients and healthy controls, and (3) synthetic data generated by the addition of increasing amounts of noise to real-world datasets. Our findings showed that RMT method was robust to the atlas used for parcellating the brain and did not discover a brain network in synthetic data when in fact a network was not present (i.e., specificity was high); RMT-identified networks in the real-world dataset had high test-retest reliability; and RMT-based method consistently discovered the same network in the presence of increasing noise in the real-world dataset.     

The role of nonlinearity in computing graph-theoretical properties of resting-state functional magnetic resonance imaging brain networks

In recent years, there has been an increasing interest in the study of large-scale brain activity interaction structure from the perspective of complex networks, based on functional magnetic resonance imaging (fMRI) measurements. To assess the strength of interaction (functional connectivity, FC) between two brain regions, the linear (Pearson) correlation coefficient of the respective time series is most commonly used. Since a potential use of nonlinear FC measures has recently been discussed in this and other fields, the question arises whether particular nonlinear FC measures would be more informative for the graph analysis than linear ones. We present a comparison of network analysis results obtained from the brain connectivity graphs capturing either full (both linear and nonlinear) or only linear connectivity using 24 sessions of human resting-state fMRI. For each session, a matrix of full connectivity between 90 anatomical parcel time series is computed using mutual information. For comparison, connectivity matrices obtained for multivariate linear Gaussian surrogate data that preserve the correlations, but remove any nonlinearity are generated. Binarizing these matrices using multiple thresholds, we generate graphs corresponding to linear and full nonlinear interaction structures. The effect of neglecting nonlinearity is then assessed by comparing the values of a range of graph-theoretical measures evaluated for both types of graphs. Statistical comparisons suggest a potential effect of nonlinearity on the local measures-clustering coefficient and betweenness centrality. Nevertheless, subsequent quantitative comparison shows that the nonlinearity effect is practically negligible when compared to the intersubject variability of the graph measures. Further, on the group-average graph level, the nonlinearity effect is unnoticeable.     

Anatomical correlates of the functional organization in the human occipitotemporal cortex

The connectivity between functionally distinct areas in the human brain is unknown because of the limitations posed by current postmortem anatomical labeling techniques. Diffusion tensor imaging (DTI) has previously been used to define large white matter tracts based on well-known anatomical landmarks in the living human brain. In the present study, we used DTI coupled with functional magnetic resonance imaging (fMRI) to assess neuronal connections between human striate and functionally defined extrastriate ventral cortical areas. Functional areas were identified with conventional fMRI mapping procedures and then used as seeding points in a DTI analysis to ascertain connectivity patterns between cortical areas, thus yielding the pattern of connections between human occipitoventral visual areas in vivo.     

Fine-scale functional connectivity in somatosensory cortex revealed by high-resolution fMRI

High-resolution functional magnetic resonance imaging (fMRI) at high field (9.4 T) has been used to measure functional connectivity between subregions within the primary somatosensory (SI) cortex of the squirrel monkey brain. The hand-face region within the SI cortex of the squirrel monkey has been previously well mapped with functional imaging and electrophysiological and anatomical methods, and the orderly topographic map of the hand region is characterized by a lateral to medial representation of individual digits in four subregions of areas 3a, 3b, 1 and 2. With submillimeter resolution, we are able to detect not only the separate islands of activation corresponding to vibrotactile stimulations of single digits but also, in subsequent acquisitions, the degree of correlation between voxels within the SI cortex in the resting state. The results suggest that connectivity patterns are very similar to stimulus-driven distributions of activity and that connectivity varies on the scale of millimeters within the same primary region. Connectivity strength is not a reflection of global larger-scale changes in blood flow and is not directly dependent on distance between regions. Preliminary electrophysiological recordings agree well with the fMRI data. In human studies at 7 T, high-resolution fMRI may also be used to identify the same subregions and assess responses to sensory as well as painful stimuli, and to measure connectivity dynamically before and after such stimulations.     

Networks behind the morphology and structural design of living systems

Technological advances in imaging techniques and biometric data acquisition have enabled us to apply methods of network science to study the morphology and structural design of organelles, organs, and tissues, as well as the coordinated interactions among them that yield a healthy physiology at the level of whole organisms. We here review research dedicated to these advances, in particular focusing on networks between cells, the topology of multicellular structures, neural interactions, fluid transportation networks, and anatomical networks. The percolation of blood vessels, structural connectivity within the brain, the porous structure of bones, and relations between different anatomical parts of the human body are just some of the examples that we explore in detail. We argue and show that the models, methods, and algorithms developed in the realm of network science are ushering in a new era of network-based inquiry into the morphology and structural design of living systems in the broadest possible terms. We also emphasize that the need and applicability of this research is likely to increase significantly in the years to come due to the rapid progress made in the development of bioartificial substitutes and tissue engineering.     

Combining independent component analysis and correlation analysis to probe interregional connectivity in fMRI task activation datasets

A new approach in studying interregional functional connectivity using functional magnetic resonance imaging (fMRI) is presented. Functional connectivity may be detected by means of cross correlating time course data from functionally related brain regions. These data exhibit high temporal coherence of low frequency fluctuations due to synchronized blood flow changes. In the past, this fMRI technique for studying functional connectivity has been applied to subjects that performed no prescribed task ("resting" state). This paper presents the results of applying the same method to task-related activation datasets. Functional connectivity analysis is first performed in areas not involved with the task. Then a method is devised to remove the effects of activation from the data using independent component analysis (ICA) and functional connectivity analysis is repeated. Functional connectivity, which is demonstrated in the "resting brain," is not affected by tasks which activate unrelated brain regions. In addition, ICA effectively removes activation from the data and may allow us to study functional connectivity even in the activated regions.     

Inferring Excitatory and Inhibitory Connections in Neuronal Networks

The comprehension of neuronal network functioning, from most basic mechanisms of signal transmission to complex patterns of memory and decision making, is at the basis of the modern research in experimental and computational neurophysiology. While mechanistic knowledge of neurons and synapses structure increased, the study of functional and effective networks is more complex, involving emergent phenomena, nonlinear responses, collective waves, correlation and causal interactions. Refined data analysis may help in inferring functional/effective interactions and connectivity from neuronal activity. The Transfer Entropy (TE) technique is, among other things, well suited to predict structural interactions between neurons, and to infer both effective and structural connectivity in small- and large-scale networks. To efficiently disentangle the excitatory and inhibitory neural activities, in the article we present a revised version of TE, split in two contributions and characterized by a suited delay time. The method is tested on in silico small neuronal networks, built to simulate the calcium activity as measured via calcium imaging in two-dimensional neuronal cultures. The inhibitory connections are well characterized, still preserving a high accuracy for excitatory connections prediction. The method could be applied to study effective and structural interactions in systems of excitable cells, both in physiological and in pathological conditions.     

Biophysical and neural basis of resting state functional connectivity: Evidence from non-human primates

Functional MRI (fMRI) has evolved from simple observations of regional changes in MRI signals caused by cortical activity induced by a task or stimulus, to task-free acquisitions of images in a resting state. Such resting state signals contain low frequency fluctuations which may be correlated between voxels, and strongly correlated regions are deemed to reflect functional connectivity within synchronized circuits. Resting state functional connectivity (rsFC) measures have been widely adopted by the neuroscience community, and are being used and interpreted as indicators of intrinsic neural circuits and their functional states in a broad range of applications, both basic and clinical. However, there has been relatively little work reported that validates whether inter-regional correlations in resting state fluctuations of fMRI (rsfMRI) signals actually measure functional connectivity between brain regions, or to establish how MRI data correlate with other metrics of functional connectivity. In this mini-review, we summarize recent studies of rsFC within mesoscopic scale cortical networks (100 μm–10 mm) within a well defined functional region of primary somatosensory cortex (S1), as well as spinal cord and brain white matter in non-human primates, in which we have measured spatial patterns of resting state correlations and validated their interpretation with electrophysiological signals and anatomic connections. Moreover, we emphasize that low frequency correlations are a general feature of neural systems, as evidenced by their presence in the spinal cord as well as white matter. These studies demonstrate the valuable role of high field MRI and invasive measurements in an animal model to inform the interpretation of human imaging studies.     

基于持久同调的复杂脑网络动态演化分析

有研究表明阿尔茨海默病（Alzheimer's disease,AD）的认知状态与动态功能连接时间特性的改变有关,持久同调指标分析方法可为AD动态脑网络的研究提供更深的见解,但是目前研究主要集中在空间演化方面,尚未有针对时变方面的脑网络演化研究．本文基于静息态功能磁共振成像（resting state-functional magnetic resonanceimaging,rs-fMRI）,对AD患者和正常被试（normal controls,NC）的静态脑网络和基于滑动窗口构建的动态脑网络进行功能连接性分析．对基于持久同调和基于图论的分析结果进行了比较,并采用k均值聚类进行了时间属性的分析．结果表明相对图论指标,持久同调的指标在AD患者和NC被试间具有更显著的差异性;而且相对于静态脑网络,基于持久同调的动态脑网络演化分析可为脑功能网络标志物的检测提供新思路．     

Improving the resolution of functional brain imaging: analyzing functional data in anatomical space

The accurate mapping of functional magnetic resonance imaging (fMRI) activations to anatomical structures is critical for fMRI studies of brain organization. In the commonly used functional space analysis method, functional images are realigned to a functional reference image and processed in low-resolution functional space. The average functional activations are then projected into high-resolution anatomical space for visualization. Here, we describe a new technique, anatomical space analysis (ASA), whereby low-resolution functional images are first coregistered and resampled directly into high-resolution anatomical space with all subsequent data processing performed in high-resolution space. A major advantage of ASA is that minor scanner sampling instabilities and small head movements can increase spatial resolution by providing multiple samples of the relationship between functional and anatomical space. Both simulations and analyses of real fMRI data show that ASA improves the precision, objectivity and reproducibility of functional brain mapping.     

Abnormal hubs of white matter networks in the frontal-parieto circuit contribute to depression discrimination via pattern classification

Previous studies had explored the diagnostic and prognostic value of the structural neuroimaging data of MDD and treated the whole brain voxels, the fractional anisotropy and the structural connectivity as classification features. To our best knowledge, no study examined the potential diagnostic value of the hubs of anatomical brain networks in MDD. The purpose of the current study was to provide an exploratory examination of the potential diagnostic and prognostic values of hubs of white matter brain networks in MDD discrimination and the corresponding impaired hub pattern via a multi-pattern analysis. We constructed white matter brain networks from 29 depressions and 30 healthy controls based on diffusion tensor imaging data, calculated nodal measures and identified hubs. Using these measures as features, two types of feature architectures were established, one only included hubs (HUB) and the other contained both hubs and non hubs. The support vector machine classifiers with Gaussian radial basis kernel were used after the feature selection. Moreover, the relative contribution of the features was estimated by means of the consensus features. Our results presented that the hubs (including the bilateral dorsolateral part of superior frontal gyrus, the left middle frontal gyrus, the bilateral middle temporal gyrus, and the bilateral inferior temporal gyrus) played an important role in distinguishing the depressions from healthy controls with the best accuracy of 83.05%. Moreover, most of the HUB consensus features located in the frontal-parieto circuit. These findings provided evidence that the hubs could be served as valuable potential diagnostic measure for MDD, and the hub-concentrated lesion distribution of MDD was primarily anchored within the frontal-parieto circuit.     

Fast reproducible identification and large-scale databasing of individual functional cognitive networks

Background  

Although cognitive processes such as reading and calculation are associated with reproducible cerebral networks, inter-individual variability is considerable. Understanding the origins of this variability will require the elaboration of large multimodal databases compiling behavioral, anatomical, genetic and functional neuroimaging data over hundreds of subjects. With this goal in mind, we designed a simple and fast acquisition procedure based on a 5-minute functional magnetic resonance imaging (fMRI) sequence that can be run as easily and as systematically as an anatomical scan, and is therefore used in every subject undergoing fMRI in our laboratory. This protocol captures the cerebral bases of auditory and visual perception, motor actions, reading, language comprehension and mental calculation at an individual level.     

Epidemic outbreaks in complex heterogeneous networks

We present a detailed analytical and numerical study for the spreading of infections with acquired immunity in complex population networks. We show that the large connectivity fluctuations usually found in these networks strengthen considerably the incidence of epidemic outbreaks. Scale-free networks, which are characterized by diverging connectivity fluctuations in the limit of a very large number of nodes, exhibit the lack of an epidemic threshold and always show a finite fraction of infected individuals. This particular weakness, observed also in models without immunity, defines a new epidemiological framework characterized by a highly heterogeneous response of the system to the introduction of infected individuals with different connectivity. The understanding of epidemics in complex networks might deliver new insights in the spread of information and diseases in biological and technological networks that often appear to be characterized by complex heterogeneous architectures. Received 20 September 2001 and Received in final form 4 February 2002     

A conditional Granger causality model approach for group analysis in functional magnetic resonance imaging

Granger causality model (GCM) derived from multivariate vector autoregressive models of data has been employed to identify effective connectivity in the human brain with functional magnetic resonance imaging (fMRI) and to reveal complex temporal and spatial dynamics underlying a variety of cognitive processes. In the most recent fMRI effective connectivity measures, pair-wise GCM has commonly been applied based on single-voxel values or average values from special brain areas at the group level. Although a few novel conditional GCM methods have been proposed to quantify the connections between brain areas, our study is the first to propose a viable standardized approach for group analysis of fMRI data with GCM. To compare the effectiveness of our approach with traditional pair-wise GCM models, we applied a well-established conditional GCM to preselected time series of brain regions resulting from general linear model (GLM) and group spatial kernel independent component analysis of an fMRI data set in the temporal domain. Data sets consisting of one task-related and one resting-state fMRI were used to investigate connections among brain areas with the conditional GCM method. With the GLM-detected brain activation regions in the emotion-related cortex during the block design paradigm, the conditional GCM method was proposed to study the causality of the habituation between the left amygdala and pregenual cingulate cortex during emotion processing. For the resting-state data set, it is possible to calculate not only the effective connectivity between networks but also the heterogeneity within a single network. Our results have further shown a particular interacting pattern of default mode network that can be characterized as both afferent and efferent influences on the medial prefrontal cortex and posterior cingulate cortex. These results suggest that the conditional GCM approach based on a linear multivariate vector autoregressive model can achieve greater accuracy in detecting network connectivity than the widely used pair-wise GCM, and this group analysis methodology can be quite useful to extend the information obtainable in fMRI.     

Natural Connectivity of Complex Networks

The concept of natural connectivity is reported as a robustness measure of complex networks. The natural connectivity has a clear physical meaning and a simple mathematical formulation. It is shown that the natural connectivity can be derived mathematically from the graph spectrum as an average eigenvalue and that it changes strictly monotonically with the addition or deletion of edges. By comparing the natural connectivity with other typical robustness measures within a scenario of edge elimination, it is demonstrated that the naturM connectivity has an acute discrimination which agrees with our intuition.     

The impact of sampling density upon cortical network analysis: regions or points

The choice of representation has a fundamental influence on the network analysis results of an empirical data set. The answers to two basic questions - how to define a node and how to define an edge between a pair of nodes - are not obvious in the network analysis of brain imaging data. We considered the first question in the case of magnetic resonance imaging (MRI)-based cortical thickness networks. We selected network nodes to represent vertices of a cortical surface mesh or cortical brain regions. The first network represents the maximal level of detail available in the analysis of cortical thickness networks, while the latter network represents the typical level of detail in the current network analysis studies. We compared the network analysis results between these two representations. The basic network measures behaved approximately as expected when the level of detail increased. However, the overall connectivity of nodes was greater in the vertex level, degree of clustering was smaller in the vertex level, and the node centralities were different between the levels. Further, many parameters of vertex-level network were more robust to the selection of the correlation threshold used to define the edges of network. We conclude that albeit many qualitative network properties were consistent between the two resolution levels, the vertex-level resolution revealed details that were not visible at the regional-level networks, and this additional detail could be useful for some applications. Finally, a similar methodology as the one used here could be used to study effects of the sampling density in other brain-imaging-based networks, for example, in resting-state functional MRI.     

Detecting functional connectivity in the resting brain: a comparison between ICA and CCA

Independent component analysis (ICA) and cross-correlation analysis (CCA) are general tools for detecting resting-state functional connectivity. In this study, we jointly evaluated these two approaches based on simulated data and in vivo functional magnetic resonance imaging data acquired from 10 resting healthy subjects. The influence of the number of independent components (maps) on the results of ICA was investigated. The influence of the selection of the seeds on the results of CCA was also examined. Our results reveal that significant differences between these two approaches exist. The performance of ICA is superior as compared with that of CCA; in addition, the performance of ICA is not significantly affected by structured noise over a relatively large range. The results of ICA could be affected by the number of independent components if this number is too small, however. Converting the spatially independent maps of ICA into z maps for thresholding tends to overestimate the false-positive rate. However, the overestimation is not very severe and may be acceptable in most cases. The results of CCA are dependent on seeds location. Seeds selected based on different criteria will significantly affect connectivity maps.     

Estimation of the effective and functional human cortical connectivity with structural equation modeling and directed transfer function applied to high-resolution EEG

Different brain imaging devices are presently available to provide images of the human functional cortical activity, based on hemodynamic, metabolic or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions are interconnected. The concept of brain connectivity plays a central role in the neuroscience, and different definitions of connectivity, functional and effective, have been adopted in literature. While the functional connectivity is defined as the temporal coherence among the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally among cortical sites. The structural equation modeling (SEM) is the most used method to estimate effective connectivity in neuroscience, and its typical application is on data related to brain hemodynamic behavior tested by functional magnetic resonance imaging (fMRI), whereas the directed transfer function (DTF) method is a frequency-domain approach based on both a multivariate autoregressive (MVAR) modeling of time series and on the concept of Granger causality.

This study presents advanced methods for the estimation of cortical connectivity by applying SEM and DTF on the cortical signals estimated from high-resolution electroencephalography (EEG) recordings, since these signals exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. To estimate correctly the cortical signals, we used a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual MRI, a distributed source model and a regularized linear inverse source estimates of cortical current density. Before the application of SEM and DTF methodology to the cortical waveforms estimated from high-resolution EEG data, we performed a simulation study, in which different main factors (signal-to-noise ratio, SNR, and simulated cortical activity duration, LENGTH) were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). The statistical analysis returned that during simulations, both SEM and DTF estimators were able to correctly estimate the imposed connectivity patterns under reasonable operative conditions, that is, when data exhibit an SNR of at least 3 and a LENGTH of at least 75 s of nonconsecutive EEG recordings at 64 Hz of sampling rate.

Hence, effective and functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high-resolution EEG techniques and linear inverse estimation with SEM or DTF methods. We conclude that the estimation of cortical connectivity can be performed not only with hemodynamic measurements, but also with EEG signals treated with advanced computational techniques.