- Research article
- Open Access
Biomechanical role of osteoporosis affects the incidence of adjacent segment disease after percutaneous transforaminal endoscopic discectomy
Journal of Orthopaedic Surgery and Researchvolume 14, Article number: 131 (2019)
Variation in the biomechanical characteristics of intervertebral discs adjacent to the segment disc after undergoing percutaneous transforaminal endoscopic discectomy (PTED) in models with normal and abnormal bone mineral density (BMD) was estimated using the finite element method.
The study investigated the change in the incidence of adjacent segment disease (ASD) after PTED in patients without and with osteoporosis.
PTED has been widely used for treating lumbar disc herniation (LDH); changes in BMD will affect biomechanical characteristics, possibly leading to changes in the incidence of ASD after PTED. However, this issue remains largely unclear.
A non-linear, lumbosacral finite element model was reconstructed based on imaging data and validated using compared values computed by the current model from published and well-validated, in vitro biomechanical experiment studies. Corresponding PTED models with normal and abnormal BMDs were also reconstructed. Shear and von Mises stresses on the annulus fibrosis, the von Mises stress on the endplates in L5–S1 segment discs, and the total deformation of current lumbosacral models were computed in different body positions by changing loading conditions, including flexion, extension, left and right lateral bending, and axial rotation.
In most loading conditions, biomechanical characteristics of the lumbosacral segment discs with normal BMDs after PTED slightly increased. However, in the PTED model with osteoporosis, most of the biomechanical characteristics dramatically increased.
Osteoporosis leads to the deterioration of biomechanical characteristics in the adjacent segment disc after PTED; this variation may also result in an increase in the incidence of ASD. However, further studies on the interactions between pathological changes are warranted.
PTED can lead to a slight deterioration in the biomechanical characteristics of adjacent segment discs, which are significantly intensified in an osteoporosis model.
Because women are more likely to suffer from osteoporosis, intensified postoperative biomechanical deterioration in patients with osteoporosis can be seen as a credible reason for the higher risk of ASD after PTED among elderly women.
The interaction of risk factors of ASD was not estimated in the current study, warranting further research.
ASD is a series of postoperative complications, including types of symptoms as low back pain (LBP) and low limb pain or numbness [1, 2]. The etiological factors and demographic characteristics of ASD have been widely discussed. The risk factors concerning demographic data, such as senility, female sex, and obesity, as well as etiological triggers such as severe stress concentration, hypermotility, and preexisting degenerative changes in adjacent segment discs have been mentioned [1, 3,4,5,6].
Biomechanical deterioration and resulting injury to spine structures can be seen as the most significant triggers for the development of disc degeneration change and which is one the most important reasons of ASD [1,2,3, 7,8,9,10]. While discussing the risk factors for ASD, demographic characteristics are always assumed to be defined by some biomechanical pathogenesis. For instance, overweight patients suffer from excessive pressure in the lumbar spine, which may accelerate disc degeneration. This can explain why patients with higher body mass indices (BMI) are more prone to ASD [1, 3, 11]. And multiple studies have revealed that the incidence of ASD is high in senile patients and that preexisting disc degeneration changes is a vital risk factor of ASD [1,2,3]. Considering that disc degeneration becomes more pronounced with increasing age, we believe that senile patients always suffer from ASD because of preexisting degenerations [5, 6].
A part of patients suffer from lumbar disc herniation (LDH) accompanied with osteoporosis, and women, especially elderly women, are more likely to suffer from osteoporosis . As mentioned earlier, women are more vulnerable to ASDs [3, 5]; therefore, logically, osteoporosis should be implicated in the occurrence of ASDs. Nevertheless, the clinical review by Wang  indicated that the risk of ASD can be altered by changing the BMD and that the most important indicator of osteoporosis was not statistically different in patients with or without ASD. However, a review by Park et al.  highlighted that osteoporosis can be deemed as an independent predictor of ASD.
Meanwhile, some published literature revealed that changes in BMD exert a vital impact on the biomechanical characteristics of the spine as well as affect the incidence of postoperative complications [13,14,15]. For biomechanical deterioration is important in the development of ASD, exploration of biomechanical mechanism is of great significance for the investigation of the pathogenesis of ASD, and different surgical methods have been computed by finite element methods (FEM) to estimate if which could decrease the risk of ASD or not [16,17,18]. However, no published studies illustrated this issue after PTED, although it became increasingly popular for treating LDH. In consideration of this, the ASD will lead to a poor prognosis for patients . Investigations concerning risks of ASD in PTED patients are essential.
We hypothesized that this change affects the incidence of ASD, and no statistically significant conclusion from some clinical report may origin from relatively small sample size or short follow-up period. In addition, to the best of our knowledge, no published English biomechanical study has focused on whether osteoporosis affects the incidence of ASD after PTED. To test this hypothesis and to clarify the abovementioned issue, a lumbosacral three-dimensional model was reconstructed and validated, and corresponding PTED models with and without osteoporosis were reconstructed to determine whether changes in BMD influence shear and the von Mises stresses on the annulus fibrosis, the von Mises stress on the endplates, and the total deformation of the lumbosacral spine. All the characteristics were assessed in the L5–S1 segment, except for the total deformation of spine model, owing to the high incidence of disc degenerative diseases in this segment .
Material and methods
Reconstruction of a complete spine model
A 3D lumbosacral model was reconstructed based on the imaging data obtained from an adult male volunteer (L3–S1, three segments of the intervertebral disc (IVD), six facet joints, and six ligaments) without any history of lumbar diseases using computed tomography (CT). The structures in the current model could not be clearly differentiated using CT and have therefore been constructed based on our anatomical observations [13, 19]. In the current study, the vertebra included a cortical shell (thickness 0.8 mm), a trabecular core, superior and inferior endplates (thickness: 0.8 mm), and posterior structures. IVD comprised the annulus fibrosis and inner nucleus, which occupied 44% of the cross-sectional area of the disc, whose location was slightly backward to the disc center, and comprised facet joints including capsule and facet cartilage surfaces (0.25 mm) and ligaments including the anterior longitudinal ligament, posterior longitudinal ligament, ligamentum flavum, intertransverse ligament, interspinous ligament, and supraspinous ligament [18, 20, 21].
Reconstruction of PTED and osteoporosis models
The L4–L5 segment, which is one of the most common sites of LDH, was selected for the simulation of PTED. PTED was simulated on the right side. To imitate tears in the annulus fibrosis, a 5-mm incision was made and to represent discectomy, the nucleus in this segment was deleted [22, 23]. Moreover, facetectomy and ligamentum flavum excision are inevitable in PTED, and therefore, sufficient space was created to insert instruments and remove the nucleus. In this, one third of each of these structures were removed [21, 22, 24, 25]. The schematic diagram of current models is shown in Fig. 1, and the magnetic resonance imaging (MRI) date of patients subjected to PTED is shown in Fig. 2.
The reconstruction of osteoporosis FEM models has been described in the scientific literature. The Young modulus in different parts of the current model has been changed, but morphological features remained the same [14, 15, 17]. The Young modulus for cortical, endplate, and posterior structures decreased to 67% of its original value. The Young modulus in cancellous bone decreased to 34% [13,14,15].
Boundary and loading conditions
Boundary and loading conditions were uniformly applied in the three models in this study. To simulate the complicated spinal motion, five different loading conditions, including flexion, extension, left and right bending, and axial torsion, were applied by setting a 10-Nm moment on the superior surface of the lumbosacral model and a 800-N vertical compressive force in the basic loading condition [9, 17].
To reduce the computed error, tetrahedral elements were selected to fill the current model because they seemed suitable for models with complex surfaces; moreover, mesh refinement was set in distortion areas. In current models, the contact type was set as “bounded”; however, for facet surfaces, it was set as “frictionless” [9, 17]. Regarding the material properties, the vertebra and facet cartilages were defined as isotropic, homogeneous materials; the Mooney–Rivlin hyper-elastic material was set in the annulus fibrosis; and the nucleus was defined as an incompressible semifluid material. Moreover, “tension-only” cable elements were used for defining the ligaments (Tables 1 and 2) [13, 18]. In addition, for simplification, the complete lumbosacral model was defined as model 1, the model that underwent PTED with normal BMD as model 2, and the model with osteoporosis as model 3.
The intradiscal pressure in the L4-L5 segment was validated by comparing with a previously published, well-validated in vitro biomechanical study, and the disc compression values were validated by comparing with another in vitro research. To calculate the intradiscal pressure, 300, 1000, and 2000 N of vertical compressive pressure was applied on the superior surface of model 1 and 1200 N of vertical compressive pressure was applied to assess disc compression [22, 25].
At 300, 1000, and 2000 N of vertical compressive force, the intradiscal pressure in our model was 0.2, 0.78, and 1.49 Mpa and that previously reportedly was 0.3 ± 0.09, 0.9 ± 0.26, and 1.85 ± 0.46 Mpa, respectively . Disc compression values in the L3–S1 segment in the current model were 1.7, 1.5, and 1.1 mm and those from a well-validated research were 1.6 ± 0.55, 1.6 ± 0.5, and 1.3 ± 0.5 mm in L3–L4, L4–L5, and L5–S1 segment disc, respectively . Thus, the values in the present study were within one standard deviation from those reported in the well-validated studies except the intradiscal pressure under 300-N compressive pressure (Fig. 3). Considering the only exception listed immediately below the standard deviation and equaling to the lower value in the original data, we believe that the current model can make a credible representation of the lumbosacral spine in actual body positions .
Biomechanical characteristics variation
Notably, the trend in biomechanical characteristics increased in model 2 in most loading conditions, except for the von Mises stress in the annulus in extension condition and when the endplates were in left lateral bending. Meanwhile, except for the shear stress of the annulus and the von Mises stress in the endplates in the flexion condition, most of the biomechanical characteristics increased by < 5%.
However, this trend dramatically changed with respect to the differences between model 3 and the other two models. Compared with that in model 2, the von Mises stress in endplates was lower in flexion and left lateral bending condition in model 3, and this value in the annulus in extension condition increased in model 3 but was lower than that in model 1. In addition, all the other biomechanical characteristics increased dramatically in model 3 as compared with those in models 1 and 2 (Fig. 4).
The risk factors of ASD have been widely investigated, and the demographic results from published literature indicate that old age, high BMI, and female sex increase the risk of incidence of ASD [1, 3, 4, 11]. Meanwhile, hypermotility of the lumbar spine, abnormal stress concentration, and preexisting degenerative changes in the adjacent discs can also trigger ASD [5, 6].
The aberrant distribution of stress and the resulting injury of disc structures have been regarded as the principal factors of disc degeneration and have been proven to be the primary risk factors for ASD [3, 7, 8]. Therefore, we believe that the high risk of the incidence of ASD in a specific population may be connected with a biomechanical cause. Moreover, our hypothesis seems reliable considering that a high incidence of ASD in senile patients may root from multiple segment degenerations in lumbar discs before surgery and that in obese patients may be connected to the greater load on the spine .
However, it is unclear as to why women are more vulnerable to ASD, considering that the biomechanical factors that are common in women and are related to ASD have not been clearly estimated. Published studies have illustrated that higher expression of estrogen receptor in the facet cartilage is connected to degenerative spondylolisthesis (a typical manifestation of lumbar instability) [27, 28] and is an important reason for ASD [1, 8, 23]. This seems to be a reasonable explanation for the relatively higher risk of ASD in women. However, considering the decisive role of biomechanical deterioration in the pathological changes of ASD [1, 3], an investigation of the relationship between biomechanical factors common to females and the pathogenesis of ASD is still necessary.
Women, particularly elderly women, are susceptible to osteoporosis. As per the epidemiological investigation report by Ha , women are approximately three times more prone to osteoporosis than men. Considering that Chinese and Korean studies focus on Asian populations, this data can be used as a reference in studies including Chinese population, indicating that the relationship between the higher risk of ASD in women may be connected with a higher incidence of osteoporosis. Published studies have reported that osteoporosis may increase the risk of disc degeneration in the lumbar spine by adversely changing the loading pattern and intradiscal pressure [29, 30]. We therefore hypothesize that osteoporosis is a key reason for the relatively higher risk of ASD in women.
LDH is one of the most common disc degenerative diseases that is widely treated by surgery, and investigation record-related risk factors of ASD have always focused on the fusion methods rather than on PTED although it has been widely used in the past decade [2, 31]. For the abovementioned reasons, we believe study concerning the biomechanical influence of the alteration of BMD for the ASD risk for PTED patients is clinically important. Moreover, to investigate this issue, a lumbosacral finite element model was reconstructed to determine the change in shear and the von Mises stresses on the annulus fibrosis, the von Mises stress on the endplates, and the total deformation of models with and without osteoporosis undergoing PTED.
Stress concentration typically results from abnormal stress distribution, which eventually leads to the formation of injury structures in IVDs . The endplates play an irreplaceable role during the process of loading conditions; therefore, damage to the endplates can result in various issues. For instance, stress concentration on the trabecular bone may result in microfractures, disrupting the nutrient pathway of the IVDs. Considering that IVD is an avascular structure, the disruption of nutrition may severely accelerate the process of disc degeneration . The von Mises stress in endplates increased in most loading conditions in the osteoporosis model, which may increase the incidence of ASD associated with changes in the endplates.
Abnormal stress concentration on the annulus can result in tears. This trend is connected with the generation of ASD [3, 7]. For example, the von Mises stress concentration is connected with the circumferential tears and increased shear stress is strongly linked to radical tears that can lead to discogenic LBP and increase the risk of LDH in adjacent segment discs, which are two typical manifestations of ASD [9, 21, 33]. Moreover, endplate injury intensifies abnormal stress distribution in the annulus and aggravates the incurred damage [7, 34]. As presented in our FEM results, the trend of stress concentration can be observed in the PTED model of osteoporosis. As mentioned in the theories presented above, the risk of annulus tears and subsequent ASD may increase in osteoporosis patients (Figs. 5 and 6).
Considering that the deformations of the spine mainly occur in IVDs rather than in vertebral bodies, we believe that the increase in total deformation in the current model is a good representation of the hypermotility of the adjacent IVDs . The hypermotility of IVDs has been proven to trigger LBP, a typical symptom of ASD [31, 35]. Moreover, a change in this index is an important and independent predictor of ASD, which could accelerate the process of disc degeneration . Therefore, increase in the total deformation in the osteoporosis model may predict an increasing trend of the incidence of ASD in patients with osteoporosis after PTED.
In a summary, it can be concluded that osteoporosis may be an independent risk factor of ASD for patients undergoing PTED, which may be the cause of the high incidence of ASD in women.
Nevertheless, some limitations in the current study still exist. For instance, several pathological changes, such as osteoporosis, osteoarthritis of facet cartilage, and multiple-segment degeneration, can be observed in PTED patients (especially elderly patients). The interaction of these pathological changes may strongly impact the risk of ASD and affect the credibility of our conclusion [9, 15, 32]. These issues require further investigation. Hence, the existing conclusion may be inaccurate or only applicable to patients without degenerative changes in other segments. Besides, the model constructed in this study omitted soft tissues around the spine. As a result, the current FEM data pertain only to ideal cases. All of these factors should be explored in future studies.
Osteoporosis may be regarded as an independence risk factor of ASD after PTED, and this association can explain the phenomenon of women being more likely to suffer from ASD.
Adjacent segment disease
Bone mineral density
Body mass indices
Finite element method
Low back pain
Lumbar disc herniation
Magnetic resonance imaging
Percutaneous transforaminal endoscopic discectomy
Wang H, Ma L, Yang D, Wang T, Liu S, Yang S, Ding WJM. Incidence and risk factors of adjacent segment disease following posterior decompression and instrumented fusion for degenerative lumbar disorders. 2017;96(5):e6032.
Liang J, Dong Y, Zhao HJJoOS, Research. Risk factors for predicting symptomatic adjacent segment degeneration requiring surgery in patients after posterior lumbar fusion. J Orthop Surg Res. 2014;9(1):97.
Park P, Garton HJ, Gala VC, Hoff JT, Mcgillicuddy JE. Adjacent segment disease after lumbar or lumbosacral fusion: review of the literature. Spine. 2004;29(17):1938–44.
Rahm MD, Hall BB. Adjacent-segment degeneration after lumbar fusion with instrumentation: a retrospective study. J Spinal Disord. 1996;9(5):392–400.
Ghasemi AA. Adjacent segment degeneration after posterior lumbar fusion: An analysis of possible risk factors. Clin Neurol Neurosurg. 2016;143(125):15–8.
Kim HJ, Kang KT, Chun HJ, Lee CK, Chang BS, Yeom JS. The influence of intrinsic disc degeneration of the adjacent segments on its stress distribution after one-level lumbar fusion. Eur Spine J. 2015;24(4):1–11.
Adams MA, Roughley PJJS. What is intervertebral disc degeneration, and what causes it? 2006;31(18):2151–61.
Adams MA, Freeman BJ, Morrison HP, Nelson IW, Dolan P. Mechanical initiation of intervertebral disc degeneration. Spine (Phila Pa 1976). 2000;25(13):1625–36.
Ruberté LM, Natarajan RN, Andersson GBJ. Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments—a finite element model study. J Biomech. 2009;42(3):341–8.
Wen-Hsien C, Shang-Chih L, Shih-Hao C, Chih-Wei W, Wen-Chi T, Yeung-Jen C, HJS J-R. Biomechanical effects of disc degeneration and hybrid fixation on the transition and adjacent lumbar segments: trade-off between junctional problem, motion preservation, and load protection. Spine. 2012;44(1):29–39.
Chien-Yu O, Tao-Chen L, Tsung-Han L, Yu-Hua HJN. Impact of body mass index on adjacent segment disease after lumbar fusion for degenerative spine disease. 2015;76(4):396–401.
Ha YC. Epidemiology of osteoporosis in Korea. J Korean Med Assoc. 2016;59(11):836–41.
Yang P, Zhang Y, Ding HW, Liu J, Ye LQ, Xiao J, Tu Q, Yang T, Wang F, Sun GG. Pedicle screw fixation with kyphoplasty decreases the fracture risk of the treated and adjacent non-treated vertebral bodies: a finite element analysis. J Huazhong Univ Sci Technolog Med Sci. 2016;36(6):887–94.
Ferguson SJ, Polikeit A, Nolte LP. The effect of cement augmentation on the load transfer in an osteoporotic functional spinal unit: finite-element analysis. Spine. 2003;28(10):991.
Zhang L, Yang G, Wu L, Yu B. The biomechanical effects of osteoporosis vertebral augmentation with cancellous bone granules or bone cement on treated and adjacent non-treated vertebral bodies: a finite element evaluation. Clin Biomech (Bristol, Avon). 2010;25(2):166–72.
Más Y, Gracia L, Ibarz E, Gabarre S, Peña D, Herrera A: Finite element simulation and clinical follow-up of lumbar spine biomechanics with dynamic fixations. PLoS One 2017, 12(11):e0188328.
Shujie T, Rebholz BJ. Does lumbar microdiscectomy affect adjacent segmental disc degeneration? A finite element study. J Surg Res. 2013;182(1):62–7.
Wen-Hsien C, Yi-Jie K, Shang-Chih L, Chih-Wei W, Shih-Hao C, Yeung-Jen C, Jiun-Ren HJS. Comparison among load-, ROM-, and displacement-controlled methods used in the lumbosacral nonlinear finite-element analysis. Spine. 2013;38(5):276–85.
Yee TJ, Terman SW, Marca FL, Park P. Comparison of adjacent segment disease after minimally invasive or open transforaminal lumbar interbody fusion. J Clin Neurosci. 2014;21(10):1796–801.
Gu G, Wang C, Gu X, Zhang H, Zhao Y, He S. Percutaneous transforaminal endoscopic discectomy for adjacent segment disease after lumbar fusion in elderly patients over 65 years old. World Neurosurg. 2018:112, e830–e836.
Sengupta DK. Clinical biomechanics of the spine. Spine. 2017;42(Suppl 7):S3.
Brinckmann P, Grootenboer H. Change of disc height, radial disc bulge, and intradiscal pressure from discectomy. An in vitro investigation on human lumbar discs. Spine. 1991;16(6):641.
Xiao L, Ni C, Shi J, Wang Z, Wang S, Zhang J, Lu A, JMSMIMJoE, Research C. Analysis of correlation between vertebral endplate change and lumbar disc degeneration. Med Sci Monit. 2017;23:4932–8.
Thomas H, Karolien BD, Michael S, Boris MJS. Endoscopic transforaminal discectomy for recurrent lumbar disc herniation: a prospective, cohort evaluation of 262 consecutive cases. Spine. 2008;33(9):973.
Renner SM, Natarajan RN, Patwardhan AG, Havey RM, Voronov LI, Guo BY, Andersson GBJ, An HS. Novel model to analyze the effect of a large compressive follower pre-load on range of motions in a lumbar spine. J Biomech. 2007;40(6):1326–32.
Luoma K, Vehmas T, Kerttula L, Grönblad M, Rinne E, JEsjopotESS, the European Spinal Deformity Society, Society tESotCSR. Chronic low back pain in relation to Modic changes, bony endplate lesions, and disc degeneration in a prospective MRI study. Eur Spine J. 2016;25(9):2873.
Zhu R, Niu WX, Zeng ZL, Tong JH, Zhen ZW, Zhou S, Yu Y, Cheng LM. The effects of muscle weakness on degenerative spondylolisthesis: a finite element study. Clin Biomech. 2017;41(1):34–8.
Kee-Yong H, Cheong-Ho C, Ki-Won K, Youn-Soo K, Kie-Ho N, Jun-Seok LJS. Expression of estrogen receptor of the facet joints in degenerative spondylolisthesis. Spine. 2005;30(5):562–6.
Zhou Z, Tian FM, Wang P, Gou Y, Zhang H, Song HP, Wang WY, Zhang LJS. Alendronate prevents intervertebral disc degeneration adjacent to a lumbar fusion in ovariectomized rats. Spine. 2015;40(20):1.
Alexander T, Stylianos Orestis S, Kleovoulos A, Nikolaos M, Sashin A, JSJOJotNASS. Assessment of stress patterns on a spinal motion segment in healthy versus osteoporotic bony models with or without disc degeneration: a finite element analysis. Spine J. 2015;15(3):S17–22.
Lee CK, Langrana NA. Lumbosacral spinal fusion. A biomechanical study. Spine. 1984;9(6):574.
Park WM, Kim K, Kim YH. Effects of degenerated intervertebral discs on intersegmental rotations, intradiscal pressures, and facet joint forces of the whole lumbar spine: Pergamon Press, Inc; 2013.
Clancy C, Quinn A, Wilson F. The aetiologies of failed back surgery syndrome: a systematic review. J Back Musculoskelet Rehabil. 2016;30(3):1–8.
Mcnally DS, Shackleford IM, Goodship AE, Mulholland RCJS. In vivo stress measurement can predict pain on discography. Spine. 1996;21(22):2580–7.
Bastian L, Lange U, Knop C, Tusch G, Blauth M. Evaluation of the mobility of adjacent segments after posterior thoracolumbar fixation: a biomechanical study. Eur Spine J. 2001;10(4):295.
We acknowledge Professor Tao Wu for the guidance of the model reconstruction and Dr. Wuguang Lu for the guidance of the manuscript preparation.
The current research was funded by the Foundation for leading talent in traditional Chinese medicine of Jiangsu province (2018 SLJ0210), Natural Science Foundation of Jiangsu (BK20151604), Project of Jiangsu province health and family planning commission (2016 BJ16026), and Project of Jiangsu province health and family planning commission (H2018025)
Availability of data and materials
All the data of the manuscript are presented.
Ethics approval and consent to participate
Ethical approval for this research was obtained by the Ethics Committees of Jiangsu ProvinceHospital on Integration of Chinese and Western Medicine.
The CT data for the reconstruction of FEM model was from a male volunteer, who knows all of our experiment, and informed consent was signed before the CT scan of his lumbar spine.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
- Adjacent segment diseases
- Percutaneous transforaminal endoscopic discectomy
- Bone mineral density
- Finite element method