A simulation-based analysis of motor unit number index (MUNIX) technique using motoneuron pool and surface electromyogram models

Model-based closed-loop control of thalamic deep brain stimulation

Introduction: Closed-loop control of deep brain stimulation (DBS) is beneficial for effective and automatic treatment of various neurological disorders like Parkinson's disease (PD) and essential tremor (ET). Manual (open-loop) DBS programming solely based on clinical observations relies on neurologists' expertise and patients' experience. Continuous stimulation in open-loop DBS may decrease battery life and cause side effects. On the contrary, a closed-loop DBS system uses a feedback biomarker/signal to track worsening (or improving) of patients' symptoms and offers several advantages compared to the open-loop DBS system. Existing closed-loop DBS control systems do not incorporate physiological mechanisms underlying DBS or symptoms, e.g., how DBS modulates dynamics of synaptic plasticity. Methods: In this work, we propose a computational framework for development of a model-based DBS controller where a neural model can describe the relationship between DBS and neural activity and a polynomial-based approximation can estimate the relationship between neural and behavioral activities. A controller is used in our model in a quasi-real-time manner to find DBS patterns that significantly reduce the worsening of symptoms. By using the proposed computational framework, these DBS patterns can be tested clinically by predicting the effect of DBS before delivering it to the patient. We applied this framework to the problem of finding optimal DBS frequencies for essential tremor given electromyography (EMG) recordings solely. Building on our recent network model of ventral intermediate nuclei (Vim), the main surgical target of the tremor, in response to DBS, we developed neural model simulation in which physiological mechanisms underlying Vim-DBS are linked to symptomatic changes in EMG signals. By using a proportional-integral-derivative (PID) controller, we showed that a closed-loop system can track EMG signals and adjust the stimulation frequency of Vim-DBS so that the power of EMG reaches a desired control target. Results and discussion: We demonstrated that the model-based DBS frequency aligns well with that used in clinical studies. Our model-based closed-loop system is adaptable to different control targets and can potentially be used for different diseases and personalized systems.

The filling factor of the sEMG signal at low contraction forces in the quadriceps muscles is influenced by the thickness of the subcutaneous layer

Introduction: It has been shown that, for male subjects, the sEMG activity at low contraction forces is normally "pulsatile", i.e., formed by a few large-amplitude MUPs, coming from the most superficial motor units. The subcutaneous layer thickness, known to be greater in females than males, influences the electrode detection volume. Here, we investigated the influence of the subcutaneous layer thickness on the type of sEMG activity (pulsatile vs. continuous) at low contraction forces. Methods: Voluntary surface EMG signals were recorded from the quadriceps muscles of healthy males and females as force was gradually increased from 0% to 40% MVC. The sEMG filling process was examined by measuring the EMG filling factor, computed from the non-central moments of the rectified sEMG signal. Results: 1) The sEMG activity at low contraction forces was "continuous" in the VL, VM and RF of females, whereas this sEMG activity was "pulsatile" in the VL and VM of males. 2) The filling factor at low contraction forces was lower in males than females for the VL (p = 0.003) and VM (p = 0.002), but not for the RF (p = 0.54). 3) The subcutaneous layer was significantly thicker in females than males for the VL (p = 0.001), VM (p = 0.001), and RF (p = 0.003). 4) A significant correlation was found in the vastus muscles between the subcutaneous layer thickness and the filling factor (p < 0.05). Discussion: The present results indicate that the sEMG activity at low contraction forces in the female quadriceps muscles is "continuous" due to the thick subcutaneous layer of these muscles, which impedes an accurate assessment of the sEMG filling process.

Compound muscle action potential (CMAP) scan examination of paretic and contralateral muscles reveals motor unit alterations after stroke

This study presents a novel compound muscle action potential (CMAP) examination of motor unit changes in paretic muscle post stroke. CMAP scan of the first dorsal interosseous (FDI) muscle was performed bilaterally in 16 chronic stroke subjects. Various parameters were derived from the CMAP scan to examine paretic muscle changes, including CMAP amplitude, D50, step index (STEPIX) and amplitude index (AMPIX). A significant decrease in CMAP amplitude and STEPIX was observed in paretic muscles compared with contralateral muscles (CMAP amplitude: paretic (9.0±0.5) mV, contralateral (11.3±0.9) mV, P=0.024; STEPIX: paretic 101.2±7.6, contralateral 121.9±6.5, P=0.020). No significant difference in D50 and AMPIX was observed between the paretic and contralateral sides (P>0.05). The findings revealed complex paretic muscle changes including motor unit degeneration, muscle fiber denervation, reinnervation and atrophy, providing useful insights to help understand neuromuscular mechanisms associated with weakness and other functional deterioration post stroke. The CMAP scan experimental protocols and the applied processing methods are noninvasive, convenient, and automated, offering practical benefits for clinical application.

Motor unit properties do not correlate between MUNIX and needle EMG in remote polio in the biceps brachii muscle

Objective

To compare the utility of MUNIX (motor unit number index) with needle EMG in characterizing motor unit (MU) properties in the biceps brachii (BB) muscle in subjects with remote polio.

Methods

Thirty subjects suffering from remote polio were investigated with MUNIX and needle EMG, all with Macro EMG and 16 of these subjects with concentric needle EMG.

Results

Both MUNIX and the needle EMG methods showed abnormal results. Fiber density (FD) was the most sensitive parameter for showing signs of reinnervation. At a group level, the methods showed neurogenic findings, but there was no correlation between the results of the MUNIX and needle EMG investigations.

Conclusions

Both MUNIX and needle EMG are valuable methods for measuring neurogenic involvement in the BB muscle. However, there was a lack of correlation between the MUNIX and needle EMG findings. The cause for this missing correlation may be multifactorial as there are several differences between the methods.

Significance

The reason for the lack of correlation between the MUNIX and needle EMG results is discussed. By combining the needle and surface recorded methods one can obtain more information on the denervation and reinnervation process compared to using just one of the methods alone.

Application Value of the Motor Unit Number Index in Patients With Kennedy Disease

The aim of this study was to evaluate the usefulness of the motor unit number index (MUNIX) technique in Kennedy disease (KD) and test the correlation between the MUNIX and other clinical parameters. The MUNIX values of the bilateral deltoid, abductor digiti minimi (ADM), quadriceps femoris (QF), and tibialis anterior (TA) were determined and compared with the course of the disease. The MUNIX sum score was calculated by adding the MUNIX values of these 8 muscles. Disability was evaluated using the spinal and bulbar muscular atrophy functional rating scale (SBMAFRS). The MUNIX scores of patients with KD were negatively correlated with the course of the disease (p < 0.05), whereas their motor unit size index (MUSIX) scores were positively correlated with the course the of disease (p < 0.05). MUNIX sum scores were correlated with SBMAFRS scores (r = 0.714, p < 0.05). MUNIX was more sensitive than compound muscle action potentials or muscle strength as an indicator of neuron loss and axonal collateral reinnervation. The MUNIX sum score is an objective and a reliable indicator of disease progression, and it is a potential choice for therapeutic clinical trials. The MUNIX can assess the functional loss of motor axons and is correlated with disability. The MUNIX sum score may be especially suitable as an objective parameter.

Model-Based Analysis of Muscle Strength and EMG-Force Relation with respect to Different Patterns of Motor Unit Loss

This study presents a model-based sensitivity analysis of the strength of voluntary muscle contraction with respect to different patterns of motor unit loss. A motor unit pool model was implemented including simulation of a motor neuron pool, muscle force, and surface electromyogram (EMG) signals. Three different patterns of motor unit loss were simulated, including (1) motor unit loss restricted to the largest ones, (2) motor unit loss restricted to the smallest ones, and (3) motor unit loss without size restriction. The model outputs including muscle force amplitude, variability, and the resultant EMG-force relation were quantified under two different motor neuron firing strategies. It was found that motor unit loss restricted to the largest ones had the most dominant impact on muscle strength and significantly changed the EMG-force relation, while loss restricted to the smallest motor units had a pronounced effect on force variability. These findings provide valuable insight toward our understanding of the neurophysiological mechanisms underlying experimental observations of muscle strength, force control, and EMG-force relation in both normal and pathological conditions.

Modeling motor-evoked potentials from neural field simulations of transcranial magnetic stimulation

OBJECTIVE

To develop a population-based biophysical model of motor-evoked potentials (MEPs) following transcranial magnetic stimulation (TMS).

METHODS

We combined an existing MEP model with population-based cortical modeling. Layer 2/3 excitatory and inhibitory neural populations, modeled with neural-field theory, are stimulated with TMS and feed layer 5 corticospinal neurons, which also couple directly but weakly to the TMS pulse. The layer 5 output controls mean motoneuron responses, which generate a series of single motor-unit action potentials that are summed to estimate a MEP.

RESULTS

A MEP waveform was generated comparable to those observed experimentally. The model captured TMS phenomena including a sigmoidal input-output curve, common paired pulse effects (short interval intracortical inhibition, intracortical facilitation, long interval intracortical inhibition) including responses to pharmacological interventions, and a cortical silent period. Changes in MEP amplitude following theta burst paradigms were observed including variability in outcome direction.

CONCLUSIONS

The model reproduces effects seen in common TMS paradigms.

SIGNIFICANCE

The model allows population-based modeling of changes in cortical dynamics due to TMS protocols to be assessed in terms of changes in MEPs, thus allowing a clear comparison between population-based modeling predictions and typical experimental outcome measures.

Motor unit number index (MUNIX) in myopathic disorders: Clinical correlations and potential pitfalls

Our aim was to study motor unit number index (MUNIX) in myopathic disorders. We studied 11 patients with myopathy, and healthy controls. We obtained MUNIX, compound muscle action potential (CMAP), motor unit size index (MUSIX) and alpha (α, power exponent from MUNIX equation) measurements from three different muscles. MUNIX and CMAP were significantly lower in one muscle. This MUNIX decrease may not be related to motor neuron loss, but rather to muscle fiber atrophy. MUSIX and α did not change and may be useful in determining whether the MUNIX decrease is indeed due to motor unit loss.

Simulation of electromyographic recordings following transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) is a technique that enables noninvasive manipulation of neural activity and holds promise in both clinical and basic research settings. The effect of TMS on the motor cortex is often measured by electromyography (EMG) recordings from a small hand muscle. However, the details of how TMS generates responses measured with EMG are not completely understood. We aim to develop a biophysically detailed computational model to study the potential mechanisms underlying the generation of EMG signals following TMS. Our model comprises a feed-forward network of cortical layer 2/3 cells, which drive morphologically detailed layer 5 corticomotoneuronal cells, which in turn project to a pool of motoneurons. EMG signals are modeled as the sum of motor unit action potentials. EMG recordings from the first dorsal interosseous muscle were performed in four subjects and compared with simulated EMG signals. Our model successfully reproduces several characteristics of the experimental data. The simulated EMG signals match experimental EMG recordings in shape and size, and change with stimulus intensity and contraction level as in experimental recordings. They exhibit cortical silent periods that are close to the biological values and reveal an interesting dependence on inhibitory synaptic transmission properties. Our model predicts several characteristics of the firing patterns of neurons along the entire pathway from cortical layer 2/3 cells down to spinal motoneurons and should be considered as a viable tool for explaining and analyzing EMG signals following TMS. NEW & NOTEWORTHY A biophysically detailed model of EMG signal generation following transcranial magnetic stimulation (TMS) is proposed. Simulated EMG signals match experimental EMG recordings in shape and amplitude. Motor-evoked potential and cortical silent period properties match experimental data. The model is a viable tool to analyze, explain, and predict EMG signals following TMS.

Motor unit number index (MUNIX) derivation from the relationship between the area and power of surface electromyogram: a computer simulation and clinical study

OBJECTIVE

The motor unit number index (MUNIX) is a technique based on the surface electromyogram (sEMG) that is gaining acceptance as a method for monitoring motor neuron loss, because it is reliable and produces less discomfort than other electrodiagnostic techniques having the same intended purpose. MUNIX assumes that the relationship between the area of sEMG obtained at increasing levels of muscle activation and the values of a variable called 'ideal case motor unit count' (ICMUC), defined as the product of the ratio between area and power of the compound muscle action potential (CMAP) by that of the sEMG, is described by a decreasing power function. Nevertheless, the reason for this comportment is unknown. The objective of this work is to investigate if the definition of MUNIX could derive from more basic properties of the sEMG.

APPROACH

The CMAP and sEMG epochs obtained at different levels of muscle activation from (1) the abductor pollicis brevis (APB) muscle of persons with and without a carpal tunnel syndrome (CTS) and (2) from a computer model of sEMG generation previously published were analysed.

MAIN RESULTS

MUNIX reflects the power relationship existing between the area and power of a sEMG. The exponent of this function was smaller in patients with motor CTS than in the rest of the subjects. The analysis of the relationship between the area and power of a sEMG could aid in distinguishing a MUNIX reduction due to a motoneuron loss from that due to a loss of muscle fibre.

SIGNIFICANCE

MUNIX is derived from the relationship between the area and power of a sEMG. This relationship changes when there is a loss of motor units (MUs), which partially explains the diagnostic sensibility of MUNIX. Although the reasons for this change are unknown, it could reflect an increase in the proportion of MUs of great amplitude.

Motor unit number index in the nasalis muscle in healthy subjects and patients with amyotrophic lateral sclerosis

INTRODUCTION

Motor unit number index (MUNIX) is a quick and feasible electrophysiological technique that estimates the number of motor neurons in limb muscles in healthy and amyotrophic lateral sclerosis (ALS) subjects. In this study we explored the feasibility, reliability, and differences of MUNIX in nasalis muscles in healthy subjects and ALS patients.

METHODS

MUNIX of the nasalis muscle of 50 healthy and 20 ALS subjects with bulbar involvement was compared. Functional impairment was evaluated by the ALS Functional Rating Scale-Revised and its bulbar subscore.

RESULTS

MUNIX was well tolerated and quickly performed. Bulbar ALS patients showed non-significant lower nasalis MUNIX values and a lower functional bulbar subscore. Intra- and interrater reliability showed high intraclass correlation coefficients (ICCs) in healthy subjects (0.87) and ALS patients (0.92).

CONCLUSION

MUNIX of the nasalis muscle is a reproducible method, but it showed no significant difference between healthy and bulbar ALS subjects and seems not to be a useful marker of disease progression in ALS. Muscle Nerve 54: 733-737, 2016.

Motor unit number estimation based on high-density surface electromyography decomposition

OBJECTIVE

To advance the motor unit number estimation (MUNE) technique using high density surface electromyography (EMG) decomposition.

METHODS

The K-means clustering convolution kernel compensation algorithm was employed to detect the single motor unit potentials (SMUPs) from high-density surface EMG recordings of the biceps brachii muscles in eight healthy subjects. Contraction forces were controlled at 10%, 20% and 30% of the maximal voluntary contraction (MVC). Achieved MUNE results and the representativeness of the SMUP pools were evaluated using a high-density weighted-average method.

RESULTS

Mean numbers of motor units were estimated as 288±132, 155±87, 107±99 and 132±61 by using the developed new MUNE at 10%, 20%, 30% and 10-30% MVCs, respectively. Over 20 SMUPs were obtained at each contraction level, and the mean residual variances were lower than 10%.

CONCLUSIONS

The new MUNE method allows a convenient and non-invasive collection of a large size of SMUP pool with great representativeness. It provides a useful tool for estimating the motor unit number of proximal muscles.

SIGNIFICANCE

The present new MUNE method successfully avoids the use of intramuscular electrodes or multiple electrical stimuli which is required in currently available MUNE techniques; as such the new MUNE method can minimize patient discomfort for MUNE tests.

Modified motor unit number index: A simulation study of the first dorsal interosseous muscle

The motor unit number index (MUNIX) technique has provided a quick and convenient approach to estimating motor unit population changes in a muscle. Reduction in motor unit action potential (MUAP) amplitude can lead to underestimation of motor unit numbers using the standard MUNIX technique. This study aims to overcome this limitation by developing a modified MUNIX (mMUNIX) technique. The mMUNIX uses a variable that is associated with the area of compound muscle action potential (CMAP) rather than an arbitrary fixed value (20 mV ms) as used in the standard MUNIX to define the output. The performance of the mMUNIX was evaluated using motoneuron pool and surface electromyography (EMG) models. With a fixed motor unit number, the mMUNIX output remained relatively constant with varying degrees of MUAP amplitude changes, while the standard MUNIX substantially underestimated the motor unit number in such cases. However, when MUAP amplitude remained unchanged, the mMUNIX showed less sensitivity than the standard MUNIX in tracking motor unit loss. The current simulation study demonstrated both the advantages and limitations of the standard and modified MUNIX techniques, which can help guide appropriate application and interpretation of MUNIX measurements.

Application of the F-Response for Estimating Motor Unit Number and Amplitude Distribution in Hand Muscles of Stroke Survivors

The F-response was used in this study to assess changes in the first dorsal interosseous (FDI) muscle after a hemispheric stroke. The number of motor units and their sizes were estimated bilaterally in 12 stroke survivors by recording both the compound muscle action potential (CMAP) and F wave responses. These F waves were induced by applying a large number of electrical stimuli to the ulnar nerve. The amplitude distribution of individual motor unit action potentials (MUAPs) was also compared between paretic and contralateral muscles. When averaged across all the subjects, a significantly lower motor unit number estimate was obtained for the paretic FDI muscle ( 88 ±13) compared with the contralateral side ( 139 ±11) ( ). Pooled surface MUAP amplitude analysis demonstrated a right-skewed distribution for both paretic (kurtosis 3.0) and contralateral (kurtosis 8.52) muscles. When normalized to each individual muscle's CMAP, the surface MUAP amplitude ranged from 0.22% to 4.94% (median 1.17%) of CMAP amplitude for the paretic muscle, and from 0.13% to 3.2% (median 0.62%) of CMAP amplitude for the contralateral muscle. A significant difference in MUAP outliers was also observed between the paretic and contralateral muscles. The findings of this study suggest significant motor unit loss and muscle structural reorganization after stroke.

Alterations in multidimensional motor unit number index of hand muscles after incomplete cervical spinal cord injury

The objective of this study was to apply a novel multidimensional motor unit number index (MD-MUNIX) technique to examine hand muscles in patients with incomplete cervical spinal cord injury (SCI). The MD-MUNIX was estimated from the compound muscle action potential (CMAP) and different levels of surface interference pattern electromyogram (EMG) at multiple directions of voluntary isometric muscle contraction. The MD-MUNIX was applied in the first dorsal interosseous (FDI), thenar and hypothenar muscles of SCI (n = 12) and healthy control (n = 12) subjects. The results showed that the SCI subjects had significantly smaller CMAP and MD-MUNIX in all the three examined muscles, compared to those derived from the healthy control subjects. The multidimensional motor unit size index (MD-MUSIX) demonstrated significantly larger values for the FDI and hypothenar muscles in SCI subjects than those from healthy control subjects, whereas the MD-MUSIX enlargement was marginally significant for the thenar muscles. The findings from the MD-MUNIX analyses provide an evidence of motor unit loss in hand muscles of cervical SCI patients, contributing to hand function deterioration.

An examination of motor unit number index in adults with cerebral palsy

Spinal motor neuron loss may be a factor contributing to weakness in central disorders. The aim of this study was to assess whether motor unit numbers are reduced in the hand musculature of adults with cerebral palsy (CP) using the motor unit number index (MUNIX) technique. In this prospective, case-control study, 10 adults with CP were matched with healthy controls. MUNIX was computed using area and power of voluntary surface hypothenar electromyographic (EMG) signals and the compound muscle action potential (CMAP) recorded with ulnar nerve stimulation. The motor unit size index (MUSIX) was calculated based on maximum CMAP amplitude and MUNIX value. Gross Motor Function Classification Scale (GMFCS) and Manual Abilities Classification Scale (MACS) levels were rated for CP subjects. MUNIX was significantly lower for CP participants (Mean 167.8 vs. 214.4, p=.022). MUNIX values did not correlate with GMFCS or MACS. MUSIX values were higher, though not significantly, for CP subjects (p=.11). MUSIX increased with increasing MACS levels (r(2)=.4017, p=.049). Thus, motor unit numbers in ulnar hand muscles may be decreased with CP. MUSIX values are associated with greater hand impairment. Therefore, peripheral motor unit loss as a component of the weakness found with CP deserves further evaluation.

Examination of hand muscle activation and motor unit indices derived from surface EMG in chronic stroke

In this study, we used muscle and motor unit indices, derived from convenient surface electromyography (EMG) measurements, for examination of paretic muscle changes post stroke. For 12 stroke subjects, compound muscle action potential and voluntary surface EMG signals were recorded from paretic and contralateral first dorsal interosseous, abductor pollicis brevis, and abductor digiti minimi muscles. Muscle activation index (AI), motor unit number index (MUNIX), and motor unit size index (MUSIX) were then calculated for each muscle. There was a significant AI reduction for all the three muscles in paretic side compared with contralateral side, providing an evidence of muscle activation deficiency after stroke. The hand MUNIX (defined by summing the values from the three muscles) was significantly reduced in paretic side compared with contralateral side, whereas the hand MUSIX was not significantly different. Furthermore, diverse changes in MUNIX and MUSIX were observed from the three muscles. A major feature of the present examinations is the primary reliance on surface EMG, which offers practical benefits because it is noninvasive, induces minimal discomfort and can be performed quickly.

Reliability of a modified motor unit number index (MUNIX) technique

INTRODUCTION

The purpose of this study was to examine the relative and absolute between-day reliability of the motor unit number index (MUNIX).

METHODS

Young, healthy adults (n=19) attended two testing sessions separated by 4-weeks where their maximal pinch-grip strength, MUNIX, and motor unit size index (MUSIX) were assessed in the abductor pollicis brevis muscle. Reliability was assessed by intraclass correlation coefficients (ICC), coefficient of variation (CV) and limits of agreement (LOA).

RESULTS

No mean differences were observed for MUNIX or MUSIX. The CV for the MUNIX and MUSIX measures were between 13.5% and 17.5%. The ICC for both measures were moderate to moderately-high (0.73-0.76), The LOA for both indicated a homoscedastic relationship.

DISCUSSION

Our findings indicate moderate to moderately-high reliability for both MUNIX and MUSIX. Future work is needed to ensure both measures are reliable in other muscles and cohorts, and further investigations are required to examine the validity of MUNIX.

Interrelationship between muscle strength, motor units, and aging

The interrelationship between muscle strength, motor unit (MU) number, and age is poorly understood, and in this study we sought to determine whether age-related differences in muscle strength are moderated by estimates of functioning MU number and size. Eighteen older adults (OA; 67 ± 1.20 years) and 24 young adults (YA; 22 ± 0.74 years) participated in this study. Maximum voluntary pinch-grip strength of the nondominant hand was determined and estimates of MU number were obtained from the abductor pollicis brevis muscle using the noninvasive motor unit number index (MUNIX) technique. The MUNIX technique was also utilized to derive a motor unit size index (MUSIX). An analysis of covariance (Age Group × MUNIX or MUSIX) was used to test heterogeneity of regression slopes, with body mass and gender serving as covariates. We observed that the slope of pinch-grip strength on the estimated number of MUs between YA and OA differed, indicated by an Age Group × MUNIX interaction (p = 0.04). Specifically, after controlling for the effect of body mass and gender, the slope in OA was significantly positive (0.14 ± 0.06 N/MUs, p = 0.03), whereas no such relationship was found in YA (-0.08 ± 0.09 N/MUs, p = 0.35). A significant Age Group × MUSIX interaction was also observed for strength (p < 0.01). In contrast to MUNIX, the slope in younger adults was significantly positive (0.48 ± 0.11 N/μV, p < 0.01), whereas no such relationship was found in older adults (-0.30 ± 0.22 N/μV, p = 0.18). These findings indicate that there is an interrelationship between muscle strength, MU numbers, and aging, which suggests that a portion of muscle weakness in seniors may be attributable to the loss of functioning motor units.