Abstract
Objective.
High-density surface electromyography (HD-sEMG) allows the reliable identification of individual motor unit (MU) action potentials. Despite the accuracy in decomposition, there is ...a large variability in the number of identified MUs across individuals and exerted forces. Here we present a systematic investigation of the anatomical and neural factors that determine this variability.
Approach
. We investigated factors of influence on HD-sEMG decomposition, such as synchronization of MU discharges, distribution of MU territories, muscle-electrode distance (MED—subcutaneous adipose tissue thickness), maximum anatomical cross-sectional area (ACSA
max
), and fiber cross-sectional area. For this purpose, we recorded HD-sEMG signals, ultrasound and magnetic resonance images, and took a muscle biopsy from the biceps brachii muscle from 30 male participants drawn from two groups to ensure variability within the factors—untrained-controls (UT = 14) and strength-trained individuals (ST = 16). Participants performed isometric ramp contractions with elbow flexors (at 15%, 35%, 50% and 70% maximum voluntary torque—MVT). We assessed the correlation between the number of accurately detected MUs by HD-sEMG decomposition and each measured parameter, for each target force level. Multiple regression analysis was then applied.
Main results.
ST subjects showed lower MED (UT = 5.1 ± 1.4 mm; ST = 3.8 ± 0.8 mm) and a greater number of identified MUs (UT: 21.3 ± 10.2 vs ST: 29.2 ± 11.8 MUs/subject across all force levels). The entire cohort showed a negative correlation between MED and the number of identified MUs at low forces (
r
= −0.6,
p
= 0.002 at 15% MVT). Moreover, the number of identified MUs was positively correlated to the distribution of MU territories (
r
= 0.56,
p
= 0.01) and ACSA
max
(
r
= 0.48,
p
= 0.03) at 15% MVT. By accounting for all anatomical parameters, we were able to partly predict the number of decomposed MUs at low but not at high forces.
Significance.
Our results confirmed the influence of subcutaneous tissue on the quality of HD-sEMG signals and demonstrated that MU spatial distribution and ACSA
max
are also relevant parameters of influence for current decomposition algorithms.
Objective. A signal-based metric for assessment of accuracy of motor unit (MU) identification from high-density surface electromyograms (EMG) is introduced. This metric, so-called ...pulse-to-noise-ratio (PNR), is computationally efficient, does not require any additional experimental costs and can be applied to every MU that is identified by the previously developed convolution kernel compensation technique. Approach. The analytical derivation of the newly introduced metric is provided, along with its extensive experimental validation on both synthetic and experimental surface EMG signals with signal-to-noise ratios ranging from 0 to 20 dB and muscle contraction forces from 5% to 70% of the maximum voluntary contraction. Main results. In all the experimental and simulated signals, the newly introduced metric correlated significantly with both sensitivity and false alarm rate in identification of MU discharges. Practically all the MUs with PNR > 30 dB exhibited sensitivity >90% and false alarm rates <2%. Therefore, a threshold of 30 dB in PNR can be used as a simple method for selecting only reliably decomposed units. Significance. The newly introduced metric is considered a robust and reliable indicator of accuracy of MU identification. The study also shows that high-density surface EMG can be reliably decomposed at contraction forces as high as 70% of the maximum.
Quercetin is one type of ergogenic aid and its effects on the neuromuscular system have recently attracted interest, but its dose-effect is not yet fully understood. The aim of this study was to ...examine the effect of different doses of quercetin ingestion on motor unit firing patterns and muscle contractile properties in humans. Thirteen young males and females conducted neuromuscular performance tests before (PRE) and 60 min after (POST) ingestions of 500 or 200 mg of quercetin glycosides (Qg500/Qg200, respectively) or placebo (PLA) on three different days. At PRE and POST, motor unit firing rates were calculated from high-density surface electromyography of the vastus lateralis muscle during 120-s isometric contraction of knee extension at 10% of maximal voluntary contraction. Electrically elicited forces in knee extensor muscles were also measured. After 60 s of voluntary contraction, motor unit firing rates, normalized by the exerted muscle force at POST, were significantly lower at POST than PRE with Qg500 and Qg200 (
p
< 0.05), but not with PLA (
p
> 0.05). Changes in motor unit firing rates normalized by the exerted force from PRE to POST were significantly greater with Qg500 than Qg200 at the end of contraction (
p
< 0.05). Under all three conditions, the electrically elicited force did not significantly change from PRE to POST (
p
> 0.05). These results suggest that both 500 and 200-mg quercetin ingestions alter motor unit firing patterns, and that quercetin’s effect is at least partially dose-dependent.
We developed and tested the methodology that supports the identification of individual motor unit (MU) firings from the Hoffman (or H) reflex recorded by surface high-density EMG (HD-EMG). Synthetic ...HD-EMG signals were constructed from simulated 10% to 90% of maximum voluntary contraction (MVC), followed by 100 simulated H-reflexes. In each H-reflex the MU firings were normally distributed with mean latency of 20 ms and standard deviations (SDLAT) ranging from 0.1 to 1.3 ms. Experimental H-reflexes were recorded from the soleus muscle of 12 men (33.6 ± 5.8 years) using HD-EMG array of <inline-formula> <tex-math notation="LaTeX">{5}\times {13} </tex-math></inline-formula> surface electrodes. Participants performed 15 to 20 s long voluntary plantarflexions with contraction levels ranging from 10% to 70% MVC. Afterwards, at least 60 H-reflexes were electrically elicited at three levels of background muscle activity: rest, 10% and 20% MVC. HD-EMGs of voluntary contractions were decomposed using the Convolution Kernel Compensation method to estimate the MU filters. When applied to HD-EMG signals with synthetic H reflexes, MU filters demonstrated high MU identification accuracy, especially for <inline-formula> <tex-math notation="LaTeX">\text {SD}_{\text {LAT}}>{0.3} </tex-math></inline-formula> ms. When applied to experimental H-reflex recordings, the MU filters identified 14.1 ± 12.1, 18.2 ± 12.1 and 20.8 ± 8.7 firings per H-reflex, with individual MU firing latencies of 35.9 ± 3.3, 35.1 ± 3.0 and 34.6 ± 3.3 ms for rest, 10% and 20% MVC background muscle activity, respectively. Standard deviation of MU latencies across experimental H-reflexes were 1.0 ± 0.8, 1.3 ± 1.1 and 1.5 ± 1.2 ms, in agreement with intramuscular EMG studies.
Quercetin is a polyphenolic flavonoid that has reported to block the binding of adenosine to A1 receptors at central nervous system and increase calcium release from the sarcoplasmic reticulum at ...skeletal muscle. The aim of the present study was to investigate the acute effect of quercetin ingestion on motor unit activation and muscle contractile properties. High-density surface electromyography during submaximal contractions and electrically elicited contraction torque in knee extensor muscles were measured before (PRE) and 60 min after (POST) quercetin glycosides or placebo ingestions in 13 young males. Individual motor units of the vastus lateralis muscle were identified from high-density surface electromyography by the Convolution Kernel Compensation technique. Firing rates of motor units recruited at 30–50% of the maximal voluntary contraction torque (MVC) were increased from PRE to POST only with quercetin (9.0 ± 2.3 to 10.5 ± 2.0 pps,
p
= 0.034). Twitch torque during doublet stimulation was decreased from PRE to POST with placebo (77.1 ± 17.1 to 73.9 ± 17.6 Nm,
p
= 0.005), but not with quercetin (
p
> 0.05). For motor units recruited at < 10% of MVC, normalized firing rate were decreased with quercetin (1.52 ± 0.33 to 1.58 ± 0.35%MVC/pps,
p
= 0.002) but increased with placebo (1.61 ± 0.32 to 1.57 ± 0.31%MVC/pps,
p
= 0.005). These results suggest that ingested quercetin has the functional roles to: mitigate reduction in the muscle contractile properties, enhance activations of relatively higher recruitment threshold motor units, and inhibit activation of relatively lower recruitment threshold motor units.
We introduce an algorithm for automatic identification of true positive (TP) and false positive (FP) spikes in the motor unit spike train, identified by blind source separation (BSS) of high-density ...surface electromyograms (HDsEMG). The algorithm selects predefined number of spikes, so called witnesses, from identified spike train. The other spikes in the spike train are called test spikes and are classified into TP or FP spikes by our algorithm. For this purpose, the algorithm constructs as many motor unit filters as there are test spikes, using the information from all the witnesses and each individual test spike. Afterwards, it applies each motor unit filter to HDsEMG to get new estimate of MU spike train for each selected test spike and calculates previously introduced Pulse-to-Noise Ratio (PNR) on preselected witnesses in this new spike train. When accumulated over all the test spikes, these PNR values exhibit bimodal distribution with the peak at lower PNR values representing FPs and the peak at higher PNR values representing TPs. Therefore, FPs and TPs can be discriminated by applying computationally efficient segmentation algorithm to corresponding PNR values. We also propose and mutually compare different witness selection strategies and show that selection of about 40 spikes with maximal amplitude in the identified spike train minimizes the selection of FPs as witnesses and maximizes the TP vs. FP discrimination power. In our tests on 20 s long experimental HDsEMG signals from biceps brachii muscle the number of FPs decreased from 23.9 ± 4.7 to 4.1 ± 4.4 when the proposed algorithm was used.
We compared non-negative matrix factorization (NMF) and convolution kernel compensation techniques for high-density electromyogram decomposition. The experimental data were recorded from nine healthy ...persons during controlled single degree of freedom (DOF) wrist flexion-extension, supination-pronation, and ulnar-radial deviation movements. We assembled the identified motor units and NMF components into three groups. Those active mostly during the first and the second movement direction per DOF were placed in the G1 and G3 groups, respectively. The remaining components were nonspecific for movement direction and were placed in the G2 group. In ulnar and radial deviation, the relative energies of identified cumulative motor unit spike trains (CSTs) and NMF components were similarly distributed among the groups. In other two movement types, the energy of NMF components in the G2 group was significantly larger than the energy of CSTs. We further performed a coherence analysis between CSTs and sums of NMF components in each group. Both decompositions demonstrated a solid match, but only at frequencies <;3 Hz. At higher frequencies, the coherence hardly exceeded the value of 0.5. Potential reasons for these discrepancies include the negative impact of motor unit action potential shapes and noise on NMF decomposition.
Nutritional supplementation in conjunction with exercise is of interest for the prevention or improvement of declines in motor performances in older adults. An understanding of the effects on both ...young and older adults contributes to its effective application. We investigated the effect of fish protein ingestion with resistance training on neural and muscular adaptations in young adults using interventions and assessments that have already been tested in older adults.
Eighteen young adults underwent 8 weeks of isometric knee extension training. During the intervention, nine participants ingested 5 g of fish protein (
= 9, Alaska pollack protein, APP), and the other nine participants ingested casein as a control (
= 9, CAS) in addition to daily meals. Before, during, and after the intervention, the isometric knee extension force, lower extremity muscle mass, and motor unit firing pattern of knee extensor muscles were measured.
Maximum voluntary contraction (MVC) was significantly increased in both APP and CAS groups from 0 weeks to 4, 6, and 8 weeks of intervention (
< 0.001), but there were no significant differences between the groups (
= 0.546-0.931). Muscle mass was not significantly changed during the intervention in either group (
= 0.250-0.698). Significant changes in motor unit firing rates (
= 0.02 and 0.029 for motor units recruited at 20-40% of MVC and at 40-60%) were observed following the intervention in the APP but not CAS (
= 0.120-0.751) group.
These results suggest that dietary fish protein ingestion changes motor unit adaptations following resistance training in young adults.
We study the changes of Motor Unit (MU) filters in MU identification from high-density surface electromyograms recorded during isokinetic dynamic contractions of biceps brachii muscle. We demonstrate ...that these changes can be predicted for limited changes of the joint angle by linearly extrapolating previously recorded changes, allowing for the linear prediction-correction paradigm of MU filter updating. We then demonstrate the efficiency of this paradigm by implementing MU filter updating by the Kalman filter and integrating it into the previously published Convolution Kernel Compensation (CKC) MU identification method. When compared to the original CKC method and the previously published cyclostationary CKC method devoted to MU identification in repeated dynamic contractions, the Kalman based MU filter prediction yielded a superior precision of MU firing tracking in dynamic contractions. In the case of relatively fast biceps brachii contractions with full elbow flexion and extension in 2s, the Kalman based MU filter prediction tracked 21.3 ± 1.8 MUs with an average sensitivity of 95.6 ± 7.0% and precision of 96.5 ± 3.5%. In the same conditions, the original CKC method identified 7.1 ± 2.0 MUs with an average sensitivity of 62.7 ± 20.1% and precision of 98.1 ± 3.9%, whereas cyclostationary CKC tracked 18.9 ± 2.0 MUs with an average sensitivity of 91.8 ± 12.2% and precision of 94.7 ± 5.1%.
Despite the progress in understanding of neural codes, the studies of the cortico-muscular coupling still largely rely on interferential electromyographic (EMG) signal or its rectification for the ...assessment of motor neuron pool behavior. This assessment is non-trivial and should be used with precaution. Direct analysis of neural codes by decomposing the EMG, also known as neural decoding, is an alternative to EMG amplitude estimation. In this study, we propose a fully-deterministic hybrid surface EMG (sEMG) decomposition approach that combines the advantages of both template-based and Blind Source Separation (BSS) decomposition approaches, a.k.a. guided source separation (GSS), to identify motor unit (MU) firing patterns. We use the single-pass density-based clustering algorithm to identify possible cluster representatives in different sEMG channels. These cluster representatives are then used as initial points of modified gradient Convolution Kernel Compensation (gCKC) algorithm. Afterwards, we use the Kalman filter to reduce the noise impact and increase convergence rate of MU filter identification by gCKC. Moreover, we designed an adaptive soft-thresholding method to identify MU firing times out of estimated MU spike trains. We tested the proposed algorithm on a set of synthetic sEMG signals with known MU firing patterns. A grid of 9 × 10 monopolar surface electrodes with 5-mm inter-electrode distances in both directions was simulated. Muscle excitation was set to 10, 30, and 50%. Colored Gaussian zero-mean noise with the signal-to-noise ratio (SNR) of 10, 20, and 30 dB, respectively, was added to 16 s long sEMG signals that were sampled at 4,096 Hz. Overall, 45 simulated signals were analyzed. Our decomposition approach was compared with gCKC algorithm. Overall, in our algorithm, the average numbers of identified MUs and Rate-of-Agreement (RoA) were 16.41 ± 4.18 MUs and 84.00 ± 0.06%, respectively, whereas the gCKC identified 12.10 ± 2.32 MUs with the average RoA of 90.78 ± 0.08%. Therefore, the proposed GSS method identified more MUs than the gCKC, with comparable performance. Its performance was dependent on the signal quality but not the signal complexity at different force levels. The proposed algorithm is a promising new offline tool in clinical neurophysiology.