Abstract

Purpose. This study was aimed at evaluating the motor cortical excitability and connectivity underlying the neural mechanism of motor deficit in acute stroke by the combination of functional magnetic resonance imaging (fMRI) and electrophysiological measures. Methods. Twenty-five patients with motor deficit after acute ischemic stroke were involved. General linear model and dynamic causal model analyses were applied to fMRI data for detecting motor-related activation and effective connectivity of the motor cortices. Motor cortical excitability was determined as a resting motor threshold (RMT) of motor evoked potential detected by transcranial magnetic stimulation (TMS). fMRI results were correlated with cortical excitability and upper extremity Fugl-Meyer assessment scores, respectively. Results. Greater fMRI activation likelihood and motor cortical excitability in the ipsilesional primary motor area (M1) region were associated with better motor performance. During hand movements, the inhibitory connectivity from the contralesional to the ipsilesional M1 was correlated with the degree of motor impairment. Furthermore, ipsilesional motor cortex excitability was correlated with an enhancement of promoting connectivity in ipsilesional M1 or a reduction of interhemispheric inhibition in contralesional M1. Conclusions. The study suggested that a dysfunction of the ipsilesional M1 and abnormal interhemispheric interactions might underlie the motor disability in acute ischemic stroke. Modifying the excitability of the motor cortex and correcting the abnormal motor network connectivity associated with the motor deficit might be the therapeutic target in early neurorehabilitation for stroke patients.

1. Introduction

Stroke is a leading cause of long-term disability worldwide [1]; approximately 80% of acute stroke survivors suffer from persistent different degrees of motor impairment [2]. Research on the neural mechanism of stroke-induced motor impairment is essential for the rehabilitation of the disease. Currently, two related hypotheses have been offered to explain motor impairment of stroke patients. In the first theory, dysfunction of ipsilesional motor areas underlies deficits in motor control after stroke [3, 4], and another posits that transcallosal inhibition (e.g., the contralesional hemisphere exerts a pathologically enhanced inhibitory influences on the ipsilesional hemisphere) thereby contributes to the reduced motor performance of stroke patients [57]. In the acute phase of stroke, changes in the brain at different levels from the molecular level to the network level are especially active and accelerate the brain’s self-repair processes [8, 9]. Earlier rehabilitation in stroke is more critical for motor functional recovery and has been considered to be associated with better outcomes [10]. Therefore, it is more significant to investigate the neural mechanisms underlying motor deficits in the early stage of stroke, which may provide important theoretical evidence for early neurorehabilitation of ischemic stroke.

Over recent years, multimodality imaging techniques, represented by functional magnetic resource imaging (fMRI), have been widely used for investigating the neural mechanism of the motor deficits in ischemic stroke. By using sophisticated analysis approaches, fMRI has revealed the abnormalities of local activation and connectivity associated with motor deficits after stroke [11]. Dynamic causal modeling (DCM) analysis based on fMRI is a promising approach to capture effective connectivity among key regions of the cortical motor system [12]. As a complementary approach, electrophysiological assessment such as the assessment of the motor evoked potentials (MEPs) by transcranial magnetic stimulation (TMS) provides information about the excitability of the corticomotor system in stroke patients [13]. Nonetheless, these numerous studies used electrophysiological or fMRI measures in isolation. Moreover, most available evidence was primarily limited to the chronic stage after stroke, in which the mechanisms underlying neural reorganization may have returned to being stable [14]. There is little specific information on the characteristics of the motor cortical network in the acute stage following stroke.

Therefore, in the present study, we combined functional MRI and electrophysiological assessments, as well as analyses of effective connectivity to investigate the neural mechanisms underlying motor deficits in acute stroke. We first hypothesized that the early movement-related cortical activation patterns were paralleled by the motor cortical excitability measured by TMS and both relates to the degree of motor impairment after stroke. Then, combining DCM with electrophysiological assessment, we further hypothesized that pathological transcallosal inhibition originating from contralesional M1 to ipsilesional M1 could also occur in the acute stage of stroke and relate to the degree of motor impairment and motor cortical excitability.

2. Materials and Methods

2.1. Patients

Twenty-five patients (8 females and 17 males; age: 53 ± 10 years) with a first ever ischemic stroke participated in this study. The detailed inclusive criteria included the following: (1) patients being in the acute stage of stroke and having a time from stroke of less than two weeks, (2) single lesion located within the middle cerebral artery territory, as shown by structural MRI, (3) symptoms of motor deficit in the unilateral upper extremity, and (4) absence of aphasia, neglect, apraxia, or epilepsy, which were evaluated according to the patients’ clinical symptoms and signs. This study was approved by the Internal Review Board of Jinling Hospital, and written informed consent was obtained from each participant.

2.2. Clinical Assessments

The clinical deficit and stroke severity were assessed using the National Institutes of Health Stroke Scale (NIHSS) [15] and the modified Rankin Scale (mRS) [16]. The following motor tests were used to assess the affected upper limb: (1) Fugl-Meyer assessment (FMA) [17] and (2) Medical Research Council (MRC) scale [18] for wrist extension of the affected hand. The FMA is a reliable and standardized motor impairment scale, which ranges from 0 (hemiplegia) to a maximum of 100 points (normal motor performance). It is divided into 66 points for the upper limb (FMA-UL) and 34 points for the lower limb (FMA-LL), scored on a 3-point ordinal scale (0 = cannot perform, 1 = performs partially, and 2 = performs fully). MRC is also used to assess the muscle strength of the hemiplegic side (ranging from 0 to 5; 5 = normal power and 0 = no movement).

2.3. Determination of Motor Cortex Excitability

Motor cortex excitability in each hemisphere was measured using single-pulse TMS [19, 20]. TMS was performed using a MagproX100 stimulator (MagVenture Company, Farum, Denmark) with a figure-of-eight coil (outer diameter of one wing 9 cm). Electromyography (EMG) recordings from the bilateral abductor pollicis brevis (APB) were detected using two pairs of silver-silver chloride surface electrodes (Alpine Biomed ApS, Skovlunde, Denmark). The RMT was defined as the minimal stimulus intensity, which could produce a MEP with a peak-to-peak amplitude of at least 50 μV in at least 5 of 10 subsequent trials. The amplitude of MEP to cortical stimulation was measured according to the peak-to-peak voltage (mV) of the EMG response. The latency of MEP was determined as the interval between the onset of the stimulation and the onset of the EMG response (ms). The presence of a TMS-elicited MEP was categorized as present or absent.

2.4. fMRI Paradigm

A block-designed motor task, for which the thumb and the forefinger of the affected and unaffected hand had to be tapped together at a rate of 1 Hz, was used as an fMRI activation paradigm. The MRI session consists of five trainings, lasting 5 minutes in total. Each training contained three blocks: a 20-second block of movement of the affected hand, a 20-second block of movement of the unaffected hand, and a 20-second rest block. Two motor task blocks were pseudorandomized and separated by a rest block. Every block was repeated 5 times in total. An investigator carried out the passive movement in the scanner room and ensured the consistency of the motor tasks. All patients underwent motor task training before MRI scanning and performed the passive movements during task-evoked fMRI scanning.

2.5. fMRI Data Acquisition

Whole-brain fMRI scans were performed using a GE MR750 3.0 Tesla Scanner (General Electric, USA). We used a gradient echo planar imaging (EPI) sequence with the following imaging parameters: ; ; ; ; ; , no gap; and slice number = 43 and 160 volumes. The slices covered the whole brain extending from the frontoparietal cortex to lower parts of the cerebellum. We also used 3D-BRAVO sequence to generate high-resolution 3D T1-weighted structural images, with these parameters: , , , , , and .

2.6. MRI Data Preprocess

Considering the effect of stroke lesion on the spatial normalization, a cost function modification was used to remove the presence of a lesion [21]. First, the lesion masks were manually traced on the individual structural 3D T1-weighted images independently by 2 radiologists (Hu JP and Zeng FY), with the guide of the DWI image. Second, a group-sample-specific brain template was generated. For each subject, the lesion-removed whole-brain mask was used as the cost function, to normalize the 3D T1-weighted image into the standard brain template of the Montreal Neurological Institute (MNI) by using a 12-parameter affine transformation with nonlinear adjustments with 78 × 7 basis functions. All individual-normalized 3D T1-weighted anatomical images and all lesion-removed masks were averaged to yield sample-specific brain templates. Third, the averaged 3D T1-weighted anatomical image was segmented using unified segment function of SPM8, with the averaged lesion-removed mask as the cost function. Fourth, the individual space 3D T1-weighted images were segmented using unified segment function. The segment issues of the third step process were used as the templates of the current segment process, and the individual lesion-removed brain mask was used as the cost function. The segment parameters contained the affine transition matrix, which would be used in the functional MRI image normalization [21, 22].

For the functional MRI data, the process was done by DPARSF (http://rfmri.org/DPARSF) and SPM8 (http://www.fil.ion.ucl.ac.uk). Prior to data analysis, images from patients with right-sided lesions were flipped to the left side at the midsagittal plane, so that the affected hemisphere corresponded to the left hemisphere in all patients. First, slice timing and realignment were performed. Translation or rotation parameters in any given data set did not exceed 1.5 mm or 1.5 degree. Second, the functional images were coregistered to the individual 3D T1-weighted images. Third, the segment parameters were applied to the coregistered functional images for the normalization of functional images, with resample of 3 × 3 × 3 mm3 voxel size. Fourth, spatial smoothing () was performed to the normalized functional images [21, 22].

2.7. fMRI Data Analysis
2.7.1. Motor Activation Detection

For each subject, the experimental conditions were modeled in the framework of the general linear model (GLM). The block designs were convolved with a canonical hemodynamic response function applied in SPM8, and the first-order temporal derivatives were also generated. The design matrix also contained the head motion parameters to remove the effects. The contrast maps of the affected and unaffected hand movement were regrouped to calculate the second-level group effect by using SPM8.

2.7.2. Dynamic Causal Modeling

Dynamic causal modeling [23] was used to assess effective connectivity within bilateral cortical motor areas activated by the fMRI motor task. Connectivity parameters represent the connectivity strength from one brain region to another. According to the previous studies [5, 12], six motor regions of interest (ROIs) from the fMRI motor task activation were selected for the DCM analysis: the primary motor area (M1), the premotor cortex (PMC), and the supplementary motor area (SMA) from both hemispheres of each subject (Table 1). A 4 mm radius sphere was applied to the ROI extraction. Principal component analysis was used to extract the ROI that represents signals. Based on previous studies on anatomical and functional connectivity in animals and human beings [4, 5, 12, 24], we supposed the intrinsic connection of these motor areas: connections between SMA and ipsilateral and contralateral M1, between SMA and ipsilateral PMC, and between PMC and ipsilateral M1 and transcallosal connections between M1-M1 and SMA-SMA. Bayes factors were computed for the “winning” model providing the best fit between accuracy and generalizability.

2.8. Statistical Analysis

Statistical analyses were performed on SPSS for Windows version 22 (IBM Corp., Armonk, NY). The comparison of neurophysiological measures in the bilateral hemisphere was analyzed by a two-tailed paired -test, with a statistical significance set at . The estimated connectivity parameters of the “winning” model in patients were tested for statistical significance by means of one-sample t-tests (). Pearson correlation analyses were used to test the relationships between fMRI results (neural activity and connectivity parameters) and motor performance, as well as cortical excitability measures. The level of significance was set at .

3. Results

Table 2 shows the demographic, imaging, and clinical characteristics of these patients.

3.1. Cortical Excitability

In this study, MEPs could not be elicited in five of the 25 patients that were excluded from neurophysiological data analysis. We found a significant decreased cortical excitability in the ipsilesional hemisphere (increased RMT, reduced MEP amplitude, and prolonged MEP latency) compared to the contralesional hemisphere (, , and , resp., paired -test) (Figures 1(a)1(c)). The RMT of the ipsilesional hemisphere significantly correlated with FMA scores of the affected limb (, ) (Figure 1(d)), suggesting that better motor performance and less intensity were needed to evoke an MEP in ipsilesional M1. However, no significant correlation was evident between the RMT of the contralesional hemisphere and FMA scores of the affected limb (, ).

3.2. fMRI Motor Activation

We first examined the task-related activation during hand movements. Movements of the unaffected hand were associated with significant BOLD activity in the contralateral M1, SMA, bilateral PMC, and the contralateral somatosensory cortex, consistent with previous studies in healthy subjects [25]. In contrast, although movements of the affected hand revealed BOLD activity in areas homologous to those just mentioned, neural activity decreased in both hemispheres, especially within M1 of the ipsilesional hemisphere (Figure 2(a)).

We next tested for correlations between brain activation early after stroke and motor impairment as well as motor cortical excitability. In the analysis of the relationship between brain activation and motor performance, we ruled out three outliers that were the patients with FMA 2-3 (patients 5, 6, and 12). A significant positive correlation was evident between the levels of BOLD activation within ipsilesional M1 and FMA (, ), with better motor performance in patients featuring higher levels of BOLD activation in ipsilesional M1 during movement of the affected hand (Figure 2(b)). In 20 patients with the presence of MEP, we also found a significant correlation between cortical excitability of ipsilesional M1 and the BOLD activation in the ipsilesional sensorimotor area for movement of the affected hand (, ) (Figure 2(c)). In contrast, the level of the BOLD signal in the contralesional hemisphere did not correlate with motor function as well as with cortical excitability of either hemisphere.

3.3. DCM

Figure 3(a) shows the winning model for endogenous neural coupling. We first analyzed the hand-movement-associated effective connectivity among bilateral motor networks (Figure 3(b)). The interhemispheric connections between ipsilesional and contralesional M1 showed negative connectivity parameters, indicating a mutual suppression of M1-M1 interactions during movement of the affected or unaffected hand (, corrected). Furthermore, we also found an additional inhibitory connectivity from contralesional SMA to M1 during movement of the affected hand, as has been described in the previous study [26].

In the next step, we performed correlation analysis of the connectivity parameters and motor function of the affected upper limb. We found that the FMA correlated negatively with inhibitory connectivity from contralesional M1 to ipsilesional M1 (, ) and positively with inhibitory connectivity from contralesional SMA to M1 (, ) during movements of the affected hand. For movements of the unaffected hand, the FMA correlated with promoting connectivity from contralesional SMA to ipsilesional M1 (, ) as well as with the M1-M1 inhibition (reduced inhibition of ipsilesional M1) (, , resp.) (Figure 4(a)). These findings were consistent with our fMRI result showing that BOLD activity within the ipsilesional M1 correlated with better motor performance (Figure 2(b)).

Finally, when correlating DCM connectivity parameters with cortical excitability, we found that the stronger the promoting connectivity from ipsilesional PMC to M1 during movement of the affected hand was, the more active the motor cortex excitability was present in the ipsilesional hemisphere (, ; Figure 4(b)). In addition, we also found a statistically significant correlation between the M1-M1 inhibition and RMT of the ipsilesional hemisphere during movement of the unaffected hand (ipsilesional M1-contralesional M1, , ; contralesional M1-ipsilesional M1, , ) (Figure 4(b)). Hence, lower motor cortex excitability (higher intensities are needed to evoke MEP) in the ipsilesional hemisphere featured more intra- or interhemispheric inhibition of ipsilesional M1 and a relative disinhibition of contralesional M1.

4. Discussion

Combining neuroimaging and electrophysiological measures, the present study examined the relationships of effective connectivity with motor performance and motor cortical excitability in acute ischemic stroke. We first found that the brain region responding to the affected hand had weaker fMRI activation and cortical excitability relative to those of the unaffected hand. The fMRI activation and cortical excitability in the ipsilesional M1 during movement of the affected hand showed a positive correlation with motor performance. Importantly, DCM analysis revealed both enhancement of positive connectivity in the ipsilesional M1 and negative connectivity in the contralesional M1 correlated with motor performance and the ipsilesional cortical excitability.

4.1. Motor Cortical Excitability and Activity

In the present study, the assessments of neural activity and motor cortical excitability provided direct multimodal evidence that the less ipsilesional M1 recruitment signifies poor motor performance in the acute phase after stroke. Our data identified that the reduced cortical excitability and activity in the ipsilesional M1 were associated with the degree of motor impairment in stroke patients. The results were in line with a recent meta-analysis indicating activity in ipsilesional M1 characterizing favorable motor recovery after stroke [27]. Furthermore, we also found a positive correlation between the brain activity in the ipsilesional sensorimotor area and RMT of ipsilesional M1, which implied that the more recruited ipsilesional sensorimotor motor cortex, the less decreased motor cortical excitability of the ipsilesional hemisphere. The parietal lobe of the cortex is responsible for somatosensation and has been shown to be involved in sensorimotor integration for hand coordination [28]. This finding may, at least in part, be related to the fact that passive motor task influenced the neural activity in the somatosensory cortex. However, due to the variance in RMT that was considerably larger in the ipsilesional motor cortex, this observation needs further investigation.

4.2. Effective Connectivity in Motor Areas

Importantly, we combined fMRI and neurophysiological measures to investigate the effective connectivity in key motor regions in relation to cortical excitability and motor performance in the acute stage of stroke. Our findings further confirmed the former studies enrolling subacute or chronic stroke patients and demonstrated that a reduction of pathological transcallosal influences (originating from contralesional M1 to ipsilesional M1) may underlie improved motor performance [3, 26].

In the concept of interhemispheric inhibition model theory, both hemispheric M1 exhibit a mutual inhibitory influence on each other [29, 30], which is reflected in the present study data by showing a mutual suppression of M1-M1 connections. We demonstrated that patients with more enhanced inhibitory connectivity to ipsilesional M1 or less intra- and interhemispheric inhibition to contralesional M1 showed stronger motor impairment and lower ipsilesional cortical excitability. In addition, our data demonstrated that the neural coupling between contralesional SMA and contralesional M1, as well as ipsilesional M1, significantly correlated with the degree of motor impairment. More recent neuroimaging and electrophysiological data provided the evidence of the role of SMA showing that a large proportion of SMA neurons exclusively respond to contralateral hand movements only [12, 31]. These findings are compatible with our data suggesting that the SMA might represent a key region promoting or inhibiting the influence on the cortical motor network during movements of the affected or unaffected hand.

4.3. Limitation

Several limitations were noteworthy. Firstly, given the cross-sectional design of the study, we can only provide preliminary evidence of correlations across different modalities, without identifying their causative interactions. Future longitudinal studies can be performed to investigate whether the factors identified in the present study (combination of RMT, BOLD activation, and effective connectivity) have a predictive value in the acute poststroke phase for long-term recovery. Secondly, this study did not involve the healthy participants as controls.

5. Conclusions

By combining fMRI and electrophysiological measures, this study revealed the aberrances of cortex excitability and connectivity in the acute stage of stroke. The findings contributed to the understanding of the pathophysiology of motor impairment in the early phase of stroke. The inhibitory connectivity to ipsilesional M1 and a relative disinhibition of contralesional M1 may constitute an important pathophysiological aspect of motor disability in acute stroke. Moreover, the present findings might have important implication for early rehabilitative therapy in stroke patients. Aberrances of the cortical excitability and connectivity in the ipsilesional M1 might be considered the therapeutic targets in a few of novel therapeutic strategies (e.g., noninvasive brain stimulation) for the early stage of stroke [3234].

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Additional Points

Highlights. The neural mechanisms of motor deficit in acute stroke remain unclear. This multimodal study demonstrated the relationships of cortical excitability, effective connectivity, and motor impairment following stroke. Our findings provided important theoretical implication for early rehabilitative therapy in stroke patients.

Disclosure

Juan Du and Jianping Hu share the first authorship.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Juan Du and Jianping Hu contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 81530038, 81501193, 81220108008, 81701299, and 81771424), Natural Science Foundation of Jiangsu Province (no. BK20141373), the special scientific research fund of Public Welfare Profession of National Health and Family Planning Commission of the People’s Republic of China (no. 201402019), the independent research project in State Key Laboratory of Analytical Chemistry for Life Science (no. 5431ZZXM1716).