Projection of current and future distribution of adaptive genetic units in an alpine ungulate.

Climate projections predict major changes in alpine environments by the end of the 21st century. To avoid climate-induced maladaptation and extinction, many animal populations will either need to move to more suitable habitats or adapt in situ to novel conditions. Since populations of a species exhibit genetic variation related to local adaptation, it is important to incorporate this variation into predictive models to help assess the ability of the species to survive climate change. Here, we evaluate how the adaptive genetic variation of a mountain ungulate-the Northern chamois (Rupicapra rupicapra)-could be impacted by future global warming. Based on genotype-environment association analyses of 429 chamois using a ddRAD sequencing approach, we identified genetic variation associated with climatic gradients across the European Alps. We then delineated adaptive genetic units and projected the optimal distribution of these adaptive groups in the future. Our results suggest the presence of local adaptation to climate in Northern chamois with similar genetic adaptive responses in geographically distant but climatically similar populations. Furthermore, our results predict that future climatic changes will modify the Northern chamois adaptive landscape considerably, with various degrees of maladaptation risk.

Whole Genome Sequencing of Escherichia coli and Enterococcus spp. in wildlife-livestock interface: a pilot study.

OBJECTIVES: This pilot study provides a multidisciplinary investigation to monitor livestock-wildlife interface. Ecological data, microbiological investigations, and whole genome sequencing were used to characterize eight bacterial isolates obtained from sympatric domestic and wild ruminants in Maiella National Park (Italy) in terms of genetic patterns of antimicrobial resistance. METHODS: Using selective culturing of fresh fecal samples of monitored and georeferenced populations of Apennine chamois, goats, red deer, and sheep, Escherichia coli, Enterococcus faecium, and Enterococcus faecalis isolates were isolated and subjected to minimum inhibitory concentration determination and whole genome sequencing. RESULTS: The analyzed isolates showed phenotypic and genotypic resistance to tetracycline and critically important antibiotics such as linezolid and carbapenems. Virulence genes related to biofilm regulation and Shiga toxins were also detected. Furthermore, serotypes related to nosocomial infections, harbouring plasmids recognized as important mobile resistance gene transmitters, were identified. CONCLUSIONS: This multidisciplinary pilot study represents a promising initial step to identify the environmental drivers and the transmission routes of antimicrobial resistance and virulence factors, providing new data on bacteria from rare and endangered species such as Apennine chamois.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification.

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

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.

Low sequence diversity of the prion protein gene (PRNP) in wild deer and goat species from Spain.

The first European cases of chronic wasting disease (CWD) in free-ranging reindeer and wild elk were confirmed in Norway in 2016 highlighting the urgent need to understand transmissible spongiform encephalopathies (TSEs) in the context of European deer species and the many individual populations throughout the European continent. The genetics of the prion protein gene (PRNP) are crucial in determining the relative susceptibility to TSEs. To establish PRNP gene sequence diversity for free-ranging ruminants in the Northeast of Spain, the open reading frame was sequenced in over 350 samples from five species: Iberian red deer (Cervus elaphus hispanicus), roe deer (Capreolus capreolus), fallow deer (Dama dama), Iberian wild goat (Capra pyrenaica hispanica) and Pyrenean chamois (Rupicapra p. pyrenaica). Three single nucleotide polymorphisms (SNPs) were found in red deer: a silent mutation at codon 136, and amino acid changes T98A and Q226E. Pyrenean chamois revealed a silent SNP at codon 38 and an allele with a single octapeptide-repeat deletion. No polymorphisms were found in roe deer, fallow deer and Iberian wild goat. This apparently low variability of the PRNP coding region sequences of four major species in Spain resembles previous findings for wild mammals, but implies that larger surveys will be necessary to find novel, low frequency PRNP gene alleles that may be utilized in CWD risk control.

A new hybrid approach for MHC genotyping: high-throughput NGS and long read MinION nanopore sequencing, with application to the non-model vertebrate Alpine chamois (Rupicapra rupicapra).

The major histocompatibility complex (MHC) acts as an interface between the immune system and infectious diseases. Accurate characterization and genotyping of the extremely variable MHC loci are challenging especially without a reference sequence. We designed a combination of long-range PCR, Illumina short-reads, and Oxford Nanopore MinION long-reads approaches to capture the genetic variation of the MHC II DRB locus in an Italian population of the Alpine chamois (Rupicapra rupicapra). We utilized long-range PCR to generate a 9 Kb fragment of the DRB locus. Amplicons from six different individuals were fragmented, tagged, and simultaneously sequenced with Illumina MiSeq. One of these amplicons was sequenced with the MinION device, which produced long reads covering the entire amplified fragment. A pipeline that combines short and long reads resolved several short tandem repeats and homopolymers and produced a de novo reference, which was then used to map and genotype the short reads from all individuals. The assembled DRB locus showed a high level of polymorphism and the presence of a recombination breakpoint. Our results suggest that an amplicon-based NGS approach coupled with single-molecule MinION nanopore sequencing can efficiently achieve both the assembly and the genotyping of complex genomic regions in multiple individuals in the absence of a reference sequence.

Ancient mitochondrial pseudogenes reveal hybridization between distant lineages in the evolution of the Rupicapra genus.

Mitochondrial pseudogenes (numts) inserted in the nuclear genome are frequently found in population studies. Its presence is commonly connected with problems and errors when they are confounded with true mitochondrial sequences. In the opposite side, numts can provide valuable phylogenetic information when they are copies of ancient mitochondrial lineages. We show that Rupicapra individuals of different geographic origin from the Cantabrian Mountains to the Apennines and the Caucasus share a nuclear COI fragment. The numt copies are monophyletic, and their pattern of differentiation shows two outstanding features: a long evolution as differentiated true mitochondrial lineage, and a recent integration and spread through the chamois populations. The COI pseudogene is much older than the present day mitochondrial clades of Rupicapra and occupies a basal position within the Rupicapra-Ammotragus-Arabitragus node. Joint analysis of this numt and a cytb pseudogene with a similar pattern of evolution places the source mitochondrial lineage as a sister branch that separated from the Ammotragus-Arabitragus lineage 6millionyearsago (Mya). The occurrence of this sequence in the nucleus of chamois suggests hybridization between highly divergent lineages. The integration event seems to be very recent, more recent than the split of the present day mtDNA lineages of Rupicapra (1.9Mya). This observation invites to think of the spread across the genus by horizontal transfer through recent male-biased dispersal.

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.

Upscaling the niche variation hypothesis from the intra- to the inter-specific level.

The "niche variation hypothesis" (NVH) predicts that populations with wider niches should display higher among-individual variability. This prediction originally stated at the intra-specific level may be extended to the inter-specific level: individuals of generalist species may differ to a greater extent than individuals of a specialist species. We tested the NVH at intra- and inter-specific levels based on a large diet database of three large herbivore feces collected in the field and analyzed using DNA metabarcoding. The three herbivores (roe deer Capreolus capreolus, chamois Rupicapra rupicapra and mouflon Ovis musimon) are highly contrasted in terms of sociality (solitary to highly gregarious) and diet. The NVH at the intraspecific level was tested by relating, for the same population, diet breadth and inter-individual variation across the four seasons. Compared to null models, our data supported the NVH both at the intra- and inter-specific levels. Inter-individual variation of the diet of solitary species was not larger than in social species, although social individuals feed together and could therefore have more similar diets. Hence, the NVH better explained diet breadth than other factors such as sociality. The expansion of the population niche of the three species was driven by resource availability, and achieved by an increase in inter-individual variation, and the level of inter-individual variability was larger in the generalist species (mouflon) than in the specialist one (roe deer). This mechanism at the base of the NVH appears at play at different levels of biological organization, from populations to communities.

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.

Sex-specific selection for MHC variability in Alpine chamois.

BACKGROUND: In mammals, males typically have shorter lives than females. This difference is thought to be due to behavioural traits which enhance competitive abilities, and hence male reproductive success, but impair survival. Furthermore, in many species males usually show higher parasite burden than females. Consequently, the intensity of selection for genetic factors which reduce susceptibility to pathogens may differ between sexes. High variability at the major histocompatibility complex (MHC) genes is believed to be advantageous for detecting and combating the range of infectious agents present in the environment. Increased heterozygosity at these immune genes is expected to be important for individual longevity. However, whether males in natural populations benefit more from MHC heterozygosity than females has rarely been investigated. We investigated this question in a long-term study of free-living Alpine chamois (Rupicapra rupicapra), a polygynous mountain ungulate. RESULTS: Here we show that male chamois survive significantly (P = 0.022) longer if heterozygous at the MHC class II DRB locus, whereas females do not. Improved survival of males was not a result of heterozygote advantage per se, as background heterozygosity (estimated across twelve microsatellite loci) did not change significantly with age. Furthermore, reproductively active males depleted their body fat reserves earlier than females leading to significantly impaired survival rates in this sex (P < 0.008). This sex-difference was even more pronounced in areas affected by scabies, a severe parasitosis, as reproductively active males were less likely to survive than females. However, we did not find evidence for a survival advantage associated with specific MHC alleles in areas affected by scabies. CONCLUSIONS: Increased MHC class II DRB heterozygosity with age in males, suggests that MHC heterozygous males survive longer than homozygotes. Reproductively active males appear to be less likely to survive than females most likely because of the energetic challenge of the winter rut, accompanied by earlier depletion of their body fat stores, and a generally higher parasite burden. This scenario renders the MHC-mediated immune response more important for males than for females, which implies a relatively stronger selection pressure on MHC genes in males than in females.

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.

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 (musat) 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 (musatW, musatC and musatE) 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.

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

The chamois is a useful species with which to investigate the combined genetic impact of habitat fragmentation, over hunting, and translocations. Genetic variation within and between chamois (genus Rupicapra) populations was analyzed in 259 individuals from 16 sampling sites located in Italy, Spain, Slovakia, and the Czech Republic. Two mitochondrial DNA markers (control region and cytochrome b) and 11 nuclear microsatellites were typed. The principal results of this study can be summarized as follows: 1) high and significant differentiation between almost all chamois populations is observed even on a microgeographical scale, probably caused by the patchy distribution of this species, sharp geographical barriers to gene flow, and drift effects related to recent bottlenecks; 2) historical translocation events have left a clear genetic signature, including interspecific hybridization in some Alpine localities; 3) the Apennine subspecies of chamois, Rupicapra pyrenaica ornata, shows a high and similar level of divergence (about 1.5 My) from the Pyrenean (Rupicapra pyrenaica pyrenaica) and the Alpine (Rupicapra rupicapra) chamois; therefore, the specific status of these taxa should be revised. These results confirm the potential of population genetic analyses to dissect and interpret complex patterns of diversity in order to define factors important to conservation and management.

Cytochrome b phylogeography of chamois (Rupicapra spp.). Population contractions, expansions and hybridizations governed the diversification of the genus.

The chamois provides an excellent model for exploring the effect of historical and evolutionary events on diversification. We investigate cytochrome b (cytb) sequences in the 10 recognized subspecies of Rupicapra classified within 2 species: Rupicapra pyrenaica, with the subspecies parva, pyrenaica, and ornata, and Rupicapra rupicapra, with cartusiana, rupicapra, tatrica, carpatica, balcanica, asiatica, and caucasica. A fragment of 349 bp of the cytb was sequenced in 189 individuals. We identified 3 cytb lineages: Clade West in Iberia and Western Alps; Clade Central in the Apennines and the Massif of Chartreuse; and Clade East present in populations to the east of the Alps. The 2 proposed species were polyphyletic; the clades West and Central are represented in both, whereas the Clade East is restricted to R. rupicapra. In contrast to the current systematic, cytb phylogenies suggest the classification of the 10 subspecies of chamois into a single species, R. rupicapra. Phylogeny and geographical distribution of the 3 lineages show the effects of limited latitudinal range expansions, contractions, and hybridizations among highly divergent lineages, along with a major role of the glacial ice sheets of the Alps and the Pyrenees as barriers to gene flow, on the diversification of extant taxa.

Real-time PCR for quantitative detection of chamois (Rupicapra rupicapra) and pyrenean ibex (Capra pyrenaica) in meat mixtures.

A real-time quantitative polymerase chain reaction (PCR) technique was developed for the quantification of chamois and pyrenean ibex DNAs in meat mixtures by using a SYBR green detection platform. Two species-specific systems and a eukaryotic endogenous system were combined in the real-time PCR approach to quantify the target species. In the specific systems, a 133 base pair (bp) fragment of the 12S rRNA gene was amplified from chamois DNA, and an 88 bp fragment from the D-loop region was amplified from pyrenean ibex DNA. In the endogenous system, universal primers amplified a 141 bp fragment on the nuclear 18S rRNA gene from eukaryotic DNA. The threshold cycle values obtained with the 18S rRNA primers were used to normalize those obtained from chamois- or pyrenean ibex-specific systems, serving as endogenous control for the total content of PCR-amplifiable DNA in the sample. Analysis of experimental raw and heat-treated binary mixtures of chamois and pyrenean ibex meat in a swine meat matrix demonstrated the suitability of the assay for the detection and quantification of the target DNAs in the range of 0.1-0.8%, depending on the species and treatment of the meat samples.

Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois.

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen-mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d(n)/d(s) ratio and the presence of trans-species polymorphisms suggest that a strong long-term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D-loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000-30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.

Diversity and evolution of the Mhc-DRB1 gene in the two endemic Iberian subspecies of Pyrenean chamois, Rupicapra pyrenaica.

Major histocompatibility complex class II locus DRB variation was investigated by single-strand conformation polymorphism analysis and sequence analysis in the two subspecies of Pyrenean chamois (Rupicapra pyrenaica) endemic to the Iberian Peninsula. Low levels of genetic variation were detected in both subspecies, with seven different alleles in R. p. pyrenaica and only three in the R. p. parva. After applying the rarefaction method to cope with the differences in sample size, the low allele number of parva was highlighted. The low allelic repertoire of the R. p. parva subspecies is most likely the result of bottlenecks caused by hunting pressure and recent parasitic infections by sarcoptic mange. A phylogenetic analysis of both Pyrenean chamois and DRB alleles from 10 different caprinid species revealed that the chamois alleles form two monophyletic groups. In comparison with other Caprinae DRB sequences, the Rupicapra alleles displayed a species-specific clustering that reflects a large temporal divergence of the chamois from other caprinids, as well as a possible difference in the selective environment for these species.