文章 | 出版前文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 《科学报告》本出版前文章提前共享,以便您更快获取新研究。它可被引用,并具有永久DOI。最终编辑版本将自动替换本文。隐形矫治器治疗中与各种正畸牙齿移动相关的外根吸收量化:一项回顾性临床研究 Safa Kurnaz · Muhammed Hilmi Buyukcavus 收稿日期:2026年5月20日 / 接受日期:2026年10月19日 © 作者 2026 摘要 本研究探讨了隐形矫治器(CA)正畸治疗(OT)期间牙齿移动与外根吸收(ERR)之间的关联。纳入48例使用ClearCorrect矫治器接受非拔牙正畸治疗的患者。评估上颌和下颌切牙及尖牙,共576颗牙齿。将正畸治疗前后的数字模型进行重叠,以量化压低、伸长、牙冠倾斜、牙冠转矩、旋转、颊舌向平移和近远中向平移。在正畸治疗前后的全景放射片上测量根长度、投影根表面积(PRSA)和ERR。使用相关性和多变量线性回归分析评估牙齿移动与根长度、PRSA和ERR之间的关联,并在适当情况下进行多重比较校正。在上颌侧切牙和尖牙中观察到显著的移动特异性关联,其中牙冠转矩与根长度和ERR变化相关,牙冠倾斜与侧切牙中的ERR相关,而伸长与尖牙中的PRSA变化相关。在下颌牙齿中未观察到显著的移动-吸收关联。观察到的吸收变化幅度通常有限。总体而言,研究结果表明,CA治疗期间的根反应可能因牙齿移动类型和牙组而异;然而,观察到的移动特异性关联应谨慎解读。关键词 隐形矫治器 · 牙髓病学 · 外根吸收 · 正畸学 · 牙齿移动 缩写 ERR 外根吸收 OT 正畸治疗 CA 隐形矫治器 ABO 美国正畸委员会 SD 标准差 CEJ 釉牙骨质界 CBCT 锥形束计算机断层扫描 PRSA 投影根表面积 FDR 错误发现率
预出版文章《Scientific Reports》 https://doi.org/10.1038/s41598-026-77321-y 引言 正畸治疗(OT)可导致一种常见副作用,即外部牙根吸收(ERR),造成根尖牙根组织不可逆丧失1。虽然OT通常可改善美观和功能,且在大多数病例中ERR无临床意义,但严重病例可损害牙齿结构并可能导致牙齿丧失2,3。多种因素可影响OT期间ERR的风险。这些因素包括遗传、年龄、性别和全身健康状况,以及治疗相关因素,如错𬌗畸形的严重程度、牙根形态和正畸力持续时间及强度。治疗方法,包括拔牙、矫治器类型和移动幅度,也可影响ERR4,5。尽管任何牙齿在OT期间均可发生ERR,但上颌和下颌切牙的风险通常较高6。虽然多种因素与ERR有关,但其促进该过程的确切机制尚未完全明确7。固定矫治器治疗长期以来一直是OT的标准。然而,隐形矫治器(CAs)因舒适性、美观性改善和治疗时间缩短而日益流行8,9。尽管证据仍存在异质性,但近期系统综述和荟萃分析提示,CAs的ERR发生率和严重程度通常低于固定矫治器10–14。间歇性施力、较轻的力和较短的治疗时间等因素可能有助于降低该风险15,16。OT施加多方向力以实现牙齿移动,可能破坏根尖处骨吸收与骨形成之间的平衡。过度或长时间的力可超过牙周组织的适应能力,导致牙根结构进行性丧失17,18。正畸力的大小和方向均可影响ERR风险19,20。压低力比伸长力更易诱发ERR21。此外,一项动物研究报告,与整体牙齿移动相比,倾斜牙齿移动与更大的ERR相关22。此外,人体研究也显示,牙根转矩增加、倾斜和旋转与更高的ERR相关23–25。然而,难以分离出导致ERR的特定牙齿移动。这是因为正畸矫治器通常产生复杂移动的组合,包括伸长、压低、平移、倾斜、转矩和旋转20。许多临床医生认为ERR是OT期间可潜在预防的并发症,并常将ERR的发生归责于正畸医生。然而,关于牙髓病学与正畸学在治疗决策中的总体关系的详细信息仍有限26。因此,进一步研究导致ERR的因素对于改善治疗结果和尽量减少不良反应至关重要。全面的研究可确定不同正畸牙齿移动和位移类型及幅度如何影响ERR量27。近期系统综述和荟萃分析主要集中于比较CAs与固定矫治器之间ERR的发生率和严重程度,通常提示CA治疗的ERR风险较低12–14。然而,这些综述也强调了现有证据的异质性、方法学局限性以及需要进一步高质量临床研究28。重要的是,既往研究在很大程度上将ERR作为总体治疗结果进行评估,而非研究CA治疗期间单个正畸牙齿移动对ERR的具体贡献。据我们所知,尚无既往临床研究全面评估仅接受CAs治疗的患者中,不同正畸牙齿移动(包括压低、伸长、平移、牙冠倾斜、牙冠转矩和旋转)与不同牙组ERR幅度之间的关联。因此,本研究旨在探讨不同正畸牙齿移动与ERR之间的关联,并量化CA治疗后不同牙组中每种移动类型相关的ERR幅度,从而提供关于移动特异性吸收模式的临床相关证据。本研究的零假设(H0)如下:不同正畸牙齿移动在不同牙组中导致的ERR量无差异。
预出版文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 方法 本回顾性临床研究经当地大学伦理委员会批准(参考编号:2024/07-36)。所有参与者均提供了知情同意。对于未满18岁的参与者,已从其法定监护人处获得同意。针对包含七个预测变量的多变量线性回归分析进行了敏感性功效分析。在样本量为48、α水平为0.05、功效为80%的情况下,分析表明可检测的最小总体效应量为f²=0.35。对2022年7月至2024年12月期间在本诊所接受治疗患者的放射影像进行了回顾性分析。根据预先设定的纳入和排除标准筛选患者。制定纳入标准旨在获得临床同质的研究人群,并尽量减少可能影响ERR的潜在混杂因素。因此,仅纳入具有适合准确识别标志点的高质量标准化全景放射影像、恒牙列完整无缺牙(第三磨牙除外)且接受CA联合OT治疗的患者。患者还须具有骨性I类错𬌗畸形、中度拥挤,且治疗计划不涉及拔牙。此外,仅纳入无副功能习惯(如磨牙症或紧咬牙)、使用同一CA系统治疗且计划牙齿移动程度相似的患者10,11,29–31。应用排除标准旨在尽量减少可能影响ERR的生物学和治疗相关因素。若患者存在牙齿或牙根异常(包括未萌出牙,无论受累牙为何),则予以排除。有牙外伤史、既往OT史、记录缺失或拔牙史的患者也被排除。其他排除标准包括全身性疾病、颅面畸形、颞下颌关节紊乱、口腔健康状况差(龋齿或牙周病)、根管充填牙、根尖孔开放以及既存牙根吸收10,11,29–31。为进一步减少与病例复杂性相关的变异性,采用美国正畸委员会(ABO)差异指数评估每个病例的难度。排除高复杂度评分的患者,以保持研究人群更为同质30,32。治疗方案 对235例患者应用纳入和排除标准后,共纳入48例接受CA治疗的患者(26例女性,22例男性)。参与者的平均年龄为19.22±2.08岁。对每位患者的上颌和下颌中切牙、侧切牙和尖牙分别进行评估。因此,每个牙组分析96颗牙,共计576颗牙。治疗组由使用ClearCorrect CA(ClearCorrect®,ClearCorrect LLC,美国德克萨斯州Round Rock)接受OT的患者组成。使用ClearCorrect的ClearPilot™软件(ClearCorrect治疗计划工具,版本5)创建虚拟治疗计划以指导治疗过程。该计划概述了一系列治疗步骤,包括邻面去釉、附件放置以及颌间弹力圈的使用。指导患者每15天更换一次矫治器,每日佩戴约22小时。牙齿移动的测量 所有数字模型测量均采用标准化工作流程进行,以确保测量的一致性和可重复性。所有治疗前后模型均使用相同的参考标志点、感兴趣区域、配准程序和测量工具。使用3Shape TRIOS 3(3Shape TRIOS A/S,丹麦哥本哈根)口内扫描仪采集数字印模,并上传至ClearPilot平台进行治疗设计。所得模型的STL文件使用3D Slicer 5.0保存。模型重叠按照Garib等33所述的基于腭皱襞的方法进行,使用稳定的腭部参考区域进行配准。标志点置于切牙乳头后部
预出版文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 边界、第二和第三腭皱襞的内侧缘,以及距第三腭皱襞内侧缘远端10 mm处。围绕这些标志点定义了感兴趣区域。切牙乳头和第三腭皱襞投影的感兴趣区域半径为5个顶点;第二和第三腭皱襞的感兴趣区域半径为30,以确保感兴趣区域与牙齿保持距离。使用Slicer的“CMF配准”工具,通过感兴趣区域匹配将治疗后上颌骨叠加到治疗前上颌骨上。随后将下颌模型与叠加后的上颌骨咬合对齐。将上颌叠加矩阵应用于下颌模型,使其与上颌骨对齐。最后,在所有模型上放置预定义标志点,以量化治疗期间实现的正畸牙齿移动幅度(图1)。在上颌和下颌模型上放置标志点以测量:(1) 水平/垂直线性位移(尖牙颊尖和切牙近中切缘),(2) 近远中旋转(咬合面观的最远中/最近中点),以及(3) 颊舌向倾斜(切牙切端中心和最深龈点)。使用3D Slicer的基准点配准向导,记录治疗前/后和初始/最终预测模型的线性位移标志点的x、y和z坐标。坐标差异量化了实现的和预测的变化。Q3DC测量了每个模型对的近中/远中标志点之间(近远中旋转)以及切端/颈部标志点之间(颊舌向倾斜)的角度。对于评估的12颗牙齿(上颌中切牙、侧切牙和尖牙,以及下颌中切牙、侧切牙和尖牙),从系统获取相关牙齿从OT开始到结束(包括精调矫治器)发生的总压低、伸长、颊舌向平移、近远中平移、牙冠倾斜、牙冠转矩和旋转量,以度数和/或毫米为单位。外部牙根吸收的测量方法 为尽量减少测量偏倚,所有放射学测量均在相同成像和测量条件下使用标准化测量方案进行。使用OT前和OT后立即获得的数字全景X线片评估所有牙齿的ERR程度。所有X线片均在自然头位下拍摄以确保一致性。这是通过指示患者直视设备内置镜子实现自我对齐来完成的。使用Castellini X射线机(Castellini X Radius Compact,意大利伊莫拉)获取图像,并通过其集成定位系统辅助头部稳定。牙冠和牙根长度测量使用Castellini X射线机的测量软件,按照先前研究中描述的经过验证的方法进行11,29,31,34。测量从标准化的治疗前和治疗后全景X线片获得,所有牙齿使用相同的解剖参考点。线性和PRSA测量所使用的参考点和程序如图2所示。图1 使用腭皱襞法对治疗前和治疗后模型进行叠加。
预出版文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y 牙冠和牙根长度的测量分多个步骤进行。首先,通过连接近中与远中釉牙骨质界(CEJ)点画一条直线来确定釉牙骨质界。然后在初始和最终X线片上分别测量牙冠长度。该测量定义为从切缘到釉牙骨质界的最长距离。对所有牙齿的牙根长度测量采用类似方法,记录从釉牙骨质界到根尖的距离11,29,31,34。ERR按如下公式计算:ERR(mm)=(正畸前牙根长度 − 正畸后牙根长度)×(正畸前牙冠长度/正畸后牙冠长度)。牙冠长度比用作比例校正因子,以补偿正畸前后全景片之间放射学放大率的差异11,29,31,34。因此,ERR以毫米表示。为补充线性测量,从全景放射照片中分析投影根表面积(PRSA)(图2)。在进行PRSA测量之前,每张全景图像在SketchAndCalc中通过在全景放射图像上定义一段已知线性距离并输入其对应的毫米长度来进行校准,从而建立像素-毫米比例尺。然后使用SketchAndCalc软件(Axiom Welldone,https://www.sketchandcalc.com/)进行根表面积测量。沿牙根的近中和远中轮廓,从釉牙骨质界开始延伸至放射学根尖,仔细描记整个根表面。为确保测量精度,沿描记的根边界以1像素间隔手动放置数字点。根据根长度和弯曲度,每个牙根平均使用25(±5)个点。软件随后自动计算封闭面积(以平方毫米为单位),即代表PRSA11,31。对正畸前和正畸后放射照片采用相同的校准和描记程序。统计分析 使用Shapiro–Wilk检验评估数据正态性。由于数据呈正态分布,进行参数统计分析。在主要分析之前评估测量可靠性。所有测量均由一名经过校准的研究者完成。在主要研究之前,使用随机选取的16名患者全景放射照片评估可靠性和可重复性。同一研究者在间隔两周后重复测量,第二名研究者独立重复相同过程。图2 全景放射照片上牙冠长度、牙根长度和PRSA的参考点及测量程序。
预出版文章《科学报告》 https://doi.org/10.1038/s41598-026-77321-y 检查者内部和检查者间一致性的组内相关系数分别为0.953(95% CI 0.925–0.981)和0.929(95% CI 0.902–0.956),表明测量重现性极佳。描述性统计以均值 ± 标准差(SD)、频数(n)和百分比(%)表示。采用Pearson相关分析评估七项正畸牙齿移动变量(压低、伸长、颊舌向平移、近远中向平移、牙冠倾斜、牙冠转矩和旋转)与根据牙齿位置变化的根长、PRSA和ERR之间的关联。为校正多重检验,采用Benjamini–Hochberg错误发现率(FDR)程序对相关性p值进行调整,FDR校正后p值 < 0.05视为具有统计学意义。对FDR校正后至少仍有一项关联保持显著的牙齿组/结局组合进行后续多变量线性回归分析。所有七项牙齿移动变量同时纳入,并校正侧别及采用患者聚类稳健标准误以考虑患者内依赖性。每个回归模型纳入来自48例患者的96颗牙齿,结果报告为非标准化系数(B)、患者聚类稳健标准误、标准化系数(β)、95%置信区间、t统计量、p值、R²、校正R²以及七项牙齿移动预测变量的联合F统计量。这些后续回归分析被视为探索性分析。所有统计分析均使用IBM SPSS Statistics for Windows 26版(IBM Corp., Armonk, NY, USA)完成。结果 表1呈现了研究参与者的描述性统计,包括年龄、性别、矫治器数量、治疗时长以及基于ABO标准的治疗复杂度评分。CA治疗导致患者前牙牙齿移动的程度见表2。图3以图形方式概述了不同正畸牙齿移动下牙根吸收参数的分布。OT后根长、PRSA和ERR的变化见表3。牙根吸收参数与牙齿移动程度之间的相关分析结果分别见表4(上颌牙)和表5(下颌牙)。经Benjamini–Hochberg FDR校正后,仅上颌牙仍存在统计学显著的相关性。多变量线性回归分析结果见表6。对于上颌侧切牙,FDR校正后牙冠转矩与根长变化呈负相关(r = − 0.504,BH-p = 0.016),与ERR呈正相关(r = 0.527,BH-p = 0.010),而牙冠倾斜与ERR呈正相关(r = 0.555,BH-p = 0.006)。在多变量回归分析中,牙冠转矩每增加1度,根长减少约0.07 mm(B = − 0.070,95% CI −0.136至− 0.003,p = 0.041),ERR增加约0.09 mm(B = 0.093,95% CI 0.025至0.160,p = 0.008)。在校正其他牙齿移动变量和侧别后,旋转也与根长变化(B = 0.032,95% CI 0.006至0.058,p = 0.018)和ERR(B = − 0.030,95% CI −0.060至− 0.001,p = 0.044)相关。特征 数量(n) 性别 女性26 男性22 特征 均值 ± SD 实足年龄(岁)19.22 ± 2.08 ABO差异指数12.06 ± 1.91 矫治器总数32 ± 12 精调次数1.81 ± 0.75 治疗时长(年)1.31 ± 0.47 表1 研究参与者的描述性统计。ABO,美国正畸委员会;n,数量;SD,标准差。
预出版文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 对于上颌尖牙,在FDR校正后,伸长与PRSA变化呈负相关(r = − 0.552,BH-p = 0.006)。在多变量回归分析中,伸长每增加1 mm,PRSA的减少量约增加6.29 mm²(B = − 6.292,95% CI −9.385至− 3.198,p < 0.001)。讨论 尽管CA的使用日益增多,但其对ERR的影响仍研究不足7。既往研究一致表明,接受CA治疗的患者中ERR的患病率低于接受固定矫治器治疗的患者10,11,30,35。近期的系统综述和伞形综述得出了类似结论,同时也强调了现有证据的异质性、方法学局限性以及需谨慎解读的必要性13,14。然而,这些研究主要比较了不同矫治器系统之间ERR的总体发生率或严重程度,而非评估特定正畸牙齿移动对不同牙组ERR的影响。表2 前牙正畸牙齿移动的描述性统计。牙齿移动 压低(均值 ± SD,最小值–最大值,mm) 伸长(均值 ± SD,最小值–最大值,mm) 平移(颊舌向)(均值 ± SD,最小值–最大值,mm) 平移(近远中向)(均值 ± SD,最小值–最大值,mm) 牙冠倾斜(近远中向)(均值 ± SD,最小值–最大值,°) 牙冠转矩(颊舌向)(均值 ± SD,最小值–最大值,°) 旋转(均值 ± SD,最小值–最大值,°) 牙位 上颌 11 0.20 ± 0.34(0.0–1.3) 0.19 ± 0.39(0.0–1.6) 0.80 ± 0.65(0.0–2.1) 0.55 ± 0.47(0.0–1.8) 1.85 ± 1.78(0.0–6.2) 4.36 ± 3.58(0.0–13.2) 9.15 ± 7.88(0.0–30.1) 12 0.28 ± 0.49(0.0–1.6) 0.37 ± 0.52(0.0–1.9) 0.82 ± 0.59(0.0–2.3) 0.64 ± 0.58(0.0–2.1) 2.95 ± 4.14(0.0–10.7) 3.34 ± 3.11(0.0–10.1) 8.29 ± 9.59(0.0–34.4) 13 0.29 ± 0.45(0.0–1.4) 0.13 ± 0.27(0.0–1.1) 0.70 ± 0.70(0.0–2.5) 0.54 ± 0.74(0.0–2.9) 2.23 ± 3.21(0.0–11.3) 3.69 ± 3.91(0.0–17.1) 6.90 ± 8.13(0.0–30.2) 21 0.21 ± 0.31(0.0–0.9) 0.25 ± 0.38(0.0–1.4) 0.90 ± 0.50(0.2–2.2) 0.60 ± 0.64(0.0–2.7) 2.91 ± 3.32(0.0–11.1) 4.38 ± 5.18(0.0–19.1) 7.12 ± 4.64(0.0–16.1) 22 0.21 ± 0.38(0.0–1.4) 0.42 ± 0.52(0.0–1.5) 0.60 ± 0.56(0.0–2.6) 0.42 ± 0.44(0.0–1.6) 1.65 ± 2.69(0.0–12.2) 3.89 ± 3.65(0.0–10.3) 9.68 ± 8.93(0.0–34.5) 23 0.14 ± 0.24(0.0–0.9) 0.22 ± 0.36(0.0–1.3) 0.62 ± 0.55(0.0–1.8) 0.57 ± 0.64(0.0–3.1) 1.26 ± 2.10(0.0–8.1) 2.09 ± 2.47(0.0–8.3) 9.21 ± 6.58(0.0–25.1) 下颌 31 0.55 ± 0.64(0.0–1.7) 0.07 ± 0.14(0.0–0.6) 0.99 ± 0.64(0.0–2.4) 0.68 ± 0.54(0.0–1.8) 2.52 ± 2.14(0.0–8.1) 4.64 ± 5.02(0.0–18.1) 9.45 ± 7.77(0.0–28.1) 32 0.64 ± 0.60(0.0–1.8) 0.03 ± 0.08(0.0–0.3) 0.92 ± 0.74(0.0–2.6) 0.88 ± 0.76(0.0–2.8) 3.85 ± 3.27(0.0–9.1) 2.69 ± 3.24(0.0–12.1) 8.88 ± 8.64(0.0–34.5) 33 0.49 ± 0.49(0.0–2.2) 0.05 ± 0.16(0.0–0.6) 0.74 ± 0.73(0.0–2.9) 0.81 ± 0.67(0.0–2.3) 2.46 ± 3.25(0.0–10.3) 2.70 ± 2.32(0.0–7.5) 9.07 ± 6.88(0.0–26.3) 41 0.59 ± 0.65(0.0–2.6) 0.07 ± 0.18(0.0–0.7) 0.89 ± 0.67(0.0–2.3) 0.87 ± 0.83(0.0–3.3) 2.52 ± 2.53(0.0–8.1) 4.70 ± 5.50(0.0–20.1) 7.45 ± 6.14(0.0–27.2) 42 0.68 ± 0.74(0.0–3.2) 0.08 ± 0.21(0.0–0.7) 1.03 ± 0.85(0.0–2.9) 0.72 ± 0.63(0.0–2.4) 2.16 ± 2.57(0.0–9.2) 3.97 ± 4.30(0.0–13.4) 6.77 ± 5.64(0.0–17.2) 43 0.49 ± 0.44(0.0–1.7) 0.04 ± 0.13(0.0–0.5) 0.82 ± 0.85(0.0–2.9) 0.62 ± 0.47(0.0–1.7) 1.98 ± 2.96(0.0–8.5) 2.35 ± 2.89(0.0–9.1) 8.69 ± 5.15(0.0–16.5) SD,标准差;Min–Max,最小值–最大值;mm,毫米;°,度。图3 各种牙齿移动下牙根吸收参数点分布的图形表示。
预出版文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y 因此,本研究旨在探讨不同牙组中特定正畸牙齿移动与牙根吸收参数之间的关联,并量化接受隐形矫治器(CA)治疗的患者中与这些移动相关的吸收性变化。本研究结果表明,牙齿移动和牙位特异性的关联主要见于上颌前牙的牙根吸收参数。特别是在上颌侧切牙的牙冠转矩和上颌尖牙的伸长方面观察到显著关联,探索性多变量分析中也观察到相关模式。正畸力的大小、持续时间和连续性是无菌性牙根吸收(ERR)的关键决定因素36,37。较大且连续的力与ERR的患病率和严重程度增加相关4,38。相反,间歇性力可能通过促进牙骨质修复来减轻吸收37,39。CA施加的力的大小和方向与托槽和弓丝产生的力不同。据认为,CA传递较轻的间歇性力,与固定矫治器相比,可能有助于降低ERR的发生率40,41。这一解释得到近期综述的支持,这些综述提示CA的间歇性力模式可能允许牙周和牙骨质修复期,尽管临床证据仍存在异质性14,28。很少有研究评估与正畸牙齿移动相关的ERR程度,且方法学上的巨大差异限制了直接比较。本研究首次专门在CA治疗背景下探讨这一现象。虽然Phillips42观察到ERR与上颌中切牙的角向或矢状根尖移动之间无相关性,但DeShields43发现这些参数之间存在显著关联。Mirabella和Årtun44提出,前牙和后牙的牙根移动均可导致ERR。Parker和Harris20一致观察到ERR与垂直根尖移动、垂直切端移动和切牙唇倾度变化之间的预测关系。Currell等45报告了压入移动与ERR之间的关联,而Tieu等46证明了与前后向根尖位移和回缩力学之间的关联。此外,ERR的分布和定位似乎因正畸牙齿移动的类型和方向而异。据报道,暴露于压缩力的区域比承受拉伸力的区域更容易发生牙根吸收38。例如,颊向倾斜移动主要与颊颈部和舌根尖部牙根表面的吸收相关47,48。颊侧牙根转矩主要导致颊根尖区和腭颈部区域的吸收23。围绕牙齿长轴的旋转移动显示出在颊面与远中面之间以及舌面与近中面之间界面处吸收的倾向性24。远中牙根倾斜显示根尖和中三分之一远中面的吸收增加,颈部三分之一近中面也有额外吸收25。伸长移动为 表3 前牙牙根长度、投影牙根表面积和外源性牙根吸收的变化。牙根长度 投影牙根表面积 外源性牙根吸收 (平均值 ± 标准差,mm) 治疗前 (平均值 ± 标准差,mm) 治疗后 (平均值 ± 标准差,mm) 差值 (平均值 ± 标准差,mm) 治疗前 (平均值 ± 标准差,mm²) 治疗后 (平均值 ± 标准差,mm²) 差值 (平均值 ± 标准差,mm²) 牙位编号 上颌 11 13.47 ± 2.45 12.88 ± 2.26 −0.59 ± 1.36 61.02 ± 18.20 49.07 ± 14.46 −11.95 ± 9.50 0.92 ± 0.69 12 12.51 ± 2.10 11.51 ± 2.11 −1.00 ± 0.94 42.24 ± 11.37 33.63 ± 8.97 −8.61 ± 6.82 1.10 ± 0.99 13 13.14 ± 2.00 11.80 ± 2.00 −1.34 ± 1.12 44.47 ± 14.22 38.69 ± 13.09 −5.77 ± 4.76 1.24 ± 1.16 21 13.53 ± 2.42 12.54 ± 2.44 −0.99 ± 0.78 57.29 ± 19.88 45.99 ± 11.70 −11.31 ± 12.93 1.05 ± 0.83 22 12.63 ± 1.77 11.55 ± 2.14 −1.08 ± 0.93 39.96 ± 11.95 34.07 ± 10.64 −5.89 ± 6.72 1.10 ± 0.98 23 13.30 ± 1.63 12.11 ± 1.36 −1.19 ± 0.92 44.29 ± 8.83 37.56 ± 7.55 −6.73 ± 4.93 1.27 ± 1.00 下颌 31 11.77 ± 1.82 10.59 ± 1.85 −1.18 ± 1.40 32.30 ± 11.31 25.97 ± 8.82 −6.33 ± 6.60 1.27 ± 1.57 32 12.41 ± 1.96 11.42 ± 2.00 −0.99 ± 1.57 35.12 ± 10.83 29.88 ± 8.48 −5.24 ± 6.26 1.24 ± 1.94 33 13.42 ± 2.18 12.73 ± 1.96 −0.69 ± 0.92 44.09 ± 12.13 39.35 ± 10.61 −4.74 ± 6.56 0.72 ± 0.99 41 11.70 ± 1.78 10.50 ± 1.84 −1.20 ± 1.16 32.32 ± 9.81 25.75 ± 7.24 −6.57 ± 5.80 1.30 ± 1.41 42 12.39 ± 2.52 11.07 ± 1.92 −1.32 ± 1.16 35.38 ± 12.46 28.83 ± 8.83 −6.55 ± 6.71 1.49 ± 1.41 43 13.79 ± 2.09 12.72 ± 2.15 −1.07 ± 0.92 47.71 ± 12.48 40.79 ± 9.62 −6.92 ± 7.39 1.11 ± 0.94 mm,毫米;mm²,平方毫米;SD,标准差。
文章出版中 科学报告 https://doi.org/10.1038/s41598-026-77321-y ∆根长度 ∆投影根表面积 外部根吸收 ∆根长度 ∆投影根表面积 外部根吸收 牙位 牙齿移动 牙位 牙齿移动 上颌中切牙 11 压低 r −0.091 0.051 0.008 21 压低 r 0.025 0.119 −0.050 p 0.537 0.731 0.956 p 0.867 0.419 0.735 伸长 r 0.051 −0.381 −0.326 伸长 r 0.060 −0.084 −0.064 p 0.733 0.008(BH− p = 0.131) 0.024(BH− p = 0.269) p 0.683 0.571 0.665 颊舌向平移 r −0.149 −0.135 0.199 颊舌向平移 r 0.095 0.267 −0.092 p 0.314 0.361 0.176 p 0.521 0.067 0.534 近远中向平移 r −0.261 0.014 0.284 近远中向平移 r −0.030 −0.208 0.059 p 0.073 0.924 0.050 p 0.842 0.156 0.691 近远中向牙冠倾斜 r 0.054 −0.423 0.254 近远中向牙冠倾斜 r −0.224 −0.245 0.209 p 0.716 0.003(BH− p = 0.084) 0.081 p 0.126 0.094 0.153 颊舌向牙冠转矩 r −0.059 −0.351 0.234 颊舌向牙冠转矩 r −0.225 −0.353 0.214 p 0.690 0.014(BH− p = 0.201) 0.109 p 0.124 0.014(BH− p = 0.201) 0.144 旋转 r 0.161 −0.071 −0.059 旋转 r −0.080 −0.273 0.067 p 0.274 0.631 0.692 p 0.589 0.060 0.653 上颌侧切牙 12 压低 r −0.202 0.141 0.136 22 压低 r −0.155 0.277 0.205 p 0.168 0.340 0.355 p 0.293 0.057 0.162 伸长 r 0.123 −0.161 −0.033 伸长 r 0.207 −0.429 −0.292 p 0.404 0.275 0.825 p 0.157 0.002 0.044 颊舌向平移 r −0.115 −0.100 0.090 颊舌向平移 r 0.300 0.190 −0.274 p 0.436 0.501 0.543 p 0.038 0.196 0.060 近远中向平移 r −0.180 0.167 0.180 近远中向平移 r 0.020 0.008 −0.001 p 0.222 0.257 0.221 p 0.894 0.956 0.996 近远中向牙冠倾斜 r 0.088 −0.035 −0.087 近远中向牙冠倾斜 r −0.456 −0.176 0.555 p 0.551 0.811 0.558 p 0.001 0.231 < 0.001(BH− p = 0.006) 颊舌向牙冠转矩 r −0.504 −0.284 0.527 颊舌向牙冠转矩 r −0.036 0.075 0.127 p < 0.001(BH− p = 0.016) 0.051 < 0.001(BH− p = 0.010) p 0.807 0.611 0.391 旋转 r 0.215 −0.277 −0.191 旋转 r 0.040 0.109 0.030 p 0.143 0.057 0.193 p 0.790 0.459 0.842 表4 上颌牙齿吸收参数与牙齿移动量之间的相关性分析。
文章出版中 科学报告 https://doi.org/10.1038/s41598-026-77321-y ∆根长度 ∆投影根表面积 外部根吸收 ∆根长度 ∆投影根表面积 外部根吸收 上颌尖牙 13 压低 r 0.155 0.090 −0.244 23 压低 r −0.029 −0.126 0.019 p 0.294 0.544 0.095 p 0.847 0.395 0.896 伸长 r −0.212 −0.552 0.161 伸长 r −0.134 −0.283 0.128 p 0.148 < 0.001(BH− p = 0.006) 0.273 p 0.362 0.052 0.388 颊舌向平移 r −0.089 −0.379 0.123 颊舌向平移 r −0.015 −0.027 0.008 p 0.548 0.008(BH− p = 0.131) 0.404 p 0.919 0.855 0.957 近远中向平移 r 0.146 0.108 −0.177 近远中向平移 r −0.164 −0.153 0.165 p 0.321 0.464 0.229 p 0.265 0.298 0.261 近远中向牙冠倾斜 r 0.063 −0.094 −0.074 近远中向牙冠倾斜 r −0.049 0.083 0.102 p 0.669 0.526 0.616 p 0.739 0.575 0.489 颊舌向牙冠转矩 r −0.020 −0.195 0.025 颊舌向牙冠转矩 r −0.299 0.154 0.320 p 0.894 0.185 0.867 p 0.039(BH− p = 0.298) 0.296 0.027(BH− p = 0.269) 旋转 r 0.220 0.014 −0.200 旋转 r −0.156 −0.185 0.161 p 0.132 0.923 0.173 p 0.291 0.208 0.275 ∆,差值;BL,颊舌向;MD,近远中向;r,Pearson相关系数;p,未校正p值;BH-p,Benjamini–Hochberg错误发现率校正p值。仅对未校正p < 0.05的相关性报告BH-p值。表4(续)
出版前文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y Δ根长度 Δ投影根表面积 外根吸收 Δ根长度 Δ投影根表面积 外根吸收 牙位 牙齿移动 牙位 牙齿移动 下颌中切牙 31 压低 r − 0.293 0.101 0.301 41 压低 r − 0.138 − 0.055 0.229 p 0.043(BH− p = 0.300) 0.496 0.038(BH− p = 0.298) p 0.350 0.712 0.117 伸长 r − 0.059 − 0.368 0.042 伸长 r − 0.044 0.085 0.016 p 0.690 0.010(BH− p = 0.157) 0.775 p 0.765 0.564 0.915 颊舌向平移 r − 0.401 − 0.334 0.395 颊舌向平移 r − 0.213 − 0.432 0.315 p 0.005(BH− p = 0.112) 0.020(BH− p = 0.254) 0.005(BH− p = 0.113) p 0.146 0.002(BH− p = 0.084) 0.029(BH− p = 0.274) 近远中向平移 r − 0.071 − 0.009 0.002 近远中向平移 r 0.292 0.081 − 0.289 p 0.631 0.954 0.990 p 0.044(BH− p = 0.300) 0.585 0.046(BH− p = 0.308) 近远中向牙尖 r − 0.303 − 0.222 0.230 近远中向牙尖 r 0.124 0.104 − 0.129 p 0.037(BH− p = 0.298) 0.129 0.116 p 0.399 0.484 0.383 颊舌向冠转矩 r − 0.004 − 0.039 0.002 颊舌向冠转矩 r 0.103 − 0.146 − 0.016 p 0.979 0.790 0.987 p 0.488 0.324 0.911 旋转 r − 0.143 − 0.209 0.185 旋转 r 0.165 − 0.262 − 0.149 p 0.332 0.153 0.208 p 0.263 0.072 0.312 下颌侧切牙 32 压低 r − 0.272 − 0.118 0.280 42 压低 r 0.192 0.139 − 0.154 p 0.062 0.423 0.054 p 0.192 0.347 0.296 伸长 r 0.175 0.248 − 0.143 伸长 r − 0.108 − 0.145 0.103 p 0.235 0.090 0.332 p 0.465 0.324 0.487 颊舌向平移 r 0.073 − 0.214 0.062 颊舌向平移 r − 0.090 − 0.037 0.031 p 0.621 0.145 0.675 p 0.542 0.800 0.836 近远中向平移 r 0.160 0.090 − 0.102 近远中向平移 r − 0.175 − 0.326 0.072 p 0.276 0.544 0.490 p 0.234 0.024(BH− p = 0.269) 0.625 近远中向牙尖 r − 0.091 − 0.125 0.212 近远中向牙尖 r − 0.054 − 0.386 0.046 p 0.537 0.396 0.148 p 0.716 0.007(BH− p = 0.130) 0.759 颊舌向冠转矩 r 0.265 0.160 − 0.140 颊舌向冠转矩 r − 0.152 − 0.156 0.072 p 0.069 0.278 0.343 p 0.303 0.290 0.626 旋转 r 0.080 0.071 − 0.046 旋转 r 0.139 0.017 − 0.098 p 0.590 0.632 0.758 p 0.348 0.907 0.506 表5 下颌牙齿吸收参数与牙齿移动量之间的相关性分析。
出版前文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y Δ根长度 Δ投影根表面积 外根吸收 Δ根长度 Δ投影根表面积 外根吸收 牙位 牙齿移动 牙位 牙齿移动 下颌尖牙 33 压低 r 0.339 0.139 − 0.321 43 压低 r 0.038 0.161 − 0.032 p 0.019(BH− p = 0.247) 0.347 0.026(BH− p = 0.269) p 0.798 0.275 0.828 伸长 r − 0.176 0.064 0.158 伸长 r − 0.002 0.078 − 0.045 p 0.232 0.668 0.282 p 0.990 0.597 0.761 颊舌向平移 r 0.304 0.302 − 0.296 颊舌向平移 r 0.062 − 0.144 − 0.046 p 0.036(BH− p = 0.298) 0.037(BH− p = 0.298) 0.041(BH− p = 0.300) p 0.675 0.329 0.757 近远中向平移 r 0.323 − 0.011 − 0.318 近远中向平移 r − 0.053 − 0.063 0.066 p 0.025(BH− p = 0.269) 0.942 0.028(BH− p = 0.269) p 0.718 0.668 0.657 近远中向牙尖 r 0.098 0.144 − 0.113 近远中向牙尖 r 0.199 − 0.061 − 0.249 p 0.509 0.330 0.445 p 0.175 0.678 0.088 颊舌向冠转矩 r 0.004 0.094 − 0.012 颊舌向冠转矩 r − 0.037 − 0.400 0.042 p 0.980 0.524 0.937 p 0.805 0.005(BH− p = 0.112) 0.778 旋转 r 0.016 − 0.102 − 0.013 旋转 r 0.089 − 0.420 − 0.098 p 0.916 0.491 0.930 p 0.548 0.003(BH− p = 0.084) 0.509 Δ,差值;BL,颊舌向;MD,近远中向;r,Pearson相关系数;p,未校正p值;BH-p,Benjamini–Hochberg错误发现率校正p值。仅对未校正p < 0.05的相关性报告BH-p值。表5(续)
预出版文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y 主要与远端根面吸收增加相关49。关于压低移动,研究结果并不完全一致。一项研究报告称,近中和远中根面吸收增加50,尤其是在根尖三分之一区域。相反,另一项研究观察到颊侧颈部根面和根尖三分之一区域的吸收更为明显51。在本研究中,移动特异性的根吸收参数关联主要在上颌前牙中观察到。尽管有这些发现,所观察到的吸收性变化的幅度总体上是有限的。因此,观察到的关联应在这些变化幅度相对较小的背景下进行解读,不应必然被视为具有临床意义的牙根损伤的指征。这可能部分反映了有意纳入中度拥挤的非拔牙病例,以维持相对标准化的治疗条件。虽然这些标准减少了临床异质性,但也可能限制研究结果对需要更广泛或更复杂正畸移动患者的普适性。在上颌侧切牙中,牙冠转矩成为最一致的移动特异性因素,显示与根长变化和ERR均有关联。在多变量分析中,牙冠转矩每增加1度,与根长减少约0.07 mm和ERR增加约0.09 mm相关。牙冠倾斜也与ERR呈正相关;然而,这种关联在多变量模型中减弱,提示其与ERR的关系可能部分与表6 正畸牙齿移动与根吸收参数之间关联的多变量线性回归分析。牙组/移动B标准误β p t95% CI 下限上限 上颌侧切牙—Δ根长(mm) 压低− 0.1150.315− 0.0530.717− 0.364− 0.7470.518 伸长0.2450.1690.1330.1541.448− 0.0950.585 颊舌向平移0.2970.1550.1820.0621.915− 0.0150.609 近远中向平移0.0500.1710.0270.7730.290− 0.2950.394 牙冠近远中向倾斜− 0.0490.029− 0.1850.094− 1.711− 0.1080.009 牙冠颊舌向转矩− 0.0700.033− 0.2490.041*− 2.105− 0.136− 0.003 旋转0.0320.0130.3130.018*2.4570.0060.058 模型摘要R 2 = 0.171校正R 2 = 0.095F(7,47) = 5.625 p < 0.001 上颌侧切牙—外部根吸收(mm) 压低− 0.0430.287− 0.0190.881− 0.150− 0.6200.534 伸长− 0.2760.173− 0.1430.117− 1.596− 0.6240.072 颊舌向平移− 0.3110.155− 0.1810.051− 2.003− 0.6230.001 近远中向平移− 0.0110.175− 0.0060.950− 0.063− 0.3640.342 牙冠近远中向倾斜0.0600.0310.2140.0571.951− 0.0020.122 牙冠颊舌向转矩0.0930.0340.3150.008*2.7610.0250.160 旋转− 0.0300.015− 0.2800.044*− 2.072− 0.060− 0.001 模型摘要R 2 = 0.190校正R 2 = 0.115F(7,47) = 8.105 p < 0.001 上颌尖牙—Δ投影根表面积(mm 2 ) 压低− 2.0201.602− 0.1510.214− 1.260− 5.2431.204 伸长− 6.2921.538− 0.415< 0.001*− 4.091− 9.385− 3.198 颊舌向平移− 1.0521.058− 0.1350.325− 0.994− 3.1811.077 近远中向平移0.7020.8520.0990.4140.824− 1.0122.417 牙冠近远中向倾斜− 0.0650.104− 0.0360.535− 0.625− 0.2740.144 牙冠颊舌向转矩0.0460.1300.0320.7240.356− 0.2150.308 旋转− 0.0200.082− 0.0310.804− 0.250− 0.1850.144 模型摘要R 2 = 0.198校正R 2 = 0.124F(7,47) = 3.096 p = 0.009 B,非标准化回归系数;SE,患者聚类稳健标准误;β,标准化回归系数;BL,颊舌向;MD,近远中向;R 2 ,决定系数;t,t统计量。每个回归模型纳入来自48名患者的96颗牙齿。所有模型均对侧别进行了校正。使用患者聚类稳健标准误以考虑同一患者双侧牙齿之间的依赖性。F统计量代表校正侧别后七个牙齿移动预测变量的联合检验。*p < 0.05。
预出版文章《Scientific Reports》https://doi.org/10.1038/s41598-026-77321-y 伴随的牙齿移动。相比之下,在多变量分析中,对其他牙齿移动变量和侧别进行调整后,旋转仍与根长变化和ERR均相关;然而,这些发现应谨慎解读。在中度拥挤的病例中,上颌侧切牙可能因空间不足而位于腭侧,其在OT期间向颊侧重新定位可能需要大幅度的根冠转矩移动52。转矩移动由于其力-力矩特性,与根面应力集中增加相关,这可能有助于解释所观察到的与ERR的关联。然而,转矩与根长变化之间的关联也应结合全景放射摄影的局限性加以解读。根朝向的颊舌向变化可能改变根的二维投影,从而促成所测根长的表观变化。因此,转矩相关发现可能既反映对正畸加载的生物学反应,也反映二维成像固有的投影相关效应。在多变量回归分析中,较大的压出与上颌尖牙PRSA的减少相关,压出每增加1 mm,估计减少约6.29 mm²。ERR受受力表面积和所施加力的大小等因素影响6,53。尽管压出通常比压入产生更小的牙周压缩环境,但本研究中观察到的关联提示,在特定生物力学条件下,压出移动仍可能与根吸收相关。这一发现也可能反映上颌尖牙相对更大的压出需求,这可能与其频繁发生的上颌超出咬合有关。这通常由乳尖牙早期脱落引起,导致间隙丧失,继而发生恒尖牙阻生或上颌超萌54。综上所述,这些发现提示,所观察到的关联可能不仅反映压出移动本身的特征,还反映上颌尖牙特定的生物力学需求和初始位置特征。正畸诱导的ERR没有确定性治疗方法;管理主要涉及早期放射学检测,并在必要时通过调整力值或缩短治疗时长来尽量减少进一步进展55。在本研究中,所有放射学评估均在OT完成后立即进行。因此,观察到的ERR反映的是早期阶段状态。正如Asscherickx等人56所证明的,正畸微种植体植入引起的根缺损通过牙骨质沉积表现出大量修复。这提示某些吸收性变化可能具有生物学修复潜力,取决于患者个体情况。因此,治疗后吸收性发现应谨慎解读。牙髓干预应以全面的临床评估为指导,考虑患者症状、功能问题和系列放射学评估,而非仅依赖治疗后图像。在疑似晚期吸收的病例中,正畸医生与牙髓医生之间的协作可有助于支持适当的治疗决策并改善长期结局。从临床角度看,在需要大量牙齿移动的病例中——尤其是被认为更易发生ERR的牙齿——临床医生可考虑增加矫治器数量以减少每个矫治器的移动幅度。这种方法可能通过减少过度或突然的力来帮助最小化ERR风险。然而,增加矫治器数量也可能延长总体治疗时长,而治疗时长本身与根吸收增加相关。因此,在为ERR较高风险患者规划矫治器方案时,临床医生应仔细平衡力值大小和治疗时长。在此类易感牙齿中,在整个OT期间定期进行放射学监测可能有助于检测ERR的早期征象,并允许及时修改力的施加。除这些临床意义外,本研究发现还有助于理解对正畸治疗的生物学反应可能因牙齿类型和计划牙齿移动的特定生物力学而异。在解读本研究发现时,应考虑若干局限性。一个重要局限性是,在常规临床OT期间多种移动同时发生,因此难以分离特定牙齿移动对ERR的个体贡献。尽管如此,本研究在密切反映常规正畸实践的条件下提供了临床相关见解。评估特定牙齿移动对吸收的孤立效应通常需要有限元模型或in
本文为《科学报告》在线预发表文章 https://doi.org/10.1038/s41598-026-77321-y 体外模拟研究。既往动物研究和有限元分析已探讨了不同牙齿移动的吸收效应,而本研究通过在CA治疗(一种广泛应用的的正畸操作)的临床背景下检查移动特异性关联,对这一证据进行了补充。此外,多变量分析显示上颌侧切牙和尖牙的根吸收参数与移动特异性存在关联,为今后研究CA治疗期间ERR的生物力学决定因素提供了依据。另外,由于下颌骨在非CBCT数据集中无法提供经过验证的稳定解剖参考用于重叠配准,下颌模型相对于重叠配准的上颌骨进行定向,这与先前发表的数字模型分析一致。尽管严格的纳入标准将明显的咬合或骨骼变化降至最低,但不能排除咬合沉降对下颌测量的轻微影响。此外,本研究所有患者均使用单一矫治器材料的ClearCorrect系统进行治疗。由于材料特性可能影响力传导,不同系统的结果可能有所差异。另一局限性是使用了全景放射片,尽管ERR是一个三维过程。锥形束计算机断层扫描(CBCT)可对这些变化提供更准确的体积测量7,但其较高的辐射剂量限制了常规使用57。相反,二维成像方法如全景放射片可能高估ERR,如Gandhi等58所报道,二维技术测得的平均吸收值比CBCT高0.2 mm。尽管根尖片和CBCT更为敏感,但它们在OT期间并非常规使用59。本研究中,全景放射片用于OT前后的评估,反映了标准正畸实践并支持回顾性分析。既往ERR研究同样使用了全景放射片11,29,31,60。Stramotas等61报道,当连续全景放射片之间咬合面角度保持相似且咬合面倾斜不超过约10°时,连续全景放射片上的线性测量可能足够准确。为尽量减少测量误差,所有放射片均由经过培训的助理使用同一全景设备以标准化体位拍摄。然而,连续检查之间头部和咬合面体位的微小差异不能完全排除,可能导致了线性测量的变异。尽管全景成像易产生尺寸失真62,但研究全程一致的成像方案和设备使用旨在尽量减少测量变异性并提高OT前后测量的可比性。PRSA相对于根长度的比例性更大降低,应结合这些测量不同的维度和解剖特征来解读。根长度是一维测量,主要反映根尖-冠向缩短,而PRSA代表整个根轮廓所包围的二维投影面积;因此,不应期望出现等效的比例变化。此外,主要影响根尖区域的变化可能在根长度上更为明显,而涉及近中和远中根轮廓的改变可能在PRSA上产生更大变化,而根尖-冠向缩短相对有限。因此,PRSA更大的比例性降低不一定表明与根长度的相应变化存在差异,而是反映了这两种测量所捕获的根形态的互补方面。然而,必须承认二维成像的固有局限性。特别是,颊舌向转矩移动可能影响根的放射投影,可能导致测得根长度的表观减少。因此,在推断真实形态学根吸收时,涉及线性根测量的转矩相关关联应谨慎解读。同样,PRSA测量也可能受到牙齿方向变化(包括转矩和轴向旋转)的影响,因为三维结构在二维图像上的投影面积可能随其方向而变化。因此,PRSA不应被解读为形态学根丢失的直接三维测量。与既往研究不同,本研究不仅量化了根尖处的ERR,还量化了近中和远中表面的ERR,使用数字面积计算程序评估投影根表面积的变化11,31。此外,分析仅限于切牙和尖牙——最易发生ERR的牙齿——尽管未来需要涉及其他牙齿类型的研究来推广这些发现。
预出版文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 本研究具有若干关键优势。第一,研究采用严格的纳入和排除标准,以提高研究人群的一致性和可比性。第二,分析同时纳入线性与投影牙根表面测量,使牙根变化能够通过互补的二维指标进行评估。第三,本研究评估了不同前牙组中广泛的的正畸牙齿移动,并在CA治疗的临床背景下考察其与多个牙根吸收参数的关联。这种综合方法为理解CA治疗期间移动特异性牙根变化模式提供了具有临床相关性的探索性数据。结论 在本研究的局限性范围内,CA治疗期间特定正畸牙齿移动与ERR相关变化之间的关联因移动类型和牙齿组而异。上颌侧切牙的牙冠转矩和上颌尖牙的伸长显示出与牙根吸收参数最为一致的移动特异性关联。这些发现应谨慎解读,并需在未来使用三维成像和更大队列的前瞻性研究中予以证实。作者贡献 S.K.:概念构思、方法学、监督、数据获取与解读、撰写初稿。M.H.B.:概念构思、方法学、数据获取与解读、分析、对稿件进行批判性修订。所有作者均给予最终批准并同意对工作的各个方面负责。资助 作者未获得本研究、 authorship 和/或本文发表方面的财务支持。数据可用性 本研究的 datasets 可应合理要求从通讯作者处获取。声明 伦理批准与参与同意 本研究按照《赫尔辛基宣言》的原则进行。已获得所有参与者(包括患者和参与本研究的儿童的法定监护人)的知情同意。批准由土耳其屈塔希亚健康科学大学临床研究伦理委员会授予(编号:2024/07-36)。利益冲突 作者声明无利益冲突。开放获取 本文根据 Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License 许可,允许任何非商业性使用、共享、分发和以任何媒介或格式复制,只要您给予原作者和来源适当署名,提供 Creative Commons 许可的链接,并指明您是否修改了许可材料。根据本许可,您无权共享源自本文或其部分的改编材料。本文中的图像或其他第三方材料包含在文章的 Creative Commons 许可中,除非在材料的信用行中另有说明。如果材料未包含在文章的 Creative Commons 许可中,且您的预期用途不被法定法规允许或超出允许使用范围,您需要直接从版权持有人处获得许可。要查看本许可的副本,请访问 http://creativecommons.org/licenses/by-nc-nd/4.0/。参考文献 1. Topkara, A., Karaman, A. I. & Kau, C. H. 正畸力引起的根尖牙根吸收:简要综述与长期观察。Eur. J. Dent. 6, 445–453 (2012)。2. Harris, E. F., Kineret, S. E. & Tolley, E. A. 接受正畸治疗患者外部根尖牙根吸收的遗传成分。Am. J. Orthod. Dentofac. Orthop. 111, 301–309 (1997)。3. Pinheiro, L. H. M. 等。正畸患者外部根尖牙根吸收风险的遗传变异:系统综述。Clin. Oral Investig. 25, 5613–5627 (2021)。
预出版文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 4. Weltman, B., Vig, K. W., Fields, H. W., Shanker, S. & Kaizar, E. E. 正畸牙齿移动相关的牙根吸收:系统综述。Am. J. Orthod. Dentofac. Orthop. 137, 462–476 (2010)。5. Årtun, J., Van’t Hullenaar, R., Doppel, D. & Kuijpers-Jagtman, A. M. 识别存在严重根尖牙根吸收风险的正畸患者。Am. J. Orthod. Dentofac. Orthop. 135, 448–455 (2009)。6. Sameshima, G. T. & Sinclair, P. M. 预测和预防牙根吸收:第一部分。诊断因素。Am. J. Orthod. Dentofac. Orthop. 119, 505–510 (2001)。7. Butsabul, P., Kanpittaya, P. & Nantanee, R. CBCT检测的隐形矫治器治疗中的牙根吸收:系统综述和荟萃分析。Int. Dent. J. 74, 1326–1336 (2024)。8. Fujiyama, K., Honjo, T., Suzuki, M., Matsuoka, S. & Deguchi, T. 接受Invisalign矫治器治疗病例的疼痛程度分析:与固定edgewise矫治器治疗的比较。Prog. Orthod. 15, 1–7 (2014)。9. Buschang, P. H., Shaw, S. G., Ross, M., Crosby, D. & Campbell, P. M. 矫治器治疗与传统edgewise托槽的时间效率比较。Angle Orthod. 84, 391–396 (2014)。10. Li, Y. 等。隐形矫治器与固定矫治器正畸治疗期间根尖牙根吸收的患病率和严重程度:锥形束计算机断层扫描研究。Prog. Orthod. 21, 1–8 (2020)。11. Kurnaz, S. & Buyukcavus, M. H. 正畸治疗期间下颌磨牙外部牙根吸收的全景评估:固定矫治器或隐形矫治器治疗的根管充填牙与活髓牙的比较。BMC Oral Health. 24, 1152 (2024)。12. Singh, S., Jain, R. K. & Balasubramaniam, A. 接受隐形矫治器与固定正畸矫治器治疗受试者之间外部根尖牙根吸收的比较评估:系统综述和荟萃分析。J. Dent. Res. Dent. Clin. Dent. Prospects. 18, 85–94 (2024)。13. Darvizeh, A. 等。隐形矫治器与固定矫治器正畸治疗后的外部根尖牙根吸收:系统综述和荟萃分析。Dent. J. 13, 580 (2025)。14. Selvaraj, M., Nivethitha, B. & Madhan, B. 隐形矫治器与固定矫治器治疗相比正畸诱导的外部根尖牙根吸收:伞形综述。J. Oral Biol. Craniofac. Res. 15, 938–947 (2025)。15. Kumasako-Haga, T., Konoo, T., Yamaguchi, K. & Hayashi, H. 8小时间歇性正畸力对破骨细胞和牙根吸收的影响。Am. J. Orthod. Dentofac. Orthop. 135, 278e1–278e9 (2009)。16. Elhaddaoui, R., Qoraich, H. S., Bahije, L. & Zaoui, F. 正畸矫治器与牙根吸收:系统综述。Int. Orthod. 15, 1–12 (2017)。17. Brezniak, N. & Wasserstein, A. 正畸诱导的炎性牙根吸收。第一部分:基础科学方面。Angle Orthod. 72, 175–179 (2002)。18. Feller, L., Khammissa, R. A., Thomadakis, G., Fourie, J. & Lemmer, J. 正畸牙齿移动中的根尖外部牙根吸收与修复:生物学事件。Biomed. Res. Int. 2016, 4864195 (2016)。19. Segal, G., Schiffman, P. & Tuncay, O. 外部根尖牙根吸收治疗相关因素的荟萃分析。Orthod. Craniofac. Res. 7, 71–78 (2004)。20. Parker, R. J. & Harris, E. F. 与上颌中切牙外部根尖牙根吸收相关的正畸牙齿移动方向。Am. J. Orthod. Dentofac. Orthop. 114, 677–683 (1998)。21. Han, G., Huang, S., Von den Hoff, J. W., Zeng, X. & Kuijpers-Jagtman, A. M. 正畸压低和伸长后的牙根吸收:一项个体内研究。Angle Orthod. 75, 912–918 (2005)。22. Nakano, T. 等。不同类型牙齿移动和力值大小对大鼠牙齿移动量和牙根吸收的影响。Angle Orthod. 84, 1079–1085 (2014)。23. Bartley, N. 等。牙根骨水泥的物理特性:第17部分。施加2.5和15度颊侧牙根转矩4周后的牙根吸收:微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 139, e353–e360 (2011)。24. Wu, A. T. 等。牙根骨水泥的物理特性:第18部分。施加轻力和重力控制旋转正畸力4周后牙根吸收的程度:微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 139, e495–e503 (2011)。25. King, A. D. 等。牙根骨水泥的物理特性:第21部分。施加2.5和15度根尖转矩4周后牙根吸收的程度:微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 140, e299–e305 (2011)。26. Hamilton, R. & Gutmann, J. 牙髓-正畸关系:综合治疗计划挑战的综述。Int. Endod. J. 32, 343–360 (1999)。27. Flatten, J., Gedrange, T., Bourauel, C., Keilig, L. & Konermann, A. 骨和牙根吸收在承受正畸力的下颌前牙生物力学行为中的作用:有限元方法。Biomedicines 12, 1959 (2024)。28. Zhang, M., Zhang, P., Koh, J.-T., Oh, M.-H. & Cho, J.-H. 矫治器与牙根吸收的评估:系统综述的概述。J. Clin. Med. 13, 1950 (2024)。29. Lee, Y. J. & Lee, T. Y. 根管充填牙与对侧活髓牙正畸治疗期间的外部牙根吸收: contributing factors的临床研究。Am. J. Orthod. Dentofac. Orthop. 149, 84–91 (2016)。30. Almagrami, I. M. 等。隐形矫治器与固定正畸矫治器之间牙槽骨变化和根尖牙根吸收的定量三维比较研究。Prog. Orthod. 24, 6 (2023)。
出版前文章《科学报告》https://doi.org/10.1038/s41598-026-77321-y 31. Kurnaz, S. & Buyukcavus, M. H. 拔牙与非拔牙正畸治疗后根管充填牙与活髓牙的外根吸收:一项分口回顾性研究。Acta Odontol. Scand. 79, 282–288 (2021)。32. Cangialosi, T. J. 等。ABO差异指数:病例复杂度的衡量指标。Am. J. Orthod. Dentofac. Orthop. 125, 270–278 (2004)。33. Garib, D. 等。利用腭皱襞重叠上颌数字化模型:老龄化是否影响可靠性?Orthod. Craniofac. Res. 22, 183–193 (2019)。34. Linge, L. & Linge, B. O. 正畸治疗期间与根尖根吸收相关的患者特征和治疗变量。Am. J. Orthod. Dentofac. Orthop. 99, 35–43 (1991)。35. Papageorgiou, S. N., Koletsi, D., Iliadi, A., Peltomaki, T. & Eliades, T. 正畸隐形矫治器与固定矫治器的治疗结果:一项系统评价与荟萃分析。Eur. J. Orthod. 42, 331–343 (2020)。36. Liou, E. J. & Chang, P. M. 使用微螺钉整体内收和压低上颌前牙的正畸患者根尖根吸收。Am. J. Orthod. Dentofac. Orthop. 137, 207–212 (2010)。37. Acar, A., Canyürek, Ü., Kocaaga, M. & Erverdi, N. 持续力与间断力应用及根吸收。Angle Orthod. 69, 159–163 (1999)。38. Roscoe, M. G., Meira, J. B. & Cattaneo, P. M. 正畸力系统与根吸收的关联:一项系统评价。Am. J. Orthod. Dentofac. Orthop. 147, 610–626 (2015)。39. Ballard, D. J., Jones, A. S., Petocz, P. & Darendeliler, M. A. 根骨水泥的物理特性:第11部分。持续与间歇控制正畸力对根吸收的影响。一项微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 136, 8e1–8e8 (2009)。40. Yi, J., Xiao, J., Li, Y., Li, X. & Zhao, Z. 无托槽隐形矫治器治疗或固定正畸治疗后非拔牙病例的外根尖吸收。J. Dent. Sci. 13, 48–53 (2018)。41. Guan, W., Lu, Y., Yu-feng, Z., San-lian, L. & Ji-wei, Z. 无托槽隐形矫治器与直丝弓矫治器治疗患者切牙根吸收的回顾性研究。Shanghai Kou Qiang Yi Xue. 26, 121–124 (2017)。42. Phillips, J. R. 正畸治疗下的根尖根吸收。Angle Orthod. 25, 1–22 (1955)。43. DeShields, R. W. 已治疗II类1分类错颌畸形中根吸收的研究。Angle Orthod. 39, 231–245 (1969)。44. Mirabella, A. D. & Årtun, J. 成年正畸患者上颌前牙根尖根吸收的危险因素。Am. J. Orthod. Dentofac. Orthop. 108, 48–55 (1995)。45. Currell, S. D., Liaw, A., Grant, P. D. B., Esterman, A. & Nimmo, A. 正畸力学疗法及其对外根吸收的影响:一项系统评价。Am. J. Orthod. Dentofac. Orthop. 155, 313–329 (2019)。46. Tieu, L. D., Saltaji, H., Normando, D. & Flores-Mir, C. II类1分类错颌畸形非手术正畸治疗后切牙影像学确定的正畸诱导外根尖吸收:一项系统评价。Prog. Orthod. 15, 1–13 (2014)。47. Chan, E. & Darendeliler, M. A. 根骨水泥的物理特性:第5部分。轻力和重力正畸力应用后根吸收陷窝的体积分析。Am. J. Orthod. Dentofac. Orthop. 127, 186–195 (2005)。48. Barbagallo, L. J., Jones, A. S., Petocz, P. & Darendeliler, M. A. 根骨水泥的物理特性:第10部分。隐形可摘热塑性矫治器与轻力和重力正畸力对前磨牙骨水泥影响的比较。一项微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 133, 218–227 (2008)。49. Montenegro, V. C. J., Jones, A., Petocz, P., Gonzales, C. & Darendeliler, M. A. 根骨水泥的物理特性:第22部分。轻力和重力伸长正畸力应用后的根吸收:一项微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 141, e1–e9 (2012)。50. Harris, D. A., Jones, A. S. & Darendeliler, M. A. 根骨水泥的物理特性:第8部分。控制性压入轻力和重力正畸力应用后根吸收陷窝的体积分析:一项微计算机断层扫描研究。Am. J. Orthod. Dentofac. Orthop. 130, 639–647 (2006)。51. Harry, M. & Sims, M. 双尖牙压入中的根吸收:一项扫描电子显微镜研究。Angle Orthod. 52, 235–258 (1982)。52. Nanda, R. S. & Tosun, Y. S. 《正畸生物力学:原理与实践》(Quintessence Publishing Co., 2010)。53. Sameshima, G. T. & Sinclair, P. M. 预测和预防根吸收:第II部分。治疗因素。Am. J. Orthod. Dentofac. Orthop. 119, 511–515 (2001)。54. Proffit, W. R., Fields, H. W. & Sarver, D. M. 《当代正畸学》第4版(Mosby Elsevier, 2007)。55. Kim, S., Hwang, S. J. & Seo, M. S. 正畸治疗对根管治疗牙根尖周透射性的影响:一项CBCT分析。BMC Oral Health. 23, 289 (2023)。56. Asscherickx, K., Vannet, B. V., Wehrbein, H. & Sabzevar, M. M. 微螺钉损伤后的根修复。Clin. Oral Implants Res. 16, 575–578 (2005)。57. Patel, S., Dawood, A., Wilson, R., Horner, K. & Mannocci, F. 使用口内X线摄影和锥形束计算机断层扫描检测和管理根吸收病变——一项体内研究。Int. Endod. J. 42, 831–838 (2009)。
待刊文章《科学报告》 https://doi.org/10.1038/s41598-026-77321-y 58. Gandhi, V. 等。无拔牙病例中隐形矫治器与预调直丝弓矫治器对外部根尖吸收的比较:系统评价与荟萃分析。《欧洲正畸学杂志》43, 15–24 (2021)。59. Yi, J. 等。锥形束计算机断层扫描与根尖片诊断外部牙根吸收的比较:系统评价与荟萃分析。《Angle正畸学杂志》87, 328–337 (2017)。60. Llamas-Carreras, J. 等。根管充填牙与对侧活髓牙正畸治疗期间外部牙根吸收的研究。《国际牙髓病学杂志》43, 654–662 (2010)。61. Stramotas, S., Geenty, J. P., Petocz, P. & Darendeliler, M. A. 体外不同位置拍摄全景放射片上线性和角度测量的准确性。《欧洲正畸学杂志》24, 43–52 (2002)。62. Lien, L. C. & Soh, G. 口腔全景片评估正畸患者牙长度的准确性。《新加坡牙科杂志》23, 68–71 (2000)。出版商注 施普林格·自然对已出版地图和机构隶属关系中的管辖权主张保持中立。作者与隶属关系 Safa Kurnaz 1 · Muhammed Hilmi Buyukcavus 2 Safa Kurnaz safa.kurnaz@ksbu.edu.tr 1 土耳其屈塔希亚健康科学大学牙学院牙髓病学系,43270 屈塔希亚,土耳其 2 土耳其安塔利亚比林大学牙学院正畸学系,安塔利亚,土耳其
查看英文原文▾
ARTICLE | Article in Press Scientific Reportshttps://doi.org/10.1038/s41598-026-77321-y Scientific Reports This Article in Press is shared early to give you faster access to new research. It is citable and carries a permanent DOI. The final edited version will replace it automatically. Quantifying external root resorption associated with various orthodontic tooth movements in clear aligner treatment: a retrospective clinical study Safa Kurnaz · Muhammed Hilmi Buyukcavus Received: 20 May 2026 / Accepted: 19 October 2026 © The Author(s) 2026 Abstract This study investigated associations between tooth movements and external root resorption (ERR) during clear aligner (CA) orthodontic treatment (OT). Forty-eight patients who underwent non-extraction OT using ClearCorrect aligners were included. Maxillary and mandibular incisors and canines were assessed, comprising 576 teeth. Pre- and post-OT digital models were superimposed to quantify intrusion, extrusion, crown tipping, crown torque, rotation, buccolingual translation, and mesiodistal translation. Root length, projected root surface area (PRSA), and ERR were measured on pre- and post-OT panoramic radiographs. Associations between tooth movements and root length, PRSA, and ERR were evaluated using correlation and multivariable linear regression analyses, with correction for multiple comparisons where appropriate. Significant movement-specific associations were observed in the maxillary lateral incisors and canines, with crown torque associated with root length and ERR changes and crown tipping associated with ERR in the lateral incisors, while extrusion was associated with PRSA changes in the canines. No significant movement–resorption associations were observed in the mandibular teeth. The magnitude of the observed resorptive changes was generally limited. Overall, the findings suggest that root responses during CA treatment may vary according to both the type of tooth movement and tooth group; however, the observed movement-specific associations should be interpreted cautiously. Keywords Clear aligner · Endodontics · External root resorption · Orthodontics · Tooth movement Abbreviations ERR External root resorption OT Orthodontic treatment CA Clear aligner ABO American Board of Orthodontics SD Standard deviation CEJ Cementoenamel junction CBCT Cone-beam computed tomography PRSA Projected root surface area FDR False discovery rate
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y Introduction Orthodontic treatment (OT) can lead to a common side effect known as external root resorption (ERR), resulting in irreversible loss of apical root tissue 1 . While OT often results in improved aesthetics and function, with clinically insignificant ERR in most cases, severe cases can compromise tooth structure and potentially lead to tooth loss 2,3 . Several factors can influence the risk of ERR during OT. These include genetics, age, sex, and systemic health conditions, as well as treatment-related factors, such as the severity of malocclusion, root morphology, and the duration and intensity of orthodontic forces. Treatment methods, including tooth extraction, appliance type, and movement magnitude, can also influence ERR 4,5 . Although ERR can occur during the OT of any tooth, the risk is generally higher for maxillary and mandibular incisors 6 . While various factors have been implicated in ERR, the precise mechanisms by which they contribute to this process are not fully understood 7 . Fixed appliance treatments have long been the standard for OT. However, clear aligners (CAs) have gained popularity due to improved comfort, aesthetics, and reduced treatment time 8,9 . Although evidence remains heterogeneous, recent systematic reviews and meta-analyses suggest that the incidence and severity of ERR are generally lower with CAs than with fixed appliances 10–14 . Factors such as intermittent force application, lighter forces, and shorter treatment times may contribute to this reduced risk 15,16 . OT applies multidirectional forces to achieve tooth movement, potentially disrupting the balance between bone resorption and formation at the root apex. Excessive or prolonged forces can exceed the adaptive capacity of periodontal tissues, leading to the progressive loss of root structure 17,18 . Both the magnitude and direction of orthodontic forces can influence the risk of ERR 19,20 . Intrusive forces are more likely to induce ERR than extrusive forces 21 . In addition, an animal study reported that tipping tooth movements were associated with greater ERR compared to bodily tooth movements 22 . Furthermore, human studies also show that increased root torque, tipping, and rotation are associated with higher ERR 23–25 . However, it is difficult to isolate the specific tooth movements that contribute to ERR. This is because orthodontic appliances typically produce a combination of complex movements, including extrusion, intrusion, translation, tipping, torque, and rotation 20 . Many clinicians perceive ERR as a potentially preventable complication during OT and often hold orthodontists responsible for its occurrence. However, detailed information regarding the overall relationships between endodontics and orthodontics in treatment decisions remains limited 26 . Therefore, further research on factors contributing to ERR is essential to improve treatment outcomes and minimize adverse effects. A thorough investigation could determine how the amount of ERR is influenced by the different types and magnitudes of orthodontic tooth movements and displacements 27 . Recent systematic reviews and meta-analyses have primarily focused on comparing the incidence and severity of ERR between CAs and fixed appliances, generally suggesting a lower risk of ERR with CA treatment 12–14 . However, these reviews have also emphasized the heterogeneity of the available evidence, methodological limitations, and the need for further high-quality clinical studies 28 . Importantly, previous studies have largely evaluated ERR as an overall treatment outcome rather than investigating the specific contribution of individual orthodontic tooth movements to ERR during CA treatment. To the best of our knowledge, no previous clinical study has comprehensively evaluated the association between different orthodontic tooth movements, including intrusion, extrusion, translation, crown tipping, crown torque, and rotation, and the magnitude of ERR across different tooth groups in patients treated exclusively with CAs. Therefore, the present study aimed to investigate the association between different orthodontic tooth movements and ERR and to quantify the magnitude of ERR associated with each movement type across different tooth groups following CA treatment, thereby providing clinically relevant evidence on movement-specific resorption patterns. The null hypothesis (H0) of this study posits the following: There is no difference in the amount of ERR in different tooth groups resulting from different orthodontic tooth movements.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y Methods This retrospective clinical study was approved by the local university ethics committee (reference number: 2024/07-36). All participants provided informed consent. For participants under the age of 18, consent was obtained from their legal guardians. A sensitivity power analysis was performed for the multivariable linear regression analyses with seven predictors. With a sample size of 48, an α level of 0.05, and 80% power, the analysis indicated a minimum detectable overall effect size of f 2 = 0.35. The radiographs of patients treated at our clinic between July 2022 and December 2024 were retrospectively analyzed. Patients were selected based on predefined inclusion and exclusion criteria. The inclusion criteria were established to obtain a clinically homogeneous study population and to minimize potential confounding factors that could influence ERR. Accordingly, only patients with high-quality standardized panoramic radiographs suitable for accurate landmark identification, complete permanent dentition without missing teeth (excluding third molars), and who underwent OT with CAs were included. Patients were also required to have skeletal Class I malocclusion, moderate crowding, and treatment plans that did not involve tooth extraction. In addition, only patients without parafunctional habits (e.g., bruxism or clenching), treated with the same CA system, and with a similar planned degree of tooth movement were included 10,11,29–31 . The exclusion criteria were applied to minimize biological and treatment-related factors that could influence ERR. Patients were excluded if they had tooth or root anomalies, including unerupted teeth, regardless of the affected tooth. Patients with a history of dental trauma, previous OT, missing records, or tooth extractions were also excluded. Additional exclusion criteria included systemic diseases, craniofacial anomalies, temporomandibular disorders, poor oral health (caries or periodontal disease), root-filled teeth, open apices, and pre-existing root resorption 10,11,29–31 . To further minimize variability related to case complexity, the American Board of Orthodontics (ABO) discrepancy index was used to assess the difficulty of each case. Patients with high complexity scores were excluded to maintain a more homogeneous study population 30,32 . Treatment protocol After applying inclusion and exclusion criteria to 235 patients, 48 patients (26 females, 22 males) treated with CAs were included. The mean age of the participants was 19.22 ± 2.08 years. Each patient’s maxillary and mandibular central, lateral, and canine teeth were evaluated individually. This resulted in an analysis of 96 teeth per tooth group, totaling 576 teeth. The treatment group consisted of patients who had undergone OT using ClearCorrect CAs (ClearCorrect ® , ClearCorrect LLC, Round Rock, TX, USA). A virtual treatment plan was created using ClearCorrect’s Clear- Pilot™ software (ClearCorrect’s Treatment Planning Tool, version 5) to guide the treatment process. The plan outlined a sequence of treatment steps, including interproximal reduction, attachment placement, and the use of intermaxillary elastics. Patients were instructed to change their aligners every 15 days and wear them for approximately 22 h daily. Measurement of tooth movement All digital model measurements were performed using a standardized workflow to ensure measurement consistency and reproducibility. The same reference landmarks, regions of interest, registration procedure, and measurement tools were used for all pre- and post-treatment models. Digital impressions were captured using a 3Shape TRIOS 3 (3Shape TRIOS A/S, Copenhagen, Denmark) intraoral scanner and uploaded to the ClearPilot platform for treatment design. STL files of the resulting models were saved using 3D Slicer 5.0. Model superimposition was performed according to the palatal rugae-based method described by Garib et al. 33 , using stable palatal reference areas for registration. Landmarks were placed on the incisive papilla’s posterior
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y border, the medial edges of the second and third palatal rugae, and 10 mm distally from the third rugae’s medial edge. Regions of interest were defined around these landmarks. Regions of interest for the incisive papilla and third rugae projections had a 5-vertex radius; second and third rugae regions of interest had a 30 radius, ensuring regions of interest remained distant from the teeth. Using Slicer’s “CMF registration” tool, the post-treatment maxilla was superimposed onto the pre-treatment maxilla via region of interest matching. Mandibular models were then aligned in occlusion with the superimposed maxilla. The maxillary superimposition matrix was applied to the mandibular models, aligning them with the maxilla. Finally, predefined landmarks were placed on all models to quantify the magnitude of orthodontic tooth movement achieved during treatment (Fig. 1). Landmarks were placed on maxillary and mandibular models to measure: (1) horizontal/vertical linear displacements (canine buccal cusp tips and incisor mesial incisal edges), (2) mesiodistal rotations (most mesial/distal points from an occlusal view), and (3) buccolingual tipping (incisor incisal center and deepest cervical point). Using 3D Slicer’s Fiducial Registration Wizard, x, y, and z coordinates for linear displacement landmarks were recorded for pre/post-treatment and initial/final predicted models. Coordinate differences quantified the achieved and predicted changes. Q3DC measured the angles between mesial/distal landmarks (mesiodistal rotations) and incisal/cervical landmarks (buccolingual tipping) for each model pair. For each of the 12 teeth evaluated (maxillary centrals, laterals and canines, and mandibular centrals, laterals and canines), the total intrusion, extrusion, buccolingual translation, mesiodistal translation, crown tip, crown torque, and rotation amounts that occurred in the relevant teeth from the beginning to the end of the OT (including refinement aligners) were taken from the system in degrees and/or millimeters. Measurement methodology for external root resorption To minimize measurement bias, all radiographic measurements were performed using a standardized measurement protocol under identical imaging and measurement conditions. The extent of ERR was assessed for all teeth using digital panoramic radiographs obtained before and immediately after OT. All radiographs were taken in natural head position to ensure consistency. This was achieved by instructing patients to look straight ahead into the unit’s built-in mirror, allowing self-alignment. The Castellini X-ray unit (Castellini X Radius Compact, Imola, Italy) was used to acquire the images and to assist with head stabilization through its integrated positioning system. Crown and root length measurements were performed using the measurement software of the Castellini X-ray unit according to a previously validated methodology described in earlier studies 11,29,31,34 . Measurements were obtained from standardized pre- and post-treatment panoramic radiographs using identical anatomical reference points for all teeth. The reference points and procedures used for linear and PRSA measurements are illustrated in Fig. 2. Fig. 1 Superimposition of pre- and post-treatment models using the palatal rugae method.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y The crown and root length measurements were performed in several steps. First, the cementoenamel junction (CEJ) was identified by drawing a straight line between the mesial and distal CEJ points. The crown length was then measured on both the initial and final radiographs. This measurement was defined as the longest distance from the incisal edge to the CEJ. A similar method was applied for root length measurements in all teeth, where the distance from the CEJ to the root apex was recorded 11,29,31,34 . ERR was calculated as follows: ERR (mm) = (pre-OT root length − post-OT root length) × (pre-OT crown length/post-OT crown length). The crown-length ratio was used as a proportional correction factor to compensate for differences in radiographic magnification between the pre- and post-OT panoramic images 11,29,31,34 . Accordingly, ERR was expressed in millimeters. To complement linear measurements, projected root surface area (PRSA) was analyzed from panoramic radiographs (Fig. 2). Before PRSA measurements, each panoramic image was calibrated in SketchAndCalc by defining a known linear distance on the radiographic image and entering its corresponding length in millimeters, thereby establishing the pixel-to-millimeter scale. Root surface measurements were then performed using SketchAnd- Calc software (Axiom Welldone, https://www.sketchandcalc.com/). The entire root surface was carefully traced, starting at the CEJ and extending to the radiographic apex, following both the mesial and distal outlines of the root. To ensure measurement precision, digital points were manually placed at one-pixel intervals along the traced root boundary. An average of 25 (± 5) points per root was used, depending on the root length and curvature. The software then automatically computed the enclosed area in square millimeters, representing the PRSA 11,31 . The same calibration and tracing procedure was applied to the pre- and post-OT radiographs. Statistical analysis Data normality was assessed using the Shapiro–Wilk test. As the data were normally distributed, parametric statistical analyses were performed. Measurement reliability was assessed before the main analyses. All measurements were performed by a single calibrated investigator. Prior to the main study, reliability and reproducibility were evaluated using panoramic radiographs from 16 randomly selected patients. The same investigator repeated the measurements after a two-week interval, and a second investigator independently repeated the same process. Fig. 2 Reference points and measurement procedures for crown length, root length, and PRSA on panoramic radiographs.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y The intraclass correlation coefficients for intra- and inter-examiner agreement were 0.953 (95% CI 0.925–0.981) and 0.929 (95% CI 0.902–0.956), respectively, indicating excellent measurement reproducibility. Descriptive statistics were expressed as mean ± standard deviation (SD), frequencies (n), and percentages (%). Pearson correlation analyses were performed to evaluate the associations between the seven orthodontic toothmovement variables (intrusion, extrusion, buccolingual translation, mesiodistal translation, crown tip, crown torque, and rotation) and changes in root length, PRSA, and ERR according to tooth position. To account for multiple testing, correlation p values were adjusted using the Benjamini–Hochberg false discovery rate (FDR) procedure, with FDR-adjusted p values < 0.05 considered statistically significant. Follow-up multivariable linear regression analyses were performed for tooth-group/outcome combinations with at least one association remaining significant after FDR correction. All seven tooth-movement variables were entered simultaneously, with adjustment for side and patient-clustered robust standard errors to account for within-patient dependence. Each regression model included 96 teeth from 48 patients, and results are reported as unstandardized coefficients (B), patient-clustered robust standard errors, standardized coefficients (β), 95% confidence intervals, t statistics, p values, R 2 , adjusted R 2 , and the joint F statistic for the seven tooth-movement predictors. These follow-up regression analyses were considered exploratory. All statistical analyses were performed using IBM SPSS Statistics for Windows, version 26 (IBM Corp., Armonk, NY, USA). Results Table 1 presents descriptive statistics of the study participants, including age, gender, number of aligners, treatment duration, and treatment complexity scores based on ABO criteria. The extent of tooth movement in the patients’ anterior teeth resulting from the CA treatment is shown in Table 2. Figure 3 provides a graphical overview of the distribution of root resorption parameters across the different orthodontic tooth movements. Changes in root length, PRSA, and ERR following OT are presented in Table 3. The results of the correlation analyses between root resorption parameters and the extent of tooth movements are presented in Table 4 for maxillary teeth and Table 5 for mandibular teeth. After Benjamini–Hochberg FDR correction, statistically significant correlations remained only in the maxillary teeth. The results of the multivariable linear regression analyses are presented in Table 6. For the maxillary lateral incisors, crown torque was negatively correlated with root length change (r = − 0.504, BH-p = 0.016) and positively correlated with ERR (r = 0.527, BH-p = 0.010), while crown tip was positively correlated with ERR (r = 0.555, BH-p = 0.006) after FDR correction. In the multivariable regression analyses, each 1-degree increase in crown torque was associated with an approximately 0.07-mm greater reduction in root length (B = − 0.070, 95% CI −0.136 to − 0.003, p = 0.041) and approximately 0.09-mm greater ERR (B = 0.093, 95% CI 0.025 to 0.160, p = 0.008). Rotation was also associated with root length change (B = 0.032, 95% CI 0.006 to 0.058, p = 0.018) and ERR (B = − 0.030, 95% CI −0.060 to − 0.001, p = 0.044) after adjustment for the other tooth-movement variables and side. Characteristics Number (n) GenderFemale26 Male22 CharacteristicsMean ± SD Chronological age (years)19.22 ± 2.08 ABO discrepancy index12.06 ± 1.91 Number of total aligners32 ± 12 Number of refinements1.81 ± 0.75 Treatment duration (years)1.31 ± 0.47 Table 1 Descriptive statistics of the study participants. ABO, American board of orthodontics; n, Number; SD, Standard deviation.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y For the maxillary canines, extrusion was negatively correlated with PRSA change (r = − 0.552, BH-p = 0.006) after FDR correction. In the multivariable regression analysis, each 1-mm increase in extrusion was associated with an approximately 6.29-mm 2 greater reduction in PRSA (B = − 6.292, 95% CI −9.385 to − 3.198, p < 0.001). Discussion Despite the increasing use of CAs, their impact on ERR remains understudied 7 . Previous studies have consistently shown that ERR is less prevalent in patients treated with CAs than in those treated with fixed appliances 10,11,30,35 . Recent systematic reviews and umbrella reviews have reached similar conclusions while also highlighting the heterogeneity of the available evidence, methodological limitations, and the need for cautious interpretation 13,14 . However, these studies primarily compared the overall incidence or severity of ERR between appliance systems rather than evaluating the influence of specific orthodontic tooth movements on ERR in different tooth groups. Table 2 Descriptive statistics of orthodontic tooth movements in anterior teeth. Tooth movementIntrusion (mean ± SD, Min–Max, mm) Extrusion (mean ± SD, Min–Max, mm) Translation (buccolingual) (mean ± SD, Min– Max, mm) Translation (mesiodistal) (mean ± SD, Min– Max, mm) Crown tip (mesiodistal) (mean ± SD, Min–Max, °) Crown torque (buccolingual) (mean ± SD, Min–Max, °) Rotation (mean ± SD, Min–Max, °) Tooth number Maxilla110.20 ± 0.34 (0.0–1.3) 0.19 ± 0.39 (0.0–1.6) 0.80 ± 0.65 (0.0–2.1) 0.55 ± 0.47 (0.0–1.8) 1.85 ± 1.78 (0.0–6.2) 4.36 ± 3.58 (0.0–13.2) 9.15 ± 7.88 (0.0–30.1) 120.28 ± 0.49 (0.0–1.6) 0.37 ± 0.52 (0.0–1.9) 0.82 ± 0.59 (0.0–2.3) 0.64 ± 0.58 (0.0–2.1) 2.95 ± 4.14 (0.0–10.7) 3.34 ± 3.11 (0.0–10.1) 8.29 ± 9.59 (0.0–34.4) 130.29 ± 0.45 (0.0–1.4) 0.13 ± 0.27 (0.0–1.1) 0.70 ± 0.70 (0.0–2.5) 0.54 ± 0.74 (0.0–2.9) 2.23 ± 3.21 (0.0–11.3) 3.69 ± 3.91 (0.0–17.1) 6.90 ± 8.13 (0.0–30.2) 210.21 ± 0.31 (0.0–0.9) 0.25 ± 0.38 (0.0–1.4) 0.90 ± 0.50 (0.2–2.2) 0.60 ± 0.64 (0.0–2.7) 2.91 ± 3.32 (0.0–11.1) 4.38 ± 5.18 (0.0–19.1) 7.12 ± 4.64 (0.0–16.1) 220.21 ± 0.38 (0.0–1.4) 0.42 ± 0.52 (0.0–1.5) 0.60 ± 0.56 (0.0–2.6) 0.42 ± 0.44 (0.0–1.6) 1.65 ± 2.69 (0.0–12.2) 3.89 ± 3.65 (0.0–10.3) 9.68 ± 8.93 (0.0–34.5) 230.14 ± 0.24 (0.0–0.9) 0.22 ± 0.36 (0.0–1.3) 0.62 ± 0.55 (0.0–1.8) 0.57 ± 0.64 (0.0–3.1) 1.26 ± 2.10 (0.0–8.1) 2.09 ± 2.47 (0.0–8.3) 9.21 ± 6.58 (0.0–25.1) Mandible310.55 ± 0.64 (0.0–1.7) 0.07 ± 0.14 (0.0–0.6) 0.99 ± 0.64 (0.0–2.4) 0.68 ± 0.54 (0.0–1.8) 2.52 ± 2.14 (0.0–8.1) 4.64 ± 5.02 (0.0–18.1) 9.45 ± 7.77 (0.0–28.1) 320.64 ± 0.60 (0.0–1.8) 0.03 ± 0.08 (0.0–0.3) 0.92 ± 0.74 (0.0–2.6) 0.88 ± 0.76 (0.0–2.8) 3.85 ± 3.27 (0.0–9.1) 2.69 ± 3.24 (0.0–12.1) 8.88 ± 8.64 (0.0–34.5) 330.49 ± 0.49 (0.0–2.2) 0.05 ± 0.16 (0.0–0.6) 0.74 ± 0.73 (0.0–2.9) 0.81 ± 0.67 (0.0–2.3) 2.46 ± 3.25 (0.0–10.3) 2.70 ± 2.32 (0.0–7.5) 9.07 ± 6.88 (0.0–26.3) 410.59 ± 0.65 (0.0–2.6) 0.07 ± 0.18 (0.0–0.7) 0.89 ± 0.67 (0.0–2.3) 0.87 ± 0.83 (0.0–3.3) 2.52 ± 2.53 (0.0–8.1) 4.70 ± 5.50 (0.0–20.1) 7.45 ± 6.14 (0.0–27.2) 420.68 ± 0.74 (0.0–3.2) 0.08 ± 0.21 (0.0–0.7) 1.03 ± 0.85 (0.0–2.9) 0.72 ± 0.63 (0.0–2.4) 2.16 ± 2.57 (0.0–9.2) 3.97 ± 4.30 (0.0–13.4) 6.77 ± 5.64 (0.0–17.2) 430.49 ± 0.44 (0.0–1.7) 0.04 ± 0.13 (0.0–0.5) 0.82 ± 0.85 (0.0–2.9) 0.62 ± 0.47 (0.0–1.7) 1.98 ± 2.96 (0.0–8.5) 2.35 ± 2.89 (0.0–9.1) 8.69 ± 5.15 (0.0–16.5) SD, standard deviation; Min–Max, minimum–maximum; mm, millimeter; °, degree. Fig. 3 Graphical representation of point distribution for root resorption parameters across various tooth movements.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y Consequently, this study aimed to investigate the associations between specific orthodontic tooth movements and root resorption parameters across different tooth groups and to quantify the resorptive changes associated with these movements in patients undergoing CA treatment. The present findings demonstrated movement- and toothspecific associations with root resorption parameters, predominantly in the maxillary anterior teeth. Significant associations were observed particularly for crown torque in the maxillary lateral incisors and extrusion in the maxillary canines, with related patterns also observed in the exploratory multivariable analyses. The magnitude, duration, and continuity of orthodontic forces are critical determinants of ERR 36,37 . Heavier and continuous forces correlate with an increased prevalence and severity of ERR 4,38 . Conversely, intermittent forces may mitigate resorption by facilitating cementogenic repair 37,39 . The magnitude and direction of forces exerted by CAs differ from those generated by brackets and archwires. CAs are believed to deliver lighter, intermittent forces, potentially contributing to a reduced incidence of ERR compared to fixed appliances 40,41 . This interpretation is supported by recent reviews suggesting that the intermittent force pattern of CAs may allow periods of periodontal and cemental repair, although the clinical evidence remains heterogeneous 14,28 . Few studies have assessed the extent of ERR associated with orthodontic tooth movement, and substantial methodological disparities limit direct comparisons. This study is the first to explore this phenomenon specifically in the context of CA treatment. While Phillips 42 observed no correlation between ERR and angular or sagittal apical movements of the maxillary central incisor, DeShields 43 identified significant associations between these parameters. Mirabella and Årtun 44 posited that both anterior and posterior root movements can contribute to ERR. Parker and Harris 20 consistently observed predictive relationships between ERR and vertical apical movement, vertical incisal movement, and incisor proclination changes. Currell et al. 45 reported an association between intrusive movements and ERR, while Tieu et al. 46 demonstrated an association with anteroposterior apical root displacement and retraction mechanics. Moreover, the distribution and localization of ERR appear to vary according to the type and direction of orthodontic tooth movement. Areas exposed to compressive forces have been reported to be more susceptible to root resorption than those subjected to tensile forces 38 . For example, buccal tipping movements were primarily associated with resorption on the buccocervical and linguoapical root surfaces 47,48 . Buccal root torque resulted in resorption predominantly on the bucco-apical and palato-cervical regions 23 . Rotational movements around the long axis of the tooth demonstrated a predilection for resorption at the interfaces between the buccal and distal surfaces, as well as between the lingual and mesial surfaces 24 . Distal root tipping showed increased resorption on distal aspects in the apical and middle thirds, with additional resorption on mesial aspects in the cervical third 25 . Extrusion movements were Table 3 Changes in root length, projected root surface area, and external root resorption in anterior teeth. Root lengthProjected root surface areaExternal root resorption (mean ± SD, mm) Pre- (mean ± SD, mm) Post- (mean ± SD, mm) Difference (mean ± SD, mm) Pre- (mean ± SD, mm 2 ) Post- (mean ± SD, mm 2 ) Difference (mean ± SD, mm 2 ) Tooth number Maxilla1113.47 ± 2.4512.88 ± 2.26− 0.59 ± 1.3661.02 ± 18.2049.07 ± 14.46− 11.95 ± 9.500.92 ± 0.69 1212.51 ± 2.1011.51 ± 2.11− 1.00 ± 0.9442.24 ± 11.3733.63 ± 8.97− 8.61 ± 6.821.10 ± 0.99 1313.14 ± 2.0011.80 ± 2.00− 1.34 ± 1.1244.47 ± 14.2238.69 ± 13.09− 5.77 ± 4.761.24 ± 1.16 2113.53 ± 2.4212.54 ± 2.44− 0.99 ± 0.7857.29 ± 19.8845.99 ± 11.70− 11.31 ± 12.931.05 ± 0.83 2212.63 ± 1.7711.55 ± 2.14− 1.08 ± 0.9339.96 ± 11.9534.07 ± 10.64− 5.89 ± 6.721.10 ± 0.98 2313.30 ± 1.6312.11 ± 1.36− 1.19 ± 0.9244.29 ± 8.8337.56 ± 7.55− 6.73 ± 4.931.27 ± 1.00 Mandible3111.77 ± 1.8210.59 ± 1.85− 1.18 ± 1.4032.30 ± 11.3125.97 ± 8.82− 6.33 ± 6.601.27 ± 1.57 3212.41 ± 1.9611.42 ± 2.00− 0.99 ± 1.5735.12 ± 10.8329.88 ± 8.48− 5.24 ± 6.261.24 ± 1.94 3313.42 ± 2.1812.73 ± 1.96− 0.69 ± 0.9244.09 ± 12.1339.35 ± 10.61− 4.74 ± 6.560.72 ± 0.99 4111.70 ± 1.7810.50 ± 1.84− 1.20 ± 1.1632.32 ± 9.8125.75 ± 7.24− 6.57 ± 5.801.30 ± 1.41 4212.39 ± 2.5211.07 ± 1.92− 1.32 ± 1.1635.38 ± 12.4628.83 ± 8.83− 6.55 ± 6.711.49 ± 1.41 4313.79 ± 2.0912.72 ± 2.15− 1.07 ± 0.9247.71 ± 12.4840.79 ± 9.62− 6.92 ± 7.391.11 ± 0.94 mm, millimeter, mm 2 , square millimeter; SD, standard deviation.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y ∆Root length ∆Projected root sur - face area External root resorption ∆Root length ∆Projected root surface area External root resorption Tooth number Tooth movement Tooth number Tooth movement Maxillary central teeth 11 Intrusion r − 0.091 0.051 0.008 21 Intrusion r 0.025 0.119 − 0.050 p 0.537 0.731 0.956 p 0.867 0.419 0.735 Extrusion r 0.051 − 0.381 − 0.326 Extrusion r 0.060 − 0.084 − 0.064 p 0.733 0.008(BH− p = 0.131) 0.024(BH− p = 0.269) p 0.683 0.571 0.665 Translation BL r − 0.149 − 0.135 0.199 Translation BL r 0.095 0.267 − 0.092 p 0.314 0.361 0.176 p 0.521 0.067 0.534 Translation MD r − 0.261 0.014 0.284 Translation MD r − 0.030 − 0.208 0.059 p 0.073 0.924 0.050 p 0.842 0.156 0.691 Crown Tip MD r 0.054 − 0.423 0.254 Crown Tip MD r − 0.224 − 0.245 0.209 p 0.716 0.003(BH− p = 0.084) 0.081 p 0.126 0.094 0.153 Crown Torque BL r − 0.059 − 0.351 0.234 Crown Torque BL r − 0.225 − 0.353 0.214 p 0.690 0.014(BH− p = 0.201) 0.109 p 0.124 0.014(BH− p = 0.201) 0.144 Rotation r 0.161 − 0.071 − 0.059 Rotation r − 0.080 − 0.273 0.067 p 0.274 0.631 0.692 p 0.589 0.060 0.653 Maxillary lateral teeth 12 Intrusion r − 0.202 0.141 0.136 22 Intrusion r − 0.155 0.277 0.205 p 0.168 0.340 0.355 p 0.293 0.057 0.162 Extrusion r 0.123 − 0.161 − 0.033 Extrusion r 0.207 − 0.429 − 0.292 p 0.404 0.275 0.825 p 0.157 0.002 0.044 Translation BL r − 0.115 − 0.100 0.090 Translation BL r 0.300 0.190 − 0.274 p 0.436 0.501 0.543 p 0.038 0.196 0.060 Translation MD r − 0.180 0.167 0.180 Translation MD r 0.020 0.008 − 0.001 p 0.222 0.257 0.221 p 0.894 0.956 0.996 Crown Tip MD r 0.088 − 0.035 − 0.087 Crown Tip MD r − 0.456 − 0.176 0.555 p 0.551 0.811 0.558 p 0.001 0.231 < 0.001(BH− p = 0.006) Crown Torque BL r − 0.504 − 0.284 0.527 Crown Torque BL r − 0.036 0.075 0.127 p < 0.001(BH− p = 0.016) 0.051 < 0.001(BH− p = 0.010) p 0.807 0.611 0.391 Rotation r 0.215 − 0.277 − 0.191 Rotation r 0.040 0.109 0.030 p 0.143 0.057 0.193 p 0.790 0.459 0.842 Table 4 Correlation analysis between resorption parameters and tooth movement amounts for maxillary teeth.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y ∆Root length ∆Projected root sur - face area External root resorption ∆Root length ∆Projected root surface area External root resorption Maxillary canine teeth 13 Intrusion r 0.155 0.090 − 0.244 23 Intrusion r − 0.029 − 0.126 0.019 p 0.294 0.544 0.095 p 0.847 0.395 0.896 Extrusion r − 0.212 − 0.552 0.161 Extrusion r − 0.134 − 0.283 0.128 p 0.148 < 0.001(BH− p = 0.006) 0.273 p 0.362 0.052 0.388 Translation BL r − 0.089 − 0.379 0.123 Translation BL r − 0.015 − 0.027 0.008 p 0.548 0.008(BH− p = 0.131) 0.404 p 0.919 0.855 0.957 Translation MD r 0.146 0.108 − 0.177 Translation MD r − 0.164 − 0.153 0.165 p 0.321 0.464 0.229 p 0.265 0.298 0.261 Crown tip MD r 0.063 − 0.094 − 0.074 Crown tip MD r − 0.049 0.083 0.102 p 0.669 0.526 0.616 p 0.739 0.575 0.489 Crown torque BL r − 0.020 − 0.195 0.025 Crown torque BL r − 0.299 0.154 0.320 p 0.894 0.185 0.867 p 0.039(BH− p = 0.298) 0.296 0.027(BH− p = 0.269) Rotation r 0.220 0.014 − 0.200 Rotation r − 0.156 − 0.185 0.161 p 0.132 0.923 0.173 p 0.291 0.208 0.275 Δ, Difference; BL, Buccolingual; MD, Mesiodistal; r, Pearson correlation coefficient; p, unadjusted p value; BH-p, Benjamini–Hochberg false discovery rate−adjusted p value. BH-p values are reported for correlations with unadjusted p < 0.05. Table 4 (continued)
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y ∆Root length ∆Projected root surface area External root resorption ∆Root length ∆Projected root surface area External root resorption Tooth number Tooth movement Tooth number Tooth movement Man - dibular Central Teeth 31 Intrusion r − 0.293 0.101 0.301 41 Intrusion r − 0.138 − 0.055 0.229 p 0.043(BH− p = 0.300) 0.496 0.038(BH− p = 0.298) p 0.350 0.712 0.117 Extrusion r − 0.059 − 0.368 0.042 Extrusion r − 0.044 0.085 0.016 p 0.690 0.010(BH− p = 0.157) 0.775 p 0.765 0.564 0.915 Translation BL r − 0.401 − 0.334 0.395 Translation BL r − 0.213 − 0.432 0.315 p 0.005(BH− p = 0.112) 0.020(BH− p = 0.254) 0.005(BH− p = 0.113) p 0.146 0.002(BH− p = 0.084) 0.029(BH− p = 0.274) Translation MD r − 0.071 − 0.009 0.002 Translation MD r 0.292 0.081 − 0.289 p 0.631 0.954 0.990 p 0.044(BH− p = 0.300) 0.585 0.046(BH− p = 0.308) Crown Tip MD r − 0.303 − 0.222 0.230 Crown Tip MD r 0.124 0.104 − 0.129 p 0.037(BH− p = 0.298) 0.129 0.116 p 0.399 0.484 0.383 Crown Torque BL r − 0.004 − 0.039 0.002 Crown Torque BL r 0.103 − 0.146 − 0.016 p 0.979 0.790 0.987 p 0.488 0.324 0.911 Rotation r − 0.143 − 0.209 0.185 Rotation r 0.165 − 0.262 − 0.149 p 0.332 0.153 0.208 p 0.263 0.072 0.312 Man - dibular lateral teeth 32 Intrusion r − 0.272 − 0.118 0.280 42 Intrusion r 0.192 0.139 − 0.154 p 0.062 0.423 0.054 p 0.192 0.347 0.296 Extrusion r 0.175 0.248 − 0.143 Extrusion r − 0.108 − 0.145 0.103 p 0.235 0.090 0.332 p 0.465 0.324 0.487 Translation BL r 0.073 − 0.214 0.062 Translation BL r − 0.090 − 0.037 0.031 p 0.621 0.145 0.675 p 0.542 0.800 0.836 Translation MD r 0.160 0.090 − 0.102 Translation MD r − 0.175 − 0.326 0.072 p 0.276 0.544 0.490 p 0.234 0.024(BH- p = 0.269) 0.625 Crown Tip MD r − 0.091 − 0.125 0.212 Crown Tip MD r − 0.054 − 0.386 0.046 p 0.537 0.396 0.148 p 0.716 0.007 (BH− p = 0.130) 0.759 Crown Torque BL r 0.265 0.160 − 0.140 Crown Torque BL r − 0.152 − 0.156 0.072 p 0.069 0.278 0.343 p 0.303 0.290 0.626 Rotation r 0.080 0.071 − 0.046 Rotation r 0.139 0.017 − 0.098 p 0.590 0.632 0.758 p 0.348 0.907 0.506 Table 5 Correlation analysis between resorption parameters and tooth movement amounts for mandibular teeth.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y ∆Root length ∆Projected root surface area External root resorption ∆Root length ∆Projected root surface area External root resorption Tooth number Tooth movement Tooth number Tooth movement Man - dibular canine teeth 33 Intrusion r 0.339 0.139 − 0.321 43 Intrusion r 0.038 0.161 − 0.032 p 0.019(BH− p = 0.247) 0.347 0.026(BH− p = 0.269) p 0.798 0.275 0.828 Extrusion r − 0.176 0.064 0.158 Extrusion r − 0.002 0.078 − 0.045 p 0.232 0.668 0.282 p 0.990 0.597 0.761 Translation BL r 0.304 0.302 − 0.296 Translation BL r 0.062 − 0.144 − 0.046 p 0.036(BH− p = 0.298) 0.037(BH− p = 0.298) 0.041(BH− p = 0.300) p 0.675 0.329 0.757 Translation MD r 0.323 − 0.011 − 0.318 Translation MD r − 0.053 − 0.063 0.066 p 0.025(BH− p = 0.269) 0.942 0.028(BH− p = 0.269) p 0.718 0.668 0.657 Crown Tip MD r 0.098 0.144 − 0.113 Crown Tip MD r 0.199 − 0.061 − 0.249 p 0.509 0.330 0.445 p 0.175 0.678 0.088 Crown Torque BL r 0.004 0.094 − 0.012 Crown Torque BL r − 0.037 − 0.400 0.042 p 0.980 0.524 0.937 p 0.805 0.005(BH− p = 0.112) 0.778 Rotation r 0.016 − 0.102 − 0.013 Rotation r 0.089 − 0.420 − 0.098 p 0.916 0.491 0.930 p 0.548 0.003(BH− p = 0.084) 0.509 Δ, Difference; BL, Buccolingual; MD, Mesiodistal; r, Pearson correlation coefficient; p, unadjusted p value; BH-p, Benjamini–Hochberg false discovery rate-adjusted p value. BH-p values are reported for correlations with unadjusted p < 0.05. Table 5 (continued)
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y primarily associated with increased resorption on the distal root surfaces 49 . Regarding intrusion movements, the findings were not entirely consistent. One study reported increased resorption on the mesial and distal root surfaces 50 , particularly in the apical thirds. Conversely, another study observed greater resorption on the buccocervical root surfaces and in the apical thirds of the roots 51 . In the present study, movement-specific associations with root resorption parameters were observed predominantly in the maxillary anterior teeth. Despite these findings, the magnitude of the observed resorptive changes was generally limited. Therefore, the observed associations should be interpreted in the context of the relatively small magnitude of these changes and should not necessarily be considered indicative of clinically relevant root damage. This may partly reflect the intentional inclusion of non-extraction cases with moderate crowding to maintain relatively standardized treatment conditions. While these criteria reduced clinical heterogeneity, they may also limit the generalizability of the findings to patients requiring more extensive or complex orthodontic movements. In the maxillary lateral incisors, crown torque emerged as the most consistent movement-specific factor, showing associations with both root length change and ERR. In the multivariable analyses, each 1-degree increase in crown torque was associated with an approximately 0.07-mm greater reduction in root length and 0.09-mm greater ERR. Crown tipping was also positively correlated with ERR; however, this association was attenuated in the multivariable model, suggesting that its relationship with ERR may partly overlap with the effects of Table 6 Multivariable linear regression analyses of the associations between orthodontic tooth movements and root resorption parameters. Tooth group/movementBStandard errorβ p t95% CI lowerUpper Maxillary lateral teeth—Δ root length (mm) Intrusion− 0.1150.315− 0.0530.717− 0.364− 0.7470.518 Extrusion0.2450.1690.1330.1541.448− 0.0950.585 Translation BL0.2970.1550.1820.0621.915− 0.0150.609 Translation MD0.0500.1710.0270.7730.290− 0.2950.394 Crown Tip MD− 0.0490.029− 0.1850.094− 1.711− 0.1080.009 Crown Torque BL− 0.0700.033− 0.2490.041*− 2.105− 0.136− 0.003 Rotation0.0320.0130.3130.018*2.4570.0060.058 Model SummaryR 2 = 0.171Adjusted R 2 = 0.095F(7,47) = 5.625 p < 0.001 Maxillary lateral teeth—external root resorption (mm) Intrusion− 0.0430.287− 0.0190.881− 0.150− 0.6200.534 Extrusion− 0.2760.173− 0.1430.117− 1.596− 0.6240.072 Translation BL− 0.3110.155− 0.1810.051− 2.003− 0.6230.001 Translation MD− 0.0110.175− 0.0060.950− 0.063− 0.3640.342 Crown Tip MD0.0600.0310.2140.0571.951− 0.0020.122 Crown Torque BL0.0930.0340.3150.008*2.7610.0250.160 Rotation− 0.0300.015− 0.2800.044*− 2.072− 0.060− 0.001 Model summaryR 2 = 0.190Adjusted R 2 = 0.115F(7,47) = 8.105 p < 0.001 Maxillary canine teeth—Δ projected root surface area (mm 2 ) Intrusion− 2.0201.602− 0.1510.214− 1.260− 5.2431.204 Extrusion− 6.2921.538− 0.415< 0.001*− 4.091− 9.385− 3.198 Translation BL− 1.0521.058− 0.1350.325− 0.994− 3.1811.077 Translation MD0.7020.8520.0990.4140.824− 1.0122.417 Crown Tip MD− 0.0650.104− 0.0360.535− 0.625− 0.2740.144 Crown Torque BL0.0460.1300.0320.7240.356− 0.2150.308 Rotation− 0.0200.082− 0.0310.804− 0.250− 0.1850.144 Model summaryR 2 = 0.198Adjusted R 2 = 0.124F(7,47) = 3.096 p = 0.009 B, Unstandardized regression coefficient; SE, patient-clustered robust standard error; β, standardized regression coefficient; BL, buccolingual; MD, mesiodistal; R 2 , coefficient of determination; t, t statistic. Each regression model included 96 teeth from 48 patients. All models were adjusted for side. Patient-clustered robust standard errors were used to account for dependence between bilateral teeth from the same patient. The F statistic represents the joint test of the seven tooth-movement predictors after adjustment for side. *p < 0.05.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y concomitant tooth movements. In contrast, rotation remained associated with both root length change and ERR in the multivariable analyses after adjustment for the other tooth-movement variables and side; however, these findings should be interpreted cautiously. In cases of moderate crowding, maxillary lateral incisors may be positioned palatally because of space limitations, and their buccal repositioning during OT may require substantial root and crown torque movements 52 . Torque movements are associated with increased stress concentration at the root surface due to their force–moment characteristics, which may contribute to the observed association with ERR. However, the association between torque and root length change should also be interpreted in light of the limitations of panoramic radiography. Buccolingual changes in root orientation may alter the two-dimensional projection of the root and thereby contribute to apparent changes in measured root length. Accordingly, the torque-related findings may reflect both biological responses to orthodontic loading and projection-related effects inherent to two-dimensional imaging. In the multivariable regression analysis, greater extrusion was associated with a decrease in PRSA in the maxillary canines, with an estimated reduction of approximately 6.29 mm² per 1-mm increase in extrusion. ERR is influenced by factors such as the surface area subjected to force and the magnitude of the applied force 6,53 . Although extrusion generally produces a less compressive periodontal environment than intrusion, the association observed in the present study suggests that extrusive movement may still be relevant to root resorption under specific biomechanical conditions. This finding may also reflect the relatively greater extrusion requirements of maxillary canines, potentially related to their frequent supraocclusion. This often results from early loss of primary canines, leading to space loss and subsequent impaction or supraeruption of permanent canines 54 . Taken together, these findings suggest that the observed association may reflect not only the characteristics of the extrusion movement itself but also the specific biomechanical demands and initial positional characteristics of maxillary canines. Orthodontically induced ERR has no definitive treatment; management primarily involves early radiographic detection and minimizing further progression by adjusting force levels or shortening treatment duration when necessary 55 . In the present study, all radiographic evaluations were performed immediately after the completion of OT. Thus, the observed ERR reflects an early-phase status. As demonstrated by Asscherickx et al. 56 , root defects caused by orthodontic mini-implant insertion exhibited substantial repair through cementum deposition. This suggests that some resorptive changes may have the potential to repair biologically, depending on individual patient conditions. Therefore, post-treatment resorptive findings should be interpreted with caution. Endodontic intervention should be guided by a thorough clinical assessment that considers patient symptoms, functional concerns, and serial radiographic evaluations, rather than relying solely on post-treatment images. In cases of suspected advanced resorption, collaboration between orthodontists and endodontists can help support appropriate treatment decisions and improve long-term outcomes. From a clinical perspective, in cases requiring extensive tooth movement—particularly in teeth considered more susceptible to ERR—clinicians may consider increasing the number of aligners to reduce the magnitude of movement per aligner. This approach may help minimize the risk of ERR by reducing excessive or abrupt forces. However, increasing the number of aligners may also prolong overall treatment duration, which itself has been associated with increased root resorption. Therefore, clinicians should carefully balance force magnitude and treatment duration when planning aligner protocols in patients at higher risk of ERR. In such susceptible teeth, regular radiographic monitoring throughout OT may help detect early signs of ERR and allow timely modification of force application. Beyond these clinical implications, the present findings also contribute to the understanding that the biological response to orthodontic treatment may vary according to both tooth type and the specific biomechanics of the planned tooth movement. Several limitations should be considered when interpreting the findings of the present study. One important limitation is the inherent difficulty of isolating the individual contributions of specific tooth movements to ERR, because multiple movements occur simultaneously during routine clinical OT. Nevertheless, the present study provides clinically relevant insights under conditions that closely reflect routine orthodontic practice. Assessing the isolated effects of specific tooth movements on resorption typically requires finite element models or in
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y vitro simulations. While previous animal studies and finite element analyses have explored the resorptive effects of different tooth movements, the present study complements this evidence by examining movement-specific associations within the clinical context of CA treatment, a widely utilized orthodontic procedure. Furthermore, the multivariable analyses showed movement-specific associations with root resorption parameters in the maxillary lateral incisors and canines, providing evidence that may inform future investigations of the biomechanical determinants of ERR during CA treatment. Additionally, as the mandible does not provide a validated stable anatomical reference for superimposition in non-CBCT datasets, mandibular models were oriented relative to the superimposed maxilla, consistent with previously published digital model analyses. Although strict inclusion criteria minimized substantial occlusal or skeletal changes, minor effects of occlusal settling on mandibular measurements cannot be excluded. In addition, all patients in this study were treated with the ClearCorrect system using a single aligner material. As material properties may influence force transmission, outcomes may vary across different systems. Another limitation is the use of panoramic radiographs, although ERR is a three-dimensional process. While cone-beam computed tomography (CBCT) offers more accurate volumetric measurements of these changes 7 , its higher radiation dose limits routine use 57 . Conversely, two-dimensional imaging methods such as panoramic radiographs may overestimate ERR, as shown by Gandhi et al. 58 , who reported a 0.2 mm higher mean resorption value with 2D techniques compared to CBCT. Although periapical radiographs and CBCT are more sensitive, they are not routinely used during OT 59 . In this study, panoramic radiographs were used for both pre- and post- OT evaluations, reflecting standard orthodontic practice and supporting retrospective analysis. Prior ERR studies have similarly utilized panoramic radiographs 11,29,31,60 . Stramotas et al. 61 reported that linear measurements on serial panoramic radiographs may be sufficiently accurate when the angulation of the occlusal plane is kept similar between serial images and the occlusal plane is not tilted by more than approximately 10°. To minimize measurement errors, all radiographs were obtained by a trained assistant using the same panoramic unit with standardized positioning. Nevertheless, minor differences in head and occlusal-plane positioning between serial examinations cannot be completely excluded and may have contributed to variation in linear measurements. Although panoramic imaging is prone to dimensional distortion 62 , the consistent imaging protocol and equipment use throughout the study were intended to minimize measurement variability and improve the comparability of pre- and post-OT measurements. The proportionally greater reductions in PRSA compared with root length should be interpreted in light of the different dimensional and anatomical characteristics of these measurements. Root length is a one-dimensional measure primarily reflecting apicocoronal shortening, whereas PRSA represents the two-dimensional projected area enclosed by the entire root contour; therefore, equivalent proportional changes are not expected. Moreover, changes predominantly affecting the apical region may be more apparent in root length, whereas alterations involving the mesial and distal root contours may produce a greater change in PRSA with comparatively limited apicocoronal shortening. Accordingly, the greater proportional reductions in PRSA do not necessarily indicate a discrepancy with the corresponding changes in root length but rather reflect the complementary aspects of root morphology captured by these two measurements. However, the inherent limitations of two-dimensional imaging must be acknowledged. In particular, buccolingual torque movements may influence the radiographic projection of the root, potentially leading to apparent reductions in measured root length. Therefore, torque-related associations involving linear root measurements should be interpreted with caution when inferring true morphological root resorption. Similarly, PRSA measurements may also be influenced by changes in tooth orientation, including torque and axial rotation, because the projected area of a three-dimensional structure on a two-dimensional image may vary with its orientation. Therefore, PRSA should not be interpreted as a direct three-dimensional measure of morphological root loss. Unlike previous studies, this investigation quantified ERR not only at the root apex but also along the mesial and distal surfaces, using a digital area calculation program to assess changes in projected root surface area 11,31 . Moreover, the analysis was limited to incisors and canines—teeth most susceptible to ERR—although future studies involving other tooth types are needed to generalize the findings.
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y This study has several key strengths. First, the study applied strict inclusion and exclusion criteria to enhance the consistency and comparability of the study population. Second, the analysis incorporated both linear and projected root surface measurements, allowing root changes to be evaluated using complementary two-dimensional metrics. Third, this study evaluated a broad range of orthodontic tooth movements across different anterior tooth groups and examined their associations with multiple root resorption parameters within the clinical context of CA treatment. This comprehensive approach provides clinically relevant exploratory data for understanding movement-specific patterns of root change during CA therapy. Conclusions Within the limitations of the present study, the associations between specific orthodontic tooth movements and ERR-related changes during CA treatment varied according to both movement type and tooth group. Crown torque in the maxillary lateral incisors and extrusion in the maxillary canines showed the most consistent movement-specific associations with root resorption parameters. These findings should be interpreted with caution and warrant confirmation in future prospective studies using three-dimensional imaging and larger cohorts. Author contributions S.K.: Conceptualization, methodology, supervision, data acquisition and interpretation, writingoriginal draft preparation. M.H.B.: Conceptualization, methodology, data acquisition and interpretation, analysis, critically revised the manuscript. All authors gave final approval and agreed to be accountable for all aspects of the work. Funding The authors received no financial support for the research, authorship, and/or publication of this article. Data availability The datasets for the current study are available from the corresponding author on reasonable request. Declarations Ethics approval and consent to participate This study was performed in line with the principles of the Declaration of Helsinki. Informed consent was obtained from all the participants including both patients and the legal guardians of children involved in the present study. Approval was granted by the Clinical Research Ethics Committee of Kutahya Health Sciences University, Kutahya, Türkiye (no: 2024/07-36). Competing interests The authors declare no competing interests. Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit h t t p : / / c r e a t i v e c o m m o n s . o r g / l i c e n s e s / b y - n c - n d / 4 . 0 / . References 1. Topkara, A., Karaman, A. I. & Kau, C. H. Apical root resorption caused by orthodontic forces: A brief review and a longterm observation. Eur. J. Dent. 6, 445–453 (2012). 2. Harris, E. F., Kineret, S. E. & Tolley, E. A. A heritable component for external apical root resorption in patients treated orthodontically. Am. J. Orthod. Dentofac. Orthop. 111, 301–309 (1997). 3. Pinheiro, L. H. M. et al. Genetic variation involved in the risk to external apical root resorption in orthodontic patients: A systematic review. Clin. Oral Investig. 25, 5613–5627 (2021).
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y 4. Weltman, B., Vig, K. W., Fields, H. W., Shanker, S. & Kaizar, E. E. Root resorption associated with orthodontic tooth movement: a systematic review. Am. J. Orthod. Dentofac. Orthop. 137, 462–476 (2010). 5. Årtun, J., Van’t Hullenaar, R., Doppel, D. & Kuijpers-Jagtman, A. M. Identification of orthodontic patients at risk of severe apical root resorption. Am. J. Orthod. Dentofac. Orthop. 135, 448–455 (2009). 6. Sameshima, G. T. & Sinclair, P. M. Predicting and preventing root resorption: Part I. Diagnostic factors. Am. J. Orthod. Dentofac. Orthop. 119, 505–510 (2001). 7. Butsabul, P., Kanpittaya, P. & Nantanee, R. Root resorption in clear aligner treatment detected by CBCT: a systematic review and meta-analysis. Int. Dent. J. 74, 1326–1336 (2024). 8. Fujiyama, K., Honjo, T., Suzuki, M., Matsuoka, S. & Deguchi, T. Analysis of pain level in cases treated with Invisalign aligner: comparison with fixed edgewise appliance therapy. Prog. Orthod. 15, 1–7 (2014). 9. Buschang, P. H., Shaw, S. G., Ross, M., Crosby, D. & Campbell, P. M. Comparative time efficiency of aligner therapy and conventional edgewise braces. Angle Orthod. 84, 391–396 (2014). 10. Li, Y. et al. Prevalence and severity of apical root resorption during orthodontic treatment with clear aligners and fixed appliances: a cone beam computed tomography study. Prog. Orthod. 21, 1–8 (2020). 11. Kurnaz, S. & Buyukcavus, M. H. Panoramic evaluation of external root resorption in mandibular molars during orthodontic treatment: a comparison between root-filled and vital teeth treated with fixed appliances or clear aligners. BMC Oral Health. 24, 1152 (2024). 12. Singh, S., Jain, R. K. & Balasubramaniam, A. CoComparative assessment of external apical rootresorption between subjects treated with clear aligners and fixed orthodontic appliances: Asystematic review and meta-analysis.. J. Dent. Res. Dent. Clin. Dent. Prospects. 18, 85–94 (2024). 13. Darvizeh, A. et al. External apical root resorption following orthodontic treatment with clearaligners versus fixed appliances: A systematic review and meta-analysis. Dent. J. 13, 580 (2025). 14. Selvaraj, M., Nivethitha, B. & Madhan, B. Orthodontically induced external apical rootresorption with clear aligners compared to fixed appliance treatment: An umbrella review. J.Oral Biol. Craniofac. Res.. 15, 938–947 (2025). 15. Kumasako-Haga, T., Konoo, T., Yamaguchi, K. & Hayashi, H. Effect of 8-hour intermittent orthodontic force on osteoclasts and root resorption. Am. J. Orthod. Dentofac. Orthop. 135, 278e1–278e9 (2009). 16. Elhaddaoui, R., Qoraich, H. S., Bahije, L. & Zaoui, F. Orthodontic aligners and root resorption: A systematic review. Int. Orthod. 15, 1–12 (2017). 17. Brezniak, N. & Wasserstein, A. Orthodontically induced inflammatory root resorption. Part I: the basic science aspects. Angle Orthod. 72, 175–179 (2002). 18. Feller, L., Khammissa, R. A., Thomadakis, G., Fourie, J. & Lemmer, J. Apical external root resorption and repair in orthodontic tooth movement: biological events. Biomed. Res. Int. 2016, 4864195 (2016). 19. Segal, G., Schiffman, P. & Tuncay, O. Meta analysis of the treatment-related factors of external apical root resorption. Orthod. Craniofac. Res. 7, 71–78 (2004). 20. Parker, R. J. & Harris, E. F. Directions of orthodontic tooth movements associated with external apical root resorption of the maxillary central incisor. Am. J. Orthod. Dentofac. Orthop. 114, 677–683 (1998). 21. Han, G., Huang, S., Von den Hoff, J. W., Zeng, X. & Kuijpers-Jagtman, A. M. Root resorption after orthodontic intrusion and extrusion: an intraindividual study. Angle Orthod. 75, 912–918 (2005). 22. Nakano, T. et al. Effects of different types of tooth movement and force magnitudes on the amount of tooth movement and root resorption in rats. Angle Orthod. 84, 1079–1085 (2014). 23. Bartley, N. et al. Physical properties of root cementum: Part 17. Root resorption after the application of 2.5 and 15 of buccal root torque for 4 weeks: A microcomputed tomography study. Am. J. Orthod. Dentofac. Orthop. 139, e353–e360 (2011). 24. Wu, A. T. et al. Physical properties of root cementum: Part 18. The extent of root resorption after the application of light and heavy controlled rotational orthodontic forces for 4 weeks: A microcomputed tomography study. Am. J. Orthod. Dentofac. Orthop. 139, e495–e503 (2011). 25. King, A. D. et al. Physical properties of root cementum: Part 21. Extent of root resorption after the application of 2.5 and 15 tips for 4 weeks: A microcomputed tomography study. Am. J. Orthod. Dentofac. Orthop. 140, e299–e305 (2011). 26. Hamilton, R. & Gutmann, J. Endodontic-orthodontic relationships: a review of integrated treatment planning challenges. Int. Endod. J. 32, 343–360 (1999). 27. Flatten, J., Gedrange, T., Bourauel, C., Keilig, L. & Konermann, A. The role of bone and root resorption on the biomechanical behavior of mandibular anterior teeth subjected to orthodontic forces: A finite element approach. Biomedicines 12, 1959 (2024). 28. Zhang, M., Zhang, P., Koh, J.-T., Oh, M.-H. & Cho, J.-H. Evaluation of aligners and rootresorption: An overview of systematic reviews. J. Clin. Med. 13, 1950 (2024). 29. Lee, Y. J. & Lee, T. Y. External root resorption during orthodontic treatment in root-filled teeth and contralateral teeth with vital pulp: A clinical study of contributing factors. Am. J. Orthod. Dentofac. Orthop. 149, 84–91 (2016). 30. Almagrami, I. M. et al. A quantitative three-dimensional comparative study of alveolar bone changes and apical root resorption between clear aligners and fixed orthodontic appliances. Prog. Orthod. 24, 6 (2023).
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y 31. Kurnaz, S. & Buyukcavus, M. H. External root resorption in root-filled and vital teeth after extraction and non-extraction orthodontic treatments: A split-mouth retrospective study. Acta Odontol. Scand. 79, 282–288 (2021). 32. Cangialosi, T. J. et al. The ABO discrepancy index: a measure of case complexity. Am. J. Orthod. Dentofac. Orthop. 125, 270–278 (2004). 33. Garib, D. et al. Superimposition of maxillary digital models using the palatal rugae: does ageing affect the reliability? Orthod. Craniofac. Res. 22, 183–193 (2019). 34. Linge, L. & Linge, B. O. Patient characteristics and treatment variables associated with apical root resorption during orthodontic treatment. Am. J. Orthod. Dentofac. Orthop. 99, 35–43 (1991). 35. Papageorgiou, S. N., Koletsi, D., Iliadi, A., Peltomaki, T. & Eliades, T. Treatment outcome with orthodontic aligners and fixed appliances: A systematic review with meta-analyses. Eur. J. Orthod. 42, 331–343 (2020). 36. Liou, E. J. & Chang, P. M. Apical root resorption in orthodontic patients with en-masse maxillary anterior retraction and intrusion with miniscrews. Am. J. Orthod. Dentofac. Orthop. 137, 207–212 (2010). 37. Acar, A., Canyürek, Ü., Kocaaga, M. & Erverdi, N. Continuous vs. discontinuous force application and root resorption. Angle Orthod. 69, 159–163 (1999). 38. Roscoe, M. G., Meira, J. B. & Cattaneo, P. M. Association of orthodontic force system and root resorption: A systematic review. Am. J. Orthod. Dentofac. Orthop. 147, 610–626 (2015). 39. Ballard, D. J., Jones, A. S., Petocz, P. & Darendeliler, M. A. Physical properties of root cementum: Part 11. Continuous vs intermittent controlled orthodontic forces on root resorption. A microcomputed-tomography study. Am. J. Orthod. Dentofac. Orthop. 136, 8e1–8e8 (2009). 40. Yi, J., Xiao, J., Li, Y., Li, X. & Zhao, Z. External apical root resorption in non-extraction cases after clear aligner therapy or fixed orthodontic treatment. J. Dent. Sci. 13, 48–53 (2018). 41. Guan, W., Lu, Y., Yu-feng, Z., San-lian, L. & Ji-wei, Z. A retrospective study on incisor root resorption in patients treated with bracketless invisible appliance and straight wire appliance. Shanghai Kou Qiang Yi Xue. 26, 121–124 (2017). 42. Phillips, J. R. Apical root resorption under orthodontic therapy. Angle Orthod. 25, 1–22 (1955). 43. DeShields, R. W. A study of root resorption in treated Class II, Division I malocclusions. Angle Orthod. 39, 231–245 (1969). 44. Mirabella, A. D. & Årtun, J. Risk factors for apical root resorption of maxillary anterior teeth in adult orthodontic patients. Am. J. Orthod. Dentofac. Orthop. 108, 48–55 (1995). 45. Currell, S. D., Liaw, A., Grant, P. D. B., Esterman, A. & Nimmo, A. Orthodontic mechanotherapies and their influence on external root resorption: A systematic review. Am. J. Orthod. Dentofac. Orthop. 155, 313–329 (2019). 46. Tieu, L. D., Saltaji, H., Normando, D. & Flores-Mir, C. Radiologically determined orthodontically induced external apical root resorption in incisors after non-surgical orthodontic treatment of class II division 1 malocclusion: A systematic review. Prog. Orthod. 15, 1–13 (2014). 47. Chan, E. & Darendeliler, M. A. Physical properties of root cementum: Part 5. Volumetric analysis of root resorption craters after application of light and heavy orthodontic forces. Am. J. Orthod. Dentofac. Orthop. 127, 186–195 (2005). 48. Barbagallo, L. J., Jones, A. S., Petocz, P. & Darendeliler, M. A. Physical properties of root cementum: Part 10. Comparison of the effects of invisible removable thermoplastic appliances with light and heavy orthodontic forces on premolar cementum. A microcomputed-tomography study. Am. J. Orthod. Dentofac. Orthop. 133, 218–227 (2008). 49. Montenegro, V. C. J., Jones, A., Petocz, P., Gonzales, C. & Darendeliler, M. A. Physical properties of root cementum: Part 22. Root resorption after the application of light and heavy extrusive orthodontic forces: a microcomputed tomography study. Am. J. Orthod. Dentofac. Orthop. 141, e1–e9 (2012). 50. Harris, D. A., Jones, A. S. & Darendeliler, M. A. Physical properties of root cementum: Part 8. Volumetric analysis of root resorption craters after application of controlled intrusive light and heavy orthodontic forces: a microcomputed tomography scan study. Am. J. Orthod. Dentofac. Orthop. 130, 639–647 (2006). 51. Harry, M. & Sims, M. Root resorption in bicuspid intrusion: A scanning electron microscope study. Angle Orthod. 52, 235–258 (1982). 52. Nanda, R. S. & Tosun, Y. S. Biomechanics in Orthodontics. Principles and Practice (Quintessence Publishing Co., 2010). 53. Sameshima, G. T. & Sinclair, P. M. Predicting and preventing root resorption: Part II. Treatment factors. Am. J. Orthod. Dentofac. Orthop. 119, 511–515 (2001). 54. Proffit, W. R., Fields, H. W. & Sarver, D. M. Contemporary Orthodontics, 4th edn (Mosby Elsevier, 2007). 55. Kim, S., Hwang, S. J. & Seo, M. S. Effect of orthodontic treatment on the periapical radiolucency of endodontically treated teeth: A CBCT analysis. BMC Oral Health. 23, 289 (2023). 56. Asscherickx, K., Vannet, B. V., Wehrbein, H. & Sabzevar, M. M. Root repair after injury from mini-screw. Clin. Oral Implants Res. 16, 575–578 (2005). 57. Patel, S., Dawood, A., Wilson, R., Horner, K. & Mannocci, F. The detection and management of root resorption lesions using intraoral radiography and cone beam computed tomography—an in vivo investigation. Int. Endod.J. 42, 831–838 (2009).
Article in Press Scientific Reports https://doi.org/10.1038/s41598-026-77321-y 58. Gandhi, V. et al. Comparison of external apical root resorption with clear aligners and pre-adjusted edgewise appliances in non-extraction cases: A systematic review and meta-analysis. Eur. J. Orthod. 43, 15–24 (2021). 59. Yi, J. et al. Cone-beam computed tomography versus periapical radiograph for diagnosing external root resorption: A systematic review and meta-analysis. Angle Orthod. 87, 328–337 (2017). 60. Llamas-Carreras, J. et al. Study of external root resorption during orthodontic treatment in root-filled teeth compared with their contralateral teeth with vital pulps. Int. Endod. J. 43, 654–662 (2010). 61. Stramotas, S., Geenty, J. P., Petocz, P. & Darendeliler, M. A. Accuracy of linear and angular measurements on panoramic radiographs taken at various positions in vitro. Eur. J. Orthod. 24, 43–52 (2002). 62. Lien, L. C. & Soh, G. Accuracy of the orthopantomogram in assessment of tooth length in orthodontic patients. Singap. Dent. J. 23, 68–71 (2000). Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Authors and Affiliations Safa Kurnaz 1 · Muhammed Hilmi Buyukcavus 2 Safa Kurnaz safa.kurnaz@ksbu.edu.tr 1 Department of Endodontics, Faculty of Dentistry, Kutahya Health Sciences University, 43270 Kutahya, Turkey 2 Department of Orthodontics, Faculty of Dentistry, Antalya Bilim University, Antalya, Turkey
兽研