Consequences for genetic diversity and population performance of introducing continental red deer into the northern distribution range

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.

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.

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.

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.

Red deer in Iberia: Molecular ecological studies in a southern refugium and inferences on European postglacial colonization history

The red deer (Cervus elaphus) is a widespread wild ungulate in Europe that has suffered strong anthropogenic impacts over their distribution during the last centuries, but also at the present time, due its economic importance as a game species. Here we focus on the evolutionary history of the red deer in Iberia, one of the three main southern refugial areas for temperate species in Europe, and addressed the hypothesis of a cryptic refugia at higher latitudes during the Last Glacial Maximum (LGM). A total of 911 individuals were sampled, genotyped for 34 microsatellites specifically developed for red deer and sequenced for a fragment of 670 bp of the mitochondrial (mtDNA) D-loop. The results were combined with published mtDNA sequences, and integrated with species distribution models and historical European paleo-distribution data, in order to further examine the alternative glacial refugial models and the influence of cryptic refugia on European postglacial colonization history. Clear genetic differentiation between Iberian and European contemporary populations was observed at nuclear and mtDNA levels, despite the mtDNA haplotypes central to the phylogenetic network are present across western Europe (including Iberia) suggesting a panmictic population in the past. Species distribution models, fossil records and genetic data support a timing of divergence between Iberian and European populations that overlap with the LGM. A notable population structure was also found within the Iberian Peninsula, although several populations displayed high levels of admixture as a consequence of recent red deer translocations. Five D-loop sub-lineages were found in Iberia that belong to the Western European mtDNA lineage, while there were four main clusters based on analysis of nuclear markers. Regarding glacial refugial models, our findings provide detailed support for the hypothesis that red deer may have persisted in cryptic northern refugia in western Europe during the LGM, most likely in southern France, southern Ireland, or in a region between them (continental shelf), and these regions were the source of individuals during the European re-colonization. This evidence heightens the importance of conserving the high mitochondrial and nuclear diversity currently observed in Iberian populations.

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.

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.

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

We analyzed more than 600 red deer (Cervus elaphus) from large parts of its European distribution range at 13 microsatellite loci, presenting the first continent-wide study of this species using nuclear markers. Populations were clearly differentiated (overall F ST = 0.166, Jost's D est = 0.385), and the BAPS clustering algorithm yielded mainly geographically limited and adjacent genetic units. When forced into only 3 genetic clusters our data set produced a very similar geographic pattern as previously found in mtDNA phylogeographic studies: a western group from Iberia to central and parts of Eastern Europe, an eastern group from the Balkans to Eastern Europe, and a third group including the threatened relict populations from Sardinia and Mesola in Italy. This result was also confirmed by a multivariate approach to analyzing our data set, a discriminant analysis of principal components. Calculations of genetic diversity and effective population sizes (linkage disequilibrium approach) yielded the lowest results for Italian (Sardinia, Mesola; N e between 2 and 8) and Scandinavian red deer, in line with known bottlenecks in these populations. Our study is the first to present comparative nuclear genetic data in red deer across Europe and may serve as a baseline for future analyses of genetic diversity and structuring in this widespread ungulate.

Complex patterns of population genetic structure of moose, Alces alces, after recent spatial expansion in Poland revealed by sex-linked markers

In recent years, human activity directly and indirectly influenced the demography of moose in Poland. The species was close to extinction, and only a few isolated populations survived after the Second World War; then, unprecedented demographic and spatial expansions had occurred, possibly generating a very complex pattern of population genetic structure at the present-day margins of the species range in Poland. Over 370 moose from seven populations were collected from Poland, and partial sequences of the mitochondrial control region (mtDNA-cr; 607 bp) were obtained. In addition, the entire mtDNA cytochrome b gene (1,140 bp) and Y-chromosome markers (1,982 bp in total) were studied in a chosen set of individuals. Twelve mtDNA haplotypes that all belonged to the European moose phylogroup were recorded. They could be divided into two distinct clades: Central Europe and the Ural Mountains. The first clade consists of three distinct groups/branches: Biebrza, Polesie, and Fennoscandia. The Biebrza group has experienced spatial and demographic expansion in the recent past. Average genetic differentiation among moose populations in Poland at mtDNA-cr was great and significant (Φ_ST = 0.407, p < 0.001). Using mtDNA-cr data, four separate groups of population were recognized using spatial analysis of molecular variance and principal coordinate analysis, including a relict population in Biebrza National Park, a reintroduced Kampinos National Park population, as well as populations that were descendants of moose that colonized Poland from the east (Lithuania, Belarus, and Ukraine) and the north (former East Prussia). Among all the sequenced Y-chromosome markers, polymorphisms were found in the DBY14 marker in three populations only; four haplotypes were recorded in total. No significant differentiation was detected for this Y-linked marker among moose populations in Poland. Our mtDNA study revealed that a variety of different factors—bottleneck, the presence of relict, autochthonous populations, translocations, limited female dispersal, and the colonization from the east and north—are responsible for the observed complex pattern of population genetic structure after demographic and spatial expansion of moose in Poland.

Evolution of population genetic structure of the British roe deer by natural and anthropogenic processes (Capreolus capreolus)

Human influence typically impacts on natural populations of conservation interest. These interactions are varied and sometimes complex, and may be negative and unintended or associated with conservation and management strategy. Understanding the details of how these interactions influence and are influenced by natural evolutionary processes is essential to the development of effective conservation strategies. In this study, we investigate a species in Britain that has experienced both negative impact through overhunting in historical times and management efforts through culls and translocations. At the same time, there are regional populations that have been less affected by human influence. We use mtDNA and nuclear microsatellite DNA markers to investigate patterns of connectivity and diversity and find multiple insular populations in Britain that probably evolved within the Holocene (when the habitat was free of ice). We identify three concurrent processes. First, surviving indigenous populations show highly provincial patterns of philopatry, maintaining and generating population structure on a small geographic scale. Second, founder populations into habitat extirpated of native populations have expanded, but remained largely insular. Third, introductions into established populations generate some admixture. We discuss the implications for the evolution of diversity of the integration of natural processes with anthropogenic influences on population size and distribution.