- Open Access
Fat infiltration of paraspinal muscles as an independent risk for bone nonunion after posterior lumbar interbody fusion
BMC Musculoskeletal Disorders volume 23, Article number: 232 (2022)
The prognosis value of paraspinal muscle degeneration on clinical outcomes has been revealed. However no study has investigated the effect of the fat infiltration (FI) of paraspinal muscles on bone nonunion after posterior lumbar interbody fusion (PLIF).
Three hundred fifty-one patients undergoing PLIF for lumbar spinal stenosis with 1-year follow-up were retrospectively identified. Patients were categorized into bone union (n = 301) and bone nonunion (n = 50) groups based on dynamic X-ray at 1-year follow-up. The relative total cross-sectional area (rTCSA) and FI of multifidus (MF) and erector spinae (ES), and the relative functional CSA (rFCSA) of psoas major (PS) were measured on preoperative magnetic resonance imaging.
The nonunion group had a significantly higher MF FI and a higher ES FI and a smaller MF rTCSA than the union group (p = 0.001, 0.038, 0.026, respectively). Binary logistic regression revealed that MF FI (p = 0.029, odds ratio [OR] = 1.04), lumbosacral fusion (p = 0.026, OR = 2193) and length of fusion (p = 0.001, OR = 1.99) were independent factors of bone nonunion. In subgroup analysis, in one or two-level fusion group, the patients with nonunion had a higher MF FI and a higher ES FI than those of the patients with union (all p < 0.05). Similarly, in lumbosacral fusion group, the patients with nonunion had a higher MF FI and a higher ES FI than those of the patients with union (all p < 0.05). The logistic regressions showed that MF FI remained an independent factor of bone nonunion both in the patients with one or two-level fusion (p = 0.003, OR = 1.074) and in the patients with lumbosacral fusion (p = 0.006, OR = 1.073).
Higher fatty degeneration was strongly associated with bone nonunion after PLIF. Surgeons should pay attention to the FI of paraspinal muscles when performing posterior surgery for patients, especially those who need short-segment fusion or to extend fusion to S1.
Posterior lumbar interbody fusion (PLIF) is a commonly surgical treatment for degenerative spinal diseases so as to stabilize the motion segment, restore lordosis and correct deformity . However, bone nonunion as a surgical complication can be observed during follow-up. Most recent studies have reported that the rate of bone nonunion ranged from 7 to 20% [1,2,3]. The process of bone union can be affected by multiple factors including increased fused level, osteoporosis and obesity [4, 5].
The predictive value of paraspinal muscle morphometry on operative complications has been investigated [6,7,8,9]. Some studies reported that decreased cross-sectional area (CSA) of paraspinal muscles was correlated to bone nonunion rate in patients with lumbar surgery [3, 7]. However, the effect of paraspinal muscles fat infiltration (FI) on bone nonunion remains indistinct. We hypothesized that the patients with higher FI were inclined to occur bone nonunion. This study aimed to examine the relationship between FI of paraspinal muscles and bone nonunion in patients with lumbar spinal stenosis (LSS) after PLIF.
Hospitalized patients undergoing PLIF for LSS between July 2011 and December 2015 were reviewed. Inclusion criteria were: (1) aged ≥50 years, (2) underwent lumbar magnetic resonance imaging (MRI) and lumbar computed tomography (CT) within 3 months before the index surgery, (3) underwent follow-up of ≥12 months. Exclusion criteria were (1) previous spinal surgery, (2) patients with bone tumor, ankylosing spondylitis, diffuse idiopathic skeletal hyperostosis, rheumatoid arthritis, tuberculosis, or secondary osteoporosis, and (3) patients with scoliosis (> 10°).
A total of 351 patients were identified. Among them, 244 patients with PLIF were LSS and 107 were LSS combined with degenerative lumbar spondylolisthesis. Degenerative lumbar spondylolisthesis was defined as displacement of 1 vertebra over subjacent vertebra using Meyerding grading system . All surgical strategies and approaches were discussed and decided before surgery. For PLIF procedures, using the posterior midline approach, meticulous exposure of the spine and posterior decompression fusion and fixation with pedicle screw was performed. After pedicle screws had been implanted, the neural decompression by laminectomy and discectomy was performed. A polyetheretherketone (PEEK) cage packed with autogenous bone was placed into the interbody space for all patients. Posterolateral fusion was also performed simultaneously. The autograft was harvested from decompression. No bone morphogenetic protein has been used in these patients.
Bone union evaluation
Segmental fusion status was evaluated by dynamic X-ray at 1-year follow-up. We defined the bone nonunion as 1) there was no continued bone fusion mass at any fusion segment; 2) any motion (greater than 3 mm or 3°) on flexion/ extension plain radiographs [11, 12]. Based on dynamic X-ray, patients were categorized into bone union (n = 301) and bone nonunion (n = 50) groups.
Bone density evaluation
In consideration of the overestimation of the BMD of the lumbar spine in patients with lumbar degenerative diseases evaluated by dual-energy X-ray absorptiometry, three-dimensional reconstructive lumbar CT (Siemens, DEFINITION, tube voltage 120 kV) were performed preoperatively to measure the bone density. The Hounsfield unit (HU) value of L1 to L4 was measured for each patient according to the method of previous studies . An oval region of interest inclusive of trabecular bone was placed in the middle-axial CT image of vertebral body (Fig. 1). The cortical bone and posterior venous plexus were excluded in the measurement. The average HU value of L1-4 was calculated.
Paraspinal muscle evaluation on MRI
All enrolled patients had undergone preoperative MRI of lumbar area with Signa HDxt 3.0 T (General Electric Company). We measured the multifidus (MF), erector spinae (ES) and psoas major (PS) bilaterally from T2-weighted images at the center of the intervertebral disc of L4-5 level. The following parameters were measured on each level by the Image J software (National Institutes of Health, Bethesda, MD, USA; Fig. 2): total cross-sectional area (TCSA) of MF, ES and intervertebral disc; FI of MF and ES was measured by the previously reported thresholding technique [14, 15]; For PS, only functional cross-sectional area (FCSA) was measured due to the ill-defined outline of intramuscular fat and soft tissue . Thresholding technique can identify two different signal intensity peaks and classify the pixel areas with lower intensity peaks as muscle tissue and the pixel areas with higher intensity peaks as intramuscular fat. Relative cross-sectional area (rCSA, the ratio of cross-sectional area of muscle to that of disc at the same level) was introduced to reduce the effect of body shape on muscular parameters [8, 16]. rCSA of both total muscle (T) and functional muscle (F) were marked as rTCSA and rFCSA.
To test the reliability, all muscular parameters of 10 patients were randomly selected and were measured by two observers independently. After 3 weeks, the same measurements were performed by each observer. The ICCs for both intra-rater and inter-rater reliability of MF rTCSA, ES rTCSA, MF FI, ES FI and PS rFCSA were > 0.8 (Supplement Table 1).
The Mann-Whitney U test or Analysis of Variance (ANOVA) test (for continuous data) and Chi-square test (for categorical data) were conducted to determine the statistical difference of clinical characteristics and paraspinal muscle degeneration between the union group and the nonunion group. Binary logistic regression model was used to identify the independent risk factors of bone nonunion. Age, sex, HU value, smoking, lumbosacral fusion, number of fusion levels, and muscular parameters that had significant differences between groups were included in the regression model. Furthermore, we have performed subgroup analyses according to fusion length and lumbosacral fusion to compare the difference of clinical characteristics and paraspinal muscle between the union and the nonunion groups. Binary logistic regressions were also performed in subgroups. Intraclass correlation coefficient was calculated to test the intra- and inter-rater reliability. Statistical significance was set at P value < 0.05. All statistical analyses were performed using SPSS 22.0 (IBM Corp).
The bone nonunion occurred in 50 (14.2%) patients. Compared with union group, nonunion group had significantly older age, higher rate of fusions extending to the sacrum, longer fusion length, lower mean HU value of L1-4 and higher rate of smoking (p < 0.001, p = 0.002, p < 0.001, p = 0.013, p = 0.036, respectively). However, the gender, BMI, whether combined with spondylolisthesis or diabetes were not significantly different between two groups (all p > 0.05) (Table 1).
Paraspinal muscle characteristics
As shown on Table 2, MF FI and ES FI at L4 were significant higher in nonunion group than in union group (p = 0.001, 0.038, respectively). Besides, those who did not achieve union had a smaller MF rTCSA than patients with union (p = 0.026). However, there were no significant difference of ES rTCSA and PS rFCSA between two groups (both p > 0.05).
Furtherly, binary logistic regression revealed that MF FI of L4 was an independent factor of bone nonunion (p = 0.029, OR = 1.04; Table 3). Lumbosacral fusion was also an independent factor of bone nonunion (p = 0.026, OR = 2.193; Table 3). Besides, the number of fusion levels had a dramatically negative impact on union status (p = 0.001, OR = 1.99; Table 3).
Subgroup analyses according to fusion length and lumbosacral fusion
In subgroup analysis, the patients were divided into 2 sets of subgroups according to the fusion length and whether lumbosacral fusion was performed respectively. In the patients with one or two-level fusion, the nonunion group had a higher MF FI and a higher ES FI than those of the union group (p = 0.001, 0.014, respectively; Table 4). However, there was no difference of muscular parameters between the two groups in the patients with three or more-level fusion (all p > 0.05; Table 4). The logistic regression showed that MF FI was still an independent factor of bone nonunion in the patients with one or two-level fusion (p = 0.003, OR = 1.074; Supplement Table 2). Additionally, in the patients with lumbosacral fusion, the nonunion group had a higher MF FI and a higher ES FI than those of union group (p < 0.001, 0.021, respectively; Table 5), whereas the nonunion group had relatively higher MF FI and ES FI than those of union group without significance (all p > 0.05; Table 5). In the patients with lumbosacral fusion, MF FI remained an independent factor of bone nonunion in the logistic regression (p = 0.006, OR = 1.073; Supplement Table 3).
Our study showed that the patients who did not achieve union had a higher MF FI and a higher ES FI than those of the patients with union. In Lee et al’s study, it was reported that the union rate decreased as fat content of extensor muscles increased, which was accordant to our findings . However, they only used a semiquantitative scale to quantify the FI and did not investigate the MF and ES separately. Considering that MF is in the deep attaching to the lumbar vertebrae while ES is more superficial spanning more sections of the spine, evaluating them separately is reasonable . Our results were also consistent to Katsu et al’s study that focused on the patients with osteoporotic fractures . They found that FI of MF and ES were both higher in insufficient union group than in union group. Several studies have revealed that increased muscle FI was correlated to poorer muscle strength [20, 21]. Previous study indicated that incremental bending moment transmitted by the internal fixation device would increase the risk of bone nonunion . Consequently, paraspinal muscles with higher FI might be less effective on reducing the bending moment. Besides, desired paraspinal musculature could provide important vascular ingrowth into the fusion site and accelerate the bone healing . It is suggested that severe muscular degeneration might impede this process.
Of note, multifactor analysis demonstrated that FI of MF rather than ES, had an effect on nonunion. Liu et al. investigated 118 LSS patients and found that the postoperative improvement of Oswestry Disability Index (ODI) was significantly less in MF FI ≥ 25% group than in MF FI < 25% group . Besides, Hong et al. found that MF FI also contributed to superior clinical outcomes including less improvement in ODI, greater postoperative pain and higher reoperation rate . As MF is the innermost and largest one of the paraspinal muscles and provides two-thirds of spinal segmental stability , MF might have a more remarkable effect on clinical outcomes compared to ES.
Our analysis showed that the nonunion group had a greater MF atrophy than that of union group, but no statistical significance was seen in multivariate analysis. Choi et al’s study reported that CSA of MF, ES and PS in the nonunion group were all smaller than those in the union group . Furthermore, they found that only PS TCSA was correlated to fusion rate in multivariate analysis. The reason why the TCSA did not show an arresting effect on bone nonunion might be that the relationship between TCSA and muscle strength was not as significant as that of FI. A study demonstrated that FI of paraspinal muscles, not CSA, remained a significant predictor of extensor strength in multivariate regression . Our study indicated that surgeons should pay more attention to the FI of paraspinal muscles rather than atrophy when evaluating the risk of bone nonunion preoperatively.
Multivariable analysis showed that the number of levels fused significantly affected the incidence of nonunion. In previous studies, the number of fusion levels was considered to be a crucial factor in achieving solid fusion in degenerative lumbar diseases [17, 22]. In view of the possible correlation between fusion length and the degree of preoperative FI, we performed a subgroup analysis. The results exhibited that the effect of FI only existed in the patients with one or two fused levels, not in the patients with three or more fused levels. We speculated that in the patients with long-segment fusion, the fusion length had a more notable effect on the bone nonunion over paraspinal muscles FI, hence a severe fatty degeneration might not notably increase the risk of nonunion. Nevertheless, in the patients with shorter fused levels, FI of paraspinal muscles began to take effect.
In addition, we found that the patients with fusion to S1 had a higher rate of bone nonunion, which was consistent with previous studies [17, 22, 27]. In subgroup analysis, the significant difference of FI between the union group and the nonunion group only existed in the patients with lumbosacral fusion, not in the patients without lumbosacral fusion. It could be interpreted by that as a great mechanical load could be applied to the sacrum in lumbosacral fusion, patients need stronger paraspinal muscles to counteract this negative effect . Once patients have a higher FI of muscles preoperatively, lumbosacral fusion will highlight the effect of muscles and then the risk of bone nonunion will increase.
Our findings indicated that surgeons should pay attention to the FI of paraspinal muscles when performing posterior surgery for patients who need short-segment fusion or to extend fusion to S1. In the above cases, preventive measures such as the use of materials to facilitate bone grafting or screws with greater fixation strength should be considered.
We recognize some limitations in the study. First, there were no postoperative MRIs to evaluate the condition of muscle injury during the operation, which might reduce the predictive value of paraspinal muscles on bone nonunion. Besides, we did not perform CT to evaluate the fusion status which may reduce the reliability of bone nonunion. While in our study, 2 observers evaluated union independently with strict criteria for defining nonunion in order to increase the accuracy. Third, the heterogeneity of cases and the small number of cases for subgroup analysis might produce bias. In addition, we have not taken the size and length of the screw into account, which might be related to bone healing . Moreover, it may be impractical to actually measure FI in clinical practice since automated measurement software is not yet available.
This is the first study focusing on the prognostic value of back muscles FI to predict bone nonunion after PLIF. We demonstrated that higher fatty degeneration of MF was an independent factor of nonunion. Furtherly, the effect of FI only existed in the patients with one or two fused levels, not in the patients with three or more fused levels. Besides, the significant difference of FI between the union group and nonunion groups only existed in the patients with lumbosacral fusion, not in the patients without lumbosacral fusion. In cases with higher MF FI during preoperative evaluation, we considered that more rigid fixation or more graft bone might be necessary.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Functional cross-sectional area
Lumbar spinal stenosis
Magnetic resonance imaging
Oswestry Disability Index
Posterior lumbar interbody fusion
Relative cross-sectional area
Total cross-sectional area
Meng B, Bunch J, Burton D, Wang J. Lumbar interbody fusion: recent advances in surgical techniques and bone healing strategies. Eur Spine J. 2021;30(1):22–33.
Mo GY, Guo HZ, Guo DQ, Tang YC, Li YX, Yuan K, et al. Augmented pedicle trajectory applied on the osteoporotic spine with lumbar degenerative disease: mid-term outcome. J Orthop Surg Res. 2019;14:170.
Formica M, Vallerga D, Zanirato A, Cavagnaro L, Basso M, Divano S, et al. Fusion rate and influence of surgery-related factors in lumbar interbody arthrodesis for degenerative spine diseases: a meta-analysis and systematic review. Musculoskelet Surg. 2020;104:1–15.
Lee C, Dorcil J, Radomisli TE. Nonunion of the spine: a review. Clin Orthop Relat Res. 2004;419:71–5.
Vaccaro AR, Chiba K, Heller JG, Patel T, Thalgott JS, Truumees E, et al. Bone grafting alternatives in spinal surgery. Spine J. 2002;2:206–15.
Choi MK, Kim SB, Park BJ, Park CK, Kim SM. Do trunk muscles affect the lumbar interbody fusion rate?: correlation of trunk muscle cross sectional area and fusion rates after posterior lumbar interbody fusion using stand-alone cage. J Korean Neurosurg Soc. 2016;59:276–81.
Choi MK, Kim SB, Park CK, Malla HP, Kim SM. Cross-sectional area of the lumbar spine trunk muscle and posterior lumbar interbody fusion rate: a retrospective study. Clin Spine Surg. 2017;30:E798-803.
Leng J, Han G, Zeng Y, Chen Z, Li W. The effect of paraspinal muscle degeneration on distal pedicle screw loosening following corrective surgery for degenerative lumbar scoliosis. Spine (Phila Pa 1976). 2020;45:590–8.
Yuan L, Zeng Y, Chen Z, Li W, Zhang X, Mai S. Degenerative lumbar scoliosis patients with proximal junctional kyphosis have lower muscularity, fatty degeneration at the lumbar area. Eur Spine J. 2020;30:1133–43.
Bydon M, Alvi MA, Goyal A. Degenerative lumbar spondylolisthesis: Definition, natural history, conservative management, and surgical treatment. Neurosurg Clin N Am. 2019;30:299–304.
Zou D, Muheremu A, Sun Z, Zhong W, Jiang S, Li W. Computed tomography Hounsfield unit-based prediction of pedicle screw loosening after surgery for degenerative lumbar spine disease. J Neurosurg Spine. 2020;32:1–6.
Larsen JM, Rimoldi RL, Capen DA, Nelson RW, Nagelberg S, Thomas JC Jr. Assessment of pseudarthrosis in pedicle screw fusion: a prospective study comparing plain radiographs, flexion/extension radiographs, CT scanning, and bone scintigraphy with operative findings. J Spinal Disord. 1996;9:117–20.
Pickhardt PJ, Pooler BD, Lauder T, del Rio AM, Bruce RJ, Binkley N. Opportunistic screening for osteoporosis using abdominal computed tomography scans obtained for other indications. Ann Intern Med. 2013;158:588–95.
Wang W, Sun Z, Li W, Chen Z. The effect of paraspinal muscle on functional status and recovery in patients with lumbar spinal stenosis. J Orthop Surg Res. 2020;15:235.
Dohzono S, Toyoda H, Takahashi S, Matsumoto T, Suzuki A, Terai H, et al. Factors associated with improvement in sagittal spinal alignment after microendoscopic laminotomy in patients with lumbar spinal canal stenosis. J Neurosurg Spine. 2016;25:39–45.
Yu B, Jiang K, Li X, Zhang J, Liu Z. Correlation of the features of the lumbar multifidus muscle with facet joint osteoarthritis. Orthopedics. 2017;40:e793–800.
Lee C-S, Chung SS, Choi SW, Yu JW, Sohn MS. Critical length of fusion requiring additional fixation to prevent nonunion of the lumbosacral junction. Spine. 2010;35:E206–11.
Meakin JR, Fulford J, Seymour R, Welsman JR, Knapp KM. The relationship between sagittal curvature and extensor muscle volume in the lumbar spine. J Anat. 2013;222:608–14.
Katsu M, Ohba T, Ebata S, Haro H. Comparative study of the paraspinal muscles after OVF between the insufficient union and sufficient union using MRI. BMC Musculoskelet Disord. 2018;19:143.
Schlaeger S, Inhuber S, Rohrmeier A, Dieckmeyer M, Freitag F, Klupp E, et al. Association of paraspinal muscle water-fat MRI-based measurements with isometric strength measurements. Eur Radiol. 2019;29:599–608.
Goodpaster BH, Park SW, Harris TB, Kritchevsky SB, Nevitt M, Schwartz AV, et al. The loss of skeletal muscle strength, mass, and quality in older adults: the health, aging and body composition study. J Gerontol A Biol Sci Med Sci. 2006;61:1059–64.
Kornblatt MD, Casey MP, Jacobs RR. Internal fixation in lumbosacral spine fusion. A biomechanical and clinical study. Clin Orthop Relat Res. 1986;(203):141–50.
Bawa M, Schimizzi AL, Leek B, Bono CM, Massie JB, Macias B, et al. Paraspinal muscle vasculature contributes to posterolateral spinal fusion. Spine (Phila Pa 1976). 2006;31:891–6.
Yang L, Yuzeng L, Yong H, Tie L, Li G, Chen X. Fat infiltration in the multifidus muscle as a predictor of prognosis after decompression and fusion in patients with single-segment degenerative lumbar spinal stenosis: an ambispective cohort study based on propensity score matching. World Neurosurg. 2019;128:E989–E1001.
Hong X, Shi R, Wang Y-T, Liu L, Bao J-P, Wu X-T. Lumbar disc herniation treated by microendoscopic discectomy: prognostic predictors of long-term postoperative outcome. Orthopade. 2018;47:993–1001.
Wilke HJ, Wolf S, Claes LE, Arand M, Wiesend A. Stability increase of the lumbar spine with different muscle groups. A biomechanical in vitro study. Spine (Phila Pa 1976). 1995;20:192–8.
Aono H, Takenaka S, Nagamoto Y, Tobimatsu H, Yamashita T, Furuya M, et al. Fusion rate and clinical outcomes in two-level posterior lumbar interbody fusion. World Neurosurg. 2018;112:e473–8.
Otsuki B, Fujibayashi S, Tanida S, Shimizu T, Murata K, Matsuda S. Possible association of pedicle screw diameter on pseudoarthrosis rate after transforaminal lumbar interbody fusion. World Neurosurg. 2021;150:e155–61.
This work was supported by the Clinical Cohort Construction Program of Peking University Third Hospital, China National Key Research and Development Program (2018YFB1307700) and Engineering Research of Bone and Joint Precision Medicine.
Ethics approval and consent to participate
This retrospective study was approved by Peking University Third Hospital Medical Science Research Ethics Committee. Peking University Third Hospital Medical Science Research Ethics Committee granted a waiver of informed consent (M2021134). All methods were performed in accordance with the Declaration of Helsinki.
Consent for publication
The authors declare no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Intra-rater and inter-rater reliability of paraspinal muscle parameters using intraclass correlation coefficient. Supplement Table 2. Independent risk factors of bone nonunion identified by logistic regression in the one or two-level fusion group. Supplement Table 3. Independent risk factors of bone nonunion identified by logistic regression in the lumbosacral fusion group.
About this article
Cite this article
Han, G., Zou, D., Liu, Z. et al. Fat infiltration of paraspinal muscles as an independent risk for bone nonunion after posterior lumbar interbody fusion. BMC Musculoskelet Disord 23, 232 (2022). https://doi.org/10.1186/s12891-022-05178-z
- Bone nonunion
- Fat infiltration
- Paraspinal muscle
- Posterior lumbar interbody fusion