Genetic Structure and Effective Population Sizes in European Red Deer (Cervus elaphus) at a Continental Scale: Insights from Microsatellite DNA

Comparative Analysis of Microsatellite and SNP Markers for Genetic Management of Red Deer

The analysis of population genetic structure and individual multilocus heterozygosity are crucial for wildlife management and conservation. Microsatellite markers have traditionally been used to assess these genetic parameters. However, single-nucleotide polymorphisms (SNPs) are becoming increasingly popular. Our goal here was to determine to what extent SNPs can provide better insights than microsatellites into the overall genetic status and population genetic processes in the species. To this end, we genotyped 210 red deer (Cervus elaphus) in the Spanish wild population with both 11 microsatellites and 31,712 SNPs. We compared parameters related to population genetic structure and individual multilocus heterozygosity obtained with both types of markers. Our results showed correlations between parameters measured using both microsatellites and SNPs, particularly those related to the level of genetic diversity and genetic differentiation. However, we found notably lower precision of microsatellites in measuring the distribution of genetic diversity among individuals. We conclude that microsatellites can be used to monitor the overall genetic status and detect broad patterns in red deer populations. Nevertheless, the greater precision of SNPs in inferring genetic structure and multilocus heterozygosity leads us to encourage scientists and wildlife managers to prioritize their use whenever possible.

Traces of Human-Mediated Selection in the Gene Pool of Red Deer Populations

In this study, we analysed the effect of human-mediated selection on the gene pool of wild and farmed red deer populations based on genotyping-by-sequencing data. The farmed red deer sample covered populations spread across seven countries and two continents (France, Germany, Hungary, Latvia, New Zealand, Poland, and Slovakia). The Slovak and Spain wild red deer populations (the latter one in a large game estate) were used as control outgroups. The gene flow intensity, relationship and admixture among populations were tested by the Bayesian approach and discriminant analysis of principal components (DAPC). The highest gene diversity (He = 0.19) and the lowest genomic inbreeding (FHOM = 0.04) found in Slovak wild population confirmed our hypothesis that artificial selection accompanied by bottlenecks has led to the increase in overall genomic homozygosity. The Bayesian approach and DAPC consistently identified three separate genetic groups. As expected, the farmed populations were clustered together, while the Slovak and Spanish populations formed two separate clusters. Identified traces of genetic admixture in the gene pool of farmed populations reflected a strong contemporary migration rate between them. This study suggests that even if the history of deer farming has been shorter than traditional livestock species, it may leave significant traces in the genome structure.

Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history

European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.

Genetic traces of dispersal and admixture in red deer (Cervus elaphus) populations from the Carpathian Basin

After the last glacial, the Carpathian Basin was repopulated from either eastward or northward colonisation routes for various species; one of these was the emblematic member of the European megafauna, the red deer, Cervus elaphus. We analysed 303 red deer individuals from the middle of the region, in seven Hungarian game reserves, at ten microsatellite loci (C01, C229, T26, T108, T123, T156, T172, T193, T501, T507), to investigate the genetic diversity of these subpopulations. We discovered high levels of genetic diversity of red deer subpopulations; allelic richness values ranging 4.99–7.01, observed heterozygosity 0.729–0.800, polymorphic information content 0.722–0.806, and Shannon’s information index 1.668–2.064. Multi-locus analyses indicated population admixtures of various degrees that corresponded to geographical location, and complex genetic structures were shown by clustering. Populations in the south-western and the north-eastern parts of the region formed two highly separated groups, and the red deer from populations in between them were highly admixed (in western Pannonia/Transdanubia, where the Danube flows into the Carpathian Basin). This pattern corresponds to the distribution of mitochondrial as well as Y-chromosome lineages. Assignment tests showed that a large fraction of individuals (29.4%) are found outside of their population of origin, indicating that the dispersal of red deer is rather common, which could be expected considering the life course of the species.

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

Genetic insights into an Apennine population of the Italian red deer

The red deer Cervus elephus has been a common species in Italy until the Middle Ages and the Renaissance, when its distribution range started to considerably decrease, due to gradual deforestation and hunting pressure. Afterwards, the red deer has been reintroduced to many regions of the world, including Italy. In the Italian Apennines, the Acquerino-Cantagallo Natural Reserve (ACQUERINO) hosts one of the largest peninsular red deer populations, originated from a series of successful reintroductions. In this study, we meant to detect the level of genetic variability of Acquerino-Cantagallo Natural Reserve deer population and to investigate the genetic relationships with the other Italian and European populations. We identified five mitochondrial DNA control region (D-loop) haplotypes, four falling in lineage A and one falling in lineage C, derived from at least two maternal lineages, confirming that ACQUERINO population should be the result of multiple reintroductions. Haplotype diversity (H = 0.50) and nucleotide (π = 0.004) diversity were low, but included into the deer range values. ACQUERINO population showed low levels of genetic diversity when compared to other European and Mediterranean populations, confirming that this expanding population may have been generated from a low number of founders.

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

Vocal phenotype of male rutting roars and genetic markers delineate East European red deer (Cervus elaphus) from Central and West European populations

This study investigates a population of red deer Cervus elaphus, founded by 10 individuals introduced in the nineteenth century from Germany to the Voronezh region of the European part of Southern Russia and then developed without further introductions. We characterize for the first time the vocal phenotype of the Voronezh red deer male rutting calls in comparison with similar data on the Pannonian (native Central European) and Iberian (native West European) red deer obtained by the authors during preceding studies. In addition, we provide for the first time the genetic data on Pannonian red deer. In Voronezh stags, the number of roars per bout (2.85 ± 1.79) was lower than in Pannonian (3.18 ± 2.17) but higher than in Iberian (2.11 ± 1.71) stags. In Voronezh stags, the duration of main (the longest within bouts) roars was longer (2.46 ± 1.14 s) than in Pannonian (1.13 ± 0.50 s) or Iberian (1.90 ± 0.50 s) stags. The maximum fundamental frequency of main roars was similar between Voronezh (175 ± 60 Hz) and Pannonian (168 ± 61 Hz) but higher in Iberian stags (223 ± 35 Hz). Mitochondrial cytochrome b gene analysis of red deer from the three study populations partially supports the bioacoustical data, of closer similarity between Voronezh and Pannonian populations. In contrast, microsatellite DNA analysis delineates Voronezh red deer from either Pannonian or Iberian red deer. We discuss that population bottlenecking might affect the acoustics of the rutting roars, in addition to genotype.

Human-driven genetic differentiation in a managed red deer population

Nineteen red deer areas in a densely populated region with a huge network of fenced motorways and the division into administrative management units (AMUs) with restricted ecological connectivity were investigated. In the season 2018/2019, a total of 1291 red deer samples (on average 68 per area) were collected and genotyped using 16 microsatellite markers. The results show a clear genetic differentiation between most of the AMUs. Fourteen AMUs may be combined into four regions with a considerable internal genetic exchange. Five areas were largely isolated or showed only a limited gene flow with neighbouring areas. Ten of the 19 AMUs had an effective population size below 100. Effective population sizes greater than 500–1000, required to maintain the evolutionary potential and a long-term adaptation potential, were not achieved by any of the studied AMUs, even when AMUs with an appreciable genetic exchange were aggregated. Substantial genetic differentiation between areas can be associated with the presence of landscape barriers hindering gene flow, but also with the maintenance of ‘red deer–free’ areas. Efforts to sustainably preserve the genetic diversity of the entire region should therefore focus on measures ensuring genetic connectivity. Opportunities for this goal arise from the establishment of game bridges over motorways and from the protection of young male stags migrating through the statutory ‘red deer–free’ areas.

Mining the red deer genome (CerEla1.0) to develop X-and Y-chromosome-linked STR markers

Microsatellites are widely applied in population and forensic genetics, wildlife studies and parentage testing in animal breeding, among others, and recently, high-throughput sequencing technologies have greatly facilitated the identification of microsatellite markers. In this study the genomic data of Cervus elaphus (CerEla1.0) was exploited, in order to identify microsatellite loci along the red deer genome and for designing the cognate primers. The bioinformatics pipeline identified 982,433 microsatellite motifs genome-wide, assorted along the chromosomes, from which 45,711 loci mapped to the X- and 1096 to the Y-chromosome. Primers were successfully designed for 170,873 loci, and validated with an independently developed autosomal tetranucleotide STR set. Ten X- and five Y-chromosome-linked microsatellites were selected and tested by two multiplex PCR setups on genomic DNA samples of 123 red deer stags. The average number of alleles per locus was 3.3, and the average gene diversity value of the markers was 0.270. The overall observed and expected heterozygosities were 0.755 and 0.832, respectively. Polymorphic Information Content (PIC) ranged between 0.469 and 0.909 per locus with a mean value of 0.813. Using the X- and Y-chromosome linked markers 19 different Y-chromosome and 72 X-chromosome lines were identified. Both the X- and the Y-haplotypes split to two distinct clades each. The Y-chromosome clades correlated strongly with the geographic origin of the haplotypes of the samples. Segregation and admixture of subpopulations were demonstrated by the use of the combination of nine autosomal and 16 sex chromosomal STRs concerning southwestern and northeastern Hungary. In conclusion, the approach demonstrated here is a very efficient method for developing microsatellite markers for species with available genomic sequence data, as well as for their use in individual identifications and in population genetics studies.

Major inconsistencies of inferred population genetic structure estimated in a large set of domestic horse breeds using microsatellites

STRUCTURE remains the most applied software aimed at recovering the true, but unknown, population structure from microsatellite or other genetic markers. About 30% of structure-based studies could not be reproduced (Molecular Ecology, 21, 2012, 4925). Here we use a large set of data from 2,323 horses from 93 domestic breeds plus the Przewalski horse, typed at 15 microsatellites, to evaluate how program settings impact the estimation of the optimal number of population clusters K opt that best describe the observed data. Domestic horses are suited as a test case as there is extensive background knowledge on the history of many breeds and extensive phylogenetic analyses. Different methods based on different genetic assumptions and statistical procedures (dapc, flock, PCoA, and structure with different run scenarios) all revealed general, broad-scale breed relationships that largely reflect known breed histories but diverged how they characterized small-scale patterns. structure failed to consistently identify K opt using the most widespread approach, the ΔK method, despite very large numbers of MCMC iterations (3,000,000) and replicates (100). The interpretation of breed structure over increasing numbers of K, without assuming a K opt, was consistent with known breed histories. The over-reliance on K opt should be replaced by a qualitative description of clustering over increasing K, which is scientifically more honest and has the advantage of being much faster and less computer intensive as lower numbers of MCMC iterations and repetitions suffice for stable results. Very large data sets are highly challenging for cluster analyses, especially when populations with complex genetic histories are investigated.

Genetic epidemiology of the Alpine ibex reservoir of persistent and virulent brucellosis outbreak

While it is now broadly accepted that inter-individual variation in the outcomes of host-pathogen interactions is at least partially genetically controlled, host immunogenetic characteristics are rarely investigated in wildlife epidemiological studies. Furthermore, most immunogenetic studies in the wild focused solely on the major histocompatibility complex (MHC) diversity despite it accounts for only a fraction of the genetic variation in pathogen resistance. Here, we investigated immunogenetic diversity of the Alpine ibex (Capra ibex) population of the Bargy massif, reservoir of a virulent outbreak of brucellosis. We analysed the polymorphism and associations with disease resistance of the MHC Class II Drb gene and several non-MHC genes (Toll-like receptor genes, Slc11A1) involved in the innate immune response to Brucella in domestic ungulates. We found a very low neutral genetic diversity and a unique MHC Drb haplotype in this population founded few decades ago from a small number of individuals. By contrast, other immunity-related genes have maintained polymorphism and some showed significant associations with the brucellosis infection status hence suggesting a predominant role of pathogen-mediated selection in their recent evolutionary trajectory. Our results highlight the need to monitor immunogenetic variation in wildlife epidemiological studies and to look beyond the MHC.

Genetic analysis of red deer (Cervus elaphus) administrative management units in a human-dominated landscape

Red deer (Cervus elaphus) throughout central Europe are influenced by different anthropogenic activities including habitat fragmentation, selective hunting and translocations. This has substantial impacts on genetic diversity and the long-term conservation of local populations of this species. Here we use genetic samples from 480 red deer individuals to assess genetic diversity and differentiation of the 12 administrative management units located in Schleswig Holstein, the northernmost federal state in Germany. We applied multiple analytical approaches and show that the history of local populations (i.e., translocations, culling of individuals outside of designated red deer zones, anthropogenic infrastructures) potentially has led to low levels of genetic diversity. Mean expected heterozygosity was below 0.6 and we observed on average 4.2 alleles across 12 microsatellite loci. Effective population sizes below the recommended level of 50 were estimated for multiple local populations. Our estimates of genetic structure and gene flow show that red deer in northern Germany are best described as a complex network of asymmetrically connected subpopulations, with high genetic exchange among some local populations and reduced connectivity of others. Genetic diversity was also correlated with population densities of neighboring management units. Based on these findings, we suggest that connectivity among existing management units should be considered in the practical management of the species, which means that some administrative management units should be managed together, while the effective isolation of other units needs to be mitigated.

Progesterone and Cortisol Levels in Blood and Hair of Wild Pregnant Red Deer (Cervus Elaphus) Hinds

The red deer (Cervus elaphus L., 1758) is one of the largest deer species in the world. Females are seasonal polyestrous, with negative photoperiod: the increase of the night peak of melatonin determines the secretion of GnRH and, therefore, LH and FSH. To date there is little information regarding the hormonal control during pregnancy for this species; this could be due to the difficulty of sampling wild subjects, while farmed animals' hormonal concentrations may not reflect the physiology of the animal in a natural state. In this study we evaluated the concentration of cortisol and progesterone, extracted from blood and hair, on 10 wild and pregnant red deer females. Belonging to the population of the Bolognese Apennines (Italy), the hinds were sampled in the January-March 2018 period, according to the regional selective hunting plan. Plasma progesterone (P4) ranged from a minimum of 1.9 to a maximum of 7.48 ng/mL; while hair P4 concentrations varied from 41.68 to 153.57 pg/mg. The plasma and hair cortisol ranges are respectively 0.4-2.97 ng/mL and 0.03-0.55 pg/mg; the only significant correlation was found between hair concentration of P4 and the date of death. The results of this preliminary study represent a small step towards a better knowledge of this species' physiology during pregnancy.

New insights into population structure of the European golden eagle (Aquila chrysaetos) revealed by microsatellite analysis

Connectivity between golden eagle (Aquila chrysaetos) populations is poorly understood. Field studies exploring natal dispersal suggest that this raptor is a philopatric species, but with the ability to roam far. However, little is known about the population structure of the species in Europe. Our study is based on 14 microsatellite loci and is complemented by new and previously published mitochondrial control region DNA data. The present dataset includes 121 eagles from Scotland, Norway, Finland, Estonia, the Mediterranean and Alpine regions. Our sampling focused on the Alpine and Mediterranean populations because both mitochondrial DNA (mtDNA) lineages found in golden eagles, the Holarctic and the Mediterranean, are known to co-occur there. Cluster analyses of nuclear DNA support a shallow split into northern and southern populations in Europe, similar to the distribution of the two mtDNA lineages, with the Holarctic lineage occurring in the north and the Mediterranean lineage predominating in the south. Additionally, Scotland shows significant differentiation and low relative migration levels that indicate isolation from the mainland populations. Alpine and Mediterranean golden eagles do not show nuclear structure corresponding to divergent mtDNA lineages. This indicates that the presence of northern Holarctic mitochondrial haplotypes in the Alps and the Mediterranean is attributable to past admixture rather than recent long-distance dispersal.

Lessons for conservation management: Monitoring temporal changes in genetic diversity of Cape mountain zebra (Equus zebra zebra)

The Cape mountain zebra (Equus zebra zebra) is a subspecies of mountain zebra endemic to South Africa. The Cape mountain zebra experienced near extinction in the early 1900's and their numbers have since recovered to more than 4,800 individuals. However, there are still threats to their long-term persistence. A previous study reported that Cape mountain zebra had low genetic diversity in three relict populations and that urgent conservation management actions were needed to mitigate the risk of further loss. As these suggestions went largely unheeded, we undertook the present study, fifteen years later to determine the impact of management on genetic diversity in three key populations. Our results show a substantial loss of heterozygosity across the Cape mountain zebra populations studied. The most severe losses occurred at De Hoop Nature Reserve where expected heterozygosity reduced by 22.85% from 0.385 to 0.297. This is alarming, as the De Hoop Nature Reserve was previously identified as the most genetically diverse population owing to its founders originating from two of the three remaining relict stocks. Furthermore, we observed a complete loss of multiple private alleles from all populations, and a related reduction in genetic structure across the subspecies. These losses could lead to inbreeding depression and reduce the evolutionary potential of the Cape mountain zebra. We recommend immediate implementation of evidence-based genetic management and monitoring to prevent further losses, which could jeopardise the long term survival of Cape mountain zebra, especially in the face of habitat and climate change and emerging diseases.

Noninvasive sampling reveals population genetic structure in the Royle's pika, Ochotona roylei, in the western Himalaya

Understanding population genetic structure of climate-sensitive herbivore species is important as it provides useful insights on how shifts in environmental conditions can alter their distribution and abundance. Herbivore responses to the environment can have a strong indirect cascading effect on community structure. This is particularly important for Royle's pika (Lagomorpha: Ochotona roylei), a herbivorous talus-dwelling species in alpine ecosystem, which forms a major prey base for many carnivores in the Himalayan arc. In this study, we used seven polymorphic microsatellite loci to detect evidence for recent changes in genetic diversity and population structure in Royle's pika across five locations sampled between 8 and 160 km apart in the western Himalaya. Using four clustering approaches, we found the presence of significant contemporary genetic structure in Royle's pika populations. The detected genetic structure could be primarily attributed to the landscape features in alpine habitat (e.g., wide lowland valleys, rivers) that may act as semipermeable barriers to gene flow and distribution of food plants, which are key determinants in spatial distribution of herbivores. Pika showed low inbreeding coefficients (F IS) and a high level of pairwise relatedness for individuals within 1 km suggesting low dispersal abilities of talus-dwelling pikas. We have found evidence of a recent population bottleneck, possibly due to effects of environmental disturbances (e.g., snow melting patterns or thermal stress). Our results reveal significant evidence of isolation by distance in genetic differentiation (F ST range = 0.04-0.19). This is the first population genetics study on Royle's pika, which helps to address evolutionary consequences of climate change which are expected to significantly affect the distribution and population dynamics in this talus-dwelling species.

Intensive Management and Natural Genetic Variation in Red Deer (Cervus elaphus)

The current magnitude of big-game hunting has outpaced the natural growth of populations, making artificial breeding necessary to rapidly boost hunted populations. In this study, we evaluated if the rapid increase of red deer (Cervus elaphus) abundance, caused by the growing popularity of big-game hunting, has impacted the natural genetic diversity of the species. We compared several genetic diversity metrics between 37 fenced populations subject to intensive management and 21 wild free-ranging populations. We also included a historically protected population from a national park as a baseline for comparisons. Contrary to expectations, our results showed no significant differences in genetic diversity between wild and fenced populations. Relatively lower genetic diversity was observed in the protected population, although differences were not significant in most cases. Bottlenecks were detected in both wild and fenced populations, as well as in the protected population. Assignment tests identified individuals that did not belong to their population of origin, indicating anthropogenic movement. We discuss the most likely processes, which could have led to the observed high levels of genetic variability and lack of differentiation between wild and fenced populations and suggest cautionary points for future conservation. We illustrate our comparative approach in red deer. However, our results and interpretations can be largely applicable to most ungulates subject to big-game hunting as most of them share a common exploitation-recovery history as well as many ecological traits.

Using genetic tools to estimate the prevalence of non-native red deer (Cervus elaphus) in a Western European population

Game species like the red deer have been subjected to anthropogenic impacts for centuries. Translocations are often carried out—sometimes illegally—not only for sporting purposes, but also to increase trophy quality, reduce inbreeding, or mitigate bottlenecks after excessive persecution. Apart from the blurring of large‐scale genetic structure, translocations without adequate quarantine measure risk introducing pathogens into potentially immunologically naïve populations. It is therefore important to understand the frequency of clandestine translocations. Identification of non‐autochthonous animals and their potential origin is often difficult and, in red deer, has been hampered by the lack of large‐scale genotypic datasets for comparison. In the present study, we make use of a recently published European‐wide microsatellite dataset to detect and quantify the presence of non‐autochthonous red deer in a large population sample (n = 1,780) from Central Europe (Belgium). Using factorial correspondence analysis, assignment tests and Bayesian clustering algorithms we arrive at an estimate of 3.7% non‐autochthonous animals (or their descendants). Some of these animals were assigned to a nearby French population and may have immigrated into Belgium naturally, but the large majority must have been introduced by humans. Our analysis pointed to the British Isles and Germany/Poland as the potential origin of many introduced deer, regions known to have been source populations for translocations in Europe and beyond. We found evidence for recreational hunters using carcasses from farmed deer to fulfill mandatory hunting quotas. Our study is the first to quantify the extent of human‐mediated introductions in a European game species at such a large scale with large and representative sample sizes.

Phylogeography of the Tyrrhenian red deer (Cervus elaphus corsicanus) resolved using ancient DNA of radiocarbon-dated subfossils

We present ancient mitochondrial DNA analyses of 31 complete cytochrome b gene sequences from subfossil red deer remains from the Tyrrhenian islands (Corsica and Sardinia) and mainland Italy in a European-wide phylogeographic framework. Tyrrhenian and North African red deer, both going back to human introductions, were previously the only red deer to harbour the mitochondrial B lineage whose origin, however, remained unknown. Our ancient Italian samples from the central part of the peninsula that were radiocarbon-dated to an age of ca. 6300 to 15 600 cal BP all showed B haplotypes, closely related or even identical to those found on Sardinia. Genetic diversity in the mainland population was considerably higher than on the islands. Together with palaeontological evidence our genetic results identify the Italian Peninsula as the ultimate origin of the B lineage and thus the Tyrrhenian and North African red deer. This is in line with previous biogeographic findings that uncovered distinct intraspecific phylogeographic lineages in Italian mammals, underlining Italy’s status as a hotspot of European mammalian diversity.

Strong population structure in a species manipulated by humans since the Neolithic: the European fallow deer (Dama dama dama)

Species that have been translocated and otherwise manipulated by humans may show patterns of population structure that reflect those interactions. At the same time, natural processes shape populations, including behavioural characteristics like dispersal potential and breeding system. In Europe, a key factor is the geography and history of climate change through the Pleistocene. During glacial maxima throughout that period, species in Europe with temperate distributions were forced south, becoming distributed among the isolated peninsulas represented by Anatolia, Italy and Iberia. Understanding modern patterns of diversity depends on understanding these historical population dynamics. Traditionally, European fallow deer (Dama dama dama) are thought to have been restricted to refugia in Anatolia and possibly Sicily and the Balkans. However, the distribution of this species was also greatly influenced by human-mediated translocations. We focus on fallow deer to better understand the relative influence of these natural and anthropogenic processes. We compared modern fallow deer putative populations across a broad geographic range using microsatellite and mitochondrial DNA loci. The results revealed highly insular populations, depauperate of genetic variation and significantly differentiated from each other. This is consistent with the expectations of drift acting on populations founded by small numbers of individuals, and reflects known founder populations in the north. However, there was also evidence for differentiation among (but not within) physically isolated regions in the south, including Iberia. In those regions we find evidence for a stronger influence from natural processes than may be expected for a species with such strong, known anthropogenic influence.