The genetic impact of translocations and habitat fragmentation in chamois (Rupicapra) spp

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 Diversity of Wisent Bison bonasus Based on STR Loci Analyzed in a Large Set of Samples

Wisent Bison bonasus is an example of a species saved from extinction and reintroduced into nature after a few decades of captive breeding. There were only twelve founders of the Lowland–Caucasian line (LC) and even fewer (seven out of twelve) of Lowland (LB) animals. The genetic diversity in studies based on pedigree or markers is very low. In this paper, we present a summary of the long-term genetic monitoring conducted for the worldwide population of European bison. We summarized the long-term genetic monitoring studies conducted on the worldwide population of wisents to date. We genotyped 2227 wisents from two genetic lines (LC and LB) and different populations at ten microsatellite loci. We found low polymorphism, with only 2.7 alleles per locus, and much lower values of observed heterozygosity (0.380 and 0.348 in the LC and LB lines, respectively) than expected heterozygosity. The difference between the lines is only noticeable in allele proportions, so the number of markers is not enough to distinguish the two genetic lines. We also present the genetic distance among four free-roaming populations that are geographically close to each other. We found that the genetic distance of one of them is larger than that of the others, which could be the effect of genetic drift.

Traces of past reintroduction in genetic diversity: The case of the Balkan chamois (Mammalia, Artiodactyla)

The translocation of wild animal species became a common practice worldwide to re-establish local populations threatened with extinction. Archaeological data confirm that chamois once lived in the Biokovo Mountain but, prior to their reintroduction in the 1960s, there was no written evidence of their recent existence in the area. The population was reintroduced in the period 1964-1969, when 48 individuals of Balkan chamois from the neighbouring mountains in Bosnia and Herzegovina were released. The main objective of this study was to determine the accuracy of the existing historical data on the origin of the Balkan chamois population from the Biokovo Mountain and to assess the genetic diversity and population structure of the source and translocated populations 56 years after reintroduction. Sixteen microsatellite loci were used to analyse the genetic structure of three source chamois populations from Prenj, Čvrsnica and Čabulja Mountains and from Biokovo Mountain. Both STRUCTURE and GENELAND analyses showed a clear separation of the reintroduced population on Biokovo from Prenj's chamois and considerable genetic similarity between the Biokovo population and the Čvrsnica-Čabulja population. This suggests that the current genetic composition of the Biokovo population does not derive exclusively from Prenj, as suggested by the available literature and personal interviews, but also from Čvrsnica and Čabulja. GENELAND analysis recognised the Balkan chamois from Prenj as a separate cluster, distinct from the populations of Čvrsnica and Čabulja. Our results thus highlight the need to implement genetic monitoring of both reintroduced and source populations of endangered Balkan chamois to inform sustainable management and conservation strategies in order to maximise the chances of population persistence.

Phylogenetic Analysis of the Complete Mitochondrial Genomes in the Ten Rupicapra Subspecies and Implications for the Existence of Multiple Glacial Refugia in Europe

Simple Summary

The successive glaciations that took place during the Pleistocene shaped de distribution of temperate species in Europe, both plants and animals. Traditionally, it has been hypothesized that during the coldest periods, species took refuge in three areas of southern Europe (Iberian, Italian and Balkan peninsulas) and then recolonized the north when temperatures rose and conditions were more favorable. In the present work, the complete mitochondrial sequences of the ten described chamois populations have been analyzed to identify which areas of the continent have served as refuges for the species during the glaciations. The results of the present work are consistent with the existence of multiple glacial refugia across Europe, revealing a much more complex picture of the effect of glaciations on the genetics of temperate species than is commonly accepted, giving us a better understanding of how past events determine the present species.

Abstract

The current distribution of populations in Europe is marked by the effects of glaciations that occurred during the Pleistocene. Temperate species were isolated in glacial refugia that were the sources of postglacial recolonization. The traditional glacial refuge areas were the Iberian, the Italian and the Balkan peninsulas. Here we revisit the evolutionary history of chamois (Rupicapra genus) to evaluate other sites in continental Europe and Anatolia that have been suggested as potential refuges. We have obtained the complete mitochondrial sequence of seven chamois, including the subspecies parva, carpatica, caucasica, and asiatica whose mitochondrial genome had not been yet reported. These, together with the other fourteen sequences already in the GenBank, represent the different geographical populations of the Rupicapra genus. The phylogenetic analysis showed the three old clades, dating from the early Pleistocene, already reported: mtW in the Iberian Peninsula, mtC in the Appenines and the Massif of Chartreuse, and mtE comprising all the population from the Alps to the east. The genomes within each of the clades mtW and mtE, showed divergence times larger than 300 thousand years. From here, it can be argued that the present-day lineages across Europe are very old and their split dates back to the middle Pleistocene.

Microsatellite based assignment reveals history of extirpated mountain ungulate

During the early 1900s, Northern chamois (Rupicapra rupicapra) populations in the northern Dinaric Mountains were extirpated. During the 1960s and 1970s there were several reintroductions of individuals from two Northern chamois subspecies (Alpine chamois, R. r. rupicapra and Balkan chamois, R. r. balcanica) from neighbouring areas in the attempt to re-establish the population. Accurate taxonomic classification, at subspecies level, of the autochthonous extirpated population was not known. To clarify which subspecies was present before reintroduction, we genotyped four male chamois skulls originating from Velebit Mountain, collected around 25 years before the population local extinction. DNA was successfully extracted from middle layer and outer sheath of horns. Assignment based on microsatellite loci, using both Bayesian clustering in STRUCTURE (with q values between 0.55 and 0.73) and DAPC (with individual membership probabilities of 0.99 and 1.00) indicated higher assessed likelihood for the Alpine subspecies.

Genetic and spatial characterization of the red fox (Vulpes vulpes) population in the area stretching between the Eastern and Dinaric Alps and its relationship with rabies and canine distemper dynamics

Information on the population dynamics of a reservoir species have been increasingly adopted to understand and eventually predict the dispersal patterns of infectious diseases throughout an area. Although potentially relevant, to date there are no studies which have investigated the genetic structure of the red fox population in relation to infectious disease dynamics. Therefore, we genetically and spatially characterised the red fox population in the area stretching between the Eastern and Dinaric Alps, which has been affected by both distemper and rabies at different time intervals. Red foxes collected from north-eastern Italy, Austria, Slovenia and Croatia between 2006–2012, were studied using a set of 21 microsatellite markers. We confirmed a weak genetic differentiation within the fox population using Bayesian clustering analyses, and we were able to differentiate the fox population into geographically segregated groups. Our finding might be due to the presence of geographical barriers that have likely influenced the distribution of the fox population, limiting in turn gene flow and spread of infectious diseases. Focusing on the Italian red fox population, we observed interesting variations in the prevalence of both diseases among distinct fox clusters, with the previously identified Italy 1 and Italy 2 rabies as well as distemper viruses preferentially affecting different sub-groups identified in the study. Knowledge of the regional-scale population structure can improve understanding of the epidemiology and spread of diseases. Our study paves the way for an integrated approach for disease control coupling pathogen, host and environmental data to inform targeted control programs in the future.

Multilocus Intron Trees Reveal Extensive Male-Biased Homogenization of Ancient Populations of Chamois (Rupicapra spp.) across Europe during Late Pleistocene

The inferred phylogenetic relationships between organisms often depend on the molecular marker studied due to the diverse evolutionary mode and unlike evolutionary histories of different parts of the genome. Previous studies have shown conflicting patterns of differentiation of mtDNA and several nuclear markers in chamois (genus Rupicapra) that indicate a complex evolutionary picture. Chamois are mountain caprine that inhabit most of the medium to high altitude mountain ranges of southern Eurasia. The most accepted taxonomical classification considers two species, R. pyrenaica (with the subspecies parva, pyrenaica and ornata) from southwestern Europe and R. rupicapra (with the subspecies cartusiana, rupicapra, tatrica, carpatica, balcanica, asiatica and caucasica) from northeastern Europe. Phylogenies of mtDNA revealed three very old clades (from the early Pleistocene, 1.9 Mya) with a clear geographical signal. Here we analyze a set of 23 autosomal introns, comprising 15,411 nucleotides, in 14 individuals covering the 10 chamois subspecies. Introns offered an evolutionary scenario that contrasts with mtDNA. The nucleotidic diversity was 0.0013± 0.0002, at the low range of what is found in other mammals even if a single species is considered. A coalescent multilocus analysis with *BEAST indicated that introns diversified 88 Kya, in the late Pleistocene, and the effective population size at the root was lower than 10,000 individuals. The dispersal of some few migrant males should have rapidly spread trough the populations of chamois, given the homogeneity of intron sequences. The striking differences between mitochondrial and nuclear markers can be attributed to strong female philopatry and extensive male dispersal. Our results highlight the need of analyzing multiple and varied genome components to capture the complex evolutionary history of organisms.

Genetic and ecological data reveal species boundaries between viviparous and oviparous lizard lineages

Identification of cryptic species is an essential aim for conservation biologists to avoid premature extinctions of 'unrecognized' species. Integrating different types of data can undoubtedly aid in resolving the issue of species delimitation. We studied here two lineages of the common lizard Zootoca vivipara that display different reproductive mode (the viviparous Z. v. vivipara and the oviparous Z. v. carniolica) and that overlap their distributional ranges in the European Alps. With the purpose of delimiting species' boundaries, we analyzed their ecological, genetic and natural history features. More than 300 samples were collected and analyzed at cytochrome b and 11 microsatellites loci for investigating genetic variation, population structure, individual relatedness and evolutionary histories of the two lineages. Additionally, we compared their ecological niches using eight ecological variables. Genetic data showed contrasting patterns of genetic structure between the two lineages, different demographic dynamics and no hybridization events. Also strong ecological differences (such as temperature) emerged between the two lineages, and niche overlap was limited. Taken together, these results indicate that Z. v. vivipara and Z. v. carniolica should be recognized as two separate species, and particular conservation consideration should be given to the oviparous lineage that tends to live in areas threatened by increasing impact of human activities. However, recent and rapid climate warming might determine an increasing risk for the persistence of the viviparous lineage, being adapted to cold environments.

IMMOBILIZING THE VULNERABLE APENNINE CHAMOIS (RUPICAPRA PYRENAICA ORNATA) WITH A LOW-DOSE XYLAZINE-KETAMINE COMBINATION, REVERSED WITH IDAZOXAN OR ATIPAMEZOLE

Little information is available on chemical capture of the vulnerable subspecies within the genus Rupicapra. Low-dosage combinations of xylazine and ketamine were tested for immobilization of captive and free-ranging Apennine chamois, Rupicapra pyrenaica ornata (85 and 66 immobilizations, respectively) in a retrospective analysis. Of the six dosage groups, all of them providing an acceptable level of immobilization, the optimal trade-off between safety and efficacy was found following administration of a mean dosage of 0.24±0.03 mg/kg xylazine and 1.07±0.15 mg/kg ketamine, resulting in 7.50±3.31 min induction time, deep sedation with no or limited reaction to handling in 96% of the chamois, minimal deviation of physiologic parameters from previously reported physiologic values for anesthetized or physically restrained chamois, and no mortality. Intravenous injection of idazoxan (0.05±0.01 mg/kg) or atipamezole (0.38±0.37 mg/kg) resulted in faster reversal than intravenous injection of tolazoline (1.05±0.15 mg/kg) in 1.3 vs. 4.1 min. When free-ranging chamois were darted with similar xylazine and ketamine dosages, induction time was 8.49±5.48 min, 88% of the animals were deeply sedated, and a single animal died from respiratory arrest (1.5% mortality). Intramuscular atipamezole provided smoother reversal than intravenous idazoxan. The results of this study suggest that xylazine/ketamine combinations, at remarkably lower dosage than previously published in Caprinae, may be safely and effectively used in chemical capture protocols of Apennine chamois, to facilitate conservation-oriented relocation and research.

The shared mitochondrial genome of Rupicapra pyrenaica ornata and Rupicapra rupicapra cartusiana: old remains of a common past

Mitochondrial DNA (mtDNA) has largely been used for species delimitation. However, mtDNA introgression across species boundaries can lead to inconsistent phylogenies. Partial sequences of the mitochondrial genome in the chamois, genus Rupicapra, show the presence of three well differentiated clades, West (mtW), Central (mtC) and East (mtE), each with a geographically restricted distribution. The complete mtDNAs of the clades mtW and mtE (main representatives of the two currently considered species R. pyrenaica and R. rupicapra respectively) have been reported. In the present study, we sequenced the clade mtC present in populations from both species inhabiting the central area of Europe: the Apennines (R. pyrenaica ornata) and the Chartreuse Mountains (R. rupicapra cartusiana). The phylogenetic comparison with the genomes of Caprini highlights the ancient presence of chamois in Europe relative to the fossil record, and the old age of the chamois clade mtC that was split from the clade mtW in the early Pleistocene. The separation of R. pyrenaica ornata and R. rupicapra cartusiana female lineages was recent, dating of the late Pleistocene. Our data represent an example of mtDNA introgression of resident females of Chartreuse Mountains into immigrant males of R. rupicapra due to male-biased migration and female phylopatry.

Evolution of the melanocortin-1 receptor gene (MC1R) in chamois (Rupicapra spp.)

The taxonomy of chamois and the effects of historical and evolutionary events on its diversification are still under discussion given that different morphological and genetic features presented partially discordant views. One of the morphological features that differentiate the two currently considered species, Rupicapra pyrenaica (southern chamois) and R. rupicapra (northern chamois) is coat color pattern. The melanocortin-1 receptor gene (MC1R) is related with differences in coloration in different mammals and was analyzed here in a sample of 25 chamois covering the 10 subspecies recognized, three in R. pyrenaica, (parva, pyrenaica and ornata) and seven in R. rupicapra (cartusiana, rupicapra, tatrica, carpatica, balcanica, asiatica and caucasica). Comparison with other caprinae showed that the MC1R gene has evolved under strong purifying selection. Three well differentiated haplotypes were identified: one shared by the seven subspecies of R. rupicapra, other common to the two Iberian chamois, both of the species R. pyrenaica, and a third haplotype, basal in the phylogenetic tree, unique to the subspecies from the Apennines, R. pyrenaica ornata. This pattern of variation, with three conspicuous clades, concurs with previous findings on microsatellites and mtDNA and argues in favor of the old classifications that distinguished the species R. ornata.

Polymorphisms at MHC class II DRB1 exon 2 locus in Pyrenean chamois (Rupicapra pyrenaica pyrenaica)

Chamois (Rupicapra spp.) are mountain ungulates from Southern and Central Europe and the Near East. A newly reported border disease virus (BDV) has affected the easternmost populations of Pyrenean chamois, leading to a dramatic population decrease that may drive to genetic variability loss. The Major Histocompatibility Complex (MHC) is a sensitive marker for genetic variation of populations: polymorphism on the MHC genes is affected both by pathogens and population dynamics and it is ecologically relevant, as depending on host-pathogen relationships and life history features. In the present study MHC class II DRB1 exon 2 variation was investigated in 81 Pyrenean chamois (Rupicapra pyrenaica pyrenaica) belonging to four populations. Haplotype analysis, population genetics statistics and network analysis were carried out, in order to analyze variability, phylogeography and genealogy, and the effects of geography and demographic trend. Twenty-nine haplotypes were identified, 26 of them newly described, with high Gene diversity (Gd). The variability observed in the easternmost populations of Pyrenean chamois showed a higher genetic diversity than that previously reported for other populations of Pyrenean and Cantabrian chamois (Rupicapra pyrenaica parva). The most frequent allele was RupyDRB*15, previously undetected, which seems to play a significant role in genotyping the variability, suggesting a possible effect of positive selection.

Y-chromosome phylogeny in the evolutionary net of chamois (genus Rupicapra)

Background

The chamois, distributed over most of the medium to high altitude mountain ranges of southern Eurasia, provides an excellent model for exploring the effects of historical and evolutionary events on diversification. Populations have been grouped into two species, Rupicapra pyrenaica from southwestern Europe and R. rupicapra from eastern Europe. The study of matrilineal mitochondrial DNA (mtDNA) and biparentally inherited microsatellites showed that the two species are paraphyletic and indicated alternate events of population contraction and dispersal-hybridization in the diversification of chamois. Here we investigate the pattern of variation of the Y-chromosome to obtain information on the patrilineal phylogenetic position of the genus Rupicapra and on the male-specific dispersal of chamois across Europe.

Results

We analyzed the Y-chromosome of 87 males covering the distribution range of the Rupicapra genus. We sequenced a fragment of the SRY gene promoter and characterized the male specific microsatellites UMN2303 and SRYM18. The SRY promoter sequences of two samples of Barbary sheep (Ammotragus lervia) were also determined and compared with the sequences of Bovidae available in the GenBank. Phylogenetic analysis of the alignment showed the clustering of Rupicapra with Capra and the Ammotragus sequence obtained in this study, different from the previously reported sequence of Ammotragus which groups with Ovis. Within Rupicapra, the combined data define 10 Y-chromosome haplotypes forming two haplogroups, which concur with taxonomic classification, instead of the three clades formed for mtDNA and nuclear microsatellites. The variation shows a west-to-east geographical cline of ancestral to derived alleles.

Conclusions

The phylogeny of the SRY-promoter shows an association between Rupicapra and Capra. The position of Ammotragus needs a reinvestigation. The study of ancestral and derived characters in the Y-chromosome suggests that, contrary to the presumed Asian origin, the paternal lineage of chamois originated in the Mediterranean, most probably in the Iberian Peninsula, and dispersed eastwards through serial funding events during the glacial-interglacial cycles of the Quaternary. The diversity of Y-chromosomes in chamois is very low. The differences in patterns of variation among Y-chromosome, mtDNA and biparental microsatellites reflect the evolutionary characteristics of the different markers as well as the effects of sex-biased dispersal and species phylogeography.

Effects of human perturbation on the genetic make-up of an island population: the case of the Sardinian wild boar

Game species are often manipulated by human beings, whose activities can deeply affect their genetic make-up and population structure. We focused on a geographically isolated wild boar population (Sardinia, Italy), which is classified, together with the Corsican population, as a separate subspecies (Sus scrofa meridionalis). Two hundred and ten wild boars collected across Sardinia were analysed with a set of 10 microsatellites and compared with 296 reference genotypes from continental wild populations and to a sample of domestic pigs. The Sardinian population showed remarkable diversity and a high proportion of private alleles, and strongly deviated from the equilibrium. A Bayesian cluster analysis of only the Sardinian sample revealed a partition into five subpopulations. However, two different Bayesian approaches to the assignment of individuals, accounting for different possible source populations, produced consistent results and proved the admixed nature of the Sardinian population. Indeed, introgressive hybridization with boars from multiple sources (Italian peninsula, central Europe, domestic stocks) was detected, although poor evidence of crossbreeding with free-ranging domestic pigs was unexpectedly found. After excluding individuals who carried exotic genes, the population re-entered Hardy-Weinberg proportions and a clear population structure with three subpopulations emerged. Therefore, the inclusion of introgressed animals in the Bayesian analysis implied an overestimation of the number of clusters. Nonetheless, two of them were consistent between analyses and corresponded to highly pure stocks, located, respectively, in north-west and south-west Sardinia. This work shows the critical importance of including adequate reference samples when studying the genetic structure of managed wild populations.

A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep

The development of genomic resources for wild species is still in its infancy. However, cross-species utilization of technologies developed for their domestic counterparts has the potential to unlock the genomes of organisms that currently lack genomic resources. Here, we apply the OvineSNP50 BeadChip, developed for domestic sheep, to two related wild ungulate species: the bighorn sheep (Ovis canadensis) and the thinhorn sheep (Ovis dalli). Over 95% of the domestic sheep markers were successfully genotyped in a sample of fifty-two bighorn sheep while over 90% were genotyped in two thinhorn sheep. Pooling the results from both species identified 868 single-nucleotide polymorphisms (SNPs), 570 were detected in bighorn sheep, while 330 SNPs were identified in thinhorn sheep. The total panel of SNPs was able to discriminate between the two species, assign population of origin for bighorn sheep and detect known relationship classes within one population of bighorn sheep. Using an informative subset of these SNPs (n=308), we examined the extent of genome-wide linkage disequilibrium (LD) within one population of bighorn sheep and found that high levels of LD persist over 4 Mb.

Integrating phylogeographic patterns of microsatellite and mtDNA divergence to infer the evolutionary history of chamois (genus Rupicapra)

Background

The chamois, distributed over most of the medium to high altitude mountain ranges of southern Eurasia, provides an excellent model for exploring the effects of historical and evolutionary events on diversification. Populations have been grouped into two species, Rupicapra pyrenaica from southwestern Europe and R. rupicapra from eastern Europe. However, a previous study of cytochrome b revealed that the two proposed species were non-monophyletic. The reconstruction of phylogenetic relationships between animal species often depends on the markers studied. To further elucidate the evolutionary history of chamois, we extended earlier studies by analysing DNA sequences of four mitochondrial regions (ND1, 12S, tRNApro and Control Region) and microsatellites (20 loci) to include all subspecies and cover its entire distribution range.

Results

We found discordant microsatellite (μsat) and mitochondrial (mt) DNA phylogenies. Mitochondrial phylogenies form three clades, West, Central and East (mtW, mtC and mtE), at variance with taxonomic classification. Our divergence age estimates indicate an initial separation into branches mtW-mtC and mtE 1.7 million years ago (mya), in the late Pliocene-early Pleistocene, quickly followed by the split of clades mtW and mtC. Clade mtW contains haplotypes from the Iberian peninsula and the western Alps, Clade mtC includes haplotypes from the Apennines and the Massif of Chartreuse and Clade mtE comprises populations to the east of the Alps. Divergence among populations within these three major clades is recent (< 0.5 mya). New microsatellite multilocus genotypes added to previously published data revealed differences between every pair of subspecies, forming three well defined groups (μsatW, μsatC and μsatE) also with a strong geographic signature. Grouping does not correspond with the mitochondrial lineages but is closer to morphology and taxonomic classification. Recent drastic reductions in population size can be noted for the subspecies ornata as an extremely low diversity.

Conclusions

The phylogeographic patterns for mtDNA and microsatellites suggest an evolutionary history with limited range contractions and expansions during the Quaternary period and reflect a major effect of the Alpine barrier on west-east differentiation. The contrasting phylogenies for mtDNA and microsatellites indicate events of hybridization among highly divergent lineages in the central area of distribution. Our study points to the importance of reticulate evolution, with periods of isolation and reduction of population size followed by expansions and hybridizations, in the diversification at the level of close species or subspecies.

ABC as a flexible framework to estimate demography over space and time: some cons, many pros

The analysis of genetic variation to estimate demographic and historical parameters and to quantitatively compare alternative scenarios recently gained a powerful and flexible approach: the Approximate Bayesian Computation (ABC). The likelihood functions does not need to be theoretically specified, but posterior distributions can be approximated by simulation even assuming very complex population models including both natural and human-induced processes. Prior information can be easily incorporated and the quality of the results can be analysed with rather limited additional effort. ABC is not a statistical analysis per se, but rather a statistical framework and any specific application is a sort of hybrid between a simulation and a data-analysis study. Complete software packages performing the necessary steps under a set of models and for specific genetic markers are already available, but the flexibility of the method is better exploited combining different programs. Many questions relevant in ecology can be addressed using ABC, but adequate amount of time should be dedicated to decide among alternative options and to evaluate the results. In this paper we will describe and critically comment on the different steps of an ABC analysis, analyse some of the published applications of ABC and provide user guidelines.