Interhemispheric and Intrahemispheric Connectivity From the Left Pars Opercularis Within the Language Network Is Modulated by Transcranial Stimulation in Healthy Subjects

Neural activity related to language can be modulated within widespread networks following learning or in response to disruption—including the experimental application of noninvasive brain stimulation. However, the spatiotemporal characteristics of such modulation remain insufficiently explored. The present study combined transcranial magnetic stimulation (TMS) and electroencephalography (EEG) to explore the modulation of activity across the language network following continuous theta-burst stimulation (cTBS) of the left pars opercularis. In 10 healthy subjects (21 ± 2 years old, four females), neuronavigated cTBS was delivered over the left pars opercularis of the frontal operculum (part of the traditional Broca’s area) at 80% of active motor threshold (AMT) stimulation intensity. Real cTBS and sham cTBS were performed in two different visits separated by at least 48 h. Before, immediately, and 10 min after cTBS, 30 single pulses of TMS were delivered to the left pars opercularis at 80% of the resting motor threshold (RMT), whereas EEG was simultaneously recorded. We examined the cTBS-induced modulation of phase locking values (PLVs) between the TMS-evoked potentials (TEPs) recorded over the pars opercularis and those recorded over its right-hemispheric homolog area, the left supramarginal area, and the left superior temporal area in different EEG frequency bands and different time windows following cTBS. cTBS to the left pars opercularis induced within the gamma band: (1) a significant increase in TEP phase synchronization between the left and right pars opercularis at an early time window (250–350 ms) following cTBS; and (2) significantly increased PLV with the left supramarginal area and the left superior temporal area at a later time window (600–700 ms). In the theta and delta band, cTBS to the left pars opercularis induced significantly increased phase synchronization of TEPs between the left pars opercularis and the posterior left hemispheric language areas at a late time window. In sham condition, there was a significant decrease in TEP phase synchronization of the high beta band between left pars opercularis and left superior temporal area at 200–300 ms. These results contribute to characterize the dynamics of the language network and may have implications in the development of noninvasive stimulation protocols to promote the language rehabilitation in aphasia patients.

Introduction

Language is processed in widely distributed neural networks including the frontal operculum (Broca’s region) in the dominant hemisphere and its right hemispheric homolog, as well as parietotemporal areas of the left hemisphere (Pulvermüller, 2002; Pulvermüller et al., 2009; Giraud and Poeppel, 2012; Poeppel et al., 2012). Measuring the information flow between those core regions is important for understanding the function of the language network (Singer, 1999; Varela et al., 2001; Buzsaki, 2004). Information flow across nodes of a network can be studied by recording synchronizations of spontaneous cortical oscillations. For instance, a study (Giraud et al., 2007) tried to link auditory structural characteristics of speech and spontaneous brain oscillations within networks important for speech perception and production. Others have examined transient event-related phase synchronizations between two brain areas (Lachaux et al., 1999; Pulvermüller, 2002) or the phase resetting of a cortical oscillation by a perceptual stimulus (Friston et al., 1997; Rizzuto and Kahana, 2001). However, sensory stimuli or cognitive task manipulations can affect in parallel several cortical areas. On the other hand, controlled perturbations applied directly to well-defined brain regions offer a superior strategy to gain insights into the casual interactions between nodes of a given large-scale brain network. Here we applied this strategy to gain insights into the dynamics across the language network.

We applied single pulses of brain magnetic resonance imaging (MRI)-guided (neuronavigated) transcranial magnetic stimulation (TMS) to the pars opercularis of the left frontal operculum and recorded TMS-evoked cortical responses with electroencephalography (EEG) to gain insights into the connectivity dynamics across interhemispheric and intrahemispheric connections of the language network. The use of neuronavigated TMS is important to ensure precise targeting of a given brain region (Gugino et al., 2001; Bashir et al., 2011). We examined the modulation of TMS-evoked cortical responses induced by continuous theta-burst stimulation (cTBS), a noninvasive stimulation protocol that has been shown not only to decrease the excitability of the targeted area (Huang et al., 2005) but also to modulate brain oscillations and synchronizations (Vernet et al., 2013). CTBS can provide longer-lasting modulatory effects with shorter stimulation time than traditional repetitive TMS, which can have important clinical implications.

We hypothesized that such an approach, using task-free and externally triggered oscillations, would provide more direct physiological evidence for connectivity mechanisms. Furthermore, synchronization and desynchronization of these task-free oscillations in response to the plasticity-inducing cTBS protocol could provide information relevant to a future noninvasive intervention for the rehabilitation of patients with aphasia (Miniussi and Thut, 2009). Specifically, we hypothesized that cTBS to the left pars opercularis would change interhemispheric synchronization of TMS-evoked potential (TEP) between the left and the right pars opercularis and intrahemispheric interactions between the left pars opercularis and the left supramarginal and left superior temporal areas.

Materials and Methods

Participants

Ten young, healthy, right-handed adults [the average Edinburgh Handedness Inventory (Oldfield, 1971); score, 83] volunteered to participate in the study (21 ± 2 years old; range, 18–24 years old; four females). None of them had a history of psychiatric or neurological conditions, and all had normal neurological and medical examinations, and Mini-Mental State Examination scores in the normal range (Miltner et al., 1999; Burgess and Ali, 2002; Buschman and Miller, 2007; de Diego-Balaguer et al., 2011). Participants visited twice at least 48 h apart for sham and real repetitive TMS (rTMS) session; half of the participants started with the sham session first. Participants were not taking any medication known to affect motor cortical excitability at the time of the study and did not have any contraindications to TMS. All tolerated the TMS without any adverse effect or complication. All gave their written informed consent to the study, which followed international guidelines for the use of TMS (Rossi et al., 2009), had been approved by the local institutional review board (Beth Israel Deaconess Medical Center, Boston, MA, USA), and was conducted in adherence to the Declaration of Helsinki.

Experimental Setup

All participants underwent a brain MRI using a 3-T GE machine (GE Healthcare, Waukesha WI, USA) with three-dimensional (3D) modified driven-equilibrium Fourier transform protocol–prepared fast spoiled gradient echo sequence (0.94 × 0.94 × 1-mm resolution; echo time = 2.9 ms; flip angle = 15°) to generate a high-resolution anatomical 3D image to guide TMS.

The stimulation setup consisted of a Nexstim stimulator (Nexstim Limited, Helsinki, Finland) for single-pulse TMS and a MagPro stimulator (MagVenture A/S, Farum, Denmark) for the cTBS intervention. The Nexstim eXimia Neuronavigation system was used to ensure that the exact location of the left pars opercularis was targeted as defined by each individual’s brain MRI (Vernet et al., 2013).

Motor-evoked potentials (MEPs) induced by single-pulse TMS were recorded using surface electromyography (EMG) with pre-gelled, disposable Ag/AgCl electrodes attached to the skin. The active electrode was placed over the first dorsal interosseus (FDI) muscle, the reference electrode over the metacarpophalangeal joint, and the ground electrode over the wrist. The EMG signal was acquired at 3 kHz, filtered (10–500 Hz), amplified, displayed, and stored for off-line analysis.

EEG was recorded with a 60-channel TMS-compatible EEG system (eXimia EEG, Nexstim Limited). This system is designed to avoid amplifier saturation after TMS pulse by using a sample-and-hold circuit that keeps the input of the amplifiers constant from 100 μs prestimulus to 2 ms poststimulus (Virtanen et al., 1999).

Experimental Session

Participants were seated in a comfortable chair, with a headrest, and with their hands supinated and resting on their laps. Participants were monitored for drowsiness and asked to keep their eyes open during TMS (Vernet et al., 2013). Relaxation of the target muscle (FDI muscle) was controlled by continuous EMG monitoring. All participants wore earplugs to protect them from possible acoustic trauma (Rossi et al., 2009) and reduce contamination of TEPs by auditory responses to the clicks produced by the discharge of the TMS coil.

The left M1 optimal scalp location was determined as the scalp location from which single-pulse TMS induced MEPs of maximum peak-to-peak amplitude in the right FDI muscle. Once this optimal spot was identified, the neuronavigation system was used to ensure consistent coil placement and orientation at the optimal spot. Resting motor threshold (RMT) was defined as the lowest stimulus intensity of the Nexstim stimulator capable of inducing MEPs of ≥50-μV peak-to-peak amplitude in at least 5 out of 10 trials. The RMT was obtained to set the stimulation intensity for subsequent single-pulse TMS. Active motor threshold (AMT) was defined as the lowest stimulus intensity of the MagPro stimulator capable of inducing visible twitches in the FDI muscle in half of the trials while the participants maintained a contraction of the FDI muscle at approximately 20% of the maximal voluntary contraction. The AMT was used to set the stimulation intensity for the cTBS protocol.

cTBS was applied to the left pars opercularis with parameters similar to those used by Huang et al. (2005): three pulses at 50 Hz, with an interval of 200 ms between the last pulse of a triplet and the first pulse of a triplet, for a total number of 600 pulses. The intensity was fixed at 80% of AMT. Because of limitations in our experimental setup, the interstimulus interval was 240 ms compared to the interstimulus interval of 200 ms in the original paradigm introduced by Huang et al. (2005). Thus, in our cTBS paradigm, the triplet repetition rate was approximately 4.17 Hz instead of 5 Hz, both frequencies being included in the theta band. Note that previous research from our group has shown that this cTBS protocol can induce clear changes of oscillations within the motor cortex (Vernet et al., 2013).

Before cTBS, a series of 30 neuronavigated single pulses of TMS was applied to the left pars opercularis to evoke TEPs (Figure 1). These pulses were applied at an intensity of 80% of RMT. After cTBS, two series of 30 single pulses were applied to the same area and at the same intensity, one immediately after cTBS (T0) and one at 10 min later (T10).

FIGURE 1

Figure 1. Post cTBS T0 TEP recorded at the FC5 electrode. cTBS, continuous theta-burst stimulation; TEP, transcranial magnetic stimulation (TMS)-evoked potential.

EEG Data Preprocessing

After recording EEG, the data were preprocessed using an open MATLAB toolbox EEGLAB¹. Data were band-pass filtered between 1 and 50 Hz. Artifacts induced by transient magnetic fields and eye blinks were eliminated using independent component analysis. The noise-free data were segmented into epochs from 300 ms before the stimulus onset to 1,000 ms after the stimulus onset. Epochs were excluded from further analysis if they contained significant physiological artifacts (amplitude exceeding ±5 μV) in the time range from 50 to 1,000 ms after the stimulus onset at any site over all electrodes.

Phase Locking Value

Phase locking value (PLV) was used to quantify the functional connectivity between pairs of regions of interest (ROIs). Phase locking value is a well-known index to measure phase synchronization between two time signals acquired from two different electrodes in the same time interval and frequency band (Pulvermüller, 2002). Phase locking value was selected as the functional connectivity measure because it can robustly estimate the interaction between two signals with a constant time delay. Indeed, systematic time delays are expected between the responses to TMS induced at different ROIs. In order to quantify the degree of phase synchrony between two signals S_x(t) and S_y(t), the instantaneous phases φ_x(t) and φ_y(t) were first evaluated using the Hilbert transform. The Hilbert transform of S(t) is defined as:

$\tilde{S}(t)=\frac{1}{\pi }P.V.{\displaystyle {\int }_{\infty }^{\infty }\frac{s(\tau )}{t-\tau }d\tau }(1)$

where $\tilde{S}(t)$ is the Hilbert transform of the time series S(t), and P.V. denotes the Cauchy principal value. The instantaneous phase φ(t) can then be estimated by:

$\phi (t)=arctan\frac{\tilde{S}(t)}{S(t)}(2)$

The PLV was evaluated using the following definition:

$\text{PLV}=\left|\langle e^{i\Delta \phi (t)}\rangle \right|(3)$

where Δφ(t) = φ_x(t) − φ_y(t), and $\left|\langle \cdot \rangle \right|$ is the averaging operator. The PLV ranges from 0 to 1, where the value close to one represents that two signals are synchronized with a constant time lag and the value close to 0 represents that the two signals are temporally independent with each other.

In the present study, PLV features were evaluated for three electrode pairs, FC5-FC6, FC5-CP3, and FC5-CP5, where FC5, FC6, CP3, and CP5 correspond to left pars opercularis, right pars opercularis, left supramarginal area, and left superior temporal area, respectively. The PLV features were evaluated for delta (1–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), low beta (13–21 Hz), high beta (21–30 Hz), and gamma (30–50 Hz) frequency bands using a sliding window with 100-ms width and 50% overlap.

Statistical Analysis

We employed a Friedman test to evaluate the difference of PLV given different conditions, baseline, T0, and T10, for different frequency bands, different time windows, and different pairs of electrodes. Then, as post hoc test, Wilcoxon signed–rank test was utilized to evaluate the difference of PLV measures between conditions. The significance level was set to p < 0.05 after multiple-comparisons correction using Bonferroni correction.

Results

Figure 2 summarizes the significant results of modulation of TEP synchronization induced by cTBS of the left pars opercularis. Please note that no other significant modulations of TEP synchronization across the studied ROIs were found. By EEG frequency bands, the significant findings were as follows:

FIGURE 2

Figure 2. cTBS was applied on left pars opercularis area with three pulses at 50 Hz, with an interval of 200 ms between the last pulse of a triplet and the first pulse of a triplet, for a total number of 600 pulses. The figure indicates changes in phase synchronization of TEP by real cTBS on left pars opercularis area of the language network. (A) Top view. (B) Left lateral view. Black text: real cTBS condition. Blue text: sham condition. cTBS, continuous theta-burst stimulation; TEP, TMS-evoked potential.

Changes in Gamma-Band Synchronization of TEPs in Real cTBS Condition

In the gamma band (30–50 Hz) at a 250- to 350-ms time window, the phase synchronization of TEPs between the left and right pars opercularis significantly increased at T0 after real cTBS compared to baseline (corrected p = 0.009). At a later time window (600–700 ms), the phase synchronization of TEPs in the gamma band between the left pars opercularis and the left supramarginal region significantly increased at both T0 and T10 compared to baseline (corrected p = 0.022). Similarly, at the same time window (600–700 ms), the phase synchronization of TEPs in the gamma band between the left pars opercularis and the left superior temporal region significantly increased at T0 (but not T10) compared to baseline (corrected p = 0.005; Figure 3A).

FIGURE 3

Figure 3. Changes of phase locking value (PLV) by real cTBS on the left pars opercularis. cTBS, continuous theta-burst stimulation. FC5, FC6, CP3, and CP5 correspond to left pars opercularis, right pars opercularis, left supramarginal area, and left superior temporal area, respectively. *p < 0.05, **p < 0.01. (A) Gamma band changes. (B) Theta band changes. (C) Delta band chnages.

Changes in Theta-Band Synchronization of TEP in Real cTBS Condition

In the theta band (4–8 Hz) in a 450- to 550-ms time window, the phase synchronization of TEPs between the left pars opercularis and the left supramarginal region significantly increased at T0 after real cTBS compared to baseline (corrected p = 0.037; Figure 3B).

Changes in Delta-Band Synchronization of TEP in Real cTBS Condition

In the delta band (1–4 Hz) in a wide time window (350–600 ms), the phase synchronization of TEPs between the left pars opercularis and the left superior temporal area significantly increased at T10 after real cTBS compared to baseline (corrected p = 0.009). However, the synchronization of TEPs between the left pars opercularis and left supramarginal in the delta band significantly decreased at an earlier time window (250–350 ms; Figure 3C).

Changes in Beta2-Band Synchronization of TEP in Sham cTBS Condition

In the high beta band (21–30 Hz) in a time window (200–300 ms), the phase synchronization of TEPs between the left pars opercularis and the left superior temporal area significantly decreased at T0 after sham cTBS compared to baseline (corrected p = 0.014; Figure 4).

FIGURE 4

Figure 4. Changes of phase-locking value (PLV) by sham cTBS. FC5 and CP5 correspond to left pars opercularis and left superior temporal area, respectively. *p < 0.05.

Discussion

cTBS to the left pars opercularis induced significant changes in phase synchronization of TEPs evoked from the left pars opercularis. Specifically, we found: (1) a significant increase in the gamma frequency between the left and right pars opercularis at an early (250–350 ms) time window and a similar increase between the left pars opercularis and the left posterior language areas at a later time period (600–700 ms); and (2) a significant increase in the theta and delta bands between the left pars opercularis and the posterior language areas. These results provide novel insights into the interhemispheric and intrahemispheric dynamics of the language network in healthy young individuals.

cTBS is a TMS protocol that has been primarily used to induce a decrease of excitability in the stimulated area. Indeed, it has been shown that cTBS applied over M1 can induce depression of the MEPs subsequently evoked by single-pulse TMS (Huang et al., 2005). However, recent studies using TMS coupled with EEG, magnetoencephalography (MEG) and/or magnetic resonance spectroscopy have revealed that cTBS also has an impact on local and remote oscillatory systems, on synchronization, and on inhibitory γ-aminobutyric acid (GABA) concentration level (Allen et al., 2014; Du et al., 2018).

Here, we found that cTBS over the left pars opercularis induced an increase of synchronization in the gamma band in early TEPs (250 to 350 ms time window) and in late TEPs (600–700 ms). GABAergic interneurons are known to be the source of the generation of gamma oscillation (Wang and Buzsáki, 1996). Our finding of increased gamma-band synchronization after cTBS might be caused by increased activity of local inhibitory interneurons subtending the inhibitory interaction between the left and right pars opercularis region. Studies using spinal epidural recordings have shown that cTBS applied over M1 suppresses the earliest I1-wave (Di Lazzaro et al., 2005), which is primarily derived from GABA-related interneuronal circuits within the human motor cortex (Ziemann et al., 1998). Thus, cTBS might have decreased cortical excitability in the left pars opercularis and increased inhibition, which led to a change in TMS-evoked EEG gamma synchronization. Alternatively, the increase of gamma coherence might be related to the gamma frequency (50 Hz) embedded in the cTBS pattern.

Interesting deficits of gamma oscillations have been reported in patients with schizophrenia. Ferrarelli et al. (2008) showed that gamma activity evoked by single-pulse TMS applied over the dorsolateral prefrontal cortex was significantly decreased within the first 100 ms after TMS in schizophrenic patients compared to healthy controls. In such patients, spontaneous gamma activity did not show any differences between patients and healthy controls, emphasizing the utility of TMS and specifically the use of TMS–EEG to reveal physiological abnormalities.

Interhemispheric interactions between the left and right pars opercularis regions have been studied as a way to promote language function recovery following stroke. Some have argued that a stroke can impair the normal interhemispheric inhibition balance between the lesioned area and its homolog in the intact hemisphere. If so, decreasing the excitability of the latter area, for example, with noninvasive brain stimulation, might decrease the excessive inhibition it exerts on the former area and thus promote recovery. Consistent with such notions, Naeser et al. (2005) observed significant and persistent improvement in picture-naming abilities in patients with nonfluent aphasia after low frequency rTMS targeting the right frontal operculum. In healthy individuals, similar effects of rTMS on interhemispheric inhibitory interactions between homolog areas of the language network have been shown in studies using positron emission tomography (Thiel et al., 2006), functional MRI (Andoh and Paus, 2011), and behavioral responses (Cappelletti et al., 2008). In the present study, we show for the first time that such interhemispheric interactions might be subtended by synchronization in the gamma band, and that appears possible to modulate it by noninvasive stimulation of the left pars opercularis.

Beyond interhemispheric interactions, our results confirm the relevance of interregional gamma band synchronization within brain networks previously shown in cognition and perception (Sarnthein et al., 1998; Rose and Büchel, 2005; Buschman and Miller, 2007) and language learning (de Diego-Balaguer et al., 2011). For instance, Miltner et al. (1999) showed that intrahemispheric coherence in the gamma band increased during the formation of percepts and memory, linguistic processing, and other behavioral and perceptual functions. Another study showed a greater functional connectivity between the frontal and parietal cortices within the gamma band during successful recollection, as opposed to merely experiencing a feeling of familiarity (Burgess and Ali, 2002). In the language network, an MEG study by Duesburg et al. (2012) reported that a covert verb generation task increased gamma band synchronization among regions activated during the task, in particular, the inferior frontal gyrus and the medial frontal gyrus. A similar result of gamma band phase locking was observed during a silent reading task (Hirata et al., 2010). Our results provide further support for the relevance of phase coherence among gamma rhythms within the language network.

It is also worth noting that local cell assemblies in the cortex are thought to form functional units through synchronous activation in the gamma range (Engel and Singer, 2001; Bastiaansen and Hagoort, 2013). Oscillations may function as nested sets, where lower-frequency activity binds higher-frequency activity (Lisman and Buzsáki, 2008). Theta bursts occur during the up-phase of delta activity (Lakatos et al., 2005), and gamma bursts occur during the up-phase of theta activity (Lisman and Buzsáki, 2008). Thus, it has been suggested that local modules of high-frequency gamma can be functionally linked in temporal synchrony with distant modules via lower frequency theta and gamma band regulation. Bastiaansen et al. (2005) showed that the left hemisphere theta activity might be related to the processing of words with extensive semantic stores. The present results, showing increases of gamma TEP synchronization in a late time window (600–700 ms) between the left pars opercularis and the left supramarginal area, may be related to local processes associated with feature binding (Gray and Singer, 1989; Tallon-Baudry et al., 1998), whereas the increases in TEP synchronization in the theta band at the 450 to 550 ms time window between the left pars opercularis area and left supramarginal area might serve as a long-range synchronization between frontal and parietal cortices (von Stein and Sarnthein, 2000; Mellem et al., 2013).

Delta oscillations, on the other hand, have been mainly associated with slow-wave sleep and anesthesia, which are states without conscious awareness. However, delta oscillations might also play a role during wakefulness. For instance, high delta (approximately 3 Hz) appears to organize gamma oscillations during memory retrieval tasks (Burgess and Ali, 2002). Thus, the dissociation of changes in PLV induced by real cTBS within the delta band might be related to different functions and also warrant further study.

We also found a significant decrease in phase synchronization of high beta band in sham condition only between left pars opercularis and left superior temporal areas at T0. Although the PLV was relatively low compared to real cTBS condition, this decrease might have been related to functional changes in this network. There was some evidence that shows decreased coherence in beta band coherence after stimulation of the motor cortex with inhibitory protocol (Mima et al., 2003; Tamura et al., 2005), but also after periods of rest (Fuggetta et al., 2007). Recent evidence also showed a decrease in beta wave at 100–200 ms after sham cTBS (Opitz et al., 2015), in congruence with our sham result. Thus, the decrease in high beta coherence we observed in the sham condition could be related either to a placebo effect of the stimulation, or to the fact that the participants were at rest during the stimulation. We cannot entirely exclude the possible influence of the previous stimulation session due to relatively short washout period (at least 48 h), although studies using a single-session cTBS protocol with usual intensity (such as the one used in the present study) have been reported to last its effects approximately only 20–60 min (Huang et al., 2005; Chung et al., 2016).

Study Limitations

The present study has a number of limitations that should be considered when interpreting the results. Most importantly, we did not explore the functional correlation of TMS-induced oscillations with language testing. Therefore, it is not possible to experimentally assess whether these oscillation patterns are related to speech production. The left pars opercularis area is well known to be part of the network for speech function. Thus, it seems reasonable to assume that it has a functional role in language, which might be subtended by oscillatory activity similar to the observed TMS-induced oscillations. Nonetheless, this hypothesis awaits future experimental evaluation. In addition, it is important to note that our number of subjects is small, and they were all quite young. Further studies are needed to assess the relevance of the present findings for older adults and elderly individuals. Similarly, we have fewer females than males and an insufficient number of participants to look at gender effects. This too would warrant further studies. Finally, the relevance of our findings to patients with a stroke that might have caused aphasia needs to be examined.

Conclusion

Our results provide novel insights into the interhemispheric and intrahemispheric dynamics of activity across nodes of the language network and illustrate the potential of using noninvasive brain stimulation for their controlled modulation. Such findings may have implications for the future development of noninvasive stimulation protocols to promote recovery of aphasia.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Beth Israel Deaconess Medical Center, Boston, MA, USA. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

W-KY, SB, and AP-L designed the study. A-KB, FI-M, ME, and MV collected the data. W-KY and J-HK analyzed the data. W-KY, MV, J-HK, A-KB, FI-M, ME, SB, and AP-L wrote the main manuscript text. All authors reviewed the manuscript.

Funding

This work was supported in part by grants from CIMIT (Dr. AP-L), the National Center for Research Resources: Harvard Clinical and Translational Science Center (UL1 RR025758). Dr. AP-L serves on the scientific advisory boards for Nexstim, Neuronix, Starlab Neuroscience, Neosync, and Novavision, and is a listed inventor on issued and pending patents on the real-time integration of TMS with EEG and MRI. Dr. MV was supported by the Fyssen Foundation (France). Dr. W-KY was supported by a National Research Foundation of Korea grant funded by the Korean Government MEST, Basic Research Promotion Fund (NRF-2013R1A1A2012562).

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.

Acknowledgments

We thank Andrea Vatulas, Cara Burzynski and Ann Connor for their administrative support.

Footnotes

1. ^ http://sccn.ucsd.edu/eeglab/

References

Allen, C. P., Dunkley, B. T., Muthukumaraswamy, S. D., Edden, R., Evans, J. C., Sumner, P., et al. (2014). Enhanced awareness followed reversible inhibition of human visual cortex: a combined TMS, MRS and MEG study. PLoS One 9:e100350. doi: 10.1371/journal.pone.0100350

PubMed Abstract | CrossRef Full Text | Google Scholar

Andoh, J., and Paus, T. (2011). Combining functional neuroimaging with off-line brain stimulation: modulation of task-related activity in language areas. J. Cogn. Neurosci. 23, 349–361. doi: 10.1162/jocn.2010.21449

PubMed Abstract | CrossRef Full Text | Google Scholar

Bashir, S., Edwards, D., and Pascual-Leone, A. (2011). Neuronavigation increases the physiologic and behavioral effects of low-frequency rtms of primary motor cortex in healthy subjects. Brain Topogr. 24, 54–64. doi: 10.1007/s10548-010-0165-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastiaansen, M., and Hagoort, P. (2013). Event-induced theta responses as a window on the dynamics of memory. Cortex 39, 967–992. doi: 10.1016/s0010-9452(08)70873-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastiaansen, M. C., van der Linden, M., ter Keurs, M., Dijkstra, T., and Hagoort, P. (2005). Theta responses are involved in lexical-semantic retrieval during language processing. J. Cogn. Neurosci. 17, 530–541. doi: 10.1162/0898929053279469

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, A. P., and Ali, L. (2002). Functional connectivity of gamma EEG activity is modulated at low frequency during conscious recollection. Int. J. Psychophysiol. 46, 91–100. doi: 10.1016/s0167-8760(02)00108-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Buschman, T., and Miller, E. (2007). Top-down versus bottom-up control of attention in the prefrontal and posterior parietal cortices. Science 315, 1860–1862. doi: 10.1126/science.1145017

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzsaki, G. (2004). Neuronal oscillations in cortical networks. Science 304, 1926–1929. doi: 10.1126/science.1099745

PubMed Abstract | CrossRef Full Text | Google Scholar

Cappelletti, M., Fregni, F., Shapiro, K., Pascual-Leone, A., and Caramazza, A. (2008). Processing nouns and verbs in the left frontal cortex: a transcranial magnetic stimulation study. J. Cogn. Neurosci. 20, 707–720. doi: 10.1162/jocn.2008.20045

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, S., Hill, A. T., Rogasch, N. C., Hoy, K. E., and Fitzgerald, P. B. (2016). Use of theta-burst stimulation in changing excitability of motor cortex: a systematic review and meta-analysis. Neurosci. Biobehav. Rev. 63, 43–64. doi: 10.1016/j.neubiorev.2016.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

de Diego-Balaguer, R., Fuentemilla, L., and Rodriguez-Fornells, A. (2011). Brain dynamics sustaining rapid rule extraction from speech. J. Cogn. Neurosci. 23, 3105–3120. doi: 10.1162/jocn.2011.21636

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Lazzaro, D. V., Pilato, F., Saturno, E., Oliviero, A., Dileone, M., Mazzone, P., et al. (2005). Theta-burst repetitive transcranial magnetic stimulation suppresses specific excitatory circuits in the human motor cortex. J. Physiol. 565, 945–950. doi: 10.1113/jphysiol.2005.087288

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, X., Rowland, L. M., Summerfelt, A., Wijtenburg, A., Chiappelli, J., Wisner, K., et al. (2018). TMS evoked N100 reflects local GABA and glutamate balance. Brain Stimul. 11, 1071–1079. doi: 10.1016/j.brs.2018.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Duesburg, S., Vinette, S. A., Cheung, M. J., and Pang, E. W. (2012). Theta-modulated gamma-band synchronization among activated regions during a verb generation task. Front. psychol. 3:195. doi: 10.3389/fpsyg.2012.00195

PubMed Abstract | CrossRef Full Text | Google Scholar

Engel, A., and Singer, W. (2001). Temporal binding and the neural correlates of sensory awareness. Trends Cogn. Sci. 5, 16–25. doi: 10.1016/s1364-6613(00)01568-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrarelli, F., Massimini, M., Peterson, M. J., Riedner, B. A., Lazar, M., Murphy, M. J., et al. (2008). Reduced evoked gamma oscillations in the frontal cortex in schizophrenia patients: a TMS/EEG study. Am. J. Psychiatry 165, 996–1005. doi: 10.1176/appi.ajp.2008.07111733

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K., Buechel, C., Fink, G., Morris, J., Rolls, E., and Dolan, R. (1997). Psychophysiological and modulatory interactions in neuroimaging. NeuroImage 6, 218–229. doi: 10.1006/nimg.1997.0291

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuggetta, G., Pavone, E. F., Fiaschi, A., and Manganotti, P. (2007). Acute modulation of cortical oscillatory activities during short trains of high-frequency repetitive transcranial magnetic stimulation of the human motor cortex: a combined EEG and TMS study. Hum. Brain Mapp. 29, 1–13. doi: 10.1002/hbm.20371

PubMed Abstract | CrossRef Full Text | Google Scholar

Giraud, A. L., Kleinschmidt, A., Poeppel, D., Lund, T. E., Frackowiak, R. S., and Laufs, H. (2007). Endogenous cortical rhythms determine cerebral specialization for speech perception and production. Neuron 56, 1127–1134. doi: 10.1016/j.neuron.2007.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Giraud, A. L., and Poeppel, D. (2012). Cortical oscillations and speech processing: emerging computational principles and operations. Nat. Neurosci. 15, 511–517. doi: 10.1038/nn.3063

PubMed Abstract | CrossRef Full Text | Google Scholar

Gray, C., and Singer, W. (1989). Stimulus-specific neuronal oscillations in orientation columns of cat visual cortex. Proc. Natl. Acad. Sci. U S A 86, 1698–1702. doi: 10.1073/pnas.86.5.1698

PubMed Abstract | CrossRef Full Text | Google Scholar

Gugino, L. D., Romero, R. J., Aglio, L., Titone, D., Ramirez, M., Pascual-Leone, A., et al. (2001). Transcranial magnetic stimulation coregistered with MRI: a comparison of a guided versus blind stimulation technique and its effect on evoked compound muscle action potentials. Clin. Neurophysiol. 112, 1781–1792. doi: 10.1016/s1388-2457(01)00633-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirata, M., Goto, T., Barnes, G., and Umekawa, Y. (2010). Language dominance and mapping based on neuromagnetic oscillatory changes: comparison with invasive procedures: clinical article. J. Neurosurg. 112, 528–538. doi: 10.3171/2009.7.jns09239

CrossRef Full Text | Google Scholar

Huang, Y. Z., Edwards, M. J., Rounis, E., Bhatia, K. P., and Rothwell, J. C. (2005). Theta burst stimulation of the human motor cortex. Neuron 45, 201–206. doi: 10.1016/j.neuron.2004.12.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Lachaux, J. P., Rodriguez, E., Martinerie, J., and Varela, F. J. (1999). Measuring phase synchrony in brain signals. Hum. Brain Mapp. 8, 194–208. doi: 10.1002/(sici)1097-0193(1999)8:4<194::aid-hbm4>3.0.co;2-c

PubMed Abstract | CrossRef Full Text | Google Scholar

Lakatos, P., Shah, A. S., Knuth, K. H., Ulbert, I., Karmos, G., and Schroeder, C. E. (2005). An oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex. J. Neurophysiol. 94, 1904–1911. doi: 10.1152/jn.00263.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisman, J., and Buzsáki, G. (2008). A neural coding scheme formed by the combined function of gamma and theta oscillations. Schizophr. Bull. 34, 974–980. doi: 10.1093/schbul/sbn060

PubMed Abstract | CrossRef Full Text | Google Scholar

Mellem, M. S., Friedman, R. B., and Medvedev, A. V. (2013). Gamma- and theta-band synchronization during semantic priming reflect local and long-range lexical-semantic networks. Brain Lang. 127, 440–451. doi: 10.1016/j.bandl.2013.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Miltner, W. H., Braun, C., Arnold, M., Witte, H., and Taub, E. (1999). Coherence of gamma-band EEG activity as a basis for associative learning. Nature 397, 434–436. doi: 10.1038/17126

PubMed Abstract | CrossRef Full Text | Google Scholar

Mima, T., Siebner, H. R., Nagamine, T., and Shibasaki, H. (2003). Low-frequency rTMS over lateral premotor cortex induces lasting changes in regional activation and functional coupling of cortical motor areas. Clin. Neurophysiol. 114, 1628–1637. doi: 10.1016/s1388-2457(03)00063-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Miniussi, C., and Thut, G. (2009). Combining TMS and EEG offers new prospects in cognitive neuroscience. Brain Topogr. 22, 249–256. doi: 10.1007/s10548-009-0083-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Naeser, M., Martin, P., Nicholas, M., Baker, E., Seekins, H., Kobayashi, M., et al. (2005). Improved picture naming in chronic aphasia after TMS to part of right Broca–s area: An open-protocol study. Brain Lang. 93, 95–105. doi: 10.1016/j.bandl.2004.08.004

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

Opitz, A., Legon, W., Mueller, J., Barbour, A., Paulus, W., and Tyler, W. J. (2015). Is sham cTBS real cTBS? The effect on EEG dynamics. Front. Hum. Neurosci. 8:1043. doi: 10.3389/fnhum.2014.01043

PubMed Abstract | CrossRef Full Text | Google Scholar

Poeppel, D., Emmorey, K., Hickok, G., and Pylkkanen, L. (2012). Towards a new neurobiology of language. J. Neurosci. 32, 14125–14131. doi: 10.1523/JNEUROSCI.3244-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Pulvermüller, F. (2002). A brain perspective on language mechanisms: from discrete neuronal ensembles to serial order. Prog. Neurobiol. 67, 85–111. doi: 10.1016/s0301-0082(02)00014-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pulvermüller, F., Kherif, F., Hauk, O., Mohr, B., and Nimmo-Smith, I. (2009). Distributed cell assemblies for general lexical and category-specific semantic processing as revealed by fMRI cluster analysis. Hum. Brain Mapp. 30, 3837–3850. doi: 10.1002/hbm.20811

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizzuto, D. S., and Kahana, M. J. (2001). An Autoassociative Neural Network Model of Paired-Associate Learning. Neural. Comput. 13, 2075–2092. doi: 10.1162/089976601750399317

PubMed Abstract | CrossRef Full Text | Google Scholar

Rose, M., and Büchel, C. (2005). Neural coupling binds visual tokens to moving stimuli. J. Neurosci. 25, 10101–10104. doi: 10.1523/jneurosci.2998-05.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossi, S., Hallett, M., Rossini, P. M., and Pascual-Leone, A. (2009). Safety, ethical considerations, and application guidelines for the use of transcranial magnetic stimulation in clinical practice and research. Clin. Neurophysiol. 120, 2008–2039. doi: 10.1016/j.clinph.2009.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarnthein, J., Petsche, H., Rappelsberger, P., Shaw, G. L., and von Stein, A. (1998). Synchronization between prefrontal and posterior association cortex during human working memory. Proc. Natl. Acad. Sci. U S A 95, 7092–7096. doi: 10.1073/pnas.95.12.7092

PubMed Abstract | CrossRef Full Text | Google Scholar

Singer, W. (1999). Neuronal synchrony: a versatile code for the definition of relations? Neuron 24, 49–65. doi: 10.1016/s0896-6273(00)80821-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Tallon-Baudry, C., Bertrand, O., Peronnet, F., and Pernier, J. (1998). Induced gamma-band activity during the delay of a visual short-term memory task in humans. J. Neurosci. 18, 4244–4254. doi: 10.1523/jneurosci.18-11-04244.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, Y., Hoshiyama, M., Nakata, H., Hiroe, N., Inui, K., Kaneoke, Y., et al. (2005). Functional relationship between human rolandic oscillations and motor cortical excitability: an MEG study. Eur. J. Neurosci. 21, 2555–2562. doi: 10.1111/j.1460-9568.2005.04096.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiel, A., Schumacher, B., Wienhard, K., Gairing, S., Kracht, L. W., Wagner, R., et al. (2006). Direct demonstration of transcallosal disinhibition in language networks. J. Cereb. Blood Flow Metab. 26, 1122–1127. doi: 10.1038/sj.jcbfm.9600350

PubMed Abstract | CrossRef Full Text | Google Scholar

Varela, F., Lachaux, J. P., Rodriguez, E., and Martinerie, J. (2001). The brainweb: Phase synchronization and large-scale integration. Nat. Rev. Neurosci. 2, 229–239. doi: 10.1038/35067550

PubMed Abstract | CrossRef Full Text | Google Scholar

Vernet, M., Bashir, S., Yoo, W. K., Perez, J. M., Najib, U., and Pascual-Leone, A. (2013). Insights on the neural basis of motor plasticity induced by theta burst stimulation from TMS-EEG. Eur. J. Neurosci. 37, 598–606. doi: 10.1111/ejn.12069

PubMed Abstract | CrossRef Full Text | Google Scholar

Virtanen, J., Ruohonen, J., Näätänen, R., and Ilmoniemi, R. (1999). Instrumentation for the measurement of electric brain responses to transcranial magnetic stimulation. Med. Biol. Eng. Comput. 37, 322–326. doi: 10.1007/bf02513307

PubMed Abstract | CrossRef Full Text | Google Scholar

von Stein, A., and Sarnthein, J. (2000). Different frequencies for different scales of cortical integration: from local gamma to long range alpha/theta synchronization. Int. J. Psychophysiol. 38, 301–313. doi: 10.1016/s0167-8760(00)00172-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. J., and Buzsáki, G. (1996). Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model. J. Neurosci. doi: 10.1523/jneurosci.16-20-06402.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Ziemann, U., Corwell, B., and Cohen, L. G. (1998). Modulation of plasticity in human motor cortex after forearm ischemic nerve block. J. Neurosci. 18, 1115–1123. doi: 10.1523/jneurosci.18-03-01115.1998

PubMed Abstract | CrossRef Full Text | Google Scholar