2.1   文献筛选流程及结果   共从数据库中检索出429篇文献，其中PubMed(n=50)、Cochrane Library (n=2)、Web of Science(n=93)、Embase (n=120)、中国知网(n=19)、万方(n=113)、维普(n=9)及中国生物医学文献数据库(n=23)。最终纳入17篇文献，其中英文文献7篇、中文文献10篇。文献筛选流程及结果见图1。"

2.2   纳入研究的基本特征和质量评价   此次研究共纳入17篇文献，总样本量为2 156例，其中轴性症状组783例，非轴性症状组1 373例，轴性症状的总发生率为34.88%。文献的基本特征见表1。"

2.3   Meta分析结果 2.3.1   颈后路椎板成形术后轴性症状的发生率   共有17项研究报告了颈后路椎板成形术后轴性症状的发生率[5，15，17，23-36]。结果显示，研究间存在显著异质性，为了进一步探讨异质性的来源，基于文献中是否按照1∶1比例纳入对照组，进行亚组分析，然而，由于1∶1对照组中研究仅为1项，异质性结果仍然较高。此外，由于1∶1对照组在选取患者时严格遵循了与轴性症状组等比例的原则，这可能导致其结果无法准确反映实际发生率。因此，最终分析中将该亚组予以排除。为评估结果的稳定性，进行了敏感性分析，结果显示剔除任意单项研究后，结果均无显著变化。分析结果显示，颈后路椎板成形术后轴性症状的发生率为34.88%(95%CI：0.301 6-0.399 0)，结果见图2。 2.3.2   年龄   共有4项研究纳入了年龄这一变量[24，30，35-36]。异质性检验结果显示，I2=89%，P < 0.000 1，提示研究间存在较高的异质性。为进一步探讨异质性的来源，采用逐一排除法对单项研究进行敏感性分析，结果显示剔除任一研究后异质性均未发生显著变化，因此，最终采用随机效应模型对效应值进行合并分析。Meta分析结果显示，年龄与颈后路椎板成形术后轴性症状的发生无显著统计学关联(OR=1.07，95%CI：0.95-1.21，P=0.25)，结果见图3。 "

2.3.3   性别   共有4项研究纳入了性别这一变量[23-24，31-32]。异质性检验结果显示，I2=44.2%，P=0.146 5，提示研究间异质性较低。为验证结果的稳健性，采用逐一排除法进行敏感性分析，结果表明剔除任一研究后异质性未出现显著变化，因此，最终采用固定效应模型对效应值进行合并分析。Meta分析结果显示，性别与颈后路椎板成形术后轴性症状的发生无显著统计学关联(OR=1.18，95%CI：0.78-1.78，P=0.68)，结果见图4。"

2.3.4   术前存在颈肩部疼痛   共有6项研究纳入了术前存在颈肩部疼痛这一变量[23，25-27，31，34]。异质性检验结果显示，I2=0%，P=0.471 5，提示研究间无显著异质性。为了进一步验证结果的稳定性，采用逐一排除法对单项研究进行敏感性分析，结果表明剔除任一研究后异质性均未发生显著变化，因此，最终采用固定效应模型对效应值进行合并。Meta分析结果显示，术前存在颈肩部疼痛是颈后路椎板成形术后轴性症状的独立危险因素(OR=2.35，95%CI：2.20-2.50，P < 0.000 1)，结果见图5。 2.3.5   术前颈椎曲度小于12°   共有3项研究纳入了术前颈椎曲度小于12°这一变量[17，28-29]。异质性检验结果显示，I2=0%，P=0.946 6，提示研究间无显著异质性。为验证分析结果的稳健性，采用逐一排除法对单项研究进行敏感性分析，结果显示剔除任一研究后异质性均未发生显著变化，因此，最终采用固定效应模型对效应值进行合并分析。Meta分析结果显示，术前颈椎曲度小于12°是颈后路椎板成形术后轴性症状的独立危险因素(OR=3.16，95%CI：1.91-5.24，P < 0.000 1)，结果见图6。 "

2.3.6   术前颈椎椎管占位率   共有3项研究纳入了术前颈椎椎管占位率这一变量[5，15，33]。异质性检验结果显示，I2=74.2%，P=0.020 8，提示研究间存在较高异质性。为探讨异质性的潜在来源，采用逐一排除法进行敏感性分析，结果显示剔除任一研究后异质性未发生显著变化，分析结果稳定，因此，最终采用随机效应模型对效应值进行合并分析。Meta分析结果表明，术前颈椎椎管占位率是颈后路椎板成形术后轴性症状的独立危险因素(OR=1.35，95%CI：1.02-1.79，P=0.034 0)，结果见图7。 2.3.7   术中发生关节突关节损伤   共有3项研究纳入了术中发生关节突关节损伤这一变量的分析[24，28-29]。异质性检验结果显示，I2=0.00%，P=0.416 1，提示研究间无显著异质性。为验证分析结果的稳健性，采用逐一排除法对单项研究进行敏感性分析，结果表明剔除任一研究后异质性未发生显著变化，分析结果仍保持稳定，因此，采用固定效应模型对效应值进行合并分析。Meta分析结果显示，术中发生关节突关节损伤是颈后路椎板成形术后轴性症状的独立危险因素(OR=2.87，95%CI：1.96-4.19，P < 0.000 1)，结果见图8。 2.3.8   术中关节突关节损伤数量   共有3项研究纳入了术中关节突关节损伤数量这一变量的分析[5，27，34]。异质性检验结果显示，I2=89.4%，P < 0.000 1，提示研究间存在较高异质性。为探讨异质性的潜在来源，进行了逐一排除法的敏感性分析，结果表明剔除任一研究后异质性未发生显著变化，分析结果稳定，因此，采用随机效应模型对效应值进行合并分析。Meta分析结果显示，术中关节突关节损伤数量与颈后路椎板成形术后轴性症状的发生无显著统计学关联(OR=1.00，95%CI：0.57-1.75，P=0.991)，结果见图9。 "

2.3.9   术中C7棘突肌肉止点损伤   共有3项研究纳入了术中C7棘突肌肉止点损伤这一变量的分析[5，28，35]。异质性检验结果显示，I2=82.1%，P=0.003 8，提示研究间存在较高异质性。为进一步探讨异质性的潜在来源，采用逐一排除法对单项研究进行敏感性分析，结果表明剔除任一研究后异质性未发生显著变化，分析结果稳定，因此，最终采用随机效应模型对效应值进行合并分析。Meta分析结果显示，术中C7棘突肌肉止点损伤与颈后路椎板成形术后轴性症状的发生无显著统计学关联(OR=1.16，95%CI：0.29-4.59，P=0.832 8)，结果见图10。 2.3.10   颈椎曲度变化   共有4项研究纳入了该变量的分析[5，27，31-32]，初步结果显示异质性较大。为明确异质性来源，采用逐一排除法进行敏感性分析，结果表明剔除杨明坤等[31]的研究后，异质性显著降低(I2=0%，P=0.665 5)，因此，选用固定效应模型对效应值进行合并分析。分析结果显示，颈椎曲度变化是颈后路椎板成形术后轴性症状的独立危险因素(OR=1.10，95%CI=1.03-1.17，P= 0.002 9)。结果见图11。 2.3.11   颈椎活动度变化   共有3项研究纳入了颈椎活动度变化这一变量的分析[5，27，34]。异质性检验结果显示，I2=0%，P=0.480 3，提示研究间无显著异质性。通过逐一排除法进行敏感性分析，结果表明剔除任一研究后异质性未发生显著变化，分析结果稳定，因此，采用固定效应模型对效应值进行合并分析。Meta分析结果显示，颈椎活动度变化是颈后路椎板成形术后轴性症状的独立危险因素(OR=1.09，95%CI：1.04-1.14，P=0.000 5)，结果见图12。 "

2.3.12   术后早期功能锻炼   共有3项研究纳入了术后早期功能锻炼这一变量的分析[5，32，36]。异质性检验结果显示，I2=88.7%，P=0.000 1，提示研究间存在较高异质性。通过逐一排除法进行敏感性分析，结果表明剔除任一研究后异质性并未显著降低，分析结果稳定，因此，采用随机效应模型对效应值进行合并分析。Meta分析结果显示，术后早期功能锻炼与颈后路椎板成形术后轴性症状的发生无显著统计学关联(OR=0.71，95%CI：0.20-2.56，P=0.597 4)，结果见图13。"

"

2.4   发表偏倚分析   在纳入的17项研究中，对颈后路椎板成形术后轴性症状的发生率进行发表偏倚分析，结果显示，Egger’s检验结果为P=0.048 6< 0.05，提示无显著的发表偏倚；而Begg’s检验结果为P=0.165 3> 0.05，提示无显著的发表偏倚。进一步结合漏斗图分析发现，除2项研究效应量分布相对离散外，其余研究的效应量分布基本对称，提示存在发表偏倚的可能性较低，结果见图14。此外，通过逐一剔除研究进行敏感性分析，发现结果无显著改变，进一步表明结果具有高度稳定性。结果见图15。"