Modulating the resting-state functional connectivity patterns of language processing areas in the human brain with anodal transcranial direct current stimulation applied over the Broca’s area

Jianwei Cao; Hanli Liu; George Alexandrakis

doi:10.1117/1.NPh.5.2.025002

27 February 2018 Modulating the resting-state functional connectivity patterns of language processing areas in the human brain with anodal transcranial direct current stimulation applied over the Broca’s area

Jianwei Cao, Hanli Liu, George Alexandrakis

Author Affiliations +

Neurophotonics, Vol. 5, Issue 2, 025002 (February 2018). https://doi.org/10.1117/1.NPh.5.2.025002

Abstract

Cortical circuit reorganization induced by anodal transcranial direct current stimulation (tDCS) over the Broca’s area of the dominant language hemisphere in 13 healthy adults was quantified by functional near-infrared spectroscopy (fNIRS). Transient cortical reorganization patterns in steady-state functional connectivity (seed-based and graph theory analysis) and temporal functional connectivity (sliding window correlation analysis) were recorded before, during, and after applying high current tDCS (1 mA, 8 min). fNIRS connectivity mapping showed that tDCS induced significantly (p < 0.05) increased functional connectivity between Broca’s area and its neighboring cortical regions while it simultaneously decreased the connectivity to remote cortical regions. Furthermore, the anodal stimulation caused significant increases to the functional connectivity variability (FCV) of remote cortical regions related to language processing. In addition to the high current tDCS, low current tDCS (0.5 mA, 2 min 40 s) was also applied to test whether the transient effects of lower stimulation current could qualitatively predict cortical connectivity alterations induced by the higher currents. Interestingly, low current tDCS could qualitatively predict the increase in clustering coefficient and FCV but not the enhancement of local connectivity. Our findings indicate the possibility of combining future studies fNIRS with tDCS at lower currents to help guide therapeutic interventions.

1. Introduction

Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique that has been applied to modulate cognitive function so as to enhance performance in healthy subjects¹^–³ and facilitate neurorehabilitation during stroke recovery.⁴ Typically, tDCS is used with the intention to alter cortical excitability by delivering weak currents (1 to 2 mA) through a pair of anode–cathode electrodes for 8 to 20 min.⁵ The primary mechanisms of the excitability shifts are depolarization of resting membrane potential by anodal stimulation and hyperpolarization of resting membrane potential by cathodal stimulation, as has been illustrated in animal studies.⁶^,⁷ In the context of human subject studies on areas of the brain controlling language processing, which is the focus of this work, it has been suggested that anodal tDCS over either Broca’s area or Wernicke’s area can improve naming accuracy or speed both in stroke-induced aphasia patients⁸^–¹⁰ and in healthy subjects.³^,¹¹^,¹²

Despite the above-described growing number of studies investigating how tDCS affects language performance through stimulation of language processing areas, there has been relatively little study of how tDCS affects cortical functional connectivity reorganization during those interventions. Meinzer et al.¹³ used resting-state functional magnetic resonance imaging (fMRI) to assess the impact of anodal tDCS on the functional connectivity networks of healthy subjects. Their results showed an increase of functional connectivity strength in language-associated regions, such as dorsolateral and medial prefrontal regions, presupplementary motor area (pre-SMA) and SMA but a decrease in more posterior and occipital regions. Marangolo et al.¹⁴ revealed that increased connectivity strength was most pronounced in the left hemispheric structures related to planning, maintenance, and execution of speech. They used resting-state fMRI to study the effects of bilateral tDCS on aphasia patients, with the anodal patch applied over the Broca’s area of the lesion-containing left brain hemisphere and the cathodal patch applied over the homologous contralateral brain region. Also, the presently few studies of tDCS-mediated changes on connectivity have focused on steady-state rather than the temporal aspects of those changes. The fact that tDCS could affect dynamic fluctuations in functional connectivity has not yet been explored even though the study of these fluctuations has been shown to elucidate fundamental properties of brain networks.¹⁵^–¹⁸

One additional aspect of tDCS studies that remains underexplored to date is the effect of electrode placement on language performance. For example, Monti et al.¹⁹ found out that anodal tDCS over the left Broca’s area did not improve naming accuracy in chronic nonfluent aphasic patients, whereas Baker et al.⁸ reported improved naming performance for stimulation over the same area. It is conceivable that due to heterogeneity of lesion effects in aphasic patients, an optimal tDCS electrode placement for language enhancement and recovery would differ between subjects. In an ideal scenario, investigators would be able to try different tDCS electrode placements rapidly to find the one that could maximize performance improvements in each patient. However, the current intensities used in tDCS interventions produce long-lasting effects on cortical activity with durations from hours to days,²⁰ which makes the rapid testing of multiple electrode placements impossible. One potential way to reduce testing time could be to use lower tDCS currents, which reduce the duration of tDCS effects on the brain hemodynamics to a few minutes. We have previously shown that use of lower tDCS currents allows multiple tDCS placements to be tested within one single session.²¹ However, the potential utility of lower current tDCS depends on whether it could produce qualitatively similar patterns of connectivity change to those seen during the high current tDCS condition.

This work is a first step toward addressing the aformentioned knowledge gaps in tDCS studies, namely (1) to study the effect of stimulation current on the steady-state and temporal connectivity changes in brain networks related to language processing and (2) to study whether lower current tDCS patterns could produce changes in connectivity patterns that are qualitatively similar to those observed during and after the application of larger tDCS currents. Here, we have studied the effect of tDCS current intensity on the cortical connectivity maps of healthy adults by use of functional near-infrared spectroscopy (fNIRS). fNIRS can detect changes in the concentration of oxyhemoglobin (HbO) and deoxyhemoglobin (Hb) resulting from neurovascular coupling²² and has been increasingly shown to be a viable alternative neuroimaging modality to fMRI despite its lower spatial resolution and ability to only map cortical activation.²³^,²⁴ Also, the compatibility of fNIRS with tDCS makes it advatangeous for studying brain reorganization induced during tDCS.²⁵ Furthermore, fNIRS provides higher temporal resolution relative to fMRI, which is valuable for investigating time-varying functional connectivity in the brain.²⁶

Several fNIRS-tDCS studies have been reported in the past years. Merzagora et al.²⁷ stated that a significant HbO increase was observed by fNIRS over a prefrontal cortex (Fp1) stimulation site for 1 mA of anodal tDCS applied for 8 min. Takai et al.²⁸ found that both anodal and cathodal stimulation of 1 mA for 20 min applied to the primary motor cortex (M1) induced a significant HbO decrease in the contralateral premotor cortex, SMA, and M1. Our group²⁵ revealed that bihemispheric tDCS with anode on the left M1 and the cathode on the right M1 changed resting-state connectivity from intrahemispheric to interhemispheric and increased flexion speed when subjects performed a wrist flexion task. However, to our knowledge, no fNIRS studies to date have reported how tDCS affects the reorganization of language area networks. Specifically, in this work, we wanted to assess the alterations of functional connectivity due to anodal tDCS over the left Broca’s area and whether a low tDCS current (0.5 mA) could qualitatively predict cortical connectivity patterns occurring after a high tDCS current (1 mA), which is a standard intervention current choice in the literature. Cortical connectivity patterns were computed using seed-based functional connectivity with the seed located on the anodally stimulated left Broca’s area and on its homologous area on the right hemisphere. Furthermore, we computed changes in time-variant functional connectivity with the same seeds and also employed graph theory analysis to assess global pattern changes with tDCS current intensity. Cortical reorganization patterns of functional networks created in response to low current tDCS were compared to those during and after high current tDCS. The findings of this work are discussed in light of their potential future utility for helping guide therapeutic interventions.

2. Methods and Materials

2.1.

Subjects

Thirteen right-hand subjects (2 females, 11 males, $mean \pm SD age = 35.4 \pm 8.4$ ) were recruited in this study. Subject handedness was determined by the Edinburgh Handedness Inventory.²⁷ All subjects were healthy and did not have a history of neurological disorders. Written informed consent was obtained from each participant before the experiments. The studies were carried out under the approval of the University of Texas at Arlington Institutional Review Board protocol (UTA #2015-0819).

2.2.

fNIRS Imaging Setup Combined with tDCS

Figure 1(a) illustrates the overall instrumentation setup. A continuous-wave fNIRS imaging system (LABNIRS, Shimadzu, Japan) was used, which utilized near-infrared light diode sources (at wavelengths of 780, 805, and 830 nm) and photomultiplier detectors. A schematic of the fNIRS source–detector geometry is shown in Fig. 1(b). The setup geometry consisted of 26 sources and 28 detectors with a separation of 3 cm, resulting in 83 source–detector channels, which were inserted into the optode holder on the subject’s head. This probe geometry covered cortical areas known to be part of the language network, including the Broca and Wernicke’s areas of the left hemisphere and their homologous locations in the right hemisphere, as well as some prefrontal cortical regions, including the frontopolar, dorsolateral prefrontal cortex, and premotor areas. fNIRS signals were sampled at a frequency of 12.35 Hz.

Fig. 1

(a) A schematic representation of the overall instrumentation setup. The protocol only with “rest” is shown on the screen. FNIRS optodes and tDCS electrodes were placed on the subject’s head, as described in Sec. 2.2. (b) Schematic of the fNIRS probe geometry with 26 sources and 28 detectors placed over a subject’s head. The separation of all source and detectors is 3 cm (red dots: sources, blue dots: detectors). (c) Placement of the fNIRS-tDCS assembly on a subject’s head. The gray arrow points to the wire connecting the cathode patch and the red arrow points to the wire connecting the anode patch.

A coregistration procedure was performed based on cranial landmark measurements on all the subjects²³ to estimate the cortical regions covered by the fNIRS source–detector geometry. A motion tracking system (Fastrak, POLHEMUS) was used to measure five reference cranial landmarks (nasion, inion, left and right preauricular points and vertex), as well as the location of all source and detector optodes. The Montreal Neurological Institute coordinates were calculated from the real-world stereotaxic coordinates of the optodes with the five cranial landmarks as reference positions.²⁹ NIRS Statistical Parametric Mapping³⁰ was used to register the optodes on the standard MRI brain template and identify the Brodmann areas of each source–detector channel. Figure 2 shows the spatially registered channels (averaged over all 13 subjects) on the standard human brain atlas.

Fig. 2

Coregistration of fNIRS source-detector channels (mid-way points between source and detector pairs) on a standard human brain atlas: (a) Sagittal view (left), (b) sagittal view (right), (c) top view, and (d) coronal view. tDCS anodal patch (red dashed square) and cathodal patch (black dashed square) are placed on the left FC5 position (a) and right Fp2 position (d).

tDCS current was applied by a battery-driven electrical stimulator (Phoresor II, IOMED Inc., Salt Lake City, Utah) connected by a pair of saline-soaked gauze covered gel electrodes ( $5 \times 5 cm$ ; IOMED Inc., Salt Lake City). The tDCS patches were placed on the head using the EEG International $10 / 20$ system³¹ as a reference. The anodal patch was placed over the FC5 position to stimulate left Broca’s area, which is centered at channel 27 in our setup, and a cathodal patch was centered over the Fp2 position¹² as control [dashed square in Figs. 1(b), 2(a), and 2(d)]. Two 0.5-cm-diameter holes were punctured on each patch to fit through optical fiber bundles that overlapped spatially with the patches.

2.3.

Protocol Design

Subjects were asked to rest for 32 min with eyes closed but awake while fNIRS measurements were performed. The first 6 min before tDCS were used as baseline, and then 0.5-mA current tDCS, henceforth referred to as low current tDCS, was applied for 2 min 40 s. Subsequently, subjects were asked to continue to rest for 5 min and 20 s, which was ample time for any hemodynamic signatures of the low current tDCS to disappear⁵^,²¹^,³² while also maintaining a practically feasible total protocol duration. After that, a higher current of 1 mA, henceforth referred to as high current tDCS, was applied for 8 min. These current intensity and stimulation duration settings were selected to be near the lower bound of values used in current clinical interventions involving tDCS in order to minimize protocol duration and improve subject comfort. The fNIRS imaging session concluded with 10 min of acquiring data while subjects rested so as to record changes in the hemodynamics immediately after the end of the 1-mA tDCS. A schematic of the timeline of the protocol is shown in Fig. 3. All subsequent fNIRS data analyses were performed separately for each one of the four stimulation conditions: no tDCS (0 to 6 min), low current tDCS (6 to 14 min), high current tDCS (14 to 22 min), and after high current tDCS (22 to 32 min). The low current tDCS session composed of the sum of 2 min 40 s of 0.5 mA followed by 5 min and 20 s of no tDCS. It was found necessary to merge those two periods into one, combining signals measured during low current tDCS and its subsequent decay to baseline hemodynamics, in order to attain sufficient fNIRS time-series data to enable all subsequent functional connectivity analyses for this stimulation condition. Subjects were not told when tDCS was applied.

Fig. 3

The tDCS protocol timeline.

2.4.

Data Preprocessing

This study used a publically available toolbox named Homer³³ to process the time-series fNIRS data. First, the data from each channel were band-pass filtered (0.009 to 0.1 Hz) to reduce low-frequency drift and high-frequency neurophysiological noise.³⁴^,³⁵ Subsequently, changes in HbO relative to the baseline were quantified by a modified Beer–Lambert’s law.³⁵ Correlation-based signal improvement³⁶ was adopted to remove motion artifacts based on negative correlation between oxygenated and deoxygenated hemoglobin dynamics.

2.5.

Data Analysis

2.5.1.

Seed-based functional connectivity

Channels 27 and 34 were chosen as seeds. Channel 27 was located over the stimulated left Broca’s area, which is the dominant language area for right-handed subjects. Channel 34 was located over the right Broca’s homologue, which is known to endow verbal communication with additional meaning by contributing to the analysis of emotional and tonal context in language.³⁷ Pearson’s correlation coefficients were calculated between the seed channel and all the other channels for the entire time interval allocated to each stimulation condition (Fig. 3). Subsequently, a Fisher’s $z$ -transformation was employed to transform the correlation coefficients to $z$ -values³⁸ for each subject. For each stimulation, session $z$ -values between all the other channels and the seed channel were averaged across subjects. Significant seed-based functional connectivity changes between stimulation sessions were identified based on the following criteria:³⁹ (1) using a paired $t$ -test, for each $z$ -value between sessions (i.e., no tDCS versus low current tDCS) at a threshold of uncorrected $p < 0.05$ and (2) $z$ -values were significantly different from zero at $p < 0.05$ significance using a one-sample $t$ -test in at least one session. The significant changes in functional connectivity were mapped by BrainNet Viewer software.⁴⁰

2.5.2.

Time-variant functional connectivity

Time-variant functional connectivity was also calculated for the same seed locations. For each individual subject, sliding-window correlation (SWC) analysis¹⁵^,⁴¹ was performed on each one of the four session data. In our SWC analysis, a 60-s time window was selected and then shifted by an increment of 1 s along the time course of each session.¹⁵ Then, the functional connectivity within each sliding-time window was calculated between the seed channel and all other channels using the Pearson correlation method. For each stimulation session, functional connectivity variability (FCV) was calculated as the standard deviation of the correlation coefficient along time.⁴² For group analyses, the FCV of each correlation coefficient for each channel was averaged across subjects for each session. Paired $t$ -tests were used to compare the FCV of each correlation coefficient between sessions, i.e., no tDCS versus low current tDCS, to determine whether the FCV between the seed channel and other channels in each stimulation session was significantly different ( $p < 0.05$ ). The channels with significant changes in FCV of the correlation coefficients with respect to the seed channels were plotted using the BrainNet Viewer software.⁴⁰

2.5.3.

Graph theory analysis

Graph theory analysis was applied to investigate the changes in topographical patterns of functional networks across the entire cortical regions mapped by fNIRS. In our study, for each subject, we defined each fNIRS channel as a node, resulting in a total of 83 nodes. The edges $e_{i j}$ in the network were determined by setting a threshold, T, to the $83 \times 83$ Pearson’s correlation matrix values $r_{i j}$ by the following formula:⁴³

e_{i j} = {\begin{matrix} 1, & if | r_{i j} | > T \\ 0, & otherwise \end{matrix} .

A certain range of [ $T_{\min}$ $T_{\max}$ ] for $T$ was chosen so that $T_{\min}$ excluded weak and potentially nonsignificant connections while selecting the Pearson’s correlation coefficients that were significant and corresponded to $p < 0.05$ in the channel-level correlation matrix.⁴⁴ $T_{\max}$ was set based on the condition that the mean node degree $\bar{k} \geq 2 \log$ (# of channels), which was equal to 8.8 for our setup.⁴⁴^,⁴⁵ This value for $T_{\max}$ can be interpreted as the node being tested having connections with no less than 8.8 other nodes on average. Furthermore, this threshold value meant that the total number of edges in the network was no less than 365 ( $= 83 \times 8.8 / 2$ ), equivalent to around 10% of the maximum number of edges possible ( $83 \times 82 / 2 = 3403$ ) in a network of 83 nodes.⁴⁶^,⁴⁷ Once the statistically significant nodes and edges were identified, graph theory-based metrics were calculated by the Gretna software for each session.⁴³

Clustering coefficient $C_{p}$ is the average of the clustering coefficients of all nodes. $C_{p} (i)$ of certain node $i$ is defined as the following:⁴⁸

C_{p} (i) = \frac{2 N_{i}}{k_{node} (i) [k_{node} (i) - 1]},

where

N_{i}

denotes the number of existing connections among the neighbors of node

i

and

k_{node}

represents the number of edges that are connected to node

i

.

Characteristic path length $L_{p}$ is defined as the average of the shortest path lengths between all pairs of nodes:⁴⁸

L_{p} = \frac{1}{N (N - 1)} \sum_{i \neq j \in G} d_{i j},

where

N

is the total node and

d_{i j}

is the shortest path length between node

i

and node

j

.

The global efficiency $E_{glob}$ and local efficiency $E_{loc}$ are defined as the following:⁴⁹

E_{glob} = \frac{1}{N (N - 1)} \sum_{i \neq j \in G} \frac{1}{d_{i j}}, E_{loc} = \frac{1}{N (N - 1)} \sum_{i \in G} E (G_{i}),

where

N

is the total node,

d_{i j}

is the shortest path length between node

i

and node

j

and

E (G_{i})

is the global efficiency of the subgraph composed of the nearest neighbors of node

i

.

Small-world properties (clustering coefficient $C_{p}$ and characteristic path length $L_{p}$ ) and efficiency parameters (global efficiency $E_{glob}$ and local efficiency $E_{loc}$ ) were calculated by setting the threshold from $T_{\min}$ to $T_{\max}$ with a step size of 0.01.⁴⁴^,⁴⁵^,⁵⁰ Then, the area under the curve (AUC) for each metric and stimulation session was computed. For group analyses, paired $t$ -tests were used to compare the AUC for each metric and any significant changes ( $p < 0.05$ ) between stimulation sessions were identified.

3. Results

3.1.

Functional Connectivity Analysis Using the Stimulated Cortical Region and Contralateral Cortical Region as the Seed

Paired $t$ -tests for low current tDCS versus high current tDCS, and low current tDCS versus after high current tDCS were conducted when the seed was on channel 27 (Broca’s area) and channel 34 (Broca’s homologue area), and no significant changes in functional connectivity strength were found for either seed between stimulation conditions. These comparisons support our hypothesis that low current tDCS creates similar patterns to those of high current tDCS. In order to understand better how current intensity affected changes in connectivity patterns we focused all further analyses on comparisons with respect to the no tDCS condition, as described in the subsections here below.

3.1.1.

Using the stimulated left Broca’s area as the seed

Figure 4 shows the detector locations over cortical areas with significant changes in connectivity strength, as deduced from paired $t$ -tests between the no tDCS condition versus each one of the other stimulation conditions (low current tDCS, high current tDCS, and after high current tDCS). Figures 4(a)–4(c) show lateral left views and Figs. 4(d)–4(f) show top views for each stimulation condition comparison. Figures 4(a) and 4(d) indicate that functional connectivity strength decreased significantly for longer distance connections with respect to the seed location at the left Broca’s area during the low current tDCS condition. Furthermore, Figs. 4(b) and 4(e) show that high current tDCS not only decreased significantly the connectivity strength with longer distance brain regions but also significantly increased the connectivity strength with nearby regions. This near-neighbor effect persisted during the after high current tDCS condition, as shown in Figs. 4(c) and 4(f), but the change in connectivity strength compared to the no tDCS condition was smaller compared to the high current tDCS condition. The increased connectivity strength between regions indicated by green and yellow ovals in Figs. 4(b) and 4(c) was most pronounced in the left dorsolateral prefrontal [green oval in Fig. 4(b), $p = 0.0010$ ], premotor, and SMA areas [yellow oval in Figs. 4(b) and 4(c), $p = 0.0031$ and $p = 0.0241$ ], which regulate the planning, maintenance, and execution of speech.⁵¹^–⁵⁴

Fig. 4

Changes in functional connectivity strength between the left Broca’s area (seed) and other cortical regions induced by different tDCS stimulation conditions. Group-level significant differences ( $p < 0.05$ ) in connectivity strength between pairs of detector locations for the no tDCS condition versus low current tDCS [(a) and (d)], versus high current tDCS [(b) and (e)], and versus after high current tDCS [(c) and (f)] stimulation conditions. Green ovals: detector pair locations showing significantly increased connectivity with the dorsolateral prefrontal area. Yellow ovals: detector pair locations showing significantly increased connectivity with the premotor and SMA areas.

3.1.2.

Using the right Broca’s homologue area as the seed

Figure 5 shows detector locations over cortical areas with significant changes in connectivity strength using the right Broca’s homologue as the seed. These connectivity changes were deduced from paired $t$ -tests between the no tDCS condition versus each one of the other stimulation conditions. Several qualitative similarities were observed in the connectivity strength changes when comparing to the left Broca’s area stimulation conditions: (1) low current tDCS induced significant decreases in connectivity strength for long-distance brain regions [Figs. 5(a) and 5(d)]; (2) high current tDCS brought in significant increases in connectivity strength for the nearby right dorsolateral prefrontal [green oval in Fig. 5(b), $p = 0.0045$ ] and premotor and SMA areas [yellow ovals in Figs. 5(b), $p = 0.0329$ , and $p = 0.0241$ , respectively, and Fig. 5(c), $p = 0.0417$ ]; (3) similar distance-dependent effects were seen for the after high current tDCS condition [Figs. 5(c) and 5(f)]. However, in contrast to the case of the left Broca’s area was the seed, high current tDCS also brought about a significant increase in connectivity strength with the right Broca’s homologue for detector location over another language-related cortical area, the superior temporal gyrus⁵⁵ [purple oval in Fig. 5(b), $p = 0.0332$ ]. This cortical region is part of Wernicke’s area and is related with prosody comprehension.³⁷ In addition, we attribute the larger number of channels observed with significant connectivity changes compared to the case of stimulating left Broca’s area to the brain’s known asymmetry with more connections stemming from the right hemisphere of right-handed adults.⁵⁶

Fig. 5

Changes in functional connectivity strength between the right Broca’s homologue (seed) and other cortical regions induced by different tDCS stimulation conditions. Group-level significant differences ( $p < 0.05$ ) in connectivity strength between pairs of detector locations for the no tDCS condition versus low current tDCS [(a) and (d)], versus high current tDCS [(b) and (e)], and versus after high current tDCS [(c) and (f)] stimulation conditions. Green ovals: detector pair locations showing significantly increased connectivity with the dorsolateral prefrontal area. Yellow ovals: detector pair locations showing significantly increased connectivity with premotor and SMA areas. Purple oval: detector pair locations showing significantly increased connectivity with the superior temporal gyrus.

In addition, Tables 1 and 2 list the Brodmann areas of channels with significantly decreased functional connectivity indicated in Figs. 4 and 5 for the three stimulation sessions compared to the no tDCS condition when the seed was on channels 27 and 34, respectively.

Table 1

Brodmann areas with significantly decreased steady-state functional connectivity when the seed was on channel 27 (left Broca’s area).

Seed: channel 27	Low current tDCS	High current tDCS	After high current tDCS
Ipsilateral hemisphere
Brodmann area 5	FC decreased ( $p = 0.0064$ )		FC decreased ( $p = 0.0401$ )
Brodmann area 8	FC decreased ( $p = 0.0416$ )
Contralateral hemisphere
Brodmann area 5		FC decreased ( $p = 0.0420$ )
Brodmann area 7		FC decreased ( $p = 0.0187$ )
Brodmann area 21	FC decreased ( $p = 0.0311$ )

Table 2

Brodmann areas with significantly decreased steady-state functional connectivity when the seed was on channel 34 (right Broca’s homologue).

Seed: channel 34	Low current tDCS	After high current tDCS
Ipsilateral hemisphere
Brodmann area 7	FC decreased ( $p = 0.0250$ , 0.0201)	FC decreased ( $p = 0.0442$ )
Brodmann area 8	FC decreased ( $p = 0.0174$ )	FC decreased ( $p = 0.0094$ )
Brodmann area 19	FC decreased ( $p = 0.0068$ )
Contralateral hemisphere
Brodmann area 3	FC decreased ( $p = 0.0194$ )	FC decreased ( $p = 0.0277$ )
Brodmann area 5		FC decreased ( $p = 0.0496$ , 0.0241)
Brodmann area 7		FC decreased ( $p = 0.0331$ )
Brodmann area 8	FC decreased ( $p = 0.0479$ )	FC decreased ( $p = 0.0094$ , 0.0087)

3.2.

Time-Variant Functional Connectivity using the Stimulated Left Broca’s Area as Seed

Figures 6(a)–6(c) illustrate the detector locations over cortical regions (blue dots) with significantly increased FCV of the time-variant functional connectivity when the left Broca’s area was the stimulated cortical region (red dots). Each row in that figure corresponds to a distinct stimulation condition, namely low current tDCS, high current tDCS, and after high current tDCS. The cortical regions with increased FCV were similar for all three tDCS conditions: left Wernicke’s area, right Wernicke’s area, and right frontopolar area. Similar results were also found when using the right Broca’s homologue as the seed. The left and right Wernicke’s areas and the left frontopolar area also showed significantly increased FCV for all three tDCS sessions. Due to their similarity to the findings for the seed over the left Broca’s area, we do not show these results for brevity. These cortical regions belong to higher order brain regions (including language processing areas and frontopolar prefrontal cortex)⁵⁷^,⁵⁸ and are remote to the seed regions. These FCV findings are in contrast to the above steady-state connectivity strength findings that mostly highlighted local connectivity strength increases.

Fig. 6

Brain cortical regions of significant increased FCV (blue dots, FCV, $p < 0.05$ ) with stimulated cortical region (red dots) induced by different tDCS stimulation conditions: (a), (d), (g) no tDCS versus low current tDCS; (b), (e), (h) no tDCS versus high current tDCS; and (c), (f), (i) no tDCS versus after high current tDCS.

Tables 3 and 4 list Brodmann areas of channels with significantly increased FCV for the three stimulation sessions compared to the no tDCS condition when the seed was on channels 27 and 34, respectively.

Table 3

Brodmann areas with significantly increased FCV when the seed was on channel 27 (left Broca’s area).

Seed: channel 27	Low current tDCS	High current tDCS	After high current tDCS
Ipsilateral hemisphere
Brodmann area 22			FCV increased ( $p = 0.0189$ )
Brodmann area 40	FCV increased ( $p = 0.0468$ )	FCV increased ( $p = 0.0358$ )	FCV increased ( $p = 0.0421$ )
Contralateral hemisphere
Brodmann area 10	FCV increased ( $p = 0.0002$ )	FCV increased ( $p = 0.0085$ )	FCV increased ( $p = 0.0001$ )
Brodmann area 22	FCV increased ( $p = 0.0321$ )	FCV increased ( $p = 0.0303$ )	FCV increased ( $p = 0.0232$ )
Brodmann area 39		FCV increased ( $p = 0.0383$ )	FCV increased ( $p = 0.0033$ )
Brodmann area 40	FCV increased ( $p = 0.0356$ , 0.0044)	FCV increased ( $p = 0.0021$ )	FCV increased ( $p = 0.0335$ , 0.0111)

Table 4

Brodmann areas with significantly increased FCV when the seed was on channel 34 (right Broca’s homologue).

Seed: channel 34	Low current tDCS	High current tDCS	After high current tDCS
Ipsilateral hemisphere
Brodmann area 39	FCV increased ( $p = 0.0422$ )	FCV increased ( $p = 0.0420$ )	FCV increased ( $p = 0.0013$ )
Brodmann area 40			FCV increased ( $p = 0.0278$ )
Contralateral hemisphere
Brodmann area 10	FCV increased ( $p = 0.0405$ )	FCV increased ( $p = 0.0151$ )	FCV increased ( $p = 0.0073$ )
Brodmann area 39	FCV increased ( $p = 0.0230$ )	FCV increased ( $p = 0.0195$ )
Brodmann area 40			FCV increased ( $p = 0.0426$ )

In addition, paired $t$ -tests were conducted for low current tDCS versus high current tDCS, and low current tDCS versus after high current tDCS. There were no significant changes on FCV when the seed was on channel 27 or channel 34. These results are consistent and our findings from comparisons with the no tDCS condition that low current tDCS creates similar FCV patterns to those of high current tDCS.

3.3.

Graph Theory Analysis

As described in Sec. 2.5.3, a certain range of functional connectivity thresholds [0.3, 0.59] with the step of 0.01 was chosen to perform graph theory analysis.

Figures 7(a) and 7(b) illustrate the AUC for the cortical network clustering coefficient ( $C_{p}$ ) and local efficiency ( $E_{loc}$ ), respectively, with the AUC being integrated as a function of threshold $T$ ⁴³ for the different stimulation conditions. It was found that the AUC for both $C_{p}$ and $E_{loc}$ had similar trends, showing an increase in going from the no tDCS to the after high current tDCS condition. However, significant changes ( $p < 0.05$ ) were found between the no tDCS condition and all other stimulation conditions for $C_{p}$ , whereas for $E_{loc}$ , the only significant change was between the no tDCS and the after high current tDCS conditions. The significant changes in the AUC of $C_{p}$ for all stimulation conditions relative to baseline reflect an increase in the number of neighboring channel connections, even for low stimulation currents and after the high current tDCS was turned off. In addition, the significant change in the AUC for $E_{loc}$ suggests that the after high current tDCS condition enhanced the formation of network clusters, which is consistent with more efficient communication between nodes within the same immediate neighborhood.⁵⁹^,⁶⁰

Fig. 7

The differences of mean areas under the curve (AUC) for (a) clustering coefficient ( $C_{p}$ ) and (b) local efficiency ( $E_{loc}$ ) computed for the no tDCS, low current tDCS, high current tDCS, and after high current tDCS conditions. The asterisk indicates significant differences ( $p < 0.05$ ). Error bars correspond to one standard deviation to the mean across subjects.

Furthermore, no significant differences were found in the AUC of $C_{p}$ for low current tDCS versus high current tDCS ( $p = 0.604$ ), and low current tDCS versus after high current tDCS ( $p = 0.177$ ). Also, no significant differences were found in the AUC of $E_{loc}$ for low current tDCS versus high current tDCS ( $p = 0.601$ ), and low current tDCS versus after high current tDCS ( $p = 0.120$ ). These results are consistent and our findings from comparisons with the no tDCS condition that low current tDCS creates similar graph theory derived metrics to those of high current tDCS.

4. Discussion

This study systematically explored by use of resting state fNIRS the impact of two different anodal tDCS exposures applied over the left Broca’s area, on functional connectivity reorganization. A lower tDCS current (low current tDCS) was applied to create transient hemodynamic changes, whereas a higher current (high current tDCS) mimicked a current dose that is typical of tDCS interventions reported in the literatures.⁶¹^,⁶² The purpose of this work was to test how tDCS modulated connectivity in the immediate neighborhood of the anodal stimulation location over left Broca’s area versus the connectivity to distant cortical regions involved in language processing. A second goal was to explore whether lower currents could reproduce qualitatively the connectivity patterns seen for typical intervention currents, which would enable the lower currents to serve as predictors of connectivity responses to the higher currents.

4.1.

tDCS Increased Steady-State Functional Connectivity in the Immediate Neighborhood of the Stimulation Location

Seed-based analysis revealed distinct differences in functional connectivity patterns between low current tDCS, high current tDCS, and within 10 min after high current tDCS. Specifically, low current tDCS was more likely to suppress the connections with long distance brain regions, while high current tDCS encouraged increased local connectivity. Interestingly, during the after high current tDCS session, a return to the suppression of long-distance connections was seen, combined with the short-distance increased connectivity seen during the immediately preceding high current tDCS session. It is hypothesized that anodal tDCS-induced subthreshold neuronal depolarization, which reduced the amount of excitatory input required, resulting in the excitation of the stimulated cortical areas near the anodal electrode.⁶^,⁷ The local increase of spontaneous activity might have decreased the neuronal signal to noise ratio, and consequently decreased the synchronization with other remote brain regions,³⁹^,⁶³ which could help explain the observed reduction in connectivity with those regions. Similar results with respect to changes in local versus long-distance connectivity were reported in prior work⁶³ for anodal tDCS, corresponding to high current tDCS in this work, for stimulation over the left primary motor area. In this work, the brain regions with the most pronounced increases in functional connectivity strength were the dorsolateral prefrontal, premotor, and SMA areas, known to be a part of the language processing and production network.⁵¹^–⁵⁴ These results were consistent with prior studies reporting on dual tDCS over the Broca’s areas,¹³^,¹⁴ where it was illustrated that tDCS resulted in significantly increased functional connectivity with the premotor and SMA regions. High current tDCS also brought about a significantly increased functional connectivity between the right Broca’s homologue and the right superior temporal gyrus, which are known to be responsible for prosody comprehension function.³⁷ For the two channels with persistently increased connectivity strength during the after high current tDCS sessions [Figs. 4(c) and 5(c)], the increase in connectivity strength during the high current tDCS session was larger than other increased connections [Figs. 4(b) and 5(b)], with the latter going back to prestimulation levels during the after high current tDCS session [Figs. 4(c) and 5(c)]. Our results suggest that anodal stimulation over the left Broca’s area prepares the task-related language areas by enhancing functional connectivity with these cortical regions. We hypothesize that enhanced baseline connectivity strength between language areas contributed to the beneficial effects of anodal tDCS stimulation previously seen for picture naming performance¹² and verbal fluency⁶⁴ in healthy subjects. Similarly, these findings might be helpful to explain the beneficial effects of anodal stimulation over the affected hemisphere in stroke aphasia patients.⁸^,⁶⁵^,⁶⁶ In addition, low current tDCS could not predict the increased functional connectivity in the immediate neighborhood of the stimulation location induced by high current tDCS. However, for the decreased functional connectivity in the remote regions of the stimulation location, seen in Tables 1 and 2 as a whole, low current tDCS did not produce significant changes for the same Brodmann areas as high current tDCS, though it showed changes for the same 5/7 Brodmann areas as the after high current tDCS condition. These results, taken together with Figs. 4 and 5, indicate that low current tDCS is at best a weak qualitative predictor of the connectivity change patterns seen after high tDCS and no further statistical comparisons were performed to quantify its predictive ability.

4.2.

tDCS Increased Functional Connectivity Variability with High-Order Cognitive Cortical Regions

In contrast to steady-state functional connectivity that focuses on connectivity strength, FCV focuses on the variation of the connectivity strength of transient states. It was observed (Fig. 6) that all of low current tDCS, high current tDCS, and after high current tDCS sessions induced a significant increase in FCV between Wernicke’s area and the right frontopolar areas after anodal stimulation of the left Broca’s area. It is known that the left Wernicke’s area is involved in comprehension or understanding of written and spoken language, while its homologous area has a role in the processing and resolution of subordinate meanings of ambiguous words.⁶⁷ Also, the right frontopolar area is associated with memory function⁶⁸ and both Wernicke’s and frontopolar areas belong to high-order language and memory processing networks.⁵⁷^,⁵⁸ Previous studies reported that higher-order brain networks show large variability in connection strength over time.¹⁶^,⁶⁹ Our results suggest that tDCS resulted in larger FCV in these regions, which we interpret as them having a higher degree of flexibility for performing higher-order language processing tasks after stimulation.¹⁶^,⁷⁰ The more variable functional connectivity in these connections potentially reflects tDCS facilitating the emergence of a large-scale network with flexible capability in functional coordination between language and memory systems. Besides, as shown in Tables 3 and 4, low current tDCS showed FCV increases in the same 7/8 Brodmann areas during high current tDCS and the same 6/10 areas during the after high current tDCS session. The majority of Brodmann areas missed (mainly Brodmann area 40 when the seed on channel 34) in the after high current tDCS session were neighboring to the areas that existed during both low current tDCS and high current tDCS (Brodmann area 39 when the seed was on channel 34). Therefore, low current tDCS presented similar results with increased FCV in high-order cognitive cortical regions, thus offering the possibility to serve as a qualitative predictor of patterns in FCV changes induced during the high current tDCS and after high current tDCS conditions.

4.3.

tDCS Enhanced the Small-World Features of the Cortical Network Globally

In contrast with seed-based functional connectivity, graph theory analysis quantified network connectivity changes induced by tDCS across the entire cortical area mapped by fNIRS. tDCS induced an increase in the density of local connections between detector channels across the entire cortical network that was quantified by the network’s small-world features. These network features were the characteristic path length ( $L_{p}$ ), the clustering coefficient ( $C_{p}$ ), global efficiency ( $E_{gb}$ ), and the local efficiency ( $E_{loc}$ ).⁴⁹^,⁷¹^,⁷² Two of these small world properties were enhanced by tDCS: $C_{p}$ became significantly higher for all stimulation sessions, whereas $E_{loc}$ became significantly higher only in the after high current tDCS session. The enhancement of small-world features after tDCS is interpreted as the cortical networks becoming more segregated⁷¹ and clustered with more locally efficient communication and higher fault tolerance.⁵⁹ The graph theory analysis findings were consistent with our seed-based functional connectivity results, where the seed areas were tightly connected with local cortical regions and desynchronized with more distant regions. In addition, changes in graph theory metrics induced by low current tDCS were qualitatively similar to the ones induced by high current tDCS, which indicates that low current tDCS could be used as a tool for studying qualitatively the patterns brought about high current tDCS.

4.4.

Limitations of the Study and Future Work

A few potential limitations should be noted for this study. First, the gender distribution was not even (2 females versus 11 males). We conducted the same analyses when excluding the two female subjects and the vast majority of results did not change, with three exceptions. First, when comparing the no tDCS versus low current tDCS conditions, there was one more channel located in the contralateral BA21 and one channel located in the ipsilateral BA8 with significantly decreased functional connectivity with the seed on left Broca’s area and right Broca’s homologue, respectively ( $p = 0.0334$ and $p = 0.0197$ , respectively). Second, the AUC of $E_{loc}$ had significant changes when comparing the no tDCS versus low current tDCS ( $p = 0.0462$ ) and the no tDCS versus high current tDCS ( $p = 0.0268$ ) conditions for the male population only. Third, when comparing the low current tDCS versus high current tDCS conditions, there was one channel located in the ipsilateral dorsolateral prefrontal BA9 with significantly increased functional connectivity when the seed was on the left Broca’s area ( $p = 0.0378$ ). However, all these changes were marginal and cannot be attributed conclusively to gender effects since we only had two female subjects. Previous MRI studies⁷³^,⁷⁴ reveal significant gender differences within small regions of interest in language networks during semantic language processing, of which the middle temporal gyrus and pars opercularis are at accessible depths for fNIRS too. To reach more certain assertions about the effect of gender on the response to tDCS, it would be important to include more female subjects and larger numbers of subjects overall in a future study. An additional limitation was that even though low current tDCS could qualitatively predict some aspects of changes in functional connectivity patterns induced by high current tDCS, those changes were not linked to meaningful behavioral changes in this study. Therefore, in a future study, the performance of language tasks, such as picture naming accuracy⁸^,¹⁹ and response time,⁹^,¹²^,¹⁹ should be included when low current tDCS and high current tDCS are applied. Finally, the lack of resemblance between connectivity changes seen during low current tDCS and high current tDCS indicates that it is not feasible to use lower currents as a rapid qualitative predictor of patterns seen at higher currents, in the resting state at least. This finding is not surprising since the effect of local connectivity increase is known to only occur at higher currents that can create depolarization of the resting membrane potential and increasing neuronal firing rates.⁵ Nevertheless, the qualitative resemblance of connectivity change patterns seen for low current tDCS versus after high current tDCS, both for steady-state and time-dependent connectivity metrics, shows some promise for the hypothesis that lower currents could be used to predict qualitatively the connectivity change patterns seen after higher current, clinically relevant, stimulation intensities. Nevertheless, further work is needed to verify the validity of this hypothesis.

5. Conclusion

Our work demonstrates the feasibility of using resting-state fNIRS to study cortical network reorganization induced by tDCS in the language processing cortical networks of healthy subjects. Seed-based functional connectivity, graph theory analysis, and time-variant functional connectivity were used to track changes in cortical functional connectivity for these networks. Seed-based connectivity and graph theory analyses revealed increased local and decreased remote functional connectivity induced by tDCS. At the same time, time-variant functional connectivity changes suggested that tDCS increased FCV of remote connections. In addition, low current tDCS produced qualitatively similar patterns of increases in clustering coefficient and FCV, but did not produce patterns of increased local connectivity observed during and immediately after high current tDCS. The qualitative similarity of connectivity patterns seen between the low current tDCS and after high current tDCS conditions hints at the possibility of using multiple electrode placements with low stimulation currents in a single protocol to gain insight into what these changes could be at higher currents for each subject. Nevertheless, the validity of this approach needs to be tested in further work involving stimulation during performance of a task so that connectivity changes could be associated with observable task performance changes.

Disclosures

All of the authors had neither relevant financial or competing interests nor other potential conflicts of interest.

References

1.

F. Fregni et al., “Anodal transcranial direct current stimulation of prefrontal cortex enhances working memory,” Exp. Brain Res., 166 (1), 23 –30 (2005). https://doi.org/10.1007/s00221-005-2334-6 EXBRAP 0014-4819 Google Scholar

2.

T. Z. Kincses et al., “Facilitation of probabilistic classification learning by transcranial direct current stimulation of the prefrontal cortex in the human,” Neuropsychologia, 42 (1), 113 –117 (2004). https://doi.org/10.1016/S0028-3932(03)00124-6 NUPSA6 0028-3932 Google Scholar

3.

R. Sparing et al., “Enhancing language performance with non-invasive brain stimulation—a transcranial direct current stimulation study in healthy humans,” Neuropsychologia, 46 (1), 261 –268 (2008). https://doi.org/10.1016/j.neuropsychologia.2007.07.009 NUPSA6 0028-3932 Google Scholar

4.

B. Elsner et al., “Transcranial direct current stimulation (tDCS) for improving function and activities of daily living in patients after stroke,” Cochrane Database Syst. Rev., 11 CD009645 (2013). https://doi.org/10.1002/14651858.CD009645.pub2 Google Scholar

5.

M. A. Nitsche and W. Paulus, “Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation,” J. Physiol., 527 (3), 633 –639 (2000). https://doi.org/10.1111/tjp.2000.527.issue-3 Google Scholar

6.

L. J. Bindman, O. Lippold and J. Redfearn, “The action of brief polarizing currents on the cerebral cortex of the rat (1) during current flow and (2) in the production of long‐lasting after‐effects,” J. Physiol., 172 (3), 369 –382 (1964). https://doi.org/10.1113/jphysiol.1964.sp007425 Google Scholar

7.

D. P. Purpura and J. G. McMurtry, “Intracellular activities and evoked potential changes during polarization of motor cortex,” J. Neurophysiol., 28 (1), 166 –185 (1965). https://doi.org/10.1152/jn.1965.28.1.166 JONEA4 0022-3077 Google Scholar

8.

J. M. Baker, C. Rorden and J. Fridriksson, “Using transcranial direct-current stimulation to treat stroke patients with aphasia,” Stroke, 41 (6), 1229 –1236 (2010). https://doi.org/10.1161/STROKEAHA.109.576785 SJCCA7 0039-2499 Google Scholar

9.

J. Fridriksson et al., “Transcranial direct current stimulation improves naming reaction time in fluent aphasia,” Stroke, 42 (3), 819 –821 (2011). https://doi.org/10.1161/STROKEAHA.110.600288 SJCCA7 0039-2499 Google Scholar

10.

V. Fiori et al., “Transcranial direct current stimulation improves word retrieval in healthy and nonfluent aphasic subjects,” J. Cognit. Neurosci., 23 (9), 2309 –2323 (2011). https://doi.org/10.1162/jocn.2010.21579 JCONEO 0898-929X Google Scholar

11.

A. Flöel et al., “Noninvasive brain stimulation improves language learning,” J. Cognit. Neurosci., 20 (8), 1415 –1422 (2008). https://doi.org/10.1162/jocn.2008.20098 JCONEO 0898-929X Google Scholar

12.

R. Holland et al., “Speech facilitation by left inferior frontal cortex stimulation,” Curr. Biol., 21 (16), 1403 –1407 (2011). https://doi.org/10.1016/j.cub.2011.07.021 CUBLE2 0960-9822 Google Scholar

13.

M. Meinzer et al., “Electrical brain stimulation improves cognitive performance by modulating functional connectivity and task-specific activation,” J. Neurosci., 32 (5), 1859 –1866 (2012). https://doi.org/10.1523/JNEUROSCI.4812-11.2012 JNRSDS 0270-6474 Google Scholar

14.

P. Marangolo et al., “Bilateral transcranial direct current stimulation language treatment enhances functional connectivity in the left hemisphere: preliminary data from aphasia,” J. Cognit. Neurosci., 28 (5), 724 –738 (2016). https://doi.org/10.1162/jocn_a_00927 JCONEO 0898-929X Google Scholar

15.

Z. Li et al., “Dynamic functional connectivity revealed by resting-state functional near-infrared spectroscopy,” Biomed. Opt. Express, 6 (7), 2337 –2352 (2015). https://doi.org/10.1364/BOE.6.002337 BOEICL 2156-7085 Google Scholar

16.

R. M. Hutchison et al., “Dynamic functional connectivity: promise, issues, and interpretations,” Neuroimage, 80 360 –378 (2013). https://doi.org/10.1016/j.neuroimage.2013.05.079 NEIMEF 1053-8119 Google Scholar

17.

E. Damaraju et al., Static and Dynamic Functional Network Connectivity During Resting State in Schizophrenia, American College of Neuropsychopharmacology, Brentwood, Tennessee (2012). Google Scholar

18.

D. T. Jones et al., “Non-stationarity in the ‘resting brain’s’ modular architecture,” PLoS One, 7 (6), e39731 (2012). https://doi.org/10.1371/journal.pone.0039731 POLNCL 1932-6203 Google Scholar

19.

A. Monti et al., “Improved naming after transcranial direct current stimulation in aphasia,” J. Neurol. Neurosurg. Psychiatry, 79 (4), 451 –453 (2008). https://doi.org/10.1136/jnnp.2007.135277 Google Scholar

20.

L. J. Bindman, O. Lippold and J. Redfearn, “Long-lasting changes in the level of the electrical activity of the cerebral cortex produced by polarizing currents,” Nature, 196 (4854), 584 –585 (1962). https://doi.org/10.1038/196584a0 Google Scholar

21.

B. Khan et al., “Enhancing motor performance improvement by personalizing non-invasive cortical stimulation with concurrent functional near-infrared spectroscopy and multi-modal motor measurements,” Proc. SPIE, 9305 93051A (2015). https://doi.org/10.1117/12.2076590 PSISDG 0277-786X Google Scholar

22.

M. A. Franceschini et al., “Hemodynamic evoked response of the sensorimotor cortex measured noninvasively with near‐infrared optical imaging,” Psychophysiology, 40 (4), 548 –560 (2003). https://doi.org/10.1111/psyp.2003.40.issue-4 PSPHAF 0048-5772 Google Scholar

23.

F. Irani et al., “Functional near infrared spectroscopy (fNIRS): an emerging neuroimaging technology with important applications for the study of brain disorders,” Clin. Neuropsychol., 21 (1), 9 –37 (2007). https://doi.org/10.1080/13854040600910018 CLNEEC 0920-1637 Google Scholar

24.

S. Cutini, S. B. Moro and S. Bisconti, “Functional near infrared optical imaging in cognitive neuroscience: an introductory review,” J. Near Infrared Spectrosc., 20 (1), 75 –92 (2012). https://doi.org/10.1255/jnirs.969 Google Scholar

25.

B. Khan et al., “Functional near-infrared spectroscopy maps cortical plasticity underlying altered motor performance induced by transcranial direct current stimulation,” J. Biomed. Opt., 18 (11), 116003 (2013). https://doi.org/10.1117/1.JBO.18.11.116003 JBOPFO 1083-3668 Google Scholar

26.

A. Villringer, “Understanding functional neuroimaging methods based on neurovascular coupling,” Optical Imaging of Brain Function and Metabolism 2, 177 –193 Springer, New York (1997). Google Scholar

27.

A. C. Merzagora et al., “Prefrontal hemodynamic changes produced by anodal direct current stimulation,” Neuroimage, 49 (3), 2304 –2310 (2010). https://doi.org/10.1016/j.neuroimage.2009.10.044 NEIMEF 1053-8119 Google Scholar

28.

H. Takai et al., “Effect of transcranial direct current stimulation over the primary motor cortex on cerebral blood flow: a time course study using near-infrared spectroscopy,” Oxygen Transport to Tissue XXXVII, 335 –341 Springer, New York (2016). Google Scholar

29.

F. Tian et al., “Test-retest assessment of cortical activation induced by repetitive transcranial magnetic stimulation with brain atlas-guided optical topography,” J. Biomed. Opt., 17 (11), 116020 (2012). https://doi.org/10.1117/1.JBO.17.11.116020 JBOPFO 1083-3668 Google Scholar

30.

J. C. Ye et al., “NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy,” Neuroimage, 44 (2), 428 –447 (2009). https://doi.org/10.1016/j.neuroimage.2008.08.036 NEIMEF 1053-8119 Google Scholar

31.

G. H. Klem et al., “The ten-twenty electrode system of the International Federation,” Electroencephalogr. Clin. Neurophysiol., 52 (3), 3 –6 (1999). Google Scholar

32.

A. Antal et al., “Transcranial direct current stimulation over the primary motor cortex during fMRI,” Neuroimage, 55 (2), 590 –596 (2011). https://doi.org/10.1016/j.neuroimage.2010.11.085 NEIMEF 1053-8119 Google Scholar

33.

T. J. Huppert et al., “HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain,” Appl. Opt., 48 (10), D280 –D298 (2009). https://doi.org/10.1364/AO.48.00D280 APOPAI 0003-6935 Google Scholar

34.

B. R. White et al., “Resting-state functional connectivity in the human brain revealed with diffuse optical tomography,” Neuroimage, 47 (1), 148 –156 (2009). https://doi.org/10.1016/j.neuroimage.2009.03.058 NEIMEF 1053-8119 Google Scholar

35.

B. Biswal et al., “Functional connectivity in the motor cortex of resting human brain using echo‐planar MRI,” Magn. Reson. Med., 34 (4), 537 –541 (1995). https://doi.org/10.1002/(ISSN)1522-2594 MRMEEN 0740-3194 Google Scholar

36.

X. Cui, S. Bray and A. L. Reiss, “Functional near infrared spectroscopy (NIRS) signal improvement based on negative correlation between oxygenated and deoxygenated hemoglobin dynamics,” Neuroimage, 49 (4), 3039 –3046 (2010). https://doi.org/10.1016/j.neuroimage.2009.11.050 NEIMEF 1053-8119 Google Scholar

37.

D. Wildgruber et al., “Identification of emotional intonation evaluated by fMRI,” Neuroimage, 24 (4), 1233 –1241 (2005). https://doi.org/10.1016/j.neuroimage.2004.10.034 NEIMEF 1053-8119 Google Scholar

38.

H. Niu et al., “Resting-state functional connectivity assessed with two diffuse optical tomographic systems,” J. Biomed. Opt., 16 (4), 046006 (2011). https://doi.org/10.1117/1.3561687 JBOPFO 1083-3668 Google Scholar

39.

R. Polanía, M. A. Nitsche and W. Paulus, “Modulating functional connectivity patterns and topological functional organization of the human brain with transcranial direct current stimulation,” Hum. Brain Mapp., 32 (8), 1236 –1249 (2011). https://doi.org/10.1002/hbm.21104 Google Scholar

40.

M. Xia, J. Wang and Y. He, “BrainNet viewer: a network visualization tool for human brain connectomics,” PLoS One, 8 (7), e68910 (2013). https://doi.org/10.1371/journal.pone.0068910 POLNCL 1932-6203 Google Scholar

41.

R. M. Hutchison et al., “Resting‐state networks show dynamic functional connectivity in awake humans and anesthetized macaques,” Hum. Brain Mapp., 34 (9), 2154 –2177 (2013). https://doi.org/10.1002/hbm.v34.9 Google Scholar

42.

A. Kucyi, T. V. Salomons and K. D. Davis, “Mind wandering away from pain dynamically engages antinociceptive and default mode brain networks,” Proc. Natl. Acad. Sci., 110 (46), 18692 –18697 (2013). https://doi.org/10.1073/pnas.1312902110 Google Scholar

43.

H. Niu et al., “Revealing topological organization of human brain functional networks with resting-state functional near infrared spectroscopy,” PLoS One, 7 (9), e45771 (2012). https://doi.org/10.1371/journal.pone.0045771 POLNCL 1932-6203 Google Scholar

44.

J.-R. Ding et al., “Topological fractionation of resting-state networks,” PLoS One, 6 (10), e26596 (2011). https://doi.org/10.1371/journal.pone.0026596 POLNCL 1932-6203 Google Scholar

45.

W. Liao et al., “Altered functional connectivity and small-world in mesial temporal lobe epilepsy,” PLoS One, 5 (1), e8525 (2010). https://doi.org/10.1371/journal.pone.0008525 POLNCL 1932-6203 Google Scholar

46.

M. D. Humphries, K. Gurney and T. J. Prescott, “The brainstem reticular formation is a small-world, not scale-free, network,” Proc. R. Soc. London B: Biol. Sci., 273 (1585), 503 –511 (2006). https://doi.org/10.1098/rspb.2005.3354 Google Scholar

47.

M. Kaiser and C. C. Hilgetag, “Nonoptimal component placement, but short processing paths, due to long-distance projections in neural systems,” PLoS Comput. Biol., 2 (7), e95 (2006). https://doi.org/10.1371/journal.pcbi.0020095 Google Scholar

48.

D. J. Watts and S. H. Strogatz, “Collective dynamics of ‘small-world’networks,” Nature, 393 (6684), 440 –442 (1998). https://doi.org/10.1038/30918 Google Scholar

49.

V. Latora and M. Marchiori, “Efficient behavior of small-world networks,” Phys. Rev. Lett., 87 (19), 198701 (2001). https://doi.org/10.1103/PhysRevLett.87.198701 PRLTAO 0031-9007 Google Scholar

50.

S. Achard and E. Bullmore, “Efficiency and cost of economical brain functional networks,” PLoS Comput. Biol., 3 (2), e17 (2007). https://doi.org/10.1371/journal.pcbi.0030017 Google Scholar

51.

A. S. Dick et al., “Co‐speech gestures influence neural activity in brain regions associated with processing semantic information,” Hum. Brain Mapp., 30 (11), 3509 –3526 (2009). https://doi.org/10.1002/hbm.v30:11 Google Scholar

52.

M. Meinzer et al., “Motor cortex preactivation by standing facilitates word retrieval in aphasia,” Neurorehabilitation Neural Repair, 25 (2), 178 –187 (2011). https://doi.org/10.1177/1545968310376577 Google Scholar

53.

M. Meinzer et al., “Transcranial direct current stimulation of the primary motor cortex improves word-retrieval in older adults,” Front. Aging Neurosci., 6 253 (2014). https://doi.org/10.3389/fnagi.2014.00253 Google Scholar

54.

A. Fertonani et al., “Naming facilitation induced by transcranial direct current stimulation,” Behav. Brain Res., 208 (2), 311 –318 (2010). https://doi.org/10.1016/j.bbr.2009.10.030 BBREDI 0166-4328 Google Scholar

55.

K. Herholz et al., “Individual metabolic anatomy of repeating words demonstrated by MRI-guided positron emission tomography,” Neurosci. Lett., 182 (1), 47 –50 (1994). https://doi.org/10.1016/0304-3940(94)90202-X NELED5 0304-3940 Google Scholar

56.

M. Daianu et al., “Left versus right hemisphere differences in brain connectivity: 4-Tesla HARDI tractography in 569 twins,” in 9th IEEE Int. Symp. on Biomedical Imaging (ISBI ’12), 526 –529 (2012). https://doi.org/10.1109/ISBI.2012.6235601 Google Scholar

57.

D. E. Stark et al., “Regional variation in interhemispheric coordination of intrinsic hemodynamic fluctuations,” J. Neurosci., 28 (51), 13754 –13764 (2008). https://doi.org/10.1523/JNEUROSCI.4544-08.2008 JNRSDS 0270-6474 Google Scholar

58.

N. Ramnani and A. M. Owen, “Anterior prefrontal cortex: insights into function from anatomy and neuroimaging,” Nat. Rev. Neurosci., 5 (3), 184 –194 (2004). https://doi.org/10.1038/nrn1343 NRNAAN 1471-003X Google Scholar

59.

F. D. V. Fallani et al., “Cortical functional connectivity networks in normal and spinal cord injured patients: evaluation by graph analysis,” Hum. Brain Mapp., 28 (12), 1334 –1346 (2007). https://doi.org/10.1002/(ISSN)1097-0193 Google Scholar

60.

E. Sivan, H. Parnas and D. Dolev, “Fault tolerance in the cardiac ganglion of the lobster,” Biol. Cybern., 81 (1), 11 –23 (1999). https://doi.org/10.1007/s004220050541 BICYAF 0340-1200 Google Scholar

61.

H. L. Filmer, P. E. Dux and J. B. Mattingley, “Applications of transcranial direct current stimulation for understanding brain function,” Trends Neurosci., 37 (12), 742 –753 (2014). https://doi.org/10.1016/j.tins.2014.08.003 TNSCDR 0166-2236 Google Scholar

62.

M. A. Nitsche et al., “Transcranial direct current stimulation: state of the art 2008,” Brain Stimul., 1 (3), 206 –223 (2008). https://doi.org/10.1016/j.brs.2008.06.004 Google Scholar

63.

R. Polanía et al., “Introducing graph theory to track for neuroplastic alterations in the resting human brain: a transcranial direct current stimulation study,” Neuroimage, 54 (3), 2287 –2296 (2011). https://doi.org/10.1016/j.neuroimage.2010.09.085 NEIMEF 1053-8119 Google Scholar

64.

Z. Cattaneo, A. Pisoni and C. Papagno, “Transcranial direct current stimulation over Broca’s region improves phonemic and semantic fluency in healthy individuals,” Neuroscience, 183 64 –70 (2011). https://doi.org/10.1016/j.neuroscience.2011.03.058 Google Scholar

65.

P. Marangolo et al., “tDCS over the left inferior frontal cortex improves speech production in aphasia,” Front. Hum. Neurosci., 7 539 (2013). https://doi.org/10.3389/fnhum.2013.00539 Google Scholar

66.

P. Marangolo et al., “Electrical stimulation over the left inferior frontal gyrus (IFG) determines long-term effects in the recovery of speech apraxia in three chronic aphasics,” Behav. Brain Res., 225 (2), 498 –504 (2011). https://doi.org/10.1016/j.bbr.2011.08.008 BBREDI 0166-4328 Google Scholar

67.

Y. Harpaz, Y. Levkovitz and M. Lavidor, “Lexical ambiguity resolution in Wernicke’s area and its right homologue,” Cortex, 45 (9), 1097 –1103 (2009). https://doi.org/10.1016/j.cortex.2009.01.002 Google Scholar

68.

A. Costa et al., “The right frontopolar cortex is involved in visual-spatial prospective memory,” PLoS One, 8 (2), e56039 (2013). https://doi.org/10.1371/journal.pone.0056039 POLNCL 1932-6203 Google Scholar

69.

J. Gonzalez-Castillo et al., Detection of Task Transitions on 45 mins Long Continuous Multi-Task Runs Using Whole Brain Connectivity, Organization for Human Brain Mapping, Beijing, China (2012). Google Scholar

70.

A. M. Hermundstad et al., “Structural foundations of resting-state and task-based functional connectivity in the human brain,” Proc. Natl. Acad. Sci. U. S. A., 110 (15), 6169 –6174 (2013). https://doi.org/10.1073/pnas.1219562110 Google Scholar

71.

M. Rubinov and O. Sporns, “Complex network measures of brain connectivity: uses and interpretations,” Neuroimage, 52 (3), 1059 –1069 (2010). https://doi.org/10.1016/j.neuroimage.2009.10.003 NEIMEF 1053-8119 Google Scholar

72.

D. J. Smit et al., “Heritability of ‘small‐world’ networks in the brain: a graph theoretical analysis of resting‐state EEG functional connectivity,” Hum. Brain Mapp., 29 (12), 1368 –1378 (2008). https://doi.org/10.1002/hbm.v29:12 Google Scholar

73.

G. S. Harrington and S. T. Farias, “Sex differences in language processing: functional MRI methodological considerations,” J. Magn. Reson. Imaging, 27 (6), 1221 –1228 (2008). https://doi.org/10.1002/jmri.21374 Google Scholar

74.

L. C. Baxter et al., “Sex differences in semantic language processing: a functional MRI study,” Brain Lang., 84 (2), 264 –272 (2003). https://doi.org/10.1016/S0093-934X(02)00549-7 Google Scholar

Biography

Jianwei Cao received her BS and MS degrees in optical engineering from Hefei University of Technology, Hefei, China, in 2008 and 2011, respectively. She worked with George Alexandrakis’s group in bioengineering from 2013 to 2017 and received her PhD degree in bioengineering from the University of Texas at Arlington, Arlington, Texas, United States, in 2017. Her research mainly focuses on functional near infrared spectroscopy, brain imaging, and biomedical signal processing.

Hanli Liu is a full professor of bioengineering at the University of Texas at Arlington. She received her MS and PhD degrees from Wake Forest University in physics, followed by postdoctoral training at the University of Pennsylvania in tissue optics. Her current expertise lies in the field of near-infrared spectroscopy of tissues, optical sensing for cancer detection, and diffuse optical tomography for functional brain imaging, all of which are related to clinical applications.

George Alexandrakis did his undergraduate studies in physics at Oxford University, United Kingdom, and his graduate studies in medical physics at McMaster University, Canada. He was a postdoctoral fellow at Massachusetts General Hospital/Harvard Medical School and at UCLA. As faculty at the University of Texas at Arlington’s Bioengineering Department since 2006, his research is focused on the development of optical imaging methods for biomedical applications.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Jianwei Cao, Hanli Liu, and George Alexandrakis "Modulating the resting-state functional connectivity patterns of language processing areas in the human brain with anodal transcranial direct current stimulation applied over the Broca’s area," Neurophotonics 5(2), 025002 (27 February 2018). https://doi.org/10.1117/1.NPh.5.2.025002

Received: 3 October 2017; Accepted: 5 February 2018; Published: 27 February 2018

Access the abstract

JOURNAL ARTICLE
12 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 12 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Brain

Sensors

Modulation

Lab on a chip

Neurophotonics

Electrodes

Shape memory alloys

1.

Introduction

2.

Methods and Materials

2.1.

Subjects

2.2.

fNIRS Imaging Setup Combined with tDCS

Fig. 1

Fig. 2

2.3.

Protocol Design

Fig. 3

2.4.

Data Preprocessing

2.5.

Data Analysis

2.5.1.

Seed-based functional connectivity

2.5.2.

Time-variant functional connectivity

2.5.3.

Graph theory analysis

3.

Results

3.1.

Functional Connectivity Analysis Using the Stimulated Cortical Region and Contralateral Cortical Region as the Seed

3.1.1.

Using the stimulated left Broca’s area as the seed

Fig. 4

3.1.2.

Using the right Broca’s homologue area as the seed

Fig. 5

Table 1

Table 2

3.2.

Time-Variant Functional Connectivity using the Stimulated Left Broca’s Area as Seed

Fig. 6

Table 3

Table 4

3.3.

Graph Theory Analysis

Fig. 7

4.

Discussion

4.1.

tDCS Increased Steady-State Functional Connectivity in the Immediate Neighborhood of the Stimulation Location

4.2.

tDCS Increased Functional Connectivity Variability with High-Order Cognitive Cortical Regions

4.3.

tDCS Enhanced the Small-World Features of the Cortical Network Globally

4.4.

Limitations of the Study and Future Work

5.

Conclusion

Disclosures

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years