Cervical Spine Disc Deformation During In Vivo Three-Dimensional Head Movements

Impact of extracellular matrix and collagen network properties on the cervical intervertebral disc response to physiological loads: A parametric study

Current intervertebral disc finite element models are hard to validate since they describe multi-physical phenomena and contain a huge number of material properties. This work aims to simplify numerical validation/identification studies by prioritizing the sensitivity of intervertebral disc behavior to mechanical properties. A 3D fiber-reinforced hyperelastic model of a C6-C7 intervertebral disc is used to carry out the parametric study. 10 parameters describing the extracellular matrix and the collagen network behaviors are included in the parametric study. The influence of varying these parameters on the disc response is estimated during physiological movements of the head, including compression, lateral bending, flexion, and axial rotation. The obtained results highlight the high sensitivity of the disc behavior to the stiffness of the annulus fibrosus extracellular matrix for all the studied loads with a relative increase in the disc apparent stiffness by 67% for compression and by 57% for axial rotation when the annulus stiffness increases from 0.4 to 2 MPa. It is also shown that varying collagen network orientation, stiffness, and stiffening in the studied configuration range have a noticeable effect on rotational motions with a relative apparent stiffness difference reaching 6.8%, 10%, and 22%, respectively, in lateral bending. However, the collagen orientation does not affect disc response to axial load.

In vivo intervertebral disc deformation: intratissue strain patterns within adjacent discs during flexion-extension

The biomechanical function of the intervertebral disc (IVD) is a critical indicator of tissue health and pathology. The mechanical responses (displacements, strain) of the IVD to physiologic movement can be spatially complex and depend on tissue architecture, consisting of distinct compositional regions and integrity; however, IVD biomechanics are predominately uncharacterized in vivo. Here, we measured voxel-level displacement and strain patterns in adjacent IVDs in vivo by coupling magnetic resonance imaging (MRI) with cyclic motion of the cervical spine. Across adjacent disc segments, cervical flexion-extension of 10° resulted in first principal and maximum shear strains approaching 10%. Intratissue spatial analysis of the cervical IVDs, not possible with conventional techniques, revealed elevated maximum shear strains located in the posterior disc (nucleus pulposus) regions. IVD structure, based on relaxometric patterns of T₂ and T_1ρ images, did not correlate spatially with functional metrics of strain. Our approach enables a comprehensive IVD biomechanical analysis of voxel-level, intratissue strain patterns in adjacent discs in vivo, which are largely independent of MRI relaxometry. The spatial mapping of IVD biomechanics in vivo provides a functional assessment of adjacent IVDs in subjects, and provides foundational biomarkers for elastography, differentiation of disease state, and evaluation of treatment efficacy.

In vivo intervertebral kinematics and disc deformations of the human cervical spine during walking

The kinematics of the cervical spine during various functional neck motions has been widely reported. However, no data has been reported on the cervical intervertebral kinematics during walking, the most frequently performed daily functional activity. In this study, we evaluated cervical kinematics and disc deformation of asymptomatic subjects during a gait cycle using a dual fluoroscopic imaging system. Our measurements showed that the vertical translation of the cervical spine (1.6 ± 0.1 Hz) occurred at twice the frequency of the gait cycle (0.8 ± 0.1 Hz). The overall ranges of motion (ROMs) of the entire (C2-T1) cervical spine were 5.0 ± 3.1° in the flexion-extension rotation, 3.4 ± 1.0° in the lateral-bending rotation, and 5.8 ± 2.1° in the axial-twisting rotation during walking. Each intervertebral disc (measured at the disc centre location) dynamically deformed in its axial direction in a range of 16.2 ± 5.7% ~ 23.7 ± 8.7% (without significant differences among different segment levels, p > 0.05), similar to the ranges of shear deformations of the same disc (p > 0.05, except for the C7-T1 disc, where p = 0.010). These data could be useful for improvements of diagnosis and treatment methods of cervical pathologies.

Dynamic functional nucleus is a potential biomarker for structural degeneration in cervical spine discs

If intervertebral disc degeneration can be identified early, preventative treatments may be initiated before symptoms become disabling and costly. Changes in disc mechanics, such as the decrease in the compressive modulus of the nucleus, are some of the earliest signs of degeneration. Therefore, in vivo changes in the disc response to compressive load may serve as a biomarker for pending or early disc degeneration. The aim of this study was to assess the potential for using in vivo dynamic disc deformation to identify pathologic structural degeneration of the intervertebral disc. A validated model-based tracking technique determined vertebral motion from biplane radiographs collected during dynamic flexion/extension and axial rotation of the cervical spine. A computational model of the subaxial intervertebral discs was developed to identify the dynamic functional nucleus of each disc, that is, the disc region that underwent little to no additional compression during dynamic movements. The size and location of the dynamic functional nucleus was determined for 10 C5/C6 spondylosis patients, 10 C5/C6/C7 spondylosis patients, and 10 asymptomatic controls. The dynamic functional nucleus size was sensitive (significantly smaller than controls in 5 of 6 measurements at the diseased disc) and specific (no difference from controls in 9 of 10 measurements at non-diseased discs) to pathologic disc degeneration. These results provide evidence to suggest that structural disc degeneration, manifested by changes in the disc response to functional loading, may be identified in vivo from dynamic imaging collected during functional movements. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 9999:1-7, 2019.

Influence of morphological variations on cervical spine segmental responses from inertial loading

OBJECTIVES

The objective of this study was to investigate the influence of morphological variations in osteoligamentous lower cervical spinal segment responses under postero-anterior inertial loading.

METHODS

A parametric finite element model of the C5-C6 spinal segment was used to generate models. Variations in the vertebral body and facet depth (anteroposterior), posterior process length, intervertebral disc height, facet articular process height and slope, segment orientation ranging from lordotic to straight, and segment size were parameterized. These variations included male-female differences. A Latin hypercube sampling method was used to select parameter values for model generation. Forces and moments associated with the inertial loading were applied to the generated model segments. The 7 parameters were grouped as local or global depending on the number of spinal components involved in the shape variation. Four output responses representing overall segmental and soft tissue responses were analyzed for each model variation: response angle of the segment, anterior longitudinal ligament stretch, anterior capsular ligament stretch, and facet joint compression in the posterior region. Pearson's correlation coefficient was used to compute the correlations of these output responses with morphological variations.

RESULTS

Fifty models were generated from the parameterized model using a Latin hypercube sampling technique. Variation in response angle among the models was 4° and was most influenced by change in the combined dimension of vertebral body and facet depth, followed by size of the segment. The maximum anterior longitudinal ligament stretch varied between 0.1 and 0.3 and was strongly influenced by the change in the segment orientation. The anterior facet joint region sustained tension, whereas the posterior region sustained compression. For the anterior capsular ligament stretch, the most influential global variation was segment orientation, whereas the most influential local variations were the facet height and facet angle parameters. In the case of posterior facet joint compression, segment orientation was again most influential, whereas among the local variations, the facet angle had the most influence.

CONCLUSION

Shape variations in the intervertebral disc influenced segmental rotation and ligament responses; however, the influence of shape variations in the facet joint was confined to capsular ligament responses. Response angle was most influenced by the vertebral body depth variations, explaining greater segmental rotations in female spines. Straighter spine segments sustained greater posterior facet joint compression, which may offer an explanation for the higher incidence of whiplash-associated disorders among females, who exhibit a straighter cervical spine. The anterior longitudinal ligament stretch was also greater in straighter segments. These findings indicate that the morphological features specific to the anatomy of the female cervical spine may predispose it to injury under inertial loading.

Ranges of Cervical Intervertebral Disc Deformation During an In Vivo Dynamic Flexion-Extension of the Neck

While abnormal loading is widely believed to cause cervical spine disc diseases, in vivo cervical disc deformation during dynamic neck motion has not been well delineated. This study investigated the range of cervical disc deformation during an in vivo functional flexion-extension of the neck. Ten asymptomatic human subjects were tested using a combined dual fluoroscopic imaging system (DFIS) and magnetic resonance imaging (MRI)-based three-dimensional (3D) modeling technique. Overall disc deformation was determined using the changes of the space geometry between upper and lower endplates of each intervertebral segment (C3/4, C4/5, C5/6, and C6/7). Five points (anterior, center, posterior, left, and right) of each disc were analyzed to examine the disc deformation distributions. The data indicated that between the functional maximum flexion and extension of the neck, the anterior points of the discs experienced large changes of distraction/compression deformation and shear deformation. The higher level discs experienced higher ranges of disc deformation. No significant difference was found in deformation ranges at posterior points of all the discs. The data indicated that the range of disc deformation is disc level dependent and the anterior region experienced larger changes of deformation than the center and posterior regions, except for the C6/7 disc. The data obtained from this study could serve as baseline knowledge for the understanding of the cervical spine disc biomechanics and for investigation of the biomechanical etiology of disc diseases. These data could also provide insights for development of motion preservation surgeries for cervical spine.