A test of the lateral semicircular canal correlation to head posture, diet and other biological traits in “ungulate” mammals

For over a century, researchers have assumed that the plane of the lateral semicircular canal of the inner ear lies parallel to the horizon when the head is at rest, and used this assumption to reconstruct head posture in extinct species. Although this hypothesis has been repeatedly questioned, it has never been tested on a large sample size and at a broad taxonomic scale in mammals. This study presents a comprehensive test of this hypothesis in over one hundred “ungulate” species. Using CT scanning and manual segmentation, the orientation of the skull was reconstructed as if the lateral semicircular canal of the bony labyrinth was aligned horizontally. This reconstructed cranial orientation was statistically compared to the actual head posture of the corresponding species using a dataset of 10,000 photographs and phylogenetic regression analysis. A statistically significant correlation between the reconstructed cranial orientation and head posture is found, although the plane of the lateral semicircular canal departs significantly from horizontal. We thus caution against the use of the lateral semicircular canal as a proxy to infer precisely the horizontal plane on dry skulls and in extinct species. Diet (browsing or grazing) and head-butting behaviour are significantly correlated to the orientation of the lateral semicircular canal, but not to the actual head posture. Head posture and the orientation of the lateral semicircular canal are both strongly correlated with phylogenetic history.

The need for a reliable and reproducible way of orienting dry skulls for cranial measurements has led to a considerable amount of literature suggesting that the plane of the lateral semicircular canal (LSC) of the bony labyrinth (the osseous capsule of the inner ear) is horizontal when the head is held in its "habitual" (i.e. not actively attained) or "alert" positions [1][2][3][4][5][6][7][8][9][10][11][12] . This is backed by the hypothesis that a horizontal orientation of the LSC would mechanically maximize the recording of rotational and linear head movements made in the horizontal plane by placing the sensory hair cells of the semicircular canal and its associated ampulla perpendicular to the horizontal plane [12][13][14][15][16] . The subsequent use of the orientation of the plane of the LSC as a proxy to infer head posture in fossil vertebrates has grown more popular among paleontologists, as it is being applied to dozens of extinct taxa such as archosaurs, including dinosaurs, and synapsids, including mammals 12,[17][18][19][20][21][22][23][24][25][26][27][28][29] . This has raised discussion on some crucial paleobiological questions, such as the evolution of bipedalism in ancient hominin 14,18 and paleodiets. As browsers are expected to hold their head higher than grazers 30 , head posture has been invoked in reconstructing ancient diet in fossil herbivorous species 20,24 . Semi-aquatic species, on the other hand, would hold their head tilted upward 27 (but see Neenan and Scheyer 31 ). In addition, head posture is directly involved in discussions about the origin of endothermy, as blood pressure to perfuse the head, and particularly the brain, directly depends on head posture and thermophysiology (species with low metabolism have a lower blood pressure than species with a high metabolism, and therefore cannot perfuse their brain if their head is held far above their heart) 32 . Head posture may thus be crucial for inferring the evolution of endothermy in birds, mammals, and their respective

Materials and methods
Sampling. As the inclusion of a statistically-significant number of taxa was essential to this study, we chose to focus primarily on ungulate-grade mammals (i.e. Paenungulata, Artiodactyla, Perissodactyla, and Tubulidentata). Ungulates are more abundant than carnivores or primates in zoos, easily identifiable, and well represented in institutional dry skull collections. They display a wide array of body sizes, a greater variety of documented head postures 47 , a wider range of expected inner ear orientations (as hypothesized from the inclination of the snout compared to that of the brain-case 48 ), and more varied degrees of adaptation to head-butting 34 than any other mammalian group. Moreover, they are the ideal target group to address if diet (browsing v. grazing) plays a significant role in the orientation of the LSC, as previously suggested in the literature 20,24 . Finally, they usually display an elongated snout, which makes it easier to compare the orientation of the head in live animals to that of the corresponding dry skulls.
Head posture in live animals. Head posture was documented by taking pictures of zoo animals in lateral view using a camera equipped with a spirit level ( Fig. 1) to ensure that pictures were taken as close to the horizontal plane as possible. The animals were photographed in 2018 and 2019 at the National Zoological Garden, Pretoria (South Africa), Johannesburg Zoo (South Africa), Montecasino Bird Garden, Fourways (South Africa), Lory Park Animal and Owl Sanctuary, Midrand (South Africa), Ménagerie du Jardin des Plantes, Paris (France), Parc Zoologique de Paris (France), Prague Zoo (Czech Republic), Chester Zoo (United Kingdom), Zoologischer Garten Berlin (Germany), Tierpark Berlin (Germany), and Zooparc of Beauval (France). The saiga antelope pictures were kindly provided by K.H. Vogel. The dataset represents about 10,000 pictures documenting the head posture of 129 species and is available here: https ://osf.io/4vpnj /?view_only=3dc98 7012f cd44a 6a64a d7d89 49ec0 1f (https ://doi.org/10.17605 /OSF.IO/4VPNJ ). The pictures were taken from outside the enclosures to avoid interaction with the animals. It was essential for this study that the animals remain calm and act naturally, so their environment was not disturbed, and the animals were not put on leash or isolated. As such, individual identification was not possible. Representatives of both sexes are mixed in the dataset as sexes could not always be determined. The typical photography set up is illustrated in Fig. 1.
To ensure that the photographed head postures were comparable between individuals and species, all pictures were taken by one of the authors only (J.B., except for the saigas, which were taken by Alexander Sliwa from the Kölner Zoo). The pictures used for this study were selected to reflect as closely as possible what will hereafter be referred to as the "neutral" head posture. The neutral posture of an "ungulate" is here defined as the angle between the main axis of the head and the horizontal when an animal's head remains still, its attention is not attracted by a moving or immobile target, and it is not foraging, drinking, or performing any other identifiable activity involving head movements (e.g. sniffing). The animal can be standing or lying down. "Neutral" head posture differs from head posture "at rest" as it encompasses ruminating animals and individuals slowly walking with their head steady (not pitching up and down while moving). Alert postures 49 were included only if the animal's attention was not directed toward an identifiable direction, and were avoided as much as possible. That is why this study focuses on zoo animals, which are accustomed to human presence. For consistency and to enable comparisons, pictures of semiaquatic "ungulates" (e.g. hippos) were taken when the animal's head was not immersed so that their head posture was not influenced by buoyancy.
The orientation of the head compared to the horizontal plane was measured by J.B. using ImageJ as the angle between the horizontal border of the picture (horizontality of which was ensured by the use of a spirit level on the camera, Fig. 1) and the main axis of the head (traced as the axis running from just above the upper lip to the middle of the occiput on the back of the head) in strict lateral view (Fig. 2a). An average neutral head posture was then calculated for each species (Table 1). The intraspecific standard deviation (measurement error) for neutral head posture is ± 1.6°.
The bony labyrinth is one of the first organs to completely ossify in mammals as its adult size and shape are reached at mid-gestation 50,51 . However, the orientation of the LSC seems to show age-related variations in some tetrapod species, including humans, which may impact their head posture 26,52,53 . As such, juveniles were excluded from the dataset.  Supplementary Table S1 for details), were used.
The bony labyrinths of each skull were segmented manually and reconstructed in 3D using the software AVIZO 9 (FEI VSG, Hillsboro OR, USA) at the virtual imaging labs of the Evolutionary Studies Institute and the Natural History Museum of Basel. The skull was reconstructed using either the Isosurface or threshold functions under the same software. The angle between the plane of the LSC and the main axis of the skull was then measured in lateral view in 2D (Fig. 2b). The plane of the LSC was determined visually in lateral view, following most previous authors 19,20,22,[24][25][26] . The main axis of the skull was traced as the axis running from just above the www.nature.com/scientificreports/ premaxilla (approximately at the level of the centre of the nasal opening) to the middle of the occiput (Fig. 2b) in order to maximize the homology with the measurements taken on living animals. This angle represents the anterior tilting of the head if the LSC is considered horizontal. This angle is hereafter referred to as "the reconstructed cranial orientation" or "reconstructed head posture". Measurements were taken bilaterally when both bony labyrinths were available and then averaged for each species (Table 1). For consistency, all measurements were taken by the same author (J.B.). None of the samples expressed strong lateral tilting of the LSC or an undulating morphology that could impede taking this measurement or affect its accuracy. The intraspecific standard deviation (measurement error) for reconstructed head posture is ± 2.1°. The complete dataset of reconstructed head postures is available in the Supplementary Table S1. This dataset was complemented by measurements made on the published pictures from Girard and Schellhorn 5,30 (see Supplementary Table S1).
As for the picture dataset, only the individuals showing reasonable signs of maturity (e.g. cranial bone fusion, erupted molars) were considered. Data processing. The dataset was analyzed using phylogenetic comparative methods to control for the non-independence of observations [54][55][56] . We used the time-calibrated phylogenetic tree of mammals of Bininda-Emonds et al. 57 because it encompasses all the species in our dataset and fossils can be easily added to it in future analyses. The tree was pruned to match the species in our dataset using function 'drop.tip' in R package ape 58 . All subsequent analyses were performed in R v3.6.3 (R Core Team, 2020). The phylogenetic signal of individual variables was estimated using Pagel's lambda 59 for continuous features (reconstructed cranial orientation, neutral head posture, body mass) using function 'phylosig' in package phytools 60 . Lambda was chosen over the other commonly used estimator K 61 because of the latter's poor performance for trees with small sample sizes and polytomies 62,63 , both of which can be found in our dataset. For binary traits (head-butting; see below), phylogenetic signal was estimated with the D-statistic 64 using function 'phylo.d' in package caper 65 . To test whether the plane of the LSC can be used as a reliable proxy to reconstruct the neutral head posture, we regressed the neutral head posture of living animals on the reconstructed cranial orientation using data in Table 1, and phylogenetic generalized least squares (PGLS) regressions 66,67 . PGLS were compiled using the 'gls' function in package nlme 68 , with correlation structures for each evolutionary model specified in ape 58 . A model selection procedure based on the corrected Akaike information criterion (AICc) was applied to the regressions using the package AICcmodavg 69 . Five evolutionary models were considered for this selection procedure (see 70 ): Brownian Motion, Pagel's Lambda, Ornstein-Uhlenbeck, Early Burst, and White Noise -i.e. non-phylogenetic, ordinary least squares (OLS) regression. All regressions were performed using raw and log-transformed data (natural logarithm). Both variables in the models are in the same unit and order of magnitude, and the models  www.nature.com/scientificreports/ built with raw data showed a higher significance and met parametric assumptions better than models built with log-transformed data. For this reason, we used the former to assess the relationship between the two variables. Because of the high degree of body mass allometry in neuroanatomical features [71][72][73] , body mass measurements for all species in the sample were taken from the literature (Supplementary Table S1) and included as a co-predictor to be tested against models built with only the reconstructed and neutral head postures as predictors in the AICc-based model selection procedures. The coefficient of determination and p-value for generalized least squares regressions cannot be compiled straightforwardly due to the autocorrelated structure of the residuals 67 . Following Paradis 55 , we compiled a pseudo-R-squared and p-value based on McFadden's formula 74 , based on a likelihood ratio test between our model and a null model. Normality and homoscedasticity of the residuals were assessed using a Shapiro-Wilk test and a Q-Q plot, and graphically using residuals v. fit plots, respectively 75 .
Finally, phylogenetic one-way Analyses of Variance (phylANOVA) 76 with False Discovery Rate posthoc corrections 77 were used to test for a difference between groups in three separate factors, for both reconstructed cranial orientation and neutral head posture ( Table 1). The first factor is diet, for which species were categorized as "browser", "grazer", "mixed" (for a mixed diet between browsing and grazing), or "other" (for omnivorous and myrmecophagous species). The second predictor is whether a species practices head-to-head combat. The head-butting category includes wrestling and ramming species (hereafter referred to simply as head-butting species) but excludes flank-butting species (e.g. giraffes). The last predictor is the habitat, which was scored between open (savannah or steppes), closed (forest or jungle), mixed (mix of open and closed habitats), rocky (for species living on steep, rocky slopes), or semi-aquatic. The scoring of all three predictors was done using the literature (see the list in Table 1). PhylANOVAs were performed using function 'phylANOVA' in phytools 60 .

Ethics declarations.
As the animals were not approached or armed, no ethical clearance was necessary for this study.

Results
All the data in the dataset for which a phylogenetic signal could be measured (neutral head posture, reconstructed cranial orientation, body mass, and head-butting) carry a strong phylogenetic signal (lambda > 0.8 for the first three variables; D = − 0.2841056 for head-butting) (See Supplementary Table S1). Species with body mass under 100 kg have an average neutral head posture of 30° and reconstructed cranial orientation of 39°, whereas species larger than 100 kg have a neutral head posture averaging 37° and average cranial orientation of 40°. This suggests an effect of body mass on head posture as was hypothesized by Köhler 47 , but not on the orientation of the LSC (Fig. 3). This is consistent with statistical analyses, which identify a very weak effect of body mass on neutral head posture (R 2 = 0.040; p-value = 0.014), and none on reconstructed cranial orientation (R 2 = 0.023; p-value = 0.054) using OLS. However, once corrected for phylogeny using PGLS, the effect of body mass on head posture (R 2 = 0.030, p-value = 0.06025) and cranial orientation (R 2 = 0.012, p-value = 0.9903) is no longer significant.
Phylogenetic regressions (Fig. 4a) identify a statistically significant (p-value = 4.519e−07) but relatively low correlation (R 2 = 0.261) between neutral head posture and the reconstructed cranial orientation. This supports that the orientation of the LSC in life is correlated to the neutral head posture in "ungulates". The equation of the linear model is: The 95% confidence interval for the slope (0.242-0.526) is significantly different from 1, which means that this model cannot be approximated to an isometric relationship (which would be expected if the LSC was held horizontally).
The model including body mass, neutral, and reconstructed head postures, and the interaction term of the three as co-predictors was selected by AICc as fitting our data best (Fig. 4b). This model shows results very similar to those of the simple regression model, being significant with a slightly stronger correlation (R 2 = 0.325; p-value = 3.235e−07). The equation of the resulting model is written: The slope confidence interval is slightly lower than that of the previous model (for reconstructed cranial orientation: 0.148-0.508), which removes it even further from an isometric relationship (Fig. 4b). The very low   www.nature.com/scientificreports/ coefficients for body mass and reconstructed orientation × body mass are not significantly different from zero (see Supplementary Table S1), and are reported here to ensure full transparency of our results. Surprisingly, for both simple and multiple regression models, the evolutionary model selected by AICc was the Early Burst (EB) model 61 , representing a rapid adaptive radiation followed by stasis. EB models are known to be rarely selected as the best evolutionary model in such selection procedures 78 .
On average, browsers tend to hold their heads less tilted anteriorly (26°) than mixed feeders (32°), and grazers (36°) in neutral posture (Fig. 5). This seems to reflect on the reconstructed cranial orientation as browsers have a higher reconstructed head posture (33°), than mixed feeders (40°) and grazers (44°) (Fig. 5). Phylogenetic ANOVAs indicate a significant difference between browsers and grazers for the reconstructed orientation of the skull (F = 7.723; p-value = 0.046), but not for the neutral head posture (F = 2.663; p-value = 0.516). Mixed feeders are statistically indiscernible from both browsers and grazers in any case.
Variations of reconstructed cranial orientation do not show any trend with habitat preference (Fig. 6), but a slight trend toward more downwardly tilted head postures seems to occur in species living in a more open habitat (Fig. 6); however, this trend is not significant (F = 1.343; p-value = 0.792). Semi-aquatic species seem to have a more posteriorly tilted LSC resulting in a higher reconstructed cranial orientation (23°) than fully terrestrial species, but this does not reflect on the neutral head posture (Fig. 6). Unfortunately, the sample size for this category was too low to effectively test if this difference was significant or simply the result of the scarcity of semi-aquatic species in the dataset. www.nature.com/scientificreports/ Head-butting is found to have a highly significant statistical effect on reconstructed head posture (F = 39.467; p-value = 0.002), with a difference of 13° between head-butting and non-head-butting species on average (Fig. 7); however, with only 1° difference on average (Fig. 7), the same is not true for the neutral head posture (F = 3.126; p-value = 0.591).

Discussion
LSC orientation is correlated with head posture, but is not horizontal. The assumption that the plane of the LSC can be used as a reliable indicator of the horizontal plane on dry skulls and thus can serve as a proxy to reconstruct the "habitual" or "alert" head posture in extinct species has been a long-held 17 , yet insufficiently tested hypothesis in mammals. Attempts to test this hypothesis have highlighted that the plane of the LSC is often tilted upward compared to the horizontal in most mammalian species 5,6,11,30,45 . A famously  www.nature.com/scientificreports/ baffling example is that of humans, in which aligning the plane of the LSC to the horizontal plane results in a "habitual" posture of the head inclined 30° down anteriorly [1][2][3][4][5][6]8,11,12,18,46 . In the current study, the steenbuck (Raphicerus campestris) is the only species in which the average neutral and reconstructed head postures are the same (Table 1), which means that, on average, the LSC is parallel to the horizontal plane when the steenbuck's head posture is neutral. Among the species for which the neutral head posture and reconstructed cranial orientation could be compared, only half of them show a difference between the averages of the two that is below 10°. As such, even though the plane of the LSC should be horizontal on theoretical grounds 13 , this is not the rule in "ungulates". The two phylogenetic regressions provided here (Fig. 4) are the first large sample size attempts to address the existence and nature of a correlation between the orientation of the plane of the LSC and the neutral head posture in mammals across a large taxonomic sampling. We find that whether corrected or not for body mass, the correlation between the reconstructed cranial orientation and the neutral head posture is significant (p-value < 0.0001); however, if the plane of the LSC was held horizontally in neutral head posture, the regression line should not differ significantly from an isometric line. Instead, both regression lines have slopes that significantly differ from 1 (Fig. 4), which means that they cannot be approximated by isometric lines. For this reason, though there is a significant correlation between the orientation of the LSC and that of the head in "ungulates", the plane of the LSC should not be considered horizontal when reconstructing ancient head posture.
According to the phylogenetic regressions, the equation that describes the relationship between the reconstructed cranial orientation and the neutral head posture is given in Eq. (1), and that between cranial orientation, head posture, and body mass is given in Eq. (2). As estimating body mass in extinct species is always contentious [79][80][81][82] , the first of these equations may seem more practical to estimate the actual head posture of a given extinct ungulate species. In both cases, the variance of the residuals is high, which might indicate a low predictive power of these models (R 2 equals 0.26 and 0.33, respectively).
The misalignment between the plane of the LSC and the horizontal is consistent with the results obtained by Marugán-Lobón et al. 43 on birds using Duijm's dataset 8 , though their approach to the study of head posture was different from the one presented here, which limits comparisons. The reason why the LSC would not be aligned with the horizontal in the neutral posture is still unclear. It may be explained by the very function of the canals and ampullae, which are meant to record head movements and play crucial roles in the vestibulo-ocular and vestibulo-collic reflexes to compensate for the movements and accelerations of the head compared to the eyes and the rest of the body 6,13,14,43 . As such, recording movements and monitoring reflexes during locomotion and head movements along the predominant axis of yaw would be the main drivers of LSC adaptations 83 , which would result in a relaxed selection on the functions performed when the head remains still, such as aligning the plane of the LSC to the horizontal plane during neutral head posture. In support of this hypothesis, a recent study by Dunbar et al. 84 on the orientation of the LSC during locomotion in horses found a 66° inclination of the head below the horizontal during slow walk, and a higher head posture of 56° and 55° during trot and www.nature.com/scientificreports/ canter. They argued that fast locomotion brought the plane of the LSC to about 5° around the horizontal in their specimens 84 . According to Zubair et al. 85 , domestic cats keep their LSC about 10° to the horizontal during locomotion, although Hullar 12 reports a tilting up to 60° of the LSC in cats during "normal activities". Primates appear to keep the plane of their LSC within a 20° range about the horizontal during locomotion 84 , which is consistent with its orientation at rest 12 . Published data about the orientation of the LSC during locomotion are scarce and some are contradictory. Additionally, they are difficult to acquire as animals rarely maintain a static head posture during locomotion as they often pitch their heads repeatedly 84,85 . In the future, such data could nevertheless enable testing whether the orientation of the LSC would be a better predictor of the head posture during locomotion rather than at rest. Another hypothesis is that an overall misalignment of all three semicircular canals would enable all semicircular canals to record a component of horizontal and vertical accelerations 43 .
The effect of phylogeny. Both the neutral and reconstructed head postures carry an important phylogenetic signal (Lambda equals 0.97 and 0.84, respectively). For paleontologists, this strong phylogenetic signal implies that the best way to predict the head posture of an extinct "ungulate" is to look at the neutral head posture of its modern relatives. In comparison, once the data are corrected for phylogeny, diet is found to have only a weak correlation with the reconstructed cranial orientation (F = 7.723; p-value = 0.046), and no significant effect on the neutral head posture (F = 2.663; p-value = 0.516). This reflects well in the dataset.
The Tylopoda is the group in which the head is the most consistently tilted upward, with an average neutral head posture of 10°. This reflects on their reconstructed cranial orientation which averages 17°. Tylopods nevertheless include grazers (genus Vicugna), mixed feeders (Camelus bactrianus), and browsers (Camelus dromedarius) that all keep their heads relatively high (Figs. 4, 8a). On the other end of the spectrum, pigs display remarkable consistency at holding their head low (average neutral head posture for Suoidea = 46°), with for example the grazing Phacochoerus and the browsing Catagonus both keeping their head 46° below the horizontal plane on average (Table 1, Fig. 4).
The species that holds their head the lowest below the horizontal belong to the Equidae (average neutral head posture = 60°; average reconstructed cranial orientation = 49°) and Alcelaphinae (average neutral head posture = 57°; average reconstructed cranial orientation = 61°) (Figs. 4, 8b). The species with the most tilted head posture is the Grevy's zebra (Equus grevyi), with a 66° tilt on average (Table 1), comparable to the extreme 67° reconstructed cranial orientation of the sauropod dinosaur Nigersaurus 20 . Tilting of the head was hypothesized to be correlated to body size in "ungulates" 47 , with small, forest-dwelling species holding their head higher than large species adapted to savannah. A similar trend is found here, with species above 100 kg having their head tilted 37° anteriorly on average whereas species below 100 kg hold their heads 30° below the horizontal on average (Fig. 3). Statistical tests on our dataset find a significant correlation between body mass and head posture (p-value = 0.014), but this correlation is no longer significant once the data are corrected for the effect of phylogeny (p-value = 0.060). This strongly suggests that the trend observed here and by Köhler 47 actually reflects www.nature.com/scientificreports/ the fact that large savannah herbivores belong to just a few clades (e.g. equids, alcelaphins, and hippotragins) whereas the small ones mostly belong to the Antilopinae. A more significant effect of body mass might nevertheless be found while including very small-bodied species (e.g. rodents, shrews) because their head posture would be more constrained by its proximity to the substrate. Similarly, species with a sprawling posture would also have to keep their head higher, as already observed in many reptiles 12,[41][42][43] .
The relationship between LSC orientation and phylogeny has been empirically anticipated as the way the LSC enters the vestibule (either directly above the posterior ampulla or at different levels within the ampulla) is distinctive between different clades of ruminant 50,51,[86][87][88] .
The selection of an Early Burst model as the best fit for the whole dataset suggests that the evolution of head posture might represent a classic example of adaptive radiation 61,78 . This may be the effect of the abundance of bovids in our dataset, which originated in the early Neogene and rapidly adapted to a wide variety of ecological niches, diet, body mass, social behavior, and habitat 48,89 . However, the presence of a strong phylogenetic signal for both head posture variables cannot be directly interpreted as evidence for such a radiation in terms of evolutionary process, which would require further analyses of diversification rates 90 . The increase of phenotypic divergence resulting in different niches for a given character in a clade does not always correspond to an adaptive radiation, even when it closely matches the underlying phylogeny. This is due to ecological interactions between distantly related species, which can result in a similar timing of evolutionary shifts for distinct clades, the detail of which is often very difficult to decipher without exhaustively sampling each of these clades 91 . Indeed, even if an early radiation in head posture diversity would be consistent with the high discrepancy between taxonomic groups for both variables (Fig. 4), the small sample size for all groups except Ruminantia prevents a straightforward discussion of specific evolutionary constraints in that context. Furthermore, other factors such as a high proportion of sympatric species in the sample may also artificially increase the fit of an EB evolutionary model 78 . A larger and broader sampling among mammals is thus likely to blur such a signal and could result in a more homogeneous distribution of residuals that would not necessarily match so closely the observed pattern of relative phylogenetic proximity.

Diet.
As grazers have to keep their head low while foraging on grass, whereas browsers have to catch leaves higher in bushes and trees, and because herbivores spend most of their time acquiring low-energetic food 49,92,93 , it is expected on an evolutionary scale that the skull, neck musculature, and vestibular apparatus of herbivores would adapt to these different feeding strategies and that it would reflect in their head posture, even at rest 20,30,45 . As such, a gradually more anteriorly tilted neutral head posture is expected as moving from browsing species to mixed feeders, and finally grazers 20,30 . Schellhorn 30 was the first to compare the reconstructed cranial orientation of Rhinocerotidae to their actual head posture (using an open-access database of photographs) and found results consistent with such a gradient.
Measuring the reconstructed cranial orientation from Schellhorn's published figures and adding our own observations of neutral head posture, we do find a more tilted reconstructed cranial orientation in the grazing Cerathotherium (average reconstructed cranial orientation = 38°), than in the mixed feeder Rhinoceros (average reconstructed cranial orientation = 34°), and the browsing Dicerorhinus and Diceros (average reconstructed cranial orientation = 31°) (Fig. 9); however, this gradient does not reflect on the neutral head posture that shows no particular trend among rhinocerotids (Fig. 9). Despite the large difference between the average head posture of browsers and grazers, the very low head posture of mixed feeder rhinocerotids casts some doubts on the validity of the correlation between head posture and diet (Fig. 9). The Cervidae constitutes a more striking example (Fig. 10). In cervids, the head is consistently kept within about 10° around the average neutral head posture (20°) regardless of diet, and browser and grazers are not distinguishable based on neutral head posture (average = 30° for browser; average = 28° for grazers) or reconstructed cranial orientation (average = 41° for browsers; average = 45° for grazers) ( Table 1).
The situation in rhinocerotids illustrates that for the whole dataset. Diet is found to have no statistical effect on neutral head posture (p-value = 0.516), and its relationship to the reconstructed cranial orientation is barely significant (p-value = 0.046). Average values for neutral and reconstructed head postures show a visible increase in tilting with a more grass-rich diet, but there is a strong overlap between dietary groups for both variables (Fig. 5). This suggests that even though diet could potentially be reconstructed in extinct species using the orientation of the LSC, caution should be taken as i) browsers and grazers could be statistically discriminated, but mixed feeders could not; ii) the correlation between reconstructed cranial orientation and diet does not seem to reflect on the neutral head posture. The reason why remains unknown; and iii) the high p-value suggests that adding more data (particularly CT data) in the future may affect this correlation.
Semi-aquatic adaptation. No significant result indicative of a correlation between habitat and head posture or reconstructed cranial orientation was found; however, semi-aquatic species show a noticeably high reconstructed posture of the skull on average (23°) compared to other species (Fig. 6). The low number of semiaquatic species in the dataset (one rhinocerotid, three tapirids, and two hippopotamids, Table 1) likely prevents this trend to be identified as significant in our sample. A high head posture is not observed in semi-aquatic species (Fig. 6), even though it would be expected of species that have to keep breathing above water level most of the time while immersed 27,30 . Recently, a semi-aquatic habit for the Triassic archosaur Proterosuchus and the therapsid Lystrosaurus has been hypothesized as these two species would have their head tilted upward anteriorly when the plane of the LSC is horizontal 27,94 (Fig. 11c). In Proterosuchus this upward tilting would be about 17°2 7 , whereas in Lystrosaurus it would be between 19° and 23°9 4 . However, an upward tilting of the neutral head posture is also habitually observed in many fully terrestrial species, such as Camelus dromedarius and was occa-   www.nature.com/scientificreports/ sionally spotted in Capra ibex ( Fig. 11a; see Supplementary Table S1). Among modern archosaurs, an upward tilting of the head is observed in some sea birds and the common starling Sturnus vulgaris 8 . An upward tilting of the beak when the dry skull is held with the plane of the LSC horizontal seems to be observed in the razorbill (Alca torda) and the Heron (Ardea cinerea) (as inferred from the figures and measurements in Duijm 8 , see Supplementary Table S1); however, the reconstructed cranial orientations in the razorbill and heron vary between 0° and -4° only (Fig. 11b). Among fossil species, the likely terrestrial sauropod Ngwevu intloko (Fig. 11c) would also have had its reconstructed cranial orientation tilted 17° above the horizontal (identified as Massospondylus carinatus in Sereno et al. 20 ; see Chapelle et al. 95 ). The peculiar orientation of the LSC in Lystrosaurus, Proterosuchus, and Massospondylus is remarkable as such a downward tilting of the LSC superior to 15° has never been found in any modern species to date, particularly not in the semi-aquatic species studied here which all have a posteriorly tilted LSC as all other "ungulates" (average reconstructed cranial orientation of semi-aquatic species = 23°) ( Fig. 6; Table 1). The heron and razorbill mentioned above would have an anterior tilting of less than 4° (Supplementary Table S1), and Duijm's 8 dataset includes mostly semi-aquatic species, which limits comparisons. Overall, a correlation between an upward tilting of the LSC and semi-aquatic lifestyle is not supported by current data.
Noteworthily, the most aquatic species of the dataset, Hippopotamus amphibius, stands out on the scatter-plots as an outlier (Fig. 4). This is due to the large difference between its strongly anteriorly tilted neutral head posture and its almost horizontal reconstructed cranial orientation (difference = 39°). Hippopotamus amphibius normally spends very little time on land 49 , and it can be hypothesized that its bony labyrinth morphology would be more adapted to life in water than on land 96 . This hypothesis is supported by the orientation of the head while swimming in H. amphibius, which is more consistent with its reconstructed cranial orientation (Fig. 12), and the fact that the difference between the neutral and reconstructed head postures in the more terrestrial Hexaprotodon liberiensis (23°) falls more within the range of variation of "ungulates" (Fig. 4). Further field observations of underwater hippopotamus head posture will be necessary to address this hypothesis.
The influence of head-butting. Taxa that engage in head-to-head combat have almost the same average neutral orientation of their head as non-head-butting taxa (34° and 33° respectively) (Fig. 7). In sharp contrast, they show a significantly more anteriorly tilted reconstructed cranial orientation (average = 46°) compared to non-head-butting taxa (average = 33°). Unlike what is observed between browsers and grazers, the values here are markedly different (Fig. 7) and the difference is highly statistically significant (p-value = 0.002). It is unlikely that the reconstructed cranial orientation reflects the posture during head-butting, as animals keep their head Figure 13. Comparison of the cranial (transparent), endocranial (pink), and LSC (green) orientations in a nonhead-butting species (Tapirus indicus, a) and a head-butting species (Connochaetes taurinus, b).
Scientific Reports | (2020) 10:19602 | https://doi.org/10.1038/s41598-020-76757-0 www.nature.com/scientificreports/ extremely low during this activity 49,92,93 , much lower than their corresponding reconstructed cranial orientation (Fig. 13b). A phenomenon of re-orientation of the braincase and basicranium in head-butting "ungulates" that would not affect head posture overall seems more probable. Such cranial flexure would result in a misalignment of the main axis of the braincase with that of the snout (Fig. 13), a condition termed cyptocephaly and commonly encountered in head-butting "ungulates" 33,48,92,[97][98][99] . This implies that natural selection on the alignment of the plane of the LSC to the horizontal was muted by another, likely more important adaptation to head-butting. This may be the necessity to re-orientate the braincase and basicranium in order to align the fighting surface of the skull, occipital condyles, and vertebral column to help dissipate the energy of the impact to the body and away from the brain 24,100,101 . Another hypothesis would be that the re-organization of the braincase serves to accommodate the development of the large cranial apparatuses found in most head-butting species (e.g. horns and antlers) 29 . In our dataset, 56% of species without horns do not head-butt, and 83% of species with horns do head-butt (Supplementary Table S1), which would partly support this hypothesis. In contrast, intraspecific variations in the presence or absence of horns appear to have no impact on neutral head posture on average ( Fig. 14; Supplementary Table S1). These two hypotheses will need further observations and biomechanical modeling to be addressed adequately.

Concluding remarks
Neutral head posture is here found to be significantly correlated to the orientation of the plane of the LSC in "ungulate" mammals, but this relationship is loose, and it appears that diet and head-butting have an effect on LSC orientation although not on neutral head posture as would be expected. This suggests an overall relaxed constraint on the alignment of the plane of the LSC to the horizontal at rest. Head posture during locomotion and/or adaptation to head-butting might play a more significant role in the orientation of the LSC than its horizontality at rest, two possibilities that will have to be addressed further.
In this contribution, some noteworthy trends between the orientation of the LSC, body mass, diet, adaptation to a semi-aquatic environment, and head-butting are pointed out, although many of these ecological components are difficult to disentangle. "Ungulates" living in closed habitats are often smaller than in open habitats, more solitary, browsers and tend to fight for mates by stabbing each other with their horns and teeth, whereas species from more open habitats usually graze in large herds and perform head-butting to ascertain dominance and attract mates 34,36,40,47 .
Finally, although this study finds that there is some interesting ecological and behavioral signal in the orientation of the LSC of ungulates that could be exploited by paleontologists, it is crucial to highlight that the phylogenetic signal was highly significant for all the variables examined here and as such, what the orientation www.nature.com/scientificreports/ of the LSC reflects the best in "ungulates" is their phylogeny more than anything else. Further understanding of the evolutionary processes associated with such a strong phylogenetic disparity will require investigating each subclade in the sample individually and a more exhaustive sample for each of them.

Data availability
The datasets analysed during the current study are available in the Supplementary