Neuromuscular neutral zones response to cyclic lumbar flexion

The in vivo lumbar spine of the anaesthetized feline was subjected to passive cyclic anterior flexion-extension at 0.25 Hz and 40 N peak load for cumulative 60 min duration. Displacement (or displacement neuromuscular neutral zones-DNNZ) and tension (or tension neuromuscular neutral zones-TNNZ) at which reflexive EMG activity from the multifidi muscles was initiated and terminated were recorded, for single-test cycles, before and for 7h after cyclic loading. Displacement and tension NNZs increased significantly after loading. The displacement NNZs decreased exponentially to near baseline by the 7th hour of rest. The tension NNZs, however, decreased to below the baseline by the 2nd to 3rd hour after loading and continued decreasing into the 7th hour. Peak EMG significantly decreased (49-57%) to below the baseline immediately after loading and then exponentially increased, exceeding the baseline by the 2nd to 3rd hour and reaching 33-59% above baseline by the 7th hour. EMG median frequency decreased after loading and then exceeded the baseline after the 3rd hour, indicating initial de-recruitment, followed by recruitment of new motor units. These findings suggest that the lumbar spine was exposed to instability for 2-3h after cyclic loading, due to concurrent laxity of the viscoelastic tissues and deficient muscular activity. A delayed neuromuscular compensation mechanism was found to exist, triggering the musculature significantly earlier and at higher magnitude than baseline, while the viscoelastic tissues were still lax. Thus, it is suggested that prolonged cyclic loading may compromise lumbar stability during the immediate 2-3h post-loading, increasing the risk of injury.     

Frequency-dependent changes in neuromuscular responses to cyclic lumbar flexion

Repetitive lifting in the workplace has been identified to be a cause of low back disorders. Epidemiologic data further supports an hypothesis that higher repetition rate (i.e. frequency) is an added risk factor. The objective of this study was to provide experimental data testing the above hypothesis. An in vivo feline model was subjected to 20-min of cyclic lumbar loading at frequencies of 0.1 Hz and 0.5 Hz while monitoring the EMG from the L-3/4-L-5/6 multifidus muscles and the creep at the L-4/5 level. Seven hours of rest were allowed after the cyclic flexion/extension was terminated. During this rest period, a single test cycle was performed every hour to assess recovery of EMG and lumbar creep. The results demonstrate that cyclic lumbar flexion elicits a transient neuromuscular disorder consisting of EMG spasms during the cyclic loading and initial and delayed muscular hyperexcitabilities during the rest period. Cyclic loading at 0.5 Hz resulted in significant (p<0.05) increase in the hyperexcitability magnitude and duration during the recovery period. It was concluded that repetitive lumbar loading at fast rates is indeed a risk factor as it induces larger creep in the lumbar viscoelastic tissues which in turn intensify the resulting neuromuscular disorder.     

Use of additive models to represent trends in a barley field trial

Even when a field experiment has been designed with care, subsequent examination of the plot values may reveal additional unforeseen trends. In this paper we examine data from a barley pathology field trial and show that additive models provide a flexible representation of environmental trends, in one or two dimensions. Such models smooth out noise in the observed data, rather than fit an equation specified in advance. This approach tends to increase the precision of treatment comparisons relative to a classical analysis of variance. We recommend the use of residual plots to explore experimental data for underlying trends, and additive models to display these trends and estimate treatment effects.     

Comparison between line and surface mesh models to represent the rotator cuff muscle geometry in musculoskeletal models

Accurate muscle geometry (muscle length and moment arm) is required to estimate muscle function when using musculoskeletal modelling. In shoulder, muscles are often modelled as a collection of independent line segments, leading to non-physiological muscles trajectory, especially for the rotator cuff muscles. To prevent this, a surface mesh model was developed and validated against 7 MRI positions in one participant. Mean moment arm errors was 11.4% for the line vs. 8.8% for the mesh model. While the model with independent lines led to some non-physiological trajectories, the mesh model gave lower misestimations of muscle lengths and moment arms.     

Fuzzy species distribution models: a way to represent plant communities spatially

Fuzzy set theory has generally been applied to smooth classification cut‐offs, with an unavoidable loss of information. In this commentary, I rely on both advantages and disadvantages of the methods proposed in Duff et al., in this issue of the Journal of Vegetation Science, to map the variability over space of vegetation classes based on fuzzy sets and species distribution models.     

Sensitivity of lumbar spine response to follower load and flexion moment: finite element study

The follower load (FL) combined with moments is commonly used to approximate flexed/extended posture of the lumbar spine in absence of muscles in biomechanical studies. There is a lack of consensus as to what magnitudes simulate better the physiological conditions. Considering the in-vivo measured values of the intradiscal pressure (IDP), intervertebral rotations (IVRs) and the disc loads, sensitivity of these spinal responses to different FL and flexion moment magnitudes was investigated using a 3D nonlinear finite element (FE) model of ligamentous lumbosacral spine. Optimal magnitudes of FL and moment that minimize deviation of the model predictions from in-vivo data were determined. Results revealed that the spinal parameters i.e. the IVRs, disc moment, and the increase in disc force and moment from neutral to flexed posture were more sensitive to moment magnitude than FL magnitude in case of flexion. The disc force and IDP were more sensitive to the FL magnitude than moment magnitude. The optimal ranges of FL and flexion moment magnitudes were 900–1100 N and 9.9–11.2 Nm, respectively. The FL magnitude had reverse effect on the IDP and disc force. Thus, magnitude for FL or flexion that minimizes the deviation of all the spinal parameters together from the in-vivo data can vary. To obtain reasonable compromise between the IDP and disc force, our findings recommend that FL of low magnitude must be combined with flexion moment of high intensity and vice versa.     

Position sense in the lumbar spine with torso flexion and loading

Proprioception plays an important role in appropriate sensation of spine position, movement, and stability. Previous research has demonstrated that position sense error in the lumbar spine is increased in flexed postures. This study investigated the change in position sense as a function of altered trunk flexion and moment loading independently. Reposition sense of lumbar angle in 17 subjects was assessed. Subjects were trained to assume specified lumbar angles using visual feedback. The ability of the subjects to reproduce this curvature without feedback was then assessed. This procedure was repeated for different torso flexion and moment loading conditions. These measurements demonstrated that position sense error increased significantly with the trunk flexion (40%, p < .05) but did not increase with moment load (p = .13). This increased error with flexion suggests a loss in the ability to appropriately sense and therefore control lumbar posture in flexed tasks. This loss in proprioceptive sense could lead to more variable lifting coordination and a loss in dynamic stability that could increase low back injury risk. This research suggests that it is advisable to avoid work in flexed postures.     

Conclusions on motor control depend on the type of model used to represent the periphery

Within the field of motor control, there is no consensus on which kinematic and kinetic aspects of movements are planned or controlled. Perturbing goal-directed movements is a frequently used tool to answer this question. To be able to draw conclusions about motor control from kinematic responses to perturbations, a model of the periphery (i.e., the skeleton, muscle–tendon complexes, and spinal reflex circuitry) is required. The purpose of the present study was to determine to what extent such conclusions depend on the level of simplification with which the dynamical properties of the periphery are modeled. For this purpose, we simulated fast goal-directed single-joint movement with four existing types of models. We tested how three types of perturbations affected movement trajectory if motor commands remained unchanged. We found that the four types of models of the periphery showed different robustness to the perturbations, leading to different predictions on how accurate motor commands need to be, i.e., how accurate the knowledge of external conditions needs to be. This means that when interpreting kinematic responses obtained in perturbation experiments the level of error correction attributed to adaptation of motor commands depends on the type of model used to describe the periphery.     

Muscle-driven and torque-driven centrodes during modeled flexion of individual lumbar spines are disparate

Lumbar spine biomechanics during the forward-bending of the upper body (flexion) are well investigated by both in vivo and in vitro experiments. In both cases, the experimentally observed relative motion of vertebral bodies can be used to calculate the instantaneous center of rotation (ICR). The timely evolution of the ICR, the centrode, is widely utilized for validating computer models and is thought to serve as a criterion for distinguishing healthy and degenerative motion patterns. While in vivo motion can be induced by physiological active structures (muscles), in vitro spinal segments have to be driven by external torque-applying equipment such as spine testers. It is implicitly assumed that muscle-driven and torque-driven centrodes are similar. Here, however, we show that centrodes qualitatively depend on the impetus. Distinction is achieved by introducing confidence regions (ellipses) that comprise centrodes of seven individual multi-body simulation models, performing flexion with and without preload. Muscle-driven centrodes were generally directed superior–anterior and tail-shaped, while torque-driven centrodes were located in a comparably narrow region close to the center of mass of the caudal vertebrae. We thus argue that centrodes resulting from different experimental conditions ought to be compared with caution. Finally, the applicability of our method regarding the analysis of clinical syndromes and the assessment of surgical methods is discussed.

     

Decreasing the required lumbar extensor moment induces earlier onset of flexion relaxation

Flexion relaxation (FR) is characterized by the lumbar erector spinae (LES) becoming myoelectrically silent near full trunk flexion. This study was designed to: (1) determine if decreasing the lumbar moment during flexion would induce FR to occur earlier; (2) characterize thoracic and abdominal muscle activity during FR. Ten male participants performed four trunk flexion/extension movement conditions; lumbar moment was altered by attaching 0, 5, 10, or 15 lb counterweights to the torso. Electromyography (EMG) was recorded from eight trunk muscles. Lumbar moment, lumbar flexion and trunk inclination angles were calculated at the critical point of LES inactivation (CPLES). Results demonstrated that counterweights decreased the lumbar moment and lumbar flexion angle at CPLES (p < 0.0001 and p = 0.0029, respectively); the hypothesis that FR occurs earlier when lumbar moment is reduced was accepted. The counterweights did not alter trunk inclination at CPLES (p = 0.1987); this is believed to result from an altered hip to spine flexion ratio when counterweights were attached. Lumbar multifidus demonstrated FR, similar to LES, while thoracic muscles remained active throughout flexion. Abdominal muscles activated at the same instant as CPLES, except in the 15 lb condition where abdominal muscles activated before CPLES resulting in a period of increased co-contraction.     

Validation and discovery from computational biology models

Simulation software is often a fundamental component in systems biology projects and provides a key aspect of the integration of experimental and analytical techniques in the search for greater understanding and prediction of biology at the systems level. It is important that the modelling and analysis software is reliable and that techniques exist for automating the analysis of the vast amounts of data which such simulation environments generate. A rigorous approach to the development of complex modelling software is needed. Such a framework is presented here together with techniques for the automated analysis of such models and a process for the automatic discovery of biological phenomena from large simulation data sets. Illustrations are taken from a major systems biology research project involving the in vitro investigation, modelling and simulation of epithelial tissue.     

Neuromuscular manifestations of viscoelastic tissue degradation following high and low risk repetitive lumbar flexion

Cumulative lumbar disorder is common in individuals engaged in long term performance of repetitive and static occupational/sports activities with the spine. The triggering source and of the disorder, the tissues involved in the failure and the biomechanical, neuromuscular, and biological processes active in the initiation and development of the disorder are not known. The hypothesis is forwarded that static and repetitive (cyclic) lumbar flexion-extension and the associated repeated stretch of the various viscoelastic tissues (ligaments, fascia, facet capsule, discs, etc.) causes micro-damage in their collagen fibers followed by an acute inflammation, triggering pain and reflexive muscle spasms/hyper-excitability. Continued exposure to activities, over time, converts the acute inflammation into a chronic one, viscoelastic tissues remodeling/degeneration, modified motor control strategy and permanent disability. Changes in lumbar stability are expected during the development of the disorder. A series of experimental data from in-vivo feline is reviewed and integrated with supporting evidence from the literature to gain a valuable insight into the multi-factorial development of the disorder. Prolonged cyclic lumbar flexion-extension at high loads, high velocities, many repetitions and short in between rest periods induced transient creep/laxity in the spine, muscle spasms and reduced stability followed, several hours later, by an acute inflammation/tissue degradation, muscular hyper-excitability and increased stability. The major findings assert that viscoelastic tissues sub-failure damage is the source and inflammation is the process which governs the mechanical and neuromuscular characteristic symptoms of the disorder. A comprehensive model of the disorder is presented. The experimental data validates the hypothesis as well as provide insights into the development of potential treatment and prevention of the disorder.     

Validation of indicators used to assess unconsciousness in veal calves at slaughter

European legislation states that after stunning regular checks should be performed to guarantee animals are unconscious between the end of the stunning process and death. When animals are killed without prior stunning these checks should be performed before the animal is released from restraint. The validity of certain indicators used to assess unconsciousness under different stunning and slaughter conditions is under debate. The aim of this study was to validate the absence of threat-, withdrawal-, corneal- and eyelid reflex as indicators to assess unconsciousness in calves subjected to different stunning and slaughter methods. Calves (201±22 kg) were randomly assigned to one of the following four treatments: (1) Captive bolt stunning followed by neck cut in an inverted position (n=25); (2) Non-stunned slaughter in an upright position (n=7); (3) Non-stunned slaughter in an inverted position (180° rotation) (n=25); (4) Non-stunned slaughter in an upright position followed by captive bolt stunning 40 s after the neck cut (n=25). Each calf was equipped with non-invasive electroencephalogram (EEG) electrodes before the slaughter procedure. All reflexes were verified once before the slaughter procedure. At the beginning of the procedure (T=0 s) calves were stunned (treatment 1) or neck cut in an upright position (treatment 2, 4) or inverted position (treatment 3). Calves of treatment 4 were captive bolt stunned 34±8 s after the neck cut. Reflexes were assessed every 20 s from T=15 s for all treatments until all reflex tests resulted in a negative response three times in a row and a flat line EEG was observed. In addition, reflexes were assessed 5 s after captive bolt stunning in calves of treatments 1 and 4. Visual assessment of changes in the amplitude and frequency of EEG traces was used to determine loss of consciousness. Timing of loss of consciousness was related to timing of loss of reflexes. After captive bolt stunning, absence of threat-, withdrawal-, corneal- and eyelid reflex indicated unconsciousness as determined by EEG recordings. After non-stunned slaughter, both threat- and withdrawal reflex were on average lost before calves were unconscious based on EEG recordings. The eyelid- and corneal reflex were on average lost after calves had lost consciousness based on EEG recordings and appeared to be distinctly conservative indicators of unconsciousness in non-stunned slaughtered calves since they were observed until 76±50 and 85±45 s (mean±SD), respectively, after EEG-based loss of consciousness.     

Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine

Low back mechanics are important to quantify to study injury, pain and disability. As in vivo forces are difficult to measure directly, modeling approaches are commonly used to estimate these forces. Validation of model estimates is critical to gain confidence in modeling results across populations of interest, such as people with lower-limb amputation. Motion capture, ground reaction force and electromyographic data were collected from ten participants without an amputation (five male/five female) and five participants with a unilateral transtibial amputation (four male/one female) during trunk-pelvis range of motion trials in flexion/extension, lateral bending and axial rotation. A musculoskeletal model with a detailed lumbar spine and the legs including 294 muscles was used to predict L4-L5 loading and muscle activations using static optimization. Model estimates of L4-L5 intervertebral joint loading were compared to measured intradiscal pressures from the literature and muscle activations were compared to electromyographic signals. Model loading estimates were only significantly different from experimental measurements during trunk extension for males without an amputation and for people with an amputation, which may suggest a greater portion of L4-L5 axial load transfer through the facet joints, as facet loads are not captured by intradiscal pressure transducers. Pressure estimates between the model and previous work were not significantly different for flexion, lateral bending or axial rotation. Timing of model-estimated muscle activations compared well with electromyographic activity of the lumbar paraspinals and upper erector spinae. Validated estimates of low back loading can increase the applicability of musculoskeletal models to clinical diagnosis and treatment.     

Development of a two-dimension manifold to represent high dimension mathematical models of the intracellular Mammalian circadian clock

A new focus for mathematical models of the circadian pacemaker involves the encapsulation within the models of detailed biological processes responsible for generating those circadian rhythms. Representing greater biological detail requires more mathematical equations, which pose a greater challenge for the analysis of such systems. Development of a method that retains the predominant dynamics while still providing biologically detailed information is advantageous. Two high-dimension mathematical models of intracellular mammalian circadian pacemakers, Leloup-Goldbeter and Forger-Peskin, with 19 and 73 differential equations, respectively, have been published. The authors projected each of these high-dimension models onto their respective manifold using proper orthogonal functions (POFs) obtained from the empirical decomposition of the model's phase space to obtain a 2-dimension model. The resulting 2-dimension model, represented by 2 differential equations, predicts most of the salient characteristics of a biological clock including approximately 24-h oscillations, entrainment to an LD cycle, phase response curves, and the amplitude recovery dynamics that emerge following amplitude suppression. The manifold representation simplifies the mathematical analysis, since only 2 variables need to be observed and analyzed to understand the behavior of the biological clock. This reduced model derived from a model based on biological variables can be used for the development and analysis of mathematical models of the coupled mammalian oscillators to understand the dynamics of the integrated circadian pacemaker.     

Biomechanical implications of lumbar spinal ligament transection

Many lumbar spine surgeries either intentionally or inadvertently damage or transect spinal ligaments. The purpose of this work was to quantify the previously unknown biomechanical consequences of isolated spinal ligament transection on the remaining spinal ligaments (stress transfer), vertebrae (bone remodelling stimulus) and intervertebral discs (disc pressure) of the lumbar spine. A finite element model of the full lumbar spine was developed and validated against experimental data and tested in the primary modes of spinal motion in the intact condition. Once a ligament was removed, stress increased in the remaining spinal ligaments and changes occurred in vertebral strain energy, but disc pressure remained similar. All major biomechanical changes occurred at the same spinal level as the transected ligament, with minor changes at adjacent levels. This work demonstrates that iatrogenic damage to spinal ligaments disturbs the load sharing within the spinal ligament network and may induce significant clinically relevant changes in the spinal motion segment.