Where does TMS Stimulate the Motor Cortex? Combining Electrophysiological Measurements and Realistic Field Estimates to Reveal the Affected Cortex Position

M2M-InvNet: Human Motor Cortex Mapping From Multi-Muscle Response Using TMS and Generative 3D Convolutional Network

Transcranial magnetic stimulation (TMS) is often applied to the motor cortex to stimulate a collection of motor evoked potentials (MEPs) in groups of peripheral muscles. The causal interface between TMS and MEP is the selective activation of neurons in the motor cortex; moving around the TMS 'spot' over the motor cortex causes different MEP responses. A question of interest is whether a collection of MEP responses can be used to identify the stimulated locations on the cortex, which could potentially be used to then place the TMS coil to produce chosen sets of MEPs. In this work we leverage our previous report on a 3D convolutional neural network (CNN) architecture that predicted MEPs from the induced electric field, to tackle an inverse imaging task in which we start with the MEPs and estimate the stimulated regions on the motor cortex. We present and evaluate five different inverse imaging CNN architectures, both conventional and generative, in terms of several measures of reconstruction accuracy. We found that one architecture, which we propose as M2M-InvNet, consistently achieved the best performance.

Verification of neuronavigated TMS accuracy using structured-light 3D scans

Objective.To investigate the reliability and accuracy of the manual three-point co-registration in neuronavigated transcranial magnetic stimulation (TMS). The effect of the error in landmark pointing on the coil placement and on the induced electric and magnetic fields was examined.Approach.The position of the TMS coil on the head was recorded by the neuronavigation system and by 3D scanning for ten healthy participants. The differences in the coil locations and orientations and the theoretical error values for electric and magnetic fields between the neuronavigated and 3D scanned coil positions were calculated. In addition, the sensitivity of the coil location on landmark accuracy was calculated.Main results.The measured distances between the neuronavigated and 3D scanned coil locations were on average 10.2 mm, ranging from 3.1 to 18.7 mm. The error in angles were on average from two to three degrees. The coil misplacement caused on average a 29% relative error in the electric field with a range from 9% to 51%. In the magnetic field, the same error was on average 33%, ranging from 10% to 58%. The misplacement of landmark points could cause a 1.8-fold error for the coil location.Significance.TMS neuronavigation with three landmark points can cause a significant error in the coil position, hampering research using highly accurate electric field calculations. Including 3D scanning to the process provides an efficient method to achieve a more accurate coil position.

Predicting the hotspot location and motor threshold prior to transcranial magnetic stimulation using electric field modelling

Objective.To investigate whether the motor threshold (MT) and the location of the motor hotspot in transcranial magnetic stimulation (TMS) can be predicted with computational models of the induced electric field.Approach.Individualized computational models were constructed from structural magnetic resonance images of ten healthy participants, and the induced electric fields were determined with the finite element method. The models were used to optimize the location and direction of the TMS coil on the scalp to produce the largest electric field at a predetermined cortical target location. The models were also used to predict how the MT changes as the magnetic coil is moved to various locations over the scalp. To validate the model predictions, the motor evoked potentials were measured from the first dorsal interosseous (FDI) muscle with TMS in the ten participants. Both computational and experimental methods were preregistered prior to the experiments.Main results.Computationally optimized hotspot locations were nearly as accurate as those obtained using manual hotspot search procedures. The mean Euclidean distance between the predicted and the measured hotspot locations was approximately 1.3 cm with a 0.8 cm bias towards the anterior direction. Exploratory analyses showed that the bias could be removed by changing the cortical target location that was used for the prediction. The results also indicated a statistically significant relationship (p< 0.001) between the calculated electric field and the MT measured at several locations on the scalp.Significance.The results show that the individual TMS hotspot can be located using computational analysis without stimulating the subject or patient even once. Adapting computational modelling would save time and effort in research and clinical use of TMS.

Transcranial magnetic stimulation mapping of the motor cortex: comparison of five estimation algorithms

Background

There are currently five different kinds of transcranial magnetic stimulation (TMS) motor mapping algorithms available, from ordinary point-based algorithms to advanced field-based algorithms. However, there have been only a limited number of comparison studies conducted, and they have not yet examined all of the currently available algorithms. This deficiency impedes the judicious selection of algorithms for application in both clinical and basic neuroscience, and hinders the potential promotion of a potential superior algorithm. Considering the influence of algorithm complexity, further investigation is needed to examine the differences between fMRI peaks and TMS cortical hotspots that were identified previously.

Methods

Twelve healthy participants underwent TMS motor mapping and a finger-tapping task during fMRI. The motor cortex TMS mapping results were estimated by five algorithms, and fMRI activation results were obtained. For each algorithm, the prediction error was defined as the distance between the measured scalp hotspot and optimized coil position, which was determined by the maximum electric field strength in the estimated motor cortex. Additionally, the study identified the minimum number of stimuli required for stable mapping. Finally, the location difference between the TMS mapping cortical hotspot and the fMRI activation peak was analyzed.

Results

The projection yielded the lowest prediction error (5.27 ± 4.24 mm) among the point-based algorithms and the association algorithm yielded the lowest (6.66 ± 3.48 mm) among field-based estimation algorithms. The projection algorithm required fewer stimuli, possibly resulting from its suitability for the grid-based mapping data collection method. The TMS cortical hotspots from all algorithms consistently deviated from the fMRI activation peak (20.52 ± 8.46 mm for five algorithms).

Conclusion

The association algorithm might be a superior choice for clinical applications and basic neuroscience research, due to its lower prediction error and higher estimation sensitivity in the deep cortical structure, especially for the sulcus. It also has potential applicability in various other TMS domains, including language area mapping and more. Otherwise, our results provide further evidence that TMS motor mapping intrinsically differs from fMRI motor mapping.

In-vivoverified anatomically aware deep learning for real-time electric field simulation

Objective.Transcranial magnetic stimulation (TMS) has emerged as a prominent non-invasive technique for modulating brain function and treating mental disorders. By generating a high-precision magnetically evoked electric field (E-field) using a TMS coil, it enables targeted stimulation of specific brain regions. However, current computational methods employed for E-field simulations necessitate extensive preprocessing and simulation time, limiting their fast applications in the determining the optimal coil placement.Approach.We present an attentional deep learning network to simulate E-fields. This network takes individual magnetic resonance images and coil configurations as inputs, firstly transforming the images into explicit brain tissues and subsequently generating the local E-field distribution near the target brain region. Main results. Relative to the previous deep-learning simulation method, the presented method reduced the mean relative error in simulated E-field strength of gray matter by 21.1%, and increased the correlation between regional E-field strengths and corresponding electrophysiological responses by 35.0% when applied into another dataset.In-vivoTMS experiments further revealed that the optimal coil placements derived from presented method exhibit comparable stimulation performance on motor evoked potentials to those obtained using computational methods. The simplified preprocessing and increased simulation efficiency result in a significant reduction in the overall time cost of traditional TMS coil placement optimization, from several hours to mere minutes.Significance.The precision and efficiency of presented simulation method hold promise for its application in determining individualized coil placements in clinical practice, paving the way for personalized TMS treatments.

Modeling the effects of transcranial magnetic stimulation on spatial attention

Objectives. Transcranial magnetic stimulation (TMS) has been widely used to modulate brain activity in healthy and diseased brains, but the underlying mechanisms are not fully understood. Previous research leveraged biophysical modeling of the induced electric field (E-field) to map causal structure-function relationships in the primary motor cortex. This study aims at transferring this localization approach to spatial attention, which helps to understand the TMS effects on cognitive functions, and may ultimately optimize stimulation schemes.Approach. Thirty right-handed healthy participants underwent a functional magnetic imaging (fMRI) experiment, and seventeen of them participated in a TMS experiment. The individual fMRI activation peak within the right inferior parietal lobule (rIPL) during a Posner-like attention task defined the center target for TMS. Thereafter, participants underwent 500 Posner task trials. During each trial, a 5-pulse burst of 10 Hz repetitive TMS (rTMS) was given over the rIPL to modulate attentional processing. The TMS-induced E-fields for every cortical target were correlated with the behavioral modulation to identify relevant cortical regions for attentional orientation and reorientation.Main results. We did not observe a robust correlation between E-field strength and behavioral outcomes, highlighting the challenges of transferring the localization method to cognitive functions with high neural response variability and complex network interactions. Nevertheless, TMS selectively inhibited attentional reorienting in five out of seventeen subjects, resulting in task-specific behavioral impairments. The BOLD-measured neuronal activity and TMS-evoked neuronal effects showed different patterns, which emphasizes the principal distinction between the neural activity being correlated with (or maybe even caused by) particular paradigms, and the activity of neural populations exerting a causal influence on the behavioral outcome.Significance. This study is the first to explore the mechanisms of TMS-induced attentional modulation through electrical field modeling. Our findings highlight the complexity of cognitive functions and provide a basis for optimizing attentional stimulation protocols.

The VLPFC-Engaged Voluntary Emotion Regulation: Combined TMS-fMRI Evidence for the Neural Circuit of Cognitive Reappraisal

A clear understanding of the neural circuit underlying emotion regulation (ER) is important for both basic and translational research. However, a lack of evidence based on combined neuroimaging and neuromodulation techniques calls into question (1) whether the change of prefrontal-subcortical activity intrinsically and causally contributes to the ER effect; and (2) whether the prefrontal control system directly modulates the subcortical affective system. Accordingly, we combined fMRI recordings with transcranial magnetic stimulation (TMS) to map the causal connections between the PFC and subcortical affective structures (amygdala and insula). A total of 117 human adult participants (57 males and 60 females) were included in the study. The results revealed that TMS-induced ventrolateral PFC (VLPFC) facilitation led to enhanced activity in the VLPFC and ventromedial PFC (VMPFC) as well as attenuated activity in the amygdala and insula during reappraisal but not during nonreappraisal (i.e., baseline). Moreover, the activated VLPFC intensified the prefrontal-subcortical couplings via the VMPFC during reappraisal only. This study provides combined TMS-fMRI evidence that downregulating negative emotion involves the prefrontal control system suppressing the subcortical affective system, with the VMPFC serving as a crucial hub within the VLPFC-subcortical network, suggesting an indirect pathway model of the ER circuit. Our findings outline potential protocols for improving ER ability by intensifying the VLPFC-VMPFC coupling in patients with mood and anxiety disorders.SIGNIFICANCE STATEMENT Using fMRI to examine the TMS effect, we uncovered that the opposite neural changes in prefrontal (enhanced) and subcortical (attenuated) regions are not a byproduct of emotion regulation (ER); instead, this prefrontal-subcortical activity per se causally contributes to the ER effect. Furthermore, using TMS to amplify the neural changes within the ER circuit, the "bridge" role of the VMPFC is highlighted under the reappraisal versus nonreappraisal contrast. This "perturb-and-measure" approach overcomes the correlational nature of fMRI data, helping us to identify brain regions that causally support reappraisal (the VLPFC and VMPFC) and those that are modulated by reappraisal (the amygdala and insula). The uncovered ER circuit is important for understanding the neural systems underlying reappraisal and valuable for translational research.

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.

Optimizing Individual Targeting of Fronto-Amygdala Network with Transcranial Magnetic Stimulation (TMS): Biophysical, Physiological and Behavioral Variations in People with Methamphetamine Use Disorder

Background

Previous studies in people with substance use disorders (SUDs) have implicated both the frontopolar cortex and amygdala in drug cue reactivity and craving, and amygdala-frontopolar coupling is considered a marker of early relapse risk. Accumulating data highlight that the frontopolar cortex can be considered a promising therapeutic target for transcranial magnetic stimulation (TMS) in SUDs. However, one-size-fits-all approaches to TMS targets resulted in substantial variation in both physiological and behavioral outcomes. Individualized TMS approaches to target cortico-subcortical circuits like amygdala-frontopolar have not yet been investigated in SUDs.

Objective

Here, we (1) defined individualized TMS target location based on functional connectivity of the amygdala-frontopolar circuit while people were exposed to drug-related cues, (2) optimized coil orientation based on maximizing electric field (EF) perpendicular to the individualized target, and (3) harmonized EF strength in targeted brain regions across a population.

Method

MRI data including structural, resting-state, and task-based fMRI data were collected from 60 participants with methamphetamine use disorders (MUDs). Craving scores based on a visual analog scale were collected immediately before and after the MRI session. We analyzed inter-subject variability in the location of TMS targets based on the maximum task-based connectivity between the left medial amygdala (with the highest functional activity among subcortical areas during drug cue exposure) and frontopolar cortex using psychophysiological interaction (PPI) analysis. Computational head models were generated for all participants and EF simulations were calculated for fixed vs. optimized coil location (Fp1/Fp2 vs. individualized maximal PPI location), orientation (AF7/AF8 vs. orientation optimization algorithm), and stimulation intensity (constant vs. adjusted intensity across the population).

Results

Left medial amygdala with the highest (mean ± SD: 0.31±0.29) functional activity during drug cue exposure was selected as the subcortical seed region. Amygdala-to-whole brain PPI analysis showed a significant cluster in the prefrontal cortex (cluster size: 2462 voxels, cluster peak in MNI space: [25 39 35]) that confirms cortico-subcortical connections. The location of the voxel with the most positive amygdala-frontopolar PPI connectivity in each participant was considered as the individualized TMS target (mean ± SD of the MNI coordinates: [12.6 64.23 -0.8] ± [13.64 3.50 11.01]). Individual amygdala-frontopolar PPI connectivity in each participant showed a significant correlation with VAS scores after cue exposure (R=0.27, p=0.03). Averaged EF strength in a sphere with r = 5mm around the individualized target location was significantly higher in the optimized (mean ± SD: 0.99 ± 0.21) compared to the fixed approach (Fp1: 0.56 ± 0.22, Fp2: 0.78 ± 0.25) with large effect sizes (Fp1: p = 1.1e-13, Hedges'g = 1.5, Fp2: p = 1.7e-5, Hedges'g = 1.26). Adjustment factor to have identical 1 V/m EF strength in a 5mm sphere around the individualized targets ranged from 0.72 to 2.3 (mean ± SD: 1.07 ± 0.29).

Conclusion

Our results show that optimizing coil orientation and stimulation intensity based on individualized TMS targets led to stronger electric fields in the targeted brain regions compared to a one-size-fits-all approach. These findings provide valuable insights for refining TMS therapy for SUDs by optimizing the modulation of cortico-subcortical circuits.

A Systematic Review and Large-Scale tES and TMS Electric Field Modeling Study Reveals How Outcome Measure Selection Alters Results in a Person- and Montage-Specific Manner

Background

Electric field (E-field) modeling is a potent tool to examine the cortical effects of transcranial magnetic and electrical stimulation (TMS and tES, respectively) and to address the high variability in efficacy observed in the literature. However, outcome measures used to report E-field magnitude vary considerably and have not yet been compared in detail.

Objectives

The goal of this two-part study, encompassing a systematic review and modeling experiment, was to provide an overview of the different outcome measures used to report the magnitude of tES and TMS E-fields, and to conduct a direct comparison of these measures across different stimulation montages.

Methods

Three electronic databases were searched for tES and/or TMS studies reporting E-field magnitude. We extracted and discussed outcome measures in studies meeting the inclusion criteria. Additionally, outcome measures were compared via models of four common tES and two TMS modalities in 100 healthy younger adults.

Results

In the systematic review, we included 118 studies using 151 outcome measures related to E-field magnitude. Structural and spherical regions of interest (ROI) analyses and percentile-based whole-brain analyses were used most often. In the modeling analyses, we found that there was an average of only 6% overlap between ROI and percentile-based whole-brain analyses in the investigated volumes within the same person. The overlap between ROI and whole-brain percentiles was montage- and person-specific, with more focal montages such as 4Ã-1 and APPS-tES, and figure-of-eight TMS showing up to 73%, 60%, and 52% overlap between ROI and percentile approaches respectively. However, even in these cases, 27% or more of the analyzed volume still differed between outcome measures in every analyses.

Conclusions

The choice of outcome measures meaningfully alters the interpretation of tES and TMS E-field models. Well-considered outcome measure selection is imperative for accurate interpretation of results, valid between-study comparisons, and depends on stimulation focality and study goals. We formulated four recommendations to increase the quality and rigor of E-field modeling outcome measures. With these data and recommendations, we hope to guide future studies towards informed outcome measure selection, and improve the comparability of studies.

Fast And Accurate Population Level Transcranial Magnetic Stimulation via Low-Rank Probabilistic Matrix Decomposition (PMD)

Transcranial magnetic stimulation (TMS) is used to study brain function and treat mental health disorders. During TMS, a coil placed on the scalp induces an E-field in the brain that modulates its activity. TMS is known to stimulate regions that are exposed to a large E-field. Clinical TMS protocols prescribe a coil placement based on scalp landmarks. There are inter-individual variations in brain anatomy that result in variations in the TMS-induced E-field at the targeted region and its outcome. These variations across individuals could in principle be minimized by developing a large database of head subjects and determining scalp landmarks that maximize E-field at the targeted brain region while minimizing its variation using computational methods. However, this approach requires repeated execution of a computational method to determine the E-field induced in the brain for a large number of subjects and coil placements. We developed a probabilistic matrix decomposition-based approach for rapidly evaluating the E-field induced during TMS for a large number of coil placements. Our approach can determine the E-field induced in over 1 Million coil placements in 9.5 hours, in contrast, to over 5 years using a bruteforce approach. After the initial set-up stage, the E-field can be predicted over the whole brain within 2-3 milliseconds and to 2% accuracy. We tested our approach in over 200 subjects and achieved an error of < 2% in most and < 3.5% in all subjects. We will present several examples of bench-marking analysis for our tool in terms of accuracy and speed across and its applicability for population level optimization of coil placement.

Highlights

A method for practical E-field informed population-level TMS coil placement strategies is developed.This algorithm enables the determination of E-field informed optimal coil placement in seconds enabling it’s use for close-loop and on-the-fly reconfiguration of TMS.After the initial set-up stage of less than 10 hours, the E-field can be predicted for any coil placement across the whole brain within in 2-3 milliseconds.

Responses of model cortical neurons to temporal interference stimulation and related transcranial alternating current stimulation modalities

Objective.Temporal interference stimulation (TIS) was proposed as a non-invasive, focal, and steerable deep brain stimulation method. However, the mechanisms underlying experimentally-observed suprathreshold TIS effects are unknown, and prior simulation studies had limitations in the representations of the TIS electric field (E-field) and cerebral neurons. We examined the E-field and neural response characteristics for TIS and related transcranial alternating current stimulation modalities.Approach.Using the uniform-field approximation, we simulated a range of stimulation parameters in biophysically realistic model cortical neurons, including different orientations, frequencies, amplitude ratios, amplitude modulation, and phase difference of the E-fields, and obtained thresholds for both activation and conduction block.Main results. For two E-fields with similar amplitudes (representative of E-field distributions at the target region), TIS generated an amplitude-modulated (AM) total E-field. Due to the phase difference of the individual E-fields, the total TIS E-field vector also exhibited rotation where the orientations of the two E-fields were not aligned (generally also at the target region). TIS activation thresholds (75-230 V m-1) were similar to those of high-frequency stimulation with or without modulation and/or rotation. For E-field dominated by the high-frequency carrier and with minimal amplitude modulation and/or rotation (typically outside the target region), TIS was less effective at activation and more effective at block. Unlike AM high-frequency stimulation, TIS generated conduction block with some orientations and amplitude ratios of individual E-fields at very high amplitudes of the total E-field (>1700 V m-1).Significance. The complex 3D properties of the TIS E-fields should be accounted for in computational and experimental studies. The mechanisms of suprathreshold cortical TIS appear to involve neural activity block and periodic activation or onset response, consistent with computational studies of peripheral axons. These phenomena occur at E-field strengths too high to be delivered tolerably through scalp electrodes and may inhibit endogenous activity in off-target regions, suggesting limited significance of suprathreshold TIS.

The neurobiology of misophonia and implications for novel, neuroscience-driven interventions

Decreased tolerance in response to specific every-day sounds (misophonia) is a serious, debilitating disorder that is gaining rapid recognition within the mental health community. Emerging research findings suggest that misophonia may have a unique neural signature. Specifically, when examining responses to misophonic trigger sounds, differences emerge at a physiological and neural level from potentially overlapping psychopathologies. While these findings are preliminary and in need of replication, they support the hypothesis that misophonia is a unique disorder. In this theoretical paper, we begin by reviewing the candidate networks that may be at play in this complex disorder (e.g., regulatory, sensory, and auditory). We then summarize current neuroimaging findings in misophonia and present areas of overlap and divergence from other mental health disorders that are hypothesized to co-occur with misophonia (e.g., obsessive compulsive disorder). Future studies needed to further our understanding of the neuroscience of misophonia will also be discussed. Next, we introduce the potential of neurostimulation as a tool to treat neural dysfunction in misophonia. We describe how neurostimulation research has led to novel interventions in psychiatric disorders, targeting regions that may also be relevant to misophonia. The paper is concluded by presenting several options for how neurostimulation interventions for misophonia could be crafted.

Dosing Transcranial Magnetic Stimulation of the Primary Motor and Dorsolateral Prefrontal Cortices With Multi-Scale Modeling

Transcranial magnetic stimulation (TMS) can depolarize cortical neurons through the intact skin and skull. The characteristics of the induced electric field (E-field) have a major impact on specific outcomes of TMS. Using multi-scale computational modeling, we explored whether the stimulation parameters derived from the primary motor cortex (M1) induce comparable macroscopic E-field strengths and subcellular/cellular responses in the dorsolateral prefrontal cortex (DLPFC). To this aim, we calculated the TMS-induced E-field in 16 anatomically realistic head models and simulated the changes in membrane voltage and intracellular calcium levels of morphologically and biophysically realistic human pyramidal cells in the M1 and DLPFC. We found that the conventional intensity selection methods (i.e., motor threshold and fixed intensities) produce variable macroscopic E-fields. Consequently, it was challenging to produce comparable subcellular/cellular responses across cortical regions with distinct folding characteristics. Prospectively, personalized stimulation intensity selection could standardize the E-fields and the subcellular/cellular responses to repetitive TMS across cortical regions and individuals. The suggested computational approach points to the shortcomings of the conventional intensity selection methods used in clinical settings. We propose that multi-scale modeling has the potential to overcome some of these limitations and broaden our understanding of the neuronal mechanisms for TMS.

Effect of neuroanatomy on corticomotor excitability during and after transcranial magnetic stimulation and intermittent theta burst stimulation

Individual neuroanatomy can influence motor responses to transcranial magnetic stimulation (TMS) and corticomotor excitability after intermittent theta burst stimulation (iTBS). The purpose of this study was to examine the relationship between individual neuroanatomy and both TMS response measured using resting motor threshold (RMT) and iTBS measured using motor evoked potentials (MEPs) targeting the biceps brachii and first dorsal interosseus (FDI). Ten nonimpaired individuals completed sham‐controlled iTBS sessions and underwent MRI, from which anatomically accurate head models were generated. Neuroanatomical parameters established through fiber tractography were fiber tract surface area (FTSA), tract fiber count (TFC), and brain scalp distance (BSD) at the point of stimulation. Cortical magnetic field induced electric field strength (EFS) was obtained using finite element simulations. A linear mixed effects model was used to assess effects of these parameters on RMT and iTBS (post‐iTBS MEPs). FDI RMT was dependent on interactions between EFS and both FTSA and TFC. Biceps RMT was dependent on interactions between EFS and and both FTSA and BSD. There was no groupwide effect of iTBS on the FDI but individual changes in corticomotor excitability scaled with RMT, EFS, BSD, and FTSA. iTBS targeting the biceps was facilitatory, and dependent on FTSA and TFC. MRI‐based measures of neuroanatomy highlight how individual anatomy affects motor system responses to different TMS paradigms and may be useful for selecting appropriate motor targets when designing TMS based therapies.

Electric Field Distribution Induced by TMS: Differences Due to Anatomical Variation

Transcranial magnetic stimulation (TMS) is a well-established technique for the diagnosis and treatment of neuropsychiatric diseases. The numerical calculation of the induced electric field (EF) distribution in the brain increases the efficacy of stimulation and improves clinical outcomes. However, unique anatomical features, which distinguish each subject, suggest that personalized models should be preferentially used. The objective of the present study was to assess how anatomy affects the EF distribution and to determine to what extent personalized models are useful for clinical studies. The head models of nineteen healthy volunteers were automatically segmented. Two versions of each head model, a homogeneous and a five-tissue anatomical, were stimulated by the model of a Hesed coil (H-coil), employing magnetic quasi-static simulations. The H-coil was placed at two standard stimulating positions per model, over the frontal and central areas. The results show small, but indisputable, variations in the EFs for the homogeneous and anatomical models. The interquartile ranges in the anatomical versions were higher compared to the homogeneous ones, indicating that individual anatomical features may affect the prediction of stimulation volumes. It is concluded that personalized models provide complementary information and should be preferably employed in the context of diagnostic and therapeutic TMS studies.

The Effect of Inter-pulse Interval on TMS Motor Evoked Potentials in Active Muscles

Objective

The time interval between transcranial magnetic stimulation (TMS) pulses affects evoked muscle responses when the targeted muscle is resting. This necessitates using sufficiently long inter-pulse intervals (IPIs). However, there is some evidence that the IPI has no effect on the responses evoked in active muscles. Thus, we tested whether voluntary contraction could remove the effect of the IPI on TMS motor evoked potentials (MEPs).

Methods

In our study, we delivered sets of 30 TMS pulses with three different IPIs (2, 5, and 10 s) to the left primary motor cortex. These measurements were performed with the resting and active right hand first dorsal interosseous muscle in healthy participants (N = 9 and N = 10). MEP amplitudes were recorded through electromyography.

Results

We found that the IPI had no significant effect on the MEP amplitudes in the active muscle (p = 0.36), whereas in the resting muscle, the IPI significantly affected the MEP amplitudes (p < 0.001), decreasing the MEP amplitude of the 2 s IPI.

Conclusions

These results show that active muscle contraction removes the effect of the IPI on the MEP amplitude. Therefore, using active muscles in TMS motor mapping enables faster delivery of TMS pulses, reducing measurement time in novel TMS motor mapping studies.

TMS intensity and focality correlation with coil orientation at three non-motor regions

Objective. The aim of this study is to define the best coil orientations for transcranial magnetic stimulation (TMS) for three clinically relevant brain areas: pre-supplementary motor area (pre-SMA), inferior frontal gyrus (IFG), and posterior parietal cortex (PPC), by means of simulations in 12 realistic head models of the electric field (E-field). Methods. We computed the E-field generated by TMS in our three volumes of interest (VOI) that were delineated based on published atlases. We then analysed the maximum intensity and spatial focality for the normal and absolute components of the E-field considering different percentile thresholds. Lastly, we correlated these results with the different anatomical properties of our VOIs. Results. Overall, the spatial focality of the E-field for the three VOIs varied depending on the orientation of the coil. Further analysis showed that differences in individual brain anatomy were related to the amount of focality achieved. In general, a larger percentage of sulcus resulted in better spatial focality. Additionally, a higher normal E-field intensity was achieved when the coil axis was placed perpendicular to the predominant orientations of the gyri of each VOI. A positive correlation between spatial focality and E-field intensity was found for PPC and IFG but not for pre-SMA. Conclusions. For a rough approximation, better coil orientations can be based on the individual’s specific brain morphology at the VOI. Moreover, TMS computational models should be employed to obtain better coil orientations in non-motor regions of interest. Significance. Finding better coil orientations in non-motor regions is a challenge in TMS and seeks to reduce interindividual variability. Our individualized TMS simulation pipeline leads to fewer inter-individual variability in the focality, likely enhancing the efficacy of the stimulation and reducing the risk of stimulating adjacent, non-targeted areas.

TMS with fast and accurate electronic control: Measuring the orientation sensitivity of corticomotor pathways

BACKGROUND

Transcranial magnetic stimulation (TMS) coils allow only a slow, mechanical adjustment of the stimulating electric field (E-field) orientation in the cerebral tissue. Fast E-field control is needed to synchronize the stimulation with the ongoing brain activity. Also, empirical models that fully describe the relationship between evoked responses and the stimulus orientation and intensity are still missing.

OBJECTIVE

We aimed to (1) develop a TMS transducer for manipulating the E-field orientation electronically with high accuracy at the neuronally meaningful millisecond-level time scale and (2) devise and validate a physiologically based model describing the orientation selectivity of neuronal excitability.

METHODS

We designed and manufactured a two-coil TMS transducer. The coil windings were computed with a minimum-energy optimization procedure, and the transducer was controlled with our custom-made electronics. The electronic E-field control was verified with a TMS characterizer. The motor evoked potential amplitude and latency of a hand muscle were mapped in 3° steps of the stimulus orientation in 16 healthy subjects for three stimulation intensities. We fitted a logistic model to the motor response amplitude.

RESULTS

The two-coil TMS transducer allows one to manipulate the pulse orientation accurately without manual coil movement. The motor response amplitude followed a logistic function of the stimulus orientation; this dependency was strongly affected by the stimulus intensity.

CONCLUSION

The developed electronic control of the E-field orientation allows exploring new stimulation paradigms and probing neuronal mechanisms. The presented model helps to disentangle the neuronal mechanisms of brain function and guide future non-invasive stimulation protocols.

Precise Modulation Strategies for Transcranial Magnetic Stimulation: Advances and Future Directions

Transcranial magnetic stimulation (TMS) is a popular modulatory technique for the noninvasive diagnosis and therapy of neurological and psychiatric diseases. Unfortunately, current modulation strategies are only modestly effective. The literature provides strong evidence that the modulatory effects of TMS vary depending on device components and stimulation protocols. These differential effects are important when designing precise modulatory strategies for clinical or research applications. Developments in TMS have been accompanied by advances in combining TMS with neuroimaging techniques, including electroencephalography, functional near-infrared spectroscopy, functional magnetic resonance imaging, and positron emission tomography. Such studies appear particularly promising as they may not only allow us to probe affected brain areas during TMS but also seem to predict underlying research directions that may enable us to precisely target and remodel impaired cortices or circuits. However, few precise modulation strategies are available, and the long-term safety and efficacy of these strategies need to be confirmed. Here, we review the literature on possible technologies for precise modulation to highlight progress along with limitations with the goal of suggesting future directions for this field.

A probabilistic transcranial magnetic stimulation localization method

Objective.Transcranial magnetic stimulation (TMS) can be used to safely and noninvasively activate brain tissue. However, the characteristic parameters of the neuronal activation have been largely unclear. In this work, we propose a novel neuronal activation model and develop a method to infer its parameters from measured motor evoked potential signals.Approach.The connection between neuronal activation due to an induced electric field and a measured motor threshold is modeled. The posterior distribution of the model parameters are inferred from measurement data using Bayes' formula. The measurements are the active motor thresholds obtained with multiple stimulating coil locations, and the parameters of the model are the location, preferred direction of activation, and threshold electric field value of the activation site. The posterior distribution is sampled using a Markov chain Monte Carlo method. We quantify the plausibility of the model by calculating the marginal likelihood of the measured thresholds. The method is validated with synthetic data and applied to motor threshold measurements from the first dorsal interosseus muscle in five healthy participants.Main results.The method produces a probability distribution for the activation location, from which a minimal volume where the activation occurs with 95% probability can be derived. For eight or nine stimulating coil locations, the smallest such a volume obtained was approximately 100 mm³. The 95% probability volume intersected the pre-central gyral crown and the anterior wall of the central sulcus, and the preferred direction was perpendicular to the central sulcus, both findings being consistent with the literature. Furthermore, it was not possible to rule out if the activation occurred either in the white or grey matter. In one participant, two distinct activations sites were found while others exhibited a unique site.Significance.The method is both generic and robust, and it lays a foundation for a framework that enables accurate analysis and characterization of TMS activation mechanisms.

Estimation of the Motor Threshold for Near-Rectangular Stimuli Using the Hodgkin-Huxley Model

The motor threshold measurement is a standard in preintervention probing in TMS experiments. We aim to predict the motor threshold for near-rectangular stimuli to efficiently determine the motor threshold size before any experiments take place. Estimating the behavior of large-scale networks requires dynamically accurate and efficient modeling. We utilized a Hodgkin-Huxley (HH) type model to evaluate motor threshold values and computationally validated its function with known true threshold data from 50 participants trials from state-of-the-art published datasets. For monophasic, bidirectional, and unidirectional rectangular stimuli in posterior-anterior or anterior-posterior directions as generated by the cTMS device, computational modeling of the HH model captured the experimentally measured population-averaged motor threshold values at high precision (maximum error ≤ 8%). The convergence of our biophysically based modeling study with experimental data in humans reveals that the effect of the stimulus shape is strongly correlated with the activation kinetics of the voltage-gated ion channels. The proposed method can reliably predict motor threshold size using the conductance-based neuronal models and could therefore be embedded in new generation neurostimulators. Advancements in neural modeling will make it possible to enhance treatment procedures by reducing the number of delivered magnetic stimuli to participants.

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.

Understanding intracortical excitability in phantom limb pain: A multivariate analysis from a multicenter randomized clinical trial

OBJECTIVES

To explore associations of intracortical excitability with clinical characteristics in a large sample of subjects with phantom limb pain (PLP).

METHODS

Ancillary study using baseline and longitudinal data from a large multicenter randomized trial that investigated the effects of non-invasive brain stimulation combined with sensorimotor training on PLP. Multivariate regression modeling analyses were used to investigate the association of intracortical excitability, measured by percentages of intracortical inhibition (ICI) and facilitation (ICF) with clinical variables.

RESULTS

Ninety-eight subjects were included. Phantom sensation of itching was positively associated with ICI changes and at baseline in the affected hemisphere (contralateral to PLP). However, in the non-affected hemisphere (ipsilateral to PLP), the phantom sensation of warmth and PLP intensity were negatively associated with ICI (both models). For the ICF, PLP intensity (baseline model only) and age (longitudinal model) were negatively associated, while time since amputation and amputation level (both for longitudinal model only) were positively associated in the affected hemisphere. Additionally, use of antidepressants led to lower ICF in the non-affected hemisphere for the baseline model while higher amputation level also led to less changes in the ICF.

CONCLUSION

Results revealed clear associations of clinical variables and cortical excitability in a large chronic pain sample. ICI and ICF changes appear not to be mainly explained by PLP intensity. Instead, other variables associated with duration of neuroplasticity changes (such as age and duration of amputation) and compensatory mechanisms (such as itching and phantom limb sensation) seem to be more important in explaining these variables.

Characterizing off-target corticospinal responses to double-cone transcranial magnetic stimulation

INTRODUCTION

The double-cone coil (D-CONE) is frequently used in transcranial magnetic stimulation (TMS) experiments that target the motor cortex (M1) lower-limb representation. Anecdotal evidence and modeling studies have shed light on the off-target effects of D-CONE TMS but the physiological extent remains undetermined.

PURPOSE

To characterize the off-target effects of D-CONE TMS based on bilateral corticospinal responses in the legs and hands.

METHODS

Thirty (N = 30) participants (9 women, age: 26 ± 5yrs) completed a stimulus-response curve procedure with D-CONE TMS applied to the dominant vastus lateralis (cVL) and motor-evoked potentials (MEPs) recorded in each active VL and resting first dorsal interosseous (FDI). As a positive control (CON), the dominant FDI was directly targeted with a figure-of-eight coil and MEPs were similarly recorded in each active FDI and resting VL. MEP_MAX, V50 and MEP latencies were compared with repeated-measures ANOVAs or mixed-effects analysis and Bonferroni-corrected pairwise comparisons.

RESULTS

Off-target responses were evident in all muscles, with similar MEP_MAX in the target (cVL) and off-target (iVL) leg (p = 0.99) and cFDI compared with CON (p = 0.99). cFDI and CON MEP_MAX were greater than iFDI (p < 0.01). A main effect of target (p < 0.001) indicated that latencies were shorter with CON but similar in all muscles with D-CONE.

DISCUSSION

Concurrent MEP recordings in bilateral upper- and lower-extremity muscles confirm that lower-limb D-CONE TMS produces substantial distance-dependent off-target effects. In addition to monitoring corticospinal responses in off-target muscles to improve targeting accuracy in real-time, future studies may incorporate off-target information into statistical models post-hoc.

Effect of skin conductivity on the electric field induced by transcranial stimulation techniques in different head models

This study aims at quantifying the effect that using different skin conductivity values has on the estimation of the electric (E)-field distribution induced by transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS) in the brain of two anatomical models. The induced E-field was calculated with numerical simulations inside MIDA and Duke models, assigning to the skin a conductivity value estimated from a multi-layered skin model and three values taken from literature. The effect of skin conductivity variations on the local E-field induced by tDCS in the brain was up to 70%. In TMS, minor local differences, in the order of 20%, were obtained in regions of interest for the onset of possible side effects. Results suggested that an accurate model of the skin is necessary in all numerical studies that aim at precisely estimating the E-field induced during TMS and tDCS applications. This also highlights the importance of further experimental studies on human skin characterization, especially at low frequencies.

Evaluating transcranial magnetic stimulation (TMS) induced electric fields in pediatric stroke

•

Numerical TMS simulations were performed over and around perinatal stroke lesions.

•

The presence of brain lesions locally affects the electric field distribution.

•

Brain lesions do not significantly change the mean electric field strength.

•

Model driven approaches can inform TMS dosing in a pediatric stroke population.

Transcranial magnetic stimulation (TMS) is an increasingly popular tool for stroke rehabilitation. Consequently, researchers have started to explore the use of TMS in pediatric stroke. However, the application of TMS in a developing brain with pathologies comes with a unique set of challenges. The effect of TMS-induced electric fields has not been explored in children with stroke lesions. Here, we used finite element method (FEM) modeling to study how the electric field strength is affected by the presence of a lesion. We created individual realistic head models from MRIs (n = 6) of children with unilateral cerebral palsy due to perinatal stroke. We conducted TMS electric field simulations for coil locations over lesioned and non-lesioned hemispheres. We found that the presence of a lesion can strongly affect the electric field distribution. On the group level, the mean electric field strength did not differ between lesioned and non-lesioned hemispheres but exhibited a greater variability in the lesioned hemisphere. Other factors such as coil-to-cortex distance have a strong influence on the TMS electric field even in the presence of lesions. Our study has important implications for the delivery of TMS in children with brain lesions with respect to TMS dosing and coil placement.

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.

A software toolkit for TMS electric-field modeling with boundary element fast multipole method: an efficient MATLAB implementation

Objective

To present and disseminate our transcranial magnetic stimulation (TMS) modeling software toolkit, including several new algorithmic developments, and to apply this software to realistic TMS modeling scenarios given a high-resolution model of the human head including cortical geometry and an accurate coil model.

Approach

The recently developed charge-based boundary element fast multipole method (BEM-FMM) is employed as an alternative to the 1st order finite element method (FEM) most commonly used today. The BEM-FMM approach provides high accuracy and unconstrained numerical field resolution close to and across cortical interfaces. Here, the previously proposed BEM-FMM algorithm has been improved in several novel ways.

Main results

The improvements resulted in a threefold increase in computational speed while maintaining the same solution accuracy. The computational code based on the MATLAB® platform is made available to all interested researchers, along with a coil model repository and examples to create custom coils, head model repository, and supporting documentation. The presented software toolkit may be useful for post-hoc analyses of navigated TMS data using high-resolution subject-specific head models as well as accurate and fast modeling for the purposes of TMS coil/hardware development.

Significance

TMS is currently the only non-invasive neurostimulation modality that enables painless and safe supra-threshold stimulation by employing electromagnetic induction to efficiently penetrate the skull. Accurate, fast, and high resolution modeling of the electric fields may significantly improve individualized targeting and dosing of TMS and therefore enhance the efficiency of existing clinical protocols as well as help establish new application domains.

Weak rTMS-induced electric fields produce neural entrainment in humans

Repetitive transcranial magnetic stimulation (rTMS) is a potent tool for modulating endogenous oscillations in humans. The current standard method for rTMS defines the stimulation intensity based on the evoked liminal response in the visual or motor system (e.g., resting motor threshold). The key limitation of the current approach is that the magnitude of the resulting electric field remains elusive. A better characterization of the electric field strength induced by a given rTMS protocol is necessary in order to improve the understanding of the neural mechanisms of rTMS. In this study we used a novel approach, in which individualized prospective computational modeling of the induced electric field guided the choice of stimulation intensity. We consistently found that rhythmic rTMS protocols increased neural synchronization in the posterior alpha frequency band when measured simultaneously with scalp electroencephalography. We observed this effect already at electric field strengths of roughly half the lowest conventional field strength, which is 80% of the resting motor threshold. We conclude that rTMS can induce immediate electrophysiological effects at much weaker electric field strengths than previously thought.

Structural and functional motor cortex asymmetry in unilateral lower limb amputation with phantom limb pain

OBJECTIVE

The role of motor cortex reorganization in the development and maintenance of phantom limb pain (PLP) is still unclear. This study aims to evaluate neurophysiological and structural motor cortex asymmetry in patients with PLP and its relationship with pain intensity.

METHODS

Cross-sectional analysis of an ongoing randomized-controlled trial. We evaluated the motor cortex asymmetry through two techniques: i) changes in cortical excitability indexed by transcranial magnetic stimulation (motor evoked potential, paired-pulse paradigms and cortical mapping), and ii) voxel-wise grey matter asymmetry analysis by brain magnetic resonance imaging.

RESULTS

We included 62 unilateral traumatic lower limb amputees with a mean PLP of 5.9 (SD = 1.79). We found, in the affected hemisphere, an anterior shift of the hand area center of gravity (23 mm, 95% CI 6 to 38, p = 0.005) and a disorganized and widespread representation. Regarding voxel-wise grey matter asymmetry analysis, data from 21 participants show a loss of grey matter volume in the motor area of the affected hemisphere. This asymmetry seems negatively associated with time since amputation. For TMS data, only the ICF ratio is negatively correlated with PLP intensity (r = -0.25, p = 0.04).

CONCLUSION

There is an asymmetrical reorganization of the motor cortex in patients with PLP, characterized by a disorganized, widespread, and shifted hand cortical representation and a loss in grey matter volume in the affected hemisphere. This reorganization seems to reduce across time since amputation. However, it is not associated with pain intensity.

SIGNIFICANCE

These findings are significant to understand the role of the motor cortex reorganization in patients with PLP, showing that the pain intensity may be related with other neurophysiological factors, not just cortical reorganization.

Group-level analysis of induced electric field in deep brain regions by different TMS coils

Deep transcranial magnetic stimulation (dTMS) is a non-invasive technique used for the treatment of depression and obsessive compulsive disorder. In this study, we computationally evaluated group-level dosage for dTMS to characterize the targeted deep brain regions to overcome the limitations of using individualized head models to characterize coil performance in a population. We used an inter-subject registration method adapted to the deep brain regions that enable projection of computed electric fields (EFs) from individual realistic head models (n  =  18) to the average space of deep brain regions. The computational results showed consistent group-level hotspots of the EF in the deep brain regions. The halo circular assembly coils induced the highest EFs in deep brain regions (up to 50% of the maximum EF in the cortex) for optimized positioning. In terms of the trade-off between field spread and penetration, the performance of the H7 coil was the best. The computational model allowed the optimization of generalized dTMS-induced EF on deep region targets despite inter-individual differences while informing and possibly minimizing unintended stimulation of superficial regions and possible mixed stimulation effects from deep and cortical areas. These results will facilitate the decision process during dTMS interventions in clinical practice.

Electric field simulations for transcranial brain stimulation using FEM: an efficient implementation and error analysis

OBJECTIVE

Transcranial magnetic stimulation (TMS) and transcranial electric stimulation (TES) modulate brain activity non-invasively by generating electric fields either by electromagnetic induction or by injecting currents via skin electrodes. Numerical simulations based on anatomically detailed head models of the TMS and TES electric fields can help us to understand and optimize the spatial stimulation pattern in the brain. However, most realistic simulations are still slow, and the role of anatomical fidelity on simulation accuracy has not been evaluated in detail so far.

APPROACH

We present and validate a new implementation of the finite element method (FEM) for TMS and TES that is based on modern algorithms and libraries. We also evaluate the convergence of the simulations and estimate errors stemming from numerical and modelling aspects.

MAIN RESULTS

Comparisons with analytical solutions for spherical phantoms validate our new FEM implementation, which is three to six times faster than previous implementations. The convergence results suggest that accurately capturing the tissue geometry in addition to choosing a sufficiently accurate numerical method is of fundamental importance for accurate simulations.

SIGNIFICANCE

The new implementation allows for a substantial increase in computational efficiency of FEM TMS and TES simulations. This is especially relevant for applications such as the systematic assessment of model uncertainty and the optimization of multi-electrode TES montages. The results of our systematic error analysis allow the user to select the best tradeoff between model resolution and simulation speed for a specific application. The new FEM code is openly available as a part of our open-source software SimNIBS 3.0.

MRI-based visualization of rTMS-induced cortical plasticity in the primary motor cortex

Repetitive transcranial magnetic stimulation (rTMS) induces changes in cortical excitability for minutes to hours after the end of intervention. However, it has not been precisely determined to what extent cortical plasticity prevails spatially in the cortex. Recent studies have shown that rTMS induces changes in “interhemispheric” functional connectivity, the resting-state functional connectivity between the stimulated region and the symmetrically corresponding region in the contralateral hemisphere. In the present study, quadripulse stimulation (QPS) was applied to the index finger representation in the left primary motor cortex (M1), while the position of the stimulation coil was constantly monitored by an online navigator. After QPS application, resting-state functional magnetic resonance imaging was performed, and the interhemispheric functional connectivity was compared with that before QPS. A cluster of connectivity changes was observed in the stimulated region in the central sulcus. The cluster was spatially extended approximately 10 mm from the center [half width at half maximum (HWHM): approximately 3 mm] and was extended approximately 20 mm long in depth (HWHM: approximately 7 mm). A localizer scan of the index finger motion confirmed that the cluster of interhemispheric connectivity changes overlapped spatially with the activation related to the index finger motion. These results indicate that cortical plasticity in M1 induced by rTMS was relatively restricted in space and suggest that rTMS can reveal functional dissociation associated with adjacent small areas by inducing neural plasticity in restricted cortical regions.

Simulation of transcranial magnetic stimulation in head model with morphologically-realistic cortical neurons

BACKGROUND

Transcranial magnetic stimulation (TMS) enables non-invasive modulation of brain activity with both clinical and research applications, but fundamental questions remain about the neural types and elements TMS activates and how stimulation parameters affect the neural response.

OBJECTIVE

To develop a multi-scale computational model to quantify the effect of TMS parameters on the direct response of individual neurons.

METHODS

We integrated morphologically-realistic neuronal models with TMS-induced electric fields computed in a finite element model of a human head to quantify the cortical response to TMS with several combinations of pulse waveforms and current directions.

RESULTS

TMS activated with lowest intensity intracortical axonal terminations in the superficial gyral crown and lip regions. Layer 5 pyramidal cells had the lowest thresholds, but layer 2/3 pyramidal cells and inhibitory basket cells were also activated at most intensities. Direct activation of layers 1 and 6 was unlikely. Neural activation was largely driven by the field magnitude, rather than the field component normal to the cortical surface. Varying the induced current direction caused a waveform-dependent shift in the activation site and provided a potential mechanism for experimentally observed differences in thresholds and latencies of muscle responses.

CONCLUSIONS

This biophysically-based simulation provides a novel method to elucidate mechanisms and inform parameter selection of TMS and other cortical stimulation modalities. It also serves as a foundation for more detailed network models of the response to TMS, which may include endogenous activity, synaptic connectivity, inputs from intrinsic and extrinsic axonal projections, and corticofugal axons in white matter.

Conditions for numerically accurate TMS electric field simulation

BACKGROUND

Computational simulations of the E-field induced by transcranial magnetic stimulation (TMS) are increasingly used to understand its mechanisms and to inform its administration. However, characterization of the accuracy of the simulation methods and the factors that affect it is lacking.

OBJECTIVE

To ensure the accuracy of TMS E-field simulations, we systematically quantify their numerical error and provide guidelines for their setup.

METHOD

We benchmark the accuracy of computational approaches that are commonly used for TMS E-field simulations, including the finite element method (FEM) with and without superconvergent patch recovery (SPR), boundary element method (BEM), finite difference method (FDM), and coil modeling methods.

RESULTS

To achieve cortical E-field error levels below 2%, the commonly used FDM and 1st order FEM require meshes with an average edge length below 0.4 mm, 1st order SPR-FEM requires edge lengths below 0.8 mm, and BEM and 2nd (or higher) order FEM require edge lengths below 2.9 mm. Coil models employing magnetic and current dipoles require at least 200 and 3000 dipoles, respectively. For thick solid-conductor coils and frequencies above 3 kHz, winding eddy currents may have to be modeled.

CONCLUSION

BEM, FDM, and FEM all converge to the same solution. Compared to the common FDM and 1st order FEM approaches, BEM and 2nd (or higher) order FEM require significantly lower mesh densities to achieve the same error level. In some cases, coil winding eddy-currents must be modeled. Both electric current dipole and magnetic dipole models of the coil current can be accurate with sufficiently fine discretization.

Determinants of Inter-Individual Variability in Corticomotor Excitability Induced by Paired Associative Stimulation

Transcranial magnetic stimulation (TMS) is a well-established tool in probing cortical plasticity in vivo. Changes in corticomotor excitability can be induced using paired associative stimulation (PAS) protocol, in which TMS over the primary motor cortex is conditioned with an electrical peripheral nerve stimulation of the contralateral hand. PAS with an inter-stimulus interval of 25 ms induces long-term potentiation (LTP)-like effects in cortical excitability. However, the response to a PAS protocol tends to vary substantially across individuals. In this study, we used univariate and multivariate data-driven methods to investigate various previously proposed determinants of inter-individual variability in PAS efficacy, such as demographic, cognitive, clinical, neurophysiological, and neuroimaging measures. Forty-one right-handed participants, comprising 22 patients with amnestic mild cognitive impairment (MCI) and 19 healthy controls (HC), underwent the PAS protocol. Prior to stimulation, demographic, genetic, clinical, as well as structural and resting-state functional MRI data were acquired. The two groups did not differ in any of the variables, except by global cognitive status. Univariate analysis showed that only 61% of all participants were classified as PAS responders, irrespective of group membership. Higher PAS response was associated with lower TMS intensity and with higher resting-state connectivity within the sensorimotor network, but only in responders, as opposed to non-responders. We also found an overall positive correlation between PAS response and structural connectivity within the corticospinal tract, which did not differ between groups. A multivariate random forest (RF) model identified age, gender, education, IQ, global cognitive status, sleep quality, alertness, TMS intensity, genetic factors, and neuroimaging measures (functional and structural connectivity, gray matter (GM) volume, and cortical thickness as poor predictors of PAS response. The model resulted in low accuracy of the RF classifier (58%; 95% CI: 42 - 74%), with a higher relative importance of brain connectivity measures compared to the other variables. We conclude that PAS variability in our sample was not well explained by factors known to influence PAS efficacy, emphasizing the need for future replication studies.

Online repetitive transcranial magnetic stimulation during working memory in younger and older adults: A randomized within-subject comparison

Working memory is the ability to perform mental operations on information that is stored in a flexible, limited capacity buffer. The ability to manipulate information in working memory is central to many aspects of human cognition, but also declines with healthy aging. Given the profound importance of such working memory manipulation abilities, there is a concerted effort towards developing approaches to improve them. The current study tested the capacity to enhance working memory manipulation with online repetitive transcranial magnetic stimulation in healthy young and older adults. Online high frequency (5Hz) repetitive transcranial magnetic stimulation was applied over the left dorsolateral prefrontal cortex to test the hypothesis that active repetitive transcranial magnetic stimulation would lead to significant improvements in memory recall accuracy compared to sham stimulation, and that these effects would be most pronounced in working memory manipulation conditions with the highest cognitive demand in both young and older adults. Repetitive transcranial magnetic stimulation was applied while participants were performing a delayed response alphabetization task with three individually-titrated levels of difficulty. The left dorsolateral prefrontal cortex was identified by combining electric field modeling to individualized functional magnetic resonance imaging activation maps and was targeted during the experiment using stereotactic neuronavigation with real-time robotic guidance, allowing optimal coil placement during the stimulation. As no accuracy differences were found between young and older adults, the results from both groups were collapsed. Subsequent analyses revealed that active stimulation significantly increased accuracy relative to sham stimulation, but only for the hardest condition. These results point towards further investigation of repetitive transcranial magnetic stimulation for memory enhancement focusing on high difficulty conditions as those most likely to exhibit benefits.

Can Transcranial Electrical Stimulation Localize Brain Function?

Transcranial electrical stimulation (TES) uses constant (TDCS) or alternating currents (TACS) to modulate brain activity. Most TES studies apply low-intensity currents through scalp electrodes (≤2 mA) using bipolar electrode arrangements, producing weak electrical fields in the brain (<1 V/m). Low-intensity TES has been employed in humans to induce changes in task performance during or after stimulation. In analogy to focal transcranial magnetic stimulation, TES-induced behavioral effects have often been taken as evidence for a causal involvement of the brain region underlying one of the two stimulation electrodes, often referred to as the active electrode. Here, we critically review the utility of bipolar low-intensity TES to localize human brain function. We summarize physiological substrates that constitute peripheral targets for TES and may mediate subliminal or overtly perceived peripheral stimulation during TES. We argue that peripheral co-stimulation may contribute to the behavioral effects of TES and should be controlled for by "sham" TES. We discuss biophysical properties of TES, which need to be considered, if one wishes to make realistic assumptions about which brain regions were preferentially targeted by TES. Using results from electric field calculations, we evaluate the validity of different strategies that have been used for selective spatial targeting. Finally, we comment on the challenge of adjusting the dose of TES considering dose-response relationships between the weak tissue currents and the physiological effects in targeted cortical areas. These considerations call for caution when attributing behavioral effects during or after low-intensity TES studies to a specific brain region and may facilitate the selection of best practices for future TES studies.

Does Resting Motor Threshold Predict Motor Hand Recovery After Stroke?

Background: Resting Motor threshold (rMT) is one of the measurement obtained by Transcranial Magnetic Stimulation (TMS) that reflects corticospinal excitability. As a functional marker of the corticospinal pathway, the question arises whether rMT is a suitable biomarker for predicting post-stroke upper limb function. To that aim, we conducted a systematic review of relevant studies that investigated the clinical significance of rMT in stroke survivors by using correlations between upper limb motor scores and rMT. Methods: Studies that reported correlations between upper limb motor function and rMT as a measure of corticospinal excitability in distal arm muscle were identified via a literature search in stroke patients. Two authors extracted the data using a home-made specific form. Subgroup analyses were carried out with patients classified with respect to time post-stroke onset (early vs. chronic stage) and stroke location (cortical, subcortical, or cortico-subcortical). Methodological quality of the study was also evaluated by a published checklist. Results: Eighteen studies with 22 groups (n = 508 stroke patients) were included in this systematic review. Mean methodological quality score was 14.75/24. rMT was often correlated with motor function or hand dexterity (n = 15/22, 68%), explaining on average 31% of the variance of the motor score. Moreover, the results did not seem impacted if patients were examined at the early or chronic stages of stroke. Two findings could not be properly interpreted: (i) the fact that the rMT is an independent predictor of motor function as several confounding factors are well-established, and, (ii) whether the stroke location impacts this prediction. Conclusion: Most of the studies found a correlation between rMT and upper limb motor function after stroke. However, it is still unclear if rMT is an independent predictor of upper limb motor function when taking into account for age, time post stroke onset and level of corticospinal tract damage as confounding factors. Clear-cut conclusions could not be drawn at that time but our results suggest that rMT could be a suitable candidate although future investigations are needed. Systematic Review Registration Number: (https://www.crd.york.ac.uk/prospero/): ID 114317.

Phase Synchronicity of μ-Rhythm Determines Efficacy of Interhemispheric Communication Between Human Motor Cortices

The theory of communication through coherence predicts that effective connectivity between nodes in a distributed oscillating neuronal network depends on their instantaneous excitability state and phase synchronicity (Fries, 2005). Here, we tested this prediction by using state-dependent millisecond-resolved real-time electroencephalography-triggered dual-coil transcranial magnetic stimulation (EEG-TMS) (Zrenner et al., 2018) to target the EEG-negative (high-excitability state) versus EEG-positive peak (low-excitability state) of the sensorimotor μ-rhythm in the left (conditioning) and right (test) motor cortex (M1) of 16 healthy human subjects (9 female, 7 male). Effective connectivity was tested by short-interval interhemispheric inhibition (SIHI); that is, the inhibitory effect of the conditioning TMS pulse given 10-12 ms before the test pulse on the test motor-evoked potential. We compared the four possible combinations of excitability states (negative peak, positive peak) and phase relations (in-phase, out-of-phase) of the μ-rhythm in the conditioning and test M1 and a random phase condition. Strongest SIHI was found when the two M1 were in phase for the high-excitability state (negative peak of the μ-rhythm), whereas the weakest SIHI occurred when they were out of phase and the conditioning M1 was in the low-excitability state (positive peak). Phase synchronicity contributed significantly to SIHI variation, with stronger SIHI in the in-phase than out-of-phase conditions. These findings are in exact accord with the predictions of the theory of communication through coherence. They open a translational route for highly effective modification of brain connections by repetitive stimulation at instants in time when nodes in the network are phase synchronized and excitable.SIGNIFICANCE STATEMENT The theory of communication through coherence predicts that effective connectivity between nodes in distributed oscillating brain networks depends on their instantaneous excitability and phase relation. We tested this hypothesis in healthy human subjects by real-time analysis of brain states by electroencephalography in combination with transcranial magnetic stimulation of left and right motor cortex. We found that short-interval interhemispheric inhibition, a marker of interhemispheric effective connectivity, was maximally expressed when the two motor cortices were in phase for a high-excitability state (the trough of the sensorimotor μ-rhythm). We conclude that findings are consistent with the theory of communication through coherence. They open a translational route to highly effectively modify brain connections by repetitive stimulation at instants in time of phase-synchronized high-excitability states.

TMSmap - Software for Quantitative Analysis of TMS Mapping Results

The use of the MRI-navigation system ensures accurate targeting of TMS. This, in turn, results in TMS motor mapping becoming a routinely used procedure in neuroscience and neurosurgery. However, currently, there is no standardized methodology for assessment of TMS motor-mapping results. Therefore, we developed TMSmap—free standalone graphical interface software for the quantitative analysis of the TMS motor mapping results (http://tmsmap.ru/). In addition to the estimation of standard parameters (such as the size of cortical muscle representation and the center of gravity location), it allows estimation of the volume of cortical representations, excitability profile of the cortical surface map, and the overlap between cortical representations. The input data for the software includes the coordinates of the coil position (or electric field maximum) and the corresponding response in each stimulation point. TMSmap has been developed for versatile assessment and comparison of TMS maps relating to different experimental interventions including, but not limited to longitudinal, pharmacological and clinical studies (e.g., stroke recovery). To illustrate the use of TMSmap we provide examples of the actual TMS motor-mapping analysis of two healthy subjects and one chronic stroke patient.

Impact of non-brain anatomy and coil orientation on inter- and intra-subject variability in TMS at midline

OBJECTIVE

To investigate inter-subject variability with respect to cerebrospinal fluid thickness and brain-scalp distance, and to investigate intra-subject variability with different coil orientations.

METHODS

Simulations of the induced electric field (E-Field) using a figure-8 coil over the vertex were conducted on 50 unique head models and varying orientations on 25 models. Metrics exploring stimulation intensity, spread, and localization were used to describe inter-subject variability and effects of non-brain anatomy.

RESULTS

Both brain-scalp distance and CSF thickness were correlated with weaker stimulation intensity and greater spread. Coil rotations show that for the dorsal portion of the stimulated brain, E-Field intensities are highest when the anterior-posterior axis of the coil is perpendicular to the longitudinal fissure, but highest for the medial portion of the stimulated brain when the coil is oriented parallel to the longitudinal fissure.

CONCLUSIONS

Normal anatomical variation in healthy individuals leads to significant differences in the site of TMS, the intensity, and the spread. These variables are generally neglected but could explain significant variability in basic and clinical studies.

SIGNIFICANCE

This is the first work to show how brain-scalp distance and cerebrospinal fluid thickness influence focality, and to show the disassociation between dorsal and medial TMS.

TMS Motor Thresholds Correlate With TDCS Electric Field Strengths in Hand Motor Area

Transcranial direct current stimulation (TDCS) modulates cortical activity and influences motor and cognitive functions in both healthy and clinical populations. However, there is large inter-individual variability in the responses to TDCS. Computational studies have suggested that inter-individual differences in cranial and brain anatomy may contribute to this variability via creating varying electric fields in the brain. This implies that the electric fields or their strength and orientation should be considered and incorporated when selecting the TDCS dose. Unfortunately, electric field modeling is difficult to perform; thus, a more-robust and practical method of estimating the strength of TDCS electric fields for experimental use is required. As recent studies have revealed a relationship between the sensitivity to TMS and motor cortical TDCS after-effects, the aim of the present study was to investigate whether the resting motor threshold (RMT), a simple measure of transcranial magnetic stimulation (TMS) sensitivity, would be useful for estimating TDCS electric field strengths in the hand area of primary motor cortex (M1). To achieve this, we measured the RMT in 28 subjects. We also obtained magnetic resonance images from each subject to build individual three-dimensional anatomic models, which were used in solving the TDCS and TMS electric fields using the finite element method (FEM). Then, we calculated the correlation between the measured RMT and the modeled TDCS electric fields. We found that the RMT correlated with the TDCS electric fields in hand M1 (R² = 0.58), but no obvious correlations were identified in regions outside M1. The found correlation was mainly due to a correlation between the TDCS and TMS electric fields, both of which were affected by individual's anatomic features. In conclusion, the RMT could provide a useful tool for estimating cortical electric fields for motor cortical TDCS.