Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe.

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Pan-European phylogeography of the European roe deer (Capreolus capreolus).

To provide the most comprehensive picture of species phylogeny and phylogeography of European roe deer (Capreolus capreolus), we analyzed mtDNA control region (610 bp) of 1469 samples of roe deer from Central and Eastern Europe and included into the analyses additional 1541 mtDNA sequences from GenBank from other regions of the continent. We detected two mtDNA lineages of the species: European and Siberian (an introgression of C. pygargus mtDNA into C. capreolus). The Siberian lineage was most frequent in the eastern part of the continent and declined toward Central Europe. The European lineage contained three clades (Central, Eastern, and Western) composed of several haplogroups, many of which were separated in space. The Western clade appeared to have a discontinuous range from Portugal to Russia. Most of the haplogroups in the Central and the Eastern clades were under expansion during the Weichselian glacial period before the Last Glacial Maximum (LGM), while the expansion time of the Western clade overlapped with the Eemian interglacial. The high genetic diversity of extant roe deer is the result of their survival during the LGM probably in a large, contiguous range spanning from the Iberian Peninsula to the Caucasus Mts and in two northern refugia.

Environmental factors shaping stable isotope signatures of modern red deer (Cervus elaphus) inhabiting various habitats.

Stable isotope analyses of bone collagen are often used in palaeoecological studies to reveal environmental conditions in the habitats of different herbivore species. However, such studies require valuable reference data, obtained from analyses of modern individuals, in habitats of well-known conditions. In this article, we present the stable carbon and nitrogen isotope composition of bone collagen from modern red deer (N = 242 individuals) dwelling in various habitats (N = 15 study sites) in Europe. We investigated which of the selected climatic and environmental factors affected the δ13C and δ15N values in bone collagen of the studied specimens. Among all analyzed factors, the percent of forest cover influenced the carbon isotopic composition most significantly, and decreasing forest cover caused an increase in δ13C values. The δ15N was positively related to the proportion of open area and (only in the coastal areas) negatively related to the distance to the seashore. Using rigorous statistical methods and a large number of samples, we confirmed that δ13C and δ15N values can be used as a proxy of past habitats of red deer.

First Insights into the Population Genetic Structure and Heterozygosity-Fitness Relationship in Roe Deer Inhabiting the Area between the Alps and Dinaric Mountains.

Across its pan-European distribution, the European roe deer (Capreolus capreolus) faces a wide diversity of environmental and climatic conditions; therefore, several factors, including intrinsic ones, shape life-history traits and cause significant variability in parameters of fitness. By utilizing microsatellite variations in 214 roe deer females collected throughout Slovenia, Central Europe, we determined the genetic variability and population structure of this species in the contact zone between the Alps and the Dinaric Mountains, i.e., over a wider area where data on the genetic outlook for this-the most common and widespread European wild ungulate-have been completely lacking so far. Throughout the country, we found moderate microsatellite diversity (Ho = 0.57-0.65) in relation to the observed heterozygosity reported for other roe deer populations in Europe. Spatial differences in genetic diversity of the species in Slovenia can be explained by population history linked to varying approaches to population management and/or different connectivity among subpopulations in topographically differentiated habitats. A country-wide pattern of genetic structure is clearly defined by separation of the populations into three groups present in the following regions: (i) Southern sub-Mediterranean and Karst regions, (ii) Central Slovenia, and (iii) the Sub-Pannonian Region in the north-east. This is also confirmed by evidencing a moderate isolation by distance, especially by separating southern samples (coastal Slovenia) from others. Levels of genetic differentiation vary among populations, which can be explained by the effect of natural geographical barriers or the presence of anthropogenic barriers such as urban areas and highways. In the subset of 172 yearling females, we analyzed the influence of genetic advantage (individual heterozygosity) and other genetic data (reflected in the structuring of the population) on body mass and reproductive ability. We found evidence that genetic factors influence the body mass of roe deer yearling females (explaining altogether 18.8% of body mass variance), and the level of individual heterozygosity alone also positively affected body mass, which is in accordance with the theory that heterozygosity is commonly positively correlated with fitness in wild populations. However, we did not uncover any effect of heterozygosity on two parameters of reproductive ability (fertility and potential reproductive outcome), indicating that several other factors, especially environmental ones, have a predominant effect on the parameters of fitness in roe deer.

Age-related effects of body mass on fertility and litter size in roe deer.

We analysed effects of females' body mass and age on reproductive capacity of European roe deer (Capreolus capreolus) in a large sample set of 1312 females (305 yearlings and 1007 adults), hunted throughout Slovenia, central Europe, in the period 2013-2015. Body mass positively affected probability of ovulation and potential litter size (number of corpora lutea), although its effect was more pronounced in yearlings than in adults. Between age groups, we found clear differences in responses of both reproductive parameters to body mass which influences primarily reproductive performance of younger, and in particular, lighter individuals: at the same body mass yearlings would at average have smaller litters than adults, and at lower body mass also young to middle-aged adults would have smaller litters than old ones. In addition, while yearlings have to reach a critical threshold body mass to attain reproductive maturity, adult females are fertile (produce ova) even at low body mass. However, at higher body mass also younger individuals shift their efforts into the reproduction, and after reaching an age-specific threshold the body mass does not have any further effects on the reproductive output of roe deer females. Increased reproductive capacity at more advanced age, combined with declining body mass suggests that old does allocate more of their resources in reproduction than in body condition.