Biomedical subjects

Dario Farina

Publications and source records attributed to Dario Farina.

At least 19 recordsLinked to original sources

PubMed2006

Neuromuscular adaptation in experimental and clinical neck pain.

The purpose of this brief review is to present evidence from experimental and clinical neck pain studies of pain-induced neuromuscular adaptations. It has been shown that clinical neck pain is associated with a substantial reorganization in the control strategies of cervical muscles during static and dynamic tasks. Experimental neck pain models allow local elicitation of nociceptive afferents, mimicking the sensory aspects of clinical pain, without major changes in muscle properties. These models may help understand the physiological mechanisms underlying the observations from clinical neck pain studies. The knowledge obtained from the interpretation of clinical findings with experimental pain models has relevance for the development of therapeutic interventions for the rehabilitation of patients with neck pain disorders.

Adaptation, Physiological↗

PubMed2006

The recruitment order of electrically activated motor neurons investigated with a novel collision technique.

OBJECTIVE: The development of a novel collision technique for assessment of the activation order of electrically activated nerve fibers, which is an important question in functional electrical therapy or for interpretation of results of motor unit number estimates. METHODS: Compound muscle action potentials were recorded with the belly-tendon configuration from the abductor digiti minimi. A novel modified Hopf's collision technique was applied on ten healthy male subjects to determine the distributions of conduction velocities (DCV) of all ulnar nerve fibers and of the fibers activated by electrical stimuli eliciting 20%, 50%, and 80% of the maximal muscle response. RESULTS: The maximum nerve conduction velocity was (means+/-SE) 64.1+/-0.85m/s. The median conduction velocity of estimated DCV was 58.9+/-0.97m/s (stimulus at 20%), 58.0+/-0.98m/s (50%), 57.2+/-0.91m/s (80%), and 56.5+/-0.84m/s (whole nerve) (all different between each other, P<0.001). CONCLUSIONS: The proposed collision technique allows the assessment of nerve conduction velocity distributions at maximal and sub-maximal stimulation levels and provided evidence for an inverse activation order of nerve fibers with electrical stimulation. SIGNIFICANCE: The excessive fatigue seen with nerve electrical stimulation can be explained by a preferential activation of large diameter nerve fibers. The motor units first activated with electrical stimulation are likely not representative of the motor unit pool in the muscle, which poses limitations in the reliability of some of the proposed methods for motor unit counting.

Action Potentials↗

PubMed2006

Estimation of average muscle fiber conduction velocity from simulated surface EMG in pinnate muscles.

The aim of this simulation study was to assess the bias in estimating muscle fiber conduction velocity (CV) from surface electromyographic (EMG) signals in muscles with one and two pinnation angles. The volume conductor was a layered medium simulating anisotropic muscle tissue and isotropic homogeneous subcutaneous tissue. The muscle tissue was homogeneous for one pinnation angle and inhomogeneous for bipinnate muscles (two fiber directions). Interference EMG signals were obtained by simulating recruitment thresholds and discharge patterns of a set of 100 and 200 motor units for the pinnate and bipinnate muscle, respectively (15 degrees pinnation angel in both cases). Without subcutaneous layer and muscle fibers with CV 4m/s, average CV estimates from the pinnate (bipinnate) muscle were 4.81+/-0.18 m/s (4.80+/-0.18 m/s) for bipolar, 4.71+/-0.19 m/s (4.71+/-0.12 m/s) for double differential, and 4.78+/-0.16 m/s (4.79+/-0.15m/s) for Laplacian recordings. When subcutaneous layer was added (thickness 1mm) in the same conditions, estimated CV values were 4.93+/-0.25 m/s (5.16+/-0.41 m/s), 4.70+/-0.21 m/s (4.83+/-0.33 m/s), and 4.89+/-0.21 m/s (4.99+/-0.39 m/s), for the three recording systems, respectively. The main factor biasing CV estimates was the propagation of action potentials in the two directions which influenced the recording due to the scatter of the projection of end-plate and tendon locations along the fiber direction, as a consequence of pinnation. The same problem arises in muscles with the line of innervation zone locations not perpendicular to fiber direction. These results indicate an important limitation in reliability of CV estimates from the interference EMG when the innervation zone and tendon locations are not distributed perpendicular to fiber direction.

Action Potentials↗

PubMed2006

Sensitivity of the cross-correlation between simulated surface EMGs for two muscles to detect motor unit synchronization.

The purpose of the study was to evaluate the use of cross-correlation analysis between simulated surface electromyograms (EMGs) of two muscles to quantify motor unit synchronization. The volume conductor simulated a cylindrical limb with two muscles and bone, fat, and skin tissues. Models of two motor neuron pools were used to simulate 120 s of surface EMG that were detected over both muscles. Short-term synchrony was established using a phenomenological model that aligned the discharge times of selected motor units within and across muscles to simulate physiological levels of motor unit synchrony. The correlation between pairs of surface EMGs was estimated as the maximum of the normalized cross-correlation function. After imposing four levels of motor unit synchrony across muscles, five parameters were varied concurrently in the two muscles to examine their influence on the correlation between the surface EMGs: 1) excitation level (5, 10, 15, and 50% of maximum); 2) muscle size (350 and 500 motor units); 3) fat thickness (1 and 4 mm); 4) skin conductivity (0.1 and 1 S/m); and 5) mean motor unit conduction velocity (2.5 and 4 m/s). Despite a constant and high level of motor unit synchronization among pairs of motor units across the two muscles, the cross-correlation index ranged from 0.08 to 0.56, with variation in the five parameters. For example, cross-correlation of EMGs from pairs of hand muscles, each having thin layers of subcutaneous fat and mean motor unit conduction velocities of 4 m/s, may be relatively insensitive to the level of synchronization across muscles. In contrast, cross-correlation of EMGs from pairs of leg muscles, with larger fat thickness, may exhibit a different sensitivity. These results indicate that cross correlation of the surface EMGs from two muscles provides a limited measure of the level of synchronization between motor units in the two muscles.

Computer Simulation↗

PubMed2006

Experimental muscle pain results in reorganization of coordination among trapezius muscle subdivisions during repetitive shoulder flexion.

The aim of the study was to examine the effect of experimental unilateral upper trapezius muscle pain on the relative activation of trapezius muscle subdivisions bilaterally during repetitive movement of the upper limb. Surface EMG signals were detected from nine healthy subjects from the upper, middle and lower divisions of trapezius during a repetitive bilateral shoulder flexion task. Measurements were performed before and after injection of 0.5 ml hypertonic (pain condition) and isotonic (control) saline into the upper division of the right trapezius muscle in two experimental sessions. On the painful side, upper trapezius showed decreased EMG amplitude (average rectified value, ARV) and lower trapezius increased ARV throughout the entire task following the injection of hypertonic saline (40.0 +/- 22.2 vs. 26.0 +/- 17.4 microV, and 12.5 +/- 7.6 vs. 25.6 +/- 14.8 microV, respectively, at the beginning of the contraction). On the side contralateral to pain, greater estimates of ARV were identified for the upper division of trapezius as the task progressed (37.4 +/- 20.2 vs. 52.7 +/- 28.4 microV, at the end of the contraction). Muscle fiber conduction velocity did not change with pain in all three divisions of the right trapezius muscle. The results suggest that local elicitation of nociceptive afferents in the upper division of the trapezius induces reorganization in the coordinated activity of the three subdivisions of the trapezius in repetitive dynamic tasks.

Adaptation, Physiological↗

PubMed2006

The change in spatial distribution of upper trapezius muscle activity is correlated to contraction duration.

The aim of the study was to confirm the hypothesis that the longer a contraction is sustained, the larger are the changes in the spatial distribution of muscle activity. For this purpose, surface electromyographic (EMG) signals were recorded with a 13 x 5 grid of electrodes from the upper trapezius muscle of 11 healthy male subjects during static contractions with shoulders 90 degrees abducted until endurance. The entropy (degree of uniformity) and center of gravity of the EMG root mean square map were computed to assess spatial inhomogeneity in muscle activation and changes over time in EMG amplitude spatial distribution. At the endurance time, entropy decreased (mean+/-SD, percent change 2.0+/-1.6%; P<0.0001) and the center of gravity moved in the cranial direction (shift 11.2+/-6.1mm; P<0.0001) with respect to the beginning of the contraction. The shift in the center of gravity was positively correlated with endurance time (R(2)=0.46, P<0.05), thus subjects with larger shift in the activity map showed longer endurance time. The percent variation in average (over the grid) root mean square was positively correlated with the shift in the center of gravity (R(2)=0.51, P<0.05). Moreover, the shift in the center of gravity was negatively correlated to both initial and final (at the endurance) entropy (R(2)=0.54 and R(2)=0.56, respectively; P<0.01 in both cases), indicating that subjects with less uniform root mean square maps had larger shift of the center of gravity over time. The spatial changes in root mean square EMG were likely due to spatially-dependent changes in motor unit activation during the sustained contraction. It was concluded that the changes in spatial muscle activity distribution play a role in the ability to maintain a static contraction.

Adult↗

PubMed2006

Experimental muscle pain changes the spatial distribution of upper trapezius muscle activity during sustained contraction.

OBJECTIVE: To investigate the effect of local excitation of nociceptive muscle afferents on the spatial distribution of muscle activity. METHODS: Surface electromyographic (EMG) signals were recorded from the upper trapezius muscle of 10 healthy volunteers with a 5 x 13 electrode grid during 90-s isometric contractions before, during, 15 and 30 min after intramuscular injection of hypertonic (painful) or isotonic (non-painful) saline. From the multi-channel EMG recordings, two-dimensional maps of root mean square and mean power frequency were obtained. The centre of gravity of the root mean square map was used to quantify global changes in the spatial distribution of muscle activity. RESULTS: During sustained contractions, average root mean square increased, average mean frequency decreased and the centre of gravity moved cranially. During experimental muscle pain, compared to before injection, the average root mean square decreased and there was a caudal shift of the centre of gravity. Fifteen minutes after the painful injection the centre of gravity returned to its original position. CONCLUSIONS: Short-term dynamic reorganization of the spatial distribution of muscle activity occurred in response to nociceptive afferent input. SIGNIFICANCE: The study furnishes an extension of the pain adaptation model indicating heterogeneous inhibition of muscle activity.

Adult↗

PubMed2006

Two-dimensional spatial distribution of surface mechanomyographical response to single motor unit activity.

In order to better understand the mechanisms of generation of mechanomyography (MMG) signals, the two-dimensional distribution of surface MMG produced by the activity of single motor units was analyzed by a novel two-dimensional recording method. Motor unit action potentials were identified from intramuscular electromyographic (EMG) signals and used to trigger the averaging of MMG signals detected over the tibialis anterior muscle of 11 volunteers with a grid of 5x3 accelerometers (20-mm inter-accelerometer distance). The intramuscular wires were inserted between the first and second accelerometer in the middle column of the grid, proximal to the innervation zone. The subjects performed three contractions with visual feedback of the intramuscular EMG signals. In each contraction, a new motor unit was recruited at the minimum stable discharge rate (mean+/-S.D., N = 11 subjects, 7.3+/-2.3 pulse/s), resulting in torque of 2.4+/-2.8% of the maximal voluntary contraction (MVC), 4.6+/-2.7% MVC, and 6.3+/-3.1% MVC (all different, P < 0.01). For 23 out of 33 detected motor units, it was possible to extract the motor unit surface acceleration map (MUAM). A negative MUAM peak (-2.7+/-2.2 mm/s2) was detected laterally and a positive MUAM peak (4.1+/-2.4 mm/s2) medially (P < 0.001). The time-to-peak was shorter in the medial part of the muscle (2.9+/-0.4 ms) than in the other locations (3.4+/-0.5 ms, P < 0.001). The double integrated signals (muscle displacement) indicated negative deflection in the lateral part and inflation close to the tibia bone. The maps of acceleration showed spatial dependency in single motor unit MMG activities. The technique provides a new insight into motor unit contractile properties.

Acceleration↗

PubMed2006

Spatial and force dependency of mechanomyographic signal features.

The aim was to investigate with a novel technique the spatial inhomogeneity in surface mechanomyographic (MMG) response to muscle contraction at varying force levels. MMG signals were detected over the dominant tibialis anterior muscle of 10 volunteers using a 5 x 3 grid of accelerometers. The subjects performed 3 s long isometric contractions at forces ranging from 0% to 100% of the maximal force (10% increment) in a randomised order. From the two-dimensional MMG recordings, maps of absolute and normalized temporal and spectral MMG descriptors were obtained. The centroid and entropy of these maps were computed to describe the spatial centre of activity and degree of homogeneity, respectively. MMG absolute amplitude did not depend on location over the muscle while normalized amplitude did and the centroid shifted with increasing force. Amplitude increased with force and its entropy decreased. Absolute and normalized spectral variables depended on location over the muscle and their centroid shifted with increasing force. In addition, the dependency of absolute and normalized spectral variables on force was affected by location. These results highlight limitations when using single-channel MMG features for the assessment of motor unit control strategies, due to a substantial effect of position on the relation between force and MMG characteristics.

Acceleration↗

PubMed2006

Optimal wavelets for biomedical signal compression.

Signal compression is gaining importance in biomedical engineering due to the potential applications in telemedicine. In this work, we propose a novel scheme of signal compression based on signal-dependent wavelets. To adapt the mother wavelet to the signal for the purpose of compression, it is necessary to define (1) a family of wavelets that depend on a set of parameters and (2) a quality criterion for wavelet selection (i.e., wavelet parameter optimization). We propose the use of an unconstrained parameterization of the wavelet for wavelet optimization. A natural performance criterion for compression is the minimization of the signal distortion rate given the desired compression rate. For coding the wavelet coefficients, we adopted the embedded zerotree wavelet coding algorithm, although any coding scheme may be used with the proposed wavelet optimization. As a representative example of application, the coding/encoding scheme was applied to surface electromyographic signals recorded from ten subjects. The distortion rate strongly depended on the mother wavelet (for example, for 50% compression rate, optimal wavelet, mean+/-SD, 5.46+/-1.01%; worst wavelet 12.76+/-2.73%). Thus, optimization significantly improved performance with respect to previous approaches based on classic wavelets. The algorithm can be applied to any signal type since the optimal wavelet is selected on a signal-by-signal basis. Examples of application to ECG and EEG signals are also reported.

Algorithms↗

PubMed2006

Correlation of average muscle fiber conduction velocity measured during cycling exercise with myosin heavy chain composition, lactate threshold, and VO2max.

The aim of the study was to investigate the correlation between myosin heavy chain (MHC) composition, lactate threshold (LT), maximal oxygen uptake VO2max, and average muscle fiber conduction velocity (MFCV) measured from surface electromyographic (EMG) signals during cycling exercise. Ten healthy male subjects participated in the study. MHC isoforms were identified from a sample of the vastus lateralis muscle and characterized as type I, IIA, and IIX. At least three days after a measure of LT and VO2max, the subjects performed a 2-min cycling exercise at 90 revolutions per minute and power output corresponding to LT, during which surface EMG signals were recorded from the vastus lateralis muscle with an adhesive electrode array. MFCV and instantaneous mean power spectral frequency of the surface EMG were estimated at the maximal instantaneous knee angular speed. Output power corresponding to LT and VO2max were correlated with percentage of MHC I (R2=0.77; and 0.42, respectively; P<0.05). MFCV was positively correlated with percentage of MHC I, power corresponding to LT and to VO2max (R2=0.84; 0.74; 0.53, respectively; P<0.05). Instantaneous mean power spectral frequency was not correlated with any of these variables or with MFCV, thus questioning the use of surface EMG spectral analysis for indirect estimation of MFCV in dynamic contractions.

Adult↗

PubMed2006

Compression of surface EMG signals with algebraic code excited linear prediction.

Despite the interest in long timescale recordings of surface electromyographic (EMG) signals, only a few studies have focused on EMG compression. In this paper we investigate a lossy coding technique for surface EMG signals that is based on the algebraic code excited linear prediction (ACELP) paradigm, widely used for speech signal coding. The algorithm was adapted to the EMG characteristics and tested on both simulated and experimental signals. The coding parameters selected led to a compression ratio of 87.3%. For simulated signals, the mean square error in signal reconstruction and the percentage error in average rectified value after compression were 11.2% and 4.90%, respectively. For experimental signals, they were 6.74% and 3.11%. The mean power spectral frequency and third-order power spectral moment were estimated with relative errors smaller than 1.23% and 8.50% for simulated signals, and 3.74% and 5.95% for experimental signals. It was concluded that the proposed coding scheme could be effectively used for high rate and low distortion compression of surface EMG signals. Moreover, the method is characterized by moderate complexity (approximately 20 million instructions/s) and an algorithmic delay smaller than 160 samples (approximately 160ms).

Adult↗

PubMed2006

Spatial and temporal changes of upper trapezius muscle fiber conduction velocity are not predicted by surface EMG spectral analysis during a dynamic upper limb task.

The purpose of the study was to examine the temporal and spatial correlation between estimates of trapezius muscle fiber conduction velocity (CV) and surface EMG instantaneous mean power spectral frequency (iMPF) during dynamic movement of the upper limb. Surface EMG signals were detected from the upper division of the trapezius muscle in 13 healthy volunteers using linear arrays of eight electrodes at three locations in the cephalad-caudal direction. Subjects were asked to tap with their hands in a cyclic manner between targets positioned mid thigh and 120 degrees of shoulder flexion, to the beat of a metronome set at 88 beats per minute for 5 min. Muscle fiber CV and iMPF were estimated for each cycle at the time instant corresponding to 90 degrees of shoulder flexion. Non-significant correlations were identified between CV and iMPF initial values (R(2)=0.03-0.01), rate of change over time (R(2)=0.10-0.004) and normalized rate of change (R(2)=0.12-0.01) at all three locations on the upper trapezius muscle. These results demonstrate that both spatial and temporal variations in trapezius muscle fiber CV are not predicted by EMG spectral analysis during dynamic movement of the upper limb. This finding suggests that spectral analysis cannot be used to infer changes in the spatial and temporal behavior of muscle fiber CV during dynamic tasks.

Adult↗

PubMed2006

Spectral moments of mechanomyographic signals recorded with accelerometer and microphone during sustained fatiguing contractions.

The aim of this study was to analyse the trends of the first three power spectral moments of the mechanomyogram (MMG) signal recorded by a microphone (MMG(MIC)) and an accelerometer (MMG(ACC)) during sustained contractions. MMG signals were recorded from the biceps brachii muscle in 14 healthy male subjects during a 3 min isometric elbow flexion at 30% of the maximal voluntary contraction. MMG absolute and normalised root mean square (RMS), mean power frequency (MNF), power spectral variance (Mc2), and skewness (mu3) were computed. For both MMG(MIC) and MMG(ACC), absolute and normalised RMS and Mc2 increased while MNF and mu3 decreased with contraction time (P<0.001). The rates of change of RMS over time were significantly correlated (P<0.001) for MMG(MIC) and MMG(ACC) but not correlated for spectral moments. The coefficient of variation of RMS was higher for MMG(MIC) than for MMG(ACC), while the opposite was observed for mu3 (P<0.05). It was concluded that higher order spectral moments of the MMG signal change during sustained contraction, indicating a complex modification of the shape of the power spectrum and not just scaling of the bandwidth. This is most likely due to the additional motor unit recruitment with fatigue and to the non-linear summation of motor unit contributions to the signal. Moreover, the characteristics of MMG signals recorded with microphones and accelerometers have important differences, which should be taken into account when comparing results from different studies.

Adult↗

PubMed2006

Amplitude cancellation reduces the size of motor unit potentials averaged from the surface EMG.

The purpose of the study was to evaluate the influence of selected physiological parameters on amplitude cancellation in the simulated surface electromyogram (EMG) and the consequences for spike-triggered averages of motor unit potentials derived from the interference and rectified EMG signals. The surface EMG was simulated from prescribed recruitment and rate coding characteristics of a motor unit population. The potentials of the motor units were detected on the skin over a hand muscle with a bipolar electrode configuration. Averages derived from the EMG signal were generated using the discharge times for each of the 24 motor units with lowest recruitment thresholds from a population of 120 across three conditions: 1) excitation level; 2) motor unit conduction velocity; and 3) motor unit synchronization. The area of the surface-detected potential was compared with potentials averaged from the interference, rectified, and no-cancellation EMGs. The no-cancellation EMG comprised motor unit potentials that were rectified before they were summed, thereby preventing cancellation between the opposite phases of the potentials. The percent decrease in area of potentials extracted from the rectified EMG was linearly related to the amount of amplitude cancellation in the interference EMG signal, with the amount of cancellation influenced by variation in excitation level and motor unit conduction velocity. Motor unit synchronization increased potentials derived from both the rectified and interference EMG signals, although cancellation limited the increase in area for both potentials. These findings document the influence of amplitude cancellation on motor unit potentials averaged from the surface EMG and the consequences for using the procedure to characterize motor unit properties.

Action Potentials↗

PubMed2006

Muscle-fiber conduction velocity during concentric and eccentric actions on a flywheel exercise device.

A gravity-independent flywheel exercise device (FWED) has been proven effective as a countermeasure to loss of strength and muscle atrophy induced by simulated microgravity. This study assessed muscle-fiber conduction velocity (CV) and surface EMG instantaneous mean power spectral frequency (iMNF) during brief bouts of fatiguing concentric (CON) and eccentric (ECC) exercise on a FWED in order to identify electromyographic (EMG) variables that can be used to provide objective indications of muscle status when exercising with a FWED. Multichannel surface EMG signals were recorded from vastus lateralis and medialis muscles of nine men during: (1) isometric, 60-s action at 50% of maximum voluntary action (MVC); (2) two isometric, linearly increasing force ramps (0-100% MVC); and (3) dynamic CON/ECC coupled actions on the FWED. Muscle-fiber CV and iMNF were computed over time during the three tasks. During ramps, CV, but not iMNF, increased with force (P < 0.001). Conduction velocity and iMNF decreased with the same normalized rate of change in constant-force actions. During CON/ECC actions, the normalized rate of change over time was larger for CV than iMNF (P < 0.05). These results suggest that, during fatiguing, dynamic, variable-force tasks, changes in CV cannot be indirectly inferred by EMG spectral analysis. This underlines the importance of measuring both CV and spectral variables for muscle assessment in dynamic tasks.

Adult↗

PubMed2006

A finite element model for describing the effect of muscle shortening on surface EMG.

A finite-element model for the generation of single fiber action potentials in a muscle undergoing various degrees of fiber shortening is developed. The muscle is assumed fusiform with muscle fibers following a curvilinear path described by a Gaussian function. Different degrees of fiber shortening are simulated by changing the parameters of the fiber path and maintaining the volume of the muscle constant. The conductivity tensor is adapted to the muscle fiber orientation. In each point of the volume conductor, the conductivity of the muscle tissue in the direction of the fiber is larger than that in the transversal direction. Thus, the conductivity tensor changes point-by-point with fiber shortening, adapting to the fiber paths. An analytical derivation of the conductivity tensor is provided. The volume conductor is then studied with a finite-element approach using the analytically derived conductivity tensor. Representative simulations of single fiber action potentials with the muscle at different degrees of shortening are presented. It is shown that the geometrical changes in the muscle, which imply changes in the conductivity tensor, determine important variations in action potential shape, thus affecting its amplitude and frequency content. The model provides a new tool for interpreting surface EMG signal features with changes in muscle geometry, as it happens during dynamic contractions.

Action Potentials↗

PubMed2006

An analytical model for surface EMG generation in volume conductors with smooth conductivity variations.

A nonspace invariant model of volume conductor for surface electromyography (EMG) signal generation is analytically investigated. The volume conductor comprises planar layers representing the muscle and subcutaneous tissues. The muscle tissue is homogeneous and anisotropic while the subcutaneous layer is inhomogeneous and isotropic. The inhomogeneity is modeled as a smooth variation in conductivity along the muscle fiber direction. This may reflect a practical situation of tissues with different conductivity properties in different locations or of transitions between tissues with different properties. The problem is studied with the regular perturbation theory, through a series expansion of the electric potential. This leads to a set of Poisson's problems, for which the source term in an equation and the boundary conditions are determined by the solution of the previous equations. This set of problems can be solved iteratively. The solution is obtained in the two-dimensional Fourier domain, with spatial angular frequencies corresponding to the longitudinal and perpendicular direction with respect to the muscle fibers, in planes parallel to the detection surface. The series expansion is truncated for the practical implementation. Representative simulations are presented. The proposed model constitutes a new approach for surface EMG signal simulation with applications related to the validation of methods for information extraction from this signal.

Action Potentials↗