Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Hybridization between Alpine Ibex and Domestic Goat in the Alps: A Sporadic and Localized Phenomenon?

Simple Summary

The Alpine ibex (Capra ibex) is a protected wild ungulate. The species is known to have a low genetic variability and to occasionally suffer from local population decline as a consequence of epidemic diseases. Another, less investigated, threat for the species’ long-term conservation is represented by the hybridization with feral domestic goats that may breed with ibex if abandoned in the mountain at the end of the summer pasture. Further reproduction and the spread of hybrids may jeopardize the genetic integrity of wild Alpine ibex. By means of an online survey and using a network of experts, we mapped unpublished evidence on observed hybridization events between Alpine ibex and feral domestic goats. The results of this study suggested that hybrids are currently distributed in several countries, and their presence is not a rare event, with some clusters including 4–20 probable hybrids. This calls for more studies clearly quantifying hybrids in Alpine ibex colonies, but also highlights the need for conservation actions aimed at controlling this phenomenon, such as the effective management of domestic herds grazing in Alpine ibex core areas and clear guidelines on hybrid management.

Abstract

The Alpine ibex (Capra ibex) is a mountain ungulate living in the European Alps. Although being currently classified as a species of Least Concern (LC) by the IUCN, a potential threat for its long-term conservation is introgression following hybridization with domestic goats (Capra hircus). Hybridization has been documented in Switzerland in captive and free ranging animals, although accurate data to assess the extent of this phenomenon in natural conditions in the Alps are lacking. Using an online survey and a network of experts, we collected and mapped unpublished evidence of hybridization events that occurred between Alpine ibex and feral domestic goats from 2000 to 2021. The results of this study showed that hybrids are distributed in most of the Alpine countries, and their presence is not a sporadic event, with some clusters including 4–20 probable hybrids. Our results illustrated the need for establishing a standardized and effective protocol to identify hybrids in the field (such as a formal description of the morphological traits characterizing hybrids), as well as clear guidelines for hybrid management. Even more importantly, this study also highlighted the need for actions aimed at avoiding hybridization, such as the effective management of domestic herds grazing in Alpine ibex core areas.

Assessing the levels of intraspecific admixture and interspecific hybridization in Iberian wild goats (Capra pyrenaica)

Iberian wild goats (Capra pyrenaica, also known as Iberian ibex, Spanish ibex, and Spanish wild goat) underwent strong genetic bottlenecks during the 19th and 20th centuries due to overhunting and habitat destruction. From the 1970s to 1990s, augmentation translocations were frequently carried out to restock Iberian wild goat populations (very often with hunting purposes), but they were not systematically planned or recorded. On the other hand, recent data suggest the occurrence of hybridization events between Iberian wild goats and domestic goats (Capra hircus). Augmentation translocations and interspecific hybridization might have contributed to increase the diversity of Iberian wild goats. With the aim of investigating this issue, we have genotyped 118 Iberian wild goats from Tortosa-Beceite, Sierra Nevada, Muela de Cortes, Gredos, Batuecas, and Ordesa and Monte Perdido by using the Goat SNP50 BeadChip (Illumina). The analysis of genotypic data indicated that Iberian wild goat populations are strongly differentiated and display low diversity. Only three Iberian wild goats out from 118 show genomic signatures of mixed ancestry, a result consistent with a scenario in which past augmentation translocations have had a limited impact on the diversity of Iberian wild goats. Besides, we have detected eight Iberian wild goats from Tortosa-Beceite with signs of domestic goat introgression. Although rare, hybridization with domestic goats could become a potential threat to the genetic integrity of Iberian wild goats; hence, measures should be taken to avoid the presence of uncontrolled herds of domestic or feral goats in mountainous areas inhabited by this iconic wild ungulate.

A robust sequencing assay of a thousand amplicons for the high-throughput population monitoring of Alpine ibex immunogenetics

Polymorphism for immune functions can explain significant variation in health and reproductive success within species. Drastic loss in genetic diversity at such loci constitutes an extinction risk and should be monitored in species of conservation concern. However, effective implementations of genome‐wide immune polymorphism sets into high‐throughput genotyping assays are scarce. Here, we report the design and validation of a microfluidics‐based amplicon sequencing assay to comprehensively capture genetic variation in Alpine ibex (Capra ibex). This species represents one of the most successful large mammal restorations recovering from a severely depressed census size and a massive loss in diversity at the major histocompatibility complex (MHC). We analysed 65 whole‐genome sequencing sets of the Alpine ibex and related species to select the most representative markers and to prevent primer binding failures. In total, we designed ~1,000 amplicons densely covering the MHC, further immunity‐related genes as well as randomly selected genome‐wide markers for the assessment of neutral population structure. Our analysis of 158 individuals shows that the genome‐wide markers perform equally well at resolving population structure as RAD‐sequencing or low‐coverage genome sequencing data sets. Immunity‐related loci show unexpectedly high degrees of genetic differentiation within the species. Such information can now be used to define highly targeted individual translocations. Our design strategy can be realistically implemented into genetic surveys of a large range of species. In conclusion, leveraging whole‐genome sequencing data sets to design targeted amplicon assays allows the simultaneous monitoring of multiple genetic risk factors and can be translated into species conservation recommendations.

The origin of domestication genes in goats

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

Assessing possible hybridization among managed Nubian ibex in North America

Hybridization among closely related species is a concern in zoo and aquarium populations where unpedigreed animals are frequently exchanged with the private sector. In this study, we examine possible hybridization in a group of Nubian ibex (Capra nubiana) imported into the Association of Zoos and Aquariums' (AZA) Species Survival Program (SSP) from a private institution. These individuals appeared smaller in stature than adult SSP Nubian ibex and were excluded from breeding recommendations over the concern that they were hybrids. Twenty-six microsatellites were used to rule out recent hybridization with domestic goats, Siberian ibex (Capra sibirica), and Alpine ibex (Capra ibex). We argue that natural phenotypic variation across the large geographic range of Nubian ibex may account for the small stature of the imported ibex, as private institutions may have historically acquired individuals from locations that differed from the SSP founders. However, the imported Nubian ibex appeared genetically differentiated from the SSP Nubian ibex and may represent a source of genetic variation for the managed population.

Hidden MHC genetic diversity in the Iberian ibex (Capra pyrenaica)

Background

Defining hidden genetic diversity within species is of great significance when attempting to maintain the evolutionary potential of natural populations and conduct appropriate management. Our hypothesis is that isolated (and eventually small) wild animal populations hide unexpected genetic diversity due to their maintenance of ancient polymorphisms or introgressions.

Results

We tested this hypothesis using the Iberian ibex (Capra pyrenaica) as an example. Previous studies based on large sample sizes taken from its principal populations have revealed that the Iberian ibex has a remarkably small MHC DRB1 diversity (only six remnant alleles) as a result of recent population bottlenecks and a marked demographic decline that has led to the extinction of two recognized subspecies. Extending on the geographic range to include non-studied isolated Iberian ibex populations, we sequenced a new MHC DRB1 in what seemed three small isolated populations in Southern Spain (n = 132). The findings indicate a higher genetic diversity than previously reported in this important gene. The newly discovered allele, MHC DRB1*7, is identical to one reported in the domestic goat C. aegagrus hircus. Whether or not this is the result of ancient polymorphisms maintained by balancing selection or, alternatively, introgressions from domestic goats through hybridization needs to be clarified in future studies. However, hybridization between Iberian ibex and domestic goats has been reported in Spain and the fact that the newly discovered allele is only present in one of the small isolated populations and not in the others suggests introgression. The new discovered allele is not expected to increase fitness in C. pyrenaica since it generates the same protein as the existing MHC DRB1*6. Analysis of a microsatellite locus (OLADRB1) near the new MHC DRB1*7 gene reveals a linkage disequilibrium between these two loci. The allele OLADRB1, 187 bp in length, was unambiguously linked to the MHC DRB1*7 allele. This enabled us to perform a DRB-STR matching method for the recently discovered MHC allele.

Conclusions

This finding is critical for the conservation of the Iberian ibex since it directly affects the identification of the units of this species that should be managed and conserved separately (Evolutionarily Significant Units).

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0616-9) contains supplementary material, which is available to authorized users.

Introgression from domestic goat generated variation at the major histocompatibility complex of Alpine ibex

The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.

The major histocompatibility complex (MHC), a crucial component of the defense against pathogens, contains the most polymorphic functional genes in vertebrate genomes. The extraordinary genetic variation is generally considered to be ancient. We investigated whether a previously neglected mechanism, introgression from related species, provides an additional source of MHC variation. We show that introgression from domestic goat dramatically increased genetic variation at the MHC of Alpine ibex, a species that had nearly gone extinct during the 18^th century, but has been restored to large numbers since. We show that Alpine ibex share one of only two alleles at a generally highly polymorphic MHC locus with domestic goats and that the chromosomal region containing the goat-type allele has a signature of recent introgression. Our finding contradicts the long-standing view that ancient trans-species polymorphism is the sole source of the extraordinary genetic variability at the MHC. Instead, we show that in Alpine ibex introgression generated genetic diversity at a MHC locus. Our study supports the view that loci favoring genetic polymorphism may be susceptible to adaptive introgression from related species and will encourage future research to identify unexpected signatures of introgression.

Brucella melitensis in France: persistence in wildlife and probable spillover from Alpine ibex to domestic animals

Bovine brucellosis is a major zoonosis, mainly caused by Brucella abortus, more rarely by Brucella melitensis. France has been bovine brucellosis officially-free since 2005 with no cases reported in domestic/wild ruminants since 2003. In 2012, bovine and autochthonous human cases due to B. melitensis biovar 3 (Bmel3) occurred in the French Alps. Epidemiological investigations implemented in wild and domestic ruminants evidenced a high seroprevalence (>45%) in Alpine ibex (Capra ibex); no cases were disclosed in other domestic or wild ruminants, except for one isolated case in a chamois (Rupicapra rupicapra). These results raised the question of a possible persistence/emergence of Brucella in wildlife. The purpose of this study was to assess genetic relationships among the Bmel3 strains historically isolated in humans, domestic and wild ruminants in Southeastern France, over two decades, by the MLVA-panel2B assay, and to propose a possible explanation for the origin of the recent bovine and human infections. Indeed, this genotyping strategy proved to be efficient for this microepidemiological investigation using an interpretation cut-off established for a fine-scale setting. The isolates, from the 2012 domestic/human outbreak harbored an identical genotype, confirming a recent and direct contamination from cattle to human. Interestingly, they clustered not only with isolates from wildlife in 2012, but also with local historical domestic isolates, in particular with the 1999 last bovine case in the same massif. Altogether, our results suggest that the recent bovine outbreak could have originated from the Alpine ibex population. This is the first report of a B. melitensis spillover from wildlife to domestic ruminants and the sustainability of the infection in Alpine ibex. However, this wild population, reintroduced in the 1970s in an almost closed massif, might be considered as a semi-domestic free-ranging herd. Anthropogenic factors could therefore account with the high observed intra-species prevalence.

Bovine viral diarrhea virus in free-ranging wild ruminants in Switzerland: low prevalence of infection despite regular interactions with domestic livestock

BACKGROUND

In the frame of an eradication program for bovine viral diarrhea (BVD) in Swiss livestock, the question was raised whether free-ranging wildlife could threaten the success of this sanitary measure. Therefore, we conducted serological and virological investigations on BVD virus (BVDV) infections in the four indigenous wild ruminant species (roe deer, red deer, Alpine chamois and Alpine ibex) from 2009 to 2011, and gathered information on interactions between wild and domestic ruminants in an alpine environment by questionnaire survey.

RESULTS

Thirty-two sera out of 1'877 (1.7%, 95% confidence interval [CI] 1.2-2.4) were seropositive for BVDV, and a BVDV1 sub genotype h virus was found in a seropositive chamois (0.05%, 95% CI 0.001-0.3). The seropositive animals originated from sub-alpine or alpine regions and significantly more seropositive red deer, chamois and ibex than roe deer were found. There were no statistically significant differences between sampling units, age classes, genders, and sampling years. The obtained prevalences were significantly lower than those documented in livestock, and most positive wild ruminants were found in proximity of domestic outbreaks. Additionally, BVDV seroprevalence in ibex was significantly lower than previously reported from Switzerland. The survey on interspecific interactions revealed that interactions expected to allow BVDV transmission, from physical contacts to non-simultaneous use of the same areas, regularly occur on pastures among all investigated ruminant species. Interactions involving cervids were more often observed with cattle than with small ruminants, chamois were observed with all three domestic species, and ibex interacted mostly with small ruminants. Interactions related to the use of anthropogenic food sources were frequently observed, especially between red deer and cattle in wintertime.

CONCLUSIONS

To our knowledge, this is the first report of BVDV RNA isolated from an Alpine chamois. Nevertheless, our results suggest that BVDV infections are only sporadic in Swiss wild ruminants, despite regular occurrence of interactions with potentially infected livestock. Overall, serological, virological and ethological data indicate that wildlife is currently an incidental spill-over host and not a reservoir for BVDV in Switzerland.

Applicability of major histocompatibility complex DRB1 alleles as markers to detect vertebrate hybridization: a case study from Iberian ibex × domestic goat in southern Spain

Background

Hybridization between closely related wild and domestic species is of great concern because it can alter the evolutionary integrity of the affected populations. The high allelic variability of Major Histocompatibility Complex (MHC) loci usually excludes them from being used in studies to detect hybridization events. However, if a) the parental species don’t share alleles, and b) one of the parental species possesses an exceptionally low number of alleles (to facilitate analysis), then even MHC loci have the potential to detect hybrids.

Results

By genotyping the exon2 of the MHC class II DRB1 locus, we were able to detect hybridization between domestic goats (Capra hircus) and free-ranging Iberian ibex (Capra pyrenaica hispanica) by molecular means.

Conclusions

This is the first documentation of a Capra pyrenaica × Capra hircus hybridization, which presented us the opportunity to test the applicability of MHC loci as new, simple, cost-effective, and time-saving approach to detect hybridization between wild species and their domesticated relatives, thus adding value to MHC genes role in animal conservation and management.