Resting-State Neural-Activity Alterations in Subacute Aphasia after Stroke

Linguistic deficits are frequent symptoms among stroke survivors. The neural mechanism of post-stroke aphasia (PSA) was incompletely understood. Recently, resting-state functional magnetic resonance imaging (rs-fMRI) was widely used among several neuropsychological disorders. However, previous rs-fMRI studies of PSA were limited to very small sample size and the absence of reproducibility with different neuroimaging indexes. The present study performed comparisons with static and dynamic amplitude of low-frequency fluctuations (ALFF) and functional connectivity (FC) based on modest sample size (40 PSA and 37 healthy controls). Compared with controls, PSA showed significantly increased static ALFF predominantly in the bilateral supplementary motor area (SMA) and right hippocampus-parahippocampus (R HIP-ParaHip) and decreased static ALFF in right cerebellum. The increased dynamic ALFF in SMA and decreased dynamic ALFF in right cerebellum were also found in PSA. The static and dynamic ALFF in right cerebellum was positively correlated with spontaneous speech. The FC between the SMA and R HIP-ParaHip was significantly stronger in patients than controls and positively correlated with ALFF in bilateral SMA. In addition, the FC between the R HIP-ParaHip and the right temporal was also enhanced in patients and negatively correlated with repetition, naming, and comprehension score. These findings revealed consistently abnormal intrinsic neural activity in SMA and cerebellum, which may underlie linguistic deficits in PSA.


Introduction
Aphasia is a disabling linguistic deficit manifested by impaired speech expression and comprehension caused by injury to the language-dominant hemisphere. Epidemiological studies have found that approximately one-third of stroke survivors suffer from aphasia [1] and that 40% still display significant speech impairment after one year [2]. Furthermore, some patients bear lifetime disability [3]. Patients with post-stroke aphasia (PSA) commonly exhibit poor social communication and low work competence, making it difficult to fully engage in society [4] and placing a substantial burden on individual, families and social services [3]. A comprehensive understanding of the neural mechanism of PSA is essential to guide the treatment for rehabilitation doctor. Evidence from human stroke motor recovery models shows that the first 90 days after stroke are a critical 'window of opportunity' for recovery [5]. This time window also applies to recovery from aphasia [6]. Recent studies showed the better efficacy of speech and language therapy (SLT) on acute and subacute aphasia [7,8], while another report suggested that intensive SLT beginning during acute aphasia could not improve efficacy than routine SLT [9]. Similarly, contradictory results were found in the efficacy of transcranial direct current stimulation on subacute PSA [10]. These inconsistencies may be the consequence of unclear neural mechanism or greater inter-subject heterogeneity. Hence, it is very important to explore the activity of local and remote language-related brain regions in subacute aphasia. Fortunately, neuroimaging studies increase the current understanding of local and whole brain activity.
Functional magnetic resonance imaging (fMRI) represents an ideal tool and is widely applied to reveal the brain mechanisms of neurological diseases [11]. Several fMRI studies have found that the recovery of aphasia mainly depend on the preservation and functional status of language-related brain regions [12,13]. Richter et al. [14] reported that the structure of the right hemisphere, traditionally considered as the non-dominant hemisphere for language, can be integrated into the language network after severe stroke.
Resting-state fMRI(rs-fMRI [15]) uses blood oxygenation level-dependent signals to measure neural activity. Compared with task-based fMRI, it is easy to complete even for patients with severe aphasia. The amplitude of low-frequency fluctuations (ALFF [16]) is a frequently used method reflecting spontaneous activity for the analysis of rs-fMRI data [17]. Previous studies have found significantly increased ALFF in the contralesional cortices within two weeks after stroke [18] that correlate with behavioral measures of language function [19]. Dynamic ALFF (dALFF [20]) measurement is a recently developed imaging application that can reflect the temporal flexibility of resting-state brain networks. Dynamics can capture uncontrolled but recurring patterns in brain networks that cannot be detected by static analysis [21]. A combination of static and dynamic analysis can complement static changes in the pathological states. Reduced dALFF [20] has been reported in the left inferior frontal gyrus triangle of disrupted language networks in acute aphasia patients.
Language functions are implemented in widely distributed neural networks, as measured with resting-state fMRI. Disrupted corresponding network connectivity may underlie language deficits of aphasia [22]. Functional connectivity (FC [23]) is also a fairly direct and simple method to measure the level of temporal co-activation based on the computation of Pearson's correlation coefficients between two anatomically separate brain regions. Previous research showed that the FC was abnormal in acute aphasia compared with health participants [18].
In this study, we combined static and dynamic ALFF to explore the dynamics of spontaneous brain activity in subacute aphasia. Then, we measured FC to assess the cooperation of remote region after stroke-induced disruption. Furthermore, we examined the correlations between neuroimaging alterations and language function scores. We hypothesized that there would be a functional abnormality of undamaged language-related brain regions. These aberrant areas may underlie linguistic deficits in subacute PSA.

Participants
We recruited 48 patients with post-stroke aphasia during subacute stage from the First Affiliated Hospital of Anhui Medical University and Anhui acupuncture hospital (Hefei, China) and 37 healthy controls (HCs) well matched for mean age, sex ratio, and years of education (see Table 1 for a summary of patient and HC demographics). Patients were recruited according to the following inclusion criteria: (i) 18-75 years old, (ii) first stroke in the left hemisphere (see Figure 1 for lesion overlap map), (iii) right-handed, (iv) native Chinese speaker, (v) aphasia persistent at day 1 after stroke, and (vi) post onset of stroke ≤3 months. Patients were excluded according to the following criteria: (i) any past or current neurological disease or family history of hereditary neurological disorders, (ii) mental illness or family history of hereditary mental disorders, (iii) severe dysarthria, (iv) a history of head injury or surgery, (v) alcohol or substance abuse, (vi) cerebral tumor or abscess, (vii) claustrophobia or implants incompatible with MRI, and (viii) failure to complete language testing or fMRI scans. In addition, excessive motion (>3 mm in translation or >3 • in rotation) during scanning and abnormal brain structure were set as exclusion criteria for analysis but no data were excluded for these reasons. Eight patients were excluded for claustrophobia or incomplete language testing or fMRI. Finally, forty patients were enrolled in our study. HCs were recruited using the same exclusion criteria and inclusion criteria unrelated to stroke (i, iii, and iv). stroke in the left hemisphere (see Figure 1 for lesion overlap map), (iii) right-handed, (iv) native Chinese speaker, (v) aphasia persistent at day 1 after stroke, and (vi) post onset of stroke ≤3 months. Patients were excluded according to the following criteria: (i) any past or current neurological disease or family history of hereditary neurological disorders, (ii) mental illness or family history of hereditary mental disorders, (iii) severe dysarthria, (iv) a history of head injury or surgery, (v) alcohol or substance abuse, (vi) cerebral tumor or abscess, (vii) claustrophobia or implants incompatible with MRI, and (viii) failure to complete language testing or fMRI scans. In addition, excessive motion (>3 mm in translation or >3° in rotation) during scanning and abnormal brain structure were set as exclusion criteria for analysis but no data were excluded for these reasons. Eight patients were excluded for claustrophobia or incomplete language testing or fMRI. Finally, forty patients were enrolled in our study. HCs were recruited using the same exclusion criteria and inclusion criteria unrelated to stroke (i, iii, and iv).

Language Assessment
Before undergoing MRI scans, each participant received neurological examinations and neuropsychological assessments. Language deficit was assessed by the same professional language therapist using the Aphasia Battery of Chinese (ABC [24]), a Chinese standardized adaptation of the Western Aphasia Battery. The multidimensional test contains four subtypes, including spontaneous speech (semi-standardized interview, assessment of the fluency, and information of speech; total score 20), auditory comprehension

Language Assessment
Before undergoing MRI scans, each participant received neurological examinations and neuropsychological assessments. Language deficit was assessed by the same professional language therapist using the Aphasia Battery of Chinese (ABC [24]), a Chinese standardized adaptation of the Western Aphasia Battery. The multidimensional test contains four subtypes, including spontaneous speech (semi-standardized interview, assessment of the fluency, and information of speech; total score 20), auditory comprehension (yes or no question, auditory picture matching task, and verbal instruction; total score 230), repetition (words and sentences; total score 100), and naming (simple objects, colors, pictures, and situations; total score 80). The combined score is used to calculate an aphasia quotient (AQ) reflecting the overall severity of language impairment. Patients with AQ below 93.8 [25] were considered as aphasic. The clinical characteristics of aphasia are summarized in Table 1. The study conformed to the tenets of the Declaration of Helsinki and was approved by the Anhui Medical University Ethics Committee (2019H009). All patients and guardians were provided informed consent before enrolling.

Lesion Mapping
We manually traced the outlines of lesions using MRIcron (http://www.mccauslandcenter. sc.edu/mricro/mricron, accessed on 20 March 2022) and then created a lesion mask for each patient in individual T1 image. After spatial normalization (described below), the superposition of all individual lesion masks was used to create a group-level lesion mask (see Figure 1). The group-level lesion mask was removed from the whole-brain mask as a segmentation mask using SPM8.  [27]. We removed the first five functional volumes of data to ensure stable longitudinal magnetization, and corrected the remaining volumes for slice timing and head motion. Images with motion > 2.5 mm in translation or 2.5 • in rotation were rejected. Individual 3D T1 images were co-registered to functional images. Next, we segmented 3D T1 images within the range of the segmentation mask. Structural T1 images were normalized to Montreal Neurological Institute (MNI) space based on the T1 image unified segmentation by 12-parameter nonlinear transformation. In addition, we used cost-function masking to discard the brain area of the group lesion, thereby avoiding divergence during spatial normalization [28]. This process must be conducted in SPM8 and has been used in other brain imaging studies including patients with lesions [29]. After spatial normalization, functional images were nuisance regressed with 24 Friston motion parameters, white matter high signals, cerebrospinal fluid signals, and global signals as regressors. Finally, images were spatially smoothed with a 4 mm full-width at half-maximum isotropic Gaussian kernel.

ALFF and dALFF Analysis
Amplitude of low-frequency fluctuation, a marker of resting-state intrinsic regional activity at the voxel level [30], was calculated using DPARSF. First, Fast Fourier Transform was used to convert the time series into a frequency power spectrum. For each voxel, we then obtained the average square root within the 0.01-0.1 Hz range. The whole-brain ALFF map of each individual was first calculated. Individual ALFF maps were standardized by Fisher's z-score conversion to allow comparison between groups. We defined a specific group mask by subtracting the patient group level lesion mask from the whole-brain template. ALFF analysis was performed within the specific group mask in patients. Considering high heterogeneity of each participant, some of the patients had smaller lesions Brain Sci. 2022, 12, 678 5 of 13 than the group level mask. Therefore, we re-defined another mask in which the voxels with less than 10 percent (4 patients) of patients were preserved to verify our results (see Supplementary Material).
Dynamic ALFF tracks the temporal changes in ALFF in the resting state using the Dynamic BC toolbox (www.restfmri.net/forum/DynamicBC, accessed on 20 March 2022) [31], which is based on the sliding-window analysis method. We segmented the full-length time series into a square sliding window of 20 TRs (48 s) as previous studies have shown that a window size of 30 to 60 s is a reasonable choice for capturing brain dynamics [32]. The step size of the sliding window was 1TR (2.4 s). This process produced 193 windows for each subject. In addition, we used two different window lengths ( Figure S1 in Supplementary Material). The ALFF map obtained by all sliding windows was standardized using Fisher's z-transformation to improve normality. Similarly, dALFF analysis was limited to the specific group mask in patients.

FC Analysis
Functional connectivity was calculated using DPARSF, SPM8 and REST. According to the inter-group comparison of ALFF and dALFF, we defined the brain regions with significantly abnormal regions between two group as regions of interest (ROIs) for FC analyses (see Figure 2). For each participant, Pearson's correlation coefficients were calculated between the mean time series of each ROI and the time series of every other brain voxel. To improve normality, correlation coefficients were transformed to z-values using Fisher's r-to-z transformation, and FC maps were constructed for each individual.
Amplitude of low-frequency fluctuation, a marker of resting-state intrinsic regional activity at the voxel level [30], was calculated using DPARSF. First, Fast Fourier Transform was used to convert the time series into a frequency power spectrum. For each voxel, we then obtained the average square root within the 0.01-0.1 Hz range. The whole-brain ALFF map of each individual was first calculated. Individual ALFF maps were standardized by Fisher's z-score conversion to allow comparison between groups. We defined a specific group mask by subtracting the patient group level lesion mask from the wholebrain template. ALFF analysis was performed within the specific group mask in patients. Considering high heterogeneity of each participant, some of the patients had smaller lesions than the group level mask. Therefore, we re-defined another mask in which the voxels with less than 10 percent (4 patients) of patients were preserved to verify our results (see Supplementary Material).
Dynamic ALFF tracks the temporal changes in ALFF in the resting state using the Dynamic BC toolbox (www.restfmri.net/forum/DynamicBC) [31], which is based on the sliding-window analysis method. We segmented the full-length time series into a square sliding window of 20 TRs (48 s) as previous studies have shown that a window size of 30 to 60 s is a reasonable choice for capturing brain dynamics [32]. The step size of the sliding window was 1TR (2.4 s). This process produced 193 windows for each subject. In addition, we used two different window lengths ( Figure S1 in Supplementary Material). The ALFF map obtained by all sliding windows was standardized using Fisher's z-transformation to improve normality. Similarly, dALFF analysis was limited to the specific group mask in patients.

FC Analysis
Functional connectivity was calculated using DPARSF, SPM8 and REST. According to the inter-group comparison of ALFF and dALFF, we defined the brain regions with significantly abnormal regions between two group as regions of interest (ROIs) for FC analyses (see Figure 2). For each participant, Pearson's correlation coefficients were calculated between the mean time series of each ROI and the time series of every other brain voxel. To improve normality, correlation coefficients were transformed to z-values using Fisher's r-to-z transformation, and FC maps were constructed for each individual. caudate nucleus was negatively correlated with repetition, naming and AQ score; Static ALFF in the right cerebellum was positively correlated with spontaneous speech score (two-tailed, no correction). The r value denotes partial correlation coefficient. Abbreviations: PSA, post-stroke aphasia; Caud, caudate nucleus; INS, insula; Cerebe, cerebellum; PCC, bilateral posterior cingulate.

Statistical Analysis
Two-sample t-tests were performed to identify the differences in standardized static ALFF and dALFF maps between patient and HC group using SPM8. The locations of the peak maxima within significant clusters were identified using REST. The ALFF and dALFF values in regions with differences were extracted. The same procedure was applied to compare FC between patients and HCs. All statistical maps were thresholded using a cluster-level family-wise error-corrected threshold of p < 0.05 (cluster-forming threshold at voxel-level p < 0.001). We performed correlation analysis to explore the relationship between aberrant ALFF and FC and clinical features. Considering high heterogeneity of each participant, the lesion volumes and time of post onset of stroke were regressed using partial correlation analysis. These analyses were performed in regions with statistically significant between-group difference.

Demographic and Clinical Features
The demographic and clinical features of the patients and HCs are summarized in Table 1. There were no significant differences in mean age (t = 0.837, df = 75, p = 0.405), sex ratio (χ 2 = 0.956, df = 1, p = 0.328), and mean educational level (t = −1.795, df = 75, p = 0.077) between the 40 aphasic patients and 37 HCs ( Table 1). The overlap of brain lesions among patients is shown in Figure 1.

ALFF Contrast
Compared to HCs, patients displayed significantly greater ALFF in the R HIP-ParaHip, insula and caudate nucleus, left ParaHip-HIP, left cerebellum, and bilateral SMA, but significantly lower ALFF in the right cerebellum and bilateral posterior cingulate(PCC) ( Table 2) (Figure 2). In the another mask, the results of ALFF analysis were similar (see Results in the Supplementary Material). Aphasia patients also exhibited significantly lower dALFF in the right cerebellum and greater dALFF in the B-SMA (Table 3) (Figure 3).
Similarly, the results of other two different window lengths are shown in Figure S1 (see Results in the Supplementary Material).  Abbreviations: dALFF, dynamic amplitude of low-frequency fluctuation.   Table 3. (B) The z values of dALFF were extracted in significant regions between PSA and HC. (C) dALFF in the right cerebellum was positively correlated with spontaneous speech score (two-tailed, no correction). The r value denotes partial correlation coefficient.

FC Contrast
Function connectivity between the B SMA and R HIP-ParaHip, bilateral dorsomedial prefrontal cortex (DMPFC), and bilateral cerebellum were significantly stronger among aphasia patients compared with HCs (Table 4) (Figure 4). Patients also demonstrated stronger FC between the R HIP-ParaHip and right superior temporal gyrus (R STG), middle temporal gyrus (R MTG), middle occipital gyrus (R MOG), and bilateral cerebellum (Table 4) ( Figure 5). Function connectivity between the B SMA and R HIP-ParaHip, bilateral dorsomedial prefrontal cortex (DMPFC), and bilateral cerebellum were significantly stronger among aphasia patients compared with HCs (Table 4) (Figure 4). Patients also demonstrated stronger FC between the R HIP-ParaHip and right superior temporal gyrus (R STG), middle temporal gyrus (R MTG), middle occipital gyrus (R MOG), and bilateral cerebellum (Table 4) ( Figure 5).  relations between FC values and ALFF values of SMA in PSA. All statistical maps were thresholded using a cluster-level family-wise error-corrected threshold of p < 0.05 (cluster-forming threshold at voxel-level p < 0.001). (A) Warm colors represent regions with increased FC values in aphasia. Numbers below axial slices represent the z-plane coordinates of the MNI space. Further details of these regions are shown in Table 4. (B) The z values of FC were extracted in significant regions between PSA and HC. (C) ALFF in the SMA was positively corrected with FC between the SMA and R HIP-ParaHip (two-tailed, no correction).  Table 4. (B) The z values of FC were extracted in significant regions between PSA and HC. (C) FC between the R HIP-ParaHip and right middle temporal gyrus (R MTG) was negatively correlated with comprehension, repetition, and naming score; FC between the R HIP-ParaHip and right superior temporal gyrus (R STG) and right middle occipital gyrus (R MOG) were negatively correlated with comprehension score in aphasia (two-tailed, no correction). The r value denotes partial correlation coefficient.

Discussion
Linguistic ability depends on a distributed network of structures that is frequently disrupted by stroke, resulting in post-stroke aphasia. In the present study, we first applied two distinct neuroimaging indexes (static and dynamic ALFF) to explore regional neural alterations in subacute PSA. Our results suggested that PSA showed significantly increased static and dynamic ALFF in the bilateral SMA and decreased static and dynamic  Table 4. (B) The z values of FC were extracted in significant regions between PSA and HC. (C) FC between the R HIP-ParaHip and right middle temporal gyrus (R MTG) was negatively correlated with comprehension, repetition, and naming score; FC between the R HIP-ParaHip and right superior temporal gyrus (R STG) and right middle occipital gyrus (R MOG) were negatively correlated with comprehension score in aphasia (two-tailed, no correction). The r value denotes partial correlation coefficient.

Discussion
Linguistic ability depends on a distributed network of structures that is frequently disrupted by stroke, resulting in post-stroke aphasia. In the present study, we first applied two distinct neuroimaging indexes (static and dynamic ALFF) to explore regional neural alterations in subacute PSA. Our results suggested that PSA showed significantly increased static and dynamic ALFF in the bilateral SMA and decreased static and dynamic ALFF in right cerebellum. The consistent findings about static and dynamic ALFF in SMA and cerebellum may indicate crucial role in the linguistic loss or recovery of PSA.
Mounting evidence indicates that the SMA is involved in language processing [33], particularly of action language [34], and speech motor control [35]. Saur et al. reported stronger activation of the SMA in patients with subacute aphasia during auditory comprehension tasks compared to controls [36]. Due to limitations of our design, we are incapable of concluding the cognitive mechanism of the enhanced ALFF in SMA. The most likely explanation may be the maladaptive and diaschisis effect after large-regions lesion. Another slim odd is that this enhanced spontaneous activity in SMA played a compensatory role in subacute PSA. This deduction is based on an increasingly important network contribute to language recovery, i.e., domain-general networks [37,38], which includes SMA as key node. Of course, it is difficulty in understanding biological interpretations with brain susceptibility imaging, including rs-fMRI. It is also notable that the SMA is closely related with apraxia of speech [35,39], which frequently co-occur with aphasia among stroke survivors. However, the apraxia of speech was not assessed in present study, which may introduce confounding effect on our results. Future studies with more powerful tools are needed to probe the cognitive mechanism of local and connective alterations in SMA among PSA.
This study also revealed reduced ALFF/dALFF in the right cerebellum of aphasia patients. Traditional viewpoints focus the role of cerebellum on motor function. Recent evidence suggest the link between the cerebellum and a variety of linguistic processing, including language learning, semantic processing, and verbal memory [40]. Task-fMRI study in healthy subjects revealed that the activated cerebellar areas during language tasks present a right lateralizing, which is consistent with crossed projections of cerebrocerebellar circuit [41,42]. The linguistic role of right cerebellum was also extensively studied in pathological conditions, including post-stroke aphasia. Li et al. [16] reported that the ALFF of right cerebellum was weaker in patients with aphasia compared to HCs. Another study found decreased ALFF in right cerebellum and weaker cerebrocerebellar connection [43]. Consistent with these findings, our results also replicated the decreased ALFF of right cerebellum in PSA. Moreover, using a relatively novel index (dALFF), we also found abnormity in the right cerebellum. The possible explanation is diaschisis effect [44] for this abnormal activity. Diaschisis refers to the neuronal dysfunction in undamaged regions distant to the lesion. We speculate that the decreased ALFF in right cerebellum may result from the poor interaction with cerebrum. Previous study found the relationship between the impaired language performance and the diaschisis effect [45] in the subacute aphasia. Similar diaschisis effect [38] has found in acute stage. Interestingly, recent studies also found that transcranial direct current stimulation to the right cerebellum improved language function [46]. These findings implied a potential therapeutic role of the right cerebellum in PSA.
Besides the abnormities in SMA and cerebellum found in present study, several temporal regions also present aberrant intrinsic activity in PSA, including increased ALFF in R HIP-ParaHip and increased FC between the R HIP-ParaHip and the R STG, R MTG and R MOG. Piai et al. reported that the hippocampus actively contributes to language comprehension [47]. Extensive studies have validated the critical role of the left temporal network [48,49] in language processing, including comprehension [50] and production [51] of speech. But temporal regions identified in present study mainly located in the right hemisphere. The role of the right hemisphere in aphasia has been debated for more than one century and remains controversial [52]. According to most studies supporting the "disinhibition hypothesis", the increased activity may be related with severe lesion of left hemisphere [53,54]. However, a small group of studies also suggested the facilitative effect of the right hemisphere (right hemisphere hypothesis [55]). Based on this hypothesis, the enhanced connections within right temporal regions may reflect a compensation. Several factors confound the observation about the role of right hemisphere [38], including the location and extension of the left lesions, the time elapsed since the onset of aphasia. Studies with more homogeneous samples may be helpful to resolve this contradiction.
It is impossible to ignore the high heterogeneous sample in present study. Furthermore, additional limitations of our study should also be considered. First, a cross-sectional rather than longitudinal design limited the biological interpretations for our findings. Second, the lesion voxels for all patients were excluded from final analysis, which may lose sight of the alterations of perilesional tissue. Third, a relatively lenient statistical threshold was used in our study to test for the correlations between neural alterations and clinical variables. Correction for multiple correlations was not performed, which may potentially increase the risk of obtaining false-positive results. Finally, we did not assess the apraxia of speech, which may have an effect on results regarding the SMA.

Conclusions
Our results revealed the abnormal local activity in SMA and right cerebellum, as well as its dynamics. Furthermore, the increased connections were also found within temporal regions predominantly in the right. Although the biological interpretations remain unclear, these findings indicated a key role of SMA and right cerebellum in the language loss or recovery among PSA patients.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/brainsci12050678/s1, Figure S1: Brain regions displaying significant differences in dALFF between patients and healthy controls; Figure S2: Brain regions displaying significant differences in static amplitude of low-frequency fluctuation (ALFF) between patients and healthy controls (HCs) and correlations between ALFF values and clinical scores in subacute aphasia patients; Figure S3: Brain regions displaying significant differences in dynamic ALFF (dALFF) between patients and HCs and correlations between dALFF values and clinical scores in subacute aphasia patients; Figure S4: Abnormalities in SMA-based FC between PSA and HC and correlations between FC values and ALFF values of SMA in PSA. Brain regions demonstrating FC differences between PSA and HC; Figure S5: Abnormalities in R HIP-based functional connectivity (FC) between PSA and HC. Brain regions demonstrating FC differences between PSA and HC; Table S1: Regions showing significant differences in static ALFF between patients and HCs; Table S2: Regions showing significant differences in dALFF between patients and HCs; Table S3: SMA-based FC abnormalities in subacute aphasia patients; Table S4: R HIP-based FC abnormalities in subacute aphasia patients.

Institutional Review Board Statement:
The study adhered to the tenets of the Declaration of Helsinki and was approved by the Ethics Committee of Anhui Medical University.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.