Brain Functional Connectivity Plasticity Within and Beyond the Sensorimotor Network in Lower-Limb Amputees
Jingna Zhang
Ye Zhang1 
Li Wang
Xuntao Yin- 1Department of Medical Imaging, College of Biomedical Engineering, Army Medical University, Chongqing, China
- 2Department of Rehabilitation, Southwest Hospital, Army Medical University, Chongqing, China
- 3Department of Radiology, Southwest Hospital, Army Medical University, Chongqing, China
Cerebral neuroplasticity after amputation has been elucidated by functional neuroimaging. However, little is known concerning how brain network-level functional reorganization of the sensorimotor system evolves following lower-limb amputation. We studied 32 unilateral lower-limb amputees (LLAs) and 32 matched healthy controls (HCs) using resting-state functional magnetic resonance imaging (rs-fMRI). A regions of interest (ROI)-wise connectivity analysis was performed with ROIs in eight brain regions in the sensorimotor network to investigate intra-network changes, and seed-based whole-brain functional connectivity (FC) with a seed in the contralateral primary sensorimotor cortex (S1M1) was used to study the FC reorganization between the sensorimotor region (S1M1) and other parts of the brain in the LLAs. The ROI-wise connectivity analysis showed that the LLAs had decreased FC, mainly between the subcortical nuclei and the contralateral S1M1 (p < 0.05, FDR corrected). Seed-based whole-brain FC analysis revealed that brain regions with decreased FC with the contralateral S1M1 extended beyond the sensorimotor network to the prefrontal and visual cortices (p < 0.05, FDR corrected). Moreover, correlation analysis showed that decreased FC between the subcortical and the cortical regions in the sensorimotor network progressively increased in relation to the time since amputation. These findings indicated a cascade of cortical reorganization at a more extensive network level following lower-limb amputation, and also showed promise for the development of a possible neurobiological marker of changes in FC related to motor function recovery in LLAs.
Introduction
With increases in traffic accidents, malignant diseases and natural disasters, the number of lower-limb amputees (LLAs) is increasing (Barbosa et al., 2016). LLAs have difficulties in motor control and coordination, especially in the early stages following amputation because of the major structural asymmetry arising from the amputation (Molina Rueda et al., 2013). With a period of rehabilitation and the application of a prosthesis, the functional status of amputees is gradually restored (Soares et al., 2009; Rusaw and Ramstrand, 2011), which could induce functional reorganization in the sensory and motor areas of the brain.
One of the most prominent consequences of limb amputation concerns the functional reorganization in the primary somatosensory and motor cortices. Several studies in non-human primates have demonstrated that the loss of an upper limb results in functional reorganization of the primary sensorimotor cortex (S1M1) (Kaas et al., 1983; Merzenich et al., 1983; Pons et al., 1991; Jain et al., 2008). In humans, such reorganization after limb amputation has also been observed, and this has been interpreted as a form of maladaptive plasticity in these cortices that is triggered by the loss of sensory inputs (MacIver et al., 2008) or voluntary control (Raffin et al., 2016). However, most of the studies dealing with cortical reorganization included only patients that have undergone upper-limb amputation. The patterns of brain reorganization might vary significantly in LLAs due to differentiation in the functions and representations between the upper and lower limbs. A remapping of the cortical topography was observed in LLAs, with an expansion of activation maps of the stump in the S1M1 of the deafferent hemisphere, spreading to neighboring regions that represent the trunk and upper limbs (Simões et al., 2012), which indicates that lower-limb amputation also induces neuroplastic changes in the S1M1.
Although emphasis has been given to the sensorimotor cortex in the functional reorganization following limb amputation, the contribution of subcortical structures cannot be overlooked (Curt et al., 2011). The cerebellum and the basal ganglia are groups of subcortical nuclei with long-established roles in motor control (Hoshi et al., 2005; Bostan and Strick, 2010; Bostan et al., 2013). Both of these structures are important subcortical structures in the motor circuit, and perform distinct functional operations such as the initiation and execution of voluntary movement (Groenewegen, 2003). However, only a few studies have examined functional and structural alterations of these subcortical nuclei following lower-limb amputation. A fMRI study found increased activation in the contralateral basal ganglia during motor imagery of the amputated toes in LLAs (Romero-Romo et al., 2010). A voxel-based morphometric (VBM) study found decreased gray matter volume in the bilateral cerebellum in LLAs (Di Vita et al., 2018). Structural and functional abnormalities of the thalamus have also been found following peripheral deafferentation. Garraghty and Kaas (1991) found functional reorganization in the thalamus of adult monkeys after peripheral nerve injury. Draganski et al., (2006) showed a significant decrease in the gray matter volume of the thalamus following limb amputation. Results from a study by Jones and Pons (1998) demonstrated that long-standing limb amputation can cause structural reorganization of the thalamus. These findings implied that these movement-related subcortical structures are possibly involved in functional reorganization in LLAs. However, how these subcortical structures work during the process of functional reorganization following lower-limb amputation has not yet been elucidated.
In recent years, functional connectivity (FC) analyses have provided invaluable approaches for studying the human brain in healthy (Smith et al., 2009, 2012; Yeo et al., 2011) and diseased groups (Napadow et al., 2010; Collignon et al., 2013; Rocca et al., 2014) on the brain-network level. Among the varied functional neuroimaging techniques, resting-state fMRI (rs-fMRI) is a promising tool for mapping FC. This method does not require the participants to accomplish complex sensorimotor task during functional neuroimaging acquisition. Resting-state FC analysis has been used in studying the network-level reorganization of FC following arm amputation and has found reduced FC between the cortex associated with the missing hand and the sensorimotor network in amputees (Makin et al., 2015). Regarding lower-limb amputation, one diffusion tensor imaging (DTI) study observed that the mean fractional anisotropy (FA) value of the fibers in region II of the corpus callosum, which connects the premotor area and supplementary motor area (SMA), was significantly reduced, demonstrating that structural connectivity between the bilateral sensorimotor cortex was reorganized in LLAs (Li et al., 2017). To our knowledge, there have been no prior studies of changes in FC in LLAs. Thus, resting-state FC analysis may offer a new way to explore functional network reorganization in the brain after lower-limb amputation.
Besides the regions within the sensorimotor network, abnormal brain reorganization has been found in other regions (Preißler et al., 2013; Jiang et al., 2015; Makin et al., 2015). Specifically, Preißler et al. (2013) found that upper-limb amputees were associated with anatomical alterations in parts of the brain region that belong to dorsal visual stream. Such plasticity was hypothesized to reflect a brain adaptation process to new movement and coordination patterns in operating hand prosthesis. Meanwhile, Jiang et al. (2015) found a significantly lower thickness in the motor-related visual cortex after lower-limb amputation which was presumably related to the degeneration of biological motion perception or tactile motion processing. Additionally, Makin et al. (2015) found that in upper-limb amputees, the missing hand cortex gradually became functionally coupling with the default mode network (DMN) and decoupling with the sensorimotor network over the time since amputation. Such network-level cortical reorganization was supposed to be related to complex perceptual experiences of phantom sensations. Thus we conclude that limb amputation may cause more extensive brain reorganization beyond the changes that occur in the sensorimotor network.
According to the research mentioned above, amputation may result in extensive structural and functional reorganization in the brain. However, very little is known about large-scale network-level reorganization of FC following lower-limb amputation, which requires further investigation. Here, we take these ideas further by characterizing intra-network (within the sensorimotor network) and inter-network (between S1M1 and other parts of the brain) changes following lower-limb amputation using resting-state FC analysis. We hypothesized that the LLAs would show abnormal FC within and beyond the sensorimotor network compared with the healthy controls (HCs); meanwhile, when considering the long duration of recovery for some amputees, we also hypothesized that FC in the LLAs may present progressive changes in relation to the time since amputation. To test these hypotheses, we compared the regions of interest (ROI)-wise connectivity within the sensorimotor network and the seed-based whole-brain FC of the sensorimotor network between the LLAs and HCs. Further, we investigated the relationship between FC and time since amputation in LLAs.
Materials and Methods
Participants
Thirty-two individuals (23 male and 9 female) with acquired unilateral lower-limb amputation (mean age ± SD: 46.03 ± 12.04, 13 patients with amputations on the right side, see Table 1) were recruited through the Prosthetic and Orthotic Clinics at the Department of Rehabilitation, Southwestern Hospital (Chongqing, China) between March 2015 and April 2016. Nineteen amputations occurred at the transfemoral level and 13 at the transtibial level. All the patients had been fitted with prostheses. Twenty-five patients underwent amputation following a traumatic injury, and the other amputations were due to tumors (three melanoma, two osteosarcoma, and two osteomyelitis). Phantom limb pain (PLP) was assessed using a 10 cm Visual Analog Scale (VAS) (0: “no pain at all” represented by the left end point and 10: “the strongest pain I can imagine” represented by the right end point of the scale (Scott and Huskisson, 1976). Exclusion criteria were the following: (1) amputation at another part of the body, (2) history of brain injury due to trauma, (3) presence of neurological or psychiatric disorders, (4) presence of stump pain, and (5) duration between amputation and MRI scanning of less than 3 months. Thirty-two age- and sex-matched HCs without neurological or psychiatric diseases and with normal brain MRI (no brain atrophy, tumor, ischemia, hemorrhage, or congenital abnormalities) were recruited from the local community. All subjects were right-handed, as assessed using the Edinburgh Handedness Inventory (Oldfield, 1971). Each participant was informed of the purpose and methods of the study and provided a written consent form prior to MRI. All experimental procedures were approved by the Ethics Committee of the Army Medical University.
TABLE 1. Demographic and clinical characteristics of the participants.
Image Acquisition
A Siemens 3.0-T Trio MRI scanner (Magnetom Trio, Siemens, Erlangen, Germany) was used to acquire the functional and structural data. Rs-fMRI data were collected via a multi-slice, gradient-echo, echo-planar sequence for a total of 240 volumes (480 s) using the following parameters: repetition time (TR) = 2000 ms; echo time (TE) = 30 ms; flip angle = 90°; matrix = 64 × 64; field of view = 200 mm; slice thickness = 3 mm, 1-mm gap, 36 slices with a voxel size = 3 mm × 3 mm × 3 mm. During the resting-state scan, subjects were asked to relax with their eyes closed and not to think of anything in particular. High-resolution three-dimensional T1-weighted structural images were also collected in the sagittal orientation via a magnetization-prepared rapid gradient-echo sequence (MP-RAGE) using the following parameters: TR = 1,900 ms, TE = 2.52 ms, flip angle = 9°, matrix = 256 × 256, thickness = 1.0 mm, 176 slices with a voxel size = 1 mm × 1 mm × 1 mm.
Data Preprocessing
Before data preprocessing, in order to obtain a ‘homogeneous’ sample of patients with a left side amputation, we flipped the imaging data from right to left along the midsagittal line for the 13 patients with an amputation on the right side. For all patients, the left side corresponded to the cerebral hemisphere ipsilateral to the side of the amputation and the right side corresponded to the cerebral hemisphere contralateral to the side of the amputation.
Preprocessing was performed using the Data Processing Assistant for Resting-State fMRI (DPARSF) toolbox1. The first 10 images of each subject were discarded to allow the signal to reach equilibrium and the participants to adapt to the scanning noise. Then, the remaining 230 images were corrected for temporal delay in slice acquisition and co-registered to the first image for correction of rigid- body head motion. For all subjects, the translation or rotation parameters did not exceed ± 2.5 mm or ± 2.5°, respectively. Since rs-fMRI data have been demonstrated to be sensitive to head motion (Yan et al., 2013a), the head motion effects were regressed using the Friston 24-parameter model (Satterthwaite et al., 2013; Yan et al., 2013a). To further reduce the effects of confounding factors, including the signals from the white matter and cerebrospinal fluid, the mean time series of all voxels across the whole brain were removed from the data via linear regression (Yan et al., 2013b). The resulting maps were spatially normalized into a standard stereotaxic space at a resolution 3 mm × 3 mm × 3 mm using an EPI template. After normalization, the images were smoothed using a full-width at half-maximum Gaussian kernel of 6 mm to increase the signal-to-noise ratio. Finally, a temporal filter (0.01–0.08 Hz) was applied to reduce low-frequency drift and high-frequency physiological noise.
Seed Selection and FC Analysis
For ROI-wise connectivity, we selected eight subcortical and cortical sensory- and motor-related brain regions based on a freely available atlas of regions defined by correlated activation patterns2 (Shirer et al., 2012), including the bilateral S1M1 with peak Montreal Neurological Institute (MNI) coordinates of (-54, -12, 34) and (58, -10, 34), a bilateral region of the SMA with a peak MNI coordinates of (6, -10, 52), a bilateral region of the cerebellum with a peak MNI coordinates of (-24, -40, -36), the bilateral basal ganglia with peak MNI coordinates of (-22, 2, -6) and (24, 4, -4), and the bilateral thalamus with peak MNI coordinates of (-8, -16, -2) and (10, -16, -2) as the ROIs. In this atlas, the basal ganglia and thalamus on each side of the brain are divided into an entire ROI. A previous study (Draganski et al., 2006) has demonstrated thalamic structural variations due to a chronic lack of afferent input in LLAs, reflecting the role of the thalamus in structural reorganization and functional restoration after injury to the central and peripheral nervous system. Therefore, we considered the thalamus as an independent ROI. In this study, the WFU_Pickatlas toolbox3 was used to define the thalamus as an independent ROI, and the remaining part in the original region from the atlas was defined as the ROI of the basal ganglia (Guan et al., 2017). All the ROIs of the subcortical and cortical sensorimotor network were visualized with BrainNet Viewer4 (Xia et al., 2013) (Figure 1). For simplicity, in the following sections, the brain region contralateral to the side of the amputation was called the contralateral region, and the brain region ipsilateral to the side of the amputation was called the ipsilateral region in LLAs. For example, the S1M1 contralateral to the side of the amputation was called the contralateral S1M1. ROI-wise connectivity was measured by extracting the residual time course data from each ROI listed above and calculating the Pearson’s correlation coefficients between each ROI pair. The resulting correlation coefficients were then Fisher-transformed into ’Z’ scores to increase normality.
FIGURE 1. Visual depiction of the regions of interest (ROIs) within the subcortical and cortical sensorimotor network that were used in the analyses described in the text. BG, basal ganglia; THA, thalamus; S1M1, the primary sensorimotor cortex; SMA, the supplementary motor area; I indicates ipsilateral to the side of the amputation in the lower-limb amputees (LLAs); C indicates contralateral to the side of the amputation in the LLAs. I corresponds to the left brain in the healthy controls (HCs), and the C corresponds to the right brain in the HCs.
To further elaborate inter-FC of the sensorimotor network, we calculated seed-based whole-brain voxel-wise FC (Jung et al., 2013; Fan et al., 2017). The contralateral S1M1 in the LLAs was selected as the seed because it is the most prominent brain region for functional reorganization following amputation (MacIver et al., 2008; Simões et al., 2012; Raffin et al., 2016). The whole-brain resting-state FC maps were obtained for the right S1M1 in HCs and the contralateral S1M1 in patients using the DPARSF toolbox.
Statistical Analysis
Analysis of Demographic and Clinical Characteristics
Two-tailed two sample t-tests and chi-square tests (only for sex) were performed to compare the demographic and clinical data from the two groups (SPSS 20.0; SPSS, Inc., Chicago, IL, United States). The significance level was set at p < 0.05.
Analyses of FC
Analyses of the ROI-wise connectivity were performed using the GRETNA toolbox5. Two-tailed two-sample t-tests were conducted to identify the specific group differences (p < 0.05, FDR corrected). Age, sex, and years of education were controlled as covariates to ensure that our results only revealed FC differences uniquely associated with amputation.
For the FC maps (seed-based whole-brain voxel-wise connectivity), two-tailed two-sample t-tests were performed to examine group-related differences between the LLAs and the HCs. Significance thresholds for t-test were set at p < 0.05 (FDR corrected). The covariates age, sex, and educational level were controlled. These statistical analyses were carried out with the SPM toolbox (spm86).
Relationships Between FC and Clinical Characteristics
We correlated clinical characteristics with brain FC, corresponding to the significant group-related difference between the LLAs and the HCs. A two-tailed partial correlation analysis was used after controlling for age, sex, educational level as confounding variables (p < 0.05).
Results
Demographic and Clinical Characteristics
Table 1 shows the demographic and clinical characterization of patients with an amputation and the HCs. No significant differences in age, sex, or education were observed between the two groups (p > 0.05). PLP was present in 12 LLAs [ranging in pain intensity from 1.1 to 6.8 (3.46 ± 2.19)].
ROI-Wise Connectivity Within the Sensorimotor Network
For ROI-wise connectivity, compared with the HCs, the LLAs exhibited significantly decreased FC between the contralateral S1M1 and the bilateral basal ganglia, bilateral thalamus, and cerebellum, between the cerebellum and the bilateral basal ganglia and bilateral thalamus, and between the SMA and the ipsilateral thalamus and contralateral basal ganglia (p < 0.05, FDR corrected) (Figure 2 and Table 2).
FIGURE 2. Significant between-group differences in the sensorimotor network in the intensity of the brain FC. The brain networks were visualized with the BrainNet Viewer. The thickness of the line represents the intensity of the between-group difference. BG, basal ganglia; THA, thalamus; S1M1, the primary sensorimotor cortex; SMA, the supplementary motor area; I indicate ipsilateral to the side of the amputation in the LLAs; C indicates contralateral to the side of the amputation in LLAs.
TABLE 2. Regions of interest (ROI)-wise FC differences within the sensorimotor network between the LLAs and the HCs.
Seed-Based Whole-Brain FC
The seed-based analysis provides an opportunity to validate and extend the main results obtained in the ROI-wise connectivity analysis, as it allows us to identify which specific loci might underlie the changes in FC between the contralateral S1M1 and the brain regions within and beyond the sensorimotor network in LLAs. When comparing connectivity maps between amputees and the HCs, significantly decreased FC was found between the contralateral S1M1 seed and the ipsilateral cerebellum and the ipsilateral thalamus in the LLAs, which is consistent with the results of the ROI-wise connectivity analysis described above. In addition, the seed-based map in Figure 3 also shows decreased FC between the contralateral S1M1 and the ipsilateral superior frontal gyrus (medial) and the ipsilateral occipital lobe in LLAs (p < 0.05, FDR corrected) (Figure 3 and Table 3).
FIGURE 3. Decreased FC in patients with lower-limb amputation with the seed region in the contralateral S1M1, compared with connectivity in the HCs (p < 0.05, FDR corrected, voxels > 50). 
: Cerebellum VI; 
: Cerebellum Crus1; 
: Superior frontal gyrus (medial); 
: Thalamus; 
: Superior occipital Gyrus, Cuneus, and Middle occipital gyrus. All these brain regions located on the cerebral hemisphere ipsilateral to the side of the amputation. The blue-to-green patches show, on axial slices, the t statistic of the comparisons between FC in the LLAs and the HCs. Slice location (z) is displayed in Montreal Neurological Institute coordinates. I indicates ipsilateral to the side of the amputation; C indicates contralateral to the side of the amputation.
TABLE 3. Brain regions showing decreased FC with the seed in the contralateral S1M1 in LLAs.
To test for possible hemisphere-specific effects, all the calculations were performed on the patients with only right- or left-sided amputations separately. The results did not reveal hemisphere specific differences and therefore all the patients were included in one calculation with a virtual left side amputation in this study. Furthermore, to detect the possible effects of the PLP, we divided all the patients into two groups (LLA with PLP and LLA without PLP). We compared the ROI-wise and the seed-based FC between two groups and failed to find any significant differences between the two groups.
Correlations Between FC and Clinical Characteristics
For ROI-wise connectivity, as shown in Figure 4, partial correlation analysis revealed that the FC between the SMA and the contralateral basal ganglia and ipsilateral thalamus and between the contralateral S1M1 and the contralateral basal ganglia and ipsilateral thalamus was positively correlated with the time since amputation in LLAs (p < 0.05, uncorrected). Only the correlations between the time since amputation and FC between the SMA and the contralateral basal ganglia and ipsilateral thalamus survived after FDR correction for multiple comparisons (p < 0.05, FDR corrected). Meanwhile, partial correlation analysis showed that the extracted FC of the cluster (the contralateral S1M1 seed to the ipsilateral thalamus) identified in the group comparison of seed-based FC was significantly correlated with the time since amputation in LLAs (r = 0.417, p = 0.017).
FIGURE 4. Scatter plots depicting the correlations between FC and time since amputation (p < 0.05). Partial correlation coefficients (r) were corrected for age, sex, and education. Scatter plots is fitted with regression line (black line). ∗The correlations survived a critical FDR threshold. BG, basal ganglia; SMA, the supplementary motor area; S1M1, the primary sensorimotor cortex; I indicates ipsilateral to the side of the amputation; C indicates contralateral to the side of the amputation.
Discussion
Accumulating evidence has indicated that amputation induces functional reorganization of the sensorimotor areas in the brain, including the S1M1, cerebellum, thalamus, and basal ganglia, which could be related to abnormal functioning of the subcortical-cortical sensorimotor system resulting from the reduction or absence of afferent information following amputation (Draganski et al., 2006; Romero-Romo et al., 2010; Simões et al., 2012; Raffin et al., 2016; Di Vita et al., 2018). In addition, structural and functional reorganization of the brain has also been found in the visual stream (Preißler et al., 2013; Jiang et al., 2015) and the DMN (Makin et al., 2015), which provides evidence of brain reorganization beyond the sensorimotor system. However, there have been few studies on patterns of reorganization in brain functional networks, especially in patients with lower-limb amputation. Moreover, previous studies have found progressive changes in brain structure and function after amputation (Jiang et al., 2015, 2016; Li et al., 2017), which suggests the need for further exploration of the relationship between imaging indicators and time since amputation.
In the current study, we investigated changes in the organization of the intra- and inter-FC of the sensorimotor network in unilateral LLAs, and the relationship between FC and the time since amputation. We identified abnormal FC within and beyond the sensorimotor network. In addition, we found positive correlations between FC among subcortical and cortical sensorimotor regions and the time since amputation.
Altered FC Within the Sensorimotor Network
As a drastic limb injury, amputation in humans induces changes in the S1M1 of the deafferent hemisphere (Ramachandran et al., 1992; Lotze et al., 1999; Giraux et al., 2001; Karl et al., 2001; Grüsser et al., 2004; Simões et al., 2012; Raffin et al., 2016). In this study, decreased FC among the contralateral S1M1 and the subcortical motor-related brain regions was found in LLAs by using an ROI-wise connectivity analysis. Our findings provide further evidence of network-level functional reorganization of the S1M1. We also found abnormal FC among the ipsilateral thalamus, the contralateral basal ganglia and the SMA in LLAs. To our knowledge, no prior study has reported the SMA function reorganization in amputees. A previous study found that the motor imagery of both normal and amputated limb can activate the SMA, and there was no significant difference in the intensity and range of activation in upper-extremity amputees (Diers et al., 2010). As an important part of the sensorimotor network, the SMA regulates the balance of the bilateral body and controls motion (Penfield and Welch, 1951; Serrien et al., 2002), playing an especially important role in controlling the movement of the lower limb (Luft et al., 2002; Sahyoun et al., 2004). Previous studies showed that SMA activation was prominent for parts of lower limb, such as the knee (Luft et al., 2002), and brain activation of the SMA was associated with gait control (Sahyoun et al., 2004). Therefore, changes in the FC of the SMA may be a unique pattern of brain reorganization in patients with lower-limb amputation.
In addition to cortical areas, the sensorimotor network consists of subcortical areas such as the cerebellum, basal ganglia and thalamus (Middleton and Strick, 2000; Dum and Strick, 2002). A previous study suggested that the cerebellum plays a role in the generation of appropriate patterns of limb movements, dynamic regulation of balance, and adaptation of posture and locomotion (Morton and Bastian, 2004). In addition, the cerebellum and the basal ganglia, as the major centers of the extra pyramidal motor system, exert their regulatory influence on the motor cortex via the thalamus to ‘shape’ the final automatic and voluntary output movement (Herrero et al., 2002; Groenewegen, 2003; Mizelle et al., 2016). All these findings reflect the essential role of these subcortical areas in the organization of limb movements. The structural and functional reorganization in these subcortical nuclei have been corroborated in neuroimaging studies of patients who have undergone lower-limb amputation. Recent work reported cerebellar gray matter modifications in LLAs (Di Vita et al., 2018), and a lower volume in the thalamus was found in subjects with limb amputation (Draganski et al., 2006). Moreover, Romero-Romo et al. (2010) found a disrupted basal ganglia-thalamus-cortex pathway in LLAs when the amputee group performed imaginary movements of amputated toes. In our study, we chose these subcortical areas as the ROIs to study the sensorimotor network in the LLAs and found decreased FC between the subcortical nuclei and between the subcortical nuclei and the SMA and the contralateral S1M1. The changes of FC in the brains of LLAs revealed the abnormal subcortical-cortical sensorimotor network after amputation, further highlighting the important role of subcortical areas in the reorganization of the sensorimotor network after lower-limb amputation.
Altered FC Between the S1M1 and the Other Parts of the Brain
The seed-based analysis also revealed aberrant connectivity within the sensorimotor network in LLAs. We found that the FC between the contralateral S1M1 and the ipsilateral cerebellum and the ipsilateral thalamus was reduced. These results confirmed the functional reorganization of the sensorimotor network in the patients with lower-limb amputation. In addition, we found reduced FC between the contralateral S1M1 and the ipsilateral superior frontal gyrus (medial) in the LLAs. The superior frontal gyrus (medial) is a part of the DMN (Greicius et al., 2003; Buckner et al., 2008; Raichle, 2015), which is involved with the representation of internalized sensations. Abnormal FC between the region of the cortex associated with a missing hand and the DMN was found in patients with arm amputation, which was suggested to correlate with complex perceptual experiences associated with phantom sensations (Makin et al., 2015). Meanwhile, we found reduced FC between the contralateral S1M1 and the ipsilateral occipital lobe, including the superior occipital gyrus, cuneus and middle occipital gyrus, which are important for visual processing. Previous fMRI studies found that amputation of the hand leads to shifts in visuospatial perception (Makin et al., 2010), as well-reductions in visuomotor object affordances (Wilf et al., 2013). Moreover, decreased cortical thickness was also reported in visual areas in LLAs (Jiang et al., 2015). Indeed, in the healthy brain, visual information of the limb is integral for motor processing during movement (Makin et al., 2014). The reduction of FC between visual and motor areas proved that in addition to visual cortical abnormalities, there was also an abnormal visual-motor information processing loop in LLAs. These findings demonstrated that functional consequences of amputation may result from reorganization at a network-level scale beyond the sensorimotor network. In this study, all the patients had been fitted with prostheses. The confounding effects of prosthesis use should be considered in the context of these results, as they may be associated with the brain functional reorganization of the brain in LLAs.
The Associations Between the Time Since Amputation and the FC
Our correlation analyses revealed statistically significant positive associations between the time since amputation and the FC in the contralateral basal ganglia, ipsilateral thalamus and SMA, and the contralateral basal ganglia, ipsilateral thalamus and contralateral S1M1 in the LLAs, which suggested that there were progressive subcortical-cortical FC changes in LLAs. Recent MRI studies have demonstrated that the brain changes over time following amputation (Jiang et al., 2015, 2016; Li et al., 2017). However, the majority of these studies mainly focused on the structural abnormalities of the LLAs. Jiang et al. (2015) found that LLAs displayed significant FA reduction in the right inferior front occipital fasciculus and a reduced thickness of the motor-related visual cortex (Jiang et al., 2016), both of which were negatively correlated with the time since amputation. Li et al. (2017) found that mean FA values of the corpus callosum were significantly smaller in amputees and that values were correlated with the time since amputation. These previous findings suggest a progressively disrupted pattern in the brain gray matter and white matter over the time following amputation. In this study, we found progressive increases in the reductions of connectivity between the brain subcortical-cortical regions in the LLAs, which could be a form of compensation for the increase in serious structural damage over time. Another possible explanation for this restoration of FC is that the functional network reorganization after amputation is more likely to be reversed/modulated by motion exercises or prosthesis use. Future rehabilitation treatments should pay attention to the progressive changes in the FC of the sensorimotor network, and the gradual progression of FC may be an important indicator for evaluating the rehabilitation of motor function in patients following lower-limb amputation.
Limitations
Despite of the important results showing brain functional network reorganization in the LLAs mentioned above, several limitations in this study should be further considered. First, in the current study, we mainly choose the sensorimotor network to study the functional reorganization process over time since amputation, which limits our ability to speculate on global brain organization patterns after amputation. Second, the sample size of the LLAs is relatively small. Studies with larger sample sizes are needed to clarify the specific functional reorganization mechanisms in patients with different types of amputations and to disentangle the influences of the amputation site, the amputation side and the use of prostheses on cortical reorganization in amputees. Third, previous studies found that PLP may affect brain functional reorganization (Flor et al., 2006) and could have an impact on perceived motor control in the amputees (Kikkert et al., 2017). However, we did not find significant differences between LLA with PLP and LLA without PLP in our study. It may be attributed to a small number of amputees with PLP and a greater number with low intensity PLP in this study. In future studies, a large number of amputees with different PLP levels are needed to further study the effects of PLP on brain functional reorganization in the LLAs. Meanwhile, phantom sensations are not assessed in our study. This information needs to be supplemented to investigate the effects of phantom sensations on brain functional reorganization in the LLAs in a subsequent study.
Conclusion
Our results revealed that lower-limb amputation leads to subcortical and cortical reorganization in the sensorimotor network. Moreover, between-group differences in the seed-based whole-brain FC analysis showed that brain regions with decreased FC with the contralateral S1M1 extended to the prefrontal [superior frontal gyrus (medial)] and occipital cortices. These findings suggested that functional reorganization after lower-limb amputation occurred at multiple scales, within and beyond the scope of the sensorimotor network. Finally, correlation analysis results showed that decreased FC between the subcortical and the cortical regions (the SMA and the contralateral S1M1) progressively increased over the time since amputation, which reflects the process of brain functional network recovery, probably related to motion exercises or prosthesis use after amputation. All of these results hold promise for developing reliable biomarkers of FC changes related to motor function recovery in LLAs.
Ethics Statement
We declare that all human studies have been approved by the Army Medical University Ethics Committee and have therefore been performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments. We declare that all patients gave informed consent prior to inclusion in this study.
Author Contributions
JZ designed and performed the experiments, analyzed the data, and wrote the manuscript. LW and LS assisted in fMRI data analysis. YZ assisted in writing and editing the manuscript. LL and PL assisted in sample collection and clinical data analysis. XY and MQ provided overall project supervision, helped to analyze and interpret the data, and assisted in writing the manuscript. All authors have read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China for Young Scholars (81601970 and 81301205), the translational medicine cooperation project of The First Hospital Affiliated to Third Military Medical University (SWH2014ZH03), and the youth development program of Third Military Medical University (2016XPY05 and 2017XQN04).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer MS and handling Editor declared their shared affiliation.
Acknowledgments
We thank all the foundations for financial support. We are grateful for the generosity of time and effort by all the participants and researchers who made this project possible.
Footnotes
- ^ http://rfmri.org/DPARSF
- ^ http://findlab.stanford.edu/functional_ROIs.html
- ^ https://www.nitrc.org/projects/wfu_pickatlas
- ^ http://www.nitrc.org/projects/bnv/
- ^ https://www.nitrc.org/projects/gretna/
- ^ http://www.fil.ion.ucl.ac.uk/spm/software/spm8/
References
Barbosa, B. M., Monteiro, R. A., Sparano, L. F., Bareiro, R. F., Passos, A. D., and Engel, E. E. (2016). Incidence and causes of lower-limb amputations in the city of Ribeirão Preto from 1985 to 2008: evaluation of the medical records from 3,274 cases. Rev Bras Epidemiol 19, 317–325. doi: 10.1590/1980-5497201600020009
Bostan, A. C., Dum, R. P., and Strick, P. L. (2013). Cerebellar networks with the cerebral cortex and basal ganglia. Trends. Cogn. Sci. 17, 241–254. doi: 10.1016/j.tics.2013.03.003
Bostan, A. C., and Strick, P. L. (2010). The cerebellum and basal ganglia are interconnected. Neuropsychol. Rev. 20, 261–270. doi: 10.1007/s11065-010-9143-9
Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. N. Y. Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011
Collignon, O., Dormal, G., Albouy, G., Vandewalle, G., Voss, P., Phillips, C., et al. (2013). Impact of blindness onset on the functional organization and the connectivity of the occipital cortex. Brain 136, 2769–2783. doi: 10.1093/brain/awt176
Curt, A., Yengue, C. N., Hilti, L. M., and Brugger, P. (2011). Supernumerary phantom limbs in spinal cord injury. Spinal Cord 49, 588–595. doi: 10.1038/sc.2010.143
Di Vita, A., Boccia, M., Palermo, L., Nemmi, F., Traballesi, M., Brunelli, S., et al. (2018). Cerebellar grey matter modifications in lower limb amputees not using prosthesis. Sci. Rep. 8:370. doi: 10.1038/s41598-017-18772-2
Diers, M., Christmann, C., Koeppe, C., Ruf, M., and Flor, H. (2010). Mirrored, imagined and executed movements differentially activate sensorimotor cortex in amputees with and without phantom limb pain. Pain 149, 296–304. doi: 10.1016/j.pain.2010.02.020
Draganski, B., Moser, T., Lummel, N., Gänssbauer, S., Bogdahn, U., Haas, F., et al. (2006). Decrease of thalamic gray matter following limb amputation. Neuroimage 31, 951–957. doi: 10.1016/j.neuroimage.2006.01.018
Dum, R. P., and Strick, P. L. (2002). Motor areas in the frontal lobe of the primate. Physiol. Behav. 77, 677–682. doi: 10.1016/S0031-9384(02)00929-0
Fan, J., Zhong, M., Zhu, X., Gan, J., Liu, W., Niu, C., et al. (2017). Resting-state functional connectivity between right anterior insula and right orbital frontal cortex correlate with insight level in obsessive-compulsive disorder. Neuroimage Clin. 15, 1–7. doi: 10.1016/j.nicl.2017.04.002
Flor, H., Nikolajsen, L., and Staehelin Jensen, T. (2006). Phantom limb pain: a case of maladaptive CNS plasticity? Nat. Rev. Neurosci. 7, 873–881. doi: 10.1038/nrn1991
Garraghty, P. E., and Kaas, J. H. (1991). Functional reorganization in adult monkey thalamus after peripheral nerve injury. Neuroreport 2, 747–750. doi: 10.1097/00001756-199112000-00004
Giraux, P., Sirigu, A., Schneider, F., and Dubernard, J. M. (2001). Cortical reorganization in motor cortex after graft of both hands. Nat. Neurosci. 4, 691–692. doi: 10.1038/89472
Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253–258. doi: 10.1073/pnas.0135058100
Groenewegen, H. J. (2003). The basal ganglia and motor control. Neural. Plast. 10, 107–120. doi: 10.1155/NP.2003.107
Grüsser, S. M., Mühlnickel, W., Schaefer, M., Villringer, K., Christmann, C., Koeppe, C., et al. (2004). Remote activation of referred phantom sensation and cortical reorganization in human upper extremity amputees. Exp. Brain. Res. 154, 97–102. doi: 10.1007/s00221-003-1649-4
Guan, X., Zeng, Q., Guo, T., Wang, J., Xuan, M., Gu, Q., et al. (2017). Disrupted functional connectivity of basal ganglia across tremor-dominant and akinetic/rigid-dominant parkinson’s disease. Front. Aging Neurosci. 9:360. doi: 10.3389/fnagi.2017.00360
Herrero, M. T., Barcia, C., and Navarro, J. M. (2002). Functional anatomy of thalamus and basal ganglia. Childs Nerv. Syst. 18, 386–404. doi: 10.1007/s00381-002-0604-1
Hoshi, E., Tremblay, L., Féger, J., Carras, P. L., and Strick, P. L. (2005). The cerebellum communicates with the basal ganglia. Nat. Neurosci. 8, 1491–1493. doi: 10.1038/nn1544
Jain, N., Qi, H. X., Collins, C. E., and Kaas, J. H. (2008). Large-scale reorganization in the somatosensory cortex and thalamus after sensory loss in macaque monkeys. J. Neurosci. 28, 11042–11060. doi: 10.1523/JNEUROSCI.2334-08.2008
Jiang, G., Li, C., Wu, J., Jiang, T., Zhang, Y., Zhao, L., et al. (2016). Progressive thinning of visual motion area in lower limb amputees. Front. Hum. Neurosci. 10:79. doi: 10.3389/fnhum.2016.00079
Jiang, G., Yin, X., Li, C., Li, L., Zhao, L., Evans, A. C., et al. (2015). The plasticity of brain gray matterand white matter following lower limb amputation. Neural Plast. 2015:823185. doi: 10.1155/2015/823185
Jones, E. G., and Pons, T. P. (1998). Thalamic and brainstem contributions to large-scale plasticity of primate somatosensory cortex. Science 282, 1121–1125. doi: 10.1126/science.282.5391.1121
Jung, W. H., Kang, D. H., Kim, E., Shin, K. S., Jang, J. H., and Kwon, J. S. (2013). Abnormal corticostriatal-limbic functional connectivity in obsessive–compulsive disorder during reward processing and resting-state. Neuroimage Clin. 3, 27–38. doi: 10.1016/j.nicl.2013.06.013
Kaas, J. H., Merzenich, M. M., and Killackey, H. P. (1983). The reorganization of somatosensory cortex following peripheral nerve damage in adult and developing mammals. Annu. Rev. Neurosci. 6, 325–356. doi: 10.1146/annurev.ne.06.030183.001545
Karl, A., Birbaumer, N., Lutzenberger, W., Cohen, L. G., and Flor, H. (2001). Reorganization of motor and somatosensory cortex in upper extremity amputees with phantom limb pain. J. Neurosci. 21, 3609–3618. doi: 10.1523/JNEUROSCI.21-10-03609.2001
Kikkert, S., Mezue, M., Henderson Slater, D., Johansen-Berg, H., Tracey, I., and Makin, T. R. (2017). Motor correlates of phantom limb pain. Cortex 95, 29–36. doi: 10.1016/j.cortex.2017.07.015
Li, Z., Li, C., Fan, L., Jiang, G., Wu, J., Jiang, T., et al. (2017). Altered microstructure rather thanmorphology in the corpuscallosum after lower limb amputation. Sci. Rep. 7:44780. doi: 10.1038/srep44780
Lotze, M., Grodd, W., Birbaumer, N., Erb, M., Huse, E., and Flor, H. (1999). Does use of a myoelectric prosthesis prevent cortical reorganization and phantom limb pain? Nat. Neurosci. 2, 501–502. doi: 10.1038/9145
Luft, A. R., Smith, G. V., Forrester, L., Whitall, J., Macko, R. F., Hauser, T. K., et al. (2002). Comparing brain activation associated with isolate dupper and lower limb movement across corresponding joints. Hum. Brain. Mapp. 7, 131–140. doi: 10.1002/hbm.10058
MacIver, K., Lloyd, D. M., Kelly, S., Roberts, N., and Nurmikko, T. (2008). Phantom limb pain, cortical reorganization and the therapeutic effect of mentalimagery. Brain 131, 2181–2191. doi: 10.1093/brain/awn124
Makin, T. R., Brozzoli, C., Cardinali, L., Holmes, N. P., and Farnè, A. (2014). Left or right? rapid visuomotor coding of hand laterality during motor decisions. Cortex 64, 289–292. doi: 10.1016/j.cortex.2014.12.004
Makin, T. R., Filippini, N., Duff, E. P., Henderson Slater, D., Tracey, I., and Johansen-Berg, H. (2015). Network-level reorganization of functional connectivity following armamputation. Neuroimage 114, 217–225. doi: 10.1016/j.neuroimage.2015.02.067
Makin, T. R., Wilf, M., Schwartz, I., and Zohary, E. (2010). Amputees “neglect” the space near their missing hand. Psychol. Sci. 21, 55–57. doi: 10.1177/0956797609354739
Merzenich, M. M., Kaas, J. H., Wall, J. T., Sur, M., Nelson, R. J., and Felleman, D. J. (1983). Progression of change following median nerve section in the cortical representation of the hand in areas 3b and 1 in adult owl and squirrel monkeys. Neuroscience 10, 639–665. doi: 10.1016/0306-4522(83)90208-7
Middleton, F. A., and Strick, P. L. (2000). Basal ganglia and cerebellar loops: motor and cognitive circuits. Brain Res. Rev. 31, 236–250. doi: 10.1016/S0165-0173(99)00040-5
Mizelle, J. C., Oparah, A., and Wheaton, L. A. (2016). Reliability of visual and somatosensory feedback in skilled movement: the role of the cerebellum. Brain Topogr. 29, 27–41. doi: 10.1007/s10548-015-0446-2
Molina Rueda, F., Alguacil Diego, I. M., Molero Sánchez, A., Carratalá Tejada, M., Rivas Montero, F. M., and Miangolarra Page, J. C. (2013). Knee and hip internal moments and upper-body kinematics in the frontal plane in unilateral transtibial amputees. Gait Posture 37, 436–439. doi: 10.1016/j.gaitpost.2012.08.019
Morton, S. M., and Bastian, A. J. (2004). Cerebellar control of balance and locomotion. Neuroscientist 10, 247–259. doi: 10.1177/1073858404263517
Napadow, V., LaCount, L., Park, K., As-Sanie, S., Clauw, D. J., and Harris, R. E. (2010). Intrinsic brain connectivity in fibromyalgia is associated with chronic pain intensity. Arthritis Rheum. 62, 2545–2555. doi: 10.1002/art.27497
Oldfield, R. (1971). The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9, 97–113. doi: 10.1016/0028-3932(71)90067-4
Penfield, W., and Welch, K. (1951). The supplementary motor area of the cerebral cortex: a clinical and experimental study. Am. Med. Assoc. Arch. Neurol. Psychiatr. 66, 289–317. doi: 10.1001/archneurpsyc.1951.02320090038004
Pons, T. P., Garraghty, P. E., Ommaya, A. K., Kaas, J. H., Taub, E., and Mishkin, M. (1991). Massive cortical reorganization after sensory deafferentation in adult macaques. Science 252, 1857–1860. doi: 10.1126/science.1843843
Preißler, S., Dietrich, C., Blume, K. R., Hofmann, G. O., Miltner, W. H., and Weiss, T. (2013). Plasticity in the Visual System is Associated with Prosthesis Use in Phantom Limb Pain. Front. Hum. Neurosci. 7:311. doi: 10.3389/fnhum.2013.00311
Raffin, E., Richard, N., Giraux, P., and Reilly, K. T. (2016). Primary motor cortex changes after amputation correlate with phantom limb pain and the ability to move the phantom limb. Neuroimage 130, 134–144. doi: 10.1016/j.neuroimage.2016.01.063
Raichle, M. E. (2015). The brain’s default mode network. Annu. Rev. Neurosci. 38, 433–447. doi: 10.1146/annurev-neuro-071013-014030
Ramachandran, V. S., Rogers-Ramachandran, D., and Stewart, M. (1992). Perceptual correlates of massive cortical reorganization. Science 258, 1159–1160. doi: 10.1126/science.1439826
Rocca, M. A., Valsasina, P., Fazio, R., Previtali, S. C., Messina, R., Falini, A., et al. (2014). Brain connectivity abnormalities extend beyond the sensorimotor network in peripheral neuropathy. Hum. Brain Mapp. 35, 513–526. doi: 10.1002/hbm.22198
Romero-Romo, J. I., Bauer, C. C., Pasaye, E. H., Gutiérrez, R. A., Favila, R., and Barrios, F. A. (2010). Abnormal functioning of the thalamocortical system underlies the conscious awareness of the phantom limb phenomenon. Neuroradiol. J. 23, 671–679. doi: 10.1177/197140091002300605
Rusaw, D., and Ramstrand, N. (2011). Motion-analysis studies of transtibial prosthesis users: a systematic review. Prosthet. Orthot. Int. 35, 8–19. doi: 10.1177/0309364610393060
Sahyoun, C., Floyer-Lea, A., Johansen-Berg, H., and Matthews, P. M. (2004). Towards an understanding of gait control:brain activation during the anticipation, preparation and execution of foot movements. Neuroimage 21, 568–575. doi: 10.1016/j.neuroimage.2003.09.065
Satterthwaite, T. D., Elliott, M. A., Gerraty, R. T., Ruparel, K., Loughead, J., and Calkins, M. E. (2013). An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. Neuroimage 64, 240–256. doi: 10.1016/j.neuroimage.2012.08.052
Scott, J., and Huskisson, E. C. (1976). Graphic representation of pain. Pain 2, 175–184. doi: 10.1016/0304-3959(76)90113-5
Serrien, D. J., Strens, L. H., Oliveiero, A., and Brown, P. (2002). Repetitive transcranial magnetic stimulation of the supplementary motor area (SMA) degrades bimanual movement control in humans. Neurosci. Lett. 328, 89–92. doi: 10.1016/S0304-3940(02)00499-8
Shirer, W. R., Ryali, S., Rykhlevskaia, E., Menon, V., and Greicius, M. D. (2012). Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb. Cortex 22, 158–165. doi: 10.1093/cercor/bhr099
Simões, E. L., Bramati, I., Rodrigues, E., Franzoi, A., Moll, J., and Lent, R. (2012). Functional expansion of sensorimotor representation and structural reorganization of callosal connections in lower limb amputees. J. Neurosci. 32, 3211–3220. doi: 10.1523/JNEUROSCI.4592-11.2012
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Smith, S. M., Miller, K. L., Moeller, S., Xu, J., Auerbach, E. J., Woolrich, M. W., et al. (2012). Temporally-independent functional modes of spontaneous brain activity. Proc. Natl. Acad. Sci. U.S.A. 109, 3131–3136. doi: 10.1073/pnas.1121329109
Soares, A. S., Yamaguti, E. Y., Mochizuki, L., Amadio, A. C., and Serra, J. C. (2009). Biomechanical parameters of gait among transtibial amputees: a review. Sao Paulo Med. J. 127, 302–309. doi: 10.1590/S1516-31802009000500010
Wilf, M., Holmes, N. P., Schwartz, I., and Makin, T. R. (2013). Dissociating between object affordances and spatial compatibility effects using early response components. Front. Psychol. 4:591. doi: 10.3389/fpsyg.2013.00591
Xia, M., Wang, J., and He, Y. (2013). BrainNet viewer: a network visualization tool for human brain connectomics. PLoS One 8:e68910. doi: 10.1371/journal.pone.0068910
Yan, C. G., Cheung, B., Kelly, C., Colcombe, S., Craddock, R. C., Di Martino, A., et al. (2013a). A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. Neuroimage 76, 183–201. doi: 10.1016/j.neuroimage.2013.03.004
Yan, C. G., Craddock, R. C., Zuo, X. N., Zang, Y. F., and Milham, M. P. (2013b). Standardizing the intrinsic brain: towards robust measurement of inter-individual variation in 1000 functional connectomes. Neuroimage 80, 246–262. doi: 10.1016/j.neuroimage.2013.04.081
Keywords: lower-limb amputation, functional reorganization, resting-state functional connectivity, the sensorimotor network, S1M1, subcortical nuclei
Citation: Zhang J, Zhang Y, Wang L, Sang L, Li L, Li P, Yin X and Qiu M (2018) Brain Functional Connectivity Plasticity Within and Beyond the Sensorimotor Network in Lower-Limb Amputees. Front. Hum. Neurosci. 12:403. doi: 10.3389/fnhum.2018.00403
Received: 29 May 2018; Accepted: 20 September 2018;
Published: 09 October 2018.
Edited by:
Peter Sörös, University of Oldenburg, GermanyReviewed by:
Melanie Spindler, University of Oldenburg, GermanyTakufumi Yanagisawa, Osaka University, Japan
Copyright © 2018 Zhang, Zhang, Wang, Sang, Li, Li, Yin and Qiu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xuntao Yin, xuntaoyin@gmail.com Mingguo Qiu, qiumingguo@yahoo.com