Factors influencing the development of artificial intelligence in orthodontics.

OBJECTIVES: Since developing AI procedures demands significant computing resources and time, the implementation of a careful experimental design is essential. The purpose of this study was to investigate factors influencing the development of AI in orthodontics. MATERIALS AND METHODS: A total of 162 AI models were developed, with various combinations of sample sizes (170, 340, 679), input variables (40, 80, 160), output variables (38, 76, 154), training sessions (100, 500, 1000), and computer specifications (new vs. old). The TabNet deep-learning algorithm was used to develop these AI models, and leave-one-out cross-validation was applied in training. The goodness-of-fit of the regression models was compared using the adjusted coefficient of determination values, and the best-fit model was selected accordingly. Multiple linear regression analyses were employed to investigate the relationship between the influencing factors. RESULTS: Increasing the number of training sessions enhanced the effectiveness of the AI models. The best-fit regression model for predicting the computational time of AI, which included logarithmic transformation of time, sample size, and training session variables, demonstrated an adjusted coefficient of determination of 0.99. CONCLUSION: The study results show that estimating the time required for AI development may be possible using logarithmic transformations of time, sample size, and training session variables, followed by applying coefficients estimated through several pilot studies with reduced sample sizes and reduced training sessions.

East Asian and Southern European craniofacial class III phenotype: two sides of the same coin?

OBJECTIVES: The skeletal class III phenotype is a heterogeneous condition in populations of different ethnicities. This study aimed to analyse the joint and ethnicity-specific clustering of morphological features in skeletal class III patients of Asian and European origins. MATERIALS AND METHODS: This cross-sectional study involved South Korean and Spanish participants who fulfilled the cephalometric, clinical, and ethnic-related selection criteria. Radiographic records were standardised, calibrated, and measured. A total of 54 skeletal variables were selected for varimax factorial analysis (VFA). Subsequently, a cluster analysis (CA) was performed (mixed method: k-means and hierarchical clustering). Method error and precision were assessed using ICC, Student's t-test, and the Dahlberg formula. RESULTS: A total of 285 Korean and Spanish participants with skeletal class III malocclusions were analysed. After performing VFA and CA, the joint sample revealed three global clusters, and ethnicity-specific analysis revealed four Korean and five Spanish clusters. Cluster_1_global was predominantly Spanish (79.2%) and male (83.01%) and was characterised by a predominantly mesobrachycephalic pattern and a larger cranial base, maxilla, and mandible. Cluster_2_global and Cluster_3_global were mainly South Korean (73.9% and 75.6%, respectively) and depicted opposite phenotypes of mandibular projection and craniofacial pattern. CONCLUSIONS: A distinct distribution of Spanish and South Korean participants was observed in the global analysis. Interethnic and interethnic differences were observed, primarily in the cranial base and maxilla size, mandible projection, and craniofacial pattern. CLINICAL RELEVANCE: Accurate phenotyping, reflecting the complexity of skeletal class III phenotype across diverse populations, is critical for improving diagnostic predictability and future personalised treatment protocols.

Time series analysis.

INTRODUCTION: This article describes a simple method of applying a time series analysis to sample data sets using a free and open statistical software program, Language R. METHODS: Records of new patients who visited 2 different university-affiliated orthodontic departments in 2 different countries were collected. Time series analysis was performed by applying Language R software. The data sets and codes were provided for tutorial and illustrative purposes. RESULTS: Using time series decomposition, the trend component and the seasonal variation were separated and visualized graphically. CONCLUSIONS: Time series analysis may be helpful to clinicians by providing a simple tool to evaluate patient characteristics and manage the practice.

Virtually-Planned Orthognathic Surgery Achieves an Accurate Condylar Position.

PURPOSE: Accuracy in orthognathic surgery with virtual planning has been reported, but detailed analysis of accuracy according to anatomic location, including the mandibular condyle, is insufficient. The purpose of this study was to compare the virtual plan and surgical outcomes and analyze the degree and distribution of errors according to each anatomic location. PATIENTS AND METHODS: This retrospective cohort study evaluated skeletal class III patients, treated with bimaxillary surgery. The primary predictor was anatomic locations that consisted of right and left condyles, maxilla, and the distal segment of the mandible. Other variables were age and gender. The primary outcome was surgical accuracy, defined as mean 3-dimensional distance error, mean absolute error, and mean error along the horizontal, vertical, and anteroposterior axes between the virtual plan and surgical outcomes. Landmarks were compared using a computational method based on affine transformation with a 1-time landmark setting. The mean errors were visualized with multidimensional scattergrams. Bivariate and regression statistics were computed. RESULTS: This study included 52 patients, 26 men and 26 women, with a mean age of 21 years and 3 months. The mean 3D distance errors for condylar landmarks, maxillary landmarks, and landmarks on the distal segment of the mandible were 1.03, 1.25, and 2.24 mm, respectively. Condylar landmarks, maxillary landmarks, and the landmarks on the distal segment of the mandible were positioned at 0.49 mm inferior, 0.28 mm anterior, and 1.25 mm inferior, respectively. The landmark errors for the distal segment of the mandible exhibited a wider distribution than those for condylar and maxillary landmarks. CONCLUSIONS: Agreement between the planned and actual outcome aided by virtual surgical planning was highest for the condyles, followed by the maxilla, and the distal segment of the mandible. It is important to consider the tendency for surgical errors in each anatomic location during operations.

Investigation of Postoperative Complications of Intrabony Cystic Lesions in the Oral and Maxillofacial Region.

PURPOSE: The purpose of this study was to identify factors that can complicate the surgical removal of intrabony cysts and any relevant correlations between them. PATIENTS AND METHODS: The medical records of 249 patients who underwent surgical removal of intrabony cysts were retrospectively reviewed. Outcome variables were postoperative complications, infection, and recurrence. Predictor variables were patient age, gender, comorbidities, anatomic location, pathologic diagnosis, preoperative infection, previous marsupialization, and bone graft methods. Logistic regression analysis was performed to identify risk factors of postoperative infection and recurrence. RESULTS: The cystic lesion was smallest in patients who did not receive a bone graft and increased steadily in those who received a xenogeneic bone graft and an autogenous bone graft, in that order. Paresthesia occurred after enucleation of the cystic lesion in 38 cases. Pathologic fractures were observed in 4 cases. There were 59 postoperative infections. The postoperative infection rate was as high as 63.6% in patients who underwent autogenous bone grafting. In contrast, infection rates were as low as 26.8 and 19.5% in those who underwent xenogeneic bone grafting and no bone grafting, respectively. Location of the cystic lesion in the maxilla or mandible affected the infection rate. When cysts were located in the anterior mandible, no postoperative infection occurred. In contrast, the infection rate was highest for cysts in the mandibular ramus, followed by those in the posterior mandible. Cystic lesion recurrence was observed in 7 cases: 5 cases of odontogenic keratocysts, 1 case of periapical cyst, and 1 case of dentigerous cyst. These findings suggest that cyst pathologic identity affects the recurrence rate. Cyst size was meaningfully correlated with recurrence rate. CONCLUSIONS: These results suggest that autogenous bone grafts increase the risk of postoperative infection compared with absence of a bone graft.

Extensive Surgical Procedures Result in Better Treatment Outcomes for Bisphosphonate-Related Osteonecrosis of the Jaw in Patients With Osteoporosis.

PURPOSE: To identify the risk factors associated with relapse or treatment failure after surgery for bisphosphonate-related osteonecrosis of the jaw (BRONJ) in patients with osteoporosis. PATIENTS AND METHODS: We performed a retrospective cohort study of BRONJ in patients with osteoporosis who had undergone surgical procedures from 2004 to 2016 at the Department of Oral and Maxillofacial Surgery, Seoul National University Dental Hospital. The predictor variables were a set of heterogeneous variables, including demographic (age, gender), anatomic (maxilla or mandible, or both, affected location), clinical (disease stage, etiology, comorbidities, history of intravenous bisphosphonate intake), time (conservative treatment before surgery, bisphosphonate treatment before the development of BRONJ, discontinuation of the drug before surgery, interval to final follow-up, interval to reoperation in the case of relapse or treatment failure), and perioperative variables (type of anesthesia, type of surgical procedures). The primary outcome variable was relapse after surgery that required reoperation (yes vs no). The descriptive and bivariate statistics were computed to assess the relationships between the study variables and the outcome. To determine the risk factors, we conducted a survival analysis using the Cox model. RESULTS: The final sample included 325 subjects with a median age of 75 years, and 97% were women. After surgery, 30% of patients did not completely recuperate and underwent repeat surgery. The interval from the first surgery to reoperation ranged from 10 days to 5.6 years. Relapse or treatment failure most often occurred immediately after surgery. The type of surgical procedure and mode of anesthesia were the most important factors in the treatment outcome. A drug holiday did not appear to influence the likelihood of relapse after surgery. CONCLUSIONS: Treatment of BRONJ in patients with osteoporosis might benefit from more careful and extensive surgical procedures rather than curettage performed with the patient under local anesthesia.

Three-Dimensional Computed Tomographic Assessment of Temporomandibular Joint Stability After Orthognathic Surgery.

PURPOSE: Ensuring that the condyle is appropriately positioned and that positional changes are minimal is considered crucial for the temporomandibular joint (TMJ) to function without symptoms after orthognathic surgery. The purposes of this study were to evaluate condylar changes after surgery and to examine the association between these changes and TMJ symptoms. PATIENTS AND METHODS: A retrospective cohort study was conducted in patients with mandibular prognathism who underwent orthognathic surgery. Linear and angular changes in the positioning of the condyle were measured by superimposing 3-dimensional computed tomograms taken before surgery and 3 months after surgery. Clinical symptoms of TMJ pain and sound were recorded at 3, 6, 9, and 12 months after surgery. Possible associations between TMJ symptoms and clinical variables, such as postoperative condylar changes, were investigated using multiple logistic regression analysis. RESULTS: Linear condylar displacement after orthognathic surgery occurred predominantly in the anterior, medial, and inferior directions, with minimal changes (<1 mm) observed. Most angular condylar changes were smaller than 4° and occurred in the inward direction in the axial plane and the posterior direction in the sagittal plane. The best predictor of postoperative TMJ signs and symptoms was the preoperative status of TMJ signs and symptoms. Neither linear nor angular condylar displacement showed a relevant influence on postoperative pain and sound. CONCLUSIONS: Within the ranges of linear (<1 mm) and angular (<4°) condylar displacement noted in this study, displacement was not associated with postoperative TMJ pain and sound.

Can measurement of apparent diffusion coefficient before treatment predict the response to uterine artery embolization for adenomyosis?

OBJECTIVES: To determine the usefulness of the apparent diffusion coefficient (ADC) in predicting response to uterine artery embolization (UAE) for symptomatic adenomyosis. METHODS: A prospective study was performed on 25 patients who underwent diffusion-weighted (DW) magnetic resonance imaging (MRI) before UAE between June 2011 and December 2012. All patients underwent 3-month follow-up MRI after UAE using non-spherical polyvinyl alcohol (PVA) particles ranging from 150 to 500 µm. Quantitative measurement of the ADC was performed for each adenomyosis. Complete response and incomplete response were defined as ≥ 90 % or < 90 %, respectively, of the non-perfusion area with adenomyosis at the follow-up MRI. The ADC values were compared between patients who achieved complete or incomplete response which was assessed according to the MRI findings after UAE. RESULTS: Nineteen patients showed complete response, and six showed incomplete response. The ADC value ranged from 0.842 to 1.346 × 10(-3) mm(2)/s (mean 1.075 ± 0.117). The mean ADC was 1.043 ± 0.237 in the complete response group and 1.176 ± 0.429 in the incomplete response group (0.012). Using a threshold of <1.147 × 10(-3) mm(2)/s, the sensitivity and specificity of the ADC to predict success after UAE were 83.3 % and 84.2 %. CONCLUSION: The ADC of adenomyosis may potentially predict a successful response to UAE for adenomyosis. KEY POINTS: • Pre-procedural MRI might help clinicians predict response of UAE in adenomyosis • ADC might help predict UAE outcomes in adenomyosis • MR predictors might be used to counsel patients with symptomatic adenomyosis.

How to test validity in orthodontic research: a mixed dentition analysis example.

INTRODUCTION: The data used to test the validity of a prediction method should be different from the data used to generate the prediction model. In this study, we explored whether an independent data set is mandatory for testing the validity of a new prediction method and how validity can be tested without independent new data. METHODS: Several validation methods were compared in an example using the data from a mixed dentition analysis with a regression model. The validation errors of real mixed dentition analysis data and simulation data were analyzed for increasingly large data sets. RESULTS: The validation results of both the real and the simulation studies demonstrated that the leave-1-out cross-validation method had the smallest errors. The largest errors occurred in the traditional simple validation method. The differences between the validation methods diminished as the sample size increased. CONCLUSIONS: The leave-1-out cross-validation method seems to be an optimal validation method for improving the prediction accuracy in a data set with limited sample sizes.

A more accurate soft-tissue prediction model for Class III 2-jaw surgeries.

INTRODUCTION: The use of bimaxillary surgeries to treat Class III malocclusions makes the results of the surgeries more complicated to estimate accurately. Therefore, our objective was to develop an accurate soft-tissue prediction model that can be universally applied to Class III surgical-orthodontic patients regardless of the type of surgical correction: maxillary or mandibular surgery with or without genioplasty. METHODS: The subjects of this study consisted of 204 mandibular setback patients who had undergone the combined surgical-orthodontic correction of severe skeletal Class III malocclusions. Among them, 133 patients had maxillary surgeries, and 81 patients received genioplasties. The prediction model included 226 independent and 64 dependent variables. Two prediction methods, the conventional ordinary least squares method and the partial least squares (PLS) method, were compared. When evaluating the prediction methods, the actual surgical outcome was the gold standard. After fitting the equations, test errors were calculated in absolute values and root mean square values through the leave-1-out cross-validation method. RESULTS: The validation result demonstrated that the multivariate PLS prediction model with 30 orthogonal components showed the best prediction quality among others. With the PLS method, the pattern of prediction errors between 1-jaw and 2-jaw surgeries did not show a significantly difference. CONCLUSIONS: The multivariate PLS prediction model based on about 30 latent variables might provide an improved algorithm in predicting surgical outcomes after 1-jaw and 2-jaw surgical corrections for Class III patients.

Comparison of individualized facial growth prediction models based on the partial least squares and artificial intelligence.

OBJECTIVES: To compare facial growth prediction models based on the partial least squares and artificial intelligence (AI). MATERIALS AND METHODS: Serial longitudinal lateral cephalograms from 410 patients who had not undergone orthodontic treatment but had taken serial cephalograms were collected from January 2002 to December 2022. On every image, 46 skeletal and 32 soft-tissue landmarks were identified manually. Growth prediction models were constructed using multivariate partial least squares regression (PLS) and a deep learning method based on the TabNet deep neural network incorporating 161 predictor, and 156 response, variables. The prediction accuracy between the two methods was compared. RESULTS: On average, AI showed less prediction error by 2.11 mm than PLS. Among the 78 landmarks, AI was more accurate in 63 landmarks, whereas PLS was more accurate in nine landmarks, including cranial base landmarks. The remaining six landmarks showed no statistical difference between the two methods. Overall, soft-tissue landmarks, landmarks in the mandible, and growth in the vertical direction showed greater prediction errors than hard-tissue landmarks, landmarks in the maxilla, and growth changes in the horizontal direction, respectively. CONCLUSIONS: PLS and AI methods seemed to be valuable tools for predicting growth. PLS accurately predicted landmarks with low variability in the cranial base. In general, however, AI outperformed, particularly for those landmarks in the maxilla and mandible. Applying AI for growth prediction might be more advantageous when uncertainty is considerable.

Clinical use of alumina-toughened zirconia abutments for implant-supported restoration: prospective cohort study of survival analysis.

OBJECTIVES: The aim of this prospective cohort study was to compute the long-term clinical survival and complication rates of alumina-toughened zirconia abutments used for implant-supported restorations and to evaluate the effects of several clinical variables on these rates. MATERIAL AND METHODS: From May 1998 to September 2010, 213 patients aged 18 years or older were recruited. The patients received 611 external hex implants and 328 implant-supported fixed restorations using alumina-toughened zirconia abutments. During the follow-up, each restoration was coded as a dental event, which included loosening or fracture of abutment screws, and abutment fracture. From the coded data, the effects of the investigated clinical variables-restored area (anterior/posterior), number of prosthodontic units (one/two units or over), prosthesis type (single-unit/multiunit without pontic/multiunit with pontic), implant system, and patient gender-on the survival of the abutments were evaluated. Survival analysis using Kaplan-Meier method and Cox proportional hazard model was applied. The 5-year survival and complication rates of the abutments were assessed. RESULTS: The number of prosthodontic units and the type of prosthesis had a significant association with complication rates (P < 0.05). Kaplan-Meier survival analysis estimated that the cumulative 5-year complication rate of the abutments used in single restorations was 19.7%. Multiunit-fixed dental prostheses without and with pontics had complication rates of 3.9% and 3.8%, respectively. The 5-year survival rate of the abutments was more than 95%, regardless of the type of prosthesis. CONCLUSIONS: Alumina-toughened zirconia abutments are likely to exhibit excellent long-term survival in clinical use for fixed restorations. Single tooth replacement with the abutment at the molar region may require special care and extra attention.

How to report reliability in orthodontic research: Part 1.

In reporting reliability, duplicate measurements are often needed to determine if measurements are sufficiently in agreement among the observers (interobserver agreement) and/or within the same observer (intraobserver agreement). Some reports are often analyzed inappropriately using paired t tests and/or correlation coefficients. The aim of this article is to highlight the statistical problems of reliability testing using paired t tests and correlation coefficients and to encourage good reliability reporting within orthodontic research. With regard to the complex issue of reliability, a simple and singular statistical approach is not available. However, some methods are better than others. A graphic technique based on the Bland-Altman plot that can be simultaneously applied for both intra- and interobserver reliability will also be discussed.

How to report reliability in orthodontic research: Part 2.

Proper statistical analysis is an absolutely essential tool for both clinicians and researchers attempting to implement evidence-based decisions. When analyzing reliability, statistical graphic representation is the best method. Other previously published error studies of 2-dimensional measurements, such as cephalometric landmarks, have inappropriately applied 1-dimensional approaches, such as linear or angular measurements. The aim of this article is to illustrate a graphic presentation method that can be applied to 2-dimensional data sets. We propose that this technique can show errors in both the x-axis and the y-axis simultaneously and should be used when reporting the reliability of a 2-dimensional data set. Our prediction error analysis of soft-tissue changes after orthognathic surgery will be presented as an example. By using different colors in each ellipse, this method can also identify any between-group differences.

Survival analysis of a miniplate and tube device designed to provide skeletal anchorage.

INTRODUCTION: The aim of this prospective cohort study was to compute the clinical survival and complication rates of a miniplate with a tube device (C-tube) used for orthodontic treatment. METHODS: From August 2003 to May 2012, 217 patients were recruited. They received 341 C-tube miniplates. Some C-tube miniplates were removed because orthodontic treatment ended. Others remained beyond the study period and were recorded as censored data. Survival was classified as a C-tube miniplate that functioned in the mouth regardless of any complications. Success was defined as survival without complications. From the data, the effects of these clinical variables on the survival of the C-tube miniplates were evaluated: sex, age, jaw, placement sites, oral hygiene, tube clearance, inflammation, miniplate shape, number of screws, and length of the fixation screws. Survival analyses using the Kaplan-Meier method and the Cox proportional hazard model were applied. RESULTS: Of the 341 miniplates, 14 failed, and 32 had complications. Two-year survival and success rates were 0.91 and 0.80, respectively. In terms of the simple ratio statistic, this was equivalent to a success rate of 96%. The status of oral hygiene maintenance and the operators' experience had significant associations with the complication rates (P <0.001). CONCLUSIONS: The C-tube miniplate has an advantage in versatility in terms of force application. When placing a miniplate, the most important factor is maintaining good tissue health by means of good oral hygiene. Even with good hygiene, the doctor's experience in performing the flap surgery was the second most important factor for success.

Effect of maxillary impaction on mandibular surgical accuracy in virtually-planned orthognathic surgery: A retrospective study.

Although surgical accuracy has been evaluated in bi-maxillary procedures, few studies have investigated the relationship between maxillary and mandibular accuracy. The present study evaluated the effect of maxillary impaction accuracy on mandibular surgical outcome. This cohort study analyzed skeletal class III patients who underwent planned maxillary impaction in bi-maxillary surgery. The primary predictor was the difference between the virtual plan and surgical outcome in the maxilla, as determined by three-dimensional (3D) and vertical differences. The secondary predictors were the planned 3D distances in the maxilla and mandible. The primary outcome was mandibular surgical accuracy, defined as the difference between the planned and actual outcomes, calculated as 3D Euclidean distance. The study included 73 patients. Increased differences between the planned and actual outcomes in the maxilla were associated with increased differences in the mandible. The post-operative position of the mandible was closer to the planned position when the position of the impacted maxilla was superior than when it was inferior to the planned position. Moving the maxilla closer to the planned position resulted in a more accurate mandibular position. These findings suggest that careful surgical procedures are needed to avoid inferior positioning of the maxilla during maxillary impaction surgery.

A more accurate method of predicting soft tissue changes after mandibular setback surgery.

PURPOSE: To propose a more accurate method to predict the soft tissue changes after orthognathic surgery. PATIENTS AND METHODS: The subjects included 69 patients who had undergone surgical correction of Class III mandibular prognathism by mandibular setback. Two multivariate methods of forming prediction equations were examined using 134 predictor and 36 soft tissue response variables: the ordinary least-squares (OLS) and the partial least-squares (PLS) methods. After fitting the equation, the bias and a mean absolute prediction error were calculated. To evaluate the predictive performance of the prediction equations, a 10-fold cross-validation method was used. RESULTS: The multivariate PLS method showed significantly better predictive performance than the conventional OLS method. The bias pattern was more favorable and the absolute prediction accuracy was significantly better with the PLS method than with the OLS method. CONCLUSIONS: The multivariate PLS method was more satisfactory than the conventional OLS method in accurately predicting the soft tissue profile change after Class III mandibular setback surgery.

A three-dimensional analysis of the perceived proportions of maxillary anterior teeth.

OBJECTIVE: The proportions of the anterior dentition, which is important for excellent esthetics, have been extensively studied, but there have been no 3-dimensional interpretations. This study was conducted to compare real tooth sizes and perceived tooth sizes between different genders and populations and to analyze the effects of 3-dimensional tooth position and alignment. MATERIALS AND METHODS: Complete dental stone casts were prepared for a total of 139 subjects (50 males and 44 females from Korea and 46 females from Japan). Using 3-dimensional scanning and reconstructions, virtual models were constructed and the widths, lengths and rotations of maxillary anterior teeth were measured. Parameters related to the arch form were measured orthographically. Descriptive statistics and ANOVA were performed to determine the differences among the three groups. A regression model was created to interpret the values of 2-dimensional perceived widths with 3-dimensional measurements and other parameters. RESULTS: This study observed differences in the average mesiodistal perceived and real dimensions of the maxillary central incisors between Japanese and Korean females, as well as differences in lateral incisor/central incisor ratios and canine/lateral incisor ratios in the perceived 2-dimensional measurements. There were no differences in individual tooth rotations between groups. The r(2) values of the regression model decreased from the central incisors to the canine. CONCLUSIONS: Several differences were found between Japanese and Korean females and the regression models that used real dimensions, rotations and arch form parameters as independent factors were not sufficient to explain the perceived widths of anterior teeth in the study samples.

Do Class III patients have a different growth spurt than the general population?

INTRODUCTION: Understanding the timing and length of the growth spurt of Class III prognathic patients is fundamental to the strategy of interceptive orthopedic orthodontics as well as to the timing of orthognathic surgery. Consequently, this study was undertaken to determine whether there are any significant differences in the stature growth pattern of Class III subjects compared with non-Class III subjects and the general population. METHODS: Twelve-year longitudinal stature growth data were collected for 402 randomly selected adolescents in the general population, 55 Class III mandibular prognathic patients, and 37 non-Class III patients. The growth data were analyzed by using the traditional linear interpolation method and nonlinear growth functions. The 6 stature growth parameters were measured: age at takeoff, stature at takeoff, velocity at takeoff, age at peak height velocity, stature at peak height velocity, and velocity at peak height velocity. Comparisons in the stature growth parameters and 15 cephalometric variables among the general population, Class III subjects, and non-Class III subjects were made with multivariate analysis. RESULTS: Patients with Class III prognathism did not have different growth parameters compared with Class II subjects or the general population. CONCLUSIONS: This study does not allow meaningful conclusions with regard to the relationship of mandibular size and stature growth pattern. The application of nonlinear growth curves vs the traditional linear interpolation method was also discussed.

Three-dimensional biometric study of palatine rugae in children with a mixed-model analysis: a 9-year longitudinal study.

INTRODUCTION: The palatine rugae have been suggested as stable reference points for superimposing 3-dimensional virtual models before and after orthodontic treatment. We investigated 3-dimensional changes in the palatine rugae of children over 9 years. METHODS: Complete dental stone casts were biennially prepared for 56 subjects (42 girls, 14 boys) aged from 6 to 14 years. Using 3-dimensional laser scanning and reconstruction software, virtual casts were constructed. Medial and lateral points of the first anterior 3 rugae were defined as the 3-dimensional landmarks. The length of each ruga and the distance between the end points of the rugae were measured in virtual 3-dimensional space. The measurement changes over time were analyzed by using the mixed-effect method for longitudinal data. RESULTS: There were slight increases in the linear measurements in the rugae areas: the lengths of the rugae and the distances between them during the observation period. However, the amounts of the increments were relatively small when compared with the initial values and individual random variability. Although age affected the linear dimensions significantly, it was not clinically significant; the rugae were relatively stable. CONCLUSIONS: The use of the palatine rugae as reference points for superimposing and evaluating changes during orthodontic treatment was thought to be possible with special cautions.