Mapping of spatial and temporal heterogeneity of plantar flexor muscle activity during isometric contraction: correlation of velocity-encoded MRI with EMG

Muscle fiber strain rates in the lower leg during ankle dorsi-/plantarflexion exercise

Static quantitative magnetic resonance imaging (MRI) provides readouts of structural changes in diseased muscle, but current approaches lack the ability to fully explain the loss of contractile function. Muscle contractile function can be assessed using various techniques including phase-contrast MRI (PC-MRI), where strain rates are quantified. However, current two-dimensional implementations are limited in capturing the complex motion of contracting muscle in the context of its three-dimensional (3D) fiber architecture. The MR acquisitions (chemical shift-encoded water-fat separation scan, spin echo-echoplanar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) wereperformed at 3 T. PC-MRI acquisitions and performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Acquisitions (3 T, chemical shift-encoded water-fat separation scan, spin echo-echo planar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) were performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Strain rates and diffusion tensors were calculated and combined to obtain strain rates along and perpendicular to the muscle fibers in seven lower leg muscles during the dynamic dorsi-/plantarflexion movement cycle. To evaluate strain rates along the proximodistal muscle axis, muscles were divided into five equal segments. t-tests were used to test if cyclic strain rate patterns (amplitude > 0) were present along and perpendicular to the muscle fibers. The effects of proximal-distal location and load were evaluated using repeated measures ANOVAs. Cyclic temporal strain rate patterns along and perpendicular to the fiber were found in all muscles involved in dorsi-/plantarflexion movement (p < 0.0017). Strain rates along and perpendicular to the fiber were heterogeneously distributed over the length of most muscles (p < 0.003). Additional loading reduced strain rates of the extensor digitorum longus and gastrocnemius lateralis muscle (p < 0.001). In conclusion, the lower leg muscles involved in cyclic dorsi-/plantarflexion exercise showed cyclic fiber strain rate patterns with amplitudes that varied between muscles and between the proximodistal segments within the majority of muscles.

Broader Estimates of Gastrocnemius Activity Generated a More Representative Cocontraction Index: A Study in Pediatric Population

The electromyography (EMG) cocontraction index (CCI) given by the antagonistic/agonistic Root Mean Square (RMS) amplitude ratio of the same muscle is a qualified biomarker used for spastic cocontraction quantification and management in cerebral palsy children. However, this normative EMG ratio is likely subject to a potential source of errors with biased estimates when measuring the gastrocnemius plantar flexors activity. Due to the uneven distribution of electrical activity within the muscle volume, cocontraction levels can be misestimated, if EMGs are obtained from the sole traditional bipolar sensor location recommended by SENIAM. This preliminary study, on 10 healthy children (mean age 10 yr), investigated whether surface EMG detected proximally and distally via two pairs of bipolar electrodes, within the medial gastrocnemius (MG), provides a significant difference in CCI estimates during non-dynamic (isometric dorsiflexion) and dynamic (swing phases of gait) conditions. Gait cycles were extracted from Inertial Measurement Unit sensors. Medial gastrocnemius activity was greater distally than proximally during plantar flexion when it acts as an agonist (~24±18%) and it was greater proximally during dorsiflexion (~23±9%) when it is acting as an antagonist. As a direct consequence, CCI estimates from the conventional sensor location were significantly different (~36%) from the CCIs computed by considering broader MG regions. This difference arose in all subjects during isometric efforts and in two of 10 healthy children during the swing phase of gait who presented cocontraction patterns ( [Formula: see text]). EMG bipolar sampling encompassing proximal and distal gastrocnemius muscle regions may reduce bias in CCI computation and provide a more representative and accurate cocontraction index that is especially important for comparisons to the diseased state.

Linking cortex and contraction-Integrating models along the corticomuscular pathway

Computational models of the neuromusculoskeletal system provide a deterministic approach to investigate input-output relationships in the human motor system. Neuromusculoskeletal models are typically used to estimate muscle activations and forces that are consistent with observed motion under healthy and pathological conditions. However, many movement pathologies originate in the brain, including stroke, cerebral palsy, and Parkinson's disease, while most neuromusculoskeletal models deal exclusively with the peripheral nervous system and do not incorporate models of the motor cortex, cerebellum, or spinal cord. An integrated understanding of motor control is necessary to reveal underlying neural-input and motor-output relationships. To facilitate the development of integrated corticomuscular motor pathway models, we provide an overview of the neuromusculoskeletal modelling landscape with a focus on integrating computational models of the motor cortex, spinal cord circuitry, α-motoneurons  and skeletal muscle in regard to their role in generating voluntary muscle contraction. Further, we highlight the challenges and opportunities associated with an integrated corticomuscular pathway model, such as challenges in defining neuron connectivities, modelling standardisation, and opportunities in applying models to study emergent behaviour. Integrated corticomuscular pathway models have applications in brain-machine-interaction, education, and our understanding of neurological disease.

Dynamics of Quadriceps Muscles during Isometric Contractions: Velocity-Encoded Phase Contrast MRI Study

Objective: To quantify the spatial heterogeneity of displacement during voluntary isometric contraction within and between the different compartments of the quadriceps. Methods: The thigh muscles of seven subjects were imaged on an MRI scanner while performing isometric knee extensions at 40% maximal voluntary contraction. A gated velocity-encoded phase contrast MRI sequence in axial orientations yielded tissue velocity-encoded dynamic images of the four different compartments of the thigh muscles (vastus lateralis (VL), vastus medialis (VM), vastus intermedius (VI), and rectus femoris (RF)) at three longitudinal locations of the proximal–distal length: 17.5% (proximal), 50% (middle), and 77.5% (distal). The displacement, which is the time integration of the measured velocity, was calculated along the three orthogonal axes using a tracking algorithm. Results: The displacement of the muscle tissues was clearly nonuniform within each axial section as well as between the three axial locations. The ensemble average of the magnitude of the total displacement as a synthetic vector of the X, Y, and Z displacements was significantly larger in the VM at the middle location (p < 0.01), and in the VI at the distal location than in the other three muscles. The ensemble average of Z-axis displacement, which was almost aligned with the line of action, was significantly larger in VI than in the other three muscles in all three locations. Displacements of more than 20 mm were observed around the central aponeuroses, such as those between VI and the other surrounding muscles. Conclusions: These results imply that the quadriceps muscles act as one functional unit in normal force generation through the central aponeuroses despite complex behavior in each of the muscles, each of which possesses different physiological characteristics and architectures.

A split-label design for simultaneous measurements of perfusion in distant slices by pulsed arterial spin labeling

Purpose

Multislice arterial spin labeling (ASL) MRI acquisitions are currently challenging in skeletal muscle because of long transit times, translating into low‐perfusion SNR in distal slices when large spatial coverage is required. However, fiber type and oxidative capacity vary along the length of healthy muscles, calling for multislice acquisitions in clinical studies. We propose a new variant of flow alternating inversion recovery (FAIR) that generates sufficient ASL signal to monitor exercise‐induced perfusion changes in muscle in two distant slices.

Methods

Label around and between two 7‐cm distant slices was created by applying the presaturation/postsaturation and selective inversion modules selectively to each slice (split‐label multislice FAIR). Images were acquired using simultaneous multislice EPI. We validated our approach in the brain to take advantage of the high resting‐state perfusion, and applied it in the lower leg muscle during and after exercise, interleaved with a single‐slice FAIR as a reference.

Results

We show that standard multislice FAIR leads to an underestimation of perfusion, while the proposed split‐label multislice approach shows good agreement with separate single‐slice FAIR acquisitions in brain, as well as in muscle following exercise.

Conclusion

Split‐label FAIR allows measuring muscle perfusion in two distant slices simultaneously without losing sensitivity in the distal slice.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Stiffness mapping of lower leg muscles during passive dorsiflexion

It is challenging to differentiate the mechanical properties of synergist muscles in vivo. Shear wave elastography can be used to quantify the shear modulus (i.e. an index of stiffness) of a specific muscle. This study assessed the passive behavior of lower leg muscles during passive dorsiflexion performed with the knee fully extended (experiment 1, n = 22) or with the knee flexed at 90° (experiment 2, n = 20). The shear modulus measurements were repeated twice during experiment 1 to assess the inter-day reliability. During both experiments, the shear modulus of the following plantar flexors was randomly measured: gastrocnemii medialis (GM) and lateralis (GL), soleus (SOL), peroneus longus (PL), and the deep muscles flexor digitorum longus (FDL), flexor hallucis longus (FHL), tibialis posterior (TP). Two antagonist muscles tibialis anterior (TA), and extensor digitorum longus (EDL) were also recorded. Measurements were performed in different proximo-distal regions for GM, GL and SOL. Inter-day reliability was adequate for all muscles (coefficient of variation < 15%), except for TP. In experiment 1, GM exhibited the highest shear modulus at 80% of the maximal range of motion (128.5 ± 27.3 kPa) and was followed by GL (67.1 ± 24.1 kPa). In experiment 2, SOL exhibited the highest shear modulus (55.1 ± 18.0 kPa). The highest values of shear modulus were found for the distal locations of both the GM (80% of participants in experiment 1) and the SOL (100% of participants in experiment 2). For both experiments, deep muscles and PL exhibited low levels of stiffness during the stretch in young asymptomatic adults, which was unknown until now. These results provide a deeper understanding of passive mechanical properties and the distribution of stiffness between and within the plantar flexor muscles during stretching between them and thus could be relevant to study the effects of aging, disease progression, and rehabilitation on stiffness.

Behavior of human gastrocnemius muscle fascicles during ramped submaximal isometric contractions

Precise estimates of muscle architecture are necessary to understand and model muscle mechanics. The primary aim of this study was to estimate continuous changes in fascicle length and pennation angle in human gastrocnemius muscles during ramped plantar flexor contractions at two ankle angles. The secondary aim was to determine whether these changes differ between proximal and distal fascicles. Fifteen healthy subjects performed ramped contractions (0–25% MVC) as ultrasound images were recorded from the medial (MG, eight sites) and lateral (LG, six sites) gastrocnemius muscle with the ankle at 90° and 120° (larger angles correspond to shorter muscle lengths). In all subjects, fascicles progressively shortened with increasing torque. MG fascicles shortened 5.8 mm (11.1%) at 90° and 4.5 mm (12.1%) at 120°, whereas LG muscle fascicles shortened 5.1 mm (8.8%) at both ankle angles. MG pennation angle increased 1.4° at 90° and 4.9° at 120°, and LG pennation angle decreased 0.3° at 90° and increased 2.6° at 120°. Muscle architecture changes were similar in proximal and distal fascicles at both ankle angles. This is the first study to describe continuous changes in fascicle length and pennation angle in the human gastrocnemius muscle during ramped isometric contractions. Very similar changes occurred in proximal and distal muscle regions. These findings are relevant to studies modeling active muscle mechanics.

FDG-PET detects nonuniform muscle activity in the lower body during human gait

INTRODUCTION

Nonuniform muscle activity has been partially explained by anatomically defined neuromuscular compartments. The purpose of this study was to investigate the uniformity of skeletal muscle activity during walking.

METHODS

Eight participants walked at a self-selected speed, and muscle activity was quantified using [¹⁸ F]-fluorodeoxyglucose positron emission tomography imaging. Seventeen muscles were divided into 10 equal length sections, and within muscle activity was compared.

RESULTS

Nonuniform activity was detected in 12 of 17 muscles (ƒ > 4.074; P < 0.046), which included both uni- and multi-articular muscles. Greater proximal activity was detected in 6 muscles (P < 0.049), and greater distal versus medial activity was found in the iliopsoas (P < 0.042).

CONCLUSIONS

Nonuniform muscle activity is likely related to recruitment of motor units located within separate neuromuscular compartments. These findings indicate that neuromuscular compartments are recruited selectively to allow for efficient energy transfer, and these patterns may be task-dependent. Muscle Nerve 54: 959-966, 2016.