Dynamic changes in tooth displacement and bone morphometry induced by orthodontic force

Scientific Reports volume 12, Article number: 13672 (2022) Cite this article

1460 Accesses

1 Citations

Metrics details

This study used a novel 3D analysis to longitudinally evaluate orthodontic tooth movement (OTM) and bone morphometry. Twelve-week-old male Wistar rats were subjected to OTM by applying a constant orthodontic force (OF) of 25cN between one of the upper first molars and a mini-screw. In vivo micro-CTs were taken before and after 10, 17, 24 and 31 days of force application, and superimposed by a novel and rigid voxel-based registration method. Then the tooth and alveolar bone segment at different time points became comparable in the same coordinate system, which facilitated the analysis of their dynamic changes in 3D. By comparison between time points and between OF and no OF sides, this study showed that the OTM rate was not constant through time, but conformed to a ‘V’ shape changing pattern. Besides, OF induced displacement of both loaded and unloaded teeth, and the latter mirrored the former in a delayed manner. In addition, bone morphometric changes synchronized with OTM rate changes, implying that a higher OTM rate was concomitant with more alveolar bone loss. The pressure and tension areas might not be in two opposite sides, but actually adjacent and connected. These findings might provide instructive evidence for both clinical, translational and basic research in orthodontics.

Shortening orthodontic treatment time is of interest to both clinicians and patients1. As a result, numerous techniques and products have been commercialized in the last years aiming to accelerate orthodontic tooth movement (OTM)2,3. However, due to the high methodological heterogeneity of OTM evaluation, there is still a lack of knowledge regarding the physiology of tooth movement4,5. Basic research in orthodontics often focuses on a microscopic level, such as on the role of specific molecules, signaling pathways or gene expression patterns involved in OTM6,7,8,9,10. However, macroscopic changes induced by orthodontic force are less studied, despite being highly clinically relevant. This can partially be due to the difficulties inherent to this research model: studies in human would involve biopsies or increased radiological exposure only for research purposes, which cannot be justified. Orthodontic animal models, mostly developed in rodents, offer important opportunities. However, research derived from them often presents contradictory results, probably due to high methodological heterogeneity, such as the use of different orthodontic anchorage, forces and appliances.

Classically, OTM has been divided into different phases according to the rate of tooth displacement. The initial phase would be followed by an arrest period (lag phase), after which an acceleration of OTM would occur, followed by a linear increase in tooth displacement11,12,13. However, studies usually measure this displacement as a distance, which cannot reflect the real 3D situation, since displacement is actually a vector. Second, the dynamic change in alveolar bone when applying a constant orthodontic force remains largely unknown. Although previous animal research showed changes in the bone mineral density (BMD) of alveolar bone subjected to orthodontic force14,15, few studies performed complete longitudinal bone morphometric assessment16,17. Whether there is a correlation between the dynamic changes in OTM rate and the changes in bone morphometric parameters still needs to be elucidated.

The aim of this study is to investigate the dynamic changes induced by orthodontic force in tooth displacement and bone morphometry in a rat OTM model over 31 days by using a novel 3D analysis and bone anchorage. These dynamic changes can provide important evidence regarding the optimal window periods to investigate microscopic molecular changes. In addition, the potential correlation between OTM rate and bone morphometry at different time points after the application of orthodontic force could directly provide instructive guidelines not only for basic research, but also for the clinicians.

Twelve young adult male Wistar rats (9 weeks old) were used in this study. Weekly surveillance was performed in order to notice any changes in weight and guarantee normal feeding. Prior to the experiment, the minimum sample size was calculated using a previous split-mouth study18,19. The primary parameter considered was the mesial displacement of the first molar after 31 days of orthodontic force application. A power analysis in the software G*Power 3.1 (Düsseldorf, Germany) suggested a minimum sample size of 12 animals per group (OF and no OF side) for a repeated measures MANOVA when assuming 99% power with α = 0.05.

The animals were housed in three cages under continuous temperature (23 °C), a 12-h schedule alternating light and dark, and standard rat maintenance diet. All animal experiments were performed in the Laboratory Animal Center and the Molecular Small Animal Imaging Center (MoSAIC) of KU Leuven, Belgium, with the approval of KU Leuven Ethical Committee for Animal Experimentation (P197/2019) and in accordance with the EU Directive 2010/63/EU and ARRIVE 2.0. Guidelines.

Tooth movement was performed by following a previously published protocol18. Briefly, on one randomly-allocated side of the maxilla, a self-drilling mini-screw (2.5 mm × 1.3 mm × 5 mm, DEWIMED, Tuttlingen, Germany) was implanted in the alveolar bone, approximately 2 mm distal to the upper incisors and with an angulation of 45° to the occlusal plane, in order to provide bone anchorage. Healing time of 3 weeks was allowed to ensure stability20. Afterwards, a constant orthodontic force (OF) of 25 cN was loaded between the mini-screw and the upper first molar by using a Sentalloy closed coil spring (Ultra-light, Dentsply GAC, Rochecorbon, France) (Fig. 1A). At the contralateral side, no OF was applied (no OF side). At this moment, the animals were 12 weeks old.

OTM Animal model. (A) Graphic representation of the animal model. A constant force of 25 cN was loaded between by a spring coil between the first molar and mini-screw in the OF side. The rats were scanned using a low-dose in vivo micro-CT at five different time points (T0–T4). (B) Overview of the image processing workflow. Micro-CT scans from different time points (T1–T4) were spatially aligned with the micro-CT on baseline (T0). A volume of interest (VOI) was delineated on the reconstructed micro-CT data to ensure an accurate selection of teeth and alveolar bone. These were segmented using a global threshold algorithm. Based on the resulting binary images an automatic 3D morphometric quantification was performed and 3D models were created.

All interventions were conducted by the same investigator, firstly under sedation of 2.5–5% isoflurane (1000 mg/g, Iso-Vet, Dechra, Skipton, UK), followed by anesthesia of 100 mg/ml IP ketamine (80 mg/kg, Nimatek, Bladel, Netherlands) and 2% xylazine (10 mg/kg, XYL-M, V.M.D, Arendonk, Belgium). After the interventions, soft diet and analgesic medication (0.05 mg/kg Buprenorphine) were supplied for 3 days. To prevent unnecessary animal suffering and potential appliance loss, daily surveillance of body weight and intraoral examination under sedation with isoflurane was performed. The rat grimace scale21 was used to evaluate whether animals were in pain. A schematic diagram showing treatments of the rats from arrival to euthanasia is shown in Supplementary Fig. 1.

Animals were longitudinally followed up with micro-CT at five time points: right before the application of OF (T0, baseline), and 10 (T1), 17 (T2), 24 (T3), and 31 (T4) days after OF. A low-dose, high-resolution in vivo micro-CT (Skyscan 1278, Bruker, Kontich, Belgium) was used for image acquisition. A high-resolution scan protocol was used at 65 kVp, 500 μA, and 180° with an angular rotation step of 0.5°, resulting in an exposure time of 50 ms. A 1 mm aluminum filter was used to eliminate the beam hardening effect. Flat field correction was performed for calibration based on an empty field of view previous to the actual scans. Animals were placed under sedation of 2.5–5% isoflurane during image acquisition. After image acquisition, the image stacks were reconstructed with NRecon software (version 1.7.1, Bruker, Kontich, Belgium). Correction for post-alignment, and ring-artifacts reduction were optimized per scan if needed. Smoothing level and beam hardening were applied with values of 0 and 10%, respectively. The dynamic image range of histogram was set from 0.003 to 0.03.

The dynamic changes of OTM were volumetrically evaluated based on the spatial displacement of the teeth from T0 to T4, using T0 as a baseline without depending on other external reference structures. Rigid voxel-based registration steps were performed before the evaluation to ensure that all teeth were in the same coordinate system and were temporospatially comparable at each time point.

First, the follow-up micro-CT scans at T1, T2, T3 and T4 were manually superimposed with the corresponding baseline scan at T0 based on the maxillary structures as a reference, followed by an optimized automatic superimposition by MTM Scaffold Strain (KU Leuven, Leuven, Belgium)22,23. The structural compatibility of maxillary anatomical reference points was inspected to verify the validity of the manual and automatic superimposition. Second, the maxilla, mini-screw, first and second molars at each time point were delineated by the same investigator and saved as volumes of interest (VOI) in CTAnalyser software (version 1.17.5, Bruker, Kontich, Belgium). Third, the selected VOIs were automatically segmented using an adaptive threshold algorithm22,23 and imported as individual 3D standard tessellation language (stl) models, which were loaded in 3-Matic (Materialise, Leuven, Belgium) for longitudinal assessment of OTM and bone morphometry. The complete workflow is shown in Fig. 1B.

To assess first and second molar displacement, six reference points on the cusps and five on the root apices were created. Their displacement over the different time points was defined as occlusal and apical movement, respectively. The occlusal plane was defined by the mesial cusp, distobuccal and distolingual points. The angle between the occlusal planes was defined as angular movement. The occlusal, apical and angular movements of the first and second molars at T0–T1, T0–T2, T0–T3, T0–T4, T1–T2, T2–T3, T3–T4 were measured on both OF and no OF side. The OTM rate was calculated by the occlusal movement divided by the time of OTM. The longitudinal measurements of the above-mentioned parameters and their dynamic changes are shown in Fig. 2.

Longitudinal assessment of OTM. (A) In the T0 model, occlusal, apical points and occlusal plane were defined. Then, the same points were extrapolated into the T1–T4 models by automated global registration. The distances or angles between them were then measured as occlusal, apical or angular movement. (B) Dynamic changes in molar displacement and alveolar bone morphometry. (Red = OF side, Blue = No OF side).

For assessment of mini-screw displacement, two reference points were created on the center of the screw head and the screw tip in the T0 model. The line connecting these two points was defined as the central axis. Then, the distances of the head and tip points from T0 to T1 were measured as the head and tip movement in the T0–T1 period, respectively. The angle between the central axes in T0 and T1 was measured as the angular movement. In this way, the displacement of mini-screw at T0–T2, T0–T3, T0–T4, T1–T2, T2–T3, T3–T4 was also measured.

Alveolar bone morphometry was quantified using the method proposed by Chatterjee et al., as shown in Fig. 324. The registered micro-CT images were transferred to CTAnalyser and a volume of interest (VOI) containing only the molar region of each hemi maxilla was selected by defining a top and a bottom slice. The top slice was defined as the slice 2 mm mesial from that where the cusp of the first molar appeared. The bottom slice was defined as that where the mesial cusp of the third molar appeared. The alveolar bone within this region was semi-automatically selected and segmented using an automatic global threshold algorithm.

Longitudinal assessment of bone morphometry. The VOIs of the alveolar bone, first and second molar were separately selected and segmented as 3D models on the registered micro-CT in a workflow. The bone morphometry was evaluated in VOIs at each time points and compared longitudinally.

To evaluate the BMD, phantoms (standard hydroxyapatite blocks) of 0.25 g/cm3 and 0.75 g/cm3 were scanned to perform BMD calibration with respect to the attenuation values. The BMD at the different time points was calculated by linear extrapolation using: Y \(-{y}_{1}= ({y}_{2}-{y}_{1})/({x}_{2}-{x}_{1}) \times\) (X \(-{x}_{1}\)), where x1 and x2 are the greyscale indices of standard hydroxyapatite 0.25 and 0.75 g/cm3, y1 and y2 are known as 0.25 and 0.75 g/cm3, and X and Y are the grey indices.

The dynamic changes of the following parameters were also evaluated in VOIs for the assessment of bone morphometry in CTAnalyser:

Bone volume fraction (BV/TV, %): the proportion of the VOI occupied by binarised solid objects in 3D within the VOI, which directly reflects the bone volume.

Bone surface density (BS/TV, mm−1): the ratio of surface area to total volume in 3D within the VOI.

Trabecular number (Tb.N, mm−1): the number of traversals across a trabecular or solid structure made per unit length on a random linear path through the VOI.

Trabecular thickness (Tb.Th, mm): essentially the thickness of the solid voxels as defined by binarisation within the VOI.

Trabecular separation (Tb.Sp, mm): essentially the thickness of the spaces among trabecular bone.

Two-way repeated measures MANOVA and post-hoc Tukey’s test were used to compare the OTM and bone morphometry between different time points (T0–T4) and between the OF and no OF side. Non-parametric analysis (Friedman test) was performed to compare the displacement of the mini-screw among different time points. The correlation of the dynamic changes between bone morphometric parameters and OTM rate was evaluated by Pearson correlation coefficients. Non-parametric statistical methods were used when normality was not confirmed using the Kolmogorov–Smirnov test. Statistical analysis was performed in GraphPad (version 8.4.3, San Diego, USA).

No relevant weight loss was detected in the animals throughout the experiments.

Significant differences in occlusal, apical and angular movement were found among the different time points (P \(<\) 0.01) for both the first (Fig. 4) and second molars (Fig. 5) in the OF side. The OTM rate of the first molar was not constant among the different time points, but conformed to a ’V’ shape changing pattern: it first decreased from T1 to T2 (P \(<\) 0.01), and then increased from T2 to T4 (P \(<\) 0.01).

First molar displacement. The OTM was compared between different time points (T0–T4) and between the OF and no OF side by two-way repeated measures MANOVA and post-hoc Tukey’s test (n = 12). *P < 0.01 versus no OF side. †P < 0.01 versus previous period.

Second molar displacement. The OTM was compared between different time points (T0–T4) and between the OF and no OF side by two-way repeated measures MANOVA and post-hoc Tukey’s test (n = 12). *P < 0.01 versus no OF side. †P < 0.01 versus previous period.

Remarkably, the first and second molars of the no OF side also underwent significant displacement (P \(<\) 0.01) after 31 days of force application on the contralateral hemi-maxilla. Although of course significantly less OTM was found compared to the OF side (P \(<\) 0.01), the changing pattern of OTM rate of the second molar at the OF side and of both molars at the no OF side was similar to that of the loaded first molar.

No mini-screw loosened or was fully lost during the follow-up period. However, mini-screw displacement was observed, shown by box plots (minimum, Q1, median, Q3, maximum) in Fig. 6. The non-parametric analysis (Friedman test) shows the mini-screw had a significantly different displacement on the head (P \(<\) 0.01) but not on the tip (P \(>\) 0.05) over time.

Mini-screw displacement. The displacement of the mini-screw was compared among different time points by Friedman test (n = 12). *P < 0.05.

Alveolar bone morphometry was quantified and compared between the OF and no OF sides and among different time points (T0–T4) by two-way repeated measures ANOVA. The dynamic changes in bone morphometry are shown in Fig. 7.

Dynamic changes in bone morphometry. The bone morphometric parameters were compared between different time points (T0–T4) and between the OF and no OF side by two-way repeated measures MANOVA and post-hoc Tukey’s test (n = 12). BMD = bone mineral density, BV/TV = bone volume fraction, BS/TV = bone surface density, Tb.Th = trabecular thickness, Tb.N = trabecular number, Tb.Sp = trabecular separation. *P < 0.01 versus no OF side. †P < 0.01 versus previous period.

When comparing the OF and no OF sides, no significant difference in bone mineral density (BMD, g cm−3) was observed (P \(>\) 0.05). However, bone volume fraction, bone surface density and trabecular number (BV/TV, %, BS/TV, mm−1 and Tb.N, mm−1) decreased significantly in the OF side compared to the no OF side (P \(<\) 0.05). In contrast, trabecular thickness (Tb.Th, mm) increased significantly (P \(<\) 0.05). Trabecular separation (Tb.Sp, mm) had no significant difference between both sides (P \(>\) 0.05).

When comparing the different time points, BMD and trabecular separation had no significant change (P \(>\) 0.05). But BS/TV and Tb.N decreased significantly (P \(<\) 0.05) while Tb.Th increased significantly (P \(<\) 0.01) with time. Additionally, the dynamic changes of BV/TV, BS/TV, Tb.Th and Tb.N conformed to a ‘V’ shape changing pattern, showing a linear correlation with the change of OTM rate (P \(<\) 0.05) in Fig. 8.

Agreement between the changes in OTM rate and bone morphometry. The OTM rate negatively correlated with bone volume and trabecular number, and positively correlated with trabecular thickness (n = 12, P \(<\) 0.05). Besides, there is no significant linear correlation between the OTM rate and BMD or trabecular separation (P \(>\) 0.05, not shown in this figure).

From digital calipers to intraoral scanners, from histology to cephalograms and from 2 to 3D, the methodology to measure in vivo OTM has been developing for years19,25,26,27,28,29,30. Although measurements have become more and more accurate, the evaluation of the dynamic changes induced by OF still remains difficult, as it is hard to find a reliable reference structure. Previous studies in rats evaluated 2D or 3D OTM in an indirect way based on external reference structures such as the contralateral maxilla (split-mouth designs)19,31, the incisors28,32, the zygomatic arch17 or the diastema formed between the first and the second molar10,15,33,34,35,36. However, these methods rely on those reference structures remaining stable during OTM. The present study uses a novel method to longitudinally evaluate dynamic OTM regardless of external reference structures18. After rigid voxel-registration steps based on the maxilla, the scans at different time points become comparable throughout time in the same coordinate system. Thus, the tooth becomes its own self-control for OTM assessment.

Another novelty of this study lies on the dynamic measurements. First, OTM was truly evaluated in 3D instead of 2D or pseudo 3D (using a reconstructed cephalogram from CT). This is a clear advantage since real OTM is a compound displacement of vertical, transversal and sagittal movement. Second, occlusal, apical and angular movement were evaluated separately, instead of performing an average, since not only the compound body movement but also tipping and torque are crucial to orthodontists. Third, OTM was calculated both relative to T0 and between successive time points, in order to register OTM evolution. Since the true 3D displacement of the same tooth did not necessarily always happen in the same direction or plane (Supplementary Fig. 2), the evolution of OTM should be evaluated not simply by subtracting the OTM towards T0 of the adjacent periods (for example T0–T1 and T0–T2), but by directly measuring the OTM between each adjacent period (T1–T2). This represents a step forward from previous cephalometric, geometric and optometric measurements.

In the present study, 3D analysis of first molar displacement allowed us to study its pattern. After an initial phase of movement, a marked decrease in OTM rate was noted 10 days after the application of force, followed by linear increase 1 week later. This matches the reports of classical articles regarding the phases of tooth movement11,12,13. Since the observation time of our study was 31 days and most related articles only include observation times up to 4 weeks, no statements can be made regarding the possible changes happening beyond this point.

The most surprising finding of our study regards the displacement of unloaded teeth. The second molar was found to significantly displace in a mesial direction in the OF side, although no OF was directly applied on it. This phenomenon confirms the unreliability of choosing the diastema formed between first and second molars as a reference point for OTM assessment. Similarly, significant movement of the first and second molars in the no OF side was also found, following the same pattern described in the OF side. However, this pattern is in contradiction to that reported by previous studies in which the OTM in the no OF side was considered to be zero37. A possible explanation is that the 2D measurements used in these studies may have missed displacement happening in planes other than sagittal.

Considering the fact that rat molars migrate distally during aging38,39, the tooth mesialization in the no OF side found in the present study cannot be due to physiologic movement. It could be an adaptive response to the change in dental occlusion and microscopic biologic environment caused by OF in the contralateral side, which is supported by the fact that the OTM in the no OF side mirrors the pattern of that in the OF side18 (Fig. 4, 5). The pattern of OTM of the second molar in the OF side also mirrored that of the loaded first molar but started with a short delay, probably due to the time consumed for occlusal change or the stretching from the periodontal fibers. Altogether, these findings show that the OF applied in one tooth can indeed influence the other teeth, even on the contralateral side, which can therefore no longer be regarded as “control” since tooth displacement can be camouflaged, leading to systematic underestimation of real OTM.

Besides tooth displacement, alveolar bone remodeling is crucial to clinicians. Although BMD has been previously studied and can reflect bone changes15,40,41, this does not allow for differential analysis of bone volume and mineralization density42, which together with bone volume fraction, bone surface density, trabecular thickness, number and separation are important aspects of bone morphometry, as they can reflect alveolar bone remodeling23,43. In the present study, significantly less bone volume (BV/TV, BS/TV) was found in the OF side compared to the no OF side, but, relevantly, there was no significant difference in BMD between sides. This implies that the bone resorption caused by the aseptic inflammation of OTM44 does coexist with bone formation, maintaining the average BMD. This is in contrast with what happens in pathological processes such as osteoporosis45 or periodontitis46,47,48, where both bone volume and BMD decrease. In correspondence to the changes in bone volume, the trabecular number was also lower in the teeth/bone on the OF side, which also presented with higher trabecular thickness. Longitudinally, both bone volume and trabecular number decreased from T0 to T4, with a constant increase in trabecular thickness as a compensation. Similarly to BMD, the trabecular separation was not significantly different between sides. Although the changes of bone morphometry compared to baseline were statistically significant, the evolution of these changes was not, probably due to the sample size of this study. To the best of our knowledge, no previous study had longitudinally assessed bone morphometry after OTM based on such a rigid voxel-based registration method. Due to this, we calculated the sample size from previous research, based on OTM instead of on bone morphometric parameters.

Longitudinal follow-up of the sample shows that the differences between OF and no OF sides are not due to accidental errors, since they were measured and visualized at each time point, and the difference was constant. Interestingly, the curves representing the changes in bone morphometry over time show a trend (Fig. 7), contrary to the ‘V’ shape changing pattern of OTM rate (Fig. 4). Linear regression analysis showed a significant correlation between the changes in OTM and BV/TV, BS/TV, Tb.N (positively) and Tb.Th (negatively). Although the correlation does not imply causation, this finding suggests that higher OTM rate was concomitant with more alveolar bone loss, which is clinically very relevant. Finding a balance between accelerating OTM and minimum alveolar bone resorption is crucial for future research and clinical practice.

It should be noted that the method of analysis is crucial to the validity of the conclusion, since results of bone morphometry depend on the VOI. Previous studies selected VOIs either in function of the pressure/tension sides16 or anatomic sites such as at the furcation49,50 or mesial root15,51 of the first molar, in different jaws and animal species52,53. However, these areas are difficult to be determined accurately. Furthermore, since OTM is not only a simple mesial body movement (also angular movement, and not necessarily only in the sagittal direction), the pressure and tension areas are not simply the mesial and distal root areas (Supplementary Fig. 3). In the present study, the whole alveolar bone around the first and second molar was selected as VOI, which should simultaneously include pressure and tension sites, since choosing only one of the two might lead to biased BMD results. On the other hand, the relatively wide VOI might reduce the statistical power, because alveolar bone not influenced by the movement might also be selected and therefore hide the change in affected bone. However, this could be compensated by the longitudinal design of the study, since the ‘non-affected’ alveolar bone is eliminated during the subtraction of the different time points.

As mentioned before, the pattern of movement of the first molar in the OF side and all the other molars was found to be similar, but delayed. It could be argued that the displacement of unloaded molars was an adaptive result to that of the loaded first molar. Since the bone morphometric changes were synchronized with the changes in OTM rate, the changes in alveolar bone could not be the cause neither the result, but concomitant with OTM. According to the classical pressure-tension theory, the mesial body movement of the tooth results from mesial alveolar bone resorption and distal bone apposition. In this case however, there should also be a delay between the critical point in the changing pattern of OTM rate and bone morphometry, which was not observed in our sample. This reinforces the idea, already presented in literature54, that the pressure and tension areas are not in two opposite sides, but actually dispersed, adjacent and widely connected in both mesial and distal sides (Supplementary Fig. 3), with bone resorption and formation happening simultaneously. If this is correct, maybe a tilting movement, rather than bodily displacement should be the goal in clinical orthodontics, since it could yield more efficient OTM. This discrete tilting, involving minimal tipping, could be induced by forces lower than the ones currently used clinically. Future research should focus on how to take the best advantage of the tension and pressure areas to optimize total tooth movement, taking into consideration factors such as the force magnitude or root resorption. The longitudinal 3D analysis presented in the present study can offer methodological support to future research.

This study longitudinally assesses OTM and bone morphometry in 3D with a novel method in rats. Tooth displacement was not constant at the different time periods, but conformed to a ‘V’ shape changing pattern. Although no OF was applied on the second molar, this tooth, as well as both molars at the no OF side, underwent significant displacement, in a similar pattern as that of the loaded first molar. The dynamic assessment of bone morphometry demonstrated that a higher OTM rate was concomitant with more alveolar bone loss, which provides instructive evidence for both clinical, translational and basic research in OTM.

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Fu, T., Liu, S., Zhao, H., Cao, M. & Zhang, R. Effectiveness and safety of minimally invasive orthodontic tooth movement acceleration: A systematic review and meta-analysis. J. Dent. Res. 98(13), 1469 (2019).

Article CAS PubMed Google Scholar

MacDonald, L. et al. Efficacy and safety of different interventions to accelerate maxillary canine retraction following premolar extraction: A systematic review and network meta-analysis. Orthod. Craniofac. Res. 24(1), 17 (2021).

Article PubMed Google Scholar

Gao, J. et al. The significance of utilizing a corticotomy on periodontal and orthodontic outcomes: A systematic review and meta-analysis. Biology (Basel). 10(8), 803 (2021).

PubMed PubMed Central Google Scholar

Miles, P. & Fisher, E. Assessment of the changes in arch perimeter and irregularity in the mandibular arch during initial alignment with the AcceleDent Aura appliance vs no appliance in adolescents: A single-blind randomized clinical trial. Am. J. Orthod. Dentofac. Orthop. 150(6), 928 (2016).

Article Google Scholar

Arqub, S. A. et al. The effect of the local administration of biological substances on the rate of orthodontic tooth movement: A systematic review of human studies. Prog. Orthod. 22(1), 5 (2021).

Article PubMed PubMed Central Google Scholar

Deng, L., Chen, Y., Guo, J., Han, X. & Guo, Y. Roles and mechanisms of YAP/TAZ in orthodontic tooth movement. J. Cell. Physiol. 236(11), 7792–7800 (2021).

Article CAS PubMed Google Scholar

Li, Y., Jacox, L. A., Little, S. H. & Ko, C.-C. Orthodontic tooth movement: The biology and clinical implications. Kaohsiung J Med Sci. 34(4), 207 (2018).

Article PubMed Google Scholar

Behm, C. et al. MMPs and TIMPs expression levels in the periodontal ligament during orthodontic tooth movement: A systematic review of in vitro and in vivo studies. Int. J. Mol. Sci. 22(13), 6967 (2021).

Article CAS PubMed PubMed Central Google Scholar

Li, Y., Zhan, Q., Bao, M., Yi, J. & Li, Y. Biomechanical and biological responses of periodontium in orthodontic tooth movement: Up-date in a new decade. Int. J. Oral Sci. 13(1), 20 (2021).

Article PubMed PubMed Central Google Scholar

Klein, Y. et al. Immunorthodontics: In vivo gene expression of orthodontic tooth movement. Sci. Rep. 10(1), 8172 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Pilon, J. J., Kuijpers-Jagtman, A. M. & Maltha, J. C. Magnitude of orthodontic forces and rate of bodily tooth movement. An experimental study. Am. J. Orthod. Dentofac. Orthop.: Off. Publ. Am. Assoc. Orthod. Const. Soc. Am. Board Orthod. 110(1), 16 (1996).

Article CAS Google Scholar

van Leeuwen, E. J., Maltha, J. C. & Kuijpers-Jagtman, A. M. Tooth movement with light continuous and discontinuous forces in beagle dogs. Eur. J. Oral Sci. 107(6), 468 (1999).

Article PubMed Google Scholar

Krishnan, V. & Davidovitch, Z. E. Cellular, molecular, and tissue-level reactions to orthodontic force. Am. J. Orthod. Dentofac. Orthoped. 129(4), 469 (2006).

Article Google Scholar

Campos, M. J. et al. The role of orthodontic tooth movement in bone and root mineral density: A study of patients submitted and not submitted to orthodontic treatment. Med Sci Monit. 18(12), Cr752 (2012).

Article PubMed Google Scholar

Wang, C., Cao, L., Yang, C. & Fan, Y. A novel method to quantify longitudinal orthodontic bone changes with in vivo micro-CT data. J Healthc. Eng. 2018, 1651097 (2018).

PubMed PubMed Central Google Scholar

Nogueira, A. V. B. et al. Experimental models of orthodontic tooth movement and their effects on periodontal tissues remodelling. Arch. Oral Biol. 130, 105216 (2021).

Article CAS PubMed Google Scholar

Fleissig, O. et al. Analytical methodology to measure periodontal bone morphometry following orthodontic tooth movement in mice. Eur. J. Orthod. 43(6), 665 (2021).

Article PubMed Google Scholar

Cadenas de Llano-Pérula, M., Zong, C., Van Dessel, J., Kuijpers-Jagtman, A. M. & Willems, G. 3D quantification of in vivo orthodontic tooth movement in rats by means of micro-computed tomography. Clin. Oral Investig. 26(5), 3911–3920 (2022).

Article PubMed Google Scholar

Kirschneck, C. et al. Differentiated analysis of orthodontic tooth movement in rats with an improved rat model and three-dimensional imaging. Ann. Anat. 195(6), 539 (2013).

Article PubMed Google Scholar

Zhang, J.-N., Lu, H.-P., Bao, X.-C., Shi, Y. & Zhang, M.-H. Evaluation of the long-term stability of micro-screws under different loading protocols: A systematic review. Braz. Oral Res. 33, e046 (2019).

Article PubMed Google Scholar

Oliver, V. et al. Psychometric assessment of the Rat Grimace Scale and development of an analgesic intervention score. PLoS ONE 9(5), e97882 (2014).

Article ADS PubMed PubMed Central CAS Google Scholar

Van Dessel, J. et al. Accuracy and reliability of different cone beam computed tomography (CBCT) devices for structural analysis of alveolar bone in comparison with multislice CT and micro-CT. Eur. J. Oral Implantol. 10(1), 95 (2017).

PubMed Google Scholar

Van Dessel, J. et al. Quantification of bone quality using different cone beam computed tomography devices: Accuracy assessment for edentulous human mandibles. Eur. J. Oral Implantol. 9(4), 411 (2016).

PubMed Google Scholar

Chatterjee, M. et al. A robust methodology for the quantitative assessment of the rat jawbone microstructure. Int. J. Oral Sci. 9(2), 87 (2017).

Article PubMed PubMed Central Google Scholar

Zhou, Y. et al. Mechanism of catechol-O-methyltransferase regulating orofacial pain induced by tooth movement. Biomed. Res. Int. 2021, 4229491 (2021).

PubMed PubMed Central Google Scholar

Ye, J., Zhang, P. Y., Guan, Z. Q., Wang, G. X. & Kou, B. Exploration of effect of Odanacatib on inhibiting orthodontic recurrence in rats and on CatK and IGF-1 mRNA. Eur. Rev. Med. Pharmacol. Sci. 23(8), 3151 (2019).

CAS PubMed Google Scholar

Zheng, J. & Yang, K. Clinical research: Low-level laser therapy in accelerating orthodontic tooth movement. BMC Oral Health 21(1), 324 (2021).

Article CAS PubMed PubMed Central Google Scholar

Ortega, A. J., Campbell, P. M., Hinton, R., Naidu, A. & Buschang, P. H. Local application of zoledronate for maximum anchorage during space closure. Am. J. Orthod. Dentofac. Orthop. 142(6), 780 (2012).

Article Google Scholar

Proff, P. et al. Local vascularization during orthodontic tooth movement in a split mouth rat model—A MRI study. Biomedicines. 8(12), 632 (2020).

Article CAS PubMed Central Google Scholar

Qi, J. et al. Effect of a DPP-4 inhibitor on orthodontic tooth movement and associated root resorption. Biomed. Res. Int. 2020, 7189084 (2020).

PubMed PubMed Central Google Scholar

Trelenberg-Stoll, V., Drescher, D., Wolf, M. & Becker, K. Automated tooth segmentation as an innovative tool to assess 3D-tooth movement and root resorption in rodents. Head Face Med. 17(1), 3 (2021).

Article CAS PubMed PubMed Central Google Scholar

Madan, M. S., Liu, Z. J., Gu, G. M. & King, G. J. Effects of human relaxin on orthodontic tooth movement and periodontal ligaments in rats. Am. J. Orthod. Dentofac. Orthop.: Off. Publ. Am. Assoc. Orthod. Const. Soc. Am. Board Orthod. 131(1), 8.e1 (2007).

Article Google Scholar

Zhou, J. et al. Dynamic evaluation of orthodontically-induced tooth movement, root resorption, and alveolar bone remodeling in rats by in vivo micro-computed tomography. Med Sci Monit. 24, 8306 (2018).

Article CAS PubMed PubMed Central Google Scholar

Xu, X., Zhou, J., Yang, F., Wei, S. & Dai, H. Using micro-computed tomography to evaluate the dynamics of orthodontically induced root resorption repair in a rat model. PLoS ONE 11(3), e0150135 (2016).

Article PubMed PubMed Central CAS Google Scholar

Wang, Y. et al. Macrophages mediate corticotomy-accelerated orthodontic tooth movement. Sci Rep. 8(1), 16788 (2018).

Article ADS PubMed PubMed Central CAS Google Scholar

Lira Dos Santos, E. J. et al. Orthodontic tooth movement alters cementocyte ultrastructure and cellular cementum proteome signature. Bone 2021, 116139 (2021).

Article CAS Google Scholar

Kirschneck, C., Bauer, M., Gubernator, J., Proff, P. & Schröder, A. Comparative assessment of mouse models for experimental orthodontic tooth movement. Sci. Rep. 10(1), 12154 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Tsuchiya, S. et al. Physiological distal drift in rat molars contributes to acellular cementum formation. Anat. Rec. (Hoboken). 296(8), 1255 (2013).

Article PubMed Google Scholar

Sicher, H. & Weinmann, J. P. Bone growth and physiologic tooth movement. Am. J. Orthod. Oral Surg. 30(3), C109 (1944).

Article Google Scholar

Pu, H. & Hua, Y. Hydrogen sulfide regulates bone remodeling and promotes orthodontic tooth movement. Mol. Med. Rep. 16(6), 9415 (2017).

Article CAS PubMed PubMed Central Google Scholar

Chang, H.-W. et al. Effects of orthodontic tooth movement on alveolar bone density. Clin. Oral Investig. 16(3), 679 (2012).

Article PubMed Google Scholar

Ott, S. M. Chapter 2—Bone Biology and Effects of Pharmaceutical Intervention on Bone Quality (2017), 29.

Trelenberg-Stoll, V., Wolf, M., Busch, C., Drescher, D. & Becker, K. Standardized assessment of bone micromorphometry around teeth following orthodontic tooth movement: A µCT split-mouth study in mice. J. Orofac. Orthop. = Fortschritte der Kieferorthopadie: Organ/official journal Deutsche Gesellschaft fur Kieferorthopadie https://doi.org/10.1007/s00056-021-00336-9 (2021).

Article Google Scholar

Chaushu, S., Klein, Y., Mandelboim, O., Barenholz, Y. & Fleissig, O. Immune changes induced by orthodontic forces: A critical review. J. Dent. Res. 2021, 220345211016285 (2021).

Google Scholar

Kim, H. J. et al. In ovariectomy-induced osteoporotic rat models, BMP-2 substantially reversed an impaired alveolar bone regeneration whereas PDGF-BB failed. Clin. Oral Investig. 25(11), 6159 (2021).

Article PubMed Google Scholar

Brasil, S. C. et al. Influence of a high-fat diet in the progression of apical periodontitis. J Endod. 47(4), 600 (2021).

Article PubMed Google Scholar

Muramatsu, R. et al. Guanabenz inhibits alveolar bone resorption in a rat model of periodontitis. J. Pharmacol. Sci. 147(3), 294 (2021).

Article CAS PubMed Google Scholar

Deng, J. et al. A novel modified-curcumin promotes resolvin-like activity and reduces bone loss in diabetes-induced experimental periodontitis. J. Inflamm Res. 14, 5337 (2021).

Article CAS PubMed PubMed Central Google Scholar

Pereira, L. J. et al. Aerobic and resistance training improve alveolar bone quality and interferes with bone-remodeling during orthodontic tooth movement in mice. Bone 138, 115496 (2020).

Article CAS PubMed Google Scholar

Lee, S.-H. et al. Effect of nicotine on orthodontic tooth movement and bone remodeling in rats. Korean J. Orthod. 51(4), 282 (2021).

Article PubMed PubMed Central Google Scholar

Li, H. et al. The expression of extracellular matrix metalloproteinase inducer (EMMPRIN) in the compression area during orthodontic relapse. Eur. J. Orthod. 42(3), 347 (2020).

Article PubMed Google Scholar

Zhou, Y., He, X. & Zhang, D. Study of bone remodeling in corticotomy-assisted orthodontic tooth movement in rats. J. Cell Biochem. 120(9), 15952 (2019).

Article CAS PubMed Google Scholar

Machibya, F. M. et al. Effects of bone regeneration materials and tooth movement timing on canine experimental orthodontic treatment. Angle Orthod. 88(2), 171 (2018).

Article PubMed Google Scholar

d’Apuzzo, F. et al. Biomarkers of periodontal tissue remodeling during orthodontic tooth movement in Mice and Men: Overview and clinical relevance. Sci. World J. 2013, 105873 (2013).

Article CAS Google Scholar

Download references

Chen Zong was supported by the China Scholarship Council (File No. 201806270252). Greetje Vande Velde acknowledged funding by KU Leuven IF (C24/17/061, STG/15/024) and the Flemish Research Foundation (FWO; 1506114N, G057721N). We would like to thank Stephanie de Vleeschauwer, Marie Roelandt and Jordi Penedo for their help during the animal experiments, and Steffen Fieuws for his assistance in statistics.

Department of Oral Health Sciences-Orthodontics, KU Leuven and Dentistry, University Hospitals Leuven, Kapucijnenvoer 7, blok A, bus 7001, 3000, Leuven, Belgium

Chen Zong, Guy Willems & Maria Cadenas de Llano-Pérula

Oral and Maxillofacial Surgery, University Hospitals Leuven and OMFS-IMPATH Research Group, Department of Imaging & Pathology, Faculty of Medicine, KU Leuven, Leuven, Belgium

Jeroen Van Dessel

Biomedical MRI/Molecular Small Animal Imaging Center (MoSAIC), Department of Imaging & Pathology, Faculty of Medicine, KU Leuven, Leuven, Belgium

Greetje Vande Velde

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

All authors contributed to the conception and development of the research. C.Z., as the principal investigator, designed and organized the whole study, analyzed data and wrote the scientific manuscript. J.V.D., G.V.V., and M.C.D.L-P. participated in data analysis and interpretation. G.W. contributed to data interpretation and writing the manuscript. All authors read and approved the final manuscript.

Correspondence to Chen Zong.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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 http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Zong, C., Van Dessel, J., Vande Velde, G. et al. Dynamic changes in tooth displacement and bone morphometry induced by orthodontic force. Sci Rep 12, 13672 (2022). https://doi.org/10.1038/s41598-022-17412-8

Download citation

Received: 26 January 2022

Accepted: 25 July 2022

Published: 11 August 2022

DOI: https://doi.org/10.1038/s41598-022-17412-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.