Corrigendum: Modulations of static and dynamic functional connectivity among brain networks by electroacupuncture in post-stroke aphasia
- 1Department of Neurology, Dongzhimen Hospital, Beijing University of Chinese Medicine, Beijing, China
- 2Key Laboratory of Chinese Internal Medicine Ministry of Education, Beijing University of Chinese Medicine, Beijing, China
- 3Institute for Brain Disorders, Beijing University of Chinese Medicine, Beijing, China
- 4Department of Rehabilitation Sciences, The Hong Kong Polytechnic University, Hong Kong, Hong Kong SAR, China
- 5The State Key Laboratory of Brain and Cognitive Sciences, The University of Hong Kong, Hong Kong, Hong Kong SAR, China
- 6Department of Psychiatry and Psychotherapy, Comprehensive Center for Clinical Neurosciences and Mental Health, Medical University of Vienna, Vienna, Austria
- 7Division of Nuclear Medicine, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
Introduction: Post-stroke aphasia (PSA) is a language disorder caused by left hemisphere stroke. Electroacupuncture (EA) is a minimally invasive therapeutic option for PSA treatment. Tongli (HT5) and Xuanzhong (GB39), two important language-associated acupoints, are frequently used in the rehabilitation of patients with PSA. Preliminary evidence indicated functional activation in distributed cortical areas upon HT5 and GB39 stimulation. However, research on the modulation of dynamic and static functional connectivity (FC) in the brain by EA in PSA is lacking.
Method: This study aimed to investigate the PSA-related effects of EA stimulation at HT5 and GB39 on neural processing. Thirty-five participants were recruited, including 19 patients with PSA and 16 healthy controls (HCs). The BOLD signal was analyzed by static independent component analysis, generalized psychophysiological interactions, and dynamic independent component analysis, considering variables such as age, sex, and years of education.
Results: The results revealed that PSA showed activated clusters in the left putamen, left postcentral gyrus (PostCG), and left angular gyrus in the salience network (SN) compared to the HC group. The interaction effect on temporal properties of networks showed higher variability of SN (F = 2.23, positive false discovery rate [pFDR] = 0.017). The interaction effect on static FC showed increased functional coupling between the right calcarine and right lingual gyrus (F = 3.16, pFDR = 0.043). For the dynamic FC, at the region level, the interaction effect showed lower variability and higher frequencies of circuit 3, with the strongest connections between the supramarginal gyrus and posterior cingulum (F = 5.42, pFDR = 0.03), middle cingulum and PostCG (F = 5.27, pFDR = 0.036), and triangle inferior frontal and lingual gyrus (F = 5.57, pFDR = 0.026). At the network level, the interaction effect showed higher variability in occipital network–language network (LN) and cerebellar network (CN) coupling, with stronger connections between the LN and CN (F = 4.29, pFDR = 0.042). Dynamic FC values between the triangle inferior frontal and lingual gyri were anticorrelated with transcribing, describing, and dictating scores in the Chinese Rehabilitation Research Center for Chinese Standard Aphasia Examination.
Discussion: These findings suggest that EA stimulation may improve language function, as it significantly modulated the nodes of regions/networks involved in the LN, SN, CN, occipital cortex, somatosensory regions, and cerebral limbic system.
Introduction
Post-stroke aphasia (PSA) is a clinical syndrome caused by a left hemisphere stroke that results in the loss of language skills and consequently has an impact on daily life (1). Even mild aphasia can have a negative effect on functional outcomes, such as mood, quality of life, and ability to return to work (2). Language is a critical cognitive skill supported by large-scale brain networks (3). In addition to the direct effect of focal lesions on important cortical regions, damage to other remote areas within language networks and non-language-specific networks also leads to the occurrence and development of aphasia (4, 5). It has been suggested that focal stroke lesions can affect language comprehension by altering the functional connectivity (FC) of multiple networks and subnetworks in PSA (6–8).
Acupuncture is an ancient Chinese treatment that has been systematically used for 2000 years (9). It is rapidly gaining recognition for its therapeutic properties in the treatment of PSA and several other neurological conditions (10–13). Electroacupuncture (EA) is a modern form of acupuncture that features a small current passing between pairs of acupuncture needles. To observe the immediate effect of acupuncture on brain activity, several recent studies have applied EA during functional magnetic resonance imaging (fMRI) (14–16). Block design is one of the most commonly used scanning methods for exploring the potential mechanisms of EA in fMRI studies (17–20). Tongli (HT5) and Xuanzhong (GB39) are important language-implicated acupoints used in Chinese medicine. They are commonly used in PSA rehabilitation therapy because they are tailored to language-processing systems. Previous fMRI studies of EA treatment demonstrated that stimulation at HT5 and GB39 resulted in activation of language regions and the somatosensory cortex in both cerebral hemispheres and that it might modulate speech function through effects on brain networks in healthy individuals (21–23).
It is worth noting that the emergence of acupuncture as a treatment stem from its function in patients rather than in healthy people (24). Modulation of information processing in the brain through acupuncture may differ between patients with PSA and healthy individuals and occurs mainly in disorder-related areas (24). A direct comparison of the effect on neural networks between patients with PSA and healthy individuals involved in the processing of EA stimulation at HT5 and GB39 is virtually absent from the present research. Moreover, although the aforementioned studies have focused on localizing a series of brain regions showing activation patterns during the processing of EA stimulation in healthy individuals, a rapidly growing body of neuroimaging literature suggests that the effect of acupuncture in patients with stroke may be attributed to modulation of disrupted patterns of the whole-brain network rather than in one or two confined brain regions (25, 26).
Therefore, this study aimed to explore and compare how EA stimulation modulates FC and temporal properties in patients with PSA and healthy individuals. To address this question, independent component analysis (ICA) and general psychophysiological interaction (gPPI) were used to reveal static and dynamic FC changes at the levels of the region of interest (ROI) and large-scale network.
Materials and methods
Participants and ethics statements
Based on previous study results, 35 right-handed participants were recruited between August 2013 and February 2021, including 19 patients with PSA (33–78 years, 13 men) and 16 demographically matched healthy volunteers (24–63 years, right-handed, eight men). Patients with PSA were recruited from the Department of Neurology at Dongzhimen Hospital with complaints of language disorders. Healthy controls (HCs) were recruited from the local community.
The Medical Research Ethics Committee of Dongzhimen Hospital (reference number: ECPJ-BDY-2015-04) approved this study, and all participants provided written consent to participate. An overview of the study design is shown in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/956931/fneur-13-956931-HTML-r1/image_m/fneur-13-956931-g001.jpg)
Figure 1. Flow diagram. PSA, post stroke aphasia; gPPI, general psychophysiological interaction; dyn-ICA, dynamic independent component analysis; static ICA, static independent component analysis; CRRCAE, Chinese Rehabilitation Research Center for Chinese Standard Aphasia Examination; ANOVA, analysis of variance; ROI, region of interest; MVPA, multi-voxel pattern analysis.
Clinical evaluation and speech and language testing
Demographic factors, sex, age, handedness, years of education (Edu), and scores on the Boston Diagnostic Aphasia Examination were reviewed for each patient. Twelve patients with motor aphasia were diagnosed by the Chinese Rehabilitation Research Center for Chinese Standard Aphasia Examination (CRRCAE), a battery of language tests designed according to the Mandarin rules. As a result of testing the reliability and validity of the CRRCAE in a prior study, it was found to have good sensitivity and reliability and may be used as a quantitative table for the diagnosis of aphasia in Mandarin speakers (27, 28). Thirty subscales including comprehension, speaking, writing, calculation, copy, and repeating comprise the scale (28).
Stimuli and scanning procedure
Acupuncture procedures and needling sensation recording
Both the PSA and HC groups received EA stimulation at the HT5 and GB39 acupoints during fMRI acquisition. Before the start of the scan, the participants positioned themselves on the fMRI scanner bed on their backs, and needles were placed at the GB39 and HT5 acupoints. According to the “Name and Location of Acupoints” (GB/T 12346-2006), two acupoints were located on both sides. HT5 was situated radially to the flexor carpi ulnaris tendon on the anteromedial side of the forearm, 33 mm proximal to the palmar wrist crease, with insertion depths ranging from 10 to 30 mm. GB39 was needled at an insertion depth of 33 mm, 100 mm above the external malleolus tip, on the anterior fibula border (Figure 2A). A professional acupuncturist performed acupuncture. All participants reported their experience (“Deqi”) with acupuncture stimulation. Deqi featured aching, pressure, heaviness, fullness, and numbness among other feelings (23–29). The acupuncturist used 0.40 × 40-mm sterile silver acupuncture needles (Guizhou, China) with the EA technique. Han's acupoint nerve stimulator (model LH-202H) was situated outside the fMRI room, with one end of the acupoint wire linked to the acupuncture needle handle of the HT5 acupoint, and the other end connected to the acupuncture needle handle of the GB39 acupoint. The EA frequency was 2 Hz, and the electric current was 2 mA. As previously mentioned, the stimulation waveform is the dilatational wave (23).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/956931/fneur-13-956931-HTML-r1/image_m/fneur-13-956931-g002.jpg)
Figure 2. Location of acupuncture points and electroacupuncture stimulation procedure. (A) Location of HT5 and GB39. (B) n-Block design of task-based fMRI. R, rest; S, stimulation.
fMRI scanning
All fMRI data were acquired using the 3 Tesla Siemens MRI scanner (Erlangen, Germany) at Dongzhimen Hospital. A high-resolution T1-weighted structural image was obtained using an isotropic multi-echo magnetization-prepared rapid acquisition sequence (repetition time: 1900 ms, echo time: 2.13 ms, field of view: 250 mm, flip angle: 90°, voxel size: 1.0 × 1.0 × 1.0 mm, slice thickness: 1.0 mm, matrix: 256 × 256, number of slices: 176). Functional images were acquired using a single shot gradient-recalled echo planar imaging sequence (31 interleaved axial slices, repetition time: 2,000 ms, echo time: 30 ms, flip angle: 90°, field of view: 225 × 225 mm2, interslice gap: 0.7 mm, matrix: 64 × 64 mm2, slice thickness: 3.5 mm, voxel size: 3.5 × 3.5 × 3.5 mm3). The fMRI experimental design was a classic resting/stimulation block design divided between a 30-s resting interval and a 30-s stimulation period; such blocks were repeated for six cycles, six acupuncture blocks, and six rest blocks out of a total of 12 blocks that were completed in one session (30) (Figure 2B).
Data processing
The Resting-State fMRI Data Analysis Toolkit (REST plus, version 1.24; http://www.restfmri.net/forum/RESTplus) (31) and task-fMRI images were preprocessed using the Statistical Parametric Mapping 12 software suite (SPM12, http://www.fil.ion.ucl.ac.uk/spm). First, we discarded the first 10 volumes for each participant to avoid the negative impact of magnetic disequilibrium on data quality. Subsequently, all functional images were preprocessed [slice acquisition time, head motion, and a 6-mm Gaussian kernel with full width at half maximum (FWHM) smoothing]. All data had head motions < 3 mm or 3°. Next, using Advanced Normalization Tools (ANTS1.9, https://www.nitrc.org/projects/ants), the functional images were spatially normalized to the Montreal Neurological Institute space in the Linux platform. Temporal filtering was also applied, which allowed a frequency band of 0.01–0.1 Hz to pass. Subsequently, each participant was subjected to linear detrending to remove any residual effect of low-frequency drift or high-frequency physiological noise.
Statistical analysis
All statistical analyses were performed using the CONN toolbox (32). We examined the effects of EA on FC using a condition main effect (group: any effect, condition: EA stimulation > rest at baseline), the PSA-related effect using a group main effect (group: PSA-HC, condition: any effect), and the interaction effect using a two-way crossed analysis of variance (ANOVA), with years of Edu, age, and sex as covariates, with a within-subjects factor (condition: EA stimulation vs. rest at baseline) and a between-subjects factor (group: patients with PSA vs. HCs). At a threshold of p < 0.05, the false discovery rate-corrected results were considered statistically significant.
Static-independent component analysis
Implemented in CONN 18b, static-ICA employs the Fast-ICA method and GICA3 back projection to estimate independent components (ICs) and subject-level spatial map estimation (33). Twenty estimated components were chosen (34). The spatial maps and time courses were normalized to z-scores for future analyses.
Spatial components
The dorsal attention network (DAN), default mode network (DMN), sensorimotor network (SMN), salience network (SN), visual network (VN), frontoparietal network, language network (LN), and cerebellar network (CN) were identified within each component using spatial sorting and a correlational spatial match-to-template approach (spatial correlation and spatial overlap of suprathreshold areas/dice coefficient) (35). Nine components were chosen from a total of 20 ICs to represent the brain networks (Figure 3A). Among them, instead of calling network 9 the VN, we chose to call it the occipital network, as there was no visual stimulus in this study. From now on, we shall refer to components as “networks” and we will identify components by referring to their label (e.g., DMN); component descriptions are listed in Supplementary Table 1. For the spatial map, we examined how the PSA-related network structure differed between the two groups using a voxel-wise F-test (group main effect), as in previous studies (36).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/956931/fneur-13-956931-HTML-r1/image_m/fneur-13-956931-g003.jpg)
Figure 3. Effect of electroacupuncture stimulation on static independent component analysis. (A) Nine spatial maps of intrinsic resting-state networks in 32 subjects. (B) Changes in temporal properties varies of SN within each group and condition. (C) Spatial map of SN for the contrast of post stroke aphasia group vs. healthy control group. (D) Changes in temporal properties varies of VN within each group and condition. DMN, default mode network; SMN, sensorimotor network; VN, visual network; SN, salience network; DAN, dorsal attention network; FPN, frontoparietal network; LN, language network; CN, cerebellar network.
Temporal components
First, for each individual, a set of nine beta coefficients for each of the two condition regressors was derived, indicating the degree to which the task regressor relates to a specific network. Then, to verify the interaction effect (group × condition) on BOLD signal fluctuations in each brain network, we performed a two-way crossed ANOVA test for the nine remaining networks using the temporal properties implemented in the CONN toolbox.
GPPI-based FC analyses
Unsmoothed (37) but preprocessed (refer to data processing) fMRI data were inputted into the CONN toolbox for additional preprocessing and gPPI analysis (38). Before conducting the gPPI analysis, we used a component-based noise correction method to remove artifacts from the fMRI data. To investigate how EA stimuli alter brain functional architecture at the ROI (Supplementary Figure 1) and large-scale network levels (static ICA results), we used gPPI analysis according to Lee Masson et al. (23–39). The regression model used in the gPPI analysis was as follows: Y = X1 × β + X2 × βG + error, where Y = time-series BOLD response of the seed, X1 = the hemodynamic response function (HRF) convolved gPPI terms (i.e., time-series BOLD response of the seed × the psychological regressor as dictated by onset and duration of EA stimulus), β = the strength of FC (i.e., the beta estimates of gPPI terms), X2 = X1 × HRF convolved psychological regressor × covariates of no interest (age, Edu, and sex regressors), and βG = the beta estimates of the HRF-convolved time-series BOLD response of the seeds. Individual results were transformed into z-scores using the Fisher z-transformation for the group-level analysis. First, we evaluated how FC was modulated by the stimuli in each group using a one-sample paired t-test. Second, we used a two-way crossed ANOVA to assess the interaction effect on FC strength (40).
Dynamic independent component analysis
Dynamic connectivity measures were used to examine and characterize the sources of FC variability. Dyn-ICA runs an ICA on the connection time series and returns ICs that most accurately represent FC modulation over time (41, 42). For first-level signal processing, Dyn-ICA matrices (circuits) are created using dynamic connectivity measures. Seeds used to build the matrices were obtained from our previous multi-voxel pattern analysis results (23) (Supplementary Figure 1) and static ICA results in this study (Figure 3A). Then, the aggregated data were divided into 12 circuits using Dyn-ICA with 30 smoothing kernels, in line with a previous study's method (43). This processing produces multiple outputs, including the individual subject-level matrices gamma and the variability and frequency of the dynamic circuit time series. For the second level, we statistically evaluated (1) the frequency of each component (temporal component time-series frequency averaged across all participants) and (2) the variability of each component (temporal component time-series standard deviation averaged across all participants), and the spatial properties (defined as gPPI interaction terms). A two-way crossover ANOVA test was used to assess the interaction effect of temporal data on the frequency and variability of dynamic circuits with Edu, age, and sex as covariates of no interest. Finally, a main effect of condition was performed to assess the EA stimulus-dependent spatial characteristics of circuits with significant temporal properties.
Correlation analysis between behavioral data and fMRI data during EA stimulation
To evaluate the correlation between brain network modulation during EA and behavioral data (CRRCAE scores) in patients with PSA, we conducted a rank correlational analysis. Please note that the correlation analyses were exploratory; thus, a p-value < 0.05 with no corrections for multiple comparisons was implemented, in line with a previous study's methods (44). Correlation analysis was performed using the R software package (The R Project for Statistical Computing, Vienna, Austria).
Results
Clinical characteristics and indices
Detailed clinical characteristics and indices are summarized in Table 1. Patients with PSA and the HCs exhibited no significant differences in age or sex, but they showed significant differences in Edu (p ≤ 0.001).
Spatial map and temporal properties of each IC
Varied interaction effect of temporal properties
According to the results of the two-way crossed ANOVA test, there was a significant effect on the degree of synchronization, including increased variability of SN [F = 2.23, positive false discovery rate (pFDR) = 0.017] (Figure 3B) and lower frequency of CN [F = −2.93, pFDR = 0.006] (Figure 3D; Supplementary Table 1).
Spatial map for the comparison of the PSA and HC groups
According to the results of the group main effect, different contributions from certain brain regions in the nine networks were detected. We found that the PSA group showed activated clusters in the left putamen, left postcentral gyrus (PostCG), and left angular gyrus (AG) in the SN compared to the HC group. No significant differences were observed in CN. Results are summarized in Figure 3C and Table 2.
Static FC during EA stimulation
Interaction effect on static FC
We found an interaction effect on functional coupling between the right calcarine and right lingual gyrus [F(14, 27) = 3.16, pFDR = 0.043] (Figure 4A; Table 3). However, no significant differences were observed at the network level.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/956931/fneur-13-956931-HTML-r1/image_m/fneur-13-956931-g004.jpg)
Figure 4. Effect of electroacupuncture (EA) stimulation on static functional connectivity (FC). (A) Interaction effect on static FC. (B) EA stimulus-dependent changes in FC for the contrast of EA stimulation vs. rest at baseline in post stroke aphasia group. (C) EA stimulus-dependent changes in FC for the contrast of EA stimulation vs. rest at baseline in healthy control group.
Changes in FC between EA stimulation and rest in each group
We also separately identified the FCs between ROIs during EA stimulation in each group. EA stimulation reduced FC between the left precentral and right middle temporal gyrus [F(4,13) = −4.03, pFDR = 0.024] and right orbital superior frontal gyrus [F(4,11) = −3.35, pFDR = 0.048] in PSA (Figure 4B). Additionally, ROI-based gPPI analysis of the HC group revealed that EA stimulation reduced FC between the left postcentral and right orbital middle frontal regions [F(3,5) = −4.01, pFDR = 0.025], compared with the rest condition (Figure 4C). The results are presented in Table 3. Network-based gPPI analysis revealed no significant differences in either HC groups or PSA groups.
Dynamic FC during EA stimulation
Results of the interaction effects at the network level
The outlook of the 12 components is shown in Supplementary Figure S2. Patients with PSA had higher variability of occipital network LN and CN coupling (Figure 5A), with stronger connections between the LN and CN (F = 4.29, pFDR = 0.042, uncorrected p = 0.0053), DMN and CN (F = 4.28, pFDR = 0.043, uncorrected p = 0.005) (Figure 5B; Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/956931/fneur-13-956931-HTML-r1/image_m/fneur-13-956931-g005.jpg)
Figure 5. Effect of electroacupuncture (EA) stimulation on dynamic independent component analysis. (A) Changes in temporal properties varies of circuit1. (B) Spatial map of circuit1 for the contrast of condition main effect. (C) Correlation between frequency of circuit1 and clinical parameter (CRRCAE score). (D,E) Changes in temporal properties varies of circuit3. (F) Spatial map of circuit3 for the contrast of condition main effect. (G–K) Correlation between clinical parameter (CRRCAE score) and dynamic FC values of triangle inferior frontal and lingual. CRRCAE, Chinese Rehabilitation Research Center for Chinese Standard Aphasia Examination.
Results of the interaction effects at the ROI level
The outlook of the 12 components is shown in Supplementary Figure S3. Patients with PSA had higher frequencies (Figure 5D) and lower variability (Figure 5E) of circuit 3, which connects the occipital, limbic, and somatosensory systems, with the strongest connections between the triangle inferior frontal and lingual gyri (F = 5.57, pFDR = 0.026, uncorrected p = 0.001), supramarginal and posterior cingulum gyri (F = 5.42, pFDR = 0.03, uncorrected p = 0.001), and middle cingulum and postcentral gyri (F = 5.27, pFDR = 0.036, uncorrected p = 0.001) (Figure 5F; Table 4). The frequency of factor 3 was directly anticorrelated with the verbal order score (uncorrected p = 0.031) (Figure 5C), and the FC values between the triangle inferior frontal and lingual gyri were directly anticorrelated with the noun transcribing score (uncorrected p = 0.026), verb transcribing score (uncorrected p = 0.022), command describing score (uncorrected p = 0.011), action description score (uncorrected p = 0.005), and dictating noun score (uncorrected p = 0.026) (Figures 5G–K).
Discussion
This study aimed to calculate the differences in EA stimulus-dependent changes in brain connectivity patterns between patients with PSA and healthy individuals. The interaction between EA stimulation and disease was examined using static ICA, gPPI and Dyn-ICA, with age, sex, and Edu as covariates, similar to a previous study's methods (45). By using static ICA and dynamic ICA approaches, we discovered commonalities with previous studies of acupuncture neuroimaging mechanisms for PSA, as well as some prominent points/networks generated from newer methodologies. The present study is the first to identify a PSA-related neuroimaging mechanism of EA stimulation at HT5 and GB39 both at the ROI and the network level. Our research showed that EA stimulation dramatically changed the nodes of the regions/networks involved in the LN, SN, CN, occipital cortex, somatosensory regions, and cerebral limbic system. Our previous research on the effects of EA on PSA was expanded by these results.
Static ICA
Recently, it has been acknowledged that stroke pathology often affects large-scale functional network structures rather than lesions (46). Based on previous studies, the therapeutic effects of acupuncture extend to the entire brain area and even to the brain network (23). Such influences cannot be calculated using a univariate approach (47) but can be solved using an ICA approach. First, ICA can find functionally independent, potentially spatially overlapping functional networks without the need for prior information of the task paradigm (48). Second, it enables researchers to assess the effects of EA on wider functional networks in patients with PSA. Third, mean activity values within networks as well as frequency and variability distributions could be computed, and the impacts of EA on different portions of a particular network in individual patients contributed to an overall change in score, making this technique more resilient to the impact of varying structural lesions. In the present study, the mean activity levels of left AG and left PostCG in the SN were higher in patients with PSA than in HCs, and the temporal distributions showed significantly higher SN variability, same as the previous study (7). There is abundant evidence in the literature that self-generated speech activates the SN, and the SN had been provided to have a directly correlated with residual language performance in patients with PSA (49). The AG is a region associated with phonological deficits (50). The fact that only individuals with motor aphasia were included may help to explain this. Thus, our findings suggest that EA causes the SN network to transmit more information in motor aphasia.
Static FC analysis
Studies on resting-state connections in healthy persons assume that positive connectivity represents integration and coordination between different brain regions, whereas negative connectivity represents separated or conflicting systems (51). First, we computed the effect of the stimulation condition against the rest of the baseline conditions in each group. Second, we calculated the interactions between the groups and conditions. Our sample's modifications were primarily connected to the somatosensory, language, and occipital cortices. These regions were previously shown to be associated with verb learning and visual word processing (52), and the integration of the visual and sensorimotor systems sustains action naming (53). Our results showed that a crucial neuronal signature underlying EA stimulation processing in PSA may be a change in communication across key regions involved in vision and naming.
Dyn-ICA analysis
These theories assume that FC is “static,” while mounting data support the idea that FC is dynamic rather than static (54). Dyn-ICA can examine time-varying and dynamic extensions of component analysis, which can be accomplished by using a classic component analysis within a sliding window setting (55), and isolate the signal from the noise and boost the sensitivity to identify individual differences (56). Previous studies have used Dyn-ICA to observe the pathological characteristics of diseases (43), and we are the first to use Dyn-ICA to study the mechanism of the effect of acupuncture. Similarly, our studies have clearly shown that Dyn-ICA increases sensitivity to variations between individuals or conditions (56). We used Dyn-ICA to incorporate temporal information, such as frequency and variability, while maintaining the spatial identification of traditional FC at the ROI and network level and arrived at two primary conclusions. First, at the ROI level, PSA-related dynamic connection showed a higher frequency of occipital cortex-limbic system-sensorimotor system connections, and their variability decreased across different dynamic factors. The structure of the FC map showed enhanced connectivity among the occipital, limbic, and sensorimotor system cortices. The higher FC values between the triangle inferior frontal and lingual gyri in the dynamic condition main effect correlated with worse transcribing, describing, and dictating scores. The same results can be found in a previous static analysis of the effect of acupuncture with sham acupuncture as a control (26), but not in a dynamic framework, to the best of our knowledge. Second, perhaps more importantly, at the network level, our results showed that PSA-related dynamic connections between the CN, DMN, and LN are more frequent during EA stimulation than rest at baseline. We found that EA at HT5 and GB39 changes the synchronization between the cerebellum and cerebrum, which is consistent with findings of a previous study (57). The DMN is a collection of areas that are more active at rest than when performing a task (58). These regions are believed to play a role in internal or low-level attentional mechanisms (59). Considerable evidence has demonstrated the essential involvement of the CN in a variety of language functions (60), and CN maps to cerebral association networks have already been reported to correlate with cognitive function (61). Our findings suggest that an aspect of the PSA treatment effect of EA stimulation may be a time-dependent connection between these networks of language function correlations.
Limitations of the present study
This study has some limitations. First, the lesioned voxels can affect fMRI preprocessing and results. However, to the best of our knowledge, there is no accepted standard procedure for the treatment of stroke lesions. We did not mask the lesion, considering that manual mask procuring is time-consuming and subjective (62). The quality of the segmentation and normalization does not seem to be affected in this study, and by adopting the ICA method, some effects of stroke lesions may be offset (56). However, we must admit that there is still more to be done in the future to address this problem. Second, pain and dysfunction in daily life caused by the disease can lead to emotions such as anxiety and depression (63), which can also affect the results of fMRI. Therefore, these factors should be included as covariates in future studies. Third, only patients with motor aphasia were included in this study. Future research is necessary to examine the effects of EA in more types of aphasia.
Conclusions
In conclusion, EA at HT5 and GB39 may improve language function by modulating FC stability among the LN, SN, CN, occipital, somatosensory, and cerebral limbic system regions.
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 human participants were reviewed and approved by the Medical Research Ethics Committee of Dongzhimen Hospital (ECPJ-BDY-2015-04) and all subjects gave written consent to participate. The patients/participants provided their written informed consent to participate in this study.
Author contributions
YG and JC designed the study. TL, ZT, CL, XL, XH, and JX gathered and managed data. HZ and BZ provided guidance on experimental design and data analysis. MX was responsible for the manuscript's conception, methods, and writing. JC, QK, GK, and SL revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Dongzhimen Hospital (hospital/grants), and grants were supported by Special Project of International Cooperation of Traditional Chinese Medicine of State Administration of Traditional Chinese Medicine (project no. 0610-2140NF020630), Special Public Welfare Industry and Scientific Research from the State Administration of Traditional Chinese Medicine (project no. 201407001-9), and National Natural Science Foundation of China (project no. 81973790).
Acknowledgments
We are appreciative of the participants' steadfast dedication and enthusiasm for our research.
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.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.956931/full#supplementary-material
Abbreviations
PSA, post stroke aphasia; HC, healthy control; EA, electroacupuncture; PPI, Psychophysiological Interaction; static ICA, static independent component analysis; dyn, ICA, dynamic independent component analysis; DMN, default mode network; SMN, sensorimotor network; CRRCAE, Chinese Rehabilitation Research Center for Chinese Standard Aphasia test; BDAE, Boston Diagnostic Aphasia Examination; VN, visual network; SN, salience network; DAN, dorsal attention network; FPN, frontoparietal network; LN, language network; CN, cerebellar network; AG, angular gyrus; SMG, supramarginal gyrus; PostCG, postcentral gyrus; ANTs, Advanced Neuroimaging Tools.
References
1. Thiel A, Zumbansen A. The pathophysiology of post-stroke aphasia: a network approach. Restor Neurol Neurosci. (2016) 34:507–18. doi: 10.3233/RNN-150632
2. Sheppard SM, Sebastian R. Diagnosing and managing post-stroke aphasia. Expert Rev Neurother. (2021) 21:221–34. doi: 10.1080/14737175.2020.1855976
3. Zhu D, Chang J, Freeman S, Tan Z, Xiao J, Gao Y, et al. Changes of functional connectivity in the left frontoparietal network following aphasic stroke. Front Behav Neurosci. (2014) 8:167. doi: 10.3389/fnbeh.2014.00167
4. Barbieri E, Mack J, Chiappetta B, Europa E, Thompson CK. Recovery of offline and online sentence processing in aphasia: language and domain-general network neuroplasticity. Cortex. (2019) 120:394–418. doi: 10.1016/j.cortex.2019.06.015
5. Klingbeil J, Wawrzyniak M, Stockert A, Saur D. Resting-state functional connectivity: an emerging method for the study of language networks in post-stroke aphasia. Brain Cogn. (2019) 131:22–33. doi: 10.1016/j.bandc.2017.08.005
6. Wang C, Qin W, Zhang J, Tian T, Li Y, Meng L, et al. Altered functional organization within and between resting-state networks in chronic subcortical infarction. J Cereb Blood Flow Metab. (2014) 34:597–605. doi: 10.1038/jcbfm.2013.238
7. Zhang C, Xia Y, Feng T, Yu K, Zhang H, Sami MU, et al. Disrupted functional connectivity within and between resting-state networks in the subacute stage of post-stroke aphasia. Front Neurosci. (2021) 15:746264. doi: 10.3389/fnins.2021.746264
8. Yang M, Li J, Yao D, Chen H. Disrupted intrinsic local synchronization in poststroke aphasia. Medicine. (2016) 95:e3101. doi: 10.1097/MD.0000000000003101
9. Liang FR, Wu X. [The developmental status and prospect of the science of acupuncture and moxibustion abroad]. Zhongguo Zhen Jiu. (2006) 26:79–82 (Chinese).
10. Hao JJ, Mittelman M. Acupuncture: past, present, and future. Global Adv Health Med. (2014) 3:6–8. doi: 10.7453/gahmj.2014.042
11. Wen J, Chen X, Yang Y, Liu J, Li E, Liu J, et al. Acupuncture medical therapy and its underlying mechanisms: a systematic review. Am J Chin Med. (2021) 49:1–23. doi: 10.1142/S0192415X21500014
12. Xiao LY, Wang XR, Yang Y, Yang JW, Cao Y, Ma SM, et al. Applications of acupuncture therapy in modulating plasticity of central nervous system. Neuromodulation. (2018) 21:762–76. doi: 10.1111/ner.12724
13. Zhang B, Han Y, Huang X, Liu Z, Li S, Chang J, et al. Acupuncture is effective in improving functional communication in post-stroke aphasia: a systematic review and meta-analysis of randomized controlled trials. Wien Klin Wochenschr. (2019) 131:221–32. doi: 10.1007/s00508-019-1478-5
14. Li G, Jack CR Jr., Yang ES. An fMRI study of somatosensory-implicated acupuncture points in stable somatosensory stroke patients. J Magn Reson Imaging. (2006) 24:1018–24. doi: 10.1002/jmri.20702
15. Li G, Yang ES. An fMRI study of acupuncture-induced brain activation of aphasia stroke patients. Compl Ther Med. (2011) 19 Suppl 1:S49–59. doi: 10.1016/j.ctim.2010.11.004
16. Chavez LM, Huang SS, MacDonald I, Lin JG, Lee YC, Chen YH. Mechanisms of acupuncture therapy in ischemic stroke rehabilitation: A literature review of basic studies. Int J Mol Sci. (2017) 18:2270. doi: 10.3390/ijms18112270
17. Jiang Y, Wang H, Liu Z, Dong Y, Dong Y, Xiang X, et al. Manipulation of and sustained effects on the human brain induced by different modalities of acupuncture: an fMRI study. PLoS ONE. (2013) 8:e66815. doi: 10.1371/journal.pone.0066815
18. Sun L, Chen YY, Fang JL, Hong Y, Wang Y, Xu K, et al. [Correlation between blood oxygen level dependent fMRI signal and GABA content in anterior cingulate cortex after acupuncture of Hegu (LI4)]. Zhen Ci Yan Jiu. (2019) 44:878–83 (Chinese). doi: 10.13702/j.1000-0607.190597
19. Jin L, Sun J, Xu Z, Yang X, Liu P, Qin W. Intersubject synchronisation analysis of brain activity associated with the instant effects of acupuncture: an fMRI study. Acupunct Med. (2018) 36:14–20. doi: 10.1136/acupmed-2016-011327
20. Zhang W, Lang S, Zheng Y, Qin X, Chen H, You Y, et al. The effects of transcranial direct current stimulation versus electroacupuncture on working memory in healthy subjects. J Altern Compl Med. (2019) 25:637–42. doi: 10.1089/acm.2018.0532
21. Liu S, Li M, Tang W, Wang G, Lv Y. An fMRI study of the effects on normal language areas when acupuncturing the Tongli (HT(5)) and Xuanzhong (GB(39)) acupoints. J Int Med Res. (2017) 45:1961–75. doi: 10.1177/0300060517720344
22. Li LL, Liu XW, Wu F, Tong DC, Ye LP, Tao HX, et al. Electroacupuncture stimulation of language-implicated acupoint Tongli (HT 5) in healthy subjects: an fMRI evaluation study. Chin J Integr Med. (2018) 24:822–9. doi: 10.1007/s11655-017-2924-8
23. Xiao J, Zhang H, Chang JL, Zhou L, Tan ZJ, Zhong HZ, et al. Effects of electro-acupuncture at Tongli (HT 5) and Xuanzhong (GB 39) acupoints from functional magnetic resonance imaging evidence. Chin J Integr Med. (2016) 22:846–54. doi: 10.1007/s11655-015-1971-2
24. Zhang J, Song L, Xu L, Fan Y, Wang T, Tian W, et al. Knowledge domain and emerging trends in ferroptosis research: a bibliometric and knowledge-map analysis. Front Oncol. (2021) 11:686726. doi: 10.3389/fonc.2021.686726
25. Han X, Jin H, Li K, Ning Y, Jiang L, Chen P, et al. Acupuncture modulates disrupted whole-brain network after ischemic stroke: evidence based on graph theory analysis. Neural Plast. (2020) 2020:8838498. doi: 10.1155/2020/8838498
26. Cai RL, Shen GM, Wang H, Guan YY. Brain functional connectivity network studies of acupuncture: a systematic review on resting-state fMRI. J Integr Med. (2018) 16:26–33. doi: 10.1016/j.joim.2017.12.002
27. Fan R, Gao Y, Zhang H, Xin X, Sang F, Tan Z, et al. Lesion distribution and early changes of right hemisphere in Chinese patients with post-stroke aphasia. Front Aging Neurosci. (2021) 13:632217. doi: 10.3389/fnagi.2021.632217
28. Chang J, Zhang H, Tan Z, Xiao J, Li S, Gao Y. Effect of electroacupuncture in patients with post-stroke motor aphasia : Neurolinguistic and neuroimaging characteristics. Wien Klin Wochenschr. (2017) 129:102–9. doi: 10.1007/s00508-016-1070-1
29. Maeda Y, Kim H, Kettner N, Kim J, Cina S, Malatesta C, et al. Rewiring the primary somatosensory cortex in carpal tunnel syndrome with acupuncture. Brain. (2017) 140:914–27. doi: 10.1093/brain/awx015
30. Gao Y, Lin Z, Tao J, Yang S, Chen R, Jiang C, et al. Evidence of timing effects on acupuncture: a functional magnetic resonance imaging study. Exp Ther Med. (2015) 9:59–64. doi: 10.3892/etm.2014.2056
31. Song XW, Dong ZY, Long XY Li SF, Zuo XN, Zhu CZ, et al. REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS ONE. (2011) 6:e25031. doi: 10.1371/journal.pone.0025031
32. Whitfield-Gabrieli S, Nieto-Castanon A. Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. (2012) 2:125–41. doi: 10.1089/brain.2012.0073
33. Calhoun VD, Adali T, Pearlson GD, Pekar JJ, A. method for making group inferences from functional MRI data using independent component analysis. Hum Brain Mapp. (2001) 14:140–51. doi: 10.1002/hbm.1048
34. Vettore M, De Marco M, Pallucca C, Bendini M, Gallucci M, Venneri A. White-matter hyperintensity load and differences in resting-state network connectivity based on mild cognitive impairment subtype. Front Aging Neurosci. (2021) 13:737359. doi: 10.3389/fnagi.2021.737359
35. Tordjman M, Madelin G, Gupta PK, Cordova C, Kurz SC, Orringer D, et al. Functional connectivity of the default mode, dorsal attention and fronto-parietal executive control networks in glial tumor patients. J Neurooncol. (2021) 152:347–55. doi: 10.1007/s11060-021-03706-w
36. Ferré P, Jarret J, Brambati S, Bellec P, Joanette Y. Functional connectivity of successful picture-naming: age-specific organization and the effect of engaging in stimulating activities. Front Aging Neurosci. (2020) 12:535770. doi: 10.3389/fnagi.2020.535770
37. Alakörkkö T, Saarimäki H, Glerean E, Saramäki J, Korhonen O. Effects of spatial smoothing on functional brain networks. Eur J Neurosci. (2017) 46:2471–80. doi: 10.1111/ejn.13717
38. Mclaren DG, Ries ML, Xu G, Johnson SC, A. generalized form of context-dependent psychophysiological interactions (gPPI): a comparison to standard approaches. Neuroimage. (2012) 61:1277–86. doi: 10.1016/j.neuroimage.2012.03.068
39. Lee Masson H, Op De Beeck H, Boets B. Reduced task-dependent modulation of functional network architecture for positive versus negative affective touch processing in autism spectrum disorders. Neuroimage. (2020) 219:117009. doi: 10.1016/j.neuroimage.2020.117009
40. Ghahremani M, Yoo J, Chung SJ, Yoo K, Ye JC, Jeong YJ. Alteration in the local and global functional connectivity of resting state networks in Parkinson's disease. J Mov Disord. (2018) 11:13. doi: 10.14802/jmd.17061
41. Calhoun VD, Miller R, Pearlson G, Adali 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
42. Grami F, De Marco G, Bodranghien F, Manto M, Habas CJ. Cerebellar transcranial direct current stimulation reconfigures brain networks involved in motor execution and mental imagery. Cerebellum. (2022) 21:665–80. doi: 10.1007/s12311-021-01322-y
43. Aracil-Bolaños I, Martínez-Horta S, González-De-Echávarri JM, Sampedro F, Pérez-Pérez J, Horta-Barba A, et al. Structure and dynamics of large-scale cognitive networks in Huntington's Disease. Mov Disord. (2022) 37:343–53. doi: 10.1002/mds.28839
44. Guerithault N, McClure SM, Ojinnaka CO, Braden BB, Bruening M. Resting-state functional connectivity differences in college students with and without food insecurity. Nutrients. (2022) 14:2064. doi: 10.3390/nu14102064
45. Laforce R Jr., Tosun D, Ghosh P, Lehmann M, Madison CM, Weiner MW, et al. Parallel ICA of FDG-PET and PiB-PET in three conditions with underlying Alzheimer's pathology. NeuroImage Clin. (2014) 4:508–16. doi: 10.1016/j.nicl.2014.03.005
46. Salvalaggio A, De Filippo De Grazia M, Zorzi M, Thiebaut De Schotten M, Corbetta M. Post-stroke deficit prediction from lesion and indirect structural and functional disconnection. Brain. (2020) 143:2173–88. doi: 10.1093/brain/awaa156
47. Darkow R, Martin A, Würtz A, Flöel A, Meinzer M. Transcranial direct current stimulation effects on neural processing in post-stroke aphasia. Hum Brain Mapp. (2017) 38:1518–31. doi: 10.1002/hbm.23469
48. Smith SM, Miller KL, Moeller S, Xu J, Auerbach EJ, Woolrich MW, et al. Temporally-independent functional modes of spontaneous brain activity. Proc Natl Acad Sci USA. (2012) 109:3131–6. doi: 10.1073/pnas.1121329109
49. Brownsett SL, Warren JE, Geranmayeh F, Woodhead Z, Leech R, Wise RJ. Cognitive control and its impact on recovery from aphasic stroke. Brain. (2014) 137:242–54. doi: 10.1093/brain/awt289
50. Na Y, Jung J, Tench CR, Auer DP, Pyun SB. Language systems from lesion-symptom mapping in aphasia: a meta-analysis of voxel-based lesion mapping studies. NeuroImage Clin. (2022) 35:103038. doi: 10.1016/j.nicl.2022.103038
51. Fair DA, Dosenbach NU, Church JA, Cohen AL, Brahmbhatt S, Miezin FM, et al. Development of distinct control networks through segregation and integration. Proc Natl Acad Sci USA. (2007) 104:13507–12. doi: 10.1073/pnas.0705843104
52. Fiori V, Kunz L, Kuhnke P, Marangolo P, Hartwigsen G. Transcranial direct current stimulation (tDCS) facilitates verb learning by altering effective connectivity in the healthy brain. Neuroimage. (2018) 181:550–9. doi: 10.1016/j.neuroimage.2018.07.040
53. Durand E, Masson-Trottier M, Sontheimer A, Ansaldo AI. Increased links between language and motor areas: a proof-of-concept study on resting-state functional connectivity following personalized observation, execution and mental imagery therapy in chronic aphasia. Brain Cogn. (2021) 148:105659. doi: 10.1016/j.bandc.2020.105659
54. Fu Z, Iraji A, Sui J, Calhoun VD. Whole-brain functional network connectivity abnormalities in affective and non-affective early phase psychosis. Front Neurosci. (2021) 15:682110. doi: 10.3389/fnins.2021.682110
55. Leonardi N, Richiardi J, Gschwind M, Simioni S, Annoni JM, Schluep M, et al. Principal components of functional connectivity: a new approach to study dynamic brain connectivity during rest. Neuroimage. (2013) 83:937–50. doi: 10.1016/j.neuroimage.2013.07.019
56. Zhang H, Tao Y, Xu H, Zou S, Deng F, Huang L, et al. Associations between childhood chronic stress and dynamic functional connectivity in drug-naïve, first-episode adolescent MDD. J Affect Disord. (2022) 299:85–92. doi: 10.1016/j.jad.2021.11.050
57. Chen J, Wang J, Huang Y, Lai X, Tang C, Yang J, et al. Modulatory effect of acupuncture at Waiguan (TE5) on the functional connectivity of the central nervous system of patients with ischemic stroke in the left basal ganglia. PLoS ONE. (2014) 9:e96777. doi: 10.1371/journal.pone.0096777
58. Van De Ven V, Esposito F, Christoffels IK. Neural network of speech monitoring overlaps with overt speech production and comprehension networks: a sequential spatial and temporal ICA study. Neuroimage. (2009) 47:1982–91. doi: 10.1016/j.neuroimage.2009.05.057
59. Buckner R, Andrews-Hanna J, Schacter DJ. The brain's default network: anatomy, function, and relevance to disease. Ann N Y Acad Sci. (2008) 1124:1–38. doi: 10.1196/annals.1440.011
60. Zheng K, Chen M, Shen Y, Xu X, Gao F, Huang G, et al. Cerebellar continuous theta burst stimulation for aphasia rehabilitation: study protocol for a randomized controlled trial. Front Aging Neurosci. (2022) 14:909733. doi: 10.3389/fnagi.2022.909733
61. Buckner RL. The cerebellum and cognitive function: 25 years of insight from anatomy and neuroimaging. Neuron. (2013) 80:807–15. doi: 10.1016/j.neuron.2013.10.044
62. Brihmat N, Boulanouar K, Darmana R, Biganzoli A, Gasq D, Castel-Lacanal E, et al. Controlling for lesions, kinematics and physiological noise: impact on fMRI results of spastic post-stroke patients. MethodsX. (2020) 7:101056. doi: 10.1016/j.mex.2020.101056
Keywords: electroacupuncture, brain networks, post stroke aphasia, functional connectivity, psychophysiological interaction analysis, independent component analysis
Citation: Xu M, Gao Y, Zhang H, Zhang B, Lyu T, Tan Z, Li C, Li X, Huang X, Kong Q, Xiao J, Kranz GS, Li S and Chang J (2022) Modulations of static and dynamic functional connectivity among brain networks by electroacupuncture in post-stroke aphasia. Front. Neurol. 13:956931. doi: 10.3389/fneur.2022.956931
Received: 30 May 2022; Accepted: 10 October 2022;
Published: 01 December 2022.
Edited by:
Jiliang Fang, Guang'anmen Hospital, China Academy of Chinese Medical Sciences, ChinaReviewed by:
Jie Jia, Fudan University, ChinaHuangan Wu, Shanghai University of Traditional Chinese Medicine, China
Copyright © 2022 Xu, Gao, Zhang, Zhang, Lyu, Tan, Li, Li, Huang, Kong, Xiao, Kranz, Li and Chang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingling Chang, ear6979@163.com; Ying Gao, gaoying973@126.com