Neurobiological Correlates of Inhibition of the Right Broca Homolog during New-Word Learning

Repetitive transcranial magnetic stimulation (rTMS) has demonstrated beneficial effects on motor learning. It would be important to obtain a similar enhancement for verbal learning. However, previous studies have mostly assessed short-term effects of rTMS on language performance and the effect on learning is largely unknown. This study examined whether an inhibition of the right Broca homolog has long-term impact on neural processes underlying the acquisition of new words in healthy individuals. Sixteen young participants trained a new-word learning paradigm with rare, mostly unknown objects and their corresponding words immediately after continuous theta burst stimulation (cTBS) or sham stimulation of right inferior frontal gyrus (IFG) in a cross-over design. Neural effects were assessed with electroencephalography (EEG) source power analyses during the naming task as well as coherence analyses at rest 1 day before and after training. Inhibition of the right Broca homolog did not affect new word learning performance at the group level. Behavioral and neural responses to cTBS were variable across participants and were associated with the magnitude of resting-state alpha-band coherence between the stimulated area and the rest of the brain before stimulation. Only participants with high intrinsic alpha-band coherence between the stimulated area and the rest of the brain before stimulation showed the expected inhibition during naming and greater learning performance. In conclusion, our study confirms that cTBS can induce lasting modulations of neural processes which are associated with learning, but the effect depends on the individual network state.


INTRODUCTION
Learning new words is an essential aspect of human language, which is not limited to childhood as adults also add new words to their mental lexicon. In addition, learning new words in adulthood is particularly relevant for the acquisition of a second language and during recovery in patients with aphasia.
The success of word learning in adults is often limited as compared to childhood, and in particular patients with aphasia often remain disabled despite high intensity speech therapy (Pedersen et al., 2004). Adjuvant interventions that enhance the outcome of learning and of speech and language therapy would therefore be highly welcome. One potential adjuvant therapy consists in neuromodulation via non-invasive cortical stimulation through repetitive transcranial magnetic stimulation (rTMS). RTMS can induce transient excitation or inhibition of cortical areas and their connected nodes beyond the duration of the stimulation . The rationale for inhibiting or exciting brain regions during word learning comes from the observation that successful acquisition of new words is associated with a dominant left hemispheric neural pattern before, during and after the training task (Breitenstein et al., 2005), while an overactivity of right language nodes may be deleterious (Hamilton et al., 2011). More importantly, good learners showed stronger activations of brain regions in the left hemisphere, in particular temporal regions, whereas poor learners activated preferentially a more diffuse bilateral network including right temporal and right inferior frontal areas (Wong et al., 2007;Mei et al., 2008). Similar observations were made in the context of word retrieval. Indeed, naming new words (as compared to familiar words) has been associated with increased activity in the left inferior frontal cortex (Broca's area), the left temporal area and the left inferior parietal lobe (Cornelissen et al., 2004;Gronholm et al., 2005;Hulten et al., 2009). In post-stroke aphasia, a large body of research has suggested that the degree of activity in the left hemisphere was more critical for naming performance than activity in the right hemisphere (Thiel et al., 1998;Cao et al., 1999;Winhuisen et al., 2005).
On the other hand, some authors argue that the recruitment of the right hemisphere is not related to bad performance but rather to effortful language processing (Raboyeau et al., 2008) or to adaptation and plasticity after stroke (Blasi et al., 2002;Heiss et al., 2013). In consequence, the contribution of the right hemisphere to verbal learning is still under debate and the optimal parameters for rTMS are unknown.
There is some evidence that language skills may be facilitated by rTMS. Positive impact has been reported on speech repetition accuracy (Restle et al., 2012) or picture naming (Topper et al., 1998;Mottaghy et al., 1999Mottaghy et al., , 2006. For instance, excitatory rTMS protocols over Wernicke's area briefly decrease the response latency in a picture naming task in healthy people (Mottaghy et al., 2006). Inhibition of the right pars triangularis of the inferior frontal gyrus (IFG) led to decreased production latencies while inhibition of the right pars opercularis seemed to increase reaction times (Naeser et al., 2011). In stroke patients, several studies have revealed that inhibitory rTMS over the right Broca homolog can improve naming performance and/or decrease naming latencies in stroke patients or improve global language recovery Naeser et al., 2011;Kindler et al., 2012;Thiel et al., 2013).
However, despite the evidence for momentary effects of rTMS on naming performance, it is largely unknown how it influences processes related to word learning. A better understanding of the long-term impact on word learning would be crucial for an application of rTMS as a complementary tool to speech therapy and language training. Such neurobiological information would likely lead to more efficient treatment methods in association with language therapies.
The present sham-controlled cross-over study therefore aimed to investigate the effect of rTMS on the ability to learn new words, as well as on the neural mechanisms related to naming new words in healthy participants. Unlike previous studies, we were not interested in short-term effects immediately after stimulation, but in long-term neural and behavioral effects on learning. Based on previous evidence for a negative influence of the right hemisphere and because of the excellent safety profile of inhibitory protocols which makes them easily applicable also in patients with recent brain lesions, we applied inhibitory rTMS over the right Broca homolog.
Neural changes were investigated with high-density electroencephalography (EEG). We first examined neural and behavioral effects at the group level in the hypothesis that rTMS modulates naming-related neural processing and neural interactions in language networks which translate into improved learning performance. Second, given that the behavioral effect of rTMS has been previously shown to be variable and dependent on the neural state of the brain before the stimulation (Silvanto and Pascual-Leone, 2008;Rizk et al., 2013;Nicolo et al., 2015;Vallence et al., 2015), we tested the hypothesis that network states before stimulation explain inter-individual variability in rTMS effects on learning and naming-related neural processing. In particular, alpha-band functional connectivity between the stimulated area and the rest of the brain was tested as a promising predictor of the response to rTMS (Rizk et al., 2013).

Subjects
Sixteen native French speakers (mean age 25 years, range 19-35, 4 men) without neurological or psychiatric disease were recruited after written informed consent. All were right-handed as determined by self-report and by the Edinburgh Handedness Scales (Oldfield, 1971) and had normal or corrected-to-normal visual acuity. All participants had comparable educational level (see Supplementary Table S1). Procedures were approved by the ethics committee of Geneva, Switzerland (project number: 10-220). All procedures were in accordance with international ethical standards on human experimentation and with the Helsinki Declaration of 1975, as revised in 2000. The participants were paid for their participation.

Overview of Experimental Design
We performed a single-blinded within-subject sham-controlled crossover study. The complete procedure is summarized in Figure 1 and took place over 5 days within the same week. Participants underwent two learning-to-name sessions on Tuesdays and Thursdays, during which they learned the name and the meaning of pictures of very rare objects. Immediately before learning, they received either inhibitory or sham stimulation over the right IFG, pars triangularis (Brodmann area 45) in counterbalanced order. Participants were blinded to the treatment arm. One day before and after each training session, i.e., on Mondays, Wednesdays, and Fridays, we obtained high-density EEG recording for pre-and post-tests of naming (event-related picture naming tasks) as well as during task-free resting states. FIGURE 1 | Experimental design. Participants (n = 16) underwent two learning sessions on Tuesdays and Thursdays. During each session, participants first received active or sham cTBS over the right Broca homolog followed by training of one list on a computer. One day before and after each training session, overt naming skills were tested for the learned and non-learned lists.

Material
We used a picture naming task with rare real objects allowing the study of new word learning in healthy subjects with low probability of baseline knowledge while avoiding artificial pseudo-words. The stimuli consisted of 100 pictures of ancient or rare objects, tools or musical instruments and their corresponding words transformed in black and white line drawings (examples in Supplementary Material). All words have very low lexical frequency and were largely unknown at baseline. Two lists of 50 items (lists A and B) were constituted and matched on lexical frequency, first phoneme, phonological neighborhood and length in syllables and phonemes [from the French database Lexique (New et al., 2004), see details in Supplementary Table  S2]. This allowed having a trained list and a control list in each learning session.

Assessment of Naming Skills
Participants were tested individually in a soundproof dark room. The presentation of stimuli was controlled by E-Prime software (Psychology Software Tools, Inc., Sharpsburg, MD, USA). Pictures were presented in constant size of 240 pixels × 245 pixels (about 4.5 • of visual angle) on a black screen (60 cm from their chest). In a familiarization phase, all 100 pictures and their corresponding names were presented on the screen one by one. Each item was presented for 3000 ms. In the picture naming task that followed, all items were presented in a pseudo-random order and were preceded by two warming-up filler trials (pictures of two familiar objects). An experimental trial began with a fixation cross presented for 500 ms. Then a blank screen preceded the appearance of the picture. Participants were requested to produce overtly the word corresponding to the picture as soon as they could. If they did not know the answer, they were asked to overtly say "no." The picture remained on screen 3000 ms and a blank screen lasting 2000 ms was displayed before the next trial. Pictures of each list were presented twice (hence resulting in 100 trials per list) in pseudo-random order and in four separate blocks. The experiment lasted about 20 min with a break after each block.
Learning performance was computed as the difference in correct answers between post-and pre-training assessments.
Production latencies were measured by means of a voice key and were digitized for further systematic latency and accuracy check with speech analysis software (CheckVocal 2.2.6, Protopapas, 2007). Two kinds of incorrect responses were analyzed: no-responses (the participants indicated that he/she does not know the word by answering "no") and errors (i.e., the participant produced a different name than the one expected for the picture, a phonologically transformed word and/or autocorrections during articulation).

Training Sessions
The learning task consisted of four runs of 50 pictures from lists A or B in pseudorandom order. The order of the two lists was counterbalanced between subjects. Pictures of the objects were presented in the black and white line drawing format and in a photograph picture format along with their definition on a computer screen. Participants were required to use the computer keyboard to write the corresponding word and then press 'enter.' If the word was unknown or incorrect, the individuals had the option of hearing the spoken word or reading the correct spelling of the word. Each picture stimulus remained on screen until the participant found or copied the correct picture-name association. Then, the next item was displayed automatically. The number of words produced correctly without help was computed after each run (four runs in total) and recorded as learning curve during training. For this learning stage, participants did not have time limits.

Continuous Theta-Burst Stimulation (cTBS)
Continuous theta burst stimulation (cTBS), a more recent form of rTMS, has been shown to induce a decrease of neural activity in a specific network (Nyffeler et al., 2006). TBS has the advantage of inducing longer aftereffects while requiring shorter stimulation time than conventional rTMS (Huang et al., 2005;Nyffeler et al., 2006).
CTBS was delivered with a biphasic waveform through a MagPro X100 system (Medtronic Functional Diagnostics). The stimulator was connected with a figure-of-eight coil (MCF-B65) with a diameter of 2 mm × 75 mm and a geometrically identical coil for the sham condition (MCF-P-B65). Sham stimulation produced the same noise as true stimulation, but induced no magnetic field within the underlying neural tissue.
The cTBS protocol consisted of continuous trains of 801 pulses, applied in 267 bursts. Each burst contained three pulses at 30 Hz with an interburst interval of 167 ms. Total duration of a train was 44 s (Kindler et al., 2012).
The coil was positioned over the right posterior inferior frontal gyrus (pIFG, pars triangularis, BA 45) using landmarks of the international 10-20 EEG system. The Münster T2T-Converter (Deppe et al., 2003) coordinates were used to calculate the scalp position relative to EEG coordinates which corresponded to this specific localization. Stimulation was applied between electrodes C4 and F8 (10-20, 1.26/0.74; Talairach space x/yr/z, 58/31/22). The coil was held tangentially to the skull with the handle pointing upward. Stimulation intensity was expressed as percentage of stimulator output and was set to 90% of the individual resting motor threshold of the small left hand muscles (Kindler et al., 2012).
Adverse events were checked after each rTMS intervention and on subsequent days by the investigators including an experimented neurologist. No negative side effects were observed except for transient discomfort during stimulation.

EEG Recordings and Preprocessing
Electroencephalography was recorded during a task-free resting state and during the picture naming tasks (event related) at each pre-learning and post-learning session. An Active-Two EEG system (Biosemi V.O.F., Amsterdam, The Netherlands) with 128 electrodes was used to digitize signals at a rate of 512 Hz (filters: DC to 104 Hz, 3 dB/octave slope).
Resting-state recordings were obtained continuously for 10 min while participants were awake and kept their eyes-closed. For event-related recordings during naming, epochs lasting from 600 ms before to 1000 ms after the picture onset were extracted. Subjects were instructed to avoid eye movements and blinking, swallowing or any movement other than required for the task.
Epochs or resting-state periods with movement artifacts, eye blinking, other noise, or signs of sleepiness were excluded by off-line visual inspection. Bad channels containing prolonged artifacts were ignored from further analyses. Artifact-free epochs were recalculated against the average reference (all non-excluded channels).

Source Reconstruction
Analyses in source space were performed using the software Matlab (The MathWorks Inc.) with the open-source toolbox NUTMEG 1  and its functional connectivity mapping (FCM) toolbox . A leadpotential with 10 mm grid spacing (894 voxels) was computed using a spherical head model with anatomical constraints (Spinelli et al., 2000) based on the segmented gray matter of the standard Montreal Neurological Institute (MNI) brain. An adaptive spatial filter (scalar minimum variance beamformer) was used to reconstruct neural oscillations in source space from surface sensors (Sekihara et al., 2004).

Time-Frequency Decomposition
Event-related power changes (erPOW; Pfurtscheller and Andrew, 1999) were assessed to investigate oscillatory modulation of sets of neurons in response to the naming task.
Adaptive spatial filter weights were calculated for each subject from all artifact-free epochs of a given condition. In order to optimize the beamformer for the respective frequency bands, we obtained separate weights for each of the following bands: 1-20, 21-30, 31-45, and 55-95 Hz (Dalal et al., 2008).
Signals were Fourier-transformed using a sliding Hanning window of 500 ms width shifted in time steps of 50 ms. Fourier coefficients were projected to source space with the adaptive spatial filter and power was computed at each time window across all trials for each of the following frequency bands: delta (2-3 Hz), theta (4-7 Hz), alpha (8-12 Hz), low beta (13-20 Hz), high beta (21-30 Hz), gamma (31-45 Hz), and high-gamma (55-95 Hz). Time-frequency power values were log transformed and baseline corrected by subtracting a pre-stimulus baseline power average from −600 to 0 ms.

Resting-State Coherence (rsCOH)
Adaptive spatial filter weights were calculated for each subject from 300 artifact-free resting-state epochs of 1 s duration in the bandwidth from 1 to 20 Hz.
Functional connectivity between two time series was quantified as the imaginary component of coherence (IC) (Nolte et al., 2004). Although conservative, this index was chosen because it is robust to volume conduction and avoids distortions due to spatial leakage of inverse solutions .
Imaginary component of coherence was computed in source space across 5 min of artifact-free epochs for four standard frequency bands: delta (1-3 Hz), theta (4-7 Hz), alpha (8-12 Hz), and low beta (13-20 Hz). In order to obtain a global measure of functional connectivity at each voxel, we computed the graph theoretical measure of node degree in weighted networks (Newman, 2004) as the summed IC with all other voxels. This measure indicates the overall importance of each voxel in the brain network (Stam and van Straaten, 2012).
Variations in functional connectivity magnitude can be due to fluctuations in signal to noise ratio between participants or conditions. To avoid this problem, node degree maps were normalized. This was achieved by subtracting the mean node degree of all voxels of the subject from the values at each individual voxel and by dividing by the standard deviation over all voxels, hence obtaining z-scores .

Statistical Analyses
For behavioral measures, repeated measures ANOVAs were computed on rates of correct productions with session (prelearning, post-learning) and stimulation condition (true, sham) as within subject factors.
The brain region that we targeted with cTBS, i.e., the right pIFG (Brodmann area 45), as well as its left homolog, were defined as regions of interest (ROIs) using the automated anatomical labeling (AAL) atlas (Tzourio-Mazoyer et al., 2002). An additional ROI was functionally defined as the area showing significant EEG activation during picture naming at pre-test 1 (p < 0.05, 5% false discovery rate, FDR).
Event-related power changes and resting-state coherence (rsCOH) at each ROI were tested against the null hypotheses of zero mean with a t-test for one sample, as well as against the null hypothesis of zero difference between before and after learning, or between true and sham stimulation, with paired t-tests. Furthermore, we tested neural changes for correlations with learning effects using a Pearson correlation analysis. Normal distribution of variables was confirmed by visual inspection of the data as well as with a Kolmogorov-Smirnov goodness-of-fit test.
To correct for the family wise error of testing ROI eventrelated power at multiple time-frequency data points, we performed a permutation test. In brief, at each of 2000 permutation loops, the values of a random combination of subjects were inverted by multiplication with −1 and subjected to the same statistical test. The largest number of significant consecutive windows at p < 0.05 (uncorrected) across all data points was recorded. The significance of the number of consecutive windows in the real dataset was then determined from its position in the empirical distribution obtained through permutations. Since we compare against the maximum value across all data points, we effectively control for multiple testing. Time-frequency windows that belonged to a window cluster that was larger than 95% of clusters obtained in permutations were considered significant at p < 0.05, corrected.

Behavioral Results
At pre-test, naming scores did not differ between A and B lists (t < 1). Figure 2A presents the improvement of naming scores from pre-learning to post-learning sessions in each stimulation condition. There was a main effect of session indicating improved naming after training [F(1,15) = 195.3, p < 0.0001], no main effect of stimulation type (F < 1) and no interaction between session and stimulation (F < 1). The between subject variability in learning was high, with learnt words varying from 13 to 70%. Regarding naming errors, whereas the amount of responses indicating that the participant did not know the word ("no" responses) were bound to decrease after training, other kinds of erroneous responses (i.e., lexical errors, mainly confusion between learnt words, and phonological errors) increased (see Figure 2B), with a trend for increased rates of errors after cTBS compared to sham [t(15) = −2.11, p = 0.052]. Figure 3 shows the training curve over the four runs during the training sessions. Accuracy improved over the four training runs as evidenced by a main effect of run [F(3,115) = 108.4, p < 0.0001]. However, there was no significant main effect of stimulation type (F < 0.2) nor an interaction between run and stimulation (F < 0.7).
Hence, inhibition of the right Broca homolog with cTBS did not influence the ability to learn new words at the group level but slightly increased the amount of errors relative to "don't know" responses.

EEG Results
We first investigated whether the neural processing of word naming was modulated by learning and stimulation by assessing event-related power changes. A voxel-wise timefrequency decomposition during the pre-training session (pretest 1) revealed a significant cluster of naming-related power modulation which survived correction for testing multiple voxels, time windows, and frequency bands (p < 0.05, FDR corrected). High-gamma (50-95 Hz) power increased between 25 and 800 ms after picture presentation in a left temporo-occipital brain area (Figure 4, peak MNI coordinates: −42 −90 20). This area was therefore defined as ROI for neural processing of naming.
Besides the increase in high-gamma power, naming during the pre-training session also induced, in the left temporo-occipital ROI, a power decrease in beta frequency bands between 200 and 800 ms after stimulus presentation. In contrast, increased power was observed in theta band between 0 and 465 ms as well as in delta band between 0 and 320 ms (p < 0.05, corrected, Figures 4-6).
Beta and high-gamma power modulations remained unchanged at the post-test 1 day after training (Figure 5). In contrast, delta and theta band power modulations were significantly altered 1 day after learning and the alteration depended significantly on the stimulation condition. Theta and delta power modulations decreased 1 day after learning preceded by sham stimulation (p < 0.05, corrected), whereas the opposite occurred after true stimulation (p < 0.05, corrected). The difference between true and sham condition was significant (p < 0.05, corrected). This effect was observed during a time window between 120 and 460 ms after stimulus presentation for the theta band and between 120 and 710 ms for the delta band. However, theta and delta band power changes were not associated with learning performance (max correlation at any time point r < 0.49, p > 0.05 uncorrected) or with naming errors (r < 0.45, p > 0.05, uncorrected) in any stimulation condition.
No learning and stimulation effects were found at the left and right inferior frontal ROIs or on rsCOH 1 day after learning at any ROI.
As stated in the Introduction, the effect of rTMS can depend on neural states before the stimulation. To test for statedependency in our study, we examined the relationship between resting-state connectivity in the alpha frequency band, eventrelated neural changes and behavioral effects induced by cTBS over the right IFG.
A significant positive correlation was observed between resting-state connectivity in the alpha band before stimulation and learning performance (r = 0.62, p = 0.01; Figure 7A) indicating larger new-word learning in participants with greater alpha-band connectivity in the right IFG. This was the case only for true stimulation and not observed in the sham condition (r = 0.23, p = 0.37).
Although cTBS had no significant effect on rsCOH at the group level, it decreased alpha-band resting-state connectivity in participants with high pre-stimulation connectivity but increased values in participants with low pre-stimulation connectivity as indicated by the high negative correlation (r = −0.74, p = 0.001; Figure 7B) between pre-stimulation and change in restingstate connectivity in the alpha band. Furthermore, connectivity changes at the alpha-band induced by stimulation were negatively correlated with learning performance (r = −0.6, p = 0.014; Figure 7C) and positively correlated with naming errors after stimulation (r = 0.51, p = 0.041). This indicated better new-word learning performance and lower lexical confusions in participants in whom resting-state connectivity between the right Broca homolog and the rest of the brain was decreased by stimulation.
Alpha-band coherence before stimulation further influenced naming-related alpha power modulations (Figures 7D-F). In participants with high resting-state alpha-band coherence before stimulation, cTBS produced an enhancement of naming-related alpha power at the right Broca homolog, while it produced no change or even decreased power in participants with low baseline coherence. We then averaged the naming-related alpha power changes of each participant across all time windows >200 ms after picture presentation and found a positive correlation with baseline alpha-band coherence (r = 0.68, p = 0.004; Figure 7G). In addition, a significant positive correlation was also found between event-related alpha power increases and learning performance (r = 0.52, p = 0.039; Figure 7H).
These pattern of correlations in the alpha band were specific to cTBS stimulation and not observed after sham condition (r < 0.32, p > 0.05). Furthermore, they were not observed in other frequency bands (max correlations with changes at right IFG coherence for beta, theta, delta bands r < 0.430, p > 0.05).
When gender was included as confounding factor, all correlations remained significant (p < 0.03) and unpaired t-tests (men vs. women) on our variables were not significant (p > 0.3). The associations were therefore not driven by gender differences.

DISCUSSION
Non-invasive brain stimulation can enhance learning through practice in healthy people or during recovery of neurological deficits (Tanaka et al., 2011;Wessel et al., 2016). This has been shown in particular for motor learning (Reis et al., 2008;Herzfeld et al., 2014). The present study examined if a similar benefit can be obtained for verbal learning, and, if yes, by which underlying neural mechanisms.
Overall, a single session of inhibition of the right Broca homolog did not facilitate new word learning nor impact the learning curve during the training phase. Our concomitant monitoring of neural effects with EEG provides insights into the reasons for this lack of benefit. Although our protocol did induce long-term modulations of naming-related oscillations in left temporo-occipital language areas, they were not sufficient for improving behavioral learning. The lack of benefit at the group level seemed to be due to endogenous variability of the neural state at baseline. The neural and behavioral response to cTBS in an individual depended on the state of intrinsic functional connectivity of the stimulated brain region. Participants with high intrinsic alpha-band coherence between the target area and the rest of the brain before stimulation were more likely to show the expected inhibition and showed greater learning performance. Participants with low intrinsic alpha-band coherence between right inferior frontal regions and the rest of the brain before stimulation showed a paradoxical increase in coherence after cTBS and showed poorer verbal learning.

Neural Processing during Naming and Its Modulation after cTBS
Time-frequency analyses revealed power modulations of oscillations in several frequency bands prior and after training. Increased power in the high-gamma band and decreased power in beta bands in the left temporo-occipital region were not modulated as a function of learning and therefore likely reflect invariant neural processes during word naming. The activation of the temporo-occipital junction, the modulations in gamma and beta bands and their time course suggest that these neural modulations correspond to visual object processing (Bookheimer et al., 1998). Beta power suppression can reflect working memory in naming (Piai et al., 2015). However, its time line (the first 800 msec post picture onset) and the combination with gamma rhythm in occipital regions rather support functional mechanisms of visual processing (Michalareas et al., 2016) and (visual) working memory (Honkanen et al., 2015). Honkanen et al. (2015) recently showed that gamma oscillations in the visual cortex are involved in the extraction and maintenance of object features and associated to the content of visual working memory. These results are also compatible with the observation that the occipital cortex remains active during a large part of the naming processes (Hassan et al., 2015).
Delta and theta oscillations were also modulated during naming and were further influenced by learning and stimulation. Theta oscillations, in particular, have been previously described in association with working memory load (Fuentemilla et al., 2010) or encoding and retrieval of episodic memory (Hasselmo and Stern, 2014). The modulations observed in the theta band in the first 500 ms post picture onset are in line with lexical processes occurring in this time-window (Indefrey, 2011) and are compatible with recent report that theta oscillations are affected by lexical frequency (Piai et al., 2014) although in a different time-window. In the present study, the difference between pre-and post-training corresponds to the difference between unknown and familiar words, which may mimic the lexical frequency effect. Delta oscillations have been associated with a large variety of cognitive processes including motivational drive (Knyazev, 2012), inhibition of sensorial afferences or attention to internal processing during mental tasks (Harmony, 2013), memory retrieval (Fernandez et al., 1998), and cortical plasticity . However, theta-delta modulations induced by inhibition of the right Broca homolog were not sufficient for improving language learning, hence suggesting that they do not reflect critical processes for new-word learning.
Since we observed a naming-related processing mainly in the left temporo-occipital junction, it might be more efficient to target this area directly. A recent study using anodal (excitatory) tDCS over Wernicke's area, which is located much more closely to the left temporo-occipital junction, has indeed reported a facilitation of new word acquisition in healthy people (Meinzer et al., 2014).
It might be tempting to speculate that more than a single stimulation session would have been necessary to obtain learning effects at the group level. However, since we do observe robust neural effects even after one session, we consider this possibility unlikely, although we cannot exclude that more sessions may have reduced inter-individual variability.

State-Dependency and Variability in Response to cTBS
The main reason for the lack of an overall behavioral effect at the group level is the large inter-individual variability in the response to cTBS. Our EEG monitoring identified one particular factor contributing to this variability. Behavioral and neural effects of cTBS were dependent on pre-stimulation alpha-band coherence between the stimulated area and the rest of the brain. Such a dependency of the response to stimulation on the neural state before stimulation has been well known for local neural activations at the stimulation site (Silvanto and Pascual-Leone, 2008;. Alpha-band coherence was also previously reported to shape the behavioral effect of cTBS at FIGURE 5 | Comparison of induced brain responses over the left temporo-occipital region. The panels in the left column show naming-related power changes 1 day after stimulation, whereas the right-hand panels represent the difference (after minus before). Significant power increase at low frequencies was observed after cTBS in comparison to sham.
Frontiers in Human Neuroscience | www.frontiersin.org FIGURE 6 | Naming-and stimulation-related power modulations at the left temporo-occipital ROI. The left column shows power modulations (mean ± standard error of mean) before (pre) and after each type of stimulation (post), relative to baseline. The right column visualizes changes induced by stimulation. Empty symbols represent significant time windows (p < 0.05, corrected). Theta and delta power increased 1 day after true cTBS, whereas the opposite was observed after sham stimulation.
FIGURE 7 | Network states modulate the effects of cTBS. Learning and EEG after-effects of cTBS depended on the magnitude of resting-state alpha coherence between the right IFG and the rest of the brain before stimulation. Participants with large pre-stimulation coherence learned more words (A) and decreased alpha-band coherence after stimulation (B), while participants with low pre-stimulation connectivity showed the opposite pattern. The reduction of coherence between the right Broca homolog and the rest of the brain was correlated with better learning performance (C). In participants with high pre-stimulation coherence, cTBS also produced the expected local inhibition during naming, as indicated by an increase of event-related alpha power (D,E; black rectangles indicate significant differences at p < 0.05). In contrast, participants with low pre-stimulation alpha-band coherence showed no power modulation after stimulation (F). Hence, resting-state coherence (rsCOH) before stimulation was associated with greater alpha power increases between 200 and 800 ms after picture presentation (G), which was in turn correlated with better learning (H). the right parietal cortex on spatial attention (Rizk et al., 2013). Here, we reproduced a state-dependency at the network level in an independent population, at a different network, and for a different behavior, hence suggesting that this represents a general response pattern of the human brain to cTBS. Moreover, we identified mechanisms by which rsCOH may influence the response to stimulation. Alpha-band coherence was correlated with increases in event-related alpha power after stimulation, hence suggesting that resting-state network interactions translate into specific local patterns of brain oscillations during naming. Event-related alpha power increases are thought to represent inhibition during tasks in order to liberate neural resources at more critical brain areas (Pfurtscheller, 1992;Klimesch et al., 2007;Jensen et al., 2012). Hence, it seems that, in participants with large pre-stimulation alpha-band coherence between the stimulated area and the rest of the brain, our stimulation protocol had the intended effect of local inhibition during learning and naming and, consequently, improved learning. However, no or even the opposite effect seems to have occurred in about half of the participants who had relatively low baseline coherence. This may suggest that cTBS did not have a strictly inhibiting effect in all participants, but may even have produced a paradoxical excitation in some cases.

Local Effects of cTBS
In comparison to previous studies on EEG effects of rTMS (Fuggetta et al., 2008;Hamidi et al., 2009;Veniero et al., 2011;Noh et al., 2012;Vernet et al., 2013), we did not observe robust effects on rsCOH or erPOW under the stimulation site at the group level. This is explained by the fact that our recordings were made 1 day after stimulation, when excitability modulations are expected to have disappeared (Grossheinrich et al., 2009;Noh et al., 2012;Rizk et al., 2013;Vernet et al., 2013). Yet, our findings suggest that the neural response to stimuli can remain modified beyond the duration of traditional aftereffects.

CONCLUSION
In order to improve the efficacy of rTMS as adjuvant language therapy in the future, it will be critical to understand the neural signatures of verbal learning that need to be modulated in learning. Furthermore, the present results indicate that there is a need to understand and control factors underlying inter-individual variability in order to apply stimulation more selectively in participants who are likely to respond. Thereby, our results support the usefulness of EEG network imaging for monitoring and predicting the aftereffects induced by noninvasive brain stimulation.

AUTHOR CONTRIBUTIONS
PN acquired data, performed EEG analyses, and wrote the paper. RF acquired data, performed behavioral analyses, and wrote the paper. ML designed the study, acquired data, performed behavioral analyses, and wrote the paper. AG designed the study, acquired data, performed EEG analyses, and wrote the paper.