BA4 labelling and registration to the functional space.

<p>An example of the quality of the cortical segmentation and labelling of BA4 is depicted, which was performed with Freesurfer in one of the participants. Panels <b>A</b> and <b>B</b> show the BA4 label over the pial and inflated brain surface respectively. Panel <b>C</b> shows the left BA4 label registered in the original T1 image. The red line shows the portion of the cortical segmentation performed by Freesurfer that corresponds to the left BA4 label (in this image left is right). Notice the tDCS electrode over the scalp of the subject is positioned over the central sulcus. Panel <b>D</b> shows the registration of the cortical segmentation and BA4 label in the native functional space of the same subject (first row: first 14 images of the native functional space. Second row: cortical segmentation registered from Freesurfer to the native functional space was overlapped. Third row: BA4 label (red) registered from Freesurfer space to the native functional space). For visualization purposes, we show the approximate location of the tDCS electrode over left BA4 in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone.0030971.s001" target="_blank">figure S1</a>.</p

Dependency of the tDCS-induced effects on baseline functional architecture.

<p>Panel A shows that the effect of nodal clustering coefficient (<i>C</i>) increase found in the cluster of <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone-0030971-g003" target="_blank">Figure 3F</a> strongly correlated with baseline <i>C</i> (P = 0.0051; R² = 0.46). Panel B shows that the positive decrease found in the characteristic path length (<i>L</i>) maps that was found after anodal tDCS (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone-0030971-g003" target="_blank">Figure 3H</a>) also has a positive correlation with the baseline <i>L</i> (P = 0.002; R² = 0.51).</p

Global network metrics.

<p>Shown are the results of the global network parameters that were calculated in the present study (mean connectivity degree (<i>K</i>) and the small-world parameters gamma, lambda and sigma <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone.0030971-Sporns2" target="_blank">[32]</a>) calculated at each threshold T (0.25–0.35 in increasing steps of 0.002) before and after each Sham (A), Cathodal (B) and Anodal (C) tDCS. The second row - second column of each panel show that the approximate number of nodes (M1 voxels) of each undirected graph was ∼470. As expected, the mean connectivity degree monotonically decreases as <i>T</i> increases. M1 has salient small-world properties i.e. <i>lambda≈1</i>, <i>gamma≫1</i>, thus <i>sigma≫1 </i><a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone.0030971-Sporns2" target="_blank">[32]</a>, <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#pone.0030971-Achard1" target="_blank">[41]</a>. No significant differences were observed for any of the network metrics before and after each of the tDCS sessions (P>0.05 paired two-tailed t-tests). Error bars represent the s.e.m.</p

Highly efficient nodes within M1.

<p>The average <i>L</i> maps for all subjects and all before-tDCS fMRI scans were averaged. Nodes that showed the highest <i>L_rand/Li</i> values (i.e. nodes that communicate more efficiently within M1) were mapped over the flattened BA4. As an exploratory threshold we used the 15% of the voxels that showed the highest <i>L_rand/Li</i> values.</p

Graph parameter statistics at the BA4 cortical surface.

<p>(A) Shown is the flattening of the left BA4 (green labelled region) obtained from the left hemisphere surface average subject, which was used to project the statistical maps. Panels B to D show the ANOVA for the interaction effects (time*stimulation) Montecarlo cluster corrected at p<0.05 for the nodal connectivity degree maps (B), clustering coefficient maps (C) and the characteristic path length maps (D). Panels E to H show post hoc paired t-tests for the following contrasts: (E) After_Cathodal – After_Sham in the clustering coefficient maps; (F) After_Anodal – Before_Anodal in the clustering coefficient maps; (G) After_Anodal – After Sham in the characteristic path length maps; (H) After_Anodal – Before_Anodal in the characteristic path length maps. Notice that the <i>L</i> maps were <i>L_rand/Li</i> normalized, which means that the values of <i>L</i> in the significant cluster are lower after stimulation (see <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0030971#s2" target="_blank">methods</a> section).</p

Peak voxel intensity and coordinates for brain regions showing performance-related activity.

<p>Peak voxel intensity and coordinates for brain regions showing performance-related activity.</p

Regions of decreased activity for hf-tRNS.

<p>Contrast <i>sham- Hfreq</i> (A) revealed changes in the left frontal cortex. Contrast <i>Lfreq-Hfreq</i> (B) revealed additional changes in right frontal cortex and precuneous.</p

Peak voxel intensity and coordinates for brain regions showing motor task-related activity.

<p>Peak voxel intensity and coordinates for brain regions showing motor task-related activity.</p

Performance related decrease of brain activity during and after stimulation.

<p>Activity decreased with time for contrast <i>run1-run2</i> in the paracingulate gyrus, superior frontal gyrus, thalamus and hippocampus. (Z>3, P<0.05, corrected).</p

Changes in tracking error relative to the first trial.

<p>Shaded area corresponds to the stimulation period. Fifth Polynomial trendlines are superimposed on the data for easier visualization.</p