Architectural changes in alveolar bone for dental decompensation before surgery in Class III patients with differing facial divergence: a CBCT study

This study aimed to investigate alveolar bone change around mandibular anterior teeth during orthodontic decompensation in patients with skeletal Class III malocclusion and different vertical facial patterns. The records of 29 consecutive Class III patients selected from those pending two-jaw orthognathic surgery were divided into low (≤ 28°), average (30°–37°), and high (≥ 39°) mandibular plane angle (MPA) groups. The DICOM files of CBCT scans and STL files of digital dental models, taken before (T1) and after (T2) presurgical orthodontic treatment, were imported into Dolphin imaging software to reconstruct dentoskeletal images. T1 and T2 images were superimposed and analyzed for bone thickness and height at the level of root apex on each mid-sagittal slice of six mandibular anterior teeth. Differences between T1 and T2 were analyzed by non-parametric tests and mixed-effect model analysis. The results showed that the measurements of alveolar bone height generally decreased after treatment, regardless of MPA. The facial divergence, incisor irregularity, tooth site, treatment time, and change in proclination were identified as the significant factors affecting alveolar bone thickness and height during treatment. The presurgical orthodontic treatment to decompensate mandibular anterior teeth should be very careful in all MPA groups.

Retroclination of the mandibular incisors is a frequent manifestation of dental compensation in adults with mandibular prognathism. In surgical-orthodontic treatment for severe Class III malocclusion, dental decompensation is crucial to maximize surgical correction and normalize stomatognathic function. Inadequate dental decompensation may limit the amount of mandible setback and correction of chin prominence in skeletal Class III patients [1][2][3] .
The craniofacial pattern plays an important role in planning orthodontic treatment mechanics. The cortical bone thickness of dentoavleolar ridge has been reported to be greater in adults with facial hypodivergence than in those with hyperdivergence [4][5][6] . In a previous study, the smallest mandibular cortical bone thickness was noted in patients with hyperdivergent facial pattern, and the mandibular bone thickness differed more than bone height among long-face, short-face, and average-face groups 7 . In another study, Class III patients with high-MPA vertical pattern exhibited the smallest alveolar bone thickness compared to other types of malocclusions with vertical and sagittal jaw discrepancy 8 . Furthermore, the vertical facial pattern has been reported as a dentofacial trait significantly affecting cancellous bone thickness in the alveolar process of the mandibular incisors 9 . Other studies demonstrated that a change in tooth position alters the thickness of labial and lingual cortical plates at the level of the root apex 10 . Significant correlation has been reported between incisor inclination and morphological contour of the alveolar bone in patients with skeletal Class III malocclusion 11,12 . Cone-beam computed tomography (CBCT) images have revealed that the alveolar process of mandibular anterior teeth with more labial inclination exhibit thicker labial spongious bone 13 . Another study reported that the labial inclination of Scientific RepoRtS | (2020) 10:14379 | https://doi.org/10.1038/s41598-020-71126-3 www.nature.com/scientificreports/ the mandibular incisors was positively correlated with labial and total alveolar bone thickness, and negatively correlated with lingual alveolar bone height 14 .
Many previous studies assessed only the condition of the alveolar bone of the mandibular incisors in orthodontic patients because of the limitations of 2D image overlapping in cephalometric analysis. Three-dimensional CBCT images allow detailed measurements of alveolar bone dimensions, including alveolar bone height and thickness, in each mandibular anterior tooth. A CBCT study conducted by Lee et al. showed that the mandibular plane angle was negatively correlated with the buccolingual alveolar bone thickness of the mandibular incisors 15 . Another CBCT study reported alveolar bone loss around the incisors in surgically treated skeletal Class III patients; the vertical alveolar bone level was more reduced in the mandibular incisors compared with maxillary incisors, particularly on the lingual side 16 . Sun et al. 14 investigated the alveolar bone condition around the mandibular incisors in patients with skeletal Class III malocclusion. That study reported that the alveolar bone height, including the labial and lingual plates, decreased significantly after presurgical orthodontic treatment. The incisor root apex was closer to the alveolar lingual plate, whereas the alveolar bone thickness at the apex level was not significantly changed.
In addition to displaying craniofacial structures, CBCT scans enable the visualization of individual teeth and the surrounding alveolar bone in a multiplane, true-to-scale image without distortion or object-dependent nonuniform magnification 11,12,[17][18][19] . However, metallic artifacts may be present in the CBCT images of orthodontic patients because of the radiographic scattering effect caused by metal brackets. Previous studies have shown that integrating a maxillofacial CBCT image with its corresponding laser-scanned digital dental models is suitable for generating a precise 3D dentoskeletal models with clinically acceptable registration error 20,21 . Several CBCT studies revealed that orthodontic decompensation alters the alveolar bone dimensions of the mandibular anterior teeth but without differentiating the influence of vertical facial divergence 14,22,23 .
Although the alveolar bone dimensions of mandibular incisors has been investigated in patients with skeletal Class III malocclusion [24][25][26][27] , no study has included alveolar bone change during presurgical orthodontic treatment in Class III patients with different vertical facial divergences. A previous study compared change in alveolar bone morphology at different treatment stages, whereas traditional cephalometric analysis limits the individual assessment of the mandibular anterior teeth 5 . Few studies have reported the influence of presurgical orthodontic treatment on the alveolar bone thickness and height for each of the mandibular incisors, and mandibular canines have not usually been included. The purpose of this study was to analyze the alveolar bone dimension of the six mandibular anterior teeth in skeletal Class III adult patients with different vertical facial patterns, and to investigate potential factors affecting changes in the alveolar bone dimensions in presurgical orthodontic treatment.

Materials and methods
patients. This retrospective study included 29 consecutive Class III adult patients (16 men and 13 women; mean age: 21.2 ± 3.7 years) who were diagnosed with skeletal Class III malocclusion (ANB < 0°). All of the patients were treated with presurgical decompensation orthodontic treatment and two-jaw orthognathic surgery at National Taiwan University Hospital from January 2017 to December 2018. Complete high-quality records, including CBCT scans and dental models, were obtained before (T1) and at the completion of pre-surgical orthodontic treatment (T2). The study protocol was approved by the institutional review board (201705106RINB). Written informed consent was obtained from all patients. Patients with other congenital anomalies, a history of maxillofacial trauma, periodontal attachment loss, or previous root resorption were excluded from the study. The patients were categorized into three groups of low (< 28°), average (28°-39°), and high (> 39°) mandibular plane angles (MPA), indicating the difference in vertical facial divergence 28 . Table 1 shows the demographic characteristics of the patients in the three groups.
Although these patients were treated in a postgraduate orthodontic clinic, all patients received 0.018 × 0.025″ slot edgewise appliances. The wire sequence in the treatment followed the common principles in biomechanics. www.nature.com/scientificreports/ At the initial stage, 0.016″ Ni-Ti wires were routinely used, followed by a sequential wire change from 0.016″ stainless wire, 0.016 × 0.022″ copper Ni-Ti wires, to 0.016 × 0.022″ stainless steel wires for complete alignment and leveling. The presurgical treatment time was 11 months on average. Complete high-quality records, including CBCT scans and dental models, were obtained before (T1) and at the completion of presurgical orthodontic treatment (T2). The dental models at T1 and T2 were scanned using 3Shape TRIOS (Cart, 3Shape A/S, Copenhagen, Denmark), and the data were exported as STL files. The CBCT images were obtained using a Morita 3D Accuitomo XYZ Slice View Tomograph (Model: MCT-1, Type: EX 1/2 F17) (J. Morita MFG. Corp., Kyoto, Japan). The machine setting parameters were 90 kVp, 5 mA, EFOV scan mode, 35 s, and an isotropic 0.25-mm voxel size. The scan data were exported in DICOM format and further processed using Dolphin Imaging Premium 3D software (Ver. 11.95; Dolphin Imaging and Management Solutions, Chatsworth, CA, USA). The 3D CBCT images were re-oriented using three reference planes based on the nasion, sella turcica, basion, porion, and orbitale landmarks, as described in our previous study 29 . Afterwards, the scanned data of dental models were combined with the CBCT data to reconstruct dentoskeletal models (Fig. 1A). To improve the integration procedures, the brackets on T2 models were removed from the regions of interest. Furthermore, the T1 and T2 dentoskeletal models were superimposed using Dolphin 3D voxel-based superimposition on the mandibular symphysis 30 . Assessment of tooth movement. To quantitatively assess orthodontic tooth movement incurred during treatment, the mesiobuccal cusp tips of the bilateral mandibular first molars and the midpoint between the incisal edges of the mandibular central incisors were used to determine the mandibular occlusal plane (MOP). The plane passing the midpoints of the left and right central incisor edges and perpendicular to the MOP was defined as the vertical reference plane. To analyze the orthodontic movement of mandibular anterior teeth, the T1 and T2 dentoskeletal models were superimposed with registration on the mandibular symphysis 31 . The "transfer orientation" algorithm in Dolphin software was used to analyze difference in tooth position at T1 and T2 by using the MOP at T1 as the reference plane 22 . Difference in the proclination of tooth axis (D_proclination) and difference in the incisal-edge position (D_crown) were analyzed separately for each individual tooth (Fig. 2). The use of the MOP at T1 helped eliminate the influence of occlusion change and alteration of the mandiblular position, which typically occurs during presurgcial orthodontic treatment. Moreover, difference in the root length at T1 and T2 was divided by the root length at T1 and denoted as DTLp (percentage of difference in root length) for evaluation of root resorption incurred during the treatment. Furthermore, several patient-related clinical variables were recorded, including space available as a result of missing tooth, interdental space, and extraction of teeth other than third molars, incisor irregularity index modified from Little's irregularity index (IRI, mild: < 4, moderate: 4-8, severe > 8) 32 , and treatment time. In addition to the variables mentioned previously, three commonly used cephalometric measurements, ANB, Wit value, L1-MP (labial inclination of central incisor) were recorded.

Statistical analysis.
To evaluate measurement error, 10 pairs of CBCT images at T1 and T2 were randomly selected for integration with their digital dental models, and repeatedly analyzed at an interval of 2 weeks. Random errors were calculated according to Dahlberg's formula. Statistical analysis was performed using SAS 9.4 (SAS Institute Inc., Cary, NC, USA). The measurements of the left and right mandibular anterior teeth were averaged and expressed as the mean and standard deviation according to the tooth site (central incisor, lateral incisor, and canine).
The Kruskal-Wallis test and post hoc Dunn's test were performed to check for differences of the measurements before treatment (T1) among the low-, normal-, and high-MPA groups. The Wilcoxon signed-rank test was then used to check for differences between the T1 and T2 measurements. Linear mixed-effect model analysis was  Table 1 shows the distribution and demographic characteristics of the patients according to their vertical facial patterns. No statistically significant difference was observed in terms of age, sex, Wit's value, ANB, L1-MP, space, presence of third molars, incisor irregularity index (IRI), and treatment time among the low-, normal-, and high-MPA groups. The mean measurement error of alveolar bone thickness ranged from 0.15 ± 0.05 mm to 0.21 ± 0.04 mm, and that of alveolar bone height was 0.47 ± 0.04 mm to 0.49 ± 0.05 mm. The random errors in the cephalometric measurements were within 0.5° (range: 0.43°-0.49°). The mean percentage of the difference in root length at T1 and T2 (DTLp) was 2.13% in the central incisors, 2.53% in the lateral incisors, and 1.36% in the canines.

Results
Alveolar bone thickness and height among patients with differing facial divergence. Table 2 shows the alveolar bone measurements of the three groups before treatment (T1). The high-MPA patients had lower alveolar bone thickness and height than the normal-and low-MPA patients did in all mandibular incisors and canines. The lowest values were generally noted in the high-MPA group, followed by the normal-angle group; the highest values were measured in the low-MPA patients ( Supplementary Fig. S1). The Kruskal-Wallis test revealed that the lingual thickness (LT), total thickness (TT), and lingual heights (LH) of alveolar bone at all tooth sites were significantly greater in the low-MPA group than in the high-MPA group. However, the measurements of buccal thickness (BT) and buccal height (BH) did not show any significant group difference at any tooth site.
changes in alveolar bone thickness and height during presurgical decompensation. Table 3 shows the differences in measurements between T1 and T2. The results of the Wilcoxon signed-rank test revealed significant differences in all of the measurements for overall changes at three different tooth sites. The measure-  Potential risk factors identified in univariate and multivariate linear mixed-effect models. Table 4 shows the results of univariate mixed-effect linear models for the respective measurements. The beta coefficient indicated the effect of an individual group relative to the reference group with a coefficient of 1. The significant factors affecting the buccal alveolar bone thickness change (D_BT) included MPA group and tooth axis proclination (D_proclination). Compared to the low-MPA group, the D_BT was − 0.746 mm (p = 0.039) and − 1.162 mm (p = 0.003) in the average-and high-MPA groups, respectively, indicating a decrease in BT in the average-and high-MPA groups compared with the low-MPA group. Moreover, D_BT increased by 0.099 mm for a 1° increase in D_proclination (p < 0.0001). In the multivariate mixed-effect linear model (Table 5), D_BT was − 0.734 mm (p = 0.008) and − 0.947 mm (p = 0.001) for the average-and high-MPA groups, respectively, compared with the low-MPA group, and increased by 0.097 mm for a 1° increase in D_proclination (p < 0.0001). The significant factors affecting lingual alveolar bone thickness (D_LT) included MPA group, treatment time, tooth site, and D_proclination (Table 4). Compared to the low-MPA group, the D_LT increased by 0.966 mm (p = 0.037) in the high-MPA group, and increased by 0.544 mm (p = 0.024) in the canine compared with the central incisors. The increase in treatment time (1-month) induced a small but statistically significant decrease in D_LT (− 0.094 mm, p = 0.044). The increase in D_proclination showed a small but statistically significant decrease in D-LT (− 0.147 mm, p < 0.0001). However, the effects of treatment time and tooth site were non-significant in Table 3. Change in alveolar bone thickness and height after treatment. D_BT, difference in buccal thickness; D_LT, difference in lingual thickness; D_TT, difference in total thickness; D_BH, difference in buccal height; D_LH, difference in lingual height, IQR, inter-quartile range; NS, non-significant in post-hoc inter-group comparison. a Wilcoxon signed-rank test. b Kruskal-Wallis test. #Significant difference between low-angle and high-angle groups. Statistically significant items are indicated in bold.

Measurement
Overall (N = 29) p value a www.nature.com/scientificreports/ the multivariate mixed-effect models (Table 5), which showed an increase of 0.72 mm in D_LT in the high-MPA group (p = 0.043) relative to the low-MPA group, and a decrease of 0.137 mm for 1° increase in D_proclination (p < 0.0001). The significant factors affecting total alveolar bone thickness (D_TT) included IRI, tooth site, and D-proclination (Table 4). The D_TT decreased in the patients with more incisor crowding (− 0.606 mm, severe versus moderate IRI) (p = 0.015). The D_TT increased (0.46 mm, p = 0.028) in the canines compared with the central incisors, whereas it decreased (− 0.041 mm, p = 0.024) with tooth proclination (D_proclination). In the multivariate mixed-effect model, the effect of D_proclination on D_TT was non-significant (Table 5), whereas D_TT showed an increase in the canines (0.417 mm, p = 0.035) relative to the central incisors, and decreased in the patients with severe incisor irregularity (− 0.576 mm, p < 0.020) compared to moderate crowding.

Median (IQR) p value b Low-angle (N = 8) Normal-angle (N = 11) High-angle (N = 10)
The significant factors affecting the changes in alveolar bone height included MPA group and TT for D_BH, and treatment time as well as presence of space for D_LH (Table 4). Compared to the low-MPA group, D_BH www.nature.com/scientificreports/ decreased by − 1.748 mm (p = 0.018) in the average-MPA group, but increased with the initial alveolar bone thickness (0.279 mm, p = 0.039). However, these two variables did not reach statistical significance in the multivariate mixed-effect models (Table 5). Regarding D-LH (Table 4), the increase in treatment time (1-month) induced a significant decrease in D_LH (− 0.318 mm, p = 0.001). The presence of space induced a decrease in D_LH (− 1.799 mm, p = 0.049). In the multivariate mixed-effect models (Table 5), only the treatment time was the significant variable affecting lingual height (D-LH, − 0.284 mm, p = 0.005).

Discussion
In patients with skeletal Class III malocclusion and mandibular prognathism, insufficient incisor decompensation can compromise surgical outcomes. However, extensive tooth movement may be accompanied by a negative effect on the alveolar bone housing of the anterior teeth (e.g., bone dehiscence and fenestration). Special attention should be focused on the alveolar bone around the mandibular anterior teeth in treating Class III patients with mandibular prognathism. Evaluating and identifying the clinical factors associated with changes in alveolar bone dimensions in presurgical orthodontic decompensation is also crucial.
Our present study revealed that in skeletal Class III patients, the alveolar bone dimensions of the mandibular anterior teeth differed among patients with different vertical facial patterns. All of the measurements investigated, including the alveolar bone thickness at the level of root apex (BT, LT, TT) and the alveolar crest height (BH, LH), were smaller in the high-MPA patients than in the low-MPA patients. Interestingly, the group differences were statistically significant on the lingual side (LT, LH), but not on the labial side (BT, BH). After presurgical decompensation, the measurements for LT, TT, BH, and LH decreased regardless of facial divergence, and these changes reached the level of statistical significance. By using linear mixed-effect model analysis, five variables including MPA group, IRI (incisor irregularity index), tooth site, treatment time, and change in proclination (D_proclination) were identified as significantly predictive of alveolar bone measurements after treatment. The significant variables for each measurement item were MPA group and D_proclination for D_BT, MPA group and D_proclination for D_LT, tooth site and IRI for D-TT, and treatment time for D_LH. These results implied that presurgical orthodontic treatment should be meticulously and efficiently provided for the mandibular anterior teeth, particularly in patients with high MPA.
In patients with limited space deficiency of the mandibular incisors, labial tipping movement can help gain space for tooth alignment and maximize anterior crossbite for mandibular setback. Therefore, non-extraction therapy is a common strategy for mandibular dentition in the presurgical orthodontic treatment of Class III patients. However, more crowding of the mandibular anterior teeth requires more labial tipping of the incisors, Table 5. Multivariate linear mixed-effect models for change of alveolar bone thickness and height after treatment (adjusted models). D_BT, difference in buccal thickness; D_LT, difference in lingual thickness; D_ TT, difference in total thickness; D_BH, difference in buccal height; D_LH, difference in lingual height. www.nature.com/scientificreports/ possibly leading to a higher risk of alveolar bone loss. Previous studies on anterior alveolar bone dimensions included patients with mild crowding in the mandibular arch 14,15,22,23 . In the current study, the severity of space deficiency of the mandibular anterior teeth was categorized according to the irregularity index modified from Little's irregularity index 17 . Multivariate mixed-effect models showed that the IRI was significantly predictive of the change in total alveolar bone thickness (D_TT). The total alveolar bone thickness after treatment was smaller in patients with severe IRI than those with moderate IRI (Table 5). Significant associations have been identified among gingival biotypes, tooth inclination, and alveolar plate thickness surrounding mandibular incisors in skeletal Class III patients 33 . Conventional cephalometric study reported that vertical facial pattern did influence presurgical decompensation of the mandibular incisors in skeletal Class III malocclusion 5 . Previous CBCT studies showed presurgical orthodontic treatment-related alveolar bone loss around the incisors in skeletal Class III patients 14,16 . Forward movement of the mandibular incisors during presurgical decompensation causes the root apex of the incisors to be closer to the lingual alveolar plate and leads to loss of vertical alveolar height 14,16 . A previous study reported significant decreases in labial height (1.43 and 1.58 mm) and lingual height (0.94 mm and 0.98 mm) for central and lateral incisors, respectively 22 . The present study also revealed significant decreases in labial height (D_BH, − 2.70, − 3.45, − 0.6 mm) and lingual height (D_LH, − 1.0, − 1.35, − 1.28 mm) for the mandibular central incisors, lateral incisors, and canine. However, no significant difference in D_BH or D_LH was noted among the low-, normal-and high-MPA groups.
Our study demonstrated the significant treatment changes in total alveolar bone thickness at the level of the root apex of the six mandibular anterior teeth. These results conflict with those of a previous study that reported no significant change in total alveolar bone thickness 14 . Another finding of decreasing bone thickness on the labial and lingual aspects is consistent with other studies 16,22 . Geometrically, uncontrolled labial tipping of mandibular anterior teeth increases tooth axis proclination, leading to an increase in BT and decrease in LT. This implies that bone thickness decreased on the side toward which a root had moved and increased on the opposite side. In other words, the incisor root apex is closer to the lingual side of alveolar bone after presurgical orthodontic treatment. Therefore, labial proclination of mandibular anterior teeth would lead to an increase in BT, whereas a decrease in BT might imply alveolar bone loss incurred during orthodontic decompensation. In the current study, a significant correlation between D_proclination and D_BT (Pearson correlation coefficient = 0.753) was noted in the low-MPA group, but not in the average-and high-MPA groups. These findings are consistent with the results of multivariate analysis (Table 5, adjusted mixed-effect models), which showed the significant effect of MPA group and tooth proclination on the variation of D_BT and D_LT. An increase in decompensation of 1° (D_proclination) for labial tipping increased BT by 0.097 mm (D_BT) and reduced LT by 0.0137 mm (D_LT). The influence on the lingual thickness was statistically significant, but may not be clinically significant.
According to the multivariate analysis, the tooth site, incisor irregularity index, and change in proclination were significant in affecting the total alveolar bone thickness after decompensation (D_TT). Surprisingly, none of the investigated variables significantly influenced the buccal bone height after treatment (D_BH), and the change in LH (D_LH) was significantly affected only by the treatment time. These results implied that some factors not investigated in the present study, such as lip pressure and gingival biotype, may influence the buccal bone height during treatment. In presurgcial orthodontics for Class III malocclusion, more crowding and lingual inclination of the mandibular incisors typically requires more dental decompensation, thus leading to an increase in treatment time. After adjusting the effect of space, the decrease in lingual alveolar bone height was 0.284 mm if the treatment time was increased by 1 month (Table 5, adjusted models).
Alveolar plate thickness is affected by a variety of factors, e.g. tooth position, tooth inclination, and gingival biotype 33,34 . Previous CBCT studies revealed the significant correlation between the changes of tooth inclination and alveolar bone thickness in orthodontic patients 10,22 . Our study revealed that the alveolar height of mandibular anterior teeth was lower on the lingual side than on the labial side before treatment, and decreased on both the labial and lingual sides after decompensation. Previous studies have reported that in skeletal Class III patients, the lingual inclination of the mandibular incisors is significantly correlated with vertical facial divergence, and is smaller in high-MPA patients than in normal-or low-MPA patients 5,15,26 . The patients included in the current study exhibited severe skeletal Class III malocclusion (ANB: − 11.8° (− 19.6° to − 1.9°), and thus their mandibular incisors exhibited obvious lingual tipping (L1-MP: 78.8° (68.6°-90.1°)). All of the patients in this study were treated with non-extraction therapy for mandibular dentition. The extraction of mandibular third molars can often provide space for distalization of the posterior teeth to help relieve crowding in the middle segment of the mandibular arch. The adjusted mixed-effect models revealed that the IRI was significantly predictive of TT after treatment (D_TT), indicating a decrease in total alveolar bone thickness induced by mandibular anterior teeth with more crowding ( Table 5). The multivariate analysis also showed that the presurgcial orthodontic treatment time was the only significant variable affecting treatment change in lingual alveolar bone height (D_LH). Surprisingly, none of the investigated variables were significantly predictive of the labial alveolar bone height in multivariate mixed-effect models (Table 5). In this study, all of the mandibular anterior teeth were moved labially with an increased labial proclination of the tooth axis. The effect of tooth site on the total alveolar thickness after treatment (D_TT) was also statistically significant, with an increase of 0.417 mm in the canines relative to the central incisors, which might have resulted from measurements of the bone thickness at the anatomic midpoint of the canines, a site lateral to the direction of decompensation, not the direct compression zone in the direction of movement (as was the case for the incisors).
The limits of alveolar housing and the risk of developing bone dehiscence have to be taken account when planning presurgical orthodontic treatment. CBCT is increasingly used as a tool for linear measurements of dental-alveolar structures, even there is a risk of overestimating fenestrations and dehiscence 35 . The presence of fenestration/ dehiscence was most prevalent in the mandibular incisors and canines 36,37 . Moreover, labial malposition of canines is frequently seen in the mandibular dentition with lingually inclined and crowded incisors. In addition to appropriate biomechanical design and periodontal health maintenance, thorough assessment of Scientific RepoRtS | (2020) 10:14379 | https://doi.org/10.1038/s41598-020-71126-3 www.nature.com/scientificreports/ pre-treatment alveolar bone features is also crucial for the prevention of further bone loss during treatment. Distalization of dentition to avoid excessive labial tipping of the mandibular anterior teeth or augmentation of the bone housing should be considered, particularly in patients with high facial divergence. Lateral ridge augmentation with corticotomy of the mandibular anterior teeth increases the labial bone thickness directly, thus allowing a greater amount of labial movement, and lessening the decrease in alveolar bone height, implying the potential of reducing the risk of dehiscence/ fenestration 38,39 . In addition to accelerating tooth movement, this procedure can protect alveolar bone housing from the risk of gingival recessions when substantial proclination of the mandibular incisors is needed. This retrospective study had some limitations. The investigation of treatment changes focused only on hard tissue, not the overlying soft tissue. The assessment of alveolar bone thickness was limited to the midsagittal plane at the level of the root apex. The periodontal health of the mandibular anterior teeth could be more comprehensively assessed by combining clinical observation of gingival recession and topographic examination of CBCT for bone fenestration and dehiscence.

conclusion
The alveolar bone thickness and height of the mandibular anterior teeth were smaller in the high-MPA patients than in the low-MPA patients. All of the lingual thickness, total thickness, labial height, and lingual height of the alveolar bone significantly decreased after presurgcial orthodontic treatmnet, regardless of the extent of facial divergence. The multivariate linear mixed-effect models revealed five factors significantly affecting alveolar bone dimensions during decompensation, including facial divergence, incisor irregularity, tooth site, treatment time, and change in proclination. Collectively, the results of this study provide additional evidence to help more comprehensive consideration for the surgical-orthodontic treatment of Class III patients. The pre-surgical orthodontic treatment to decompensate mandibular anterior teeth should be approached very carefully in all MPA groups.