Bringing transcranial mapping into shape: Sulcus-aligned mapping captures motor somatotopy in human primary motor hand area

Corticospinal excitability during motor preparation of upper extremity reaches reflects flexor muscle synergies: A novel principal component-based motor evoked potential analyses

Background

Previous research has shown that noninvasive brain stimulation can be used to study how the central nervous system (CNS) prepares the execution of a motor task. However, these previous studies have been limited to a single muscle or single degree of freedom movements (e.g., wrist flexion). It is currently unclear if the findings of these studies generalize to multi-joint movements involving multiple muscles, which may be influenced by kinematic redundancy and muscle synergies.

Objective

The objective of this study was to characterize corticospinal excitability during motor preparation in the cortex prior to functional upper extremity reaches.

Methods

20 participants without neurological impairments volunteered for this study. During the experiment, the participants reached for a cup in response to a visual "Go Cue". Prior to movement onset, we used transcranial magnetic stimulation (TMS) to stimulate the motor cortex and measured the changes in motor evoked potentials (MEPs) in several upper extremity muscles. We varied each participant's initial arm posture and used a novel synergy-based MEP analysis to examine the effect of muscle coordination on MEPs. Additionally, we varied the timing of the stimulation between the Go Cue and movement onset to examine the time course of motor preparation.

Results

We found that synergies with strong proximal muscle (shoulder and elbow) components emerged as the stimulation was delivered closer to movement onset, regardless of arm posture, but MEPs in the distal (wrist and finger) muscles were not facilitated. We also found that synergies varied with arm posture in a manner that reflected the muscle coordination of the reach.

Conclusions

We believe that these findings provide useful insight into the way the CNS plans motor skills.

Transcranial magnetic stimulation over supramarginal gyrus stimulates primary motor cortex directly and impairs manual dexterity: implications for TMS focality

Based on human motor cortex, the effective spatial resolution of transcranial magnetic stimulation (TMS) is often described as 5-20 mm, because small changes in TMS coil position can have large effects on motor-evoked potentials (MEPs). MEPs are often studied at rest, with muscles relaxed. During muscle contraction and movement, corticospinal excitability is higher, thresholds for effective stimulation are lower, and MEPs can be evoked from larger regions of scalp, so the effective spatial resolution of TMS is larger. We found that TMS over the supramarginal gyrus (SMG) impaired manual dexterity in the grooved pegboard task. It also resulted in short-latency MEPs in hand muscles, despite the coil being 55 mm away from the motor cortex hand area (M1). MEPs might be evoked by either a specific corticospinal connection from SMG or a remote but direct electromagnetic stimulation of M1. To distinguish these alternatives, we mapped MEPs across the scalp during rest, isotonic contraction, and manual dexterity tasks and ran electric field simulations to model the expected M1 activation from 27 scalp locations and four coil orientations. We also systematically reviewed studies using TMS during movement. Across five experiments, TMS over SMG reliably evoked MEPs during hand movement. These MEPs were consistent with direct M1 stimulation and substantially decreased corticospinal thresholds during natural movement. Systematic review suggested that 54 published experiments may have suffered from similar motor activation confounds. Our results have implications for the assumed spatial resolution of TMS, and especially when TMS is presented within 55 mm of the motor cortex.NEW & NOTEWORTHY Transcranial magnetic stimulation (TMS) is often described as having an effective spatial resolution of ∼10 mm, because of the limited area of the scalp on which TMS produces motor-evoked potentials (MEPs) in resting muscles. We find that during natural hand movement TMS evokes MEPs from a much larger scalp area, in particular when stimulating over the supramarginal gyrus 55 mm away. Our results show that TMS can be effective at much larger distances than generally assumed.

The Influence of Resistance Training Experience on the Efficacy of Motor Imagery for Acutely Increasing Corticospinal Excitability

Both motor imagery and resistance-training enhance motor function and corticospinal excitability. We tested the hypothesis that young participants with significant resistance-training experience would show heightened corticospinal excitability during a single session of motor imagery training. Fifty-six participants (mean ± SD age = 22 ± 2 years) were divided into resistance-trained and untrained groups. Forty-one upper-body resistance trained (21 males, 20 females; mean ± SD relative one repetition maximum bench press = 0.922 ± 0.317 kg/kg) and 15 untrained (4 males, 11 females; mean ± SD relative one repetition maximum bench press = 0.566 ± 0.175 kg/kg) participants visited the laboratory on three separate occasions. The first visit served as the familiarization session. During visits 2 and 3, participants engaged in a hand/wrist motor imagery protocol or rested quietly (control condition) in a randomized order. Before and after the interventions, single-pulse transcranial magnetic stimulation (TMS) over the motor cortex was used to measure resting motor-evoked potential amplitude of the first dorsal interosseous muscle. Our main finding was that motor imagery acutely increased corticospinal excitability by ~64% (marginal means pre = 784.1 µV, post = 1246.6 µV; p < 0.001, d = 0.487). However, there was no evidence that the increase in corticospinal excitability was influenced by resistance-training experience. We suspect that our results may have been influenced by the specific nature of the motor imagery task. Our findings have important implications for motor imagery prescription and suggest that motor imagery training may be equally beneficial for both resistance-trained and untrained populations. This study was prospectively registered at ClinicalTrials.gov (Identifier: NCT03889548).

A fast direct solver for surface-based whole-head modeling of transcranial magnetic stimulation

When modeling transcranial magnetic stimulation (TMS) in the brain, a fast and accurate electric field solver can support interactive neuronavigation tasks as well as comprehensive biophysical modeling. We formulate, test, and disseminate a direct (i.e., non-iterative) TMS solver that can accurately determine global TMS fields for any coil type everywhere in a high-resolution MRI-based surface model with ~ 200,000 or more arbitrarily selected observation points within approximately 5 s, with the solution time itself of 3 s. The solver is based on the boundary element fast multipole method (BEM-FMM), which incorporates the latest mathematical advancement in the theory of fast multipole methods-an FMM-based LU decomposition. This decomposition is specific to the head model and needs to be computed only once per subject. Moreover, the solver offers unlimited spatial numerical resolution. Despite the fast execution times, the present direct solution is numerically accurate for the default model resolution. In contrast, the widely used brain modeling software SimNIBS employs a first-order finite element method that necessitates additional mesh refinement, resulting in increased computational cost. However, excellent agreement between the two methods is observed for various practical test cases following mesh refinement, including a biophysical modeling task. The method can be readily applied to a wide range of TMS analyses involving multiple coil positions and orientations, including image-guided neuronavigation. It can even accommodate continuous variations in coil geometry, such as flexible H-type TMS coils. The FMM-LU direct solver is freely available to academic users.

Reassessing referral of touch following peripheral deafferentation: The role of contextual bias

Some amputees have been famously reported to perceive facial touch as arising from their phantom hand. These referred sensations have since been replicated across multiple neurological disorders and were classically interpreted as a perceptual correlate of cortical plasticity. Common to all these and related studies is that participants might have been influenced in their self-reports by the experimental design or related contextual biases. Here, we investigated whether referred sensations reports might be confounded by demand characteristics (e.g., compliance, expectation, suggestion). Unilateral upper-limb amputees (N = 18), congenital one-handers (N = 19), and two-handers (N = 22) were repeatedly stimulated with computer-controlled vibrations on 10 body-parts and asked to report the occurrence of any concurrent sensations on their hand(s). To further manipulate expectations, we gave participants the suggestion that some of these vibrations had a higher probability to evoke referred sensations. We also assessed similarity between (phantom) hand and face representation in primary somatosensory cortex (S1), using functional Magnetic Resonance Imaging (fMRI) multivariate representational similarity analysis. We replicated robust reports of referred sensations in amputees towards their phantom hand; however, the frequency and distribution of reported referred sensations were similar across groups. Moreover, referred sensations were evoked by stimulation of multiple body-parts and similarly reported on both the intact and phantom hand in amputees. Face-to-phantom-hand representational similarity was not different in amputees' missing hand region, compared with controls. These findings weaken the interpretation of referred sensations as a perceptual correlate of S1 plasticity and reveal the need to account for contextual biases when evaluating anomalous perceptual phenomena.

Neural Plasticity Changes Induced by Motor Robotic Rehabilitation in Stroke Patients: The Contribution of Functional Neuroimaging

Robotic rehabilitation is one of the most advanced treatments helping people with stroke to faster recovery from motor deficits. The clinical impact of this type of treatment has been widely defined and established using clinical scales. The neurofunctional indicators of motor recovery following conventional rehabilitation treatments have already been identified by previous meta-analytic investigations. However, a clear definition of the neural correlates associated with robotic neurorehabilitation treatment has never been performed. This systematic review assesses the neurofunctional correlates (fMRI, fNIRS) of cutting-edge robotic therapies in enhancing motor recovery of stroke populations in accordance with PRISMA standards. A total of 7, of the initial yield of 150 articles, have been included in this review. Lessons from these studies suggest that neural plasticity within the ipsilateral primary motor cortex, the contralateral sensorimotor cortex, and the premotor cortices are more sensitive to compensation strategies reflecting upper and lower limbs' motor recovery despite the high heterogeneity in robotic devices, clinical status, and neuroimaging procedures. Unfortunately, the paucity of RCT studies prevents us from understanding the neurobiological differences induced by robotic devices with respect to traditional rehabilitation approaches. Despite this technology dating to the early 1990s, there is a need to translate more functional neuroimaging markers in clinical settings since they provide a unique opportunity to examine, in-depth, the brain plasticity changes induced by robotic rehabilitation.

A fast direct solver for surface-based whole-head modeling of transcranial magnetic stimulation

Background

When modeling transcranial magnetic stimulation (TMS) in the brain, a fast and accurate electric field solver can support interactive neuronavigation tasks as well as comprehensive biophysical modeling.

Objective

We formulate, test, and disseminate a direct (i.e., non-iterative) TMS solver that can accurately determine global TMS fields for any coil type everywhere in a high-resolution MRI-based surface model with ~200,000 or more arbitrarily selected observation points within approximately 5 sec, with the solution time itself of 3 sec.

Method

The solver is based on the boundary element fast multipole method (BEM-FMM), which incorporates the latest mathematical advancement in the theory of fast multipole methods - an FMM-based LU decomposition. This decomposition is specific to the head model and needs to be computed only once per subject. Moreover, the solver offers unlimited spatial numerical resolution.

Results

Despite the fast execution times, the present direct solution is numerically accurate for the default model resolution. In contrast, the widely used brain modeling software SimNIBS employs a first-order finite element method that necessitates additional mesh refinement, resulting in increased computational cost. However, excellent agreement between the two methods is observed for various practical test cases following mesh refinement, including a biophysical modeling task.

Conclusion

The method can be readily applied to a wide range of TMS analyses involving multiple coil positions and orientations, including image-guided neuronavigation. It can even accommodate continuous variations in coil geometry, such as flexible H-type TMS coils. The FMM-LU direct solver is freely available to academic users.

Multi-scale and cross-dimensional TMS mapping: A proof of principle in patients with Parkinson's disease and deep brain stimulation

Introduction

Transcranial magnetic stimulation (TMS) mapping has become a critical tool for exploratory studies of the human corticomotor (M1) organization. Here, we propose to gather existing cutting-edge TMS-EMG and TMS-EEG approaches into a combined multi-dimensional TMS mapping that considers local and whole-brain excitability changes as well as state and time-specific changes in cortical activity. We applied this multi-dimensional TMS mapping approach to patients with Parkinson's disease (PD) with Deep brain stimulation (DBS) of the sub-thalamic nucleus (STN) ON and OFF. Our goal was to identifying one or several TMS mapping-derived markers that could provide unprecedent new insights onto the mechanisms of DBS in movement disorders.

Methods

Six PD patients (1 female, mean age: 62.5 yo [59-65]) implanted with DBS-STN for 1 year, underwent a robotized sulcus-shaped TMS motor mapping to measure changes in muscle-specific corticomotor representations and a movement initiation task to probe state-dependent modulations of corticospinal excitability in the ON (using clinically relevant DBS parameters) and OFF DBS states. Cortical excitability and evoked dynamics of three cortical areas involved in the neural control of voluntary movements (M1, pre-supplementary motor area - preSMA and inferior frontal gyrus - IFG) were then mapped using TMS-EEG coupling in the ON and OFF state. Lastly, we investigated the timing and nature of the STN-to-M1 inputs using a paired pulse DBS-TMS-EEG protocol.

Results

In our sample of patients, DBS appeared to induce fast within-area somatotopic re-arrangements of motor finger representations in M1, as revealed by mediolateral shifts of corticomuscle representations. STN-DBS improved reaction times while up-regulating corticospinal excitability, especially during endogenous motor preparation. Evoked dynamics revealed marked increases in inhibitory circuits in the IFG and M1 with DBS ON. Finally, inhibitory conditioning effects of STN single pulses on corticomotor activity were found at timings relevant for the activation of inhibitory GABAergic receptors (4 and 20 ms).

Conclusion

Taken together, these results suggest a predominant role of some markers in explaining beneficial DBS effects, such as a context-dependent modulation of corticospinal excitability and the recruitment of distinct inhibitory circuits, involving long-range projections from higher level motor centers and local GABAergic neuronal populations. These combined measures might help to identify discriminative features of DBS mechanisms towards deep clinical phenotyping of DBS effects in Parkinson's Disease and in other pathological conditions.

Reduced connectivity of primary auditory and motor cortices during exposure to auditory white noise

Auditory white noise (WN) is widely used in daily life for inducing sleep, and in neuroscience to mask unwanted environmental noise and cues. However, WN was recently reported to influence corticospinal excitability and behavioral performance. Here, we expand previous preliminary findings on the influence of WN exposure on cortical functioning, and we hypothesize that it may modulate cortical connectivity. We tested our hypothesis by performing magnetoencephalography in 20 healthy subjects. WN reduces cortical connectivity of the primary auditory and motor regions with very distant cortical areas, showing a right lateralized connectivity reduction for primary motor cortex. The present results, together with previous finding concerning WN impact on corticospinal excitability and behavioral performance, further support the role of WN as a modulator of cortical function. This suggest avoiding its unrestricted use as a masking tool, while purposely designed and controlled WN application could be exploited to harness brain function and to treat neuropsychiatric conditions.

Revisiting the focality of non-invasive brain stimulation - implications for studies of human cognition

Non-invasive brain stimulation techniques are popular tools to investigate brain function in health and disease. Although transcranial magnetic stimulation (TMS) is widely used in cognitive neuroscience research to probe causal structure-function relationships, studies often yield inconclusive results. To improve the effectiveness of TMS studies, we argue that the cognitive neuroscience community needs to revise the stimulation focality principle - the spatial resolution with which TMS can differentially stimulate cortical regions. In the motor domain, TMS can differentiate between cortical muscle representations of adjacent fingers. However, this high degree of spatial specificity cannot be obtained in all cortical regions due to the influences of cortical folding patterns on the TMS-induced electric field. The region-dependent focality of TMS should be assessed a priori to estimate the experimental feasibility. Post-hoc simulations allow modeling of the relationship between cortical stimulation exposure and behavioral modulation by integrating data across stimulation sites or subjects.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

No evidence for a difference in lateralization and distinctiveness level of transcranial magnetic stimulation-derived cortical motor representations over the adult lifespan

This study aimed to investigate the presence and patterns of age-related differences in TMS-based measures of lateralization and distinctiveness of the cortical motor representations of two different hand muscles. In a sample of seventy-three right-handed healthy participants over the adult lifespan, the first dorsal interosseus (FDI) and abductor digiti minimi (ADM) cortical motor representations of both hemispheres were acquired using transcranial magnetic stimulation (TMS). In addition, dexterity and maximum force levels were measured. Lateralization quotients were calculated for homolog behavioral and TMS measures, whereas the distinctiveness between the FDI and ADM representation within one hemisphere was quantified by the center of gravity (CoG) distance and cosine similarity. The presence and patterns of age-related changes were examined using linear, polynomial, and piecewise linear regression. No age-related differences could be identified for the lateralization quotient of behavior or cortical motor representations of both intrinsic hand muscles. Furthermore, no evidence for a change in the distinctiveness of the FDI and ADM representation with advancing age was found. In conclusion this work showed that lateralization and distinctiveness of cortical motor representations, as determined by means of TMS-based measures, remain stable over the adult lifespan.

Transcranial magnetic stimulation of the brain: What is stimulated? - A consensus and critical position paper

Transcranial (electro)magnetic stimulation (TMS) is currently the method of choice to non-invasively induce neural activity in the human brain. A single transcranial stimulus induces a time-varying electric field in the brain that may evoke action potentials in cortical neurons. The spatial relationship between the locally induced electric field and the stimulated neurons determines axonal depolarization. The induced electric field is influenced by the conductive properties of the tissue compartments and is strongest in the superficial parts of the targeted cortical gyri and underlying white matter. TMS likely targets axons of both excitatory and inhibitory neurons. The propensity of individual axons to fire an action potential in response to TMS depends on their geometry, myelination and spatial relation to the imposed electric field and the physiological state of the neuron. The latter is determined by its transsynaptic dendritic and somatic inputs, intrinsic membrane potential and firing rate. Modeling work suggests that the primary target of TMS is axonal terminals in the crown top and lip regions of cortical gyri. The induced electric field may additionally excite bends of myelinated axons in the juxtacortical white matter below the gyral crown. Neuronal excitation spreads ortho- and antidromically along the stimulated axons and causes secondary excitation of connected neuronal populations within local intracortical microcircuits in the target area. Axonal and transsynaptic spread of excitation also occurs along cortico-cortical and cortico-subcortical connections, impacting on neuronal activity in the targeted network. Both local and remote neural excitation depend critically on the functional state of the stimulated target area and network. TMS also causes substantial direct co-stimulation of the peripheral nervous system. Peripheral co-excitation propagates centrally in auditory and somatosensory networks, but also produces brain responses in other networks subserving multisensory integration, orienting or arousal. The complexity of the response to TMS warrants cautious interpretation of its physiological and behavioural consequences, and a deeper understanding of the mechanistic underpinnings of TMS will be critical for advancing it as a scientific and therapeutic tool.

Stimulation with acoustic white noise enhances motor excitability and sensorimotor integration

Auditory white noise (WN) is widely used in neuroscience to mask unwanted environmental noise and cues, e.g. TMS clicks. However, to date there is no research on the influence of WN on corticospinal excitability and potentially associated sensorimotor integration itself. Here we tested the hypothesis, if WN induces M1 excitability changes and improves sensorimotor performance. M1 excitability (spTMS, SICI, ICF, I/O curve) and sensorimotor reaction-time performance were quantified before, during and after WN stimulation in a set of experiments performed in a cohort of 61 healthy subjects. WN enhanced M1 corticospinal excitability, not just during exposure, but also during silence periods intermingled with WN, and up to several minutes after the end of exposure. Two independent behavioural experiments highlighted that WN improved multimodal sensorimotor performance. The enduring excitability modulation combined with the effects on behaviour suggest that WN might induce neural plasticity. WN is thus a relevant modulator of corticospinal function; its neurobiological effects should not be neglected and could in fact be exploited in research applications.

Brain Reorganization and Neural Plasticity in Elite Athletes With Physical Impairments

Use-dependent and impairment-specific brain plasticity are hypothesized to interact and enhance neural reorganization in the central nervous system (CNS) of athletes with physical impairments. Paralympic brain studies are helpful in achieving a fundamental understanding of the underlying neural mechanism related to CNS reorganization after physical therapy or athletic training. Information learned from these individuals also provides new insights into sports- and rehabilitation-related neuroscience.

Exploring the spatial resolution of TMS-EEG coupling on the sensorimotor region

The use of TMS-EEG coupling as a neuroimaging tool for the functional exploration of the human brain recently gained strong interest. If this tool directly inherits the fine temporal resolution from EEG, its spatial counterpart remains unknown. In this study, we explored the spatial resolution of TMS-EEG coupling by evaluating the minimal distance between two stimulated cortical sites that would significantly evoke different response dynamics. TMS evoked responses were mapped on the sensorimotor region in twenty participants. The stimulation grid was composed of nine targets separated between 10 and 15 mm on average. The dynamical signatures of TMS evoked activity were extracted and compared between sites using both local and remote linear regression scores and spatial generalized mixed models. We found a significant effect of the distance between stimulated sites on their dynamical signatures, neighboring sites showing differentiable response dynamics. Besides, common dynamical signatures were also found between sites up to 25-30 mm from each other. This overlap in dynamical properties decreased with distance and was stronger between sites within the same Brodmann area. Our results suggest that the spatial resolution of TMS-EEG coupling might be at least as high as 10 mm. Furthermore, our results reveal an anisotropic spatial resolution that was higher across than within the same Brodmann areas, in accordance with the TMS induced E-field modeling. Common cytoarchitectonic leading to shared dynamical properties within the same Brodmann area could also explain this anisotropy. Overall, these findings suggest that TMS-EEG benefits from the spatial resolution of TMS, which makes it an accurate technique for meso-scale brain mapping.

A Randomized, Sham-Controlled Trial of Repetitive Transcranial Magnetic Stimulation Targeting M1 and S2 in Central Poststroke Pain: A Pilot Trial

OBJECTIVES

Central poststroke pain (CPSP), a neuropathic pain condition, is difficult to treat. Repetitive transcranial magnetic stimulation (rTMS) targeted to the primary motor cortex (M1) can alleviate the condition, but not all patients respond. We aimed to assess a promising alternative rTMS target, the secondary somatosensory cortex (S2), for CPSP treatment.

MATERIALS AND METHODS

This prospective, randomized, double-blind, sham-controlled three-arm crossover trial assessed navigated rTMS (nrTMS) targeted to M1 and S2 (10 sessions, 5050 pulses per session at 10 Hz). Participants were evaluated for pain, depression, anxiety, health-related quality of life, upper limb function, and three plasticity-related gene polymorphisms including Dopamine D2 Receptor (DRD2). We monitored pain intensity and interference before and during stimulations and at one month. A conditioned pain modulation test was performed using the cold pressor test. This assessed the efficacy of the descending inhibitory system, which may transmit TMS effects in pain control.

RESULTS

We prescreened 73 patients, screened 29, and included 21, of whom 17 completed the trial. NrTMS targeted to S2 resulted in long-term (from baseline to one-month follow-up) pain intensity reduction of ≥30% in 18% (3/17) of participants. All stimulations showed a short-term effect on pain (17-20% pain relief), with no difference between M1, S2, or sham stimulations, indicating a strong placebo effect. Only nrTMS targeted to S2 resulted in a significant long-term pain intensity reduction (15% pain relief). The cold pressor test reduced CPSP pain intensity significantly (p = 0.001), indicating functioning descending inhibitory controls. The homozygous DRD2 T/T genotype is associated with the M1 stimulation response.

CONCLUSIONS

S2 is a promising nrTMS target in the treatment of CPSP. The DRD2 T/T genotype might be a biomarker for M1 nrTMS response, but this needs confirmation from a larger study.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {cm}^2$$\end{document}cm2 along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

Multi-locus transcranial magnetic stimulation system for electronically targeted brain stimulation

BACKGROUND

Transcranial magnetic stimulation (TMS) allows non-invasive stimulation of the cortex. In multi-locus TMS (mTMS), the stimulating electric field (E-field) is controlled electronically without coil movement by adjusting currents in the coils of a transducer.

OBJECTIVE

To develop an mTMS system that allows adjusting the location and orientation of the E-field maximum within a cortical region.

METHODS

We designed and manufactured a planar 5-coil mTMS transducer to allow controlling the maximum of the induced E-field within a cortical region approximately 30 mm in diameter. We developed electronics with a design consisting of independently controlled H-bridge circuits to drive up to six TMS coils. To control the hardware, we programmed software that runs on a field-programmable gate array and a computer. To induce the desired E-field in the cortex, we developed an optimization method to calculate the currents needed in the coils. We characterized the mTMS system and conducted a proof-of-concept motor-mapping experiment on a healthy volunteer. In the motor mapping, we kept the transducer placement fixed while electronically shifting the E-field maximum on the precentral gyrus and measuring electromyography from the contralateral hand.

RESULTS

The transducer consists of an oval coil, two figure-of-eight coils, and two four-leaf-clover coils stacked on top of each other. The technical characterization indicated that the mTMS system performs as designed. The measured motor evoked potential amplitudes varied consistently as a function of the location of the E-field maximum.

CONCLUSION

The developed mTMS system enables electronically targeted brain stimulation within a cortical region.

Sensorimotor integration within the primary motor cortex by selective nerve fascicle stimulation

KEY POINTS

Cortical integration of sensory inputs is crucial for dexterous movement. Short-latency somatosensory afferent inhibition of motor cortical output is typically produced by peripheral whole-nerve stimulation. We exploited intraneural multichannel electrodes used to provide sensory feedback for prosthesis control to assess whether and how selective intraneural sensory stimulation affects sensorimotor cortical circuits in humans. The activation of the primary somatosensory cortex (S1) was explored by recording scalp somatosensory evoked potentials. Sensorimotor integration was tested by measuring the inhibitory effect of the afferent stimulation on the output of the primary motor cortex (M1) generated by transcranial magnetic stimulation. We demonstrate in humans that selective intraneural sensory stimulation elicits a measurable activation of S1 and that it inhibits the output of M1 at the same time range of whole-nerve superficial stimulation.

ABSTRACT

The integration of sensory inputs in the motor cortex is crucial for dexterous movement. We recently demonstrated that a closed-loop control based on the feedback provided through intraneural multi-channel electrodes implanted in the median and ulnar nerves of a participant with upper limb amputation improved manipulation skills and increased prosthesis embodiment. Here we assessed, in the same participant, whether and how selective intraneural sensory stimulation also elicits a measurable cortical activation and affects sensorimotor cortical circuits. After estimating the activation of the primary somatosensory cortex evoked by intraneural stimulation, sensorimotor integration was investigated by testing the inhibition of primary motor cortex (M1) output to transcranial magnetic stimulation, after both intraneural and perineural stimulation. Selective sensory intraneural stimulation evoked a low-amplitude, 16 ms-latency, parietal response in the same area of the earliest component evoked by whole-nerve stimulation, compatible with fast-conducting afferent fiber activation. For the first time, we show that the same intraneural stimulation was also capable of decreasing M1 output, at the same time range of the short-latency afferent inhibition effect of whole-nerve superficial stimulation. The inhibition generated by the stimulation of channels activating only sensory fibers was stronger than the one due to intraneural or perineural stimulation of channels activating mixed fibers. We demonstrate in a human subject that the cortical sensorimotor integration inhibiting M1 output previously described after the experimental whole-nerve stimulation is present also with a more ecological selective sensory fiber stimulation. Abstract Figure: Double-sided filament electrodes (ds-FILE), bearing 16 active sites, and perineural Cuff electrodes were implanted in the median and ulnar nerve of the arm in a hand amputee (upper left panel, single nerve represented). Selectivity of stimulation (1), evoked activity in the somatosensory cortex (2), and sensorimotor integration (3) were investigated. TMS: transcranial magnetic stimulation. This article is protected by copyright. All rights reserved.

Face-hand sensorimotor interactions revealed by afferent inhibition

Reorganization of the sensorimotor cortex following permanent (e.g., amputation) or temporary (e.g., local anaesthesia) deafferentation of the hand has revealed large-scale plastic changes between the hand and face representations that are accompanied by perceptual correlates. The physiological mechanisms underlying this reorganization remain poorly understood. The aim of this study was to investigate sensorimotor interactions between the face and hand using an afferent inhibition transcranial magnetic stimulation protocol in which the motor evoked potential elicited by the magnetic pulse is inhibited when it is preceded by an afferent stimulus. We hypothesized that if face and hand representations in the sensorimotor cortex are functionally coupled, then electrocutaneous stimulation of the face would inhibit hand muscle motor responses. In two separate experiments, we delivered an electrocutaneous stimulus to either the skin over the right upper lip (Experiment 1) or right cheek (Experiment 2) and recorded muscular activity from the right first dorsal interosseous. Both lip and cheek stimulation inhibited right first dorsal interosseous motor evoked potentials. To investigate the specificity of this effect, we conducted two additional experiments in which electrocutaneous stimulation was applied to either the right forearm (Experiment 3) or right upper arm (Experiment 4). Forearm and upper arm stimulation also significantly inhibited the right first dorsal interosseous motor evoked potentials, but this inhibition was less robust than the inhibition associated with face stimulation. These findings provide the first evidence for face-to-hand afferent inhibition.

Cortical Thickness of Brain Areas Beyond Stroke Lesions and Sensory-Motor Recovery: A Systematic Review

Background: The clinical outcome of patients suffering from stroke is dependent on multiple factors. The features of the lesion itself play an important role but clinical recovery is remarkably influenced by the plasticity mechanisms triggered by the stroke and occurring at a distance from the lesion. The latter translate into functional and structural changes of which cortical thickness might be easy to quantify one of the main players. However, studies on the changes of cortical thickness in brain areas beyond stroke lesion and their relationship to sensory-motor recovery are sparse. Objectives: To evaluate the effects of cerebral stroke on cortical thickness (CT) beyond the stroke lesion and its association with sensory-motor recovery. Materials and Methods: Five electronic databases (PubMed, Embase, Web of Science, Scopus and the Cochrane Library) were searched. Methodological quality of the included studies was assessed with the Newcastle-Ottawa Scale for non-randomized controlled trials and the Risk of Bias Cochrane tool for randomized controlled trials. Results: The search strategy retrieved 821 records, 12 studies were included and risk of bias assessed. In most of the included studies, cortical thinning was seen at the ipsilesional motor area (M1). Cortical thinning can occur beyond the stroke lesion, typically in regions anatomically connected because of anterograde degeneration. Nonetheless, studies also reported cortical thickening of regions of the unaffected hemisphere, likely related to compensatory plasticity. Some studies revealed a significant correlation between changes in cortical thickness of M1 or somatosensory (S1) cortical areas and motor function recovery. Discussion and Conclusions: Following a stroke, changes in cortical thickness occur both in regions directly connected to the stroke lesion and in contralateral hemisphere areas as well as in the cerebellum. The underlying mechanisms leading to these changes in cortical thickness are still to be fully understood and further research in the field is needed. Systematic Review Registration: https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42020200539; PROSPERO 2020, identifier: CRD42020200539.

Electric Field Comparison for TMS Using Different Neuroimaging Segmentation Methods

Computational electromagnetic modeling is a powerful technique to evaluate the effects of electrical stimulation of the human brain. The results of these simulations can vary depending on the specific segmentation of the head and brain generated from the patient images. Using an existing boundary element fast multipole method (BEM-FMM) electromagnetic solver, this work evaluates the electric field differences modeled using two neuroimaging segmentation methods. A transcranial magnetic stimulation (TMS) coil targeting both the primary motor cortex and the dorsolateral prefrontal cortex (DLPFC) was simulated. Average field differences along a 100 mm line from the coil were small (2% for motor cortex, 3% for DLPFC) and the average field differences in the regions directly surrounding the target stimulation point were 5% for the motor cortex and 2% for DLPFC. More studies evaluating different coils and other segmentation options may further improve the computational modeling for robust TMS treatment.Clinical relevance- Patient-specific computational modeling will provide more information to clinicians for improved localization and targeting of neuromodulation therapies.

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

Efficient high-resolution TMS mapping of the human motor cortex by nonlinear regression

Transcranial magnetic stimulation (TMS) is a powerful tool to investigate causal structure-function relationships in the human brain. However, a precise delineation of the effectively stimulated neuronal populations is notoriously impeded by the widespread and complex distribution of the induced electric field. Here, we propose a method that allows rapid and feasible cortical localization at the individual subject level. The functional relationship between electric field and behavioral effect is quantified by combining experimental data with numerically modeled fields to identify the cortical origin of the modulated effect. Motor evoked potentials (MEPs) from three finger muscles were recorded for a set of random stimulations around the primary motor area. All induced electric fields were nonlinearly regressed against the elicited MEPs to identify their cortical origin. We could distinguish cortical muscle representation with high spatial resolution and localized them primarily on the crowns and rims of the precentral gyrus. A post-hoc analysis revealed exponential convergence of the method with the number of stimulations, yielding a minimum of about 180 random stimulations to obtain stable results. Establishing a functional link between the modulated effect and the underlying mode of action, the induced electric field, is a fundamental step to fully exploit the potential of TMS. In contrast to previous approaches, the presented protocol is particularly easy to implement, fast to apply, and very robust due to the random coil positioning and therefore is suitable for practical and clinical applications.

Influence of different transcranial magnetic stimulation current directions on the corticomotor control of lumbar erector spinae muscles during a static task

Different directions of transcranial magnetic stimulation (TMS) can activate different neuronal circuits. Whereas posteroanterior current (PA-TMS) depolarizes mainly interneurons in primary motor cortex (M1), an anteroposterior current (AP-TMS) has been suggested to activate different M1 circuits and perhaps axons from the premotor regions. Although M1 is also involved in the control of axial muscles, no study has explored whether different current directions activate different M1 circuits that may have distinct functional roles. The aim of the study was to compare the effect of different current directions (PA- and AP-TMS) on the corticomotor control and spatial cortical organization of the lumbar erector spinae muscle (LES). Thirty-four healthy participants were recruited for two independent experiments, and LES motor-evoked potentials (MEPs) were recorded. In experiment 1 (n = 17), active motor threshold (AMT), MEP latencies, recruitment curve (90% to 160% AMT), and excitatory and inhibitory intracortical mechanisms by paired-pulse TMS (80% followed by 120% AMT stimuli at 2-, 3-, 10-, and 15-ms interstimulus intervals) were tested with a double-cone (n = 12) and a figure-of-eight (n = 5) coil. In experiment 2 (n = 17), LES cortical representations were tested with PA- and AP-TMS. AMT was higher for AP- compared with PA-TMS (P = 0.002). Longer latencies with AP-TMS were present compared with PA-TMS (P = 0.017). AP-TMS produced more inhibition compared with PA-TMS at 2 ms and 3 ms (P = 0.010), but no difference was observed for longer intervals. No difference was found for recruitment curve and mapping. These findings suggest that PA- and AP-TMS may activate different cortical circuits controlling low back muscles, as proposed for hand muscles.NEW & NOTEWORTHY For the first time, anteroposterior and posteroanterior induced electric currents in the brain were compared when targeting back muscle representation with transcranial magnetic stimulation. The use of the anteroposterior current resulted in later response latency, larger inhibition probed by paired-pulse stimulation, and higher motor threshold. These important differences between current directions suggest that each of the current directions may recruit specific cortical circuits involved in the control of back muscles, similar to that for hand muscles.

TMS Motor Mapping Methodology and Reliability: A Structured Review

Motor cortical representation can be probed non-invasively using a transcranial magnetic stimulation (TMS) technique known as motor mapping. The mapping technique can influence features of the maps because of several controllable elements. Here we review the literature on six key motor mapping parameters, as well as their influence on outcome measures and discuss factors impacting their selection. 132 of 1,587 distinct records were examined in detail and synthesized to form the basis of our review. A summary of mapping parameters, their impact on outcome measures and feasibility considerations are reported to support the design and interpretation of TMS mapping studies.

Degree of improving TMS focality through a geometrically stable solution of an inverse TMS problem

The Transcranial Magnetic Stimulation (TMS) inverse problem (TMS-IP) investigated in this study aims to focus the TMS induced electric field close to a specified target point defined on the gray matter interface in the M1_HAND area while otherwise minimizing it. The goal of the study is to numerically evaluate the degree of improvement of the TMS-IP solutions relative to the well-known sulcus-aligned mapping (a projection approach with the 90^∘ local sulcal angle). In total, 1536 individual TMS-IP solutions have been analyzed for multiple target points and multiple subjects using the boundary element fast multipole method (BEM-FMM) as the forward solver. Our results show that the optimal TMS inverse-problem solutions improve the focality - reduce the size of the field "hot spot" and its deviation from the target - by approximately 21-33% on average for all considered subjects, all observation points, two distinct coil types, two segmentation types, two intracortical observation surfaces under study, and three tested values of the field threshold. The inverse-problem solutions with the maximized focality simultaneously improve the TMS mapping resolution (differentiation between neighbor targets separated by approximately 10 mm) although this improvement is quite modest. Coil position/orientation and conductivity uncertainties have been included into consideration as the corresponding de-focalization factors. The present results will change when the levels of uncertainties change. Our results also indicate that the accuracy of the head segmentation critically influences the expected TMS-IP performance.

Mapping of Motor Function with Neuronavigated Transcranial Magnetic Stimulation: A Review on Clinical Application in Brain Tumors and Methods for Ensuring Feasible Accuracy

Navigated transcranial magnetic stimulation (nTMS) has developed into a reliable non-invasive clinical and scientific tool over the past decade. Specifically, it has undergone several validating clinical trials that demonstrated high agreement with intraoperative direct electrical stimulation (DES), which paved the way for increasing application for the purpose of motor mapping in patients harboring motor-eloquent intracranial neoplasms. Based on this clinical use case of the technique, in this article we review the evidence for the feasibility of motor mapping and derived models (risk stratification and prediction, nTMS-based fiber tracking, improvement of clinical outcome, and assessment of functional plasticity), and provide collected sets of evidence for the applicability of quantitative mapping with nTMS. In addition, we provide evidence-based demonstrations on factors that ensure methodological feasibility and accuracy of the motor mapping procedure. We demonstrate that selection of the stimulation intensity (SI) for nTMS and spatial density of stimuli are crucial factors for applying motor mapping accurately, while also demonstrating the effect on the motor maps. We conclude that while the application of nTMS motor mapping has been impressively spread over the past decade, there are still variations in the applied protocols and parameters, which could be optimized for the purpose of reliable quantitative mapping.

Mapping of multiple muscles with transcranial magnetic stimulation: absolute and relative test-retest reliability

The spatial accuracy of transcranial magnetic stimulation (TMS) may be as small as a few millimeters. Despite such great potential, navigated TMS (nTMS) mapping is still underused for the assessment of motor plasticity, particularly in clinical settings. Here, we investigate the within-limb somatotopy gradient as well as absolute and relative reliability of three hand muscle cortical representations (MCRs) using a comprehensive grid-based sulcus-informed nTMS motor mapping. We enrolled 22 young healthy male volunteers. Two nTMS mapping sessions were separated by 5-10 days. Motor evoked potentials were obtained from abductor pollicis brevis (APB), abductor digiti minimi, and extensor digitorum communis. In addition to individual MRI-based analysis, we studied normalized MNI MCRs. For the reliability assessment, we calculated intraclass correlation and the smallest detectable change. Our results revealed a somatotopy gradient reflected by APB MCR having the most lateral location. Reliability analysis showed that the commonly used metrics of MCRs, such as areas, volumes, centers of gravity (COGs), and hotspots had a high relative and low absolute reliability for all three muscles. For within-limb TMS somatotopy, the most common metrics such as the shifts between MCR COGs and hotspots had poor relative reliability. However, overlaps between different muscle MCRs were highly reliable. We, thus, provide novel evidence that inter-muscle MCR interaction can be reliably traced using MCR overlaps while shifts between the COGs and hotspots of different MCRs are not suitable for this purpose. Our results have implications for the interpretation of nTMS motor mapping results in healthy subjects and patients with neurological conditions.

Simultaneous Recording of Motor Evoked Potentials in Hand, Wrist and Arm Muscles to Assess Corticospinal Divergence

The purpose of this study was to further develop methods to assess corticospinal divergence and muscle coupling using transcranial magnetic stimulation (TMS). Ten healthy right-handed adults participated (7 females, age 34.0 ± 12.9 years). Monophasic single pulses were delivered to 14 sites over the right primary motor cortex at 40, 60, 80 and 100% of maximum stimulator output (MSO), using MRI-based neuronavigation. Motor evoked potentials (MEPs) were recorded simultaneously from 9 muscles of the contralateral hand, wrist and arm. For each intensity, corticospinal divergence was quantified by the average number of muscles that responded to TMS per cortical site, coactivation across muscle pairs as reflected by overlap of cortical representations, and correlation of MEP amplitudes across muscle pairs. TMS to each muscle's most responsive site elicited submaximal MEPs in most other muscles. The number of responsive muscles per cortical site and the extent of coactivation increased with increasing intensity (ANOVA, p < 0.001). In contrast, correlations of MEP amplitudes did not differ across the 60, 80 and 100% MSO intensities (ANOVA, p = 0.34), but did differ across muscle pairs (ANOVA, p < 0.001). Post hoc analysis identified 4 sets of muscle pairs (Tukey homogenous subsets, p < 0.05). Correlations were highest for pairs involving two hand muscles and lowest for pairs that included an upper arm muscle. Correlation of MEP amplitudes may quantify varying levels of muscle coupling. In future studies, this approach may be a biomarker to reveal altered coupling induced by neural injury, neural repair and/or motor learning.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.

Dosimetry Analysis in Non-brain Tissues During TMS Exposure of Broca's and M1 Areas

For human protection, the internal electric field is used as a dosimetric quantity for electromagnetic fields lower than 5–10 MHz. According to international standards, in this frequency range, electrostimulation is the main adverse effect against which protection is needed. One of the topics to be investigated is the quantification of the internal electric field threshold levels of perception and pain. Pain has been reported as a side effect during transcranial magnetic stimulation (TMS), especially during stimulation of the Broca’s (speech) area of the brain. In this study, we designed an experiment to conduct a dosimetry analysis to quantify the internal electric field corresponding to perception and pain thresholds when targeting the Broca’s and M1 areas from magnetic stimulator exposure. Dosimetry analysis was conducted using a multi-scale analysis in an individualized head model to investigate electrostimulation in an axonal model. The main finding is that the stimulation on the primary motor cortex has higher perception and pain thresholds when compared to Broca’s area. Also, TMS-induced electric field applied to Broca’s area exhibited dependence on the coil orientation at lower electric field threshold which was found to be related to the location and thickness of pain fibers. The derived dosimetry quantities provide a scientific rationale for the development of human protection guidelines and the estimation of possible side effects of magnetic stimulation in clinical applications.

Fast computational optimization of TMS coil placement for individualized electric field targeting

Background:

During transcranial magnetic stimulation (TMS) a coil placed on the scalp is used to non-invasively modulate activity of targeted brain networks via a magnetically induced electric field (E-field). Ideally, the E-field induced during TMS is concentrated on a targeted cortical region of interest (ROI). Determination of the coil position and orientation that best achieve this objective presently requires a large computational effort.

Objective:

To improve the accuracy of TMS we have developed a fast computational auxiliary dipole method (ADM) for determining the optimum coil position and orientation. The optimum coil placement maximizes the E-field along a predetermined direction or, alternatively, the overall E-field magnitude in the targeted ROI. Furthermore, ADM can assess E-field uncertainty resulting from precision limitations of TMS coil placement protocols.

Method:

ADM leverages the electromagnetic reciprocity principle to compute rapidly the TMS induced E-field in the ROI by using the E-field generated by a virtual constant current source residing in the ROI. The framework starts by solving for the conduction currents resulting from this ROI current source. Then, it rapidly determines the average E-field induced in the ROI for each coil position by using the conduction currents and a fast-multipole method. To further speed-up the computations, the coil is approximated using auxiliary dipoles enabling it to represent all coil orientations for a given coil position with less than 600 dipoles.

Results:

Using ADM, the E-fields generated in an MRI-derived head model when the coil is placed at 5900 different scalp positions and 360 coil orientations per position (over 2.1 million unique configurations) can be determined in under 15 min on a standard laptop computer. This enables rapid extraction of the optimum coil position and orientation as well as the E-field variation resulting from coil positioning uncertainty. ADM is implemented in SimNIBS 3.2.

Conclusion:

ADM enables the rapid determination of coil placement that maximizes E-field delivery to a specific brain target. This method can find the optimum coil placement in under 15 min enabling its routine use for TMS. Furthermore, it enables the fast quantification of uncertainty in the induced E-field due to limited precision of TMS coil placement protocols, enabling minimization and statistical analysis of the E-field dose variability.

Focality of the Induced E-Field Is a Contributing Factor in the Choice of TMS Parameters: Evidence from a 3D Computational Model of the Human Brain

Transcranial magnetic stimulation (TMS) is a promising, non-invasive approach in the diagnosis and treatment of several neurological conditions. However, the specific results in the cortex of the magnitude and spatial distribution of the secondary electrical field (E-field) resulting from TMS at different stimulation sites/orientations and varied TMS parameters are not clearly understood. The objective of this study is to identify the impact of TMS stimulation site and coil orientation on the induced E-field, including spatial distribution and the volume of activation in the cortex across brain areas, and hence demonstrate the need for customized optimization, using a three-dimensional finite element model (FEM). A considerable difference was noted in E-field values and distribution at different brain areas. We observed that the volume of activated cortex varied from 3000 to 7000 mm3 between the selected nine clinically relevant coil locations. Coil orientation also changed the induced E-field by a maximum of 10%, and we noted the least optimal values at the standard coil orientation pointing to the nose. The volume of gray matter activated varied by 10% on average between stimulation sites in homologous brain areas in the two hemispheres of the brain. This FEM simulation model clearly demonstrates the importance of TMS parameters for optimal results in clinically relevant brain areas. The results show that TMS parameters cannot be interchangeably used between individuals, hemispheres, and brain areas. The focality of the TMS induced E-field along with its optimal magnitude should be considered as critical TMS parameters that should be individually optimized.

Large-scale analysis of interindividual variability in theta-burst stimulation data: Results from the 'Big TMS Data Collaboration'

BACKGROUND

Many studies have attempted to identify the sources of interindividual variability in response to theta-burst stimulation (TBS). However, these studies have been limited by small sample sizes, leading to conflicting results.

OBJECTIVE/HYPOTHESIS

This study brought together over 60 TMS researchers to form the 'Big TMS Data Collaboration', and create the largest known sample of individual participant TBS data to date. The goal was to enable a more comprehensive evaluation of factors driving TBS response variability.

METHODS

118 corresponding authors of TMS studies were emailed and asked to provide deidentified individual TMS data. Mixed-effects regression investigated a range of individual and study level variables for their contribution to iTBS and cTBS response variability.

RESULTS

430 healthy participants' TBS data was pooled across 22 studies (mean age = 41.9; range = 17-82; females = 217). Baseline MEP amplitude, age, target muscle, and time of day significantly predicted iTBS-induced plasticity. Baseline MEP amplitude and timepoint after TBS significantly predicted cTBS-induced plasticity.

CONCLUSIONS

This is the largest known study of interindividual variability in TBS. Our findings indicate that a significant portion of variability can be attributed to the methods used to measure the modulatory effects of TBS. We provide specific methodological recommendations in order to control and mitigate these sources of variability.

The Clinical Utility of Transcranial Magnetic Stimulation in Determining Hemispheric Dominance for Language: A Magnetoencephalography Comparison Study

PURPOSE

Transcranial magnetic stimulation (TMS) has recently emerged as a noninvasive alternative to the intracarotid sodium amytal (Wada) procedure for establishing hemispheric dominance (HD) for language. The accuracy of HD determined by TMS was examined by comparing against the HD derived by magnetoencephalography (MEG), a prominent clinical technique with excellent concordance with the Wada procedure.

METHODS

Sixty-seven patients (54 patients ≤18 years) underwent language mapping with TMS and MEG as part of clinical epilepsy and tumor presurgical assessment. Language was mapped in MEG during an auditory word recognition paradigm, and a laterality index was calculated using the number of dipoles and their spatial extent in the two hemispheres. Transcranial magnetic stimulation language mapping was performed as patients performed a naming task, and TMS-induced speech disruptions were recorded during 5-Hz TMS applied to anterior and posterior language cortices. Transcranial magnetic stimulation laterality index was estimated using the number and type of speech disruption in the language regions of each hemisphere.

RESULTS

Transcranial magnetic stimulation and MEG estimates of HD were concordant in 42 (63%) patients, resulting in a sensitivity of 74% and a specificity of 72%. The overall accuracy of TMS was 73%, equivalent to an odds ratio of 7.35.

CONCLUSIONS

In this first large-scale comparative study in a clinical population, we demonstrate that TMS is a safe and reliable noninvasive tool in determining HD for language. Improving the accuracy of TMS by optimizing TMS parameters and improving task choice will further facilitate the use of TMS to characterize language function, especially in pediatrics.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Concordance Between Transcranial Magnetic Stimulation and Functional Magnetic Resonance Imaging (MRI) Derived Localization of Language in a Clinical Cohort

Transcranial magnetic stimulation (TMS) is a newer noninvasive language mapping tool that is safe and well-tolerated by children. We examined the accuracy of TMS-derived language maps in a clinical cohort by comparing it against functional magnetic resonance imaging (MRI)-derived language map. The number of TMS-induced speech disruptions and the volume of activation during functional MRI tasks were localized to Brodmann areas for each modality in 40 patients with epilepsy or brain tumor. We examined the concordance between TMS- and functional MRI-derived language maps by deriving statistical performance metrics for TMS including sensitivity, specificity, accuracy, and diagnostic odds ratio. Brodmann areas 6, 44, and 9 in the frontal lobe and 22 and 40 in the temporal lobe were the most commonly identified language areas by both modalities. Overall accuracy of TMS compared to functional MRI in localizing language cortex was 71%, with a diagnostic odds ratio of 1.27 and higher sensitivity when identifying left hemisphere regions. TMS was more accurate in determining the dominant hemisphere for language with a diagnostic odds ratio of 6. This study is the first to examine the accuracy of the whole brain language map derived by TMS in the largest cohort examined to date. While this comparison against functional MRI confirmed that TMS reliably localizes cortical areas that are not essential for speech function, it demonstrated only slight concordance between TMS- and functional MRI-derived language areas. That the localization of specific language cortices by TMS demonstrated low accuracy reveals a potential need to use concordant tasks between the modalities and other avenues for further optimization of TMS parameters.

Alterations of hand sensorimotor function and cortical motor representations over the adult lifespan

Using a cross sectional design, we aimed to identify the effect of aging on sensorimotor function and cortical motor representations of two intrinsic hand muscles, as well as the course and timing of those changes. Furthermore, the link between cortical motor representations, sensorimotor function, and intracortical inhibition and facilitation was investigated. Seventy-seven participants over the full adult lifespan were enrolled. For the first dorsal interosseus (FDI) and abductor digiti minimi (ADM) muscle, cortical motor representations, GABAA-mediated short-interval intracortical inhibition (SICI), and glutamate-mediated intracortical facilitation (ICF) were assessed using transcranial magnetic stimulation over the dominant primary motor cortex. Additionally, participants' dexterity and force were measured. Linear, polynomial, and piecewise linear regression analyses were conducted to identify the course and timing of age-related differences. Our results demonstrated variation in sensorimotor function over the lifespan, with a marked decline starting around the mid-thirties. Furthermore, an age-related reduction in cortical motor representation volume and maximal MEP of the FDI, but not for ADM, was observed, occurring mainly until the mid-forties. Area of the cortical motor representation did not change with advancing age. Furthermore, cortical motor representations, sensorimotor function, and measures of intracortical inhibition and facilitation were not interrelated.

Use-Dependent Plasticity in Human Primary Motor Hand Area: Synergistic Interplay Between Training and Immobilization

Training and immobilization are powerful drivers of use-dependent plasticity in human primary motor hand area (M1_HAND). In young right-handed volunteers, corticomotor representations of the left first dorsal interosseus and abductor digiti minimi muscles were mapped with neuronavigated transcranial magnetic stimulation (TMS) to elucidate how finger-specific training and immobilization interact within M1_HAND. A first group of volunteers trained to track a moving target on a smartphone with the left index or little finger for one week. Linear sulcus shape-informed TMS mapping revealed that the tracking skill acquired with the trained finger was transferred to the nontrained finger of the same hand. The cortical representations of the trained and nontrained finger muscle converged in proportion with skill transfer. In a second group, the index or little finger were immobilized for one week. Immobilization alone attenuated the corticomotor representation and pre-existing tracking skill of the immobilized finger. In a third group, the detrimental effects of finger immobilization were blocked by concurrent training of the nonimmobilized finger. Conversely, immobilization of the nontrained fingers accelerated learning in the adjacent trained finger during the first 2 days of training. Together, the results provide novel insight into use-dependent cortical plasticity, revealing synergistic rather than competitive interaction patterns within M1_HAND.

Cortical reorganization of lower-limb motor representations in an elite archery athlete with congenital amputation of both arms

•

We investigated cortical reorganization in an amputated archer who used his feet.

•

Lower-limb motor representations were examined using fMRI and TMS mapping.

•

M1 areas innervating lower-limb muscles were larger in the amputated athlete.

•

The toe and knee representations were expanded towards the lateral part of the M1.

•

Paralympic athletes have a unique and dynamic M1 plasticity.

Despite their disabilities, top Paralympic athletes have better motor skills than able-bodied athletes. However, the neural underpinnings of these better motor skills remain unclear. We investigated the reorganization of the primary motor cortex (M1) in a Paralympic athlete with congenital amputation of both arms who holds the world record for the farthest accurate shot in archery (Amputee Archer: AA). We recorded brain activity during contraction of right toe, ankle, knee, and hip joint muscles in the AA and 12 able-bodied control subjects using functional magnetic resonance imaging. The results revealed that M1 activation was more widespread in the AA compared with control subjects during all tasks, and shifted towards the lateral part of the M1 during contraction of toe and knee muscles. We also conducted a motor mapping experiment using navigated transcranial magnetic stimulation. The M1 area receiving stimulation elicited motor-evoked potentials from the toe, lower-leg, and thigh muscles, which were larger in the AA compared with 12 control subjects. Furthermore, the AA's motor maps were shifted towards the lateral side of M1. These results suggest an expansion of lower-limb M1 representation towards the lateral side of M1, including the trunk and upper-limb representations, and an expansion of the area of corticomotor neurons innervating the lower limb muscles in the AA. This unique M1 reorganization could underpin the AA's excellent archery performance in the absence of upper limbs. The current results suggest that Paralympic athletes may exhibit extreme M1 plasticity, which could arise through a combination of rigorous long-term motor training and compensatory M1 reorganization for missing body parts.

The effect of coil placement and orientation on the assessment of focal excitability in motor mapping with navigated transcranial magnetic stimulation

BACKGROUND

Navigated transcranial magnetic stimulation (nTMS) is used for mapping muscle representations in the primary motor cortex. We used sulcus-aligned mapping and electric field (E-field) modeling to investigate the excitability of the motor hand area for further understanding the methodological limitations of nTMS.

NEW METHOD

We studied 10 healthy volunteers to locate the cortical target eliciting the largest responses (the hotspot) in the first dorsal interosseous (FDI) muscle. Six additional targets were placed along the central sulcus at 5-mm distances. Resting motor thresholds (rMTs) and optimal coil orientations were determined at all targets, and a conventional motor mapping was conducted. The cortical E-fields, induced by stimulating the targets with rMT intensities and optimal coil orientations, were modeled in a realistic head geometry to estimate the activated cortical sites.

RESULTS

The rMTs increased with increasing distance from the hotspot (p < 0.001). The greatest motor-evoked potential (MEP) amplitudes occurred with the coil perpendicular to the sulcus and with the coil pointing towards the hotspot or the center of gravity of the motor map. The E-field strengths at the hotspot (99±26 V/m) remained above previously estimated thresholds for activation.

COMPARISON WITH EXISTING METHODS

Depending on the target location, optimal coil orientations may deviate significantly from the conventional perpendicular-to-sulcus angle, which is often assumed optimal. These orientations seem to maintain the E-field stable in the hand knob, regardless of the sulcal shape near the stimulated target.

CONCLUSIONS

The coil orientation is crucial for the accuracy of motor mapping, and the apparent motor map may extend due to remote hotspot activation.

Conditioning transcranial magnetic stimulation of ventral premotor cortex shortens simple reaction time

The ventral premotor cortex (PMv) is a key area of the sensorimotor control loop, it subtends complex motor sequences, especially when the hand is involved. However, its specific contribution to simple motor response to sensory cue is still not completely clear. To investigate the role of PMv, we used transcranial magnetic stimulation (TMS) to interfere with its function during a simple reaction time (SRT) task. We ran two experiments where participants were required to respond as fast as possible to a median nerve stimulation (go-signal), while sub-M1-threshold single pulse TMS was delivered either on left (contralateral) PMv or right (ipsilateral to sensory stimulus and motor response) PMv, 5-65 ms after the go-signal. TMS delivered on either PMv up to 25 ms after the go-signal shortened reaction time. This is the time window compatible with the arrive of sensory afferences, as if conditioning before sensory afferences arrive lower the threshold needed to release the pre-planned motor program to the primary motor cortex. This is in line with a putative PMv function of buffer of pre-planned motor program not strictly lateralized in one hemisphere.

Neural effects of transcranial magnetic stimulation at the single-cell level

Transcranial magnetic stimulation (TMS) can non-invasively modulate neural activity in humans. Despite three decades of research, the spatial extent of the cortical area activated by TMS is still controversial. Moreover, how TMS interacts with task-related activity during motor behavior is unknown. Here, we applied single-pulse TMS over macaque parietal cortex while recording single-unit activity at various distances from the center of stimulation during grasping. The spatial extent of TMS-induced activation is remarkably restricted, affecting the spiking activity of single neurons in an area of cortex measuring less than 2 mm in diameter. In task-related neurons, TMS evokes a transient excitation followed by reduced activity, paralleled by a significantly longer grasping time. Furthermore, TMS-induced activity and task-related activity do not summate in single neurons. These results furnish crucial experimental evidence for the neural effects of TMS at the single-cell level and uncover the neural underpinnings of behavioral effects of TMS.

Transcranial Magnetic Stimulation (TMS) can modulate human brain activity, but the extent of the cortical area activated by TMS is unclear. Here, the authors show that TMS affects monkey single neuron activity in an area less than 2 mm diameter, while TMS-induced activity and task-related activity do not summate.