Validation of a noninvasive technique to precisely measure in vivo three-dimensional cervical spine movement

In vivo kinematic study of lumbar center of rotation under different loads

BACKGROUND

Dual fluoroscopic imaging system (DFIS) was employed to identify the Center of Rotation(COR) in the lower lumbar spine and determine its relationship with weight bearing.

METHODS

In this study, twenty participants were recruited. A 3D model of each participant's lumbar spine was created using CT images, and their relative positions were determined through DFIS. By integrating CT imaging with DFIS, the kinematic data of the participants' spines during movement were captured. The lower lumbar spine's COR was calculated using the method of perpendicular bisectors.

RESULTS

While flexing and extending, the Center of Rotation (COR) initially moved downward with increasing load, followed by upward movement as the load further increased. During flexion and extension, the COR coordinates of L3-4 at 0 kg, 5 kg and 10 kg are(0.3549 ± 0.2176,0.0177 ± 0.1317),(0.0598 ± 0.2095,-0.1806 ± 0.1719),(0.1427 ± 0.1440,-0.0911 ± 0.2722); The center of rotation coordinates of L4-5 at 0 kg, 5 kg and 10 kg are(0.0566 ± 0.2693,-0.0727 ± 0.2132),(0.0964 ± 0.2671,-0.2037 ± 0.2299),(0.1648 ± 0.1520,-0.0049 ± 0.1641). The anterior-posterior position of the COR shifted posteriorly with increasing weight-bearing. During lateral bending, the center of rotation coordinates of L3-4 at 0 kg, 5 kg and 10 kg are(0.0745 ± 0.1229,0.0966 ± 0.3403) (-0.0438 ± 0.1281,0.1161 ± 0.1584), (-0.0464 ± 0.1517,0.1320 ± 0.2730); The center of rotation coordinates of L4-5 at 0 kg, 5 kg and 10 kg are(-0.0314 ± 0.1411,-0.0355 ± 0.2088), (-0.0764 ± 0.3135,0.0105 ± 0.3230),(-0.0376 ± 0.1701,0.0285 ± 0.2395). Throughout the lateral bending exercises, the upper and lower COR positions increased as the load increased, while the left and right COR positions remained unaffected by the load increment. The COR height differed between flexion and lateral bending. We observed variations in the COR position of the lumbar spine during lateral bending and flexion-extension movements. This enhanced our comprehension of coupled motion patterns within the lumbar spine.

CONCLUSIONS

Position of the lumbar spine COR changes with variations in the load. During different movements, the COR location of the lower lumbar spine varied. This finding suggests the presence of distinct motion patterns in the lower lumbar spine. As the load increases, the lumbar COR position changes significantly. Abnormal movement patterns of the lower lumbar spine under different loads may be one of the factors that accelerate lumbar disc degeneration.

Inverse kinematics in cervical spine models: Effects of scaling and model degrees of freedom for extension and flexion movements

Intervertebral kinematics can affect model-predicted loads and strains in the spine; therefore knowledge of expected vertebral kinematics error is important for understanding the limitations of model predictions. This study addressed how different kinematic models of the neck affect the prediction of intervertebral kinematics from markers on the head and trunk. Eight subjects executed head and neck extension-flexion motion with simultaneous motion capture and biplanar dynamic stereo-radiography (DSX) of vertebrae C1-C7. A generic head and neck model in OpenSim was scaled by marker data, and three versions of the model were used with an inverse kinematics solver. The models differed according to the number of independent degrees of freedom (DOF) between the head and trunk: 3 rotational DOF with constraints defining intervertebral kinematics as a function of overall head-trunk motion; 24DOF with 3 independent rotational DOF at each level, skull-T1; 48DOF with 3 rotational and 3 translational DOF at each level. Marker tracking error was lower for scaled models compared to generic models and decreased as model DOF increased. The largest mean absolute error (MAE) was found in extension-flexion angle and anterior-posterior translation at C1-C2, and superior-inferior translation at C2-C3. Model scaling and complexity did not have a statistically significant effect on most error metrics when corrected for multiple comparisons, but ranges of motion were significantly different from DSX in some cases. This study evaluated model kinematics in comparison to gold standard radiographic data and provides important information about intervertebral kinematics error that are foundational to model validity.

The deflection of fatigued neck

The human neck is a unique mechanical structure, highly flexible but fatigue prone. The rising prevalence of neck pain and chronic injuries has been attributed to increasing exposure to fatigue loading in activities such as prolonged sedentary work and overuse of electronic devices. However, a causal relationship between fatigue and musculoskeletal mechanical changes remains elusive. This work aimed to establish this relationship through a unique experiment design, inspired by a cantilever beam mechanical model of the neck, and an orchestrated deployment of advanced motion-force measurement technologies including dynamic stereo-radiographic imaging. As a group of 24 subjects performed sustained-till-exhaustion neck exertions in varied positions-neutral, extended, and flexed, their cervical spine musculoskeletal responses were measured. Data verified the occurrence of fatigue and revealed fatigue-induced neck deflection which increased cervical lordosis or kyphosis by 4-5° to 11°, depending on the neck position. This finding and its interpretations render a renewed understanding of muscle fatigue from a more unified motor control perspective as well as profound implications on neck pain and injury prevention.

Dynamic in vivo 3D atlantooccipital kinematics during multiplanar physiologic motions

Normal biomechanics of the upper cervical spine, particularly at the atlantooccipital joint, remain poorly characterized. The purpose of this study was to determine the intervertebral kinematics of the atlantooccipital joint (occiput-C1) during three-dimensional in vivo physiologic movements. Twenty healthy young adults performed dynamic flexion/extension, axial rotation, and lateral bending while biplane radiographs were collected at 30 images per second. Motion at occiput-C1 was tracked using a validated volumetric model-based tracking process that matched subject-specific CT-based bone models to the radiographs. The occiput-C1 total range of motion (ROM) and helical axis of motion (HAM) was calculated for each movement. During flexion/extension, the occiput-C1 moved almost exclusively in-plane (ROM: 17.9 ± 6.9°) with high variability in kinematic waveforms (6.3°) compared to the in-plane variability during axial rotation (1.4°) and lateral bending (0.9°) movements. During axial rotation, there was small in-plane motion (ROM: 4.2 ± 2.5°) compared to out-of-plane flexion/extension (ROM: 12.7 ± 5.4°). During lateral bending, motion occurred in-plane (ROM: 9.0 ± 3.1°) and in the plane of flexion/extension (ROM: 7.3 ± 2.7°). The average occiput-C1 axis of rotation intersected the sagittal and coronal planes 7 mm to 18 mm superior to the occipital condyles. The occiput-C1 axis of rotation pointed 60° from the sagittal plane during axial rotation but only 10° from the sagittal plane during head lateral bending. These novel results are foundational for future work on upper cervical spine kinematics.

In vivo cervical vertebrae kinematic studies based on dual fluoroscopic imaging system measurement: A narrative review

Understanding the motion characteristics of cervical spine through biomechanical analysis aids in the identification of abnormal joint movements. This knowledge is essential for the prevention, diagnosis, and treatment of related disorders. However, the anatomical structure of the cervical spine is complex, and traditional medical imaging techniques have certain limitations. Capturing the movement characteristics of various parts of the cervical spine in vivo during motion is challenging. The dual fluoroscopic imaging system (DFIS) is able to quantify the motion and motion patterns of individual segments. In recent years, DFIS has achieved accurate non-invasive measurements of dynamic joint movements in humans. This review assesses the research findings of DFIS about the cervical spine in healthy and pathological individuals. Relevant study search was conducted up to October 2023 in Web of Science, PubMed, and EBSCO databases. After the search, a total of 30 studies were ultimately included. Among them, 13 studies focused on healthy cervical spines, while 17 studies focused on pathological cervical spines. These studies mainly centered on exploring the vertebral bodies and associated structures of the cervical spine, including intervertebral discs, intervertebral foramina, and zygapophyseal joints. Further research could utilize DFIS to investigate cervical spine motion in different populations and under pathological conditions.

A scoping review of human skeletal kinematics research using biplane radiography

Biplane radiography has emerged as the gold standard for accurately measuring in vivo skeletal kinematics during physiological loading. The purpose of this scoping review was to map the extent, range, and nature of biplane radiography research on humans from 2004 through 2022. A literature search was performed using the terms biplane radiography, dual fluoroscopy, dynamic stereo X-ray, and biplane videoradiography. All articles referenced in included publications were also assessed for inclusion. A secondary search was then performed using the names of the most frequently appearing principal investigators among included papers. A total of 379 manuscripts were identified and included. The first studies published in 2004 focused on the native knee, followed by studies of the ankle joint complex in 2006, the shoulder in 2007, and the spine in 2008. Nearly half (180, 47.5%) of all manuscripts investigated knee kinematics. The average number of publications increased from 21.6 per year from 2012 to 2017 to 34.6 per year from 2017 to 2022. The average number of participants per study was 16, with a range from 1 to 101. A total of 90.2% of studies featured cohorts of 30 or less. The most prolific research groups for each joint were: Mass General Hospital (lumbar spine and knee), Henry Ford Hospital (shoulder), the University of Utah (ankle and hip), The University of Pittsburgh (cervical spine), and Brown University (hand/wrist/elbow). Future advancements in biplane radiography research are dependent upon increased availability of these imaging systems, standardization of data collection protocols, and the development of automated approaches to expedite data processing.

The Effects of Cervical Orthoses on Head and Intervertebral Range of Motion

STUDY DESIGN

Prospective Cohort.

OBJECTIVE

Quantify and compare the effectiveness of cervical orthoses in restricting intervertebral kinematics during multiplanar motions.

SUMMARY OF BACKGROUND DATA

Previous studies evaluating the efficacy of cervical orthoses measured global head motion and did not evaluate individual cervical motion segment mobility. Prior studies focused only on the flexion/extension motion.

METHODS

Twenty adults without neck pain participated. Vertebral motion from the occiput through T1 was imaged using dynamic biplane radiography. Intervertebral motion was measured using an automated registration process with validated accuracy better than 1 degree. Participants performed independent trials of maximal flexion/extension, axial rotation, and lateral bending in a randomized order of unbraced, soft collar (foam), hard collar (Aspen), and cervical thoracic orthosis (CTO) (Aspen) conditions. Repeated-measures ANOVA was used to identify differences in the range of motion (ROM) among brace conditions for each motion.

RESULTS

Compared with no collar, the soft collar reduced flexion/extension ROM from occiput/C1 through C4/C5, and reduced axial rotation ROM at C1/C2 and from C3/C4 through C5/C6. The soft collar did not reduce motion at any motion segment during lateral bending. Compared with the soft collar, the hard collar reduced intervertebral motion at every motion segment during all motions, except for occiput/C1 during axial rotation and C1/C2 during lateral bending. The CTO reduced motion compared with the hard collar only at C6/C7 during flexion/extension and lateral bending.

CONCLUSIONS

The soft collar was ineffective as a restraint to intervertebral motion during lateral bending, but it did reduce intervertebral motion during flexion/extension and axial rotation. The hard collar reduced intervertebral motion compared with the soft collar across all motion directions. The CTO provided a minimal reduction in intervertebral motion compared with the hard collar. The utility in using a CTO rather than a hard collar is questionable, given the cost and little or no additional motion restriction.

Changes in intervertebral sagittal alignment of the cervical spine from supine to upright

Cervical sagittal alignment is a critical component of successful surgical outcomes. Unrecognized differences in intervertebral alignment between supine and upright positions may affect clinical outcomes; however, these differences have not been quantified. Sixty-four patients scheduled to undergo one or two-level cervical arthrodesis for symptomatic pathology from C4-C5 to C6-C7, and forty-seven controls were recruited. Upright sagittal alignment was obtained through biplane radiographic imaging and measured using a validated process with accuracy better than 1° in rotation. Supine alignment was obtained from computed tomography scans. Coordinate systems used to measure supine and upright alignment were identical. Distances between adjacent bony endplates were measured to calculate disc height in each position. For both patients and controls, the C1-C2, C2-C3, and C3-C4 motion segments were in more lordosis when upright as compared with supine (all p < 0.001). However, the C4-C5, C5-C6, and C6-C7 motion segments were in less lordosis when upright as compared with supine (all p ≤ 0.004). There was an interaction between group and position at the C1-C2 (p = 0.002) and C2-C3 (p = 0.001) motion segments, with the controls demonstrating a greater increase in lordosis at both motion segments when moving from supine to upright. The results indicate that cervical motion segment alignment changes between supine and upright positioning, those changes differ among motion segments, and cervical pathology affects the magnitude of these changes. Clinical Significance: Surgeons should be mindful of the differences in alignment between supine and upright imaging and the implications they may have on clinical outcomes.

Intervertebral kinematics during neck motion 6.5 years after fusion and artificial disc replacement

BACKGROUND

Arthroplasty with artificial disc replacement for surgical treatment of cervical spine degeneration was introduced with the notion that motion-preserving approaches would prevent development of adjacent segment disease. Though clinical outcomes favor arthroplasty over the commonly used anterior cervical discectomy with fusion approach, clinical studies confirming the biomechanical basis of these results are lacking. The aim of this study was to compare intervertebral kinematics between arthroplasty and fusion patients 6.5 years post-surgery during physiological motion of the neck.

METHODS

Using a biplane dynamic X-ray system, computed tomography imaging and model based tracking algorithms, three dimensional intervertebral kinematics were measured during neck axial rotation and extension in 14 patients treated for cervical radiculopathy with fusion (n = 8) or arthroplasty (n = 6). The measurements were performed at 2-year (baseline) and 6.5 year post-surgical time points, with the main interest being in the interaction between surgery types and time points. 3 translations and 3 rotations were investigated for the index (C5C6), and upper- (C4C5) and lower adjacent levels (C6C7).

FINDINGS

Surgery-time interaction was significant for axial rotation (P < 0.04) and flexion-extension rotation (P < 0.005) in C4C5 during neck axial rotation, left-right translation (P < 0.04) in C5C6 and anterior-posterior translation in C6C7 (P < 0.04) during neck extension. In contrast with the expectations, axial rotation and flexion-extension decreased in C4C5 during neck rotation and anterior-posterior translation decreased in C6C7 during neck extension for fusion.

INTERPRETATION

The findings do not support the notion that adjacent segment motion increases after fusion.

In Vivo Evidence of Early Instability and Late Stabilization in Motion Segments Immediately Superior to Anterior Cervical Arthrodesis

STUDY DESIGN

Prospective cohort study.

OBJECTIVE

The aim was to identify patient factors that affect adjacent segment kinematics after anterior cervical discectomy and fusion (ACDF) as measured by biplane radiography.

SUMMARY OF BACKGROUND DATA

The etiology of adjacent segment disease (ASD) may be multifactorial. Previous studies have investigated associations between patient factors and ASD, although few attempted to link patient factors with mechanical changes in the spine that may explain ASD development. Previous studies manually measured intervertebral motion from static flexion/extension radiographs, however, manual measurements are unreliable, and those studies failed to measure intervertebral motion during rotation.

METHODS

Patients had continuous cervical spine flexion/extension and axial rotation movements captured at 30 images per second in a dynamic biplane radiography system preoperatively and 1 year after ACDF. Digitally reconstructed radiographs generated from subject-specific computed tomography scans were matched to the biplane radiographs using a validated tracking process. Dynamic kinematics and preoperative disc height were calculated from this tracking process. Preoperative magnetic resonance imagings were evaluated for disc bulge. Patient age, sex, body mass index, smoking status, diabetes, psychiatric history, presence of an inciting event, and length of symptoms were collected. Multivariate linear regression was performed to identify patient factors associated with 1-year postoperative changes in adjacent segment kinematics.

RESULTS

Sixty-three patients completed preoperative and postoperative testing. Superior adjacent segment disc height and disc bulge predicted the change in superior adjacent segment range of motion after surgery. Inferior adjacent segment disc bulge, smoking history, and the use of psychiatric medications predicted the change in inferior adjacent segment flexion/extension range of motion after surgery.

CONCLUSIONS

Preexisting adjacent segment disc degeneration, as indicated by disc height and disc bulge, was associated with reduced adjacent segment motion after ACDF, while lack of preexisting adjacent disc degeneration was associated with increased adjacent segment motion after ACDF. These findings provide in vivo evidence supporting early instability and late stabilization in the pathophysiology of disc degeneration.

Within-subject effects of standardized prosthetic socket modifications on physical function and patient-reported outcomes

Background

Among the challenges of living with lower limb loss is the increased risk of long-term health problems that can be either attributed directly to the amputation surgery and/or prosthetic rehabilitation or indirectly to a disability-induced sedentary lifestyle. These problems are exacerbated by poorly fit prosthetic sockets. There is a knowledge gap regarding how the socket design affects in-socket mechanics and how in-socket mechanics affect patient-reported comfort and function. The objectives of this study are (1) to gain a better understanding of how in-socket mechanics of the residual limb in transfemoral amputees are related to patient-reported comfort and function, (2) to identify clinical tests that can streamline the socket design process, and (3) to evaluate the efficacy and cost of a novel, quantitatively informed socket optimization process.

Methods

Users of transfemoral prostheses will be asked to walk on a treadmill wearing their current socket plus 8 different check sockets with designed changes in different structural measurements that are likely to induce changes in residual limb motion, skin strain, and pressure distribution within the socket. Dynamic biplane radiography and pressure sensors will be used to measure in-socket residual limb mechanics. Patient-reported outcomes will also be collected after wearing each socket. The effects of in-socket mechanics on both physical function and patient-reported outcomes (aim 1) will be assessed using a generalized linear model. Partial correlation analysis will be used to examine the association between research-grade measurements and readily available clinical measurements (aim 2). In order to compare the new quantitative design method to the standard of care, patient-reported outcomes and cost will be compared between the two methods, utilizing the Wilcoxon-Mann-Whitney non-parametric test (aim 3).

Discussion

Knowledge on how prosthetic socket modifications affect residual bone and skin biomechanics itself can be applied to devise future socket designs, and the methodology can be used to investigate and improve such designs, past and present. Apart from saving time and costs, this may result in better prosthetic socket fit for a large patient population, thus increasing their mobility, participation, and overall health-related quality of life.

Trial registration

ClinicalTrials.gov NCT05041998. Date of registration: Sept 13, 2021.

Validation and application of a novel in vivo cervical spine kinematics analysis technique

To validate the accuracy of Cone beam computed tomography (CBCT) cervical spine modeling with three dimensional (3D)-3D registration for in vivo measurements of cervical spine kinematics. CBCT model accuracy was validated by superimposition with computed tomography (CT) models in 10 healthy young adults, and then cervical vertebrae were registered in six end positions of functional movements, versus a neutral position, in 5 healthy young adults. Registration errors and six degrees of freedom (6-DOF) kinematics were calculated and reported. Relative to CT models, mean deviations of the CBCT models were < 0.6 mm. Mean registration errors between end positions and the reference neutral position were < 0.7 mm. During flexion-extension (F-E), the translation in the three directions was small, mostly < 1 mm, with coupled LB and AR both < 1°. During lateral bending (LB), the bending was distributed roughly evenly, with coupled axial rotation (AR) opposite to the LB at C1-C2, and minimal coupled F-E. During AR, most of the rotation occurred in the C1-C2 segment (29.93 ± 7.19° in left twist and 31.38 ± 8.49° in right twist) and coupled LB was observed in the direction opposite to that of the AR. Model matching demonstrated submillimeter accuracy in cervical spine kinematics data. The presently evaluated low-radiation-dose CBCT technique can be used to measure 3D spine kinematics in vivo across functional F-E, AR, and LB positions, which has been especially challenging for the upper cervical spine.

In vivo 3-Dimensional Kinematics Study of the Healthy Cervical Spine Based on CBCT Combined with 3D-3D Registration Technology

STUDY DESIGN

A cervical biomechanical study.

OBJECTIVE

We sought to demonstrate the three-dimensional (3D) intervertebral motion characteristics of the cervical spine in healthy volunteers using cone beam computed tomography (CBCT) combined with 3D-3D registration technology.

SUMMARY OF BACKGROUND DATA

No previous studies have used CBCT combined with 3D-3D registration technology to successfully documented in vivo 3D intervertebral six-degrees-of-freedom (6-DOF) motions of the cervical spine.

METHODS

Twenty healthy subjects underwent cervical (C1-C7) CBCT scans in seven functional positions. Segmented 3D vertebral body models were established according to the cervical CBCT images. A 3D-to-3D registration was then performed for each vertebral body in the different positions to calculate the 3D segmental motion characteristics in vivo.

RESULTS

During flexion-extension, the range-of-motion (ROM) of C1-C2 and C4-C5 was significantly greater than the other segments. The average coupled axial rotation and lateral bending of each segment were between 0.6° and 3.2°. The average coupling translations in all directions were between 0.2 and 2.1 mm. During axial rotation, the ROM of C1-C2 was 65.8 ± 5.9°, which accounted for approximately 70% of all axial rotation. The motion and displacement of C1-C2 coupled lateral bending were 11.4 ± 5.2° and 8.3 ± 1.9 mm, respectively. During lateral bending, the ROM of C3-C4 was significantly greater than C1-C2, C5-C6, and C6-C7. The coupled axial rotation of C1-C2 was 34.4 ± 8.1°, and the coupled lateral translation was 3.8 ± 0.5 mm. The coupled superoinferior and anteroposterior translation of each cervical segment were between 0.1 and 0.6 mm.

CONCLUSION

CBCT combined with 3D-3D registration was used to accurately measure and record the ROMs of lateral bending, axial rotation, and flexion-extension in cervical vertebrae under physiological-load conditions. Our findings may contribute to the diagnosis of cervical spinal disease, the development of new surgical techniques, and the restoration of normal, cervical segmental movement.Level of Evidence: 3.

Surgery-related Factors Do Not Affect Short-term Adjacent Segment Kinematics After Anterior Cervical Arthrodesis

STUDY DESIGN

Prospective cohort study.

OBJECTIVE

The aim of this study was to identify surgical factors that affect adjacent segment kinematics after anterior cervical discectomy and fusion (ACDF) as measured by biplane radiography.

SUMMARY OF BACKGROUND DATA

Previous studies investigated the effect of surgical factors on spine kinematics as a potential etiology for adjacent segment disease (ASD). Those studies used static flexion-extension radiographs to evaluate range of motion. However, measurements from static radiographs are known to be unreliable. Furthermore, those studies were unable to evaluate the effect of ACDF on adjacent segment axial rotation.

METHODS

Patients had continuous cervical spine flexion/exten- sion and axial rotation movements captured at 30 images per second in a dynamic biplane radiography system preoperatively and 1 year after ACDF. Digitally reconstructed radiographs generated from subject-specific CT scans were matched to biplane radiographs using a previously validated tracking process. Dynamic kinematics, postoperative segmental kyphosis, and disc distraction were calculated from this tracking process. Plate-to-disc distance was measured on postoperative radiographs. Graft type was collected from the medical record. Multivariate linear regression was performed to identify surgical factors associated with 1-year post-surgery changes in adjacent segment kinematics. A secondary analysis was also performed to compare adjacent segment kinematics between each of the surgical factors and previously defined thresholds believed to be associated with adjacent segment degeneration.

RESULTS

Fifty-nine patients completed preoperative and postoperative testing. No association was found between any of the surgical factors and change in adjacent segment flexion/exten- sion or axial rotation range of motion (all P > 0.09). The secondary analysis also did not identify differences between adjacent segment kinematics and surgical factors (all P > 0.07).

CONCLUSION

Following ACDF for cervical spondylosis, factors related to surgical technique were not associated with short-term changes in adjacent segment kinematics that reflect the hypermobility hypothesized to lead to the development of ASD.Level of Evidence: 2.

Biomechanics of the Canine Elbow Joint

The canine elbow joint is a complex joint, whose musculoskeletal anatomy is well investigated. During the last 30 years kinematic analysis has gained importance in veterinary research and kinematics of the healthy and medial coronoid disease affected canine elbow joint are progressively investigated. Video-kinematographic analysis represents the most commonly used technique and multiple studies have investigated the range of motion, angular velocity, duration of swing and stance phase, stride length and other kinematic parameters, mostly in the sagittal plane only. However, this technique is more error-prone and data gained by video-kinematography represent the kinematics of the whole limb including the soft tissue envelope. A more precise evaluation of the in vivo bone and joint movement can only been achieved using fluoroscopic kinematography. Based on recent studies significant differences in the motion pattern between healthy joints and elbows with medial coronoid disease could be detected. Thereby not only adaptive changes, caused by pain and lameness, could be described, but primary changes in the micromotion of the joint forming bones could be found, which potentially represent new factors in the pathogenesis of medial coronoid disease. This chapter gives a review of current literature on elbow joint kinematics, with particular focus onto pathologic biomechanics in dysplastic canine elbows.

Residual Motion and Graft Type Do Not Influence Patient-reported Outcomes Following One- or Two-level Anterior Cervical Discectomy and Fusion

STUDY DESIGN

Prospective cohort.

OBJECTIVE

The aim of this study was to determine the effect of graft type on residual motion and the relationship among residual motion, smoking, and patient-reported outcome (PRO) scores following anterior cervical discectomy and fusion (ACDF).

SUMMARY OF BACKGROUND DATA

Although most patients develop solid fusion based on static imaging following ACDF, dynamic imaging has revealed that many patients continue to have residual motion at the arthrodesis.

METHODS

Forty-eight participants performed dynamic neck flexion/extension and axial rotation within a biplane radiography system 1 year following ACDF (21 one-level, 27 two-level). PRO scores included the Short Form-36, Neck Disability Index, and Cervical Spine Outcomes Questionnaire. An automated model-based tracking process matched subject-specific bone models to the biplane radiographs with sub-millimeter accuracy. Residual motion was measured across the entire arthrodesis site for both one- and two-level fusions in patients who received either allograft or autograft. Patients were divided into "pseudarthrosis" (>3° of flexion/extension residual motion) and "solid fusion" groups. Residual motion and PROs were compared between groups using Student t tests.

RESULTS

Patients who received allograft showed more total flexion/extension residual motion (4.1° vs. 2.8°, P = 0.12), although this failed to reach significance. No differences were noted in PROs based on graft type (all P > 0.08) or the presence of pseudarthrosis (all P > 0.13). No differences were noted in residual motion between smokers and nonsmokers (all P > 0.15); however, smokers who received allograft reported worse outcomes than nonsmokers who received allograft and smokers who received autograft.

CONCLUSION

Allograft may result in slightly more residual motion at the arthrodesis site 1 year after ACDF. However, there is minimal evidence that PROs are adversely affected by slightly increased residual motion, suggesting that the current definition of pseudarthrosis correlates poorly with clinically significant findings. Additionally, autograft appears to result in superior outcomes in patients who smoke.Level of Evidence: 2.

Assessing the biofidelity of in vitro biomechanical testing of the human cervical spine

In vitro biomechanical studies of the osteoligamentous spine are widely used to characterize normal biomechanics, identify injury mechanisms, and assess the effects of degeneration and surgical instrumentation on spine mechanics. The objective of this study was to determine how well four standards in vitro loading paradigms replicate in vivo kinematics with regards to the instantaneous center of rotation and arthrokinematics in relation to disc deformation. In vivo data were previously collected from 20 asymptomatic participants (45.5 ± 5.8 years) who performed full range of motion neck flexion-extension (FE) within a biplane x-ray system. Intervertebral kinematics were determined with sub-millimeter precision using a validated model-based tracking process. Ten cadaveric spines (51.8 ± 7.3 years) were tested in FE within a robotic testing system. Each specimen was tested under four loading conditions: pure moment, axial loading, follower loading, and combined loading. The in vivo and in vitro bone motion data were directly compared. The average in vitro instant center of rotation was significantly more anterior in all four loading paradigms for all levels. In general, the anterior and posterior disc heights were larger in the in vitro models than in vivo. However, after adjusting for gender, the observed differences in disc height were not statistically significant. This data suggests that in vitro biomechanical testing alone may fail to replicate in vivo conditions, with significant implications for novel motion preservation devices such as cervical disc arthroplasty implants.

Inter, intra-examiner reliability and validity of inertial sensors to measure the active cervical range of motion in patients with primary headache

We analyzed the inter- and intra-examiner reliability of Werium inertial sensors and the cervical range of motion (CROM) instrument for the measurement of active CROM (AcROM) in patients with primary headache. Another objective is to analyze the validity of the inertial sensors (Werium). The literature has reported symptomatology features in patients diagnosed with primary headache similar to that of patients with cervicogenic headache. The International Classification of Headache (ICHD-III) established the presence of reduced AcROM as a diagnostic criterion for cervicogenic headache. Several instruments are used for this measurement, with limitations in their applicability in daily clinical practice. A prospective longitudinal repeated measures study was conducted to assess the intra- and inter-rater reliability and validity of Werium inertial sensors in 20 adults with chronic primary headache. For the inter-rater analysis, the intraclass correlation coefficient (ICC) values were above 0.75 for all movements, indicating a good level of reliability. For the intra-rater results, the ICC values obtained by the Werium inertial sensors for all cervical movements were good for rater A (ICC >0.80) and rater B (ICC >0.84). For the validity, the ICCs obtained by the Werium inertial sensors compared with the CROM instrument for all cervical movements were moderate for both raters (ICC > 0.70, respectively). Values obtained in the standard error of measurement, minimum detectable change at 90% and limits of agreement also indicated good agreement. Werium inertial sensors have shown good to excellent reliability results, both intra- and inter-examiner (ICC > 0.75). Likewise, when the sensors were compared with another validated instrument (CROM device) they obtained high reliability results (ICC > 0.70). These results plus its relatively low price and ease of use allow us to recommend it in daily clinical practice to measure AcROM in patients with chronic primary headache.

The effects of age, pathology, and fusion on cervical neural foramen area

Cervical radiculopathy is a relatively common neurological disorder, often resulting from mechanical compression of the nerve root within the neural foramen. Anterior cervical discectomy and fusion (ACDF) is a common treatment for radicular symptoms that do not resolve after conservative treatment. One mechanism by which ACDF is believed to resolve symptoms is by replacing degenerated disc tissue with bone graft to increase the neural foramen area, however in vivo evidence demonstrating this is lacking. The aim of this study was to evaluate the effects of age, pathology, and fusion on bony neural foramen area. Participants included 30 young adult controls (<35 years old), 23 middle-aged controls (36 to 60 years old), and 36 cervical arthrodesis patients tested before and after ACDF surgery. Participants' cervical spines were imaged in the neutral, full flexion, and full extension positions while seated within a biplane radiography system. A validated model-based tracking technique determined three-dimensional vertebral position and orientation and automated software identified the neural foramen area in each head position. The neural foramen area decreased throughout the entire sub-axial cervical spine with age and pathology, however, no changes in neural foramen area were observed due solely to replacing degenerated disc tissue with bone graft. The neural foramen area was not associated with disc height in young adult controls, but moderate to strong associations were observed in middle-aged controls. The results provide evidence to inform the debate regarding localized versus systemic spinal degeneration and provide novel insight into the mechanism of pain relief after ACDF.

Pedicle Screws Challenged: Lumbar Cortical Density and Thickness Are Greater in the Posterior Elements Than in the Pedicles

STUDY DESIGN

Controlled laboratory study.

OBJECTIVE

To measure the total bone mineral density (BMD), cortical volume, and cortical thickness in seven different anatomical regions of the lumbar spine.

METHODS

Using computed tomography (CT) images, 3 cadaveric spines were digitally isolated by applying filters for cortical and cancellous bone. Each spine model was separated into 5 lumbar vertebrae, followed by segmentation of each vertebra into 7 anatomical regions of interest using 3-dimensional software modeling. The average Hounsfield units (HU) was determined for each region and converted to BMD with calibration phantoms of known BMD. These BMD measurements were further analyzed by the total volume, cortical volume, and cancellous volume. The cortical thickness was also measured. A similar analysis was performed by vertebral segment. St Mary's Medical Center's Institutional Review Board approved this study. No external funding was received for this work.

RESULTS

The lamina and inferior articular process contained the highest total BMD, thickest cortical shell, and largest percent volumes of cortical bone. The vertebral body demonstrated the lowest BMD. The BMDs of the L4 and L5 segments were lower; however, there were no statistically significant differences in BMD between the L1-L5 vertebral segments.

CONCLUSION

Extrapedicular regions of the lumbar vertebrae, including the lamina and inferior articular process, contain denser bone than the pedicles. Since screw pullout strength relies greatly on bone density, the lamina and inferior articular processes may offer stronger fixation of the lumbar spine.

Dynamic foraminal dimensions during neck motion 6.5 years after fusion and artificial disc replacement

Objective

To compare changes in foraminal motion at two time points post-surgery between artificial disc replacement (ADR) and anterior cervical discectomy and fusion (ACDF).

Methods

Eight ACDF and 6 ADR patients (all single-level C5-6) were tested at 2 years (T1) and 6.5 years (T2) post-surgery. The minimum foraminal height (FH.Min) and width (FW.Min) achieved during neck axial rotation and extension, and the range of these dimensions during motion (FH.Rn and FW.Rn, respectively) were measured using a biplane dynamic x-ray system, CT imaging and model-based tracking while patients performed neck axial rotation and extension tasks. Two-way mixed ANOVA was employed for analysis.

Results

In neck extension, significant interactions were found between year post-surgery and type of surgery for FW.Rn at C5-6 (p<0.006) and C6-7 (p<0.005), and for FH.Rn at C6-7 (p<0.01). Post-hoc analysis indicated decreases over time in FW.Rn for ACDF (p<0.01) and increases in FH.Rn for ADR (p<0.03) at the C6-7 adjacent level. At index level, FW.Rn was comparable between ACDF and ADR at T1, but was smaller for ACDF than for ADR at T2 (p<0.002). In axial rotation, differences were found between T1 and T2 but did not depend on type of surgery (p>0.7).

Conclusions

Changes were observed in the range of foraminal geometry at adjacent levels from 2 years to 6.5 years post-surgery that were different between ACDF and ADR for neck extension. These changes are contrary to the notion that motion at adjacent levels continue to increase following ACDF as compared to ADR over the long term.

Validation of an automated shape-matching algorithm for biplane radiographic spine osteokinematics and radiostereometric analysis error quantification

Biplane radiography and associated shape-matching provides non-invasive, dynamic, 3D osteo- and arthrokinematic analysis. Due to the complexity of data acquisition, each system should be validated for the anatomy of interest. The purpose of this study was to assess our system’s acquisition methods and validate a custom, automated 2D/3D shape-matching algorithm relative to radiostereometric analysis (RSA) for the cervical and lumbar spine. Additionally, two sources of RSA error were examined via a Monte Carlo simulation: 1) static bead centroid identification and 2) dynamic bead tracking error. Tantalum beads were implanted into a cadaver for RSA and cervical and lumbar spine flexion and lateral bending were passively simulated. A bead centroid identification reliability analysis was performed and a vertebral validation block was used to determine bead tracking accuracy. Our system’s overall root mean square error (RMSE) for the cervical spine ranged between 0.21–0.49mm and 0.42–1.80° and the lumbar spine ranged between 0.35–1.17mm and 0.49–1.06°. The RMSE associated with RSA ranged between 0.14–0.69mm and 0.96–2.33° for bead centroid identification and 0.25–1.19mm and 1.69–4.06° for dynamic bead tracking. The results of this study demonstrate our system’s ability to accurately quantify segmental spine motion. Additionally, RSA errors should be considered when interpreting biplane validation results.

In Vivo Ankle Kinematics Revealed Through Biplane Radiography: Current Concepts, Recent Literature, and Future Directions

PURPOSE OF REVIEW

Lateral ligament repair, specifically the modified Broström-Gould (BG) procedure, has been described for patients with chronic ankle instability (CAI) after failure of nonoperative management. However, there is minimal data about native in vivo ankle bone kinematics and how repairs such as the BG procedure affect the kinematics. The objective of this review is to appraise existing literature that used biplane radiography to measure in vivo kinematics of the ankle in healthy, CAI, and BG populations.

RECENT FINDINGS

Results showed that the tibiotalar joint contributes more to dorsi/plantarflexion, the subtalar joint contributes more to inversion/eversion and internal/external rotation, and that both joints are capable of complex three-dimensional (3D) motion. Preliminary data suggests that demanding activities (as opposed to walking) are necessary to elicit kinematic differences between healthy and CAI populations. Results also indicate that the BG procedure restores static kinematics and range of motion. All but one of the studies identified in this review collected static, quasi-stance, or partial gait capture data. The strength of our current knowledge is low given the small sample sizes, exploratory nature of previous work, and lack of rigorous experimental design in previous studies. Future directions include development of an improved protocol for establishing coordinate systems in the ankle bones, continued development of a database of normal kinematics during a variety of activities, and large-scale, longitudinal studies of CAI and BG patients.

Measurement of three-dimensional cervical segmental kinematics: Reliability of whole vertebrae and facet-based approaches

BACKGROUND

Previous studies have used orientation and translation of whole-vertebrae to describe three-dimensional cervical segmental kinematics. Describing kinematics using facet joint movement may be more relevant to pathology and effects of interventions but has not been investigated in the cervical spine. This study compared the reliability of two different methods (whole-vertebrae vs facet joint) to evaluate cervical kinematics.

METHODS

Two healthy adults each had six cervical (C1 to T1) magnetic resonance imaging scans, two each in neutral and left and right rotation. A semi-automated method of segmentation and alignment determined the relative orientation and translation of each whole-vertebrae and translation of each facet joint. Intra-rater and inter-rater reliability was determined using limits of agreement (LOA) with 95% confidence intervals and intraclass correlation coefficients (ICC_3,1 for intra- and ICC_2,1 for inter-rater).

RESULTS

The LOA for intra-rater evaluation of facet movement was superior to whole vertebra translation. Both methods showed excellent intra-rater ICC_3,1 (0.80-0.99) and inter-rater ICC_2,1 (0.79-0.85) for all variables except for Euler angle for flexion/extension which was good (0.65). Intra-and inter-rater ICCs were better for facet movement than all measures of whole of vertebrae movement except Euler angles of axial rotation where no difference was detected.

CONCLUSIONS

Measurement of three-dimensional segmental kinematics using either the facet joint or the whole-vertebrae method demonstrated excellent and comparable reliability. These findings support the use of the facet joint method as an option for describing and investigating cervical segmental kinematics.

Locating the Instant Center of Rotation in the Subaxial Cervical Spine with Biplanar Fluoroscopy during In Vivo Dynamic Flexion-Extension

Background

Recently, biplanar fluoroscopy is used to evaluate the cervical kinematics, especially to locate the instant center of rotation (ICR) during in vivo motion. This study aims to ascertain the ICR at each cervical segment in the sagittal plane during dynamic motion and assess the differences from previous studies.

Methods

While three healthy subjects were performing full flexion-extension, two oblique views aligned horizontally and angled at approximately 55° were obtained by biplanar fluoroscopy. The minimum degree to detect significant movement in a helical axis model was set at 2°, and anterior-posterior and superior-inferior locations of each ICR were defined. To evaluate the possible distribution area and overlapping area of the ICR with disc space, we drew a circle by using the calculated distance between each coordination and the mean coordination of ICR as the radius.

Results

During flexion-extension motion, the mean superior-inferior location of the ICR became progressively more superior, except the C5–6 segment (p = 0.015), and the mean anterior-posterior location of the ICR became progressively more anterior without exception from C2–3 to C6–7 segments, but anterior-posterior ICR locations were not significantly different among segments. The overlapping area with the distribution circle of ICR was mainly located in the posterior half in the C3–4 segment, but the overlapping area was about 80% of the total disc space in C4–5 and C6–7 segments. The overlapping was more noticeable in the lower cervical segments after exclusion of the outlier data of the C5–6 segment in subject 1.

Conclusions

The ICR in the cervical spine showed a trend of moving progressively more superiorly and anteriorly and the disc space overlapping the distribution circle of ICR increased along the lower motion segments except the C5–6 segment. These findings could provide a good basis for level-specific cervical arthroplasty designs.

Long-Duration Tracking of Cervical-Spine Kinematics With Ultrasound

Cervical-spine (C-spine) pathoanatomy is commonly evaluated by plane radiographs, computed tomography (CT), or magnetic resonance imaging (MRI); however, these modalities are unable to directly measure the dynamic mechanical properties of the functional spinal units (FSU) comprising the C-spine that account for its functional performance. We have developed an ultrasound-based technique that provides a non-invasive, real-time, quantitative, in vivo assessment of C-spine kinematics and FSU viscoelastic properties. The fidelity of the derived measurements is predicated on accurate tracking of vertebral motion over a prolonged time duration. The purpose of this work was to present a bundle adjustment method that enables accurate tracking of the relative motion of contiguous cervical vertebrae from ultrasound radio-frequency data. The tracking method was validated using both a plastic anatomical model of a cervical vertebra undergoing prescribed displacements and also human cadaveric C-spine specimens subjected to physiologically relevant loading configurations. While the velocity of motion and thickness of the surrounding soft tissue envelope affected accuracy, using the bundle adjustment method, B-mode ultrasound was capable of accurately tracking vertebral motion under clinically relevant physiologic conditions. Therefore, B-mode ultrasound can be used to evaluate in vivo real-time C-spine kinematics and FSU mechanical properties in environments where radiographs, CT, or MRI cannot be used.

Normal intervertebral segment rotation of the subaxial cervical spine: An in vivo study of dynamic neck motions

BACKGROUND

Accurate knowledge of the intervertebral center of rotation (COR) and its corresponding range of motion (ROM) can help understand development of cervical pathology and guide surgical treatment.

METHODS

Ten asymptomatic subjects were imaged using MRI and dual fluoroscopic imaging techniques during dynamic extension-flexion-extension (EFE) and axial left-right-left (LRL) rotation. The intervertebral segment CORs and ROMs were measured from C34 to C67, as the correlations between two variables were analyzed as well.

RESULTS

During the EFE motion, the CORs were located at 32.4 ± 20.6%, -2.4 ± 11.7%, 21.8 ± 12.5% and 32.3 ± 25.5% posteriorly, and the corresponding ROMs were 13.8 ± 4.3°, 15.1 ± 5.1°, 14.4 ± 7.0° and 9.2 ± 4.3° from C34 to C67. The ROM of C67 was significantly smaller than other segments. The ROMs were not shown to significantly correlate to COR locations (r = -0.243, p = 0.132). During the LRL rotation cycle, the average CORs were at 85.6 ± 18.2%, 32.3 ± 25.3%, 15.7 ± 12.3% and 82.4 ± 31.3% posteriorly, and the corresponding ROMs were 3.5 ± 1.7°, 6.9 ± 3.8°, 9.6 ± 4.1° and 2.6 ± 2.5° from C34 to C67. The ROMs of C34 and C67 was significantly smaller than those of C45 and C56. A more posterior COR was associated with a less ROM during the neck rotation (r = -0.583, p < 0.001). The ROMs during EFE were significantly larger than those during LRL in each intervertebral level.

CONCLUSION

The CORs and ROMs of the subaxial cervical intervertebral segments were segment level- and neck motion-dependent during the in-vivo neck motions.

THE TRANSLATIONAL POTENTIAL OF THIS ARTICLE

Our study indicates that the subaxial cervical intervertebral CORs and ROMs were segment level- and neck motion-dependent. This may help to improve the artificial disc design as well as surgical technique by which the neck functional motion is restored following the cervical arthroplasty.

Unloader knee brace increases medial compartment joint space during gait in knee osteoarthritis patients

PURPOSE

The purpose of the present study was to investigate the effect of the unloader brace on medial compartment dynamic joint space (DJS) during gait, while simultaneously recording ground reaction force (GRF) in varus knee osteoarthritis (OA) patients using a highly accurate biplane radiography system which allowed continuous measurement of DJS from heel strike through the midstance phase of gait. The hypothesis was that DJS in the medial compartment would be greater with the unloader brace than without the brace during gait.

METHODS

After 2 weeks of daily use of the unloader brace, ten varus knee OA patients (age 52 ± 8 years) walked with and without the brace on an instrumented treadmill, while biplane radiographs of the OA knees were acquired at 100 Hz. Medial compartment DJS was determined from heel strike to terminal stance (0-40% of the gait cycle) using a validated volumetric model-based tracking process that matches subject-specific 3D bone models to the biplane radiographs. The GRF during gait was collected at 1000 Hz. Repeated-measures ANOVA was used to explore differences in medial compartment DJS and GRF between the unbraced and braced conditions. A patient-reported subjective questionnaire related to the brace use was collected at the time of the test.

RESULTS

Medial compartment DJS was significantly greater with the unloader brace than without the brace during gait (P = 0.005). The average difference was 0.3 mm (95% confidence interval 0.1-0.4 mm). No significant difference was observed in terms of vertical GRF between the two conditions. The questionnaire showed participants felt reduced pain when wearing the brace.

CONCLUSION

The unloader brace has the significant effect of increasing medial compartment DJS during gait, which supports the underlying premise that the unloader brace reduces pain by increasing medial joint space during dynamic loading activities.

LEVEL OF EVIDENCE

III.

In Vivo Measurement of the Human Lumbar Spine Using Magnetic Resonance Imaging to Ultrasound Registration

OBJECTIVE

This study aimed to refine a magnetic resonance imaging (MRI)-ultrasound registration (ie, alignment) technique to make noninvasive, nonionizing, 3-dimensional measurement of the lumbar segmental motion in vivo.

METHODS

Five healthy participants participated in this validation study. We scanned the lumbar region of each participant 5 times using an ultrasound probe while he or she kept a prone lying posture on a plinth. Participant-specific models of L1-L5 were constructed from magnetic resonance (MR) images and aligned with the 3-dimensional ultrasound dataset of each scan using 4 variants of MRI-ultrasound registration approach (simplified intensity-based registration [1] with and [2] without including the transverse processes and their surrounding soft tissues [denoted as TP complex]; and hierarchical intensity-based registration [3] with and [4] without including the TP complex). The robustness and precision of these registration approaches were compared.

RESULTS

Although all registration approaches converged to a similar solution, excluding the TP complex improved the percentage of successful registration from 92% to 100%. There was no significant difference in the precision among the 4 MRI-ultrasound registration variants. For the simplified intensity-based registration without including the TP complex, average precision at each degree of freedom was 1.33° (flexion-extension), 2.48° (lateral bending), 1.32° (axial rotation), 2.15 mm (left/right), 1.08 mm (anterior-posterior), and 1.16 (superior-inferior), respectively.

CONCLUSION

Given that using simplified intensity-based MRI-ultrasound registration can substantially streamline the registration process and excluding the TP complex would improve the robustness of the registration, we conclude that this combination is the method of choice for in vivo human applications.

Optimization of compressive loading parameters to mimic in vivo cervical spine kinematics in vitro

The human cervical spine supports substantial compressive load in vivo. However, the traditional in vitro testing methods rarely include compressive loads, especially in investigations of multi-segment cervical spine constructs. Previously, a systematic comparison was performed between the standard pure moment with no compressive loading and published compressive loading techniques (follower load - FL, axial load - AL, and combined load - CL). The systematic comparison was structured a priori using a statistical design of experiments and the desirability function approach, which was chosen based on the goal of determining the optimal compressive loading parameters necessary to mimic the segmental contribution patterns exhibited in vivo. The optimized set of compressive loading parameters resulted in in vitro segmental rotations that were within one standard deviation and 10% of average percent error of the in vivo mean throughout the entire motion path. As hypothesized, the values for the optimized independent variables of FL and AL varied dynamically throughout the motion path. FL was not necessary at the extremes of the flexion-extension (FE) motion path but peaked through the neutral position, whereas, a large negative value of AL was necessary in extension and increased linearly to a large positive value in flexion. Although further validation is required, the long-term goal is to develop a "physiologic" in vitro testing method, which will be valuable for evaluating adjacent segment effect following spinal fusion surgery, disc arthroplasty instrumentation testing and design, as well as mechanobiology experiments where correct kinematics and arthrokinematics are critical.

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.

Asymmetry in healthy adult knee kinematics revealed through biplane radiography of the full gait cycle

One commonly used criterion in evaluating a patients' response to knee surgery or rehabilitation is bilateral symmetry. However, the natural symmetry in uninjured healthy adult knee kinematics remains relatively unknown, making it challenging to determine if clinical treatment has adequately restored bilateral symmetry. The primary purpose of this study was to determine the typical side-to-side differences in 6 degree of freedom (DOF) knee kinematics over the entire gait cycle in healthy adults using biplane radiography. Six DOF tibiofemoral kinematics were measured during treadmill walking in 19 participants using a validated volumetric model-based tracking process that matched subject-specific bone models to biplane radiographs collected at 100 images/s. Average absolute side-to-side differences in knee kinematics at foot strike were 1.3 mm or less in translation and 3.8° or less in rotation. Peak side-to-side differences in knee kinematics occurred during the swing phase and were up to 2.2 mm in translation and 7.1° in rotation. Dominant versus non-dominant leg differences were 0.8 mm and 2.8° or less at foot strike and reached maximum values of 0.8 mm and 7.2° over the full gait cycle. Statement of Clinical Significance: This study quantifies the inherent asymmetry of knee kinematics in healthy individuals over the entire gait cycle. The values of asymmetry presented here may serve as a guide for evaluating functional outcomes and restoration of so-called "normal" kinematics after injury and clinical intervention. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

Intervertebral kinematics of the cervical spine before, during, and after high-velocity low-amplitude manipulation

BACKGROUND CONTEXT

Neck pain is one of the most commonly reported symptoms in primary care settings, and a major contributor to health-care costs. Cervical manipulation is a common and clinically effective intervention for neck pain. However, the in vivo biomechanics of manipulation are unknown due to previous challenges with accurately measuring intervertebral kinematics in vivo during the manipulation.

PURPOSE

The objectives were to characterize manual forces and facet joint gapping during cervical spine manipulation and to assess changes in clinical and functional outcomes after manipulation. It was hypothesized that patient-reported pain would decrease and intervertebral range of motion (ROM) would increase after manipulation.

STUDY DESIGN/SETTING

Laboratory-based prospective observational study.

PATIENT SAMPLE

12 patients with acute mechanical neck pain (4 men and 8 women; average age 40 ± 15 years).

OUTCOME MEASURES

Amount and rate of cervical facet joint gapping during manipulation, amount and rate of force applied during manipulation, change in active intervertebral ROM from before to after manipulation, and numeric pain rating scale (NPRS) to measure change in pain after manipulation.

METHODS

Initially, all participants completed a NPRS (0-10). Participants then performed full ROM flexion-extension, rotation, and lateral bending while seated within a custom biplane radiography system. Synchronized biplane radiographs were collected at 30 images/s for 3 seconds during each movement trial. Next, synchronized, 2.0-milliseconds duration pulsed biplane radiographs were collected at 160 images/s for 0.8 seconds during the manipulation. The manipulation was performed by a licensed chiropractor using an articular pillar push technique. For the final five participants, two pressure sensors placed on the thumb of the chiropractor (Novel pliance system) recorded pressure at 160 Hz. After manipulation, all participants repeated the full ROM movement testing and once again completed the NPRS. A validated volumetric model-based tracking process that matched subject-specific bone models (from computed tomography) to the biplane radiographs was used to track bone motion with submillimeter accuracy. Facet joint gapping was calculated as the average distance between adjacent articular facet surfaces. Pre- to postmanipulation changes were assessed using the Wilcoxon signed-rank test.

RESULTS

The facet gap increased 0.9 ± 0.40 mm during manipulation. The average rate of facet gapping was 6.2 ± 3.9 mm/s. The peak force and rate of force application during manipulation were 65 ± 4 N and 440 ± 58 N/s. Pain score improved from 3.7 ± 1.2 before manipulation to 2.0 ± 1.4 after manipulation (p <. 001). Intervertebral ROM increased after manipulation by 1.2° (p = .006), 2.1° (p = .01), and 3.9° (p = .003) at the C4/C5, C5/C6, and C6/C7 motion segments, respectively, during flexion-extension; by 1.5° (p = .028), 1.9° (p = .005), and 1.3° (p = .050) at the C3/C4, C4/C5, and C5/C6 motion segments, respectively, during rotation; and by 1.3° (p = .034) and 1.1° (p = .050) at the C4/C5 and C5/C6 motion segments, respectively, during lateral bending. Global head ROM relative to the torso increased after manipulation by 8º (p = .023), 10º (p = .002), and 13º (p = .019) during lateral bending, axial rotation and flexion-extension, respectively, after manipulation.

CONCLUSIONS

This study is the first to measure facet gapping during cervical manipulation on live humans. The results demonstrate that target and adjacent motion segments undergo facet joint gapping during manipulation and that intervertebral ROM is increased in all three planes of motion after manipulation. The results suggest that clinical and functional improvement after manipulation may occur as a result of small increases in intervertebral ROM across multiple motion segments. This study demonstrates the feasibility of characterizing in real time the manual inputs and biological responses that comprise cervical manipulation, including clinician-applied force, facet gapping, and increased intervertebral ROM. This provides a basis for future clinical trials to identify the mechanisms behind manipulation and to optimize the mechanical factors that reliably and sufficiently impact the key mechanisms behind manipulation.

Trapeziometacarpal stabilization through dorsoradial ligament reconstruction: An early post-surgery in vivo biomechanical analyses

Ligament reconstruction can provide pain relief in patients with a painful, unstable, pre-arthritic trapeziometacarpal (TMC) joint. Imbrication of the dorsoradial ligament (DRL) has been proposed as a minimal invasive stabilization technique. It requires less invasive surgery than an Eaton-Littler technique and shows promising long-term clinical outcome. We used dynamic CT to objectively review the effects of the imbrication. Four patients with pain and laxity at the TMC joint, but without radiographic signs of osteoarthritis, were recruited. Dynamic CT scans were made during active thumb abduction-adduction, flexion-extension, and two functional grip tasks using a radiolucent jig. Scans of the patients were acquired before and 3 to 6 months after DRL reconstruction. Motion of each bone in the articular chain of the thumb was quantified. In addition, we mapped changes in the contact patterns between the articular facets during the entire thumb motion. After DRL imbrication, we found no overall decrease in MC1 movement in three out of four patients. Furthermore, no increase in TMC joint congruency, defined as proximity area size, was found for three out of four patients. Pre- and post-operative differences in congruency across different tasks were patient-dependent and relatively small. We demonstrated that, from a biomechanical perspective, there is high variability in post-operative outcome between patients that undergo identical surgical procedures performed by the same surgeon. A post-operative decrease in range of motion, increase in joint congruency or decrease in proximity area shift during thumb motion is not omnipresent. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:2851-2864, 2018.

In vivo validation of patellofemoral kinematics during overground gait and stair ascent

BACKGROUND

The patellofemoral (PF) joint is a common site for non-specific anterior knee pain. The pathophysiology of patellofemoral pain may be related to abnormal motion of the patella relative to the femur, leading to increased stress at the patellofemoral joint. Patellofemoral motion cannot be accurately measured using conventional motion capture.

RESEARCH QUESTION

The aim of this study was to determine the accuracy of a biplane radiography system for measuring in vivo PF motion during walking and stair ascent.

METHODS

Four subjects had three 1.0 mm diameter tantalum beads implanted into the patella. Participants performed three trials each of over ground walking and stair ascent while biplane radiographs were collected at 100 Hz. Patella motion was tracked using radiostereophotogrammetric analysis (RSA) as a "gold standard", and compared to a volumetric CT model-based tracking algorithm that matched digitally reconstructed radiographs to the original biplane radiographs.

RESULTS

The average RMS difference between the RSA and model-based tracking was 0.41 mm and 1.97° when there was no obstruction from the contralateral leg. These differences increased by 34% and 40%, respectively, when the patella was at least partially obstructed by the contralateral leg. The average RMS difference in patellofemoral joint space between tracking methods was 0.9 mm or less.

SIGNIFICANCE

Previous validations of biplane radiographic systems have estimated tracking accuracy by moving cadaveric knees through simulated motions. These validations were unable to replicate in vivo kinematics, including patella motion due to muscle activation, and failed to assess the imaging and tracking challenges related to contralateral limb obstruction. By replicating the muscle contraction, movement velocity, joint range of motion, and obstruction of the patella by the contralateral limb, the present study provides a realistic estimate of patellofemoral tracking accuracy for future in vivo studies.

Model-based tracking of the bones of the foot: A biplane fluoroscopy validation study

Measuring foot kinematics using optical motion capture is technically challenging due to the depth of the talus, small bone size, and soft tissue artifact. We present a validation of our biplane X-ray system, demonstrating its accuracy in tracking the foot bones directly. Using an experimental linear/rotary stage we imaged pairs of tali, calcanei, and first metatarsals, with embedded beads, through 30 poses. Model- and bead-based algorithms were employed for semi-automatic tracking. Translational and rotational poses were compared to the experimental stage (a reference standard) to determine registration performance. For each bone, 10 frames per pose were analyzed. Model-based: The resulting overall translational bias of the six bones was 0.058 mm with a precision of ± 0.049 mm. The overall rotational bias of the six bones was 0.42° with a precision of ± 0.41°. Bead-based: the overall translational bias was 0.037 mm with a precision of ± 0.032 mm and for rotation was 0.29° with a precision of ± 0.26°. We validated the accuracy of our system to determine the spatial position and orientation of isolated foot bones, including the talus, calcaneus, and first metatarsal over a range of quasi-static poses. Although the accuracy of dynamic motion was not assessed, use of an experimental stage establishes a reference standard.

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.

The Pedicles Are Not the Densest Regions of the Lumbar Vertebrae: Implications for Bone Quality Assessment and Surgical Treatment Strategy

STUDY DESIGN

Cadaver study.

OBJECTIVE

To determine the bone density of lumbar vertebral anatomic subregions. Bone mineral density (BMD) is a major factor in osseous fixation construct strength. The standard region for implant fixation of the spine is the pedicle; however, other regions may be more viable options with higher bone quality.

METHODS

Using computed tomography images, the spine was digitally isolated by applying a filter for adult bone. The spine model was separated into 5 lumbar vertebrae, followed by segmentation of each vertebra into 7 regions and determination of average Hounsfield units (HU). HU was converted to BMD with calibration phantoms of known BMD.

RESULTS

Overall mean BMD in vertebral regions ranged from 172 to 393 mg/cm³ with the highest and lowest BMD in the lamina and vertebral body, respectively. Vertebral regions formed 3 distinct groups (P < .03). The vertebral body and transverse processes represent one group with significantly lower BMD than other regions. Spinous process, pedicles, and superior articular processes represent a second group with moderate BMD. Finally, inferior articular process (IAP) and lamina represent a third group with significantly higher BMD than other regions.

CONCLUSIONS

Standard lumbar fusion currently uses the vertebral body and pedicles as primary locations for fixation despite their relatively low BMD. Utilization of posterior elements, especially the lamina and IAP, may be advantageous as a supplement to modern constructs or the primary site for fixation, possibly mitigating construct failures due to loosening or pullout.

Validation of 2 noninvasive, markerless reconstruction techniques in biplane high-speed fluoroscopy for 3-dimensional research of bovine distal limb kinematics

Lameness severely impairs cattle's locomotion, and it is among the most important threats to animal welfare, performance, and productivity in the modern dairy industry. However, insight into the pathological alterations of claw biomechanics leading to lameness and an understanding of the biomechanics behind development of claw lesions causing lameness are limited. Biplane high-speed fluoroscopic kinematography is a new approach for the analysis of skeletal motion. Biplane high-speed videos in combination with bone scans can be used for 3-dimensional (3D) animations of bones moving in 3D space. The gold standard, marker-based animation, requires implantation of radio-opaque markers into bones, which impairs the practicability for lameness research in live animals. Therefore, the purpose of this study was to evaluate the comparative accuracy of 2 noninvasive, markerless animation techniques (semi-automatic and manual) in 3D animation of the bovine distal limb. Tantalum markers were implanted into each of the distal, middle, and proximal phalanges of 5 isolated bovine distal forelimbs, and biplane high-speed x-ray videos of each limb were recorded to capture the simulation of one step. The limbs were scanned by computed tomography to create bone models of the 6 digital bones, and 3D animation of the bones' movements were subsequently reconstructed using the marker-based, the semi-automatic, and the manual animation techniques. Manual animation translational bias and precision varied from 0.63 ± 0.26 mm to 0.80 ± 0.49 mm, and rotational bias and precision ranged from 2.41 ± 1.43° to 6.75 ± 4.67°. Semi-automatic translational values for bias and precision ranged from 1.26 ± 1.28 mm to 2.75 ± 2.17 mm, and rotational values varied from 3.81 ± 2.78° to 11.7 ± 8.11°. In our study, we demonstrated the successful application of biplane high-speed fluoroscopic kinematography to gait analysis of bovine distal limb. Using the manual animation technique, kinematics can be measured with sub-millimeter accuracy without the need for invasive marker implantation.

Dynamic in vivo 3D atlantoaxial spine kinematics during upright rotation

Diagnosing dysfunctional atlantoaxial motion is challenging given limitations of current diagnostic imaging techniques. Three-dimensional imaging during upright functional motion may be useful in identifying dynamic instability not apparent on static imaging. Abnormal atlantoaxial motion has been linked to numerous pathologies including whiplash, cervicogenic headaches, C2 fractures, and rheumatoid arthritis. However, normal C1/C2 rotational kinematics under dynamic physiologic loading have not been previously reported owing to imaging difficulties. The objective of this study was to determine dynamic three-dimensional in vivo C1/C2 kinematics during upright axial rotation. Twenty young healthy adults performed full head rotation while seated within a biplane X-ray system while radiographs were collected at 30 images per second. Six degree-of-freedom kinematics were determined for C1 and C2 via a validated volumetric model-based tracking process. The maximum global head rotation (to one side) was 73.6±8.3°, whereas maximum C1 rotation relative to C2 was 36.8±6.7°. The relationship between C1/C2 rotation and head rotation was linear through midrange motion (±20° head rotation from neutral) in a nearly 1:1 ratio. Coupled rotation between C1 and C2 included 4.5±3.1° of flexion and 6.4±8.2° of extension, and 9.8±3.8° of contralateral bending. Translational motion of C1 relative to C2 was 7.8±1.5mm ipsilaterally, 2.2±1.2mm inferiorly, and 3.3±1.0mm posteriorly. We believe this is the first study describing 3D dynamic atlantoaxial kinematics under true physiologic conditions in healthy subjects. C1/C2 rotation accounts for approximately half of total head axial rotation. Additionally, C1 undergoes coupled flexion/extension and contralateral bending, in addition to inferior, lateral and posterior translation.

Validation of static and dynamic radiostereometric analysis of the knee joint using bone models from CT data

Objectives

Static radiostereometric analysis (RSA) using implanted markers is considered the most accurate system for the evaluation of prosthesis migration. By using CT bone models instead of markers, combined with a dynamic RSA system, a non-invasive measurement of joint movement is enabled. This method is more accurate than current 3D skin marker-based tracking systems. The purpose of this study was to evaluate the accuracy of the CT model method for measuring knee joint kinematics in static and dynamic RSA using the marker method as the benchmark.

Methods

Bone models were created from CT scans, and tantalum beads were implanted into the tibia and femur of eight human cadaver knees. Each specimen was secured in a fixture, static and dynamic stereoradiographs were recorded, and the bone models and marker models were fitted to the stereoradiographs.

Results

Results showed a mean difference between the two methods in all six degrees of freedom for static RSA to be within -0.10 mm/° and 0.08 mm/° with a 95% limit of agreement (LoA) ranging from ± 0.49 to 1.26. Dynamic RSA had a slightly larger range in mean difference of -0.23 mm/° to 0.16 mm/° with LoA ranging from ± 0.75 to 1.50.

Conclusions

In a laboratory-controlled setting, the CT model method combined with dynamic RSA may be an alternative to previous marker-based methods for kinematic analyses.

Cite this article: K. Stentz-Olesen, E. T. Nielsen, S. De Raedt, P. B. Jørgensen, O. G. Sørensen, B. L. Kaptein, M. S. Andersen, M. Stilling. Validation of static and dynamic radiostereometric analysis of the knee joint using bone models from CT data. Bone Joint Res 2017;6:376–384. DOI: 10.1302/2046-3758.66.BJR-2016-0113.R3.

Three-dimensional intervertebral range of motion in the cervical spine: Does the method of calculation matter?

Intervertebral range of motion (ROM) is commonly calculated using ordered rotations or projection angles. Ordered rotations are sequence-dependent, and projection angles are dependent upon on which orientation vectors are projected. This study assessed the effect of calculation method on intervertebral ROM in the subaxial cervical spine (C3-C7) during in vivo dynamic, three-dimensional, functional movement. Biplane radiographs were collected at 30 images per second while 29 participants performed full ROM flexion/extension, axial rotation and lateral bending movements of their cervical spine. In vivo bone motion was tracked with sub-millimeter accuracy using a validated volumetric model-based tracking technique. Intervertebral rotations were calculated using six Cardan angle sequences and two projection angle combinations. Within-subject comparisons revealed significant differences in intervertebral ROM among calculation methods (all p<0.002). Group mean ROM differences were small, but significantly different among calculation methods (p<0.001). A resampling technique demonstrated that as group size increases, the differences between calculation methods decreases substantially. It is concluded that the method used to calculate intervertebral rotations of the sub-axial cervical spine can significantly affect within-subject and between group comparisons of intervertebral ROM.

Longitudinal Study of the Six Degrees of Freedom Cervical Spine Range of Motion During Dynamic Flexion, Extension, and Rotation After Single-level Anterior Arthrodesis

STUDY DESIGN

A longitudinal study using biplane radiography to measure in vivo intervertebral range of motion (ROM) during dynamic flexion/extension, and rotation.

OBJECTIVE

To longitudinally compare intervertebral maximal ROM and midrange motion in asymptomatic control subjects and single-level arthrodesis patients.

SUMMARY OF BACKGROUND DATA

In vitro studies consistently report that adjacent segment maximal ROM increases superior and inferior to cervical arthrodesis. Previous in vivo results have been conflicting, indicating that maximal ROM may or may not increase superior and/or inferior to the arthrodesis. There are no previous reports of midrange motion in arthrodesis patients and similar-aged controls.

METHODS

Eight single-level (C5/C6) anterior arthrodesis patients (tested 7 ± 1 months and 28 ± 6 months postsurgery) and six asymptomatic control subjects (tested twice, 58 ± 6 months apart) performed dynamic full ROM flexion/extension and axial rotation whereas biplane radiographs were collected at 30 images per second. A previously validated tracking process determined three-dimensional vertebral position from each pair of radiographs with submillimeter accuracy. The intervertebral maximal ROM and midrange motion in flexion/extension, rotation, lateral bending, and anterior-posterior translation were compared between test dates and between groups.

RESULTS

Adjacent segment maximal ROM did not increase over time during flexion/extension, or rotation movements. Adjacent segment maximal rotational ROM was not significantly greater in arthrodesis patients than in corresponding motion segments of similar-aged controls. C4/C5 adjacent segment rotation during the midrange of head motion and maximal anterior-posterior translation were significantly greater in arthrodesis patients than in the corresponding motion segment in controls on the second test date.

CONCLUSION

C5/C6 arthrodesis appears to significantly affect midrange, but not end-range, adjacent segment motions. The effects of arthrodesis on adjacent segment motion may be best evaluated by longitudinal studies that compare maximal and midrange adjacent segment motion to corresponding motion segments of similar-aged controls to determine if the adjacent segment motion is truly excessive.

LEVEL OF EVIDENCE

3.