Differences in high-definition transcranial direct current stimulation over the motor hotspot versus the premotor cortex on motor network excitability

The effectiveness of transcranial direct current stimulation (tDCS) placed over the motor hotspot (thought to represent the primary motor cortex (M1)) to modulate motor network excitability is highly variable. The premotor cortex—particularly the dorsal premotor cortex (PMd)—may be a promising alternative target to reliably modulate motor excitability, as it influences motor control across multiple pathways, one independent of M1 and one with direct connections to M1. This double-blind, placebo-controlled preliminary study aimed to differentially excite motor and premotor regions using high-definition tDCS (HD-tDCS) with concurrent functional magnetic resonance imaging (fMRI). HD-tDCS applied over either the motor hotspot or the premotor cortex demonstrated high inter-individual variability in changes on cortical motor excitability. However, HD-tDCS over the premotor cortex led to a higher number of responders and greater changes in local fMRI-based complexity than HD-tDCS over the motor hotspot. Furthermore, an analysis of individual motor hotspot anatomical locations revealed that, in more than half of the participants, the motor hotspot is not located over anatomical M1 boundaries, despite using a canonical definition of the motor hotspot. This heterogeneity in stimulation site may contribute to the variability of tDCS results. Altogether, these preliminary findings provide new considerations to enhance tDCS reliability.

some individuals will show some level of change in cortical excitability, even without stimulation. We also found considerable inter-individual variability in MEP changes following HD-tDCS across all the 3 groups (Fig. 1). The number of participants with increased MEP amplitudes post-HD-tDCS showed a trending difference across the 3 groups (χ 2 (2) = 5.497, p = 0.06). Finally, as the boundary between responders and nonresponders using canonical definitions is very narrow 12,37 , we also explored different, more robust methods to classify responders/non-responders, which resulted in similar results. This additional analysis is summarized in Supplementary Data 1, Supplementary Table 2. Group level changes in MEP amplitude: Then, we examined whether HD-tDCS induced significant changes in cortical motor excitability across the 3 groups by comparing the MEP amplitudes before and after HD-tDCS (Table 1, Fig. 2).
A two-way RM-ANOVA was performed on MEP amplitude, with "Time" (pre-HD-tDCS, post-HD-tDCS) and "Group" (Motor Hotspot, Premotor, Sham) as factors. There were no significant interaction or group effects Resting-state fMRI. Resting-state functional connectivity. Next, we checked if HD-tDCS induced changes in rsfMRI motor network connectivity across the 3 groups by examining whether there were overall changes in the functional organization of the motor network. We defined four regions of interest (ROIs) as the left M1, right M1, left PMd and right PMd, using a meta-analysis from Hardwick and collaborators 38 to define the location of the ROIs in the left hemisphere and then flipped them to the right hemisphere. We compared cross-correlation rankings for the six ROI-to-ROI pairs (1)   www.nature.com/scientificreports www.nature.com/scientificreports/ The global organization of the motor network did not change after HD-tDCS in any of the groups, as demonstrated by a strong rank correlation within each group between pre and post-HD-tDCS measurements (Motor Hotspot: n = 15, r = 0.88, p = 0.03, Premotor: n = 15, r = 0.88, p = 0.03, Sham: n = 16, r = 0.94, p = 0.02). Moreover, there were also no significant differences in correlation ranks (Fisher's Z test 39 ) between groups (Motor Hotspot > Premotor: z = 0.0, p = 1; Motor Hotspot > Sham: z = 0.78, p = 0.43, Premotor > Sham: z = 0.78, p = 0.43). This suggests that there was no global reorganization of the connectivity strengths in the motor network following HD-tDCS over either of the targets. We also found no significant changes post-HD-tDCS in the connectivity each ROI-to-ROI pair in any of the three groups (ANOVA, n = 46, thresholded at qFDR = 0.05; CONN toolbox). Additional exploratory analyses focusing on ROI-to-ROI comparisons and ICA network comparisons at a less stringent threshold (cluster-level corrected p < 0.01) are reported in Supplementary  Multi-scale entropy of resting-state fMRI. Then, we checked if HD-tDCS induced changes in motor network complexity by examining HD-tDCS-associated changes in MSE across the 3 groups (Fig. 3).
RM-ANOVAs with "Time" (pre-HD-tDCS, post-HD-tDCS) and "Group" ( Comparison of stimulation site location. Finally, given the high inter-individual variability of the results, we performed an additional post-hoc exploratory analysis to examine if the individual functional location of the left motor hotspot and the left premotor cortex matched the conventional/anatomical definitions for these regions (Fig. 4), and whether this could explain variability in HD-tDCS effects. www.nature.com/scientificreports www.nature.com/scientificreports/ In the Motor Hotspot group, the average minimal Hausdorff 's distance between individual's cortex-projected stimulated region (i.e., the motor hotspot) and the standard M1 ROI was 3.79 ± 4.64 mm, and with the standard PMd ROI, it was 3.4 6 ± 4.89 mm. The motor hotspot site overlapped with the standard anatomical M1 ROI in 6 participants and with the standard anatomical PMd ROI in 6 participants. In the remaining 3 participants, the stimulation site was localized between the standard M1 and PMd ROIs.
In the Premotor group, the average Hausdorff 's distance between the anatomical location of the individual stimulation (i.e., premotor cortex stimulation) site and the standard PMd was 4.93 ± 4.80 mm, and with the standard M1, it was 15.44 ± 4.78 mm. The premotor cortex stimulation site overlapped with the standard PMd ROI in 5 participants and between the standard M1 and PMd ROIs in 2 participants. In the other 8 participants, the anatomical location of the stimulation site was localized slightly anterior to the standard PMd ROI, but still within the anatomical (atlas-based) boundaries of PMd.
Given the wide variability in the anatomical location of the stimulation site in the Motor Hotspot group, we reassigned the participants to groups based on the anatomical location of their stimulation site (see Supplementary Data 4, Supplementary Fig. 3). We confirmed our previous results, finding again a significant  These results suggest that there is large inter-individual variability in the left motor hotspot and the premotor cortex anatomical localization using canonical methods employed in a majority of tDCS studies (e.g., TMS-defined motor hotspot [6][7][8][9][10]12,[40][41][42][43][44][45][46] and PMd defined as 2.5 cm anterior to the motor hotspot [47][48][49][50] ). In particular, many individuals in the Motor Hotspot group had stimulation sites that were bordering or overlapping with the standard PMd ROIs, or somewhere between the standard M1 and PMd ROIs. This heterogeneity in the stimulation site could contribute to the variability of responsiveness to HD-tDCS.
We, therefore, compared the location of the stimulation sites between responders and non-responders. As shown in Fig. 5, more than half of the responders, including participants in both the Motor Hotspot and Premotor groups, received stimulation over the same area inside PMd anatomical boundaries. We found that 10 responders had an overlapping stimulation site within anatomical PMd region, regardless of group allocation (see Supplementary Fig. 4), while the anatomical location of the stimulation site in non-responders across groups was more widespread.

Discussion
In this preliminary study, we performed a systematic, double-blind investigation using concurrent HD-tDCS and fMRI to examine changes in brain network connectivity and complexity after HD-tDCS stimulation of either the left motor hotspot (assumed to reflect the hand knob location in M1 3 ) or the left premotor cortex. To our knowledge, this is the first known group study using concurrent HD-tDCS and fMRI to compare the left motor hotspot and left premotor cortex stimulation. We examined changes in cortical motor excitability using both TMS and fMRI after HD-tDCS. Finally, we explored inter-individual variability in the stimulation sites.
Overall, we observed that HD-tDCS applied over either the left motor hotspot or left premotor cortex produced highly variable changes across individuals. Due to the large inter-individual variability in response to tDCS, the magnitude of changes in cortical excitability was not significantly different between the 3 groups. However, despite this, we found that slightly more participants showed an increased MEP amplitude following HD-tDCS over left premotor cortex than HD-tDCS over the left motor hotspot or sham, and when comparing intragroup excitability changes, we found that MEP amplitude increased significantly after HD-tDCS only in the Premotor group. The slight increase of cortical excitability in the Premotor group was also accompanied by modulations in motor network complexity.
Finally, we showed that there was significant variability in the location of the anatomical stimulation sites compared to the expected stimulation sites based on canonical electrode placement -a factor which may have introduced variability into previous brain stimulation studies. In addition, the highest and most consistent number of responders was observed when stimulation occurred over left anatomical PMd. However, there was still large inter-individual variability in the stimulated anatomical location of responders, and some individuals in the non-responder group also were stimulated within PMd, suggesting that there is likely not a single anatomical target across all participants that would effectively increase excitability.
Altogether, these preliminary results suggest that HD-tDCS applied over premotor cortex may lead to small but more reliable changes in cortical excitability than over the motor hotspot. As previously observed, there was large variability in the MEP measurement even in the Sham group, suggesting that this measure is highly variable over time and across individuals. However, in line with our initial hypothesis, we observed a non-statistically significant trend of a more reliable increase of MEP amplitude in the Premotor group, resulting in 75% responders compared to only 50% responders in the Motor Hotspot group (and 31% in the Sham group). Moreover, only the Premotor group showed a slightly significant increase in MEP amplitude post-tDCS. These results support previous findings 12 that the motor hotspot may not be the best target to consistently enhance cortical excitability.
These findings may also help to explain previous findings with a standard tDCS montage. Using a standard tDCS montage, the distribution of the electrical current is quite diffuse and thus, when M1 is targeted, at least a part of PMd 27 is also likely stimulated 51 . This has made distinguishing the effects associated with tDCS over M1 versus PMd difficult. The present preliminary work aimed to differentially target M1 versus PMd by using a highly focal montage in which the peak of the stimulation was concentrated over one small center electrode (constrained by a ring). The preliminary work presented in this paper suggests a more reliable impact of PMd than M1 on motor excitability changes. This is in line with previous observations that standard tDCS montages over the motor hotspot actually show strong involvement of PMd, which may impact the observed effects 52 . Future research should examine the potential role of the PMd in both standard and high-definition stimulation montages.
In addition, HD-tDCS over both the motor hotspot and premotor cortex led to no differences in motor network functional connectivity. Although previous work has suggested that tDCS may induce changes in functional connectivity of the motor network 53,54 , our analyses showed no significant changes. Our exploratory analyses, at a lower threshold, presented in Supplementary Data 2 and 3, Supplementary Figs. 1 and 2, and Supplementary Discussion, demonstrate a trend towards modifications in functional connectivity after HD-tDCS, but these modifications are only significant at lower statistical threshold. It is possible that the short duration of the HD-tDCS session (7 minutes) and low amplitude of current (1 mA) might not have been sufficient to induce reliable and significantly different changes in the functional connectivity of the motor network across individuals, which has been previously observed after 10 to 20 minutes of stimulation 53,54 .
On the other hand, we did find changes in multiscale entropy in left and right PMd following HD-tDCS over left premotor cortex. MSE, a method for measuring the complexity of finite length time series, quantifies the regularity of a time series by calculating the complexity of a signal at different time scales 55 . This measure reflects the dynamics of the neural networks 33 . Previous work has shown that increasing brain complexity is associated with brain adaptation 34,35 . Importantly, MSE is thought to measure fine-grained changes in brain function following tDCS 34,36,56 . Despite large variability in the baseline MSE measurements, HD-tDCS over premotor cortex, as hypothesized, led to significant and more reliable changes in MSE than HD-tDCS over the motor hotspot or than Sham HD-tDCS. The MSE changes found in both hemispheres could reflect the expected HD-tDCS associated changes in cortical excitability (i.e., an enhancement of the excitability in the stimulated hemisphere associated with a decrease of the excitability in the contralateral hemisphere 57,58 ). In line with previous studies 35 , MSE and functional connectivity seemed to be complementary metrics: MSE may be more sensitive and able to capture local changes in the motor network whereas traditional rs-fMRI functional connectivity would capture more robust changes between distant brain areas.
In the present study, we also observed high inter-individual variability in HD-tDCS-induced changes in cortical excitability. This variability is similar to what has been recently reported, with up to 50% of participants showing no response to tDCS 11,46,59 . Here, we found a modest enhancement of the number of participants who responded to HD-tDCS when applied over premotor cortex. However, when we reallocated individuals based on their anatomical stimulation location, we still found a consistent increase in MEP amplitude only for the Premotor group (for those that received HD-tDCS over premotor cortex). We also consistently found the PMd to be the most-stimulated anatomical location in responders. This suggests that stimulating the premotor cortex to modulate the motor network might help reduce inter-individual variability compared to stimulating the motor hotspot, as long as it is actually the PMd that is being stimulated. These preliminary findings reflect an initial attempt to explore the potentially greater impact of PMd tDCS over M1 tDCS on cortical excitability; however, this work requires larger samples to confirm these findings. Finally, our results suggest there is no single brain area that would ensure a positive response to tDCS across all participants. Although more than half of the participants who responded to tDCS received stimulation over an overlapping PMd spot, for the other half, the stimulation site was somewhere between the anatomical M1 and PMd. Moreover, in the participants who did not respond to tDCS, the stimulation site was all over the motor network, including in PMd, with no strong overlap across participants. Again, this suggests that more than the absolute location, differences in brain structure and function 60 should be taken into account to determine the responsiveness of each individual to tDCS. Future studies should determine which stimulation parameters (neural target, nature of stimulation, duration, etc.) are the most effective to induce changes in motor network excitability and function at the individual level instead of following a one-size-fits-all method. www.nature.com/scientificreports www.nature.com/scientificreports/ The current study also has several limitations. First, the size of our groups (n = 15-16 individuals per group) is relatively small, especially given the large variability within each group. Thus, confirming these results in a larger sample is required. A second limitation is that the current Sham group included half motor hotspot and half the premotor cortex electrode positions, which showed no significant differences with one another, and which we, therefore, combined into a single Sham group. However, to provide proper control for each experimental group, a matched full Sham group for each experimental group should be used in the future. Third, we used a short 7-minute stimulation duration at 1 mA current, which may not have been enough to induce a strong effect, particularly in functional connectivity changes. Moreover, the cortical excitability measurements post-stimulation were measured around 15 minutes after the end of the stimulation (less than 10 min after the end of the third resting-state fMRI scan). This time delay could have impacted the cortical excitability measures. However, as noted previously, this stimulation duration is in line with previous studies showing effects of 5-7 minutes of anodal tDCS on cortex excitability (for up to 30 min after the end of the stimulation) [61][62][63] . In addition, this study demonstrated that stimulating a brain region other than the motor hotspot (in this case, the premotor cortex) even for a short period can lead to an enhancement of cortical excitability and complexity, suggesting the duration was long enough to elicit some changes. However, future work should examine both cortical excitability and whole-brain network changes following longer durations or higher intensities of HD-tDCS. Forth, when turned on, the MRI-compatible HD-tDCS introduces noise into the BOLD signal with some undesirable, hard to remove, artifacts, that are even been observed in a cadaver 64 . Therefore, we were not able to analyze rs-fMRI data during brain stimulation-only immediately before and after. However, it would be very useful in future research to examine also BOLD signal changes during brain stimulation. Fifth, individual anatomy could have influenced the tDCS observed effect, as our anatomical findings suggest. Future studies that more carefully target stimulation based on individual anatomy, using each individual's own anatomical brain scan with neuronavigation, would help address this issue. In addition, the present study does not include electrical field modeling, which is beyond the scope of the current investigation. However, in future studies, the impact of the anatomical properties of the stimulated area in the non-normalized brain of each individual could be used to explore the deeper impacts of the anatomical stimulation location on motor network excitability. Finally, we used a parallel study design, as opposed to a cross-over design. Each participant received HD-tDCS over only one neural target. Consequently, we cannot confirm that the participants who responded to HD-tDCS over premotor cortex would have not also responded to HD-tDCS over the motor hotspot, or vice-versa. Future studies involving a larger sample size and a cross-over design are needed to answer this specific question.

Conclusion
To our knowledge, this preliminary work is the first multi-participant study employing HD-tDCS with fMRI in a double-blind design with multiple stimulation sites. Using this design, we investigated the neurophysiological effects of tDCS on the motor network immediately after stimulation. We compared HD-tDCS over the motor hotspot versus the premotor cortex and showed that the premotor cortex, and especially PMd, may be a potential alternative neural target to modulate motor network excitability and neural complexity with less variability compared to HD-tDCS over the motor hotspot. These results also suggest that multiscale entropy may be a sensitive measure of changes in neuronal excitability following noninvasive brain stimulation within a single region and may be able to detect changes not observed with standard functional connectivity measures. Finally, we demonstrate that conventional methods for defining tDCS targets may be influenced by large inter-individual variability in the anatomy of the motor hotspot location. Future work should carefully consider and account for the anatomical, versus functional, stimulation locations. Altogether, these preliminary results open new considerations for improving the efficacy of tDCS on modulating the motor network.

Material and Methods
Population. The protocol was approved by the University of Southern California Institutional Review Board.
All subjects gave written informed consent in accordance with the 1964 Declaration of Helsinki. The inclusion criteria were the following: (1) being a healthy volunteer aged 18-50 years, (2) being right-handed. Exclusion criteria were any contra-indications to MRI or brain stimulation.

Study design.
The study involved a single three-hour visit which included the following 3 sessions (see Fig. 6): (1) a pre-MRI TMS session, including the localization of the left motor hotspot and the left premotor cortex as well as the recording of neurophysiological measurements, (2) a concurrent HD-tDCS/MRI session, including the placement of the HD-tDCS cap, acquisition of structural MRI sequences, a pre-HD-tDCS rs-fMRI scan, the double-blinded HD-tDCS stimulation with concurrent acquisition of rs-fMRI scan, as well as a post-HD-tDCS rs-fMRI scan, and (3) a post-MRI TMS session to measure changes in neurophysiological measurements following HD-tDCS. To our knowledge, this is the first known randomized, double-blinded group study using a concurrent HD-tDCS and fMRI paradigm in humans. The safety of these concurrent designs was established in a pilot study using both phantoms and a single participant 65 . They also examined the tolerability of their specific double ring HD-tDCS montage in a behavioral session (without fMRI) in 30 participants.
Pre-MRI TMS session. Briefly, we used single-and paired-pulse transcranial magnetic stimulation (TMS) (Magstim 200² device with BiStim module; Magstim Inc., UK) with a figure-of-eight coil, along with the Brainsight neuronavigation system (Rogue Resolutions Ltd, UK) and surface electrodes to record muscle activity (electromyography, or EMG; Delsys Trigno wireless sensors, Delsys Inc, MA). To localize the motor hotspot, we recorded EMG at the right first dorsal interosseous (FDI). We then recorded pre-HD-tDCS cortical excitability, including each individual's RMT and single-pulse motor-evoked potentials (MEPs) with a test stimulus at the intensity needed to generate an MEP of 0.75 mV peak-to-peak amplitude. Notably, this test stimulus is in www.nature.com/scientificreports www.nature.com/scientificreports/ accordance with previous studies showing a range of test stimuli between 0.2 mV and 4mV 66,67 . We also measured short-interval intracortical inhibition (SICI) using a paired-pulse TMS paradigm (see Brain Stimulation Intervention section below for details on the TMS paradigms). MEPs were recording using Signal 6.05 software (Cambridge Electronic Design Ltd., UK).
MRI session. Immediately following the first TMS session, participants were fitted with the HD-tDCS cap and began an MRI session including the acquisition of an anatomical image (T1w MPRAGE), and three 7-minute baseline rs-fMRI. Sham or HD-tDCS were applied during the second rs-fMRI. An additional experimenter, who was unblinded to the stimulation condition (real or sham), turned the tDCS stimulation device on and off. In this way, the primary experimenter and participant remained blinded to the stimulation condition. During the 3 rs-fMRI runs, participants were asked to look at a fixation cross and to not focus on any specific thoughts. The instructions were repeated at the beginning of each rs-fMRI run. Before the second rs-fMRI run, we informed the participants that the HD-tDCS stimulation was about to start, and we checked with them after the ramp-up to see if they were fine with the associated scalp sensations. Before the third rs-fMRI run, we let them know that the HD-tDCS had ended and that it would be the final rs-fMRI run.
Post-MRI TMS session. Immediately following the MRI session, the participants took part in a second TMS session to quantify the HD-tDCS-induced changes on cortical excitability. The cortical excitability measurements of this second TMS session started less than 10 minutes after the end of the third resting-state fMRI scan.
We defined the "responders" subgroup as participants in the Motor Hotspot or Premotor group who showed an increase in MEP amplitude following HD-tDCS [(averaged MEP amplitude over the 20 trials post-tDCS)/ (averaged MEP amplitude over the 20 trials pre-tDCS) > 1)]. The "non-responders" subgroup refers to participants in the Motor Hotspot or Premotor group who demonstrated a decrease or no change in MEP amplitude following HD-tDCS [(averaged MEP amplitude over the 20 trials post-tDCS)/(averaged MEP amplitude over the 20 trials pre-tDCS) ≤ 1)]. We also examined additional definitions of responders (see Supplementary Materials).

Brain stimulation intervention. TMS paradigm.
We defined the left motor hotspot using a canonical definition of the spot at which the maximal MEP peak-to-peak amplitude was recorded in the FDI muscle using a single TMS pulse at a given suprathreshold intensity (starting at 50% of the maximum stimulator output) [40][41][42][43][44][45] . We initially positioned the coil at the C3 location and moved around until finding a location that induced the largest MEP amplitude for a fixed intensity [68][69][70] . At this location, we defined the RMT as the intensity of the stimulator that induced an MEP of at least 50 µV peak-to-peak amplitude in 5 out of 10 trials 71 and the TS 0.75mV peak-to-peak amplitude in 5 out of 10 trials 66,67 .
Then, we acquired a MEP-SICI paradigm which consisted of 40 trials in a randomized order with a 7 s ± 10% inter-trial duration. This paradigm included 20 single-pulse trials using the TS 0.75mV (the peak-to-peak amplitude of these trials will be referred to as MEP amplitude), and 20 paired-pulse trials using a conditioning stimulus (CT) at 80% of the RMT 2.5 ms before the TS 0.75mV (HD-tDCS did not impact the amplitude of the SICI trials, for more details see Supplementary Data 5).
For 6 participants (2 in each group), we were unable to reach the intensity needed to elicit an MEP of 0.75 mV, or even 0.5 mV. In these cases, we used the RMT as the test stimulus to acquire MEP amplitude. For one additional participant in the Motor Hotspot group, we were only able to localize the motor hotspot but did not acquire any other neurophysiological measurements due to the participant's reported discomfort. This participant was included only in the rs-fMRI analyses.
HD-tDCS experiment. We used the Soterix MRI-compatible HD-tDCS system (4 × 1 ring configuration) to deliver tDCS (1 mA) for 7 minutes during an fMRI session (see Fig. 6). We used 5 small gel-filled electrodes (base diameter: 2.4 cm, gel-skin contact area of ~4.5 cm 2 ) positioned in adjacent electrode locations within the cap, and we placed www.nature.com/scientificreports www.nature.com/scientificreports/ the center electrode on the cap over the location of the neural target (as defined below). Each of the 4 cathodal ring electrodes delivered a -0.25 mA current leading to a 1 mA anodal stimulation over the center electrode.
Localization of the HD-tDCS targets. For the Motor Hotspot group, we positioned the center electrode over the left motor hotspot. For the Premotor group, we positioned the center electrode 2.5 cm anteriorly to the motor hotspot, which is thought to represent PMd [72][73][74][75][76] . The use of this MRI-compatible device also gave us the unique opportunity to obtain an anatomical MRI with the exact position of the electrodes on the participant's own head. We therefore also report the location of individual stimulation sites across participants in the Results section.
Resting-state fMRI processing. We analyzed the rs-fMRI data using SPM12, the CONN toolbox (functional connectivity toolbox 77 ; http://www.nitrc.org/projects/conn) and the Complexity Toolbox (http://loft-lab.org/ index-5.html). The preprocessing steps included slice-timing correction, motion realignment, noise correction using white matter, CSF and motion parameters as regressors, and band-pass filtering (0.01-0.1 Hz). We also performed co-registration between the functional scans and 3D-T1w MPRAGE scans of each participant. Finally, we normalized the functional scans to MNI space and smoothed them using a Gaussian filter of 6 mm.
To measure functional changes in the motor network, we compared the functional connectivity between regions of interest (ROI-to-ROI analyses) and analyzed the neural complexity by computing the MSE of the BOLD time series in each ROI (defined below) before and after HD-tDCS. As mentioned in the results, we defined four ROIs as the left M1, right M1, left PMd and right PMd, using a meta-analysis from Hardwick and collaborators 38 to define the location of the ROIs in the left hemisphere and then flipped them to the right hemisphere. The exact coordinates of the center of each ROI, which was defined as a sphere with a radius of 10 mm, ROI-to-ROI analyses. We extracted the cross-correlations of the 4 ROIs (left/right M1, left/right PMd) and ordered them for each group and each condition based on the correlation intensity. To compare the organization of the motor network following each stimulation session, we compared these orders using Spearman's rank correlations coefficient.
We also used RM-ANOVAs (one RM-ANOVA for each ROI as a seed region) in the CONN toolbox to compare combinations of pairs between the 4 ROIs (the 6 functional connectivity pairs were: left M1-right M1, left PMd-right PMd, left M1-left PMd, left M1-right PMd, left PMd-right M1, right M1-right PMd) using "Time" (pre-HD-tDCS, post-HD-tDCS) and "Group" (Motor Hotspot, Premotor and Sham) as factors.
Multiscale-entropy. MSE was used to evaluate changes in neural complexity following HD-tDCS in each of the 4 ROIs (left M1, right M1, left PMd, right PMd). The MSE metric quantifies the regularity or complexity of biological signals across a range of coarse-grained timescales or temporal frequencies. Complex systems with 1/f power spectra exhibit constant entropy over various timescales (due to their fractal properties), whereas random noise shows a marked decrease in entropy at longer time scales (as random fluctuations are smoothed out) 35 . We used a pattern matching threshold (r) of 0.5 and a pattern length (m) of 2 35,78,79 . In total, we obtained 20 coarse-sampled scales. Based on the band-pass filter applied, we kept only 9 scales for the analyses (scales 12 to 20: frequency range [1/(TR*scale)]: 0.06-0.1 Hz).

Statistical analyses.
We used R (https://www.r-project.org/about.html) to perform all the statistical analyses, except for the analyses performed using the CONN toolbox (rs-fMRI functional connectivity). We used the CONN toolbox for ROI-to-ROI analysis, and the results are presented at a qFalse discovery rate = 0.05 (corrected for multiple comparisons).
We performed a repeated measures RM-ANOVAs to compare changes in MEP amplitude across groups, and separate RM-ANOVAS to compare changes in MSE for each of the 4 motor ROIs. We used factors of "Time" (pre-HD-tDCS, post-HD-tDCS) and "Group" (Motor Hotspot, Premotor and Sham) in all RM-ANOVAs. To evaluate HD-tDCS-induced changes within each group, we also performed a priori analyses of pre/post-HD-tDCS differences in MEP amplitude for each of the 3 groups using two-way paired t-tests. We report the test statistics and p-values, along with the effect size, using partial eta squared (ƞ p 2 ) for the ANOVA main effects, and Cohen's d for the t-tests. Any post-hoc tests are corrected for multiple comparisons using the Tukey HSD. In addition, given the large inter-individual variability in the data, we also reported the Bayes Factor (BF) for the group level changes over time in MEP amplitude and MSE analyses using the R package "BayesFactor". The BF is used to report a ratio of the strength of the experimental effect compared to the null hypothesis and has recently been used with highly www.nature.com/scientificreports www.nature.com/scientificreports/ variable experimental data, including in brain stimulation studies 80,81 . In our data, we used the anovaBF to determine which model explained the best the variation observed in the results. We tested 3 models: impact of "Time", impact of "Group" and impact of the interaction Group*Time. For each anovaBF, the model with the highest BF value was the be the best combination of models among these. The strength of the impact of this model on the data was determined by the size of the BF value. Briefly, a BF of 0-3 should be considered as weak evidence, a BF of 3-10 as substantial evidence, a BF of 10-30 as strong evidence and a BF > 30 as very strong evidence 82 . We also reported the Jeffrey-Zellner-Siow Prio BF for the t-test analyses to evaluate if in regards to Bayesian statistics the results will favor the alternative, against the null, hypothesis 83 .
Anatomical location of the stimulation site. We used the previously defined ROIs to define the standard M1 and PMd locations (see Fig. 7). The standard M1 ROI matched the location of the hand knob 3 , which is commonly thought to represent the anatomical location of the motor hand area. The standard PMd ROI was within the anatomical PMd boundaries (but due to the large extent of the anatomical PMd the functional ROI only covers a small portion of it).
In all the participants, the stimulation ROI was defined by a cortex-projection of the center electrode of the 4*1 ring. We projected the scalp coordinates of the center of this electrode on the MNI-normalized scans to the cortex. The distance of each voxel of the grey matter from "P" was calculated, enabling the creation of a spherical ROI centered at "P" (10 mm diameter). This resulted in each participant with a cortex-projected stimulated region. We quantified the minimal distance of each individual's cortex-projected stimulated region from the standard M1/PMd ROIs by computing the minimum Hausdorff 's distance (i.e., the minimum distance between the two ROI boundaries) 84 .

Data availability
The unthresholded statistical maps used to generate the motor network (i.e., ICA analyses) are stored on the NEUROVAULT repository (https://neurovault.org/collections/UTZRKARR/). The raw data supporting these analyses will be made available following reasonable requests to the corresponding author. These data are not publicly available due to IRB restrictions.