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

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.

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

Simple Summary

We determined the genetic variability, population structure, and influence of genetic factors on two parameters of fitness (body mass and reproductive ability) in roe deer females in the contact zone between the Alps and the Dinaric Mountains by utilizing microsatellite variations in 214 individuals collected throughout Slovenia, Central Europe. Spatial differences in the genetic diversity of the species can be explained by population history, different approaches to population management and/or different connectivity among subpopulations. The population genetic structure confirms the high side fidelity of roe deer, but also shows the existence of admixtures of genes among different areas. We found evidence that genetic factors, including individual heterozygosity, influence body mass, confirming that heterozygosity positively affects fitness in wild populations. However, as the effect of genetic factors is usually masked or overruled by the influence of environmental factors, i.e., availability of resources, data on the joint influence of external and intrinsic factors on fitness and other life-history traits are needed to better predict the population dynamics of targeted species, which would enable sustainable, science-based population management.

Abstract

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.

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (N_e ), despite very large census sizes in modern C. capreolus populations and show that low N_e has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.

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.

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 diversity and genetic structure of the Siberian roe deer (Capreolus pygargus) populations from Asia

Background

The roe deer, Capreolus sp., is one of the most widespread meso-mammals of Palearctic distribution, and includes two species, the European roe deer, C. capreolus inhabiting mainly Europe, and the Siberian roe deer, C. pygargus, distributed throughout continental Asia. Although there are a number of genetic studies concerning European roe deer, the Siberian roe deer has been studied less, and none of these studies use microsatellite markers. Natural processes have led to genetic structuring in wild populations. To understand how these factors have affected genetic structure and connectivity of Siberian roe deer, we investigated variability at 12 microsatellite loci for Siberian roe deer from ten localities in Asia.

Results

Moderate levels of genetic diversity (H_E = 0.522 to 0.628) were found in all populations except in Jeju Island, South Korea, where the diversity was lowest (H_E = 0.386). Western populations showed relatively low genetic diversity and higher degrees of genetic differentiation compared with eastern populations (mean Ar = 3.54 (east), 2.81 (west), mean F_ST = 0.122). Bayesian-based clustering analysis revealed the existence of three genetically distinct groups (clusters) for Siberian roe deer, which comprise of the Southeastern group (Mainland Korea, Russian Far East, Trans-Baikal region and Northern part of Mongolia), Northwestern group (Western Siberia and Ural in Russia) and Jeju Island population. Genetic analyses including AMOVA (F_RT = 0.200), Barrier and PCA also supported genetic differentiation among regions separated primarily by major mountain ridges, suggesting that mountains played a role in the genetic differentiation of Siberian roe deer. On the other hand, genetic evidence also suggests an ongoing migration that may facilitate genetic admixture at the border areas between two groups.

Conclusions

Our results reveal an apparent pattern of genetic differentiation among populations inhabiting Asia, showing moderate levels of genetic diversity with an east-west gradient. The results suggest at least three distinct management units of roe deer in continental Asia, although genetic admixture is evident in some border areas. The insights obtained from this study shed light on management of Siberian roe deer in Asia and may be applied in conservation of local populations of Siberian roe deer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0244-6) contains supplementary material, which is available to authorized users.

Weak population structure in European roe deer (Capreolus capreolus) and evidence of introgressive hybridization with Siberian roe deer (C. pygargus) in northeastern Poland

We investigated contemporary and historical influences on the pattern of genetic diversity of European roe deer (Capreolus capreolus). The study was conducted in northeastern Poland, a zone where vast areas of primeval forests are conserved and where the European roe deer was never driven to extinction. A total of 319 unique samples collected in three sampling areas were genotyped at 16 microsatellites and one fragment (610 bp) of mitochondrial DNA (mtDNA) control region. Genetic diversity was high, and a low degree of genetic differentiation among sampling areas was observed with both microsatellites and mtDNA. No evidence of genetic differentiation between roe deer inhabiting open fields and forested areas was found, indicating that the ability of the species to exploit these contrasting environments might be the result of its phenotypic plasticity. Half of the studied individuals carried an mtDNA haplotype that did not belong to C. capreolus, but to a related species that does not occur naturally in the area, the Siberian roe deer (C. pygargus). No differentiation between individuals with Siberian and European mtDNA haplotypes was detected at microsatellite loci. Introgression of mtDNA of Siberian roe deer into the genome of European roe deer has recently been detected in eastern Europe. Such introgression might be caused by human-mediated translocations of Siberian roe deer within the range of European roe deer or by natural hybridization between these species in the past.

Influence of Holocene environmental change and anthropogenic impact on the diversity and distribution of roe deer

Extant patterns of population structure and levels of diversity are a consequence of factors that vary in both space and time. Our objective in this study is to investigate a species that has responded to both natural and anthropogenic changes in ways that have shaped modern populations and provide insight into the key processes. The roe deer (Capreolus capreolus) is one of the two species of deer native to Britain. During the last glacial maximum (LGM), the British habitat was largely under ice and there was a land bridge to mainland Europe. As the Earth warmed during the early Holocene, the land bridge was lost. Subsequent hunting on the British mainland left the southern region extirpated of roe deer, whereas a refugial population remained in the north. Later reintroductions from Europe led to population expansion, especially in southern United Kingdom. Here, we combine data from ancient and modern DNA to track population dynamics and patterns of connectivity, and test hypotheses about the influence of natural and anthropogenic environmental change. We find that past expansion and divergence events coincided with a warming environment and the subsequent closure of the land bridge between Europe and the United Kingdom. We also find turnover in British roe deer haplotypes between the late-Holocene and modern day that have likely resulted from recent human disturbance activities such as habitat perturbation, overhunting and restocking.