Assessing functional connectivity in the human brain by fMRI

Functional magnetic resonance imaging (fMRI) is widely used to detect and delineate regions of the brain that change their level of activation in response to specific stimuli and tasks. Simple activation maps depict only the average level of engagement of different regions within distributed systems. FMRI potentially can reveal additional information about the degree to which components of large-scale neural systems are functionally coupled together to achieve specific tasks. In order to better understand how brain regions contribute to functionally connected circuits, it is necessary to record activation maps either as a function of different conditions, at different times or in different subjects. Data obtained under different conditions may then be analyzed by a variety of techniques to infer correlations and couplings between nodes in networks. Several multivariate statistical methods have been adapted and applied to analyze variations within such data. An approach of particular interest that is suited to studies of connectivity within single subjects makes use of acquisitions of runs of MRI images obtained while the brain is in a so-called steady state, either at rest (i.e., without any specific stimulus or task) or in a condition of continuous activation. Interregional correlations between fluctuations of MRI signal potentially reveal functional connectivity. Recent studies have established that interregional correlations between different components of circuits in each of the visual, language, motor and working memory systems can be detected in the resting state. Correlations at baseline are changed during the performance of a continuous task. In this review, various methods available for assessing connectivity are described and evaluated.     

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.     

A serial functional connectivity MRI study in healthy individuals assessing the variability of connectivity measures: reduced interhemispheric connectivity in the motor network during continuous performance

To date, little data is available on the reproducibility of functional connectivity MRI (fcMRI) studies. Here, we tested the variability and reproducibility of both the functional connectivity itself and different statistical methods to analyze this phenomenon. In the main part of our study, we repeatedly examined two healthy subjects in 10 sessions over 6 months with fcMRI. Cortical areas involved in motor function were examined under two different cognitive states: during continuous performance (CP) of a flexion/extension task of the fingers of the right hand and while subjects were at rest. Connectivity to left primary motor cortex (lSM1) was calculated by correlation analysis. The resulting correlation coefficients were transformed to z-scores of the standard normal distribution. For each subject, multisession statistical analyses were carried out with the z-score maps of the resting state (RS) and the CP experiments. First, voxel based t tests between the two groups of fcMRI experiments were performed. Second, ROI analyses were carried out for contralateral right SM1 and for supplementary motor area (SMA). For both ROI, mean and maximum z-score were calculated for each experiment. Also, the fraction of significantly (P<.05) correlated voxels (FCV) in each ROI was calculated. To evaluate the differences between the RS and the CP condition, paired t tests were performed for the mean and maximum z-scores, and Wilcoxon signed ranks tests for matched pairs were carried out for the FCV. All statistical methods and connectivity measures under investigation yielded a distinct loss in left–right SM1 connectivity under the CP condition. For SMA, interindividual differences were apparent. We therefore repeated the fcMRI experiments and the ROI analyses in a group of seven healthy subjects (including the two subjects of the main study). In this substudy, we were able to verify the reduction of left–right SM1 connectivity during unilateral performance. Still, the direction of SMA to lSM1 connectivity change during the CP condition remained undefined as four subjects showed a connectivity increase and three showed a decrease. In summary, we were able to demonstrate a distinct reduction in left–right SM1 synchrony in the CP condition compared to the RS both in the longitudinal and in the multisubject study. This effect was reproducible with all statistical methods and all measures of connectivity under investigation. We conclude that despite intra- and interindividual variability, serial and cross-sectional assessment of functional connectivity reveals stable and reliable results.     

Disrupted functional connectivity in social anxiety disorder: a resting-state fMRI study

Dysfunction of the corticolimbic circuitry has been highlighted in social anxiety disorder (SAD) during social stimuli. However, few studies have investigated functional connectivity in SAD during the resting state, which may improve our understanding of SAD pathophysiology. The aim of this study was to investigate whether whole-brain functional connectivity might be aberrant in SAD patients, and if so, whether these changes are related to the measured clinical severity. Seventeen SAD patients and 19 healthy controls participated in resting-state functional magnetic resonance imaging. The brain was first divided into 90 paired brain regions and functional connectivity was then estimated by temporal correlation between each of these regions. Furthermore, connections that were significantly disrupted in SAD patients were correlated with clinical severity measured using the Liebowitz Social Anxiety Scale. Compared with healthy controls, SAD patients showed decreased positive connections within the frontal lobe and decreased negative connections between the frontal and occipital lobes. In particular, the weaker negative connections between the frontal lobe, which mainly involved the right median prefrontal cortex, and the occipital lobe had a significant positive correlation with the severity of SAD symptoms. The results support the hypothesis that some abnormalities of functional connectivity exist in SAD patients, which relate to the frontal cortex and occipital cortex. In addition, decreased functional connectivity between the frontal and occipital lobes and within the frontal lobe might be related to abnormal information processing and reflect disturbed neural organization resulting in defective social cognition, which could represent an early imaging biomarker for SAD.     

Occipital-parietal interactions during shifts of exogenous visuospatial attention: trial-dependent changes of effective connectivity

We studied neural interactions between brain areas involved in exogenous (stimulus-driven) control of visuospatial attention. With event-related functional magnetic resonance imaging (fMRI), we investigated changes of connectivity during shifts of spatial attention from an attended location to a previously unattended target location. Using a 3-T scanner, fMRI data were acquired from three healthy volunteers. According to a central visual cue, participants directed endogenous spatial attention to the left or the right visual hemifield for blocks of 56 s. Peripheral visual targets were presented unpredictably in either the attended hemifield (valid trials, 80%) or in the unattended hemifield (invalid trials, 20%) and participants performed a two-alternative forced-choice discrimination task with the target, irrespective of cue validity. In accordance with previous results, we found that the temporal–parietal junction (TPJ) mediates the shift of spatial attention toward stimuli presented at the unattended side (i.e., invalid trials). We critically studied the interaction between occipital areas responding to the visual stimuli and other brain regions in order to find regions functionally coupled with the occipital cortex during invalid trials. We found that the coupling between occipital areas processing visual stimuli and the TPJ selectively increased during invalid trials. Our results highlight how changes of connectivity between brain areas can describe attentional processes such as stimulus-driven shifts of spatial attention.     

Fractional amplitude of low-frequency fluctuation changes in functional dyspepsia: A resting-state fMRI study

Recently, there is an increasing interest in the study of the role of brain dysfunction in the pathogenesis of symptoms of functional dyspepsia (FD). More specifically, abnormal brain activities in patients with FD during the resting state have been proven by several positron emission tomography (PET) studies. Resting-state functional magnetic resonance imaging (fMRI) is also a valuable tool in investigating spontaneous brain activity abnormalities in pathological conditions. In the present study, we examined the amplitude of low-frequency fluctuations (ALFF) and fractional (f)ALFF changes in patients with FD by using fMRI. Twenty-nine patients with FD and sixteen healthy controls participated in this study. Between-group differences in ALFF/fALFF were examined using a permutation-based nonparametric test after accounting for the gender and age effects. The results revealed a significant between-group difference in fALFF but not in ALFF in multiple brain regions including the right insula, brainstem and cerebellum. Seed-based resting-state functional connectivity analysis revealed that FD patients have increased correlations between the right cerebellum and multiple brain regions including the bilateral brainstem, bilateral cerebellum, bilateral thalamus, left para-/hippocampus, left pallidum and left putamen. Furthermore, fLAFF values in the right insula were positively correlated with the severity of the disease. These findings have provided further evidence of spontaneous brain activity abnormalities in FD patients which might contribute to our understanding of the pathophysiology of the disease.     

Electric stimulation fMRI of the perforant pathway to the rat hippocampus

The hippocampal formation is a brain system that is implicated in learning and memory. The major input to the hippocampus arrives from the entorhinal cortex (EC) to the dentate gyrus (DG) through the perforant path. In the present work, we have investigated the functional properties of this connection by concomitantly applying electrophysiological techniques, deep-brain electric microstimulation and functional magnetic resonance imaging in anesthetized rats. We systematically delivered different current intensities at diverse stimulation frequencies to the perforant path while recording electrophysiological and blood-oxygenation-level-dependent (BOLD) signals. We observed a linear relationship between the current intensity used to stimulate the hippocampal formation and the amplitude and extension of the induced BOLD response. In addition, we found a frequency-dependent spatial pattern of activation. With stimulation protocols and train frequencies used for kindling, the activity strongly spreads ipsilaterally through the hippocampus, DG, subiculum and EC.     

Spontaneous low-frequency blood oxygenation level-dependent fluctuations and functional connectivity analysis of the 'resting' brain

Functional magnetic resonance imaging techniques using the blood oxygenation level-dependent (BOLD) contrast are widely used to map human brain function by relating local hemodynamic responses to neuronal stimuli compared to control conditions. There is increasing interest in spontaneous cerebral BOLD fluctuations that are prominent in the low-frequency range (<0.1 Hz) and show intriguing spatio-temporal correlations in functional networks. The nature of these signal fluctuations remains unclear, but there is accumulating evidence for a neural basis opening exciting new avenues to study human brain function and its connectivity at rest. Moreover, an increasing number of patient studies report disease-dependent variation in the amplitude and spatial coherence of low-frequency BOLD fluctuations (LFBF) that may afford greater diagnostic sensitivity and easier clinical applicability than standard fMRI. The main disadvantage of this emerging tool relates to physiological (respiratory, cardiac and vasomotion) and motion confounds that are challenging to disentangle requiring thorough preprocessing. Technical aspects of functional connectivity fMRI analysis and the neuroscientific potential of spontaneous LFBF in the default mode and other resting-state networks have been recently reviewed. This review will give an update on the current knowledge of the nature of LFBF, their relation to physiological confounds and potential for clinical diagnostic and pharmacological studies.     

Images-based suppression of unwanted global signals in resting-state functional connectivity studies

Correlated fluctuations of low-frequency fMRI signal have been suggested to reflect functional connectivity among the involved regions. However, large-scale correlations are especially prone to spurious global modulations induced by coherent physiological noise. Cardiac and respiratory rhythms are the most offending component, and a tailored preprocessing is needed in order to reduce their impact. Several approaches have been proposed in the literature, generally based on the use of physiological recordings acquired during the functional scans, or on the extraction of the relevant information directly from the images. In this paper, the performances of the denoising approach based on general linear fitting of global signals of noninterest extracted from the functional scans were assessed. Results suggested that this approach is sufficiently accurate for the preprocessing of functional connectivity data.     

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.     

FMRI connectivity analysis of acupuncture effects on the whole brain network in mild cognitive impairment patients

The increased risk for the elderly with mild cognitive impairment (MCI) to progress to Alzheimer's disease makes it an appropriate condition for investigation. While the use of acupuncture as a complementary therapeutic method for treating MCI is popular in certain parts of the world, the underlying mechanism is still elusive. We sought to investigate the acupuncture effects on the functional connectivity throughout the entire brain in MCI patients compared to healthy controls (HC). The functional magnetic resonance imaging experiment was performed with two different paradigms, namely, deep acupuncture (DA) and superficial acupuncture (SA), at acupoint KI3. We first identified regions showing abnormal functional connectivity in the MCI group compared to HC during the resting state and subsequently tested whether these regions could be modulated by acupuncture. Then, we made the comparison of MCI vs. HC to test whether there were any specific modulatory patterns in the poststimulus resting brain between the two groups. Finally, we made the comparisons of DA vs. SA in each group to test the effect of acupuncture with different needling depths. We found the temporal regions (hippocampus, thalamus, fusiform gyrus) showing abnormal functional connectivity during the resting state. These regions are implicated in memory encoding and retrieving. Furthermore, we found significant changes in functional connectivity related with the abnormal regions in MCI patients following acupuncture. Compared to HC, the correlations related with the temporal regions were enhanced in the poststimulus resting brain in MCI patients. Compared to SA, significantly increased correlations related with the temporal regions were found for the DA condition. The enhanced correlations in the memory-related brain regions following acupuncture may be related to the purported therapeutically beneficial effects of acupuncture for the treatment of MCI. The heterogeneous modulatory patterns between DA and SA may suggest that deep muscle insertion of acupuncture is necessary to achieve the appreciable clinical effect.     

Functional connectivity in the rat brain: a complex network approach

Functional connectivity analyses of fMRI data can provide a wealth of information on the brain functional organization and have been widely applied to the study of the human brain. More recently, these methods have been extended to preclinical species, thus providing a powerful translational tool. Here, we review methods and findings of functional connectivity studies in the rat. More specifically, we focus on correlation analysis of pharmacological MRI (phMRI) responses, an approach that has enabled mapping the patterns of connectivity underlying major neurotransmitter systems in vivo. We also review the use of novel statistical approaches based on a network representation of the functional connectivity and their application to the study of the rat brain functional architecture.     

Real-time MR artifacts filtering during continuous EEG/fMRI acquisition

The purpose of this study was the development of a real-time filtering procedure of MRI artifacts in order to monitor the EEG activity during continuous EEG/fMRI acquisition. The development of a combined EEG and fMRI technique has increased in the past few years. Preliminary “spike-triggered” applications have been possible because in this method, EEG knowledge was only necessary to identify a trigger signal to start a delayed fMRI acquisition. In this way, the two methods were used together but in an interleaved manner. In real simultaneous applications, like event-related fMRI study, artifacts induced by MRI events on EEG traces represent a substantial obstacle for a right analysis. Up until now, the methods proposed to solve this problem are mainly based on procedures to remove post-processing artifacts without the possibility to control electrophysiological behavior of the patient during fMRI scan. Moreover, these methods are not characterized by a strong “prior knowledge” of the artifact, which is an imperative condition to avoid any loss of information on the physiological signals recovered after filtering. In this work, we present a new method to perform simultaneous EEG/fMRI study with real-time artifacts filtering characterized by a procedure based on a preliminary analytical study of EPI sequence parameters-related EEG-artifact shapes. Standard EEG equipment was modified in order to work properly during ultra-fast MRI acquisitions. Changes included: high-performance acquisition device; electrodes/cap/wires/cables materials and geometric design; shielding box for EEG signal receiver; optical fiber link; and software. The effects of the RF pulse and time-varying magnetic fields were minimized by using a correct head cap wires-locked environment montage and then removed during EEG/fMRI acquisition with a subtraction algorithm that takes in account the most significant EPI sequence parameters. The on-line method also allows a further post-processing utilization.     

CORSICA: correction of structured noise in fMRI by automatic identification of ICA components

When applied to functional magnetic resonance imaging (fMRI) data, spatial independent component analysis (sICA), a data-driven technique that addresses the blind source separation problem, seems able to extract components specifically related to physiological noise and brain movements. These components should be removed from the data to achieve structured noise reduction and improve any subsequent detection and analysis of signal fluctuations related to neural activity. We propose a new automatic method called CORSICA (CORrection of Structured noise using spatial Independent Component Analysis) to identify the components related to physiological noise, using prior information on the spatial localization of the main physiological fluctuations in fMRI data. As opposed to existing spectral priors, which may be subject to aliasing effects for long-TR data sets (typically acquired with TR >1 s), such spatial priors can be applied to fMRI data, regardless of the TR of the acquisitions. By comparing the proposed automatic selection to a manual selection performed visually by a human operator, we first show that CORSICA is able to identify the noise-related components for long-TR data with a high sensitivity and a specificity of 1. On short-TR data sets, we validate that the proposed method of noise reduction allows a substantial improvement of the signal-to-noise ratio evaluated at the cardiac and respiratory frequencies, even in the gray matter, while preserving the main fluctuations related to neural activity.     

On the connectivity index for lattices of nonintegral dimensionality

We define the connectivity indexc for an infinite graph by the requirement that to disconnect a subset of at leastV points from the rest of the graph requires the deletion of a minimum ofS(V) bonds whereS(V) ∼V ^(c−1)/c for largeV. For ad-dimensional hypercubical lattice withd integral,c=d. We construct explicit examples of lattices with nonintegral connectivity indexc, 1<c<∞. It is argued that the connectivity index is an important parameter determining the critical behaviour of Hamiltonians on these lattices.     

静息态脑电信号动态功能连接分析

静息态脑功能连接分析是近年来脑研究的一个热点问题, 对于某些脑疾病的诊断及成因理解具有重要意义. 已有的脑功能连接研究基本上都假设功能连接网络在一段时间内是稳定不变的, 但越来越多的证据表明它应该是随时间动态变化的. 对25名被试睁眼和闭眼状态的64电极脑电生理信号, 采用独立成分分析、滑动时间窗、低分辨率脑电断层溯源、图论等方法和技术进行动态功能连接分析, 展现了睁眼和闭眼两种基线状态下视觉网络、默认网络等功能连接网络随时间的动态变化, 并对动态连接矩阵进行主成分分析得到了在整个时间段内具有代表意义的功能连接模式. 该结论支持和补充了传统稳态脑功能连接的研究, 也将为相关实验设计以及脑电信号临床研究提供基线选择依据.     

Effects of spatial smoothing on fMRI group inferences

The analysis of functional magnetic resonance imaging (fMRI) data involves multiple stages of data pre-processing before the activation can be statistically detected. Spatial smoothing is a very common pre-processing step in the analysis of functional brain imaging data. This study presents a broad perspective on the influence of spatial smoothing on fMRI group activation results. The data obtained from 20 volunteers during a visual oddball task were used for this study. Spatial smoothing using an isotropic gaussian filter kernel with full width at half maximum (FWHM) sizes 2 to 30 mm with a step of 2 mm was applied in two levels — smoothing of fMRI data and/or smoothing of single-subject contrast files prior to general linear model random-effects group analysis generating statistical parametric maps. Five regions of interest were defined, and several parameters (coordinates of nearest local maxima, t value, corrected threshold, effect size, residual values, etc.) were evaluated to examine the effects of spatial smoothing. The optimal filter size for group analysis is discussed according to various criteria. For our experiment, the optimal FWHM is about 8 mm. We can conclude that for robust experiments and an adequate number of subjects in the study, the optimal FWHM for single-subject inference is similar to that for group inference (about 8 mm, according to spatial resolution). For less robust experiments and fewer subjects in the study, a higher FWHM would be optimal for group inference than for single-subject inferences.     

Evaluation of preprocessing steps to compensate for magnetic field distortions due to body movements in BOLD fMRI

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is currently the dominant technique for non-invasive investigation of brain functions. One of the challenges with BOLD fMRI, particularly at high fields, is compensation for the effects of spatiotemporally varying magnetic field inhomogeneities (ΔB₀) caused by normal subject respiration and, in some studies, movement of the subject during the scan to perform tasks related to the functional paradigm. The presence of ΔB₀ during data acquisition distorts reconstructed images and introduces extraneous fluctuations in the fMRI time series that decrease the BOLD contrast-to-noise ratio. Optimization of the fMRI data-processing pipeline to compensate for geometric distortions is of paramount importance to ensure high quality of fMRI data. To investigate ΔB₀ caused by subject movement, echo-planar imaging scans were collected with and without concurrent motion of a phantom arm. The phantom arm was constructed and moved by the experimenter to emulate forearm motions while subjects remained still and observed a visual stimulation paradigm. These data were then subjected to eight different combinations of preprocessing steps. The best preprocessing pipeline included navigator correction, a complex phase regressor and spatial smoothing. The synergy between navigator correction and phase regression reduced geometric distortions better than either step in isolation and preconditioned the data to make them more amenable to the benefits of spatial smoothing. The combination of these steps provided a 10% increase in t-statistics compared to only navigator correction and spatial smoothing and reduced the noise and false activations in regions where no legitimate effects would occur.     

fMRI of the auditory system: understanding the neural basis of auditory gestalt

Functional magnetic resonance imaging (fMRI) has rapidly become the most widely used imaging method for studying brain functions in humans. This is a result of its extreme flexibility of use and of the astonishingly detailed spatial and temporal information it provides. Nevertheless, until very recently, the study of the auditory system has progressed at a considerably slower pace compared to other functional systems. Several factors have limited fMRI research in the auditory field, including some intrinsic features of auditory functional anatomy and some peculiar interactions between fMRI technique and audition. A well known difficulty arises from the high intensity acoustic noise produced by gradient switching in echo-planar imaging (EPI), as well as in other fMRI sequences more similar to conventional MR sequences. The acoustic noise interacts in an unpredictable way with the experimental stimuli both from a perceptual point of view and in the evoked hemodynamics. To overcome this problem, different approaches have been proposed recently that generally require careful tailoring of the experimental design and the fMRI methodology to the specific requirements posed by the auditory research. The novel methodological approaches can make the fMRI exploration of auditory processing much easier and more reliable, and thus may permit filling the gap with other fields of neuroscience research. As a result, some fundamental neural underpinnings of audition are being clarified, and the way sound stimuli are integrated in the auditory gestalt are beginning to be understood.