Next Article in Journal
Simple Prediction of Metabolic Equivalents of Daily Activities Using Heart Rate Monitor without Calibration of Individuals
Next Article in Special Issue
Influence of a Water-Based Exercise Program in the Rate of Spontaneous Birth: A Randomized Clinical Trial
Previous Article in Journal
Psychological Safety in Ghana: Empirical Analyses of Antecedents and Consequences
Previous Article in Special Issue
Investigation the EMG Activities of Lower Limb Muscles When Doing Squatting Exercise in Water and on Land
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 17
Issue 1
10.3390/ijerph17010215
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Priming Effects of Water Immersion on Paired Associative Stimulation-Induced Neural Plasticity in the Primary Motor Cortex

by
Daisuke Sato
1,2,*,
Koya Yamashiro
1,2,
Yudai Yamazaki
2,3,
Koyuki Ikarashi
2,3,
Hideaki Onishi
2,
Yasuhiro Baba
1 and
Atsuo Maruyama
4
1
Department of Health and Sports, Niigata University of Health and Welfare, Shimamicho 1398, Kita-ku, Niigata City, Niigata 950-3198, Japan
2
Institute for Human Movement and Medical Science, Niigata University of Health and Welfare, Shimamicho 1398, Kita-ku, Niigata City, Niigata 950-3198, Japan
3
Graduate School, Niigata University of Health and Welfare, Shimamicho 1398, Kita-ku, Niigata City, Niigata 950-3198, Japan
4
Department of Rehabilitation Medicine, Kagoshima University, Sakuragaoka 8-35-1, Kagoshima City, Kagoshima 890-8520, Japan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(1), 215; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17010215
Submission received: 17 November 2019 / Revised: 21 December 2019 / Accepted: 23 December 2019 / Published: 27 December 2019
(This article belongs to the Special Issue Immersing in an Enriched Environment: Aquatic Therapy in a Modern Perspective)
Browse Figures
Review Reports Versions Notes

Abstract

:
We aimed to verify whether indirect-wave (I-wave) recruitment and cortical inhibition can regulate or predict the plastic response to paired associative stimulation with an inter-stimulus interval of 25 ms (PAS25), and also whether water immersion (WI) can facilitate the subsequent PAS25-induced plasticity. To address the first question, we applied transcranial magnetic stimulation (TMS) to the M1 hand area, while alternating the direction of the induced current between posterior-to-anterior and anterior-to-posterior to activate two independent synaptic inputs to the corticospinal neurons. Moreover, we used a paired stimulation paradigm to evaluate the short-latency afferent inhibition (SAI) and short-interval intracortical inhibition (SICI). To address the second question, we examined the motor evoked potential (MEP) amplitudes before and after PAS25, with and without WI, and used the SAI, SICI, and MEP recruitment curves to determine the mechanism underlying priming by WI on PAS25. We demonstrated that SAI, with an inter-stimulus interval of 25 ms, might serve as a predictor of the response to PAS25, whereas I-wave recruitment evaluated by the MEP latency difference was not predictive of the PAS25 response, and found that 15 min WI prior to PAS25 facilitated long-term potentiation (LTP)-like plasticity due to a homeostatic increase in cholinergic activity.

1. Introduction

The use-dependent changes in synaptic strength (plasticity) are essential for learning and memory [1]. Among the multiple types of synaptic plasticity described thus far [2,3], Hebbian long-term potentiation (LTP) and depression (LTD) are the most thoroughly studied as potential neuroplastic mechanisms for learning. LTP and LTD are mutually dependent as induction thresholds, and are influenced by the prior history of neuronal activity and functional state of the synapse, termed as metaplasticity [4]. Metaplasticity represents a neuroprotective mechanism that stabilizes synaptic weights in neuronal networks while maintaining the capacity for synaptic plasticity by varying the induction thresholds according to the integrated postsynaptic activity [5].
Various non-invasive transcranial magnetic stimulation (TMS) protocols can be used to explore the neurophysiological mechanisms underlying synaptic plasticity in the human cortex. In particular, paired associative stimulation (PAS) [6] involves the application of an electrical stimulus to the median nerve at the wrist, followed by TMS over the primary motor cortex (M1) after 25 ms (PAS25). Repeated pairings at this interval increase the corticospinal excitability for 30–60 min, which is measured as motor-evoked potentials (MEPs). However, recent studies show that PAS25-induced plasticity is not observed in 25%−64% of healthy young volunteers [7,8,9,10,11]. Numerous factors have been proposed to explain this non-responsive population, such as in the review article by Ridding and Ziemann [12]. As non-responsiveness affects both research design and the extent of potential therapeutic applications, we investigated the factors that may predict or alter the plastic response to PAS25.
A recent study reported that the individual response variability to theta burst stimulation (TBS) was strongly influenced by the recruited interneuron networks [13]. Thus, we first questioned whether an analogous recruitment process accounts for individual variability in the PAS25 response; in particular, we examined whether the variation in potentiating indirect-wave (I-wave) recruitment influences the PAS25-induced plasticity. Furthermore, γ-aminobutyric acid (GABA)ergic interneurons have been reported to modulate the plastic response to PAS25 [14,15]. In fact, cortical GABAergic activity negatively interacts with short-latency afferent inhibition (SAI) [16,17]. As the PAS25 protocol is actually a repeated SAI protocol, it is likely that PAS25-induced potentiation may be associated with a progressively decreasing inhibition by the conditioning afferent stimuli on later I-wave recruitment [15]. Thus, we predicted that individual SAI with an inter-stimulus interval (ISI) of 25 ms may also be associated with the plastic response to PAS25, because both SAI with an ISI of 25 ms and PAS25 recruit the same sensorimotor networks.
The history of prior cortical activity may be another possible factor influencing the PAS25 response [8,18,19]. Water immersion (WI) is one simple way of modulating the cortical activity in the sensorimotor system [20,21]. Our previous study showed that WI decreased SAI [21] through a form of cortical inhibition that originates from the central cholinergic modulation of inhibitory circuits different from those underlying short-interval intracortical inhibition (SICI) [22]. Thus, WI may suppress the specific inhibitory circuits in the M1 that are modulated by cholinergic activity. However, effects of WI on cholinergic activity have remained unclear. Considering the characteristics of cholinergic activity, which contribute to the homeostatic modulation of sleep and wake states [23], cholinergic activity may invert after WI due to a homeostatic response. In that case, WI would promote PAS25 because the effects of PAS25 could be strengthened by cholinergic activity [24].
Here, we tested whether WI primed PAS25 by altering the cholinergic modulation of inhibitory circuits and M1 excitability. The elucidation of the effects of WI priming on the neural plasticity of M1 could increase the reliability of PAS25 induction for research on human cortical sensorimotor integration and synaptic plasticity, and could facilitate the development of improved aquatic therapies for neurological patients.

2. Materials and Methods

2.1. Subjects

A total of 27 right-handed subjects (defined by the Edinburgh handedness questionnaire; Oldfield, 1971) participated in this study. We enrolled 18 participants for experiments 1 and 2 and nine for experiment 3. Informed consent was obtained from all subjects. The present study was conducted in accordance with the Declaration of Helsinki and approved by the ethics committee of Niigata University of Health and Welfare (approval number 17656-160404).

2.2. TMS

TMS was performed by using two Magstim 200 stimulators (Magstim, Dyfed, United Kingdom) connected via a Y-cable to a figure-of-eight-shaped coil, with an internal wing diameter of 7 cm. The coil was held with the handle pointing backward and 45° laterally to the interhemispheric line to induce an anteriorly directed current in the brain, and was optimally positioned to obtain MEPs in the abductor pollicis brevis (APB) muscle. We verified the target position of the coil relative to the brain anatomy by using a frameless TMS navigation system (Brainsight, Rogue Resolution, United Kingdom). The position of the coil was fixed using this system throughout all the experiments.
The resting motor threshold (RMT) was defined as the minimum stimulation intensity over the motor hotspot required to evoke an MEP of ≥50 μV in 5 out of 10 trials [25]. The active motor threshold (AMT) was defined as the lowest intensity required to evoke an MEP of 200 μV in more than 5 of 10 consecutive trials, whereas the subjects maintained approximately 10% contraction of the target muscle. To maintain 10% contraction, as performed by previous studies [26], a rectified running average EMG with an averaging window of 175 ms was used to provide visual feedback on the monitor to the participants.

2.3. Electromyographic Recording

Surface electromyographic (EMG) recordings were acquired in a belly-to-tendon montage from the APB and the first dorsal interosseous (FDI) muscle of the right hand. The raw signal was amplified and band-pass filtered between 5 Hz and 2 kHz (AB-601G, Nihon Kohden, Tokyo, Japan), and transferred through a Micro 1401 Laboratory Interface (Cambridge Electronic Design, Cambridge, United Kingdom) to a personal computer for further offline analysis. All electrodes were covered with a transparent film (Tegaderm Hydrocolloid Dressing, 3M Japan, Tokyo, Japan) for waterproofing, as described previously [27]. The waterproofing was applied before starting the experiment and removed after the final assessment in all experiments.

2.4. PAS Session

The PAS session comprised the application 180 electrical stimuli to the right median nerve at the wrist, paired with a single posterior-to-anterior (PA)-directed TMS over the hotspot of the left APB muscle. Square-wave pulses (duration, 200 µs) at three times the sensory threshold were applied through a bipolar electrode (cathodal proximal). TMS was delivered through a figure-of-eight-shaped coil connected to a Magstim 200 magnetic stimulator that was held in the same position, by using the brain navigation system described above. Stimulus pairs were delivered at 0.2 Hz.
Stimulation was applied at an intensity adjusted to evoke an MEP in the relaxed APB muscle of approximately 1 mV peak-to-peak at baseline (TS1mV_base). The effects of PAS with an ISI of 25 ms between the peripheral and TMS stimulus (PAS25) on the MEP amplitude were tested. This paradigm has been shown to induce a long-lasting MEP increase [6,28]. Subjects were instructed to look at their stimulated hand and count the peripheral electrical stimuli they perceived; they were then asked the actual count by the investigator three or four times during the PAS session [28]. The MEPs evoked in the APB were stored for off-line analysis.

2.5. WI Intervention

The WI intervention lasted for 15 min. The subjects wore only swimwear and were seated on a comfortable reclining armchair with a mounted headrest. The subjects assumed the same body position under all interventions, and three belts were secured around the thigh, abdomen, and chest to prevent movement. Both the right and left hands were also fixed in the same relaxed position on the armrests by using a belt to prevent muscle contractions in water. We monitored the EMG for any contraction in the APB and FDI muscle during WI. Participants were instructed to focus their gaze at the wall facing them throughout the experiments to divert their attention from their right hand. For each intervention, the ambient air temperature was maintained at 30 ± 1 °C, and water temperature was maintained at 34 ± 1 °C. The ambient air and water temperatures were controlled to minimize the effect of changing skin temperature. The tank in which the subjects were seated was filled with water up to the axillary level.

2.6. Experimental Design and Parameters

Figure 1 shows the overall design of experiments 1, 2, and 3, including all the interventions and parameters. All trial blocks were conducted in random order for each subject, with an interval of at least 1 week. All experiments were performed in the afternoon to control for circadian influences on motor excitability and plasticity [29].

2.6.1. Experiment 1

First, RMT and AMTs with posterior-to-anterior (PA), anterior-to-posterior (AP), and lateral-to-medial (LM) currents—defined as AMTpa, AMTap, and AMTlm—were measured; thereafter, MEPs during contraction were measured to evaluate the MEP onset latencies. All measurements were performed at the hotspot determined for PA currents, as previous experiments have shown that the direction of the current does not significantly influence the position of the hotspot [30,31].
We measured the onset latency of MEPs using the PA, AP, and LM currents during mild contractions (<10% of the maximum voluntary contraction) of the target muscle (Figure 2), according to the method of Hamada, Murase, Hasan, Balaratnam, and Rothwell [13]. Previous studies found that PA currents preferentially activate early I-waves, whereas AP currents mainly activate later I-waves (e.g., I3) [31,32,33]. Therefore, we decided to use the latency difference between the LM- and AP-evoked MEP onset as a measure of late I-wave recruitment (the longer the latency difference, the more efficient the later I-wave recruitment).
M1 excitability was assessed by measuring the peak-to-peak MEP amplitude from the relaxed right APB muscle in response to a single-pulse TMS with PA current over the contralateral APB hotspot. The test stimulus (TS) intensity was fixed to produce MEPs of approximately 1.0 mV in the right APB muscle at baseline (TS1mV_base). The TS was applied every 5 s. During the experiments, EMG activity of the APB muscle was monitored using an oscilloscope. Trials contaminated with voluntary EMG activity were discarded from the analysis. M1 excitability assessments consisted of baseline measurements, one measurement between the priming session and PAS25 session (“between”), and three “post” measurements after the PAS25 session: immediately, 15 min, and 30 min after PAS25 (post0, post15, and post30, respectively).

2.6.2. Experiment 2

Paired TMS pulses were administered through the same stimulating coil over the left motor cortex to assess SICI. The test stimulus intensity was adjusted to elicit an unconditioned test MEP in the relaxed left APB of approximately 1 mV peak-to-peak amplitude (TS1mV_adjusted). Test stimuli were applied every 5 s. SICI was tested with a subthreshold conditioning stimulus (CSSICI) at 90% of the AMT, applied 3 ms before the TS, as described in previous studies [34,35].
SAI was studied by using a previously described protocol [36,37,38]. Conditioning electrical pulses (CSSAI duration, 200 µs) were applied through a bipolar electrode to the right median nerve at the wrist (cathode proximal). The intensity of the conditioning stimulus was set at approximately three times the sensory threshold. The intensity of the TMS test pulse over the left motor cortex was adjusted to evoke an unconditioned MEP in the relaxed APB of approximately 1 mV peak-to-peak. The conditioning stimulus to the median nerve preceded the TMS test pulse by inter-stimulus intervals (ISIs) set according to the individual latency of the median nerve somatosensory evoked potential N20 component. To record the somatosensory evoked potentials, the active electrode was attached at a site 3 cm posterior to C3 (according to the international 10–20 system), and the reference was set at Fz. A total of 500 responses were averaged to identify the latency of the N20 peak. The ISIs corresponding to the N20 latency plus 2 and 4 ms, and an ISI of 25 ms were investigated [36]. Test stimuli were applied every 5 s. During the experiments, the EMG activity of the APB muscle was monitored using an oscilloscope. Trials contaminated with voluntary EMG activity were discarded.
After SEP measurements, RMT, AMT, TS1mV_base, and CSSICI were established, and the baseline measurement block was started. Each measurement block consisted of five different stimuli: test alone, SICI with an ISI of 3 ms, SAI with an ISI of N20 latency plus 2 or 4 ms, and an ISI of 25 ms. The order was randomized by a computer, and 12 trials of each type were recorded per block. All assessments consisted of baseline measurements and six post measurements after intervention (non-WI and WI). Four post measurements were acquired after intervention with a short break of 1.5 min (post0, post1, post2, and post3), and two additional post measurements were acquired 15 and 30 min after post3 (post4, post5). The post1 and -2 measurements were acquired to examine the cholinergic and GABAergic neural activity, respectively, during the PAS25 session in experiment 1, because cholinergic activity has been shown to enhance PAS25-induced plasticity [24]. However, there is no evidence for changes in GABAergic and cholinergic neuronal activity immediately after WI [39].
During the experiments, the EMG activity of the APB muscle was monitored using an oscilloscope. Trials contaminated with voluntary EMG activities were discarded. The mean amplitudes of the conditioned MEPs in response to these stimulus paradigms were expressed as a percentage of the mean amplitude of the corresponding unconditioned test MEP. These data were averaged across all stimulus paradigms to obtain a grand mean single value of SICI3ms, SAIN20+2, SAIN20+4, and SAI25ms.

2.6.3. Experiment 3

After evaluating RMT, we measured the MEPs evoked by single TMS pulses of increasing stimulus intensity (50%, 80%, 90%, 100%, 110%, 120%, 130%, and 150% of the RMT) in each participant to construct individual MEP recruitment curves. A total of eight pulses were delivered for each stimulus intensity in random order. Stimuli were delivered every 5 s. To avoid startle and reflex responses, we excluded the first MEP for each trial from the analysis. The assessment consisted of baseline measurements during WI and five post-WI measurements—three starting 4 min after WI (post0, post1, and post2) and two additional post assessments 15 and 30 min after post2 (post3 and post4).

2.7. Data Analysis and Statistics

The participant characteristics and both TMS and electrical stimulation (ES) parameters for experiments 1, 2, and 3 are presented in Table 1. The intraclass correlation coefficients (ICCs) were calculated to test the consistency of the latency values in experiment 1 and the inhibited parameters in experiment 2 between the no priming and priming WI trials.

2.7.1. Experiment 1

MEP amplitudes were entered into a two-way repeated-measures ANOVA (rmANOVA) with “trial” (no priming and priming WI) and “time” (baseline, between, post0, post15, and post30) as within-subject factors. The MEP amplitudes normalized to baseline were entered into an rmANOVA, with “trial” (no priming and priming WI) and “time” (baseline, between, post0, post15, and post30) as within-subject factors.
The PAS25 effects were also assessed by the grand averaging of the normalized MEP amplitudes measured at post0–post30, to evaluate the correlations between all the measurement values and the PAS25 effects. Pearson’s correlation coefficients were calculated to measure the strengths of these correlations. Moreover, the responder and non-responder groups were defined operationally according to the presence of grand average PAS25 responses of below and above 1. We primarily aimed to elucidate factors that were predictive of or that regulated M1 plasticity in response to PAS25, and hence, we performed rmANOVA with “trial” and “time” as within-subject factors separately in responders and non-responders.

2.7.2. Experiment 2

SICI3ms, SAIN20+2, SAIN20+4, and SAI25ms were entered into a two-way rmANOVA, with “trial” (CON and WI) and “time” (baseline, post0, post1, post2, post3, post4, and post5) as within-subject factors. We calculated the Pearson’s correlation coefficients between the grand average normalized MEP amplitudes measured at post1−3 in experiment 1 and each inhibitory measure in experiment 2, in order to evaluate the strengths of the associations with PAS25. Because each experiment was conducted on a separate day, the ICCs were calculated to test the consistency of the inhibitory measures in experiment 2 between the CON and WI trials.

2.7.3. Experiment 3

The MEP amplitudes were analyzed by two-way rmANOVA, with “intensity” (from 50% to 150% RMT) and “time” (baseline, during, post0, post1, post2, post3, and post4) as the within-subject factors.
In all analyses using rmANOVA, the Greenhouse–Geisser correction was used, if necessary, to correct for non-sphericity, whereas Tukey’s post hoc tests were used for pair-wise comparisons. A p-value of <0.05 was considered significant. Data were analyzed using the Statistical Software Package (IBM SPSS Version 18, USA). All data are expressed as mean ± standard error of the mean (SEM).

3. Results

All 18 participants completed 2 sessions in experiments 1 and 2, and 9 completed Experiment 3.

3.1. Latency Difference among Different Coil Orientations

As explained in the Methods section, D-wave latency was estimated following the onset of large MEPs evoked by LM stimulation. We estimated I-wave circuit activation during PAS25 by measuring the onset latency of the near-threshold MEPs evoked by AP stimulation, relative to D-wave activation. Finally, we also measured the onset latencies of the near-threshold MEPs evoked by PA stimulation. Table 1 shows the MEP latencies for each TMS stimulus direction and the latency differences.
The interclass correlation coefficients (ICC) calculated to test consistency between latency values measured in the same individual on different days (no priming and priming WI, Figure 3) were 0.921 for PA latency, 0.926 for AP latency, 0.905 for LM latency, 0.852 for the PA–LM latency difference, and 0.874 for the AP–LM latency difference (all p < 0.001), indicating that the data spread was almost entirely due to inter-individual differences in the TMS response.

3.2. PAS25-Induced Plasticity in M1

Figure 3 shows the time course of PAS25-induced plasticity in the no priming trials (Figure 4A) and priming WI trials (Figure 4B) for all 18 participants, as well as the group average in both the trials (Figure 4C).
Two-way rmANOVA of the MEP amplitudes revealed a significant interaction between “trial” and “time” (F [4, 68] = 9.417, p < 0.01), as well as the main effects of both “trial” (F [1, 17] = 13.539, p < 0.01) and “time” (F [4, 68] = 7.256, p < 0.01). Post hoc comparisons revealed significant differences in the values at post15 and post30, compared to those at baseline and the “between” time point only in the WI priming trials (p < 0.05, Bonferroni-corrected). Moreover, there were significant differences in the values between the WI and non-WI (CON) trials at post0, post15, and post30 (p < 0.05, Bonferroni-corrected).
Two-way rmANOVA of the MEP amplitudes normalized to baseline revealed a significant interaction between “trial” and “time” (F [4, 68] = 10.784, p < 0.01), as well as the main effects of “trial” (F [1, 17] = 62.882, p < 0.01) and “time” (F [4, 68] = 6.635, p < 0.01). Post hoc comparisons indicated significant differences in the values at post15 and post30, compared to those at baseline and the “between” time point only in the WI priming trials (p < 0.05, Bonferroni-corrected). There were also significant differences between no priming and WI priming trials at post0, post15, and post30 (p < 0.05, Bonferroni-corrected).
Figure 4D,E shows the group average time courses of PAS25-induced plasticity for responders and non-responders. Two-way rmANOVA of MEP amplitudes normalized to baseline revealed a significant interaction between “trial” and “time” in both responders (F [4, 36] = 5.059, p < 0.01) and non-responders (F [4, 28] = 6.672, p < 0.01). Moreover, the main effects of “trial” and “time” in responders (F [1, 9] = 49.075, p < 0.01 and F [4, 36] = 10.778, p < 0.01) and non-responders (F [1, 7] = 53.766, p < 0.01 and F [4, 28] = 0.355, p = 0.852) were observed. In addition, the post hoc comparisons revealed significant differences in the values at post0, post15, and post30, compared to those at baseline and the “between” time point in the WI priming trials in responders (p < 0.05, Bonferroni-corrected), and even in the values at post15 and post30, compared to those at baseline and the “between” time point in WI priming trials in non-responders (p < 0.05, Bonferroni-corrected). Moreover, there were significant differences in the values between the WI and CON trials at post30 in the responders (p < 0.05, Bonferroni-corrected) and at post15 and post30 in the non-responders (p < 0.05, Bonferroni-corrected).
We tested whether the PAS25 response of each individual in the no priming trials was correlated with any of the baseline physiological measures collected. The PAS25 response in the CON trials was significantly related to RMT (r = −0.472, p < 0.05), AMTpa (r = −0.582, p < 0.05) and TS1mV_base (r = −0.473, p < 0.05), but not to the other baseline physiological measures (Table 2).

3.3. SAI and SICI

Figure 5 shows the time courses of SICI3ms (Figure 5A), SAIN20+2 (Figure 5B), SAIN20+4 (Figure 5C), and SAI25ms (Figure 5D) in both the WI and CON trials. Two-way rmANOVA of the SICI3ms revealed no interaction between “trial” and “time” (F [6, 102] = 0.796, p = 0.576), and no main effects of “trial” (F [1, 17] = 0.026, p = 0.874) and “time” (F [6, 102] = 0.687, p = 0.661) either.
Two-way rmANOVA on the SAIN20+2 revealed a significant interaction between “trial” and “time” (F [6, 102] = 3.255, p < 0.01), as well as the main effects of “trial” (F [1, 17] = 1.251, p = 0.279) and “time” (F [6, 102] = 4.120, p < 0.01). Post hoc comparisons revealed significant differences in the values at post1, compared to those at baseline, as well as in the values at post0 and at post2, compared to those at baseline and post5 only in the WI trials (p < 0.05, Bonferroni-corrected); however, there were no significant pair-wise differences in the CON trials. Moreover, a significant difference was observed in the values between the CON and WI trials at post1 (p < 0.05, Bonferroni-corrected).
Two-way rmANOVA of the SAIN20+4 revealed a significant interaction between “trial” and “time” (F [6, 102] = 2.898, p < 0.05), as well as the main effects of “trial” (F [1, 17] = 4.938, p < 0.05) and “time” (F [6, 102] = 5.606, p < 0.01). Post hoc comparisons revealed significant differences in the values at post1, compared to those at baseline and post0, as well as in the values at post2, compared to those at baseline, post0, post3, post4, and post5 only in the WI trials (p < 0.05, Bonferroni-corrected). Alternatively, there were no significant pair-wise differences in the values in the CON trials. Nevertheless, significant differences were found in the values between the CON and WI trials at post1 and post2 (p < 0.05, Bonferroni-corrected).
Two-way rmANOVA of the SAI25ms revealed a significant interaction between “trial” and “time” (F [6, 102] = 2.307, p < 0.05), as well as the main effect of “trial” (F [1, 17] = 6.200, p < 0.05), but no main effect of “time” (F [6, 102] = 1.714, p = 0.125). Moreover, post hoc comparisons revealed significant differences in the values at post1 and post2, compared to those at baseline only in the WI trials (p < 0.05, Bonferroni-corrected); however, there were no significant differences in the values in the CON trials. In addition, there were significant differences in the values between the CON and WI trials at post1 and post2 (p < 0.05, Bonferroni-corrected).
The ICC values calculated to test the consistency between the inhibitory values measured in the same individual on different days were 0.881 for SICI3ms, 0.893 for SAIN20+2, 0.924 for SAIN20+4, and 0.969 for SAI25ms (all p < 0.001) (Figure 6). We tested whether the PAS25 responses in each individual obtained in experiment 1 were correlated with any of the inhibitory measures collected in experiment 2, and which found a significant correlation between the PAS25 response and SAI25ms (r = −0.534, p < 0.05) (Figure 7, Table 2).

3.4. MEP Recruitment Curve

Figure 8 illustrates the MEP recruitment curves before, during, and after the WI intervention. Two-way rmANOVA revealed no interaction between “intensity” and “time” (F [3.67, 29.33] = 1.08, p = 0.36), and no main effect of “time” (F [2.40, 19.16] = 0.89, p = 0.45), although there was a main effect of “intensity” (F [1.14, 9.12] = 32.94, p < 0.00).

4. Discussion

We present two major findings. First, SAI with an ISI of 25 ms appeared to be predictive of PAS25-induced plasticity, whereas I-wave recruitment did not. Second, the PAS25 response was enhanced by prior WI. This priming effect was observed in 17 participants, irrespective of whether they were responders or non-responders to PAS25 alone, as there was little inter-individual variability in the homeostatic SAI increase after WI. These results support our initial hypothesis that WI priming enhances the PAS25-induced plasticity in the primary motor cortex.
WI can alter numerous physiological parameters and has recently been used to improve the activities of daily living [40,41] and cognitive function in the elderly subjects [42]. Moreover, WI can influence neural activity, such as activation in the primary somatosensory area and multimodal sensory processing [21]. However, to our knowledge, no study had examined the WI priming effect on neural plasticity, despite the observed decrease in cholinergic activity during WI [39]. Priming the activity of specific interneurons synapsing in M1 was found to alter the neural plasticity without modifying M1 excitability [43,44], thus supporting the notion that heterosynaptic plasticity may account for the effects of priming. Experiments 2 and 3 indicated that WI did not alter corticospinal excitability or intracortical excitability, consistent with previous studies [39,45]. In contrast, the cholinergic activity, as evaluated by SAI, significantly increased for at least 8.5 min after WI and returned to baseline at 21.5 min post-WI. Thus, the results of experiment 1 suggest that this cholinergic facilitation phase after WI may enhance the M1 plasticity induced by PAS25.
Consistent with previous studies examining synaptic plasticity in the human motor cortex induced by TMS, we observed a large inter-individual variability in the PAS25 response. In particular, there was no significant response to PAS25 in the entire cohort, and only 55% of participants exhibited MEP facilitation after PAS25, which is within the proportion range noted in previous studies (e.g., 22% [11], 64% [10], and 68% [8]). This inter-individual variability was associated with RMT and AMTpa, but no other initial thresholds, including the indices of I-wave recruitment. The effect of attention [28,46,47] was unlikely because the participants were instructed to focus on the median nerve stimuli during the PAS25 protocol. Circadian effects [29] were also unlikely because all the experiments were performed in the afternoon. Finally, prior and ongoing voluntary muscle contractions [8,48] were carefully monitored and corrupted trials were exempted from the analysis. Therefore, the inter-individual variability presumably reflected the intrinsic differences in the PAS25 response between participants.
Furthermore, responders exhibited lower RMT, AMTpa, and TS1mV_base than non-responders. Muller-Dahlhaus et al. [9] examined the age-dependence of the PAS25 response and found that RMT and TS1mV_base were predictive of PAS25 facilitation in elderly. However, the present study enrolled only young adults. Another possible explanation for the lower thresholds in responders is genetic polymorphisms associated with synaptic plasticity. For example, a single-nucleotide missense polymorphism in the gene encoding brain-derived neurotrophic factor (BDNF Val66Met) modified the experience-dependent motor cortical plasticity [49] and the PAS25 response [50]. This notion is also supported by our findings that responders showed greater M1 excitability than non-responders, and that the threshold correlated negatively with the PAS25 response. Motor thresholds and TS1mV_base may thus reflect a genetically determined endophenotype [51] determining the magnitude and direction of stimulation-induced M1 plasticity in a given individual.
We found no evidence that I-wave recruitment is related to PAS25 response, although a relationship between intermittent theta burst stimulation (iTBS)-induced neural plasticity and I-wave recruitment was reported by using the same methods [13]. A previous study showed that PAS25 led to a pronounced increase in the excitability of the cortical circuits generating later I-waves, whereas the earliest I-wave remained unaffected [52]. We hypothesized that the individual variability in the PAS25 response depends on the relative recruitment of late versus early I-waves. However, our results show that the recruitment of later I-waves is unlikely to predict the PAS25 response. This may be attributed to the distinct mechanisms for PAS25- and iTBS-evoked plasticity. PAS25 is a form of spike timing-dependent plasticity (STDP) [53], whereas TBS responses are related to activity-dependent changes in the synaptic strength between cortical neurons [54]. Moreover, the GABAergic activity involved in SICI was not affected by PAS25, but was decreased after iTBS [55].
The inter-individual variability in the PAS25 response was significantly related to SAI25ms, although it was not correlated with SICI3ms, SAIN20+2, or SAIN20+4. Thus, a higher SAI induced by a longer ISI (25 ms) predicted a greater PAS25 response, possibly because of an optimal interval between the pairing of conditioning afferent stimulation with TMS stimulation. SAI is induced by the increased excitability of GABAergic neurons through the activation of thalamocortical cholinergic projections [56]. The PAS25 induction protocol consists of repeated stimulus pairings with an ISI of 25 ms every 5 s, which increases the late I-wave amplitude without changing the I1 wave [52]. This LTP-like phenomena is reminiscent of STDP and may be produced either by a change in the efficacy of excitatory synaptic connections between the thalamocortical cholinergic projections activated by peripheral nerve stimulation or between P2−P3 pyramidal neurons activated by TMS, or by a change in the excitability of inhibitory neurons directly activated by afferent inputs and projecting onto P5 cells [56]. Considering these shared mechanisms, the present correlation between PAS25 and SAI25ms appears reasonable. Strigaro et al. [57] showed that repeated pairing stimulation with the same ISI is important for the induction of the PAS plasticity response. Unfortunately, we cannot conclude that SAI25ms is a predictor of the PAS25 response due to the lack of data showing whether the PAS changes at other ISIs (N20 + 2 and N20 + 4) can be predicted by SAI responses to the same ISI.
Priming WI enhanced the PAS25 response in almost all participants, although there was no significant MEP facilitation by PAS25 alone due to the high inter-individual variability. Metaplasticity could explain these results. Animal experiments have shown that prior synaptic activity can influence the subsequent response to a stimulation protocol designed to induce LTP or LTD [4], and priming effects have been investigated in humans using non-invasive brain stimulation [18,19,48,58]. The LTP-like response to a facilitatory PAS protocol was decreased when preceded by similar facilitatory PAS, but enhanced when preceded by a PAS protocol that alone induced a LTD-like MEP decrease [19]. These results are consistent with the influence of homeostatic–homosynaptic metaplasticity. However, considering that RMT, MEP size, and SICI did not change 30 min after WI in the present study, the occurrence of metaplasticity among M1 pyramidal neurons and GABAAergic interneurons is less likely than a change in the cholinergic activity.
The significantly increased SAI in experiment 2 indicated the facilitation of cholinergic activity between sensory and motor areas after WI. WI reportedly decreases SAI through widespread somatosensory inputs from the entire body surface [39]. Therefore, we speculate that a homeostatic response of cholinergic neurons between the sensory and motor areas occurs after WI, which could facilitate the PAS25 response. Cholinergic activity is a powerful regulator of synaptic plasticity. In animal studies, cholinergic blockade reduced LTP in the hippocampus, piriform cortex, and neocortex [59,60,61,62,63]; in contrast, in humans, use-dependent plasticity in the motor cortex was facilitated by an acetylcholinesterase inhibitor and blocked by a cholinergic antagonist [64,65]. The regulation of plasticity by cholinergic activity may facilitate the detection of incoming afferent inputs and decrease intrinsic feedback excitability, thus enhancing the encoding of relevant associated information [66,67]. Kuo et al. [24] reported that acetylcholine (Ach) enhanced the synapse-specific cortical excitability increase induced by PAS25 and consolidated the PAS10-induced reduction in motor cortical excitability, whereas it prevented the global excitatory aftereffects produced by anodal tDCS. They described that cholinergic nervous activity improved the efficacy of PAS by enhancing the signal-to-noise ratio, thereby facilitating the processing of meaningful (associative) inputs, and by suppressing non-meaningful/irrelevant asynchronous inputs [68,69,70]. Therefore, the enhanced plasticity following WI may be related to increased cholinergic activity associated with WI-induced homeostatic aftereffects.
In the field of neurological rehabilitation, not all patients can achieve expected results by multiple rehabilitation rounds due to several reasons. However, it should be noted that M1 plasticity plays essential roles in motor learning and memory, which are integral to neurological rehabilitation. Similar to results in previous studies [10,71], the present results show a high inter-individual variability for PAS25 and M1 plasticity. On the basis of these, the present finding that WI priming will promote LTP-like plasticity induction may apply to home-based rehabilitation. WI prior to starting rehabilitation will promote M1 plasticity induction, and this may help to facilitate the effects of rehabilitation. To date, beyond this theory, no direct evidence has been found to show that WI priming facilitates the effects of rehabilitation, hence, further studies are required in this direction.
The present results suggest that WI prior to PAS25 facilitates LTP-like plasticity due to a homeostatic increase in the cholinergic activity. The I-wave recruitment evaluated by the MEP latency difference was not related to the PAS25 plasticity response, in contrast to the TBS response [13].

Author Contributions

D.S. designed whole experiments, wrote the manuscript and draft the figure and tables. K.Y. and A.M. supervised this study. D.S., Y.Y., K.I. and Y.B. conducted experiments. K.Y., H.O. and A.M. reviewed and edited the manuscript and figures. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by JSPS KAKENHI grant number 15K12712 and 18H03134.

Conflicts of Interest

The authors declare no competing financial interests.

References

  1. Martin, S.J.; Grimwood, P.D.; Morris, R.G. Synaptic plasticity and memory: An evaluation of the hypothesis. Annu. Rev. Neurosci. 2000, 23, 649–711. [Google Scholar] [CrossRef] [Green Version]
  2. Sjostrom, P.J.; Turrigiano, G.G.; Nelson, S.B. Rate, timing, and cooperativity jointly determine cortical synaptic plasticity. Neuron 2001, 32, 1149–1164. [Google Scholar] [CrossRef] [Green Version]
  3. Nelson, S.B.; Turrigiano, G.G. Strength through diversity. Neuron 2008, 60, 477–482. [Google Scholar] [CrossRef] [Green Version]
  4. Abraham, W.C.; Bear, M.F. Metaplasticity: The plasticity of synaptic plasticity. Trends Neurosci. 1996, 19, 126–130. [Google Scholar] [CrossRef]
  5. Muller-Dahlhaus, F.; Ziemann, U. Metaplasticity in human cortex. Neuroscientist 2015, 21, 185–202. [Google Scholar] [CrossRef] [PubMed]
  6. Stefan, K.; Kunesch, E.; Cohen, L.G.; Benecke, R.; Classen, J. Induction of plasticity in the human motor cortex by paired associative stimulation. Brain 2000, 123, 572–584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Fratello, F.; Veniero, D.; Curcio, G.; Ferrara, M.; Marzano, C.; Moroni, F.; Pellicciari, M.C.; Bertini, M.; Rossini, P.M.; De Gennaro, L. Modulation of corticospinal excitability by paired associative stimulation: Reproducibility of effects and intraindividual reliability. Clin. Neurophysiol. 2006, 117, 2667–2674. [Google Scholar] [CrossRef] [PubMed]
  8. Stefan, K.; Wycislo, M.; Gentner, R.; Schramm, A.; Naumann, M.; Reiners, K.; Classen, J. Temporary occlusion of associative motor cortical plasticity by prior dynamic motor training. Cereb. Cortex 2006, 16, 376–385. [Google Scholar] [CrossRef] [Green Version]
  9. Muller-Dahlhaus, J.F.; Orekhov, Y.; Liu, Y.; Ziemann, U. Interindividual variability and age-dependency of motor cortical plasticity induced by paired associative stimulation. Exp. Brain Res. 2008, 187, 467–475. [Google Scholar] [CrossRef]
  10. Hamada, M.; Galea, J.M.; Di Lazzaro, V.; Mazzone, P.; Ziemann, U.; Rothwell, J.C. Two distinct interneuron circuits in human motor cortex are linked to different subsets of physiological and behavioral plasticity. J. Neurosci. 2014, 34, 12837–12849. [Google Scholar] [CrossRef]
  11. Lopez-Alonso, V.; Cheeran, B.; Rio-Rodriguez, D.; Fernandez-Del-Olmo, M. Inter-individual variability in response to non-invasive brain stimulation paradigms. Brain Stimul. 2014, 7, 372–380. [Google Scholar] [CrossRef] [PubMed]
  12. Ridding, M.C.; Ziemann, U. Determinants of the induction of cortical plasticity by non-invasive brain stimulation in healthy subjects. J. Physiol. 2010, 588, 2291–2304. [Google Scholar] [CrossRef] [PubMed]
  13. Hamada, M.; Murase, N.; Hasan, A.; Balaratnam, M.; Rothwell, J.C. The role of interneuron networks in driving human motor cortical plasticity. Cereb. Cortex 2013, 23, 1593–1605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Murase, N.; Cengiz, B.; Rothwell, J.C. Inter-individual variation in the after-effect of paired associative stimulation can be predicted from short-interval intracortical inhibition with the threshold tracking method. Brain Stimul. 2015, 8, 105–113. [Google Scholar] [CrossRef] [PubMed]
  15. Suzuki, M.; Kirimoto, H.; Sugawara, K.; Watanabe, M.; Shimizu, S.; Ishizaka, I.; Yamada, S.; Matsunaga, A.; Fukuda, M.; Onishi, H. Induction of cortical plasticity for reciprocal muscles by paired associative stimulation. Brain Behav. 2014, 4, 822–832. [Google Scholar] [CrossRef] [PubMed]
  16. Udupa, K.; Ni, Z.; Gunraj, C.; Chen, R. Interactions between short latency afferent inhibition and long interval intracortical inhibition. Exp. Brain Res. 2009, 199, 177–183. [Google Scholar] [CrossRef]
  17. Udupa, K.; Ni, Z.; Gunraj, C.; Chen, R. Effects of short-latency afferent inhibition on short-interval intracortical inhibition. J. Neurophysiol. 2014, 111, 1350–1361. [Google Scholar] [CrossRef] [Green Version]
  18. Iyer, M.B.; Schleper, N.; Wassermann, E.M. Priming stimulation enhances the depressant effect of low-frequency repetitive transcranial magnetic stimulation. J. Neurosci. 2003, 23, 10867–10872. [Google Scholar] [CrossRef]
  19. Muller, J.F.; Orekhov, Y.; Liu, Y.; Ziemann, U. Homeostatic plasticity in human motor cortex demonstrated by two consecutive sessions of paired associative stimulation. Eur. J. Neurosci. 2007, 25, 3461–3468. [Google Scholar] [CrossRef]
  20. Sato, D.; Onishi, H.; Yamashiro, K.; Iwabe, T.; Shimoyama, Y.; Maruyama, A. Water immersion to the femur level affects cerebral cortical activity in humans: Functional near-infrared spectroscopy study. Brain Topogr. 2012, 25, 220–227. [Google Scholar] [CrossRef]
  21. Sato, D.; Yamashiro, K.; Onishi, H.; Shimoyama, Y.; Yoshida, T.; Maruyama, A. The effect of water immersion on short-latency somatosensory evoked potentials in human. BMC Neurosci. 2012, 13, 13. [Google Scholar] [CrossRef] [Green Version]
  22. Di Lazzaro, V.; Oliviero, A.; Profice, P.; Pennisi, M.A.; Di Giovanni, S.; Zito, G.; Tonali, P.; Rothwell, J.C. Muscarinic receptor blockade has differential effects on the excitability of intracortical circuits in the human motor cortex. Exp. Brain Res. 2000, 135, 455–461. [Google Scholar] [CrossRef]
  23. Brown, R.E.; Basheer, R.; McKenna, J.T.; Strecker, R.E.; McCarley, R.W. Control of sleep and wakefulness. Physiol. Rev. 2012, 92, 1087–1187. [Google Scholar] [CrossRef] [Green Version]
  24. Kuo, M.F.; Grosch, J.; Fregni, F.; Paulus, W.; Nitsche, M.A. Focusing effect of acetylcholine on neuroplasticity in the human motor cortex. J. Neurosci. 2007, 27, 14442–14447. [Google Scholar] [CrossRef]
  25. Rossi, S.; Hallett, M.; Rossini, P.M.; Pascual-Leone, A. Safety, ethical considerations, and application guidelines for the use of transcranial magnetic stimulation in clinical practice and research. Clin. Neurophysiol. 2009, 120, 2008–2039. [Google Scholar] [CrossRef] [Green Version]
  26. Buharin, V.E.; Shinohara, M. Corticospinal excitability for flexor carpi radialis decreases with baroreceptor unloading during intentional co-contraction with opposing forearm muscles. Exp. Brain Res. 2019, 237, 1947–1958. [Google Scholar] [CrossRef]
  27. Sato, D.; Yamashiro, K.; Onishi, H.; Yasuhiro, B.; Shimoyama, Y.; Maruyama, A. Whole-hand water flow stimulation increases motor cortical excitability: A study of transcranial magnetic stimulation and movement-related cortical potentials. J. Neurophysiol. 2015, 113, 822–833. [Google Scholar] [CrossRef] [Green Version]
  28. Stefan, K.; Wycislo, M.; Classen, J. Modulation of associative human motor cortical plasticity by attention. J. Neurophysiol. 2004, 92, 66–72. [Google Scholar] [CrossRef]
  29. Sale, M.V.; Ridding, M.C.; Nordstrom, M.A. Cortisol inhibits neuroplasticity induction in human motor cortex. J. Neurosci. 2008, 28, 8285–8293. [Google Scholar] [CrossRef] [Green Version]
  30. Arai, N.; Okabe, S.; Furubayashi, T.; Terao, Y.; Yuasa, K.; Ugawa, Y. Comparison between short train, monophasic and biphasic repetitive transcranial magnetic stimulation (rTMS) of the human motor cortex. Clin. Neurophysiol. 2005, 116, 605–613. [Google Scholar] [CrossRef]
  31. Sakai, K.; Ugawa, Y.; Terao, Y.; Hanajima, R.; Furubayashi, T.; Kanazawa, I. Preferential activation of different I waves by transcranial magnetic stimulation with a figure-of-eight-shaped coil. Exp. Brain Res. 1997, 113, 24–32. [Google Scholar] [CrossRef]
  32. Day, B.L.; Rothwell, J.C.; Thompson, P.D.; Dick, J.P.; Cowan, J.M.; Berardelli, A.; Marsden, C.D. Motor cortex stimulation in intact man. 2. Multiple descending volleys. Brain 1987, 110, 1191–1209. [Google Scholar] [CrossRef]
  33. Rothwell, J.C.; Thompson, P.D.; Day, B.L.; Dick, J.P.; Kachi, T.; Cowan, J.M.; Marsden, C.D. Motor cortex stimulation in intact man. 1. General characteristics of EMG responses in different muscles. Brain 1987, 110, 1173–1190. [Google Scholar] [CrossRef]
  34. Kujirai, T.; Caramia, M.D.; Rothwell, J.C.; Day, B.L.; Thompson, P.D.; Ferbert, A.; Wroe, S.; Asselman, P.; Marsden, C.D. Corticocortical inhibition in human motor cortex. J. Physiol. 1993, 471, 501–519. [Google Scholar] [CrossRef]
  35. Ziemann, U.; Rothwell, J.C.; Ridding, M.C. Interaction between intracortical inhibition and facilitation in human motor cortex. J. Physiol. 1996, 496, 873–881. [Google Scholar] [CrossRef]
  36. Tokimura, H.; Di Lazzaro, V.; Tokimura, Y.; Oliviero, A.; Profice, P.; Insola, A.; Mazzone, P.; Tonali, P.; Rothwell, J.C. Short latency inhibition of human hand motor cortex by somatosensory input from the hand. J. Physiol. 2000, 523, 503–513. [Google Scholar] [CrossRef]
  37. Di Lazzaro, V.; Oliviero, A.; Saturno, E.; Dileone, M.; Pilato, F.; Nardone, R.; Ranieri, F.; Musumeci, G.; Fiorilla, T.; Tonali, P. Effects of lorazepam on short latency afferent inhibition and short latency intracortical inhibition in humans. J. Physiol. 2005, 564, 661–668. [Google Scholar] [CrossRef]
  38. Di Lazzaro, V.; Pilato, F.; Dileone, M.; Profice, P.; Ranieri, F.; Ricci, V.; Bria, P.; Tonali, P.A.; Ziemann, U. Segregating two inhibitory circuits in human motor cortex at the level of GABAA receptor subtypes: A TMS study. Clin. Neurophysiol. 2007, 118, 2207–2214. [Google Scholar] [CrossRef]
  39. Sato, D.; Yamashiro, K.; Yoshida, T.; Onishi, H.; Shimoyama, Y.; Maruyama, A. Effects of water immersion on short- and long-latency afferent inhibition, short-interval intracortical inhibition, and intracortical facilitation. Clin. Neurophysiol. 2013, 124, 1846–1852. [Google Scholar] [CrossRef]
  40. Sato, D.; Kaneda, K.; Wakabayashi, H.; Nomura, T. Comparison of 2-year effects of once and twice weekly water exercise on activities of daily living ability of community dwelling frail elderly. Arch. Gerontol. Geriatr. 2009, 49, 123–128. [Google Scholar] [CrossRef]
  41. Sato, D.; Kaneda, K.; Wakabayashi, H.; Shimoyama, Y.; Baba, Y.; Nomura, T. Comparison of once and twice weekly water exercise on various bodily functions in community-dwelling frail elderly requiring nursing care. Arch. Gerontol. Geriatr. 2011, 52, 331–335. [Google Scholar] [CrossRef]
  42. Sato, D.; Seko, C.; Hashitomi, T.; Sengoku, Y.; Nomura, T. Differential effects of water-based exercise on the cognitive function in independent elderly adults. Aging Clin. Exp. Res. 2015, 27, 149–159. [Google Scholar] [CrossRef]
  43. Hamada, M.; Hanajima, R.; Terao, Y.; Okabe, S.; Nakatani-Enomoto, S.; Furubayashi, T.; Matsumoto, H.; Shirota, Y.; Ohminami, S.; Ugawa, Y. Primary motor cortical metaplasticity induced by priming over the supplementary motor area. J. Physiol. 2009, 587, 4845–4862. [Google Scholar] [CrossRef]
  44. Popa, T.; Velayudhan, B.; Hubsch, C.; Pradeep, S.; Roze, E.; Vidailhet, M.; Meunier, S.; Kishore, A. Cerebellar processing of sensory inputs primes motor cortex plasticity. Cereb. Cortex 2013, 23, 305–314. [Google Scholar] [CrossRef] [Green Version]
  45. Sato, D.; Yamashiro, K.; Onishi, H.; Baba, Y.; Nakazawa, S.; Shimoyama, Y.; Maruyama, A. Whole-body water flow stimulation to the lower limbs modulates excitability of primary motor cortical regions innervating the hands: A transcranial magnetic stimulation study. PLoS ONE 2014, 9, 102472. [Google Scholar] [CrossRef]
  46. Kamke, M.R.; Hall, M.G.; Lye, H.F.; Sale, M.V.; Fenlon, L.R.; Carroll, T.J.; Riek, S.; Mattingley, J.B. Visual attentional load influences plasticity in the human motor cortex. J. Neurosci. 2012, 32, 7001–7008. [Google Scholar] [CrossRef]
  47. Kamke, M.R.; Ryan, A.E.; Sale, M.V.; Campbell, M.E.; Riek, S.; Carroll, T.J.; Mattingley, J.B. Visual spatial attention has opposite effects on bidirectional plasticity in the human motor cortex. J. Neurosci. 2014, 34, 1475–1480. [Google Scholar] [CrossRef] [Green Version]
  48. Ziemann, U.; Ilic, T.V.; Pauli, C.; Meintzschel, F.; Ruge, D. Learning modifies subsequent induction of long-term potentiation-like and long-term depression-like plasticity in human motor cortex. J. Neurosci. 2004, 24, 1666–1672. [Google Scholar] [CrossRef]
  49. Kleim, J.A.; Chan, S.; Pringle, E.; Schallert, K.; Procaccio, V.; Jimenez, R.; Cramer, S.C. BDNF val66met polymorphism is associated with modified experience-dependent plasticity in human motor cortex. Nat. Neurosci. 2006, 9, 735–737. [Google Scholar] [CrossRef] [Green Version]
  50. Cheeran, B.; Talelli, P.; Mori, F.; Koch, G.; Suppa, A.; Edwards, M.; Houlden, H.; Bhatia, K.; Greenwood, R.; Rothwell, J.C. A common polymorphism in the brain-derived neurotrophic factor gene (BDNF) modulates human cortical plasticity and the response to rTMS. J. Physiol. 2008, 586, 5717–5725. [Google Scholar] [CrossRef]
  51. Wassermann, E.M. Variation in the response to transcranial magnetic brain stimulation in the general population. Clin. Neurophysiol. 2002, 113, 1165–1171. [Google Scholar] [CrossRef]
  52. Di Lazzaro, V.; Dileone, M.; Pilato, F.; Profice, P.; Oliviero, A.; Mazzone, P.; Insola, A.; Capone, F.; Ranieri, F.; Tonali, P.A. Associative motor cortex plasticity: Direct evidence in humans. Cereb. Cortex 2009, 19, 2326–2330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Muller-Dahlhaus, F.; Ziemann, U.; Classen, J. Plasticity resembling spike-timing dependent synaptic plasticity: The evidence in human cortex. Front. Synaptic Neurosci. 2010, 2, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hallett, M. Transcranial magnetic stimulation: A primer. Neuron 2007, 55, 187–199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Huang, Y.Z.; Edwards, M.J.; Rounis, E.; Bhatia, K.P.; Rothwell, J.C. Theta burst stimulation of the human motor cortex. Neuron 2005, 45, 201–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Di Lazzaro, V.; Profice, P.; Ranieri, F.; Capone, F.; Dileone, M.; Oliviero, A.; Pilato, F. I-wave origin and modulation. Brain Stimul. 2012, 5, 512–525. [Google Scholar] [CrossRef] [PubMed]
  57. Strigaro, G.; Hamada, M.; Murase, N.; Cantello, R.; Rothwell, J.C. Interaction between different interneuron networks involved in human associative plasticity. Brain Stimul. 2014, 7, 658–664. [Google Scholar] [CrossRef]
  58. Todd, G.; Flavel, S.C.; Ridding, M.C. Priming theta-burst repetitive transcranial magnetic stimulation with low- and high-frequency stimulation. Exp. Brain Res. 2009, 195, 307–315. [Google Scholar] [CrossRef]
  59. Hirotsu, I.; Hori, N.; Katsuda, N.; Ishihara, T. Effect of anticholinergic drug on long-term potentiation in rat hippocampal slices. Brain Res. 1989, 482, 194–197. [Google Scholar] [CrossRef]
  60. Blitzer, R.D.; Gil, O.; Landau, E.M. Cholinergic stimulation enhances long-term potentiation in the CA1 region of rat hippocampus. Neurosci. Lett. 1990, 119, 207–210. [Google Scholar] [CrossRef]
  61. Huerta, P.T.; Lisman, J.E. Heightened synaptic plasticity of hippocampal CA1 neurons during a cholinergically induced rhythmic state. Nature 1993, 364, 723–725. [Google Scholar] [CrossRef] [PubMed]
  62. Brocher, S.; Artola, A.; Singer, W. Agonists of cholinergic and noradrenergic receptors facilitate synergistically the induction of long-term potentiation in slices of rat visual cortex. Brain Res. 1992, 573, 27–36. [Google Scholar] [CrossRef]
  63. Hasselmo, M.E.; Barkai, E. Cholinergic modulation of activity-dependent synaptic plasticity in the piriform cortex and associative memory function in a network biophysical simulation. J. Neurosci. 1995, 15, 6592–6604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Meintzschel, F.; Ziemann, U. Modification of practice-dependent plasticity in human motor cortex by neuromodulators. Cereb. Cortex 2006, 16, 1106–1115. [Google Scholar] [CrossRef] [PubMed]
  65. Sawaki, L.; Boroojerdi, B.; Kaelin-Lang, A.; Burstein, A.H.; Butefisch, C.M.; Kopylev, L.; Davis, B.; Cohen, L.G. Cholinergic influences on use-dependent plasticity. J. Neurophysiol. 2002, 87, 166–171. [Google Scholar] [CrossRef]
  66. Blokland, A.; Honig, W.; Raaijmakers, W.G. Effects of intra-hippocampal scopolamine injections in a repeated spatial acquisition task in the rat. Psychopharmacology 1992, 109, 373–376. [Google Scholar] [CrossRef]
  67. Winters, B.D.; Bussey, T.J. Removal of cholinergic input to perirhinal cortex disrupts object recognition but not spatial working memory in the rat. Eur. J. Neurosci. 2005, 21, 2263–2270. [Google Scholar] [CrossRef]
  68. Patil, M.M.; Linster, C.; Lubenov, E.; Hasselmo, M.E. Cholinergic agonist carbachol enables associative long-term potentiation in piriform cortex slices. J. Neurophysiol. 1998, 80, 2467–2474. [Google Scholar] [CrossRef]
  69. Rasmusson, D.D.; Dykes, R.W. Long-term enhancement of evoked potentials in cat somatosensory cortex produced by co-activation of the basal forebrain and cutaneous receptors. Exp. Brain Res. 1988, 70, 276–286. [Google Scholar] [CrossRef]
  70. Tremblay, N.; Warren, R.A.; Dykes, R.W. Electrophysiological studies of acetylcholine and the role of the basal forebrain in the somatosensory cortex of the cat. II. Cortical neurons excited by somatic stimuli. J. Neurophysiol. 1990, 64, 1212–1222. [Google Scholar] [CrossRef] [Green Version]
  71. Hordacre, B.; Goldsworthy, M.R.; Vallence, A.M.; Darvishi, S.; Moezzi, B.; Hamada, M.; Rothwell, J.C.; Ridding, M.C. Variability in neural excitability and plasticity induction in the human cortex: A brain stimulation study. Brain Stimul. 2017, 10, 588–595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijerph 17 00215 g001 550
Figure 1. The experimental protocol. (A) Experiment 1 examined the priming effect of water immersion (WI) on paired associative stimulation with ISI of 25 ms (PAS25) by comparison between no priming and priming WI trials. After measuring onset of motor-evoked potential (MEP) latency induced by single pulse transcranial magnetic stimulation (TMS) with posterior-to-anterior (PA), anterior-to-posterior (AP), and lateral-to-medial (LM) current during a little contraction, MEPTEST and resting motor threshold (RMT) were measured at baseline, before PAS25 (between), and at 0 min, 15 min, and 30 min after PAS25 (post0, post15, and post30). (B) Experiment 2 examined the effects of WI on short-interval intracortical inhibition (SICI) and short-latency afferent inhibitions (SAIs) by comparison between non-WI (CON) and WI trials. After short latency evoked potential (SEP) measurement to determine stimulus paradigm of SAIs, SICI and SAIs were measured before (baseline) and after intervention (post0, post1, post2, post3, post4, and post5). (C) Experiment 3 investigated the changes in MEP recruitment curve before (baseline), during (during), and after WI (post0, post1, post2, post3, and post4).
Figure 1. The experimental protocol. (A) Experiment 1 examined the priming effect of water immersion (WI) on paired associative stimulation with ISI of 25 ms (PAS25) by comparison between no priming and priming WI trials. After measuring onset of motor-evoked potential (MEP) latency induced by single pulse transcranial magnetic stimulation (TMS) with posterior-to-anterior (PA), anterior-to-posterior (AP), and lateral-to-medial (LM) current during a little contraction, MEPTEST and resting motor threshold (RMT) were measured at baseline, before PAS25 (between), and at 0 min, 15 min, and 30 min after PAS25 (post0, post15, and post30). (B) Experiment 2 examined the effects of WI on short-interval intracortical inhibition (SICI) and short-latency afferent inhibitions (SAIs) by comparison between non-WI (CON) and WI trials. After short latency evoked potential (SEP) measurement to determine stimulus paradigm of SAIs, SICI and SAIs were measured before (baseline) and after intervention (post0, post1, post2, post3, post4, and post5). (C) Experiment 3 investigated the changes in MEP recruitment curve before (baseline), during (during), and after WI (post0, post1, post2, post3, and post4).
Ijerph 17 00215 g001
Ijerph 17 00215 g002 550
Figure 2. Representative MEP waveform and measurement of MEP latency. (A) Schematic representation of the coil orientations and typical examples of MEPs during contraction by each stimulus. (B) Arrow indicates the timing of TMS and arrowhead indicates the onset of MEPs. The PA–LM latency difference was 1.3 ms in this case, whereas AP–LM was 3 ms, compatible with known latency differences between D- and I1-waves or I3-waves [32].
Figure 2. Representative MEP waveform and measurement of MEP latency. (A) Schematic representation of the coil orientations and typical examples of MEPs during contraction by each stimulus. (B) Arrow indicates the timing of TMS and arrowhead indicates the onset of MEPs. The PA–LM latency difference was 1.3 ms in this case, whereas AP–LM was 3 ms, compatible with known latency differences between D- and I1-waves or I3-waves [32].
Ijerph 17 00215 g002
Ijerph 17 00215 g003 550
Figure 3. Test-retest reliability of each MEP latency (A–C) and latency difference (D,E). Black lines show regression lines in each parameter. The x- and y-axis show each value in no priming and priming WI trials. Intraclass correlation coefficient (ICC) values (detail in results) mean higher test–retest reliability in all parameters.
Figure 3. Test-retest reliability of each MEP latency (A–C) and latency difference (D,E). Black lines show regression lines in each parameter. The x- and y-axis show each value in no priming and priming WI trials. Intraclass correlation coefficient (ICC) values (detail in results) mean higher test–retest reliability in all parameters.
Ijerph 17 00215 g003
Ijerph 17 00215 g004 550
Figure 4. Time course of PAS25. (A,B) show individual responses to PAS25 in no priming and priming WI trials, respectively. Though there was high inter-individual variability for PAS25 response, (C) shows significant group difference of PAS25 response—increased in WI trial and unchanged in CON trial. (D,E) show PAS25 response in responders and non-responders. In responders, PAS25 response was facilitated in WI trial compared to CON trial. In non-responders, expected PAS25 responses were present in WI trial, whereas there were opposite responses or no response for PAS25 in CON trial. The x- and y-axes show measurement time points and MEP amplitude normalized to baseline, respectively. Asterisk (*) represents significant difference compared to baseline and “between” time points (p < 0.05). Dagger (†) shows significant difference between no priming and priming WI trials (p < 0.05).
Figure 4. Time course of PAS25. (A,B) show individual responses to PAS25 in no priming and priming WI trials, respectively. Though there was high inter-individual variability for PAS25 response, (C) shows significant group difference of PAS25 response—increased in WI trial and unchanged in CON trial. (D,E) show PAS25 response in responders and non-responders. In responders, PAS25 response was facilitated in WI trial compared to CON trial. In non-responders, expected PAS25 responses were present in WI trial, whereas there were opposite responses or no response for PAS25 in CON trial. The x- and y-axes show measurement time points and MEP amplitude normalized to baseline, respectively. Asterisk (*) represents significant difference compared to baseline and “between” time points (p < 0.05). Dagger (†) shows significant difference between no priming and priming WI trials (p < 0.05).
Ijerph 17 00215 g004
Ijerph 17 00215 g005 550
Figure 5. Time course of SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D) in both WI and CON trials. (A) shows the change in SICI3ms before and after WI, which, as shown, did not change throughout the experiment. (B–D) present the change in SAIN20+2, SAIN20+4, and SAI25ms before and after WI, which all showed temporal increase in SAI after WI irrespective of ISI. The x- and y-axes show measurement time points and each value, respectively. Asterisk (*) represents significant difference among each time points (p < 0.05). Dagger (†) shows significant difference between non-WI and WI trials (p < 0.05).
Figure 5. Time course of SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D) in both WI and CON trials. (A) shows the change in SICI3ms before and after WI, which, as shown, did not change throughout the experiment. (B–D) present the change in SAIN20+2, SAIN20+4, and SAI25ms before and after WI, which all showed temporal increase in SAI after WI irrespective of ISI. The x- and y-axes show measurement time points and each value, respectively. Asterisk (*) represents significant difference among each time points (p < 0.05). Dagger (†) shows significant difference between non-WI and WI trials (p < 0.05).
Ijerph 17 00215 g005
Ijerph 17 00215 g006 550
Figure 6. Test–retest reliability of SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D). Black lines show regression lines in each parameter. The x- and y-axes show each value in CON and WI trials. ICC values (detail in results) mean higher test–retest reliability in all parameters.
Figure 6. Test–retest reliability of SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D). Black lines show regression lines in each parameter. The x- and y-axes show each value in CON and WI trials. ICC values (detail in results) mean higher test–retest reliability in all parameters.
Ijerph 17 00215 g006
Ijerph 17 00215 g007 550
Figure 7. The relationships between PAS25 response and SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D). Black lines show regression lines in each parameter. The x- and y-axes show PAS25 response and each value in CON trials, respectively. As shown in (A–C), there was no correlation between PAS25 response and each value (SICI3ms, SAIN20+2, SAIN20+4). On the other hand, (D) shows a significant negative correlation between PAS25 response and SAI25ms, which indicates that SAI25ms would be a predictor for PAS25 response.
Figure 7. The relationships between PAS25 response and SICI3ms (A), SAIN20+2 (B), SAIN20+4 (C), and SAI25ms (D). Black lines show regression lines in each parameter. The x- and y-axes show PAS25 response and each value in CON trials, respectively. As shown in (A–C), there was no correlation between PAS25 response and each value (SICI3ms, SAIN20+2, SAIN20+4). On the other hand, (D) shows a significant negative correlation between PAS25 response and SAI25ms, which indicates that SAI25ms would be a predictor for PAS25 response.
Ijerph 17 00215 g007
Ijerph 17 00215 g008 550
Figure 8. MEP recruitment curves before, during, and after WI. The x- and y-axes show stimulus intensity of TMS (% resting MT) and MEP amplitude, respectively. MEP recruitment curves were measured before, during, and after WI. There was no change in MEP recruitment curve, which indicates that WI did not change corticospinal excitability.
Figure 8. MEP recruitment curves before, during, and after WI. The x- and y-axes show stimulus intensity of TMS (% resting MT) and MEP amplitude, respectively. MEP recruitment curves were measured before, during, and after WI. There was no change in MEP recruitment curve, which indicates that WI did not change corticospinal excitability.
Ijerph 17 00215 g008
Table
Table 1. Participant characteristics and parameters of TMS and ES in experiments 1 to 3.
Table 1. Participant characteristics and parameters of TMS and ES in experiments 1 to 3.
Experiment 1
Participants’ Characteristics
Number18 (female 3, male 15)
Age (years old)21.37 ± 0.18
TMS and ES Parameter
No PrimingPriming WI
RMT (%)42.28 ± 1.9340.72 ± 1.68
AMTpa (%)34.83 ± 1.5532.61 ± 1.29
AMTap (%)49.89 ± 2.2447.50 ± 2.00
AMTlm (%)37.56 ± 1.6736.50 ± 1.26
TS1mV_base52.72 ± 2.3752.83 ± 2.23
ST (mA)3.63 ± 0.283.73 ± 0.34
MEP Latency
PA (ms)22.10 ± 0.2621.96 ± 0.25
AP (ms)24.17 ± 0.3923.98 ± 0.38
LM (ms)20.77 ± 0.2620.59 ± 0.25
Latency Differences
PA–LM (ms)1.33 ± 0.171.38 ± 0.19
AP–LM (ms)3.40 ± 0.323.39 ± 0.31
AP–PA (ms)2.08 ± 0.252.02 ± 0.28
Experiment 2
Participants’ Characteristics
Number18 (female 3, male 15)
Age (years old)21.37 ± 0.18
N20 (msec)18.62 ± 0.16
TMS and ES Parameter
CON TrialWI Trial
RMT (%)45.29 ± 1.6544.00 ± 1.02
AMTpa (%)35.12 ± 1.2633.61 ± 0.83
TS1mV_base58.88 ± 1.9957.44 ± 1.92
ST (mA)3.35 ± 0.213.67 ± 0.20
CSSICI (%)31.46 ± 1.2230.48 ± 0.71
CSSAI (mA)10.04 ± 0.6410.97 ± 0.55
Experiment 3
Participants’ Characteristics
Number9 (female 1, male 8)
Age (y.o.)21.89 ± 0.31
RMT (%)42.71 ± 1.77
ES: electrical stimulation; AMT: active motor threshold; TS: test stimulus; PA: posterior-to-anterior; AP: anterior-to-posterior; LM: lateral-to-medial; ST: sensory threshold; CS: conditioning stiumulus.
Table
Table 2. Relationship between PAS25 response and each baseline value in experiments 1 and 2.
Table 2. Relationship between PAS25 response and each baseline value in experiments 1 and 2.
Experiment 1
rp
RMT−0.4720.048 *
AMTpa−0.5820.011 *
AMTap−0.3510.154
AMTlm−0.3480.157
TS1mV_base−0.4370.047 *
ST0.2880.246
PA latency0.1690.502
AP latency0.0400.876
LM latency0.0950.706
PA–LM0.1270.615
AP–LM−0.0190.940
Experiment 2
rp
SICI3ms−0.240.333
SAIN20+2−0.150.566
SAIN20+4−0.230.368
SAI25ms−0.530.022 *
* Significant correlation with PAS25 response.

Share and Cite

MDPI and ACS Style

Sato, D.; Yamashiro, K.; Yamazaki, Y.; Ikarashi, K.; Onishi, H.; Baba, Y.; Maruyama, A. Priming Effects of Water Immersion on Paired Associative Stimulation-Induced Neural Plasticity in the Primary Motor Cortex. Int. J. Environ. Res. Public Health 2020, 17, 215. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17010215

AMA Style

Sato D, Yamashiro K, Yamazaki Y, Ikarashi K, Onishi H, Baba Y, Maruyama A. Priming Effects of Water Immersion on Paired Associative Stimulation-Induced Neural Plasticity in the Primary Motor Cortex. International Journal of Environmental Research and Public Health. 2020; 17(1):215. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17010215

Chicago/Turabian Style

Sato, Daisuke, Koya Yamashiro, Yudai Yamazaki, Koyuki Ikarashi, Hideaki Onishi, Yasuhiro Baba, and Atsuo Maruyama. 2020. "Priming Effects of Water Immersion on Paired Associative Stimulation-Induced Neural Plasticity in the Primary Motor Cortex" International Journal of Environmental Research and Public Health 17, no. 1: 215. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17010215

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop