Brain Response Induced with Paired Associative Stimulation Is Related to Repetition Suppression of Motor Evoked Potential

Repetition suppression (RS), i.e., the reduction of neuronal activity upon repetition of an external stimulus, can be demonstrated in the motor system using transcranial magnetic stimulation (TMS). We evaluated the RS in relation to the neuroplastic changes induced by paired associative stimulation (PAS). An RS paradigm, consisting of 20 trains of four identical suprathreshold TMS pulses 1 s apart, was assessed for motor-evoked potentials (MEPs) in 16 healthy subjects, before and following (at 0, 10, and 20 min) a common PAS protocol. For analysis, we divided RS into two components: (1) the ratio of the second MEP amplitude to the first one in RS trains, i.e., the “dynamic” component, and (2) the mean of the second to fourth MEP amplitudes, i.e., the “stable” component. Following PAS, five subjects showed change in the dynamic RS component. However, nearly all the individuals (n = 14) exhibited change in the stable component (p < 0.05). The stable component was similar between subjects showing increased MEPs and those showing decreased MEPs at this level (p = 0.254). The results suggest the tendency of the brain towards a stable state, probably free from the ongoing dynamics, following PAS.


Introduction
Owing to its dynamicity, the brain responds to an intense, novel stimulus with enhanced, transient neural activity. This rapid response, referred to as a startle, is considered to play a critical function in promoting survival [1]. However, exposure to a higher number of identical sensory stimuli yields attenuation of neural activity in the responding network, a phenomenon known as repetition suppression (RS) [2]. RS has been well-characterized across several brain regions, employing various stimulus categories and modalities [3][4][5][6]. In the motor system, RS has been demonstrated as a decrement in the amplitude of motor-evoked potentials (MEPs) when transcranial magnetic stimulation (TMS) is applied to an optimal motor cortex location [7,8]. Although it has been suggested that the attenuation observed in RS may serve to provide an energy-efficient neuronal information processing [9], the exact mechanisms underlying RS have remained elusive. RS was initially portrayed merely as an expression of bottom-up mechanisms [2,3,10,11]. However, more recent theories have emphasized the role of top-down mechanisms within a predictive coding scheme, relying on iterative The stimulation procedure was initiated by locating the optimal motor representation of the right abductor pollicis brevis (APB) muscle, i.e., APB "hotspot", with the corresponding optimized coil orientation. The hotspot was defined as the cortical site repeatedly eliciting the greatest peak-to-peak MEP responses compared to adjacent stimulation sites. Once the hotspot was determined, the resting motor threshold (rMT) was identified at this cortical site using a system-integrated iterative threshold assessment tool [29]. In the RS paradigm, trials of four TMS stimuli were applied over the APB hotspot, with an inter-stimulus interval (ISI) of 1 s, at an intensity of 120% rMT. The RS paradigm, comprising 20 trials of four single biphasic TMS pulses, was employed with an inter-train interval (ITI) of 17 s [30], before (RS-baseline) and immediately (0 min), 10 min, and 20 min after the PAS intervention ( Figure 1).
We recorded MEPs via an integrated electromyography (EMG) system (Nexstim Plc) at a sampling frequency of 3 kHz. A pair of disposable Ag-Cl electrodes was utilized, with the active electrode over the belly of the APB muscle while the reference electrode was placed over the joint distal to the active electrode ( Figure 1). The MEP data were processed offline in MATLAB (version R2017b, MathWorks Inc., Natick, MA, USA), and only the MEPs with no preceding muscle activation and peak-to-peak amplitude greater than 50 µV were included as responses.
Brain Sci. 2020, 10, x FOR PEER REVIEW 3 of 11 The stimulation procedure was initiated by locating the optimal motor representation of the right abductor pollicis brevis (APB) muscle, i.e., APB "hotspot", with the corresponding optimized coil orientation. The hotspot was defined as the cortical site repeatedly eliciting the greatest peak-topeak MEP responses compared to adjacent stimulation sites. Once the hotspot was determined, the resting motor threshold (rMT) was identified at this cortical site using a system-integrated iterative threshold assessment tool [29]. In the RS paradigm, trials of four TMS stimuli were applied over the APB hotspot, with an inter-stimulus interval (ISI) of 1 s, at an intensity of 120% rMT. The RS paradigm, comprising 20 trials of four single biphasic TMS pulses, was employed with an inter-train interval (ITI) of 17 s [30], before (RS-baseline) and immediately (0 min), 10 min, and 20 min after the PAS intervention ( Figure 1).
We recorded MEPs via an integrated electromyography (EMG) system (Nexstim Plc) at a sampling frequency of 3 kHz. A pair of disposable Ag-Cl electrodes was utilized, with the active electrode over the belly of the APB muscle while the reference electrode was placed over the joint distal to the active electrode ( Figure 1). The MEP data were processed offline in MATLAB (version R2017b, MathWorks Inc., Natick, MA, USA), and only the MEPs with no preceding muscle activation and peak-to-peak amplitude greater than 50 µV were included as responses. The RS paradigm was divided into two components for the analysis: "dynamic RS", i.e., the ratio of the second MEP to the first one, and "stable RS", i.e., the mean of the second, third, and fourth MEPs. (c) RS applied before (baseline) and after PAS intervention (at 0 min, 10 min, and 20 min). In the PAS intervention, electrical stimulation of the median nerve-innervated APB muscle was delivered prior to TMS at an ISI of 25 ms to generate plasticity.

Paired Associative Stimulation (PAS)
PAS consisting of 180 single stimuli was applied over the right median nerve at an intensity of 300% of the sensory threshold (ST) [31]. A bipolar stimulation electrode was placed over the median nerve, and the ST was measured by adjusting the stimulation current until the subject indicated sensation of the stimulus. The pairing with TMS at the APB hotspot was implemented with a selfbuilt triggering and delayer device. To generate a plasticity effect, median nerve stimulation at a The RS paradigm was divided into two components for the analysis: "dynamic RS", i.e., the ratio of the second MEP to the first one, and "stable RS", i.e., the mean of the second, third, and fourth MEPs. (c) RS applied before (baseline) and after PAS intervention (at 0 min, 10 min, and 20 min). In the PAS intervention, electrical stimulation of the median nerve-innervated APB muscle was delivered prior to TMS at an ISI of 25 ms to generate plasticity.

Paired Associative Stimulation (PAS)
PAS consisting of 180 single stimuli was applied over the right median nerve at an intensity of 300% of the sensory threshold (ST) [31]. A bipolar stimulation electrode was placed over the median nerve, and the ST was measured by adjusting the stimulation current until the subject indicated sensation of the stimulus. The pairing with TMS at the APB hotspot was implemented with a self-built triggering and delayer device. To generate a plasticity effect, median nerve stimulation at a frequency of 0.2 Hz was delivered 25 ms prior to TMS [24], with the TMS pulses delivered at 120% of rMT. The median nerve stimulation was conducted using a constant-current electrical stimulator (Digitimer model DS7A, Digitimer, Welwyn Garden City, Herts, UK), using a rectangular pulse form (0.2 ms, maximum voltage of 300 V).

Statistical Analysis
The MEP amplitudes of each subject were first averaged based on their ordinal position in a trial, i.e., the first, second, third, and fourth. To evaluate the dynamic component of RS, the average of the second stimulus MEP amplitudes was divided by the average of the first stimulus MEP amplitudes. Further, to assess the stable component of RS, the mean of the averaged responses was computed over the second, third, and fourth stimuli per subject.
Considering the inherent heterogeneity of the neurophysiological characteristics, the analysis for identifying significant PAS-effects was initially performed at the individual level using the non-parametric Wilcoxon signed rank test. Individuals with a statistically significant increase in MEP amplitudes at a stable level at 0 min were identified as those showing LTP-like plasticity as an immediate response to PAS (as a higher MEP amplitude is considered as an index of elevated cortical excitability), and clustered as the "LTP-like group". In addition, individuals with decreased MEP amplitudes at a stable level were considered as those exhibiting LTD-like plasticity as an immediate effect to PAS and clustered as the "LTD-like group".
To test the change of the dynamic RS and stable RS over a time course of 20 min, the Friedman test was employed. Post hoc comparisons were performed using the Wilcoxon signed rank test.
A comparison of the two clusters prior to and following PAS was made using the Mann-Whitney U test. Statistical analysis was conducted using SPSS (v. 25.0, SPSS Inc., IBM Company, Armonk, NY, USA) and MATLAB (version R2017b, MathWorks Inc., Natick, MA, USA), and p < 0.05 indicated statistical significance.

Results
Eleven subjects showed no significant change of dynamic RS following PAS (p > 0.1). Only one subject exhibited significantly milder dynamic RS (a lower drop from the second MEP to the first one) (p < 0.05), and four subjects showed significantly stronger dynamic RS (a higher drop from the second to the first MEP) (p < 0.05) (Figure 2a).
Fourteen subjects exhibited a significant change at stable RS ( Figure 2b). The stable RS levels were significantly higher in six subjects immediately following PAS compared to those before PAS (p < 0.05). This heightened post-intervention MEP amplitude was assumed to be linked to LTP-like plasticity. However, eight subjects demonstrated significantly diminished MEP amplitudes at stable RS (i.e., the LTD-like group) (p < 0.05), and two subjects showed no significant change (p > 0.1). A non-parametric Friedman test revealed no change in the trend over the time of measurement following PAS (0, 10, and 20 min). One subject demonstrated delayed LTP-like plasticity at 20 min after exhibiting no effect at earlier time points.
Furthermore, a Mann-Whitney U test revealed that the dynamic and stable RS at the baseline was significantly higher in the LTD-like group compared to the LTP-like group (p < 0.05). Following PAS, no significant difference in these two components was observed between the two groups (p > 0.1) ( Figure A1).
The low between-group and high within-group homogeneities were observed at stable RS prior and following the PAS, respectively (Figure 3). , across the subjects from two clusters, within the trials before (baseline) and at 0, 10, and 20 min after the induction of short-term plasticity with PAS. Although dynamic RS was mild in the LTD-like group at baseline, a trend towards a stronger suppression and low variability at this level was observed in the LTD-like group following PAS. However, this trend in dynamic RS was towards recovery in the LTP-like group, as this component became milder with time. Similar trends, i.e., sustaining of the suppression and recovery from it, were also observed in stable RS after PAS in the LTD-like group and the LTP-like group, respectively. An asterisk indicates significant differences for pairwise comparisons between time points (p < 0.01 for *** and p < 0.05 for **). , across the subjects from two clusters, within the trials before (baseline) and at 0, 10, and 20 min after the induction of short-term plasticity with PAS. Although dynamic RS was mild in the LTD-like group at baseline, a trend towards a stronger suppression and low variability at this level was observed in the LTD-like group following PAS. However, this trend in dynamic RS was towards recovery in the LTP-like group, as this component became milder with time. Similar trends, i.e., sustaining of the suppression and recovery from it, were also observed in stable RS after PAS in the LTD-like group and the LTP-like group, respectively. An asterisk indicates significant differences for pairwise comparisons between time points (p < 0.01 for *** and p < 0.05 for **).

Discussion
Our study investigated two distinct components of RS as a measure of neuroplasticity: (1) immediate changes in motor response upon the first repetition ("dynamic RS") and (2) the suppressed level of RS ("stable RS"). Surprisingly, induction of plasticity with PAS with a 25 ms ISI resulted in different trends whereof one was rather LTD-like. Irrespective of such a discrepancy, the brain demonstrated an overall tendency towards a common level in stable RS following PAS intervention ( Figure A2).
Minimizing the surprise encountered in the face of a novel stimulus is the principle behind the free energy principle [12]. According to this principle, to maintain its integrity, any adaptive biological system, like the brain, seeks to minimize its free energy [32]. It has been proposed that minimizing the free energy rests on either changing the top-down predictions, which are the conceptual internal models, or the bottom-up predicted sensory inputs [32]. In this regard, the

Discussion
Our study investigated two distinct components of RS as a measure of neuroplasticity: (1) immediate changes in motor response upon the first repetition ("dynamic RS") and (2) the suppressed level of RS ("stable RS"). Surprisingly, induction of plasticity with PAS with a 25 ms ISI resulted in different trends whereof one was rather LTD-like. Irrespective of such a discrepancy, the brain demonstrated an overall tendency towards a common level in stable RS following PAS intervention ( Figure A2).
Minimizing the surprise encountered in the face of a novel stimulus is the principle behind the free energy principle [12]. According to this principle, to maintain its integrity, any adaptive biological system, like the brain, seeks to minimize its free energy [32]. It has been proposed that minimizing the free energy rests on either changing the top-down predictions, which are the conceptual internal models, or the bottom-up predicted sensory inputs [32]. In this regard, the dynamic RS depicts an update of the predictions in response to a twice repeated stimulus, ensuring an efficient sensory processing in a known environment. In other terms, the attenuation of the stimulus-evoked motor response upon the first repetition of the stimulus reduces the prediction errors originating from the incoming sensory information. In this respect, the suppressed level of the MEPs during the stable, suppressed part of RS may reflect a level of cortical excitability that is relatively free from ongoing dynamics in cortical excitability, which exhibits as a characteristically high intra-individual variance in the MEPs and may affect the sensitivity of MEPs to reveal longitudinal changes in excitability due to long-term neuromodulation and -plasticity.
RS has been demonstrated to last over short timescales in the visual and auditory systems, indicating a memory trace of the recently viewed or heard stimulus [33]. This short-term storage of information is reflected in our findings in the stable RS. A potential explanation for the observed stable RS might go back to the existence of a short-term internal representation of the perceived involuntary movement ("automatic memory"). Evidence consistent with this postulate is the lack of RS while an ITI of less than 3 s was employed in a TMS study, with the RS being more pronounced with longer ITIs [30]. Apart from the initial motor response, the subsequent responses elicited by TMS are modulated by sensory feedback, i.e., their magnitudes are controlled by the sensory inputs onto the motor neurons. The brain embodies a dynamic interconnected hierarchal processing organization that enables the reciprocal influence of current and past information. A plethora of positive and negative-feedback connections at both the cellular and network levels is central to sustaining the encoded sensory information on a timescale of seconds [34,35]. Hence, to maintain the automatic memory over a short timescale in RS, a negative feedback probably needs to be provided via recruiting inhibitory pathways to sustain the underlying neural activity. These pathways include the intracortical sensory areas and subcortical areas, among which the basal ganglia and thalamus play a key role. It has been demonstrated via RS that this stable state cannot be achieved in patients with progressive myoclonus type 1, who have impaired neuroplasticity in the thalamo-cortical connections [13,36].
Both dynamic and stable RS might reflect alterations in synaptic efficacy. The persistent changes in synaptic efficacy serve as a window into the formation of synaptic plastic changes, a candidate mechanism through which PAS works [26]. If the neuronal network is provided with only a positive feedback loop, that is, the spiking activity of the presynaptic neuron is correlated with the spiking activity of the postsynaptic neuron, its stability gets disrupted. In fact, this unidirectional process reduces the threshold for the presynaptic neuron to stimulate the postsynaptic neuron, thus precluding the stability and reversibility of the system. To counteract this instability and to tune the neuronal activity within a functional dynamic boundary, the brain employs an array of homeostatic mechanisms [37]. Homeostatic plasticity provides the necessary negative feedback loop to prevent the neural circuits from hyper-or hypo-activity.
A well-established proposed mechanism for homeostatic plasticity is the Bienenstock, Cooper, and Munro (BCM) model [38]. This model assumes a bidirectional synaptic plasticity, where the threshold for LTP/LTD induction varies as a function of the dynamic state of the brain. Considering this model, the more excitable the corticospinal pathway is, the more capacity for inhibition may be required. This can in part explain the reversal of the LTP-like plasticity effect to LTD-like plasticity in individuals showing higher pre-intervention MEP amplitudes (baseline). The degree of the modifications of neuronal plasticity depends on updating the synaptic efficacy. Thus, assessing the RS in the mentioned terms, i.e., are dynamic RS and stable RS, can provide us with information on how the alteration in synaptic efficacy following PAS can be reflected in RS.
A few limitations need to be acknowledged in this study. First, we applied the PAS paradigm using a fixed ISI of 25 ms. Inter-individual variability in responses has been reported for PAS due to non-optimized timing of the peripheral stimulus [39]. The potential decrease in this variability might have been achieved by employing an individualized ISI [40]. However, we did not measure individual sensory evoked potential to optimize PAS for LTP-like effects. This was by design to enable more inter-individual latency variance in the induced PAS effect, and to make the sessions shorter for the subjects. Second, in spite of having a sample size within the range of other studies in this field, the number of subjects was still small to account for generalization in large populations or in patients. We consider this a successful proof-of-concept study, but for application in patient groups, a larger-scale trial is required considering more inter-individual heterogeneities. Thirdly, we identified the PAS effect from the suppressed responses of RS (stable RS) to avoid the dynamicity of causing variance in the identification of the plastic effects, as we observed in the case of the first responses in the RS trials. This is not common practice with PAS. However, since no previous studies have been conducted with PAS in relation to RS, we had no point of reference.  Figure A1, Comparison of LTD-like group and LTP-like group prior (baseline) and following PAS (at 0, 10, and 20 min). Figure A2, RS at 0, 10, and 20 min after PAS in LTD-like and LTP-like groups.

Appendix A
Brain Sci. 2020, 10, x FOR PEER REVIEW 8 of 11 this field, the number of subjects was still small to account for generalization in large populations or in patients. We consider this a successful proof-of-concept study, but for application in patient groups, a larger-scale trial is required considering more inter-individual heterogeneities. Thirdly, we identified the PAS effect from the suppressed responses of RS (stable RS) to avoid the dynamicity of causing variance in the identification of the plastic effects, as we observed in the case of the first responses in the RS trials. This is not common practice with PAS. However, since no previous studies have been conducted with PAS in relation to RS, we had no point of reference.  Effects of LTP/LTD-like plasticity were measured with respect to "baseline" pre-intervention neural activity. Irrespective of the dynamic state of RS, the brain shows a tendency towards a suppressed Figure A1. (a) Normalized MEP amplitude (mean ± standard error) at dynamic level of RS and (b) MEP amplitude (mean ± standard error) at stable level of RS, across the subjects from the LTP/LTD-like groups within the trials before (baseline) and at 0, 10, and 20 min after the induction of short-term plasticity with PAS. Subjects in the LTD-like group exhibited significantly higher pre-intervention MEP amplitudes (baseline) in dynamic RS (p = 0.005), and stable RS (p = 0.007). No significant difference was observed between the two clusters following PAS in either component. *** indicates significant differences for pairwise comparisons (p < 0.01).
Brain Sci. 2020, 10, x FOR PEER REVIEW 9 of 11 (b) Figure A1. (a) Normalized MEP amplitude (mean ± standard error) at dynamic level of RS and (b) MEP amplitude (mean ± standard error) at stable level of RS, across the subjects from the LTP/LTDlike groups within the trials before (baseline) and at 0, 10, and 20 min after the induction of short-term plasticity with PAS. Subjects in the LTD-like group exhibited significantly higher pre-intervention MEP amplitudes (baseline) in dynamic RS (p = 0.005), and stable RS (p = 0.007). No significant difference was observed between the two clusters following PAS in either component. *** indicates significant differences for pairwise comparisons (p < 0.01). Figure A2. RS at 0, 10, and 20 min after the induction of short-term plasticity with PAS (post-PAS). Effects of LTP/LTD-like plasticity were measured with respect to "baseline" pre-intervention neural activity. Irrespective of the dynamic state of RS, the brain shows a tendency towards a suppressed Figure A2. RS at 0, 10, and 20 min after the induction of short-term plasticity with PAS (post-PAS). Effects of LTP/LTD-like plasticity were measured with respect to "baseline" pre-intervention neural activity. Irrespective of the dynamic state of RS, the brain shows a tendency towards a suppressed amplitude level of the second to the fourth MEPs within trials upon repetition of the stimulus. This static status was maintained within a narrow range in two groups after PAS.