Developmental Changes in Morphology of the Middle and Posterior External Cranial Base in Modern Homo sapiens

The basicranium has been described as phylogenetically informative, developmentally stable, and minimally affected by external factors and consequently plays an important role in cranial size and shape in subadult humans. Here basicranial variation of subadults from several modern human populations was investigated and the impact of genetic relatedness on basicranial morphological similarities was investigated. Three-dimensional landmark data were digitized from subadult basicrania from seven populations. Published molecular data on short tandem repeats were statistically compared to morphological data from three ontogenetic stages. Basicranial and temporal bone morphology both reflect genetic distances in childhood and adolescence (5–18 years), but not in infancy (<5 years). The occipital bone reflects genetic distances only in adolescence (13–18 years). The sphenoid bone does not reflect genetic distances at any ontogenetic stage but was the most diagnostic region evaluated, resulting in high rates of correct classification among populations. These results suggest that the ontogenetic processes driving basicranial development are complex and cannot be succinctly summarized across populations or basicranial regions. However, the fact that certain regions reflect genetic distances suggests that the morphology of these regions may be useful in reconstructing population history in specimens for which direct DNA evidence is unavailable, such as archaeological sites.


Introduction
Cranial morphology is frequently studied with the purpose of identifying and interpreting the extensive range of variation that exists among modern human populations and in the hominin fossil record (e.g., [1][2][3][4][5][6][7][8][9][10]). The cranium is a valuable structure for studying the genetic and ontogenetic relationships among Homo sapiens based on geographic provenance, population affinities, and dietary proclivities [1][2][3][4][5][6][7]. The vast majority of studies on the relationship between cranial morphology and genetic relatedness, however, have focused on adult crania. To date, little data exist on the ontogenetic trajectories of the developing cranium, especially the basicranium, as they relate to molecular distances among populations. Determining how basicranial shape develops in infants, juveniles, and young adults has implications for bioarchaeology and paleoanthropology, in which skeletal specimens are often discovered that contain no DNA and thus no direct genetic evidence of ancestry. In particular, there would be great utility for future bioarchaeological studies to determine whether the basicranium can be used to classify isolated subadult specimens in populations and draw conclusions regarding the ancestry or population history from basicranial shape of a subadult specimen or group of specimens.
Previous studies have demonstrated that the basicranium is a phylogenetically informative region among human and nonhuman primate adults [11][12][13][14][15]. Human populations differ significantly in their adult basicranial morphology, and these differences reflect the genetic relatedness among populations [1,3,4,6]. The basicranium begins to ossify early in the prenatal period, at approximately 11 to 12 weeks [16][17][18], experiences minimal strain [19][20][21][22][23], and is less susceptible to external factors, such as environment and diet, than other 2 BioMed Research International regions of the cranium [2,3,15]. It is also developmentally stable [11,13,14,24,25], because it emerges from a cartilaginous template early in fetal life, making this region less susceptible to nongenetic forces during ontogeny than the externally sensitive intramembranous bones of the facial skeleton [19,20]. Additionally, the osseous morphology of the basicranium experiences a lower degree of masticatory stresses than other cranial modules, such as the splanchnocranium (face). The higher strain regions are consequently more subject to biomechanical stress and therefore exhibit higher levels of variability [26]. As such, it has been argued that low strain cranial regions, such as the basicranium, should be more stable and more phylogenetically informative than higher strain regions, such as the face [26]. However, it should be noted that this prediction has not held up to empirical testing in several hominoid (ape) species [5,[26][27][28].
Nonetheless, since the cranium starts to develop before birth [16][17][18] and continues to develop until adulthood, it is fundamental to understand the ontogenetic processes, including a comparison of these processes in diverse groups of humans. This can lead to a thorough understanding of how the basicranium develops and whether it reflects among population genetic distances at various stages of subadult ontogeny.
Previous studies have compared the shape of the developing facial skeleton among groups from different geographic locations and genetic backgrounds to determine whether similar ontogenetic processes characterize divergent groups [7][8][9][10]. A wide range of variation in human craniofacial form exists and can be relatively easily altered through minor shifts in the ontogenetic process [9]. There have been a limited number of previous studies comparing cranial ontogenetic patterns among human populations, and most have been limited to the temporal bone [7,[29][30][31]. Thus, the morphological patterning of human cranial ontogeny is still relatively poorly understood.
To date, there have been extensive studies on the morphology of the temporal bone and its applications to phylogeny [32][33][34], evolution, and ontogeny [7]. The temporal bone demonstrates genetic and geographic patterning that is consistent with a predominantly neutral evolutionary history, shaped primarily by mutation, genetic drift, and gene flow [1][2][3][4][5][6][7]34]. However, despite the abundance of studies on temporal bone shape, there is comparatively little information on the ontogenetic trajectories that result in morphological variation at various subadult stages of development.
Modern human populations can be readily distinguished from one another based on one particular component of the basicranium and temporal bone shape [2][3][4]34]. More recently, Smith and colleagues [7] compared subadult and adult temporal bone morphologies in modern Homo sapiens and found that significant differences among populations originate early in ontogeny. Specifically, individuals from different populations can be differentiated based on temporal bone shape prior to the eruption of the first molar [7]. Subsequent developmental stages progress in a relatively similar manner with essentially parallel ontogenetic trajectories, however, maintaining those original differences into adulthood [7]. These findings are limited to the temporal bone but encourage further analysis of other cranial regions, including other larger regions, such as the basicranium, in relation to genetics, geography, and population differences. Several researchers have attempted to determine whether the basicranium is a more reliable region for reconstructing genetic distances than other regions of the cranium [2,3,32]. Smith [3] found that the basicranium was significantly more highly correlated with molecular distances than the cranial vault. Morphology of the temporal bone and the upper face were also correlated with a molecular distance matrix, demonstrating that these bones are phylogenetically informative. von Cramon-Taubadel concluded, however, that the basicranium is not significantly more congruent with genetic data than the intramembraneously ossifying cranial modules (i.e., cranial vault and face) [5,36]. In sum, the basicranium "is a good place to look for reliable characters. . .that describe developmental processes or events" [13, page 159], but its phylogenetic utility during various stage of development is minimally understood.
In light of the previous work and unresolved research questions outlined above, the aims of this study are to determine at what point during the ontogenetic process population-specific basicranial morphologies emerge and to compare basicranial ontogeny among human populations. Specifically, we test the following hypotheses: (H1) Human populations differ significantly in the shape of the basicranium and its various components irrespective of ontogenetic stage.
(H2) Differences among human populations in basicranial morphology are significantly correlated with their genetic distances throughout ontogeny.

Data Collection.
Three-dimensional data on basicranial morphology were collected by one of us (DHD). The data were collected from skulls of seven modern human populations at various ontogenetic stages (Table 1) from the American Museum of Natural History, New York, NY, and The National Museum of Natural History at the Smithsonian Institute, Washington, DC. These populations were chosen based on a single primary criterion: having sufficient cranial ontogenetic series available in museum collections. This criterion is a limiting factor, because museum collections typically consist of primarily adult skulls, since fewer individuals die as subadults and end up housed in collections. Even for individuals that die young, the smaller size and more delicate nature of their crania lead to taphonomic processes being more destructive to subadult skulls, and even those that end up in collections are more likely to be too damaged and fragmentary to measure. Thus, the populations chosen for inclusion in this study were primarily determined by the availability of well-preserved subadult cranial material. Additionally, for the sake of comparability, these populations were identical to the seven populations included in a previous study that we conducted on human temporal bone ontogeny [7]. Prior to traveling to the museums, an intraobserver error study was conducted to ensure the accuracy of the data collection. Two adult crania were digitized ten times each. To minimize the effect of investigator fatigue, the data were collected on two separate days. A paired samples t-test was performed, comparing the Procrustes residuals from the different trials on each skull. The paired t-test did not yield significant differences ( value = 0.96 to 0.99), indicating that the individual landmarks were consistently measured without significant error.
Forty-four landmarks were collected from the basicranium and subcategorized into temporal, occipital, and sphenoidal regions (Tables 2(a)-2(c); Figure 1). These cranial regions have been historically underrepresented in traditional craniometric studies (e.g., [37]). Thus, in order to sufficiently capture the morphology of these regions, it was necessary to include a few empirical, instrumentally determined landmarks in addition to the traditional craniometric points (Tables 2(a)-2(c)). Due to the fact that subadult specimens are more vulnerable to taphonomic processes and often fragmentary, not all subadult specimens here were complete for all cranial regions evaluated, and thus not all specimens were included in all analyses. These landmark data were collected using a MicroScribe G2 digitizer (Immersion Corp., San Jose). Following Smith et al. [7], subadult specimens were assigned a developmental age estimate based on established dental eruption standards [37]: Age category 1 (AC1) = M1s not yet erupted; Age category 2 (AC2) = M1s erupted but not M2s; Age category 3 (AC3) = M2s erupted but not M3s [33]. In humans, these stages correspond roughly to chronological ages of <5 years, 5-12 years, and 13-18 years [37][38][39]. A sample of adult specimens from each population was also digitized; however, since it is not possible to accurately estimate ages for adult crania, these specimens were included in only a subset of the analyses (explained in further detail below).
Data on individual genotypes for 783 short tandem repeats (STRs) were compiled from Ramachandran et al. [40,41] and Rosenberg et al. [35,42] for molecular population representatives matching as closely as possible the populations from which the morphological data were collected. This practice of matching morphological and molecular populations has been successfully employed by previous studies [1-7, 34, 36]. Some of our morphological populations (i.e., the Alaskan and Egyptian samples) were not included in the Rosenberg studies from which the molecular data were derived, and thus, STR data were not available for these populations. In these cases, we attempted to choose a molecular population representative that was as similar as possible to the morphological sample. For the Alaskan sample, we chose the Yakut (a Siberian native population), and for the Egyptian sample we chose the Mozabite (an Algerian Berber group). We recognize that these imperfectly matched molecular representatives introduce a certain degree of incompatibility between the two data types; however, this mismatch renders any correlations that we obtain between the morphological and molecular distances to be a conservative, minimum estimate of the relationships between these types of data.
STRs are composed of back-to-back repeating segments of two to six nucleotides and are found at many locations within the genome [38,39]. They are often variable in their number of repeats, and the lengths of the STRs differ among different populations, in patterns that reflect their ancestry and relatedness. STRs are appropriate for assessing human population distances because they evolve primarily neutrally and are homologous [35,[40][41][42][43]. They are also well-typed for a large number of human populations from geographically diverse locations, cultural backgrounds, and linguistic traditions [36,37].

Analytical Methods.
The morphological data were analyzed using a Generalized Procrustes Analysis (GPA), in which the digitized points were rotated and translated and specimens were scaled to the same size, such that the only remaining differences among them were directly attributable to shape. The GPA was followed by a Principal Components Analysis (PCA) using MorphoJ [44]. Adult specimens were  Most posterolateral point on the occipital condyle 8 Most posterior point on the occipital condyle along the margin of the foramen magnum 9 Mid-point of the occipital condyle (inferior aspect) 10 Opisthion (posterior most point on the foramen magnum) 11 Mid-point on the median nuchal line between the external occipital protuberance and foramen magnum 12 Anteromedial point on the hypoglossal canal 13 Asterion (temporal, occipital, and parietal meet) * * Landmarks are repeated with regard to the overlapping cranial regions.
Asterion (where temporal, occipital, and parietal meet) * 14 Parietal notch (not depicted) 15 Mastoidale (center of the inferior point on the mastoid process) 16 Most lateral point on the margin of the stylomastoid foramen 17 Most lateral point on the vagina of the styloid process (whether process is present or absent) 18 Most posterolateral point on the jugular fossa 19 Most posterolateral point on the margin of the carotid canal entrance 20 Point on anterior margin of tympanic element that is closest to carotid canal 21 Most posterolateral point on the external acoustic meatus 22 Most inferior point on the external acoustic meatus 23 Point on lateral margin of zygomatic process of the temporal bone at the position of the postglenoid process 24 Point of inflection where the braincase curves laterally into the supraglenoid gutter, in coronal plane of mandibular fossa (not depicted) 25 Point on the anterior of the lateral margin of the articular surface of the articular eminence 26 Most inferior point on the postglenoid process 27 Deepest point within the mandibular fossa (instrumentally determined) 28 Mid-point of the articular eminence 29 Most anterior point on the articular surface of the articular eminence 30 Auriculare (most lateral point on the temporal) 31 Suture between temporal and zygomatic bones on inferior aspect of zygomatic process 32 Most inferior point at the sphenotemporal suture closer to the midline (on the sphenoid if disconnected) 33 Most lateral point on the greater wing of the sphenoid (intersection between sphenoid, temporal, and parietal bone) * (not depicted) 34 Most frontolateral point on the greater wing of the sphenoid (intersection between sphenoid, temporal, and frontal bone) * (not depicted) 36 Most posterior, inferior point on the sphenotemporal suture * 37 Apex of the petrous part of the temporal bone * Landmarks are repeated with regard to the overlapping cranial regions.

Sphenoid 35
Most anterior inferior point on the sphenozygomatic suture (sphenozygomatic) (not depicted) 36 Most posterior, inferior point on the sphenotemporal suture * 38 Most lateral point of the foramen spinosum 39 Most lateral point on the margin of foramen ovale 40 Most anterolateral point of the lateral pterygoid plate 41 Most inferior part of the pterygoid hamulus (not depicted) 42 The most anteromedial point of the sphenoidal region on the sphenovomer suture 43 Most posterior point where the vomer meets the medial pterygoid plate 44 Point on the sphenoid in the midline in contact with the vomer (vomer notch) 1 Most anterior point on the basioccipital in the midline (sphenobasion, on the occipital if not connected) * 2 Most lateral point on the basioccipital * * Landmarks are repeated with regard to the overlapping cranial regions. Each set of resulting pairwise population distances was entered into distance matrices for each cranial dataset. The significance of these population differences was assessed using a permutation test of 1000 replicates for each region. In order to assess the degree to which basicranial morphology can be utilized to correctly classify individuals of various ages into the population from which they derived, a discriminant function analysis (DFA) was conducted using the Principal Component (PC) scores from the PCA. The DFA was used to determine whether groups could be classified reliably or if there was excessive morphological overlap. This analysis was conducted for the morphology of the basicranium and then for each of its major components, the temporal, occipital, and sphenoidal regions. These tests were conducted with cross-validation using SPSS version 11.0.1 (SPSS, Chicago, IL). This analysis indicated how well basicranial shape discriminates among populations at the three subadult stages of ontogeny (AC1, AC2, AC3).
A matrix of Slatkin's molecular distances among the molecular representative of these populations was calculated from the published molecular data using Arlequin 3.11 [45] ( Table 3). The matrix of pairwise population Procrustes distances for each cranial region was then statistically compared to the molecular distance matrix using Mantel tests [46] to statistically assess and quantify the correlation among matrices. This analysis indicated whether the morphology of the basicranium, and the regions contained therein, reflected genetic distances at each developmental stage.
A series of descriptive analyses were conducted to evaluate the ontogenetic trajectory in the sample. Principal Component scores were regressed against centroid size and biological age (as determined using dental eruption of each specimen, following Ubelaker [37]) to determine the PCs that were significantly correlated with size or age and therefore indicative of age-related changes in basicranial shape. The regression of centroid size and age against PC scores reveals how shape changes with age (ontogeny) and size (allometry) [8]. Finally, a wireframe was constructed by connecting points in shape-space and morphing them along each major PC axis in Morphologika 2 [47], for the PCs that were found to explain >5% of the variance. This process allowed us to visualize shape changes along these PC axes in the different populations and age categories and describe how basicranial shape changes across all populations and age categories.

Differentiation and Classification of Populations.
Procrustes distances based on basicranial morphology of populations sampled in this study were found to be statistically significant for all combined-age subadult samples (Table 4). When separated into age categories, the majority of the Procrustes distances were significantly different among populations. There were no consistent patterns of population differences among the age categories or regions of the basicranium. The Procrustes distances showed a higher number of pairwise significant differences among populations as the age of the individuals and sample sizes increased. Hence, at AC1, given the smaller sample size for some regions (especially the unfused occipital bone), most populations were not significantly different from each other (Table 5). In AC2, more significant differences were found (Table 6), increasing in AC3 (Table 7). Therefore, (H1)-that human populations differ in basicranial shape at every ontogenetic stage-was not supported. DFAs were performed on all regions of the basicranium: temporal, occipital, and sphenoidal. The DFA for the entire basicranium resulted in cross-validated classification rates ranging from 13.3-34.8% (Table 8). The individual cranial regions yielded higher average correct classification rates: temporal = 41.3%, occipital = 42.6%, and sphenoid = 50.0% (Tables 6-8).
The Egyptian population was one of the most correctly classified in the basicranial (33.3%), temporal bone (55%), and occipital bone (35.7%) data sets (Tables 8-11). The Egyptian population was the most correctly classified population in all regions of the basicranium except the sphenoid bone (Tables 8-11). As predicted, the three closely related Native American populations (Mexican, Utahan, and Peruvian) were frequently classified as each other. Consequently, the least correctly classified group overall was the Utah population.

Correspondence between Genetic
Distances and Basicranial Morphology. The morphological Procrustes distances based on the different age categories and regions of the basicranium were statistically compared to the molecular ST distances based on STR data using Mantel tests. In the combined subadult sample, the morphology of the basicranium and occipital bone was significantly correlated with molecular distances (Table 12). Following the analysis with all populations, populations were systemically removed from the analyses to determine whether the results were unduly influenced by the inclusion of any particular population. When the Egyptian population was excluded from the analyses, the correlations between morphology and molecular distance increased to the point of significance in several additional subgroups (Table 13). Results from both analyses are reported in further detail below.
When the subadult samples for each cranial region were divided into separate age categories, more distinct patterns emerged (Tables 12 and 13). Morphology in AC3 individuals was significantly correlated with genetic distances for the basicranium, temporal, and occipital regions (Table 13). In AC2 specimens, morphology of the basicranium and temporal bone was correlated with the genetic matrix (Table 13). There were no significant correlations between the morphology of any cranial region and genetic distances for AC1 (Table 13); however, the sample sizes of these datasets were much smaller than the other age categories, rendering these results more tenuous.
Overall, the above patterns show that the shape of the basicranium, temporal bone, and occipital bone, for all populations excluding the Egyptian population, reflects genetic distances for the combined subadult sample and AC3 (13-18 years of age), thus supporting (H2) (basicranial differences are correlated with genetic distances) for these data sets. The temporal bone and basicranium also reflect genetic distances during AC2 (5-12 years of age) when the Egyptian population is excluded. The morphology of the sphenoid bone does not reflect genetic distances for subadults of any age category; however, its morphology can be used to discriminate among populations with slightly higher rates than the other two cranial regions. Therefore, (H2) is not supported for the sphenoid, AC1 (for any morphological region), or the occipital bone for AC2.

Age. The regression of biological age with Principal
Components scores (Tables 14 and 15) for the entire basicranium revealed significant correlations between age and PC1 (14.83% of the variance), PC2 (9.31% variance), and PC4 (6.53% variance) (Figures 2 and 3). While these PCs each describe a relatively low amount of variance, this pattern is typical in PCAs following a Procrustes superimposition, because the Procrustes analysis essentially removes the effect of size. PC1 was found to be correlated with age for all populations for morphology of both the entire basicranium and the isolated temporal bone (Table 15). Additionally, PC4 was correlated with biological age for sphenoid morphology in all individual populations except the Polynesians, and PC6 was correlated with age for the occipital bone in all individual populations except the Peruvian sample (Table 15).

Centroid
Size. The regression analysis of centroid size with PC scores for the entire basicranium revealed significant correlations between size and PC1 (14.83% of the variance), PC2 (9.31% of the variance), PC3 (7.57% of the variance), and PC4 (6.53% of the variance) (Tables 14 and 16). The regression plot of PC1 scores versus log centroid size is illustrated in            revealed by the age versus PC1 comparisons (Table 15). This reflects the fact that size and age are related, as expected. PC3 was significantly correlated with centroid size for all individual populations in the occipital bone, a pattern not seen in all of the various age categories. Similar to the age comparisons, the centroid size of the sphenoid bone and PC4 were significantly correlated for all populations individually. Mean PC scores for each age group in each population are provided for comparative purposes (Tables 17 and 18).

Shape Space Exploration.
In order to further explore how populations and age categories differed with regard to basicranial morphology, we conducted a shape space exploration in which wireframes of each cranial data set were morphed along the major PC axes to visualize how the shape varied along each PC. In basicranial shape, subadults were commonly located on the +PC1 axis (Figure 3(a)). A positive PC1 score is associated a mediolaterally narrow basicranium, with an elongated external acoustic meatus (EAM). The external occipital protuberance (EOP) also moves relatively further away from the occipital condyle with an increasingly positive PC1 score. Given that most of the subadult values appeared on the +PC1 side, this indicates that the subadult basicranium is longer and narrower but widens into adulthood, while relatively decreasing in length. A negative PC1 score is characterized by the reverse morphological conditions (i.e., mediolaterally wider basicranial, shortened EAM, closer proximity of the EOP, and occipital condyle).
Associated with a positive PC2 score (Figure 3(a)), the EAM becomes relatively larger, but the overall basicranial length decreases. The hypoglossal canal and articular eminence move relatively further posteriorly. The length between the occipital condyle and the EOP decreases as well, thus shortening the basicranium. The EAM length relatively increases.

Discussion
To date, extensive research has been conducted on the temporal bone and how it reflects genetic relationships in nonhuman primates (e.g., [13,14,[31][32][33]48]). In modern humans, it has been determined that the temporal bone reflects genetic relationships in both subadults and adults and that many populations significantly differ in their ontogenetic trajectories [7]. The current study focused on a larger area of the skull, the basicranium and its constituent regions, and their ontogenetic patterns.

Hypotheses.
The hypotheses that formed the basis for the study were found to be only partially supported. In actuality, the patterns of morphological variation turned out to be more complex than the relatively simply stated hypotheses. The patterns of the individual basicranial bones differed from each other and across the various age categories. Thus, the ontogenetic processes driving basicranial development cannot be succinctly summarized across all regions of the basicranium.
(H1) Human populations differ in the shape of the basicranium and its various components irrespective of ontogenetic stage.
This hypothesis was not supported by our results. All populations were found to differ significantly in basicranial morphology in the combined subadult samples. However, in the separate age cateories there was a trend toward increasing differentiation with age. In AC1, several population were not significantly different, by AC2 most populations differed, and by AC3 all populations were significantly different. Interestingly, of all the individual cranial bones, the sphenoid was revealed to be the most distinct among populations and therefore the most reliable region for population classification.
(H2) Differences among human populations in basicranial morphology are significantly correlated with their genetic distances throughout ontogeny.
This hypothesis was partially supported for several subsets of our data. Our results show that morphology of the basicranium, occipital, and temporal bones each significantly Table 15: Principal Components of the basicranium and its regions significantly correlated with age in the regression analysis.

Classification among Populations.
The discriminant function analyses revealed that the Egyptian population was the most correctly classified population for the basicranium, occipital, and temporal bone in subadults. The Egyptian population was most likely classified correctly because it is the most geographically and genetically distant of the populations in this study. This finding is consistent with the previous description of Egyptian cranial morphology as unique and distinctive [49]. Unsurprisingly, the closely related Native American groups of Utah and Mexico were often classified as each other in basicranial, temporal, and occipital bone for all subadults. This trend has been supported by previous studies [6,7] and supports the hypothesis that closely related populations should share a similar basicranial shape. Relethford [6] showed that the Peruvian population was phenotypically very similar to other Native American populations in his study (Arikara, Greenland Inuit and Santa Cruz) despite its remote geographic location compared to the other Native American populations. Consistent with the results of this study, Smith and colleagues [7] also found that the temporal bone of the subjects from Utah and Mexico was most often classified as another Native American population. In fact, the low levels of correct classification for these populations suggest that their morphological variation overlaps considerably in the younger age categories. These findings, taken together, show that the basicranial morphology can generally be used as an indicator of genetic relatedness among children and adolescents but considerable overlap may still exist for populations which are close in genetic proximity, such as Native American populations.
Based on the findings of this research, the basicranium, temporal, and occipital bone reflect genetic distances in childhood and adolescence, but this study suggests that these differences are not seen in infancy. AC1 was not correlated with molecular distances among populations for any of the cranial regions. Two possible scenarios may explain these results in infancy. First, infant basicranial morphology may be similar among populations and the observable differences in later ontogenetic stages have not yet developed. However, an alternative explanation is that as a result of the fact that many basicranial bones are not fully fused in infants, our sample sizes in the AC1 age category may have been insufficient to reveal subtle differences among populations in infancy.
These findings contrast with some of the results found previously for the temporal bone [7,31], which suggest that the relationships in temporal bone originate early in ontogeny and these differences reach adulthood via different ontogenetic trajectories. Terhune et al. [31] also found that subadult great apes and humans demonstrate differences in temporal bone shape early in ontogeny. Contrary to Smith et al. [7] and Terhune et al. [31] for the temporal bone, and Viðarsdóttir et al. [8] for the face, the present study did not find that early subadult morphology (i.e., AC1) reflects adult population-specific differences. This finding is surprising given that the temporal bone has been found previously to reflect genetic distances in adult humans [2][3][4][5][6][7]. However, these seemingly contrasting results may reflect the choice of different landmarks used in the present study. Alternatively, these differences may be explained by the relatively small sample sizes necessitated by the limited subadult cranial material available in museum collections. Interestingly, while morphology of the sphenoid was not found to reflect the molecular distance matrix at any ontogenetic stage evaluated (Table 9), it did perform quite well in the DFAs (Table 8). Thus, the morphology of this bone can be inferred to be distinct among the populations sampled here; however, those differences do not reflect genetic relatedness. In other words, more closely related populations do not share similar sphenoid morphology, but instead, each population is unique in its sphenoid shape.

Cranial Morphology and Genetic Relationships.
Overall patterns revealed by this study show that the shape of the basicranium, the temporal bone, and occipital bone, for all populations excluding the Egyptian population, reflects genetic distances in subadults. Removing the Egyptian population from the analyses yielded significant results in the temporal bone, occipital bone, and age categories AC2 and AC3. It cannot be definitively ascertained why the Egyptian population deviated from the patterns of the rest of the included samples, but one possibility is the relatively high degree of mismatch between the morphological Egyptian sample and its molecular population representative, the Mozabite. The Mozabite people live in Algeria and speak a Berber language. As one of the few northern Saharan populations to have been extensively studied for neutral molecular loci, they are the closest well-typed molecular representative for the Egyptians, but certainly not a perfect match.
Patterns on the basicranial PC wireframe showed that the Egyptian crania appeared on the +PC2 scale, indicating that the basicranium in the Egyptian population is relatively shorter anteroposteriorly compared to the other populations. Most of the Egyptian population clustered towards an increasingly long occipital bone (-PC2, +PC3, and -PC4) and increasing large temporal bone. It appears that the Egyptian cranium might have a relatively smaller overall basicranium but that the temporal and occipital bone increased in length, with the sphenoid decreasing in length to compromise. A study by A. C. Berry and R. J. Berry [49] found that the Egyptian population was characterized by distinctive cranial morphology, which had persisted for centuries and was dissimilar to the crania of other populations in their study. Based on the current findings, the Egyptian population showed a relatively decreasing anteroposterior length but increasingly elongated temporal and occipital bones compared to other population samples in this study. Given its genetic and geographic distance and cranial stability, the Egyptian population could have likely acted as a morphological and genetic outlier, and only with its exclusion did the results begin to show significant differences among populations in various regions of the basicranium and age categories.
One potentially limiting factor of this study is the degree of mismatch between the molecular and morphological population representatives. The morphological samples were chosen based on the availability of subadult cranial material in museum collections. Thus, we were not able to include a wide geographically distribution of populations. In North American natural history collections, there is a natural bias towards Native American specimens. Consequently, the present study included a few Native American samples. Similarly, the molecular representative for each population was not a perfect match with these morphological populations. However, there is precedent for this approach. As noted by Roseman (2004), such mismatch is not necessarily problematic, but the correlations obtained from this type of analysis should be interpreted as minimal approximations of their actual value. Thus, the significant correlations obtained here should be interpreted as minimum estimates of a real biological relationship between morphology of the cranial region and genetic relatedness.
The landmark dataset employed in the present research expands the smaller landmark sets used by previous studies. Harvati and Weaver [2] and Smith et al. [34] each used a much smaller and less anatomically distributed set of landmarks, 13 and 15, respectively, while our study included 24. There were also some differences in landmark coverage and dispersion between this and prior studies. The temporal bone landmarks in this study were widely distributed to include points that overlap with the sphenoidal and occipital regions of the basicranium, covering a large area of all components of the basicranium. We used landmarks distributed from the asterion in the occipital region and auriculare in the temporal  Figure 3: Plot of the first four principal components of the basicranial landmark configuration that were significantly correlated with developmental age and centroid size: (a) PC1 versus PC2 (representing 14.83% and 9.31% of the variance, resp.); (b) PC3 versus PC4 (7.57% and 6.53% variance, resp.). Each quadrant of the plot contains a wireframe indicating the morphology typical for the specimens in that quadrant. region, to the suture between the temporal and zygomatic bones, for a total of 44 landmarks. Harvati and Weaver [2] also used several landmarks on the petrotympanic crest, which this study did not include. Smith et al. [34] included landmarks on the mandibular fossa, porion, and entoglenoid process, which were not used in this study.
Occipital bone morphology is valuable for assessing relatedness in adolescents starting around the age of 13 years old. The sphenoid bone, on the other hand, does not reflect genetic relatedness at any ontogenetic stage; however, despite this, the morphology is distinct enough among populations to allow unidentified specimens to be classified by population