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

Article metrics

View details

952

Views

85

Downloads

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

Sequential brain MRI images from Z=7 to Z=47 show overlapping areas of stroke lesion. A sagittal view on the right features blue slice markers. A color bar below ranges from dark red (0) to white (32), and the images are labeled with L for left and R for right.

Overlap map of stroke lesions across all participants. Color intensity indicates the number of patients with overlapping lesions at each voxel. Axial slices are displayed along the z-axis from slice 7 to slice 47 in Montreal Neurological Institute (MNI) space (MNI, Montreal Neurological Institute; R, right; L, left).

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

Brain maps showing functional connectivity evolving over time, represented by color gradients from blue to red. Each model is labeled W with subscript indices, connected by a red arrow indicating time flow, T. A vertical color scale is on the right.

DFC maps across different sliding windows for a representative participant. Each map illustrates the time-varying connectivity between the left primary motor cortex (M1) and other brain regions. The color bar denotes Fisher’s Z-transformed correlation coefficients (T, time; W, window).

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

Flowchart of a convolutional neural network architecture for processing 3D dFC maps. The sequence includes Conv3 layers with varying filters, MaxPooling, and dropout layers, culminating in fully connected (FC) layers and outputs labeled SSP and MSP. Components are color-coded for clarity, including detailed inset modules illustrating operations like Batch Normalization and ReLU activation.

VGG-Net-based 3D-CNN architecture. SSP, severely stroke patients; MSP, mild stroke patients; dFC, dynamic functional connectivity.

Table 1

Layer Feature map Stride Kernel Activation structure
3D Input 61 × 73 × 61
Convolution 64 1x1x1 3x3x3 Conv
Convolution 64 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Maxpool 2x2x2 2x2x2
Convolution 128 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 128 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Maxpool 2x2x2 2x2x2
Convolution 256 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 256 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 256 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Maxpool 2x2x2 2x2x2
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Maxpool 2x2x2 2x2x2
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Convolution 512 1x1x1 3x3x3 Batchnorm + ReLU + Conv
Maxpool 1x1x1 2x2x2
Fully Connected 4,096 Dropout Rate o.7 ReLU
Fully Connected 4,096 Dropout Rate o.7 ReLU
Output
Fully Connected 2 ReLU
Softmax
Classification Layer Argmax

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)a 54.03 ± 11.95 51.28 ± 10.73 0.392
Sex (male: female) 8/32 5/24 0.513
Lesion volume (ml)a 6.7 ± 4.73 4.97 ± 4.70 0.137
Duration of illness (month)a 8.01 ± 6.52 9.40 ± 7.10 0.410
Lesion side (left: right) 22/18 15/14 0.490
UL-FMAa 8.85 ± 4.49 31.4 ± 7.84 0.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

Fold Train ACC (%) Test ACC (%) p-value AUC Spec (%) Sen (%) F-score (%) BAC (%)
1 99.72 99.45 <0.001 0.9994 99.4 99.51 99.34 99.45
2 99.96 99.65 <0.001 0.9999 99.88 99.34 99.59 99.61
3 99.72 99.69 <0.001 0.9999 99.76 99.67 99.67 99.71
4 100 99.93 <0.001 1 99.88 100 99.91 99.94
5 100 100 <0.001 1 100 100 100 100
6 99.99 99.86 <0.001 1 99.88 99.84 99.84 99.86
7 99.95 100 <0.001 1 100 100 100 100
8 99.92 99.51 <0.001 0.9999 99.76 99.17 99.42 99.47
9 100 99.93 <0.001 1 100 99.83 99.92 99.92
10 99.97 100 <0.001 1 100 100 100 100
Mean ± SD 99.23 ± 0.10 99.80 ± 0.20 0.9999 ± 0.0002 99.86 ± 0.18 99.74 ± 0.29 99.78 ± 0.24 99.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

Index State Group comparisons Brain regions (AAL) Peak MNI coordinates Cluster voxels Peak T values
X Y Z
A 2 MSP > SSP Angular. R 60 −66 21 52 2.9
B 2 MSP > SSP MOF. R 6 33 −12 161 3.37
C 2 MSP > SSP SFG. R 30 60 6 86 3.01
D 2 MSP > SSP SFG. R 18 33 39 58 2.64
E 2 MSP > SSP SPL. L −36 −66 57 103 2.92
F 3 MSP > SSP SMA. L −3 0 69 50 2.89
G 3 MSP > SSP STP. L −42 −6 −15 127 3.136
H 3 MSP > SSP Thalamus. L 3 −6 15 126 3.15
I 3 MSP < SSP Angular. L −39 −63 42 471 3.75
J 3 MSP < SSP Cerebellum 3 −45 −24 77 2.88
K 3 MSP < SSP Rectus. R 9 39 −12 77 3.31
L 3 MSP < SSP MTG. R 60 −63 18 129 3.52
M 3 MSP < SSP Pre. L −3 −63 36 102 3.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

Brain scan images in a grid layout, each showing different sections with highlighted areas in red and blue. The images are labeled with states and regions such as Angular, MOF, SFG, SPL, SMA, STP, Thalamus, Rectus, MTG, and Pre, along with peak MNI coordinates. A color scale ranges from blue (-3.75) to red (3.44) indicating T-value intensity.

Discriminative features of selected states. (A–E) correspond to state 2, and (F–M) correspond to state 3, respectively. 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. Red indicates regions where functional connectivity was higher in MSP than in SSP; blue indicates the opposite (MSP, mild stroke patients; SSP, severe stroke patients).

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.

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

  • 2.

    Nichols-Larsen DS Clark P Zeringue A Greenspan A Blanton S . Factors influencing stroke survivors’ quality of life during subacute recovery. Stroke. (2005) 36:1480–4. doi: 10.1161/01.STR.0000170706.13595.4f

  • 3.

    Wang L Yu C Chen H Qin W He Y Fan F et al . Dynamic functional reorganization of the motor execution network after stroke. Brain. (2010) 133:1224–38. doi: 10.1093/brain/awq043

  • 4.

    Zhao Z Wu J Fan M Yin D Tang C Gong J et 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.

    Griffis JC Metcalf NV Corbetta M Shulman GL . Structural disconnections explain brain network dysfunction after stroke. Cell Rep. (2019) 28:2527–40. doi: 10.1016/j.celrep.2019.07.100

  • 6.

    Kwakkel G Veerbeek JM van Wegen EE Wolf SL . Constraint-induced movement therapy after stroke. Lancet Neurol. (2015) 14:224–34. doi: 10.1016/S1474-4422(14)70160-7

  • 7.

    Hirakawa Y Takeda K Tanabe S Koyama S Motoya I Sakurai H et 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.

    Gladstone DJ Danells CJ Black SE . 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.

    Meyer AF Jaasko L Leyman I Olsson S Steglind S . The post stroke hemiplegic patient. I. A method for evaluation of physical performance. Scand J Rehabil Med. (1975) 7:13–31.

  • 10.

    Gabel M Gilad-Bachrach R Renshaw E Schuster A , eds. Full body gait analysis with Kinect. 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society; (2012): IEEE.

  • 11.

    Kurillo G Chen A Bajcsy R Han JJ . Evaluation of upper extremity reachable workspace using Kinect camera. Technol Health Care. (2013) 21:641–56. doi: 10.3233/THC-130764

  • 12.

    Park C-h Chang WH Ohn SH Kim ST Bang OY Pascual-Leone A et 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-Butel AG Lin G Shiner CT McNulty PA . 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.

    Lu J Liu H Zhang M Wang D Cao Y Ma Q et 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 Meer MP Otte WM van der Marel K Nijboer CH Kavelaars A van der Sprenkel JWB et 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.

    Zhao Z Wang X Fan M Yin D Sun L Jia J et 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.

    Cassidy JM Mark JI Cramer SC . Functional connectivity drives stroke recovery: shifting the paradigm from correlation to causation. Brain. (2022) 145:1211–28. doi: 10.1093/brain/awab469

  • 18.

    Carter AR Astafiev SV Lang CE Connor LT Rengachary J Strube MJ et al . Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke. Ann Neurol. (2010) 67:365–75. doi: 10.1002/ana.21905

  • 19.

    Allen EA Damaraju E Plis SM Erhardt EB Eichele T Calhoun VD . Tracking whole-brain connectivity dynamics in the resting state. Cereb Cortex. (2014) 24:663–76. doi: 10.1093/cercor/bhs352

  • 20.

    Calhoun VD Miller R Pearlson G Adalı 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.

    Li R Liao W Yu Y Chen H Guo X Tang YL et 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.

    Chen J Sun D Shi Y Jin W Wang Y Xi Q et 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.

    Damaraju E Allen EA Belger A Ford JM McEwen S Mathalon D et 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.

    Duc NT Ryu S Qureshi MNI Choi M Lee KH Lee B . 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.

    Shamrat FJM Akter S Azam S Karim A Ghosh P Tasnim Z et 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.

    Qureshi MNI Oh J Lee B . 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.

    Cao M Yang M Qin C Zhu X Chen Y Wang J et 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.

    Qureshi MNI Ryu S Song J Lee KH Lee B . 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.

    Alorf A Khan MUG . 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.

    Ji S Xu W Yang M Yu K . 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.

    Woytowicz EJ Rietschel JC Goodman RN Conroy SS Sorkin JD Whitall J et 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.

    Saad ZS Gotts SJ Murphy K Chen G Jo HJ Martin A et 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.

    Shen K Hutchison RM Bezgin G Everling S McIntosh AR . Network structure shapes spontaneous functional connectivity dynamics. J Neurosci. (2015) 35:5579–88. doi: 10.1523/JNEUROSCI.4903-14.2015

  • 34.

    Andersen SM Rapcsak SZ Beeson PM . 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.

    Yan C-G Wang X-D Zuo X-N Zang Y-F . DPABI: data processing and analysis for (resting-state) brain imaging. Neuroinformatics. (2016) 14:339–51. doi: 10.1007/s12021-016-9299-4

  • 36.

    Hindriks R Adhikari MH Murayama Y Ganzetti M Mantini D Logothetis NK et 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.

    Yin D Song F Xu D Peterson BS Sun L Men W et 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.

    Yin D Song F Xu D Sun L Men W Zang L et 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 Ü Pearlson GD Kiehl KA Wang YM Michael AM Calhoun VD . 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.

    Cohen J Cohen P West SG Aiken LS . Applied multiple regression/correlation analysis for the behavioral sciences. New York: Routledge (2013).

  • 41.

    Kingma D. P. Ba J. Adam: a method for stochastic optimization. [Epubh ahead of preprint] (2014). doi: 10.48550/arXiv.1412.6980

  • 42.

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

  • 43.

    He K Zhang X Ren S Sun J , eds. Delving deep into rectifiers: surpassing human-level performance on imagenet classification. Proceedings of the IEEE International Conference on Computer Vision; (2015).

  • 44.

    Aggarwal CC Hinneburg A Keim DA , 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.

    Forkert ND Verleger T Cheng B Thomalla G Hilgetag CC Fiehler J . 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.

    Rehme AK Volz LJ Feis DL Eickhoff SB Fink GR Grefkes C . 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.

    Rondina JM Filippone M Girolami M Ward NS . Decoding post-stroke motor function from structural brain imaging. NeuroImage. (2016) 12:372–80. doi: 10.1016/j.nicl.2016.07.014

  • 48.

    Tozlu C Edwards D Boes A Labar D Tsagaris KZ Silverstein J et 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.

    Favre I Zeffiro TA Detante O Krainik A Hommel M Jaillard A . 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.

    Marshall RS Zarahn E Alon L Minzer B Lazar RM Krakauer JW . Early imaging correlates of subsequent motor recovery after stroke. Ann Neurol. (2009) 65:596–602. doi: 10.1002/ana.21636

  • 51.

    Krizhevsky A Sutskever I Hinton GE . 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.

    Preti MG Bolton TA Van De Ville D . 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