Abstract
Word retrieval in bilingual speakers partly depends on executive control systems in the left prefrontal cortex – including dorsolateral prefrontal cortex (DLPFC). We tested the hypothesis that DLPFC modulates word production of words specifically in a second language (L2) by measuring the effects of anodal transcranial direct current stimulation (anodal-tDCS) over the DLPFC on picture naming and word translation and on event-related potentials (ERPs) and their sources. Twenty-six bilingual participants with “unbalanced” proficiency in two languages were given 20 min of 1.5 mA anodal or sham tDCS (double-blind stimulation design, counterbalanced stimulation order, 1-week intersession delay). The participants then performed the following tasks: verbal and non-verbal fluency during anodal-tDCS stimulation and first and second language (L1 and L2) picture naming and translation [forward (L1 → L2) and backward (L2 → L1)] immediately after stimulation. The electroencephalogram (EEG) was recorded during picture naming and translation. On the behavioral level, anodal-tDCS had an influence on non-verbal fluency but neither on verbal fluency, nor on picture naming and translation. EEG measures revealed significant interactions between Language and Stimulation on picture naming around 380 ms post-stimulus onset and Translation direction and Stimulation on translation around 530 ms post-stimulus onset. These effects suggest that L2 phonological retrieval and phoneme encoding are spatially and temporally segregated in the brain. We conclude that anodal-tDCS stimulation has an effect at a neural level on phonological processes and, critically, that DLPFC-mediated activation is a constraint on language production specifically in L2.
Introduction
Learning a second language (L2) appears to be effortless early in life [early bilinguals; age of acquisition (AoA) of L2 before the age of 6–7] but becomes more demanding with age (; ). Bilingualism can be defined as the fluent production of more than one language although proficiency in each language can vary. At the neural level, the dominant view is that language representations share the same brain structures, although degree of overlap between languages depends on multiple factors such as AoA, proficiency, degree of immersion in L2, and linguistic (e.g., morphosyntactic) similarity between languages (; ).
Language production is dependent on executive functions (EFs). EFs refer to multiple meta-cognitive abilities that are required to optimize performance when multiple cognitive processes have to be coordinated and comprise processes such as mental shifting, updating, and inhibitory control (IC) of prepotent responses (Miyake et al., 2000). Bilingual language production is particularly dependent on EF, and according to all models of bilingual language processing, verbal fluency involves an interaction between language proficiency and EF. Specifically, EFs are required to control language production in the first and second languages (Postma, 2000; ; ; Metuki et al., 2012) because language representations must be monitored both within the language being spoken and between languages to select the appropriate target vocabulary and syntactical construction. In the model of , IC is recruited to avoid interference when the dominant language (L1) competes for access while speaking the non-dominant second language (L2; ; Sunderman and Kroll, 2006).
Bilateral prefrontal cortices, and most notably the dorsolateral prefrontal cortex (DLPFC) and the dorsal anterior cingulate cortex (dACC; ; ), are assumed to be crucially involved in EFs. DLPFC is mostly involved in selecting a target response and suppression of inappropriate response and the manipulation of information held in working memory (Rowe et al., 2000), while ACC is involved in conflict monitoring and error detection (). For example, patients with lesions in prefrontal cortex show reduced verbal fluency, which in the most extreme cases can resemble mutism (; ). Thus, assume that the brain regions activated during EFs are also recruited during language production for fluent bilingual speakers (Price, 2000, 2010).
Language production in bilingual speakers have been measured at the word (lexeme) level on picture naming and word translation tasks which put differential requirements on EF: relatively low in picture naming and relatively high in word translation. Models of picture naming assume that word retrieval requires a number of processing stages: visual attentional control, access to semantic representations (from 0 to 200 ms after image presentation; Noppeney et al., 2004; Whitney et al., 2012), lexeme selection (to 275 ms post-picture presentation), speech monitoring and inhibition of lexical competitors followed by phonological code retrieval (at around 350 ms), syllabification (at around 450 ms), and phonetic encoding (motor programming of individual syllables) at around 600 ms which precedes articulation (; ).
The Revised Hierarchical Model (RHM) assumes a separate lexicon for each language (; ). These separate lexica are linked to a single semantic or conceptual system, which contains word meaning representations. The RHM assumes asymmetry in the strength of the connections between words and concepts in two languages, thereby postulating that associations between lexical nodes and concepts are weaker in a less proficient language (Stein et al., 2009). Therefore, in low proficiency bilinguals, access to L2 is assumed to require mediation via translation of L1 equivalents, which in turn leads to a processing cost and therefore slower responses. According to the RHM, forward translation (from L1 to L2) will be slower compared to backward translation (from L2 to L1) for bilinguals who are less proficient due to the later acquisition of L2 after early childhood or when L2 is the less dominant language (; Sholl et al., 1995); importantly, this model does not claim the importance of EF for the difference in L1 and L2 processing.
Abutalebi and Green’s IC model (; ) offers an alternative explanation for the translation asymmetry effects. According to their model, forward translation requires inhibition of L1 lexical nodes to allow spoken production of L2 words. For bilinguals who have “unbalanced” proficiency in two languages, the lexical nodes in L2 are assumed to be less active than lexical nodes in L1. This requires additional attentional resources for processing the L1 lexical nodes, i.e., for redundant nodes to be suppressed. Consequently, forward translation is slower when compared to backward translation because the tasks require differential inhibitory demands (Price, 2000; Sunderman and Kroll, 2006; ).
Behavioral and fMRI studies report significant associations between neural activity in the prefrontal cortex and performance on picture naming and word translation tasks (; Price, 2000; ). However, it is not clear whether the DLPFC is necessary for word retrieval in these tasks. FMRI studies report correlational information. Most relevant to neurocognitive models of bilingualism is the putative role of the DLPFC in cognitive control within a language (e.g., naming in L1 or L2) vs. between languages (e.g., translation). Current evidence for a causal role of the DLPFC in cognitive control during language processing in bilingual speakers is restricted to lesion studies, which can only confirm whether the DLPFC is sufficient for language control and not whether it is a necessary mechanism. Transcranial direct current stimulation (tDCS) is a safe, reliable, transient, and effective method to study the necessary and sufficient conditions for different cognitive functions such as language processing (Nitsche et al., 2008). tDCS is assumed to induce changes in resting potentials of neurons and these appear to be polarity specific, that is, anodal-tDCS increases and cathodal-tDCS decreases cortical excitability (Nitsche et al., 2003). Although there is a general trend for anodal-tDCS to have a facilitating effect on cognitive performance, there are mixed results on effects of anodal-tDCS over DLPFC on picture naming in monolingual speakers (). These controversies resonate with meta-analyses on the effect of tDCS over DLPFC on working memory tasks, namely, improvement in reaction times but not accuracy (; ; ). In addition, a meta-analysis investigating the effect of tDCS over DLPFC on different cognitive functions (EFs, language, and memory) showed no effect of anodal-tDCS on performance (). In their review, argue that the lack of tDCS-induced behavioral effects can be due to inter-subject differences and that the lack of detailed data (especially for null-effects) makes it difficult to explain the absence of effects. They also argue that tDCS is a weak form of stimulation that might not be able to modulate the functioning of the healthy brain. Anodal-tDCS has been used over the DLPFC to modulate domain-general and language-specific cognitive control (; ; Metuki et al., 2012; ; ; ) but it has not yet been used in bilingual language processing.
According to the IC model, the DLPFC constrains performance in language production. The goal of this study is to investigate whether anodal-tDCS over the left DLPFC improves language production in bilingual speakers. Following the argument by , stating the relatively weak modulation of the healthy brain, to overcome the limitations of discrete behavioral measures, we recorded the electroencephalogram (EEG) during bilingual word production (picture naming in L1 and L2 and backward and forward translation) to assess putative effects of anodal-tDCS over DLPFC on the brain processes underlying L1 and L2 production. We tested the specific hypothesis that increasing brain activity of cognitive control areas could modify networks implicated in lexical search and selection during word production in L2 but not in L1. In addition to picture naming and translation tasks, we used a non-verbal fluency task as a measure of domain general EF, and a phonemic verbal fluency task as a measure of lexical access and EF (Shao et al., 2014). Based on the RHM, the relationship between lexica and the conceptual system depends on language proficiency which is weaker in the less proficient language (; ). According to the IC model (), less proficient bilinguals rely more on strategic control during language production in L2 (). As L2 proficiency increases, the lexico-semantic processing of L2 target words is less dependent on the activation of executive control and therefore the prefrontal cortex (Stein et al., 2009). Critically, the IC model assumes that cognitive control in language production is not domain (language) specific. Following these assumptions (; ), we therefore hypothesized that anodal-tDCS stimulation over the left DLPFC would enhance behavioral performance on both verbal and non-verbal fluency tasks.
The RHM and IC models offer complementary explanations why picture naming is more difficult in L2 than L1 and why forward translation is more difficult than backward translation. But only the IC model assumes a role of EF in this process. Hence according to the IC but not the RHM model, anodal-tDCS should improve performance that relies more on EF. Therefore, improvements should be reflected maximally in improved response accuracy on verbal fluency tasks in L2, and response accuracy and reaction time in picture naming in L2 and word translation during forward translation (which is more demanding on the control system and requires activating lexical nodes in L2 while inhibiting lexical nodes in dominant L1). At the level of EEG activity, stimulation should impact on the time windows related to executive functioning. These effects should be observed in early pre-lexical attentional processes including language switching, language selection, and target lexical selection (>250 ms post-stimulus in naming and ∼400 ms post-stimulus in translation) according to the results reported by .
Materials and Methods
Participants
We recruited native French speakers who acquired English (L2) as a second language after 7 years of age. All were all right-handed according to the Edinburgh Handedness Inventory (Oldfield, 1971) with no history of neurological, psychological, or other health problems. We recruited initially a total of 28 participants. However, for verbal and non-verbal fluency task, only 26 participants (mean age: 24.3 ± 6.4, N = 9 males) were included (one participant was excluded due to very low L2 proficiency and another was excluded due to missing data error). For the picture naming and translation tasks, 24 participants (mean age: 24.5 ± 6.5, n = 8 males) were included (three of the initially recruited participants were excluded because of missing data and one due to very low L2 proficiency).
Second Language Proficiency
Age of Acquisition and Immersion
Age of acquisition, immersion, and self-evaluation of language use was assessed using a custom-made questionnaire (Radman et al., 2016). We assessed immersion in French and English by asking the following information: the AoA, how long they lived in a region where the dominant spoken language was English or French, which language they spoke with their family members, in school, in present activities (watching TV/listening to radio, reading books, arithmetic), and if the language was acquired in school or out of school only. In the same questionnaire, participants were asked to indicate in percentage terms how well they would estimate their English reading, speaking, comprehension, and writing skills (Tables 1, 2). As expected, greater usage and immersion of French was reported compared to English.
Table 1
| Variable | Mean | SD |
|---|---|---|
| Age of acquisition (years) | 12.62 | 2.6 |
| L2 Self-evaluation (%) | ||
| Speaking | 60.5 | 15.5 |
| Comprehension | 70.5 | 13.8 |
| Reading | 75.5 | 12.2 |
| Writing | 57.6 | 20.6 |
| L2 Vocabulary Tests | ||
| DIALANG L2 vocabulary score (minimum = 0; maximum = 1000) | 656.6 | 156.6 |
Age of acquisition and L2 proficiency of participants (n = 24): using a questionnaire, participants were asked to indicate in percentage terms how well they would estimate their English reading, speaking, comprehension, and writing skills.
The vocabulary subtest from the computer-based DIALANG language diagnosis system was administered to evaluate vocabulary knowledge in English.
Table 2
| Language usage | Mean | SD | |
|---|---|---|---|
| L2 immersion and use | Use at work/studies (%) | 19.8 | 17.6 |
| TV/radio (%) | 27.2 | 23.2 | |
| With friends (%) | 11 | 12 | |
| Reading books (%) | 25 | 24 | |
| L1 immersion and use | Use at work/studies (%) | 80.21 | 17.6 |
| TV/radio (%) | 73.95 | 23.36 | |
| With friends (%) | 89.58 | 12.32 | |
| Reading books (%) | 75 | 23.93 |
L1 and L2 immersion and usage of participants (n = 24).
Computer-Based Vocabulary Knowledge Evaluation in English
The vocabulary subtest was taken from the computer-based DIALANG language diagnosis system and administered to evaluate vocabulary knowledge in English (Zhang and Thompson, 2004). In this test, participants indicate for each of 75 stimuli whether it is a correct word in English or a highly word-like pseudo-word. The mean DIALANG score in English vocabulary was 656 ± 159 in our participants, confirming an intermediate level of English vocabulary. Information on L2 skills is summarized in Table 1.
Stimuli and Task Procedure
Online Tasks
Verbal fluency task
Participants were asked to produce words starting with a specific letter within 1 min (we used letter “P” in L1 and letter “S” in L2). The selected letters have been used in previous studies with healthy and clinical populations in both languages (Paradis and Libben, 1987; Tombaugh et al., 1999). “P” and “S” are considered to be “easy letters” (). Participants were instructed to produce any word except for proper nouns, repetitions, and verb conjugations. Scoring was based on the number of words produced.
Non-verbal fluency task
Participants were given a fixed configuration of dots and asked to generate as many novel designs as possible in 120 s while avoiding repetitions and other rule-breaks (). Blind scoring was based on the number of unique designs produced and the percentage of perseverative errors [(perseverative errors/total unique designs) × 100; ].
Offline Tasks
Picture naming
Stimuli comprised two series of 70 pictures to be named in French (list 1) and English (list 2). Items were selected from the Snodgrass image corpus (Snodgrass and Vanderwart, 1980). Selected names were non-cognates. Word frequency, name agreement, and image agreement were matched across lists based on data reported by . French words were longer than English words (mean letter count 6.24 vs. 5.16, t = 3.47, p < 0.001). Image familiarity was matched across list. Image complexity [rating of the detail/intricacies of the line drawing (Snodgrass and Vanderwart, 1980)] was higher for the French list (average complexity index was 3.07 for French and 2.71 for English lists, t = 2.49, p = 0.01). All images consisted of line drawings of approximately the same size (no larger than 540 × 400 pixels) with a consistent white background. Participants were instructed to name pictures as quickly as possible. Stimuli were presented for 2000 ms at the center of a 15′ LCD screen with a refresh rate of 60 Hz. Each picture was preceded by a fixation cross for 300 ms. A white color screen of 3000 ms followed the picture. A short training was performed for the tasks during tDCS stimulation. The training consisted of series of five trials for naming and translation (different from the stimuli of the study), which were the same for all participants. The training was performed in both sessions (Figure 1B).
FIGURE 1
Word translation
Stimuli consisted of two lists of five- to eight-letter concrete nouns; List 1 contained 50 French words to be translated in English (forward translation) and List 2 contained 50 English words to be translated in French (backward translation). Items were selected from Wordgen () and CELEX () databases. None of the lists contained cognate words. The lists were matched for word length, frequency, and neighborhood size based on Wordgen and CELEX databases. Each word was presented on a 15′ LCD screen with a refresh rate of 60 Hz for 2000 ms preceded by a fixation cross of 200 ms and followed by a blank screen of 2000–3000 ms. Stimuli in both languages were presented in the same block in a random order. Participants were instructed to translate the presented word as fast as possible in the other language. No cue was given about the direction of translation (Figure 1C).
For the picture naming and translation tasks, stimulus presentation and responses were recorded using a voice activated key controlled using E-Prime 2.0 (Psychology Tools, Inc., Pittsburgh, PA, United States). In these tasks, only the first oral responses given after stimulus onset were considered. Synonymous words (in the picture naming task according to Snodgrass normative data and in word translation task according to Merriam-Webster’s dictionary) were also considered as correct responses.
Study Design
Each participant participated in two sessions presented 1 week apart. In each session, participants received 20 min of anodal or sham tDCS. Because the effect of tDCS is known to reach the maximum at around 4 min after the start of stimulation (Nitsche et al., 2003, 2007, 2008), we performed verbal phonemic (in L1 and L2) and non-verbal fluency tasks (a modified five points test, as a control task) in a controlled order (1 – verbal fluency in L1, 2 – non-verbal fluency, and 3 – verbal fluency in L2) starting 6 min before the end of stimulation (online). After anodal-tDCS stimulation, participants were asked to perform picture naming and word translation tasks while EEG was recorded simultaneously (Figure 1). The order of stimulation type was counterbalanced between the participants.
After preparing the tDCS electrode montage [active electrode over F3 and reference electrode over the right supraorbital area, cf. Section “Transcranial Direct Current Stimulation (tDCS)”], an EEG cap was placed on the head, over the tDCS electrodes. In order to minimize the time interval between the end of stimulation and the beginning of the tasks, the EEG electrodes were already attached before the stimulation (except for the EEG electrodes of the frontal part because of the overlap between tDCS and EEG electrodes). Each participant then received anodal or sham tDCS. After the end of stimulation, the tDCS electrodes were removed and the EEG montage was completed. Participants were then asked to perform three tasks while EEG was recorded. The tasks started with two blocks of picture naming [one block in L1 and one in L2 (counterbalanced order across participants)] and one block of word translation (mixed forward and backward) always in the last part of the session. The ensemble of three language tasks was designed to last maximally 25 min to ensure the stability of tDCS after effect (Nitsche et al., 2008; Figure 1A).
The protocol was first evaluated in a pilot study to confirm the feasibility of EEG recording during the stability of tDCS effects (i.e., considering the time limit of tDCS after-effect and the effect of the wet area of the tDCS on subsequent EEG recording). Our pilot study confirmed that the planned EEG recording after tDCS could indeed be conducted during the expected tDCS after-effect. EEG recording was not affected by the wet area as during EEG data preprocessing, only limited bridges between electrodes was found as a result of the wet area (see section “EEG Acquisition and Preprocessing”).
At the beginning of the first session, participants filled in the questionnaires on their health status, handedness, and their second language proficiency (AoA, immersion, and self-evaluation). They also performed a vocabulary subtest using the DIALANG language diagnosis system (for more details, see section “Second Language Proficiency”). At the end of each session, participants were asked to fill in a questionnaire on tDCS side effects that they might have noticed during that session.
This study was carried out in accordance with the recommendations of “Swiss Ethics Committees on research involving humans” with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by Cantonal Ethics Committee of Vaud. This study was also registered as a clinical trial in clinical trial registry of the U.S. National Institutes of Health1 (ID: NCT02289521).
Transcranial Direct Current Stimulation (tDCS)
Region of Interest
Stimulation was performed using DC-STIMULATOR PLUS (Eldith, neuroConn GmbH, Germany) through a pair of rubber iso-potential electrodes fully enclosed in saline-soaked sponge pocket over the left DLPFC using tDCS. We used isotonic (0.9%) NaCl saline solution. The active electrode was placed over F3 of a 10–20 EEG system. F3 has been constantly used to target DLPFC (Nitsche et al., 2009). To find the F3 accurately, we used the system introduced by and confirmed the point with an EEG cap. The return electrode was placed over the right supraorbital area.
Electrode Size
To increase the focality of stimulation under the active electrode and decrease the effect of stimulation under the return electrode, we selected different sizes of electrodes (Nitsche et al., 2007, 2008): the smaller electrode (7 cm × 5 cm) as active and the larger electrode (7 cm × 10 cm) as return electrode. The sponge–electrode–sponge thickness was 2–1.9–2 mm.
Figure 1D shows an estimate of the electric field induced by our tDCS montage in the anodal condition. The estimation of the distribution of the electric field was generated in SimNIBS 2.0.1 (Thielscher et al., 2015). The model is based on the following conductivity values for its key anatomical components [SimNIBS default values, as in, e.g., Thielscher et al., 2015 or : scalp (σ = 0.465 S/m), bone (σ = 0.010 S/m), cerebrospinal fluid (σ = 1.654 S/m), gray matter (σ = 0.275 S/m), and white matter (σ = 0.126 S/m)]. The volume mesh and visualization were generated through Gmsh ().
Duration and Current Density
A constant current of 1.5 mA was applied for 20 min (fade in: 15 s and fade out: 20 s) for anodal-tDCS resulting in a current density of 0.042 mA/cm2 under the active electrode. This stimulation duration is suggested to result in a shift in cortical excitability up to 60 min (Nitsche et al., 2008). In the sham condition, there was a short direct current of 30 s (8 s fade in and 5 s fade out) at the beginning of the stimulation. There was no direct current stimulation after the first 43 s, except for small pulses of 3 ms every 550 ms emitted by the stimulator for impedance checking. The average current over time is not more than 2 μA. This way, the participant is less able to detect any difference between sham and real stimulation (Neurocare Group, 2013).
Analysis of Behavioral Data
Online Tasks
Scores on the verbal fluency task were analyzed using a 2 × 2 within subject design with factors Language (L1; L2) and Stimulation (anodal; sham) in a two-way repeated measures ANOVA. Additionally, we have performed this analysis while controlling for the DIALANG score (as the measure of L2 proficiency), implemented in the model as a covariate. For the non-verbal fluency task, a paired t-test compared scores of the number of unique designs as well as perseverative errors between anodal and sham tDCS.
Offline Tasks
Response accuracy and voice onset time (VOT) on the picture naming task were implemented in a 2 × 2 within subject design with factors Language (L1; L2) and Stimulation (anodal; sham; using two-way repeated measures ANOVA). We have also performed this analysis while controlling for the DIALANG score (as the measure of L2 proficiency). The behavioral data from the translation task (response accuracy and VOT) were subjected to a 2 × 2 within subject design with factors Translation direction (Forward; Backward) and Stimulation (anodal; sham). We have also performed this analysis while controlling for the DIALANG score (as the measure of L2 proficiency). Behavioral analyses were performed using SPSS 21.
Comparison of tDCS Effects Across Online and Offline tasks
have put forth two explanations for the lack of behavioral effects of tDCS in cognitive tasks. They argue that (i) tDCS is not strong enough to produce an effect or that (ii) inter-individual differences obliterate the effects. If the latter is true, such inter-individual differences should be consistent across tasks. We hence assessed the consistency of tDCS across online and offline tasks and their relation to the DIALANG vocabulary score. To that end, we computed the effects of tDCS in the verbal and non-verbal fluency tasks and on the VOT and error rates for naming in L1 and L2 and backward and forward translation by subtracting the respective values in the sham from those of the anodal stimulation condition. Next, we computed correlation coefficients between those values and the DIALANG vocabulary score and corrected the so-obtained p-values using the False Discovery Rate ().
EEG Acquisition and Preprocessing
Electroencephalogram was sampled continuously at 1024 Hz from 64 preamplified Ag/AgCl electrodes using an ActiveTwo system (Biosemi, Inc., Amsterdam, Netherlands). Electrodes were placed in an elastic cap according to the extended 10–10 system, the reference electrode was placed at the vertex (“Cz”), and impedances were kept below 10 kΩ. Electrode preparation was done with SignaGel (Parker Laboratories, Inc., Fairfield, NJ, United States), in order to use the same electrolyte (NaCl) for EEG and tCDS electrodes. This was done to avoid the buildup of low frequency battery potentials that can arise from mixing different electrolytes which can affect EEG recordings. Offline, the data were recomputed to average reference and digitally band-pass filtered between 0.18 and 40 Hz using second-order Butterworth filter with a -12 db/octave roll-off. The filter was computed linearly with two passes (one forward and one backward), eliminating the phase shift, and with poles calculated each time to the desired cut-off frequency. This was done using epochs spanning -500 to +500 ms around the selected epochs. Additionally, a notch filter (50 Hz) was used to remove AC noise.
All offline analyses were performed using the Cartool software by Denis Brunet2 and the STEN toolbox developed by Jean-François Knebel.3 To avoid possible topographic distortions caused by electrodes bridges (especially because of a possibly wet skin over frontal area caused by tDCS electrodes), the eBridge toolbox in EEGLAB was used to find any possible bridges in the raw EEG files ().
Epochs from 100 ms before the presentation of the stimuli to 600 ms (for picture naming task) and to 800 ms (for word translation task) after the stimulus onset were extracted for each condition. Data were baseline corrected over the whole epoch at single epoch level. Event-related potentials (ERPs) were calculated by averaging the extracted epochs for each participant and condition separately (i.e., in picture naming task: picture naming in L1 and L2 following anodal and sham stimulation and in word translation task: forward and backward translation following anodal and sham stimulation). Only trials with correct responses were included. Epochs with eye-blink or other artifacts (as determined by amplitude changes exceeding 80 μV on at least one electrode during an epoch) were rejected before averaging (). Since longer epochs were extracted for the translation task, more epochs were contaminated with eye-blinks. This resulted in too many trials being removed because of eye-blink. Therefore, for the translation task, independent component analysis (ICA) was implemented using the extended runica algorithm in EEGLAB to remove eye-blinks. This led to increase the number of trials included in the ERP. Before group averaging, electrodes exhibiting substantial artifacts as well as bridged electrodes from each participant were interpolated using a three-dimensional spline algorithm before statistical analyses (mean 6.25% interpolated electrodes; Perrin et al., 1987). There was on average 1.25 (±1.3) EEG sensor bridges per participant. Nine subjects showed no EEG sensor bridges in both the two sessions.
Electrical Neuroimaging Data Analysis
Analysis of ERP Waveforms
Electrode and time-point wise analysis of the ERPs was conducted by comparing the ERPs to L1/Anodal, L1/Sham, L2/Anodal, and L2/Sham for picture naming using a 2 × 2 within-subjects ANOVA with factors Language and Stimulation for word translation at each scalp electrode as a function of peri-stimulus time. The same analysis was performed for ERPs to Forward/Anodal; Forward/Sham; Backward/Anodal; and Backward/Sham conditions for word translation with factors Translation direction (Forward; Backward) × Stimulation (Anodal; Sham). This data-driven analysis allows us to identify peri-stimulus time windows showing sustained effects without being limited to a priori selected time windows (e.g., specific ERP components). Correction was made for temporal and spatial auto-correlation through the application of a >11 contiguous data points (10 ms) on at least 10% of electrodes criterion for the persistence of significant effects (). The analysis was performed using the STEN toolbox.
For the periods of interaction between the studied factors, to identify the direction of the effects, sources estimations were computed and statistically analyzed. In fact, in contrast to the source analyses, waveform analyses depend on the selection of the electrode. The direction of the effects found at waveform level is not necessarily similar to the results obtained from statistical analyses on source level. We therefore interpret the results based on the statistical analyses at source level. The details on intracranial source estimation and analyses come in the next part.
Analysis of Intracranial Sources and Source Differences
The estimation of intracranial generators for a given scalp topography is an ill-posed problem, because a given scalp topography can in principle be generated by any combination of intracranial sources (inverse problem). This inverse problem can be overcome by implementing known biophysical constraints of the generation and propagation of intracranial sources. We used a distributed linear inverse solution based on local autoregressive average (LAURA) regularization approach (; Michel et al., 2004). LAURA selects the source configuration that best mimics the biophysical behavior of the electric fields; it confines the solution space in the gray matter of the brain and takes into account how the activity in a given area diminishes with the distance from the scalp and thereby assumes smoothness between neighboring sources. The solution space was based on a simplified realistic head model (SMAC; Spinelli et al., 2000) and comprised 3005 points, selected from a 6 mm × 6 mm × 6 mm grid, homogeneously distributed within the gray matter of the average brain of the Montreal Neurological Institute (MNI) that was used as template source space for all subjects. The lead field (or the forward solution) was then solved with an analytical solution with a three-shell spherical head model (brain, skull, and scalp). Numerous experimental and clinical studies have demonstrated that LAURA yields reliable source estimates (; ; Plomp et al., 2010; ).
Notwithstanding this careful implementation to solve the inverse problem, one is confronted with the problem of thresholding of intracranial sources. There is not and there cannot be a predefined criterion according to which an estimated source can considered as “active.” One way of overcoming this problem of considering absolute estimates of intracranial current density is to perform statistical comparisons between conditions in the source space (; ).
The sources of the ERP were calculated for each subject and each condition previously averaged over the period of interest (i.e., the period showing the significant ERP waveform modulation) and then statistically compared using the same within-subject design as the waveform analysis. To control for multiple comparisons, only significant clusters with a minimal size of 14 consecutive points (KE) were retained. This spatial criterion was determined with the AlphaSim program.4 There was a false positive probability of p < 0.005 for observing a cluster of >14 contiguous nodes (see also ; for the same approach).
Results
Tolerance of tDCS Stimulation
We have evaluated the effects of the stimulation based on a questionnaire proposed by . All of the participants tolerated the stimulation well. Participants reported no different effect of scalp tingling sensation at the stimulation site following both anodal and sham tDCS (Mann–Whitney U-test, p > 0.05). No participant reported adverse effects or asked to interrupt the experiment.
Behavioral Results
Online Tasks
Verbal fluency
Figure 2A depicts the results of the verbal fluency task. We found a main effect of Language [F(1,25) = 87.1, p < 0.001, = 0.78], with more words generated in L1 than L2. However, there was no effect of Stimulation [F(1,25) = 1.6, p = 0.21, = 0.061] nor interaction between Language and Stimulation [F(1,25) = 0.13, = 0.005]. After controlling for L2 proficiency, the main effect of Language [F(1,25) = 5.17, p = 0.03, = 0.17 remained significant with more words were generated in L1 than L2]. However, there was no effect of Stimulation [F(1,25) = 1.27, p = 0.27, = 0.05] nor interaction between Language and Stimulation [F(1,25) = 1.23, p = 0.78, = 0.049].
FIGURE 2
Non-verbal fluency
Figure 2B depicts the results of the non-verbal fluency task: there was a higher number of unique designs and a lower percentage of perseverative errors after anodal than sham tDCS [t(25) = 2.76, p = 0.011 and t(25) = 2.70, p = 0.012, respectively].
Offline Tasks
Picture naming
Figure 2C depicts the behavioral results of the picture naming task. We found a main effect of Language [F(1,23) = 220.75, p < 0.001, = 0.90] with better performance in L1 than L2. However, there was no effect of Stimulation [F(1,23) = 0.46, = 0.02] nor an interaction between Language and Stimulation [F(1,23) = 0.03, = 0.001; Figure 2C, 1]. After controlling for the effect of L2 proficiency, the effect of Language [F(1,23) = 24.75, p < 0.001, = 0.52] remained significant with better performance in L1 than L2, but there was no effect of Stimulation [F(1,23) = 2.45, = 0.1] nor an interaction between Language and Stimulation [F(1,23) = 0.02, = 0.001].
For naming latency (VOT), there was an effect of Language [F(1,23) = 70.4, p < 0.001, = 0.75] with faster latencies in L1 than L2, but no effect of Stimulation [F(1,23) < 0.001, = 0.001] nor an interaction between Language and Stimulation [F(1,23) = 0.46, = 0.02; Figure 2C, 2]. After controlling for the effect of L2 proficiency, the main effect of Language [F(1,23) = 10.99, p = 0.003, = 0.33] remained significant, but there was no effect of Stimulation [F(1,23) = 0.04, = 0.002] nor an interaction between Language and Stimulation [F(1,23) = 0.172, = 0.008].
Word translation
Figure 2D depicts the behavioral results of the Translation task. For translation performance, there was a main effect of Translation direction [F(1,23) = 31.2, p < 0.001, = 0.58] with backward translation better than forward translation. However, there was no effect of Stimulation [F(1,23) = 0.59, = 0.02] nor an interaction between Translation and Stimulation [F(1,23) = 2.13, p = 0.16, = 0.08; Figure 2D, 1]. After controlling for the effect of L2 proficiency, the main effect of Translation direction no longer reached significance [F(1,23) = 3.82, p = 0.063, = 0.15], and neither the main effect of stimulation, nor the interaction between these factors were significant [F(1,23) = 0.19, = 0.009 and F(1,23) = 0.07, = 0.003, respectively].
For translation latency (VOT), there was an effect of Translation direction [F(1,23) = 32.7, p < 0.001, = 0.58] with backward translation faster than forward. However, there were no main effects of Stimulation [F(1,23) = 0.047, = 0.002] nor significant interaction between Language and Stimulation [F(1,23) = 0.87, = 0.03; Figure 2D, 2]. After controlling for L2 proficiency, the main effect of translation direction failed to reach significance [F(1,23) = 3.78, p = 0.06, = 0.147], and there was no main effect of Stimulation nor a significant interaction between Translation direction and Stimulation [F(1,23) = 1.08, = 0.047 and F(1,23) = 0.046, = 0.002, respectively]. We should add that because the data on accuracy in picture naming and translation did not follow normal distribution in some conditions, supplementary analyses were performed on rank transformed picture naming and translation data. This information can be found in Supplementary Material, Presentation 1.
Consistency of tDCS Effect Across Tasks
Even though there was no overt effect of Stimulation except in the online non-verbal fluency task, tDCS had consistent effects across all tasks, both online and offline: in some subjects, it had facilitating effects and in others it had inhibitory effects, and it did so consistently across all tasks. We computed correlation coefficients between the behavioral measures in the six tasks we employed and the DIALANG score and adjusted the significance level using the false discovery rate; the adjusted p-value is 0.0110. These results are depicted in Figure 3 (Supplementary Figure 1 shows the effects in all tasks for each subject). tDCS had consistent effects in all online tasks, the only ones that did not correlate significantly were verbal fluency in L1 and L2 (see also Supplementary Figure 1). Offline, all measures of VOT (naming in L1 and L2, forward and backward translation) were highly and significantly correlated (Supplementary Figure 1), i.e., subjects were consistently faster or slower in pronouncing words in L1 and L2 after stimulation. A similar effect was obtained for errors: errors committed for naming in L2 were significantly correlated with errors in both forward and backward translation (Supplementary Figure 1), and VOT in forward translation was significantly correlated with errors in naming in L2 and forward translation. Finally, the effects of tDCS in the online non-verbal fluency task was significantly correlated with the online verbal fluency task and with errors committed in all offline tasks as well as with the VOT in naming in L2 and forward translation (the two most difficult conditions, in which words had to be produced in L2).
FIGURE 3
It is important to note here that the DIALANG vocabulary score did not correlate with any of the behavioral measures, i.e., the consistent effects of tDCS on word production cannot be attributed to language performance measured with the DIALANG vocabulary score.
Electrical Neuroimaging Results
ERP Waveforms
Picture naming
A time-wise 2 × 2 voltage waveform analysis on all electrodes with the factors Language and Stimulation showed main effects of Language at 110–150 ms on left central and parietal electrodes and 235–380 ms post-stimulus presentation on fronto-central electrodes. There was also an interaction between Language and Stimulation at 365–410 ms on left and posterior electrodes [F(1,23) = 15.8, p < 0.001; >11TF]. There was no effect of stimulation on picture naming. The topographic maps for group mean of the ERPs in all conditions showed a different map for picture naming performance in L2 following anodal stimulation when compared to all other conditions (Figures 4A,B).
FIGURE 4
Word translation
A time-wise 2 × 2 voltage waveform analysis on all electrodes with Translation direction and Stimulation showed a main effect of Translation direction at 350–600 ms on bilateral fronto-parietal electrodes and 675–750 ms post-stimulus presentation on left and right parietal electrodes and an interaction between Translation direction and Stimulation at 445–610 ms on left frontal electrodes [F(1,23) = 16.6, p < 0.001; >11TF]. There was no main effect of stimulation. Visual inspection of the topographic maps for group mean of the ERPs in different conditions suggests a different map for backward translation after anodal stimulation compared to all other conditions (Figures 5A,B).
FIGURE 5
Electrical Source Estimation
Repeated measures ANOVA of distributed source estimates for all conditions were performed for each of the 3005 solution points for the period of significant effects identified in voltage waveform analyses, i.e., 365–410 ms for picture naming and 445–610 ms for translation.
Picture naming
Analysis revealed that the ERP modulation at around 365–410 ms post-stimulus presentation originated from an interaction between Language and Stimulation [F(1,23) = 4.5, p < 0.05, KE > 14 nodes] in the left inferior parietal lobule (IPL; BA 40). This interaction was driven by an increase of activity within the left IPL during picture naming in L2 after anodal stimulation relative to sham stimulation (). There were no effects of anodal-tDCS on picture naming in L1 (depicted in Figure 4C).
Word translation
The ERP modulation at 445–610 ms post-stimulus presentation originated from interactions between Translation direction and Stimulation [F(1,23) = 4.3, p < 0.05, KE > 14 nodes] in the right anterior medial prefrontal cortex (AMPFC; BA10). This interaction was driven by an increase of activity within the right AMPFC during backward translation after anodal stimulation relative to sham stimulation (p < 0.01) and the increase of activity within this area following anodal stimulation in backward but not forward translation (p = 0.01). There were no effects of anodal-tDCS on forward translation (depicted in Figure 5C).
Discussion
We investigated the role of EFs in bilingual language production by activating DLPFC by means of anodal-tDCS. To that end, we investigated behavioral performance and ERPs and their sources during picture naming in L1 and L2 and backward and forward translation during the steady tDCS aftereffect. In addition, we investigated verbal and non-verbal fluency during anodal-tDCS stimulation of DLPFC.
According to the neurocognitive model of bilingual language processing (), it was expected that stimulation of DLPFC using anodal-tDCS would improve fluency measures on a range of tasks. The only effect of anodal-tDCS we found was an improvement in non-verbal fluency during stimulation. Contrary to our expectations, we found no behavioral improvements on verbal fluency during stimulation and no effect on picture naming or word translation tasks following translation. This is in line with evidence summarized in a recent review that shows no effect of tDCS on cognitive performance () in healthy individuals. The authors of this study argue that the effects of tDCS are either to weak or obliterated by inter-individual differences. We tackled this issue by assessing the consistency of the effects of tDCS across the six behavioral tasks we used in our study. Even though tDCS does not improve word production in L1 and L2, we can show that it has very consistent effects across all behavioral measures, both online and offline. In a subset of subjects, it has consistently facilitating effects (more drawings produced in the non-verbal fluency task, more words produced in L1 and L2, faster VOT, and less errors in naming in L2 and both forward and backward translation), whereas it has consistently impeding effects in others (less drawings produced in the non-verbal fluency task, less words produced in L1 and L2, slower VOT, and more errors in naming in L2 and both forward and backward translation). Importantly, tDCS only affected L2 production, but not L1 production (apart from VOT). Inter-individual differences in the consistent effect of tDCS have obliterated group effects in the present study. Surprisingly, these inter-individual differences cannot be attributed to L2 performance, at least not to L2 performance measured with the vocabulary subtest of the DIALANG test. An important lesson to be learned from this is that tDCS can have consistently facilitating or inhibitory effects across multiple tasks that are obliterated by inter-individual differences. A solution to this issue is to attempt to elucidate the nature of those consistent effects.
Our EEG source imaging results draw a more differentiated picture of the effects of anodal-tDCS on bilingual language production: we found evidence that anodal-tDCS had a larger effect on L2 picture naming than L1 picture naming at an epoch related to executive functioning between 365 and 410 ms post-stimulus presentation. This was reflected in increased activity in the left IPL (BA 40) after anodal relative to sham stimulation in L2. This is in line with a recent study highlighting the importance of this area in L2 acquisition and performance (). In a 12-week intensive language training, activity in this area prior to training predicted success of the training, and activity in area after training correlated with L2 performance.
We also found evidence that anodal-tDCS had a larger effect on backward translation than forward translation at an epoch that is related to phonological retrieval between 445 and 610 ms post-stimulus presentation. This was reflected by increased activity in the right AMPFC (BA 10) in anodal relative to sham stimulation in backward translation. In broad terms, these results are consistent with the results reported by who found larger N400 amplitudes for backward translation.
In the picture naming task, we found a main effect of language on response accuracy and VOT, but no main effect of stimulation and no interaction between language and stimulation. The effect of language is compatible with previous studies (faster and more accurate in L1 than L2; e.g., ; ). The EEG results in the picture naming task help to refine this view; even though we do not find an interaction between stimulation and language, we do find this interaction in the ERPs and their concomitant source differences: the effect of anodal-tDCS was bigger than sham stimulation in L2, and no such effects were found in L1. The present findings therefore suggest to us that anodal-tDCS has a specific effect on access to word forms (lexemes) and word form encoding in L2. Note that L2 is assumed to rely more on executive processes that may be reflected here. By means of EEG source localization (LAURA), we determined that the facilitation of L2 picture naming at around 365–410 ms post-stimulus onset is driven by increased activity within the left IPL. Even though we did not find a behavioral effect of anodal-tDCS over DLPFC on picture naming in L2, we can show that stimulating DLPFC recruits left IPL exclusively for L2 production after stimulation in a time window (365–410 ms) important for form encoding (phonological code retrieval, phonological encoding, and syllabification; ; ). While the left IPL is involved in identification and discrimination of syllables (), suggest a specific role for left IPL in bilingual language control, including switching and inhibitory functions. IPL is also known for the integration of semantic information in both language production and comprehension tasks ().
In the translation task, contrary to our prediction, we did not find a main effect of stimulation nor did we find an interaction between stimulation and translation direction. Like in the picture naming task, our EEG source imaging results shed a more differentiated light onto the underlying brain processes: we found an interaction between translation direction and stimulation in the time window 445–610 ms post-stimulus onset in the ERPs and their concomitant sources. We found increased activity in AMPFC after anodal than sham stimulation in the forward translation condition. AMPFC is recruited more in attention demanding processes as a part of a network guiding cognitive performance (Zysset et al., 2003). used PET and showed that the pattern of brain activation in complex language tasks (metaphor comprehension) engages the analog regions in the right hemisphere including right prefrontal cortex.
Our ERP results in the translation task are in line with the results reported by who found that neural processing differed according to the direction of translation at around 200 and 400 ms following stimulus presentation with larger P2 amplitudes for forward translation and larger N400 amplitudes for backward translation. They explained their results as due to the inversed language effect (Meuter and Allport, 1999), which has been interpreted in terms of reduced access to L1 or the inhibition of L1. This may be the result of adaptation to mixed-language situations as a consequence of sustained language control. Translation is a mixed-language situation where access to L1 words is inhibited to facilitate L2 processing, which will result in a higher activation threshold for L1 than for L2. An alternative explanation however is more effortful conceptual processing in backward compared to forward translation. For example, Phillips et al. (2006) report phonological mismatch negativity (PMN) which occurs around 250–350 ms post-stimulus presentation and N400 components amplitude changes in between-language compared to within-language repetitions with English–French bilinguals. The authors argued that these effects (i.e., larger PMN and N400 to the forward translation and larger PMN but minimal N400 activity to the backward translation) indicate a more strongly engaged conceptual processing in backward than forward translation.
The slower RTs in forward relative to backward the translation can be explained in terms of the RHM model () by assuming that forward translation in less proficient L2 speakers is mediated via the L1 translation equivalent, which then leads to slower responses. In terms of the IC Model (), the results confirm that forward translation is slower and more error prone and associated with inhibition of dominant L1 lexical nodes – and this will be more difficult.
Our behavioral data confirmed longer latencies in the translation compared to the picture naming tasks (around 200 ms longer) suggesting additional processing. This might be related to target lexical selection after semantic processing or control necessary in switching from the input language to the target language. Moreover, retrieval of corresponding phonological codes may be more complex in translation compared to picture naming possibly leading to more processing steps. Thus, although the stimulation effect is seen in time windows that are different to picture naming results, the time period associated with interactions between translation direction and stimulation may also reflect phonological retrieval and form encoding.
Anodal stimulation over DLPFC improved non-verbal fluency performance, i.e., it improved the number of unique designs generated and decreased perseverative errors. Non-verbal fluency is a task that is sensitive to changing demands on executive functioning (). Non-verbal fluency relies on problem-solving to produce unique designs. Robins Wahlin et al. (2015) reported that non-verbal fluency tasks are more sensitive compared to verbal fluency tasks to evaluate cognitive decline in prodromal Huntington disease. The authors suggested that phonemic fluency requires retrieving items stored in the individual’s lexicon, while non-verbal fluency requires imagination of unknown designs not existing in any lexicon. They assume that the less constrained nature of the task may require more EFs. Interestingly, positive effects of neuromodulation on non-verbal fluency are not reported to date. The presence of the effect of anodal-tDCS on non-verbal fluency task can be also explained by performing this task online during tDCS. Previous studies have indeed shown that this approach may be more efficient in inducing behavioral effects as it controls for the changes in strength and consistency of tDCS aftereffects (Romei et al., 2016). One reason for null effects on the verbal fluency may a relatively weak contribution of DLPFC in language-related EF should also be considered. suggested a preferential role of basal ganglia and dorso-lateral temporal areas in phonemic verbal fluency. Because basal ganglia are connected to cortical structures such as DLPFC via the ganglia–thalamo-cortical loops, fluency impairment following disruption in ganglia–thalamo-cortical loops has been demonstrated (). Our results may be explained by neural specificity and type of fluency tested (left-phoneme generation and right-design generation, e.g., Ruff et al., 1994); left prefrontal cortex has also been shown to influence figure fluency performance.
The observed pattern of the behavioral findings is not in line with the IC model. That is, our behavioral results do not support the role of (DLPFC-mediated) EF in the processing of L2. It is notable that according to the IC model, several brain areas including DLPFC, ACC, and basal ganglia mediate EF in processing language in bilingual brain. Therefore, we are not able to totally rule out the IC model as we have only considered DLPFC mediated executive functioning. However, because the RHM model does not assume the importance of the EF in bilingual language processing, our behavioral findings could be explained by the RHM model.
The absence of behavioral results can be due the weak relative contribution of DLPFC-mediated control processes involved in language production processes even in L2. The lack of sensitivity of discrete behavioral measures to track subtle tDCS-induced changes in brain activity could explain the absence of behavioral results. The selection of stimulation zone can also explain the absence of behavioral improvement in the translation task. Notably, stimulating areas with more language-specific executive control effects [e.g., anterior cingulate cortex (ACC), which has been suggested to be involved in conflict monitoring and response inhibition, especially in translation tasks] could have been more advantageous. However, ACC is rather difficult to be reached and stimulated using transcranial stimulation techniques. Alternatively, the absence of behavioral results in presence of ERP findings can be explained also by applying different neurocognitive strategies to perform the same task (e.g., picture naming). A further alternative explanation would be that, tDCS stimulation in the present study did not produce enough change at neuronal level to become behaviorally measurable (e.g., ). However, the applied stimulation montage and parameters have been selected following the recommendations of the literature on tDCS stimulation over DLPFC area.
To our knowledge, this is the first study that used the combination of anodal-tDCS over DLPFC and electrical source imaging to assess the role of EF in bilingual language production (picture naming and translation). While we can reliably reproduce behavioral effects of language [picture naming is faster and more accurate in L1 than L2, and backward translation (L2 → L1) is faster and more accurate than forward translation (L1 → L2)], we did not find a main effect of stimulation of DLPFC on bilingual language production nor an interaction between language and stimulation. At first sight, this absence of a stimulation effect is in line with many behavioral studies that failed to find effects of tDCS on cognitive performance (). However, absence of evidence does not automatically imply evidence of absence. We added a continuous direct measure of brain activity with a high temporal resolution to the discrete behavioral measures which allowed us to elucidate the effects of stimulating DLPFC on bilingual language production: in the picture naming task, we found an interaction between stimulation and language and could show that IPL, an area crucially involved in L2 acquisition and performance is more strongly recruited in L2 after anodal than sham stimulation in a time window related to lexical form encoding. Likewise, in the translation task, we found an interaction between stimulation and translation direction that was absent in the behavioral data: we can show that AMPFC is recruited more strongly during backward than forward translation after anodal than sham stimulation. Continuous measures of brain activity such as EEG and ERPs provide valuable information about how the brain accomplishes a task, and it is important to note that such brain mechanisms do not always have a concomitant behavioral consequence: for example, during shifting and maintenance of visual–spatial attention, alpha power is retinotopically suppressed and enhanced which facilitates the suppression of irrelevant distractors (Rihs et al., 2007, 2009). The fact that this important mechanism underlying the governing of attentional selection has no direct behavioral consequence (alpha power is not related to reaction times and accuracy) does not diminish its relevance. In the language domain, it has been shown that high- and low-working memory span readers show no difference in the comprehension of simple sentences with a thematic violation. However, they differently use animacy and world knowledge to solve the thematic violation which manifests by an N400 for low-span readers and a P600 for high-span readers. Hence, even though there is no behavioral difference in comprehension between the two groups, ERP results provide important information of how comprehension is differentially achieved (Nakano et al., 2010).
In the present study, discrete behavioral measures such as reaction time and accuracy might not be sensitive enough to track the subtle differences induced by anodal-tDCS on L2 production, but our ERP and ERP source imaging results provide more detailed information of how tDCS can affect word production in a second language.
We note that our study has several limitations that encourage further investigation. First, combining tDCS and EEG is and remains a technical challenge and requires solving the trade-off between the rather short-lived tDCS after-effect (∼30 min) and proper electrode preparation. Because it is not possible to completely remove the saline solution used to prepare the tDCS electrodes and to then properly prepare the 64 scalp EEG electrodes, we decided to prepare the EEG montage before the tDCS montage and to then put the remaining EEG electrodes in place. This is not ideal but under the given circumstances the best solution. We used the same electrolyte (NaCl) for both tDCS and EEG electrode preparation to avoid battery potentials and we checked for potential electrode bridges. We can show that the P100 component evoked by the visual stimuli has the canonical waveform, scalp distribution and intracranial sources, which indicates that the tDCS montage did not fundamentally alter the electric field (Supplementary Figure 2). Second, the stimulation would have stronger effect if the target area was selected among the areas which are more specifically involved in language-related EF (e.g., ACC and IPL), although those areas are more difficult to target using tDCS. Third, because of the non-focal nature of tDCS stimulation, one could argue that the observed effect could be due to the effect of stimulation on other cortical areas rather than DLPFC (e.g., the area under the return electrode). This could be investigated using control tDCS montages to compare with used montage in the present study.
Overall, we can learn two positive lessons from the negative behavioral result in our study on anodal-tDCS on word production in the native and a second language: first, it is not the case that anodal-tDCS does not affect word production at all; it rather has strong facilitatory effects in some subjects and strong inhibitory effects in others which is obliterated by considering the effect of stimulation across all subjects. More importantly so, we can show that it has very consistent effects in the six behavioral tasks (two online: verbal and non-verbal fluency, four offline: naming in L1 and L2 and translating into L1 and L2) using two behavioral measures (reaction times and accuracy) in our study. The challenge in future studies is to elucidate the source of such consistent effects. In our case, it would be obvious to attribute these differences to L2 performance, which is not the case: the scores on the DIALANG vocabulary comprehension test do not predict the production of words in a foreign language. Second, using a continuous measure of brain activity with high temporal resolution such as EEG mapping and source analysis appears to be a promising new avenue to study the subtle effects of tDCS on brain and behavior. Before we can make the claim that tDCS has no overt behavioral effect in a task, one should consider the consistency of tDCS across tasks and try to elucidate the nature of inter-individual differences that obliterate its effects on the group level.
Statements
Author contributions
NR contributed to conception, design, acquisition of data, data analyses, and interpretation of the results and writing of the manuscript. JB contributed to data analyses, interpretation of the results and writing the manuscript. BW contributed to the interpretation of the results and writing the manuscript. KB and J-MA contributed to conception, design, and interpretation of the results. LS contributed to conception of the study and interpretation of the results.
Acknowledgments
This work was supported by the Swiss National Science Foundation (#325130_156937 to J-MA and #32003B_156854 and #316030_144998 to LS). The STEN toolbox (http://doi.org/10.5281/zenodo.1164038) has been programmed by Jean-François Knebel and Michael Notter, from the Laboratory for Investigative Neurophysiology (the LINE), Lausanne, Switzerland, and is supported by the Center for Biomedical Imaging (CIBM) of Geneva and Lausanne and by National Center of Competence in Research project “SYNAPSY – The Synaptic Bases of Mental Disease”; project no. 51AU40_125759. The Cartool software (cartoolcommunity.unige.ch) has been programmed by Denis Brunet, from the Functional Brain Mapping Laboratory (FBMLab), Geneva, Switzerland, and is supported by the Center for Biomedical Imaging (CIBM) of Geneva and Lausanne.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2018.00490/full#supplementary-material
References
1
AbutalebiJ.AnnoniJ. M.ZimineI.PegnaA. J.SeghierM. L.Lee-JahnkeH.et al (2008). Language control and lexical competition in bilinguals: an event-related FMRI study.Cereb. Cortex181496–1505. 10.1093/cercor/bhm182
2
AbutalebiJ.GreenD. (2007). Bilingual language production: the neurocognition of language representation and control.J. Neurolinguistics20242–275. 10.1016/j.jneuroling.2006.10.003
3
AbutalebiJ.GreenD. (2016). Neuroimaging of language control in bilinguals: neural adaptation and reserve.Biling. Lang. Cogn.19689–698. 10.1017/S1366728916000225
4
AlarioF.-X.FerrandL. (1999). A set of 400 pictures standardized for French: norms for name agreement, image agreement, familiarity, visual complexity, image variability, and age of acquisition.Behav. Res. Methods Instrum. Comput.31531–552. 10.3758/BF03200732
5
AlschulerD. M.TenkeC. E.BruderG. E.KayserJ. (2014). Identifying electrode bridging from electrical distance distributions: a survey of publicly-available EEG data using a new method.Clin. Neurophysiol.125484–490. 10.1016/j.clinph.2013.08.024
6
BaayenR. H.PiepenbrockR.GulikersL. (1996). The CELEX Lexical Database (Release 2) [CD-ROM].Philadelphia, PA: Data Consortium, Univ. Pennsylvania.
7
BarbeauE.ChaiX.ChenJ.SolesJ.BerkenJ.BaumS.et al (2017). The role of the left inferior parietal lobule in second language learning: an intensive language training fMRI study.Neuropsychologia98169–176. 10.1016/j.neuropsychologia.2016.10.003
8
BeamW.BorckardtJ. J.ReevesS. T.GeorgeM. S. (2009). An efficient and accurate new method for locating the F3 position for prefrontal TMS applications.Brain Stimul.250–54. 10.1016/j.brs.2008.09.006
9
BenjaminiY.HochbergY. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing.J. R. Stat. Soc.57289–300.
10
BorkowskiJ. G.BentonA. L.SpreenO. (1967). Word fluency and brain damage.Neuropsychologia5135–140. 10.1016/0028-3932(67)90015-2
11
BottiniG.CorcoranR.SterziR.PaulesuE.SchenoneP.ScarpaP.et al (1994). The role of the right hemisphere in the interpretation of figurative aspects of language. A positron emission tomography activation study.Brain1171241–1253. 10.1093/brain/117.6.1241
12
BotvinickM.CohenJ. D.CarterC. S. (2004). Conflict monitoring and anterior cingulate cortex: an update.Trends Cogn. Sci.8539–546. 10.1016/j.tics.2004.10.003
13
BritzJ.PittsM. A. (2011). Perceptual reversals during binocular rivalry: ERP components and their concomitant source differences.Psychophysiology481490–1499. 10.1111/j.1469-8986.2011.01222.x
14
BrodbeckV.SpinelliL.LascanoA. M.WissmeierM.VargasM.-I.VulliemozS.et al (2011). Electroencephalographic source imaging: a prospective study of 152 operated epileptic patients.Brain1342887–2897. 10.1093/brain/awr243
15
BrunoniA. R.AmaderaJ.BerbelB.VolzM. S.RizzerioB. G.FregniF. (2011). A systematic review on reporting and assessment of adverse effects associated with transcranial direct current stimulation.Int. J. Neuropsychopharmacol.141133–1145. 10.1017/S1461145710001690
16
BrunoniA. R.VanderhasseltM. A. (2014). Working memory improvement with non-invasive brain stimulation of the dorsolateral prefrontal cortex: a systematic review and meta-analysis.Brain Cogn.861–9. 10.1016/j.bandc.2014.01.008
17
CattelaniR.Dal SassoF.CorsiniD.PosteraroL. (2011). The modified five-point test: normative data for a sample of Italian healthy adults aged 16-60.Neurol. Sci.32595–601. 10.1007/s10072-011-0489-4
18
ChouiterL.HolmbergJ.ManuelA. L.ColomboF.ClarkeS.AnnoniJ. M.et al (2016). Partly segregated cortico-subcortical pathways support phonologic and semantic verbal fluency: a lesion study.Neuroscience329275–283. 10.1016/j.neuroscience.2016.05.029
19
ChristoffelsI. K.GanushchakL.KoesterD. (2013). Language conflict in translation: an ERP study of translation production.J. Cogn. Psychol.25646–664. 10.1080/20445911.2013.821127
20
CoelhoC. (1995). Discourse production deficits following traumatic brain injury a critical review of the recent literature.Aphasiology9409–429. 10.1080/02687039508248707
21
CunilleraT.BrignaniD.CucurellD.FuentemillaL.MiniussiC. (2016). The right inferior frontal cortex in response inhibition: a tDCS-ERP co-registration study.Neuroimage14066–75. 10.1016/j.neuroimage.2015.11.044
22
de LaraG.AlekseichukI.TuriZ.LehrA.AntalA.PaulusW. (2017). Perturbation of theta-gamma coupling at the temporal lobe hinders verbal declarative memory.Brain Stimul.1131023–31029. 10.1016/j.brs.2017.12.007
23
De LuciaM.ClarkeS.MurrayM. M. (2010). A temporal hierarchy for conspecific vocalization discrimination in humans.J. Neurosci.3011210–11221. 10.1523/JNEUROSCI.2239-10.2010
24
DedonckerJ.BrunoniA. R.BaekenC.VanderhasseltM. A. (2016). A systematic review and meta-analysis of the effects of transcranial direct current stimulation (tDCS) over the dorsolateral prefrontal cortex in healthy and neuropsychiatric samples: influence of stimulation parameters.Brain Stimul.9501–517. 10.1016/j.brs.2016.04.006
25
DewingK.BattyeG. (1971). Attention deployment and nonverbal fluency.J. Pers. Soc. Psychol.17214–218. 10.1037/h0030307
26
DunnO. J. (1961). Multiple comparisons among means.J. Am. Stat. Assoc.5652–64. 10.1080/01621459.1961.10482090
27
DuyckW.DesmetT.VerbekeL.BrysbaertM. (2004). WordGen: a tool for word selection and nonword generation in Dutch, English, German, and French.Behav. Res. Methods Instrum. Comput.36488–499. 10.3758/BF03195595
28
FabbroF. (2001). The bilingual brain: bilingual aphasia.Brain Lang.79201–210. 10.1006/brln.2001.2480
29
FernandezA. L.MoroniM. A.CarranzaJ. M.FabbroN.LebowitzB. K. (2009). Reliability of the five-point test.Clin. Neuropsychol.23501–509. 10.1080/13854040802279675
30
FertonaniA.RosiniS.CotelliM.RossiniP. M.MiniussiC. (2010). Naming facilitation induced by transcranial direct current stimulation.Behav. Brain Res.208311–318. 10.1016/j.bbr.2009.10.030
31
GeranmayehF.LeechR.WiseR. J. (2015). Semantic retrieval during overt picture description: left anterior temporal or the parietal lobe?Neuropsychologia76125–135. 10.1016/j.neuropsychologia.2014.12.012
32
GeuzaineC.RemacleJ. (2009). Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities.Int. J. Numer. Methods Eng.791309–1331. 10.1002/nme.2579
33
GoralM.LevyE. S.KastlR. (2007). Cross-language treatment generalisation: a case of trilingual aphasia.Aphasiology103203–204. 10.1016/j.bandl.2007.07.116
34
Grave de Peralta MenendezR.Gonzalez AndinoS.LantzG.MichelC. M.LandisT. (2001). Noninvasive localization of electromagnetic epileptic activity. I. Method descriptions and simulations.Brain Topogr.14131–137. 10.1023/A:1012944913650
35
GreenD. (1998). Mental control of the bilingual lexico-semantic system.Biling. Lang. Cogn.167–81. 10.1017/S1366728998000133
36
GrosjeanF. (2010). Bilingual: Life and Reality.Cambridge, MA: Harvard University Press. 10.4159/9780674056459
37
GuthrieD.BuchwaldJ. S. (1991). Significance testing of difference potentials.Psychophysiology28240–244. 10.1111/j.1469-8986.1991.tb00417.x
38
HernandezA.LiP.MacWhinneyB. (2005). The emergence of competing modules in bilingualism.Trends Cogn. Sci.9220–225. 10.1016/j.tics.2005.03.003
39
HickokG.PoeppelD. (2007). The cortical organization of speech processing.Nat. Rev. Neurosci.8393–402. 10.1038/nrn2113
40
HillA. T.FitzgeraldP. B.HoyK. E. (2016). Effects of anodal transcranial direct current stimulation on working memory: a systematic review and meta-analysis of findings from healthy and neuropsychiatric populations.Brain Stimul.9197–208. 10.1016/j.brs.2015.10.006
41
HorvathJ. C.ForteJ. D.CarterO. (2015). Quantitative review finds no evidence of cognitive effects in healthy populations from single-session transcranial direct current stimulation (tDCS).Brain Stimul.8535–550. 10.1016/j.brs.2015.01.400
42
HusseyE. K.WardN.ChristiansonK.KramerA. F. (2015). Language and memory improvements following tDCS of left lateral prefrontal cortex.PLoS One10:e0141417. 10.1371/journal.pone.0141417
43
IndefreyP. (2011). The spatial and temporal signatures of word production components: a critical update.Front. Psychol.2:255. 10.3389/fpsyg.2011.00255
44
IndefreyP.LeveltW. J. (2004). The spatial and temporal signatures of word production components.Cognition92101–144. 10.1016/j.cognition.2002.06.001
45
JamesC. E.BritzJ.VuilleumierP.HauertC. A.MichelC. M. (2008). Early neuronal responses in right limbic structures mediate harmony incongruity processing in musical experts.Neuroimage421597–1608. 10.1016/j.neuroimage.2008.06.025
46
JeonS. Y.HanS. J. (2012). Improvement of the working memory and naming by transcranial direct current stimulation.Ann. Rehabil. Med.36585–595. 10.5535/arm.2012.36.5.585
47
KhatebA.AbutalebiJ.MichelC. M.PegnaA. J.Lee-JahnkeH.AnnoniJ. M. (2007). Language selection in bilinguals: a spatio-temporal analysis of electric brain activity.Int. J. Psychophysiol.65201–213. 10.1016/j.ijpsycho.2007.04.008
48
KhatebA.MichelC. M.PegnaA. J.O’DochartaighS. D.LandisT.AnnoniJ. M. (2003). Processing of semantic categorical and associative relations: an ERP mapping study.Int. J. Psychophysiol.4941–55. 10.1016/S0167-8760(03)00076-X
49
KleinD.MilnerB.ZatorreR. J.MeyerE.EvansA. C. (1995). The neural substrates underlying word generation: a bilingual functional-imaging study.Proc. Natl. Acad. Sci. U.S.A.922899–2903. 10.1073/pnas.92.7.2899
50
KnebelJ. F.MurrayM. M. (2012). Towards a resolution of conflicting models of illusory contour processing in humans.Neuroimage592808–2817. 10.1016/j.neuroimage.2011.09.031
51
KoechlinE.OdyC.KouneiherF. (2003). The architecture of cognitive control in the human prefrontal cortex.Science3021181–1185. 10.1126/science.1088545
52
KolbB.TaylorL. (1981). Affective behavior in patients with localized cortical excisions: role of lesion site and side.Science21489–91. 10.1126/science.7280683
53
KrollJ. F.StewartE. (1994). Category interference in translation and picture naming: evidence for asymmetric connections between bilingual memory representations.J. Mem. Lang.33149–174. 10.1016/j.jecp.2008.10.004
54
KrollJ. F.van HellJ. G.TokowiczN.GreenD. W. (2010). The revised hierarchical model: a critical review and assessment.Biling13373–381. 10.1017/S136672891000009X
55
LeeT. W.GirolamiM.SejnowskiT. J. (1999). Independent component analysis using an extended infomax algorithm for mixed subgaussian and supergaussian sources.Neural Comput.11417–441. 10.1162/089976699300016719
56
MacDonaldA. W.IIICohenJ. D.StengerV. A.CarterC. S. (2000). Dissociating the role of the dorsolateral prefrontal and anterior cingulate cortex in cognitive control.Science2881835–1838. 10.1126/science.288.5472.1835
57
MetukiN.SelaT.LavidorM. (2012). Enhancing cognitive control components of insight problems solving by anodal tDCS of the left dorsolateral prefrontal cortex.Brain Stimul.5110–115. 10.1016/j.brs.2012.03.002
58
MeuterR. F. I.AllportA. (1999). Bilingual language switching in naming: asymmetrical costs of language selection.J. Mem. Lang.4025–40. 10.1037/a0034060
59
MichelC. M.MurrayM. M.LantzG.GonzalezS.SpinelliL.Grave de PeraltaR. (2004). EEG source imaging.Clin. Neurophysiol.1152195–2222. 10.1016/j.clinph.2004.06.001
60
MiyakeA.FriedmanN. P.EmersonM. J.WitzkiA. H.HowerterA.WagerT. D. (2000). The unity and diversity of executive functions and their contributions to complex ‘Frontal Lobe’ tasks: a latent variable analysis.Cogn. Psychol.4149–100. 10.1006/cogp.1999.0734
61
NakanoH.SaronC.SwaabT. Y. (2010). Speech and Span: working memory capacity impacts the use of animacy but not of world knowledge during spoken sentence comprehension.J. Cogn. Neurosci.222886–2898. 10.1162/jocn.2009.21400
62
Neurocare Group (2013). DC-STIMULATOR PLUS, User Manual.Sydney, NSW: Neurocare Group.
63
NitscheM. A.BoggioP. S.FregniF.Pascual-LeoneA. (2009). Treatment of depression with transcranial direct current stimulation (tDCS): a Review.Exp. Neurol.21914–19. 10.1016/j.expneurol.2009.03.038
64
NitscheM. A.CohenL. G.WassermannE. M.PrioriA.LangN.AntalA.et al (2008). Transcranial direct current stimulation: state of the art 2008.Brain Stimul.1206–223. 10.1016/j.brs.2008.06.004
65
NitscheM. A.DoemkesS.KarakoseT.AntalA.LiebetanzD.LangN.et al (2007). Shaping the effects of transcranial direct current stimulation of the human motor cortex.J. Neurophysiol.973109–3117. 10.1152/jn.01312.2006
66
NitscheM. A.LiebetanzD.LangN.AntalA.TergauF.PaulusW. (2003). Safety criteria for transcranial direct current stimulation (tDCS) in humans.Clin. Neurophysiol.1142220–2223. 10.1016/S1388-2457(03)00235-9
67
NoppeneyU.PhillipsJ.PriceC. (2004). The neural areas that control the retrieval and selection of semantics.Neuropsychologia421269–1280. 10.1016/j.neuropsychologia.2003.12.014
68
OldfieldR. C. (1971). The assessment and analysis of handedness: the edinburgh inventory.Neuropsychologia997–113. 10.1016/0028-3932(71)90067-4
69
ParadisM.LibbenG. (1987). The Assessment of Bilingual Aphasia.Hillsdale, NJ: Lawrence Erlbaum Associates.
70
PerrinF.PernierJ.BertrandO.GiardM. H.EchallierJ. F. (1987). Mapping of scalp potentials by surface spline interpolation.Electroencephalogr. Clin. Neurophysiol.6675–81. 10.1016/0013-4694(87)90141-6
71
PhillipsN. A.KleinD.MercierJ.de BoyssonC. (2006). ERP measures of auditory word repetition and translation priming in bilinguals.Brain Res.1125116–131. 10.1016/j.brainres.2006.10.002
72
PlompG.MichelC. M.HerzogM. H. (2010). Electrical source dynamics in three functional localizer paradigms.Neuroimage53257–267. 10.1016/j.neuroimage.2010.06.037
73
PostmaA. (2000). Detection of errors during speech production: a review of speech monitoring models.Cognition7797–132. 10.1016/S0010-0277(00)00090-1
74
PriceC. J. (2000). The anatomy of language: contributions from functional neuroimaging.J. Anat.197335–359. 10.1046/j.1469-7580.2000.19730335.x
75
PriceC. J. (2010). The anatomy of language: a review of 100 fMRI studies published in 2009.Ann. N. Y. Acad. Sci.119162–88. 10.1111/j.1749-6632.2010.05444.x
76
RadmanN.MouthonM.Di PietroM.GaytanidisC.LeemannB.AbutalebiJ.et al (2016). The role of the cognitive control system in recovery from bilingual aphasia: a multiple single-case fmri study.Neural Plast.2016:8797086. 10.1155/2016/8797086
77
RihsT.MichelC. M.ThutG. (2007). Mechanisms of selective inhibition in visual spatial attention are indexed by alpha-band EEG synchronization.Eur. J. Neurosci.25603–610. 10.1111/j.1460-9568.2007.05278.x
78
RihsT.MichelC. M.ThutG. (2009). A bias for posterior α-band power suppression versus enhancement during shifting versus maintenance of spatial attention.Neuroimage44190–199. 10.1016/j.neuroimage.2008.08.022
79
Robins WahlinT.-B.LuszczM.WahlinA.ByrneG. (2015). Non-Verbal and verbal fluency in prodromal huntington’s disease.Dement. Geriatr. Cogn. Disord.5517–529. 10.1159/000441942
80
RomeiV.ThutG.SilvantoJ. (2016). Information-based approaches of noninvasive transcranial brain stimulation.Trends Neurosci.39782–795. 10.1016/j.tins.2016.09.001
81
RoweJ.ToniI.JosephsO.FrackowiakR. S.PassinghamR. (2000). The prefrontal cortex: response selection or maintenance within working memory?Science2881650–1660. 10.1126/science.288.5471.1656
82
RuffR. M.AllenC. C.FarrowC. E.NiemannH.WylieT. (1994). Figural fluency: differential impairment in patients with left versus right frontal lobe lesions.Arch. Clin. Neuropsychol.941–55. 10.1093/arclin/9.1.41
83
ShaoZ.JanseE.VisserK.MeyerA. S. (2014). What do verbal fluency tasks measure? Predictors of verbal fluency performance in older adults.Front. Psychol.5:772. 10.3389/fpsyg.2014.00772
84
ShollA.SankaranarayananA.KrollJ. F. (1995). Transfer between picture naming and translation: a test of asymmetries in bilingual memory.Psychol. Sci.645–49. 10.1111/j.1467-9280.1995.tb00303.x
85
SnodgrassJ. G.VanderwartM. (1980). A standardized set of 260 pictures: norms for name agreement, image agreement, familiarity, and visual complexity.J. Exp. Psychol. Hum. Learn.6174–215. 10.1037/0278-7393.6.2.174
86
SpinelliL.Gonzalez AndinoS.LantzG.SeeckM.MichelC. M. (2000). Electromagnetic inverse solutions in anatomically constrained spherical head models.Brain Topogr.13115–125. 10.1023/A:1026607118642
87
SteinM.FederspielA.KoenigT.WirthM.LehmannC.WiestR.et al (2009). Reduced frontal activation with increasing 2nd language proficiency.Neuropsychologia472712–2720. 10.1016/j.neuropsychologia.2009.05.023
88
SundermanG.KrollJ. F. (2006). First language activation during second language lexical processing: an investigation of lexical form, meaning, and grammatical class.Stud. Second Lang. Acquis.28387–422. 10.1017/S0272263106060177
89
ThielscherA.AntunesA.SaturninoG. (2015). “Field modeling for transcranial magnetic stimulation: a useful tool to understand the physiological effects of TMS?” inProceedings of the 37th Annual International Conference of the EEE Enginering in Medicine and Biology Society (Milan) 222–225. 10.1109/EMBC.2015.7318340
90
TombaughT. N.KozakJ.ReesL. (1999). Normative data stratified by age and education for two measures of verbal fluency: FAS and animal naming.Arch. Clin. Neuropsychol.14167–177.
91
WhitneyC.KirkM.O’SullivanJ.Lambon RalphM. A.JefferiesE. (2012). Executive semantic processing is underpinned by a large-scale neural network: revealing the contribution of left prefrontal, posterior temporal, and parietal cortex to controlled retrieval and selection using TMS.J. Cogn. Neurosci.24133–147. 10.1162/jocn_a_00123
92
ZhangS.ThompsonN. (2004). DIALANG:A diagnostic language assessment system (review).Can. Mod. Lang. Rev.61290–293. 10.1353/cml.2005.0011
93
ZyssetS.HuberO.SamsonA.FerstlE. C.von CramonD. Y. (2003). Functional specialization within the anterior medial prefrontal cortex: a functional magnetic resonance imaging study with human subjects.Neurosci. Lett.335183–186. 10.1016/S0304-3940(02)01196-5
Summary
Keywords
tDCS, EEG/ERP, dorsolateral prefrontal cortex, bilingualism, language production
Citation
Radman N, Britz J, Buetler K, Weekes BS, Spierer L and Annoni J-M (2018) Dorsolateral Prefrontal Transcranial Direct Current Stimulation Modulates Language Processing but Does Not Facilitate Overt Second Language Word Production. Front. Neurosci. 12:490. doi: 10.3389/fnins.2018.00490
Received
15 December 2017
Accepted
29 June 2018
Published
25 July 2018
Volume
12 - 2018
Edited by
Roberto Cecere, University of Glasgow, United Kingdom
Reviewed by
Juha Silvanto, University of Westminster, United Kingdom; James Dowsett, Ludwig-Maximilians-Universität München, Germany
Updates
Copyright
© 2018 Radman, Britz, Buetler, Weekes, Spierer and Annoni.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Narges Radman, radmann@ipm.ir; narges.radman@gmail.com
This article was submitted to Perception Science, a section of the journal Frontiers in Neuroscience
Disclaimer
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.