We retrospectively identified 2564 inpatients with LDH from L4-S1 and 2590 outpatients during the period of March 2011 and September 2021. After age and gender matched, there were 516 subjects in each group (LDH group and Control group). Between the two groups, there were no differences in demographic and clinical characteristics (Table 1). Both the inter-observer(kappa = 0.85, P༜0.001) and intra-observer (kappa = 0.96, P༜0.001) variability in this research demonstrated excellent consistency.
Table 1
Demographic data and clinical characteristics
| LDH group(n = 516) | Control group(n = 516) | P Value |
Age(years) | 52.22 ± 11.72 | 52.22 ± 11.72 | 1 |
Gender, n (%) | | | 1 |
Male | 308(59.7) | 308(59.7) | |
Female | 208(40.3) | 208(40.3) | |
BMI (kg/m2) | 23.88 ± 3.47 | 24.34 ± 6.45 | 0.18 |
Comorbidity | | | |
Hypertension | 85 | 72 | 0.29 |
Diabetes | 28 | 19 | 0.23 |
Osteoporosis | 116 | 92 | 0.07 |
Others | 86 | 78 | 0.55 |
Level of disc herniation | | | |
L4/5 | 263 | - | - |
L5/S1 | 148 | - | - |
L4/S1 | 105 | - | - |
BMI, Body Mass Index |
Prevalence of VAR in lumbar spine
The total prevalence of VAR in the entire lumbar spine was 65.7% in LDH group and 46.7% in Control group, which showed a significant difference (P༜0.001). When comparing the prevalence of VAR at each lumbar level, LDH group was significantly higher than Control group (L1, 49.8% vs 34.5%,P༜0.001; L2, 54.7% vs 39.5%,P༜0.001; L3, 37.2% vs 23.3%,P༜0.001; L4, 10.7% vs 3.5%,P༜0.001; L5, 3.5% vs 0.8%,P = 0.03). And, it showed that either in close or distant neighboring vertebrae, the prevalence of VAR was significantly higher in LDH group than Control group.
According to our findings, the prevalence of VAR was not significantly different between males and females in the LDH group (Total, 65.6% vs 65.9%, P = 0.947; L1, 49.7% vs 50.0%, P = 0.942; L2, 55.5% vs 53.4%, P = 0.63; L3, 36.7% vs 38.0%, P = 0.766; L4, 10.7% vs 10.6%, P = 0.96; L5, 3.3% vs 3.9%, P = 0.716). Similarly, there was no evidence of a gender difference in the Control group (Total, 45.8% vs 48.1%, P = 0.608; L1, 33.4% vs 36.1%, P = 0.54; L2, 39.3% vs 40.4%, P = 0.802; L3, 21.8% vs 25.5%, P = 0.326; L4, 3.9% vs 2.9%, P = 0.539; L5, 0% vs 1.9%, P = 0.608).
Distribution characteristic of VAR in lumbar spine
By measuring of the entire lumbar vertebral rotation, we found that in LDH group, the L2 vertebra had the highest prevalence of VAR (49.5%), followed by L1(45.1%), and from L3 to L5, the prevalence decreased gradually(33.6%,8.9%,3.1%). The Control group also noticed a similar tendency (L2,39.8%; L1,34.6%; L3,23.2%; L4,3.1%; L5,0.8%) (Fig. 1a).
The subgroup analysis showed that when the disc herniation occurred at L4-L5 segment, the difference of VAR is significant in each lumbar vertebra from L1 to L5(P༜0.05) ( Fig. 1b). When the disc herniation occurred at L5-S1, the difference is evident in L2, L3, L4, L5, but not in L1 (Fig. 1c). When the disc herniation occurred at L4-S1, the difference is significant in L1, L2 and L3, but not in L4 and L5 (Fig. 1d). It appeared that disc herniation was not only associated with increased prevalence of VAR in close neighboring vertebrae, but also in distant neighboring vertebrae.
Table 2 displayed a thorough comparison of VAR in various segments of disc herniation. There is no significant difference of VAR in either entire lumbar spine or close neighboring vertebrae among the subgroups. Yet, distant neighboring vertebrae showed a significant difference of VAR. L5-S1 subgroup displayed lower prevalence of VAR compared with L4-L5 subgroup (35.41% vs 47.53%, P = 0.015) or L4-S1 subgroup (35.41% vs 52.07%, P = 0.006); while the difference between L4-L5 subgroup and L4-S1 subgroup were not statistical (47.53% vs 52.07%, P = 0.4).
Table 2
Subgroup analysis of lumbar vertebral axial rotation at different disc herniation segment
| L4/5 (n = 263) | L5/S1 (n = 148) | L4/S1 (n = 105) | P value | P value (L4/5 vs L5/S1) | P value (L5/S1 vs L4/S1) | P value (L4/5 vs L4/S1) |
Total ✝ | 409(31.10%) | 216(29.19%) | 176(33.52%) | 0.766 | 0.653 | 0.468 | 0.689 |
Close neighboring vertebrae✝ | 34(6.46%) | 8 (5.41%) | 12 (5.71%) | 0.901 | 0.666 | 0.916 | 0.789 |
Distant neighboring vertebrae✝ | 375(47.53%) | 208(35.41%) | 164(52.07%) | 0.012* | 0.015# | 0.006# | 0.400 |
✝,the number and percent of vertebrae axial rotation; *,P value༜0.05; The pairwise comparison among these three subgroups was performed by Tukey’s post hoc test, and corrected α was 0.017 ; #, P value༜0.017 |
The risk factors for VAR in patients with LDH
According to existence of VAR or not, the LDH group was divided into two groups (Rotation group and Non-rotation group). It indicated that each vertebra's Nash-Moe index in the Non-rotation group was 0. While in Rotation group, at least one vertebra was more than Nash-Moe grade 1. The results of univariate analysis revealed that mean age of Rotation group was significantly greater than Non-rotation group(P༜0.001).There was no significant difference in gender, BMI, status of smoking and drinking, and clinical medical comorbidities, with the exception of osteoporosis (P = 0.003). Moreover, there was no statistical difference in VAS back and symptom duration between two groups(Table 3). Age and osteoporosis were suggested to be two potential predicative variables for VAR. After stepwise forward logistic regression, the results showed that age is the only risk factor for VAR (OR = 1.245, 95%CI [1.111,1.394], P༜0.001).༈Table 4༉. In addition, ROC curve indicated that age is associated with incidence of VAR in patients with LDH(AUC = 0.589, cutoff point was 43.5 years)༈Fig. 2༉.
Table 3
Difference of characteristics in LDH group with and without VAR
| Rotation group(n = 265) | Non-rotation group(n = 251) | P value |
Age, years✝ | 46.7(14.9) | 42.3(13.7) | ༜0.001* |
Gender, n (%) | | | 0.615 |
Male | 133(50.2) | 131(52.2) | |
Female | 132(49.8) | 120(47.8) | |
BMI, kg/m2✝ | 24.1(3.5) | 23.7(3.4) | 0.138 |
Hypertension | 50 | 35 | 0.154 |
Diabetes | 16 | 12 | 0.565 |
Osteoporosis | 74 | 42 | 0.003* |
Smoking | 47 | 37 | 0.404 |
Drinking | 17 | 16 | 0.985 |
VAS Back☨ | 4(3) | 4(5) | 0.076 |
Symptom duration☨ | 36(76) | 24(54) | 0.128 |
✝, mean and standard deviation; BMI, Body Mass Index; ☨,median and interquartile range; *, P༜0.05 |
Table 4
Multivariate logistic regression analysis of risk factors for vertebral axial rotation in patients with lumbar disc herniation
| OR | 95% CI | P value |
lower | upper |
Age, per year | 1.022 | 1.011 | 1.034 | ༜0.001 |
OR, Odds ratio; CI, Confidence interval |