ORIGINAL RESEARCH article

Front. Neurol., 13 November 2025

Sec. Artificial Intelligence in Neurology

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1666991

Assessment of functional decline in stroke patients using 3D deep learning and dynamic functional connectivity based on resting-state fMRI

  • 1. The First Hospital of Xinjiang Production and Construction Group, Aksu, Xinjiang, China

  • 2. Department of Biomedical Engineering, Children's Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health, Zhejiang University, Hangzhou, China

  • 3. Department of Rehabilitation, Huashan Hospital, Fudan University, Shanghai, China

Abstract

Introduction:

This study aimed to develop an automated approach for assessing upper limb (UL) motor impairment severity in stroke patients using a deep learning framework applied to resting-state functional magnetic resonance imaging (rs-fMRI).

Methods:

Dynamic functional connectivity (dFC) was computed with the ipsilesional primary motor cortex (M1) as a seed and extracted from rs-fMRI data of 69 stroke patients. These dFC features were used to train a three-dimensional convolutional neural network (3D-CNN) for automatic classification of UL motor impairment severity. Patients were divided into two groups according to UL Fugl-Meyer Assessment (UL-FMA) scores: mild-to-moderate impairment (UL-FMA > 20; n = 29, maximum = 66) and severe impairment (0 ≤ UL-FMA ≤ 20; n = 40). UL-FMA scores served as labels for supervised learning.

Results:

The model achieved a balanced accuracy of 99.8% ± 0.2%, with a specificity of 99.9% ± 0.2% and a sensitivity of 99.7% ± 0.3%. Several brain regions—including the angular gyrus, medial orbitofrontal cortex, dorsolateral superior frontal gyrus, superior parietal lobule, supplementary motor area, thalamus, cerebellum, and middle temporal gyrus—were linked to UL motor impairment severity.

Discussion:

These findings demonstrate that a 3D deep learning framework based on dFC features from rs-fMRI enables highly accurate and objective classification of UL motor impairment in stroke patients. This approach may provide a valuable alternative to manual UL-FMA scoring, particularly in clinical settings with limited access to experienced evaluators.

1 Introduction

Stroke is a leading cause of long-term disability, with upper limb (UL) motor dysfunction among the most prevalent and debilitating outcomes (1). Such impairments substantially reduce quality of life (2). Although the mechanisms linking stroke lesions to persistent motor deficits remain incompletely understood, neuroimaging studies consistently implicate both structural and functional brain alterations (3–5). Given the heterogeneity of stroke, treatment strategies vary widely and depend critically on the severity of impairment. Accurate assessment of UL motor dysfunction is therefore essential for guiding clinical decision-making, tailoring rehabilitation, and setting realistic recovery expectations for patients and caregivers (6, 7).

The Fugl-Meyer Assessment (FMA) is widely regarded as a reliable and validated tool for evaluating post-stroke motor function (8). It evaluates five domains—motor function, sensation, balance, joint range of motion, and joint pain—with the motor domain most commonly applied to quantify impairment and recovery (9). The UL-FMA specifically examines motor function in the shoulder-arm, wrist, hand, and coordination/speed subsections, progressing from proximal to distal and from synergistic to isolated voluntary movements (9). Each of the 33 items is scored from 0 (absent) to 2 (normal), yielding a maximum of 66 points that reflects full functional capacity (10). Despite its clinical value, the UL-FMA has several limitations: it requires trained therapists, depends on patient cooperation, and is vulnerable to inter-rater variability (9–11). It may also exhibit ceiling effects in higher-functioning patients and is challenging to administer in resource-limited settings or with uncooperative individuals (12, 13). These limitations highlight the need for an objective, automated method of assessing UL motor impairment severity.

Resting-state functional magnetic resonance imaging (rs-fMRI) offers a complementary neuroimaging approach to evaluate post-stroke motor dysfunction, particularly through analyses of functional connectivity (FC) within motor-related regions. Prior work has demonstrated post-stroke FC reorganization involving the ipsilesional primary M1 and regions such as the thalamus, SMA, middle frontal gyrus, and cerebellum (3, 14–17). Interhemispheric M1 connectivity has also been positively correlated with motor performance during the subacute phase of stroke (12, 18), suggesting that ipsilesional M1 FC may represent a biomarker of motor impairment. Traditional rs-fMRI analyses, however, often assume that FC remains stable during the scanning. Emerging evidence indicates that brain connectivity is inherently dynamic (19–21). Dynamic FC (DFC) captures temporal fluctuations in FC, providing a richer characterization of brain activity than static approaches (21). For example, dFC between the ipsilesional M1 and contralesional precentral gyrus has been negatively correlated with FMA scores in stroke patients (22). Using sliding time-window methods, dFC generates time-resolved connectivity matrices that can be clustered to identify distinct connectivity states and their temporal properties (23).

Deep learning has become a powerful tool in medical image analysis, offering automated feature extraction and classification capabilities that surpass traditional approaches (24, 25). In stroke research, deep learning holds promise for minimizing reliance on expert evaluation in the classification of motor impairment severity (26, 27). Unlike conventional algorithms that require manual feature engineering and operate on one-dimensional data, deep learning models—particularly three-dimensional convolutional neural networks (3D-CNNs)—can directly process volumetric neuroimaging data (28, 29). 3D-CNNs have shown strong performance in various medical imaging tasks, excelling in classification and pattern recognition within supervised frameworks (24, 26, 30). Although their application to brain imaging remains relatively limited due to computational demands, dFC-derived features from the ipsilesional M1 can reduce data dimensionality while retaining spatiotemporal patterns relevant to stroke severity.

In this study, we propose a novel framework for automated dichotomous classification of UL motor impairment (mild/moderate vs. severe) in stroke patients. We hypothesize that dFC features derived from rs-fMRI, particularly those involving the ipsilesional M1, can serve as effective biomarkers of motor impairment severity. By leveraging a 3D-CNN architecture, our goal is to enable accurate, objective classification and provide a viable alternative to conventional clinical assessment tools such as the UL-FMA.

2 Methods

2.1 Subjects

A total of 69 patients in the chronic phase of subcortical stroke were recruited from Huashan Hospital and underwent neuroimaging at the Shanghai Key Laboratory of Magnetic Resonance. Inclusion criteria were: (1) first-ever subcortical stroke; (2) duration of illness > 3 months; (3) Montreal Cognitive Assessment (MoCA) score > 23, indicating preserved global cognition; and (4) right-handedness. Exclusion criteria were: (1) contraindications to MRI; (2) bilateral stroke lesions; and (3) comorbid neurological or psychiatric disorders unrelated to stroke. Cognitive function was assessed using the MoCA, administered by an experienced therapist at admission. UL motor function was evaluated using the UL-FMA scale. The study was conducted in accordance with the Declaration of Helsinki, approved by the Review Board of Ethics Committee of Huashan Hospital, and registered at the Chinese Clinical Trial Registry (ChiCTR-TRC-08003005). Written informed consent was obtained from all participants.

To investigate resting-state functional connectivity alterations associated with UL motor impairment severity, patients were categorized into two groups based on their UL-FMA scores: those scoring > 20 were classified as having mild-to-moderate impairment (Mild/Moderate Stroke Patients, MSP; n = 29), and those with scores ≤ 20 as having severe impairment (Severe Stroke Patients, SSP; n = 40) (9, 13, 31). Demographic comparisons between groups were performed using SPSS version 23.0 for Windows.

2.2 MRI data acquisition

All imaging was conducted on a 3.0-T Siemens MRI scanner (Erlangen, Germany) at the Shanghai Key Laboratory of Magnetic Resonance. High-resolution T1-weighted structural images were acquired using a magnetization-prepared rapid gradient echo (MPRAGE) sequence with 192 sagittal slices and the following parameters: repetition time (TR) = 1900 ms; echo time (TE) = 3.42 ms; inversion time (TI) = 900 ms; flip angle = 9°; field of view (FOV) = 240 × 240 mm2; matrix = 256 × 256; slice thickness/gap = 1/0.5 mm. T2-weighted images for lesion localization were obtained with a turbo spin-echo sequence (30 axial slices; TR = 6,000 ms; TE = 93 ms; FOV = 220 × 220 mm2; flip angle = 120°; matrix = 320 × 320; slice thickness/gap = 5/0 mm). Rs-fMRI data were collected using an echo-planar imaging (EPI) sequence (30 axial slices; TR = 2000 ms; TE = 30 ms; flip angle = 90°; FOV = 220 × 220 mm2; matrix = 64 × 64; slice thickness/gap = 4/0.8 mm; 240 volumes). The total acquisition time was 8 min 6 s, including a 6-s dummy scan to ensure magnetization equilibrium. During rs-fMRI scanning, participants were instructed to remain still, keep their eyes closed, stay awake, and avoid focused thought.

2.3 Lesion analysis

Lesion volumes were manually delineated on T1- and T2-weighted images by two experienced radiologists using MRIcron software1. Lesions were outlined slice by slice for each patient, with T1-weighted images used for primary outlining and T2-weighted images used for validation. To standardize lesion laterality, all images were mirrored so that lesions were represented on the left hemisphere (Figure 1). The resulting lesion masks were used for visualization, lesion volume calculation, and exclusion of intra-lesion voxels during dFC analysis. Notably, lesion volume was not significantly correlated with UL motor function, as measured by UL-FMA scores (Pearson’s r = −0.208, p = 0.087), suggesting that stroke size alone did not account for impairment severity in this cohort.

Figure 1

2.4 fMRI data preprocessing

Image preprocessing and analysis were performed using Statistical Parametric Mapping (SPM122) and the Data Processing Assistant for RS-fMRI (DPABI3). For patients with right-hemisphere lesions, images were flipped along the midsagittal plane to standardize orientation, designating the left hemisphere as ipsilesional and the right as contralesional for all participants. Preprocessing steps were as follows: (1) The first 10 volumes of each functional scan were discarded to allow for magnetization equilibrium and participant adaptation. The remaining 230 volumes were corrected for slice timing, and head motion. No subject was excluded, as none exceeded the motion threshold of 2.0 mm translation or 2.0° rotation; (2) Linear trends, mean white matter and cerebrospinal fluid (CSF) signals, and 24 motion parameters (translations, rotations, their temporal derivatives, and quadratic terms) were regressed out of each voxel time series. Global signal regression was not applied to preserve correlation matrix structure, particularly for sliding-window analyses (32, 33); (3) Functional images were spatially normalized to Montreal Neurological Institute (MNI) space using each patient’s lesion mask and the DARTEL (Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra) algorithm (34). Images were resampled to 3 × 3 × 3 mm3 voxels using the DARTEL-derived deformation fields; (4) Time-series data were spatially smoothed with a 6 mm full-width at half-maximum (FWHM) Gaussian kernel, linearly detrended, and band-pass filtered (0.01–0.1 Hz) to reduce physiological and low-frequency noise.

2.5 Dynamic functional connectivity

DFC analysis was conducted using the Temporal Dynamic Analysis (TDA) module within the DPABI toolbox (35). A sliding-window approach was applied to capture temporal fluctuations in connectivity between the left primary M1 and all other brain voxels (36). The left M1 seed was defined at MNI coordinates (−38, −22, 56), consistent with prior studies implicating this region in UL motor function after stroke (3, 12, 16, 37, 38). Each sliding window covered 22 TRs (44 s) and shifted by 1 TR (2 s), enabling fine-grained tracking of temporal dynamics (19). This window length was selected according to the lower limit of the frequency spectrum (<0.5/flower = 50 s; here, flower = 0.01 Hz), ensuring reliable estimation while preserving sensitivity to short-term fluctuations (39). A 22-TR window therefore represents a balance between temporal resolution and estimation stability, consistent with prior literature (16, 19, 39). From the 230 rs-fMRI volumes, each participant contributed 209 overlapping sliding windows. Within each window, Pearson correlation coefficients were calculated between the left M1 time series and all other brain voxels, generating a dynamic series of whole-brain correlation maps (Figure 2). Fisher’s r-to-z transformation was then applied to improve normality of the correlation values (40).

Figure 2

2.6 Feature extraction

The 3D volumetric functional connectivity maps were used as inputs to a 3D-CNN to classify participants into either the MSP or SSP group. Consistent with prior neuroimaging studies, no standard data augmentation techniques were applied, as introducing synthetic data could bias the training process. Rs-fMRI signals are particularly sensitive; even a 3 mm spatial shift can substantially alter their interpretation. For this reason, conventional augmentation strategies are generally discouraged in the neuroimaging community (24, 26, 28).

2.7 Deep learning and 3D-CNN framework

A 3D-CNN–based deep learning framework was developed for participant classification. The model was implemented in MATLAB 2019b with GPU acceleration on an NVIDIA Quadro RTX 5000. Training employed the Adam optimizer (41) with an initial learning rate of 0.001, which was reduced by half every 10 epochs. Each epoch consisted of 32-sample mini-batches over 50 epochs in total. An epsilon value of 0.001 was applied, and cross-entropy loss was used as the cost function. To mitigate overfitting, 10-fold cross-validation was performed to estimate mean classification accuracy at the dFC level. For each fold, 10% of the dFC maps were reserved for testing, while the remaining data were split into training (80%) and validation (10%) sets. The network was based on a modified VGG-Net architecture (42), incorporating batch normalization in convolutional layers and Rectified Linear Unit (ReLU) activation functions. A dropout rate of 0.7 was applied to the fully connected layers to enhance generalization. Weight initialization followed the method proposed by He et al. (43), using a zero-mean Gaussian distribution with a standard deviation scaled to network depth. Hyperparameters—including learning rate, epsilon, dropout rate, batch size, and number of epochs—were optimized based on prior work (24, 26). Specifically, epsilon was tuned in the range [0.1:0.05:1], learning rate across a logarithmic scale [1, 0.1, 0.01, 0.001, 0.0001, 0.00001], and dropout rate between [0.1:0.05:1]. Batch size was optimized using the maximum available GPU memory, while the number of epochs was tuned within [10:1:200]. The complete 3D-CNN architecture is shown in Figure 3 and detailed in Table 1.

Figure 3

Table 1

LayerFeature mapStrideKernelActivation structure
3D Input61 × 73 × 61
Convolution641x1x13x3x3Conv
Convolution641x1x13x3x3Batchnorm + ReLU + Conv
Maxpool2x2x22x2x2
Convolution1281x1x13x3x3Batchnorm + ReLU + Conv
Convolution1281x1x13x3x3Batchnorm + ReLU + Conv
Maxpool2x2x22x2x2
Convolution2561x1x13x3x3Batchnorm + ReLU + Conv
Convolution2561x1x13x3x3Batchnorm + ReLU + Conv
Convolution2561x1x13x3x3Batchnorm + ReLU + Conv
Maxpool2x2x22x2x2
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Maxpool2x2x22x2x2
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Convolution5121x1x13x3x3Batchnorm + ReLU + Conv
Maxpool1x1x12x2x2
Fully Connected4,096Dropout Rate o.7ReLU
Fully Connected4,096Dropout Rate o.7ReLU
Output
Fully Connected2ReLU
Softmax
Classification LayerArgmax

Details of the three-dimensional convolutional neural network (3D-CNN) architecture.

2.8 Dynamic states and clustering

To examine the contribution of common functional connectivity to deep learning classification and to better understand how deep neural networks iteratively optimize classification weights (24), we analyzed dFC between the ipsilesional M1 and other brain regions across individual sliding windows. Because multiple pairwise comparisons can introduce statistical bias, we applied a k-means clustering algorithm to dFC maps from all participants to identify distinct connectivity states at rest and to assess intra-individual variability. Following previous studies (44), L1 distance was chosen as the similarity metric due to its robustness in high-dimensional data. The clustering procedure consisted of two main steps. First, exemplar datasets were selected based on local maxima in functional connectivity variance across windows, consistent with prior work (19, 23). Second, rather than relying on conventional criteria such as the elbow method or Silhouette index—which may not guarantee that each participant is represented across multiple states—we determined the optimal number of clusters (k) by evaluating the participation rate (PR) of subjects across candidate solutions (19, 23). Cluster solutions were examined for k values ranging from 2 to 8. A three-cluster solution (k = 3) was selected because it yielded an adequate distribution of subject participation (state 1: PR = 1/69; state 2: PR = 68/69; state 3: PR = 69/69). States 2 and 3 were retained for further analysis, while state 1 was excluded due to insufficient subject representation, which precluded statistically meaningful comparisons.

2.9 Significance testing

To evaluate the statistical significance of our results, we performed a permutation test on the 3D-CNN classification accuracies, followed by a two-sample t-test to assess between-group differences. For the classification analysis, test data labels were randomly shuffled 1,000 times within each of the 10 cross-validation folds to estimate the probability of achieving an accuracy greater than that obtained with the true labels. Between-group differences in dFC (MSP vs. SSP) were tested using two-sample t-tests, with age, sex, and duration of illness included as covariates. A significance threshold of uncorrected p < 0.05 was adopted.

3 Results

3.1 Demographic characteristics

No significant differences were observed between the MSP and SSP groups in age (MSP: 54.3 ± 11.95 years; SSP: 51.28 ± 10.73 years; p = 0.392, two-sample t-test), sex distribution (MSP: 32 males; SSP: 24 males; p = 0.513, χ2 test), duration of illness (MSP: 8.01 ± 6.5 months; SSP: 9.40 ± 7.10 months; p = 0.410, two-sample t-test), or lesion volume (MSP: 6.7 ± 4.73 mL; SSP: 4.97 ± 4.70 mL; p = 0.513, two-sample t-test). As expected, UL-FMA scores were significantly different, with higher scores in the MSP group compared to the SSP group (31.4 ± 7.84 vs. 8.85 ± 4.49; p < 0.0001). Detailed demographic and clinical characteristics for all participants, including healthy controls (HCs), are presented in Table 2.

Table 2

SSP (n = 40)
Mean ± std
MSP (n = 29)
Mean ± std
SSP vs. MSP
p-value
Age (years)a54.03 ± 11.9551.28 ± 10.730.392
Sex (male: female)8/325/240.513
Lesion volume (ml)a6.7 ± 4.734.97 ± 4.700.137
Duration of illness (month)a8.01 ± 6.529.40 ± 7.100.410
Lesion side (left: right)22/1815/140.490
UL-FMAa8.85 ± 4.4931.4 ± 7.840.0001

Demographic, clinical and structural MRI data of the participants.

aIndependent t-test. UL-FMA, up limb Fugl-Meyer assessment; SSP, severely stroke patients; MSP, mild/moderate stroke patients.

3.2 Classification

The 3D-CNN achieved high classification performance under 10-fold cross-validation. To mitigate potential bias from class imbalance, balanced accuracy was reported alongside other evaluation metrics. The mean performance across folds was as follows: training accuracy = 99.23%, test accuracy = 99.80%, specificity = 99.86%, sensitivity = 99.74%, F-score = 99.78%, and balanced accuracy = 99.80% (Table 3).

Table 3

FoldTrain ACC (%)Test ACC (%)p-valueAUCSpec (%)Sen (%)F-score (%)BAC (%)
199.7299.45<0.0010.999499.499.5199.3499.45
299.9699.65<0.0010.999999.8899.3499.5999.61
399.7299.69<0.0010.999999.7699.6799.6799.71
410099.93<0.001199.8810099.9199.94
5100100<0.0011100100100100
699.9999.86<0.001199.8899.8499.8499.86
799.95100<0.0011100100100100
899.9299.51<0.0010.999999.7699.1799.4299.47
910099.93<0.001110099.8399.9299.92
1099.97100<0.0011100100100100
Mean ± SD99.23 ± 0.1099.80 ± 0.200.9999 ± 0.000299.86 ± 0.1899.74 ± 0.2999.78 ± 0.2499.80 ± 0.21

Classification accuracy using 10-fold cross-validation.

Test ACC, test accuracy; Train ACC, train accuracy; Spec, specificity; Sen, sensitivity; BAC, balanced accuracy.

Statistical significance was confirmed using permutation testing. A threshold of p = 0.001 was applied across all folds, and permutation tests consistently yielded p < 0.001, indicating that the observed classification performance was highly unlikely to arise by chance. Fold-wise permutation-derived p-values are presented in Table 3.

3.3 Clinical significance

Discriminative brain regions contributing to the deep learning framework were identified using two-sample t-tests across dynamic states 2 and 3. Significant regions included the angular gyrus (Angular), MOF, SFG, SPL, SMA, superior temporal pole (STP), thalamus, cerebellum, gyrus rectus, middle temporal gyrus (MTG), and precuneus. Uncorrected t-values indicated group differences in functional connectivity: indices A-E corresponded to state 2 (t = 2.64–3.37), and indices F-M corresponded to state 3 (t = 2.88–3.75) (Table 4 and Figure 4). In Figure 4, red shading denotes regions where functional connectivity was higher in the MSP group compared to the SSP group, while blue shading indicates higher connectivity in SSP relative to MSP.

Table 4

IndexStateGroup comparisonsBrain regions (AAL)Peak MNI coordinatesCluster voxelsPeak T values
XYZ
A2MSP > SSPAngular. R60−6621522.9
B2MSP > SSPMOF. R633−121613.37
C2MSP > SSPSFG. R30606863.01
D2MSP > SSPSFG. R183339582.64
E2MSP > SSPSPL. L−36−66571032.92
F3MSP > SSPSMA. L−3069502.89
G3MSP > SSPSTP. L−42−6−151273.136
H3MSP > SSPThalamus. L3−6151263.15
I3MSP < SSPAngular. L−39−63424713.75
J3MSP < SSPCerebellum3−45−24772.88
K3MSP < SSPRectus. R939−12773.31
L3MSP < SSPMTG. R60−63181293.52
M3MSP < SSPPre. L−3−63361023.25

Statistical analysis of each state.

L, left; R, right; Angular, angular gyrus; MOF, medial orbitofrontal cortex; SFG, dorsolateral superior frontal gyrus; SPL, superior parietal lobule; SMA, supplementary motor area; STP, superior temporal pole; Thalamus, thalamus; Cerebellum, cerebellum; Rectus, gyrus rectus; MTG, middle temporal gyrus; Pre, precuneus; MNI, Montreal Neurological Institute. Two-tailed two-sample t-test, p < 0.05 (|T| > 2.0), cluster size ≥ 50, uncorrected.

Figure 4

4 Discussion

Using the 3D-CNN framework, we achieved high classification performance, with a mean balanced accuracy of 99.80%, specificity of 99.86%, and sensitivity of 99.74%. Accurate assessment of UL motor impairment in stroke is critical for guiding treatment decisions, highlighting the clinical relevance of automated classification methods. To identify features driving classification, we examined state-specific alterations in resting-state dFC of the ipsilesional primary M1. Comparisons between MSP and SSP groups across two connectivity states revealed several discriminative regions—including the angular gyrus, MOF, SFG, SPL, SMA, STP, thalamus, cerebellum, gyrus rectus, and MTG—which were subsequently integrated into the deep learning framework. These findings support our hypothesis that ipsilesional M1-derived dFC can serve as a functional biomarker for classifying UL motor impairment severity using 3D-CNN.

Previous studies have demonstrated the feasibility of using machine learning with multimodal neuroimaging data, including structural MRI and clinical assessments, for automated classification and outcome prediction in stroke (24, 26, 28, 29, 45–48). Several studies successfully predicted UL motor impairment using these approaches (45–48). While structural MRI can distinguish MSP from SSP, it is often limited by inter-individual variability and may fail to capture subtle changes, particularly when lesion volume, shape, or corticospinal tract damage exceeds certain thresholds (47). In contrast, rs-fMRI provides a dynamic and individualized assessment, with shorter acquisition times that enhance clinical efficiency.

Prior research has shown that interventions such as repetitive transcranial magnetic stimulation can improve FMA scores in patients with severe UL impairment by increasing excitability in the ipsilesional M1 (7). Additionally, dFC between the ipsilesional M1 and contralesional precentral and middle frontal gyri has been negatively correlated with FMA scores (22), supporting the potential of M1-derived dFC as a predictor of motor impairment severity. To our knowledge, this study is the first to use dFC maps in combination with a 3D-CNN framework to classify stroke patients into mild/moderate versus severe UL-FMA groups, representing a novel contribution to automated stroke severity assessment.

In our 10-fold cross-validation, the 3D-CNN achieved a mean test accuracy of 99.80%. Although deep neural networks iteratively optimize classification weights (24–26, 28), ranking the contribution of individual dFC features remains challenging without feature selection. To address this, we analyzed dFC differences between MSP and SSP across states using k-means clustering. Statistical analysis revealed significant connectivity differences between the ipsilesional M1 and 13 other regions (p < 0.05, cluster size ≥ 50, uncorrected). Prior work has implicated the contralesional cerebellum (46), SMA (49), and contralesional parietal cortex (50) in motor recovery or stroke severity. Our findings are largely consistent with these reports, while also revealing increased connectivity in the ipsilesional SMA and SPL, which may reflect specific resting-state configurations rather than averaged connectivity patterns (22). Notably, significant connectivity differences were observed in the contralesional angular gyrus, MOF, SFG, gyrus rectus, MTG, as well as the ipsilesional STP, thalamus, angular gyrus, precuneus, and cerebellum, consistent with widespread functional plasticity following stroke (12, 23, 37). These regions likely contributed substantially to classifier performance, although additional regions may also play roles through subtle, weighted patterns learned by the network.

Consistent with previous findings, reduced connectivity between the ipsilesional M1 and contralesional hemisphere is observed after unilateral motor network damage (12). Conversely, increased connectivity within ipsilesional regions may reflect compensatory plasticity. In state 2, MSP patients exhibited higher connectivity in the contralesional angular gyrus, MOF, SFG, and ipsilesional SPL, suggesting better preservation of interhemispheric integration relative to SSP. State 3 revealed more complex patterns: MSP patients showed higher connectivity in the ipsilesional SMA, STP, and thalamus, but lower connectivity in the ipsilesional angular gyrus, precuneus, and contralesional gyrus rectus, MTG, and cerebellum. These results suggest greater network reorganization and motor adaptation in MSP during this state (4), demonstrating a meaningful relationship between dFC patterns and UL-FMA scores (51).

Our study introduces several innovations by integrating rs-fMRI-derived dFC with deep learning for stroke severity classification. However, it has limitations. Despite augmenting sample diversity through the sliding window approach (209 time points per subject), the overall dataset remains modest, limiting the application of end-to-end learning approaches, static functional connectivity-based 3D-CNN classification, and regression models for predicting actual UL-FMA scores. External validation is also lacking. Future studies should expand sample sizes and replicate findings in independent cohorts. Additionally, deep learning models inherently reduce interpretability. Although k-means clustering helped identify salient dFC differences between MSP and SSP, it does not fully reveal how individual features contribute to network decisions. Explainable deep learning techniques should be employed in future work to visualize influential 3D dFC features. Furthermore, our study focused exclusively on dFC; complementary approaches-such as joint time-frequency analysis and dynamic graph theory-may provide additional insights into time-varying connectivity patterns (20, 52). Future research should integrate these methods to build multi-modal representations and improve classification accuracy. Although the 3D-CNN provides an end-to-end approach without additional feature extraction, comparisons with conventional machine learning remain important once feature extraction methods are fully developed.

In conclusion, rs-fMRI-derived dFC combined with 3D-CNN offers a reliable, objective tool for assessing UL motor impairment severity in stroke. This approach may complement clinical scales such as the UL-FMA, particularly in settings with limited access to experienced rehabilitation specialists. In the future, our framework may provide a practical means to explore functional connectivity patterns associated with motor impairment severity under specific resting-state configurations.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Ethics Committee of Huashan Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YG: Validation, Conceptualization, Funding acquisition, Formal analysis, Writing – original draft, Visualization. GX: Formal analysis, Visualization, Methodology, Investigation, Data curation, Software, Writing – review & editing. JP: Software, Data curation, Investigation, Resources, Writing – review & editing. CH: Methodology, Software, Investigation, Validation, Writing – review & editing. SW: Visualization, Software, Validation, Writing – review & editing. MW: Validation, Investigation, Writing – review & editing, Software. HW: Writing – review & editing, Supervision, Investigation, Methodology, Resources, Validation. ZZ: Conceptualization, Investigation, Resources, Writing – original draft, Formal analysis, Supervision, Project administration, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by President’s Fund of Tarim University for Key Guided Medical Special Projects (grant number: TDZKYXZD202502).

Conflict of interest

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.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    DobkinBH. Strategies for stroke rehabilitation. Lancet Neurol. (2004) 3:528–36. doi: 10.1016/S1474-4422(04)00851-8

  • 2.

    Nichols-LarsenDSClarkPZeringueAGreenspanABlantonS. Factors influencing stroke survivors’ quality of life during subacute recovery. Stroke. (2005) 36:1480–4. doi: 10.1161/01.STR.0000170706.13595.4f

  • 3.

    WangLYuCChenHQinWHeYFanFet al. Dynamic functional reorganization of the motor execution network after stroke. Brain. (2010) 133:1224–38. doi: 10.1093/brain/awq043

  • 4.

    ZhaoZWuJFanMYinDTangCGongJet al. Altered intra-and inter-network functional coupling of resting-state networks associated with motor dysfunction in stroke. Hum Brain Mapp. (2018) 39:3388–97. doi: 10.1002/hbm.24183

  • 5.

    GriffisJCMetcalfNVCorbettaMShulmanGL. Structural disconnections explain brain network dysfunction after stroke. Cell Rep. (2019) 28:2527–40. doi: 10.1016/j.celrep.2019.07.100

  • 6.

    KwakkelGVeerbeekJMvan WegenEEWolfSL. Constraint-induced movement therapy after stroke. Lancet Neurol. (2015) 14:224–34. doi: 10.1016/S1474-4422(14)70160-7

  • 7.

    HirakawaYTakedaKTanabeSKoyamaSMotoyaISakuraiHet al. Effect of intensive motor training with repetitive transcranial magnetic stimulation on upper limb motor function in chronic post-stroke patients with severe upper limb motor impairment. Top Stroke Rehabil. (2018) 25:321–5. doi: 10.1080/10749357.2018.1466971

  • 8.

    GladstoneDJDanellsCJBlackSE. The Fugl-Meyer assessment of motor recovery after stroke: a critical review of its measurement properties. Neurorehabil Neural Repair. (2002) 16:232–40. doi: 10.1177/154596802401105171

  • 9.

    MeyerAFJaaskoLLeymanIOlssonSSteglindS. The post stroke hemiplegic patient. I. A method for evaluation of physical performance. Scand J Rehabil Med. (1975) 7:13–31.

  • 10.

    GabelMGilad-BachrachRRenshawESchusterA, eds. Full body gait analysis with Kinect. 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society; (2012): IEEE.

  • 11.

    KurilloGChenABajcsyRHanJJ. Evaluation of upper extremity reachable workspace using Kinect camera. Technol Health Care. (2013) 21:641–56. doi: 10.3233/THC-130764

  • 12.

    ParkC-hChangWHOhnSHKimSTBangOYPascual-LeoneAet al. Longitudinal changes of resting-state functional connectivity during motor recovery after stroke. Stroke. (2011) 42:1357–62. doi: 10.1161/STROKEAHA.110.596155

  • 13.

    Thompson-ButelAGLinGShinerCTMcNultyPA. Comparison of three tools to measure improvements in upper-limb function with poststroke therapy. Neurorehabil Neural Repair. (2015) 29:341–8. doi: 10.1177/1545968314547766

  • 14.

    LuJLiuHZhangMWangDCaoYMaQet al. Focal pontine lesions provide evidence that intrinsic functional connectivity reflects polysynaptic anatomical pathways. J Neurosci. (2011) 31:15065–71. doi: 10.1523/JNEUROSCI.2364-11.2011

  • 15.

    van MeerMPOtteWMvan der MarelKNijboerCHKavelaarsAvan der SprenkelJWBet al. Extent of bilateral neuronal network reorganization and functional recovery in relation to stroke severity. J Neurosci. (2012) 32:4495–507. doi: 10.1523/JNEUROSCI.3662-11.2012

  • 16.

    ZhaoZWangXFanMYinDSunLJiaJet al. Altered effective connectivity of the primary motor cortex in stroke: a resting-state fMRI study with granger causality analysis. PLoS One. (2016) 11:e0166210. doi: 10.1371/journal.pone.0166210

  • 17.

    CassidyJMMarkJICramerSC. Functional connectivity drives stroke recovery: shifting the paradigm from correlation to causation. Brain. (2022) 145:1211–28. doi: 10.1093/brain/awab469

  • 18.

    CarterARAstafievSVLangCEConnorLTRengacharyJStrubeMJet al. Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke. Ann Neurol. (2010) 67:365–75. doi: 10.1002/ana.21905

  • 19.

    AllenEADamarajuEPlisSMErhardtEBEicheleTCalhounVD. Tracking whole-brain connectivity dynamics in the resting state. Cereb Cortex. (2014) 24:663–76. doi: 10.1093/cercor/bhs352

  • 20.

    CalhounVDMillerRPearlsonGAdalıT. The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery. Neuron. (2014) 84:262–74. doi: 10.1016/j.neuron.2014.10.015

  • 21.

    LiRLiaoWYuYChenHGuoXTangYLet al. Differential patterns of dynamic functional connectivity variability of striato–cortical circuitry in children with benign epilepsy with centrotemporal spikes. Hum Brain Mapp. (2018) 39:1207–17. doi: 10.1002/hbm.23910

  • 22.

    ChenJSunDShiYJinWWangYXiQet al. Alterations of static functional connectivity and dynamic functional connectivity in motor execution regions after stroke. Neurosci Lett. (2018) 686:112–21. doi: 10.1016/j.neulet.2018.09.008

  • 23.

    DamarajuEAllenEABelgerAFordJMMcEwenSMathalonDet al. Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia. NeuroImage. (2014) 5:298–308. doi: 10.1016/j.nicl.2014.07.003

  • 24.

    DucNTRyuSQureshiMNIChoiMLeeKHLeeB. 3D-deep learning based automatic diagnosis of Alzheimer’s disease with joint MMSE prediction using resting-state fMRI. Neuroinformatics. (2020) 18:71–86. doi: 10.1007/s12021-019-09419-w

  • 25.

    ShamratFJMAkterSAzamSKarimAGhoshPTasnimZet al. AlzheimerNet: an effective deep learning based proposition for alzheimer’s disease stages classification from functional brain changes in magnetic resonance images. IEEE Access. (2023) 11:16376–95. doi: 10.1109/ACCESS.2023.3244952

  • 26.

    QureshiMNIOhJLeeB. 3D-CNN based discrimination of schizophrenia using resting-state fMRI. Artif Intell Med. (2019) 98:10–7. doi: 10.1016/j.artmed.2019.06.003

  • 27.

    CaoMYangMQinCZhuXChenYWangJet al. Using DeepGCN to identify the autism spectrum disorder from multi-site resting-state data. Biomed Signal Process Control. (2021) 70:103015. doi: 10.1016/j.bspc.2021.103015

  • 28.

    QureshiMNIRyuSSongJLeeKHLeeB. Evaluation of functional decline in Alzheimer’s dementia using 3D deep learning and group ICA for rs-fMRI measurements. Front Aging Neurosci. (2019) 11:8. doi: 10.3389/fnagi.2019.00008

  • 29.

    AlorfAKhanMUG. Multi-label classification of Alzheimer's disease stages from resting-state fMRI-based correlation connectivity data and deep learning. Comput Biol Med. (2022) 151:106240. doi: 10.1016/j.compbiomed.2022.106240

  • 30.

    JiSXuWYangMYuK. 3D convolutional neural networks for human action recognition. IEEE Trans Pattern Anal Mach Intell. (2012) 35:221–31. Available at: https://icml.cc/Conferences/2010/papers/100.pdf

  • 31.

    WoytowiczEJRietschelJCGoodmanRNConroySSSorkinJDWhitallJet al. Determining levels of upper extremity movement impairment by applying a cluster analysis to the Fugl-Meyer assessment of the upper extremity in chronic stroke. Arch Phys Med Rehabil. (2017) 98:456–62. doi: 10.1016/j.apmr.2016.06.023

  • 32.

    SaadZSGottsSJMurphyKChenGJoHJMartinAet al. Trouble at rest: how correlation patterns and group differences become distorted after global signal regression. Brain Connect. (2012) 2:25–32. doi: 10.1089/brain.2012.0080

  • 33.

    ShenKHutchisonRMBezginGEverlingSMcIntoshAR. Network structure shapes spontaneous functional connectivity dynamics. J Neurosci. (2015) 35:5579–88. doi: 10.1523/JNEUROSCI.4903-14.2015

  • 34.

    AndersenSMRapcsakSZBeesonPM. Cost function masking during normalization of brains with focal lesions: still a necessity?NeuroImage. (2010) 53:78–84. doi: 10.1016/j.neuroimage.2010.06.003

  • 35.

    YanC-GWangX-DZuoX-NZangY-F. DPABI: data processing and analysis for (resting-state) brain imaging. Neuroinformatics. (2016) 14:339–51. doi: 10.1007/s12021-016-9299-4

  • 36.

    HindriksRAdhikariMHMurayamaYGanzettiMMantiniDLogothetisNKet al. Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI?NeuroImage. (2016) 127:242–56. doi: 10.1016/j.neuroimage.2015.11.055

  • 37.

    YinDSongFXuDPetersonBSSunLMenWet al. Patterns in cortical connectivity for determining outcomes in hand function after subcortical stroke. PLoS One. (2012) 7:e52727. doi: 10.1371/journal.pone.0052727

  • 38.

    YinDSongFXuDSunLMenWZangLet al. Altered topological properties of the cortical motor-related network in patients with subcortical stroke revealed by graph theoretical analysis. Hum Brain Mapp. (2014) 35:3343–59. doi: 10.1002/hbm.22406

  • 39.

    SakoğluÜPearlsonGDKiehlKAWangYMMichaelAMCalhounVD. A method for evaluating dynamic functional network connectivity and task-modulation: application to schizophrenia. MAGMA. (2010) 23:351–66. doi: 10.1007/s10334-010-0197-8

  • 40.

    CohenJCohenPWestSGAikenLS. Applied multiple regression/correlation analysis for the behavioral sciences. New York: Routledge (2013).

  • 41.

    KingmaD. P.BaJ.Adam: a method for stochastic optimization. [Epubh ahead of preprint] (2014). doi: 10.48550/arXiv.1412.6980

  • 42.

    SimonyanKZissermanA. Very deep convolutional networks for large-scale image recognition. [Epubh ahead of preprint] (2014). doi: 10.48550/arXiv.1409.1556

  • 43.

    HeKZhangXRenSSunJ, eds. Delving deep into rectifiers: surpassing human-level performance on imagenet classification. Proceedings of the IEEE International Conference on Computer Vision; (2015).

  • 44.

    AggarwalCCHinneburgAKeimDA, eds. On the surprising behavior of distance metrics in high dimensional space. Database Theory—ICDT 2001: 8th International Conference London, UK, January 4–6, 2001 Proceedings 8; (2001): Springer.

  • 45.

    ForkertNDVerlegerTChengBThomallaGHilgetagCCFiehlerJ. Multiclass support vector machine-based lesion mapping predicts functional outcome in ischemic stroke patients. PLoS One. (2015) 10:e0129569. doi: 10.1371/journal.pone.0129569

  • 46.

    RehmeAKVolzLJFeisDLEickhoffSBFinkGRGrefkesC. Individual prediction of chronic motor outcome in the acute post-stroke stage: behavioral parameters versus functional imaging. Hum Brain Mapp. (2015) 36:4553–65. doi: 10.1002/hbm.22936

  • 47.

    RondinaJMFilipponeMGirolamiMWardNS. Decoding post-stroke motor function from structural brain imaging. NeuroImage. (2016) 12:372–80. doi: 10.1016/j.nicl.2016.07.014

  • 48.

    TozluCEdwardsDBoesALabarDTsagarisKZSilversteinJet al. Machine learning methods predict individual upper-limb motor impairment following therapy in chronic stroke. Neurorehabil Neural Repair. (2020) 34:428–39. doi: 10.1177/1545968320909796

  • 49.

    FavreIZeffiroTADetanteOKrainikAHommelMJaillardA. Upper limb recovery after stroke is associated with ipsilesional primary motor cortical activity: a meta-analysis. Stroke. (2014) 45:1077–83. doi: 10.1161/STROKEAHA.113.003168

  • 50.

    MarshallRSZarahnEAlonLMinzerBLazarRMKrakauerJW. Early imaging correlates of subsequent motor recovery after stroke. Ann Neurol. (2009) 65:596–602. doi: 10.1002/ana.21636

  • 51.

    KrizhevskyASutskeverIHintonGE. Imagenet classification with deep convolutional neural networks. Adv Neural Inf Process Syst. (2012) 25. Available at: https://proceedings.neurips.cc/paper/2012/file/c399862d3b9d6b76c8436e924a68c45b-Paper.pdf

  • 52.

    PretiMGBoltonTAVan De VilleD. The dynamic functional connectome: state-of-the-art and perspectives. NeuroImage. (2017) 160:41–54. doi: 10.1016/j.neuroimage.2016.12.061

Summary

Keywords

stroke, resting-state functional magnetic resonance imaging, dynamic functional connectivity, ipsilesional primary motor cortex, three-dimensional convolutional neural network

Citation

Gao Y, Xu G, Peng J, Han C, Wu S, Wang M, Wang H and Zhao Z (2025) Assessment of functional decline in stroke patients using 3D deep learning and dynamic functional connectivity based on resting-state fMRI. Front. Neurol. 16:1666991. doi: 10.3389/fneur.2025.1666991

Received

16 July 2025

Accepted

20 October 2025

Published

13 November 2025

Volume

16 - 2025

Edited by

Chuanming Li, Chongqing University Central Hospital, China

Reviewed by

Gaoxing Zheng, University of Texas Southwestern Medical Center, United States

Xingwang Yong, Mass General Brigham, United States

Updates

Copyright

*Correspondence: Zhiyong Zhao, Hewei Wang,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics