Greece: a Balkan subrefuge for a remnant red deer (cervus elaphus) population

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.

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.

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.

Genetic analyses of brown hare (Lepus europaeus) support limited migration and translocation of Greek populations

Numerous studies have shown that the phylogeography of many species, including European brown hare, has been affected by the climatic oscillations of the Pleistocene. During this period the Balkans acted as a major refugium offering habitable conditions for many species. However, few studies have focused on the specific role of the Greek peninsula in the phylogeographic history of species in this southernmost margin of Balkans. We, therefore analyzed a 528 bp fragment of the D-loop region of mitochondrial DNA in 154 wild brown hare individuals from unsampled areas from both mainland and island Greece and compared it to 310 available brown hare sequences (including 110 Greek samples). Newly identified haplotypes show characteristic distribution in specific Greek areas reinforcing the theory that Greece can be considered as a subrefuge within Balkans for a number of species, with several “refugia within refugia” spots, holding significant genetic diversity. No haplotypes from wild Greek individuals clustered with the Central and Western Europe group revealing a minimal contribution of this area to the colonization of central Europe. One hundred and ten reared brown hares were also analyzed to elucidate the impact of introductions on local populations. Most of these samples presented close genetic affinity with haplotypes from Central and Western Europe indicating that farms in Greece use breeders from those areas. Therefore, despite human translocation of individuals, the genetic structure of brown hare has mostly been influenced by paleoclimatic conditions and minimally by human actions.

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.

Microevolution of the noble crayfish (Astacus astacus) in the Southern Balkan Peninsula

Background

The noble crayfish (Astacus astacus) displays a complex historical and contemporary genetic status in Europe. The species divergence has been shaped by geological events (i.e. Pleistocene glaciations) and humanly induced impacts (i.e. translocations, pollution, etc.) on its populations due to species commercial value and its niche degradation. Until now, limited genetic information has been procured for the Balkan area and especially for the southernmost distribution of this species (i.e. Greece). It is well known that the rich habitat diversity of the Balkan Peninsula offers suitable conditions for genetically diversified populations. Thus, the present manuscript revisits the phylogenetic relationships of the noble crayfish in Europe and identifies the genetic make-up and the biogeographical patterns of the species in its southern range limit.

Results

Mitochondrial markers (i.e. COI and 16S) were used in order to elucidate the genetic structure and diversity of the noble crayfish in Europe. Two of the six European haplotypic lineages, were found exclusively in Greece. These two lineages exhibited greater haplotypic richness when compared with the rest four (of “Central European” origin) while they showed high genetic diversity. Divergence time analysis identified that the majority of this divergence was captured through Pleistocene, suggesting a southern glacial refugium (Greece, southern Balkans). Furthermore, six microsatellite markers were used in order to define the factors affecting the genetic structure and demographic history of the species in Greece. The population structure analysis revealed six to nine genetic clusters and eight putative genetic barriers. Evidence of bottleneck effects in the last ~5000 years (due to climatic and geological events and human activities) is also afforded. Findings from several other research fields (e.g. life sciences, geology or even archaeology) have been utilized to perceive the genetic make-up of the noble crayfish.

Conclusions

The southernmost part of Balkans has played a major role as a glacial refugium for A. astacus. Such refugia have served as centres of expansion to northern regions. Recent history of the noble crayfish in southern Balkans reveals the influence of environmental (climate, geology and/or topology) and anthropogenic factors.

Electronic supplementary material

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

Chromosomal and mitochondrial diversity in Melitaea didyma complex (Lepidoptera, Nymphalidae): eleven deeply diverged DNA barcode groups in one non-monophyletic species?

It is generally accepted that cases of species' polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, and the detected cases reflect misidentifications or/and methodological errors. Here we studied the problem of species' non-monophyly through chromosomal and molecular analysis of butterfly taxa close to Melitaea didyma (Esper, 1779) (Lepidoptera, Nymphalidae). We found absence or low interspecific chromosome number variation and presence of intraspecific variation, therefore we conclude that in this group, chromosome numbers have relatively low value as taxonomic markers. Despite low karyotype variability, the group was found to have unexpectedly high mitochondrial haplotype diversity. These haplotypes were clustered in 23 highly diverged haplogroups. Twelve of these haplogroups are associated with nine traditionally recognized and morphologically distinct species Melitaea chitralensis Moore, 1901, Melitaea deserticola Oberthür, 1909, Melitaea didymoides Eversmann, 1847, Melitaea gina Higgins, 1941, Melitaea interrupta Colenati, 1846, Melitaea latonigena Eversmann, 1847, Melitaea mixta Evans, 1912, Melitaea saxatilis Christoph, 1873 and Melitaea sutschana Staudinger, 1892. The rest of the haplogroups (11 lineages) belong to a well-known west-palaearctic species Melitaea didyma. The last species is particularly unusual in the haplotypes we obtained. First, it is clearly polyphyletic with respect to COI gene. Second, the differentiation in COI gene between these mostly allopatric (but in few cases sympatric) eleven lineages is extremely high (up to 7.4%), i.e. much deeper than the "standard" DNA barcode species threshold (2.7-3%). This level of divergence normally could correspond not even to different species, but to different genera. Despite this divergence, the bearers of these haplogroups were found to be morphologically indistinguishable and, most importantly, to share absolutely the same ecological niches, i.e. demonstrating the pattern which is hardly compatible with hypothesis of multiple cryptic species. Most likely such a profound irregularity in barcodes is caused by reasons other than speciation and represents an extraordinary example of intra-species barcode variability. Given the deep level of genetic differentiation between the lineages, we assume that there was a long period (up to 5.0 My) of allopatric differentiation when the lineages were separated by geographic or/and ecological barriers and evolved in late Pliocene and Pleistocene refugia of north Africa, the Iberian and Balkan Peninsulas, the Middle East and Central Asia. We discuss the refugia-within-refugia concept as a mechanism explaining the presence of additional diverged minor haplogroups within the areas of the major haplogroups. We also provide the first record of Melitaea gina in Azerbaijan and the record of Melitaea didyma turkestanica as a new taxon for Russia and Europe.

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.

Genetic diversity of Thaumetopoea pityocampa in Greece: the role of Quaternary changes in Aegean Sea

Thaumetopoea pityocampa, the winter pine processionary moth, is one of the most important pests of pine trees in the Mediterranean region. To learn more about the refugial areas within the greater refugium of Greece, samples from 15 local populations were collected. Analysis of the Cytochrome Oxidase I region identified 15 haplotypes. One strongly supported clade was found, separating the Aegean island of Lesvos from the mainland populations. Mdiv analysis showed that this clade diverged from the general clade 155,000 y.a., suggesting the sea level changes during the Quaternary Period as a possible reason for its seclusion. Additionally, comparisons with similar studies in this region revealed a possible colonization of the eastern Aegean islands from the Turkish mainland rather than from Greece. However, no strong isolation-by-distance events were detected among the mainland populations, which could be attributed to the joint effect of regular gene flow and the lack of insurmountable geographic barriers. Finally, regarding population structure, Bayesian analysis as well as neutrality tests pointed towards an ongoing population expansion which verifies the potential invasiveness of this pest species, something that will have to be dealt with under the influence of climate change.

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

Due to a restriction of the distributional range of European red deer (Cervus elaphus L.) during the Quaternary and subsequent recolonization of Europe from different refugia, a clear phylogeographical pattern in genetic structure has been revealed using mitochondrial DNA markers. In Central Europe, 2 distinct, eastern and western, lineages of European red deer are present; however, admixture between them has not yet been studied in detail. We used mitochondrial DNA (control region and cytochrome b gene) sequences and 22 microsatellite loci from 522 individuals to investigate the genetic diversity of red deer in what might be expected to be an intermediate zone. We discovered a high number of unique mtDNA haplotypes belonging to each lineage and high levels of genetic diversity (cyt b H = 0.867, D-loop H = 0.914). The same structuring of red deer populations was also revealed by microsatellite analysis, with results from both analyses thus suggesting a suture zone between the 2 lineages. Despite the fact that postglacial recolonization of Central Europe by red deer occurred more than 10000 years ago, the degree of admixture between the 2 lineages is relatively small, with only 10.8% admixed individuals detected. Direct translocations of animals by humans have slightly blurred the pattern in this region; however, this blurring was more apparent when using maternally inherited markers than nuclear markers.