DELMEP: a deep learning algorithm for automated annotation of motor evoked potential latencies.

The analysis of motor evoked potentials (MEPs) generated by transcranial magnetic stimulation (TMS) is crucial in research and clinical medical practice. MEPs are characterized by their latency and the treatment of a single patient may require the characterization of thousands of MEPs. Given the difficulty of developing reliable and accurate algorithms, currently the assessment of MEPs is performed with visual inspection and manual annotation by a medical expert; making it a time-consuming, inaccurate, and error-prone process. In this study, we developed DELMEP, a deep learning-based algorithm to automate the estimation of MEP latency. Our algorithm resulted in a mean absolute error of about 0.5 ms and an accuracy that was practically independent of the MEP amplitude. The low computational cost of the DELMEP algorithm allows employing it in on-the-fly characterization of MEPs for brain-state-dependent and closed-loop brain stimulation protocols. Moreover, its learning ability makes it a particularly promising option for artificial-intelligence-based personalized clinical applications.

Advanced Pipeline for Designing Multi-Locus TMS Coils With Current Density Constraints.

OBJECTIVE: This work aims for a method to design manufacturable windings for transcranial magnetic stimulation (TMS) coils with fine control over the induced electric field (E-field) distributions. Such TMS coils are required for multi-locus TMS (mTMS). METHODS: We introduce a new mTMS coil design workflow with increased flexibility in target E-field definition and faster computations compared to our previous method. We also incorporate custom current density and E-field fidelity constraints to ensure that the target E-fields are accurately reproduced with feasible winding densities in the resulting coil designs. We validated the method by designing, manufacturing, and characterizing a 2-coil mTMS transducer for focal rat brain stimulation. RESULTS: Applying the constraints reduced the computed maximum surface current densities from 15.4 and 6.6 kA/mm to the target value 4.7 kA/mm, yielding winding paths suitable for a 1.5-mm-diameter wire with 7-kA maximum currents while still replicating the target E-fields with the predefined 2.8% maximum error in the FOV. The optimization time was reduced by two thirds compared to our previous method. CONCLUSION: The developed method allowed us to design a manufacturable, focal 2-coil mTMS transducer for rat TMS impossible to attain with our previous design workflow. SIGNIFICANCE: The presented workflow enables considerably faster design and manufacturing of previously unattainable mTMS transducers with increased control over the induced E-field distribution and winding density, opening new possibilities for brain research and clinical TMS.

TMS-Induced Modulation of EEG Functional Connectivity Is Affected by the E-Field Orientation.

Coregistration of transcranial magnetic stimulation (TMS) and electroencephalography (EEG) allows non-invasive probing of brain circuits: TMS induces brain activation due to the generation of a properly oriented focused electric field (E-field) using a coil placed on a selected position over the scalp, while EEG captures the effects of the stimulation on brain electrical activity. Moreover, the combination of these techniques allows the investigation of several brain properties, including brain functional connectivity. The choice of E-field parameters, such as intensity, orientation, and position, is crucial for eliciting cortex-specific effects. Here, we evaluated whether and how the spatial pattern, i.e., topography and strength of functional connectivity, is modulated by the stimulus orientation. We systematically altered the E-field orientation when stimulating the left pre-supplementary motor area and showed an increase of functional connectivity in areas associated with the primary motor cortex and an E-field orientation-specific modulation of functional connectivity intensity.

Accuracy and precision of navigated transcranial magnetic stimulation.

Objective.Transcranial magnetic stimulation (TMS) induces an electric field (E-field) in the cortex. To facilitate stimulation targeting, image-guided neuronavigation systems have been introduced. Such systems track the placement of the coil with respect to the head and visualize the estimated cortical stimulation location on an anatomical brain image in real time. The accuracy and precision of the neuronavigation is affected by multiple factors. Our aim was to analyze how different factors in TMS neuronavigation affect the accuracy and precision of the coil-head coregistration and the estimated E-field.Approach.By performing simulations, we estimated navigation errors due to distortions in magnetic resonance images (MRIs), head-to-MRI registration (landmark- and surface-based registrations), localization and movement of the head tracker, and localization of the coil tracker. We analyzed the effect of these errors on coil and head coregistration and on the induced E-field as determined with simplistic and realistic head models.Main results.Average total coregistration accuracies were in the range of 2.2-3.6 mm and 1°; precision values were about half of the accuracy values. The coregistration errors were mainly due to head-to-MRI registration with average accuracies 1.5-1.9 mm/0.2-0.4° and precisions 0.5-0.8 mm/0.1-0.2° better with surface-based registration. The other major source of error was the movement of the head tracker with average accuracy of 1.5 mm and precision of 1.1 mm. When assessed within an E-field method, the average accuracies of the peak E-field location, orientation, and magnitude ranged between 1.5 and 5.0 mm, 0.9 and 4.8°, and 4.4 and 8.5% across the E-field models studied. The largest errors were obtained with the landmark-based registration. When computing another accuracy measure with the most realistic E-field model as a reference, the accuracies tended to improve from about 10 mm/15°/25% to about 2 mm/2°/5% when increasing realism of the E-field model.Significance.The results of this comprehensive analysis help TMS operators to recognize the main sources of error in TMS navigation and that the coregistration errors and their effect in the E-field estimation depend on the methods applied. To ensure reliable TMS navigation, we recommend surface-based head-to-MRI registration and realistic models for E-field computations.

Transcranial magnetic stimulation set-up for small animals.

Transcranial magnetic stimulation (TMS) is widely applied on humans for research and clinical purposes. TMS studies on small animals, e.g., rodents, can provide valuable knowledge of the underlying neurophysiological mechanisms. Administering TMS on small animals is, however, prone to technical difficulties, mainly due to their small head size. In this study, we aimed to develop an energy-efficient coil and a compatible experimental set-up for administering TMS on rodents. We applied a convex optimization process to develop a minimum-energy coil for TMS on rats. As the coil windings of the optimized coil extend to a wide region, we designed and manufactured a holder on which the rat lies upside down, with its head supported by the coil. We used the set-up to record TMS-electromyography, with electromyography recorded from limb muscles with intramuscular electrodes. The upside-down placement of the rat allowed the operator to easily navigate the TMS without the coil blocking their field of view. With this paradigm, we obtained consistent motor evoked potentials from all tested animals.

Closed-loop optimization of transcranial magnetic stimulation with electroencephalography feedback.

BACKGROUND: Transcranial magnetic stimulation (TMS) is widely used in brain research and treatment of various brain dysfunctions. However, the optimal way to target stimulation and administer TMS therapies, for example, where and in which electric field direction the stimuli should be given, is yet to be determined. OBJECTIVE: To develop an automated closed-loop system for adjusting TMS parameters (in this work, the stimulus orientation) online based on TMS-evoked brain activity measured with electroencephalography (EEG). METHODS: We developed an automated closed-loop TMS-EEG set-up. In this set-up, the stimulus parameters are electronically adjusted with multi-locus TMS. As a proof of concept, we developed an algorithm that automatically optimizes the stimulation orientation based on single-trial EEG responses. We applied the algorithm to determine the electric field orientation that maximizes the amplitude of the TMS-EEG responses. The validation of the algorithm was performed with six healthy volunteers, repeating the search twenty times for each subject. RESULTS: The validation demonstrated that the closed-loop control worked as desired despite the large variation in the single-trial EEG responses. We were often able to get close to the orientation that maximizes the EEG amplitude with only a few tens of pulses. CONCLUSION: Optimizing stimulation with EEG feedback in a closed-loop manner is feasible and enables effective coupling to brain activity.

Quantifying arousal and awareness in altered states of consciousness using interpretable deep learning.

Consciousness can be defined by two components: arousal (wakefulness) and awareness (subjective experience). However, neurophysiological consciousness metrics able to disentangle between these components have not been reported. Here, we propose an explainable consciousness indicator (ECI) using deep learning to disentangle the components of consciousness. We employ electroencephalographic (EEG) responses to transcranial magnetic stimulation under various conditions, including sleep (n = 6), general anesthesia (n = 16), and severe brain injury (n = 34). We also test our framework using resting-state EEG under general anesthesia (n = 15) and severe brain injury (n = 34). ECI simultaneously quantifies arousal and awareness under physiological, pharmacological, and pathological conditions. Particularly, ketamine-induced anesthesia and rapid eye movement sleep with low arousal and high awareness are clearly distinguished from other states. In addition, parietal regions appear most relevant for quantifying arousal and awareness. This indicator provides insights into the neural correlates of altered states of consciousness.

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.

Local brain-state dependency of effective connectivity: a pilot TMS-EEG study.

Background: Spontaneous cortical oscillations have been shown to modulate cortical responses to transcranial magnetic stimulation (TMS). However, whether these oscillations influence cortical effective connectivity is largely unknown. We conducted a pilot study to set the basis for addressing how spontaneous oscillations affect cortical effective connectivity measured through TMS-evoked potentials (TEPs). Methods: We applied TMS to the left primary motor cortex and right pre-supplementary motor area of three subjects while recording EEG. We classified trials off-line into positive- and negative-phase classes according to the mu and beta rhythms. We calculated differences in the global mean-field amplitude (GMFA) and compared the cortical spreading of the TMS-evoked activity between the two classes. Results: Phase affected the GMFA in four out of 12 datasets (3 subjects × 2 stimulation sites × 2 frequency bands). Two of the observed significant intervals were before 50 ms, two between 50 and 100 ms, and one after 100 ms post-stimulus. Source estimates showed complex spatial differences between the classes in the cortical spreading of the TMS-evoked activity. Conclusions: TMS-evoked effective connectivity seems to depend on the phase of local cortical oscillations at the stimulated site. This work paves the way to design future closed-loop stimulation paradigms.

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.

Trade-off between stimulation focality and the number of coils in multi-locus transcranial magnetic stimulation.

Objective. Coils designed for transcranial magnetic stimulation (TMS) must incorporate trade-offs between the required electrical power or energy, focality and depth penetration of the induced electric field (E-field), coil size, and mechanical properties of the coil, as all of them cannot be optimally met at the same time. In multi-locus TMS (mTMS), a transducer consisting of several coils allows electronically targeted stimulation of the cortex without physically moving a coil. In this study, we aimed to investigate the relationship between the number of coils in an mTMS transducer, the focality of the induced E-field, and the extent of the cortical region within which the location and orientation of the maximum of the induced E-field can be controlled.Approach.We applied convex optimization to design planar and spherically curved mTMS transducers of different E-field focalities and analyzed their properties. We characterized the trade-off between the focality of the induced E-field and the extent of the cortical region that can be stimulated with an mTMS transducer with a given number of coils.Main results.At the expense of the E-field focality, one can, with the same number of coils, design an mTMS transducer that can control the location and orientation of the peak of the induced E-field within a wider cortical region.Significance. With E-fields of moderate focality, the problem of electronically targeted TMS becomes considerably easier compared with highly focal E-fields; this may speed up the development of mTMS and the emergence of new clinical and research applications.

Effect of stimulus orientation and intensity on short-interval intracortical inhibition (SICI) and facilitation (SICF): A multi-channel transcranial magnetic stimulation study.

Besides stimulus intensities and interstimulus intervals (ISI), the electric field (E-field) orientation is known to affect both short-interval intracortical inhibition (SICI) and facilitation (SICF) in paired-pulse transcranial magnetic stimulation (TMS). However, it has yet to be established how distinct orientations of the conditioning (CS) and test stimuli (TS) affect the SICI and SICF generation. With the use of a multi-channel TMS transducer that provides electronic control of the stimulus orientation and intensity, we aimed to investigate how changes in the CS and TS orientation affect the strength of SICI and SICF. We hypothesized that the CS orientation would play a major role for SICF than for SICI, whereas the CS intensity would be more critical for SICI than for SICF. In eight healthy subjects, we tested two ISIs (1.5 and 2.7 ms), two CS and TS orientations (anteromedial (AM) and posteromedial (PM)), and four CS intensities (50, 70, 90, and 110% of the resting motor threshold (RMT)). The TS intensity was fixed at 110% RMT. The intensities were adjusted to the corresponding RMT in the AM and PM orientations. SICI and SICF were observed in all tested CS and TS orientations. SICI depended on the CS intensity in a U-shaped manner in any combination of the CS and TS orientations. With 70% and 90% RMT CS intensities, stronger PM-oriented CS induced stronger inhibition than weaker AM-oriented CS. Similar SICF was observed for any CS orientation. Neither SICI nor SICF depended on the TS orientation. We demonstrated that SICI and SICF could be elicited by the CS perpendicular to the TS, which indicates that these stimuli affected either overlapping or strongly connected neuronal populations. We concluded that SICI is primarily sensitive to the CS intensity and that CS intensity adjustment resulted in similar SICF for different CS orientations.

Individual head models for estimating the TMS-induced electric field in rat brain.

In transcranial magnetic stimulation (TMS), the initial cortical activation due to stimulation is determined by the state of the brain and the magnitude, waveform, and direction of the induced electric field (E-field) in the cortex. The E-field distribution depends on the conductivity geometry of the head. The effects of deviations from a spherically symmetric conductivity profile have been studied in detail in humans. In small mammals, such as rats, these effects are more pronounced due to their less spherical head, proportionally much thicker neck region, and overall much smaller size compared to the TMS coils. In this study, we describe a simple method for building individual realistically shaped head models for rats from high-resolution X-ray tomography images. We computed the TMS-induced E-field with the boundary element method and assessed the effect of head-model simplifications on the estimated E-field. The deviations from spherical symmetry have large, non-trivial effects on the E-field distribution: for some coil orientations, the strongest stimulation is in the brainstem even when the coil is over the motor cortex. With modelling prior to an experiment, such problematic coil orientations can be avoided for more accurate targeting.

Automated search of stimulation targets with closed-loop transcranial magnetic stimulation.

Transcranial magnetic stimulation (TMS) protocols often include a manual search of an optimal location and orientation of the coil or peak stimulating electric field to elicit motor responses in a target muscle. This target search is laborious, and the result is user-dependent. Here, we present a closed-loop search method that utilizes automatic electronic adjustment of the stimulation based on the previous responses. The electronic adjustment is achieved by multi-locus TMS, and the adaptive guiding of the stimulation is based on the principles of Bayesian optimization to minimize the number of stimuli (and time) needed in the search. We compared our target-search method with other methods, such as systematic sampling in a predefined cortical grid. Validation experiments on five healthy volunteers and further offline simulations showed that our adaptively guided search method needs only a relatively small number of stimuli to provide outcomes with good accuracy and precision. The automated method enables fast and user-independent optimization of stimulation parameters in research and clinical applications of TMS.

Interhemispheric symmetry of µ-rhythm phase-dependency of corticospinal excitability.

Oscillatory activity in the µ-frequency band (8-13 Hz) determines excitability in sensorimotor cortex. In humans, the primary motor cortex (M1) in the two hemispheres shows significant anatomical, connectional, and electrophysiological differences associated with motor dominance. It is currently unclear whether the µ-oscillation phase effects on corticospinal excitability demonstrated previously for the motor-dominant M1 are also different between motor-dominant and motor-non-dominant M1 or, alternatively, are similar to reflect a ubiquitous physiological trait of the motor system at rest. Here, we applied single-pulse transcranial magnetic stimulation to the hand representations of the motor-dominant and the motor-non-dominant M1 of 51 healthy right-handed volunteers when electroencephalography indicated a certain µ-oscillation phase (positive peak, negative peak, or random). We determined resting motor threshold (RMT) as a marker of corticospinal excitability in the three µ-phase conditions. RMT differed significantly depending on the pre-stimulus phase of the µ-oscillation in both M1, with highest RMT in the positive-peak condition, and lowest RMT in the negative-peak condition. µ-phase-dependency of RMT correlated directly between the two M1, and interhemispheric differences in µ-phase-dependency were absent. In conclusion, µ-phase-dependency of corticospinal excitability appears to be a ubiquitous physiological trait of the motor system at rest, without hemispheric dominance.

The shaky ground truth of real-time phase estimation.

Instantaneous phase of brain oscillations in electroencephalography (EEG) is a measure of brain state that is relevant to neuronal processing and modulates evoked responses. However, determining phase at the time of a stimulus with standard signal processing methods is not possible due to the stimulus artifact masking the future part of the signal. Here, we quantify the degree to which signal-to-noise ratio and instantaneous amplitude of the signal affect the variance of phase estimation error and the precision with which "ground truth" phase is even defined, using both the variance of equivalent estimators and realistic simulated EEG data with known synthetic phase. Necessary experimental conditions are specified in which pre-stimulus phase estimation is meaningfully possible based on instantaneous amplitude and signal-to-noise ratio of the oscillation of interest. An open source toolbox is made available for causal (using pre-stimulus signal only) phase estimation along with a EEG dataset consisting of recordings from 140 participants and a best practices workflow for algorithm optimization and benchmarking. As an illustration, post-hoc sorting of open-loop transcranial magnetic stimulation (TMS) trials according to pre-stimulus sensorimotor µ-rhythm phase is performed to demonstrate modulation of corticospinal excitability, as indexed by the amplitude of motor evoked potentials.

Short-interval intracortical inhibition in human primary motor cortex: A multi-locus transcranial magnetic stimulation study.

Short-interval intracortical inhibition (SICI) has been studied with paired-pulse transcranial magnetic stimulation (TMS) by administering two pulses at a millisecond-scale interstimulus interval (ISI) to a single cortical target. It has, however, been difficult to study the interaction of nearby cortical targets with paired-pulse TMS. To overcome this limitation, we have developed a multi-locus TMS (mTMS) device, which allows controlling the stimulus location electronically. Here, we applied mTMS to study SICI in primary motor cortex with paired pulses targeted to adjacent locations, aiming to quantify the extent of the cortical region producing SICI in the location of a test stimulus. We varied the location and timing of the conditioning stimulus with respect to a test stimulus targeted to the cortical hotspot of the abductor pollicis brevis (APB) in order to study their effects on motor evoked potentials. We further applied a two-coil protocol with the conditioning stimulus given by an oval coil only to the surroundings of the APB hotspot, to which a subsequent test stimulus was administered with a figure-of-eight coil. The strongest SICI occurred at ISIs below 1 ms and at ISIs around 2.5 ms. These ISIs increased when the conditioning stimulus receded from the APB hotspot. Our two-coil paired-pulse TMS study suggests that SICI at ISIs of 0.5 and 2.5 ms originate from different mechanisms or neuronal elements.

Brain State-dependent Brain Stimulation with Real-time Electroencephalography-Triggered Transcranial Magnetic Stimulation.

The effect of a stimulus to the brain depends not only on the parameters of the stimulus but also on the dynamics of brain activity at the time of the stimulation. The combination of electroencephalography (EEG) and transcranial magnetic stimulation (TMS) in a real-time brain state-dependent stimulation system allows the study of relations of dynamics of brain activity, cortical excitability, and plasticity induction. Here, we demonstrate a newly developed method to synchronize the timing of brain stimulation with the phase of ongoing EEG oscillations using a real-time data analysis system. This real-time EEG-triggered TMS of the human motor cortex, when TMS is synchronized with the surface EEG negative peak of the sensorimotor µ-alpha (8-14 Hz) rhythm, has shown differential corticospinal excitability and plasticity effects. The utilization of this method suggests that real-time information about the instantaneous brain state can be used for efficacious plasticity induction. Additionally, this approach enables personalized EEG-synchronized brain stimulation which may lead to the development of more effective therapeutic brain stimulation protocols.

Connectivity differences between consciousness and unconsciousness in non-rapid eye movement sleep: a TMS-EEG study.

The neuronal connectivity patterns that differentiate consciousness from unconsciousness remain unclear. Previous studies have demonstrated that effective connectivity, as assessed by transcranial magnetic stimulation combined with electroencephalography (TMS-EEG), breaks down during the loss of consciousness. This study investigated changes in EEG connectivity associated with consciousness during non-rapid eye movement (NREM) sleep following parietal TMS. Compared with unconsciousness, conscious experiences during NREM sleep were associated with reduced phase-locking at low frequencies (<4 Hz). Transitivity and clustering coefficient in the delta and theta bands were also significantly lower during consciousness compared to unconsciousness, with differences in the clustering coefficient observed in scalp electrodes over parietal-occipital regions. There were no significant differences in Granger-causality patterns in frontal-to-parietal or parietal-to-frontal connectivity between reported unconsciousness and reported consciousness. Together these results suggest that alterations in spectral and spatial characteristics of network properties in posterior brain areas, in particular decreased local (segregated) connectivity at low frequencies, is a potential indicator of consciousness during sleep.

The effect of experimental pain on short-interval intracortical inhibition with multi-locus transcranial magnetic stimulation.

Chronic neuropathic pain is known to alter the primary motor cortex (M1) function. Less is known about the normal, physiological effects of experimental neurogenic pain on M1. The objective of this study is to determine how short-interval intracortical inhibition (SICI) is altered in the M1 representation area of a muscle exposed to experimental pain compared to SICI of another muscle not exposed to pain. The cortical representation areas of the right abductor pollicis brevis (APB) and biceps brachii (BB) muscles of 11 subjects were stimulated with a multi-locus transcranial magnetic stimulation device while the resulting motor-evoked potentials (MEPs) were recorded with electromyography. Single- and paired-pulse TMS was administered in seven conditions, including one with the right hand placed in cold water. The stimulation intensity for the conditioning pulses in the paired-pulse examination was 80% of the resting motor threshold (RMT) of the stimulated site and 120% of RMT for both the test and single pulses. The paired-pulse MEP amplitudes were normalized with the mean amplitude of the single-pulse MEPs of the same condition and muscle. SICI was compared between conditions. After the cold pain, the normalized paired-pulse MEP amplitudes decreased in APB, but not in BB, indicating that SICI was potentially increased only in the cortical area of the muscle subjected to pain. These data suggest that SICI is increased in the M1 representation area of a hand muscle shortly after exposure to pain has ended, which implies that short-lasting pain can alter the inhibitory balance in M1.