The relationship between concave angle of vertebral endplate and lumbar intervertebral disc degeneration

Augmenting Screw Technique to Prevent TLIF Cage Subsidence: A Biomechanical In Vitro Study

(1) Cage subsidence in spine surgery is a frequent clinical challenge. This study aimed to assess a novel screw augmentation technique for Transforaminal Lumbar Interbody Fusion in cadavers of reduced bone mineral density (BMD). (2) Forty human lumbar vertebrae (BMD 84.2 ± 24.4 mgHA/cm3, range 51-119 mgHA/cm3) were assigned to two groups: augmenting screw group and control group. The augmentation technique comprised placement of two additional subcortical screws. Ten constructs per group were loaded with a quasi-static load-to-failure protocol and other ten were cyclically loaded. Failure modes were documented. (3) During the quasi-static load-to-failure testing, the augmenting screw technique showed a significantly higher failure load (1426.0 ± 863.6 N) versus the conventional technique in the control group (682.2 ± 174.5 N, p = 0.032). Cyclic loading revealed higher number of cycles and corresponding load until reaching 5 mm subsidence and significantly higher number of cycles and corresponding load until reaching 10 mm subsidence for the augmenting screw technique (9645 ± 3050; 1164.5 ± 305.0 N) versus the conventional technique in the control group (5395 ± 2340; 739.5 ± 234.0 N, p < 0.05). Failure modes were different and showed bending of the augmenting screws, followed by cut-out. (4) The investigated augmenting screw technique demonstrated higher failure loads and cycles to failure against cage subsidence compared to conventional cage placement. Failure modes were different between the two techniques and may lead to a different kind of complications.

Development of Modic Changes After Percutaneous Endoscopic Transforaminal Lumbar Discectomy: From Risk Analysis to Prediction Modeling

Objective

This study examines the occurrence of Modic changes (MC) within the first year following percutaneous endoscopic transforaminal lumbar discectomy (PETD) and investigates associated risk factors.

Methods

This study adopted a retrospective cohort design. Between January 2019 and June 2023, 538 patients diagnosed with single-level lumbar disc herniation and treated with PETD were included. The patients were divided into a training set and a validation set based on their surgery dates. Preoperative radiographic parameters and perioperative indicators were evaluated. Univariate analysis examined risk factors for postoperative MC. Gender-specific subgroups were analyzed. Binary logistic regression developed a predictive model for postoperative MC, assessed using ROC, calibration, and decision curves.

Results

The incidence of MC at one year after PETD was 24.8%. Logistic regression identified 8 significant risk factors for MC after PELD: longer symptom duration, proximity of herniated segment to sacrum, severe disc degeneration, reduced disc height, greater vertebral endplate concavity angle, segmental instability, and lumbar-sacral fusion. Menopause and herniation type were identified as female-specific risk factors. In males, total cholesterol levels were additionally found to be a risk factor for postoperative MC. The male and female subgroup models exhibited satisfactory performance across ROC analysis, calibration plots, and decision curve analysis. Specifically, for male patients, the area under the curve (AUC) was 0.831 for the training set and 0.820 for the validation set; for female patients, the AUC was 0.911 for the training set and 0.868 for the validation set. A nomogram was developed to visualize the model.

Conclusion

This study explored the relevant risk factors of MC after PETD and visualized the prediction model by nomogram, which is beneficial to optimize the surgical scheme of PETD to improve the clinical efficacy.

Association between vertebral endplate defects and patient-reported symptoms: an immunohistochemical study investigating the COX-2/PGE-2/EP-4 axis

BACKGROUND CONTEXT

Vertebral endplate defects are often implicated in degenerative disc disorders, yet their connection to patient-reported symptoms remains unclear. COX-2 and PGE-2 are known for their roles in inflammation and pain, with EP-4 receptor involvement in pain signaling. Examining their expression in vertebral endplate tissues may provide insights into pathomechanism of low back pain.

PURPOSE

To investigate the association between endplate defects and patient-reported symptoms and to further clarify the role of the COX-2/PGE-2/EP-4 axis in the pathogenesis of chronic low back pain.

STUDY DESIGN/SETTING

Retrospective study.

PATIENT SAMPLE

A total of 71 patients who had undergone single-level L4/5 or L5/S1 modified laminectomy decompression preserving proximal upper laminae and transforaminal lumbar interbody fusion surgery were included in this study, including 18 patients diagnosed with lumbar disc herniation, 19 with lumbar disc herniation accompanied by degenerative lumbar spinal stenosis, and 34 with degenerative spondylolisthesis.

OUTCOME MEASURES

Demographic data, Pfirrmann grade, Modic changes, endplate defect score, visual analog scale (VAS) for back and leg pain, and Oswestry Disability Index (ODI) before surgery, 3-month and 6-month follow-up, and the percentage of immune-positive cells (COX-2, PGE-2, and EP-4) in endplate tissue sections.

METHODS

Patients were divided into defect and nondefect groups according to endplate morphology on lumbar MR. All intraoperative endplate specimens were immediately fixed in 10% formaldehyde, and then embedded in paraffin 3 days later for tissue sections. The outcome measures were compared between the defect group and nondefect group. Data were analyzed using independent t-tests and χ² tests. Pearson's rank correlation test was used to assess correlations between patient-reported symptoms and the percentage of immune-positive cells in the groups. Multivariable logistic regression models using the forward stepwise likelihood ratio method were used to identify the factors that were independently associated with endplate defects.

RESULTS

The age of Defect group was significantly higher than that of nondefect group (52.5±7.7 vs 57.2±9.1. p=.024). There were no significant differences in gender, diagnosis, BMI, comorbidities, or surgical level between the two groups. Modic changes (Type Ⅱ/Type Ⅲ) were more common in patients of Defect group than nondefect group (38.5% vs 11.1%, p<.001), and so was disc degeneration (Pfirrmann grade Ⅳ/Ⅴ) (69.2% vs 33.3%, p<.001). Defect group had significantly higher VAS-Back (6.5±2.0 vs 4.9±1.6, p<.001) and ODI scores (62.9±10.7 vs 45.2±14.8, p<.001) than nondefect group, while there was no significant differences between the two groups during the 3 and 6-month follow-up after surgery. Histologically, Defect group was characterized by upregulation of COX-2, PGE-2, and EP-4 in endplate tissue sections. Both in defect and nondefect groups, VAS-Back showed moderate positive correlations with the expressions of COX-2 (r=0.643; r=0.558, p both<.001), PGE-2 (r=0.611; r=0.640, p both<.001), and EP-4 (r=0.643; r=0.563, p both<.001). Multivariate regression analyses reveled that percentage of COX-2-positive cells was associated with endplate defects (OR=1.509, 95%CI [1.048-2.171], p=.027), as well as percentage of PGE-2-positive (OR=1.291, 95%CI [1.106-1.508], p=.001) and EP-4-positive cells (OR=1.284, 95%CI [1.048∼2.171], p=.003).

CONCLUSIONS

Patients with endplate defects had worse quality of life, more severe disc degeneration and Modic changes, and up-regulated COX-2/PGE-2/EP-4 axis expression in cartilage endplates in patients with defected endplates. Inflammatory factors may significantly contribute to the onset and progression of chronic low back pain in patients with endplate defects, consequently impacting patient-reported symptoms.

Correlation between sagittal morphology of lower lumbar end plate and degenerative changes in patients with lumbar disc herniation

Objective

As an important anatomic factor in the process of lumbar disc herniation (LDH), the correlation between end plate sagittal morphology and intervertebral disc degeneration (IDD) is unclear. Moreover, research on imaging data of lumbar end plate in patients with LDH is still insufficient. Our study aimed to observe the morphological change of the lower lumbar end plate (L3-S1) in patients with LDH on magnetic resonance imaging (MRI) and analyze its correlation with the degree of IDD.

Materials and Methods

A total of 116 patients were included in the study. Based on their MRI, we divided end plates into three types (concave, flat, and irregular), assigned intervertebral discs with Grade I-V given 1-5 points successively according to the Pfirrmann system, and determined whether there was Modic change of each end plate. The correlation between the morphology of the end plate and the degree of IDD was analyzed.

Results

There was an excellent interobserver agreement for each item we analyzed (interclass correlation coefficient >0.75). Concave end plate appeared most frequently (187, 53.7%) and was mainly distributed in L3/4 and L4/5, whereas irregular end plate was the least common type (54, 15.5%) and mainly concentrated in L5/S1. The IDD degree of the corresponding disc increased gradually from concave (3.27 ± 0.81) to irregular end plates (4.25 ± 0.79) (P < 0.05). Irregular end plates were more likely to have Modic changes than concave and flat end plates (P < 0.05).

Conclusion

The sagittal morphology of the lower lumbar end plate is related to modic changes and degree of IDD (based on the Pfirrmann grading system) in patients with LDH, and the concave end plate mostly reflects a lower degree of lumbar disc degeneration, which has substantial clinical significance.

Comparison of the biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients: A finite element analysis

BACKGROUND

Some older patients who suffered from both conditions (disc degeneration and osteoporosis) have higher surgical risks and longer postoperative recovery times. Understanding the relation between disc degeneration and osteoporosis is fundamental to know the mechanisms of orthopedic disorders and improve clinical treatment. However, there is a lack of finite element (FE) studies to predict the combined effects of disc degeneration and osteoporosis. So the aim of the present study is to explore the differences of biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients.

METHODS

A normal lumbar spine finite element model (FEM) was developed based on the geometric information of a healthy male subject (age 35 years; height 178 cm; weight 65 kg). This normal lumbar spine FEM was modified to build three lumbar spine degeneration models simulating mild, moderate and severe grades of disc degeneration at the L4-L5 segment. Then the degenerative lumbar spine models for osteoporotic patients were constructed on the basis of the above-mentioned degeneration models. Firstly, the normal model (flexion: 8 Nm; extension: 6 Nm; lateral bending: 6 Nm; torsion: 4 Nm) and degenerative models (10 Nm) were calibrated under pure moment load, respectively. Secondly, under a 400 N follower load, the 7.5 Nm moments of different directions were applied on all models to simulate different motion postures. Finally, under the above loading conditions, we calculated and analyzed the range of motion (ROM), Mises stress in cortical (MSC1), Mises stress in endplate (MSE), Mises stress in cancellous (MSC2), and Mises stress in post (MSP).

RESULTS

Compared with disc degeneration patients without osteoporosis, the ROM, MSC1, and MSE of osteoporosis patients with various disc degeneration decreased in all postures, while the MSC2 and MSP increased. With increase in the degree of disc degeneration, the reduction proportions of ROM and MSE in osteoporotic patients gradually increased, while the reduction percentages in MSC1 of osteoporotic patients gradually decreased. The increase percentages of MSC2 in osteoporotic patients gradually increased. Given the progressive changes of disc degeneration, the changes in MSP in osteoporosis patients were uneven.

CONCLUSION

In summary, the effect of disc degeneration on flexibility in the two kinds of patients (osteoporosis and non-osteoporosis patients) was nearly same. By comparing the remaining biomechanical parameters (MSC1, MSE, MSC2, and MSP), we found that degenerated intervertebral discs caused changes in loading patterns of osteoporosis patients. Disc degeneration reduced the Mises stress in the cortical and endplate, which increased the Mises stress in the cancellous and post. That is to say, in order to cope with the changes in bone stresses caused by disc degeneration and osteoporosis, clinicians should be more careful in choosing the surgical option for osteoporotic patients with disc degeneration.

Spinopelvic alignment and lumbar vertebral shape in children: associations with structural spinal abnormalities and body composition in the generation R study

PURPOSE

To investigate the spinopelvic alignment and vertebral shape in children, and associations with body composition and structural spinal abnormalities on magnetic resonance imaging (MRI).

METHODS

We performed a cross-sectional study embedded in the Generation R Study, a prospective population-based birth cohort. Pelvic incidence and vertebral concavity ratios for each lumbar level were determined on sagittal MRI images in 9-year-old children, and structural spinal abnormalities were scored semi-quantitatively. The BMI-SD score was calculated, and body composition was assessed using DXA scans. Associations of pelvic incidence and vertebral concavity ratios with structural abnormalities and body composition measures were assessed using (multilevel) regression analyses.

RESULTS

This study included 522 participants (47.7% boys), aged 9.9 years (IQR 9.7-10.0). The mean pelvic incidence was 36.6° (SD 8.0). Vertebral concavity ratios ranged from 0.87 to 0.90, with significantly lower ratios for boys compared to girls. Associations were found for a larger pelvic incidence with decreased disc height [OR 1.03 (95% CI 1.02-1.05)], and a pelvic incidence in the lowest tertile with less disc bulging [OR 0.73 (95% CI 0.56-0.95)]. Increased vertebral concavity ratio was associated with decreased disc height [OR 14.16 (95% CI 1.28-157.13)]. Finally, increased fat-free mass index was associated with a smaller pelvic incidence [adjusted OR 0.85 (95% CI 0.07-1.63)].

CONCLUSION

The mean pelvic incidence of 9-year-old children is 36.6° on supine MRI images, and a slightly concave shape of the lumbar vertebrae is seen. Spinopelvic alignment is associated with structural spinal abnormalities, and might itself be influenced by the children's body composition.

Diagnostic Role of Magnetic Resonance Imaging in Low Back Pain Caused by Vertebral Endplate Degeneration

Low back pain (LBP) is a common health issue worldwide with a huge economic burden on healthcare systems. In the United States alone, the cost is estimated to be $100 billion each year. Intervertebral disc degeneration is considered one of the primary causes of LBP. Moreover, the critical role of the vertebral endplates in disc degeneration and LBP is becoming apparent. Endplate abnormalities are closely correlated with disc degeneration and pain in the lumbar spine. Imaging modalities such as plain film radiography, computed tomography, and fluoroscopy are helpful but not very effective in detecting the causes behind LBP. Magnetic resonance imaging (MRI) can be used to acquire high-quality three-dimensional images of the lumbar spine without using ionizing radiation. Therefore, it is increasingly being used to diagnose spinal disorders. However, according to the American College of Radiology, current referral and justification guidelines for MRI are not sufficiently clear to guide clinical practice. This review aimed to evaluate the role of MRI in diagnosing LBP by considering the correlative contributions of vertebral endplates. The findings of the review indicate that MRI allows for fine evaluations of endplate morphology, endplate defects, diffusion and perfusion properties of the endplate, and Modic changes. Changes in these characteristics of the endplate were found to be closely correlated with disc degeneration and LBP. The collective evidence from the literature suggests that MRI may be the imaging modality of choice for patients suffering from LBP. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 3.

Nomograms for Predicting Recurrent Herniation in PETD with Preoperative Radiological Factors

Purpose

To investigate the preoperative radiological risk factors for recurrent lumbar disc herniation (rLDH) within 1 year after percutaneous endoscopic transforaminal discectomy (PETD).

Patients and Methods

A retrospective case–control study was conducted. Between January 2013 and November 2019, a total of 1210 patients with single-level L4/5 LDH who underwent PETD were enrolled in the present study. In total, 62 rLDH patients were diagnosed and collected based on the clinical and radiological manifestations, and 224 non-rLDH controls were selected from the remaining 1148 patients. Preoperative radiological parameters were collected and measured. An age threshold to distinguish patients into subgroups was established using the Youden index. The relationships between the risk factors and rLDH were evaluated by univariate and multivariate analyses in two subgroups. Predictive models were established based on logistic analysis. The area under the curve (AUC) of receiver operating characteristic (ROC) curves, calibration plots, and decision curve analysis (DCA) were used to assess the predictive models.

Results

In the present study, logistic analysis identified six significant predictors associated with rLDH in the young group: superior endplate concave angle (ECA), sacral slope, Modic changes, sagittal range of motion (sROM), extension intervertebral angle (IVA), and lumbar lordosis. Four significant predictors were identified in the elderly group: disc height index (DHI), retrolisthesis (posterior spondylolisthesis), sROM, and extension IVA. Validation of both models demonstrated excellent model discrimination (AUC= 0.940 and 0.946, respectively). DCA also showed excellent clinical utility and benefits.

Conclusion

The nomograms that we constructed could accurately predict individual patient recurrence risk. Individualized measures should be taken for patients of different ages with the above risk factors, and tailored postoperative surveillance of patients who underwent PETD can be planned.

The effect of cervical intervertebral disc degeneration on the motion path of instantaneous center of rotation at degenerated and adjacent segments: A finite element analysis

BACKGROUND

The motion path of instantaneous center of rotation (ICR) is a crucial kinematic parameter to dynamically characterize cervical spine intervertebral patterns of motion; however, few studies have evaluated the effect of cervical disc degeneration (CDD) on ICR motion path. The purpose of this study was to investigate the effect of CDD on the ICR motion path of degenerated and adjacent segments.

METHOD

A validated nonlinear three-dimensional finite element (FE) model of a healthy adult cervical spine was used. Progressive degeneration was simulated with six FE models by modifying intervertebral disc height and material properties, anterior osteophyte size, and degree of endplate sclerosis at the C5-C6 level. All models were subjected to a pure moment of 1 Nm and a compressive follower load of 73.6 N to simulate physical motion. ICR motion paths were compared among different models.

RESULTS

The normal FE model results were consistent with those of previous studies. In degenerative models, average ICR motion paths shifted significantly anterior at the degenerated segment (β = 0.27 mm; 95% CI: 0.22, 0.32) and posterior at the proximal adjacent segment (β = -0.09 mm; 95% CI: -0.15, -0.02) than those of the normal model.

CONCLUSION

CDD significantly affected ICR motion paths at the degenerated and proximal adjacent segments. The changes at adjacent segments may be a result of compensatory mechanisms to maintain the balance of the cervical spine. Surgical treatment planning should take into account the restoration of ICR motion path to normal. These findings could provide a basis for prosthesis design and clinical practice.

A New Hope in Spinal Degenerative Diseases: Piezo1

As a newly discovered mechanosensitive ion channel protein, the piezo1 protein participates in the transmission of mechanical signals on the cell membrane and plays a vital role in mammalian biomechanics. Piezo1 has attracted widespread attention since it was discovered in 2010. In recent years, studies on piezo1 have gradually increased and deepened. In addition to the discovery that piezo1 is expressed in the respiratory, cardiovascular, gastrointestinal, and urinary systems, it is also stably expressed in cells such as mesenchymal stem cells (MSCs), osteoblasts, osteoclasts, chondrocytes, and nucleus pulposus cells that constitute vertebral bodies and intervertebral discs. They can all receive external mechanical stimulation through the piezo1 protein channel to affect cell proliferation, differentiation, migration, and apoptosis to promote the occurrence and development of lumbar degenerative diseases. Through reviewing the relevant literature of piezo1 in the abovementioned cells, this paper discusses the effect of piezo1 protein expression under mechanical stress stimuli on spinal degenerative disease, providing the molecular basis for the pathological mechanism of spinal degenerative disease and also a new basis, ideas, and methods for the prevention and treatment of this degenerative disease.

Study of global dynamics of COVID-19 via a new mathematical model

The theme of this paper focuses on the mathematical modeling and transmission mechanism of the new Coronavirus shortly noted as (COVID-19), endangering the lives of people and causing a great menace to the world recently. We used a new type epidemic model composed on four compartments that is susceptible, exposed, infected and recovered (SEIR), which describes the dynamics of COVID-19 under convex incidence rate. We simulate the results by using nonstandard finite difference method (NSFDS) which is a powerful numerical tool. We describe the new model on some random data and then by the available data of a particular regions of Subcontinents.

Biomechanical Effect of L4 -L5 Intervertebral Disc Degeneration on the Lower Lumbar Spine: A Finite Element Study

OBJECTIVE

To ascertain the biomechanical effects of a degenerated L₄ -L₅ segment on the lower lumbar spine through a comprehensive simulation of disc degeneration.

METHODS

A three-dimensional nonlinear finite element model of a normal L₃ -S₁ lumbar spine was constructed and validated. This normal model was then modified such that three degenerated models with different degrees of degeneration (mild, moderate, or severe) at the L₄ -L₅ level were constructed. While experiencing a follower compressive load (500 N), hybrid moment loads were applied to all models to determine range of motion (ROM), intradiscal pressure (IDP), maximum von Mises stress in the annulus, maximum shear stress in the annulus, and facet joint force.

RESULTS

As the degree of disc degeneration increased, the ROM of the L₄ -L₅ degenerated segment declined dramatically in all postures (flexion: 5.79°-1.91°; extension: 5.53°-2.62°; right lateral bending: 4.47°-1.46°; left lateral bending: 4.86°-1.61°; right axial rotation: 2.69°-0.74°; left axial rotation: 2.69°-0.74°), while the ROM in adjacent segments increased (1.88°-8.19°). The largest percent decrease in motion of the L₄ -L₅ segment due to disc degeneration was in right axial rotation (75%), left axial rotation (69%), flexion (67%), right lateral bending (67%), left lateral bending right (67%), and extension (53%). The change in the trend of the IDP was the same as that of the ROM. Specifically, the IDP decreased (flexion: 0.592-0.09 MPa; extension: 0.678-0.334 MPa; right lateral bending: 0.498-0.205 MPa; left lateral bending: 0.523-0.272 MPa; right axial rotation: 0.535-0.246 MPa; left axial rotation: 0.53-0.266 MPa) in the L₄ -L₅ segment, while the IDP in adjacent segments increased (0.511-0.789 MPa). The maximum von Mises stress and maximum shear stress of the annulus in whole lumbar spine segments increased (L₄ -L₅ segment: 0.413-2.626 MPa and 0.412-2.783 MPa, respectively; adjacent segment of L₄ -L₅ : 0.356-1.493 MPa and 0.359-1.718 MPa, respectively) as degeneration of the disc progressively increased. There was no apparent regularity in facet joint force in the degenerated segment as the degree of disc degeneration increased. Nevertheless, facet joint forces in adjacent healthy segments increased as the degree of disc degeneration increased (extension: 49.7-295.3 N; lateral bending: 3.5-171.2 N; axial rotation: 140.2-258.8 N).

CONCLUSION

Degenerated discs caused changes in the motion and loading pattern of the degenerated segments and adjacent normal segments. The abnormal load and motion in the degenerated models risked accelerating degeneration in the adjacent normal segments. In addition, accurate simulation of degenerated facet joints is essential for predicting changes in facet joint loads following disc degeneration.

Spondylolysis and spinal adaptations for bipedalism: The overshoot hypothesis

Background and objectives

The study reported here focused on the aetiology of spondylolysis, a vertebral pathology usually caused by a fatigue fracture. The goal was to test the Overshoot Hypothesis, which proposes that people develop spondylolysis because their vertebral shape is at the highly derived end of the range of variation within Homo sapiens.

Methodology

We recorded 3D data on the final lumbar vertebrae of H. sapiens and three great ape species, and performed three analyses. First, we compared H. sapiens vertebrae with and without spondylolysis. Second, we compared H. sapiens vertebrae with and without spondylolysis to great ape vertebrae. Lastly, we compared H. sapiens vertebrae with and without spondylolysis to great ape vertebrae and to vertebrae of H. sapiens with Schmorl’s nodes, which previous studies have shown tend to be located at the ancestral end of the range of H. sapiens shape variation.

Results

We found that H. sapiens vertebrae with spondylolysis are significantly different in shape from healthy H. sapiens vertebrae. We also found that H. sapiens vertebrae with spondylolysis are more distant from great ape vertebrae than are healthy H. sapiens vertebrae. Lastly, we found that H. sapiens vertebrae with spondylolysis are at the opposite end of the range of shape variation than vertebrae with Schmorl’s nodes.

Conclusions

Our findings indicate that H. sapiens vertebrae with spondylolysis tend to exhibit highly derived traits and therefore support the Overshoot Hypothesis. Spondylolysis, it appears, is linked to our lineage’s evolutionary history, especially its shift from quadrupedalism to bipedalism.

Lay summary: Spondylolysis is a relatively common vertebral pathology usually caused by a fatigue fracture. There is reason to think that it might be connected with our lineage’s evolutionary shift from walking on all fours to walking on two legs. We tested this idea by comparing human vertebrae with and without spondylolysis to the vertebrae of great apes. Our results support the hypothesis. They suggest that people who experience spondylolysis have vertebrae with what are effectively exaggerated adaptations for bipedalism.

On modeling of coronavirus-19 disease under Mittag-Leffler power law

This paper investigates a new model on coronavirus-19 disease (COVID-19) with three compartments including susceptible, infected, and recovered class under Mittag-Leffler type derivative. The mentioned derivative has been introduced by Atangana, Baleanu, and Caputo abbreviated as ( ABC ) . Upon utilizing fixed point theory, we first prove the existence of at least one solution for the considered model and its uniqueness. Also, some results about stability of Ulam-Hyers type are also established. By applying a numerical technique called fractional Adams-Bashforth (AB) method, we develop a scheme for the approximate solutions to the considered model. Using some real available data, we perform the concerned numerical simulation corresponding to different values of fractional order.

Three-Dimensional Morphometric Analysis of Lumbar Vertebral End Plate Anatomy

BACKGROUND

Information on the three-dimensional (3D) shape of vertebral end plates is lacking. Previous studies have analyzed two-dimensional shape; however, 3D data are important because they may help improve our understanding of how differences in shape are related to age, gender, race, size, and other parameters, which may subsequently help improve device design for interbody prosthesis.

OBJECTIVE

To study the 3D shape of lumbar vertebral end plates from normal adult lumbar spines and correlate them with age, gender, spinal/end plate level, end plate surface area, concave depth, and size.

METHODS

An in vivo analysis was undertaken of lumbar vertebral end plate 3D shape. A total of 136 patients' computed tomography scans were used to create 3D models of the lumbar spine for each patient, which were subsequently analyzed.

RESULTS

The shape of the superior end plates is different compared with inferior end plates. Across the lumbar spine (L1-S1), the shape of inferior end plates is similar; however, the shape of the superior end plate varies between spinal levels significantly. There was no clear relationship between age and principal component (PC) shapes but there was a strong correlation between end plate concave depth and end plate PC shape scores.

CONCLUSIONS

Future interbody (disc replacement and fusion) device designs could use the findings that inferior end plate shape is similar throughout the length of the lumbar spine, whereas superior end plate shape changes. Further, future implants could be level-specific because the present study shows that end plate shape varies through the length of the lumbar spine.

Effect of disc degeneration on the mechanical behavior of the human lumbar spine: a probabilistic finite element study

BACKGROUND CONTEXT

Intervertebral disc degeneration has been subject to numerous in vivo and in vitro investigations and numerical studies during recent decades, reporting partially contradictory findings. However, most of the previous studies were limited in the number of specimens investigated and, therefore, could not consider the vast variety of the specimen geometries, which are likely to strongly influence the mechanical behavior of the spine.

PURPOSE

To complement the understanding of the mechanical consequences of disc degeneration, whereas considering natural variations in the major spinal geometrical parameters.

DESIGN/SETTING

A probabilistic finite element study.

METHODS

A parametric finite element model of a human L4-L5 motion segment considering 40 geometrical parameters was developed. One thousand individual geometries comprising four degeneration grades were generated in a probabilistic manner, and the influence of the severity of disc degeneration on the mechanical response of the motion segment to different loading conditions was statistically evaluated.

RESULTS

Variations in the individual structural parameters resulted in marked variations in all evaluated parameters within each degeneration grade. Nevertheless, the effect of degeneration in almost all evaluated response values was statistically significant. With degeneration, the intradiscal pressure progressively decreased. At the same time, the facet loads increased and the ligament tension was reduced. The initially nonlinear load-deformation relationships became linear whereas the segment stiffness increased.

CONCLUSIONS

Results indicate significant stiffening of the motion segment with progressing degeneration and gradually increasing loading of the facets from nondegenerated to moderately degenerated conditions along with a significant reduction of the ligament tension in flexion.

Asymmetry between the superior and inferior endplates is a risk factor for lumbar disc degeneration

Endplate pathology plays an important role in the development of lumbar disc degeneration. Previous research paid little attention to differences between the superior and inferior endplates as a possible risk factor for disc degeneration. The purpose of this study was to test the hypothesis that asymmetry between the superior and inferior endplates is a risk factor for the development of lumbar disc degeneration. A total of 134 patients with lumbar disc herniation (LDH) and 100 healthy adults ("Controls") underwent magnetic resonance imaging scans. Each disc was categorized as non-degenerated (Pfirrmann grades I-II) or degenerated (Pfirrmann grades III-V) and get the following three groups: "Degenerated LDH" discs (n = 145), "Non-degenerated LDH" discs (n = 525) and "Non-degenerated Control" discs (n = 500). On mid-sagittal image, the lumbar endplate morphology could be categorized into three types: Flat, concave, and irregular. Superior and inferior endplates of a given disc were "symmetric" if both were of the same type, and "asymmetric" if they were of different types. The proportion of asymmetric endplates at L4-5 was higher in the "Degenerated LDH" discs group (47%) than in the "Non-degenerated LDH" discs group (21%) or "Non-degenerated Control" discs group (7%) (p < 0.05). At L5-S1 the proportions were 73%, 55%, and 38% (p < 0.05). Asymmetry of superior and inferior endplates in the mid-sagittal plane is a risk factor for lumbar disc degeneration. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:2469-2475, 2018.

Leptin regulates disc cartilage endplate degeneration and ossification through activation of the MAPK-ERK signalling pathway in vivo and in vitro

Recent findings demonstrate that leptin plays a significant role in chondrocyte and osteoblast differentiation. However, the mechanisms by which leptin acts on cartilage endplate (CEP) cells to give rise to calcification are still unclear. The aim of this study was to evaluate the effects of leptin that induced mineralization of CEP cells in vitro and in vivo. We constructed a rat model of lumbar disc degeneration and determined that leptin was highly expressed in the presence of CEP calcification. Rat CEP cells treated with or without leptin were used for in vitro analysis using RT‐PCR and Western blotting to examine the expression of osteocalcin (OCN) and runt‐related transcription factor 2 (Runx2). Both OCN and Runx2 expression levels were significantly increased in a dose‐ and time‐dependent manner. Leptin activated ERK1/2 and STAT3 phosphorylation in a time‐dependent manner. Inhibition of phosphorylated ERK1/2 using targeted siRNA suppressed leptin‐induced OCN and Runx2 expression and blocked the formation of mineralized nodules in CEP cells. We further demonstrated that exogenous leptin induced matrix mineralization of CEP cells in vivo. We suggest that leptin promotes the osteoblastic differentiation of CEP cells via the MAPK/ERK signal transduction pathway and may be used to investigate the mechanisms of disc degeneration.

Structural vertebral endplate nomenclature and etiology: a study by the ISSLS Spinal Phenotype Focus Group

PURPOSE

Vertebral endplate abnormalities may be associated with disc degeneration and, perhaps, pain generation. However, consensus definitions for endplate findings on spine MRI do not exist, posing a challenge to compare findings between studies and ethnic groups. The following survey was created to characterize the variability among the global spine community regarding endplate structural findings with respect to nomenclature and etiology.

METHODS

A working group within the International Society for the Study of the Lumbar Spine (ISSLS) Spinal Phenotype Focus Group was established to assess the endplate phenotype. A survey which consisted of 13 T2-weighted sagittal MRIs of the human lumbar spine illustrating the superior and inferior endplates was constructed based on discussion and agreement by the working group. A list of nomenclature and etiological terms with historical precedence was generated. Participants were asked to describe the endplates of each image and select from 14 possible nomenclatures and 10 etiological terms along with the option of free text response. The survey was entered into RedCap and was circulated throughout the ISSLS membership for data capture. Participants' demographics were also noted.

RESULTS

The survey was completed by 55 participants (87% males; 85% above 45 years of age, 39 clinicians, and 16 researchers). Sixty-eight percent of researchers and seventy-four percent of clinicians reported more than 16 and 20 years of research and clinical experience. Considerable variation existed in selection of nomenclature, etiology, and degree of severity of the endplate structural findings (reliability coefficients for single measures in each case were 0.3, 0.08, and 0.2, respectively). Sixty-seven percent regarded Modic changes as being a structural endplate finding. Approximately 84 and 80% of clinicians and researchers, respectively, agreed that a standardized endplate nomenclature and understanding the etiology is clinically important and needed.

CONCLUSIONS

This study found that variations exist with respect to endplate nomenclature and etiology between clinicians and basic scientists, and paves the way for a consensus process to formalize the definitions.

Association between Measures of Vertebral Endplate Morphology and Lumbar Intervertebral Disc Degeneration

Purpose

The aim of this study was to evaluate the association between vertebral endplate morphology and the degree of lumbar intervertebral disc degeneration via magnetic resonance imaging (MRI).

Methods

In total, 150 patients who met the inclusion criteria and were 20–60 years of age were retrospectively evaluated. Patients were evaluated for the presence of intervertebral disc degeneration or herniation, and the degree of degeneration was assessed at all lumbar levels. Vertebral endplate morphology was evaluated based on the endplate sagittal diameter, endplate sagittal concave angle (ECA), and endplate sagittal concave depth (ECD) on sagittal MRI. The association between intervertebral disc degeneration or herniation and endplate morphological measurements was analysed.

Results

In MRI, superior endplates ( ie, inferior endplates of the superior vertebra) were concave and inferior endplates ( ie, superior endplates of the inferior vertebra) were flat at all disc levels. A decrease in ECD and an increase in ECA were detected at all lumbar levels as disc degeneration increased ( P < .05). At the L4-L5 and L5-S1 levels, a decrease in ECD and an increase in ECA were detected in the group with herniated lumbar discs ( P < .05). There was no association between lumbar disc degeneration or herniation and endplate sagittal diameter at lumbar intervertebral levels ( P > .05). At all levels, ECD of women was significantly lesser than that of men and ECA of women was significantly greater than that of men ( P < .05).

Conclusions

There is an association between vertebral endplate morphology and lumbar intervertebral disc degeneration. Vertebral endplates at the degenerated disc level become flat; the severity of this flattening is correlated with the degree of disc degeneration.

Validity and reliability of computerized measurement of lumbar intervertebral disc height and volume from magnetic resonance images

BACKGROUND CONTEXT

Magnetic resonance (MR) examinations of morphologic characteristics of intervertebral discs (IVDs) have been used extensively for biomechanical studies and clinical investigations of the lumbar spine. Traditionally, the morphologic measurements have been performed using time- and expertise-intensive manual segmentation techniques not well suited for analyses of large-scale studies..

PURPOSE

The purpose of this study is to introduce and validate a semiautomated method for measuring IVD height and mean sagittal area (and volume) from MR images to determine if it can replace the manual assessment and enable analyses of large MR cohorts.

STUDY DESIGN/SETTING

This study compares semiautomated and manual measurements and assesses their reliability and agreement using data from repeated MR examinations.

METHODS

Seven healthy asymptomatic males underwent 1.5-T MR examinations of the lumbar spine involving sagittal T2-weighted fast spin-echo images obtained at baseline, pre-exercise, and postexercise conditions. Measures of the mean height and the mean sagittal area of lumbar IVDs (L1-L2 to L4-L5) were compared for two segmentation approaches: a conventional manual method (10-15 minutes to process one IVD) and a specifically developed semiautomated method (requiring only a few mouse clicks to process each subject).

RESULTS

Both methods showed strong test-retest reproducibility evaluated on baseline and pre-exercise examinations with strong intraclass correlations for the semiautomated and manual methods for mean IVD height (intraclass correlation coefficient [ICC]=0.99, 0.98) and mean IVD area (ICC=0.98, 0.99), respectively. A bias (average deviation) of 0.38 mm (4.1%, 95% confidence interval 0.18-0.59 mm) was observed between the manual and semiautomated methods for the IVD height, whereas there was no statistically significant difference for the mean IVD area (0.1%±3.5%). The semiautomated and manual methods both detected significant exercise-induced changes in IVD height (0.20 and 0.28 mm) and mean IVD area (5.7 and 8.3 mm(2)), respectively.

CONCLUSIONS

The presented semiautomated method provides an alternative to time- and expertise-intensive manual procedures for analysis of larger, cross-sectional, interventional, and longitudinal MR studies for morphometric analyses of lumbar IVDs.

Is post-contrast MRI a valuable method for the study of the nutrition of the intervertebral disc?

Decreased nutrition has been proposed as a potential mechanism leading to intervertebral disc degeneration. A method to investigate it in vivo is the MRI evaluation of the transport of a paramagnetic contrast agent, which is assumed to diffuse through the endplate to the disc using the same mechanisms as the cell nutrients. However, previous numerical studies questioned the value of this method as a model to investigate disc nutrition. To assess its validity, a parametric osmoporoelastic finite element model of a lumbar intervertebral disc incorporating diffusion and convection of a solute (representing the contrast agent) was developed. A Taguchi sensitivity analysis was performed in order to assess the relevance of various parameters which influence the solute transport. Subsequently, a full-factorial sensitivity analysis was used to investigate specifically the diffusion coefficients of the contrast agent. The most important parameters in determining the results were the disc height, the diffusion coefficients and the pharmacokinetic of the contrast agent. However, diffusion coefficients values as measured in in vitro studies would lead to insubstantial enhancement of the MRI signal. Thus, transport mechanisms other than pure diffusion should be active in in vivo transport of the contrast agent. In conclusion, the study showed that post-contrast MRI may not be suited for a quantitative analysis, but only for a qualitative examination aimed for example to detect endplate lesions. Open questions remain on the use of post-contrast MRI for the investigation of the relevance of reduced nutrition as a trigger to disc degeneration.

Sagittal geometry of the middle and lower cervical endplates

PURPOSE

Construct subsidence is a relatively common complication following anterior cervical fusion. Its occurrence has been revealed to be closely related to endplate-implant contact interface. But current literature focusing on the anatomy of cervical endplate is very scarce. The purpose of this morphometric study was to analyse the sagittal geometry, especially the concavity and slope, of vertebral endplates from C3 to C7 by employing data from CT scans.

METHODS

Reformatted CT scans of 97 individuals were analyzed and endplate concavity depth, endplate concavity apex location, as well as endplate slope were measured in midsagittal plane. Those specific parameters were compared among different age and gender groups. Meanwhile, comparison between superior and inferior endplate of each vertebra was also performed.

RESULTS

Age and gender did not influence endplate concavity depth, endplate concavity apex location, or endplate slope significantly (P > 0.05). Endplate concavity depths of superior endplates (range 0.9-1.2 mm) were significantly smaller than those of inferior endplates (range 2.1-2.7 mm). Endplate concavity apex was always located in the posterior half of the endplate, with the superior one ranged from 56 to 67% and the inferior one 52 to 57%. Average endplate slopes of superior endplates were between 4.5° and 9.0°, and average inferior endplate slopes ranged from 4.5° to 7.5°. Among all measured segments, C5 had the largest endplate slope values, while C7 the least.

CONCLUSIONS

Superior endplate is more flat than its inferior counterpart in middle and lower cervical spine, and the concavity apex is always located in the posterior half of the endplate. Endplate slope is correlated with cervical curvature, greater slope implying more significant lordosis. These sagittal endplate geometrical parameters should be taken into consideration when investigating implant subsidence following anterior cervical fusion.

Association of endothelin-1 expression and cartilaginous endplate degeneration in humans

Background

Inflammatory cytokines are involved in intervertebral disc (IVD) degeneration. Endothelin-1 (ET-1), a 21-amino-acid cytokine implicated with cartilage degradation, is secreted by vascular endothelial cells and also by many other cell types. The expression of ET-1 in human IVD cartilage endplate (CEP) and its role in disc degeneration have not been explored.

Methods and Findings

The expression of ET-1 in degenerated CEP was analyzed by immunohistochemical staining and Western blotting; ET-1 was demonstrated in cartilaginous endplate cells (CECs) by immunofluorescent staining. The ET-1 mRNA expression and protein production by CECs stimulated by tumor necrosis factor alpha (TNF-α), a pro-inflammatory cytokine, were determined by real-time PCR analysis and Western blotting, respectively. The matrix metalloprotease-1 (MMP-1), MMP-13 and tissue inhibitor of metalloproteases-1 (TIMP-1) levels in the supernatant of cultured CECs treated with ET-1 were determined using enzyme-linked immunosorbent assays. Nitric oxide (NO) release and nitric oxide synthase (NOS) activity were measured using a spectrophotometric assay. The apoptosis of CECs by ET-1 was measured by an Annexin V-FITC detection assay. The production of ET-1 in degenerated cartilage endplate was significantly higher than normal CEP. The results showed that ET-1 was expressed by CECs and modulated by TNF-α in a dose-dependent manner. ET-1 increased production of MMP-1 and MMP-13, decreased TIMP-1 production, and induced NO and NOS release by cultured CECs. The direct stimulation of CECs by ET-1 did not promote cell apoptosis.

Conclusion

The study results suggest that ET-1 played a pivotal role in human CEP degeneration, and may be a new target for development of therapies for this condition.