Changes in brain functional network connectivity after stroke

Wei Li, Yapeng Li, Wenzhen Zhu, Xi Chen

1 Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan, Hubei Province, China

2 Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, Hubei Province, China

3 Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei Province, China

Changes in brain functional network connectivity after stroke

Wei Li1,2, Yapeng Li1,2, Wenzhen Zhu3, Xi Chen1,2

1 Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan, Hubei Province, China

2 Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, Hubei Province, China

3 Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei Province, China

Studies have shown that functional network connection models can be used to study brain network changes in patients with schizophrenia. In this study, we inferred that these models could also be used to explore functional network connectivity changes in stroke patients. We used independent component analysis to find the motor areas of stroke patients, which is a novel way to determine these areas. In this study, we collected functional magnetic resonance imaging datasets from healthy controls and right-handed stroke patients following their fi rst ever stroke. Using independent component analysis, six spatially independent components highly correlated to the experimental paradigm were extracted. Then, the functional network connectivity of both patients and controls was established to observe the differences between them. The results showed that there were 11 connections in the model in the stroke patients, while there were only four connections in the healthy controls. Further analysis found that some damaged connections may be compensated for by new indirect connections or circuits produced after stroke. These connections may have a direct correlation with the degree of stroke rehabilitation. Our fi ndings suggest that functional network connectivity in stroke patients is more complex than that in healthy controls, and that there is a compensation loop in the functional network following stroke. This implies that functional network reorganization plays a very important role in the process of rehabilitation after stroke.

nerve regeneration; brain injury; stroke; motor areas; functional magnetic resonance imaging; brain network; independent component analysis; functional network connectivity; neural plasticity; NSFC grant; neural regeneration

Funding: This study was supported by the National Natural Science Foundation of China, No. 60905024.

Li W, Li YP, Zhu WZ, Chen X. Changes in brain functional network connectivity after stroke. Neural Regen Res. 2014;9(1):51-60.

Introduction

Many stroke patients still have long-term functional de fi cits, especially in motor function[1-4]. There are many investigators who have attempted to explore the mechanism of functional recovery after stroke[5-7], yet until now there have been no breakthroughs in this area[8]. In this study, we ask the question: are there changes in the connectivity of motor areas following stroke, and what kinds of changes happen?

In the past two decades, functional magnetic resonance imaging (fMRI), as one of the most recently developed forms of neuroimaging, has been used as a reliable method of clinical diagnosis and assessment for stroke[9-11]. Owing to the relatively low invasiveness and high spatial resolution of fMRI, scientists have started to use this method of neuroimaging to focus on cerebral connectivity[12-14].

A growing number of studies have tried to use connectivity analysis to understand the changes in brain networks between stroke patients and healthy controls as well as the mechanism of neural plasticity after stroke[15-17]. Using dynamic causal models, Mintzopoulos et al.[18]illustrated that stroke patients showed increased coupling from the SMA to M1 and the SMA to the cerebellum. The intrinsic neural coupling between M1 and the cerebellum has also been found to decrease. The results suggest that alterations in brain connectivity between motor areas could compensate for the abnormal function of M1 in chronic stroke patients[18].

With graph-theoretical approaches to the fMRI data of stroke patients, Wang et al.[19]explored functional changes in and the reorganization of the brain network after stroke. The study illustrated that the brain network may evolve to a random mode during the process of rehabilitation after stroke[19]. Grefkes and Gereon[20]pointed out that connectivity analysis could be a useful way to develop new hypothesis-driven treatment methods to promote motor function recovery for stroke patients[20]. Westlake and Nagarajan[21]also found that cerebral network modeling is a powerful tool to explore cerebral reorganization and to predict the prog-nosis of stroke.

Previous studies on brain connectivity for stroke focused on functional connectivity and effective connectivity among motor areas[15,18]. Functional connectivity is defined as the correlation between spatially remote neurophysiological events, and can be used to study the relationship among individual voxels or regions with fMRI activation images[22].

Effective connectivity can be defined as the influence that one neural system exerts over another. The various approaches of effective connectivity, such as dynamic causal modeling, accommodate the nonlinear and dynamic performance of the interactions of neural units[23], and can describe the task-dependent in fl uence that the activation of one area exerts over another to infer the connectivity between selected areas within a brain network[23-24]. However, these methods always require a priori knowledge about neurophysiology and the experimental paradigm, or given assumptions about the speci fi c processing of fMRI data.

Can we use a data-driven approach to study the connectivity between brain areas? Independent component analysis is a useful method of studying functional connectivity with little reliance on previous knowledge[25-27]. It can decompose fMRI data into time courses of spatially independent components to distinguish task-related or non-task-related signal components[28-31]. A signi fi cant advantage of independent component analysis is that it is a model-free and data-driven approach. Thus, independent component analysis is a good selection for the study of brain network connectivity when we lack a priori knowledge to use other techniques.

McKeown et al.[25]showed that independent component analysis is an useful way to decompose fMRI datasets into mutually independent components of the space, from both the experimental fMRI (artificially constructed fMRI datasets) datasets and the real datasets (actually collected fMRI datasets)[25]. Many studies have used independent component analysis to decompose such datasets and produce meaningful results. For example, Demirci et al.[12]used independent component analysis to pretreat the fMRI datasets, after which a Granger causality model was established among the independent component analysis components. Their results suggested that there was a signi fi cant difference in the models of healthy controls and patients with schizophrenia. These results have important theoretical and clinical signi fi cance for the diagnosis and treatment of schizophrenia[12]. Other studies have also used independent component analysis to investigate fMRI datasets, and have achieved good results[31].

Functional network connectivity, into which independent component analysis is integrated, has been regarded as a strong tool to study the functional interactions among selected independent components that represent different activated regions without an a priori model. It can also describe the constrained maximal time-lagged correlation for every pair of component combinations. Therefore, functional network connectivity approaches have been applied in both the study of brain network changes in schizophrenia[26]and the temporal sequence of the activated hemispheric network during semantic processing[32]. Using this method, Jafri et al.[26]reported that there was more functional connectivity and a longer lag time between the independent component analysis components in patients with schizophrenia than in healthy controls.

Assaf et al.[32]used a similar method to study functional network connectivity among brain areas in the left and right hemispheres during the processing of semantic memory. The results indicated an early activation in the right hemisphere that was closely followed by an area of activation in the left hemisphere. However, there is as yet no research that uses functional network connectivity to investigate possible alterations in the brain network after stroke. A major question is whether a new theory of neurological rehabilitation can be discovered by exploring the brains of stroke patients using functional network connectivity.

In this study, acute stroke patients and healthy controls were recruited to take part in blood oxygenation level dependent (BOLD) fMRI experiments while performing alternating unilateral finger-to-thumb opposition movements. Group independent component analysis was used to obtain spatially independent components of brain activation with the time courses of decomposed components[25,27]. Then, to study brain connectivity changes after stroke, functional network connectivity was constructed from the selected components.

Based on the previous study in stroke patients versus healthy controls, it has been hypothesized that the functional network connectivity of stroke patients is more complex than that of healthy controls; moreover, we assumed that functional compensation could be achieved through the other neural circuits in the stroke patients. The purpose of this study was to explore the changes in brain network connectivity after stroke.

Results

Quantitative analysis of participants

The fMRI datasets of eight right-handed patients with acute stroke and eight healthy controls were collected in this study. The independent component analysis method was used to decompose the datasets to independent components after preprocessing each dataset.

In the preprocessing of fMRI datasets, two patients and two healthy controls were excluded for exceeding the head movement parameters. Therefore, in the final analysis, six stroke patients and six healthy controls were used.

Baseline data analysis of subjects

There were six stroke patients, five males and one female, with a mean age of 41.5 (range 17–65) years. Detailed clinical and demographic data for all patients are shown in Table 1. There were six healthy controls, four males and two females, with a mean age of 53 (range 40–58) years.

There were no significant differences in age, gender, and years of education between stroke patients and healthy controls (Table 2).

Selection of independent components

To obtain the connectivity differences between stroke pa-tients and healthy controls, six independent components highly correlated to the experimental paradigm were chosen. The components chosen (A–F) are shown in Figure 1. The Brodmann areas (BA)[33]and MNI coordinates corresponding to the activated regions are also shown in Table 3. The temporal multiple regression correlation between each selected component and the experimental paradigm was computed by temporal sorting analysis as follows. Component A: R = 0.22; component B: R = 0.22; component C: R = 0.215; component D: R = 0.200; component E: R = 0.186; component F: R = 0.185.

Table 1 Clinical and demographic data of stroke patients

Establishment of the functional network connectivity model

Functional network connectivity for healthy controls and stroke patients is shown in Figure 2 and Figure 3. The links indicate that there was a signi fi cant correlation between two different components (P ＜ 0.05). From the figures, there were four links in healthy controls, and 11 links in stroke patients. The arrow direction describes the relationship of the lag between the two linked components. For example, A B means that component A preceded component B according to the group lag average. Through graphical analysis, we foundsigni fi cant correlations for D→ E and C →F for both healthy controls and stroke patients. For healthy controls only, we found connectivity for C →D and E →C, while we found connectivity in stroke patients only for D→ B, B →C, E →A, A →F, F→ D, A→ C, B→ A, F→ B and B →E.

Table 2 Baseline data of stroke patients and healthy controls

From the results, we found compensatory connection loop in the functional network of stroke patients. For example, we found that the correlations C→ D and E →C only existed in healthy controls. There were no direct connections in stroke patients between components C and D, or components E and C.

However, there was an indirect connection between components C and D in stroke patients through component F with the path: C →F →D. There was also an indirect connection in stroke patients between components E and C through component A with the path: E →A→ C.

Discussion

In this paper, functional network connectivity was used to study functional connectivity during rehabilitation following stroke and was used to explore the differences between stroke patients and healthy controls.

Previous studies have shown that independent component analysis is one of the most effective methods for fMRI dataset analysis[29,34-35]. Independent component analysis and functional network connectivity can be used to distinguish the difference between patients with schizophrenia and controls[26], and can also be used to study the temporal sequence of hemispheric network activation during semantic processing[32]. We studied changes in the functional network connectivity of stroke patients during movements using independent component analysis and the functional network connectivity method. Our results illustrated that the functional network connectivity of stroke patients was more complex than that of healthy controls.

In addition, a compensatory connection loop existed in the functional network of stroke patients. These results illustrate that the changes in functional network connectivity can be explored by the dual methods of independent component analysis and functional network connectivity.

The independent component analysis method is a data-driven method for preprocessing fMRI datasets. Many fMRI datasets have been decomposed in this way in a large number of recent studies. Independent component analysis can decompose the fMRI datasets into mutually independent components in the spatial domain, and a single component of time synchronization. This method has been used in different studies and has obtained good results[30-31]. In this study, the fMRI datasets of stroke patients and healthy controls were fi rst preprocessed using the independent componentanalysis method.

Table 3 The details of activated regions for each component

Figure 2 Functional network connectivity model for the healthy controls.

In this study, the functional regions involved in component A include the supplementary motor cortex and premotor cortex. These regions have been shown to play a crucial role in motor control[22,36-37]. Component B mainly represents primary somatosensory cortex[38]. These regions of component C are generally associated with the recognition of human beings[39]. They are also considered to have a relationship with the motor cortex. Component D may be related to sensory systems[40].

Component E has been shown to be involved in directing attention in space both when an individual makes movements and while imaging or preparing the movements[41]. From the selected components, we can fi nd that many of the areas of movement were found by the independent component analysis method, although we did not obey the previous seed-voxel method. The results show the feasibility of selecting components with this method, and provide a theoretical basis for future research.

Previous studies on stroke have demonstrated that cortical activation is increased in stroke patients in comparison with healthy controls in the contralesional primary sensorimotor cortex during the performance of movement[42]. It also has been suggested that the disorganization of the motor network after stroke is increased, notably including more bilateral involvement of BA4p[43]. This result indicates that the neurovascular coupling may be affected by infarct lesions in patients[44-46].

Figure 3 Functional network connectivity model for the stroke patients.

Wang et al.[19]found from analysis of the whole brain network that the connectivity of patients was much more complex than that of healthy controls. The result was consistent with our hypothesis. The same results were also found in a study by Jafri et al.[26]. Thus, the optimal structure of the brain may change during disease, which is also consistent with a previous study of brain injury[47].

The network connections in the healthy controls were simple and of good group consistency, which may imply a simple but efficient network. For stroke patients, we could interpret their complex connectivity as evidence that plenty of new connections were produced to compensate for the damaged or injured connections and nerves[48-49]. Although the network connections in stroke patients may be less stable and ef fi cient than those in controls, this solution evolved after the stroke and does work[50-51].

Previous studies have suggested that the functional networks of patients with brain disease were more complex than those of healthy controls[52]and our study also demonstrated this fact, which could help us to comprehend the mechanism of neurological rehabilitation after stroke. These results have very important theoretical and clinical significance for stroke.

The analysis of functional network connectivity could also be used to obtain the functional connectivity of various other brain regions and the temporal relationships among the identi fi ed components[26]. Theoretically, the time lag between the two components could be due to temporal differencesin blood supply to different brain regions[32,53-54]. Results from this study showed that the signi fi cant correlations D E and C F existed in both healthy controls and stroke patients, which suggests that these connections, including correlations and temporal relationships, existed after stroke in patients, and still worked during performance of movements.

From our results, we also found some connections which were found in healthy controls disappeared in the stroke patients, while these components can be connected through other components. This shows that some cerebral regions play an important role as a bridge to connect damaged connections, such as component F or component D. Thus, it suggests that although some connections were damaged after stroke, indirect connections were produced and functional compensation may be achieved through such new indirect connections or neural circuits[55-56].

Generally, the studies illustrated that the human brain has a strong robustness[57], and that the function of the injured brain regions could be compensated for by their peripheral regions or nerves that have similar functions. This may be due to the long-term evolution of the brain. In this study, a loop existed in the functional network connectivity of stroke patients although small components were selected. This conclusion further demonstrates the robustness of the brain, and confirms the reliability of the previous studies. At the same time, it illustrates that functional compensation plays a very important role in functional rehabilitation for stroke patients.

It should be mentioned here that there are several limitations to this study. The sample set is not large enough, and the behavioral data have not been included in the network connectivity analysis. In future research, it may be meaningful to combine the behavioral data and functional network connectivity analysis to explore the possible correlations between them after stroke.

In conclusion, differences in functional network connectivity between stroke patients and healthy controls were explored in this paper. Six spatially independent components highly correlated to the experimental paradigm were extracted using independent component analysis to fi nd the nodes of the network. The results of network modeling illustrated that functional network connectivity in stroke patients was much more complex than that in healthy controls. This result suggests that many new connections are produced to compensate for the damaged connections and nerves after stroke.

Furthermore, it also has been found that some damaged connections may be compensated for by new indirect connections or circuits produced after stroke, which implies that functional network reorganization plays a very important role in the process of rehabilitation after stroke.

Subjects and Methods

Design

A case-control study.

Time and setting

This study was performed at the Laboratory of Nuclear Magnetic Resonance, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, China between April 2010 and November 2010.

Subjects

Each subject was recognized as having had a stroke using a T1 image. The patients were recruited by the Department of Rehabilitation, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, China. All participants were stroke patients who received treatment in Tongji Hospital and came from Hubei Province, China.

All stroke patients fulfilled the following inclusion criteria: (1) first-ever ischemic stroke; (2) within 2 weeks after stroke; (3) showing motor de fi cits with acute unilateral loss of hand strength (grade 4 on the Medical Research Council (MRC) scale (0–5, the higher the score, the greater the muscle strength, a score of 5 = normal))[58]; (4) good hand motor function.

Exclusion criteria for stroke patients were as follows: (1) language or cognitive de fi cits that would impair cooperation in fMRI examination; (2) significant somatosensory (light touch or proprioception) deficits in the stroke-affected hand; (3) mirror movements; and (4) contraindication to magnetic resonance imaging.

In this study, healthy controls (4 male; age range: 40–58 years; mean: 53 years) with a similar age and education to the patients with stroke (5 male; age range: 17–65 years; mean: 41.5 years) were selected. Healthy controls were recruited by the Department of Rehabilitation, Tongji Hospital, and all of them came from Hubei Province. All healthy controls fulfilled the following inclusion criteria: (1) No previous history of brain injury; (2) no history of mental illness; (3) normal motor function. Exclusion criteria for the healthy controls were as follows: (1) contraindication to MRI; (2) the subjects could not complete the experiment.

This study was approve by the Ethics Committee of Tongji Hospital, Tongji Medical college, Huazhong Unversity of Science and Technology, China, all subjects gave their informed consent before being submitted to fMRI examination.

The clinical characteristics of the stroke patients are summarized in Table 2. All six patients had a subcortical stroke encompassing the left basal ganglia or the left frontal cortex (Table 1), and the illustration of the lesion location is shown in red in Figure 4. For all stroke patients, the localization of infarct was almost the same.

Methods

Experimental paradigm

All the subjects were instructed to execute alternating unilateral fi nger-to-thumb opposition movements at a frequency of 1 Hz in a block-design fMRI paradigm.

The task occurred in 20-second blocks of movements alternated with 20-second intervals as rest periods: Rest –Movements (Left) – Rest – Movements (Right) – Rest –Movements (Left) – Rest – Movements (Right) – Rest –Movements (Left) – Rest – Movements (Right) – Rest. It lasted for a total of 260 seconds as shown in Figure 5.

Figure 4 Illustration of lesion location in red for each patient.

During the fMRI procedure, the subjects kept their eyes open, and their head was kept motionless. Prior to scanning, the subjects were instructed and trained until they fully understood the task and were able to adequately follow the cues and instructions.

Parameters of fMRI and data acquisition

MRI scans were acquired on a 3.0T GE Signaxs scanner (Signa HDxt, GE Healthcare, Fair fi eld, CT, USA) with a custom-built head coil. A high resolution T1-weighted spoiled grass gradient recalled (SPGR) inversion recovery three-dimensional MRI sequence was performed for each subject with the following parameters: inversion time = 400 ms; repetition time = 6.5 ms; echo time = 2.1 ms; fl ip angle = 15 degrees; fi eld of view = 25.6 cm; 132 slices in coronal plane; matrix = 256 × 256; number of excitations = 1; acquired resolution = 1.5 mm × 0.9 mm × 1.1 mm.

BOLD signal was collected with a T2-weighted gradient echo spiral in-out pulse sequence[59]with the following parameters: repetition time = 2,000 ms; echo time = 30 ms; fl ip angle = 90 degrees; 1 interleave; fi eld of view = 24 cm; 64 × 64 matrix. A total of 32 axial slices (5.0 mm thickness, 0 mm skip), which were parallel to the AC-PC line and covered the whole brain, were obtained with a temporal resolution of 2 seconds. For the whole task, 120 images were obtained for a total duration of 4 minutes. Structural and functional scans were acquired in the same scan session.

Preprocessing of fMRI data

The fMRI datasets were preprocessed using SPM8 software (www.fil.ion.ucl.ac.uk/spm). First, the dicom dataset was converted into ⋆.img/⋆.hdr document. Then, all image volumes were realigned to the mean volume. The subjects whose head displacements were more than 2 mm in the X, Y, Z direction or whose head rotation exceeded 1° were excluded (two patients were excluded for this reason and are not shown in Table 1).

Using the unified segmentation approach[60], functional images were normalized to the MNI template (voxel 3 mm × 3 mm × 3 mm). In order to decrease spatial noise, volumes were smoothed by a 4-mm full-width half maximum Gaussian kernel.

Independent component selection

The procedure of component identification and selection described by Jafri et al.[26]was followed in this study. First, the preprocessed data of all twelve participants (six stroke patients and six healthy controls) were decomposed using the group spatial independent component analysis. Then, 36 independent components were obtained using the GIFT toolbox (http://icatb.sourceforge. net/). The rules of component selection should include the principles described as follows.

Figure 5 The movement paradigm of the subjects in the magnetic fi eld.

First, only those components highly correlated to gray matter would be chosen, as the activation of those components always represent signi fi cant hemodynamic change[27]. Second, different from prior studies in which M1, S1 and PMd regions of interest were extracted from the template[2], only those components highly correlated to our experimental paradigm would be chosen. In standard general linear model analysis, the canonical hemodynamic response function obtained from SPM8 software was convolved with the experimental paradigm time course to get the BOLD time courses. Then, GIFT was used to obtain the association coef fi cients between the time courses of components and our experimental paradigm. GIFT is a software package that includes different independent component analysis methods.

In this study, the information algorithm was used in the preprocessing of the fMRI datasets. The independent components whose correlation coef fi cients were greater than 0.15 between the activation of independent component and the task were chosen.

Establishing the functional network connectivity model

In this study, six spatially independent component analysis components were selected. However, according to their temporal properties, there were correlations between them. A constrained maximal lagged correlation between every pairwise combination of time courses of the independent component analysis components selected was computed by a functional network connectivity toolbox (http://mialab. mrn.org/software).

Theoretically, there should be a maximum of 15 connections among the six components. The statistical analysis was performed using SPSS 13.0 software (SPSS, Chicago, IL, USA). The connections between any two components were fi rst calculated, and then an independent sample t-test was used to test the signi fi cance of any connection. If the P value was ＜ 0.05, the connection was established. The group average of the lag was evaluated for two groups (stroke patients and healthy controls). For more detail on this method, refer to Jafri et al.[26].

Acknowledgments:We would like to thank all the patients and volunteers for their participation in this study. We also express our thanks to Qin YY for data acquisition and Xu J and Luo X for their insightful comments.

Author contributions:Li W conceived and designed this study, conducted data analysis, revised the manuscript and guided the study. Li YP conducted data analysis and wrote the manuscript. Zhu WZ collected experimental data. Chen X participated in statistical analysis and manuscript revision. All authors approved the final version of the paper.

Con fl icts of interest:None declared.

Peer review:The present findings confirm that the functional network connectivity of motor areas in stroke patients is more complex than that in healthy controls, and that some damaged connections can be compensated by new indirect connections or circuits produced after stroke. These findings provide theoretical and clinical significance for brain rehabilitation and nerve regeneration after stroke.

[1] Calautti C, Naccarato M, Jones PS, et al. The relationship between motor de fi cit and hemisphere activation balance after stroke: A 3T fMRI study. NeuroImage. 2007;34(1):322-331.

[2] Richards LG, Stewart KC, Woodbury ML, et al. Movement-dependent stroke recovery: a systematic review and meta-analysis of TMS and fMRI evidence. Neuropsychologia. 2008;46(1):3-11.

[3] Kokotilo KJ, Eng JJ, McKeown MJ, et al. Greater activation of secondary motor areas is related to less arm use after stroke. Neurorehabil Neural Repair. 2010;24(1):78-80.

[4] Calautti C, Jones PS, Naccarato M, et al. The relationship between motor deficit and primary motor cortex hemispheric activation balance after stroke: longitudinal fMRI study. J Neurol Neurosurg Psychiatry. 2007;81(7):788-792.

[5] Ward NS, Brown MM, Thompson AJ, et al. Neural correlates of motor recovery after stroke: a longitudinal fMRI study. Brain. 2003;126(11):2476-2496.

[6] Askim T, Indredavik B, Håberg A. Internally and externally paced finger movements differ in reorganization after acute ischemic stroke. Arch Phys Med Rehabil. 2010; 91(10):1529-1536.

[7] Rehme AK, Fink GR, Cramon DY, et al. The role of the contralesional motor cortex for motor recovery in the early days after stroke assessed with longitudinal fMRI. Cerebral Cortex. 2010;21(4):2-13.

[8] Enzinger C, Dawes H, Berg HJ, et al. Brain activity changes associated with treadmill training after stroke. Stroke. 2009;40(7):2460-2467.

[9] Higashida RT, Furlan AJ, Roberts H et al. Trial design and reporting standards for intra-arterial cerebral thrombolysis for acute ischemic stroke. Stroke. 2003;34(8):e109-137.

[10] Förster A, Griebe M, Ottomeyer C, et al. Cerebral network disruption as a possible mechanism for impaired recovery after acute pontine stroke. Cerebrovasc Dis. 2011;31(5):499-505.

[11] Arsava EM. The role of MRI as a prognostic tool in ischemic stroke. J Neurochem. 2012;123(2):22-28.

[12] Demirci O, Stevens MC, Andreasen NC, et al. Investigation of relationships between fMRI brain networks in the spectral domain using ICA and Granger causality reveals distinct differences between schizophrenia patients and healthy controls. Neuroimage. 2009;46(2):419-431.

[13] Eickhoff SB, Dafotakis M, Grefkes C, et al. Central adaptation following heterotopic hand replantation probed by fMRI and efective connectivity analysis. Exp Neurol. 2008;212(1):132-144.

[14] Grefkes C, Nowak DA, Wang LE, et al. Modulating cortical connectivity in stroke patients by rTMS assessed with fMRI and dynamic causal modeling. NeuroImage. 2010;50(1):233-242.

[15] Murase N, Duque J, Mazzocchio R, et al. In fl uence of interhemispheric interactions on motor function in chronic stroke. Ann Neurol. 2004;55(3):400-409.

[16] Hummel F, Celnik P, Giraux P, et al. Effects of non-invasive cortical stimulation on skilled motor function in chronic stroke. Brain. 2005;128(3):490-409.

[17] He BJ, Snyder AZ, Vincent JL, et al. Breakdown of functional connectivity in frontoparietal networks underlies behavioral de fi cits in spatial neglect. Neuron. 2007;53(6):905-118.

[18] Mintzopoulos D, Astrakas LG, Khanicheh A, et al. Connectivity alterations assessed by combining fMRI and MR-compatible hand robots in chronic stroke. Neuroimage. 2009;47 Suppl 2:T90-T97.

[19] Wang L, Yu CS, Chen H, et al. Dynamic functional reorganization of the motor execution network after stroke. Brain. 2010;133(4):1224-1238.

[20] Grefkes C and Gereon RF. Reorganization of cerebral networks after stroke: new insights from neuroimaging with connectivity approaches. Brain. 2011;134(5):1264-1276.

[21] Westlake KP, Nagarajan SS. Functional connectivity in relation to motor performance and recovery after stroke. Front Syst Neurosci. 2011;5:8.

[22] Biswal B, Yetkin FZ, Haughton VM, et al. Functional connectivity in the motor cortex of resting human brain using echoplanar MRI. J Magn Reson Med. 1995;34(4):537- 541.

[23] Friston KJ, Harrison L, Penny W. Dynamic causal modelling. Neuroimage. 2003;19(4):1273-1302.

[24] Friston KJ, Frith CD, Liddle PF, et al. Functional connectivity: the principal component analysis of large (PET) data sets. J Cereb Blood Flow Metab. 1993;13(1):5-14.

[25] McKeown MJ, Makeig S, Brown GG, et al. Analysis of fmri data by blind separation into independent spatial components. Hum Brain Mapp. 1998;6(3):160-188.

[26] Jafri MJ, Pearlson GD, Stevens M, et al. A method for functional network connectivity among spatially independent resting-state components in schizophrenia. NeuroImage. 2008;39:1666-1681.

[27] Calhoun VD, Adali T, Pearlson GD, et al. Spatial and temporal independent component analysis of functional MRI data containing a pair of task-related waveforms. Hum Brain Mapp. 2001;13:43-53.

[28] Turner GH, Twieg DB. Study of temporal stationarity and spatial consistency of fMRI noise using independent component analysis. IEEE Trans Med Imag. 2005;24(6):712-718.

[29] Calhoun VD, Adali T, Stevens MC, et al. Semi-blind ICA of fMRI: a method for utilizing hypothesis-derived time courses in a spatial ICA analysis. Neuroimage. 2005;25(2):527-538.

[30] Correaa N, AdalVa T, Calhoun VD. Performance of blind source separation algorithms for fMRI analysis using a group ICA method. Magn Reson Imaging. 2007;25(5):684-694.

[31] Chen H, Yao D. Discussion on the choice of separated components in fMRI data analysis by spatial independent component analysis. Magn Reson Imaging. 2004;22(6): 827-833.

[32] Assaf M, Jagannathan K, Calhoun VD, et al. Temporal sequence of hemispheric network activation during semantic processing: a functional network connectivity analysis. Brain Cogn. 2009;70(2): 238-246.

[33] Brodmann K. Vergleichende Lokalisationslehre der Grosshirnrinde. Leipzig: Johann Ambrosius Bart, Germany. 1909.

[34] Sui J, Adali T, Pearlson GD, et al. An ICA-based method for the identification of optimal FMRI features and components using combined group-discriminative techniques. Neuroimage. 2009;46(1):73-86.

[35] Calhoun VD, Liu J, Adali T. A review of group ICA for fMRI data and ICA for joint inference of imaging, genetic, and ERP data. Neuroimage. 2009;45(S1):63-72.

[36] Ruby P, Decety J. Effect of subjective perspective taking during simulation of action: a PET investigation of agency. Nat Neurosci. 2001;4(5):546-550.

[37] Corbetta D, Sirtori V, Moja L, et al. Constraint-induced movement therapy in stroke patients: systematic review and meta-analysis. Eur J Phys Rehabil Med. 2010;46(4):537-544.

[38] Wasaka T, Kakigi R. The effect of unpredicted visual feedback on activation in the secondary somatosensory cortex during movement execution. BMC Neurosci. 2012;5(13):138.

[39] Tomasi D, Volkow ND. Ultrafast method for mapping local functional connectivity hubs in the human brain. Conf Proc IEEE Eng Med Biol Soc. 2010;2010:4274-4277.

[40] Simpson GV, Weber DL, Dale CL, et al. Dynamic activation of frontal, parietal, and sensory regions underlying anticipatory visual spatial attention. J Neurosci. 2011;31(39):13880-13889.

[41] Cavanna AE, Trimble MR. The precuneus: a review of its functional anatomy and behavioural correlates. Brain. 2006;129(3):564-583.

[42] Schaechter JD, Perdue KL. Enhanced cortical activation in the contralesional hemisphere of chronic stroke patients in response to motor skill challenge. Cerebral Cortex. 2008;18(3):638-647.

[43] Sharma N, Simmons LH, Jones PS, et al. Motor imagery after subcortical stroke a functional magn reson imaging study. Stroke. 2009;40(4):1315-1324.

[44] Gerloff C, Bushara K, Sailer A, et al. Multimodal imaging of brain reorganization in motor areas of thecontralesional hemisphere of well recovered patients after capsularstroke. Brain. 2006; 129(3):791-808.

[45] Nomura EM, Gratton C, Visser RM, et al. Double dissociation of two cognitive control networks in patientswith focal brain lesions. Proc Natl Acad Sci USA. 2010;107(26):12017-12022.

[46] Murata Y, Sakatani K, Hoshino T, et al. Effects of cerebral ischemia on evoked cerebral blood oxygenation responses and BOLD contrast functional MRI in stroke patients. Stroke. 2006;37(10):2514-2520.

[47] Stam C, Jones B, Nolte G, et al. Small-World networks and functional connectivity in Alzheimer’s disease. Cerebral Cortex. 2007;17(1):92-99.

[48] Warren JE, Crinion JT, Lambon Ralph MA, et al. Anterior temporal lobe connectivity correlates with functional outcome after aphasic stroke. Brain. 2009;132(12):3428-3442.

[49] van Meer MP, van der MK, Wang K, et al. Recovery of sensorimotor function after experimental stroke correlates with restoration of resting-state interhemispheric functional connectivity. J Neurosci. 2010;30(11):3964-3972.

[50] Kollen B, Kwakkel G, Lindeman E. Functional recovery after stroke: a review of current developments in stroke rehabilitation research. Rev Recent Clin Trials. 2006;1(1):75-80.

[51] Voytek B, Davis M, Yago E, et al. Dynamic neuroplasticity after human prefrontal cortex damage. Neuron. 2010;68(3):401-408

[52] Rowe JB. Connectivity analysis is essential to understand neurological disorders. Front Syst Neurosci. 2010;4:144.

[53] Calvisi ML, Szeri AJ, Liley DT, et al. Theoretical study of BOLD response to sinusoidal input. Conf Proc IEEE Eng Med Biol Soc. 2004;1:659-662.

[54] Crijns SP, Raaymakers BW, Lagendijk JJ. Real-time correction of magnetic field inhomogeneity-induced image distortions for MRI-guided conventional and proton radiotherapy. Phys Med Biol. 2011;56(1):289-297.

[55] de Vico FF, Astol fi L, Cincotti F, et al. Evaluation of the brain network organization from EEG signals: a preliminary evidence in stroke patient. Anat Rec (Hoboken). 2009;292(12):2023-2031.

[56] Johansen-Berg H, Scholz J, Stagg CJ. Relevance of structural brain connectivity to learning and recovery from stroke. Front Syst Neurosci. 2010;2(4):146.

[57] Bullmore E, Sporns O. The economy of brain network organization. Nat Rev Neurosci. 2012;13(5):336-349.

[58] Medical Research Council. Aids to the examination of the peripheral nervous system, Memorandum no. 45, Her Majesty’s Stationery Of fi ce, London. 1981.

[59] Glover GH, Law CS. Spiral-in/out BOLD fMRI for increased SNR and reduced susceptibility artifacts. Magn Reson Med. 2001; 46(3):515-522.

[60] Ashburner J, Friston KJ. Unified segmentation. Neuroimage. 2005;26(3):839-851.

Copyedited by Jackson C, Stow A, Yin JZ, Jiang GH, Wang J, Li CH, Song LP, Liu WJ, Zhao M

10.4103/1673-5374.125330

Wei Li, Associate professor, Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan 430074, Hubei Province, China, Liwei 0828@mail.hust.edu.cn.

http://www.nrronline.org/

Accepted: 2013-07-16