Length data pointed out statistically significant differences among groups (p < 0.01); the Medieval population, usually, showed higher mean values in comparison to other groups. Data referred to heights revealed statistical significance for Nasion-prosthion height (NPH) and Nasal height (NLH) variables, with lower mean values being observed in the Domus population. As regards the breadth measurements, statistically significant differences were observed for the variables Minimum cranial breadth (WCB), Biorbital breadth (EKB), Orbital breadth (OBB), and Nasal breadth (NLB). A progressive increment of WCB was detected through historical ages. The NLB had an unusual figure: Nuragici, Medieval and Modern groups showed similar values and only Domus significantly differed from the others (p < 0.01), showing the lowest values (Table 1).
Table 1
Comparison of craniometric variables with Howell’s Code25 among the Sardinian population: mean and standard deviation (in brackets). Means differences assessed by Anova one-way analysis.
Craniometric Variable | Domus | Nuragici | Medieval | Modern | p-value |
Basion-nasion length (BNL) | 105.2 (5.2) | 106.8 (4.2) | 114.7 (6.5) | 110.4 (6.1) | < 0.01 |
Basion-prosthion length (BPL) | 99.2 (6.5) | 102.2 (6.7) | 109.3 (6.3) | 107.3 (5.8) | < 0.01 |
Maxillary-alveolar length (MAL) | 45.5 (4.7) | 48.1 (4.5) | 53.9 (4.6) | 51.2 (4.0) | < 0.01 |
Nasion-prosthion height (NPH) | 67.2 (3.5) | 66.9 (1.9) | 72.5 (4.5) | 68.7 (5.1) | < 0.01 |
Orbital height (OBH) | 31.8 (1.4) | 31.2 (1.7) | 32.7 (3.1) | 33.4 (2.8) | 0.30 |
Nasal height (NLH) | 51.5 (2.8) | 52.3 (1.0) | 55.1 (3.6) | 52.6 (3.4) | 0.01 |
Maximum cranial breadth (XCB) | 143.0 (7.0) | 140.8 (4.4) | 145.9 (7.9) | 147.6 (5.9) | 0.15 |
Minimum cranial breadth (WCB) | 104.0 (3.3) | 105.2 (1.1) | 111,8 (5.7) | 110.1 (5.9) | < 0.01 |
Maximum frontal breadth (XFB) | 118.9 (5.4) | 115.3 (2.3) | 120.4 (6.8) | 118.2 (7.0) | 0.39 |
Biorbital breadth (EKB) | 102.9 (3.6) | 101.0 (2.1) | 111.9 (5.5) | 106.8 (10.0) | < 0.01 |
Bizygomatic breadth (ZYB) | 136.0 (7.8) | 130.9 (4.9) | 138.2 (6.7) | 139.6 (5.7) | 0.06 |
Bimaxillary breadth (ZMB) | 96.3 (3.1) | 96.0 (1,2) | 98.1 (5.3) | 98.5 (3.4) | 0.46 |
Orbital breadth (OBB) | 40.2 (2.0) | 40.6 (2.8) | 43.3 (2.3) | 41.5 (1.9) | < 0.01 |
Nasal breadth (NLB) | 23.5 (1.9) | 26.1 (0.7) | 26.9 (2.3) | 25.9 (2.3) | < 0.01 |
Palate breadth external (MAB) | 64.7 (3.9) | 63.7 (3.8) | 65.7 (5.4) | 67.0 (6.7) | 0.58 |
Anterior-posterior nasal spine | 48.4 (5.4) | 47.1 (7.4) | 53.9 (5.2) | 57.2 (5.0) | < 0.01 |
Intraorbital breadth | 27.6 (2.9) | 26.2 (3.0) | 29.4 (3.2) | 28.4 (1.9) | 0.08 |
Alveolar spinal height | 17.5 (1.5) | 17.0 (1.8) | 18.5 (3.3) | 19.6 (1.8) | 0.07 |
Minimum alveolar height | 42.4 (2.4) | 42.2 (0.5) | 43.6 (7.4) | 44.0 (3.9) | 0.85 |
Prosthion total angle | 89.6 (1.8) | 91.0 (5.0) | 93.7 (5.8) | 92.0 (4.3) | 0.16 |
Prosthion nasal angle | 91.5 (2.6) | 93.0 (4.8) | 91.6 (7.2) | 93.4 (3.7) | 0.57 |
The Domus sample showed high values of distance to the Medieval population (544.7) and to the Modern one (330.6). The Modern population showed no similarity to prehistoric populations nor to Bolton’s standard. Similar results were obtained using the bootstrap technique on the original data set, without Bolton’s standard, allowing us to build up the means and standard errors after 500 iterations for each variable (Table 2).
Table 2
Squared Euclidean Dissimilarity Coefficient Matrix obtained for the populations
| Medieval | Bolton1 | Modern | Nuragici |
Bolton1 | 603.8 | | | |
Modern | 117.2 | 293.4 | | |
Nuragici | 545.0 | 253.8 | 379.9 | |
Domus | 544.7 | 279.8 | 330.6 | 87.3 |
1 Bolton values were used as a reference |
Similar results of UPGMA were obtained using the bootstrap technique on the original data set, without Bolton’s standard, allowing us to build up the means and standard errors after 500 iterations for each variable (Table 2,3).
Table 3
Squared Euclidean Dissimilarity Coefficient Matrix after 500 bootstrapped iterations below diagonal and relative standard errors, above diagonal, for the populations.
| Medieval | Modern | Nuragici | Domus |
Medieval | - | .0003 | .0011 | .0005 |
Modern | 115.9 | - | .0009 | .0004 |
Nuragici | 544.9 | 379.9 | - | .0005 |
Domus | 544.7 | 330.6 | 87.3 | - |
The similarity among populations is immediately appreciable from the dendrogram (Fig. 1).
A heed emerging from the view tree is the clear differentiation in two clusters separated by Bolton’s standard; Nuragici and Domus at the top and Medieval and Modern at the bottom.
Figure 2 displays the D-plot of the first two principal coordinates derived from Principal Component Analysis (PCA). The first two eigenvalues collectively account for 80.9% of the total variance (59.1% and 21.8%, respectively).
Clearly, Domus and Nuragici populations tend to form separate clusters. The goodness of fit was high (≥0.75) for Bolton, Medieval and Nuragici, and medium-high (0.5–0.75) for the other groups.
The length and the breadth of the neurocranium, Basion-nasion length (BNL), and Basion-prosthion length (BPL), highly contribute to the first principal component. The variance along the second axis was principally generated by Maxillary-alveolar length (MAL), Orbital height (OBH), Maximum cranial breadth (XCB) and Palate breadth external (MAB) positive sense, Maximum frontal breadth (XFB), EKB and OBB, negative sense. BNL, BPL and XFB clearly separate between groups (Domus and Nuragici versus Modern and Medieval), the former with a large frontal, biorbital and orbital breadth, the latter with a longer and narrower skull pattern (BNL and BPL).
The first component places the Domus close to the Nuragici population (negative sense) and separates them from Modern basically due to BNL, MAL, OBH, XCB, WCB, EKB, ZMB, and MAB.