两种类型滑脱的脊柱骨盆矢状位参数与骶骨参数差异及相关性分析

Difference and correlation analysis of spinal pelvic sagittal parameters and sacral parameters between two types of spondylolisthesis

摘要

背景 既往研究发现骶骨参数与脊柱骨盆矢状位参数存在相关性，而对不同类型滑脱患者骶骨形态参数研究较少。 目的 分析不同类型滑脱患者脊柱骨盆矢状位参数及骶骨形态参数的特征及其相关性，探讨其在不同类型滑脱发生机制中的作用。 方法 回顾性分析2019年1月— 2023年12月就诊于解放军总医院第一医学中心的腰椎滑脱患者的脊柱站立位全长片，选取同期因体格检查就诊的健康成年人作为对照。测量脊柱骨盆矢状位参数，包括胸椎后凸角(thoracic kyphosis，TK)，腰椎前凸角(lumbar lordosis，LL)，骨盆入射角(pelvic incidence，PI)，骨盆倾斜角(pelvic tilt，PT)，骶骨倾斜角(sacral slope，SS)，矢状面轴向距离(sagittal vertical axis，SVA)，滑脱角(slip angle，SA)。测量骶骨形态参数，包括骶1-2倾斜角(s1-2 tilt，S1-2T)，骶骨后凸(sacral kyphosis，SK)，S5Co1倾斜线水平角(the horizontal angle of S5Co1，S5Co1HA)，比较两组间参数有无差异。组内脊柱骨盆矢状位参数与骶骨参数进行相关性分析，并建立线性回归方程。 结果 滑脱组91例，年龄(61.45±11.56)岁。其中退行性滑脱57例，男26例，女31例，年龄(65.49±8.29)岁；峡部裂性滑脱34例，男18例，女16例，年龄(54.68±13.13)岁；对照组35例，男16例，女19例，年龄(64.77±10.47)岁。3组间性别分布差异无统计学意义(P＞0.05)，年龄分布差异有统计学意义(P＜0.01)。退行性滑脱组、峡部裂型滑脱组TK(29.78°±13.71°和30.56°±13.61° vs 38.94°± 10.51°)、低于对照组，PI(54.06°±10.99°和56.10°±10.76° vs 49.33°±8.92°)、PT(22.84°±8.68°和21.45°±9.31° vs 12.32°±5.86°)均高于对照组(P＜0.05)。退行性滑脱组的LL(40.42°±20.44° vs 49.25°±11.17°)、SS(31.42°±10.59° vs 37.01°±8.11°)低于对照组(P＜0.05)，SA显著低于峡部裂型滑脱组(P＜0.05)。退行性滑脱组及峡部裂型滑脱组的SK(57.23°±13.29°和56.01°±11.69° vs 50.35°±9.77°)高于对照组(P＜0.05)，S5Co1HA(80.70°±9.68°和81.45°±11.56° vs 94.45°±12.02°)低于对照组(P＜0.05)。退行性滑脱组的S1-2T(17.09°±7.08° vs 13.15°±9.08°)高于峡部裂型滑脱组(P＜0.05)。多元线性回归方程：退行性滑脱组PI=0.383SK-0.753(S1-2T)+44.989(R=0.727)，峡部裂型滑脱组PI=0.529SK-0.876(S1-2T)+38.001(R=0.815)。 结论 退变性与峡部裂性腰椎滑脱患者表现出相似的骶骨形态学改变特征，包括骶骨后凸角显著增大和S5Co1倾斜线水平角明显减小。这些特征性改变与脊柱-骨盆矢状位失衡存在显著关联，其中骶骨后凸角与骶1-2倾斜角的组合可有效预测骨盆入射角。研究结果提示，骶骨形态参数可作为评估腰椎滑脱患者脊柱-骨盆平衡状态的重要补充指标，这为临床制定个性化治疗方案提供了新的解剖学参考。

Abstract

Background Previous studies have found correlations between sacral parameters and spinal-pelvic sagittal parameters, but there is limited research on the sacral morphology parameters in patients with different types of spondylolistheses. Objective To analyze the characteristics and correlations of spinal pelvic sagittal parameters and sacral morphology parameters in patients with different types of spondylolistheses, and explore their roles in the mechanisms of different types of spondylolistheses. Methods A retrospective analysis was conducted on the full-length standing spine radiographs of the patients with lumbar spondylolisthesis treated in the Department of Orthopedics of the First Medical Center of PLA General Hospital from January 2019 to December 2023, the healthy adults who visited for physical examination during the same period were selected as controls. Spinal-pelvic sagittal parameters were measured, including thoracic kyphosis (TK), lumbar lordosis (LL), pelvic incidence (PI), pelvic tilt (PT), sacral slope (SS), sagittal vertical axis (SVA), and slip angle (SA). Sacral morphological parameters were measured, including S1-2 tilt (S1-2T), sacral kyphosis (SK), and the horizontal angle of S5Co1 (S5Co1HA). Independent sample t-tests was used to compare parameters between groups. Correlation analysis was performed between spinal pelvic sagittal parameters and sacral parameters within groups, and linear regression equations were obtained. Results A total of 91 patients with spondylolisthesis were included in the study group, with a mean age of (61.45 ± 11.56) years. Among them, 57 patients had degenerative spondylolisthesis (26 males and 31 females, mean age of 65.49 ± 8.29 years), and 34 patients had isthmic spondylolisthesis (18 males and 16 females, mean age of 54.68 ± 13.13 years). The control group included 35 healthy adults (16 males and 19 females, mean age of 64.77 ± 10.47 years,). There was no statistical difference in gender distribution among the three groups (P＞0.05), while the difference in age distribution was significant (P＜0.01). In the degenerative spondylolisthesis group and the isthmic spondylolisthesis group, the TK (29.78°±13.71° and 30.56°±13.61° vs 38.94° ± 10.51°) was lower than those in the healthy control group, PI (54.06°±10.99° and 56.10°±10.76° vs 49.33°±8.92°) and PT (22.84°±8.68° and 21.45°±9.31° vs 12.32°±5.86°) were all higher than those in the healthy control group (P＜0.05). In the degenerative spondylolisthesis group, the LL (40.42°±20.44° vs 49.25°±11.17°) and SS (31.42°± 10.59° vs 37.01°±8.11°) were lower than those in the healthy control group (P＜0.05). The SA in the degenerative spondylolisthesis group was significantly lower than that in the isthmic spondylolisthesis group (P＜0.05). The SK (57.23°±13.29° and 56.01°±11.69° vs 50.35°±9.77°) in both the degenerative spondylolisthesis group and the isthmic spondylolisthesis group was higher than that in the healthy control group (P＜0.05), while the S5Co1HA (80.70°±9.68° and 81.45°±11.56° vs 94.45°±12.02°) was lower than that in the healthy control group (P＜0.05). The S1-2T (17.09°±7.08° vs 13.15°±9.08°) in the degenerative spondylolisthesis group was higher than that in the isthmic spondylolisthesis group (P＜0.05). The multiple linear regression equation for the degenerative spondylolisthesis group was PI=0.383SK-0.753 (S1-2T) +44.989 (R=0.727), and the multiple linear regression equation for the isthmic spondylolisthesis group was PI=0.529SK-0.876 (S1-2T) +38.001 (R=0.815). Conclusion This study demonstrates that patients with both degenerative and isthmic lumbar spondylolisthesis exhibit similar characteristic changes in sacral morphology, including a significantly increased sacral kyphosis angle and a markedly reduced S5Co1 inclination line horizontal angle. These characteristic alterations show significant correlations with spino-pelvic sagittal imbalance, with the combination of sacral kyphosis angle and S1-2 inclination angle effectively predicting pelvic incidence. The findings suggest that sacral morphological parameters can serve as important complementary indicators for assessing spino-pelvic balance in patients with lumbar spondylolisthesis, providing new anatomical basis for developing personalized treatment strategies in clinical practice.

Graphical abstract

关键词

腰椎滑脱 / 骶骨形态参数 / 脊柱-骨盆矢状位参数 / 影像学测量 / 骨盆入射角 / 成年人

Key words

sacrum morphology parameters / lumbar spondylolisthesis / spinal-pelvic sagittal parameters / radiographic measurements / pelvic incidence / adult

腰椎滑脱在人群中的患病率约为6%，由此带来的下腰痛及神经根压迫症状严重影响生活质量。1976年，Wiltse等^[1]根据病因将腰椎滑脱分为5个主要类别：发育不良型、峡部裂型、退行性、创伤性和病理性。其中以峡部裂型和退行性滑脱最为常见。随着外科技术的不断发展，对腰椎滑脱研究的关注点从局部解剖逐渐发展到整体矢状位平衡^[1-5]。滑脱的发生使得脊柱和骨盆的整体对位情况发生变化，一旦超过骨盆代偿时则会出现矢状位的失平衡。既往研究发现在正常人群中骶骨形态参数与脊柱矢状位参数存在相关性，然而对不同类型滑脱的骶骨形态参数了解较少。本研究旨在通过对比分析退行性与峡部裂型腰椎滑脱患者及健康人群的影像学参数，探讨两类滑脱患者的骶骨形态与脊柱-骨盆矢状位参数的特征性改变，并深入探究其内在相关性，以期为揭示不同滑脱类型的内在机制及临床精准诊疗提供影像学依据。

1 对象和方法

1.1 研究对象

收集解放军总医院第一医学中心2019年1月— 2023年12月就诊的腰椎滑脱患者为研究对象。纳入标准：(1)明确诊断为退行性腰椎滑脱或峡部裂型滑脱并伴有腰腿痛症状；(2)具有完整的影像学资料并且清晰可辨认。排除标准：(1)伴有脊柱肿瘤、脊柱骨折、脊柱畸形等其他疾病；(2)既往因脊柱相关疾病行手术治疗；(3)发育不良型、创伤性、病理性腰椎滑脱；(4)多节段腰椎滑脱。

选取同期在解放军总医院第一医学中心健康体检的35例健康成人作为对照组。纳入标准：(1)脊柱站立位全长片影像无异常；(2)无脊柱肿瘤、脊柱骨折、脊柱畸形等其他疾病；(3)年龄为45 ~ 80岁。排除标准：(1)脊柱手术史；(2)下肢关节手术史。

本研究经解放军总医院第一医学中心医学伦理委员会批准(批号：2023-446-02)。

1.2 影像学参数测量

所有患者均拍摄脊柱站立位全长片，测量软件采用Surgimap Spine(Nemaris Inc.，New York)。测量工作由一名脊柱外科医师完成，测量2次，间隔1周后再次测量，取平均值。测量的参数包括两类。

1.2.1 脊柱-骨盆参数

(1)胸椎后凸角(thoracic kyphosis，TK)：胸5椎体上终板切线与胸12椎体下终板切线的夹角；(2)腰椎前凸角(lumbar lordosis，LL)：腰1椎体上终板的切线与骶1椎体上终板的切线的夹角；(3)骨盆入射角(pelvic incidence，PI)：髋轴中点与骶骨平台中点连线与骶骨平台中垂线的夹角；(4)骨盆倾斜角(pelvic tilt，PT)：髋轴中点与骶骨平台中点连线和铅垂线的夹角；(5)骶骨倾斜角(sacral slope，SS)：骶骨平台与水平线的夹角；(6)矢状面轴向距离(sagittal vertical axis，SVA)：C7铅垂线与S1后上角垂线之间的水平距离；(7)滑脱角(slip angle，SA)：滑脱椎体下缘与尾端椎体上缘的夹角。

1.2.2 骶骨参数

(1)骶1-2倾斜角(S1-2 tilt，S1-2T)：Ferguson法标记S1-2的倾斜线，测量其于与S1终板垂线的夹角；(2)骶骨后凸(sacral kyphosis，SK)：S1-2倾斜线与骶5尾1(S5Co1)倾斜线的夹角；(3)S5Co1倾斜线水平角(the horizontal angle of S5Co1，S5Co1HA)：S5Co1倾斜线与水平线夹角。见图1。

1.3 分组及研究方法

根据Wiltse分型将滑脱患者分为退行性滑脱组和峡部裂型滑脱组，比较退行性滑脱组、峡部裂型滑脱组和对照组脊柱-骨盆矢状位参数、骶骨形态参数的差异，观察各组内脊柱-骨盆矢状位参数与骶骨形态参数是否存在相关性。

1.4 统计学分析

采用SPSS 26.0软件进行统计分析。计数资料采用χ²检验，符合正态分布的计量资料以\overline{x}\pm s表示，多组间参数比较采用单因素方差分析；部分偏态资料以中位数M(IQR)表示，两组比较采用Wilcoxon秩检验。组内骶骨参数与脊柱骨盆参数的相关关系采用Pearson检验，并对部分显著指标建立线性回归方程以量化解析其相关关系。P＜0.05为差异有统计学意义。

2 结果

2.1 三组基线资料比较

滑脱组91例，年龄(61.45±11.56)岁。其中退行性滑脱57例，男26例，女31例，年龄(65.49±8.29)岁。峡部裂性滑脱34例，男18例，女16例，年龄(54.68±13.13)岁。Meyerding分型分为：I度滑脱64例，Ⅱ度滑脱24例，Ⅲ度滑脱1例，Ⅳ度滑脱2例。对照组35例，男16例，女19例，年龄(64.77±10.47)岁。三组间性别分布无统计学差异(P＞0.05)，年龄分布存在统计学差异(P＜0.01)。见表1。

2.2 三组脊柱-骨盆矢状位参数比较

退行性滑脱组、峡部裂型滑脱组TK(29.78°± 13.71°、30.56°±13.61°)低于对照组、PI(54.06°± 10.99°、56.10°±10.76°)、PT(22.84°±8.68°、21.45°±

9.31°)均高于对照组。退行性滑脱组的LL(40.42°± 20.44°)、SS(31.42°±10.59°)低于对照组。退行性滑脱组SA(5.70°±4.95°)明显低于峡部裂型滑脱组(P均＜0.05)。见表1。

2.3 三组骶骨参数比较

退行性滑脱组及峡部裂型滑脱组的SK(57.23°±13.29°及56.01°±11.69°)高于对照组，S5Co1HA(80.70°±9.68°及81.45°±11.56°)低于对照组。峡部裂滑脱组的S1-2T(13.15°±9.08°)显著低于退行性滑脱组(P均＜0.05)。见表1。

2.4 两滑脱组脊柱-骨盆矢状位参数与骶骨参数的相关性

经Pearson相关分析，各组内骶骨参数(SK、S1-2T)与矢状位参数(PI)有较好相关性(P均＜0.05)，见表2。进一步行多元线性回归，建立内脊柱-骨盆矢状位参数(PI)与骶骨参数(SK、S1-2T)的线性回归模型(方程)：退行性滑脱组PI=0.383×SK-0.753×(S1-2T)+44.989(R=0.727)，峡部裂型滑脱组PI=0.529×SK-0.876×(S1-2T)+38.001(R=0.815)。见表3。

3 讨论

随着对脊柱整体平衡认识的不断加深，脊柱骨盆矢状位形态在腰椎滑脱的治疗过程中逐渐得到重视^[6-9]。Roussouly等^[10]提出通过外科手术矫正可以改善腰椎滑脱患者的脊柱平衡状态。Elmorsy等^[11]进行了一项纳入40例低中级别退行性或峡部裂型滑脱患者的前瞻性研究，结果显示通过手术恢复正常的矢状位平衡参数可以获得更好的临床结果。因此了解腰椎滑脱的脊柱骨盆参数特征对于外科医师十分重要，本研究发现腰椎滑脱患者的脊柱-骨盆矢状位形态以及骶骨形态与健康人群存在差异。滑脱患者的TK、S5Co1HA显著低于对照组，PI、PT、SK显著高于对照组。退行性滑脱组的LL、SS低于对照组。在两组滑脱患者之间。退行性滑脱组的SA低于峡部裂滑脱组、S1-2T高于峡部裂滑脱组。

Whitesides等^[12]发现，患有峡部裂型滑脱的成人PI平均为64.5°，显著高于正常青少年的41.8°。Duval-Beaupere等^[13]也发现PI在峡部裂型滑脱发生机制中起着重要作用。Funao等^[14]通过对比退行性滑脱与单纯腰椎管狭窄患者的矢状位参数发现退行性滑脱患者具有更大的PI，并且与滑脱程度呈正相关。在本研究中，退行性滑脱及峡部裂型滑脱患者的PI分别为54.06°±10.99°和56.10°±10.76°，显著大于正常人群PI值，并且两种类型滑脱患者的PI值无明显差异，提示两组滑脱患者具有相似的骨盆形态。尽管滑脱组与对照组年龄存在显著差异，但既往研究已经证明成年后PI不会发生改变。Fredrickson等^[15]比较在不同个体的矢状位参数，PI增大时SS的增加明显超过了PT，较大的PI使得骶骨上终板倾斜角度增大^[16]，腰椎向前下方的剪切力增大，脊柱后方结构及椎间盘受力增加^[17]，关节囊和韧带随着年龄的增长而逐渐变得松弛，这反过来又导致关节突关节对前方剪切力的抵抗力大幅下降，由此增大了腰椎滑脱的风险^[18]。这可能是腰椎滑脱发生机制的一种解释。

除此之外本研究发现两组滑脱患者的TK小于对照组，PT值显著大于对照组，退行性滑脱组SS与LL显著小于对照组。在正常人群中，PT的波动范围较窄，PI增大时SS的变化较PT更为显著。然而本研究中滑脱患者具有较大的PI，但SS却没有随之增大，并且在退行性滑脱患者中SS显著小于对照组，这表明在退行性滑脱和峡部裂型滑脱患者中，调动了骨盆的代偿机制。在滑脱发生后，上位椎体相对前移，矢状面整体重心前移，骨盆作为矢状面平衡最重要的调节器，通过其后倾，PT增大，SS减小进一步维持矢状面平衡^[1]。在Roussouly对人群脊柱骨盆矢状位分型中，PI较大对应较大的TK和LL。本研究中两组滑脱患者TK显著小于对照组，提示为维持脊柱矢状面平衡，减少重心前倾，胸椎后凸代偿性减小。同时LL小于对照组，提示在腰椎滑脱发生过程中，LL为适应骨盆的后倾，同样发生了代偿性变化。

本研究中SA在退行性滑脱和峡部裂型滑脱患者之间存在统计学差异。退行性腰椎滑脱患者椎间盘退变更加显著，椎间隙高度明显下降。椎间盘前缘高度的丢失会使退变节段的腰前凸角减小，同时滑脱的发生会使SA进一步减小。而峡部裂型滑脱没有了后方结构的限制，在滑脱发生后，局部SA角度增大，以维持正常的矢状面平衡。这可能是退行性滑脱SA小于峡部裂型滑脱的原因之一。SA角度在维持下腰椎区域性前凸和矢状位整体平衡有重要作用。Lee等^[19]在一项研究中发现，腰椎融合术后SA减小的患者较SA增加的患者矢状面恶化的风险更高，这提示我们在对退行性滑脱患者制订手术方案时应更加注意对SA角度的纠正，维持下腰椎的局部前凸，以获得更好的手术效果^[20]。

骶骨作为脊柱最下方的延续，同时也是骨盆的重要组成部分。骨盆的前倾与后倾通过骶骨向上传递，调节人体的矢状位形态。腰椎滑脱的实质为下腰椎与骶骨之间的解剖异常^[21]，从而发生整体矢状位对线变化。既往研究观察到了骶骨形态在腰椎滑脱发生发展中的作用^[22]。有研究认为，STA较小的人群中峡部裂型滑脱患病率更高^[12]。Zeliha^[23]在退行性滑脱患者中发现了相同的表现，认为较小的STA增加了滑脱的倾向。然而Kazuhiro等^[24]发现晚期峡部裂型滑脱患者的STA偏低，而在早期和进展期峡部裂型滑脱患者中并未观察到相同的特征，推测较低的STA来自滑脱发生后剪切力对骶骨终板的重塑。Wang等^[25]在一项纳入131例发育性滑脱症患者的研究中发现，滑脱患者的STA显著低于对照组，并且在两组人群中STA与PI都存在相关性，同时还发现了滑脱患者具有更大的SK，推测SK较高是椎体滑移的易感因素。Hresko等^[21]的研究中发现滑脱患者的S1-S2和S2-S3之间的后凸角度更大。可见骶骨形态与滑脱的发生发展存在相关性。宋凯等^[26]在正常成年人群中发现骶骨形态与脊柱骨盆矢状位形态存在显著相关性，并提出了一系列新的骶骨形态参数。在本研究中，对退行性滑脱及峡部裂型滑脱患者的骶骨形态进行了比较分析，发现退行性滑脱和峡部裂型滑脱患者的SK分别为57.23°±13.29°和56.01°±11.69°，显著大于正常人群的50.35°± 9.77°，具有统计学差异，而两组滑脱患者的SK无明显差异。既往研究证明，不仅是胸椎及腰椎的曲度，骶骨曲度的形成也与人类直立行走有明显的关联，并且骶骨曲度与腰椎前凸有一定相关性^[27]。通过研究青少年骶骨形态的特征发现，在发育期SK与LL之间也存在相关性，而PI与LL无显著相关性^[28]。因此可以推测滑脱患者具有较大的SK，对应较大的LL，造成下腰部受力增加，可能是腰椎滑脱发生的因素之一。S5Co1HA是一个位置参数^[26]，在两组滑脱患者中的S5Co1HA值都显著减小，提示骨盆发生了明显的后旋代偿。

峡部裂型滑脱组的S1-2T小于退行性滑脱组和对照组。与SS不同的是，S1-2T是一个解剖参数，S1-2T越小代表骶1上终板相对更加垂直于地面。在峡部裂型滑脱中，骶骨参数的特征是SK较大，而S1-2T较小，即代表较为弯曲的脊柱形态以及高耸的骶骨后上角。由此在脊柱伸展时，后方关节突关节会受到L4及S1关节突的钳夹和撞击。这是对峡部裂型滑脱发生的一种推测。

既往研究证明了在正常人群中骶骨参数与脊柱骨盆矢状位参数之间存在的相关性^[26]。本研究在滑脱患者中也发现了这种规律。在滑脱患者中SK与S1-2T两个骶骨形态参数与PI均有相关性，通过多元线性回归分析，在退行性滑脱组PI=0.383SK-0.753(S1-2T)+44.989(R=0.727)，峡部裂型滑脱组PI=0.529SK-0.876(S1-2T)+38.001(R=0.815)。骶骨矢状面形态是每个个体特有的解剖变量，在成年后不易发生改变^[21]。骨盆入射角作为脊柱矢状位形态测量的经典参数，在评估脊柱平衡、手术方案制订方面具有重要意义。本研究发现了骶骨参数与骨盆入射角的相关性，可在股骨头或髋臼无法良好识别时作为一种补充测量方法，通过上述方程进行预测PI等参数，可以帮助临床医师制定更精确的治疗计划，提高手术成功率，改善患者的长期预后^[29-31]。

本研究存在一定的局限性。首先，由于峡部裂型滑脱患者发病年龄小于退行性滑脱患者，所选观察对象年龄不在同一基线水平，不能排除年龄因素对研究结果造成影响。其次，所纳入对象中大部分为低度滑脱患者，可能形成选择偏倚，对骶骨形态与滑脱程度的相关性有待进一步分析。

综上所述，在腰椎滑脱患者与正常人群之间，骶骨形态与脊柱骨盆矢状位参数存在显著差异，除SA，S1-2T外，在退行性滑脱与峡部裂型滑脱之间未观察到显著差异。滑脱患者的骶骨形态参数与脊柱

基金资助

国家自然科学基金项目(82372499)