Next Article in Journal
The Burden of Respiratory Alterations during Sleep on Comorbidities in Obstructive Sleep Apnoea (OSA)
Previous Article in Journal
Update on Non-Pharmacological Interventions for Treatment of Post-Traumatic Headache
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Coil Orientation on the Stimulation of the Pre–Supplementary Motor Area: A Combined TMS and EEG Study

1
Non−Invasive Brain Stimulation Unit, IRCCS Santa Lucia Foundation, 00179 Rome, Italy
2
Department of Clinical and Movement Neurosciences, UCL Queen Square Institute of Neurology, University College London, London WC1N 3BG, UK
3
Department of Human Neurosciences, Sapienza University of Rome, 00185 Rome, Italy
4
IRCCS Neuromed, 86077 Pozzilli, Italy
5
Department of Bioengineering, Imperial College, London SW7 2AZ, UK
6
BSICoS Group, I3A Institute, University of Zaragoza, IIS Aragón, 50009 Zaragoza, Spain
7
Neurology Unit, Department of Clinical and Experimental Sciences, University of Brescia, 25121 Brescia, Italy
8
The Royal’s Institute of Mental Health Research, University of Ottawa, Ottawa, ON K1N 6N5, Canada
9
Département de Psychoéducation et Psychologie, Université du Québec en Outaouais, Gatineau, QC J8X 3X7, Canada
10
Department of Cellular and Molecular Medicine, Faculty of Medicine, University of Ottawa, Ottawa, ON K1N 6N5, Canada
11
Department of Medical Sciences and Public Health, University of Cagliari, 09124 Cagliari, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to the present work.
Brain Sci. 2022, 12(10), 1358; https://doi.org/10.3390/brainsci12101358
Submission received: 28 August 2022 / Revised: 26 September 2022 / Accepted: 3 October 2022 / Published: 6 October 2022
(This article belongs to the Section Neurotechnology and Neuroimaging)

Abstract

:
Studies using transcranial magnetic stimulation (TMS) have demonstrated the importance of direction and intensity of the applied current when the primary motor cortex (M1) is targeted. By varying these, it is possible to stimulate different subsets of neural elements, as demonstrated by modulation of motor evoked potentials (MEPs) and motor behaviour. The latter involves premotor areas as well, and among them, the presupplementary motor area (pre–SMA) has recently received significant attention in the study of motor inhibition. It is possible that, similar to M1, different neuronal populations can be activated by varying the direction and intensity of TMS; however, the absence of a direct electrophysiological outcome has limited this investigation. The problem can be solved by quantifying direct cortical responses by means of combined TMS and electroencephalography (TMS–EEG). We investigated the effect of variable coil orientations (0°, 90°, 180° and 270°) and stimulation intensities (100%, 120% and 140% of resting motor threshold) on local mean field potential (LMFP), transcranial evoked potential (TEP) peaks and TMS–related spectral perturbation (TRSP) from pre–SMA stimulation. As a result, early and late LMFP and peaks were larger, with the coil handle pointing posteriorly (0°) and laterally (90°). This was true also for TRSP in the β–γ range, but, surprisingly, θ–α TRSP was larger with the coil pointing at 180°. A 90° orientation activated the right M1, as shown by MEPs elicitation, thus limiting the spatial specificity of the stimulation. These results suggest that coil orientation and stimulation intensity are critical when stimulating the pre–SMA.

1. Introduction

Transcranial magnetic stimulation (TMS) is a noninvasive brain stimulation technique that has been extensively used in the past decades to study the basic physiology of the cerebral cortex, as well as human behaviour in health and disease. Most studies so far have been performed on the primary motor area (M1); a single TMS pulse applied here activates excitatory inputs to corticospinal neurons. This ultimately leads to the generation of descending volleys, which give rise to a compound muscle action potential called motor evoked potential (MEP) [1]. The current direction induced in the brain by the TMS pulse has a critical role in this process. It is well known that the latency of the MEP is longer if the current applied has an anterior to posterior (AP) direction compared with posterior to anterior (PA). Possibly, this is due to the AP and PA currents activating different sets of inputs to corticospinal neurons [2,3,4,5,6,7]. Besides coil orientation (CO), this differential activation also depends on the stimulation intensity (SI) used, i.e., a low SI is more selective in activating subsets of inputs to corticospinal neurons [5]. Even more importantly, this specific activation seems to be behaviourally relevant: it has been suggested that the suppression of a set of motor cortical neurons, obtained by a specific current direction of TMS, may have an impact on movement preparation [8,9].
Motor planning is a complex phenomenon also involving higher–order motor areas, including the presupplementary motor area (pre–SMA). In humans, it is located in the dorsomedial frontal cortex, anterior to the leg representation of the primary motor cortex (M1) [10,11]. Several lines of evidence have pointed towards a role of the pre–SMA in complex motor behaviour, including self–initiated activity, generation of action sequences and motor learning [12]. This has been confirmed in human studies using TMS; in particular, the human pre–SMA has received attention in the context of switching and stopping behaviour [13,14,15,16,17]. However, in most studies, parameters such as CO and SI have been selected without a clear rationale. It is possible that, as for M1, different neuronal populations can be activated in the pre–SMA by varying CO and SI; however, investigating this outside M1 is difficult due to the absence of a direct readout such as MEP. This problem can be solved with the use of combined TMS and electroencephalography (TMS–EEG), which allows to record postsynaptic potentials generated by the magnetic pulse in the form of transcranial evoked potentials (TEPs) [7,18,19,20] or oscillations [21,22], and, thus, to have a readout from cortical areas outside M1.
In the present study, we assessed TMS–evoked EEG responses from the pre–SMA in a cohort of young, healthy individuals, with the aim of investigating signal changes due to different COs and SIs. To assess cortical excitability, we analysed TEPs in terms of discrete peaks and local mean field potential (LMFP), a reference–free measure commonly used to measure local excitability of a specific area, from a cluster of nearby electrodes [6,23,24]. To assess cortical oscillations, we computed TMS–related spectral perturbation (TRSP), a measure reflecting the power of TMS–evoked response in the frequency domain [21,22,24]. These variables were computed locally to the SMA since we were interested in the response of this area obtained with three different SIs (100%, 120% and 140%) of resting motor threshold (RMT) and four different COs (0°, 90°, 180° and 270°, starting with the coil handle pointing posteriorly on the transverse plane and proceeding counterclockwise). Overall, we found that both CO and SI are critical in determining TMS–EEG responses, and this likely suggests that, by varying these parameters, it is possible to preferentially stimulate different neural elements within the pre–SMA.

2. Materials and Methods

2.1. Subjects and Experimental Sessions

Sixteen healthy subjects (7 female, age 29.5 ± 4.6), all right–handed [25], were enrolled in the study. They had no history of neurological or psychiatric diseases and were not taking drugs active at the central nervous system level at the time of the experiment. Subjects gave their written consent to participate prior to the experimental sessions. All procedures were performed in accordance with the Declaration of Helsinki. Approval to conduct the experiments was obtained from the human subjects review board of the University College London. The experiment consisted of two sessions; participants underwent a total of 12 blocks of TMS–EEG recording, in which four different COs (0°, 90°, 180° and 270°, starting with the coil handle pointing posteriorly on the transverse plane and proceeding counterclockwise) and, for each of them, three different SIs (100%, 120% and 140% RMT) were tested (Figure 1, panel A). Each recording block consisted of 100 TMS pulses, and the order of the blocks was randomised across the two sessions and within the same session. Single–pulse TMS was applied over the right pre–SMA during EEG recording with an interpulse interval of 4 ± 10% s (3.6–4.4 s); EMG was recorded bilaterally from the first dorsal interosseous (FDI) muscle to check whether stimulation was able to induce MEPs.

2.2. Electric Field Modelling

To ensure that the chosen stimulation intensities could generate an electric field (E–field) in the cortex sufficient for a reliable TEP, i.e., at least 40 V/m [22,26], we computed the induced E–field over the TMS target with SimNIBS v3.2, an open–source simulation package that integrates segmentation of MRI scans, mesh generation and FEM E–field estimate [27]. The E–field was computed for the twelve stimulation conditions resulting from the combination of four COs (0°, 90°, 180°, 360°) and three SIs (100%, 120%, 140% RMT). The E–fields were estimated based on the MNI standard brain (ERNIE) provided in SimNIBS software as an anatomical reference [28].

2.3. TMS, Electromyographic Recording and Analysis

EMG activity was recorded through a pair of Ag/AgCl electrodes placed over the right and left FDI muscle in a belly–tendon montage. EMG signal was amplified and filtered (gain 1000x; bandwidth 5 Hz–2 kHz) with a Digitimer D360 (Digitimer Ltd., Welwyn Garden City, UK), then digitally converted with a CED 1401 analogue–to–digital laboratory interface (Cambridge Electronic Design Ltd., Milton, Cambridge, UK). Single–pulse TMS was performed using a Magstim 200 stimulator with a 70 mm figure of eight coil (Magstim Co Ltd., Whitland, UK), which produces stimuli with a monophasic waveform and a pulse width of ~80 µs.
The RMT was measured on the FDI hotspot of the right M1, which was defined as the site where TMS evoked the largest MEP in the left FDI muscle. The RMT was calculated as the lowest magnetic stimulator intensity able to evoke an MEP of at least 50 µV in 5 out of 10 consecutive trials in the relaxed FDI [29,30]. Peak–to–peak amplitudes of MEP were calculated by summing the absolute minimum and maximum EMG values in a region between 10 and 60 ms after the TMS artefact. The right pre–SMA cortical site where TMS was applied was identified based on MNI coordinates from previous studies (x = 20; y = 6; z = 62) [31]. These coordinates were used to initially locate the area and to maintain the coil in the correct position throughout the stimulation blocks by using a Brainsight navigation system (Brainbox Ltd., Cardiff, UK) coupled with a Polaris Spectra optical measurement system (Northern Digital Inc, Waterloo, Canada). The chosen location is more lateral and posterior compared with the pre–SMA identified elsewhere [32], and it has been specifically linked to correct stopping behaviour during a stop signal task [31]. Note that this area is far from craniofacial muscles, whose activation by TMS could potentially affect the EEG signal [23,33]. An estimated individualised MRI scan in the MNI space was used for all the participants. Previous studies demonstrated that the mean accuracy of the estimated MRI scans is comparable with the spatial resolution of TMS [34]. To increase the reproducibility of our data between the two sessions, in the first experiment, the position of the recording electrodes was digitised in the MNI space, and the coordinates were used to ensure their accurate placement in the second experiment.

2.4. Electroencephalographic Recording and Analysis

EEG was recorded using a TMS–compatible amplifier (TruScan EEG, Deymed Diagnostic s.r.o, Hronov, Czech Republic). The system minimises the TMS–induced artefact by removing the AC coupling from 0.5 ms before to 1 ms after the TMS pulse. Signals were recorded from 62 TMS–compatible Ag/AgCl passive electrodes mounted on a cap produced by the same manufacturer, according to the 10–20 international EEG system, including: Fp1, Fp2, F7, F3, Fz, F4, F8, T7, C3, Cz, C4, T8, P7, P3, Pz, P4, P8, O1, O2, FC5, FC1, FC2, FC6, AF7, CP5, CP1, CP2, CP6, AF8, Oz, FPz, AF3, FC3, AF4, C6, Iz, FC4, FT8, F5, C2, F1, AFz, C5, F2, TP7, F6, C1, FCz, FT7, CP3, CPz, CP4, TP8, P5, P1, P2, P6, PO7, PO3, POz, PO4 and PO8.
Recordings were referenced online to linked mastoids, and the ground electrode was placed above the nasion. Skin impedances were kept below 5 kΩ, and the sampling frequency during recording was 3000 Hz. In order to mask the TMS–induced noise and avoid possible auditory evoked potentials, participants wore earplugs continuously playing a white noise mixed with specific time–varying frequencies of the TMS [35,36]. Additionally, a 0.5 cm foam layer was placed underneath the coil to minimise bone conduction of the TMS click and scalp sensation caused by coil vibration.
Offline EEG preprocessing was performed with EEGLAB 14.1.1 [37] with the addition of some functions included in the TMS−EEG signal analyser (TESA) toolbox [38] and in Fieldtrip open−source MATLAB toolbox [39], all running in MATLAB environment (Version 2016b, MathWorks Inc., Natick, MA, USA).
EEG signal recorded in all blocks was epoched (−1.3 to 1.3 s) and demeaned using a baseline from −1000 to −10 ms. Epochs were visually inspected, and those with excessively noisy EEG were excluded. The TMS artefact was removed from −5 to 10 ms around the trigger and interpolated by means of a cubic function. A first round of independent component decomposition analysis (ICA) was run using a fastICA algorithm. Only the 15 components explaining the largest variance were plotted in a time window ranging from −200 to 500 ms, and those reflecting residual scalp muscle or voltage decay artefacts were eliminated after visual inspection based on time, frequency, scalp distribution and amplitude criteria [40,41]. After this, a band–pass (1–100 Hz) and a band–stop (48–52 Hz) fourth−order Butterworth filter were applied. The signal was further epoched (−1 to 1 s) to reduce possible edge artefacts, and a second round of fastICA was performed to remove residual artefacts (e.g., eyeblinks, continuous muscle activity, etc.). Lastly, a common average reference was applied.
Since the aim of the study was to investigate the local effects related to the variation in CO and SI, we calculated several TMS–EEG measures in a cluster of electrodes surrounding the stimulation site (Fz, F2, FCz, FC2). To assess the local cortical activation induced by TMS in the time domain, we computed the LMFP as the square root of squared TEPs averaged across the four channels of interest, as performed in previous studies [23,24,42,43]. The time–domain signal was further analysed, considering discrete TEP peaks. After computing the grand average signal across all subjects and conditions, we identified five main waves by visual inspection, peaking at 26 (PI), 44 (PII), 62 (PIII), 118 (PIV) and 198 (PV) ms. Maximum (for positive) and minimum (for negative) amplitude values were extracted for each subject, condition and peak within the following time windows: 22–30 ms (PI), 39–49 ms (PII), 58–66 ms (PIII), 108–128 ms (PIV) and 178–228 ms (PV) (Figure 2).
For the time–frequency analysis of TEP, spectral estimations of the EEG epochs were obtained for frequencies between 1 and 60 Hz (1 Hz resolution) and times in the interval from −500 to 500 ms. A sliding window (5 ms steps), linearly increasing its length across frequencies (1 cycle length for 1 Hz up to 7 cycles for 60 Hz), was used to extract amplitude and power values of all time–frequency bins. These values were estimated using the multitapers method as implemented in fieldtrip’s ft_freqanalysis function. For these estimations, Hanning tapers were used, and the amount of spectral smoothing factor was set to 0.1 times the frequency analysed in each bin. Then, TRSP was computed as follows:
T R S P ( f , t ) = 1 n k = 1 n | F k ( f , t ) | 2
where, for n trials, the spectral power estimate F was computed at trial k, at frequency f and time t [37]. Both the LMFP and the TRSP were measured in two time windows (early: 10–70 ms; late: 70–250 ms); these should reflect more local vs. more distributed brain activation, respectively [36]. Another rationale for the indicated time windows was to compare the LMFP and TRSP, the latter showing separate response clusters in the two intervals (Figure 1, panels C and D). TSRP values were averaged from 10 to 70 ms for γ (31–48 Hz) and β (14–30 Hz) frequency bands and from 70 to 250 ms for α (8–13 Hz) and θ (5–8 Hz) frequency bands.

2.5. Statistical Analysis

To measure differences in LMFP induced by TMS delivered with different CO and SI, we performed two separate two–way repeated measures ANOVAs on early and late LMFP, respectively. Factors of analysis were “CO” (0°, 90°, 180° and 270°) and “SI” (100%, 120% and 140% RMT). Five ANOVAs with the same structure were used to assess the effects of CO and SI on TEP peaks. Possible differences in TRSP were assessed again by means of repeated measures ANOVAs with the same factors and levels as before. This time, values were averaged both across frequencies and time windows (10–70 ms for β–γ range and 70–250 ms for θ–α range). To check whether TMS over the pre−SMA effectively stimulated M1, we also performed a three–way ANOVA on MEP amplitude, using “side” (left, right) “SI” (100%, 120% and 140% RMT) and “CO” (0°, 90°, 180° and 270°) as factors of analysis. Before undergoing ANOVAs, normal distribution of data was assessed by means of Shapiro–Wilk’s test. All p-values < 0.05 were considered significant. Greenhouse–Geisser correction was used when necessary to correct for non–sphericity (i.e., Mauchly’s test < 0.05). To correct for multiple comparisons, Bonferroni’s correction was used for main effects, interactions and post hoc analyses following the ANOVAs. Statistical analyses were performed with IBM SPSS v24 (Armonk, NY, USA: IBM Corp).

3. Results

The test sessions were well tolerated, and no participants reported any side effects. Results are expressed as average ± standard deviation (SD) if not otherwise specified. Average RMT was 55.1 ± 9.4 of the maximum stimulator output (MSO). TMS pulses induced a pattern of negative and positive deflections consistent with previous literature [44,45] (Figure 1, panel B; Figure 2). The TRSP showed a prominent and early increase in power peaking in the γ frequency range, compatible with activation of medial motor areas [26], and a later increase in the θ–α frequency bands (Figure 1, panel D).

3.1. Electric Field Modelling

Table 1 and Figure 3 report the E–fields computed for the 12 conditions explored. Our results showed that, even with the lowest intensity of stimulation (100% RMT), the estimated E–fields were well above the threshold of 40 V/m to evoke a reliable TMS–evoked EEG response [22,26]. To note that the E–fields obtained at 0° and 90° CO were higher than those at 180° and 270°; this might, to an extent, have contributed to our results (see below).

3.2. Local Mean Field Potential and TEP Peaks

Globally, LMFP induced by the 0° and 90° CO was larger than 180° and 270°. The ANOVA on early LMFP showed a significant main effect of “CO” (F3,39 = 5.736, p = 0.02), “SI” (F2,26 = 9.731, p < 0.001) and a significant “CO × SI” interaction (F6,78 = 3.12, p = 0.04). Post hoc comparisons showed a consistently larger LMFP at 0° and 90° CO compared with 180° and 270°. Although this was true for all the tested SIs, the effect reached statistical significance only at 140% RMT (all p values < 0.01). A similar pattern was found in the ANOVA on late LMFP, where a significant main effect of “CO” (F3,39 = 3.506, p = 0.024), “SI” (F2,26 = 11.804, p < 0.001) and a significant “CO × SI” interaction (F6,78 = 4.633, p < 0.001) were found. This time, only LMFP evoked by 0° was larger than 180° and 270°, and this again occurred only when an SI of 140% RMT was used (Figure 4).
The results of the ANOVAs on TEP peaks are summarised in Table 2 and Figure 5. Overall, these analyses confirmed the trend observed for the LMFP, i.e., larger amplitude values for 0° and 90° CO compared with 180° and 270°. Again, this was clearer with increasing stimulation intensities, except for PI and PIV, in which some statistically significant differences were observed for 100% RMT as well (Figure 5).

3.3. TMS–Related Spectral Perturbation

TRSP showed a trend similar to LMFP and TEP peaks in the early time window, in which power in the β–γ bands was considered. The ANOVA on TRSP at β–γ frequencies showed a significant main effect of “CO” (F3,39 = 8.144, p < 0.001), “SI” (F3,39 = 53.502, p < 0.001) and a significant “CO × SI” interaction (F3,39 = 3.816, p = 0.02). Post hoc comparisons showed that 0° induced a higher TRSP than 180° and 270° at 100% SI (p = 0.023 and 0.015, respectively). At higher intensities (120% and 140% RMT), both 0° and 90° COs induced a higher TRSP compared with 180° and 270° (all p values < 0.01). Interestingly, results were very different in the TRSP measured in the θ–α bands. In this case, the ANOVA showed a significant main effect of “CO” (F3,39 = 7.467, p = 0.01), “SI” (F3,39 = 27.433, p < 0.001) and a significant “CO × SI” interaction (F3,39 = 5.318, p < 0.001). This time, a higher TRSP was induced by stimulation at 180° compared with the other three COs, and this was true for 120% and 140% RMT SIs (all p values < 0.01) (Figure 4).

3.4. Motor Evoked Potentials

No stimulation condition elicited MEPs in the right FDI, whereas clear MEPs were recorded in the left FDI with a 90° CO, both at 120% and 140% SIs. The related ANOVA showed a significant main effect of “side” (F1,15 = 10.309, p = 0.006), “CO” (F3,45 = 11.947, p < 0.001), “SI” (F2,30 = 10.414, p < 0.001), as well as significant interactions of “side × CO” (F3,45 = 11.579, p < 0.001), “side × SI” (F2,30 = 7.979, p = 0.002), “CO × SI” (F2,30 = 10.474, p < 0.001) and “side × CO × SI” (F2,30 = 9.129, p < 0.001). Post hoc comparisons demonstrated that, when using 90° CO, both at 120% and 140% SIs, MEP size was larger compared with all other conditions (all p values < 0.05) (Figure 4).

4. Discussion

In this study, we demonstrated that TMS coil direction, already known to influence responses from M1, also plays a major role when the pre–SMA is stimulated. Specifically, early time–domain, late time–domain, and early time/frequency–domain EEG responses to TMS were generally larger when the coil handle was pointing posteriorly and 90° laterally. Albeit this result might be partly explained by a larger E–field induced in these two conditions (Figure 3), it is likely that local neuronal dynamics contributed as well, since late TRSP was greater when the coil was pointing at 180°. Additionally, the present results suggest that stimulating the pre–SMA with a 90° orientation activates the right M1, thus limiting the spatial specificity of the stimulation.
Despite a large number of studies, TMS CO has not been sufficiently addressed when targeting the pre–SMA, both with single–pulse and repetitive TMS. So far, different TMS studies have used lateral [16,17] or posterior [15,46,47,48] CO. In some cases, CO was not specified [13,14], and in others, a cone coil was used [49,50,51]. Several studies specifically addressed the issue of TMS–EEG responses linked to different COs when M1 was stimulated. Results were mixed, with at least one reporting differences in TEP peaks [52]. The LMFP has been reported to be less sensitive to coil orientation, provided that stimulation intensity is adjusted by RMT when M1 is stimulated [6,7]. The effects of CO in other brain areas have been less studied. Casarotto and colleagues, for example [33], targeted a cortical spot in Brodmann’s area 6, possibly within the SMA; they found that CO affects EEG responses induced by TMS. However, in this study, as well as others investigating medial premotor areas [53,54], a biphasic stimulator was used, thus limiting inferences about the effects of CO on TEP.
The first important finding of the present study was that the early LMFP and corresponding TEP peaks (I–III), calculated between 10 and 70 ms after the TMS pulse, are larger with COs of 0° and 90° compared with 180° and 270°. Statistical significance was reached only at the maximal SI used (140% RMT) for the LMFP, while it also occurred for lower stimulation intensities for some peaks. This is in apparent contrast with previous investigations in M1, in which differences in MEP latency and modulation were obtained with a low stimulation intensity [5,55]. In this regard, it is worth noting that neurons in the pre–SMA have a higher threshold for electrical stimulation than those in M1 [56]. Additionally, they vary considerably in terms of sensitivity to somatosensory inputs and threshold to evoke movements [57]. It is thus possible that a higher TMS intensity recruits neurons with a higher threshold. Alternatively, since the spatial spread of the TMS–induced electric field depends on stimulation intensity, activation of cortical sites adjacent to the stimulated one may have contributed to CO–specific effects at high stimulation intensities. In this regard, since early TEPs mostly reflect local cortical activation [36], the spread of electric field within the same area may result in larger TEPs via in−phase summation of homogenous neural activity, whereas concurrent activation of different areas may have resulted in smaller TEPs via out–of–phase cancellation.
Interestingly, the observed effect was clearer on the β–γ TRSP compared with the LMFP. TRSP measured with a 0° CO was higher than 180° and 270° at all SIs. Similarly, TRSP at 90° was higher than 180° and 270° when a SI of 120% and 140% RMT were used. A possible explanation is that the response of the stimulated neurons is only partly phase–locked; thus, a measure such as the TRSP, which also takes into account nonphase–locked activity, is likely to be more sensitive to variations in the effectiveness of TMS. Such a notion has already been suggested to explain a higher sensitivity of TRSP compared with LMFP in detecting cortical excitability changes induced by continuous theta–burst stimulation [24]. The exact nature of the neural elements stimulated in the present study is difficult to ascertain. Previous data suggested that early components of TEP, which have frequency content in the beta range, could be a reflection of either excitatory postsynaptic potentials mediated by NMDA receptors [58] or inhibitory postsynaptic potentials due to GABAa receptors activation [59], the two hypotheses being not mutually exclusive.
When looking at longer latency responses, LMFP and TEP peaks showed a pattern similar to the earlier potentials. In this case, only 0° CO gave rise to significantly larger LMFP than 180° and 270°, and only at 140% RMT SI; by contrast, TEP PIV was modulated by CO even with lower stimulation intensities. Surprisingly, the pattern of θ–α TRSP was very different, i.e., it was larger with a CO of 180° compared with all the other COs, and the effect was significant at 120% and 140% RMT. It is thus possible that a CO of 180° stimulates neural elements which give rise to a more desynchronised global response; this might explain the discrepancy between time–domain signals and TRSP. Again, we can only speculate about the cellular mechanisms involved in these late responses to TEP. Some lines of evidence suggest that they might be mostly generated by local or interhemispheric inhibitory circuits involving GABAb receptors [60].
A further comment is needed with regards to the comparison between 0° and 90° COs. Albeit with slight differences in terms of post hoc comparisons with other COs, both gave rise to similar LMFP, TEP peaks and TRSP. However, a CO of 90° also evoked clear MEPs in the contralateral but not the ipsilateral FDI when SIs of 120% and 140% RMT were used. The most likely explanation is that when the coil handle points laterally, the magnetic field generated stimulates M1, which is slightly lateral to the pre–SMA. If this is the case, part of the signals observed when using a CO of 90° might be due to activation of M1, and thus it might be overestimated. Even if the measured cortical signal was not contaminated by activity generated in M1, the latter is nonetheless activated as indicated by MEP generation; thus, a CO of 90° might not be selective enough for stimulation of the pre–SMA. Lastly, part of the TMS–EEG signals observed in this case may be due to reafferent activity due to muscle twitch caused by MEP [43,61]. M1 activation might be partly dependent on the pre–SMA coordinates used here. We chose a site based on the peak of fMRI activation found by Sharp and colleagues in the context of motor response inhibition [31]. Compared to pre–SMA coordinates used in other studies [14], this site is more posterior and lateral, close to the border of the superior frontal gyrus and, thus, closer to M1; this might have facilitated M1 activation during our experiments.
It has recently been proposed that part of the TEP might be due to non–neural sources, including auditory and somatosensory evoked potentials, the latter being generated by activation of craniofacial muscles and cutaneous nerve fibres under the area of TMS stimulation [62,63,64]. However, as noted in the experimental procedure section, we took great care in minimising the contribution of non–neural sources by masking the TMS click with an appropriate noise and by placing a thin foam layer underneath the coil. It is also worth noticing that, due to the stimulated area being medial, scalp muscle activation was minimised [23,33] and that the remaining sources of the EEG signal not caused by direct cortical activation (e.g., afferent volleys generated by activation of somatosensory fibres in the skin) would not differ across experimental blocks. Additionally, since we compared conditions that were homogeneous in terms of stimulation intensity and scalp position, we are confident that our conclusion applies even if part of the responses were due to sensory input.

5. Conclusions

To sum up, the present findings suggest that, by varying CO and SI, it might be possible to target different neural populations in the pre–SMA, with different properties in terms of stimulation threshold and response synchronisation to the TMS pulse. It is known that neurons within the pre–SMA subserve different functions. These include modality–specific changes of activity in the context of reaching movements [65] and encoding information for the numerical order of components in sequences of movements [57,66]. Additionally, some neurons in the pre–SMA might be involved in non–motor tasks, such as comprehension of accelerated speech [67], and they may show variable sensitivity to somatosensory stimuli and threshold to evoke movements [57]. This evidence considered, further research is warranted to understand whether the different TMS–EEG responses observed in the present work can be used to characterise the subset of neurons in the pre–SMA involved in diverse behaviour. This will have importance in the context of physiological studies, especially involving motor planning, and in clinical studies in which the pre–SMA has received attention as a potential therapeutic target with repetitive TMS [68].

Author Contributions

Conceptualisation, E.P.C., J.I., A.B., V.R., S.T., A.L., J.C.R. and L.R.; methodology, E.P.C., J.I., A.B., V.R., S.T. and L.R.; software, G.L., J.I., V.R. and L.R.; validation, G.L., J.I., A.B., V.R., S.T., A.L. and L.R.; formal analysis, G.L., J.I., V.R. and L.R.; investigation, E.P.C., A.B., V.R., S.T., A.L. and L.R.; resources, J.C.R.; data curation, E.P.C., G.L., V.R., A.L. and L.R.; writing—original draft preparation, E.P.C., G.L., J.I., A.B., V.R., S.T., A.L. and L.R.; writing—review and editing, E.P.C., G.L., J.I., A.B., V.R., S.T., A.L., J.C.R. and L.R.; visualisation, V.R. and L.R.; supervision, J.C.R. and L.R.; project administration, L.R.; funding acquisition, J.C.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the human subjects review board of the University College London (protocol code 10037/001).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are available upon request to the corresponding author.

Acknowledgments

E.P.C was supported by the European Commission with a Marie Sklodowska–Curie Individual Fellowship H2020–MSCA–IF Grant No. 798957. We would like to acknowledge the help of Dr Lucia Mencarelli for the E–field modelling.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Day, B.L.; Dressler, D.; de Noordhout, A.M.; Marsden, C.D.; Nakashima, K.; Rothwell, J.C.; Thompson, P.D. Electric and magnetic stimulation of human motor cortex: Surface EMG and single motor unit responses. J. Physiol. 1989, 412, 449–473. [Google Scholar] [CrossRef] [PubMed]
  2. Di Lazzaro, V.; Rothwell, J.C. Corticospinal activity evoked and modulated by non–invasive stimulation of the intact human motor cortex. J. Physiol. 2014, 592, 4115–4128. [Google Scholar] [CrossRef] [PubMed]
  3. Hanajima, R.; Ugawa, Y.; Terao, Y.; Sakai, K.; Furubayashi, T.; Machii, K.; Kanazawa, I. Paired–pulse magnetic stimulation of the human motor cortex: Differences among I waves. J. Physiol. 1998, 509, 607–618. [Google Scholar] [CrossRef] [PubMed]
  4. D’Ostilio, K.; Goetz, S.M.; Hannah, R.; Ciocca, M.; Chieffo, R.; Chen, J.A.; Peterchev, A.V.; Rothwell, J.C. Effect of coil orientation on strength–duration time constant and I–wave activation with controllable pulse parameter transcranial magnetic stimulation. Clin. Neurophysiol. 2016, 127, 675–683. [Google Scholar] [CrossRef] [Green Version]
  5. Hannah, R.; Rothwell, J.C. Pulse Duration as Well as Current Direction Determines the Specificity of Transcranial Magnetic Stimulation of Motor Cortex during Contraction. Brain Stimul. 2017, 10, 106–115. [Google Scholar] [CrossRef] [Green Version]
  6. Rocchi, L.; Spampinato, D.A.; Pezzopane, V.; Orth, M.; Bisiacchi, P.S.; Rothwell, J.C.; Casula, E.P. Cerebellar noninvasive neuromodulation influences the reactivity of the contralateral primary motor cortex and surrounding areas: A TMS–EMG–EEG study. Cerebellum 2022. [Google Scholar] [CrossRef]
  7. Casula, E.P.; Rocchi, L.; Hannah, R.; Rothwell, J.C. Effects of pulse width, waveform and current direction in the cortex: A combined cTMS–EEG study. Brain Stimul. 2018, 11, 1063–1070. [Google Scholar] [CrossRef] [Green Version]
  8. Hannah, R.; Cavanagh, S.E.; Tremblay, S.; Simeoni, S.; Rothwell, J.C. Selective Suppression of Local Interneuron Circuits in Human Motor Cortex Contributes to Movement Preparation. J. Neurosci. 2018, 38, 1264–1276. [Google Scholar] [CrossRef] [Green Version]
  9. Rawji, V.; Modi, S.; Rocchi, L.; Jahanshahi, M.; Rothwell, J.C. Proactive inhibition is marked by differences in the pattern of motor cortex activity during movement preparation and execution. J. Neurophysiol. 2022, 127, 819–828. [Google Scholar] [CrossRef]
  10. Picard, N.; Strick, P.L. Motor areas of the medial wall: A review of their location and functional activation. Cereb. Cortex 1996, 6, 342–353. [Google Scholar] [CrossRef]
  11. Picard, N.; Strick, P.L. Imaging the premotor areas. Curr. Opin. Neurobiol. 2001, 11, 663–672. [Google Scholar] [CrossRef]
  12. Nachev, P.; Kennard, C.; Husain, M. Functional role of the supplementary and pre–supplementary motor areas. Nat. Rev. Neurosci. 2008, 9, 856–869. [Google Scholar] [CrossRef] [PubMed]
  13. Rushworth, M.F.; Hadland, K.A.; Paus, T.; Sipila, P.K. Role of the human medial frontal cortex in task switching: A combined fMRI and TMS study. J. Neurophysiol. 2002, 87, 2577–2592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lee, H.W.; Lu, M.S.; Chen, C.Y.; Muggleton, N.G.; Hsu, T.Y.; Juan, C.H. Roles of the pre–SMA and rIFG in conditional stopping revealed by transcranial magnetic stimulation. Behav. Brain Res. 2016, 296, 459–467. [Google Scholar] [CrossRef]
  15. Obeso, I.; Robles, N.; Marron, E.M.; Redolar-Ripoll, D. Dissociating the Role of the pre–SMA in Response Inhibition and Switching: A Combined Online and Offline TMS Approach. Front. Hum. Neurosci. 2013, 7, 150. [Google Scholar] [CrossRef] [Green Version]
  16. Obeso, I.; Cho, S.S.; Antonelli, F.; Houle, S.; Jahanshahi, M.; Ko, J.H.; Strafella, A.P. Stimulation of the pre–SMA influences cerebral blood flow in frontal areas involved with inhibitory control of action. Brain Stimul. 2013, 6, 769–776. [Google Scholar] [CrossRef]
  17. Chen, C.Y.; Muggleton, N.G.; Tzeng, O.J.; Hung, D.L.; Juan, C.H. Control of prepotent responses by the superior medial frontal cortex. Neuroimage 2009, 44, 537–545. [Google Scholar] [CrossRef]
  18. Ilmoniemi, R.J.; Kicic, D. Methodology for combined TMS and EEG. Brain Topogr. 2010, 22, 233–248. [Google Scholar] [CrossRef] [Green Version]
  19. Casula, E.P.; Mayer, I.M.S.; Desikan, M.; Tabrizi, S.J.; Rothwell, J.C.; Orth, M. Motor cortex synchronization influences the rhythm of motor performance in premanifest huntington’s disease. Mov. Disord. 2018, 33, 440–448. [Google Scholar] [CrossRef]
  20. Leodori, G.; De Bartolo, M.I.; Guerra, A.; Fabbrini, A.; Rocchi, L.; Latorre, A.; Paparella, G.; Belvisi, D.; Conte, A.; Bhatia, K.P.; et al. Motor Cortical Network Excitability in Parkinson’s Disease. Mov. Disord. 2022, 37, 734–744. [Google Scholar] [CrossRef]
  21. Biondi, A.; Rocchi, L.; Santoro, V.; Rossini, P.G.; Beatch, G.N.; Richardson, M.P.; Premoli, I. Spontaneous and TMS–related EEG changes as new biomarkers to measure anti–epileptic drug effects. Sci. Rep. 2022, 12, 1919. [Google Scholar] [CrossRef] [PubMed]
  22. Casula, E.P.; Pellicciari, M.C.; Bonnì, S.; Borghi, I.; Maiella, M.; Assogna, M.; Minei, M.; Motta, C.; D’Acunto, A.; Porrazzini, F.; et al. Decreased Frontal Gamma Activity in Alzheimer Disease Patients. Ann. Neurol. 2022, 92, 464–475. [Google Scholar] [CrossRef] [PubMed]
  23. Casarotto, S.; Canali, P.; Rosanova, M.; Pigorini, A.; Fecchio, M.; Mariotti, M.; Lucca, A.; Colombo, C.; Benedetti, F.; Massimini, M. Assessing the effects of electroconvulsive therapy on cortical excitability by means of transcranial magnetic stimulation and electroencephalography. Brain Topogr. 2013, 26, 326–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rocchi, L.; Ibanez, J.; Benussi, A.; Hannah, R.; Rawji, V.; Casula, E.; Rothwell, J. Variability and Predictors of Response to Continuous Theta Burst Stimulation: A TMS–EEG Study. Front. Neurosci. 2018, 12, 400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Oldfield, R.C. The assessment and analysis of handedness: The Edinburgh inventory. Neuropsychologia 1971, 9, 97–113. [Google Scholar] [CrossRef]
  26. Rosanova, M.; Casali, A.; Bellina, V.; Resta, F.; Mariotti, M.; Massimini, M. Natural frequencies of human corticothalamic circuits. J. Neurosci. 2009, 29, 7679–7685. [Google Scholar] [CrossRef]
  27. Thielscher, A.; Antunes, A.; Saturnino, G.B. Field modeling for transcranial magnetic stimulation: A useful tool to understand the physiological effects of TMS? In Proceedings of the 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Milan, Italy, 25–29 August 2015; Volume 2015, pp. 222–225. [Google Scholar] [CrossRef]
  28. Windhoff, M.; Opitz, A.; Thielscher, A. Electric field calculations in brain stimulation based on finite elements: An optimized processing pipeline for the generation and usage of accurate individual head models. Hum. Brain Mapp. 2013, 34, 923–935. [Google Scholar] [CrossRef]
  29. Rocchi, L.; Erro, R.; Antelmi, E.; Berardelli, A.; Tinazzi, M.; Liguori, R.; Bhatia, K.; Rothwell, J. High frequency somatosensory stimulation increases sensori–motor inhibition and leads to perceptual improvement in healthy subjects. Clin. Neurophysiol. 2017, 128, 1015–1025. [Google Scholar] [CrossRef] [Green Version]
  30. Rossini, P.M.; Barker, A.T.; Berardelli, A.; Caramia, M.D.; Caruso, G.; Cracco, R.Q.; Dimitrijevic, M.R.; Hallett, M.; Katayama, Y.; Lucking, C.H.; et al. Non–invasive electrical and magnetic stimulation of the brain, spinal cord and roots: Basic principles and procedures for routine clinical application. Report of an IFCN committee. Electroencephalogr. Clin. Neurophysiol. 1994, 91, 79–92. [Google Scholar] [CrossRef]
  31. Sharp, D.J.; Bonnelle, V.; De Boissezon, X.; Beckmann, C.F.; James, S.G.; Patel, M.C.; Mehta, M.A. Distinct frontal systems for response inhibition, attentional capture, and error processing. Proc. Natl. Acad. Sci. USA 2010, 107, 6106–6111. [Google Scholar] [CrossRef] [Green Version]
  32. Li, C.S.; Huang, C.; Constable, R.T.; Sinha, R. Imaging response inhibition in a stop–signal task: Neural correlates independent of signal monitoring and post–response processing. J. Neurosci. 2006, 26, 186–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Casarotto, S.; Lauro, L.J.R.; Bellina, V.; Casali, A.G.; Rosanova, M.; Pigorini, A.; Defendi, S.; Mariotti, M.; Massimini, M. EEG responses to TMS are sensitive to changes in the perturbation parameters and repeatable over time. PLoS ONE 2010, 5, e10281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Herwig, U.; Schonfeldt-Lecuona, C.; Wunderlich, A.P.; von Tiesenhausen, C.; Thielscher, A.; Walter, H.; Spitzer, M. The navigation of transcranial magnetic stimulation. Psychiatry Res. 2001, 108, 123–131. [Google Scholar] [CrossRef]
  35. Mancuso, M.; Sveva, V.; Cruciani, A.; Brown, K.; Ibáñez, J.; Rawji, V.; Casula, E.; Premoli, I.; D’Ambrosio, S.; Rothwell, J.; et al. Transcranial Evoked Potentials Can Be Reliably Recorded with Active Electrodes. Brain Sci. 2021, 11, 145. [Google Scholar] [CrossRef]
  36. Rocchi, L.; Di Santo, A.; Brown, K.; Ibáñez, J.; Casula, E.; Rawji, V.; Di Lazzaro, V.; Koch, G.; Rothwell, J. Disentangling EEG responses to TMS due to cortical and peripheral activations. Brain Stimul. 2021, 14, 4–18. [Google Scholar] [CrossRef]
  37. Delorme, A.; Makeig, S. EEGLAB: An open source toolbox for analysis of single–trial EEG dynamics including independent component analysis. J. Neurosci. Methods 2004, 134, 9–21. [Google Scholar] [CrossRef] [Green Version]
  38. Rogasch, N.C.; Sullivan, C.; Thomson, R.H.; Rose, N.S.; Bailey, N.W.; Fitzgerald, P.B.; Farzan, F.; Hernandez-Pavon, J.C. Analysing concurrent transcranial magnetic stimulation and electroencephalographic data: A review and introduction to the open–source TESA software. Neuroimage 2017, 147, 934–951. [Google Scholar] [CrossRef]
  39. Oostenveld, R.; Fries, P.; Maris, E.; Schoffelen, J.M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011, 2011, 156869. [Google Scholar] [CrossRef]
  40. Rogasch, N.C.; Thomson, R.H.; Daskalakis, Z.J.; Fitzgerald, P.B. Short–latency artifacts associated with concurrent TMS–EEG. Brain Stimul. 2013, 6, 868–876. [Google Scholar] [CrossRef]
  41. Rogasch, N.C.; Thomson, R.H.; Farzan, F.; Fitzgibbon, B.M.; Bailey, N.W.; Hernandez-Pavon, J.C.; Daskalakis, Z.J.; Fitzgerald, P.B. Removing artefacts from TMS–EEG recordings using independent component analysis: Importance for assessing prefrontal and motor cortex network properties. Neuroimage 2014, 101, 425–439. [Google Scholar] [CrossRef]
  42. Pellicciari, M.C.; Brignani, D.; Miniussi, C. Excitability modulation of the motor system induced by transcranial direct current stimulation: A multimodal approach. Neuroimage 2013, 83, 569–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Fecchio, M.; Pigorini, A.; Comanducci, A.; Sarasso, S.; Casarotto, S.; Premoli, I.; Derchi, C.C.; Mazza, A.; Russo, S.; Resta, F.; et al. The spectral features of EEG responses to transcranial magnetic stimulation of the primary motor cortex depend on the amplitude of the motor evoked potentials. PLoS ONE 2017, 12, e0184910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hill, A.T.; Rogasch, N.C.; Fitzgerald, P.B.; Hoy, K.E. TMS–EEG: A window into the neurophysiological effects of transcranial electrical stimulation in non–motor brain regions. Neurosci. Biobehav. Rev. 2016, 64, 175–184. [Google Scholar] [CrossRef] [PubMed]
  45. Leodori, G.; Rocchi, L.; Mancuso, M.; De Bartolo, M.I.; Baione, V.; Costanzo, M.; Belvisi, D.; Conte, A.; Defazio, G.; Berardelli, A. The effect of stimulation frequency on transcranial evoked potentials. Transl. Neurosci. 2022, 13, 211–217. [Google Scholar] [CrossRef]
  46. Drummond, N.M.; Cressman, E.K.; Carlsen, A.N. Offline continuous theta burst stimulation over right inferior frontal gyrus and pre–supplementary motor area impairs inhibition during a go/no–go task. Neuropsychologia 2017, 99, 360–367. [Google Scholar] [CrossRef]
  47. Obeso, I.; Wilkinson, L.; Teo, J.T.; Talelli, P.; Rothwell, J.C.; Jahanshahi, M. Theta burst magnetic stimulation over the pre–supplementary motor area improves motor inhibition. Brain Stimul. 2017, 10, 944–951. [Google Scholar] [CrossRef] [Green Version]
  48. Allen, C.; Singh, K.D.; Verbruggen, F.; Chambers, C.D. Evidence for parallel activation of the pre–supplementary motor area and inferior frontal cortex during response inhibition: A combined MEG and TMS study. R. Soc. Open Sci. 2018, 5, 171369. [Google Scholar] [CrossRef] [Green Version]
  49. Georgiev, D.; Rocchi, L.; Tocco, P.; Speekenbrink, M.; Rothwell, J.C.; Jahanshahi, M. Continuous Theta Burst Stimulation Over the Dorsolateral Prefrontal Cortex and the Pre–SMA Alter Drift Rate and Response Thresholds Respectively During Perceptual Decision–Making. Brain Stimul. 2016, 9, 601–608. [Google Scholar] [CrossRef]
  50. Mendez, J.C.; Rocchi, L.; Jahanshahi, M.; Rothwell, J.; Merchant, H. Probing the timing network: A continuous theta burst stimulation study of temporal categorization. Neuroscience 2017, 356, 167–175. [Google Scholar] [CrossRef]
  51. Conte, A.; Rocchi, L.; Nardella, A.; Dispenza, S.; Scontrini, A.; Khan, N.; Berardelli, A. Theta–burst stimulation–induced plasticity over primary somatosensory cortex changes somatosensory temporal discrimination in healthy humans. PLoS ONE 2012, 7, e32979. [Google Scholar] [CrossRef]
  52. Bonato, C.; Miniussi, C.; Rossini, P.M. Transcranial magnetic stimulation and cortical evoked potentials: A TMS/EEG co–registration study. Clin. Neurophysiol. 2006, 117, 1699–1707. [Google Scholar] [CrossRef] [PubMed]
  53. Salo, K.S.; Vaalto, S.M.; Mutanen, T.P.; Stenroos, M.; Ilmoniemi, R.J. Individual activation patterns after the stimulation of different motor areas–a TMS–EEG study. Brain Connect. 2018, 8, 420–428. [Google Scholar] [CrossRef] [PubMed]
  54. Casarotto, S.; Turco, F.; Comanducci, A.; Perretti, A.; Marotta, G.; Pezzoli, G.; Rosanova, M.; Isaias, I.U. Excitability of the supplementary motor area in Parkinson’s disease depends on subcortical damage. Brain Stimul. 2018, 12, 152–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hannah, R.; Rocchi, L.; Tremblay, S.; Rothwell, J.C. Controllable Pulse Parameter TMS and TMS–EEG As Novel Approaches to Improve Neural Targeting with rTMS in Human Cerebral Cortex. Front. Neural Circuits 2016, 10, 97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Escola, L.; Michelet, T.; Macia, F.; Guehl, D.; Bioulac, B.; Burbaud, P. Disruption of information processing in the supplementary motor area of the MPTP–treated monkey: A clue to the pathophysiology of akinesia? Brain 2003, 126, 95–114. [Google Scholar] [CrossRef]
  57. Clower, W.T.; Alexander, G.E. Movement sequence–related activity reflecting numerical order of components in supplementary and presupplementary motor areas. J. Neurophysiol. 1998, 80, 1562–1566. [Google Scholar] [CrossRef]
  58. Rogasch, N.C.; Fitzgerald, P.B. Assessing cortical network properties using TMS–EEG. Hum. Brain Mapp. 2013, 34, 1652–1669. [Google Scholar] [CrossRef]
  59. Premoli, I.; Castellanos, N.; Rivolta, D.; Belardinelli, P.; Bajo, R.; Zipser, C.; Espenhahn, S.; Heidegger, T.; Muller-Dahlhaus, F.; Ziemann, U. TMS–EEG signatures of GABAergic neurotransmission in the human cortex. J. Neurosci. 2014, 34, 5603–5612. [Google Scholar] [CrossRef] [Green Version]
  60. Ziemann, U.; Reis, J.; Schwenkreis, P.; Rosanova, M.; Strafella, A.; Badawy, R.; Muller-Dahlhaus, F. TMS and drugs revisited 2014. Clin. Neurophysiol. 2015, 126, 1847–1868. [Google Scholar] [CrossRef]
  61. Petrichella, S.; Johnson, N.; He, B. The influence of corticospinal activity on TMS–evoked activity and connectivity in healthy subjects: A TMS–EEG study. PLoS ONE 2017, 12, e0174879. [Google Scholar] [CrossRef] [Green Version]
  62. Gordon, P.C.; Desideri, D.; Belardinelli, P.; Zrenner, C.; Ziemann, U. Comparison of cortical EEG responses to realistic sham versus real TMS of human motor cortex. Brain Stimul. 2018, 11, 1322–1330. [Google Scholar] [CrossRef] [PubMed]
  63. Belardinelli, P.; Biabani, M.; Blumberger, D.M.; Bortoletto, M.; Casarotto, S.; David, O.; Desideri, D.; Etkin, A.; Ferrarelli, F.; Fitzgerald, P.B.; et al. Reproducibility in TMS–EEG studies: A call for data sharing, standard procedures and effective experimental control. Brain Stimul. 2019, 12, 787–790. [Google Scholar] [CrossRef] [Green Version]
  64. Conde, V.; Tomasevic, L.; Akopian, I.; Stanek, K.; Saturnino, G.B.; Thielscher, A.; Bergmann, T.O.; Siebner, H.R. The non–transcranial TMS–evoked potential is an inherent source of ambiguity in TMS–EEG studies. Neuroimage 2019, 185, 300–312. [Google Scholar] [CrossRef]
  65. Matsuzaka, Y.; Tanji, J. Changing directions of forthcoming arm movements: Neuronal activity in the presupplementary and supplementary motor area of monkey cerebral cortex. J. Neurophysiol. 1996, 76, 2327–2342. [Google Scholar] [CrossRef] [PubMed]
  66. Nakajima, T.; Hosaka, R.; Mushiake, H.; Tanji, J. Covert representation of second–next movement in the pre–supplementary motor area of monkeys. J. Neurophysiol. 2009, 101, 1883–1889. [Google Scholar] [CrossRef] [PubMed]
  67. Dietrich, S.; Hertrich, I.; Muller-Dahlhaus, F.; Ackermann, H.; Belardinelli, P.; Desideri, D.; Seibold, V.C.; Ziemann, U. Reduced Performance During a Sentence Repetition Task by Continuous Theta–Burst Magnetic Stimulation of the Pre–supplementary Motor Area. Front. Neurosci. 2018, 12, 361. [Google Scholar] [CrossRef]
  68. Badran, B.W.; Glusman, C.E.; Austelle, C.W.; Jenkins, S.; DeVries, W.H.; Galbraith, V.; Thomas, T.; Adams, T.G., Jr.; George, M.S.; Revuelta, G.J. A Double–Blind, Sham–Controlled Pilot Trial of Pre–Supplementary Motor Area (Pre–SMA) 1 Hz rTMS to Treat Essential Tremor. Brain Stimul. 2016, 9, 945–947. [Google Scholar] [CrossRef]
Figure 1. Panel (A): Experimental protocol. Single–pulse TMS was applied on the right pre–SMA using four different coil orientations (0°, 90°, 180°, 270°), starting with the coil handle pointing posteriorly on the transverse plane and proceeding counterclockwise). For each coil orientation, three stimulation intensities were used (100%, 120% and 140% of the RMT). Panel (B): Example of TEP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition. Each line represents a signal from one electrode; all 62 recording electrodes are plotted. Panel (C): Example of LMFP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition. Panel (D): Example of TRSP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition.
Figure 1. Panel (A): Experimental protocol. Single–pulse TMS was applied on the right pre–SMA using four different coil orientations (0°, 90°, 180°, 270°), starting with the coil handle pointing posteriorly on the transverse plane and proceeding counterclockwise). For each coil orientation, three stimulation intensities were used (100%, 120% and 140% of the RMT). Panel (B): Example of TEP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition. Each line represents a signal from one electrode; all 62 recording electrodes are plotted. Panel (C): Example of LMFP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition. Panel (D): Example of TRSP obtained by averaging signals from all subjects in the 0° CO and 140% SI condition.
Brainsci 12 01358 g001
Figure 2. Grande average TEP across subjects and conditions. Each red line depicts the TEP from one electrode. The thick black line indicates the TEP averaged across the four electrodes from which the LMFP was calculated (Fz, F2, FCz, FC2). The yellow panels indicate the time windows from which maximum/minimum values of TEP peak amplitudes were extracted (see text for details).
Figure 2. Grande average TEP across subjects and conditions. Each red line depicts the TEP from one electrode. The thick black line indicates the TEP averaged across the four electrodes from which the LMFP was calculated (Fz, F2, FCz, FC2). The yellow panels indicate the time windows from which maximum/minimum values of TEP peak amplitudes were extracted (see text for details).
Brainsci 12 01358 g002
Figure 3. Pictorial representation of the E–field induced by TMS in the 12 conditions explored (see text for details).
Figure 3. Pictorial representation of the E–field induced by TMS in the 12 conditions explored (see text for details).
Brainsci 12 01358 g003
Figure 4. Summary of results on LMFP and TRSP. Panel (A): Early LMFP (10–70 ms) was larger at 0° and 90° compared with 180° and 270°. Although this was true for all the tested SIs, the effect reached statistical significance only at 140% RMT. Panel (B): Late LMFP (70–250 ms) evoked by 0° CO was larger than 180° and 270°, and again, this occurred only when an SI of 140% RMT was used. Panel (C): For TRSP in the β–γ frequency range, 0° CO induced a higher TRSP than 180° and 270° at 100% SI. At higher intensities, both 0° and 90° COs induced a higher TRSP compared with 180° and 270°. Panel (D): TRSP in the θ–α frequency range showed higher values for 180° compared with the other CO, both at 120% and 140% RMT SIs. Panel (E): MEP recorded from the left FDI were larger with 90° CO when using 120% and 140% RMT SIs compared with all other conditions. Panel (F): MEP recorded from the right FDI showed no difference across different CO and SI. Note: TRSP is expressed as the relative change compared with baseline. Error bars indicate the standard error of the mean. Asterisks indicate statistical significance (* p < 0.05).
Figure 4. Summary of results on LMFP and TRSP. Panel (A): Early LMFP (10–70 ms) was larger at 0° and 90° compared with 180° and 270°. Although this was true for all the tested SIs, the effect reached statistical significance only at 140% RMT. Panel (B): Late LMFP (70–250 ms) evoked by 0° CO was larger than 180° and 270°, and again, this occurred only when an SI of 140% RMT was used. Panel (C): For TRSP in the β–γ frequency range, 0° CO induced a higher TRSP than 180° and 270° at 100% SI. At higher intensities, both 0° and 90° COs induced a higher TRSP compared with 180° and 270°. Panel (D): TRSP in the θ–α frequency range showed higher values for 180° compared with the other CO, both at 120% and 140% RMT SIs. Panel (E): MEP recorded from the left FDI were larger with 90° CO when using 120% and 140% RMT SIs compared with all other conditions. Panel (F): MEP recorded from the right FDI showed no difference across different CO and SI. Note: TRSP is expressed as the relative change compared with baseline. Error bars indicate the standard error of the mean. Asterisks indicate statistical significance (* p < 0.05).
Brainsci 12 01358 g004
Figure 5. Summary of results on TEP peaks: panel (A), PI; panel (B), PII; panel (C), PIII; panel (D), PIV; panel (E), PV. Error bars indicate standard error of the mean. Asterisks indicate statistical significance (* p < 0.05).
Figure 5. Summary of results on TEP peaks: panel (A), PI; panel (B), PII; panel (C), PIII; panel (D), PIV; panel (E), PV. Error bars indicate standard error of the mean. Asterisks indicate statistical significance (* p < 0.05).
Brainsci 12 01358 g005
Table 1. E–field values for the 12 conditions explored (measured in V/m).
Table 1. E–field values for the 12 conditions explored (measured in V/m).
RMT90°180°270°
100% RMT110108.890.590.6
120% RMT132.2130.9108.8108.9
140% RMT154.4152.9127.1127.2
Table 2. Summary statistics of the ANOVAs on TEP peaks.
Table 2. Summary statistics of the ANOVAs on TEP peaks.
COSICO × SI
F3,45pF2,30pF6,90p
PI6.7630.0011.0310.3690.2890.941
PII1.7090.1792.6460.0873.5200.004
PIII7.697<0.0010.4750.6262.3550.037
PIV4.1310.0118.8990.0010.2830.944
PV7.316<0.0014.5040.0190.7430.617
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Casula, E.P.; Leodori, G.; Ibáñez, J.; Benussi, A.; Rawji, V.; Tremblay, S.; Latorre, A.; Rothwell, J.C.; Rocchi, L. The Effect of Coil Orientation on the Stimulation of the Pre–Supplementary Motor Area: A Combined TMS and EEG Study. Brain Sci. 2022, 12, 1358. https://doi.org/10.3390/brainsci12101358

AMA Style

Casula EP, Leodori G, Ibáñez J, Benussi A, Rawji V, Tremblay S, Latorre A, Rothwell JC, Rocchi L. The Effect of Coil Orientation on the Stimulation of the Pre–Supplementary Motor Area: A Combined TMS and EEG Study. Brain Sciences. 2022; 12(10):1358. https://doi.org/10.3390/brainsci12101358

Chicago/Turabian Style

Casula, Elias P., Giorgio Leodori, Jaime Ibáñez, Alberto Benussi, Vishal Rawji, Sara Tremblay, Anna Latorre, John C. Rothwell, and Lorenzo Rocchi. 2022. "The Effect of Coil Orientation on the Stimulation of the Pre–Supplementary Motor Area: A Combined TMS and EEG Study" Brain Sciences 12, no. 10: 1358. https://doi.org/10.3390/brainsci12101358

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