Abstract
Learning of new skills may occur through Hebbian associative changes in the synaptic strength of cortical connections [spike-timing-dependent plasticity (STDP)], but how the precise temporal relationship of the presynaptic and postsynaptic inputs determines the STDP effects in humans is poorly understood. We used a novel paired associative stimulation protocol to repeatedly activate the short-latency connection between the posterior parietal cortex and the primary motor cortex (M1) of the left-dominant hemisphere. In different experiments, we systematically varied the temporal relationships between the stimuli and the preferential activation of different M1 neuronal populations by applying transcranial magnetic stimulation over M1 with different coil orientations and in different states of cortical excitability (rest vs muscular contraction). We found evidence for the existence of both Hebbian and anti-Hebbian STDP in human long-range connections. The induction of bidirectional long-term potentiation or depression in M1 depended not only on the relative timing between the stimuli but, crucially, on the stimulation of specific neuronal populations and the activity state of the cortex. Our findings demonstrate that these mechanisms are not fixed but susceptible to rapid adaptations. This sudden transition from anti-Hebbian to Hebbian plasticity likely involves local dynamics of interaction with different populations of postsynaptic neurons.
Introduction
The temporal proximity of novel encoding stimuli is a key element in determining the formation of a new memory. A classic example is Pavlov's dog, which learned to associate a tone with food when these stimuli were repeatedly associated. These ideas have been formalized by the Hebbian rule, proposing that neurons near-simultaneously activated increase their efficacy if the synapse consistently assists the postsynaptic target neuron to generate action potentials (APs) (Markram et al., 2011). In slice preparations, repeated coincidence of EPSPs and postsynaptic APs results in long-term potentiation (LTP) at both hippocampal and neocortical synapses (Markram et al., 1997; Sjöström et al., 2001), thus satisfying the Hebbian postulate. Typically, LTP is induced when presynaptic activity occurs just before postsynaptic spiking in the target cell (pre-pairing → post-pairing). Conversely, long-term depression (LTD) is usually induced when the postsynaptic cell fires before the presynaptic input (post-pairing → pre-pairing). This type of plasticity is referred to as spike-timing-dependent plasticity (STDP) (Magee and Johnston, 1997; Markram et al., 2011). However, these temporal rules are not irremovable. Recent findings showed that, at synapses more distal from the soma, the timing requirements for pre-pairing/post-pairing shift such that the magnitude, and eventually the sign, of synaptic modifications during STDP at distal synapses are profoundly different from that found at proximal inputs (anti-Hebbian STDP) (Letzkus et al., 2006; Sjöström and Häusser, 2006).
A main objection is that these complex mechanisms have been investigated extensively in brain slices and do not take in account spontaneous synaptic activity that runs in vivo in humans (Feldman, 2012). Although there is some indirect evidence that STDP is a key factor in governing some mechanisms of associative learning in humans (McMahon and Leopold, 2012), its existence is still matter of debate. Moreover, the experiments performed in brain slices cannot take in account the modules of large-scale connectivity that are crucial in ruling cognition (Varela et al., 2001). Here, by using bifocal transcranial magnetic stimulation (TMS) (Rizzo et al., 2009, 2011; Arai et al., 2011; Buch et al., 2011), we provide novel strong evidence supporting the existence of concomitant antithetic forms of STDP mechanisms in large-scale connections in humans. We used a novel paired associative stimulation (PAS) protocol (Stefan et al., 2000) to repeatedly activate the short-latency connection between the posterior parietal cortex (PPC) and the primary motor cortex (M1) of the left-dominant hemisphere (Koch et al., 2007) (see Fig. 1A). In different experiments, we systematically varied (1) the timing of the stimuli between PPC and M1 stimulation, (2) the current direction for M1 stimulation, and (3) the hand muscle activity. The temporal relationship between the stimuli was varied to test for the existence of Hebbian bidirectional plasticity, current direction over M1 was changed to investigate the role of different M1 neuronal populations, and muscular activity (rest vs muscular contraction) was varied to explore the role of different states of M1 cortical excitability. We show that PAS protocols can induce both Hebbian and anti-Hebbian STD-like plasticity that are fully compatible with SPTD as defined at the cellular level. These forms of STDP do not merely depend on the timing of the stimuli but indeed on the dynamics of interaction with specific cortical circuits in the postsynaptic neurons.
Materials and Methods
In total, 29 healthy volunteers (15 males and 14 females; ranging from 21 to 38 years) took part in this study (Table 1). All subjects were right-handed based on the Edinburgh Handedness Inventory (Oldfield, 1971). The experimental procedures used here were approved by the local ethics committee and were performed in accordance with the Declaration of Helsinki. Written informed consent was obtained from all recruited subjects before study initiation.
Experiment 1: PAS between the PPC and the M1 with posterior–anterior orientation.
We used bifocal TMS to repeatedly activate the connection between the PPC and M1 of the left-dominant hemisphere (Koch et al., 2007, 2008). In experiment 1A (n = 10), left PPC TMS preceded or followed the M1PA (PPC–M1PAPAS) stimulation by 5 ms. A third condition in which the interstimulus interval (ISI) varied randomly in the same temporal range was included as a control. We then tracked the time course of the aftereffects of PAS protocol by measuring motor-evoked potentials (MEPs) amplitude as an index of M1 excitability. Notably, we assumed that the inputs from PPC to M1 may discharge the same neurons that are activated by the M1 TMS itself (Koch et al., 2007).
We used a bifocal stimulation technique based on two high-power Magstim 200 machines (Fig. 1A). The magnetic stimuli had a nearly monophasic pulse configuration with a rise time of ∼100 μs, decaying back to zero over ∼0.8 ms. To measure the MEPs, electromyographic (EMG) traces were recorded from right first dorsal interosseous (FDI) muscles using 9-mm-diameter, Ag–AgCl surface cup electrodes. The active electrode was placed over the belly muscle, whereas the reference electrode was located over the metacarpophalangeal joint of the index finger. Responses were amplified using a Digitimer D360 amplifier through filters set at 20 Hz and 2 kHz with a sampling rate of 5 kHz and then recorded by a computer using SIGNAL software (Cambridge Electronic Devices). For M1 TMS, the coil was positioned over the hand motor area of left M1, defined as the point at which stimulation evoked the largest MEPs from the contralateral FDI muscle. The stimulator for M1 was connected to a small custom-made figure-of-eight-shaped coil (50 mm external diameter). The intensity of M1 TMS was adjusted to evoke an MEP of ∼1 mV peak to peak in the relaxed contralateral FDI muscles. The coil was positioned to induce a posterior–anterior (PA) directed current, being held posterolaterally at an angle of ∼45° to the midline.
To best activate the ipsilateral PPC–M1 connection, the conditioning stimulus (CS) was applied over the left PPC at an intensity of 90% of the ipsilateral resting motor threshold (RMT) (Fig. 1A). RMT was defined as the lowest intensity that evoked five small responses (∼50 μV) in the contralateral FDI muscle, in a series of 10 stimuli when the subject kept the FDI muscles relaxed in both hands (Rossini et al., 1994). The RMT was assessed at the beginning of each experimental session for each stimulating coil (PPC coil and M1 coil) separately. The magnetic pulse over PPC was applied using a second small custom-made figure-of-eight-shaped coil (50 mm external diameter). We used a neuronavigation system (Softaxic; E.M.S.) to precisely position the coil over the PPC sites, using individual T1-weighted magnetic resonance imaging volumes as anatomical reference; this technique has been described previously in detail (Koch and Rothwell, 2009; Koch et al., 2010). The stimulation points were determined before the experiment and were then marked on the adherent plastic cap worn by the subject. The individual coordinates of each stimulation site were then normalized a posteriori into the Montreal Neurological Institute (MNI) coordinate system and averaged (Fig. 1A). To target PPC, the coil was positioned over the angular gyrus, close to a posterior part of the adjoining caudal intraparietal sulcus. The coordinates of the left PPC were stored and used for neuronavigation in consecutive sessions to ensure comparable stimulation conditions among sessions. The center of the coils was positioned tangentially to the skull with the handle pointing downward and rotated medially by 15° (Fig. 1A). Indeed, we previously combined TMS with diffusion tensor imaging-based tractography to provide evidence for the anatomical bases of the PPC–M1 interaction, revealing that PPC is strongly connected with M1 by bundles of the superior longitudinal fasciculus (Koch et al., 2010, 2011).
For the PAS protocol, 100 pairs of stimuli were continuously delivered at a rate of 0.2 Hz for ∼8.3 min. The PPC stimulus preceded (PPC–M1PAPAS+5 ms) or followed (PPC–M1PAPAS−5 ms) the M1PA TMS by 5 ms (Fig. 1B). In a third control condition (PPC–M1PAPAS Random), the ISI varied randomly for each pair of stimuli in the same temporal range (ISIs: −5, −4, −3, −2, −1, +1, +2, +3, +4, +5 ms). Each subject underwent the three PAS conditions (PPC–M1PAPAS− 5 ms, PPC–M1PAPAS+5 ms, and PPC–M1PAPAS Random) in 3 different days. These sessions were performed at least 1 week apart. In each session, 20 MEPs were collected and averaged at baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS. For measuring MEPs, the coil was positioned over the M1 cortical site with the same PA orientation as for the PAS using a standard 7 cm figure-of-eight coil connected with a Magstim 200 stimulator.
In experiment 1B, we tested in eight subjects the temporal limits of the PAS interactions. The experimental procedure was the same as experiment 1A, but for the PAS protocol, the PPC preceded or followed the M1PA pulse by 50 or 20 ms (PPC preceding M1PA TMS: PPC–M1PAPAS+50 ms, PPC–M1PAPAS+20 ms; PPC after M1 TMS: PPC–M1PAPAS−50 ms, PPC–M1PAPAS−20 ms). We administered four blocks of stimulation performed in 4 different days. These sessions were performed at least 1 week apart. As in experiment 1, 20 MEPs were collected and averaged at baseline. Then, over the same hotspot, 20 MEPs were recorded at 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS and averaged.
Experiment 2: PAS between the PPC and the M1 with anterior–posterior orientation.
In experiment 2, we tested whether the effects of PAS could change depending on the preferential activation of different neuronal populations, by varying the direction of the induced current over M1 during the PAS protocol. It is known that the neuronal populations activated by TMS in M1 depend on the direction of the induced current in the brain (Day et al., 1989; Di Lazzaro et al., 2006). By simply rotating the stimulating coil over the central sulcus from a PA to an anterior–posterior (AP) direction, it is possible to obtain MEPs with different latencies. This occurs because these different coil orientations activate different set of neurons (Ni et al., 2011), likely related to the activity of the so-called I-waves. These circuits represent synchronous activity of cortico-spinal axons originating from trans-synaptic activation of cortico-spinal cells that occur with different regular delays of ∼1.5 ms (Day et al., 1989; Di Lazzaro et al., 2006) and are known to be mediated by different neuronal populations (Ni et al., 2011).
To confirm that varying the coil orientation would activate MEPs with different latencies in our subjects, we first measured onset latency of MEPs using PA, AP, and lateral–medial (LM) currents during mild contraction of the target muscle (∼10% of the maximum voluntary contraction). For each orientation (PA, AP, and LM), the hotspot was identified as the position at which the maximum and more stable MEP response was achieved. The hotspot was redefined for each coil orientation. The intensity of stimulation was adjusted relatively on the active motor threshold (AMT), defined as the lowest intensity to evoke an MEP of 200 μV in more than 5 of 10 consecutive trials while subjects maintained ∼10% contraction of the target muscle (Rothwell, 1997). AMT was assessed for each coil orientation. For PA and AP, stimulus intensity was set at 110% AMT; for LM, it was set at 150% AMT (Hamada et al., 2012). MEP latencies with different coil orientations were computed by visual inspection in each trial.
We then applied the PAS procedure as in experiment 1A, but the M1 coil was positioned to induce current in the AP direction (PPC–M1APPAS). M1AP TMS was applied over the hotspot of the FDI by rotating the coil handle of 180° from the PA direction. The intensity of M1 TMS was adjusted to evoke an MEP of ∼1 mV peak to peak in the relaxed contralateral FDI muscles. The CS was applied over the left PPC at an intensity of 90% of the ipsilateral RMT. The PPC stimulus preceded (PPC–M1APPAS+5 ms) or followed (PPC–M1APPAS−5 ms) the M1AP TMS by 5 ms. Each subject underwent the two PAS conditions in 2 different days, performed at least 1 week apart. Twenty MEPs were collected and averaged at baseline and at 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS and averaged.
In the experiment 2A (n = 10), pre-PAS and post-PAS MEPs measures were collected over M1 with PA current orientation. Then we performed a control experiment [experiment 2B (n = 7)] in which MEPs were collected with both PA and AP current orientation. The PPC stimulus preceded (PPC–M1APPAS+5 ms) the M1AP TMS by 5 ms.
In experiments 2A and 2B, the coil positioned over M1 was oriented to induce current in AP direction. This orientation requires higher intensity to be applied to stimulate M1 if compared with the intensity needed with the PA current direction (experiment 1). Such a difference could complicate the interpretation of the results obtained with the two coil orientations because of differences in the magnitude and spatial extent of M1 stimulation rather than differential effects on specific intracortical I-wave circuits.
Thus, we performed another control experiment (experiment 2C; n = 6) in which we increased the intensity of M1 stimulation for the PPC–M1PAPAS protocol to make the two experiments comparable in terms of intensity of M1 TMS, regardless of the coil orientation used (AP or PA) (Table 1). We decided to increase M1 TMS in the PPC–M1PAPAS protocol instead of lowering M1 TMS for the PPC–M1AP PAS given that, in this latter case, we would have caused a subthreshold activation of M1 AP circuits and this could have been an additional confound for the comparison between the two experiments. The PPC stimulus preceded (PPC–M1PAPAS+5 ms) the M1PA TMS by 5 ms.
Experiment 3: state dependency of PAS.
In experiment 3A, performed in a sample of eight subjects, we tested whether the effects of PPC–M1PAPAS protocol could change depending on the state of M1 excitability. Hence PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5 ms protocols were applied with subjects performing a mild muscular contraction of the index finger and the thumb finger of the right hand at 10% of maximal voluntary contraction throughout the entire duration of the PAS protocol. The level of muscular contraction was visually monitored by the experimenter on the basis of the EMG amplitude. M1PA TMS was applied over the hotspot of the FDI with intensity sufficient to evoke a 1 mV MEP as in the experiment 1. Twenty MEPs were collected and averaged at baseline. Then, 20 MEPs were recorded over the same hotspots at 0–5, 6–10, 11–15, 16–20, and 21–25 min after each PAS intervention and averaged.
In experiment 3A, the intensity of M1 TMS was substantially reduced compared with experiment 1 (Table 1). Therefore, we performed a control experiment in six subjects (experiment 3B) in which PPC–M1PAPAS was applied during a mild muscular contraction as in experiment 3A, but the intensity of M1 TMS was increased to be matched with that used during PPC–M1PAPAS at rest (experiment 1). The PPC stimulus preceded (PPC–M1PAPAS+5 ms) the M1PA TMS by 5 ms.
Experiment 4: topographical M1 specificity of PAS effects.
In this experiment, performed in a sample of eight subjects, we tested the effects of PPC–M1PAPAS protocols on different hand muscles to investigate whether the effects of the PAS were specific for the stimulated portion of cortical area within M1. MEPs were recorded from surface electrodes over the contralateral abductor digiti minimi (ADM) and FDI muscles on the right hand. We first determined the location of the hotspot for evoking MEPs in the FDI and ADM muscles and marked these points on the scalp. The PPC–M1PAPAS+5 ms protocol was applied, with the M1PA TMS delivered over the hotspot of the FDI as in experiment 1. Twenty MEPs were collected and averaged at baseline for each muscle. Then, over the same hotspots, 20 MEPs were recorded at 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS and averaged.
Experiment 5: effects of PAS on the intracortical circuits of M1.
In this experiment, performed in a sample of nine subjects, we tested the effects of PPC–M1PAPAS protocols (PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5 ms) on the short intracortical inhibition (SICI) and intracortical facilitation (ICF) circuits (experiment 5A). We used a 7 cm figure-of-eight coil connected with two Magstim 200 stimulators to apply paired TMS over the left motor cortex. The coil was placed at the optimal position for eliciting MEPs from the left FDI muscle. We evaluated the changes occurring in intracortical circuits testing both SICI at 2 and 3 ms ISIs and ICF at 15 ms ISI (Kujirai et al., 1993) at baseline and after 0–5 min after PAS. CS was set at 80% AMT while the intensity of the test stimulus (TS) was adjusted to evoke an MEP of ∼1 mV peak to peak in the relaxed FDI. The amplitude of the conditioned MEP at each ISI was expressed as a percentage of the mean peak-to-peak amplitude size of the unconditioned TS in that block. For the PAS protocol, the conditions PPC–M1PAPAS−5 ms and PPC–M1PAPAS+5 ms were administered in 2 different days at least 1 week apart.
In experiment 5B, performed in a sample of seven subjects, we tested the effects of PPC–M1PAPAS protocols on the short interval intracortical facilitation (SICF) and the short latency afferent inhibition (SLAI) circuits of M1. We evaluated the changes occurring in SICF and SLAI in two blocks performed at baseline and after PAS. The order of presentation of the blocks was pseudorandomized across subjects. Recordings started 1 min for the first block and 10 min after PAS for the second block. For the PAS protocol, the conditions PPC–M1PAPAS−5 ms and PPC–M1PAPAS+5 ms were administered in 2 different days at least 1 week apart.
For SICF circuit, the TS given alone or the TS followed by the CS at various ISIs were intermixed randomly in one block (Ziemann et al., 1998a). In each block, seven conditions were randomly intermingled: TS alone (MEP) and CS + TS (conditioned MEP for each six different ISIs: 1, 1.3, 2.1, 2.5, 3.3, and 4.1 ms). Before each block, the intensity of TS was adjusted to evoke an MEP of ∼1 mV peak to peak in the relaxed left FDI. CS intensity was set at 90% RMT. Ten responses were collected for paired conditioned MEP for each ISI and for TS alone with a total number of 70 trials in each block. The intertrial interval was set at 5 s (±10%), for a total duration of ∼6 min. Measurements were made on each individual trial. The mean peak-to-peak amplitude of the conditioned MEP at each ISI was expressed as a percentage of the mean peak-to-peak amplitude size of the unconditioned test pulse in that block. SLAI measures were studied and recorded using the technique that has been described previously (Sailer et al., 2003). CSs consisted of single pulses (200 μs) of electrical stimulation applied to the right median nerve at the wrist through bipolar electrodes (proximal cathode). The intensity of the CS was set at intensity just above motor threshold to evoke a visible twitch of the thenar muscles. The intensity of the cortical magnetic TS was adjusted to evoke a muscle response of ∼1 mV peak-to-peak amplitude in the relaxed right FDI. The CS delivered to the peripheral nerve preceded the magnetic TS by different ISIs (16, 20, 24, and 28 ms). Ten responses were collected for each ISI and for TS alone with a total number of 50 trials in each block. The intertrial interval was set to 5 s (±10%), for a total duration of ∼5 min. Measurements were collected for each individual trial. The mean peak-to-peak amplitude of the conditioned MEP at each ISI was expressed as a percentage of the mean peak-to-peak amplitude of the unconditioned single-pulse magnetic stimulus in that block.
Data analyses.
Data were analyzed using SPSS for Windows version 11.0 on the mean MEPs amplitude in each condition. In experiment 1A, a two-way repeated-measures ANOVA was performed with protocol (PPC–M1PAPAS+5 ms, PPC–M1PAPAS−5 ms, and PPC–M1PAPAS Random) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factors. In experiment 1B, a two-way repeated-measures ANOVA was performed with protocol (PPC–M1PAPAS+50 ms, PPC–M1PAPAS+20 ms, PPC–M1PAPAS−20 ms, and PPC–M1PAPAS−50ms) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factors.
In experiment 2A, a two-way repeated-measures ANOVA was performed with protocol (PPC–M1PAPAS+5 ms, PPC–M1PAPAS−5 ms, PPC–M1AP+5 ms, PPC–M1AP−5 ms) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factors. In experiment 2B, a two-way repeated-measures ANOVA was performed with coil orientation (PA and AP) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factors. In experiment 2C, a one-way repeated-measures ANOVA was performed with time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factor.
In experiment 3, a two-way repeated-measures ANOVA was performed with protocol (PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5 ms) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factors. In experiment 3B, a one-way repeated-measure ANOVA was performed with time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PAS) as within-subjects main factor. In experiment 4, a two-way repeated-measure ANOVA was performed with hand muscle (FDI and ADM) and time (baseline and 0–5, 6–10, 11–15, 16–20, and 21–25 min after PPC–M1PAPAS−5 ms) as within-subjects main factors. In experiment 5A, a three-way repeated-measures ANOVA was performed with protocol (PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5 ms), time (baseline and 0–5 min after PAS), and ISI (2, 3, and 15 ms) as within-subjects main factors. In experiment 5B, two different three-way repeated-measures ANOVA were performed: one for SLAI protocol with protocol (PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5ms), time (baseline and after PAS), and ISI (16, 20, 24, and 28 ms) as within-subjects main factors for and another for SICF protocol with protocol (PPC–M1PAPAS+5 ms and PPC–M1PAPAS−5 ms), time (baseline and after PAS), and ISI (1.0, 1.3, 2.1, 2.5, 3.3, and 4.1 ms) as within-subjects main factors.
The Greenhouse–Geisser correction was used for nonspherical data. When a significant main effect was reached, paired t tests with Bonferroni's correction were used to characterize the different effects of the specific ISIs. Mauchley's test examined for sphericity. For all statistical analyses, a p value <0.05 was considered to be significant.
Results
In experiment 1A, 100 stimuli of PPC stimulation were paired with M1 TMS with PA orientation (PPC–M1PA). We found evidence that this PAS protocol induces an anti-Hebbian bidirectional cortico-cortical associative plasticity in the PPC–M1PA network as revealed by the ANOVA showing a main effect of protocol (F(2,18) = 9.12; p = 0.0018) as well as a protocol × time interaction (F(10,90) = 3.26; p = 0.0012). When PPC precedes M1PA stimulation (PPC–M1PAPAS+5 ms), there is a long-lasting decrease of the excitability of M1, indicating an LTD-like effect; conversely, when PPC follows M1PA stimulation (PPC–M1PAPAS−5 ms), there is a long-lasting increase of the excitability of M1, indicating an LTP-like effect. Post hoc analyses showed that for PPC–M1PAPAS+5 ms there is an inhibition of MEP amplitude compared with baseline at 11–15 min (p = 0.010) and 16–20 min (p = 0.009); for PPC–M1PAPAS−5ms, there is a facilitation of MEP amplitude compared with baseline at 0–5 min (p = 0.010), 6–10 min (p = 0.009), 16–20 min (p = 0.011), and 21–25 min (p = 0.007). No significant change is observed when the delay between the two stimuli varies randomly between −5 and +5 ms (Fig. 2A).
The same effects are still evident when PPC precedes or follows M1PA by 20 ms but not by 50 ms, as revealed by experiment 1B (ANOVA showing a main effect of protocol (F(3,21) = 8.30; p = 0.0007), time (F(5,35) = 4.45; p = 0.0030), as well as a protocol × time interaction (F(15,105) = 3.15; p = 0.0002). Post hoc analyses showed that for PPC–M1PAPAS+20 ms there is an inhibition of MEP amplitude compared with baseline at 6–10 min (p = 0.03) and 11–15 min (p = 0.012); for PAS−20 ms, there is a facilitation of MEP amplitude compared with baseline at 0–5 min (p = 0.003), 6–10 min (p = 0.010), 11–15 min (p = 0.008), 16–20 min (p = 0.012), and 21–25 min (p = 0.011) (Fig. 2B). Therefore, experiments 1A and 1B demonstrate that, in the PPC–M1PA network, LTD is obtained when the presynaptic input precedes the postsynaptic activation (positive paring), whereas LTP is evoked in the opposite condition (negative paring), fully satisfying the rules of anti-Hebbian STDP (Fig. 2C). For each experimental condition (PPC–M1PAPAS−5 ms, PPC–M1PAPAS−20 ms, PPC–M1PAPAS+5 ms, and PPC–M1PAPAS+20 ms), we also analyzed eventual changes occurring during the plasticity-inducing protocols. In these cases, there were no significant differences in the amplitudes of the MEP conditioned by PPC stimulation during the administration of 100 paired stimulations as shown by ANOVA reporting no effects for both the condition (F(3,21) = 1.09; p = 0.37) and time (F(9,63) = 0.76; p = 0.64) main factors and for the condition × time interaction (F(27,189) = 0.88; p = 0.63) (Fig. 2D).
In experiment 2A, we tested whether the effects of PAS could change depending on the preferential activation of different groups of interneurons. Hence, during PAS administration, we varied the orientation of the stimulating coil over M1 to reverse the direction of induced current from PA to AP direction (PPC–M1AP PAS). To verify that this procedure was able to activate different neuronal populations in our sample, we first measured MEP latency at baseline, before any PAS protocol. As expected, we found that different current orientations result in MEPs with different latencies (Fig. 3A). The MEP latencies were 19.5 ± 1.3 ms for LM, 20.5 ± 2.1 ms for PA, and 23.4 ± 2.4 ms for AP current direction (one-way ANOVA main effect, F(2,18) = 50.77; p < 0.00001; comparisons between LM and PA, p < 0.01; between PA and AP, p < 0.001; between LM and AP, p < 0.001). Notably, the mean conditioned MEP amplitude during cortico-cortical PAS in this experiment using AP orientation did not differ from that used for PPC–M1 cortico-cortical PAS with PA orientation in experiment 1A. We then found that the PPC–M1APPAS protocol induced opposite effects as those obtained in experiment 1A, as revealed by the ANOVA showing main effects of protocol (F(1,9) = 6.81; p = 0.02) and time (F(5,45) = 3.72; p = 0.006), as well as a protocol × time interaction (F(5,45) = 3.59; p = 0.008). When PPC precedes M1AP stimulation (PPC–M1APPAS+5 ms), there is a long-lasting increase of the excitability of M1, indicating an LTP-like effect; conversely, when PPC follows M1AP stimulation (PPC–M1APPAS−5 ms), there is a trend toward a long-lasting decrease of the excitability of M1. Post hoc analyses showed that, for PPC–M1APPAS+5 ms, there is a facilitation of MEPs amplitude compared with baseline at 0–5 min (p = 0.007), 6–10 min (p = 0.009), 11–15 min (p = 0.010), and 16–20 min (p = 0.009); for PPC–M1APPAS−5 ms, there is a trend toward a significant inhibition of MEP amplitude compared with baseline at 11–15 min (p = 0.06) (Fig. 3B). Therefore, when we activated different M1 neuronal populations, by changing the direction of the induced current in M1 (PPC–M1APPAS), we observed an opposite temporal relationship: LTP was evoked when the presynaptic input preceded the postsynaptic activation, whereas for negative timing, there was a trend toward LTD. Hence, STDP followed the temporal rules of classical Hebbian STDP.
In these experiments, although PAS was performed with AP orientation, MEPs were evaluated with the standard PA orientation as in experiments 1A and 1B. Therefore, in a control experiment (experiment 2B), we verified whether the PPC–M1APPAS effects would differ if MEPs were tested with different orientations. We found that PPC–M1APPAS+5 ms protocol induces similar LTP-like effects in both the AP and PA MEPs conditions as revealed by the ANOVA showing main effects of coil orientation (F(1,6) = 12.20; p = 0.017) and time (F(5,25) = 4.76; p = 0.003) but not a coil orientation × time interaction (F(5,25) = 1.14; p = 0.36). Post hoc analyses revealed that the LTP-like effects tested with PA MEPs were more durable, showing a facilitation of PA MEP amplitude compared with baseline at 0–5 min (p = 0.007), 6–10 min (p = 0.009), 11–15 min (p = 0.010), and 16–20 min (p = 0.009). Conversely, the LTP-like effects tested with AP MEPs were shorter, being significant only at 0–5 min (p = 0.04) and 6–10 min (p = 0.006) (Fig. 3C).
We also analyzed eventual changes occurring during the plasticity-inducing protocols. In these cases, there were no significant differences in the amplitudes of the MEP conditioned by PPC stimulation during the administration of 100 paired stimulations (Fig. 3D).
Experiment 2C showed that the different Hebbian effects obtained in experiments 1 and 2 by changing the orientation of M1 TMS during the PAS protocol cannot be ascribed to differences in the M1 intensity of stimulation. When the M1 TMS in the PPC–M1PAPAS+5 ms protocol was increased to be comparable with that used in experiment 2A, we found the expected LTD-like effect as in experiment 1 (ANOVA showing main effects of time, F(5,25) = 2.72; p = 0.043) (Table 2).
Experiment 3A revealed that the direction of cortical plasticity is also strictly dependent on the M1 state of excitability. When subjects were asked to keep a slight muscular contraction during the PPC–M1PAPAS protocol, the effects described in experiment 1 with subjects at rest were fully reversed. The ANOVA showed a main effect of protocol (F(1,7) = 11.79; p = 0.018) and a protocol × time interaction (F(5,25) = 7.51; p = 0.002). Although with subjects at rest PAS+5 ms induces LTD-like STDP, during contraction cortical plasticity turns toward LTP, even applying the same PPC–M1PAPAS protocol. Post hoc analysis showed that, for PPC–M1PAPAS+5 ms, there is a facilitation of MEP amplitude, instead of the expected inhibition observed at rest, at 0–5 min (p = 0.010), 6–10 min (p = 0.011), and 16–20 min (p = 0.008) compared with baseline (paired t tests). Similarly, for PAS−5 ms, there is an inhibition of MEP amplitude instead of the expected LTP found at rest at 0–5 min (p = 0.004), 6–10 min (p = 0.002), 11–15 min (p = 0.003), and 16–20 min (p = 0.0007) compared with baseline (Fig. 4). Thus, a transition from anti-Hebbian to classical Hebbian rules is found when the PPC–M1PAPAS protocols are applied with subjects performing a slight muscular contraction.
Experiment 3B showed that the different Hebbian effects obtained in experiments 1 and 3 by changing the state of muscular activity during the PAS protocol cannot be ascribed to differences in the M1 intensity of stimulation. When the M1 TMS in the PPC–M1PAPAS+5 ms protocol with muscular contraction was increased to be comparable with that used in experiment 1, we found the same LTP-like aftereffects as in experiment 3A (ANOVA showing main effects of time, F(5,25) = 2.67; p = 0.045) (Table 2).
Experiment 4 showed that the aftereffects of PAS are highly specific to the M1 cortical area that is paired with PPC stimulation. In fact, although the LTP effect is observed when PPC–M1PAPAS occurs between the PPC and the cortical area of the FDI that is stimulated during the PAS protocol, no change is observed for the motor cortical representation of another close hand muscle in M1, the ADM [ANOVA showing a main effect of muscle (F(1,7) = 17.70; p = 0.008), time (F(5,25) = 3.65; p = 0.012), as well as a muscle × time interaction (F(5,25) = 4.78; p = 0.003). Also in this experiment, for PPC–M1PAPAS−5 ms, there is a facilitation of MEP amplitude compared with baseline at 0–5 min (p = 0.007), 6–10 min (p = 0.011), 11–15 min (p = 0.009), and 16–20 min (p = 0.012) (Fig. 5).
Experiment 5A performed with PPC–M1PAPAS demonstrated the plastic aftereffects are associated with specific changes in M1 intracortical circuits [ANOVA showing a main effect of ISI (F(2,16) = 102.96; p = 0.000001) as well as a protocol × time interaction (F(1,8) = 15.32; p = 0.004) and a protocol × time × ISI interaction (F(2,16) = 3.32; p = 0.041)]. Namely, PPC–M1PAPAS+5 ms induces a decrease of ICF at 15 ms (p = 0.002), whereas the opposite effect is observed for the PPC–M1PAPAS−5 ms protocol (p = 0.031), suggesting that long-lasting changes in M1 likely involve the specific set of interneurons that mediate ICF (Fig. 6A,B).
Experiment 5B showed that SLAI and SICF circuits were not modified by the PAS protocols tested. For SLAI measurements ANOVA showed only a main effect of ISI (F(3,18) = 15.13; p = 0.00003) but not any other significant main factor or interaction (Fig. 7A,B). For SICF measurements, ANOVA showed only a main effect of ISI (F(5,25) = 7.82; p = 0.0001) but not any other significant main factor or interaction (Fig. 7C,D).
Discussion
We provide here novel strong evidence revealing the existence of both anti-Hebbian and Hebbian STDP in human long-range connections. Our findings reinforce the notion that mechanisms of Hebbian-like plasticity can be an important property of cortico-cortical connections (Rizzo et al., 2009; Arai et al., 2011; Buch et al., 2011) and demonstrate that these mechanisms are not fixed but susceptible to rapid adaptations. As discussed below, this sudden transition from anti-Hebbian to Hebbian plasticity likely involves local dynamics of interaction with different populations of postsynaptic neurons.
Our findings reinforce that the opposite Hebbian and anti-Hebbian STDP effects described here in humans for the first time resemble some recent findings obtained in animal models, describing a switch from Hebbian to anti-Hebbian STDP as the distance of the synaptic contact from the soma increases (Froemke et al., 2005; Sjöström and Häusser, 2006). It was shown that the temporal window of STDP is different at proximal and distal synapses, leading to different degrees of potentiation or depression (Kampa et al., 2007) and that the site of dendritic interaction with the afferent input is a critical element that locally controls the direction of plastic changes within the target cortical neuron (Sjöström and Häusser, 2006). Depending on the relative distance between the site of synaptic interaction with dendritic arbor and the soma of the pyramidal cell, striking different LTP or LTD effects can be induced even when the presynaptic and postsynaptic stimuli are delivered with an identical temporal order. If the presynaptic input precedes the backpropagating APs (bAPs) by few milliseconds (positive timing) and targets dendrites that are localized close to the soma in layer 5, this scheme induces LTP; conversely, if the same experimental procedure is performed targeting the more distal dendrites that are localized in layer 2/3, an opposite LTD phenomenon emerges (Sjöström and Häusser, 2006). Similarly, when the presynaptic input follows the postsynaptic bAP by few milliseconds (negative timing), there is LTP for distal but LTD for proximal inputs (Froemke et al., 2005; Kampa et al., 2007). Notably, the sign of distal plasticity switches from negative (LTD) to positive (LTP) once a sufficient number of distal inputs cooperates, when more inputs are recruited by increasing the stimulation strength (Sjöström and Häusser, 2006).
In our study, the opposite PAS aftereffects observed by varying the orientation of the stimulating coil likely involve the interactions with different M1 neuronal populations. In the M1, TMS with AP and PA orientations is known to induce opposite current flows that activate a set of neurons likely located in different cortical sites (Ni et al., 2011). These opposite AP and PA current flows are supposed to induce upward- or downward-directed current along the dendritic tree of the pyramidal neurons in the anterior bank of the central sulcus (Sommer et al., 2012). It has been suggested that the AP orientation would favor the activation of more apical neuronal populations and be more prone to induce stronger bAPs. In contrast, the PA orientation would induce a current flow directed downward from layer 1 to layer 6 that would be more depolarizing at deeper sites (Sommer et al., 2012). This background could explain why PPC–M1APPAS induces STDP following the temporal rules of classical STDP. AP TMS would more promptly induce stronger bAPs, favoring LTP at more apical layers, in which the PPC inputs arrive (Koch et al., 2007). Conversely, anti-Hebbian STDP would emerge for the PPC–M1PAPAS because PA TMS provokes weaker distal depolarization, thereby inducing more easily LTD.
It is also possible that these distinct effects could also involve the activity of I-wave interneurons that could mediate the distal afferent input from PPC. According to the canonical model of cortical circuitry, early I-waves are likely to reflect monosynaptic input to cortico-spinal neurons from proximal interneurons, whereas more complex oligo-synaptic circuits, including inhibitory interneurons, are involved in late I-wave generation (Day et al., 1989; Sakai et al., 1997; Di Lazzaro et al., 2006, 2012). These interneurons could also have contributed locally to the generation of different forms of STDP. However, the activation of early or late I-wave circuits could have simply caused a temporal shift in the reversal point of the learning rule rather than flipping the sign of the temporal learning rule. Moreover, we did not observe any change in the SICF circuits after PPC–M1PAPAS protocols. Therefore, the possible contribution of the I-wave circuits remains controversial. Another possibility is that AP and PA activated circuits could have involved different synaptic relation with GABAA interneurons that have also be invoked in mediating SPTD in humans (Elahi et al., 2012). Therefore, changes in the synaptic efficacy between different neuronal species (pyramidal vs inhibitory neurons) could also produce such effects.
Fascinatingly, the STDP effects were also modulated when the PPC–M1PAPAS protocols were applied with subjects performing a slight muscular contraction. Also in this case, the direction of induced plasticity flipped completely, now following the timing rules of standard Hebbian STDP. It may well be that, by changing the activation state of M1 during muscular contraction, the same circuits activated during PPC–M1APPAS could become more engaged, leading to a stronger distal dendritic depolarization (Sjöström and Häusser, 2006). The discharge of pyramidal cortical neurons in the course of muscular contraction could increase the coincidence between distally evoked EPSPs with APs satisfying mechanisms of cooperativity among inputs. Interestingly, this state dependency property of cortical plasticity resembles also recent works testing the PAS interactions between other cortical areas with M1 (Arai et al., 2011; Buch et al., 2011) and other protocols, such as theta burst stimulation (Huang et al., 2005; Gentner et al., 2008).
In this regard, it has been suggested that the transition from LTD toward LTP (and vice versa) in these experiments involves the modulation of Ca2+ dynamics at the level of the dendritic synaptic contact. Interestingly, Ca2+ dynamics regulate both the polarity changes of naive synapses exposed to a plasticity-inducing stimulation protocol and the physiological events induced by isometric contractions (Wankerl et al., 2010).
Whereas LTP-like plasticity was associated with increased ICF activity, LTD was related with an opposite effect. Although ICF has been related to glutamatergic activity (Ziemann et al., 1998b), the mechanisms of ICF are less clear than those related to SICI. It has been proposed that ICF might result from the recruitment of circuits related to the activation of long-range connections originating from remote areas (Ziemann et al., 1998b). Given that the current PAS protocol is based on the activation of such long-range connections, it is conceivable that the ICF changes could involve the modulation of the same pathways activated by PPC–M1 PAS. Indeed, the specificity of the LTP/LTD-like effect observed in our study was further highlighted by the findings that PPC–M1 PAS induced critical changes only at precise cortical regions within M1. Only MEPs evoked from the motor hotspot stimulated by PAS (recorded from the FDI) were increased after the protocol, whereas no effect at all was found from a close cortical spot in which the excitability of the ADM was tested. These findings resemble the input specificity described in classical PAS protocols. However, the input from PPC likely reaches both the FDI and ADM cortical spots, and we were not able to control for other cortical hotspots that do not receive the cortical projection from PPC.
We did not directly control for any possible spinal effect induced by PPC TMS during the PAS protocols. However, we demonstrated previously that PPC TMS at 90% RMT does not change the amplitude of the H reflex (Koch et al., 2007), making unlikely this possibility. We also found that the strength of PPC–M1 interactions did not vary in different protocols. Although this might be surprising, this lack of change in the strength of the stimulated conditioning pathway is consistent with classical PAS in which SLAI representing motor cortical inhibition mediated by an ascending central cholinergic projection (Di Lazzaro et al., 2000) does not change (Stefan et al., 2000; Hamada et al., 2012).
The current results obtained with PPC–M1PAPAS at rest have an opposite timing to that characterizing the classical median nerve PAS interactions (Stefan et al., 2000). One possibility is that cortico-cortical inputs (mediating the current PPC–M1 interactions) could terminate in different layers than thalamic ones (which mediate the median nerve PAS effects). In fact, thalamo-cortical neurons synapse on the dendrites of the pyramidal neurons in several cortical layers (McFarland and Haber, 2002). Apart from the more widespread, diffuse terminals to layer 1, there are also strong topographically specific thalamo-cortical projections to deeper layers L2/3 and L5 (McFarland and Haber, 2002).
An intrinsic critique of the current study is that we did not have a direct recording of APs at the postsynaptic neurons, and therefore the current data could be considered as an indirect demonstration of the underlying STDP mechanisms. However, MEPs can be considered as a global measure of cortical excitability, which reflects the summation of thousands of APs. Indeed, the striking analogies between our current findings with the extensive amount of data obtained in animal slices make us very confident in claiming that what we observed were STPD processes occurring at a system level. In particular, our data are compatible with previous animal studies using paired recordings in dissociated neuronal cultures, reporting an ∼40-ms-long coincidence window, with an astoundingly rapid 1 ms transition between LTP and LTD for near-perfect coincidence between presynaptic and postsynaptic cell activity (Bi and Poo, 1998).
In conclusion, we demonstrate that the time-locked activation of parallel distinct human cortico-cortical connections may lead to antithetic forms of STDP. Additional studies would clarify in the future the behavioral role played by the different forms of Hebbian plasticity described here. For instance, it is possible that these antithetic cortical interactions could play a specific role once new motor memories are encoded or when unnecessary information has to be erased. This interplay could be relevant during processes of adaptation and de-adaptation in visuomotor tasks (Tanaka et al., 2009). From a clinical perspective, the current data could be important to investigate plasticity of specific damaged cortical circuits in patients with motor disorders or to promote novel strategies for stroke-related cortical reorganization (Ridding and Rothwell, 2007; Koch et al., 2012).
Footnotes
This work was supported by Italian Ministry of Health Grants RF08.18, MR08.7, and 09.219 (G.K.). We are grateful to Dr. Silvia Picazio and Dr. Sonia Bonnì for their help in data collection. Many thanks to Prof. John Rothwell, Dr. Masashi Hamada, Dr. Binith Cheeran, and Prof. Stefano Ferraina for their constructive comments on a previous version of this manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Giacomo Koch, Laboratory of Clinical and Behavioral Neurology, Santa Lucia Foundation, Institute for Inpatient Treatment and Scientific Studies, Via Ardeatina, 306, I-00179 Rome, Italy. g.koch{at}hsantalucia.it