#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Were ancient foxes far more carnivorous than recent ones?—Carnassial morphological evidence


Authors: Elwira Szuma aff001;  Mietje Germonpré aff002
Authors place of work: Mammal Research Institute Polish Academy of Sciences, Białowieża, Poland aff001;  Operational Direction “Earth and History of Life”, Royal Belgian Institute of Natural Sciences, Brussels, Belgium aff002
Published in the journal: PLoS ONE 15(1)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0227001

Summary

Crown shape variation of the first lower molar in the arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) was analyzed using five groups of morphotypes. Carnassial morphologies were compared between the species and between spatially and temporally distant populations: one Late Pleistocene (n = 45) and seven modern populations of the arctic fox (n = 259), and one Late Pleistocene (n = 35) and eight modern populations of the red fox (n = 606). The dentition of Holocene red foxes had larger morphotype variability than that of arctic foxes. The lower carnassials of the red fox kept have some primitive characters (additional cusps and stylids, complex shape of transverse cristid), whereas the first lower molars of the arctic fox have undergone crown shape simplification, with the occlusal part of the tooth undergoing a more pronounced adaptation to a more carnivorous diet. From the Late Pleistocene of Belgium to the present days, the arctic fox’s crown shape has been simplified and some primitive characters have disappeared. In the red fox chronological changes in the morphology of the lower carnassials were not clearly identified. The phyletic tree based on morphotype carnassial characteristics indicated the distinctiveness of both foxes: in the arctic fox line, the ancient population from Belgium and recent Greenland made separate branches, whereas in the red foxes the ancient population from Belgium was most similar to modern red foxes from Belgium and Italy.

Keywords:

Mammals – Teeth – Foxes – Dentition – Molars – Paleobiology – Pleistocene epoch – Belgium

Introduction

The carnassials–the upper fourth premolars (P4) and the opposite first lower molars (M1) are specialized and very important in the dental systems of carnivorous mammals. This pair of teeth occupies a central position in the tooth-row, which gives them advantage by maximizing bite force and gape. The main function of the carnassials is slicing and grinding of food; thus, the shape of the occlusal crown surface is intimately associated with dietary preferences in carnivores.

Contemporary studies of the types of food to which teeth and skulls are adapted use both traditional and advanced approaches. Traditional analyses of tooth function are based on measuring the size, size proportions, and angles of teeth [1, 2]. Recently, two-dimensional geometric morphometric [35] or three-dimensional modeling [6, 7] have often been used for analyzing crown complexity. The established approach relies on defining tooth shape into morphotype units [810]. This approach requires thorough knowledge of the dental morphology of the target species. Up till now morphotype analyses have been carried out during research of the dental characteristics of modern mammal species such as the red fox Vulpes vulpes [1115]; arctic fox Vulpes lagopus [1416]; corsac Vulpes corsac [14, 15]; raccoon dog Nyctereutes procyonoides [17]; European pine marten Martes martes, stone marten Martes foina and sable Martes zibellina [18, 19]; Japanese marten Martes melampus, American pine marten Martes americana, yellow-throated marten Martes flavigula, and fisher Martes pennanti [18]; badger Meles meles [20, 21]; brown bear Ursus arctos and arctic bear Ursus maritimus [22]; arctic ground squirrel Citellus parryi [23]; hares Lepus sp. [2427]; and some vole species Microtinae [10, 2832]. Such detailed analyses of the tooth morphologies of recent mammals are useful for studying fossil dental remains. Knowledge of tooth structure morphologies and evolutionary, ecological, behavioral and physiological constraints of the dental patterns of modern species enables us to picture the habitats and adaptations of extinct mammal species of past epochs.

Morphotype studies of recent materials are useful, especially as most identifications and classifications of extinct mammals are based on fossil remains such as molars and premolars, and pieces of mandibles with teeth [15, 3336]. Morphotype fractions like another measure–orientation patch count (OPC) can be used to quantify variation in dental complexity for many groups of mammals [6, 12]. Variation in the morphological dental patterns of mammals can be analyzed for various levels of systematic and ecological units. In each mammal species, subspecies, and population, specific dental patterns are a result of the multilateral effects of different factors, such as evolutionary history, food habits, climate, intra and inter-specific competition, and population density. Some shifts in tooth structure complexity between mammal taxa, subspecies, and geographically or chronologically separated populations of the same species are direct phenotypic responses to changes in the natural environment [13, 37]. Analyses of tooth occlusal surfaces (the complicated patterns of their crests and cusps) in microchiroptera, insectivore, carnivore, primate, and marsupial species indicate that more conspicuous and stronger occlusal tooth crown structures are related to harder foods; whereas, lower cusps and more complicated enamel crest patterns are related to consuming much softer-objects [38, 39]. Over the last century, some shifts have been registered in the morphotype frequencies of the dental pattern of the red fox in Poland [37]. The observations reflect the red fox’s broadening food niche and more optimal utilization of anthropogenic habitats. Subsequent studies on European red foxes have confirmed these observations. Over the last 50 years skull and body size in recent red foxes have grown in different populations [40, 41]. Both the temporal shift in dental morphology to a more generalist pattern and the increase in skull and body size of the red fox are a result of higher food availability and clear commensalism with man in the last decades of the last century.

Morphotype analysis of tooth shape in mammals can directly show dietary preferences and also reliably characterize their natural habitat. Morphological dental patterns constructed at the population, subspecies or species levels offer the opportunity to study phyletic or phylogeographic relationships and distances. Szuma [13], based on morphotype dental patterns of the red fox, showed variation between populations in carnivore specialization level and indicated the most probable place of diversification of the V. vulpes line: the area between Asia Minor and the Middle East. Morphotype characterization of the red fox throughout its whole native range uncovered the phylogeographic relations of the species. Analogical analysis was also carried out for the arctic fox across almost the whole of its Arctic range [16].

The aim of this research was to characterize the morphotype variation of the lower carnassials in the red and arctic foxes during the Late Pleistocene of Belgium. The morphotype patterns of the first lower molars in ancient samples were compared with morphological patterns found previously in extant populations of both species [13, 16]. Up till now there have been no studies on the phylogenetic relationships between Late Pleistocene arctic foxes in Europe. We expected to uncover reliable phyletic relationships between recent and fossil populations based on the crown shape of the lower carnassials. In our view the dental carnivorous specialization of the arctic and red foxes has weakened over the epochs following the Pleistocene. Over the last century the dietary opportunism of the European red fox has increased. Moreover, the dental pattern of Polish red foxes has been found to have become less carnivorous [37]. Based on our results, we discuss the evolutionary changes in the dental patterns of both foxes from the Pleistocene to present.

Material and methods

In this study we used fossil material of the red (V. vulpes Linnaeus, 1758) and arctic foxes (V. lagopus Linnaeus, 1758) from the Late Pleistocene of Belgium and recent material of the two species from geographically distinct populations. The fossil sample of the arctic fox contained 45 lower carnassials from four Palaeolithic sites (Trou de Chaleux, Goyet, Trou des Nutons, Trou Magritte; S1 Table). The fossil material of the red fox consisted of 35 isolated first lower molars from four Upper Palaeolithic sites (Trou de Chaleux, Goyet, Trou du Frontal, Trou des Nutons) in Belgium (S2 Table). Some carnassials could be from one specimen (the same fox): left and right carnassials; we studied isolated teeth so we could not assess specimen number. We analyzed left and right teeth together, as all specimens were equally valuable for drawing good conclusions. Morphological dental characters generally indicated bilateral symmetry; however there were also many cases of unilateral expression of morphotype characters. Determining how the collective analysis of left and right teeth impacted the study outcomes is difficult. However, the number of ancient samples of the red arctic fox seems to have been large enough to minimize the effect of the method.

All these fossil sites are in the Ardennes Massif and contain Pleistocene fauna. Additionally, a sample of 61 skulls of recent red foxes from Belgium was from the collections of the Royal Belgian Institute of Natural Sciences (RBINS).

The Trou de Chaleux cave is situated on the right bank of the River Lesse, a tributary of the Meuse River [42]. Inside this cave a major bone bed with a tremendous number of skeletal remains was discovered by Dupont [43]. Several AMS (Accelerator Mass Spectrometry) dates are available with calibrated ages ranging from 15 840 cal BP to 14 130 cal BP (S3 Table). The calibrations were calculated using the Oxcal 3.4 program (https://c14.arch.ox.ac.uk/oxcal/OxCal.html). Arctic fox remains represent 11.5% of the number of identified specimens (NISP) assigned to the Pleistocene fauna, whereas red fox remains account for 1.9% [44].

The Goyet cave forms part of a large karstic system on the right bank of the Samson valley, a tributary of the Meuse River. In this cave numerous Pleistocene mammal bones and human remains, and large quantities of Middle and Upper Palaeolithic artefacts have been found [4547]. The calibrated ages range from 44 990 cal BP to 14 460 cal BP for Horizon 1 and from 40 210 cal BP to 14 120 cal BP for Horizon 2 (S3 Table). Foxes constitute 3.8% of the NISP from Bed 1 [44, 48] and 3.7% of the NISP from Bed 2 ([44]; Germonpré unpublished data).

The Trou des Nutons cave is situated on the right bank of the Lesse River. The main bone bed yielded Magdalenian artefacts, and the Pleistocene mammal assemblage includes horse (Equus sp.) and reindeer (Rangiffer tarandus) bones [49]. Nevertheless, the AMS dates indicate that this assemblage consists of different components with calibrated ages ranging from 26 230 cal BP to 14 240 cal BP ([49, 50]; S3 Table). The bones of the arctic and red foxes constitute resp. 27.8% and 11.5% of the NISP of the Pleistocene assemblage (Germonpré unpublished data).

The Trou du Frontal cave is located on the right bank of the Lesse River near the Trou des Nutons cave. The mammal assemblage contains a mixture of remains dating from the Late Glacial and Postglacial. The calibrated ages of cut-marked horse bones range from 14 780 cal BP to 15 730 cal BP (S3 Table). Arctic fox remains form 7.1% of the NISP of the Pleistocene assemblage in this cave, whereas red fox remains form more, at 10.7% ([44]; Germonpré unpublished data).

The Trou Magritte cave is about 3 km upstream of the confluence of the Lesse and Meuse Rivers. The Palaeolithic assemblages have been assigned to the Aurignacian and Middle Palaeolithic [51]. The mammal assemblage is primarily composed of remains of horses, reindeer and woolly rhinoceros (Coelodonta antiquitatis), some of which bear anthropogenic traces while others show evidence of having been gnawed by cave hyenas [52, 53]. The remains of the foxes constitute 3.1% of the NISP of the Pleistocene assemblage [54]. An AMS date of a mammoth bone indicates that the cave was in use between 46 640 cal BP and 42 970 cal BP (S3 Table).

To follow evolutionary changes in the morphological structure of the lower carnassials in the arctic and red foxes, recent materials were incorporated into the analyses. The contemporary materials included eight sampled populations (545 recent M1) of the red fox from the Palearctic region and seven sampled populations (259 recent M1) of the arctic fox from across its whole arctic range (Fig 1). A list of samples, abbreviations of population names, systematic statuses, collection data and collection names can be found in Table 1. The modern fox samples had been used before in our earlier studies [13, 16, 53, 55]. The number of red foxes in each sample was higher than 50 specimens, whereas the number of arctic foxes in samples ranged from n = 11 (Greenland) to n = 58 (Taymyr). Shape analysis of M1 for each specimen was carried out according to morphotypes defined and illustrated previously for the red fox ([13]; Fig 2). Five groups of morphotypes (K, L, P, R, and S) were used to depict M1 shape variation in the foxes. Morphotypes within each group were identified by Arabian numbers, with higher numbers indicating more complex morphologies. For group P, Szuma [12] classified 5 morphotypes based on numerous variants of crown shape. Definitions of the dental morphotypes used to analyze dental polymorphism in the arctic fox are listed in S4 Table.

Geographic distribution of the recent and fossil samples of the arctic <i>Vulpes lagopus</i> and red foxes <i>Vulpes vulpes</i> used in the study.
Fig. 1. Geographic distribution of the recent and fossil samples of the arctic Vulpes lagopus and red foxes Vulpes vulpes used in the study.
Black rhomb–recent populations of the arctic fox measured by Szuma (2011), empty rhomb–fossil populations of the arctic fox measured by Szuma for the study, grey circle–recent populations of the red fox measured by Szuma (2011); empty circle–fossil population of the red fox measured by Szuma for the study. Information about the samples can be found in Table 1.
Graphic depiction of morphological variation of the first lower molar (M<sub>1</sub>) in the red and arctic foxes.
Fig. 2. Graphic depiction of morphological variation of the first lower molar (M1) in the red and arctic foxes.
Group K–M1 in lateral view, lingual side; group L–M1 in lateral view, buccal side; group P, R, S–distal part of M1 (talonid) in occlusal view. Definitions of the morphotypes, and morphotype variants are listed in S4 Table.
Tab. 1. List of samples of the red and arctic foxes used in the study with specimen numbers, collection times, and sources of the specimens.
List of samples of the red and arctic foxes used in the study with specimen numbers, collection times, and sources of the specimens.

For the analysis of fossil materials, both left and right carnassials were used. Previous analyses of frequency distributions of morphotypes in female and male Polish red foxes showed no significant differences between sexes [11]; thus, the sexes in recent samples of the arctic and red foxes were combined. This is very convenient for comparing the ancient and recent fox samples, and does not add error to the statistical analysis.

Frequency distributions of the five morphotype groups were analyzed in ten red and eight arctic fox populations. The significance of interpopulation variation in relation to the dental morphotypes in the foxes was determined with the chi-square tests χ2. The tooth dimensions: length of M1 (LM1), width of M1 (WM1) and carnivorous index (LM1/WM1) in relation to the particular morphotypes were verified using MANOVA. One of the authors (ES) has analyzed the sizes and shapes of the lower carnassials of foxes in various other studies [1113, 37, 5557].

Before we carried out the multidimensional analyses, we converted the percentage share of each morphotype using angular transformation in conformity with the equation y = arc sinx. Angular transformation is needed if the number n on which the individual percentages are based is very variable [58]. The transformations were carried out in MS Excel 97.

Based on the transformed morphotype frequencies of shape variation of M1, we carried out multidimensional scaling (MDS). An agglomeration method was used to determine levels of morphological similarity among recent and fossil populations of the arctic and red foxes. In the agglomeration analysis the city distance (Manhattan) and Ward’s method of linking were used. All statistical analyses were performed using Statistica.PL version 9.0.

Results

Morphotype frequency in the arctic fox

The distribution of group K morphotypes in fossil arctic foxes from Belgium (pop. 1) was very close to the recent fox population from the Yamal Peninsula (pop. 5, Fig 3)–they had a low share of the morphotype K1 (17% in pop. 1, and 24% in pop. 5), and a presence of the K3 morphotype. The lowest share of morphotype K1 (6%) was observed in arctic foxes from Yakutskaya oblast (pop. 8), whereas the highest was in Beringia (71%, pop. 2), and Greenland (100%, pop. 4). In group L, for the arctic fox we found morphotypes L1 and L2. The frequency of the morphotype L2 in the ancient arctic foxes (pop. 1) was 27%, whereas among the recent populations of the arctic fox the L2 share ranged from 14% in the Yamal Peninsula (pop. 5) to 0% in Beringia (pop. 2). In the ancient sample (pop. 1) we identified four group P morphotypes (P1–P4), and the cumulative share of morphotypes P2, P3, and P4 (55%) was larger than in recent populations of the species. Among the recent foxes only the Baffin Land population (pop. 3) had a large share of the more complex morphotypes (P3, P4)– 53%. In all other populations of the species the more complicated morphotypes of the transverse cristid of M1 (see Fig 1) had a share lower than 28% (Fig 3).

Share of morphotypes in fossil and recent populations of the arctic fox.
Fig. 3. Share of morphotypes in fossil and recent populations of the arctic fox.
Numerals in each part of the pie chart are the numbers of specimens with the morphotypes. Illustrations of the tooth variants and morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

The distal part of the talonid of the lower carnassials in the arctic fox showed a simplistic occlusal structure. In three populations a small hypocingulid (morphotype R2) appeared at the levels of 2% (pop. 2), 4% − Kozhevnikovo (pop. 7), and 3% (pop. 1) (Fig 3). We did not observe morphotype S2 anywhere across the whole range of the species. In recent and fossil populations of the arctic fox morphotype S1 occurred with a 100% frequency.

Morphotype frequency in the red fox

The percentage shares of the lower carnassial morphotypes relating to shape of M1 differ between the red and arctic foxes. In the red fox, morphotype K2 oscillated from 64% in Novosibirskaya oblast (pop. 15) to 96% in Iran (pop. 14; Fig 4). In contrast, the share of K1 ranged from 3% (pop. 14) to 34% (pop. 15). In the ancient red foxes (pop. 9), morphotype K1 reached 13%, K2 84%, and K3 just 3%. Group L was represented by all three morphotypes; however, L1 varied from 94% in Belgium (pop. 10) to 57% in Chukotka (pop. 17). In the ancient sample (pop. 9) morphotype L1 constituted 85%, L2 − 11%, and L3 − 4% (Fig 4). Shares of group P morphotypes in the studied samples were characterized by large levels of differentiation. Morphotype P1 oscillated from 0% in Italy (pop. 12) to 30% in pop. 15. The frequency of morphotype P3 in the red fox ranged from 27% (pop. 15) to 47% in Yakutskaya oblast (pop. 16). In red fox populations, the most frequent morphotype was P4 and its values oscillated from 37% (pop. 14) to 61% (pop. 12). The lower carnassials in the red fox were characterized by the presence of two morphotypes, R1 and R2. In population 14, individuals with morphotype R2 comprised 14% of the whole population, whereas in Kamchatka (pop. 18) and Chukotka (pop. 17) red foxes had morphotype R2 frequencies of 81% and 89%, respectively (Fig 4). Among the fossil red foxes (pop. 9) 50% of individuals possessed the morphotype R1 and the other 50% − R2. Morphotype S2 occurred with lower frequencies than S1 in all populations of the species. The highest levels of the morphotype S2 were found in red foxes from Chukotka (pop. 17; 28%), Kamchatka (pop. 18; 23%), ancient Belgium (pop. 9; 23%), and Italy (pop. 12; 21%), whereas the lowest levels of S2 were in Iran (pop. 14; 2%), Egypt (pop. 13; 3%), and Belgium (pop. 10; 6%).

Share of morphotypes in fossil and recent populations of the red fox.
Fig. 4. Share of morphotypes in fossil and recent populations of the red fox.
Numerals in each part of the pie chart are the numbers of specimens with the morphotypes. Illustrations of the tooth variants and morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

Morphotypes versus carnassial size

Group K

The arctic fox: The median length of the lower carnassials was highest for the K3 morphotype group, and was lower for K2 and K1; however, the differences were not statistically significant (p>0.05; Fig 5). The width of the carnassial crown was also highest in samples with morphotype K3 and lower in the other two (p = 0.005; Fig 5). The carnivorous index of arctic fox samples with morphotypes K1, K2 or K3 did not vary significantly (p>0.05; Fig 5). The MANOVA for the arctic fox was only carried out for a few samples because of low n in the samples and morphotype fractions. However, there were statistically significant differences in LM1/WM1 between three populations with various group K morphotypes (p = 0.007; Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various K group morphotypes in the arctic fox with and without population division.
Fig. 5. Comparison of median values for length, width, and size proportion of M1 between various K group morphotypes in the arctic fox with and without population division.
Illustrations of group K morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.
Tab. 2. Size variation of the first lower molar (M1) in the arctic and red foxes according to population division (POP), morphotype differentiation (morphotype groups: K, L, P, R, S), and interactions between both factors (POP*K; POP*L; POP*P; POP*R; POP*S).
Size variation of the first lower molar (M<sub>1</sub>) in the arctic and red foxes according to population division (POP), morphotype differentiation (morphotype groups: K, L, P, R, S), and interactions between both factors (POP*K; POP*L; POP*P; POP*R; POP*S).
Bold letters indicate statistical significances (p<0.05) of the relationships.

The red fox: The lower carnassials in red foxes with morphotype K1 had the highest median crown length (Fig 6). There was a statistically significant difference in the median values for LM1 between foxes with morphotypes K1, K2, and K3 (p = 0.001). Morphotype K3 was characterized by a moderate median carnassial crown length. Similarly, the median width of the lower carnassials was highest in group K1, moderate in K3, and lowest in K2 (Fig 6). With regards to the carnivorous index (LM1/WM1), red foxes with the various K morphotypes did not show any significant differences (p>0.05) in the central distribution measure (Fig 6). However, when taking into account the population division of the samples with the various K morphotypes, there were statistically significant differences in all the M1 parameters (Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group K morphotypes in the red fox with and without population division.
Fig. 6. Comparison of median values for length, width, and size proportion of M1 between various group K morphotypes in the red fox with and without population division.
Illustrations of group K morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

Group L

The arctic fox: Samples with morphotypes L1 and L2 had similar median values of WM1 and carnivorous index (LM1/WM1) (p>0.05; Fig 7). Medians for crown length of the lower carnassials were higher in morphotype L2 than in L1 (p<0.05). There were statistically significant differences (p<0.001) between some populations with the various L morphotypes, as calculated with the LM1/WM1 index. When considering the population division of samples with L morphotypes, there was statistically significant variation in WM1 and LM1/WM1 between groups of arctic foxes (Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group L morphotypes in the arctic fox with and without population division.
Fig. 7. Comparison of median values for length, width, and size proportion of M1 between various group L morphotypes in the arctic fox with and without population division.
Illustrations of group L morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

The red fox: Medians for LM1, WM1, and LM1/WM1 showed similar values in L1, L2, and L3 morphotype (p>0.05; Fig 8). However, when considering population and L group divisions there were statistically significant differences in the means of WM1 and LM1/WM1 between samples (p<0.05; Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group L morphotypes in the red fox with and without population division.
Fig. 8. Comparison of median values for length, width, and size proportion of M1 between various group L morphotypes in the red fox with and without population division.
Illustrations of group L morphotypes are shown in Fig 2 Information about the samples can be found in Table 1.

Group P

The arctic fox: Group P morphotypes were similar (p>0.05) in all metrical parameters and had similar carnivorous indices (Fig 9). When considering population and morphotype division (P1–P4), a significant cumulative effect of both factors on LM1/WM1 between sample variation (p<0.05; Table 2) was found.

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group P morphotypes in the arctic fox with and without population division.
Fig. 9. Comparison of median values for length, width, and size proportion of M1 between various group P morphotypes in the arctic fox with and without population division.
Illustrations of group P morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

The red fox: Group P morphotypes had similar (p>0.05) medians of metrical parameters and carnivorous indices (Fig 10). However, when considering differentiation between populations, significant variation (p<0.005) was observed. The MANOVA procedure confirmed there was statistically significant variation for all parameters of M1 between populations with various P morphotypes (Table 2). The cumulative effect of population and morphotype division translated into significant variation in the size and size proportion of M1 between samples.

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group P morphotypes in the red fox with and without population division.
Fig. 10. Comparison of median values for length, width, and size proportion of M1 between various group P morphotypes in the red fox with and without population division.
Illustrations of P group morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

Group R

The arctic fox: Median length and width were significantly larger in the arctic foxes with the R2 morphotype than those with R1 (Fig 11). Samples with morphotypes R1 and R2 had similar (p>0.05) carnivorous indices (Table 2). The cumulative effect of both factors is difficult to verify because the df is only 2 (Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various group R morphotypes in the red fox with and without population division.
Fig. 11. Comparison of median values for length, width, and size proportion of M1 between various group R morphotypes in the red fox with and without population division.
Illustrations of group R morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

The red fox: Specimens with morphotypes R1 and R2 had significant differences (p<0.001) in the median values of length and width of M1 (Fig 11). Only the median of carnivorous index was similar between the morphotypes (p>0.05; Fig 11). When considering population and R morphotype divisions, all metrical parameters of M1 showed highly significant variation (p<0.001) between samples. The cumulative effect of the population and R division was visible in WM1 (Table 2).

Group S

The arctic fox: The S1 morphotype occurred with 100% frequency.

The red fox: Specimens with morphotype S2 had higher median values of metrical characters than those with morphotype S1 (p<0.001; Fig 12); the LM1/WM1 index was also higher in those with morphotype S2 than those with morphotype S1, based on the level p = 0.05 (Fig 12). The MANOVA procedure found statistically significant (p<0.001) differentiation in LM1, WM1, and LM1/WM1 in population and morphotype separated red foxes. The cumulative effect of the population and S division was visible in WM1 (p<0.05; Table 2).

Comparison of median values for length, width, and size proportion of M<sub>1</sub> between various S group morphotypes in the red fox with and without population division.
Fig. 12. Comparison of median values for length, width, and size proportion of M1 between various S group morphotypes in the red fox with and without population division.
Illustrations of group S morphotypes are shown in Fig 2. Information about the samples can be found in Table 1.

Morphological distances between populations

In the arctic fox, the MDS result, based on the angularly transformed frequencies of M1 morphotypes, suggested that the ancient sample (pop. 1) was distant compared with modern samples of the species. The Beringia and Greenland populations of the arctic fox (pops 2, 4; Fig 13) were very different from both fossil and all modern samples.

Distribution of the populations of the arctic and red foxes in a two-dimensional plane based on the morphotype variations of the first lower molar.
Fig. 13. Distribution of the populations of the arctic and red foxes in a two-dimensional plane based on the morphotype variations of the first lower molar.
Black rhomb–recent populations of the arctic fox, empty rhomb–fossil populations of the arctic fox, grey circle–recent populations of the red fox, empty circle–fossil population of the red fox. Information about the samples can be found in Table 1.

The red fox had smaller morphological distances between fossil and modern foxes than the arctic fox. The fossil Belgian red foxes (pop. 9) were closest to recent foxes from Belgium (pop. 10), Italy (pop. 12), England (pop. 11), and Egypt (pop. 13), and most distant to populations from Kamchatka and Chukotka (pops 18, 17; Fig 13).

The cluster tree based on the angularly transformed frequencies of M1 morphotypes separated the red and arctic foxes very well (Fig 14).

Agglomeration of the recent and fossil populations of the arctic and red foxes constructed by the Ward method and Manhattan distances using the percentage shares of the morphotypes of the first lower molar.
Fig. 14. Agglomeration of the recent and fossil populations of the arctic and red foxes constructed by the Ward method and Manhattan distances using the percentage shares of the morphotypes of the first lower molar.

Discussion

Lower carnassial evolution in the arctic fox

The current knowledge on the evolutionary history of Vulpes suggests that the genus appeared ca. 7 Ma BP [59]; the oldest fox remains have been found in the late Miocene deposits of Africa, Chad–Vulpes riffautae [59]. The arctic fox probably evolved from an ancestral line of a cold adapted ancestor–an early Pliocene (3.60–5.08 Myr ago) fox, Vulpes qiuzhudingi, from the Tibetan Plateau [60, 61]. This early arctic fox from the Himalaya (holotype: IVPP V18923) was characterized by very sharp, hypercarnivorous teeth. In this specimen, the talonid of the first lower molar was narrow, with a tendency to be distally narrower, and its second lower premolar had a relatively narrow crown. The shape of the carnassials had a primitive character, exemplified by the presence of additional small conules (a postmetaconulid and entoconulid) in the concavity between the metaconid and entoconid–morphotype K3. The Pleistocene arctic foxes from Belgium had a low share of K1 morphotypes (17%; an absence of small conule/-s in the concavity between the metaconid and entoconid). The more complicated shapes K2, K3, and L2 had a larger share in the fossil arctic foxes than in the recent ones. Moreover, the occlusal surface of the talonid on M1 in the fossil arctic foxes from the Late Pleistocene of Belgium had more complicated enamel cristids in the central part of the talonid basin than the recent arctic foxes. The comparison of the morphotype characterization of the lower carnassials between the ancient arctic foxes from Belgium, the ancient holotype from the Tibetan Plateau, and individuals from geographically distant recent populations revealed that the arctic fox has undergone distinctive, parallel evolutionary changes in terms of tooth shape, simplification of crown shape and disappearance of some primitive characters, while also undergoing changes in the size and size proportions of the carnassials [62]).

The newly described, oldest arctic fox V. qiuzhudingi was bigger than both the late Pleistocene foxes from Belgium, and the recent arctic foxes from across their whole range. This new species of arctic fox displayed great predatory adaptations in its dental system because of its much smaller talonids and the presence of several primitive characters in the lower carnassials. The large gaps between the Early Pliocene arctic fox from the Tibetan Plateau and Late Pleistocene arctic fox from Belgium, and then between the latter and recent arctic fox could suggest intensive evolutionary processes have been taking place since the first occurrence of the arctic fox progenitor. The heavy hypercarnivorous specialization of the new species, V. qiuzhudingi, was linked with its adaptation to the very cold conditions present in the Himalayan region at that time, and the presence of sympatric cold-adapted mammals (large body size, long hair, snow sweeping structures)–a potential food base for foxes. A new hypothesis was recently formulated suggesting that the paleo-elevation of the Tibetan Plateau was the ‘starting line’ for most Quaternary Ice Age fauna [61]. The oldest representatives of the woolly rhino Coelodonta thibetana have been identified from the fossil beds of the Tibetan Plateau [62]. Deng et al. [61, 63] indicated that C. thibetana was the baseline form in the woolly rhino lineage history. Scavenging of the then Tibetan foxes V. qiuzhudingi on carcasses of big, cold-adapted forms of mammals required the relevant body and tooth construction properties. The crown length of the lower carnassials of the ancestral form of the arctic fox reached 16.6 mm, whereas among the recent arctic foxes the specimen with the largest length of this molar– 15.51 mm–was recorded in the Yamal Peninsula population, with the whole sample having a mean of 13.74 mm [55]. In the fossil arctic foxes from the Late Pleistocene of Belgium the largest length of M1 was 15.49 mm, whereas the mean of LM1 was 14.21 mm. In the evolutionary history of the arctic fox we can observe a progressive reduction in tooth size and an undermining of carnivorous specialization in the dental system. Considering, the great gap between the progenitor of the arctic fox line–V. qiuzhudingi, the next step of the V. lagopus from the Late Pleistocene, and then later the recent arctic fox–the magnitude of the changes in the molar’s structure and size fit the length of the time intervals between the studied evolutionary stages. We can observe a spectacular change in the molar’s morphology from the Pliocene to the Late Pleistocene (about 4 Myr BP) and then from the Late Pleistocene to the recent Holocene (about 42–12 kya BP).

Lower carnassial evolution in the red fox

The recent red fox is the largest of the extant representatives of the Vulpes genus and occupies the largest natural range with a very wide spectrum of habitats. The species has a very malleable genotype that has enabled it adapt to rapid changes in recent times. With great success the red fox has inhabited a new ecological niche–landscapes heavily transformed or occupied by man. The fox’s great plasticity means that the species is very expansive, leading to it occupying even the northernmost areas of the globe [64, 65]. The northern expansion of the red fox in Sweden over the last century is a consequence of the population decrease of the grey wolf and European lynx, and increasing ungulate densities. Thus the red fox has expanded into the native range of the arctic fox and competes with it effectively [66]. This evolutionary success of the red fox is most probably a result of its very flexible phenotype, e.g., inter alia, its tooth morphology allows it to use many different types of food; the red fox has 15 morphotypes of crown structure of the lower carnassial, whereas the arctic fox has just 12. The main indicator of the lower carnassials of the red fox is the larger significance of the talonid, i.e. the presence of an entostylid, and a more developed transverse cristid than the arctic fox. Among the worldwide population of the red fox, we have previously observed 34 variants in the form of the transverse cristid in the talonid basin [14], whereas in the arctic fox only 14 [16].

The ancient sample of the red fox from the Late Pleistocene of Belgium had smaller differences in the lower carnassial morphology compared to modern samples of the species, than those observed between the arctic fox samples. Only for group R and S morphotypes we found significant differentiation in the size of the first lower molar (p<0.05). We observed that the red fox specimens with additional cusps present in the most distal part of the lower carnassials (R2, S2 morphotypes) tend to have longer or wider teeth, e.g. the Late Pleistocene of Belgium, Chukotka, and Kamchatka red foxes. It seems that the presence of the additional distal cusps on M1 could have been affected by the threshold mechanism. For some characters to be expressed during ontogeny (i.e for some structures and their elements to occur), they need to reach or exceed an appropriate threshold value (size, concentration) during the early stages of development [21, 67]. Additional cusps on the anterior or distal parts of the tooth crown are specific to ancestral, primitive carnivorous mammals [68]. Over the course of red fox’s evolutionary history, the small cusps on the talonid have occurred at a relatively stable frequency.

The baseline and most primitive set of dental characters in the red fox have been found in the southern regions of Asia (India, Turkmenistan, Morocco, Amurskaya oblast), indicating that the V. vulpes line most probably radiated in Asia Minor, Middle Asia, and India [13]. For example, the frequency of occurrence of the entoconulid (K2, K3 morphotypes) is highest in recent red foxes from Morocco, Saudi Arabia, Iran, Turkmenistan, Pakistan, India, and Futsing province in China [13]. In a specimen of V. riffatuae from the Late Miocene of Chad, the entoconulid was not found due to heavy wearing of the occlusal surface of the M1 [59]. However, among other ancient foxes, e.g. in Vulpes beihaiensis, a well shaped entoconulid was confirmed [69], while in Vulpes praecorsac an entoconulid was found on one of the three first lower molars [68]. The presence of entoconulids in some specimens of ancient foxes (V. beihaiensis, V. praecorsac, V. vulpes: pop. 9) and comparable frequencies of additional cusps in the occlusal structures of M1 in ancient and recent red foxes indicated there was no substantial shift in the morphology of the tooth from the Pleistocene to the recent time. Some changes in tooth morphology toward hypocarnivory were observed in the Polish red foxes over seventy years of the 20th century [70]. Microevolutionary changes in the occlusal structure of the premolars (i.e. P4, P3, P4) and the incisor (I3) were found, whereas no relevant trends were detected in the other cheek teeth of the species [70]. The evolutionary history of the Carnivore group has revealed that less specialized forms (e.g. Miacids, recent Canids) retained a greater ability to react to environmental changes than heavily specialized ones (e.g. ancient Hyaenodontidae, Nimravidae). Over the course of evolution, the dental pattern of the red fox has oscillated around a certain norm, so that at any given point in the fox’s phylogeny, it could react effectively to new factors and environmental changes. Red fox teeth in modern populations are an example of a prompt evolutionary response of tooth morphology to dietary changes [70]. This morphological shift has confirmed that the Vulpini group has great abilities to adapt to fast-changing environments.

Phylogeny of the arctic and red foxes

The morphological patterns of the lower carnassials showed very clear separation between populations representing the arctic and red fox lineages. The percentage shares of particular crown shape morphotypes evidently separated the fox species, and their geographically and chronologically distant populations. In the phenogram we constructed, the Beringian and Greenland recent arctic foxes made a separate branch. The Siberian arctic foxes (Kozhevnikovo, Taymyr, Yamal, Yakutsk) and the population of Baffin Land were the most similar morphologically, whereas the population of ancient arctic foxes from Belgium was more different to recent samples. Another scenario was observed between red fox populations, where the ancient red foxes from the Late Pleistocene of Belgium were morphologically very close to the recent ones from Italy and Belgium. The most morphologically distant to other populations were the recent red foxes from Kamchatka and Chukotka.

The phylogenetic tree based on the morphological characteristics of M1 for the arctic fox to some extent agreed with a tree based on metric parameters [62]. Both studies indicate fossil arctic foxes to be distant to recent populations of the species. This evolutionary line of the arctic fox occupied Western Europe during cold periods of the Pleistocene, and after the retreat of the last glaciers arctic foxes from Belgium did not go north. Dalén et al. [71] based on ancient DNA analyses suggested that European fossil populations of the arctic fox most probably died out with the end of the Pleistocene epoch. The morphological analysis of the lower carnassials indicated that the arctic foxes from the Late Pleistocene of Belgium shared the same fate as other foxes populating glacier free-areas in western and southern Europe during the Late Pleistocene.

A comparison of the cladograms for the red fox populations based on morphotype variation and metric characteristics of the lower carnassials showed that they uniformly indicated that the fossil red foxes from the Pleistocene of Belgium are various to recent foxes. In the size of M1 and proportion index (LM1/WM1) ancient Belgian foxes were similar to recent foxes from Siberia, i.e. Novosibirskaya oblast, Yakutskaya oblast, Chukotka, Kamchatka [62]. This was clear because in carnivores very high correlations between length of the M1 and body mass (r25 = 0.97, p<0.001), according to Flower [72] has been reported before. The red fox follows Bergmann’s rule, so Siberian foxes are much bigger than red foxes from the southern areas of their range [73].

Based on the morphotype patterns of the lower carnassials we revealed that ancient Belgian red foxes were most similar to recent Belgian and Italy populations of the species. These three red fox populations were most similar to each other morphologically and the closest geographically. By comparing the relationships in our phenogram with the phylogenetic tree of the relations between populations across the whole range of the red fox based on dental polymorphism [16], we can conclude that the western and southern European red foxes make up a distinct line from the central and east European populations. This observation suggests the Late Pleistocene red foxes of Belgium and recent Belgian and Italian populations share a common gene pool, and during the Last Glacial Maximum red foxes from western Europe survived in a southern–Italian refugium. The recent Polish population of red foxes was most probably populated by red foxes from some eastern refugium.

Conclusions

Morphotype analysis of the lower carnassials in ancient and extant populations of the arctic and red fox revealed distinctive evolutionary trends in the two species. In the arctic fox the progressive changes in dental form tended towards hypercarnivory, whereas in the red fox the most primitive characters have been retained from the Pleistocene to recent times. Such evolutionary tendencies for both foxes are a result of their varying life strategies and evolutionary trajectories in the Late Pleistocene, and Holocene. The arctic fox occupies the most extreme habitats of the northern hemisphere, and over the last decades its range has been declining in response to pressure from the red fox.

Supporting information

S1 Table [doc]
List of the fossil first lower molars of the arctic fox () from RBINS, Brussels, used in the study with information about site excavation and left/right position of the molar in the tooth-row.

S2 Table [doc]
List of the fossil first lower molars of the red fox () from RBINS, Brussels, used in the study with information about site excavation and left/right position of the molar in the tooth-row.

S3 Table [doc]
AMS dates and their calibration () of mammals present at the Belgian sites analyzed in this study (anthrop.: Anthropological).

S4 Table [doc]
List of dental morphotypes used to assess of dental polymorphism in the red fox () and the arctic fox ().


Zdroje

1. Binder WJ, Van Valkenburgh B. Development of bite strength and feeding behavior in juvenile spotted hyenas (Crocuta crocuta). J Zool Lond. 2000;252:273–283.

2. Evans AR, Sanson GD. The tooth of perfection: functional and spatial constraints on mammalian tooth shape. Biol J Linn Soc. 2003;78:173–191.

3. Meloro C, Raia P, Piras P, Barbera C, O’Higgins P. The shape of the mandibular corpus in large fissiped carnivores: allometry, function and phylogeny. Zool J Linn Soc. 2008;154:832–845.

4. Polly D. Development and phenotypic correlations: the evolution of tooth shape in Sorex araneus. Evol Dev. 2005;7:29–41. doi: 10.1111/j.1525-142X.2005.05004.x 15642087

5. White J. Geometric morphometric investigation of molar shape diversity in modern lemurs and lorises. Anat Rec. 2009;292:701–719.

6. Smits PD, Evans AR. Functional constraints on tooth morphology in carnivorous mammals. BMC Evol Biol. 2012;12:146. doi: 10.1186/1471-2148-12-146 22899809

7. Tarquini SD, Chemisquy MA, Prevosti FJ. Evolution of carnassial in living mammalian Carnivores (Carnivora, Didelphimorphia, Dasyuromorphia): diet, phylogeny, and allometry. J Mamm Evol. 2018, Aug 23. doi: 10.1007/s10914-018-9448-7

8. Nadachowski A. Late Quaternary rodents of Poland, with special reference to morphotype dentition analysis of voles. Warszawa: Państwowe Wydawnictwo Naukowe; 1982.

9. Pozdnyakov AA. Morphotypic variation in molar masticatory surface pattern in gray voles of the “maximowiczi” group (Rodentia, Arvicolidae): Experience in quantitative statistical analysis. Zool Zh. 1993;72:114–125.

10. Markova EA. Assessment of tooth complexity in Arvicolines (Rodentia). A morphotype ranking approach. Biol Bull. 2014;41:589–600.

11. Szuma E. Dental polymorphism in a population of the red fox (Vulpes vulpes) from Poland. J Zool, Lond. 2002;256:243–253.

12. Szuma E. Evolutionary implications of morphological variation in the lower carnassial of red fox Vulpes vulpes. Acta Theriol. 2004;49:433–447.

13. Szuma E. Geography of dental polymorphism in the red fox Vulpes vulpes and its evolutionary implications. J Zool, Lond. 2007;90:61–84.

14. Gimranov DO. Species diagnostics of the corsac (Vulpes corsac), fox (Vulpes vulpes) and arctic fox (Vulpes lagopus, Carnivora, Canidae) using the upper teeth. Zool Zh. 2017;96:684–697.

15. Gimranov DO, Kosintsev PA, Gasilin VV. Species diagnostic of corsac (Vulpes corsac), fox (Vulpes vulpes), and arctic fox (Vulpes lagopus) according to ontological characteristics of mandible teeth. Zool Zh. 2015;94:1338–1350.

16. Szuma E. Ecological and evolutionary determinants of dental polymorphism in the arctic fox Vulpes (Alopex) lagopus. Ann Zool Fenn. 2011;48:191–213.

17. Korablev NP, Szuma E, Korablev PN, Zinoviev AV. Dental polymorphism of the raccoon dog in indigenous and invasive populations: internal and external causation. Mamm Res. 2017;62:163–177. doi: 10.1007/s13364-016-0293-x

18. Wolsan M. Dental polymorphism in the genus Martes (Carnivora: Mustelidae) and its evolutionary significance. Acta Theriol. 1989;34:545–593.

19. Gimranov DO, Kosintsev PA. Differentiation of three Martes species (M. martes, M. zibellina, M. foina) by tooth morphotypes. CR Palevol. 2015;14:647–656.

20. Baryshnikov GF, Potapova OR. Variability of the dental system in badgers (Meles, Carnivora) of the USSR fauna. Zool Zh. 1990;69:84–97. Russian.

21. Szuma E. Quasi-continuous variation of the first premolars in the Polish population of the badger Meles meles. Acta Theriol. 1994;39:201–208.

22. Gimranov DO. Morphotypic characteristics of the fourth premolars of the brown (Ursus arctos) and polar (Ursus maritimus) bears (Carnivora, Ursidae). Zool Zh. 2018;97:205–223.

23. Serdyuk VA. Variability of teeth structure in Citellus parryi rich. and possible ways of migration of this species in the north-east of USSR. Zool Zh. 1979;58:1692–1702.

24. Suchentrunk F. Epigenetic dental asymmetry of Israeli hares: Developmental stability along an environmental gradient. Isr J Zool. 2000;46:103–118.

25. Suchentrunk F. Phylogenetic relationships between Indian and Burmese hares (Lepus nigricollis and L. peguensis) inferred from epigenetic dental characters. Mamm Biol. 2004;69:28–45.

26. Suchentrunk F, Flux JEC. Minor dental traits in East African cape hares and savanna hares (Lepus capensis and Lepus victoriae): A study of intra- and interspecific variability. J Zool, Lond. 1996;238:495–511.

27. Suchentrunk F, Willing R, Hartl GB. Nonmetrical polymorphism of the first lower premolar (P-3) in Austrian brown hares (Lepus europaeus)–a study on regional differentiation. J Zool, Lond. 1994;232:79–91.

28. Maleeva AG. Evolutsiya gryzunov i istoriya formirovaniya ikh sovremennoy fauny. In: Sbornik nauchnykh trudov, Leningrad 1976:48–57. Russian.

29. Markova EA, Borodin AV, Gileva EA. Odontological features in Microtus arvalis s.str. and M. rossiaemeridionalis from the Ural region and their taxonomic value. Zool Zh. 2003;82:1086–1094.

30. Markova EA, Malygin V, Montuire S, Nadachowski A, Quéré J-P, Ochman K. Dental Variation in sibling species Microtus arvalis and M. rossiaemeridionalis (Arvicolinae, Rodentia): Between-species comparisons and geography of morphotype dental patterns. J Mam Evol. 2010;17:121–139.

31. Nadachowski A. Systematic, geographic variation and evolution of snow voles (Chionomys) based on dental characters. Acta Theriol. 1991;36:1–45.

32. Pozdnyakov AA. The structure of morphological variability of common voles (Microtus: Rodentia, Arvicolidae) in the light of epigenetic theory of evolution. Usp sovrem biol. 2007;127:416–424.

33. Berthaume MA, Schroer K. Extant ape dental topography and its implications for reconstructing the emergence of early Homo. J Hum Evol. 2017;112:15–29. doi: 10.1016/j.jhevol.2017.09.001 29037413

34. Ercoli MD, Candela AM, Rasia LL, Ramírez MA. Dental shape variation of Neogene Pachyrukhinae (Mammalia, Notoungulata, Hegetotheriidae): systematics and evolutionary implications for the late Miocene Paedotherium species. J Syst Paleontol. 2018;16:1073–1095.

35. Gomez Cano AR, Kimura Y, Blanco F, Menéndez I, Álvarez-Sierra MA, Fernández MH. Ecomorphological characterization of murines and non-arvicoline cricetids (Rodentia) from south-western Europe since the latest Middle Miocene to the Mio-Pliocene boundary (MN 7/8_MN13). PeerJ 2017;5:e3646. doi: 10.7717/peerj.3646 28966888

36. Zazhigin VS, Voyta LL. A new middle Miocene crocidosoricine shrew from the Mongolian Shargain Gobi Desert. Acta Palaeontol Pol. 2018;63:171–187.

37. Szuma E. Microevolutionary trends in the dentition of the red fox (Vulpes vulpes). J Zool Syst Evol Res. 2003;41:47–56.

38. Kay RF. Molar structure and diet in extant Cercopithecidae. In: Butler PM, Joysey KA (eds) Development, function and evolution of teeth. New York: Academic Press; 1978.

39. Strait SG. Molar morphology and food texture among small-bodied insectivorous mammals. J Mammal. 1993;74:391–402.

40. Yom-Tov Y, Yom-Tov S, Baagøe HJ. Increase of skull size in the red fox (Vulpes vulpes) and Eurasian badger (Meles meles) in Denmark during the 20th century: an effect of improved diet? Evol Ecol Res. 2003;5:1037–1048.

41. Yom-Tov Y, Yom-Tov S, Zachos FE. Temporal and geographical variation in skull size of the red fox (Vulpes vulpes) and the Eurasian badger (Meles meles) in Austria. J Linn Soc, Lond. 2013;108:579–585.

42. Otte M, Baden-Powell D. Le Magdalénien du Trou de Chaleux (Hulsonniaux-Belgique). ERAUL 1994;60:241–253.

43. Dupont E. L’homme pendant les âges de la pierre dans les environs de Dinant-sur-Meuse. 2nd ed. Bruxelles: C Muquardt; 1873.

44. Germonpré M, Sablin MV. Systematics and osteometry of Late Glacial foxes from Belgium. Bull Inst r sci nat Belg. 2004;74:175–188.

45. Germonpré M. A reconstruction of the spatial distribution of the faunal remains from Goyet, Belgium. Notae praehistoricae. 2001;21:57–65.

46. Posth C, Renaud G, Mittnik A, Drucker D, Rougier H, Cupillard Ch, et al. Pleistocene mitochondrial genomes suggest a single major dispersal of Non-Africans and Late Glacial Population turnover in Europe. Curr Biol. 2016;26:827–833. doi: 10.1016/j.cub.2016.01.037 26853362

47. Rougier H, Crevecoeur I, Beauval C, Posth C, Flas D, Wiẞing Ch, et al. Neandertal canibalism and Neandertal bones used as tools in Northern Europe. Sci Rep. 2016;6:29005. doi: 10.1038/srep29005 27381450

48. Germonpré M. Preliminary results on the mammals of the Magdalenian upper horizon of Goyet (Belgium). Notae Praehistoricae. 1996;16:75–85.

49. Charles R. Late Magdalenian chronology and faunal exploitation in the North-Western Ardennes. In: BAR Int Ser. 1998;737:1–246.

50. Germonpré M, Sablin MV, Stevens RE, Hedges REM, Hofreiter M, Stiller M, et al. Fossil dogs and wolves from Paleolithic sites in Belgium, the Ukraine and Russia: osteometry, ancient DNA and stable isotopes. J Archaeol Sci. 2009;36:473–490.

51. Otte M. Le paléolithique supérieur ancien en Belgique. Monographies d’Archéologie Nationale 5. Bruxelles: Musées Royaux d’art et d’histoire; 1979.

52. Jimenez E-L. Palaeoecology and subsistence strategies in Belgium and Northwestern Europe during the MIS 3 through the reassessment of forgotten collections: A methodological approach. PIA 2016;26(1):Art.3. doi: 10.5334/pia-486

53. Jimenez E-L, Smolderen A, Jadin I, Germonpré M. Exhumation de la collection faunique d’Édouard Dupont provenant du Trou Magrite (Pont-à-Lesse) Quelles données et quelles perspectives pour une collection do XIXe siècle? Notae Praehistoricae 2016;36:167–190. Franch.

54. Jimenez E-L. Modalités d’occupation du territoire et relations humains-grands carnivores durant le Pléistocène supérieur. Approche archéozoologique, taphonomique et paléoécologique du Bassin mosan belge dans son contexte nord-ouest européen. [dissertation], Histoire de l’art et Archéologie, Université Libre de Bruxelles; 2017. Franch.

55. Szuma E. Geographic variation of tooth and skull sizes in the arctic fox Vulpes (Alopex) lagopus. Ann Zool Fenn. 2008;45:185–199.

56. Szuma E. Evolutionary and climatic factors affecting tooth size in the red fox Vulpes vulpes in the Holarctic. Acta Theriol. 2008;53:289–332.

57. Szuma E. Variation and correlation patterns in the dentition of the red fox from Poland. Ann Zool Fenn. 2000;37:113–127.

58. Gilbert N. Biometrical interpretation. Making sense of statistics in Biology. 2nd ed. Oxford, New York: Oxford University Press; 1989.

59. Bonis de L, Peigné S, Likius A, Mackaye HT, Vignaud P, Brunet M. The oldest African fox (Vulpes riffautae n. sp., Canidae, Carnivora) recovered in the late Miocene deposits of the Djurab desert, Chad. Naturwissenschaften. 2007;94:575–580. doi: 10.1007/s00114-007-0230-6 17361401

60. Wang X, Tseng ZJ, Li Q, Takeuchi GT, Xie G. From ‘third pole’ to north pole: a Himalayan origin for the arctic fox. Proc Royal Soc. B 2014, Jul 22. 281(1787). doi: 10.1098/rspb.2014.0893 24920475

61. Deng T, Wang X, Wu F, Wang Y, Li Q, Wang S, et al. Review: implications of vertebrate fossils for paleo-elevations of the Tibetan Plateau. Glob Planet Change. 2019;174:58–69.

62. Szuma E, Germonpré M. Size of the lower carnassial in the arctic and the red fox from Late Pleistocene in Belgium compared to other ancient and extant populations. Mammal Res. doi: 10.1007/s13364-019-00459-w

63. Deng T, Wang XM, Fortelius M, Li Q, Wang Y, Tseng ZJ, et al. Out of Tibet: Pliocene woolly rhino suggests highplateau origin of Ice Age megaherbivores. Science. 2011;333:1285–1288. doi: 10.1126/science.1206594 21885780

64. Elmhagen B, Rushton SP. Trophic control of mesopredators in terrestrial ecosystems: top-down or bottom-up? Ecol Lett. 2007;10:197–206. doi: 10.1111/j.1461-0248.2006.01010.x 17305803

65. Selås V, Vik JO. Possible impact of snow depth and ungulate carcasses on red fox (Vulpes vulpes) populations in Norway, 1897–1976. J Zool. 2006;268:299–308.

66. Hersteinsson P, Macdonald DW. Interspecific competition and the geographical-distribution of red and arctic foxes Vulpes vulpes and Alopex lagopus. Oikos. 1990;64:505–515.

67. Roff DA, Stirling G, Fairbairn DJ. The evolution of threshold traits: A quantitative genetic analysis of the physiological and life-history correlates of wing dimorphism in the sand cricket. Evolution. 1997;51:1910–1919. doi: 10.1111/j.1558-5646.1997.tb05113.x 28565097

68. Wang X, Tedford RH. Dogs. Their fossil relatives abd evolutionary history. New York: Columbia University Press; 2010.

69. Qiu Z, Tedford RH. A Pliocene species of Vulpes from Yushe, Shanxi. Vertebrat PalAsiatic. 1990;10:245–258.

70. Szuma E. Microevolutionary trends in the dentition of the Red fox (Vulpes vulpes). J Zool Syst Evol Res. 2002,40:1–10.

71. Dalén L, Nyström V, Valdiosera C, Germonpré M, Sablin M, Turner E, et al. Ancient DNA reveals lack of postglacial habitat tracking in the arctic fox. PNAS. 2007;104:6726–6729. doi: 10.1073/pnas.0701341104 17420452

72. Flower LDH. New body mass estimates of British Pleistocene wolves: Palaeoenvironmental implications and competitive interactions. Quat Sci Rev. 2016;149:230–247

73. Meiri S, Dayan T, Simberloff D. Carnivores, biases and Bergmann’s rule. Biol J Linn Soc. 2004;81:579–588.


Článok vyšiel v časopise

PLOS One


2020 Číslo 1
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Získaná hemofilie - Povědomí o nemoci a její diagnostika
nový kurz

Eozinofilní granulomatóza s polyangiitidou
Autori: doc. MUDr. Martina Doubková, Ph.D.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#