Continuous theta burst stimulation-induced suppression of the right fronto-thalamic-cerebellar circuit accompanies improvement in language performance in poststroke aphasia: A resting-state fMRI study

Introduction

Aphasia, one of the most common impairments following stroke, is usually caused by lesions affecting the network of cortical and subcortical structures. Poststroke aphasia (PSA), characterized by deficits in spontaneous speech, auditory comprehension, repetition, naming, etc., is associated with prolonged hospitalization, decreased quality of life, and increased mortality (Lazar and Boehme, 2017). Thus, the treatment of PSA needs to be further explored.

In addition to traditional speech and language therapy (SLT), repetitive transcranial magnetic stimulation (rTMS), one of the most promising noninvasive brain stimulation (NIBS) modalities, has emerged as an ideal method to improve the language function of patients with PSA (Ren et al., 2014; Kapoor, 2017; Gholami et al., 2022). rTMS can modulate the excitability and activity of targeted cortical regions, thereby promoting the reorganization of functional networks (Kapoor, 2017; Yoo et al., 2020; Jargow et al., 2021). Continuous theta burst stimulation (cTBS) is a novel rTMS method that reduces excitability. It has shorter stimulation times and better therapeutic effects than conventional paradigms (Bai et al., 2022; Chen et al., 2022). A recent study showed that cTBS over the right pars triangularis of individuals with chronic PSA enhanced word retrieval by facilitating phonological access (Harvey et al., 2019). Yoo et al. (2020) delivered neuronavigated cTBS over the left pars opercularis in healthy subjects, resulting in significantly increased phase synchronization of TMS-evoked potentials (TEPs) in language areas between the two hemispheres and within the left hemisphere. These results may have implications for the development of a cTBS protocol to assist language rehabilitation in patients with PSA.

Given interhemispheric competition, inhibiting the healthy hemisphere will promote the recovery of the lesioned hemisphere. Inhibitory low-frequency rTMS targeting the pars triangularis, an area in the right posterior inferior frontal gyrus (pIFG), has positive effects in treating aphasia (Ren et al., 2014; Kapoor, 2017; Hu et al., 2018). The right posterior superior temporal gyrus (pSTG), homotopic to Wernicke’s area in the left hemisphere, is a critical area for semantic processing and is thus another optimal target for rTMS stimulation. However, few studies have targeted the pSTG in aphasic patients. A randomized sham-controlled study found significant improvement in auditory comprehension and repetition in patients with subacute global aphasia following 3 weeks of low-frequency rTMS over the right pSTG coupled with SLT (Ren et al., 2019). Kawamura et al. administered low-frequency rTMS over the right pSTG to a patient with severe aphasia due to injury of the left IFG and demonstrated subsequent improvement in language perception (Kawamura et al., 2019). Moreover, significant activations of the left and right temporal gyrus, left posterior cingulate gyrus, and left thalamus during a language task following rTMS over the pSTG were observed by functional magnetic resonance imaging (fMRI; Kawamura et al., 2019). A pilot study by Versace demonstrated that excitatory intermittent TBS (iTBS) over Wernicke’s area led to transient facilitation of auditory comprehension in patients with chronic fluent aphasia (Versace et al., 2019). However, these studies have several limitations, including the lack of exploration of underlying mechanisms and insufficient sample size.

However, research is lacking on whether cTBS of the contralesional pSTG improves language outcomes in patients with PSA, and the underlying mechanisms are still unknown. Resting-state fMRI (rs-fMRI) has proven to be a feasible tool for the exploration of networkwide lesion effects with the potential to provide insight into network reconfiguration following stroke and is increasingly being utilized to investigate PSA (Ilkhani et al., 2018; Allendorfer et al., 2021a,b). Therefore, the present study was conducted to determine whether inhibitory cTBS over the right pSTG (coupled with SLT) benefits patients with subacute aphasia and to explore the remodeling mechanism of neural plasticity. Specifically, the current study aimed to explore changes in intrinsic neural activity by examining the outcomes of rs-fMRI, including whole-brain analysis of fractional amplitude of low-frequency fluctuations (fALFF) and functional connectivity (FC).

Materials and methods

Participants

A total of 34 participants with aphasia were recruited from Wuxi Mental Health Center (Wuxi Central Rehabilitation Hospital). Details about the study design and data collection are shown in Figure 1. The inclusion criteria were as follows: (1) provided written informed consent; (2) aged 45 to 75 years; (3) were right-handed (as assessed by the Edinburgh Handedness Inventory; Oldfield, 1971), native Chinese speakers without professional vocal or instrumental training; (4) experienced their first middle cerebral artery (MCA) stroke on the left side, with an MRI-identified lesion site; (5) were 4–12 weeks poststroke; (6) had aphasia confirmed by the Chinese version of Western Aphasia Battery (WAB) scores; and (7) had at least elementary education and normal or corrected-to-normal vision and hearing. The exclusion criteria were as follows: (1) history of epilepsy, drug abuse, or neuropsychiatric diseases; (2) speech impairment or severe dysarthria before the onset of this disease; (3) contraindications to TMS or fMRI (e.g., skull defects or skin lesions at the site of stimulation, intracranial implants, pacemakers, or implanted drug pumps); (4) other neurological diseases; and (5) history of neurosurgical treatment.

FIGURE 1

Figure 1. CONSORT diagram of patient flow through the study.

The present randomized study was approved by the Ethics Committee of Wuxi Mental Health Center (Wuxi Central Rehabilitation Hospital) (No. WXMHCIR2021LLky054) and registered in the Chinese Clinical Trial Registry (No. ChiCTR210052962).

Procedure

A prospective, randomized, and sham-controlled study was performed. The random allocation sequence was generated by a completely randomized digital table. A blinded investigator randomly assigned patients to the cTBS group and the sham group using numbers in opaque envelopes. The patients in both groups (the cTBS group: n = 15; the sham group: n = 14) underwent 15 sessions (5 times a week, from Monday to Friday) of real or sham cTBS coupled with 30 min of SLT over 3 weeks. The rehabilitation programs focused on specific training that stimulated various aspects of the language system, such as semantic, syntactic, phonological and motor speech performance.

Language assessment was obtained before and after 3 weeks of real or sham cTBS treatment with the WAB, which consists of the AQ score and subscales of performance on spontaneous speech, auditory verbal comprehension, repetition, and naming. The AQ score reflects the severity of aphasia and can be used to assess the improvement in or worsening of aphasia symptoms. The highest AQ score is 100; the normal range is 98.4–99.6. An AQ score < 93.8 indicates aphasia.

Image acquisition

MRI scanning was performed with a 3.0 T MAGNETOM Skyra scanner (Siemens, Germany) the day before and after the intervention. A total of 240 brain volumes were obtained using an echo-planar imaging (EPI) pulse fMRI sequence with the following parameters: repetition time (TR) = 2,000 ms; echo time (TE) = 30 ms; 35 slices, field of view (FOV) = 224 mm × 224 mm; slice thickness = 3.5 mm; layer spacing = 0.7 mm; flip angle (FA) = 90°; acquisition matrix = 64 × 64; and voxel size = 3.5 mm × 3.5 mm × 3.5 mm. The parameters of the high-resolution three-dimensional (3D) structural T1-weighted images were as follows: 192 sagittal slices; slice thickness = 1 mm; TR = 6.6 ms; TE = 3.1 ms; and FA = 12°. Each subject was instructed to remain still, close their eyes, and not think of anything in particular before scanning.

Lesion analysis

The lesions were manually segmented into individual structural images obtained from the T1-weighted sequence using ITK-SNAP tools (v. 3.2) (Yushkevich et al., 2016) and registered to Montreal Neurological Institute (MNI) space using SMP12. The sum of all the individual lesions was used to construct a voxelwise map of the lesions.

Preprocessing of resting-state functional MRI (rs-fMRI) data

rs-fMRI data were preprocessed on the MATLAB platform (MathWorks, Natick, MA, United States) with DPABI software (v. 6.1) (Yan et al., 2016). The main steps were as follows: (1) conversion of the raw data in DICOM format to NIFTI format using MRIConvert software (v. 2.1.0, Lewis Center, Eugene, OR, United States); (2) removal of the first 10 slices and slice-timing correction of the remaining 230 slices; (3) correction of head motion using Friston 24 parameters; (4) spatial normalization of the functional images to MNI space and resampling to a voxel size of 3 mm × 3 mm × 3 mm; (5) application of a bandpass filter (0.01–0.08 Hz) to each voxel to reduce the influence of low-frequency fluctuations and high-frequency noise (Lv et al., 2018); and (6) smoothing of the functional data with a 6-mm full-width at half-maximum (FWHM) Gaussian kernel. An enantiomorphic approach was performed to replace the lesioned brain tissue with the contralateral mirrored scans during spatial normalization (Nachev et al., 2008).

Calculation of the fractional amplitude of low-frequency fluctuations (fALFF)

fALFF indicates the amplitude of spontaneous fluctuation in blood oxygen level-dependent (BOLD) signal intensity in a specific brain area in fMRI, reflecting the level of spontaneous activity of a voxel in the resting state (Zou et al., 2008). The time series data were converted into the frequency spectrum using a fast Fourier transform (FFT), and the power spectrum was obtained. In each voxel, the square root of the signal over 0.01–0.08 Hz was calculated. The fALFF values of each voxel were normalized through transformation to z scores.

Seed-based functional connectivity analyses

Seed-based FC analyses were used to investigate functional changes in the right pSTG and right pars triangularis compared to all the other voxels in the brain. Two 10-mm spherical regions of interest (ROIs) were created based on the Harvard-Oxford Cortical Atlas in the DPABI toolbox (Whitfield-Gabrieli and Nieto-Castanon, 2012; Yan et al., 2016). The MNI coordinates of the right pSTG and right pars triangularis were x = 61, y = −24, z = 1 and x = 48, y = 30, z = 6, respectively (Amad et al., 2017). The time courses of all voxels within a 5-mm sphere around the center coordinate were averaged. Connectivity maps were calculated by calculating the Fisher z-transformed correlation coefficient between the average time course and the time course of each voxel in the whole-brain mask for each subject.

cTBS protocol

T1-weighted structural imaging data were acquired before cTBS intervention for precise targeting of regions in individuals. The individual stimulation sites around the right pSTG (MNI coordinates x = 61, y = −24, z = 1) were identified in each patient. Individual structural images were loaded into the visor2^™ TMS neuronavigation system (ANT Neuro, Netherlands) to enable fMRI-guided cTBS intervention. The cTBS intervention was conducted with a Magneiro^® (Vishee Medical Technology Co., Ltd., Nanjing, China) with a 70-mm figure-8-shaped stimulation coil. Before intervention, a single-pulse TMS was applied over the right primary motor cortex to determine the active motor threshold (AMT). The AMT was defined as the lowest stimulus intensity that elicited a motor evoked potential (MEP) >200 μV in at least 5–10 consecutive trials, with a 10% maximal voluntary contraction of the first interosseous muscle of the left hand (Koch et al., 2020). Each subject received real or sham cTBS over the right pSTG at 80% of the AMT (Nee and D’Esposito, 2017). The standard cTBS protocol consisted of a burst of three 50-Hz pulses, presented repeatedly at a frequency of 5 Hz, resulting in a total of 600 pulses applied over 40 s. The coil was placed tangentially to the skull with the handle pointing backward and laterally angled at 45° from the sagittal axis during real stimulation. In contrast, the stimulation coil was rotated 90° to place one side of the coil on the target during sham stimulation.

Statistical analysis

Statistical analysis was performed using SPSS for Windows, version 24.0 (SPSS, Inc., Chicago, IL, United States). The sample size estimation was determined according to a previous randomized controlled study (Thiel et al., 2013). One-way analyses of variance (ANOVAs) were used for continuous data, and the χ² test was used for categorical data in group comparisons of baseline demographic characteristics. Between-group differences in WAB score changes were assessed using independent-sample t tests and a 2-factor repeated-measures ANOVA.

The overall pre-and posttreatment maps of fALFF and FC values were performed using a 2-factor mixed-effects ANOVA, with time (pre-vs. post-treatment) as a within-subjects factor and treatment (cTBS vs. sham) as a between-subjects factor. Framewise displacement (FD) was included as a covariate. Multiple comparisons across voxels were corrected by a Gaussian random field (GRF) (voxelwise p < 0.001, clusterwise p < 0.05, two-tailed). Linear regression analysis was used to determine the relationship between changes in WAB scores and fALFF values.

Results

Participant characteristics

Thirty-four eligible subjects participated in the study after application of the inclusion and exclusion criteria. Four patients (two from the sham group and two from the cTBS group) withdrew from the trial after providing informed consent and after being allocated to the groups. One patient (from the cTBS group) dropped out because of complications. Thus, a total of 29 participants completed the study (cTBS group: n = 15; sham group: n = 14) (Figure 1).

The baseline demographic characteristics of all 29 subjects were shown in Table 1. No significant differences between the two groups were observed at baseline in terms of age (p = 0.47), sex (p = 0.76), time since stroke (p = 0.59), or WAB scores (p > 0.05).

TABLE 1

Table 1. Demographic and clinical characteristics of the participants.

Language measures

A repeated-measures ANOVA revealed a significant group × time interaction effect on WAB-AQ scores [F(2,27) = 7.687, p = 0.010, η² = 0.222] (Table 2). There was a significant main effect of time [F(2,27) = 79.172, p < 0.001, η² = 0.746], but no significant main effect of treatment [F(2,27) = 3.559, p = 0.070, η² = 0.116]. Independent-samples t tests showed that the change in WAB-AQ scores was significantly higher in the cTBS group than in the sham group (t = 2.773, 95% CI = 1.966 to 12.162, p = 0.010) (Figure 2).

TABLE 2

Table 2. Statistical data of WAB score changes with pre-and post-treatment between the cTBS group and the sham group.

Regarding scores on the WAB subscales, there was a significant group × time interaction effect on auditory comprehension [F(2,27) = 5.870, p = 0.022, η² = 0.179] and repetition scores [F(2,27) = 4.927, p = 0.035, η² = 0.154] (Table 2). There was no significant interaction effect on spontaneous speech or naming scores (p > 0.10). The changes in auditory comprehension scores (t = 2.423, 95% CI = 0.159 to 1.926, p = 0.022) and repetition scores (t = 2.220, 95% CI = 0.092 to 2.345, p = 0.035) were significantly higher in the cTBS group than in the sham group (Figure 2).

FIGURE 2

Figure 2. The change in the aphasia quotient (AQ) score of the Western Aphasia Battery (WAB) was significantly greater in the continuous theta burst stimulation (cTBS) group than in the sham group (A). cTBS significantly improved comprehension (COMP) and repetition (REP) performance but not spontaneous speech (SPO) or naming (NAM) performance (B). *Significant at p < 0.05. ** Significant at p ≤ 0.01.

fALFF changes after cTBS

Figure 3 presents the lesion overlap map and the distribution of lesions for all patients.

FIGURE 3

Figure 3. Lesion overlap map for 29 patients with PSA to illustrate the distribution of lesions. The color scale in the spectrograms represents the frequency of injury in a given location.

Figure 4 and Table 3 indicate the clusters with significant between-group differences in fALFF before and after cTBS intervention. A 2-factor mixed-effects ANOVA showed that the great significant difference in group and time interaction occurred in brain regions including the pars triangularis of the inferior frontal gyrus (IFG), right middle frontal gyrus, right thalamus, left orbitofrontal cortex and left cerebellar Crus I (corrected by GRF, voxelwise p < 0.001, clusterwise p < 0.05, two-tailed). Specifically, the intrinsic activity indexed by the fALFF values in the right pars triangularis of IFG, right middle frontal gyrus, right thalamus, and left cerebellar Crus I decreased significantly from pre-to post-treatment in the cTBS group compared to those in the sham group. In contrast, the fALFF value of the left orbitofrontal cortex increased significantly from pre-to post-treatment in the cTBS group compared to the sham group.

FIGURE 4

Figure 4. (A) F map of differences in spontaneous activity pre-and post-treatment between the cTBS group and sham group. Brilliant colors (red to bright yellow) indicate a greater significant difference over pre-and post-treatment between the cTBS and sham groups (group × time interaction). (B–F) Fractional amplitude of low-frequency fluctuations (fALFF) decreased significantly after cTBS in the right pars triangularis, left cerebellar crus I, right middle frontal gyrus, and right thalamus and increased significantly after cTBS in the left orbitofrontal cortex. We used a 2-factor mixed-effects ANOVA. ** indicates p < 0.001 in (B–F). Bars show means, and error bars represent the standard error of the mean (SEM). (G) The regression analysis showed a significant linear correlation between changes in fALFF in the right pars triangularis and changes in the WAB comprehension score. IFG, inferior frontal gyrus; WAB, Western Aphasia Battery.

TABLE 3

Table 3. Brain regions with significantly different values of fALFF in the intragroup (pre-and post-treatment) and intergroup (cTBS and sham groups) comparisons with group × time interactions (GRF, voxelwise p < 0.001, clusterwise p < 0.05, two-tailed).

A linear regression analysis was performed to test whether intrinsic regional activity, represented by changes in fALFF, correlated with clinical recovery. The regression analysis revealed that the change in WAB comprehension scores were significantly associated with the change in the fALFF value of the right IFG pars triangularis in both groups (R² = 0.651, p = 0.009).

Functional connectivity

We defined the right IFG pars triangularis as the seed region based on the results of the linear regression analysis. A 2-factor mixed-effects ANOVA showed that the cTBS group demonstrated increased FC between the right IFG pars triangularis and left paracingulate gyrus but reduced FC of these regions in the sham group (Figure 5). The right pSTG, the stimulation site, was set as the other seed region. The FC of the right pSTG with the right angular gyrus and posterior cingulate gyrus was significantly increased after treatment in the cTBS group but was reduced after treatment in the sham group (Figure 6).

FIGURE 5

Figure 5. cTBS-induced changes in intrinsic functional connectivity (FC) between the right pars triangularis and the left paracingulate gyrus (A). The cTBS group demonstrated increased FC between the right pars triangularis and left paracingulate gyrus after treatment, but FC between these regions was reduced after treatment in the sham group. Bars show means, and error bars represent the standard error of the mean (SEM). We used a 2-factor mixed-effects ANOVA. ** indicates p < 0.001 in (B).

FIGURE 6

Figure 6. cTBS-induced changes in intrinsic functional connectivity (FC) of the right posterior superior temporal gyrus (pSTG) with the right angular gyrus and the posterior cingulate gyrus (A). The cTBS group demonstrated increased FC of the right pSTG with the right angular gyrus and posterior cingulate gyrus after treatment; however, in the sham group, the FC between these areas was reduced after treatment. Bars show means, and error bars represent the standard error of the mean (SEM). We used a 2-factor mixed-effects ANOVA. ** indicates p < 0.001 in (B,C).

Discussion

Our prior work documented the effectiveness of low-frequency rTMS (LF-rTMS) over the right pSTG for improving language performance in stroke patients with aphasia (Ren et al., 2019). The present study revealed that cTBS over the right pSTG (homotopic to Wernicke’s area) coupled with SLT improved language performance, especially auditory comprehension and repetition, in patients with PSA, as demonstrated by the increase in WAB scores on the AQ, auditory comprehension and repetition. These findings expand potential stimulation targets for effective intervention beyond anterior language areas, in which most studies targeted areas homotopic to Broca’s area, to posterior temporal language areas in the right hemisphere. Notably, the observed improvements in our study were paralleled by suppression of spontaneous neural activity in the right fronto-thalamic-cerebellar circuit and an increase in FC of the right temporoparietal lobe, as measured by rs-fMRI. Therefore, this study may provide a fruitful target for improving the language function of auditory comprehension as well as valuable insight into how cTBS reorganizes language networks in aphasic patients.

We found that the spontaneous neural activity of the right pars triangularis was significantly suppressed after patients received inhibitory cTBS over the right pSTG; this decrease was closely related to the improvement in auditory comprehension. Therefore, cTBS of the right pSTG may downregulate the activity of the right pars triangularis through structural or functional connections between the two major nodes, such as the arcuate fasciculus. The right pars triangularis, the homotopic area to Broca’s area frequently used as a stimulation target, was inhibited in most of the previous studies addressing interhemispheric imbalance by decreasing activity in the right hemisphere or increasing activity in the left hemisphere (Thiel et al., 2013; Heikkinen et al., 2018; Lee et al., 2022; Zumbansen et al., 2022). Our fALFF results provided compelling evidence for cTBS-induced reduction in spontaneous neural activity not only in the right frontal language regions but also in subcortical regions, including the right thalamus and the left cerebellum, in aphasic patients compared with those who received sham treatment. This result is consistent with that of a positron emission tomography (PET) study, which showed hypometabolism in the right hemisphere and left cerebellum after stimulating Wernicke’s area in a patient with crossed aphasia in dextral areas (Lu et al., 2014). Additionally, the previous clinical trial illustrated that clusters in the right thalamus exhibited suppressed fALFF in patients who received inhibitory LF-rTMS of the right pars triangularis (Lee et al., 2022), which was in accordance with our findings.

We noticed that a similar synchronized decline in fALFF was present in the frontal and subcortical areas belonging to the contralesional cerebro-cerebellar circuit, including the right IFG pars triangularis, right middle frontal gyrus, right thalamus and left cerebellar hemisphere. This finding demonstrates that the inhibitory stimulation of cTBS over the right pSTG downregulated the intrinsic neural activity in the right frontal-thalamic-cerebellar circuit compared with that of the sham group, which might be due to the major neuroplastic alteration of the contralesional cerebro-cerebellar system not only in response to stroke and left cerebral hemispheric glioma involving the language network (Wessel and Hummel, 2018; Zhang et al., 2018) but also following the inhibitory stimulation of cTBS over the contralesional hemisphere cortex. The cerebro-cerebellar circuit, which transports information between the contralateral cerebellum and the ipsilateral cerebral cortex via the thalamus, is properly organized in a crossing pattern (Habas et al., 2009; Stoodley and Schmahmann, 2010; Buckner, 2013; De Smet et al., 2013). A number of anatomical, clinical, and functional imaging studies have confirmed that the left cerebro-thalamic-cerebellar circuit is involved in language processing with strong lateralization (Marien et al., 2001; Baillieux et al., 2010; Stoodley and Schmahmann, 2010; Mariën and Beaton, 2014; Schwartze and Kotz, 2016; Mariën and Borgatti, 2018). Additionally, a previous task-state fMRI study in which aphasic patients exhibited stronger lateralization to the right cerebellar hemisphere during a semantic decision/tone decision task after 3 weeks of iTBS over the ipsilesional hemisphere (Szaflarski et al., 2021). These findings highlight the potentially critical role of the cerebro-cerebellar circuit in the remodeling process following TBS intervention in the cerebral language network. Moreover, the spontaneous neural activity in the left orbitofrontal cortex significantly increased after 3 weeks of cTBS treatment in aphasic patients. The orbitofrontal cortex has long been recognized as being involved in affection and decision-making (Rolls and Grabenhorst, 2008). Its action during language processing remains unknown. However, the orbitofrontal cortex has been found to be coactivated with a network of prefrontal regions and the bilateral thalamus, which are involved in cognitive functions including language and memory (Zald et al., 2014). Accordingly, we can infer that the increased neural activity in the left orbitofrontal cortex may be viewed as the recruitment of intra-hemispheric language-network plasticity following the cTBS intervention.

The static FC of the right pars triangularis with the left paracingulate gyrus and of the right pSTG with the right angular gyrus and posterior cingulate gyrus was increased after treatment in the cTBS group but decreased after the sham intervention. These increased FC areas belong to the domain-general cognitive network, which implies that the domain-general cognitive network may be associated with recovery from aphasia (Brownsett et al., 2014). Cognitive deficits in aphasic patients are difficult to measure because of challenges with verbal communication (Arheix-Parras et al., 2021). The anterior paracingulate gyrus participates in speech processing through its connection to the prefrontal cortex, as confirmed by a previous fMRI study using word generation tasks (Crosson et al., 1999). Moreover, the angular gyrus may be included in the “language association areas,” as shown by the activation of the bilateral angular gyri (BA39) and the other language-related cortices following language tasks in a meta-analysis of fMRI data (Rosselli et al., 2015). The current study revealed that cTBS enhanced the involvement of the right angular cortex and the posterior cingulate gyrus in language processing, as reflected by the increased FC and improved clinical-assessment scores. This result is consistent with previous findings regarding structural and functional connectivity and supports the integrative role of the angular gyrus in language functions (Margulies and Petrides, 2013; Sroka et al., 2015; Ardila et al., 2016; Yourganov et al., 2016). A connectome-based analysis reported that the connectivity of the temporoparietal junction (TPJ), a multimodal region critical for language tasks, could predict repetition scores, as with, for example, the connectivity of the TPJ with the angular gyrus and pSTG (Yourganov et al., 2016).

However, this study also has some noteworthy limitations. We did not conduct subgroup analysis of aphasia subtypes because of the relatively small sample size. In future studies, the aphasia subtype and severity of language impairment should be considered during random assignment and treatment. Additionally, long-term linguistic improvement and neural reorganization could not be observed because outcomes were measured immediately after the intervention. Long-term follow-up research could supplement these findings and provide a more comprehensive understanding of brain remodeling.

Conclusion

To the best of our knowledge, this is the first randomized, sham-controlled clinical trial to investigate the effectiveness and underlying neural mechanism of cTBS over the right pSTG to treat subacute stroke patients with aphasia. Our data provide preliminary evidence that inhibitory cTBS of the right pSTG combined with SLT boosts the recovery of language performance, especially auditory comprehension and repetition. Moreover, the rs-fMRI data showed that the recovery of language function was associated with inhibition of the right fronto-thalamic-cerebellar circuit and increases in the neuroplasticity of the right temporoparietal lobe. Future studies with large samples are needed to elucidate the precise mechanisms underlying language recovery in the fronto-thalamic-cerebellar circuit. However, our results further optimize the TMS protocol and clarify the role of neural remodeling in the right hemisphere during aphasia recovery.

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 Ethics Committee of Wuxi Mental Health Center (Wuxi Central Rehabilitation Hospital; No. WXMHCIR2021LLky054). The patients/participants provided their written informed consent to participate in this study.

Author contributions

CR and GZ designed the study. CR, KZ, XX, and YJ drafted the protocol. KZ coordinated the trial. HF, FG, and LB were responsible for data acquisition and clinical evaluation. GH and BS administered the interventions. CR, KZ, and GZ revised the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by the Natural Science Foundation of Jiangsu Province (no. BK20201138), the Social Development Project of Jiangsu Province (no. BE2022700), the Wuxi Taihu Talent Project (no. WXTTP2020008), and the Top Talent Support Program for Young and Middle-Aged People of Wuxi Health Committee and General Project from Wuxi Health Commission (no. MS201911).

Acknowledgments

The authors thank the study participants and Hui Zheng from Shanghai Jiaotong University for his comprehensive assistance with data analysis.

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.

References

Allendorfer, J. B., Nenert, R., Nair, S., Vannest, J., and Szaflarski, J. P. (2021a). Functional magnetic resonance imaging of language following constraint-induced aphasia therapy primed with intermittent theta burst stimulation in 13 patients with post-stroke aphasia. Med. Sci. Monit. 27:e930100. doi: 10.12659/MSM.930100

CrossRef Full Text | Google Scholar

Allendorfer, J. B., Nenert, R., Vannest, J., and Szaflarski, J. P. (2021b). A pilot randomized controlled trial of intermittent theta burst stimulation as stand-alone treatment for post-stroke aphasia: effects on language and verbal functional magnetic resonance imaging (fMRI). Med. Sci. Monit. 27:e934818. doi: 10.12659/MSM.934818

CrossRef Full Text | Google Scholar

Amad, A., Seidman, J., Draper, S. B., Bruchhage, M. M. K., Lowry, R. G., Wheeler, J., et al. (2017). Motor learning induces plasticity in the resting brain-drumming up a connection. Cereb. Cortex 27, 2010–2021. doi: 10.1093/cercor/bhw048

PubMed Abstract | CrossRef Full Text | Google Scholar

Ardila, A., Bernal, B., and Rosselli, M. (2016). How localized are language brain areas? A review of Brodmann areas involvement in oral language. Arch. Clin. Neuropsychol. 31, 112–122. doi: 10.1093/arclin/acv081

PubMed Abstract | CrossRef Full Text | Google Scholar

Arheix-Parras, S., Barrios, C., Python, G., Cogné, M., Sibon, I., Engelhardt, M., et al. (2021). A systematic review of repetitive transcranial magnetic stimulation in aphasia rehabilitation: leads for future studies. Neurosci. Biobehav. Rev. 127, 212–241. doi: 10.1016/j.neubiorev.2021.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, Z., Zhang, J., and Fong, K. N. K. (2022). Effects of transcranial magnetic stimulation in modulating cortical excitability in patients with stroke: a systematic review and meta-analysis. J. Neuroeng. Rehabil. 19:24. doi: 10.1186/s12984-022-00999-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Baillieux, H., De Smet, H. J., Dobbeleir, A., Paquier, P. F., De Deyn, P. P., and Mariën, P. (2010). Cognitive and affective disturbances following focal cerebellar damage in adults: a neuropsychological and SPECT study. Cortex 46, 869–879. doi: 10.1016/j.cortex.2009.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Brownsett, S. L., Warren, J. E., Geranmayeh, F., Woodhead, Z., Leech, R., and Wise, R. J. (2014). Cognitive control and its impact on recovery from aphasic stroke. Brain 137, 242–254. doi: 10.1093/brain/awt289

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L. (2013). The cerebellum and cognitive function: 25 years of insight from anatomy and neuroimaging. Neuron 80, 807–815. doi: 10.1016/j.neuron.2013.10.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Lin, T., Wu, M., Cai, G., Ding, Q., Xu, J., et al. (2022). Effects of repetitive transcranial magnetic stimulation on upper-limb and finger function in stroke patients: a systematic review and meta-analysis of randomized controlled trials. Front. Neurol. 13:940467. doi: 10.3389/fneur.2022.940467

PubMed Abstract | CrossRef Full Text | Google Scholar

Crosson, B., Sadek, J. R., Bobholz, J. A., Gökçay, D., Mohr, C. M., Leonard, C. M., et al. (1999). Activity in the paracingulate and cingulate sulci during word generation: an fMRI study of functional anatomy. Cereb. Cortex 9, 307–316. doi: 10.1093/cercor/9.4.307

PubMed Abstract | CrossRef Full Text | Google Scholar

De Smet, H. J., Paquier, P., Verhoeven, J., and Mariën, P. (2013). The cerebellum: its role in language and related cognitive and affective functions. Brain Lang. 127, 334–342. doi: 10.1016/j.bandl.2012.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Gholami, M., Pourbaghi, N., and Taghvatalab, S. (2022). Evaluation of rTMS in patients with poststroke aphasia: a systematic review and focused meta-analysis. Neurol. Sci. 43, 4685–4694. doi: 10.1007/s10072-022-06092-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Habas, C., Kamdar, N., Nguyen, D., Prater, K., Beckmann, C. F., Menon, V., et al. (2009). Distinct cerebellar contributions to intrinsic connectivity networks. J. Neurosci. 29, 8586–8594. doi: 10.1523/JNEUROSCI.1868-09.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Harvey, D. Y., Mass, J. A., Shah-Basak, P. P., Wurzman, R., Faseyitan, O., Sacchetti, D. L., et al. (2019). Continuous theta burst stimulation over right pars triangularis facilitates naming abilities in chronic post-stroke aphasia by enhancing phonological access. Brain Lang. 192, 25–34. doi: 10.1016/j.bandl.2019.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Heikkinen, P. H., Pulvermuller, F., Makela, J. P., Ilmoniemi, R. J., Lioumis, P., Kujala, T., et al. (2018). Combining rTMS with intensive language-action therapy in chronic aphasia: a randomized controlled trial. Front. Neurosci. 12:1036. doi: 10.3389/fnins.2018.01036

CrossRef Full Text | Google Scholar

Hu, X. Y., Zhang, T., Rajah, G. B., Stone, C., Liu, L. X., He, J. J., et al. (2018). Effects of different frequencies of repetitive transcranial magnetic stimulation in stroke patients with non-fluent aphasia: a randomized, sham-controlled study. Neurol. Res. 40, 459–465. doi: 10.1080/01616412.2018.1453980

PubMed Abstract | CrossRef Full Text | Google Scholar

Ilkhani, M., Shojaie Baghini, H., Kiamarzi, G., Meysamie, A., and Ebrahimi, P. (2018). The effect of low-frequency repetitive transcranial magnetic stimulation (rTMS) on the treatment of aphasia caused by cerebrovascular accident (CVA). Med. J. Islam Repub. Iran 32:25. doi: 10.14196/mjiri.32.25

CrossRef Full Text | Google Scholar

Jargow, J., Zwosta, K., Korb, F. M., Ruge, H., and Wolfensteller, U. (2021). Low-frequency TMS results in condition-related dynamic activation changes of stimulated and contralateral inferior parietal lobule. Front. Hum. Neurosci. 15:684367. doi: 10.3389/fnhum.2021.684367

PubMed Abstract | CrossRef Full Text | Google Scholar

Kapoor, A. (2017). Repetitive transcranial magnetic stimulation therapy for post-stroke non-fluent aphasia: a critical review. Top. Stroke Rehabil. 24, 547–553. doi: 10.1080/10749357.2017.1331417

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawamura, M., Takahashi, N., and Kobayashi, Y. (2019). Effect of repetitive transcranial magnetic stimulation on the right superior temporal gyrus for severe aphasia caused by damage to the left inferior frontal gyrus. Case Rep. Neurol. 11, 189–198. doi: 10.1159/000500669

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, G., Esposito, R., Motta, C., Casula, E. P., Di Lorenzo, F., Bonnì, S., et al. (2020). Improving visuo-motor learning with cerebellar theta burst stimulation: behavioral and neurophysiological evidence. NeuroImage 208:116424. doi: 10.1016/j.neuroimage.2019.116424

PubMed Abstract | CrossRef Full Text | Google Scholar

Lazar, R. M., and Boehme, A. K. (2017). Aphasia as a predictor of stroke outcome. Curr. Neurol. Neurosci. Rep. 17:83. doi: 10.1007/s11910-017-0797-z

CrossRef Full Text | Google Scholar

Lee, I. T., Huang, C. C., Hsu, P. C., Lin, C. P., and Tsai, P. Y. (2022). Resting-state network changes following transcranial magnetic stimulation in patients with aphasia-a randomized controlled study. Neuromodulation 25, 528–537. doi: 10.1016/j.neurom.2021.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, H., Wu, H., Cheng, H., Wei, D., Wang, X., Fan, Y., et al. (2014). Improvement of white matter and functional connectivity abnormalities by repetitive transcranial magnetic stimulation in crossed aphasia in dextral. Int. J. Clin. Exp. Med. 7, 3659–3668.

Google Scholar

Lv, H., Wang, Z., Tong, E., Williams, L. M., Zaharchuk, G., Zeineh, M., et al. (2018). Resting-state functional MRI: everything that nonexperts have always wanted to know. AJNR Am. J. Neuroradiol. 39, 1390–1399. doi: 10.3174/ajnr.A5527

PubMed Abstract | CrossRef Full Text | Google Scholar

Margulies, D. S., and Petrides, M. (2013). Distinct parietal and temporal connectivity profiles of ventrolateral frontal areas involved in language production. J. Neurosci. 33, 16846–16852. doi: 10.1523/JNEUROSCI.2259-13.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariën, P., and Beaton, A. (2014). The enigmatic linguistic cerebellum: clinical relevance and unanswered questions on nonmotor speech and language deficits in cerebellar disorders. Cerebellum Ataxias 1:12. doi: 10.1186/2053-8871-1-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariën, P., and Borgatti, R. (2018). Language and the cerebellum. Handb. Clin. Neurol. 154, 181–202. doi: 10.1016/B978-0-444-63956-1.00011-4

CrossRef Full Text | Google Scholar

Marien, P., Engelborghs, S., Fabbro, F., and De Deyn, P. P. (2001). The lateralized linguistic cerebellum: a review and a new hypothesis. Brain Lang. 79, 580–600. doi: 10.1006/brln.2001.2569

PubMed Abstract | CrossRef Full Text | Google Scholar

Nachev, P., Coulthard, E., Jäger, H. R., Kennard, C., and Husain, M. (2008). Enantiomorphic normalization of focally lesioned brains. NeuroImage 39, 1215–1226. doi: 10.1016/j.neuroimage.2007.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Nee, D. E., and D'Esposito, M. (2017). Causal evidence for lateral prefrontal cortex dynamics supporting cognitive control. eLife 6:e28040. doi: 10.7554/eLife.28040

PubMed Abstract | CrossRef Full Text | Google Scholar

Oldfield, R. C. (1971). The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 9, 97–113. doi: 10.1016/0028-3932(71)90067-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, C. L., Zhang, G. F., Xia, N., Jin, C. H., Zhang, X. H., Hao, J. F., et al. (2014). Effect of low-frequency rTMS on aphasia in stroke patients: a meta-analysis of randomized controlled trials. PLoS One 9:e102557. doi: 10.1371/journal.pone.0102557

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, C., Zhang, G., Xu, X., Hao, J., Fang, H., Chen, P., et al. (2019). The effect of rTMS over the different targets on language recovery in stroke patients with global aphasia: a randomized sham-controlled study. Biomed. Res. Int. 2019:4589056. doi: 10.1155/2019/4589056

CrossRef Full Text | Google Scholar

Rolls, E. T., and Grabenhorst, F. (2008). The orbitofrontal cortex and beyond: from affect to decision-making. Prog. Neurobiol. 86, 216–244. doi: 10.1016/j.pneurobio.2008.09.001

CrossRef Full Text | Google Scholar

Rosselli, M., Ardila, A., and Bernal, B. (2015). Angular gyrus connectivity model for language: a functional neuroimaging meta-analysis. Rev. Neurol. 60, 495–503. doi: 10.33588/rn.6011.2015039

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartze, M., and Kotz, S. A. (2016). Contributions of cerebellar event-based temporal processing and preparatory function to speech perception. Brain Lang. 161, 28–32. doi: 10.1016/j.bandl.2015.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sroka, M. C., Vannest, J., Maloney, T. C., Horowitz-Kraus, T., Byars, A. W., Holland, S. K., et al. (2015). Relationship between receptive vocabulary and the neural substrates for story processing in preschoolers. Brain Imaging Behav. 9, 43–55. doi: 10.1007/s11682-014-9342-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoodley, C. J., and Schmahmann, J. D. (2010). Evidence for topographic organization in the cerebellum of motor control versus cognitive and affective processing. Cortex 46, 831–844. doi: 10.1016/j.cortex.2009.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Szaflarski, J. P., Nenert, R., Allendorfer, J. B., Martin, A. N., Amara, A. W., Griffis, J. C., et al. (2021). Intermittent theta burst stimulation (iTBS) for treatment of chronic post-stroke aphasia: results of a pilot randomized, double-blind. Sham-controlled trial. Med. Sci. Monit. 27:e931468. doi: 10.12659/MSM.931468

CrossRef Full Text | Google Scholar

Thiel, A., Hartmann, A., Rubi-Fessen, I., Anglade, C., Kracht, L., Weiduschat, N., et al. (2013). Effects of noninvasive brain stimulation on language networks and recovery in early poststroke aphasia. Stroke 44, 2240–2246. doi: 10.1161/STROKEAHA.111.000574

PubMed Abstract | CrossRef Full Text | Google Scholar

Versace, V., Schwenker, K., Langthaler, P. B., Golaszewski, S., Sebastianelli, L., Brigo, F., et al. (2019). Facilitation of auditory comprehension after theta burst stimulation of Wernicke's area in stroke patients: a pilot study. Front. Neurol. 10:1319. doi: 10.3389/fneur.2019.01319

CrossRef Full Text | Google Scholar

Wessel, M. J., and Hummel, F. C. (2018). Non-invasive cerebellar stimulation: a promising approach for stroke recovery? Cerebellum 17, 359–371. doi: 10.1007/s12311-017-0906-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Whitfield-Gabrieli, S., and Nieto-Castanon, A. (2012). Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. 2, 125–141. doi: 10.1089/brain.2012.0073

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoo, W. K., Vernet, M., Kim, J. H., Brem, A. K., Bashir, S., Ifert-Miller, F., et al. (2020). Interhemispheric and intrahemispheric connectivity from the left pars opercularis within the language network is modulated by transcranial stimulation in healthy subjects. Front. Hum. Neurosci. 14:63. doi: 10.3389/fnhum.2020.00063

PubMed Abstract | CrossRef Full Text | Google Scholar

Yourganov, G., Fridriksson, J., Rorden, C., Gleichgerrcht, E., and Bonilha, L. (2016). Multivariate connectome-based symptom mapping in post-stroke patients: networks supporting language and speech. J. Neurosci. 36, 6668–6679. doi: 10.1523/JNEUROSCI.4396-15.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Yushkevich, P. A., Yang, G., and Gerig, G. (2016). ITK-SNAP: an interactive tool for semi-automatic segmentation of multi-modality biomedical images. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2016, 3342–3345. doi: 10.1109/EMBC.2016.7591443

PubMed Abstract | CrossRef Full Text | Google Scholar

Zald, D. H., McHugo, M., Ray, K. L., Glahn, D. C., Eickhoff, S. B., and Laird, A. R. (2014). Meta-analytic connectivity modeling reveals differential functional connectivity of the medial and lateral orbitofrontal cortex. Cereb. Cortex 24, 232–248. doi: 10.1093/cercor/bhs308

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, N., Xia, M., Qiu, T., Wang, X., Lin, C. P., Guo, Q., et al. (2018). Reorganization of cerebro-cerebellar circuit in patients with left hemispheric gliomas involving language network: a combined structural and resting-state functional MRI study. Hum. Brain Mapp. 39, 4802–4819. doi: 10.1002/hbm.24324

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Q. H., Zhu, C. Z., Yang, Y., Zuo, X. N., Long, X. Y., Cao, Q. J., et al. (2008). An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: fractional ALFF. J. Neurosci. Methods 172, 137–141. doi: 10.1016/j.jneumeth.2008.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zumbansen, A., Kneifel, H., Lazzouni, L., Ophey, A., Black, S. E., Chen, J. L., et al. (2022). Differential effects of speech and language therapy and rTMS in chronic versus subacute post-stroke aphasia: results of the NORTHSTAR-CA trial. Neurorehabil. Neural Repair 36, 306–316. doi: 10.1177/15459683211065448

PubMed Abstract | CrossRef Full Text | Google Scholar