Surgical instrumentation for the in vivo determination of lumbar spinal segment stiffness and viscoelasticity

The definition of spinal instability is still controversial. For this reason, it is essential to better understand the difference in biomechanical behaviour between healthy and degenerated human spinal segments in vivo. A novel computer-assisted instrument was developed with the objective to characterize the biomechanical parameters of the spinal segment. Investigation of the viscoelastic properties as well as the dynamic spinal stiffness was performed during a minimally invasive procedure (microdiscectomy) on five patients. Measurements were performed intraoperatively and the protocol consisted of a dynamic part, where spinal stiffness was computed, and a static part, where force relaxation of the segment under constant elongation was studied. The repeatability of the measurement procedure was demonstrated with five replicated tests. The spinal segment tissues were found to have viscoelastic properties. Preliminary tests confirmed a decrease in stiffness after decompression surgery. Patients with non-relaxed muscles showed higher stiffness and relaxation rate compared to patients with relaxed muscles, which can be explained by the contraction and relaxation reflex of muscles under fast and then static elongation. The results show the usefulness of the biomechanical characterization of the human lumbar spinal segment to improve the understanding of the contribution of individual anatomical structures to spinal stability.

Early stage disc degeneration does not have an appreciable affect on stiffness and load transfer following vertebroplasty and kyphoplasty.

Vertebroplasty and kyphoplasty have been reported to alter the mechanical behavior of the treated and adjacent-level segments, and have been suggested to increase the risk for adjacent-level fractures. The intervertebral disc (IVD) plays an important role in the mechanical behavior of vertebral motion segments. Comparisons between normal and degenerative IVD motion segments following cement augmentation have yet to be reported. A microstructural finite element model of a degenerative IVD motion segment was constructed from micro-CT images. Microdamage within the vertebral body trabecular structure was used to simulate a slightly (I = 83.5% of intact stiffness), moderately (II = 57.8% of intact stiffness), and severely (III = 16.0% of intact stiffness) damaged motion segment. Six variable geometry single-segment cement repair strategies (models A-F) were studied at each damage level (I-III). IVD and bone stresses, and motion segment stiffness, were compared with the intact and baseline damage models (untreated), as well as, previous findings using normal IVD models with the same repair strategies. Overall, small differences were observed in motion segment stiffness and average stresses between the degenerative and normal disc repair models. We did however observe a reduction in endplate bulge and a redistribution in the microstructural tissue level stresses across both endplates and in the treated segment following early stage IVD degeneration. The cement augmentation strategy placing bone cement along the periphery of the vertebra (model E) proved to be the most advantageous in treating the degenerative IVD models by showing larger reductions in the average bone stresses (vertebral and endplate) as compared to the normal IVD models. Furthermore, only this repair strategy, and the complete cement fill strategy (model F), were able to restore the slightly damaged (I) motion segment stiffness above pre-damaged (intact) levels. Early stage IVD degeneration does not have an appreciable effect in motion segment stiffness and average stresses in the treated and adjacent-level segments following vertebroplasty and kyphoplasty. Placing bone cement in the periphery of the damaged vertebra in a degenerative IVD motion segment, minimizes load transfer, and may reduce the likelihood of adjacent-level fractures.

Effects of tensioning the lumbar fasciae on segmental stiffness during flexion and extension - Young investigator award winner

Study Design. Biomechanical study of unembalmed human lumbar segments. Objective. To investigate the effects of tensioning the lumbar fasciae ( transversus abdominis [TrA]) aponeurosis) on segment stiffness during flexion and extension. Summary of Background Data. Animal and human studies suggest that TrA may influence intersegmental movement via tension in the middle and posterior layers of lumbar fasciae ( MLF, PLF). Methods. Compressive flexion and extension moments were applied to 17 lumbar segments from 9 unembalmed cadavers with 20 N lateral tension of the TrA aponeurosis during: 1) static tests: load was compared when fascial tension was applied during static compressive loads into flexion-extension; 2) cyclic loading tests: load, axial displacement, and stiffness were compared during repeated compressive loading cycles into flexion-extension. After testing, the PLF was incised to determine the tension transmitted by each layer. Results. At all segments and loads (< 200 N), fascial tension increased resistance to flexion loads by similar to 9.5 N. In 15 of 17, fascial tension decreased resistance to extension by similar to 6.6 N. Fascial tension during cyclic flexion loading decreased axial displacement by 26% at the onset of loading (0 - 2 N) and 2% at 450 N ( 13 of 17). During extension loading, fascial tension increased displacement at the onset of loading ( 10 of 17) by similar to 23% and slightly (1%) decreased displacement at 450 N. Segment stiffness was increased by 6 N/mm in flexion (44% at 25 N) and decreased by 2 N/mm (8% at 25 N) in extension. More than 85% of tension was transmitted through the MLF. Conclusions. Tension on the lumbar fasciae simulating moderate contraction of TrA affects segmental stiffness, particularly toward the neutral zone.

Numerical simulations of the spine for patient-specific models of severe scoliosis

In silico methods, such as musculoskeletal modelling, may aid the selection of the optimal surgical treatment for highly complex pathologies such as scoliosis. Many musculoskeletal models use a generic, simplified representation of the intervertebral joints, which are fundamental to the flexibility of the spine. Therefore, to model and simulate the spine, a suitable representation of the intervertebral joint is crucial. The aim of this PhD was to characterise specimen-specific models of the intervertebral joint for multi-body models from experimental datasets. First, the project investigated the characterisation of a specimen-specific lumped parameter model of the intervertebral joint from an experimental dataset of a four-vertebra lumbar spine segment. Specimen-specific stiffnesses were determined with an optimisation method. The sensitivity of the parameters to the joint pose was investigate. Results showed the stiffnesses and predicted motions were highly depended on both the joint pose. Following the first study, the method was reapplied to another dataset that included six complete lumbar spine segments under three different loading conditions. Specimen-specific uniform stiffnesses across joint levels and level-dependent stiffnesses were calculated by optimisation. Specimen-specific stiffness show high inter-specimen variability and were also specific to the loading condition. Level-dependent stiffnesses are necessary for accurate kinematic predictions and should be determined independently of one another. Secondly, a framework to create subject-specific musculoskeletal models of individuals with severe scoliosis was developed. This resulted in a robust codified pipeline for creating subject-specific, severely scoliotic spine models from CT data. In conclusion, this thesis showed that specimen-specific intervertebral joint stiffnesses were highly sensitive to joint pose definition and the importance of level-dependent optimisation. Further, an open-source codified pipeline to create patient-specific scoliotic spine models from CT data was released. These studies and this pipeline can facilitate the specimen-specific characterisation of the scoliotic intervertebral joint and its incorporation into scoliotic musculoskeletal spine models.

The role of B-mode ultrasonography in the musculoskeletal anatomical evaluation of the cervical region of the dog spine

This study characterized the normal musculoskeletal anatomy of the cervical segment of the spine of dogs by means of B-mode ultrasonography. The objective was to establish the role of B-mode ultrasonography for the anatomical evaluation of the cervical spine segment in dogs, by comparing the ultrasonographic findings with images by computed tomography and magnetic resonance imaging. The ultrasound examination, in transverse and median sagittal sections, allowed to identify a part of the epaxial cervical musculature, the bone surface of the cervical vertebrae and parts of the spinal cord through restricted areas with natural acoustic windows, such as between the atlanto-occipital joint, axis and atlas, and axis and the third cervical vertebra. The images, on transverse and sagittal planes, by low-field magnetic resonance imaging, were superior for the anatomical identification of the structures, due to higher contrast between the different tissues in this modality. Computed tomography showed superiority for bone detailing when compared with ultrasonography. As for magnetic resonance imaging, in addition to the muscles and cervical vertebrae, it is possible to identify the cerebrospinal fluid and differentiate between the nucleus pulposus and annulus fibrosus of the intervertebral discs. Although not the scope of this study, with knowledge of the ultrasonographic anatomy of this region, it is believed that some lesions can be identified, yet in a limited manner, when compared with the information obtained mainly with magnetic resonance imaging. The ultrasound examination presented lower morphology diagnostic value compared with the other modalities.

The role of B-mode ultrasonography in the musculoskeletal anatomical evaluation of the cervical region of the dog spine

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Patient-specific spinal stiffness in AIS: a preoperative and noninvasive method.

INTRODUCTION The clinical tests currently used to assess spinal biomechanics preoperatively are unable to assess true mechanical spinal stiffness. They rely on spinal displacement without considering the force required to deform a patient's spine. We propose a preoperative method for noninvasively quantifying the three-dimensional patient-specific stiffness of the spines of adolescent idiopathic scoliosis patients. METHODS The technique combines a novel clinical test with numerical optimization of a finite element model of the patient's spine. RESULTS A pilot study conducted on five patients showed that the model was able to provide accurate 3D reconstruction of the spine's midline and predict the spine's stiffness for each patient in flexion, bending, and rotation. Statistically significant variation of spinal stiffness was observed between the patients. CONCLUSION This result confirms that spinal biomechanics is patient-specific, which should be taken into consideration to individualize surgical treatment.

2D-3D Rigid-Body Registration of X-Ray Fluoroscopy and CT Images

The registration of pre-operative volumetric datasets to intra- operative two-dimensional images provides an improved way of verifying patient position and medical instrument loca- tion. In applications from orthopedics to neurosurgery, it has a great value in maintaining up-to-date information about changes due to intervention. We propose a mutual information- based registration algorithm to establish the proper align- ment. For optimization purposes, we compare the perfor- mance of the non-gradient Powell method and two slightly di erent versions of a stochastic gradient ascent strategy: one using a sparsely sampled histogramming approach and the other Parzen windowing to carry out probability density approximation. Our main contribution lies in adopting the stochastic ap- proximation scheme successfully applied in 3D-3D registra- tion problems to the 2D-3D scenario, which obviates the need for the generation of full DRRs at each iteration of pose op- timization. This facilitates a considerable savings in compu- tation expense. We also introduce a new probability density estimator for image intensities via sparse histogramming, de- rive gradient estimates for the density measures required by the maximization procedure and introduce the framework for a multiresolution strategy to the problem. Registration results are presented on uoroscopy and CT datasets of a plastic pelvis and a real skull, and on a high-resolution CT- derived simulated dataset of a real skull, a plastic skull, a plastic pelvis and a plastic lumbar spine segment.

Prevalencia de sintomatología musculoesquelética en trabajadores de una empresa avícola de Cundinamarca en el año 2013

RESUMEN Objetivos Determinar la prevalencia de sintomatología musculoesquelética en trabajadores en una empresa avícola de Cundinamarca en el año 2013. Metodología Se eligió un diseño observacional, descriptivo, transversal mediante un método de recolección de datos a través de una herramienta estandarizada el “Cuestionario Nórdico musculoesquelético”(1), se adicionaron datos sociodemográficos como edad, género, cargo laboral y antigüedad laboral. El análisis del presente estudio se realizó a través de herramienta informática estadística “EPIINFO versión 3.5.4”(2). Resultados Del total de la población estudiada en la empresa avícola de Cundinamarca que fue de 229 trabajadores, se excluyeron 65 trabajadores por llevar menos de un año de labor en la empresa; quedando 164 trabajadores de las que correspondieron al género femenino 36.0% (n=58) y del género masculino 64.0% (n=105). El promedio de edad fue de 38.1 ± 9.4 años con un mínimo de 18 años y un máximo de 56 años. La distribución de acuerdo a nivel educativo principalmente fueron bachilleres con el 48.8% (n=80) y primaria fueron 34.8% (n=57). La prevalencia de síntomas osteomusculares en los últimos 12 meses por segmentos corporales fue más frecuente en manos y muñecas con un 65.2%, seguido del segmento de columna lumbar con un 48.0%; por número de segmentos comprometidos se encontró una prevalencia de 4 o más segmentos del 52.2%. Conclusión La prevalencia global de sintomatología musculoesquelética de nuestro estudio aunque coincidente con la literatura, evidencia a nivel de segmentos diferencias importantes como en el de manos muñecas, donde la séptima encuesta nacional de condiciones de trabajo del Instituto Nacional de Higiene en el Trabajo de España (INHT) reporta el 10.8%(3). En nuestro estudio se encontró una mayor prevalencia de síntomas en 4 o más segmentos comprometidos, comparado con el estudio Prevalence of multisite musculoskeletal symptoms: a French cross-sectional working population-based study que reporta una prevalencia para 4 o más segmentos del 27.0% al 35.0%(12), lo que nos hace pensar en la coexistencia de múltiples factores de riesgos ergonómicos.

Caracterização da rigidez póstero-anterior da coluna lombar in vivo: uma análise biomecânica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Vibration as a method for muscular stimulation

Human reactions to vibration have been extensively investigated in the past. Vibration, as well as whole-body vibration (WBV), has been commonly considered as an occupational hazard for its detrimental effects on human condition and comfort. Although long term exposure to vibrations may produce undesirable side-effects, a great part of the literature is dedicated to the positive effects of WBV when used as method for muscular stimulation and as an exercise intervention. Whole body vibration training (WBVT) aims to mechanically activate muscles by eliciting neuromuscular activity (muscle reflexes) via the use of vibrations delivered to the whole body. The most mentioned mechanism to explain the neuromuscular outcomes of vibration is the elicited neuromuscular activation. Local tendon vibrations induce activity of the muscle spindle Ia fibers, mediated by monosynaptic and polysynaptic pathways: a reflex muscle contraction known as the Tonic Vibration Reflex (TVR) arises in response to such vibratory stimulus. In WBVT mechanical vibrations, in a range from 10 to 80 Hz and peak to peak displacements from 1 to 10 mm, are usually transmitted to the patient body by the use of oscillating platforms. Vibrations are then transferred from the platform to a specific muscle group through the subject body. To customize WBV treatments, surface electromyography (SEMG) signals are often used to reveal the best stimulation frequency for each subject. Use of SEMG concise parameters, such as root mean square values of the recordings, is also a common practice; frequently a preliminary session can take place in order to discover the more appropriate stimulation frequency. Soft tissues act as wobbling masses vibrating in a damped manner in response to mechanical excitation; Muscle Tuning hypothesis suggest that neuromuscular system works to damp the soft tissue oscillation that occurs in response to vibrations; muscles alters their activity to dampen the vibrations, preventing any resonance phenomenon. Muscle response to vibration is however a complex phenomenon as it depends on different parameters, like muscle-tension, muscle or segment-stiffness, amplitude and frequency of the mechanical vibration. Additionally, while in the TVR study the applied vibratory stimulus and the muscle conditions are completely characterised (a known vibration source is applied directly to a stretched/shortened muscle or tendon), in WBV study only the stimulus applied to a distal part of the body is known. Moreover, mechanical response changes in relation to the posture. The transmissibility of vibratory stimulus along the body segment strongly depends on the position held by the subject. The aim of this work was the investigation on the effects that the use of vibrations, in particular the effects of whole body vibrations, may have on muscular activity. A new approach to discover the more appropriate stimulus frequency, by the use of accelerometers, was also explored. Different subjects, not affected by any known neurological or musculoskeletal disorders, were voluntarily involved in the study and gave their informed, written consent to participate. The device used to deliver vibration to the subjects was a vibrating platform. Vibrations impressed by the platform were exclusively vertical; platform displacement was sinusoidal with an intensity (peak-to-peak displacement) set to 1.2 mm and with a frequency ranging from 10 to 80 Hz. All the subjects familiarized with the device and the proper positioning. Two different posture were explored in this study: position 1 - hack squat; position 2 - subject standing on toes with heels raised. SEMG signals from the Rectus Femoris (RF), Vastus Lateralis (VL) and Vastus medialis (VM) were recorded. SEMG signals were amplified using a multi-channel, isolated biomedical signal amplifier The gain was set to 1000 V/V and a band pass filter (-3dB frequency 10 - 500 Hz) was applied; no notch filters were used to suppress line interference. Tiny and lightweight (less than 10 g) three-axial MEMS accelerometers (Freescale semiconductors) were used to measure accelerations of onto patient’s skin, at EMG electrodes level. Accelerations signals provided information related to individuals’ RF, Biceps Femoris (BF) and Gastrocnemius Lateralis (GL) muscle belly oscillation; they were pre-processed in order to exclude influence of gravity. As demonstrated by our results, vibrations generate peculiar, not negligible motion artifact on skin electrodes. Artifact amplitude is generally unpredictable; it appeared in all the quadriceps muscles analysed, but in different amounts. Artifact harmonics extend throughout the EMG spectrum, making classic high-pass filters ineffective; however, their contribution was easy to filter out from the raw EMG signal with a series of sharp notch filters centred at the vibration frequency and its superior harmonics (1.5 Hz wide). However, use of these simple filters prevents the revelation of EMG power potential variation in the mentioned filtered bands. Moreover our experience suggests that the possibility of reducing motion artefact, by using particular electrodes and by accurately preparing the subject’s skin, is not easily viable; even though some small improvements were obtained, it was not possible to substantially decrease the artifact. Anyway, getting rid of those artifacts lead to some true EMG signal loss. Nevertheless, our preliminary results suggest that the use of notch filters at vibration frequency and its harmonics is suitable for motion artifacts filtering. In RF SEMG recordings during vibratory stimulation only a little EMG power increment should be contained in the mentioned filtered bands due to synchronous electromyographic activity of the muscle. Moreover, it is better to remove the artifact that, in our experience, was found to be more than 40% of the total signal power. In summary, many variables have to be taken into account: in addition to amplitude, frequency and duration of vibration treatment, other fundamental variables were found to be subject anatomy, individual physiological condition and subject’s positioning on the platform. Studies on WBV treatments that include surface EMG analysis to asses muscular activity during vibratory stimulation should take into account the presence of motion artifacts. Appropriate filtering of artifacts, to reveal the actual effect on muscle contraction elicited by vibration stimulus, is mandatory. However as a result of our preliminary study, a simple multi-band notch filtering may help to reduce randomness of the results. Muscle tuning hypothesis seemed to be confirmed. Our results suggested that the effects of WBV are linked to the actual muscle motion (displacement). The greater was the muscle belly displacement the higher was found the muscle activity. The maximum muscle activity has been found in correspondence with the local mechanical resonance, suggesting a more effective stimulation at the specific system resonance frequency. Holding the hypothesis that muscle activation is proportional to muscle displacement, treatment optimization could be obtained by simply monitoring local acceleration (resonance). However, our study revealed some short term effects of vibratory stimulus; prolonged studies should be assembled in order to consider the long term effectiveness of these results. Since local stimulus depends on the kinematic chain involved, WBV muscle stimulation has to take into account the transmissibility of the stimulus along the body segment in order to ensure that vibratory stimulation effectively reaches the target muscle. Combination of local resonance and muscle response should also be further investigated to prevent hazards to individuals undergoing WBV treatments.

Combining 3D tracking and surgical instrumentation to determine the stiffness of spinal motion segments: a validation study

The spine is a complex structure that provides motion in three directions: flexion and extension, lateral bending and axial rotation. So far, the investigation of the mechanical and kinematic behavior of the basic unit of the spine, a motion segment, is predominantly a domain of in vitro experiments on spinal loading simulators. Most existing approaches to measure spinal stiffness intraoperatively in an in vivo environment use a distractor. However, these concepts usually assume a planar loading and motion. The objective of our study was to develop and validate an apparatus, that allows to perform intraoperative in vivo measurements to determine both the applied force and the resulting motion in three dimensional space. The proposed setup combines force measurement with an instrumented distractor and motion tracking with an optoelectronic system. As the orientation of the applied force and the three dimensional motion is known, not only force-displacement, but also moment-angle relations could be determined. The validation was performed using three cadaveric lumbar ovine spines. The lateral bending stiffness of two motion segments per specimen was determined with the proposed concept and compared with the stiffness acquired on a spinal loading simulator which was considered to be gold standard. The mean values of the stiffness computed with the proposed concept were within a range of ±15% compared to data obtained with the spinal loading simulator under applied loads of less than 5 Nm.