Cyril Etienne
email: cyril.etienne@cri-paris.org

Christophe Mallet

Raphael Cornette

Alexandra Houssaye

Raphaël Cornette

Influence of mass on tarsus shape variation: a morphometrical investigation among Rhinocerotidae (Mammalia: Perissodactyla)

Keywords:

come

Introduction

In vertebrate locomotion, bone is a rigid organ of paramount importance. It provides support for the body as well as attachment points for the muscles, via the tendons [START_REF] Hildebrand | Analysis of vertebrate structure[END_REF][START_REF] Biewener | Biomechanics of mammalian terrestrial locomotion[END_REF]. Bone shape varies with a diversity of factors, one of the main ones being the size of the animal [START_REF] Hildebrand | Functional vertebrate morphology[END_REF][START_REF] Biewener | Mammalian terrestrial locomotion and size[END_REF][START_REF] Polly | Limbs in mammalian evolution[END_REF][START_REF] Biewener | Physical and Biological Properties and Principles[END_REF]. Evolutionary convergences are usually observed when similar selective pressures are applied to the same structure independently in different groups. Accordingly, in a given clade, an increasing mass generally results in e.g. a more vertical orientation of the pelvis [START_REF] Polly | Limbs in mammalian evolution[END_REF], an increasing diameter of the femur [START_REF] Alexander | Body support, scaling, and allometry[END_REF], and micro-anatomical changes such as a high thickness of cortical bone [START_REF] Houssaye | Hyperspecialization in some South American endemic ungulates revealed by long bone microstructure[END_REF]. Postural and locomotor factors such as facultative bipedalism or cursoriality are also important factors influencing the shape of the skeleton, [START_REF] Hildebrand | Analysis of vertebrate structure[END_REF][START_REF] Polly | Limbs in mammalian evolution[END_REF]. Analysing the relationship between mass and the shape of the bones, while considering also factors such as posture and locomotion, would allow a better understanding of the way animals with different body plans adapt to an increasing mass.

Rhinocerotidae [START_REF] Gray | On the natural arrangement of vertebrose animals[END_REF] seem to be an excellent group to study morphological variations in bone related to mass and to a varying body plan. Today comprising five species and four genera, rhinocerotids are found only in tropical regions [START_REF] Dinerstein | Family Rhinocerotidae (Rhinoceroses)[END_REF]. They were much more diverse during the Cenozoic, appearing in the middle of the Eocene and comprising more than a hundred species [START_REF] Cerdeño | Diversity and evolutionary trends of the family Rhinocerotidae (Perissodactyla)[END_REF]. Rhinocerotids have been found in Eurasia, Africa and North America (as far south as Panama; [START_REF] Macfadden | North American Miocene land mammals from Panama[END_REF]), and existed in a diversity of habitats, e.g. cold steppes, dense forests or swamps [START_REF] Prothero | The history of the Rhinocerotoidea[END_REF][START_REF] Mörs | Biostratigraphy and paleoecology of continental Tertiary vertebrate faunas in the Lower Rhine Embayment (NW-Germany)[END_REF][START_REF] Prothero | Biogeography and diversity patterns[END_REF]. All of them were relatively heavy animals as compared to the average body mass of mammals [START_REF] Gardezi | Diversity in relation to body size in mammals: a comparative study[END_REF], ranging from approximately 150 kg for the lightest taxa of the Eocene to five tons for Elasmotherium [START_REF] Cerdeño | Diversity and evolutionary trends of the family Rhinocerotidae (Perissodactyla)[END_REF][START_REF] Becker | Paléoécologie et paléoclimats de la Molasse du Jura (Oligo-Miocène) : apport des Rhinocerotoidea (Mammalia) et des minéraux argileux[END_REF]Antoine in press).

They showed several independent extreme increases of body mass, up to more than two tons, during their evolutionary history (e.g. in Rhinocerotina, Teleoceratina or Elasmotheriinae; [START_REF] Cerdeño | Diversity and evolutionary trends of the family Rhinocerotidae (Perissodactyla)[END_REF]. They also vary in terms of body plan: some taxa are massive and sturdy (e.g.

Coelodonta, Ceratotherium), some are extremely short-legged (e.g. Teleoceras, Brachypotherium, Prosantorhinus in the Teleoceratina tribe), others have been described as gracile and cursorial (e.g. Protaceratherium, Hispanotherium, [START_REF] Cerdeño | Diversity and evolutionary trends of the family Rhinocerotidae (Perissodactyla)[END_REF]. Rhinocerotids should therefore be a very interesting group to study bone shape variations, and analyse its relationship with different masses and morphologies. Rhinocerotidae are part of the order Perissodactyla [START_REF] Owen | Description of teeth and portions of jaws of two extinct Anthracotherioid quadrupeds (Hyopotamus vectianus and Hyop. bovinus) discovered by the Marchioness of Hastings in the Eocene deposits on the N.W. coast of the Isle of Wight: with an attempt to develop Cuvier's idea of the classification of pachyderms by the number of their toes[END_REF]. If they nowadays only include the five species of Rhinocerotidae, the seven species of Equidae and the four species of Tapiridae, they were much more diversified during the Cenozoic, in terms of number of species and families (Prothero & Schoch [1989], see figure 1). They included notably the Paraceratheriidae (sometimes considered a subfamily of Hyracodontidae, see [START_REF] Wang | Earliest known unequivocal rhinocerotoid sheds new light on the origin of giant rhinos and phylogeny of early rhinocerotoids[END_REF], which included some of the heaviest land mammals that ever lived [START_REF] Prothero | Rhinoceros giants: the paleobiology of indricotheres[END_REF]. Perissodactyla also included the intriguing Chalicotheriidae, among which some species (the subfamily Chalicotheriinae) were facultative bipeds with a gorilla-like stance, with very short hindlimbs and walking on the knuckles of their forelimbs, whereas others (the Schizotheriinae) had front and hindlimbs of approximately equal length [START_REF] Coombs | Large mammalian clawed herbivores: a comparative study[END_REF][START_REF] Semprebon | Potential bark and fruit browsing as revealed by stereomicrowear analysis of the peculiar clawed herbivores known as chalicotheres (Perissodactyla, Chalicotherioidea)[END_REF]. This order therefore encompasses a great diversity in terms of mass and body plan that can be related to their bone shape, and compared with the diversity observed among Rhinocerotidae.

Our study focused on two bones of the tarsus: the astragalus and the calcaneus. These bones are at the junction between the hind zeugopodium and the autopodium, and are essential to the movement of the foot and consequently of the entire animal. The astragalus serves as the pivot, or fulcrum, and the calcaneus as the lever arm of the foot [START_REF] Carrano | Morphological indicators of foot posture in mammals: a statistical and biomechanical analysis[END_REF]. These bones have been extensively studied, from taxonomic and phylogenetic points of view (see [START_REF] Gladman | A calcaneus attributable to the primitive late Eocene anthropoid Proteopithecus sylviae: phenetic affinities and phylogenetic implications[END_REF][START_REF] Guérin | Les rhinocéros (Mammalia, Perissodactyla) du Miocène terminal au Pléistocène supérieur en Europe occidentale : comparaison avec les espèces actuelles[END_REF][START_REF] Missiaen | Asian gliriform origin for arctostylopid mammals[END_REF][START_REF] Stains | Use of the calcaneum in studies of taxonomy and food habits[END_REF], but also in a morphofunctional context, with multiple studies trying to link their shape to the animal's mass (see [START_REF] Dagosto | Estimating the body size of Eocene primates: a comparison of results from dental and postcranial variables[END_REF][START_REF] Martinez | The astragalus of Paleogene artiodactyls: comparative morphology, variability and prediction of body mass[END_REF][START_REF] Tsubamoto | Estimating body mass from the astragalus in mammals[END_REF], habitat [START_REF] Degusta | A method for inferring paleohabitats from the functional morphology of bovid astragali[END_REF][START_REF] Plummer | Habitat preference of extant African bovids based on astragalus morphology: operationalizing ecomorphology for palaeoenvironmental reconstruction[END_REF][START_REF] Curran | Expanding ecomorphological methods: geometric morphometric analysis of Cervidae post-crania[END_REF][START_REF] Barr | Functional morphology of the bovid astragalus in relation to habitat: controlling phylogenetic signal in ecomorphology[END_REF], and mode of locomotion [START_REF] Nakatsukasa | Functional morphology of the postcranium and locomotor behavior of Neosaimiri fieldsi, a Saimiri-like Middle Miocene platyrrhine[END_REF][START_REF] Panciroli | Correlates between calcaneal morphology and locomotion in extant and extinct carnivorous mammals[END_REF]. These studies all found a link between the mass of a species and the shape of its astragalus and of its calcaneus, generally represented by linear measurements or ratios. These studies concerned a great variety of mammals, e.g. bovids, cervids, carnivorans, primates, but none ever specifically studied the relationship between mass and shape on both ankle bones in perissodactyls.

In the present study, we investigated the variation of the shape of the astragalus and calcaneus across a diversity of extant and fossil Rhinocerotidae and additional Perissodactyla.

Our primary aim was to identify shape variations associated with an increase of mass. We expect that mass will have a strong influence on those bones since they are extremely important for the support and movement in mammals. We expect that bones of large animals will be more robust and more resistant, with wider and flatter articular facets to help dissipate forces. We also expect that adaptations will vary according to the general body plan and mode of locomotion of the animal. We studied variations of bone shape, and tested the influence of the bone size and of the mean mass of each species. We first focused on shape variations across Rhinocerotidae, and then across all Perissodactyla sampled in order to compare the variations observed among Rhinocerotidae to more diverse forms (e.g. Paraceratherium), and thus better interpret the drivers acting on this variation.

Material and methods

Material

We studied 112 astragali and 94 calcanei belonging to 43 different species across five different families of Perissodactyla, with varying masses, morphologies or locomotor modes (Appendix 1; Fig. 1). Taxa were chosen in order to encompass as much as possible of the variation within Rhinocerotidae. A few specimens of other families of Perissodactyla with particular characteristics (e.g. an extremely high mass for Paraceratherium, cursoriality for horses, shortened hindlimbs for Chalicotheriinae) were included to provide a comparison point for the shape variations observed in Rhinocerotidae. Rhinocerotidae constitutes our main subsample, with 40 specimens of each bone type for living species. Our sample also includes 43 astragali and 31 calcanei of fossil Rhinocerotidae, including small, cursorial genera (Protaceratherium, Pleuroceros), and at least three lineages presenting independent rises of body mass above two tons (in Elasmotherium, Brachypotherium, Coelodonta, plus the living Ceratotherium and Rhinoceros unicornis). We studied 10 astragali and calcanei of extant [START_REF] Antoine | Phylogénie et évolution des Elasmotheriina (Mammalia, Rhinocerotidae)[END_REF][START_REF] Antoine | A revision of Aceratherium blanfordi Lydekker, 1884 (Mammalia: Rhinocerotidae) from the Early Miocene of Pakistan: postcranials as a key[END_REF], Antoine et al. in press, Holbrook and Lapergola, 2011[START_REF] Piras | Functional and phylogenetic constraints in Rhinocerotinae craniodental morphology[END_REF][START_REF] Steiner | Molecular phylogeny and evolution of the Perissodactyla: phylogeny of the perissodactyls[END_REF]. Occurrence dates were estimated using [START_REF] Antoine | Aegyrcitherium beonensis, nouvel élasmothère (Mammalia, Perissodactyla) du gisement miocène (MN 4b) de Montréal-du-Gers (Gers, France). Position phylogénétique au sein des Elasmotheriini[END_REF][START_REF] Piras | Functional and phylogenetic constraints in Rhinocerotinae craniodental morphology[END_REF][START_REF] Antoine | A revision of Aceratherium blanfordi Lydekker, 1884 (Mammalia: Rhinocerotidae) from the Early Miocene of Pakistan: postcranials as a key[END_REF][START_REF] Geraads | A new rhinoceros, Victoriaceros kenyensis gen. et sp. nov., and other Perissodactyla from the middle Miocene of Maboko, Kenya[END_REF][START_REF] Guérin | Anisodon grande (Perissodactyla, Chalicotheriidae) de Sansan[END_REF][START_REF] Prothero | Rhinoceros giants: the paleobiology of indricotheres[END_REF], as well as data recorded on http://fossilworks.org. 1: Elasmotheriinae; 2: Rhinocerotinae; 3: Aceratheriina; 4: Teleoceratina; 5: Rhinocerotina; 6: Schizotheriinae; 7: Chalicotheriinae. Species in bold are living species. The relationships between the five extant rhino species are still debated, specifically the position of Dicerorhinus sumatrensis and its close fossil relatives is uncertain, placed either as sister taxa to Ceratotherium and Diceros, as sister taxa to Rhinoceros, or as sister taxa to a group formed of Rhinoceros, Ceratotherium and Diceros [START_REF] Willerslev | Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution[END_REF][START_REF] Gaudry | Molecular phylogenetics of the rhinoceros clade and evolution of UCP1 transcriptional regulatory elements across the mammalian phylogeny[END_REF]. This phylogenetic uncertainty is here represented by a polytomy.

Mass estimations used

Mass data were retrieved from the literature (Table 1). Methods of mass reconstructions are detailed in the references; they usually relied on regression equations and measurements on the molars or on proximal limb segments. None of them used measurements on the astragalus or the calcaneus. When only a range of masses was available with no average, the mean of the minimal and maximal mass was used.

Data acquisition

The specimens were digitized using either an Artec Eva surface scanner and the Artec Studio Professional v12.1.5.1 software (Artec 3D, 2018), or a Nikon D5500 camera (automatic mode, without flash, focal length 50 mm, aperture f/1.8) and the photogrammetry software Agisoft PhotoScan v1.4.0 (Agisoft LLC, 2017). The 3D meshes were then exported, decimated down to 60,000 faces and mirrored to have only right side astragali and calcanei, using MeshLab v2016.12 [START_REF] Cignoni | MeshLab: an open-source mesh processing tool[END_REF]. In two cases (The astragalus of Brachypotherium snowi NHM-PAL-PV-M-29279 and the calcaneus of Hispanotherium beonense MHNT-2015-0-837), the specimens were slightly damaged where a curve would pass and had to be partially reconstructed using Geomagic (3D Systems, 2017).

Geometric morphometrics

Bone shape was modelled using anatomical landmarks and semi-landmarks sliding on curves (Gunz & Mitteroecker, 2013). Landmarks were all placed by the same operator (C.E.).

Given that there can be marked differences in bone shape between rhinocerotids and the other perissodactyls, we split the analysis in two. Two sets of landmarks and curves were therefore defined: one for all the rhinocerotids, and another for all the perissodactyls (see Figs. 2,3;

Appendix 3 for descriptions of the landmarks and curves), with fewer landmarks and curves but able to encompass a broader number of taxa. The second set is mostly a subset of the first one, only two curves had to be redefined. Landmarks and curves were digitized on the meshes using the IDAV Landmark software package [START_REF] Wiley | Landmark. Institute for data analysis and visualization[END_REF]. All the analyses and statistical tests were run using R (R Development Core Team, 2005) and RStudio [START_REF] Rstudio | RStudio[END_REF]. The curves were resampled using the algorithm provided in [START_REF] Botton-Divet | Morphological analysis of long bones in semi-aquatic mustelids and their terrestrial relatives[END_REF]. Then, as the algorithm can result in some semi-landmarks being slightly above or below the mesh surface, the semi-landmarks were reprojected on the meshes using the closemeshKD function of the Morpho R package [START_REF] Schlager | Morpho: a toolbox providing methods for dataacquisition, visualisation and statistical methods related to Geometric Morphometrics and shape analysis[END_REF], which uses the coordinates of each semi-landmark to calculate its closest match on the surface of the mesh.

Landmarks were superimposed using a Generalized Procrustes Analysis (GPA), which translates, scales and rotates each set of landmarks in order to remove the information of size, position and angle, and to minimize the sum of square distances between landmark configurations [START_REF] Bookstein | Morphometric tools for landmark data: geometry and biology[END_REF]. The curve semi-landmarks were slid along the curves in order to minimize the bending energy of a Thin-Plate-Spline as described in [START_REF] Gunz | Semilandmarks in three dimensions[END_REF]. The bending energy is a scalar quantity that roughly represents the amount of local shape deformation between a reference set of landmarks (chosen arbitrarily among our sample). More technically, it is the integral of the squared second derivatives of the deformation (see [START_REF] Gunz | Semilandmarks: a method for quantifying curves and surfaces[END_REF][START_REF] Mitteroecker | Advances in geometric morphometrics[END_REF]. To assess the repeatability of the landmarks, and prior to placing the landmarks on the whole sample, we placed each landmark ten times on each of three specimens of Diceros bicornis, alternating between each specimen. The three specimens were assessed by sight to be the three morphologically closest. These 30 landmark sets were then superimposed using a Generalized Procrustes Analysis and visualized using a Principal Components Analysis (PCA), to check that landmark error per specimen was smaller than inter individual variation (Appendix 4).

Statistical analyses

After the GPA, the aligned landmarks coordinates were used in a PCA, in order to reduce dimensionality of our data and assess the shape variation patterns in our sample. Neighbour-Joining trees were generated using a Euclidian distance matrix based on the PC-scores, in order to visualize the phenotypic similarities between each specimen or group in a multivariate manner, instead of one axis at a time, which is useful if each axis explains a small percentage of variance. PC-scores were used instead of Procrustes coordinates to reduce number of dimensions and thus lower the computing power required.

We tested the influence of the centroid size of each bone on its shape. Centroid size is defined as the square root of the sum of the square of the distance of each point to the centroid of the landmark set; it is most commonly used to assess the variations of shape that are due to variations of size, or allometry [START_REF] Mitteroecker | A brief review of shape, form, and allometry in geometric morphometrics, with applications to human facial morphology[END_REF][START_REF] Klingenberg | Size, shape, and form: concepts of allometry in geometric morphometrics[END_REF]. Logarithms were used for the centroid size values, as recommended by [START_REF] Bookstein | Morphometric tools for landmark data: geometry and biology[END_REF] and [START_REF] Klingenberg | Multivariate allometry[END_REF]. Procrustes coordinates were correlated against centroid size using a multivariate regression. The allometry-free residuals from the tests were used to create allometry-free shapes for each individual, allowing analyses where the influence of size is entirely removed (see e.g. [START_REF] Evin | Phenotypic diversification and island evolution of pipistrelle bats (Pipistrellus pipistrellus group) in the Mediterranean region inferred from geometric morphometrics and molecular phylogenetics[END_REF][START_REF] Perrard | Differences in caste dimorphism among three hornet species (Hymenoptera: Vespidae): forewing size, shape and allometry[END_REF].

The centroid size of both bones is statistically linked to the mean mass of the species in our sample (see Appendix 5; p < 0.0001, R² between 0.46 and 0.82 depending on the bone and the landmark set). The R-square is however different from 1, therefore mass could have an influence on the shape of the bones that is independent from its centroid size. For example, two astragali or calcanei with the same centroid size belonging to species with different mean masses would exhibit divergent shapes. This was tested using a multivariate regression of the allometry-free shapes generated earlier on the logarithm of the cubic root of the mean mass of the species. Given that we could not find mass estimation for 14 sampled species, they have been removed from this analysis.

In order to assess what the effect of shared evolutionary history of different species is on the shape of the astragalus and calcaneus, the degree of phylogenetic signal in the morphological data was also assessed, using a multivariate K statistic (K-mult) based on the PC-scores. It compares the observed rate of morphological change to the expected change under a Brownian motion (see [START_REF] Adams | A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data[END_REF][START_REF] Blomberg | Testing for phylogenetic signal in comparative data: behavioral traits are more labile[END_REF]. The phylogeny used is provided in figure 1. Results are provided in Appendix 6.

Thin Plate Splines were used to visualize the results of our analyses: for each set of landmarks on the calcaneus and astragalus, the mean-shape generated by the GPA was mapped onto the specimen closest to the mean value. Then, this mean-shaped model was deformed towards the shape resulting from our analyses, for instance the shape extremes of each PCA axis.

Results

Rhinocerotidae

Astragalus

Morphological variations

The Neighbour-Joining tree (Fig. 4) generally shows a greater morphological proximity between members of the same species than between members of different species, which indicates that interspecific variation is generally greater than intraspecific variation. The Teleoceratina, the short-legged rhinocerotids, are clearly grouped and separated from the others, except Diaceratherium. Among them, Teleoceras has a very long branch indicating a very derived morphology for this individual. Dicerorhinus, Ceratotherium simum and Rhinoceros all form homogenous groups, but Diceros has two specimens that are separated from the others.

This might be due to them possibly belonging to different subspecies or a different sex, but this is unknown for these specimens. All the fossil dicerorhinins (Dihoplus, Stephanorhinus and Coelodonta) are grouped together, but Dicerorhinus sumatrensis is separated from them.

Aceratheriina, Elasmotheriinae, Diaceratherium, Protaceratherium, Pleuroceros, and

Plesiaceratherium tend to all group together, with a few exceptions. The two specimens of Iranotherium morgani are clearly separated from one another in the tree; this is also the case for the two Ceratotherium neumayri specimens. Above the tribe level, there are no clusters that seem to follow the phylogeny.

A low percentage of variance is explained by each axis of the PCA on Rhinocerotidae astragali (63.1% for the first ten axes). Only the principal components that are correlated to centroid size, or that highlight variations of shape that could be due to differences in terms of animal general body plan, are described. Thus, only PC1, PC2 and PC4 are described here. PC1

highlights differences in astragalus shape between species with different morphologies, PC2

and PC4 are weakly but significantly correlated with centroid size. PC2 is positively correlated with centroid size (p < 0.001, R² = 0.17), and PC4 is negatively correlated with centroid size (p < 0.001, R² = 0.12). Vector representations of the deformations along the principal components are provided in Appendix 7. Teleoceratina, on the negative side, from the other Rhinocerotidae. The axis is characterized in its negative extremity by a great proximo-distal compression of the bone; a flatter and symmetrical trochlea with medio-laterally wider and lower ridges; and a trochlea with a more proximal than anterior orientation. As for the articular facets, in the negative side of the PC1 the proximal facet with the calcaneus is distally elongate; the medial facet for the calcaneus is proximo-distally compressed, twice as broad (medio-laterally) as high (proximo-distally), and not fused with the distal facet; the distal facet for the calcaneus is medio-laterally very short;

and the facet for the navicular is broader than on the positive extremity of the axis, covering most of the distal face of the astragalus.

Impact of allometry and mass

The centroid size has a significant but weak effect on the shape of the astragalus (p < 0.01 and R² = 0.04, according to multivariate regression of the logarithm of the centroid size on the Procrustes coordinates). A large astragalus (Fig. 6A) is characterized by a medio-laterally wider and triangle-shaped proximal facet with the calcaneus; medio-laterally wider and fused medial and distal facets with the calcaneus; and an articular facet with the navicular positioned less laterally offset, more directly underneath the rest of the bone.

Once the influence of the centroid size is removed, there is only a weak influence of the mass of the species on the shape of the astragalus (p < 0.05, R² = 0.03). The shape variations are minimal: in an astragalus pertaining to a heavy species (Fig. 6B), the facets for both malleoli are enlarged, the crescent they form being wider; the proximal facet for the calcaneus is slightly more triangle-shaped; the medial facet is slightly wider medio-laterally; the facet for the cuboid and the distal facet for the calcaneus are anteriorly extended.

Calcaneus

Morphological variations

The Neighbour-Joining tree on the calcaneal morphology (Fig. 7) indicates, as for the astragalus, a tendency for individuals of the same species to be grouped together. Teleoceratina are grouped with Aceratheriina, and Teleoceras fossiger has again the longest branch of all species. Elasmotherium sibiricum has also a particularly long branch. There seem to be fewer clusters of species belonging to the same higher-rank taxon than for the astragalus. The genus Rhinoceros, notably, sees its three species scattered. Ceratotherium neumayri is close to Chilotherium or Elasmotherium sibiricum, and most of extinct dicerorhinins are grouped with Rhinoceros palaeindicus, Hispanotherium beonense and Plesiaceratherium. Again, for taxa of higher rank than tribes, there are not clusters following the phylogeny.

Impact of allometry and mass

The centroid size has a significant but weak effect on the calcaneus shape (p < 0.01 and R² = 0.04; multivariate regression). A large calcaneus (Fig. 9A) has a medio-laterally wider proximal part of the proximal facet for the astragalus; a wider medial facet for the astragalus, expanding distally and merging with the distal facet for the astragalus; an elongated distal facet for the astragalus; and a sustentaculum tali oriented more distally, whereas it is oriented anterodistally in small calcanei.

The mass has a slightly stronger influence on the allometry-free shapes of the calcaneus than on those of the astragalus (p < 0.001, R² = 0.06). A calcaneus belonging to a heavy species (Fig. 9B) has a more robust tuber calcanei than a calcaneus belonging to a light species; the proximal facet for the astragalus is more triangular, and slightly extended medially; the distal facet for the astragalus is extended distally; and the facet for the cuboid is piriform, and wider proximally.

Perissodactyla Astragalus

The Neighbour-Joining tree on the astragalar morphology (Fig. 10

Impact of allometry and mass

The centroid size of the astragalus has a significant effect on its shape (p < 0.001 and R² = 0.14 for the astragalus, multivariate regression). A large-sized astragalus (Fig. 12A) is characterized by an overall flat bone, twice as wide medio-laterally as high proximo-distally; medio-laterally wide and low trochlear ridges; a trochlea oriented proximally; a medio-laterally wide and triangle-shaped proximal facet for the calcaneus; a round-shaped medial facet for the calcaneus; a wide facet for the navicular, flat overall and positioned below the body of the astragalus. A small sized astragalus is as wide as it is high, has higher trochlear ridges and a deeper trochlear groove; a more square-shaped proximal facet for the calcaneus, with a small latero-distal extension; a rectangle-shaped medial facet for the calcaneus, higher proximodistally than wide medio-laterally; and a smaller facet for the navicular, not directly below the body of the astragalus but medially offset.

Species mass has a statistically significant but very weak effect on allometry-free astragalus shape (p < 0.05, R² = 0.027). An astragalus pertaining to a heavier species (Fig. 12B) has a flatter trochlea with lower ridges; a triangle-shaped proximal facet for the calcaneus; and a medial facet for the calcaneus located more laterally.

Calcaneus

The Neighbour-Joining tree on the calcaneus morphology (Fig. 13) also shows a relative

Impact of allometry and mass

The centroid size of the calcaneus has a significant influence on its shape (p < 0.001 and R² = 0.11, multivariate regression). A large-sized calcaneus (Fig. 15A) is characterized by an extreme proximo-distal compression, the tuber calcanei being very robust; a much wider proximal facet for the astragalus, extended in all directions, especially in its proximal half; a distally-oriented sustentaculum tali and medial facet for the astragalus; and a wider facet for the cuboid, triangle-shaped and latero-posteriorly extended. A small-sized calcaneus has a very thin tuber calcanei as compared to large-sized ones; a relatively much smaller overall proximal facet for the astragalus; an anteriorly oriented sustentaculum tali and medial facet for the astragalus; and a relatively smaller facet for the cuboid.

There is a statistically significant influence of species mass on allometry-free calcaneus shape (p < 0.001, R² = 0.09). Shape differences are clear (Fig. 15B), unlike those observed for that same analysis on Rhinocerotidae alone. In our sample, a calcaneus belonging to a heavier species is, on average, characterized by a stouter tuber calcanei; a wider overall proximal facet for the astragalus; a slightly wider medial facet for the astragalus, oriented distally along with the whole sustentaculum tali; and a wider facet for the cuboid, expanding more proximally.

Discussion

First and foremost, it is worth noting that the percentage of variance explained by the first axes of the PCA is usually low (around 66% for the first five axes for the Perissodactyla dataset, 40% for the Rhinocerotidae dataset). The first four or five axes describe the variations between species or families, but the following axes usually distinguish one or two individuals from other specimens of the same taxon. There is indeed a great intraspecific variation in the species studied, even if it remains inferior to inter specific variations (Figs. 4,7). For example, Dicerorhinus sumatrensis ZSM-1908-571 presents an astragalus with a wide medio-distal extension of the medial facet for the calcaneus, an extension absent in all the other specimens.

Some specimens of Ceratotherium simum present a calcaneus with a fusion of the medial and distal facets for the astragalus, whereas in others they are separated by a deep groove. This could explain the low PC-score values. Indeed, if there are many different variations observed in the sample, these cannot be described on one PC alone and thus the percentage of variation explained by the first axes lowers. On the other hand, if there were a factor clearly driving a continuum of variations in all our sample, we would see a higher percentage of variance for the first axis. That factor is often size (see e.g. [START_REF] Bonnan | What lies beneath: sub-articular long bone shape scaling in eutherian mammals and saurischian dinosaurs suggests different locomotor adaptations for gigantism[END_REF][START_REF] Cardini | Why the long face? Kangaroos and wallabies follow the same 'rule' of cranial evolutionary allometry (CREA) as placentals[END_REF][START_REF] Knigge | Three-dimensional geometric morphometric analysis of talar morphology in extant gorilla taxa from highland and lowland habitats[END_REF]; here it seems clear that size does not have a strong influence on the shape of the bones, especially in our Rhinocerotidae dataset. We already noted qualitatively this intraspecific variability between individuals of the same species of Rhinocerotidae when digitizing the bones, and it has also been observed by [START_REF] Guérin | Les rhinocéros (Mammalia, Perissodactyla) du Miocène terminal au Pléistocène supérieur en Europe occidentale : comparaison avec les espèces actuelles[END_REF] on various bones of the tarsus of extant rhinocerotids, by [START_REF] Harrison | Extreme carpal variability in Teleoceras (Rhinocerotidae, Mammalia)[END_REF] on the carpus bones of Teleoceras, and by [START_REF] Heissig | Les Rhinocerotidae (Perissodactyla) de Sansan[END_REF] on several limb bones, including the astragalus, of aceratheres. Variations in the age of the specimens, especially for individuals for which we have only an astragalus and no calcaneus or long bones associated, could account for some intraspecific variation. It is difficult to determine the age of individuals using only the astragalus, considering there are no epiphysis on this bone. Additionally, the high number of species, mostly with only one or two individuals, could also result in a greater diversity of morphological variations in our sample and thus lower the variance explained by the first axes.

Finally, it seems that Rhinocerotids are a relatively conserved group in terms of the morphology of their astragalus and calcaneus. This could mean that PC scores are more driven by small, individual-specific or species-specific variations than by large-scale variations such as variations linked to size or mass.

General influence of mass

The centroid size of both the astragalus and calcaneus has an effect on their shape (Figs. 6,9). Considering that the centroid size of the bone is linked to the mean mass of the species it belongs to (Appendix 5), especially on the Perissodactyla dataset, this means that mass has an influence on the shape of the bones in our sample of Perissodactyla. The percentage of variance explained by the centroid size or the mass, however, is lower than we originally expected. We could indeed expect mass to have a very strong influence on the shape of limb bones, explaining at least 50% of the total variance [START_REF] Hildebrand | Functional vertebrate morphology[END_REF][START_REF] Biewener | Mammalian terrestrial locomotion and size[END_REF][START_REF] Campione | A universal scaling relationship between body mass and proximal limb bone dimensions in quadrupedal terrestrial tetrapods[END_REF]. It seems that, especially when studying only rhinocerotids' astragali and calcanei, which do not vary much in terms of shape, size has no overwhelming influence on the shape of the bone but is instead one factor among others (e.g. possibly habitat, phylogeny and intraspecific variations). Another study on the limb long bones of extant rhinocerotids found a relatively low influence of centroid size on the shape of the bones, although higher than what is observed here (R-square between 10 and 18%; [START_REF] Mallet | Interspecific variation in the limb long bones among modern rhinoceroses-extent and drivers[END_REF]. It is therefore possible that long bones are more affected by size than astragalus and calcaneus.

When centroid size increases, both the astragalus and the calcaneus present an increase in the size of the articular facets. Moreover, in the analysis on Rhinocerotidae alone, the distal and medial facets of each bone for the articulation with the other are fused in specimens with a high centroid size. It could be suggested that wider facets result in a more intricate association between the bones, making the talocalcaneal complex more suited to dissipate compressive forces during limb loading, during plantarflexion or dorsiflexion of the foot (i.e. flexion or extension of the ankle). In the large astragali belonging to Perissodactyla, the trochlea is oriented proximally, directly facing the tibia and fibula, and has a lower angular extent (Fig. 12). One can assume that this orientation permits a complete unfolding of the crurotarsal joint, putting the foot in the exact same axis as the rest of the limb. This results in a general columnar posture for the limb, as is characteristic of graviportal animals. This columnar posture would help resist twisting, bending and compression forces, and reduce the possibility of dorsiflexion of the autopodium, reducing running speed [START_REF] Hildebrand | Analysis of vertebrate structure[END_REF]. It has already been said that Paraceratherium had columnar limbs [START_REF] Osborn | Baluchitherium grangeri, a giant hornless rhinoceros from Mongolia[END_REF][START_REF] Prothero | Rhinoceros giants: the paleobiology of indricotheres[END_REF], and our results corroborate that statement, for the hind autopodium-zeugopodium at least. It can be assumed that the flatter trochlea observed in large astragali, associated with a proximal-distal compression of the bone, fulfils the same role of resistance to twisting and compression. A deeper trochlea would provide more stability for the crurotarsal joint [START_REF] Polly | Limbs in mammalian evolution[END_REF], but lead to thinner and therefore more fragile ridges of the trochlea, unable to resist the high forces expected on the ankle of a very heavy animal. This flattening is also observed in Brontotheriidae [START_REF] Osborn | The Titanotheres of ancient Wyoming[END_REF], Elephantidae [START_REF] Csuti | The elephant's foot: prevention and care of foot conditions in captive Asian and African elephants[END_REF], and sauropod dinosaurs [START_REF] Bonnan | Pes anatomy in sauropod dinosaurs: implications for functional morphology, evolution, and phylogeny[END_REF], which supports our hypothesis.

For large-sized calcanei belonging to Perissodactyla, beyond the increasing size of the articular surfaces, the main characteristic is that the tuber calcanei is more robust, thicker both medio-laterally and antero-posteriorly, and shorter proximo-distally compared to the total length of the calcaneus (Fig. 15). However, this is only clearly observable when studying our Perissodactyla dataset. The tuber calcanei is a lever arm for the plantarflexion of the foot; with two of the muscles inserting on it being responsible for plantarflexion: the gastrocnemius and the soleus [START_REF] Beddard | On the anatomy of Rhinoceros sumatrensis[END_REF][START_REF] White | Chapter 16 -Foot: tarsals, metatarsals, & phalanges[END_REF]. A more robust and shorter tuber would presumably lead to a lower mechanical advantage, requiring a weaker pull from the muscles, which would be easier to resist and reduce bending stress. It could also be a consequence of the proximo-distal shortening of the foot generally observed in heavier species of our sample (Rhinocerotidae, Chalicotheriidae).An animal with an elongated foot would need a longer tuber calcaneus to keep the mechanical advantage constant; conversely, an animal with a short foot would not need a very long lever arm, presuming the mechanical advantage is indeed constant [START_REF] Biewener | Mammalian terrestrial locomotion and size[END_REF]. The correlation of the mean mass of the species on allometry-free shapes of the calcaneus corroborates this result. For two calcanei of the same size but belonging to species of different masses, the one belonging to the heaviest species will have a more robust tuber and wider articular facets (Fig. 15B). This is observable in our analysis in C. simum and R. unicornis: C. simum is heavier and has a slightly more robust calcaneus, but on average, the centroid size of its calcaneus is smaller than that of R. unicornis (Figs. 5,6).

This increased robustness of the body of the calcaneus is found in other mammal families of high body mass (i.e. more than about two tons), such as Elephantidae [START_REF] Chen | On the hindfoot bones of Mammuthus trogontherii from Shanshenmiaozui in Nihewan Basin, China[END_REF], and also in fossils such as Pyrotherium [START_REF] Shockey | Pyrotherium macfaddeni, sp nov (Late Oligocene, Bolivia) and the pedal morphology of pyrotheres[END_REF]. Interestingly, this is not the case in Hippopotamidae which have a rather elongate calcaneus (see e.g. Fig. 6 and 8 in [START_REF] Fisher | Hind limb myology of the common hippopotamus, Hippopotamus amphibius (Artiodactyla: Hippopotamidae)[END_REF]). Hippopotamuses have a body plan close to Teleoceratina, with very short limbs, which also present an elongate calcaneus. Possibly the forces exerted on the calcaneus are less intense for animals with short legs; comparisons with e.g. Suidae and Amynodontidae could yield insights in this regard.

Particular cases linked to body plan and locomotion

As expected, some observed variations in bone shape appear to be linked to the diverse body plans and modes of locomotion of the taxa studied. For the astragalus, equids are characterized by the great depth of their trochlea, a common characteristic in cursorial mammals allowing a stabilization of the crurotarsal joint by restricting movement to a parasagittal plan [START_REF] Polly | Limbs in mammalian evolution[END_REF].

The trochlea is also moderately deep in Rhinocerotidae (except most Teleoceratina) and

Tapiridae, but not in our Chalicotheriinae, animals that most likely could not gallop [START_REF] Coombs | Large mammalian clawed herbivores: a comparative study[END_REF]. Teleoceratina specimens also possess very shallow trochleas. Considering their similarity in terms of body-plan with hippopotamuses, which cannot gallop [START_REF] Lewison | Family Hippopotamidae (Hippopotamuses)[END_REF], it is likely that they could not gallop either, and the shape of their astragalus is consistent with this. Paraceratherium and the Chalicotheriinae possess the flattest trochlea of all of our specimens, but still with clearly distinguishable ridges (Figs. 11,12), unlike elephants for instance (see e.g. [START_REF] Scarborough | Insular adaptations in the astragalus-calcaneus of Sicilian and Maltese dwarf elephants[END_REF]. Equids also display a greater angular extent of their trochlea, presumably allowing a greater flexion and extension of the ankle. The most cursorial species (i.e. equids and, to a lower extent, tapirs) possess, on their astragalus, curved facets for the navicular and a curved proximal facet for the calcaneus, whereas those facets are mostly flat in Paraceratherium and Chalicotherium (Fig. 11). Perhaps the curved facets help to lock the talocalcaneal and talonavicular joints and provide stability for the ankle. The flat facets of heavy species could help homogeneously dissipate the forces in the foot and facilitate the formation of robust ridges.

A particular shape variation linked to body proportions is the proximo-distal compression of the astragalus across most of our Teleoceratina (Fig. 5). Diaceratherium is the only Teleoceratina from our sample presenting an astragalus similar to that of other rhinocerotids in this regard. Teleoceratina had extremely short, columnar limbs, like modern hippopotamuses.

This compression of the astragalus could be linked to the general shortening of the limbs, each segment being proximo-distally shortened, including the basipodium. Interestingly, the astragalus is not compressed in Diaceratherium intermedium, a Teleoceratina that is Chilotherium is characterized by short legs for its members [START_REF] Geraads | Rhinocerotidae (Mammalia) from the late Miocene of Bulgaria[END_REF], but no studies have been done specifically on D. intermedium. If this species was indeed shortlegged, the shortening of the limbs would predate the flattening of the astragalus in our sample.

The compression of the astragalus does not seem to be dependent on the size of the animal in our Teleoceratina. This condition is observed in both small (e.g. Prosantorhinus; <800kg) and large (Brachypotherium; >2000kg) Teleoceratina [START_REF] Cerdeño | Diversity and evolutionary trends of the family Rhinocerotidae (Perissodactyla)[END_REF][START_REF] Becker | Paléoécologie et paléoclimats de la Molasse du Jura (Oligo-Miocène) : apport des Rhinocerotoidea (Mammalia) et des minéraux argileux[END_REF]. It is worth noting that Paraceratherium presents the same flattening of the astragalus as our Teleoceratina.

Both groups are indeed very close regarding the morphology of their astragalus (Figs. 10,11).

Paraceratherium is, however, very different from Teleoceratina in that it is very high-legged.

It seems that different constraints, i.e. the very high mass of Paraceratherium and the short legs and lighter mass of Teleoceratina, can produce a similar result in terms of morphology. A study incorporating Amynodontidae (rhinocerotoids with some members, like Metamynodon, being short-legged like Teleoceratina, [START_REF] Wall | The phylogenetic history and adaptive radiation of Amynodontidae[END_REF], could also yield more insights on the matter.

Teleoceratina astragali also differ from those of other rhinocerotids by the distal elongation of their proximal facet with the calcaneus. The facet almost reaches the distal side of the bone, whereas it reaches only halfway in other rhinocerotids (Fig. 5). The facet might need to remain relatively long in order to maintain cohesion between the astragalus and the calcaneus. Thus, when the bone is proximally reduced, the facet retains the same length and occupies relatively more space on the posterior face. Teleoceras being an extremely variable genus in terms of bone morphology [START_REF] Harrison | Extreme carpal variability in Teleoceras (Rhinocerotidae, Mammalia)[END_REF]), a study with more individuals could yield insights on more subtle shape variations.

Chalicotheriinae also present a proximo-distally compressed astragalus. They differ from Teleoceratina and Paraceratherium in that their trochlea is oriented more anteriorly, and has a greater angular extent. This seems logical when looking at the angle of the crurotarsal articulation: the angle is clearly superior to 90°, almost reaching a flat angle, in Teleoceras and Paraceratherium, giving the limb a columnar posture. It is however approximatively equal to 90° in Chalicotherium, whose hindlimb is much more crouched (see e.g. Fig. 7B in [START_REF] Coombs | Large mammalian clawed herbivores: a comparative study[END_REF]. The extremely flattened astragalus of the Chalicotheriinae is not found in Moropus, which shows an astragalus closer to a rhinocerotid. The trochlea in particular is deeper in Moropus, whereas it is shallow in Chalicotherium and Anisodon. The extreme proximal-distal compression of Chalicotheriinae astragali could be a consequence of the reduction in the length of the hind limb, with each part of the limb being reduced, just as in Teleoceratina. This shortening could also be linked to the heavier mass carried by shorter hindlimbs, whereas body mass would be more evenly spread on fore and hindlimbs if they were of equal length. It could also be a consequence of their posture. Chalicotheriinae are indeed described as bipedal browsers. It is assumed that they could adopt an erected posture on their hindlimbs and use their forelimbs to grasp at branches and twigs [START_REF] Zapfe | Chalicotherium grande (BLAINV.) aus der miozänen Spaltenfüllung von Neudorf an der March (Děvinská Nová Ves)[END_REF][START_REF] Coombs | Large mammalian clawed herbivores: a comparative study[END_REF][START_REF] Schulz-Kornas | Feeding ecology of the Chalicotheriidae (Mammalia, Perissodactyla, Ancylopoda). Results from dental micro-and mesowear analyses[END_REF]. Most of their weight would therefore be supported by the hindlimbs, which would be in accordance with a stronger, flatter astragalus as observed in Paraceratherium. Moropus, and presumably the Schizotheriinae in general, have hindlimbs and forelimbs of approximately the same length, and were postulated to use bipedal browsing less frequently [START_REF] Coombs | Chalicotheres (Perissodactyla) as large terrestrial mammals[END_REF]1989).

This would reduce the advantage of a flatter astragalus. Further studies are needed to confirm or refute these hypotheses, with more individuals belonging to more genera, especially for Schizotheriinae (e.g. Ancylotherium or Metaschizotherium).

Another bone presenting a shape much different than what would be expected if mass were the single driving factor is Paraceratherium's calcaneus. Paraceratherium is by far the heaviest species of our sample, almost twice as heavy as Elasmotherium (Table 1). However, its calcaneus is elongated when compared to Elasmotherium's, which has the most robust calcaneus (Figs. 15,16). This could be a consequence of its general body plan: Paraceratherium had higher legs than all the rhinocerotids. One might thus suppose that longer legs, and thus, a longer autopodium like observed in Paraceratherium [START_REF] Prothero | Rhinoceros giants: the paleobiology of indricotheres[END_REF], lead to an elongation of the tuber calcanei to keep the mechanical advantage of the lever system of the foot constant. have suggested it could at least attain a trot. The fact that their astragalus retains clear ridges and that their calcaneus is quite elongated, characteristics reminiscent of cursorial animals [START_REF] Polly | Limbs in mammalian evolution[END_REF][START_REF] Bassarova | The calcaneum-on the heels of marsupial locomotion[END_REF], is consistent with this suggestion. Elephants possess a completely flat trochlea and a short calcaneus, and are unable of trotting or galloping.

Conclusion

Overall, it seems that mass has an influence on the shape of the astragalus and calcaneus in Rhinocerotidae and in our sample of Perissodactyla. However, that influence is lower than we initially thought, especially among Rhinocerotidae alone. This seems to indicate that Rhinocerotidae is a relatively conserved group in terms of the morphology of those bones, and that other factors, such as the phylogeny or intraspecific variations, have more influence. An ecomorphological study could help determine if habitat could have a role, but would require reliable habitat assignations for the fossil species. Nonetheless, bones belonging to heavier Rhinocerotidae present larger articular facets, presumably to help better dissipate the larger forces involved in the locomotion of heavier animals. The calcaneus is also more robust. In our sample of Perissodactyla, a stronger influence of mass is noted, with again heavier facets and stronger bones overall. We observe a flattened trochlea of the astragalus that would limit the risk of breaking, as compared to lighter animals which have a deeper trochlea with thin ridges for the stability of the crurotarsal joint. Although these features can thus be explained morphofunctionally, the phylogenetic signal is significant and could also explain variations between the families. A larger study encompassing large and small species of Perissodactyla will be necessary to determine more specifically what drives the shape of these bones in this order. Moreover, comparison between rhinocerotids and other perissodactyls reveal that body plan has a clear influence on the shape of the bones. Short-legged Teleoceratina display a flattened astragalus and an elongate calcaneus. Chalicotheres with short hindlimbs also display a flattened astragalus compared to chalicotheres with hindlimbs as long as their forelimbs, perhaps linked to the increased mass supported by the hindlimbs. Finally, Paraceratherium, which is extremely heavy and relatively high-legged compared to other rhinocerotoids, displays a flat astragalus as expected but a relatively elongate tuber calcanei, perhaps linked to either its elongate metapodials or to its phylogenetic history.

  Tapiridae, and eight astragali and six calcanei of Equidae. Three astragali and one calcaneus belong to Paraceratherium (Paraceratheriidae), and eight astragali and six calcanei belong to Chalicotheriidae. To our knowledge, all the bones belonged to adult specimens. A description of the bones including the nomenclature used for the main anatomical features of the bones is provided in Appendix 2. Specimens come from the collections of the Muséum National d'Histoire Naturelle (MNHN, Paris, France), the Muséum d'Histoire Naturelle de Toulouse (MHNT, Toulouse, France), the Claude Bernard University (UCBL, Lyon, France), the Natural History Museum (NHM, London, United Kingdom), the Powell-Cotton Museum (BICPC, Birchington-on-Sea, United Kingdom), the Naturhistorisches Museum Wien (NMW, Vienna, Austria), the Zoologische Staatssammlung München (ZSM, Munich, Germany), and the Bayerische Staatssammlung für Paläontologie und Historiche Geologie (BSPHM, Munich, Germany; Appendix 1).

Fig. 1 .

 1 Fig. 1. Composite phylogeny including the species sampled, modified from[START_REF] Antoine | Phylogénie et évolution des Elasmotheriina (Mammalia, Rhinocerotidae)[END_REF][START_REF] Antoine | A revision of Aceratherium blanfordi Lydekker, 1884 (Mammalia: Rhinocerotidae) from the Early Miocene of Pakistan: postcranials as a key[END_REF], Antoine et al. in press, Holbrook and Lapergola, 2011[START_REF] Piras | Functional and phylogenetic constraints in Rhinocerotinae craniodental morphology[END_REF][START_REF] Steiner | Molecular phylogeny and evolution of the Perissodactyla: phylogeny of the perissodactyls[END_REF]. Occurrence dates were estimated using[START_REF] Antoine | Aegyrcitherium beonensis, nouvel élasmothère (Mammalia, Perissodactyla) du gisement miocène (MN 4b) de Montréal-du-Gers (Gers, France). Position phylogénétique au sein des Elasmotheriini[END_REF][START_REF] Piras | Functional and phylogenetic constraints in Rhinocerotinae craniodental morphology[END_REF][START_REF] Antoine | A revision of Aceratherium blanfordi Lydekker, 1884 (Mammalia: Rhinocerotidae) from the Early Miocene of Pakistan: postcranials as a key[END_REF][START_REF] Geraads | A new rhinoceros, Victoriaceros kenyensis gen. et sp. nov., and other Perissodactyla from the middle Miocene of Maboko, Kenya[END_REF][START_REF] Guérin | Anisodon grande (Perissodactyla, Chalicotheriidae) de Sansan[END_REF][START_REF] Prothero | Rhinoceros giants: the paleobiology of indricotheres[END_REF], as well as data recorded on http://fossilworks.org. 1: Elasmotheriinae; 2: Rhinocerotinae; 3: Aceratheriina; 4: Teleoceratina; 5: Rhinocerotina; 6: Schizotheriinae; 7: Chalicotheriinae. Species in bold are living species. The relationships between the five extant rhino species are still debated, specifically the position of Dicerorhinus sumatrensis and its close fossil relatives is uncertain, placed either as sister taxa to Ceratotherium and Diceros, as sister taxa to Rhinoceros, or as sister taxa to a group formed of Rhinoceros, Ceratotherium and Diceros[START_REF] Willerslev | Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution[END_REF][START_REF] Gaudry | Molecular phylogenetics of the rhinoceros clade and evolution of UCP1 transcriptional regulatory elements across the mammalian phylogeny[END_REF]. This phylogenetic uncertainty is here represented by a polytomy.

Fig. 2 .

 2 Fig. 2. Representation of the landmarks and the curves placed on the astragalus of Rhinoceros unicornis MNHN-ZM-AC-1960-59. A: anterior, B: medial, C: posterior and D: distal views. White dots denote the 22 anatomical landmarks, dotted black lines the 9 curves. Description of the landmarks and curves is provided in Appendix 3.

Fig. 3 .

 3 Fig. 3. Representation of the landmarks and curves placed on the calcaneus of Ceratotherium simum MNHN-ZM-MO-2005-297. A: medial, B: anterior, C: latero-anterior and D: distal views. White dots denote the 13 anatomical landmarks, dotted black lines the 4 curves. Description of the landmarks and curves is provided in Appendix 3.

Fig. 4 .

 4 Fig. 4. Neighbour-Joining tree generated from a matrix of the Euclidian distance between every specimen, on the astragali of Rhinocerotidae. Legend as in Fig. 5.

Fig. 5 .

 5 Fig. 5. Results of the PCA performed on the astragalus of Rhinocerotidae. Left: repartition of the Rhinocerotidae astragali studied across PC1, PC2 and PC4. Right: TPS deformation of a mean shape towards the maximal and minimal value of each axis. The view is first posterior then anterior. Red dots denote landmarks, blue dots denote curve semi-landmarks. Vector representations of the deformations are provided in Appendix 7A.

Fig. 6 .

 6 Fig. 6. A: Regression of the Common Allometric Component on the logarithm of the centroid size, with representations of the shapes corresponding to the theoretical maximum and minimum of allometry, on Rhinocerotidae astragali. B: Regression of the Common Allometric Component of allometry-free shapes, on the logarithm of the cubic root of the mean mass of the species, with representations of the shapes corresponding to the theoretical maximum and minimum of mass, on Rhinocerotidae astragali. Legend as in Fig. 5. Vector representations are available in Appendix 7B.

Fig. 7 .

 7 Fig. 7. Neighbour-Joining tree generated from a matrix of the Euclidian distance between every specimen, on the calcanei of Rhinocerotidae. Legend as in Fig. 8.

Fig. 8 .

 8 Fig. 8. Results of the PCA performed on the calcanei of Rhinocerotidae. Left: repartition of the Rhinocerotidae calcanei studied across the first two PCA axes. Right: TPS deformation of a mean shape towards the maximal and minimal values of each axis. The view is first anteromedial then postero-medial. Red dots denote landmarks, blue dots denote curve semilandmarks. Vector representations of the deformations are provided in Appendix 7C.

Fig. 9 .

 9 Fig. 9. A: Regression of the Common Allometric Component on the logarithm of the centroid size, with representations of the shapes corresponding to the theoretical maximum and minimum of allometry, on Rhinocerotidae calcanei. B: Regression of the Common Allometric Component of allometry-free shapes, on the logarithm of the cubic root of the mean mass of the species, with representations of the shapes corresponding to the theoretical maximum and minimum of mass, on Rhinocerotidae calcanei. Legend as in Fig. 8. Vector representations are available in Appendix 7D.

  ) shows a relative clustering of the families. However, there are exceptions. Moropus is not grouped with the other chalicotheres but is closer in morphology to the Rhinocerotidae, although it has a long branch, which indicates a particular morphology. Kalobatippus, a three-toed anchitheriine equid from the Oligocene, has a morphology closer to tapirs and rhinocerotids than to a modern one-toed equinine or hipparionine equid. The Teleoceratina are found relatively close to Paraceratherium, and to a lesser extent to the Chalicotheriinae, as compared to other rhinocerotids. Teleoceras itself is the closest to the Chalicotheriinae, sharing with them an extremely proximo-distally flattened astragalus. Families are not grouped together according to phylogenetic proximity.

Fig. 10 .

 10 Fig. 10. Neighbour-Joining tree generated from a matrix of the Euclidian distance between every specimen, on the astragali of Perissodactyla. Extant species are represented as dots and extinct species, as squares. Teleoceratina, Kalobatippus and Moropus belong to Rhinocerotidae, Equidae and Chalicotheriidae, respectively, but are highlighted with regard to their particular positions in the tree. Legend as in Fig. 11.

Fig. 11 .

 11 Fig. 11. Results of the PCA performed on the astragalus of Perissodactyla. Left: repartition of the Perissodactyla astragali studied across the first three PCA axes. Right: TPS deformation of a mean shape towards the maximal and minimal value of each axis. The view is first posterior then anterior. Red dots denote landmarks, blue dots denote curve semi-landmarks. Vector representations of the deformations are provided in Appendix 7E.

Fig. 12 .

 12 Fig. 12. A: Regression of the Common Allometric Component on the logarithm of the centroid size, with representations of the shapes corresponding to the theoretical maximum and minimum of allometry, on Perissodactyla astragali. B: Regression of the Common Allometric Component of allometry-free shapes, on the logarithm of the cubic root of the mean mass of the species, with representations of the shapes corresponding to the theoretical maximum and

  homogeneity of the families. The closest specimen to the Paraceratherium specimen is Teleoceras fossiger. Contrary to what was observed on the astragalus, Moropus is grouped with the others Chalicotheriidae, although it is not as close to them as they are to each other. Again, families are not grouped following their phylogenetic relationships.

Fig. 13 .Fig. 14 .

 1314 Fig. 13. Neighbour-Joining tree generated from a matrix of the Euclidian distance between every specimen, on the calcanei of Perissodactyla. Teleoceratina and Moropus belong to Rhinocerotidae and Chalicotheriidae, respectively, but are highlighted with regard to their particular positions in the tree. Legend as in Fig. 14.

Fig. 15 .

 15 Fig. 15. A: Regression of the Common Allometric Component on the logarithm of the centroid size, with representations of the shape corresponding to the theoretical maximum and minimum

  phylogenetically the sister-group to the other Teleoceratina from our sample (Figs. 1, 5). It is unclear if D. intermedium was short-legged like the others Teleoceratina. The species was placed in the genus Chilotherium for a long time, before being reassigned by Antoine et al. in press.

Antoine

  et al. (2004) have indeed observed similarities between the calcaneus of Paraceratherium and that of a Giraffa. However, Teleoceratina have very short legs and a rather elongated calcaneus as well compared to other rhinocerotids (Fig.6), and elephants have high legs but a very short calcaneus. Others individuals of Paraceratherium and Elasmotherium are needed to confirm these results, as well as smaller members of the Paraceratheriidae family (e.g. Pappaceras and Juxia;[START_REF] Wang | Earliest known unequivocal rhinocerotoid sheds new light on the origin of giant rhinos and phylogeny of early rhinocerotoids[END_REF]. A study including other families of Rhinocerotoidea, like the small cursorial Hyracondontidae and the short-legged Amynodontidae could provide a better understanding of the question. Ultimately, comparing heavy and stocky mammals, such as Mixotoxodon (Notoungulata, Meridiungulata) or Hippopotamus amphibus (Hippopotamidae), with heavy and slender mammals like Titanotylopus (Camelidae) or Giraffa camelopardalis (Giraffidae) could also help understand the adaptations in the basipodium of Paraceratherium. However, in these extremely disparate taxa, one must be careful for a phylogenetical signal which could mask the changes of shape linked to mass. It is unclear what gait Paraceratherium was capable of adopting besides walking.[START_REF] Paul | Forelimb posture in neoceratopsian dinosaurs: implications for gait and locomotion[END_REF]

Fig. 16 .

 16 Fig. 16. Anterior views of the calcanei of Paraceratherium bugtiense NHM-PAL-PV-M-100418 (Paraceratheriidae) and Elasmotherium sibiricum NHM-PAL-PV-M-12429 (Rhinocerotidae).

Table 1 .

 1 List of the masses used for the species studied here. †: fossil taxon. *: two separate mass estimates were used for Equus caballus, one for average-sized horses and one for draught horses, given the wide morphological differences between the two.

		Iranotherium morgani †	Missing data
	Paraceratheriidae † Paraceratherium bugtiense †	7400	Fortelius & Kappelman, 1993
	Tapiridae	Tapirus pinchaque	150-200	Medici, 2011
	Tapiridae	Tapirus terrestris	220	Medici, 2011
	Tapiridae	Tapirus indicus	280-400	Medici, 2011
	Chalicotheriidae †	Chalicotherium sp. †	924	Costeur, 2004
	Chalicotheriidae †	Anisodon grande †	1500	Guérin, 2012
	Chalicotheriidae †	Moropus sp. †	1179	Damuth, 1990
	Equidae	Equus zebra	240-380	Rubenstein, 2011
	Equidae	Equus quagga	175-320	Rubenstein, 2011
	Equidae	Equus przewalski	200-300	Rubenstein, 2011
	Equidae	Equus caballus	380-600	Bongianni, 1988

Acknowledgements

This work acknowledges financial support from ERC-2016-STG GRAVIBONE allocated to

A. Houssaye. We acknowledge the very helpful comments of the three reviewers, P.-O. Antoine of the University of Montpellier, J. MacLaren of the University of Liège and M. Mihlbachler of the New York Institute of Technology. We wish to thank the curators of all the collections where we digitized specimens: J. Lesur, C. Bens, A. Verguin & G. Billet (Muséum National d'Histoire Naturelle, Paris, France), E. Robert (Université Claude Bernard, Lyon, France), Y. Laurent (Muséum d'Histoire Naturelle de Toulouse, Toulouse, France), C. West, R. Jennings & M. Cobb (Powell-Cotton Museum, Birchington-on-Sea, United Kingdom), P. Brewer, R. Pappa & R. Portela Miguez (Natural History Museum, London, United Kingdom), A.H. van Heteren (Zoologische Staatssammlung München, Munich, Germany), G. Rößner (Bayerische Staatssammlung fur Palaöntologie und historiche Geologie, Munich, Germany), and F. Zachos, A. Bibl & U. Göhlich (Naturhistorisches Museum, Vienna, Austria).

Supporting Information

Additional Supporting Information may be found in the online version of this article at the publisher's web-site.

Appendix 1. List of all the specimens studied.