Admixture of Eastern and Western European Red Deer Lineages as a Result of Postglacial Recolonization of the Czech Republic (Central Europe)

Red Deer Resequencing Reveals the Importance of Sex Chromosomes for Reconstructing Late Quaternary Events

Sex chromosomes differ in their inheritance properties from autosomes and hence may encode complementary information about past demographic events. We compiled and analyzed a range-wide resequencing data set of the red deer (Cervus elaphus), one of the few Eurasian herbivores of the Late Pleistocene megafauna still found throughout much of its historic range. Our analyses of 144 whole genomes reveal striking discrepancies between the population clusters suggested by autosomal and X-chromosomal data. We postulate that the genetic legacy of Late Glacial population structure is better captured and preserved by the X chromosome than by autosomes, for two reasons. First, X chromosomes have a lower Ne and hence lose genetic variation faster during isolation in glacial refugia, causing increased population differentiation. Second, following postglacial recolonization and secondary contact, immigrant males pass on their X chromosomes to female offspring only, which effectively halves the migration rate when gene flow is male mediated. Our study illustrates how a comparison between autosomal and sex chromosomal phylogeographic signals unravels past demographic processes that otherwise would remain hidden.

Genome-wide SNP assessment of contemporary European red deer genetic structure highlights the distinction of peripheral populations and the main admixture zones in Europe

Genome-wide technologies open up new possibilities to clarify questions on genetic structure and phylogeographic history of taxa previously studied with microsatellite loci and mitochondrial sequences. Here, we used 736 individual red deer (Cervus elaphus) samples genotyped at 35,701 single nucleotide polymorphism loci (SNPs) to assess the population structure of the species throughout Europe. The results identified 28 populations, with higher degrees of genetic distinction in peripheral compared to mainland populations. Iberian red deer show high genetic differentiation, with lineages in Western and Central Iberia maintaining their distinctiveness, which supports separate refugial ranges within Iberia along with little recent connection between Iberian and the remaining Western European populations. The Norwegian population exhibited the lowest variability and the largest allele frequency differences from mainland European populations, compatible with a history of bottlenecks and drift during post-glacial colonization from southern refugia. Scottish populations showed high genetic distance from the mainland but high levels of diversity. Hybrid zones were found between Eastern and Western European lineages in Central Europe as well as in the Pyrenees, where red deer from France are in close contact with Iberian red deer. Anthropogenic restocking has promoted the Pyrenean contact zone, admixture events in populations on the Isle of Rum and in the Netherlands, and at least partly the admixture of the two main lineages in central-eastern Europe. Our analysis enabled detailed resolution of population structure of a large mammal widely distributed throughout Europe and contributes to resolving the evolutionary history, which can also inform conservation and management policies.

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.

Preservation of genetic diversity in a highly fragmented population of the gray-sided vole Myodes rufocanus in an intensive farming region

Individual dispersal plays an important role in preserving genetic diversity in density-fluctuating populations of arvicoline rodents. When habitats are fragmented and dispersal between habitats is severely constrained, genetic diversity can be lost. Here, I investigated whether genetic diversity in the gray-sided vole Myodes rufocanus was preserved in an intensive farming region in Japan, where voles inhabited isolated windbreak forests along the borders of plowed lands. Genetic structure was examined in 673 vole samples (330 in spring and 343 in fall) collected at 34 windbreak forests located 0.35-20 km apart. A part of the control region (425 bp) of mitochondrial DNA (mtDNA) was sequenced in 673 voles, yielding 76 haplotypes. Genetic differentiation of maternally inherited mtDNA among trapping sites was markedly lower in males than in females in both seasons, indicating strong male-biased dispersal. Genotypes at six microsatellite DNA loci were determined in 494 voles (245 in spring and 249 in fall) from 18 trapping sites, and loci harbored 16-24 alleles. The mean number of alleles per locus (allelic diversity) at trapping sites was positively correlated with the number of examined individuals (density) in both seasons, and the relationship was very similar to that of a previous study performed in much less fragmented populations. The genetic differentiation of microsatellite DNA among trapping sites decreased considerably from spring to fall. In a STRUCTURE analysis with a most probable cluster number of two, closer trapping sites showed more similar mean values of cluster admixture proportions. The present findings indicate that gene flow among isolated windbreak forests, which occurred mainly by dispersal of males, was not restrained in this intensive farming region. Furthermore, the results suggest that genetic diversity in the study population was preserved as well as in less fragmented populations.

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.

Uncovering the Fecal Bacterial Communities of Sympatric Sika Deer (Cervus nippon) and Wapiti (Cervus canadensis)

Simple Summary

There are many microbial communities in the digestive tracts of animals, and the complex gut microbiome constitutes an intricate ecosystem and intestinal microbial community which has co-adapted with its host. The intestinal microecology plays an important role in the host’s maintenance of normal physical activities, such as substance metabolism, energy transmission, signal transduction, and the immune system. This study used high-throughput sequencing technology to sequence the fecal microbiota of sika deer (Cervus nippon) and wapiti (Cervus canadensis) in order to explore the composition of, and the similarity between, the fecal microbiota structures of sika deer and wapiti in the similar living environment. The species composition, relative abundance of fecal microbiota, alpha diversity, and differences in beta diversity were analyzed. The maintenance of the composition of the gut microbiota and a balanced intestinal environment through the diet plays a key role in maintaining the host’s health. The results demonstrate that the fecal microbiota of sika deer and wapiti share a similar fecal microbiota structure, but there was some evidence showing that the gut microbiota of these two animals exhibit a clear divergence at the species level.

Abstract

Microbial symbiotic associations may be beneficial, neutral, or harmful to the host. Symbionts exploit the host space and nutrition or use hosts as carriers to spread to other environments. In order to investigate the fecal bacterial communities of wild sika deer (Cervus nippon) and wapiti (Cervus canadensis), this study aimed to sequence and explore the composition of, and similarity between, the fecal microbiota of sika deer and wapiti using high-throughput sequencing. The composition and relative abundance of fecal microbiota, alpha diversity, and differences in beta diversity between the two species were analyzed. We found that no pathogenic bacteria were present in large quantities in the hosts. The dominant bacterial phyla found in the two deer species were similar and included Firmicutes, Bacteroidetes, Proteobacteria, and Spirochaetes. Moreover, the deer also shared similar dominant genera, including the Rikenellaceae RC9 gut group, Ruminococcaceae_UCG-010, Ruminococcaceae_UCG-005, and Bacteroides. These results demonstrate that the sika deer and wapiti share a similar fecal microbiotal structure, probably due to their common diet and living environment, but there was some evidence of a difference at the species level. These analyses provide new insights into the health status of deer populations outside protected environments and offer a scientific framework for monitoring the health conditions of sika deer and wapiti.

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.

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.

Variation in the prion protein gene (PRNP) sequence of wild deer in Great Britain and mainland Europe

Susceptibility to prion diseases is largely determined by the sequence of the prion protein gene (PRNP), which encodes the prion protein (PrP). The recent emergence of chronic wasting disease (CWD) in Europe has highlighted the need to investigate PRNP gene diversity in European deer species, to better predict their susceptibility to CWD. Here we report a large genetic survey of six British deer species, including red (Cervus elaphus), sika (Cervus nippon), roe (Capreolus capreolus), fallow (Dama dama), muntjac (Muntiacus reevesii), and Chinese water deer (Hydropotes inermis), which establishes PRNP haplotype and genotype frequencies. Two smaller data sets from red deer in Norway and the Czech Republic are also included for comparison. Overall red deer show the most PRNP variation, with non-synonymous/coding polymorphisms at codons 98, 168, 226 and 247, which vary markedly in frequency between different regions. Polymorphisms P168S and I247L were only found in Scottish and Czech populations, respectively. T98A was found in all populations except Norway and the south of England. Significant regional differences in genotype frequencies were observed within both British and European red deer populations. Other deer species showed less variation, particularly roe and fallow deer, in which identical PRNP gene sequences were found in all individuals analysed. Based on comparison with PRNP sequences of North American cervids affected by CWD and limited experimental challenge data, these results suggest that a high proportion of wild deer in Great Britain may be susceptible to CWD.

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.

Introgression of exotic Cervus (nippon and canadensis) into red deer (Cervus elaphus) populations in Scotland and the English Lake District

Since the mid-19th century, multiple introductions of Japanese sika deer (Cervus nippon nippon) and North American wapiti (C. canadensis) have taken place in the British Isles. While wapiti have generally been unsuccessful, sika have been very successful, especially in Scotland where they now overlap at least 40% of the range of native red deer (C. elaphus). Hybridization between these two species and red deer has been demonstrated in captivity and in the wild. Using a panel of 22 microsatellite loci that are highly diagnostic between red deer and sika, and moderately diagnostic between red deer and wapiti, we investigated the extent of introgression between these species in 2,943 deer sampled from around Scotland and from the English Lake District using the Bayesian clustering software STRUCTURE. We also used a diagnostic mitochondrial marker for red deer and sika. Our survey extends previous studies indicating little introgression of wapiti nuclear alleles into red deer, in particular in Northern Scotland, Kintyre, and the Lake District. We found a new area of extensive sika introgression in South Kintyre. In the North Highlands, we show for the first time geographically scattered evidence of past hybridization followed by extensive backcrossing, including one red-like individual with sika introgression, two sika-like individuals with red deer introgression, and six individuals that were apparently pure sika at the nuclear markers assessed but which carried red deer mitochondria. However, there has not been a collapse of assortative mating in this region. Similarly, in the English Lake District red deer, we found only traces of past sika introgression. No sika alleles were detected in the Central Highlands or the Hebridean red deer refugia. We make suggestions for management to prevent further spread of sika alleles into red deer and vice versa.

Genetic diversity, genetic structure and diet of ancient and contemporary red deer (Cervus elaphus L.) from north-eastern France

In north-eastern France, red deer (Cervus elaphus L.) populations were rebuilt from a few hundred individuals, which have subsisted in remote valleys of the Vosges mountains, and to a lesser extent from individuals escaped from private enclosures; at present times, this species occupies large areas, mainly in the Vosges Mountains. In this study, we examined the population dynamics of red deer in the Vosges Mountains using ancient and contemporary mitochondrial DNA (mtDNA) from 140 samples (23 ancient + 117 modern) spanning the last 7'000 years. In addition, we reconstructed the feeding habits and the habitat of red deer since the beginning of agriculture applying isotopic analyses in order to establish a basis for current environmental management strategies. We show that past and present red deer in the Vosges Mountains belong to mtDNA haplogroup A, suggesting that they originated from the Iberian refugium after the last glacial maximum (LGM). Palaeogenetic analysis of ancient bone material revealed the presence of two distinct haplotypes with different temporal distributions. Individuals belonging to the two haplotype groups apparently occupied two different habitats over at least 7'000 years. AM6 correlates with an ecological type that feeds in densely forested mountain landscapes, while AM235 correlates with feeding in lowland landscapes, composed of a mixture of meadows and riverine, herb-rich woodlands. Our results suggest that red deer of north-eastern France was able to adapt, over the long term, to these different habitat types, possibly due to efficient ethological barriers. Modern haplotype patterns support the historical record that red deer has been exposed to strong anthropogenic influences as a major game species.

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.

Purifying Selection, Density Blocking and Unnoticed Mitochondrial DNA Diversity in the Red Deer, Cervus elaphus

The trajectories of postglacial range expansions, the occurrence of lineage patches and the formation and maintenance of secondary contact between lineages may mostly reflect neutral demographic processes, including density blocking, that may leave long-lasting genetic signatures. However, a few studies have recently shown that climate may also play a role. We used red deer, a large, mobile herbivore that is assumed to be sensitive to climate change, to test hypotheses of possible selection on the mitochondrial DNA cytochrome b gene (mtDNA cytb) and competitive and/or density-blocking (using mtDNA control region). We searched for a possible link between the phylogeographic structure and abiotic climatic variables. Finally, we tested for isolation by distance and isolation by environment and assessed the impact of human-mediated translocations on the genetic structure of red deer. Our analysis of 30 red deer populations in Poland using the mtDNA control region (N = 357) and cytochrome b (N = 50) markers not only confirmed the presence of the Western and South-Eastern lineages of the species but also indicated the presence of a previously unnoticed, rare relic haplotype that grouped together C. e. italicus from Italy (the Mesola deer). No significant signs of positive selection were detected for the mtDNA cytb gene in the studied red deer. However, a significant signal for purifying selection was found in our study that may explain the narrowness of the contact zone because gene flow between the Western and South-Eastern lineages should drive relatively strong mito-nuclear incompatibilities. MtDNA control region differentiation among red deer populations in Poland correlated with different abiotic climatic variables. Strikingly, the southernmost ice sheet limits during the Elsterian was the most important factor, and it explained the largest amount of variation. However, neither isolation by distance (IBD) nor isolation by environment (IBE) were recorded, and a very limited impact of human translocations was evident. The above-mentioned results suggest that in contemporary red deer populations in Poland, the phylogeographic pattern is well preserved, and long-term processes (density and/or competitive blocking) still play a major role.

Genetic structure and diversity of the selfing model grass Brachypodium stacei (Poaceae) in Western Mediterranean: out of the Iberian Peninsula and into the islands

Annual Mediterranean species of the genus Brachypodium are promising model plants for energy crops since their selfing nature and short-life cycles are an advantage in breeding programs. The false brome, B. distachyon, has already been sequenced and new genomic initiatives have triggered the de-novo genome sequencing of its close relatives such as B. stacei, a species that was until recently mistaken for B. distachyon. However, the success of these initiatives hinges on detailed knowledge about the distribution of genetic variation within and among populations for the effective use of germplasm in a breeding program. Understanding population genetic diversity and genetic structure is also an important prerequisite for designing effective experimental populations for genomic wide studies. However, population genetic data are still limited in B. stacei. We therefore selected and amplified 10 nuclear microsatellite markers to depict patterns of population structure and genetic variation among 181 individuals from 19 populations of B. stacei occurring in its predominant range, the western Mediterranean area: mainland Iberian Peninsula, continental Balearic Islands and oceanic Canary Islands. Our genetic results support the occurrence of a predominant selfing system with extremely high levels of homozygosity across the analyzed populations. Despite the low level of genetic variation found, two different genetic clusters were retrieved, one clustering all SE Iberian mainland populations and the island of Minorca and another one grouping all S Iberian mainland populations, the Canary Islands and all Majorcan populations except one that clustered with the former group. These results, together with a high sharing of alleles (89%) suggest different colonization routes from the mainland Iberian Peninsula into the islands. A recent colonization scenario could explain the relatively low levels of genetic diversity and low number of alleles found in the Canary Islands populations while older colonization events are hypothesized to explain the high genetic diversity values found in the Majorcan populations. Our study provides widely applicable information about geographical patterns of genetic variation in B. stacei. Among others, the genetic pattern and the existence of local alleles will need to be adequately reflected in the germplasm collection of B. stacei for efficient genome wide association studies.

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.