A strong genetic footprint of the re-introduction history of Alpine ibex (Capra ibex ibex)

Genetic mixing in conservation translocations increases diversity of a keystone threatened species, Bettongia lesueur

Translocation programmes are increasingly being informed by genetic data to monitor and enhance conservation outcomes for both natural and established populations. These data provide a window into contemporary patterns of genetic diversity, structure and relatedness that can guide managers in how to best source animals for their translocation programmes. The inclusion of historical samples, where possible, strengthens monitoring by allowing assessment of changes in genetic diversity over time and by providing a benchmark for future improvements in diversity via management practices. Here, we used reduced representation sequencing (ddRADseq) data to report on the current genetic health of three remnant and seven translocated boodie (Bettongia lesueur) populations, now extinct on the Australian mainland. In addition, we used exon capture data from seven historical mainland specimens and a subset of contemporary samples to compare pre-decline and current diversity. Both data sets showed the significant impact of population founder source (whether multiple or single) on the genetic diversity of translocated populations. Populations founded by animals from multiple sources showed significantly higher genetic diversity than the natural remnant and single-source translocation populations, and we show that by mixing the most divergent populations, exon capture heterozygosity was restored to levels close to that observed in pre-decline mainland samples. Relatedness estimates were surprisingly low across all contemporary populations and there was limited evidence of inbreeding. Our results show that a strategy of genetic mixing has led to successful conservation outcomes for the species in terms of increasing genetic diversity and provides strong rationale for mixing as a management strategy.

A centenary tale: population genetic insights into the introduction history of the oriental fire-bellied toad (Bombina orientalis) in Beijing

Background

The successful establishment of a species population following a single introduction of a few individuals to a non-native area has been limited. Nevertheless, the oriental fire-bellied toad (Bombina orientalis) population in Beijing is purportedly descended from a single introduction of about 200 individuals translocated from Yantai, Shandong Province, China, in 1927.

Results

To resolve the introduction process and to understand the genetic consequences since that introduction approximately 90 years ago, we investigated the population’s genetic diversity and structure using 261 toads from Beijing and two native Shandong populations and inferred the species’ introduction history using simulation-based approaches. Analysis of mitochondrial DNA (mtDNA) sequences showed the two haplotypes found in Beijing nested within Yantai haplotypes, thus corroborating the historical record of the translocation source. The mtDNA and 11 nuclear microsatellite markers revealed both considerably lower genetic diversity in Beijing than in the source population and strong genetic differentiation between them. Although the current census population in Beijing may be in the range of a few thousand, the effective population size was estimated at only 20–57. Simulations also suggest that this population may have descended from 40–60 founders.

Conclusions

The Beijing population’s genetic patterns were consistent with the consequences of a severe bottleneck during introduction followed by genetic drift. The introduction trajectory constructed for this B. orientalis population reveals the genetic footprints of a small population sustained in isolation for nearly a century. Our results provide an intriguing example of establishment success from limited founders and may inform ex situ conservation efforts as well as the management of biological invasions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02072-z.

Ancient mitochondrial and modern whole genomes unravel massive genetic diversity loss during near extinction of Alpine ibex

Population bottlenecks can have dramatic consequences for the health and long-term survival of a species. Understanding of historic population size and standing genetic variation prior to a contraction allows estimating the impact of a bottleneck on the species genetic diversity. Although historic population sizes can be modelled based on extant genomics, uncertainty is high for the last 10-20 millenia. Hence, integrating ancient genomes provides a powerful complement to retrace the evolution of genetic diversity through population fluctuations. Here, we recover 15 high-quality mitogenomes of the once nearly extinct Alpine ibex spanning 8601 BP to 1919 CE and combine these with 60 published modern whole genomes. Coalescent demography simulations based on modern whole genomes indicate population fluctuations coinciding with the last major glaciation period. Using our ancient and historic mitogenomes, we investigate the more recent demographic history of the species and show that mitochondrial haplotype diversity was reduced to a fifth of the pre-bottleneck diversity with several highly differentiated mitochondrial lineages having co-existed historically. The main collapse of mitochondrial diversity coincides with elevated human population growth during the last 1-2 kya. After recovery, one lineage was spread and nearly fixed across the Alps due to recolonization efforts. Our study highlights that a combined approach integrating genomic data of ancient, historic and extant populations unravels major long-term population fluctuations from the emergence of a species through its near extinction up to the recent past.

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.

Phylogeography and population genetic structure of the European roe deer in Switzerland following recent recolonization

In the early 1800s, the European roe deer (Capreolus capreolus) was probably extirpated from Switzerland, due to overhunting and deforestation. After a federal law was enacted in 1875 to protect lactating females and young, and limiting the hunting season, the roe deer successfully recovered and recolonized Switzerland. In this study, we use mitochondrial DNA and nuclear DNA markers to investigate the recolonization and assess contemporary genetic structure in relation to broad topographic features, in order to understand underlying ecological processes, inform future roe deer management strategies, and explore the opportunity for development of forensic traceability tools. The results concerning the recolonization origin support natural, multidirectional immigration from neighboring countries. We further demonstrate that there is evidence of weak genetic differentiation within Switzerland among topographic regions. Finally, we conclude that the genetic data support the recognition of a single roe deer management unit within Switzerland, within which there is a potential for broad‐scale geographic origin assignment using nuclear markers to support law enforcement.

OCCURRENCE OF FOOTROT IN FREE-RANGING ALPINE IBEX (CAPRA IBEX) COLONIES IN SWITZERLAND

Footrot is a worldwide economically important, debilitating disease caused by Dichelobacter nodosus. In sheep (Ovis aries), it is characterized by lesions of varying severity, depending on the strain, whereas Alpine ibex (Capra ibex) seem to develop severe lesions, whatever the strain. Healthy carriers occur in livestock but are rare in wild ruminants. Using a triangulation approach (retrospective questionnaire survey, necropsy database screening, and pathogen prevalence estimation in selected ibex colonies with and without footrot), we aimed at evaluating the importance of footrot in the ibex population, identifying potential risk factors for disease occurrence in this species, and defining the epidemiological role of ibex. Our study revealed that footrot occurs throughout the entire ibex territory (34% of the Swiss ibex colonies affected) but only as a sporadic disease (mostly one case per disease event), although the situation differed among footrot-positive colonies because half of them had experienced outbreak recurrences. Risk factor analysis for the occurrence of footrot in ibex colonies suggested an absence of an effect of meteorologic conditions, region, contacts with sheep or cattle (known to be very common healthy carriers of D. nodosus) and existing local disease control program. We found a significant effect only of contacts with sheep having footrot. Pathogen prevalence was very low in all investigated colonies. In conclusion, our results support previous data suggesting that ibex are susceptible spillover hosts, likely infected mainly by sympatric sheep displaying clinical signs.

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.

Disentangling adaptation from drift in bottlenecked and reintroduced populations of Alpine ibex

Identifying local adaptation in bottlenecked species is essential for conservation management. Selection detection methods have an important role in species management plans, assessments of adaptive capacity, and looking for responses to climate change. Yet, the allele frequency changes exploited in selection detection methods are similar to those caused by the strong neutral genetic drift expected during a bottleneck. Consequently, it is often unclear what accuracy selection detection methods have across bottlenecked populations. In this study, simulations were used to explore if signals of selection could be confidently distinguished from genetic drift across 23 bottlenecked and reintroduced populations of Alpine ibex (Capra ibex). The meticulously recorded demographic history of the Alpine ibex was used to generate comprehensive simulated SNP data. The simulated SNPs were then used to benchmark the confidence we could place in outliers identified in empirical Alpine ibex RADseq derived SNP data. Within the simulated data set, the false positive rates were high for all selection detection methods (F_ST outlier scans and Genetic‐Environment Association analyses) but fell substantially when two or more methods were combined. True positive rates were consistently low and became negligible with increased stringency. Despite finding many outlier loci in the empirical Alpine ibex SNPs, none could be distinguished from genetic drift‐driven false positives. Unfortunately, the low true positive rate also prevents the exclusion of recent local adaptation within the Alpine ibex. The baselines and stringent approach outlined here should be applied to other bottlenecked species to ensure the risk of false positive, or negative, signals of selection are accounted for in conservation management plans.

Bronchopulmonary Nematodes in Alpine Ibex: Shedding of First Stage Larvae Analyzed at the Individual Host Level

Pneumonia is the most frequent cause of death for Alpine ibex (Capra ibex) in Gran Paradiso National Park, (Italy). The etiology of this form of pneumonia is currently unknown and the identification of the primary etiological agent remains difficult due to biological and logistic constraints. Uncovering individual differences in Protostrongylid prevalence and intensity is important to further investigate the epidemiology of respiratory diseases and their relationship to heterozygosity and inbreeding in a once almost extinct population like C. ibex. In a group of 21 individually recognizable adult male we monitored monthly prevalence and intensity of Protostrongylid first-stage larvae using Baerman's technique from June to September 2019. First-stage larvae of 5 genera were detected. Muellerius (P = 100%, CI_95% = 84–100) and Protostrongylus (P = 86%, CI_95%:71–100) were two dominant genera according to Bush's importance index. Neostrongylus (P = 38%,CI_95%: 17–59), Cystocaulus (P = 33%,CI_95% = 13–53) were classified as co-dominant genera while Dictyocaulus filaria (P = 0.05%, CI_95% = 0.04–0.13) was detected, for the first time in Alpine ibex, in one subject. Protostrongylidae larval excretion varied significantly over time, with minimum L1 excretion in July. Individual median larval intensity ranged from 4.4 lpg to 82.2 lpg with Poulin's discrepancy index showing highly aggregated distribution patterns for Muellerius spp. (D = 0.283, CI_95% = 0.760–0.895) and Protostrongylus spp. (D = 0.635, CI_95% = 0.580–0.705). Presented data provide the necessary base point to further investigate how lungworm infection account for the different rates of progression of pneumonia in C. ibex. Individual aggregation of larval intensity must be further evaluated to determine whether these differences mirror different levels of parasitic infection related to individual differences in immune response, hormonal-states or genetic fitness.

Evaluating wildlife translocations using genomics: A bighorn sheep case study

Wildlife restoration often involves translocation efforts to reintroduce species and supplement small, fragmented populations. We examined the genomic consequences of bighorn sheep (Ovis canadensis) translocations and population isolation to enhance understanding of evolutionary processes that affect population genetics and inform future restoration strategies. We conducted a population genomic analysis of 511 bighorn sheep from 17 areas, including native and reintroduced populations that received 0-10 translocations. Using the Illumina High Density Ovine array, we generated datasets of 6,155 to 33,289 single nucleotide polymorphisms and completed clustering, population tree, and kinship analyses. Our analyses determined that natural gene flow did not occur between most populations, including two pairs of native herds that had past connectivity. We synthesized genomic evidence across analyses to evaluate 24 different translocation events and detected eight successful reintroductions (i.e., lack of signal for recolonization from nearby populations) and five successful augmentations (i.e., reproductive success of translocated individuals) based on genetic similarity with the source populations. A single native population founded six of the reintroduced herds, suggesting that environmental conditions did not need to match for populations to persist following reintroduction. Augmentations consisting of 18-57 animals including males and females succeeded, whereas augmentations of two males did not result in a detectable genetic signature. Our results provide insight on genomic distinctiveness of native and reintroduced herds, information on the relative success of reintroduction and augmentation efforts and their associated attributes, and guidance to enhance genetic contribution of augmentations and reintroductions to aid in bighorn sheep restoration.

Genetic monitoring of the greater stick-nest rat meta-population for strategic supplementation planning

Translocation is an increasingly common component of species conservation efforts. However, translocated populations often suffer from loss of genetic diversity and increased inbreeding, and thus may require active management to establish gene flow across isolated populations. Assisted gene flow can be laborious and costly, so recipient and source populations should be carefully chosen to maximise genetic diversity outcomes. The greater stick-nest rat (GSNR, Leporillus conditor), a threatened Australian rodent, has been the focus of a translocation program since 1985, resulting in five extant translocated populations (St Peter Island, Reevesby Island, Arid Recovery, Salutation Island and Mt Gibson), all derived from a remnant wild population on the East and West Franklin Islands. We evaluated the genetic diversity in all extant GSNR populations using a large single nucleotide polymorphism dataset with the explicit purpose of informing future translocation planning. Our results show varying levels of genetic divergence, inbreeding and loss of genetic diversity in all translocated populations relative to the remnant source on the Franklin Islands. All translocated populations would benefit from supplementation to increase genetic diversity, but two—Salutation Island and Mt Gibson—are of highest priority. We recommend a targeted admixture approach, in which animals for supplementation are sourced from populations that have low relatedness to the recipient population. Subject to assessment of contemporary genetic diversity, St Peter Island and Arid Recovery are the most appropriate source populations for genetic supplementation. Our study demonstrates an effective use of genetic surveys for data-driven management of threatened species.

Genetic epidemiology of the Alpine ibex reservoir of persistent and virulent brucellosis outbreak

While it is now broadly accepted that inter-individual variation in the outcomes of host-pathogen interactions is at least partially genetically controlled, host immunogenetic characteristics are rarely investigated in wildlife epidemiological studies. Furthermore, most immunogenetic studies in the wild focused solely on the major histocompatibility complex (MHC) diversity despite it accounts for only a fraction of the genetic variation in pathogen resistance. Here, we investigated immunogenetic diversity of the Alpine ibex (Capra ibex) population of the Bargy massif, reservoir of a virulent outbreak of brucellosis. We analysed the polymorphism and associations with disease resistance of the MHC Class II Drb gene and several non-MHC genes (Toll-like receptor genes, Slc11A1) involved in the innate immune response to Brucella in domestic ungulates. We found a very low neutral genetic diversity and a unique MHC Drb haplotype in this population founded few decades ago from a small number of individuals. By contrast, other immunity-related genes have maintained polymorphism and some showed significant associations with the brucellosis infection status hence suggesting a predominant role of pathogen-mediated selection in their recent evolutionary trajectory. Our results highlight the need to monitor immunogenetic variation in wildlife epidemiological studies and to look beyond the MHC.

Purging of highly deleterious mutations through severe bottlenecks in Alpine ibex

Human activity has caused dramatic population declines in many wild species. The resulting bottlenecks have a profound impact on the genetic makeup of a species with unknown consequences for health. A key genetic factor for species survival is the evolution of deleterious mutation load, but how bottleneck strength and mutation load interact lacks empirical evidence. We analyze 60 complete genomes of six ibex species and the domestic goat. We show that historic bottlenecks rather than the current conservation status predict levels of genome-wide variation. By analyzing the exceptionally well-characterized population bottlenecks of the once nearly extinct Alpine ibex, we find genomic evidence of concurrent purging of highly deleterious mutations but accumulation of mildly deleterious mutations. This suggests that recolonization bottlenecks induced both relaxed selection and purging, thus reshaping the landscape of deleterious mutation load. Our findings highlight that even populations of ~1000 individuals can accumulate mildly deleterious mutations. Conservation efforts should focus on preventing population declines below such levels to ensure long-term survival of species.

Although there is extensive theory predicting the effects of population bottlenecks on mutation load, there is little empirical evidence from recent bottlenecks. Here, Grossen et al. compare the consequences of population bottlenecks in six ibex species for genome-wide variation and mutation load.

Inbreeding reduces long-term growth of Alpine ibex populations

Many studies document negative inbreeding effects on individuals, and conservation efforts to preserve rare species routinely employ strategies to reduce inbreeding. Despite this, there are few clear examples in nature of inbreeding decreasing the growth rates of populations, and the extent of population-level effects of inbreeding in the wild remains controversial. Here, we take advantage of a long-term dataset of 26 reintroduced Alpine ibex (Capra ibex ibex) populations spanning nearly 100 years to show that inbreeding substantially reduced per capita population growth rates, particularly for populations in harsher environments. Populations with high average inbreeding (F ≈ 0.2) had population growth rates reduced by 71% compared with populations with no inbreeding. Our results show that inbreeding can have long-term demographic consequences even when environmental variation is large and deleterious alleles may have been purged during bottlenecks. Thus, efforts to guard against inbreeding effects in populations of endangered species have not been misplaced.

The genetic legacy of 50 years of desert bighorn sheep translocations

Conservation biologists have increasingly used translocations to mitigate population declines and restore locally extirpated populations. Genetic data can guide the selection of source populations for translocations and help evaluate restoration success. Bighorn sheep (Ovis canadensis) are a managed big game species that suffered widespread population extirpations across western North America throughout the early 1900s. Subsequent translocation programs have successfully re-established many formally extirpated bighorn herds, but most of these programs pre-date genetically informed management practices. The state of Nevada presents a particularly well-documented case of decline followed by restoration of extirpated herds. Desert bighorn sheep (O. c. nelsoni) populations declined to less than 3,000 individuals restricted to remnant herds in the Mojave Desert and a few locations in the Great Basin Desert. Beginning in 1968, the Nevada Department of Wildlife translocated ~2,000 individuals from remnant populations to restore previously extirpated areas, possibly establishing herds with mixed ancestries. Here, we examined genetic diversity and structure among remnant herds and the genetic consequences of translocation from these herds using a genotyping-by-sequencing approach to genotype 17,095 loci in 303 desert bighorn sheep. We found a signal of population genetic structure among remnant Mojave Desert populations, even across geographically proximate mountain ranges. Further, we found evidence of a genetically distinct, potential relict herd from a previously hypothesized Great Basin lineage of desert bighorn sheep. The genetic structure of source herds was clearly reflected in translocated populations. In most cases, herds retained genetic evidence of multiple translocation events and subsequent admixture when founded from multiple remnant source herds. Our results add to a growing literature on how population genomic data can be used to guide and monitor restoration programs.

A diet rich in C3 plants reveals the sensitivity of an alpine mammal to climate change

Plant-herbivore interactions provide critical insights into the mechanisms that govern the spatiotemporal distributions of organisms. These interactions are crucial to understanding the impacts of climate change, which are likely to have an effect on the population dynamics of alpine herbivores. The Royle's pika (Ochotona roylei, hereafter pika) is a lagomorph found in the western Himalaya and is dependent on alpine plants that are at risk from climate change. As the main prey of many carnivores in the region, the pika plays a crucial role in trophic interactions. We examined topographical features, plant genera presence and seasonal dynamics as drivers of the plant richness in the pika's diet across an elevational gradient (2,600-4,450 m). We identified 79 plant genera in the faecal pellets of pikas, of which 89% were forbs, >60% were endemic to the Himalaya, and 97.5% of the diet plant genera identified followed the C3 photosynthetic pathway. We found that, during the premonsoon season, the number of genera in the pika's diet decreased with increasing elevation. We demonstrate that a large area of talus supports greater plant diversity and, not surprisingly, results in higher species richness in the pika's diet. However, in talus habitat with deep crevices, pikas consumed fewer plant genera suggesting they may be foraging suboptimally due to predation risk. The continued increase in global temperature is expected to have an effect on the distribution dynamics of C3 plants and consequently influence pika diet and distribution, resulting in a significant negative cascading effect on the Himalayan ecosystem.

Human-Mediated Dispersal and the Rewiring of Spatial Networks

Humans fundamentally affect dispersal, directly by transporting individuals and indirectly by altering landscapes and natural vectors. This human-mediated dispersal (HMD) modifies long-distance dispersal, changes dispersal paths, and overall benefits certain species or genotypes while disadvantaging others. HMD is leading to radical changes in the structure and functioning of spatial networks, which are likely to intensify as human activities increase in scope and extent. Here, we provide an overview to guide research into HMD and the resulting rewiring of spatial networks, making predictions about the ecological and evolutionary consequences and how these vary according to spatial scale and the traits of species. Future research should consider HMD holistically, assessing the range of direct and indirect processes to understand the complex impacts on eco-evolutionary dynamics.

Heterozygosity-fitness correlation at the major histocompatibility complex despite low variation in Alpine ibex (Capra ibex)

Crucial for the long-term survival of wild populations is their ability to fight diseases. Disease outbreaks can lead to severe population size reductions, which makes endangered and reintroduced species especially vulnerable. In vertebrates, the major histocompatibility complex (MHC) plays an important role in determining the immune response. Species that went through severe bottlenecks often show very low levels of genetic diversity at the MHC. Due to the known link between the MHC and immune response, such species are expected to be at particular risk in case of disease outbreaks. However, so far, only few studies have shown that low MHC diversity is correlated with increased disease susceptibility in species after severe bottlenecks. We investigated genetic variation at the MHC and its correlations with disease resistance and other fitness-related traits in Alpine ibex (Capra ibex), a wild goat species that underwent a strong bottleneck in the last century and that is known to have extremely low genetic variability, both genome-wide and at the MHC. We studied MHC variation in male ibex of Gran Paradiso National Park, the population used as a source for all postbottleneck reintroductions. We found that individual MHC heterozygosity (based on six microsatellites) was not correlated with genome-wide neutral heterozygosity. MHC heterozygosity, but not genome-wide heterozygosity, was positively correlated with resistance to infectious keratoconjunctivitis and with body mass. Our results show that genetic variation at the MHC plays an important role in disease resistance and, hence, should be taken into account for successfully managing species conservation.

Introduction history overrides social factors in explaining genetic structure of females in Mediterranean mouflon

Fine‐scale spatial genetic structure of populations results from social and spatial behaviors of individuals such as sex‐biased dispersal and philopatry. However, the demographic history of a given population can override such socio‐spatial factors in shaping genetic variability when bottlenecks or founder events occurred in the population. Here, we investigated whether socio‐spatial organization determines the fine‐scale genetic structure for both sexes in a Mediterranean mouflon (Ovis gmelini musimon × Ovis sp.) population in southern France 60 years after its introduction. Based on multilocus genotypes at 16 loci of microsatellite DNA (n = 230 individuals), we identified three genetic groups for females and two for males, and concurrently defined the same number of socio‐spatial units using both GPS‐collared individuals (n = 121) and visual resightings of marked individuals (n = 378). The socio‐spatial and genetic structures did not match, indicating that the former was not the main driver of the latter for both sexes. Beyond this structural mismatch, we found significant, yet low, genetic differentiation among female socio‐spatial groups, and no genetic differentiation in males, with this suggesting female philopatry and male‐biased gene flow, respectively. Despite spatial disconnection, females from the north of the study area were genetically closer to females from the south, as indicated by the spatial analysis of the genetic variability, and this pattern was in accordance with the common genetic origin of their founders. To conclude, more than 14 generations later, genetic signatures of first introduction are not only still detectable among females, but they also represent the main factor shaping their present‐time genetic structure.

Population genomics analyses of European ibex species show lower diversity and higher inbreeding in reintroduced populations

Restoration of lost species ranges to their native distribution is key for the survival of endangered species. However, reintroductions often fail and long‐term genetic consequences are poorly understood. Alpine ibex (Capra ibex) are wild goats that recovered from <100 individuals to ~50,000 within a century by population reintroductions. We analyzed the population genomic consequences of the Alpine ibex reintroduction strategy. We genotyped 101,822 genomewide single nucleotide polymorphism loci in 173 Alpine ibex, the closely related Iberian ibex (Capra pyrenaica) and domestic goat (Capra hircus). The source population of all Alpine ibex maintained genetic diversity comparable to Iberian ibex, which experienced less severe bottlenecks. All reintroduced Alpine ibex populations had individually and combined lower levels of genetic diversity than the source population. The reintroduction strategy consisted of primary reintroductions from captive breeding and secondary reintroductions from established populations. This stepwise reintroduction strategy left a strong genomic footprint of population differentiation, which increased with subsequent rounds of reintroductions. Furthermore, analyses of genomewide runs of homozygosity showed recent inbreeding primarily in individuals of reintroduced populations. We showed that despite the rapid census recovery, Alpine ibex carry a persistent genomic signature of their reintroduction history. We discuss how genomic monitoring can serve as an early indicator of inbreeding.

Molecular Analyses Reveal Unexpected Genetic Structure in Iberian Ibex Populations

BACKGROUND

Genetic differentiation in historically connected populations could be the result of genetic drift or adaptation, two processes that imply a need for differing strategies in population management. The aim of our study was to use neutral genetic markers to characterize C. pyrenaica populations genetically and examine results in terms of (i) demographic history, (ii) subspecific classification and (iii) the implications for the management of Iberian ibex.

METHODOLOGY/PRINCIPAL FINDINGS

We used 30 neutral microsatellite markers from 333 Iberian ibex to explore genetic diversity in the three main Iberian ibex populations in Spain corresponding to the two persisting subspecies (victoria and hispanica). Our molecular analyses detected recent genetic bottlenecks in all the studied populations, a finding that coincides with the documented demographic decline in C. pyrenaica in recent decades. Genetic divergence between the two C. pyrenaica subspecies (hispanica and victoriae) was substantial (FST between 0.39 and 0.47). Unexpectedly, we found similarly high genetic differentiation between two populations (Sierra Nevada and Maestrazgo) belonging to the subspecies hispanica. The genetic pattern identified in our study could be the result of strong genetic drift due to the severe genetic bottlenecks in the studied populations, caused in turn by the progressive destruction of natural habitat, disease epidemics and/or uncontrolled hunting.

CONCLUSIONS

Previous Capra pyrenaica conservation decision-making was based on the clear distinction between the two subspecies (victoriae and hispanica); yet our paper raises questions about the usefulness for conservation plans of the distinction between these subspecies.

Microsatellite diversity of the Nordic type of goats in relation to breed conservation: how relevant is pure ancestry?

In the last decades, several endangered breeds of livestock species have been re-established effectively. However, the successful revival of the Dutch and Danish Landrace goats involved crossing with exotic breeds and the ancestry of the current populations is therefore not clear. We have generated genotypes for 27 FAO-recommended microsatellites of these landraces and three phenotypically similar Nordic-type landraces and compared these breeds with central European, Mediterranean and south-west Asian goats. We found decreasing levels of genetic diversity with increasing distance from the south-west Asian domestication site with a south-east-to-north-west cline that is clearly steeper than the Mediterranean east-to-west cline. In terms of genetic diversity, the Dutch Landrace comes next to the isolated Icelandic breed, which has an extremely low diversity. The Norwegian coastal goat and the Finnish and Icelandic landraces are clearly related. It appears that by a combination of mixed origin and a population bottleneck, the Dutch and Danish Land-races are separated from the other breeds. However, the current Dutch and Danish populations with the multicoloured and long-horned appearance effectively substitute for the original breed, illustrating that for conservation of cultural heritage, the phenotype of a breed is more relevant than pure ancestry and the genetic diversity of the original breed. More in general, we propose that for conservation, the retention of genetic diversity of an original breed and of the visual phenotype by which the breed is recognized and defined needs to be considered separately.

Utility of several microsatellite markers for the genetic characterisation of three ex situ populations of threatened caprine taxa (<i>Capra aegagrus</i>, <i>C. cylindricornis</i> and <i>C. falconeri</i>)

Abstract. Caprines belong to the most endangered group of mammals and artiodactyls suffering from many negative human impacts. Fortunately, many of them are protected and managed by national and international legislation and in situ and ex situ conservation actions. Although many microsatellite markers have been developed for wild and domestic caprines, they remain uninvestigated in respect of their utility for some taxa. We examined the utility of the International Society for Animal Genetics microsatellite set for genetic characterisations of three wild and one domestic Capra species from captive or semi-captive ex situ populations in Europe. Our data suggest the utility of this microsatellite set for detecting shared and species-specific alleles, characterising the genetic variability, and determining phylogenetic relationships and intraspecific structures in investigated taxa. We detected a depleted genetic variability in Capra falconeri and Capra cylindricornis in European ex situ populations; unrelated individuals are therefore needed for improving genetic variability parameters, as they are for the extralimital population of Capra aegagrus in the Vřísek game reserve (Czech Republic), for which we identified no genetic introgression from the domestic goat and great dissimilarity with some analysed individuals from European zoos. Current results here indicate some difficulties with the historical evidence, for example with respect to the origin and purity of particular individuals under breeding programmes.

Comparative genetic diversity in a sample of pony breeds from the U.K. and North America: a case study in the conservation of global genetic resources

Most species exist as subdivided ex situ daughter population(s) derived from a single original group of individuals. Such subdivision occurs for many reasons both natural and manmade. Traditional British and Irish pony breeds were introduced to North America (U.S.A. and Canada) within the last 150 years, and subsequently equivalent breed societies were established. We have analyzed selected U.K. and North American equivalent pony populations as a case study for understanding the relationship between putative source and derived subpopulations. Diversity was measured using mitochondrial DNA and a panel of microsatellite markers. Genetic signatures differed between the North American subpopulations according to historical management processes. Founder effect and stochastic drift was apparent, particularly pronounced in some breeds, with evidence of admixture of imported mares of different North American breeds. This demonstrates the importance of analysis of subpopulations to facilitate understanding the genetic effects of past management practices and to lead to informed future conservation strategies.

Direct and indirect causal effects of heterozygosity on fitness-related traits in Alpine ibex

Heterozygosity-fitness correlations (HFCs) are a useful tool to investigate the effects of inbreeding in wild populations, but are not informative in distinguishing between direct and indirect effects of heterozygosity on fitness-related traits. We tested HFCs in male Alpine ibex (Capra ibex) in a free-ranging population (which suffered a severe bottleneck at the end of the eighteenth century) and used confirmatory path analysis to disentangle the causal relationships between heterozygosity and fitness-related traits. We tested HFCs in 149 male individuals born between 1985 and 2009. We found that standardized multi-locus heterozygosity (MLH), calculated from 37 microsatellite loci, was related to body mass and horn growth, which are known to be important fitness-related traits, and to faecal egg counts (FECs) of nematode eggs, a proxy of parasite resistance. Then, using confirmatory path analysis, we were able to show that the effect of MLH on horn growth was not direct but mediated by body mass and FEC. HFCs do not necessarily imply direct genetic effects on fitness-related traits, which instead can be mediated by other traits in complex and unexpected ways.

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.

Retrospective and prospective perspectives on zoonotic brucellosis

Members of the genus Brucella are pathogenic bacteria exceedingly well adapted to their hosts. The bacterium is transmitted by direct contact within the same host species or accidentally to secondary hosts, such as humans. Human brucellosis is strongly linked to the management of domesticated animals and ingestion of their products. Since the domestication of ungulates and dogs in the Fertile Crescent and Asia in 12000 and 33000 ya, respectively, a steady supply of well adapted emergent Brucella pathogens causing zoonotic disease has been provided. Likewise, anthropogenic modification of wild life may have also impacted host susceptibility and Brucella selection. Domestication and human influence on wild life animals are not neutral phenomena. Consequently, Brucella organisms have followed their hosts’ fate and have been selected under conditions that favor high transmission rate. The “arm race” between Brucella and their preferred hosts has been driven by genetic adaptation of the bacterium confronted with the evolving immune defenses of the host. Management conditions, such as clustering, selection, culling, and vaccination of Brucella preferred hosts have profound influences in the outcome of brucellosis and in the selection of Brucella organisms. Countries that have controlled brucellosis systematically used reliable smooth live vaccines, consistent immunization protocols, adequate diagnostic tests, broad vaccination coverage and sustained removal of the infected animals. To ignore and misuse tools and strategies already available for the control of brucellosis may promote the emergence of new Brucella variants. The unrestricted use of low-efficacy vaccines may promote a “false sense of security” and works towards selection of Brucella with higher virulence and transmission potential.

The impact of translocations on neutral and functional genetic diversity within and among populations of the Seychelles warbler

Translocations are an increasingly common tool in conservation. The maintenance of genetic diversity through translocation is critical for both the short- and long-term persistence of populations and species. However, the relative spatio-temporal impacts of translocations on neutral and functional genetic diversity, and how this affects genetic structure among the conserved populations overall, have received little investigation. We compared the impact of translocating different numbers of founders on both microsatellite and major histocompatibility complex (MHC) class I diversity over a 23-year period in the Seychelles warbler (Acrocephalus sechellensis). We found low and stable microsatellite and MHC diversity in the source population and evidence for only a limited loss of either type of diversity in the four new populations. However, we found evidence of significant, but low to moderate, genetic differentiation between populations, with those populations established with fewer founders clustering separately. Stochastic genetic capture (as opposed to subsequent drift) was the main determinant of translocated population diversity. Furthermore, a strong correlation between microsatellite and MHC differentiation suggested that neutral processes outweighed selection in shaping MHC diversity in the new populations. These data provide important insights into how to optimize the use of translocation as a conservation tool.

The number of markers and samples needed for detecting bottlenecks under realistic scenarios, with and without recovery: a simulation-based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection method sexist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. (2012) analysed the power of two approaches, the M-ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100),which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when post bottle neck recovery has occurred. For the realistic situations examined,we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals,(ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii)recovery after a moderate amount of time or gradual change in size reduces power, by up to one-half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated,using simulation modelling, as a routine component of molecular ecology studies.

European springtime temperature synchronises ibex horn growth across the eastern Swiss Alps

Direct effects of climate change on animal physiology, and indirect impacts from disruption of seasonal synchrony and breakdown of trophic interactions are particularly severe in Arctic and Alpine ecosystems. Unravelling biotic from abiotic drivers, however, remains challenging because high-resolution animal population data are often limited in space and time. Here, we show that variation in annual horn growth (an indirect proxy for individual performance) of 8043 male Alpine ibex (Capra ibex) over the past four decades is well synchronised among eight disjunct colonies in the eastern Swiss Alps. Elevated March to May temperatures, causing premature melting of Alpine snowcover, earlier plant phenology and subsequent improvement of ibex food resources, fuelled annual horn growth. These results reveal dependency of local trophic interactions on large-scale climate dynamics, and provide evidence that declining herbivore performance is not a universal response to global warming even for high-altitude populations that are also harvested.

Genetic structure of introduced populations: 120-year-old DNA footprint of historic introduction in an insular small mammal population

Wildlife populations have been introduced to new areas by people for centuries, but this human-mediated movement can disrupt natural patterns of genetic structure by altering patterns of gene flow. Insular populations are particularly prone to these influences due to limited opportunities for natural dispersal onto islands. Consequently, understanding how genetic patterns develop in island populations is important, particularly given that islands are frequently havens for protected wildlife. We examined the evolutionary origins and extent of genetic structure within the introduced island population of red squirrels (Sciurus vulgaris) on the Channel Island of Jersey using mitochondrial DNA (mtDNA) control region sequence and nuclear microsatellite genotypes. Our findings reveal two different genetic origins and a genetic architecture reflective of the introductions 120 years ago. Genetic structure is marked within the maternally inherited mtDNA, indicating slow dispersal of female squirrels. However, nuclear markers detected only weak genetic structure, indicating substantially greater male dispersal. Data from both mitochondrial and nuclear markers support historic records that squirrels from England were introduced to the west of the island and those from mainland Europe to the east. Although some level of dispersal and introgression across the island between the two introductions is evident, there has not yet been sufficient gene flow to erase this historic genetic "footprint." We also investigated if inbreeding has contributed to high observed levels of disease, but found no association. Genetic footprints of introductions can persist for considerable periods of time and beyond traditional timeframes of wildlife management.

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

Human influence typically impacts on natural populations of conservation interest. These interactions are varied and sometimes complex, and may be negative and unintended or associated with conservation and management strategy. Understanding the details of how these interactions influence and are influenced by natural evolutionary processes is essential to the development of effective conservation strategies. In this study, we investigate a species in Britain that has experienced both negative impact through overhunting in historical times and management efforts through culls and translocations. At the same time, there are regional populations that have been less affected by human influence. We use mtDNA and nuclear microsatellite DNA markers to investigate patterns of connectivity and diversity and find multiple insular populations in Britain that probably evolved within the Holocene (when the habitat was free of ice). We identify three concurrent processes. First, surviving indigenous populations show highly provincial patterns of philopatry, maintaining and generating population structure on a small geographic scale. Second, founder populations into habitat extirpated of native populations have expanded, but remained largely insular. Third, introductions into established populations generate some admixture. We discuss the implications for the evolution of diversity of the integration of natural processes with anthropogenic influences on population size and distribution.

A novel approach for choosing summary statistics in approximate Bayesian computation

The choice of summary statistics is a crucial step in approximate Bayesian computation (ABC). Since statistics are often not sufficient, this choice involves a trade-off between loss of information and reduction of dimensionality. The latter may increase the efficiency of ABC. Here, we propose an approach for choosing summary statistics based on boosting, a technique from the machine-learning literature. We consider different types of boosting and compare them to partial least-squares regression as an alternative. To mitigate the lack of sufficiency, we also propose an approach for choosing summary statistics locally, in the putative neighborhood of the true parameter value. We study a demographic model motivated by the reintroduction of Alpine ibex (Capra ibex) into the Swiss Alps. The parameters of interest are the mean and standard deviation across microsatellites of the scaled ancestral mutation rate (θ(anc) = 4N(e)u) and the proportion of males obtaining access to matings per breeding season (ω). By simulation, we assess the properties of the posterior distribution obtained with the various methods. According to our criteria, ABC with summary statistics chosen locally via boosting with the L(2)-loss performs best. Applying that method to the ibex data, we estimate θ(anc)≈ 1.288 and find that most of the variation across loci of the ancestral mutation rate u is between 7.7 × 10(-4) and 3.5 × 10(-3) per locus per generation. The proportion of males with access to matings is estimated as ω≈ 0.21, which is in good agreement with recent independent estimates.

Loss of genetic diversity means loss of geological information: the endangered Japanese crayfish exhibits remarkable historical footprints

Intra-specific genetic diversity is important not only because it influences population persistence and evolutionary potential, but also because it contains past geological, climatic and environmental information. In this paper, we show unusually clear genetic structure of the endangered Japanese crayfish that, as a sedentary species, provides many insights into lesser-known past environments in northern Japan. Over the native range, most populations consisted of unique 16S mtDNA haplotypes, resulting in significant genetic divergence (overall F_ST = 0.96). Owing to the simple and clear structure, a new graphic approach unraveled a detailed evolutionary history; regional crayfish populations were comprised of two distinct lineages that had experienced contrasting demographic processes (i.e. rapid expansion vs. slow stepwise range expansion) following differential drainage topologies and past climate events. Nuclear DNA sequences also showed deep separation between the lineages. Current ocean barriers to dispersal did not significantly affect the genetic structure of the freshwater crayfish, indicating the formation of relatively recent land bridges. This study provides one of the best examples of how phylogeographic analysis can unravel a detailed evolutionary history of a species and how this history contributes to the understanding of the past environment in the region. Ongoing local extinctions of the crayfish lead not only to loss of biodiversity but also to the loss of a significant information regarding past geological and climatic events.