Marietta Krüsi,Demetrios J. Halazonetis,Theodore Eliades,Vasiliki Koretsi
Covariance patterns between ramus morphology and the rest of the face: A geometric morphometric study
Mandibular ramus,Face,Morphometrics,Covariance
Objective The growth and development of the mandible strongly depend on modeling changes occurring at its ramus. Here, we investigated covariance patterns between the morphology of the ramus and the rest of the face. Methods Lateral cephalograms of 159 adults (55 males and 104 females) with no history of orthodontic treatment were collected. Geometric morphometrics with sliding semi-landmarks was used. The covariance between the ramus and face was investigated using a two-block partial least squares analysis (PLS). Sexual dimorphism and allometry were also assessed. Results Differences in the divergence of the face and anteroposterior relationship of the jaws accounted for 24.1% and 21.6% of shape variation in the sample, respectively. Shape variation was greater in the sagittal plane for males than for females (30.7% vs. 17.4%), whereas variation in the vertical plane was similar for both sexes (23.7% for males and 25.4% for females). Size-related allometric differences between the sexes accounted for the shape variation to a maximum of 6% regarding the face. Regarding the covariation between the shapes of the ramus and the rest of the face, wider and shorter rami were associated with a decreased lower anterior facial height as well as a prognathic mandible and maxilla (PLS 1, 45.5% of the covariance). Additionally, a more posteriorly inclined ramus in the lower region was correlated with a Class II pattern and flat mandibular plane. Conclusions The width, height, and inclination of the ramus were correlated with facial shape changes in the vertical and sagittal planes.
Covariance patterns between ramus morphology and the rest of the face: A geometric morphometric study

The growth and development of the mandible strongly depend on modeling changes occurring at its ramus. Here, we investigated covariance patterns between the morphology of the ramus and the rest of the face.

Lateral cephalograms of 159 adults (55 males and 104 females) with no history of orthodontic treatment were collected. Geometric morphometrics with sliding semi-landmarks was used. The covariance between the ramus and face was investigated using a two-block partial least squares analysis (PLS). Sexual dimorphism and allometry were also assessed.

Differences in the divergence of the face and anteroposterior relationship of the jaws accounted for 24.1% and 21.6% of shape variation in the sample, respectively. Shape variation was greater in the sagittal plane for males than for females (30.7% vs. 17.4%), whereas variation in the vertical plane was similar for both sexes (23.7% for males and 25.4% for females). Size-related allometric differences between the sexes accounted for the shape variation to a maximum of 6% regarding the face. Regarding the covariation between the shapes of the ramus and the rest of the face, wider and shorter rami were associated with a decreased lower anterior facial height as well as a prognathic mandible and maxilla (PLS 1, 45.5% of the covariance). Additionally, a more posteriorly inclined ramus in the lower region was correlated with a Class II pattern and flat mandibular plane.

The width, height, and inclination of the ramus were correlated with facial shape changes in the vertical and sagittal planes.

Growth of the craniofacial complex is an extremely complicated process involving structural and spatial changes in all its units in different developmental periods and to different extents, rates, and directions. Historically, the mandibular condyle has been considered a primary growth center, which, under the influence of intrinsic factors, can control mandibular growth and displace the mandible inferiorly and anteriorly, thus defining the sagittal relationship of the jaws. In contrast, the theory of a functional matrix accounts for the mandibular condyle as a locus where only secondary, compensatory periosteal growth occurs, and mandibular growth is attributed to the morphologic effects of both capsular and periosteal matrices. According to studies on animal and human mandibles, it is common knowledge that the mandible grows forward by bone deposition at the posterior border of the ramus rather than the anterior mandibular region. Simultaneously, the anterior border of the ramus is resorptive; therefore, the body of the mandible elongates. Following the aforementioned modeling of the ramus, coupled with modeling changes at the condylar neck, the mandible is translated downwards and forward relative to the cranial base by growing larger in an upward and backward direction. According to Enlow’s counterpart principle of craniofacial growth, the modeling process at the mandibular ramus results in the lengthening of the mandibular corpus and simultaneous anterior displacement of the mandible, both of which help to keep pace with the amount of anterior maxillary displacement due to growth. Simultaneously, since the expansion of the middle cranial fossa has unequally displaced forward the nasomaxillary complex and the mandible, different amounts of bone deposition and resorption at the posterior and anterior borders of the ramus lead to increased ramus width, which in turn displaces the mandible forward in a coordinated position with the maxilla. Traditional cephalometry is not the most appropriate method for quantifying and/or evaluating the shape of the craniofacial structures since the collection of landmarks, angles, ratios, and linear distances only provide partial and localized information based on the assumption that reference planes are stable. Instead, shape can be investigated with geometric morphometrics by statistically analyzing the coordinates of both traditional landmarks and semi-landmarks, the latter of which can enable the quantification of homologous curves instead of just points. Geometric morphometrics is applicable in orthodontics and can reveal substantial shape information by circumventing the important limiting conventions of traditional cephalometry. The ramus might play a vital, primary role in the development of the mandible and its anteroposterior relationship with the maxilla. Although orthodontic treatment plans often aim to modify the growth of the mandible, orthodontists may underestimate the ramus as an important mandibular module, either by downgrading its role compared to condylar growth or because of limitations of conventional cephalometry. If ramal morphology covaries with the anteroposterior position of the mandible, it could become an essential diagnostic element leading to better-informed treatment plans, which directly benefit patients. Therefore, this study aimed to investigate the covariation between the shape of the ramus and the rest of the face using geometric morphometrics in an adult population. According to Enlow’s counterpart theory as well as the growth mechanism of the mandible at its ramus, wider rami might be correlated with more prognathic mandibles.

The study sample was obtained from the Clinic of Orthodontics and Pediatric Dentistry archives at the Center of Dental Medicine, University of Zurich, Switzerland, between 2003 and 2016. Patients signed an informed consent form before treatment initiation, and the research project was approved by the Cantonal Ethics Committee Zurich (BASEC-Nr. 2020-03012). Since morphological covariation between the ramus and the rest of the face was to be investigated with the RV coefficient, the sample size estimation was based on the findings by Fruciano et al., who reported relatively stable RV values with sample sizes greater than 100. Consequently, at least 100 individuals had to be included in this study. The sample was collected based on the following eligibility criteria: (1) adult patients aged > 18 years, (2) cephalometric radiographs of good quality in habitual occlusion with no functional shifts, (3) no previous orthodontic treatment, and (4) no craniofacial malformations, systemic diseases, syndromes, or myofunctional disorders. Lateral cephalograms were scanned at a resolution of 150 dpi. They were then imported to Viewbox version 4.1.0.10 (dHAL software, Kifissia, Greece) for landmark acquisition. In total, 66 landmarks were digitized: 15 landmarks of types I and II and 51 semi-landmarks, which were organized into five curves and two datasets. The curves corresponded to the outlines of the posterior and anterior borders of the ramus, mandible, maxilla, and frontonasal area. These curves were organized into three datasets and depicted the morphology of the ramus (two curves: posterior and anterior borders of ramus), rest of the face (three curves: mandible, maxilla, and frontonasal area), and whole configuration (all curves; Figure 1). After landmark digitization, semi-landmarks were automatically placed along their adjusted curves in equidistant positions. Semi-landmarks were allowed to slide along their respective curves to minimize bending energy, with the average configuration as the reference shape. Each time, a new average configuration was derived from sliding and set as the new reference shape for subsequent sliding. This process was repeated three times. Intra- and inter-rater measurement errors were evaluated with a Procrustes ANOVA on a set of 30 lateral cephalograms, each against the biological variation of the whole configuration both in shape (Procrustes coordinates) and size (centroid size). Intra-rater assessment was performed at two-week intervals, and inter-rater assessment was performed by two raters. Outliers were visually checked in all three datasets using boxplots for size and Procrustes or squared Mahalanobis distances for shape. The subsequent analyses and graphs were produced in MorphoJ.

Information unrelated to shape was removed by a generalized partial least squares Procrustes superimposition (GLS) with orthogonal projection on the tangent plane. GLS was applied to the raw coordinates of the whole configuration of all individuals in the sample to acquire the Procrustes coordinates and the average shape of the sample in the shape space. A principal component analysis (PCA) was then performed to investigate patterns of shape variation in the study sample. Principal components (PCs) that explained most of the shape variance and were biologically meaningful were considered.

Sexual dimorphism regarding shape was separately assessed in all three datasets using Procrustes distance and permutation tests on the mean values between males and females. Signs of separation between sexes were additionally sought in the shape space by visually inspecting the PC plots of the whole configuration. Sexual dimorphism in size was tested by comparing the logarithm of centroid size (InCS) between males and females for all three datasets with an independent t-test. Furthermore, since the male-to-female ratio of the sample was approximately 1:2, static allometry was expected because males are, on average, larger than females. We assessed its significance with a multivariate regression analysis of shape (dependent variables: Procrustes coordinates) on size (independent variable: InCS) for all the datasets. To further control for allometry, new shape variables corrected for the effect of size were obtained for each dataset by exporting the residuals of the aforementioned respective multivariate regressions.

Two-block partial least squares (PLS) analysis was performed to investigate any covariation patterns between the morphology of the ramus and the rest of the face. The dataset whole configuration was divided into two blocks: ramus and the rest of the face. Two separate Procrustes fits were applied to evaluate joint changes in shape within the ramus and the rest of the face and, consequently, investigate the covariation between the ramus and the rest of the face. The strength of the covariation between the datasets was quantified using the RV coefficient, which can take values between zero and one. PLS analyses were carried out with (Procrustes coordinates) and without (regression residuals) the allometric effect. Furthermore, to better control for the potential influence on the outcome due to the unbalanced male-to-female ratio, PLS analyses were repeated with only females (104 females, 35.5 (standard deviation [SD]) = 11.7 years old).

The sample consisted of 159 individuals (55 males and 104 females) with a mean age of 34.4 years (SD = 11.4). Males were 32.4 (SD = 10.8) and females 35.5 (SD = 11.7) years, respectively. The Mann–Whitney U test did not identify any statistically significant age differences between males and females (U = 2,404, p = 0.1). The baseline cephalometric data for the entire sample are provided in Supplementary Table 1. No outliers were detected in terms of size and shape as far as all the datasets were concerned. The intra- and inter-rater measurement errors explained less than 1% of the total shape variance in the whole configuration for both shape and size.

The first four PCs accounted for 66.1% of the total shape variance of the whole configuration (PC1 24.1%, PC2 21.6%, PC3 11.5%, and PC4 8.9%). PC1 depicted shape differences corresponding to skeletal hyperdivergence and decreased facial depth, whereas a hypodivergent skeletal pattern with increased facial depth was evident in the opposite direction. PC2 described the shape differences in the sagittal plane between Class II and III skeletal patterns (Figure 2). Alveolar shape variability towards bimaxillary protrusion with a steeper mandibular plane or bimaxillary retrusion with a flatter mandibular plane was seen in PC3, whereas PC4 showed variability in the width of the ramus; wider rami were combined with bimaxillary protrusion and a flatter mandibular plane, whereas narrower rami were combined with bimaxillary retrusion and a steeper mandibular plane (Supplementary Figure 1).

Shape differences between males and females were detected in all datasets (9,999 permutations, whole configuration p = 0.0001, ramus p = 0.0102, rest of the face p = 0.0001; Table 1). Although males and females largely overlapped in the PC plots, some separation hints along PC1, PC2, and PC3 were evident (Figure 2 and Supplementary Figure 2). To visualize sexual dimorphism in shape, we performed additional PCA separately for males and females. Most of the variance was explained by shape differences in the sagittal plane in males (30.7%), whereas females mostly presented shape variation in the vertical plane (25.4%). However, PC1 and PC2 described the same shape variation in reversed order for males and females and jointly accounted for 54.4% and 42.8% of the total variance in males and females, respectively. PC3 explained only approximately 10% of the variance in both sexes, which showed changes in the gonial area for males and alveolar changes in females (Supplementary Figure 2). For all datasets, significant size differences between the sexes were detected (p < 0.0001; Table 2). However, although regressions were highly statistically significant for all datasets (p ≤ 0.0028; Table 3), shape variance explained by allometry only ranged from 2.2% (ramus) to 6.0% (the rest of the face).

PLS analyses of all individuals with and without the allometric effect as well as only females produced similar results. Among the groups, the total covariance explained by PLS 1 ranged from 41.9% in the female group to 45.5% in all individuals with the allometric effect included. Results for PLS 2 only ranged from 25.1% in all individuals with the allometric effect removed to 25.5% in the female group and in the group of all individuals with the allometric effect included (Table 4). By inspecting the associated plots and outline drawings, slight shape changes between females and all individuals, with and without the allometric effect, were observed for PLS 1 (Figure 3 and Supplementary Figure 3). These slight shape changes between females and the other groups mostly referred to the sagittal plane, especially the anteroposterior position of the maxilla. In all individuals with the allometric effect included (Figure 3), 45.5% (p < 0.0001) and 25.5% (p < 0.0001) of the total covariance between the ramus and the rest of the face was explained by PLS 1 and PLS 2, respectively. PLS 1 described covariance mainly between the width and height of the ramus, sagittal and vertical positions of the jaws, and thickness of the symphysis. In PLS 1, a wider and shorter ramus with an increased posterior curvature in the lower area (gonial area) was related to a decreased lower anterior facial height, a slightly prognathic mandible with a thicker symphysis, and a prognathic maxilla. PLS 2 depicted covariance pertaining to the position of the lower area of the ramus and changes in the sagittal relationship of the jaws and the inclination of the mandibular plane. A more posteriorly positioned ramus with an increased posterior curvature in the lower area (gonial area) was related to the Class II pattern and flatter mandibular plane. In females, PLS 1 explained 41.9% (p < 0.0001) of the covariance between the ramus and the rest of the face (Table 4). It depicted the covariance between the width of the ramus and the vertical dimension of the face, as well as the sagittal position of the mandible and the thickness of the symphysis. Wider rami in the upper posterior and lower anterior rims were related to a decreased lower anterior facial height and a prognathic mandible with an increased thickness of symphysis (Supplementary Figure 3).

The present study included an extensive range of skeletal patterns in adults. These morphological variations were revealed by PCA analyses, where PC1 described morphological variation in the vertical plane as the divergence of the maxillary and mandibular planes, while variation in the sagittal plane corresponding to Class II or III skeletal patterns was described by PC2. These results are in accordance with those of other orthodontic populations and, collectively, may suggest that facial variation originating from the vertical plane is predominant. Indeed, vertical variation, which is age- and size-related, may be the most significant source of facial variation in humans. Specifically, the overall facial pattern is established early, and skull proportions are determined by growth in the vertical plane, which involves posterior rotation of the mandible as far as anterior facial height variation in humans is concerned. However, if conventional superimpositions on cranial base structures are used, growth rotations are generally masked because of modeling within the jaws, which might impede diagnostic accuracy with conventional cephalometrics in the vertical plane. In this study, the adult orthodontic population was investigated using geometric morphometrics. Such growth rotations are expected to have been completed, along with growth in other planes in adults, which could minimize confounding arising from growth in different structures and developmental times. Furthermore, geometric morphometrics based on sliding semi-landmarks enabled the investigation of curves, which would have been impossible with conventional cephalometric measurements, irrespective of the reference planes. It is important to note that another common area of facial morphological variation is the gonial angle, which, in orthodontic populations, consistently contributes 8–10% to facial skeletal variation. This was also confirmed in the present investigation, since PC3 and PC4, except for alveolar variation, also depicted variations in the gonial area (PC3) and ramus width (PC4), amounting to 11.5% and 8.9%, respectively (Supplementary Figure 1). Ramus width is highly variable among humans, and variation is primarily located at the posterior border. According to Enlow’s counterpart theory, the mandibular ramus is the specific structural counterpart of the middle cranial fossa and pharynx; its horizontal growth can bridge the pharyngeal space and span the middle cranial fossa, resulting in proper intermaxillary fit. However, the lower correlation between the midline cranial base and mandibular ramus, compared to the correlation between the lateral basicranial elements and ramus, indicates that the lateral basicranial structures are more relevant than the midline cranial base. Indeed, in a morphometric study with a non-orthodontic adult population, the correlation between the petrosal part of the middle cranial fossa and the posterior border of the ramus was significantly higher than that between the midline cranial base and the posterior border of the ramus or the basicranium and the anterior border of the ramus. However, common developmental timing appears to strongly affect such integrations since the higher integration between the midline cranial base and the face in children shifts in adulthood; adults demonstrate higher integration between the lateral parts of the basicranium and the face. Since horizontal growth of the ramus might establish the anteroposterior position of the mandible and, consequently, its relationship with the rest of the face, this study investigated the integration of the ramus with the rest of the face. The curves used to capture the ramal shape were drawn along the posterior and anterior ramal margins, excluding the sigmoid notch region, which is biologically irrelevant to the anteroposterior position of the mandible (Figure 1). The shape patterns corresponding to the most covariance-explained (45.5%, p < 0.0001) pertained to the width and height of the ramus in relation to the anteroposterior position of the jaws, the thickness of the symphysis and divergence of the face. Specifically, the wider and shorter the ramus, the more prognathic the jaws were, and the thicker the symphysis and the more hypodivergent the faces were (Figure 3). These findings follow the well-discernible facial patterns in the vertical plane. Compared to hyperdivergent individuals, rami have been reported to be wider in hypodivergent individuals. Furthermore, compared with facial hyperdivergence, more elongated and anteriorly positioned maxilla and mandible, including increased symphyseal dimensions, appear in facial hypodivergence. Finally, the present results corroborate other morphometric findings stating that a narrower and more elongated ramus correlates with facial hyperdivergence, whereas a wider ramus correlates with a hypodivergent pattern. PLS 2, which explained 25.5% (p < 0.0001) of the covariance between the face and ramus, referred to the position of both posterior and anterior rims of the ramus at its caudal region, the sagittal jaw relationship, and the steepness of the mandibular plane (Figure 3). The more posteriorly the ramus was positioned in its lower region, the more Class II the jaw relationship and the flatter the mandibular plane. The posterior position of the ramus and the closed gonial angle leading to a flat mandibular plane may have led to a mandibular retrusive effect. Regarding the differences in the covariance patterns between the ramus and the rest of the face in only females, PLS 1 (41.9%, p < 0.0001) mainly described changes in the anteroposterior dimensions of the ramus in relation to facial divergence and the mandibular anteroposterior position. The wider the ramus, the more hypodivergent the face and the more prognathic the mandible. The covariance patterns described in this study were similar between males and females. However, sexual dimorphism of the mandible has been reported for both size and shape, and the ramus is a highly dimorphic area.

From our findings, an increased ramal width correlates with a decreased lower anterior facial height and prognathic mandible and maxilla. In contrast, a posterior orientation of the ramus in its lower region correlates with a retrognathic mandible and a flat mandibular plane within a Class II skeletal pattern. The width, height, and inclination of the ramus are related to facial changes in both the sagittal and vertical planes.