Introduction

Sensory inputs are necessary for the successful execution and acquisition of skillful voluntary movements1,2,3,4,5,6,7. The primary somatosensory (S1) cortex is a central structure involved in motor learning and motor performance1,8,9,10,11,12. S1 has strong reciprocal connections with the primary motor cortex (M1)1 and one-third of the cortico-spinal tract fibers that control movement originate from S113. Physiological evidence from animal models indicates that a lesion to S1 in monkeys impaired learning of new motor behaviors3. In humans, the extent of contralesional M1 and S1 activity correlated with the severity of post-stroke motor deficit14. Disrupting somatosensation by applying inhibitory 1 Hz repetitive transcranial magnetic stimulation (rTMS) over S1 in healthy individuals prior to skilled motor practice decreased motor skill acquisition12.

Motor learning can be enhanced via the use of non-invasive brain stimulation methods over M1 that modulate neuroplasticity15,16. One non-invasive method that may improve motor performance in healthy individuals and those with neurological disorders is transcranial direct current stimulation (tDCS)16,17,18,19. This is a safe, easy to administer and painless stimulation method that delivers weak direct currents (usually 0.5–2 mA) through surface electrodes placed on the skull. It alters spontaneous brain activity and excitability by the subthreshold modulation of neuronal membranes in a polarity dependent manner20. It is commonly assumed that anodal stimulation (anodal electrode is placed over the region of interest) increases cortical excitability (inducing greater Motor Evoked Potentials (MEP), a commonly used measure of cortical excitability), while cathodal stimulation (cathodal electrode is placed over the region of interest) decreases it15,20; however, this is an overly simplistic assumption, and the effects of tDCS are likely to be much more complex. A non-linear dose–response relationship was found in neurophysiological15,21,22,23,24 and behavioral measures25,26. In order to avoid the assumption that anodal stimulation necessarily reflects excitatory stimulation, instead of the term anodal tDCS, we have described it in this paper as active tDCS with an anodal electrode placed over the region of interest.

In the studies which found active tDCS to be an effective means of improving motor performance in healthy participants, as well as patients suffering from neurological diseases such as stroke and Parkinson’s disease16,17,18,19, the anodal electrode was usually placed over M1. Considering the importance of S1 in motor learning and motor performance, stimulating S1 via tDCS may also be relevant in neurorehabilitation, especially in individuals with stroke who have a large amount of damage in M1. As recovery after stroke is attributed to plastic reorganization in the central nervous system27,28, enhancing the recruitment of relevant regions of interest (among them S1) for motor performance may have important clinical significance. Yet, to the best of our knowledge, only one study—that of Faraji et al.29—investigated the effects of tDCS over S1 on motor performance, and it was conducted in rats. They found that both unilateral and bilateral anodal somatosensory stimulation significantly improved reaching performance. It should be noted that the few tDCS studies in which healthy participants received S1 stimulation focused on sensation measures30,31,32,33,34, and indeed, some of them found improved sensation perception30,31,33. Also, in animals who received S1 stimulation, measures other than motor behavioral measures were mainly reported35,36,37, e.g., the application of tDCS with anode and cathode electrodes over S1 induced polarity-specific bidirectional changes in the N1 component of the sensory-evoked potentials and associated gamma oscillations35.

As far as we are aware, the present study is the first attempt to determine the effects of tDCS over S1 on motor performance in adults (it should be noted that some TMS studies already showed that stimulating S1 enhanced motor learning in healthy participants and individuals with chronic stroke8,11,38). Due to the relevance of S1 to sensation perception, we also assessed tactile sensation and proprioception as secondary outcome measures. The motor task consisted of sequential point-to-point reaching movements on a graphics tablet, a version of a similar, previously used task25,39,40. The sensation tasks consisted of (1) two-point orientation discrimination (TPOD) of fingers, in which the participants were required to indicate the orientation of the prongs41 and (2) proprioception of reaching movements, in which the participants were asked to actively reproduce reaching movements without vision after being passively guided by the examiner (a digitized version of the previously used paper and pencil task—the Brief Kinesthesia Test42). As stimulation of S1 via conventional tDCS, which uses large pad electrodes, delivers current to diverse brain regions rather than to the targeted region of interest only, we used high definition (HD) tDCS with optimized electrode configurations for maximal focal stimulation to S1 or M1. Improved spatial focality of tDCS can be achieved using HD-tDCS43,44,45,46. In comparison to conventional large pad tDCS, HD-tDCS (4 × 1 ring electrode configurations) demonstrated a peak induced electric field magnitude at the sulcus and adjacent gyri directly beneath the active electrode44.

In addition, we took into consideration the effect of sex, and investigated a possible interaction between sex and brain stimulation because it may mediate the effects of tDCS/HD-tDCS on cortical induced electric field current47,48, intracortical excitability49 and different behaviors50,51,52,53,54,55 such as visually guided reaching movements52 and social cognition skills50. In some studies, tDCS significantly affected behavior of women only and in others—in men50,51,54. Only in women, tDCS, in which the anodal electrode was placed over the medial prefrontal cortex, led to slower reaction times (vs. sham) for the social intention attribution task50, and tDCS, in which the anodal electrode was placed over the orbitofrontal cortex, increased the performance of risky decision making54. In men, tDCS, in which the anodal electrode was placed over the dorsolateral prefrontal cortex, enhanced practice-related changes in accuracy of response execution, but in women—prevented practice-related changes51. Concerning motor behavior, Gorbet and Stains52 found that reaction time slowed in men but not in women following inhibitory continuous theta burst stimulation over the contralateral dorsal premotor cortex during a reaching task.

Proposed explanations for sex being a factor that mediates the effect of brain stimulation on cortical excitability and behavior include hormonal levels, neurotransmitter balances and cortical bone structure56. Fluctuations in estradiol and progesterone across the menstrual cycle in women affect cortical excitability56,57,58. It appears that estradiol facilitates cortical excitability, likely through glutamatergic mechanisms, while progesterone metabolites dampen cortical excitability, likely through gamma-aminobutyric acid (GABA)59. Alterations in GABA following active tDCS are related to motor learning capacity60,61. Effects of hormonal fluctuations on manual dexterity are controversial62,63,64. Although Maki et al.63 reported that manual dexterity improved in the mid-luteal phase compared to the follicular phase, Keenan et al.62 reported that dexterity did not change throughout the menstrual cycle. In addition, evidence regarding the amount of induced electrical current at the cortex of men versus women are inconsistent. Whereas Russell et al.47 found that men would receive ~ 45% more current at the cortex than women, perhaps because of the more cancellous bone composition in men compared to women, induced electric field was found to be higher in female head models than male head models across several metrics48.

We hypothesized that HD-tDCS over S1 would be more effective in improving (1) motor performance than sham tDCS and (2) sensation perception than HD-tDCS over M1 and sham tDCS, and that HD-tDCS over M1 would be more effective in improving (1) sensation perception than sham tDCS and (2) motor performance than HD-tDCS over S1 and sham tDCS. We expected to find sex-related differences in the effects of active HD-tDCS over S1 and M1 on motor performance. We could not hypothesize regarding the specific influence of sex on the behavioral response due to the lack of specific relevant tDCS literature in the motor domain, the mixed results regarding the interaction between sex and tDCS with respect to other behavioral domains50,51,52,53,54,55 and the mixed evidence of menstrual cycle phase-related motor performance differences in women62,63,64.

Methods

Study design

This was a single-blind, parallel, randomized, sham-controlled study. Data were collected in a brain and motor behavior laboratory based at Ariel University, Israel. Participants were randomly assigned with a 1:1 ratio, using a random number generator in WINPEPI (by researcher SFT), to one of three groups: (1) HD-tDCS over the S1 (S1 group); (2) HD-tDCS over the M1 (M1 group); and (3) sham HD-tDCS (sham group). All participants were blinded to group allocation. To ensure blinding of participants, the stimulator monitor was hidden from the participants, and the sham stimulation increased and decreased in a ramp-like fashion (see HD-tDCS section). The researcher (YS) who administered the HD-tDCS application and measured the outcomes received allocation information via coded email from another researcher (SFT). Blinding of group allocation was maintained during the data analysis [for a similar approach see25]. The trial was prospectively registered at the ClinicalTrials.gov registry on 06/11/2020 with the trial registration number NCT04618614. The protocol is available on the following website:

https://register.clinicaltrials.gov/prs/app/action/SelectProtocol?sid=S000AC9A&selectaction=Edit&uid=U0005AKF&ts=2&cx=-nam3zd.

Participants

The sample size for this study was determined based on a power analysis calculation that was conducted using G*Power version 3.1.9.7. Power analysis yielded a total sample size of 42 individuals (14 individuals per group) for the detection of a significant interaction with an assumed effect size of 0.25 and a power of 80%. To account for potential data loss, we aimed for a sample size of 15 individuals per group (in total 45 participants). Forty-five participants (24 women, 21 men; aged 24 ± 2 years) participated in the study between November 2020 and January 2021. Participants were included if they were aged between 20 and 35, were right-hand dominant and were healthy according to their report. They were excluded if they took psychiatric medications, had a history of drug abuse or dependence, had any psychiatric or neurological disorder, had a history of seizures, had metal implants in their head or had musculoskeletal deficits interfering with task performance (proper reaching performance with the left arm in sitting) [for a similar approach see25]. Participants signed an informed consent form prior to participating in the study. All the procedures were approved by Ariel University Institutional Review Board (approval number: AU-HEA-SFT-20190326-B) and were performed in accordance with relevant guidelines and regulations. Participants were paid $20 for their participation.

HD-tDCS

The stimulation was administered noninvasively using an M x N 9-channel high-definition transcranial electrical current stimulator from Soterix Medical (New York, NY). Five sintered Ag/AgCl electrodes were attached to plastic holders, filled with conductive gel, and embedded in a HD cap, according to the extended 10–20 method of electrode placing. We administered a single session of 15 min of active stimulation at 1 mA targeting the right S1 (primary somatosensory cortex; postcentral gyrus, based on Talairach labels) by positioning electrodes at the following sites: CP4 (0.97 mA), CP6 (0.03 mA), Cz (− 0.71 mA), CPz (− 0.12 mA), and O2 (− 0.17 mA) or right M1 (primary motor cortex, Brodmann area 4; precentral gyrus, based on Talairach labels) by positioning electrodes at the following sites: C4 (0.82 mA), Fz (− 0.43 mA), F1 (0.18 mA), F6 (− 0.24 mA), and FT8 (− 0.33 mA). HD-Targets brain modelling software (Soterix Medical, New York, NY) was used to determine the tDCS montage for maximal focal stimulation of the right S1 and right M1 (Figs. 1a, 2). In the S1 and M1 groups, the current increased in a ramp-like fashion over the course of the first 30 s, and decreased in a ramp-like fashion over the course of the last 30 s. In the sham group, the position of the electrodes was similar to the position in the S1 group. Once the current reached 1 mA over the first 30 s, it was ramped back down over 30 s. In the last min of the simulation an identical ramp up and ramp down occurred [for a similar approach see65,66,67]. Participants were asked to report any adverse effects and to rank their discomfort from 1 to 10 following one min of stimulation. If the participant felt that he does not exactly know how to define his feeling, he was asked if he means tingling or burning or itching or headache.

Figure 1
figure 1

Course of study. (a) High-definition transcranial direct current stimulation (HD-tDCS) montage for maximal focal stimulation of the (1) right primary somatosensory cortex (S1) and (2) right primary motor cortex (M1) using the HD-Targets modelling software (Soterix Medical, New York, NY). The location and current intensity value of each stimulating electrode are shown. Red denotes anodal stimulation while blue denotes cathodal stimulation. (b) Experimental stimuli. (1) TPOD = Two-point orientation discrimination. The prongs were oriented randomly either across or down (trial-by-trial) to the left index and little fingers. (2) General setup of the proprioception task. The non-dominant left arm of the participant that held the stylus was placed on a custom-made device. The participants were required to reproduce a reaching movement from the starting position to target 2, with their eyes closed, after being passively moved by the device. This was repeated 10 times. (3) General setup of the motor task. Note that in this figure, for clarity, the targets and numbers are shown 3 times their relative size compared to those shown in the experiment.

Figure 2
figure 2

Current flow modeling during 1 mA High-definition transcranial direct current stimulation (HD-tDCS) using the HD-Target software (Soterix Medical, New York, NY). Current-flow models of (a) the right primary somatosensory cortex (S1) and (b) the right primary motor cortex (M1) are shown on 2D and 3D reconstructions of the cortical surface. Skin, skull, and cerebrospinal fluid (CSF) masks are suppressed to reveal the underlying gray matter mask. A head model derived from the MNI 152 dataset was used.

Motor sequence learning task

In all participants, the non-dominant left arm was tested. The non-dominant arm was tested to challenge the motor performance of healthy participants and to allow more room for improvement in the motor performance (aiming to avoid ceiling effects). This approach is acceptable in motor learning studies [e.g.,68,69] and tDCS studies which investigate upper limb motor performance/learning16,70. After placing the tDCS cap on the head, the participants performed a sequential point-to-point movement task on the graphics tablet, a version of a similar, previously used task39,40. The stimuli consisted of a starting point and five targets around a semicircle (Fig. 1b1). The participants were instructed that the targets would change their colors following the sequence: 4-1-3-2-5, and to perform the task as fast and accurately as possible. A detailed description of the task is provided in a previous study25.

Initially, the participants were required to perform 3 sequences without errors to familiarize themselves with the task. Then, they performed the pretest which consisted of one block of 6 sequences. Two min after starting the appropriate stimulation, simultaneously with the stimulation, they performed 3 blocks of 6 sequences, i.e., 18 sequences. Brief breaks (30 s) were given between each block, while within a block breaks were not given. After finishing the tDCS stimulation, the participants performed a posttest, which was identical to the pretest (i.e., one block of 6 sequences). The stimuli were visible throughout each trial, the duration of each movement was subject-dependent (it started after the stylus was placed at the start position for 500 ms and ended when the participants reached the target and remained there for 500 ms).

Three outcome measures for the motor task were used. Movement time (s) of the reaching movements was defined as the time from movement onset (first time the tangential velocity was greater than 5% of the peak tangential velocity) until the end of the movement (the last time the tangential velocity was greater than 5% of the peak tangential velocity). Reaction time (s) was defined as the time between when the target appeared in green, and movement onset (as defined above). Endpoint error (cm) was defined as the straight-line distance from the stylus location at the end of the movement to the center of the target (cm). Improved motor performance was indicated by a shorter movement time, a shorter reaction time and a smaller endpoint error.

Sensation tasks

TPOD: This task involved presentation of a two-pronged instrument (spacing 5, 4, 3 and 2 mm) to the palmar side of the distal pads of the left index and little fingers (Fig. 1b2). These prongs were oriented randomly either across or down (trial-by-trial), with respect to the proximal–distal finger axis, using descending spacing order (see41 for advantages of this task over the traditional two-point discrimination task). Initially, the participants were required to indicate each orientation of the prongs to the palmar pad of the middle finger with their eyes open and closed to familiarize themselves with the task. Then, participants verbally indicated the orientation at each trial with their eyes closed. Proportion of correct trials were averaged as a measure of tactile spatial acuity41.

Proprioception: The participants performed a point-to-point movement task on the same graphics tablet that was used for the motor task. The task is a version of a similar, previously used paper and pencil task—the Brief Kinesthesia Test42, in which the participant was asked to actively reproduce reaching movements without vision after being passively guided by the examiner. The passive movement of the participant’s arm was performed here using a custom-made device, which was connected to the tablet (Fig. 1b3). The linear actuator of the device was designed as a toothed-belt-driven carriage, sliding on two circular steel rods of diameter 10 mm. It was actuated by a NEMA-17 stepper motor, controlled by an Arduino Uno controller and a DRV8825 driver. The resolution of the movement is 50 pulses to mm, i.e., it provides a high accuracy of movement. At the beginning of every actuation, a homing was done by a microswitch located at the starting position. As a preliminary step for using the linear actuator, we tested it to make sure it could withstand the required loads. It can carry more than 20 kg without slipping or velocity reduction.

Initially, to familiarize the participants with the task, they were required to hold the stylus, without the device, in their left hand and reproduce a reaching movement from the starting position to target 3 after being passively moved by the examiner’s hand, with their eyes open and eyes closed. Then, with eyes closed and blindfolded, the assessment with the device began. The non-dominant left arm of the participant that held the stylus was placed on the device. After holding the stylus at the starting point, the device passively moved the arm, from the starting point to target 2, 17 cm distant from the starting point, over two sec. After remaining at target 2 for 500 ms, the tablet’s screen changed its color from black to white, and the examiner lifted the participant’s arm from the device and returned it to the starting position to allow the participant to actively reproduce the reaching movement, without the device. The participants were instructed to reach toward the target and stay there. After 500 ms, the tablet’s screen changed its color from black to white, and the examiner lifted the participant’s arm and returned it to the starting position for the next passive reaching movement. The pretest and the posttest consisted of 10 passive and 10 active reaching movements. Two outcome measures for the proprioception were used: movement time (s) of the reaching movements and endpoint error (cm) (see the section of motor sequence learning task).

Trials that were not properly recorded by the tablet for the motor and proprioception tasks were discarded.

Statistical analysis

Age and sex were compared between groups (S1, M1, sham) using Kruskal–Wallis (as age was not normally distributed) and chi-squared tests, respectively. Differences between groups with respect to each of the outcome measures in the pretest were investigated using one-way analysis of variance (ANOVA) with Bonferroni correction for multiple comparisons. The effects of stimulation, sex and time on the outcome measures were investigated using a mixed design ANOVA with time (pretest, posttest) as the within-subject factor, and group (S1, M1, sham) and sex (male, female) as the between-subject factors with Bonferroni correction for multiple comparisons. Since reaction time values were not normally distributed in each group, they were log-transformed before this analysis (the original values are presented for clarity). The differences between groups with respect to the frequency of adverse effects were investigated using a chi-squared test. The differences between groups with respect to the discomfort from adverse effects were investigated using Kruskal–Wallis with Bonferroni correction for multiple comparisons [for a similar approach see43]. In addition, due to an imbalanced number of males and females within the groups, and low number of participants in each group, a mixed model ANCOVA was conducted with time as the within-subject factor, group as the between-subject factors and sex as covariate with Bonferroni correction for multiple comparisons. All tests were done using SPSS (version 26.0) with initial significance levels of p < 0.05.

Results

The flowchart illustrating the process of the study is shown in Fig. 3. Fifty-four participants underwent the pre-enrollment screening evaluation. Of those, nine did not meet inclusion criteria. Age (S1 group: 23.5 ± 2.6 years; M1 group: 23.9 ± 2.6 years; sham group: 23.3 ± 1.8 years) and sex (S1 group: six women; M1 group: eight women; sham group: seven women) did not differ between groups (p > 0.717, for all).

Figure 3
figure 3

Trial flowchart. HD-tDCS S1/M1 = High-definition transcranial direct current stimulation of primary somatosensory cortex (S1)/primary motor cortex (M1).

Motor sequence learning task

Mean values of movement time (s), reaction time (s) and endpoint error (cm) by group and time are shown in Table 1. All these outcome measures did not show significant differences between groups in the pretest (p > 0.508, for all).

Table 1 Means, standard deviations and confidence intervals of reaction time, movement time and endpoint error of the motor task for stimulation groups in time points.

Effects on reaction time (s)

A main effect of Sex (F(1,39) = 5.206; p = 0.028; partial η2 = 0.12; observed power = 0.61) showed that, overall, reaction time was faster in women (0.38 ± 0.27 s) compared to men (0.46 ± 0.12 s). Group x Time interaction varied across Sex as shown by a second-order interaction between Sex x Group x Time (F(2,39) = 4.419; p = 0.019; partial η2 = 0.19; observed power = 0.73). Only in men, the interaction Group x Time was significant (F(2,18) = 3.914; p = 0.039; partial η2 = 0.30; observed power = 0.63) such that only in M1 group, reaction time decreased significantly in posttest (0.44 ± 0.11 s) compared to pretest (0.47 ± 0.13 s) (F(1,7) = 15.109; p = 0.006; partial η2 = 0.68; observed power = 0.91) (Fig. 4). No other significant effects were observed.

Figure 4
figure 4

Mean reaction time (s) of reaching movements in pretest and posttest in the three groups in men (a) versus women (b). Error bars show standard deviation. Asterisks denote a significant difference.

Effects on movement time (s)

A main effect of Time (F(1,39) = 13.550; p = 0.001; partial η2 = 0.26; observed power = 0.95) showed that, overall, movement time decreased significantly in posttest (0.77 ± 0.22 s) compared to the pretest (0.83 ± 0.28 s). This effect was, however, modulated by Sex, as was shown by the interaction of Sex x Time (F(1,39) = 12.460, p = 0.001; partial η2 = 0.24; observed power = 0.93). Only in women, movement time decreased significantly in posttest (0.80 ± 0.21 s) compared to the pretest (0.92 ± 0.27 s), across groups (F(1,21) = 21.210, p < 0.001; partial η2 = 0.51; observed power = 0.99) (Fig. 5). No other significant effects were observed.

Figure 5
figure 5

Mean movement time (s) of reaching movements in women and men in pretest and posttest, across stimulation groups (collapse across stimulation groups is presented because the interaction Time × Sex × Group was not significant). Error bars show standard deviation. Asterisks denote a significant difference.

Effects on endpoint error (cm)

A main effect of Sex (F(1,39) = 14.005; p = 0.001; partial η2 = 0.26; observed power = 0.95) showed that, overall, endpoint error was smaller in women (0.33 ± 0.10 cm) compared to men (0.44 ± 0.13 cm). A main effect of Time (F(1,39) = 4.483; p = 0.041; partial η2 = 0.10; observed power = 0.54) showed that, overall, endpoint error was smaller in posttest (0.37 ± 0.10 cm) compared to pretest (0.40 ± 0.15 cm). This effect was, however, modulated by Sex, as was shown by the interaction of Sex x Time (F(1,39) = 6.335; p = 0.016; partial η2 = 0.14; observed power = 0.69). Only in men, the endpoint error was significantly smaller in posttest (0.40 ± 0.09 cm) compared to pretest (0.48 ± 0.15 cm) (F(1,18) = 9.938; p = 0.006; partial η2 = 9.94; observed power = 0.85) (Fig. 6). No other significant effects were observed.

Figure 6
figure 6

Mean endpoint error (cm) of reaching movements in women and men in pretest and posttest, across stimulation groups (collapse across stimulation groups is presented because the interaction Time × Sex × Group was not significant). Error bars show standard deviation. Asterisks denote a significant difference.

Sensation tasks

Mean values of proportion of correct trials in the TPOD task, and movement time (s) and endpoint error (cm) in the proprioception task by group and time are shown in Table 2. All these outcome measures did not differ between groups in the pretest (p > 0.360, for all).

Table 2 Means, standard deviations and confidence intervals of TPOD, movement time and endpoint error of the sensation tasks for stimulation groups in time points.

Two-point orientation discrimination

A main effect of Time (F(1,39) = 8.149; p = 0.007; partial η2 = 0.17; observed power = 0.80) showed that, overall, the percent of correct trials was higher in posttest (72.28 ± 14.63%) compared to pretest (67.23 ± 13.50%). The interaction of Group x Time reached border-line significance (F(2,39) = 2.765; p = 0.075; partial η2 = 0.12; observed power = 0.51). Our interest was focused on clarifying whether the proportion of correct trials differed between groups at timepoints and whether the proportion of correct trials differed between time points within each group. Therefore, despite the borderline significance, the interaction was further investigated. Only in the S1 group, the proportion of correct trials increased significantly in posttest (81.25 ± 11.81%) compared to pretest (70.31 ± 12.90%) (F(1,13) = 8.493; p = 0.012; partial η2 = 0.40; observed power = 0.77) (Fig. 7). The effect size of the change in the correct trials from pretest to posttest was high for S1 (D = 0.7) whereas it was low for M1 group (D = 0.2) and sham group (D = 0.3). In addition, only in posttest, the proportion of correct trials differed between groups (F(2,42) = 5.339; p = 0.009) such that it was significantly higher in the S1 group compared to the M1 group (65.83 ± 16.34%; pBonferroni = 0.009), and just tended to be higher compared to the sham group (70.00 ± 11.38%; p = 0.026; pBonferroni = 0.079). No other significant effects were observed.

Figure 7
figure 7

Mean percentage of correct answers in the two-point orientation discrimination (TPOD) test in each group at pretest and posttest. Error bars show standard deviation. Asterisks denote a significant difference.

Proprioception

Effects on movement time (s)

A main effect of Time (F(1,39) = 26.688; p < 0.001; partial η2 = 0.41; observed power = 0.99) showed that, overall, movement time was shorter in posttest (1.97 ± 0.63 s) compared to pretest (2.24 ± 0.70 s). No other significant effects were observed.

Effects on endpoint error (s)

No significant effects were found.

Adverse effects

The stimulation was well tolerated by the participants, and no sessions were aborted due to adverse effects. The occurrence of adverse effects in the S1, M1 and sham groups are displayed in Table 3. Frequency of adverse effects and discomfort from the adverse did not differ between the groups.

Table 3 Frequency and discomfort of adverse effects.

The mixed design ANOVA effects for motor and sensation tasks are displayed in detail in Supplementary Table 1. The ANCOVA results for motor and sensation tasks are displayed in detail in Supplementary Table 2.

Discussion

To the best of our knowledge, this is the first study that evaluated the effects of tDCS over S1 on motor performance in adults. In addition, we considered the sex-mediated effects of the stimulation on the response. By means of HD-tDCS applied with optimized electrode configurations for maximal focal active stimulation, a considerably nuanced determination of stimulation site (S1, M1, and sham stimulation) effects on motor performance was achieved. We also investigated the effects on sensation perception, as a secondary outcome measure. We found that reaction time, movement time, and endpoint error of a sequence of reaching movements did not differ between S1, M1, and sham groups. Nonetheless, reaction time improved from pretest to posttest in the men from the M1 group only. Reaching movement time improved from pretest to posttest in women only, in a similar manner in all stimulation groups, whereas endpoint error improved in men only. Regarding sensation perception, percent of correct trials for the TPOD was significantly higher in the S1 group compared to the M1 group in the posttest (but not in the pretest), and improved from pretest to posttest in the S1 group only. The reaching movement time of the proprioception task improved from pretest to posttest, in a similar manner in all groups, and endpoint error did not change.

Our finding that HD-tDCS over S1 did not improve motor performance as compared to sham stimulation in healthy participants is not in line with our first hypothesis. We had assumed positive effects in motor performance for the S1 stimulation based on the neuroanatomical1 and behavioral links8,11,29,38 between the sensory and motor systems. Improved motor performance was found following tDCS over sensory areas in rats29, and following TMS over S1 in both healthy participants38 and individuals with stroke8,11. Exciting ipsilesional S1 using 5 Hz rTMS with skilled motor practice enhanced motor learning in individuals with chronic stroke8. Meehan et al.11 found that continuous theta burst stimulation over contralesional S1 or M1 in individuals with chronic stroke enhanced performance of serial targeting tasks performed by the hemiparetic upper limb. They even found that stimulation over contralesional S1 led to significant changes in the functional ability of the hemiparetic upper limb, indexed by time to complete the Wolf Motor Function Test, as compared to stimulation over the contralesional M1, or control stimulation. The different mechanisms underlying the effects of tDCS versus TMS may partly explain the more potent effects of TMS than those of tDCS. In contrast to TMS, which directly induces action potentials71, tDCS alters spontaneous brain activity and excitability by modulating the neuronal membranes in a polarity dependent manner20. It is also possible that the stimulation protocol of the S1 group in the current study is not optimal for enhancing tDCS effects on motor performance. Variable and even inverted neurophysiological15,21,22,23,24 and behavioral effects25,26 were found in response to different tDCS protocols, depending on several parameters such as stimulation intensity and duration. For example, conventional anodal stimulation at 2 mA led to faster reaction times than cathodal stimulation, but only before 13 min of stimulation had elapsed, whereas after 13 min, the reaction times under cathodal stimulation were faster72. However, the optimal protocol is yet to be determined. On the other hand, it should be noted that, as with our results, in some TMS studies, stimulating S1 in healthy participants73 and those with stroke74 did not improve motor performance either. The mixed results in the non-invasive neuromodulation studies may be explained by the different protocols (such as rTMS vs. theta burst TMS), the outcome measures used, and variable responses in humans to non-invasive brain stimulation techniques, e.g., due to differences in the individual morphology of the cerebrospinal fluid and brain21,75.

Whereas HD-tDCS over S1 in the current experimental setup did not improve motor performance, it did improve the percent of correct trials for the TPOD from pretest to posttest; moreover, the percent of correct trials was significantly higher in the S1 group compared to the M1 group in the posttest (but not in the pretest). The percent of correct trials in the S1 group was not significantly higher compared to the sham group (p = 0.026; pBonferroni = 0.079). However, the effect sizes of the change between pretest and posttest were high in the S1 group (D = 0.7) but low in the sham and M1 groups (D = 0.3 and D = 0.2, respectively). These results are reasonably in line with our second hypothesis. Given S1’s essential role in somatosensory processing, it is reasonable to assume that stimulating S1 would produce more effects than the sham and M1 stimulations on sensation perception. The modulatory excitatory and inhibitory effects that are exerted by M1 on S1 via anatomical projections1,76,77 may have also contributed to sensorimotor integration and to the augmented effects in sensation perception specifically after S1 stimulation because simultaneously to the S1 stimulation, the participants performed 18 sequences of the motor task. Our results are also consistent with previous findings in which healthy participants30,33 and those with stroke31 improved their performance of the tactile spatial discrimination task following active tDCS with an anodal electrode placed over the S1. However, it should be noted that, as mentioned in the methods and results section, a mixed model ANCOVA was also conducted (with sex as covariate) due to an imbalanced number of males and females within the groups, and low number of participants in each group. In the ANCOVA results (Supplementary Table 2), the percent of correct trials only tended to be higher in the S1 group compared to the M1 group in the posttest (p = 0.068). Therefore, further studies with larger sample sizes are required for validation of the current results.

Even though tDCS can modulate the proprioceptive afferent system by changing the excitability of projections to propriospinal neurons78,79, HD-tDCS over S1 did not improve proprioception in the current study. Muffel et al.32 found that active tDCS with an anodal electrode placed over S1 elicited opposing effects on proprioceptive accuracy as a function of age. While young adults (27.0 ± 2.4 years) showed a slight improvement, older adults (69.4 ± 4.9 years) exhibited a decline in performance during a-tDCS. This study, similar to our own, used 1 mA tDCS for 15 min; however, in contrast to our study, Muffel et al.32 used a different experimental design (crossover), other assessments of somatosensory perception (an arm position matching task in a robotic environment), and a different tDCS montage (conventional). Taken together, these differences could explain the contradictory results. It is also important to note that the current simulation used in the Muffel et al. study32 showed a non-focal distribution of the induced electric field throughout the sensorimotor system. In light of the interhemispheric interactions between S1 and the secondary somatosensory cortex (S2) in humans80,81, and the involvement of S2 in sensation perception82, the larger spread of the current in the Muffel et al. study32 could have led to the improved proprioception in young adults. Indeed, sensation perception improved following tDCS over S2 compared to sham stimulation in individuals with stroke31. Our finding that reaching movement time for the proprioception task improved from pretest to posttest similarly in all groups can be related to the initial fast learning phase during the session83. Interestingly, such an initial fast learning phase did not occur in any group for the endpoint error of the reaching movement in the proprioception task. This may reflect the complexity of performing a reaching movement toward the center of the target with eyes closed. It is possible that additional repetitions would have improved the endpoint error.

Sex modulated the effects of active HD-tDCS over M1 on reaction time, which improved from pretest to posttest only in the men in the M1 group. Similarly, Gorbet and Stains52 found that reaction time (but not movement time) slowed in men but not in women following inhibitory continuous theta burst stimulation (using TMS) over the contralateral dorsal premotor cortex during visually guided reaching movements, a task similar to the one used in the current study. In the Gorbet and Stains study52 and in ours, the specific influence of non-invasive brain stimulation on reaction time in men may be related to previous findings, which pointed to more bilateral patterns of activity in women and unilateral activity in men in different brain regions during motor tasks84 and others, such as mental rotation85. If men rely more on the contralateral hemisphere, they might be more responsive than women to stimulation in that hemisphere. The finding of sex differences in reaction time but not movement time in the Gorbet and Stains study52 led them to suggest that sex differences in visually guided reaching movements are probably more strongly associated with movement planning than with motor execution. This suggestion may explain our finding that only reaction time and not movement time improved following HD-tDCS in men.

Additional reasons for sex-related differences in the responses to tDCS may be related to hormonal differences and bone density. Hormonal levels fluctuate significantly more in women than in men56,86. As it appears that progesterone drives the increase of cortical inhibition and estradiol enhances excitability56,57,58, it is possible that the phase of menstrual cycle of the women that participated in the current study affected their responsiveness to the tDCS. However, we did not control for the menstrual phase. With regard to bone density, a computational study indicated that the induced electric field is higher in female head models, on average, than male head models across several metrics48. By contrast, in other tDCS studies, men were found to receive more current than women47,87,88. In addition, men were found to have a more cancellous parietal bone and females a denser parietal bone. Indeed, studies have indicated sex-related anatomical differences in head structures89,90,91. Our finding supports the evidence that sex moderates the effects of tDCS51,52,53,54,55. It seems that the high rates of inter- and intra-individual variability with regard to the effects of tDCS on motor performance92 may also be related to sex56. tDCS studies may find more meaningful results if they are analyzed according to sex.

The improvement in perception discrimination but not motor performance following HD-tDCS over S1, and the improvement in reaction time but not sensation perception in the men from the M1 group may point to a selective HD-tDCS influence that is dependent on the stimulation site. Despite reciprocal projections between S1 and M1, it seems that this input is relatively weak with regards to producing behavioral effects. Indeed, Kinnischtzkeet al.77 found that M1 provides input to nearly all S1 pyramidal neurons in the mouse, yet most of that input on its own was unlikely to make S1 neurons fire. It seems that using HD-tDCS44,46,53 together with the brain modelling software to determine the tDCS montage for maximal focal stimulation of the right S1 and M1 improved the spatial focality of the current and induced stimulation site effects. On the other hand, since the measures of the motor reaching sequence task did not differ between groups, and the TPOD that differed between the S1 and M1 group only tended to do so between the S1 and sham group, it should be acknowledged that this interpretation relates to moderate effects that were mainly found within groups and not between groups.

We also found sex-mediated effects that are not related specifically to the stimulation type. For the motor task, reaching movement time improved from pretest to posttest in a similar manner in all stimulation groups, only in women, whereas reaching endpoint error improvement was found only in men. This finding is in line with previous evidence regarding differences in motor ability between the sexes93,94,95,96,97. Similar to the specific improvement found here in reaching reaction time in women, females also perform the 9-hole peg test, a measure of finger dexterity, faster than males94,95. There is a biological basis for sex differences98, which also encompasses the motor domain, for example, female rats have 20% fewer dopaminergic neurons which are involved in motor control99. The similar improvement in movement time and endpoint error from pretest to posttest across groups, including sham stimulation, in women and men, respectively, may reflect motor learning over time and a placebo effect. The motor learning over time can be related to the initial fast-learning phase within session83. The possible placebo effect could be related to the finding that the frequency of adverse effects and discomfort from the adverse did not differ between the groups.

The study has several limitations. First, despite improved spatial focality of current using HD-tDCS and brain modelling software to determine the tDCS montage for maximal focal stimulation of S1 and M1, one cannot assume that stimulation was exclusively delivered to S1 or M1, and the electric fields (V/m) could have differed between groups because of differences in the individual participants’ anatomical features100. Responses to the HD-tDCS may also be related to variable anatomical features. Second, as the current study design was a single-blind randomized controlled study, the experimenter was not blinded to group allocation. It should be noted, though, that besides the TPOD task which was based on the experimenters’ scoring (based on a structured protocol), the scoring of the motor and proprioception task was automatically computed by the MATLAB software. Third, because changes in the levels of estradiol and progesterone during the follicular and luteal phases affect cortical excitation/inhibition and behavior50,51,52,53,54,55,56,57,58, it is possible that not recruiting women in the same phase of the menstrual cycle increased the variability of the tDCS related effects on motor performance and sensation perception. For example, it is possible that women who are tested in the second half of the follicular phase which is characterized by high level of cortical excitation86 may respond better to the tDCS. Fourth, we did not take into consideration the smoking habits of the participants. This potentially confounding factor could have affected the results because nicotine may affect cortical excitability101 and tDCS effects on MEPs102. Fifth, the montage of the sham group was similar to the montage of the S1 group only. It may have been useful to add a sham group with a montage that is similar to the M1 group or alternatively to randomize the montage of the sham group to be similar to either the S1 group or M1 group. It should be noted, though, that the frequency of adverse effects and discomfort from the adverse did not differ between the groups. Lastly, despite sample size calculation, the value of power and effect size (partial η2) was smaller than 0.8 and 0.14, respectively (these values are equivalent to acceptable power and high effect size), for some of the results (Supplementary Table 1), and number of males and females was not equal in each of the groups. Therefore, the findings reported here should be reproduced in larger cohorts.

Future tDCS studies should consider menstrual phases of women to control for the hormonal fluctuations and their effect on cortical excitation and behavioral effects56,86. Men lack the cyclic fluctuation of sex hormones56; however, the relative contribution of testosterone may also be taken into consideration. Hormonal levels and bone composition of each individual should be measured. As sex differences are founded in biological variants inherent in men and women103, it is important to take into account the participants’ biological sex when analyzing effects of tDCS.

Conclusions

A montage of 15-min active HD-tDCS at 1 mA with focal stimulation of the primary sensory cortex did not affect motor performance but did improve tactile discrimination perception in young healthy participants. Such stimulation over M1 improved reaction time in men only but did not affect movement time and endpoint error of reaching movements or sensation perception. These findings suggest that despite the links between the sensory and motor systems, the effect of excitatory stimulation of S1 and M1, at least in the current montages, may be more specific for influencing sensation perception and motor performance, respectively, i.e., stimulation focality matters. In addition, the results demonstrate sex-mediating effects of HD-tDCS on motor performance. The influence of stimulation site and sex should be taken into consideration in clinical tDCS studies that aim to improve upper limb sensorimotor functioning in individuals with stroke.