Dynamic Structural and Functional Reorganizations Following Motor Stroke.

BACKGROUND The combined effects of bilateral corticospinal tract (CST) reorganization and interhemispheric functional connectivity (FC) reorganization on motor recovery of upper and lower limbs after stroke remain unknown. MATERIAL AND METHODS A total of 34 patients underwent magnetic resonance imaging (MRI) examination at weeks 1, 4, and 12 after stroke, with a control group of 34 healthy subjects receiving 1 MRI examination. Interhemispheric FC in the somatomotor network (SMN) was calculated using the resting-state functional MRI (rs-fMRI). Fractional anisotropy (FA) of bilateral CST was recorded as a measure of reorganization obtained from diffusion tensor imaging (DTI). After intergroup comparisons, multiple linear regression analysis was used to explore the effects of altered FA and interhemispheric FC on motor recovery. RESULTS Interhemispheric FC restoration mostly occurred within 4 weeks after stroke, and FA in ipsilesional remained CST consistently elevated within 12 weeks. Multivariate linear regression analysis showed that the increase in both interhemispheric FC and ipsilesional CST-FA were significantly correlated with greater motor recovery from week 1 to week 4 following stroke. Moreover, only increased FA of ipsilesional CST was significantly correlated with greater motor recovery during weeks 4 to 12 after stroke compared to interhemispheric FC. CONCLUSIONS Our results show dynamic structural and functional reorganizations following motor stroke, and structure reorganization may be more related to motor recovery at the late subacute phase. These results may play a role in guiding neurological rehabilitation.

Background

Although the mortality rate has declined in the past 30 years, stroke remains the most important cause of adult disability [1]. Most stroke patients generally experience a certain degree of motor recovery within the first several months following stroke, even without treatment. Several studies have shown that reorganization of structure and functional connection in the motor system significantly contributes to the restoration of motor functions after stroke [2-4].

The corticospinal tract (CST) is the major descending pathway and plays a key role in motor control. In some research, changes in microstructural integrity were found using fractional anisotropy (FA) extracted from diffusion tensor imaging (DTI) [5-7]. The elevated FA values were assumed to reflect extensive white matter reorganization and improvement in the structural integrity of nerve fibers during stroke recovery [8-10]. Some previous studies have found that FA in perilesional or contralesional CST is positively correlated with motor performance in murine models or in patients at the chronic stage after stroke [8,9,11,12].

Numerous studies have investigated the altered interhemispheric functional connectivity (FC) in the somatomotor network (SMN) using resting-state functional MRI (rs-fMRI) after stroke, which has been shown to be related to motor performance [13-17]. In fact, almost all patients exhibit a decrease in interhemispheric FC at the initial stage after stroke. Subsequently, increased interhemispheric FC shows a positive association with motor recovery in the subacute or chronic phase of stroke [18-20].

Although reorganization in CST or interhemispheric FC is highly relevant to motor recovery after stroke, their combined effects on motor recovery over time remain unknown. In the current study, we combined rs-fMRI and DTI techniques to explore longitudinal changes in interhemispheric FC and bilateral CST across a continuum from acute to subacute stages of stroke and to examine the combined effect of these 2 processes on motor recovery over time.

Material and Methods

Participants

This study was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University, and all participants signed written informed consent before their enrolment. Stroke patients were prospectively recruited at the Department of Neurology of Guangxi Medical University between January 2019 and December 2019. Patients were included in this study if they: (1) were between 18 and 80 years old; (2) had their first-ever unilateral ischemic stroke confirmed by MRI scan, and the lesion locations were restricted to the internal capsule and neighboring regions; (3) exhibited clinically detectable motor deficits based on neurological examination; (4) had sufficient cognition to participate in the study and take the Mini-Mental State Examination (MMSE) ≥21. Patients were excluded if they had: (1) previous motor impairments; (2) a history of neurological or psychiatric disease; (3) MRI contraindications, such as claustrophobia, among others; and (4) heart and renal failure, malignancies, as well as other intercurrent severe illnesses that preclude active participation in research.

Age- and sex-matched healthy controls were also recruited if they: (1) exhibited a normal neurological examination; (2) denied having suffered from a neurological or mental illness; and (3) did not have MRI contraindications.

Clinical Evaluation

We recorded participants’ baseline characteristics, including age and sex, lesions, and vascular risk factors, as well as parameters from neurological examination using the National Institutes of Health Stroke Scale (NIHSS) and Mini-Mental State Examination (MMSE) at admission. Patients were discharged from the neurology ward after a period of stabilization. Subsequently, they underwent conventional stroke rehabilitation treatment at the outpatient department or rehabilitation wards. Each patient’s synergistic motor control of the paretic extremities was evaluated using the upper- (scale 0–66) and lower- (scale 0–34) extremity section of the Fugl-Meyer Assessment (FMA), 2 h before each MRI scan. The total possible FMA score was 0 to 100. In addition, the degree of disability in performing daily living activities was assessed based on the Barthel Index (BI), which uses a scale ranging from 0 to 100.

MRI Data Acquisition

We obtained rs-fMRI and DTI data from all enrolled patients on weeks 1, 4, and 12 after stroke onset, hereafter denoted W1, W4, and W12, respectively, with a control group of 34 healthy subjects receiving 1 MRI examination. Images were obtained using Siemens Prisma 3.0 Tesla MR scanners (Siemens, Erlangen, Germany) in the Radiology Department at the First Affiliated Hospital of Guangxi Medical University.

A gradient echo-planar imaging (EPI) sequence was used for rs-fMRI scanning (TR/TE=2000/30 ms, flip angle = 90° field of view (FOV) = 240×240 mm, number of slices=38, slice thickness=3 mm, gap=1 mm, and number of time points=240, total scan time=6 min 20 s). These 38 slices covered the entire cerebral cortex. High-resolution structural scans were obtained with T1-weighted brain volume sequence (TR/TE =8/3.5 ms, flip angle=11°, FOV=256×256 mm, slice thickness=1 mm, no gap, voxel resolution=1.0×1.0×1.0 mm). A single-shot echo-planar imaging (EPI) sequence which consisted of 30 directions (TR/TE=9600/92 ms, FOV =256×256 mm, flip angle=90°, slice thickness=3 mm, no gap, and b value: 0 and 1000 s/mm2) was used for DTI scanning. During the scanning, all subjects were instructed to close their eyes but not sleep. In addition, they were asked not to think about anything.

Lesion Analysis

MRIcron software () was used to outline the lesion overlap images for all stroke patients. More details about the processing of lesion overlap maps have been reported in a previous publication [21]. Maps of the lesion overlap are shown in Figure 1.

Figure 1

Lesion overlapped across the 34 stroke patients. Color bar indicates the number of subjects having lesions in each voxel.

DTI Data Pre-processing

DTI data were preprocessed using the FMRIB’s Software Library (FSL) version 5.0.9 (). Briefly, pre-processing included eddy current correction, head motion correction, and brain extraction. After pre-processing, diffusion indices were calculated for each voxel, and diffusion tensors were reconstructed using the linear least-square fitting algorithm [22]. Then whole-brain fiber tractography was performed using Diffusion Toolkit (). The CST tracts were tracked using the regions of interest (ROIs) approach implemented in the Diffusion Toolkit (). Specifically, 3 ROIs, selected based on previous studies [19,23], were placed in the precentral gyrus, the posterior limb of the internal capsule (PLIC), and the cerebral peduncle on the individual axial FA map. Exact anatomic positions were defined according to a neuroanatomical atlas [24]. Fibers passing through all ROIs were defined as the CST. We only extracted FA values of CST from the PLIC using the Trackvis software, as the CST fibers in PLIC have a relatively coherent arrangement. Further details on image reconstruction can be found in a previous publication [19].

rs-fMRI data pre-processing

Pre-processing of rs-fMRI data was performed using MATLAB version 9.1 () and Statistical Parametric Mapping software (SPM8, ) [25]. Briefly, we removed the first 5 functional volumes and initial pre-processing for the remaining volumes included compensation for slice-dependent time shifts and correction of head motion. Participants with a maximum displacement of less than 2 mm and a maximum rotation of less than 2.0° were selected for this study. The pre-processing steps included: resampling to the standard space of the Montreal Neurological Institute (MNI) (voxel size=3×3×3 mm); selection of a band pass frequency filter (0.01–0.08 Hz) to eliminate unwanted frequencies; and application of the 8mm gaussian kernel for spatial smoothing. Then, linear regression analysis was performed to remove several nuisance variables from the smoothened and filtered data including head motion parameters, as well as signals of the whole brain, white matter, and cerebrospinal fluids. The residual BOLD time course was also recorded for FC analysis. The ROIs were selected in the somatomotor network (SMN) along the left and right primary motor cortex (M1) as previously defined [26]. (left: MNI coordinates, −11, −26 and 78; radius=6 mm; right: MNI coordinates, 11, −26, and 78; radius=6mm). The average time series for the left and right seed areas were extracted and their linear correlation coefficients (r) were calculated. Fisher r-z transformation was then performed to convert the r to a normally distributed z score. The z scores were subsequently used as the FC of each time point for statistical analysis.

Statistical Analysis

All statistical analyses were performed using SPSS version 23.0. Differences in FC, FA, BI, and FMA indicators at week 1 between the stroke group and healthy controls were analyzed by two-sample independent t tests. Data of each parameter at 3 timepoints were compared using a repeated-measures ANOVA in the stroke group. Linear correlation analysis was also implemented to assess the relationship between percentage changes in FA and those in FC in weeks 1–4 and 4–12. A multivariate linear regression method was then applied to explore the effects of percentage changes in FA and interhemispheric FC on percentage changes in FMA in weeks 1–4 and 4–12. With reference to our previous study [27], the percentage change was calculated as: {[(W4-W1)/W1]×100%} in weeks 1–4, and {[(W12-W4)/W1]×100%} in weeks 4–12. Values of P<0.05 were considered statistically significant. Because the initial motor impairment (initial FMA score) is a comprehensive parameter reflecting characteristics of the lesions, the initial FMA score was included in multivariable models based on the literature [12]. The multicollinearity was tested using SPSS when conducting multiple linear regression analyses. If multicollinearity was found among covariates (Variance Inflation Factor, VIF >5), variables with VIF >5 were excluded.

Results

Clinical Data

The process for patient selection and the resulting data are summarized in Figure 2, whereas the general clinical baseline characteristics across patient groups and healthy controls are outlined in Table 1. The lesion was on the left side in 18 cases, and the right side in 16. The sizes of lesions ranged from 0.8 to 54.3 cm3, median 9.4 cm3. Overall, most patients exhibited mild to moderate stroke, with a median NIHSS score of 7 (Inter Quartile Range 4–13). All stroke patients had received conventional rehabilitation with or without acupuncture for at least 4 weeks, with 15 patients in rehabilitation wards and 19 in the outpatient department. Twenty patients completed conventional rehabilitation for 8 weeks. FMA and BI scores in stroke patients further showed a progressive increase across the 3 time points (Figure 3C, Table 2).

Figure 2

Patient selection.

Table 1

Baseline characteristics of patients and controls.

Characteristics	Patients, n=34	Controls, n=34	P value
Age (years)	63.4 [52.1–73.8]	61.2 [49.3–70.5]	0.532*
Female N (%)	14 (41.2%)	16 (47.1%)	0.808**
Vascular risk factors N (%)
Hypertension	23 (67.6%)	18 (52.9%)	0.215**
Diabetes	19 (55.9%)	13 (38.2%)	0.145**
Hypercholesterolemia	20 (58.8%)	15 (44.1%)	0.225**
NIHSS admission	7 [4–13]	–	–
MMSE	≥21	–	–
Lesion size (cm3)	9.4 [3.2–35.1]	–	–
Lesion side (left) N (%)	18 (52.9%)	–	–

NIHSS – National Institutes of Health Stroke Scale; MMSE – Mini-Mental State Examination. Continuous variables presented as median [P25–P75].

Assessed by Mann-Whitney U test;

assessed by chi-square test.

Figure 3

Changes over time. FA – fractional anisotropy, FMA – Fugl-Meyer Assessment, FC – functional connectivity. FA values increased longitudinally in ipsilesional CST (A). Interhemispheric FC restoration happened mostly within 4 weeks after stroke (B). The Fugl-Meyer Assessment (FMA) scores of stroke patients improved gradually (C).

Table 2

Clinical and imaging indicators of subjects at different time points.

	Controls (n=34)	Patients W1 (n=34)	Patients W4 (n=34)	Patients W12 (n=34)
FMA (0–100)	–	51.315±9.209	73.355±8.854**	89.703±7.870**,***
Barthel index (0–100)	–	54.732±12.937	77.708±8.634**	86.368±14.409**,***
Interhemispheric FC	1.01±0.043	0.587±0.119*	0.929±0.077*,**	0.969±0.093*,**,***
FA in ipsilesional CST	0.547±0.058	0.448±0.029*	0.492±0.033*,**	0.544±0.045**,***
FA in contralesional CST	0.547±0.058	0.538±0.055	0.544±0.048	0.548±0.040

FA – fractional anisotropy; FMA – Fugl-Meyer Assessment; FC – functional connectivity. The FA of healthy controls were derived from averages of bilateral CST measurements.

Compared with the control group, P<0.05 obtained by two-sample independent t tests;

compared with week 1, P<0.05 obtained by RM-ANOVA;

compared with week 4, P<0.05 obtained by RM-ANOVA.

CST Reorganization

FA values in the ipsilesional remaining CST showed a decline during the first week and then increased longitudinally over the next 12 weeks after stroke (Figures 3A and 4, Table 2). Results further revealed that mean percentage changes of FA values in the ipsilesional remaining CST across 1–4 weeks (22.68±2.50%) and 4–12 weeks (23.01±1.90%) were not significantly different (P=0.524). In addition, contralesional CST did not show a significant change from week 1 to week 12 after stroke (Table 2).

Figure 4

FA pseudocolor pictures of one stroke patient at W1, W4, W12, levels of the internal capsule. Arrows indicate the increased FA in the perilesional CST.

Reorganization of the Interhemispheric FC

A significantly lower interhemispheric FC was found in stroke patients relative to healthy controls in the first week, after which a longitudinal increase through week 12 was evident. Interhemispheric FC restoration occurred mostly within 4 weeks after stroke, with only a minor recovery in the following 8 weeks (Figure 3B, Table 2). The difference in mean percentage changes of interhemispheric FC between weeks 1–4 (58.26±10.56%) and weeks 4–12 (6.81±3.98%) was statistically significant (P<0.001).

Correlation Analysis

We found a positive correlation between percentage changes of FA in ipsilesional remaining CST and those of interhemispheric FC in patients during weeks 1 to 4 after stroke (r=0.695, P<0.001, Figure 5A). However, no correlation was observed between them from weeks 4 to 12 (r=0.306, P=0.079, Figure 5B).

Figure 5

Significantly positive correlation between percentage changes of FA in ipsilesional remaining CST and those of interhemispheric FC during weeks 1 to 4 after stroke (A). No correlation between them from weeks 4 to 12 (B). FA – fractional anisotropy, FC – functional connectivity, CST – corticospinal tract.

Multivariate Linear Regression

The initial FMA score was removed from the multivariable model due to multicollinearity. Multivariate linear regression analysis showed that the increase in both interhemispheric FC and ipsilesional CST-FA were significantly correlated with the greater change of FMA between weeks 1 and 4 following stroke. Moreover, only increased FA of ipsilesional CST was significantly correlated with the greater change of FMA at weeks 4–12 after stroke compared to interhemispheric FC (Table 3).

Table 3

Multiple linear regression of factors influencing motor recovery.

Model 1		Model 2
R2	0.701	R2	0.309
Total (P value)	<0.001	Total (P value)	<0.001
Percentage changes of FA(W1–W4)Beta (P value)	0.304 (0.011)	Percentage changes of FA(W4–W12)Beta (P value)	0.557 (0.001)
Percentage changes of interhemispheric FC(W1–W4)Beta (P value)	0.647 (<0.001)	Percentage changes of interhemispheric FC(W4–W12)Beta (P value)	0.086 (0.571)

Beta – standardized regression coefficients. Model 1 with percentage changes of FMA(W1–W4) as dependent variable. Model 2 with percentage changes of FMA(W4–W12) as dependent variable.

Discussion

The key findings of the current study included: (1) Ipsilesional remaining CST and interhemispheric FC underwent reorganization within 12 weeks after stroke and showed different time courses; and (2) The combined effects of these 2 processes on motor recovery varied with time. These findings can provide a better understanding of motor recovery mechanisms. Taken together, the results of this study may play a role in guiding neurological rehabilitation.

Previous studies have shown that stroke is accompanied by microstructural changes in the white matter, even in undamaged areas [8,28,29]. In ipsilesional white matter areas, DTI studies using stroke animal models have reported significantly increased FA values in the perilesional white matter, with postmortem histological sections indicating the presence of high-density axons and myelin in this area [8,30,31]. The present study showed an increase in longitudinal FA in ipsilesional remaining CST regions, consistent with the findings of Lu et al [32], who placed the ROIs within the CST but outside of the lesion to extract the FA values. These findings indicate ipsilesional remaining CST undergoes remodeling after stroke. However, it is not clear whether contralesional corticospinal pathways undergo reorganization after stroke. A previous study using severely impaired long-term stroke patients showed higher FA values in contralesional CST compared to normal controls, and this was associated with favorable motor outcomes [9]. However, the results of the current study indicated that FA values in the contralesional CST were relatively stable from weeks 1 to 12 after stroke. In the present study, our patient population mainly exhibited mild to moderate stroke, unlike in previous studies that showed contralesional CST reorganization based on subjects that mostly had a severe stroke. This variation in stroke severity possibly explains the variation in the range of CST reorganization.

Several longitudinal studies have demonstrated that interhemispheric FC in the SMN first decreases in the early stages after stroke, then increases in the following weeks or months [15,18,33]. Some studies that selected subjects similar to those in the current research have explored the dynamic evolution of interhemispheric FC after stroke [15,33,34]. For example, Golestani et al [15] used resting-state fMRI scans on patients across 3 time points (<24 h, 7 days, and 90 days) after stroke, and found that interhemispheric FC reestablished normality after 1 week in some patients who had a mild stroke. In another study, mild to moderate stroke patients who underwent rs-fMRI at week 5 and week 26 after stroke exhibited no change in interhemispheric FC, although their motor performance scores improved during this period [34]. In the current study, we selected participants with mild to moderate stroke and found that their interhemispheric FC improved longitudinally. Notably, the interhemispheric FC restoration occurred mostly within 4 weeks after stroke. These findings suggest that the first few weeks after stroke are the best time window for interhemispheric FC reorganization for patients with mild to moderate stroke.

A significant correlation between increased interhemispheric FC and motor recovery has been demonstrated by numerous studies [17,20,32,33,35,36]. Among those, only a few studies have reported a strong correlation between elevated FA values in CST and improved motor function within 10 weeks [35] or 3 months [32]. As post-stroke motor recovery may reflect a complex interaction of structural and functional reorganization in brain motor networks, we used multiple linear regression to explore the combined effects of structural and functional reorganization on motor recovery. Our results showed that the reorganization in both interhemispheric FC and ipsilesional CST were significantly correlated with greater motor recovery during weeks 1 to 4 following the stroke. Moreover, ipsilesional CST reorganization was more significantly correlated with greater motor recovery at weeks 4–12 after stroke compared to interhemispheric FC. These results may be related to the different time courses of structural and functional reorganization in the brain motor network. Recently, Lin et al [12] reported this difference in the time course of reorganization beyond 3 months after stroke. Their results showed that CST reorganization may occur from 3 to 12 months after stroke and is related to improved motor function, while interhemispheric FC does not change during this period. The present study focused on the longitudinal changes in bilateral CST and interhemispheric FC within 12 weeks after stroke, revealing that interhemispheric FC restoration mainly occurred within 4 weeks after stroke, with only a minor recovery observed during the following 8 weeks. However, increased FA values in ipsilesional remaining CST were comparable across these 2 stages. Our study complements those of Lin et al

Several studies have revealed a close relationship between changes in CST and upstream interhemispheric FC in SMN [19,32,35,37-39]. For example, CST impairment was found to be positively correlated to decreased somatomotor interhemispheric FC within 4 weeks after subcortical stroke [37]. A study by Van et al [35] evaluated the links between brain reorganization and neurological function recovery from days 3 to 70 after stroke in mice. They found a significant positive correlation between increased interhemispheric FC in SMN and elevated FA in the ipsilesional CST. In the current study, our results revealed a positive correlation between CST reorganization and increased interhemispheric FC during weeks 1 to 4 following stroke. However, there was no correlation between these 2 processes at weeks 4–12. Therefore, although CST reorganization is closely related to that of interhemispheric FC, it is likely that changes unrelated to the reorganization of FC also occur.

Our study had some limitations. Firstly, for the study sample was a narrow range of mild to moderate stroke patients, so the generalizability of results is probably limited. Secondly, participants’ lesions were located in subcortical regions, implying that our results cannot be generalized to individuals with a stroke occurring in the cortex. Finally, many independent variables might affect the long-term motor outcomes, such as lesion volume, lesion location, different rehabilitation programs, and the duration of rehabilitation, which could significantly affect the outcome. Future studies in a more general stroke population are needed to validate our findings.

Conclusions

Our results show the dynamic structural and functional reorganizations following motor stroke, and suggest that structure reorganization may be more related to motor recovery at the late subacute phase. These results may play a role in guiding neurological rehabilitation.