Next Article in Journal
Predicting Adherence to Walking from Anxiety, Depression, Disease Impact, Catastrophizing, and Cognitive Fusion in Patients with Fibromyalgia: Does Pain Severity Matter?
Previous Article in Journal
Managing Oral Health in the Context of Antimicrobial Resistance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model

1
Department of Orthopedic Surgery, Tazawako Hospital, Akita 014-1201, Japan
2
Department of Orthopedic Surgery, Graduate School of Medicine, Akita University, Akita 010-8543, Japan
3
Department of Rehabilitation Medicine, Akita University Hospital, Akita 010-8543, Japan
4
Department of Mechanical Engineering, Faculty of Engineering Science, Akita University, Akita 010-8502, Japan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(24), 16452; https://doi.org/10.3390/ijerph192416452
Submission received: 12 November 2022 / Revised: 3 December 2022 / Accepted: 6 December 2022 / Published: 8 December 2022
(This article belongs to the Section Aging)

Abstract

:
The intervertebral disc loading based on compensated standing posture in patients with adult spinal deformity remains unclear. We analyzed the relationship between sagittal alignment and disc compression force ( F m ). In 14 elderly women, the alignment of the sagittal spinopelvic and lower extremities was measured. F m was calculated using the Anybody Modeling System. Patients were divided into low sagittal vertical axis (SVA) and high SVA groups. Comparisons between the two groups were performed and the relationship between the F m and each parameter was examined using Spearman’s correlation coefficient (r). The mean lumbar F m in the high SVA group was 67.6%; significantly higher than that in the low SVA group (p = 0.046). There was a negative correlation between cervical F m with T1 slope (r = −0.589, p = 0.034) and lumbar F m with lumbar lordosis (r = −0.566, p = 0.035). Lumbar F m was positively correlated with center of gravity-SVA (r = 0.615, p = 0.029), T1 slope (r = 0.613, p = 0.026), and SVA (r = 0.612, p = 0.020). The results suggested sagittal malalignment increased the load on the thoracolumbar and lower lumbar discs and was associated with cervical disc loading.

1. Introduction

Adult spinal deformity (ASD) is a complex disease characterized by various types of lumbar or thoracolumbar deformities [1,2]. There is growing concern surrounding this disorder as it is associated with chronic back pain, imbalance, and risk of falling and decreased health-related quality of life (QoL) [3,4,5]. Sagittal imbalance results in increased energy expenditure during the standing position, and the patients compensate to maintain an upright posture and horizontal gaze through thoracic hypokyphosis, pelvic retroversion, hip extension, knee flexion, and cervical spine extension [4,6,7,8]. Recently, cervical alignment has been focused to gain adequate horizontal gaze after corrective surgery for ASD [9]. Such regional and global sagittal malalignment has been shown to correlate with pain and disability and possibly cause deformity progression due to loss of disc height and increased load on the anterior part of the vertebral body [5,10,11,12]. However, the mechanical behavior of intervertebral disc loading in all regions of the spine according to postural changes remains unclear. Biomechanical analysis is necessary to clarify the relationship between changes in sagittal alignment and disc loading.
Common parameters used to define the spinopelvic sagittal alignment of ASD include the sagittal vertical axis (SVA), pelvic tilt (PT), and pelvic incidence (PI) minus lumbar lordosis (LL) [1,2,5,10]. In addition, due to the influence of the center of gravity (COG) on spinal balance and load distribution, biomechanical analyses should include an assessment of cervical alignment and the cranial center of mass [13,14,15,16]. Several studies have evaluated intradiscal pressure in the standing posture in the thoracolumbar or lumbar region by inserting a pressure sensor inside the intervertebral discs [17,18,19,20]. However, this procedure is invasive and only reported below the lower thoracic spine (below T6). Therefore, the musculoskeletal model (MSM) is widely used to evaluate the intervertebral disc load of the entire spinal column [21,22,23].
Among MSMs, the Anybody Modeling System (AMS) is increasingly being exploited and updated worldwide and has been validated as a full-body model for the estimation of lumbar spinal loads [24]. Furthermore, it enables detailed analysis of the cervical and thoracolumbar regions [25,26,27,28]. Taking into account the cervical spine compensation mechanism in patients with ASD and the influence of the COG due to the position of the cranial, evaluation of the intervertebral disc load using an MSM should adapt the sagittal spinal alignment, including the cervical spine.
Although there have been studies on AMS with altered thoracolumbar and pelvic alignment, there are no studies on altered cervical spine alignment [29,30,31]. Therefore, this study aimed to investigate the relationship between sagittal alignment, including cervical to lumbar spine, pelvis, and lower extremity in the standing position, and disc compression force ( F m ) from the cervical to the lumbar spine analyzed by patient-specific adapted AMS that also considered cervical spine alignment in elderly women.

2. Materials and Methods

2.1. Patient-Specific Data

After institutional review board approval (approval number: 2022-4), a retrospective review was conducted to identify patients who had whole spine and lateral lower extremity radiograph. The inclusion criteria were ambulatory patients diagnosed with primary osteoporosis and no history of heavy work. Patients with a history of spine surgery, two or more vertebral fractures, hip or knee arthroplasty, or complaints of severe back pain were excluded. Finally, 14 women (78.8 ± 7.0 years) were enrolled. Written informed consent was obtained from all patients.

2.2. Image Acquisition and Anatomical Parameters Extraction

Lateral radiographs of the whole spine and lower extremities were obtained with both hands placed on the clavicle in a relaxed standing position. The spinopelvic lower-extremity alignment parameters were measured by a single author. The spinal parameters included in this analysis were the occipito–C2 angle (O–C2 angle; the angle between McGregor’s line and the inferior endplate line of C2) [32], C2–C7 angle (the angle between the lines parallel to the inferior endplate of C2 and C7 vertebral bodies) [33], C2–C7 SVA (the horizontal distance between the posterosuperior corner of the C7 vertebral body and a plumb line drawn from the centroid of C2) [34], COG–SVA (the distance between a plumb line dropped from the COG, drawn from the anterior portion of the external auditory canal, to the posterior superior endplate of S1) [35], T1 slope (the angle between a horizontal line and the superior endplate of T1) [36], thoracic kyphosis (TK; Cobb angle from the upper endplate of T4 to the lower endplate of T12), LL (Cobb angle from the upper endplate of L1 to the lower endplate of S1), and SVA (horizontal distance from the C7 plumb line that originates in the middle of the C7 vertebral body to the posterior superior endplate of S1). The pelvic parameters included in this analysis were PT (the angle between the line connecting the midpoint of the sacral plate to the bi-coxo-femoral axis and the vertical plane), sacral slope (SS; the angle between the sacral plate and the horizontal plane), and PI (the angle between the line perpendicular to the sacral plate and the line connecting the midpoint of the sacral plate to the bi-coxo-femoral axis). The lower-extremity parameters included in this analysis were the femur obliquity angle (FOA; the angle between the femoral shaft and the vertical line) [37] and knee flexion angle (KFA; the angle between the line from the hip axis to the midpoint of the bilateral notches of the femoral condyles and the line from the notch to the midpoint of the distal tibial joint surfaces) [38] (Figure 1).
The C2–C7 angle and the T1 slope range from 20–35° and 13–25° in healthy subjects, respectively [36,39]. The radiological criteria for ASD, according to the International Spine Study Group, include a frontal Cobb angle > 20°, SVA > 50 mm, TK > 60°, and/or PT > 25°. In healthy subjects, the KFA in standing position was <6° [38], and the knee joint flexion was defined as >6°.

2.3. Musculoskeletal Model

Full-body MSM was obtained using AMS (AMS. V. 6.0.5.4379; Anybody Technology, Alborg, Denmark). We previously developed an MSM in which the thoracic vertebrae were divided into 12 parts, and we modified muscle attachment points and pathways in the trunk and lower extremities based on magnetic resonance imaging data [27,28]. Furthermore, a cervical spine model, including seven vertebrae with three degrees of freedom spherical joints from T1 to C2, one degree of joint between C2 and the skull, and 136 muscles, was adopted [25] (Figure 2).

2.4. Spinopelvic and Lower-Extremity Sagittal Alignment Input and Simulation Process

The input for the entire spinal sagittal alignment was based on a previously reported method [29]. On lateral whole-spine radiographs, the center of the vertebral body was defined as the intersection of the diagonals of the quadrilateral formed by each vertebra (Figure 3).
The angle between the centers of each vertebra from C2 to S1 was then measured, and the respective values were input into the AMS. The pelvic and lower-extremity alignment inputs were similar to those in the previous method, and the posture of the model was corrected based on the patient-specific PT, SS, FOA, and KFA calculated on the radiographs [31]. This method was applied to all patients, and each patient’s body weight and height were input into the AMS (Figure 4).

2.5. Model Outputs

Using AMS, the intervertebral disc load at each vertebral level in the static standing posture was calculated using inverse dynamics analysis while considering known inertia and external forces. In the AMS, the disc load was calculated using the equation below:
F m = f m + i = 1 N f i
Here, F m , f m , i, N, and f i represented the compression force [N] at the intervertebral disc, the apparent force on the intervertebral joint, the muscle number, the total number of muscles, and the estimated muscle force of the simulation, respectively (Figure 5).
F m was estimated for each intervertebral segment from C2/C3 to L5/S1 and calculated by standardization according to each individual’s weight. The mean F m for each region was defined as follows: cervical F m from C2/C3 to C7/T1, upper thoracic F m from T1/2 to T3/4, lower thoracic F m from T4/5 to T12/L1, and lumbar F m from L1/2 to L5/S1.
Subjects were divided into two groups to compare the differences in F m at the cervical, upper thoracic, lower thoracic, and lumbar spine and each disc from C2/3 to L5/S1 between the low SVA (≤50 mm) and high SVA (>50 mm) groups (n = 7 per group). The cut-off values were based on the radiological criteria of the International Spine Study Group and Schwab’s realignment objectives for ASD [40].

2.6. Statistical Analysis

All continuous variables were expressed as mean ± standard deviation (SD). Comparisons between the two groups were performed using the Mann–Whitney U test for continuous variables. The relationship between F m and each parameter was examined using Spearman’s correlation coefficient (r). The results were verified using R version 3.6.2 (R Foundation for Statistical Computing, Vienna, Austria). Statistical significance was set at p < 0.05.

3. Results

Patient characteristics are shown in Table 1.
Compensatory changes, such as pelvic retroversion and knee flexion, were observed, and sagittal decompensation was observed as high COG–SVA and SVA (Table 1). The calculated F m (%BW) was 22.3 ± 7.3 at the cervical region, 30.8 ± 6.0 at the upper thoracic, 74.5 ± 11.8 at the lower thoracic, and 106.4 ± 44.5 at the lumbar.
A comparison of the two groups is presented in Table 2. The high SVA group had older age (p = 0.046), increased COG–SVA (p = 0.001), lower LL (p = 0.001), lower SS (p = 0.041) and increased SVA (p = 0.001). F m was significantly increased by 67.6% in the high SVA group (p = 0.046), but only in the lumber region (Table 2).
Figure 6 showed F m between the two groups at each disc level.
F m was higher in the high SVA group at the T4/5 (p = 0.008), T6/7 (p = 0.033), T12/L1 (p = 0.039), L1/2 (p = 0.044), and L5/S1 (p = 0.039) levels. The peak F m was in T12/L1 (88.4 ± 12.4) for the low SVA group and in L5/S1 (150.1 ± 55.9) for the high SVA group.
Table 3 showed the Spearman’s correlation coefficients between the mean F m at each region and spinopelvic sagittal alignment parameter.
There was a negative correlation between cervical F m with T1 slope (r = −0.589, p = 0.034) and lumbar F m with LL (r = −0.566, p = 0.035). Lumbar F m was positively correlated with COG–SVA (r = 0.615, p = 0.029), T1 slope (r = 0.613, p = 0.026), and SVA (r = 0.612, p = 0.020). There were no correlated parameters in the thoracic F m .
Figure 7 showed two cases of low SVA and high SVA.
Figure 7A showed a lateral radiograph of an 80-year-old female with a C2–C7 angle of 7°, C2–C7 SVA of 19.1 mm, T1 slope of 12°, and SVA of 51 mm. Figure 7B showed the MSM input with the patient’s alignment parameters. Figure 7C showed an 86-year-old female with a C2–C7 angle of 33°, a C2–C7 SVA of 19.2 mm, a T1 slope of 41°, and an SVA of 139 mm. Figure 7D showed the MSM. For the first patient, cervical F m was 41.7 (%BW) and lumbar F m was 86.3 (%BW). For the second patient, cervical F m was 19.8 (%BW) and lumbar F m was 230.8 (%BW). In patients with a high T1 slope and high SVA, the cervical F m was lower but the lumbar F m   was higher.

4. Discussion

In this study, the relationship between regional and global sagittal alignment parameters and disc loading was investigated using a patient-specific adapted MSM of the whole sagittal spinopelvic lower-extremity alignment for elderly women in standing posture. In the group of patients with SVA > 50 mm, the lumbar F m increased by 67.6%, and the maximum value of F m shifted to L5/S1. Cervical F m was negatively correlated with the T1 slope, lumbar F m was negatively correlated with LL, and SVA was positively correlated with the T1 slope. To the best of our knowledge, this is the first study to investigate the relationship between sagittal spinopelvic alignment and disc loading in patients with ASD. Our findings suggested that sagittal malalignment due to progressive spinal deformity resulted in further deformity due to increasing loading on thoracolumbar and lower lumbar discs; it was also associated with cervical disc loading.
Our main finding was the association between lumbar F m and sagittal alignment, with COG–SVA and SVA showing similar correlation coefficients for lumbar F m . Clinically, it has been noted that there is an association between decreased LL with increased disc degeneration and increased lumbar degenerative disease with high SVA [41,42]. Measured intradiscal pressure increased with trunk forward bending and MSM analysis also reported decreased LL and trunk forward bending, which were associated with increased disc loading and vertebral compression loads in the thoracolumbar and lumbar regions [19,20,43,44,45,46,47,48]. Studies using AMS have examined the relationship between high SVA, lumbar kyphosis, and increased L4/5 and L5/S1 disc compression forces. The findings from our study were consistent with those from previous reports [30]. Increased lumbar disc load might be caused by increased erector spinae muscle activation due to inefficient standing posture by a forward shift of the gravity line and could lead to disc degeneration [30,41,48,49,50]. Furthermore, SVA has been used for global alignment because the C7 plumb line is considered a virtual COG. This parameter is easy to measure and is related to clinical outcomes [1,8,50]. However, recent studies revealed that COG is at the cranial center of mass in both normal subjects and patients with ASD [8,51]. Therefore, we evaluated the results using COG. In our study, SVA was correlated with lumbar F m to the same degree as COG–SVA, implying that SVA measurements were useful for investigating lumbar spinal disc loading. However, since lumbar F m has a large range, it is necessary to identify other factors that may affect it.
The negative correlation between T1 slope and cervical F m suggested that the T1 slope might be useful not only as an indicator of overall sagittal spinal balance, but also in elucidating the pathology of the cervical spine. The cervical spine is a complex region with the widest range of motion and is interrelated with the thoracic and lumbar vertebrae. Compensatory mechanisms act to keep the skull above the pelvis [52]. The T1 slope correlates with SVA and is useful in evaluating global sagittal balance. T1 slope has been reported to influence cervical lordosis [36,53,54]. The correlation between the C2–C7 angle and the T1 slope could be explained by the action of cervical spine extension with large cervical lordosis in ASD patients with decreased LL and high TK [51,53,55]. Although cervical kyphosis was reportedly associated with decreased health-related QoL, there are few studies on cervical spine disease in ASD patients, and no biomechanical studies have evaluated cervical disc loading and sagittal alignment. Therefore, the effect of increased lordosis remains unclear [13,53,54,56,57,58]. Conversely, a previous study analyzed the intervertebral disc load of all vertebrae, including the cervical spine, in physiological movements using AMS and reported that the cervical disc compression force increased in flexion and decreased in extension, while the disc compression force of the upper and middle thoracic spine showed only minor changes in flexion and extension [21]. As shown in Figure 7, patients with a high T1 slope extended the cervical spine and had a low cervical F m , which could explain the negative correlation between T1 slope and cervical F m . This result was consistent with a clinical study that indicated that the degree of cervical disc herniation and cervical spinal cord compression were negatively correlated with cervical lordosis and that the degree of disc herniation and disc height might recover as cervical lordosis does [59]. Our results implied that compensatory mechanisms for cervical extension in patients with ASD did not cause an increase in cervical disc. However, in this study, cervical F m did not correlate with the C2–C7 angle, nor did thoracic F m correlate with any alignment parameters, although disc compression force in T4/5 and T6/7 was significantly higher in the high SVA group. This result should take into consideration the clinical heterogeneity of ASD patients, who have various standing postures. In this regard, most of our participants had decreased LL, while few had severe cervical or thoracic deformities.
This study had several limitations. First, sample size was small (n = 14). Second, coronal alignment was not measured or adapted to the model, and the actual muscle activity and bone morphology of the individual was not considered in the analysis. Third, only a small number of patients had local deformities of the cervical and thoracic spine. Therefore, the analysis of the association between cervical and thoracic spine alignment and disc loading needs to be validated in other patient groups. However, in this study, we focused on determining the relationship between sagittal alignment and disc compression force in patients with ASD with various standing postures. The novelty of this study was that the analysis was performed considering COG and cervical alignment in ASD patients. Applying cervical spine alignment can improve the alignment of the model, as shown in Figure 4. In future large and age-matched samples, modified models’ analysis with individual muscle activity and bone morphology are needed. Furthermore, studies focusing on local deformities would be useful in revealing the relationship between spinal deformities and entire spinal disc loading.

5. Conclusions

In the group of patients with SVA > 50 mm, the lumbar disc compression force increased by 67.6%, and the maximum value of the disc compression force shifted to L5/S1. Cervical disc compression force was negatively correlated with T1 slope. Lumbar disc compression force was negatively correlated with LL and positively correlated with global alignment. Our findings suggested that sagittal malalignment due to spinal deformity resulted in further deformity due to increasing loading on thoracolumbar and lower lumbar discs; it was also associated with cervical disc loading. Future studies are needed to assess the relationship between regional deformities in the cervical thoracic spine and disc loading.

Author Contributions

Conceptualization, T.M.; methodology, T.M., Y.K., H.K., D.K., K.S., R.K. and T.I.; software, T.I.; formal analysis, T.M.; investigation, T.M.; resources, T.M.; data curation, T.M.; writing—original draft preparation, T.M.; writing—review and editing, Y.K. and M.H.; visualization, T.M.; supervision, N.M.; project administration, N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Tazawako Hospital (approval number: 2022-4; date of approval: 11 July 2022).

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schwab, F.; Ungar, B.; Blondel, B. Scoliosis research society—schwab adult spinal deformity classification. Spine 2012, 37, 1077–1082. [Google Scholar] [CrossRef] [PubMed]
  2. Bess, S.; Schwab, F.; Lafage, V.; Shaffrey, C.I.; Ames, C.P. Classifications for adult spinal deformity and use of the scoliosis research society -schwab adult spinal deformity classification. Neurosurg. Clin. N. Am. 2013, 24, 185–193. [Google Scholar] [CrossRef] [PubMed]
  3. Miyakoshi, N.; Itoi, E.; Kobayashi, M.; Kodama, H. Impact of postural deformities and spinal mobility on quality of life in postmenopausal osteoporosis. Osteoporos Int. 2003, 14, 1007–1012. [Google Scholar] [CrossRef] [PubMed]
  4. Bess, S.; Protopsaltis, T.S.; Lafage, V.; Lafage, R.; Ames, C.P.; Errico, T.; Smith, J. Clinical and Radiographic Evaluation of Adult Spinal Deformity [Internet]. 2015. Available online: www.clinicalspinesurgery.com (accessed on 12 July 2022).
  5. Glassman, S.D.; Berven, S.; Bridwell, K.; Horton, W.; Dimar, J.R. Correlation of radiographic parameters and clinical symptoms in adult scoliosis. Spine 2005, 30, 682–688. [Google Scholar] [CrossRef] [Green Version]
  6. Barrey, C.; Roussouly, P.; Perrin, G. Sagittal balance disorders in severe degenerative spine. Can we identify the compensatory mechanisms? Eur. Spine J. 2011, 20, 626–633. [Google Scholar] [CrossRef] [Green Version]
  7. Roussouly, P.; Huec, J. Compensatory mechanisms contributing to keep the sagittal balance of the spine. Eur. Spine J. 2013, 22, 834–841. [Google Scholar]
  8. Yagi, M.; Takeda, K.; Machida, M.; Asazuma, T. Discordance of gravity line and C7PL in patient with adult spinal deformity-Factors affecting the occiput-trunk sagittal discordance. Spine J. 2015, 15, 213–221. [Google Scholar] [CrossRef]
  9. Day, L.M.; Ramchandran, S.; Jalai, C.M.; Diebo, B.G.; Liabaud, B.; Lafage, R.; Protopsaltis, T.; Passias, P.G.; Schwab, F.J.; Bess, S.; et al. Thoracolumbar realignment surgery results in simultaneous reciprocal changes in lower extremities and cervical spine. Spine 2017, 42, 799–807. [Google Scholar] [CrossRef]
  10. Schwab, F.J.; Lafage, V.; Shaffrey, C.I.; Smith, J.S.; Moal, B.; Klineberg, E.O.; Ames, C.P.; Hostin, R.; Fu, K.-M.G.; Kebaish, K.M.; et al. The SRS-schwab adult spinal deformity classification: Assessment and clinical correlations based on a prospective operative and nonoperative cohort. Neurosurgery 2013, 73, 559–568. [Google Scholar] [CrossRef]
  11. Diebo, B.G.; Shah, N.V.; Boachie-Adjei, O.; Zhu, F.; Rothenfluh, D.A.; Paulino, C.B.; Schwab, F.J.; Lafage, V. Adult spinal deformity. Lancet 2019, 394, 160–172. [Google Scholar] [CrossRef]
  12. Sparrey, C.J.; Bailey, J.F.; Safaee, M.; Clark, A.J.; Lafage, V.; Schwab, F.; Smith, J.S.; Ames, C.P. Etiology of lumbar lordosis and its pathophysiology: A review of the evolution of lumbar lordosis, and the mechanics and biology of lumbar degeneration. Neurosurg. Focus 2014, 36, E1. [Google Scholar] [CrossRef] [PubMed]
  13. Sugrue, P.A.; McClendon, J., Jr.; Smith, T.R.; Halpin, R.J.; Nasr, F.F.; O’Shaughnessy, B.A.; Koski, T.R. Redefining global spinal balance. Spine 2013, 38, 484–489. [Google Scholar] [CrossRef] [PubMed]
  14. Tang, J.A.; Scheer, J.K.; Smith, J.S.; Deviren, V.; Bess, S.; Hart, R.A.; Lafage, V.; Shaffrey, C.I.; Schwab, F.; Ames, C.P. The impact of standing regional cervical sagittal alignment on outcomes in posterior cervical fusion surgery. Neurosurgery 2012, 71, 662–669. [Google Scholar] [CrossRef] [PubMed]
  15. Schwab, F.; Lafage, V.; Boyce, R.; Skalli, W.; Farcy, J. Gravity line analysis in adult volunteers age-related correlation with spinal parameters, pelvic parameters, and foot position. Spine 2006, 31, 959–967. [Google Scholar] [CrossRef] [PubMed]
  16. Lafage, V.; Schwab, F.; Skalli, W.; Hawkinson, N.; Gagey, P.-M.; Ondra, S.; Farcy, J.-P. Standing balance and sagittal plane spinal deformity. Spine 2008, 33, 1572–1578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Polga, D.J.; Beaubien, B.P.; Kallemeier, P.M.; Schellhas, K.P.; Lew, W.D.; Buttermann, G.R.; Wood, K.B. Measurement of in vivo intradiscal pressure in healthy thoracic intervertebral discs. Spine 2004, 29, 1320–1324. [Google Scholar] [CrossRef]
  18. Schultz, A.; Andersson, G.; Ortengren, R.; Haderspeck, K.; Nachemson, A. Loads on the lumbar spine. Validation of a biomechanical analysis by measurements of intradiscal pressures and myoelectric signals. J. Bone Jt. Surg. Am. 1982, 64, 713–720. [Google Scholar] [CrossRef]
  19. Nachemson, A.; Elfström, G. Intravital dynamic pressure measurements in lumbar discs. A study of common movements, maneuvers and exercises. Scand. J. Rehabil. Med. Suppl. 1970, 1, 1–40. [Google Scholar]
  20. Takahashi, I.; Kikuchi, S.I.; Sato, K.; Sato, N. Mechanical load of the lumbar spine during forward bending motion of the trunk-a biomechanical study. Spine 2006, 31, 18–23. [Google Scholar] [CrossRef]
  21. Bayoglu, R.; Galibarov, P.E.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente spine model: A thorough investigation of the spinal loads in a complete and coherent musculoskeletal model of the human spine. Med. Eng. Phys. 2019, 68, 35–45. [Google Scholar] [CrossRef]
  22. Bayoglu, R.; Geeraedts, L.; Groenen, K.H.J.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the thoracic and cervical regions of the human spine. J. Biomech. 2017, 58, 52–63. [Google Scholar] [CrossRef] [PubMed]
  23. Arshad, R.; Schmidt, H.; El-Rich, M.; Moglo, K. Sensitivity of the cervical disc loads, translations, intradiscal pressure, and muscle activity due to segmental mass, disc stiffness, and muscle strength in an upright neutral posture. Front. Bioeng. Biotechnol. 2022, 27, 10. [Google Scholar] [CrossRef] [PubMed]
  24. Bassani, T.; Stucovitz, E.; Qian, Z.; Briguglio, M.; Galbusera, F. Validation of the AnyBody full body musculoskeletal model in computing lumbar spine loads at L4L5 level. J. Biomech. 2017, 58, 89–96. [Google Scholar] [CrossRef] [PubMed]
  25. De Zee, M.; Falla, D.; Farina, D.; Rasmussen, J. A detailed rigid-body cervical spine model based on inverse dynamics. J. Biomech. 2007, 40, S284. [Google Scholar] [CrossRef]
  26. Ignasiak, D.; Dendorfer, S.; Ferguson, S.J. Thoracolumbar spine model with articulated ribcage for the prediction of dynamic spinal loading. J. Biomech. 2016, 49, 959–966. [Google Scholar] [CrossRef] [PubMed]
  27. Higuchi, R.; Komatsu, A.; Iida, J.; Iwami, T.; Shimada, Y. Construction and validation under dynamic conditions of a novel thoracolumbar spine model with defined muscle paths using the wrapping method. J. Biomech. Sci. Eng. 2019, 14, 18-00432. [Google Scholar] [CrossRef] [Green Version]
  28. Suzuki, K.; Komatsu, A.; Hatakeyama, K.; Iida, J.; Iwami, T.; Shimada, Y. Validation of the hip contact force for a novel lower extremity muscloskeletal model. In Proceedings of the 2019 International Symposium on Micro-NanoMechatronics and Human Science (MHS), Nagoya, Japan, 1–4 December 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 1–5. Available online: https://ieeexplore.ieee.org/document/9249302/ (accessed on 23 November 2020).
  29. Iida, J.; Miyakoshi, N.; Hongo, M.; Iwami, T.; Higuchi, R.; Komatsu, A.; Matsunaga, T.; Shimada, Y. Calculation of intervertebral disc pressure in the thoracic and lumbar spine in elderly women with kyphosis using a novel musculoskeletal model with isolated thoracic vertebrae and rib cage. Open J. Orthop. 2019, 9, 241–253. [Google Scholar] [CrossRef] [Green Version]
  30. Bassani, T.; Casaroli, G.; Galbusera, F. Dependence of lumbar loads on spinopelvic sagittal alignment: An evaluation based on musculoskeletal modeling. PLoS ONE 2019, 14, e0207997. [Google Scholar] [CrossRef] [Green Version]
  31. Miura, T.; Miyakoshi, N.; Saito, K.; Kijima, H.; Iida, J.; Hatakeyama, K.; Suzuki, K.; Komatsu, A.; Iwami, T.; Matsunaga, T.; et al. Association between global sagittal malalignment and increasing hip joint contact force, analyzed by a novel musculoskeletal modeling system. PLoS ONE 2021, 16, e0259049. [Google Scholar] [CrossRef]
  32. Ota, M.; Neo, M.; Aoyama, T.; Ishizaki, T.; Fujibayashi, S.; Takemoto, M.; Nakayama, T.; Nakamura, T. Impact of the O-C2 angle on the oropharyngeal space in normal patients. Spine 2011, 36, E720–E726. [Google Scholar] [CrossRef]
  33. Ohara, A.; Miyamoto, K.; Naganawa, T.; Matsumoto, K.; Shimizu, K. Reliabilities of and correlations among five standard methods of assessing the sagittal alignment of the cervical spine. Spine 2006, 31, 2585–2591. [Google Scholar] [CrossRef] [PubMed]
  34. Ames, C.P.; Smith, J.S.; Eastlack, R.; Blaskiewicz, D.J.; Shaffrey, C.I.; Schwab, F.; Bess, S.; Kim, H.J.; Mundis, G.M.; Klineberg, E.; et al. Reliability assessment of a novel cervical spine deformity classification system. J. Neurosurg. Spine 2015, 23, 673–683. [Google Scholar] [CrossRef] [PubMed]
  35. Mizutani, J.; Verma, K.; Endo, K.; Ishii, K.; Abumi, K.; Yagi, M.; Hosogane, N.; Yang, J.; Tay, B.; Deviren, V.; et al. Global spinal alignment in cervical kyphotic deformity: The importance of head position and thoracolumbar alignment in the compensatory mechanism. Neurosurgery 2018, 82, 686–694. [Google Scholar] [CrossRef] [PubMed]
  36. Knott, P.T.; Mardjetko, S.M.; Techy, F. The use of the T1 sagittal angle in predicting overall sagittal balance of the spine. Spine J. 2010, 10, 994–998. [Google Scholar] [CrossRef] [PubMed]
  37. Huec, J.C.; Leijssen, P.; Duarte, M.; Aunoble, S. Thoracolumbar imbalance analysis for osteotomy planification using a new method: FBI technique. Eur. Spine J. 2011, 20, 669–680. [Google Scholar] [CrossRef] [Green Version]
  38. Hasegawa, K.; Okamoto, M.; Hatsushikano, S. Normative values of spino-pelvic sagittal alignment, balance, age, and health-related quality of life in a cohort of healthy adult subjects. Eur. Spine J. 2016, 25, 3675–3686. [Google Scholar] [CrossRef] [Green Version]
  39. Martini, M.L.; Neifert, S.N.; Chapman, E.K.; Mroz, T.E.; Rasouli, J.J. Cervical spine alignment in the sagittal axis: A review of the best validated measures in clinical practice. Glob. Spine J. 2021, 11, 1307–1312. [Google Scholar] [CrossRef]
  40. Schwab, F.; Patel, A.; Ungar, B.; Farcy, J.P.; Lafage, V. Adult spinal deformity—postoperative standing imbalance: How much can you tolerate? An overview of key parameters in assessing alignment and planning corrective surgery. Spine 2010, 35, 2224–2231. [Google Scholar] [CrossRef]
  41. Ergun, T.; Lakadamyali, H.; Şahin, M.Ş. The relation between sagittal morphology of the lumbosacral spine and the degree of lumbar intervertebral disc degeneration. Acta Orthop. Traumatol. Turc. 2010, 44, 293–299. [Google Scholar] [CrossRef]
  42. Barrey, C.; Jund, J.; Noseda, O.; Roussouly, P. Sagittal balance of the pelvis-spine complex and lumbar degenerative diseases. A comparative study about 85 cases. Eur. Spine J. 2007, 16, 1459–1467. [Google Scholar] [CrossRef] [Green Version]
  43. Sato, K.; Kikuchi, S.; Yonezawa, T. In vivo intradiscal pressure measurement in healthy individuals and in patients with ongoing back problems. Spine 1999, 24, 2468–2474. [Google Scholar] [CrossRef] [PubMed]
  44. Wilke, H.J.; Neef, P.; Caimi, M.; Hoogland, T.; Claes, L.E. New in vivo measurements of pressures in the intervertebral disc in daily life. Spine 1999, 24, 755–762. [Google Scholar] [CrossRef]
  45. Keller, T.S.; Colloca, C.J.; Harrison, D.E.; Harrison, D.D.; Janik, T.J. Influence of spine morphology on intervertebral disc loads and stresses in asymptomatic adults: Implications for the ideal spine. Spine J. 2005, 5, 297–309. [Google Scholar] [CrossRef] [PubMed]
  46. Harrison, D.E.; Colloca, C.J.; Harrison, D.D.; Janik, T.J.; Haas, J.W.; Keller, T.S. Anterior thoracic posture increases thoracolumbar disc loading. Eur. Spine J. 2005, 14, 234–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Bruno, A.G.; Burkhart, K.; Allaire, B.; Anderson, D.E.; Bouxsein, M.L. Spinal loading patterns from biomechanical modeling explain the high incidence of vertebral fractures in the thoracolumbar region. J. Bone Miner. Res. 2017, 32, 1282–1290. [Google Scholar] [CrossRef] [Green Version]
  48. Bruno, A.G.; Anderson, D.E.; D’Agostino, J.; Bouxsein, M.L. The effect of thoracic kyphosis and sagittal plane alignment on vertebral compressive loading. J. Bone Miner. Res. 2012, 27, 2144–2151. [Google Scholar] [CrossRef] [Green Version]
  49. Dubousset, J. Three-Dimensional Analysis of the Scoliotic Deformity; Raven Press: New York, NY, USA, 1994; pp. 479–483. Available online: https://cir.nii.ac.jp/crid/1570572700207103744 (accessed on 20 June 2022).
  50. Glassman, S.D.; Bridwell, K.; Dimar, J.R.; Horton, W.; Berven, S.; Schwab, F. The impact of positive sagittal balance in adult spinal deformity. Spine 2005, 30, 2024–2029. [Google Scholar] [CrossRef]
  51. Hasegawa, K.; Okamoto, M.; Hatsushikano, S.; Shimoda, H.; Ono, M.; Homma, T.; Watanabe, K. Standing sagittal alignment of the whole axial skeleton with reference to the gravity line in humans. J. Anat. 2017, 230, 619–630. [Google Scholar] [CrossRef]
  52. Manoharan, S.R.; Joshi, D.; Owen, M.; Theiss, S.M.; Deinlein, D. Relationship of cervical sagittal vertical alignment after sagittal balance correction in adult spinal deformity: A retrospective radiographic study. Int. J. Spine Surg. 2018, 12, 269. [Google Scholar] [CrossRef]
  53. Weng, C.; Wang, J.; Tuchman, A.; Wang, J.; Fu, C.; Hsieh, P.C.; Buser, Z.; Wang, J.C. Influence of T1 slope on the cervical sagittal balance in degenerative cervical spine. Spine 2016, 41, 185–190. [Google Scholar] [CrossRef]
  54. Staub, B.N.; Lafage, R.; Kim, H.J.; Shaffrey, C.I.; Mundis, G.M.; Hostin, R.; Burton, D.; Lenke, L.; Gupta, M.C.; Ames, C.; et al. Cervical mismatch: The normative value of T1 slope minus cervical lordosis and its ability to predict ideal cervical lordosis. J. Neurosurg. Spine 2018, 30, 31–37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kim, H.J.; Lenke, L.G.; Oshima, Y.; Chuntarapas, T.; Mesfin, A.; Hershman, S.; Fogelson, J.L.; Riew, K.D. Cervical lordosis actually increases with aging and progressive degeneration in spinal deformity patients. Spine Deform. 2014, 2, 410–414. [Google Scholar] [CrossRef]
  56. Fujimori, T.; Le, H.; Schairer, W.; Inoue, S.; Iwasaki, M.; Oda, T.; Hu, S.S. The relationship between cervical degeneration and global spinal alignment in patients with adult spinal deformity. Clin. Spine Surg. A Spine Publ. 2017, 30, E423–E429. [Google Scholar] [CrossRef] [PubMed]
  57. Schairer, W.W.; Carrer, A.; Lu, M.; Hu, S.S. The increased prevalence of cervical spondylosis in patients with adult thoracolumbar spinal deformity. J. Spinal Disord. Tech. 2014, 27, E305–E308. [Google Scholar] [CrossRef] [PubMed]
  58. Okada, E.; Matsumoto, M.; Ichihara, D.; Chiba, K.; Toyama, Y.; Fujiwara, H.; Momoshima, S.; Nishiwaki, Y.; Hashimoto, T.; Ogawa, J.; et al. Does the sagittal alignment of the cervical spine have an impact on disk degeneration? Minimum 10-year follow-up of asymptomatic volunteers. Eur. Spine J. 2009, 18, 1644–1651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Gao, K.; Zhang, J.; Lai, J.; Liu, W.; Lyu, H.; Wu, Y.; Lin, Z.; Cao, Y. Correlation between cervical lordosis and cervical disc herniation in young patients with neck pain. Medicine 2019, 98, e16545. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Sagittal alignment parameters in the standing position. (A,B) Spinal parameters, O–C2 angle (occipito–C2 angle), C2–C7 angle, C2–C7 sagittal vertical axis (C2–C7 SVA), the center of gravity of the head-posterior superior endplate of S1 (COG–SVA), T1 slope, thoracic kyphosis (TK), lumbar lordosis (LL), and SVA, were measured on a lateral whole-spine radiograph. (C) Pelvic parameters, PT (pelvic tilt), sacral slope (SS), and pelvic incidence (PI), were measured. (D) Lower-extremity parameters, the femur obliquity angle (FOA) and the knee flexion angle (KFA), were measured.
Figure 1. Sagittal alignment parameters in the standing position. (A,B) Spinal parameters, O–C2 angle (occipito–C2 angle), C2–C7 angle, C2–C7 sagittal vertical axis (C2–C7 SVA), the center of gravity of the head-posterior superior endplate of S1 (COG–SVA), T1 slope, thoracic kyphosis (TK), lumbar lordosis (LL), and SVA, were measured on a lateral whole-spine radiograph. (C) Pelvic parameters, PT (pelvic tilt), sacral slope (SS), and pelvic incidence (PI), were measured. (D) Lower-extremity parameters, the femur obliquity angle (FOA) and the knee flexion angle (KFA), were measured.
Ijerph 19 16452 g001
Figure 2. The constructed musculoskeletal model using the Anybody Modeling System. The constructed full-body model altered the muscle pathways and divided the cervical and thoracic spine.
Figure 2. The constructed musculoskeletal model using the Anybody Modeling System. The constructed full-body model altered the muscle pathways and divided the cervical and thoracic spine.
Ijerph 19 16452 g002
Figure 3. Defining the center of the vertebral body and modified musculoskeletal model. The musculoskeletal model was modified according to the sagittal spinal alignment of the patient.
Figure 3. Defining the center of the vertebral body and modified musculoskeletal model. The musculoskeletal model was modified according to the sagittal spinal alignment of the patient.
Ijerph 19 16452 g003
Figure 4. A patient’s full body lateral radiograph, and musculoskeletal model (MSM) images before and after cervical spine alignment correction. (A) Full-body lateral radiograph of a patient. (B) A modified MSM with altered thoracolumbar, pelvic, and lower-extremity alignment and unaltered cervical spine alignment. The model’s cranial location was not matched to the patient alignment. (C) In this study, a modified MSM with an altered whole spine, including cervical spine alignment, was used.
Figure 4. A patient’s full body lateral radiograph, and musculoskeletal model (MSM) images before and after cervical spine alignment correction. (A) Full-body lateral radiograph of a patient. (B) A modified MSM with altered thoracolumbar, pelvic, and lower-extremity alignment and unaltered cervical spine alignment. The model’s cranial location was not matched to the patient alignment. (C) In this study, a modified MSM with an altered whole spine, including cervical spine alignment, was used.
Ijerph 19 16452 g004
Figure 5. The magnified sagittal view of the lumbar spine in the musculoskeletal model. The green arrow indicates the penetration force exerted by an adjacent body segment. The blue arrow indicates the muscle forces acting on the vertebral body. The red arrow indicates the disc compression force, which was defined as the force perpendicular to the horizontal axis in the sagittal plane.
Figure 5. The magnified sagittal view of the lumbar spine in the musculoskeletal model. The green arrow indicates the penetration force exerted by an adjacent body segment. The blue arrow indicates the muscle forces acting on the vertebral body. The red arrow indicates the disc compression force, which was defined as the force perpendicular to the horizontal axis in the sagittal plane.
Ijerph 19 16452 g005
Figure 6. Compression force of the disc between the low SVA group and the high SVA group in each disc level. The light grey box plots indicate the disc compression force in each intervertebral body in the low SVA group, and the dark grey box plots indicate this for the high SVA group.
Figure 6. Compression force of the disc between the low SVA group and the high SVA group in each disc level. The light grey box plots indicate the disc compression force in each intervertebral body in the low SVA group, and the dark grey box plots indicate this for the high SVA group.
Ijerph 19 16452 g006
Figure 7. Example cases of low SVA and high SVA patients. (A) Lateral radiographs of the whole spine and lower extremity in the standing position in an 80-year-old female. (B) The musculoskeletal model with input of the patient’s alignment parameters. (C) Lateral radiographs of an 86-year-old female. (D) The musculoskeletal model of the patient.
Figure 7. Example cases of low SVA and high SVA patients. (A) Lateral radiographs of the whole spine and lower extremity in the standing position in an 80-year-old female. (B) The musculoskeletal model with input of the patient’s alignment parameters. (C) Lateral radiographs of an 86-year-old female. (D) The musculoskeletal model of the patient.
Ijerph 19 16452 g007
Table 1. Clinical information of all patients.
Table 1. Clinical information of all patients.
VariableValue
Total patients (n)14
Age (years)78.8 ± 7.0 (68–92)
BMI (kg/m2)22.9 ± 4.3 (16.4–32.8)
O–C2 angle (°)12.3 ± 7.6 (0–28)
C2–C7 angle (°)21.2 ± 11.0 (6–42)
C2–C7 SVA (mm)21.0 ± 8.0 (7.6–38.7)
COG–SVA (mm)86.3 ± 55.6 (15–224.7)
T1 slope (°)29.4 ± 11.6 (12–58)
TK (°)36.4 ± 15.2 (12–67)
LL (°)38.2 ± 15.8 (15–69)
PT (°)24.5 ± 8.7 (12–47)
SS (°)27.8 ± 6.7 (13–38)
PI (°)52.9 ± 9.5 (30–65)
PI−LL (°)10.9 ± 19.5 (−23–41)
SVA (mm)66.1 ± 52.5 (4–162)
FOA (°)5.8 ± 3.0 (0–13)
KFA (°)8.5 ± 4.5 (0.3–18)
Values were expressed as the mean ± standard deviation (range). BMI, body mass index; O–C2 angle, occipito–C2 angle; COG, center of gravity; SVA, sagittal vertical axis; TK, thoracic kyphosis; LL, lumbar lordosis; PT, pelvic tilt; SS, sacral slope; PI, pelvic incidence; FOA, femur obliquity angle; KFA, knee flexion angle.
Table 2. Comparison of the low SVA and high SVA groups.
Table 2. Comparison of the low SVA and high SVA groups.
VariablesLow SVA Group (n = 7)High SVA Group (n = 7)p-Value
Age (years)74.0 ± 6.782.1 ± 6.00.046 *
BMI (kg/m2)25.4 ± 3.522.0 ± 3.80.142
O–C2 angle (°)12.7 ± 7.912.8 ± 7.30.980
C2–C7 angle (°)18.4 ± 12.321.3 ± 10.70.685
C2–C7 SVA (mm)17.1 ± 8.124.7 ± 7.00.089
COG–SVA (mm)45.5 ± 15.5134.0 ± 47.4<0.001 *
T1 slope (°)24.2 ± 5.735.5 ± 13.70.070
TK (°)33.7 ± 8.435.7 ± 19.00.820
LL (°)50.6 ± 11.727.9 ± 12.5<0.001 *
PT (°)23.6 ± 7.024.7 ± 9.70.819
SS (°)32.2 ± 4.524.6 ± 6.80.041 *
PI (°)55.9 ± 7.650.0 ± 10.90.338
PI−LL (°)5.3 ± 13.015.0 ± 23.40.396
SVA (mm)21.1 ± 13.7109.9 ± 39.50.001 *
FOA (°)4.3 ± 3.16.7 ± 2.50.166
KFA (°)7.4 ± 4.69.4 ± 4.60.463
Cervical F m (%BW)23.9 ± 6.524.3 ± 7.90.912
Upper thoracic F m (%BW)29.3 ± 4.533.6 ± 6.90.225
Lower thoracic F m (%BW)70.7 ± 10.482.5 ± 9.30.067
Lumbar F m (%BW)79.2 ± 8.9132.8 ± 52.30.046 *
Values were expressed as the mean ± standard deviation; * indicates significant findings. BMI, body mass index; O–C2 angle, occipito–C2 angle; COG, center of gravity; SVA, sagittal vertical axis; TK, thoracic kyphosis; LL, lumbar lordosis; PT, pelvic tilt; SS, sacral slope; PI, pelvic incidence; FOA, femur obliquity angle; KFA, knee flexion angle; F m , disc compression force.
Table 3. Spearman’s correlation coefficients between the mean disc compression force at each part and spinopelvic sagittal alignment parameter in all patients.
Table 3. Spearman’s correlation coefficients between the mean disc compression force at each part and spinopelvic sagittal alignment parameter in all patients.
Variables Cervical   F m Upper   Thoracic   F m Lower   Thoracic   F m Lumbar   F m
rp-Valuer p-Valuerp-Valuerp-Value
O–C2 angle−0.152n.s.0.089n.s.−0.149n.s.0.116n.s.
C2–C7 angle−0.456n.s.−0.203n.s.−0.022n.s.0.247n.s.
C2–C7 SVA−0.108n.s.0.305n.s.0.165n.s.0.451n.s.
COG–SVA−0.121n.s.0.385n.s.0.280n.s.0.6150.029 *
T1 slope−0.5890.034 *0.028n.s.0.245n.s.0.6130.026 *
TK −0.203n.s.0.090n.s.0.269n.s.0.176n.s.
LL0.007n.s.−0.363n.s.−0.416n.s.−0.5660.035 *
PT−0.168n.s.0.289n.s.0.187n.s.0.119n.s.
SS0.304n.s.0.255n.s.0.092n.s.−0.139n.s.
PI−0.033n.s.0.426n.s.0.072n.s.0.028n.s.
PI−LL−0.165n.s.0.455n.s.0.455n.s.0.332n.s.
SVA −0.180n.s.0.317n.s.0.310n.s.0.6120.020 *
* Indicates significant findings; n.s. indicates not significant.   F m , disc compression force; r, Spearman’s correlation coefficient; O–C2 angle, occipito–C2 angle; COG, center of gravity; SVA, sagittal vertical axis; TK, thoracic kyphosis; LL, lumbar lordosis; PT, pelvic tilt; SS, sacral slope; PI, pelvic incidence.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miura, T.; Hongo, M.; Kasukawa, Y.; Kijima, H.; Kudo, D.; Saito, K.; Kimura, R.; Iwami, T.; Miyakoshi, N. Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model. Int. J. Environ. Res. Public Health 2022, 19, 16452. https://doi.org/10.3390/ijerph192416452

AMA Style

Miura T, Hongo M, Kasukawa Y, Kijima H, Kudo D, Saito K, Kimura R, Iwami T, Miyakoshi N. Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model. International Journal of Environmental Research and Public Health. 2022; 19(24):16452. https://doi.org/10.3390/ijerph192416452

Chicago/Turabian Style

Miura, Takanori, Michio Hongo, Yuji Kasukawa, Hiroaki Kijima, Daisuke Kudo, Kimio Saito, Ryota Kimura, Takehiro Iwami, and Naohisa Miyakoshi. 2022. "Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model" International Journal of Environmental Research and Public Health 19, no. 24: 16452. https://doi.org/10.3390/ijerph192416452

APA Style

Miura, T., Hongo, M., Kasukawa, Y., Kijima, H., Kudo, D., Saito, K., Kimura, R., Iwami, T., & Miyakoshi, N. (2022). Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model. International Journal of Environmental Research and Public Health, 19(24), 16452. https://doi.org/10.3390/ijerph192416452

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop