Skip to main content

Evidence of genetic structure in the wide-ranging bearded vulture (Gypaetus barbatus (Linnaeus, 1758))



The bearded vulture is sparsely distributed across a wide geographic range that extends over three continents (Africa, Europe and Asia). Restriction to high-altitude mountainous habitats, low breeding rates, lack of food and a heightened level of persecution have left many local populations severely diminished or extinct. Understanding the genetic connectivity and population structure of this threatened vulture species is critical for accurately assessing their conservation status, and for appropriately managing local populations through captive breeding programmes or translocations. Previous genetic assessments of the species were mainly focused on the European and Asian populations and included limited representation of the geographically isolated southern African population. A single mitochondrial study, which focused on the African populations of the bearded vulture, detected limited genetic differentiation between populations in Ethiopia and southern Africa, with reduced haplotype diversity in the southern Africa population. In this study, we extend the previous genetic assessments of the species by examining the phylogeography and genetic connectivity of global G. barbatus populations using a panel of 14 microsatellite loci.


Analyses revealed spatially correlated genetic differentiation between regional populations and low levels of gene flow between these population fragments. In contrast to the mitochondrial data, the microsatellite data support the management of genetically different populations as separate entities.


Low genetic diversity and geographic isolation are known to adversely affect the evolutionary potential of a species in the long-term. The high inbreeding found in the southern African G. barbatus and, to a lesser extent, the northern African populations highlights the need for conservation programmes to effectively manage populations of this species and maintain extant genetic diversity.


Historically, bearded vultures, Gypaetus barbatus, were considerably more abundant, albeit always sparsely distributed across a wide area ranging from the Palearctic, through Afrotropic and into Indomalay regions [1]. In the past century however, populations of this species have become locally extinct, or survive in highly habitat-specific isolated refuges across their former range [2,3,4]. A population in the Pyrenees, the islands of Corsica and Crete, and a reintroduced population in the Alps currently comprise the European bearded vulture population; while southern Africa, Ethiopia and Morocco are the remaining habitats for the African bearded vultures [1, 2].

In Europe, the population declined in the nineteenth century due to food shortage and direct persecution. This led to local extinction in the Alps by the early twentieth century [5]. The Pyrenees is home to the largest population of bearded vultures in Europe, but this population has also declined severely [4, 6]. The populations on Sardinia (Italy) and in the Balkans became extinct in the 1960s and 2000s due to unintentional poisoning [7] and decreased food availability [8]. Vulture populations still occur on Corsica (France) and Crete (Greece) but these populations are small (less than 10 breeding pairs) and continued existence is precarious [8]. In 1986, the first birds from captive breeding were released back into the wild in the Alps. Currently (2019), the estimated population size in the Alps is 280–340 birds and 52 breeding pairs. The main threats for European bearded vultures remain non-intentional poisoning and collision with energy infrastructure [4, 9].

In Africa, habitat loss, electrocution and collision with energy infrastructure, unintentional poisoning through veterinary drug-use, feeding on poison baits aimed at predator control and intentional poisoning for the use of vulture body parts in the illegal wildlife and the traditional medicine trade, are the primary culprits for the population decline in many vulture species [10,11,12,13,14,15,16]. These threats have contributed to the fragmentation of the African bearded vulture population into several isolated pockets [10, 12, 15] (Fig. 1).

Fig. 1

The bearded vulture’s Gypaetus barbatus global distribution (both extinct and extant) across the Palearctic, Afrotropic and Indomalay regions. Numbered circles correspond to broad sampling locality of individuals included in this study (Additional file 1: Table S4). Map created using ArcGIS and Adobe Illustrator

Morocco is home to fewer than 20 individuals and the species is considered extinct in Tunisia [17]. Ethiopia is believed to hold a few hundred pairs, however, reliable data are not available and these values are therefore thought to be considerably overestimated [18, 19]. A 2011-survey identified only three nest-sites in Kenya, and roughly double that in Tanzania [20]. The entire southern hemisphere population is restricted to the Highlands of Lesotho and south-eastern South Africa, where they roost and breed on high altitude cliffs of the Maloti-Drakensberg mountain border between South Africa (the Free State, KwaZulu-Natal and Eastern Cape Provinces) and Lesotho [2].

Due to the reduced effective population size, small disjunct populations experience amplified genetic drift, inbreeding and a greater loss of genetic diversity [21,22,23]. Loss of genetic variability directly impacts long-term persistence [24,25,26,27].

Two subspecies of bearded vulture have been recognized based on plumage characteristics and geography [28]—Gypaetus barbatus barbatus is distributed across Europe and Asia, extending into north-west Africa (mainly in High Atlas, Morocco), and Gypaetus barbatus meridionalis is found in north-east Africa, east Africa and extends into southern Africa. Studies based on maternally inherited mitochondrial data found evidence for two distinct lineages of bearded vultures which did not correspond geographically to the morphologically described subspecies. The first lineage described by Godoy et al. (2004) occurs predominantly in western Europe, and the other lineage extends over Africa, eastern Europe and central Asia. This study detected high levels of inter-population differentiation, with admixture between the two lineages occurring in Greece and the Alps [1]. Godoy et al. (2004) hypothesised the species to be genetically structured as a consequence of the extinction of the central European population. The two lineages underwent allopatric differentiation (Mediterranean vs. Africa/Asian), and range expansion in the African/Asian cohort which resulted in secondary contact zones in central Europe and northern Africa [1]. These contemporary, remnant populations were found to be ecologically interchangeable, even between vastly segregated populations. A further study including the southern African population of bearded vultures based on mitochondrial sequence data, found that, despite vast geographical distances separating populations, there is little genetic differentiation between the isolated southern African population and the northern hemisphere bearded vulture populations [29]. Most genetic studies have recorded low levels of genetic diversity in bearded vulture populations [1, 29, 30].

The genetic results have implications for captive breeding programmes and open the possibility for the translocation of individuals into shrinking populations. Unfortunately, previous studies have focussed their analyses almost entirely on the European population (with limited representation from southern Africa). The single study that included good representation of the isolated southern African population was based on mitochondrial DNA [29]. In this study we re-examine the population structure of the bearded vulture using improved sampling across the global distribution of the species and examining genetic connectivity among geographically isolated populations using a suite of biparentally inherited nuclear microsatellite markers. We test the hypothesis that the bearded vulture is a single, widely distributed species that shows shallow genetic structure despite the isolated nature of many populations. Estimating connectivity of disjunct populations will provide conservation authorities with the necessary information to form and implement the appropriate conservation measures to better protect local populations of the bearded vulture.


A total of 236 individuals were genotyped with the 14 microsatellite markers. All loci were polymorphic with the number of alleles per individual ranging from 6 (Gf11A4) to 17 (BV6, BV12). A summary of locus-by-locus measures of genetic diversity are provided in Additional file 1: Table S1. The mean null allele frequency was 18.9%. Paired t-tests, however, did not indicate a significant difference between uncorrected and corrected FST values (p-value > 0.05), suggesting that null alleles have a very limited effect on the genetic structuring analyses in these samples of bearded vultures. In addition, many of the analyses we conducted are not strongly influenced by the presence of null alleles [31]. Linkage disequilibrium analysis detected significant (p < 0.001) associations between some pairs of loci (Additional file 1: Table S4). Processes such as genetic drift, gene flow among populations with dissimilar allele frequencies and small sample size, can lead to erroneous detection of genetic association where there is none. Given that these microsatellites have been used in previous studies, all 14 microsatellite loci were used in subsequent analyses.

Significant departures (p < 0.05) from HWE were obtained for all loci in the southern African population. Loci 1, 3 and 4 were in HWE for the European, northern African and Eastern population complexes, respectively. The remaining loci did not conform to HWE predictions. This was not unexpected as these small isolated populations naturally deviate from the assumptions of HWE.

Genetic diversity

The average number of alleles and effective alleles ranged from 5.93–7.00 and 2.81–3.23 respectively (Table 1). All populations tended towards being inbred (FIS > 1, Table 1), and showed higher levels of homozygosity than would be expected under HWE (HO < HE; Table 1).

Table 1 Genetic diversity estimates for 236 bearded vultures Gypaetus barbatus

Population structure

The analysis using the Admixture parameter detected three genetic clusters (∆K = 3, Fig. 2) as the optimal partitioning scenario using both the Puechmaille and Evanno methods. The Bayesian STRUCTURE cluster plot revealed low levels of genetic admixture among populations from different geographic regions, and genetic clustering was closely linked to geography. The pattern of clustering does not support the delimitation of the two subspecies. The birds from Ethiopia and Yemen had alleles at similar frequencies to those seen in birds from north-west Africa, whereas the southern African population was distinct from other populations. For this reason, all northern African populations were grouped together in subsequent analyses.

Fig. 2

Probabilistic population structure under the Admixed model for 236 bearded vultures Gypaetus barbatus collected from across the global distribution of the species. Individuals are sorted into four broad geographically delineated regions, and are represented by vertical lines partitioned into coloured segments which represent the estimated membership coefficients (Q-value) of each cluster represented in their genetics. ∆K = 3 was estimated as the optimal partitioning strategy

The lowest admixture was evident in the European populations (Q = 0.97), followed by the southern African population (Q = 0.92) and the Asian populations (Q = 0.88). The highest admixture was detected in the northern African populations (Q = 0.63, Table 2). Results using the No Admixture parameter, returned a similar clustering pattern although in this case ∆K = 4. Barplots of K = 2 to 4 are shown in Fig. 2.

Table 2 Comparison of membership coefficients (Q-values) for both Admixture and No Admixture models used in STRUCTURE

In all analyses, the isolated southern African bearded vulture population is genetically distinct from all other bearded vulture populations. Genotypes collected from this region appear to belong to unique genetic clusters and are characterized by unique allele frequencies. The genetic distinctiveness of the southern African population is also reflected in the number of unique alleles found in this population (36 unique alleles). The European population complex contains only 18 unique alleles and the Asian and northern African complexes have fewer than 10 unique alleles each. Only 31 alleles are shared by all populations. Results from the Mantel test indicate that there was no positive correlation between genetic and geographical distance (r2 = 0.049; P < 0.001).

Population assignment tests conducted in GenAlEx support the STRUCTURE results by showing the separation of the southern African population from the other African, European and Asian populations (Fig. 3). The latter three population complexes showed significant overlap and could not be separated in this analysis.

Fig. 3

Population assignment of 236 bearded vultures Gypaetus barbatus across 14 microsatellite loci

Pairwise FST comparisons only provide moderate support for differentiation of the southern African population from the northern African populations (FST = 0.12, Table 3) but indicate that the southern African population is genetically distinct from the European and Asian populations (FST = 0.18 and 0.16, respectively, Table 3).

Table 3 Pairwise FST among the four regional populations of bearded vulture Gypaetus barbatus

AMOVA results, conducted when grouping individuals into the four regions, showed that 49% (p < 0.001) of variation was attributed to within individual differences, while differences among individuals and among populations explained 39% (p < 0.001) and 12% (p < 0.001) of the diversity respectively.

Population connectivity

The BayesAss 3.0 analysis indicates that the highest migration seems to have occurred from the European populations to the Asian populations (28.26%; 83 shared alleles). This supports the hypothesis that, before the extinction of bearded vultures in the Alps, there was an interconnected European population from Spain to the Alps, Balkans to Turkey and Asia.

Minimal migration has occurred within Africa in the recent past. Less migration took place from northern Africa to southern Africa, than the reverse (7.05% and 15.20%, Table 4). The northern African populations contain a greater percentage (15.64%, Table 4) of the European alleles than the southern African population (< 1%, Table 4), suggesting stronger linkage among northern hemisphere populations. However, alleles dominant in the Asian populations do occur, albeit at lower frequencies in both the southern and northern African populations (< 2%, Table 4). While the genetic admixture seen in the northern African bearded vulture populations may be a reflection of active migration, it is unlikely that the southern African population has experienced any recent genetic admixing with the Asian populations. Thus, this pattern may merely reflect incomplete lineage sorting.

Table 4 Pair-wise migration rates between the four regional population complexes; bold values on the diagonal show the proportion of population which are permanent residents

The clustering patterns which emerge from the EDENetworks also indicates that the four regional population complexes are separate and that there has been limited recent connectivity. EDENetwork analysis reveals stronger linkage and overlap between Asian and European populations. The southern African population shares a weak linkage with the northern African populations but is not connected genetically with any of the other regional populations (Fig. 4). The northern African populations were connected (albeit weakly) with the Asian and European populations.

Fig. 4

Network constructed in EDENetworks grouping 236 bearded vulture Gypaetus barbatus samples by region. Populations (nodes) are linked by edges which are weighted in proportion to the strength of the genetic linkage


The population structure of endangered species reflects population collapses, bottlenecking events, fragmentation and isolation of populations [24, 32, 33]. A key indicator of population health is genetic variability, as healthy levels of variation allow for populations to respond to natural selection [27, 34, 35]. The associated loss of genetic variability in small isolated populations, through genetic drift and inbreeding, is believed to adversely affect the adaptive potential of remnant populations and thus threatens their evolutionary persistence [24, 36, 37]. Past genetic studies conducted on the wide-ranging but sparsely distributed bearded vulture have either focussed primarily on the European and Asian populations, and/or have used mitochondrial markers [1, 29]. These papers report a lack of mitochondrial genetic diversity within the species. Based on the lack of adequate mitochondrial differentiation between the geographically disjunct populations in these studies, and apparent ecological and genetic interchangeability, it was recommended that the southern and northern hemisphere populations of bearded vultures be managed as a single unit [1, 29]. Genetic structuring is often prevalent in populations which are widely distributed and display strong natal philopatric behaviour [30, 38, 39]. Given that the bearded vulture populations are separated by vast geographic distances, and philopatric behaviour has been documented in the species, it would be reasonable to expect genetic structuring [30, 40]. In the research presented here, population structuring was investigated by using biparentally inherited microsatellite fragment data.

The pattern of population differentiation found here provides clarification of the population structure and migration of bearded vultures. Although there is debate over the number of subspecies of bearded vultures [41,42,43], our results advocate for the management of the South African bearded vulture population as a separate entity, G. b. meridionalis (Keyserling & Blasius, 1840). There is debate over the splitting of the northern lineage into either one, Gypaetus barbatus barbatus [43, 44], or 3 subspecies, Gypaetus barbatus aureus (Hablizl, 1783), Gypaetus barbatus hemachalanus (Hutton, 1838) and Gypaetus barbatus barbatus (Linnaeus, 1758) [41, 42]. Based on our findings, we advocate recognition of G. b. barbatus in Eurasia and northern Africa and G. b. meridionalis in eastern and southern Africa. The diversification of the northern and southern bearded vulture lineages is hypothesised to be linked to the expanding Sahara Desert which acted as a physical barrier during the last glacial maximum [1, 45]. When the Sahara receded, bearded vultures were able to disperse into East Africa and this formed a ‘steppingstone’ allowing for admixture between the two subspecies. The highest level of admixture in this study was found in the northern African populations, and this could reflect past dispersal events, as well as the subsequent isolation in the recent past. Although there appears to be divergence between the SW Arabian and southern African population, further analysis based on a larger sample size from Yemen would be required to confirm this.


Based on this species’ philopatric dispersal behaviour, genetic differentiation between regional populations is expected to increase with geographic distance following an isolation-by-distance model [46, 47], especially in the absence of translocations or supplementations to the population. Here analyses found minimal migration between the populations. The highest emigration rates were out of Europe and little to no migration was detected to, or from, the southern African population. This supports the hypothesis that the southern African population is indeed isolated, both geographically and genetically.

Southern African bearded vulture diversity

Apart from past persecution, the southern African population of bearded vultures continues to suffer threats from several sources, adding to their continued decline [20, 48]. Past studies on the global and local bearded vulture populations have raised concern for their management based on the reduced genetic variability [1, 30, 48, 49]. Given that the southern African bearded vulture population is geographically and genetically isolated, it was expected that this population should be at risk genetically. Indeed, we found evidence for this and found the southern African bearded vulture population to harbour less genetic diversity than would be expected under Hardy–Weinberg, although this genetic diversity is still higher than that seen in other African vulture species [50]. The drastic population decline and reduced genetic variability of the bearded vulture places the southern African population alongside other ecologically and behaviourally similar vulture species impacted by “The African Vulture Crisis”. This population collapse, coupled with the isolated nature of the southern African population, raises concern for the long-term population health and persistence of the species in southern African. Genetic diversity levels will continue to decline if the population bottleneck is indefinitely sustained. The genetic variation seen in contemporary populations reflects the diversity in the previous generations. Should conservation seek to preserve this population’s gene-assemblage, active management would need to firstly minimise external threats, and secondly mitigate factors leading to population structuring, as this amplifies genetic drift.

This study highlights the utility of integrating genetic data to inform conservation strategies as it elucidates information that is not readily evident in the field. Our results suggest that the southern African bearded vulture population, being distinct from other such populations, may represent a reservoir of genetic variation that should be given separate conservation status.



A total of 236 bearded vulture specimens were analysed in this study. These include samples from across the global distribution of the species (Fig. 1), and include historical as well as recent samples (1805–2012). Bearded vultures have a long generation time, and this temporal sampling represents approximately 10 generations. Unfortunately, the few historical samples preclude analysis of temporal changes in population structure. These samples were divided into four regional groupings for analyses: a southern African group, including 52 samples from a single population spanning South Africa and Lesotho; a northern African group (n = 39) including samples from Ethiopia, Yemen, Morocco and Algeria; a European group (n = 107) including samples from Albania, Greece, Turkey, Crete, Corsica, Sardinia, Switzerland, Austria, the Pyrenees, France, Spain; and an Asian group (n = 38) including samples from Caucasus, central Asia, Turkestan, Kyrgyzstan, Russia, India, and China. Details of specimens included in the study are provided in Additional file 1: Tables S1 and S5.

DNA extraction

The NucleoSpin® Tissue kit (Macherey–Nagel) was used for all DNA extractions. The standard protocol for extracting genetic material from blood was followed for blood samples. The feather and archival tissue extractions were done using the standard protocol for muscle tissue, with the following modifications to increase DNA yield: incubation of the ~ 5 mm feather tip/skin snip and proteinase K for 48 h, lysate was then incubated in B3 (buffer) for 45 min (70 °C), the final volume of pre-warmed Buffer BE was decreased to 80 μl, incubation at 70 °C for 20 min followed by centrifuging and then reapplication of the solution onto the membrane, followed by incubation at 70 °C for an additional 5 min and a final centrifugation step. DNA concentrations were determined using the ThermoScientific NanoDrop 2000 spectrophotometer (Inqaba Biotec, South Africa). All DNA extracts were stored at -20 °C.

Microsatellite amplification

Fourteen microsatellite loci (Table 5) were amplified in the current study. These microsatellite primers were developed specifically for European populations of bearded vulture [3] and for other Gyps species [51]. Microsatellite loci were amplified in four multiplex reactions, and one singleplex (annealing temperature for BV 17 excluded it from multiplexing) (Table 5) using KAPA2G Fast multiplex PCR Kit (KAPA Biosystems, Wilmington, MA, USA). PCR reactions consisted of 5 μl KAPA2G Fast Multiplex mix, 0.1- 0.2 μM of each primer, 0.5–3.5 μl of template DNA, and purified water was added to each reaction to make up a final reaction volume of 10 μl. The thermocycler cycling parameters were: 95 °C for 3 min as the initial denaturation step, 30 cycles at 95 °C for 15 s, 58–60 °C (Ta) for 30 s, 72 °C for 1 min, with a final elongation step at 72 °C for 10 min. The samples were held at 4 °C once the cycle had completed running. The amplified products were sent to the Central Analytical Facilities (Stellenbosch University, South Africa) for fragment analysis. To ensure correct genotype scoring, all genotypes from archival material were genotyped multiple times and 20% of newly collected samples were reamplified and compared to ensure consistency. Genotyping error was < 0.2%. Genotypes were scored using the software package GeneMarker® v2.4.0 (Soft Genetics).

Table 5 Microsatellite marker multiplex combinations, and associated information for primers used across populations of bearded vulture Gypaetus barbatus

Microsatellite analysis

Null alleles in data can be problematic for calculations of genetic variability as they can decrease within-population variance and so inflate FST values [52]. The presence of null alleles was tested for using FreeNA [52]. Uncorrected FST values were compared to FST values corrected using the excluding null alleles (ENA) method using a paired t-test. The presence of scoring errors, such as incorrect assigning of genotypes and fractions in genotypes, was assessed using GenAlEx v6.5 [53] and FSTAT [54]. These packages were also used to calculate the average number of alleles at each locus (NA), allelic richness (AR), inbreeding co-efficient (FIS) and gene diversity. Effective alleles at each locus (NE), observed (HO) and expected (HE) heterozygosity were calculated in GenAlEx. FSTAT was used to determine allelic richness and gene diversity. Tests for deviation from Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium were performed in Genepop v4.2 [55] and FSTAT. The fixation indices, inbreeding co-efficient (FIS) and population differentiation (FST) [56, 57] were calculated in FSTAT [58]. FIS and FST explain the deviations seen from HWE. Positive FIS values suggest higher inbreeding values (due to higher observed homozygosity) than would be expected under HWE, and vice versa. FST shows the genetic differentiation between two populations which no longer conform to the HWE assumptions.

Population structure

Population structure across the global distribution of the bearded vulture was assessed using a Bayesian clustering method implemented in the software program STRUCTURE version 2.3.4 [59, 60]. The analysis was performed applying the correlated allele frequencies model [60] using sampling localities as prior information (LOCPRIOR parameter).

Analyses were run using both the Admixture and No Admixture model parameters. If there is reason to consider each population as completely discrete, the No Admixture model is appropriate. In contrast, the Admixture model allows for a large proportion of sampled individuals to have ancestry from multiple lineages, which implies a shallower population structure. The historically wide geographic distribution of the bearded vulture favours the Admixture model, and previous molecular work suggests that the African, European and Asian populations are genetically closely linked [1, 29]. The geographically isolated nature of these populations, however, suggests that each regional population may be acting as discrete units in recent generations. For this reason, both models were implemented and the output from both models presented.

All STRUCTURE analyses were performed with 100 000 Markov-Chain Monte Carlo (MCMC) replicates with a burnin of 10 000. The number of iterations for all assignments was 10 with K ranging from 1 to 10. Once all the simulations were completed, STRUCTURE Selector [61] was used to determine the optimum ∆K value for each assignment using the Puechmaille method [62]. Detecting the optimal number of genetic clusters was also conducted in STRUCTURE Harvester [63] using the Evanno method [64]. The online software package Pophelper [65] was used to compile bar plots of the optimum K- value (∆K).

Population differentiation (FST) values were calculated in FreeNA using the ENA-corrected method. Population assignment analysis was conducted in GenAlEx. Pairwise FST values were calculated using two grouping schemes. Individuals were grouped by country of origin, and then because some countries have very low sample sizes, individuals were also grouped by region (southern Africa, northern Africa, Europe and Asia; Additional file 1: Table S1). Analysis of Molecular Variance (AMOVA) was performed after grouping individuals into the four geographic regions. Isolation-by-distance was tested for by assessing the correlation between genetic distance and geographic distance using a paired Mantel test [66] in GenAlEx.

Regional connectivity

Regional connectivity was visualized through construction of a genetic network in EDENetworks v2.18 [67]. A genetic network with individuals grouped by region was constructed with edges or links between nodes weighted by pairwise FST genetic distances. The automatic thresholding option was used to determine the optimal percolation threshold. The final clustering pattern was recalculated 10 times to test for alternate placement of nodes or connections.

Recent migration among geographically segregated populations (southern Africa, northern Africa, Europe and Asia) was estimated using the Bayesian algorithm implemented in BayesAss 3.0 [68]. The analysis was conducted using the MCMC method with a burnin of 20 000 steps and a sampling frequency of 60 000 000 iterations as suggested by the authors of the programme [68]. The delta values for each parameter were adjusted to achieve a 20–60% acceptance rate [68]. A final delta value of 0.40 for allele frequency, migration and inbreeding rates maintained acceptance levels within limits. Migration rates below 0.10 were used to indicate demographic independent populations [69].

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the Zenodo repository,





Deoxyribonucleic acid


Population differentiation


Hardy–Weinberg equilibrium


Inbreeding co-efficient

HO :

Observed heterozygosity

HE :

Expected heterozygosity


Membership co-efficient


Analysis of molecular variance


Excluding null alleles


Polymerase chain reaction


Markov-Chain Monte Carlo


  1. 1.

    Godoy JA, Negro JJ, Hiraldo F, Donazar JA. Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L.) as revealed by mitochondrial DNA. Mol Ecol. 2004;13(2):371–90.

  2. 2.

    Brown C. Population dynamics of the bearded vulture Gypaetus barbatus in southern Africa. Afr J Ecol. 1997;35(1):53–63.

    CAS  Article  Google Scholar 

  3. 3.

    Gautschi B, Tenzer I, Müller JP, Schmid B. Isolation and characterization of microsatellite loci in the bearded vulture (Gypaetus barbatus) and cross-amplification in three Old World vulture species. Mol Ecol. 2000;9(12):2193–5.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Margalida A, Heredia R, Razin M, Hernández M. Sources of variation in mortality of the Bearded Vulture Gypaetus barbatus in Europe. Bird Conserv Int. 2008;18(1):1–10.

    Article  Google Scholar 

  5. 5.

    Mingozzi T, Estève R. Analysis of a historical extirpation of the Bearded vulture Gypaetus barbatus (L.) in the Western Alps (France-Italy): former distribution and causes of extirpation. Biol Conserv. 1997;79(2–3):155–71.

  6. 6.

    Hernández M, Margalida A. Assessing the risk of lead exposure for the conservation of the endangered Pyrenean bearded vulture (Gypaetus barbatus) population. Environ Res. 2009;109(7):837–42.

    PubMed  Article  CAS  Google Scholar 

  7. 7.

    Parvanov D, Stoynov E, Vangelova N, Peshev H, Grozdanov A, Delov V, et al. Vulture mortality resulting from illegal poisoning in the southern Balkan Peninsula. Environ Sci Pollut Res Int. 2018;25(2):1706–12.

    PubMed  Article  Google Scholar 

  8. 8.

    Tewes E, Terrasse M, Sánchez J, Fremuth W, Frey H. Action Plan for the Recovery and Conservation of Vultures on the Balkan Peninsula: activities and projects during 2002 and 2003. Raptors Worldwide WWGBP/MME. 2004.

  9. 9.

    Plaza PI, Martinez-Lopez E, Lambertucci SA. The perfect threat: pesticides and vultures. Sci Total Environ. 2019;687:1207–18.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Brown C. An investigation into the decline of the Bearded Vulture Gypaetus barbatus in southern Africa. Biol Conserv. 1991;57(3):315–37.

    Article  Google Scholar 

  11. 11.

    Ledger JA, Annegarn HJ. Electrocution hazards to the cape vulture Gyps coprotheres in South Africa. Biol Conserv. 1981;20(1):15–24.

    Article  Google Scholar 

  12. 12.

    Mundy PJ, Butchart D, Ledger J, Piper S. The vultures of Africa. London: Academic Press; 1992. p. 464.

    Google Scholar 

  13. 13.

    Ogada D, Shaw P, Beyers RL, Buij R, Murn C, Thiollay JM, et al. Another continental vulture crisis: Africa’s vultures collapsing toward extinction. Conservation Letters. 2016;9(2):89–97.

    Article  Google Scholar 

  14. 14.

    Ogada DL. The power of poison: pesticide poisoning of Africa’s wildlife. Ann N Y Acad Sci. 2014;1322(1):1–20.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Ogada DL, Keesing F, Virani MZ. Dropping dead: causes and consequences of vulture population declines worldwide. Ann N Y Acad Sci. 2012;1249(1):57–71.

    PubMed  Article  Google Scholar 

  16. 16.

    Krüger SC, Amar A. Lead exposure in the critically endangered bearded vulture (Gypaetus barbatus) population in southern Africa. J Raptor Res. 2018;52(4):491–9.

    Article  Google Scholar 

  17. 17.

    Godino A, Paz JL, Simón MÁ. Naturalistas espanoles localizan en Marruecos cinco quebrantahuesos. Quercus. 2003;205:46–7.

    Google Scholar 

  18. 18.

    Brown L. The status, population structure and breeding dates of the African Lammergeier Gypaetus barbatus meridionalis. Raptor Res. 1977;11(3):49–58.

    Google Scholar 

  19. 19.

    Ash J, Atkins J. Birds of Ethiopia and Eritrea: an atlas of distribution. New York: Bloomsbury Publishing; 2010.

    Google Scholar 

  20. 20.

    Krüger SC, Allan DG, Jenkins AR, Amar A. Trends in territory occupancy, distribution and density of the Bearded Vulture Gypaetus barbatus meridionalis in southern Africa. Bird Conserv Int. 2014;24(2):162–77.

    Article  Google Scholar 

  21. 21.

    Allendorf FW. Genetic drift and the loss of alleles versus heterozygosity. Zoo Biol. 1986;5(2):181–90.

    Article  Google Scholar 

  22. 22.

    Frankham R. Inbreeding and extinction: a threshold effect. Conserv Biol. 1995;9(4):792–9.

    Article  Google Scholar 

  23. 23.

    Lacy RC. Loss of genetic diversity from managed populations: interacting effects of drift, mutation, immigration, selection, and population subdivision. Conserv Biol. 1987;1(2):143–58.

    Article  Google Scholar 

  24. 24.

    Vucetich JA, Waite TA. Erosion of heterozygosity in fluctuating populations. Conserv Biol. 1999;13(4):860–8.

    Article  Google Scholar 

  25. 25.

    Frankham R. Genetics and extinction. Biol Conserv. 2005;126(2):131–40.

    Article  Google Scholar 

  26. 26.

    Tollington S, Jones CG, Greenwood A, Tatayah V, Raisin C, Burke T, et al. Long-term, fine-scale temporal patterns of genetic diversity in the restored Mauritius parakeet reveal genetic impacts of management and associated demographic effects on reintroduction programmes. Biol Conserv. 2013;161:28–38.

    Article  Google Scholar 

  27. 27.

    Frankham R, Bradshaw CJ, Brook BW. Genetics in conservation management: revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. Biol Conserv. 2014;170:56–63.

    Article  Google Scholar 

  28. 28.

    Orta J, Juana E, Marks JS, Sharpe CJ, Garcia EFJ. Bearded Vulture (Gypaetus barbatus), version 1.0. In: Hoyo J, Elliott A, Sargatal J, Christie DA, Juana E, editors. Birds of the World. Ithaca: Cornell Lab of Ornithology; 2020.

  29. 29.

    Krüger SC, Wesche PL, Jansen van Vuuren B. Reduced genetic diversity in Bearded Vultures Gypaetus barbatus in Southern Africa. Ibis. 2015;157(1):162–6.

  30. 30.

    Gautschi B, Jacob G, Negro JJ, Godoy JA, Müller JP, Schmid B. Analysis of relatedness and determination of the source of founders in the captive bearded vulture, Gypaetus barbatus, population. Conserv Genet. 2003;4(4):479–90.

    CAS  Article  Google Scholar 

  31. 31.

    Carlsson J. Effects of microsatellite null alleles on assignment testing. J Hered. 2008;99(6):616–23.

    CAS  Article  Google Scholar 

  32. 32.

    Wilcox BA, Murphy DD. Conservation strategy: the effects of fragmentation on extinction. Am Nat. 1985;125(6):879–87.

    Article  Google Scholar 

  33. 33.

    Lande R, Barrowclough GF. Effective population size, genetic variation, and their use in population management. In: Soulé M, editor. Viable Populations for Conservation. 87. Cambridge: Cambridge University Press; 1987. p. 87–124.

  34. 34.

    Gilpin M. Forty-eight parrots and the origins of population viability analysis. Conserv Biol. 1996;10(6):1491–3.

    Article  Google Scholar 

  35. 35.

    Shaffer M. Minimum viable populations: coping with uncertainty. In: Soulé M, editor. Viable populations for conservation. Cambridge: Cambridge University Press; 1987. p. 69–86.

    Google Scholar 

  36. 36.

    Bretagnolle V, Inchausti P, Seguin J-F, Thibault J-C. Evaluation of the extinction risk and of conservation alternatives for a very small insular population: the bearded vulture Gypaetus barbatus in Corsica. Biol Conserv. 2004;120(1):19–30.

    Article  Google Scholar 

  37. 37.

    Soulé ME. The fitness and viability of populations. In: Soulé ME, editor. Conservation biology: the science of scarcity and diversity. Sinauer Associates: Sunderland, Massachusetts; 1986. p. 13–116.

  38. 38.

    Haig SM, Gratto-Trevor CL, Mullins TD, Colwell MA. Population identification of western hemisphere shorebirds throughout the annual cycle. Mol Ecol. 1997;6(5):413–27.

    CAS  Article  Google Scholar 

  39. 39.

    Wenink PW, Baker AJ, Tilanus M. Hypervariable-control-region sequences reveal global population structuring in a long-distance migrant shorebird, the Dunlin (Calidris alpina). Proc Natl Acad Sci. 1993;90(1):94–8.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    López-López P, Zuberogoitia Í, Alcántara M, Gil JA. Philopatry, natal dispersal, first settlement and age of first breeding of Bearded Vultures Gypaetus barbatus in central Pyrenees. Bird Study. 2013;60(4):555–60.

    Article  Google Scholar 

  41. 41.

    IOC World Bird List (v 10.2). 2020. Accessed 01 Dec 2020

  42. 42.

    Dickinson E, Christidis L. The Howard and Moore complete checklist of the birds of the world. Berlin: Springer; 2014. p. 461.

    Google Scholar 

  43. 43.

    Hiraldo F, Delibes M, Calderon J. Comments on the taxonomy of the Bearded Vulture Gypaetus barbatus (Linnaeus, 1758). Bonner Zoologische Beiträge. 1984;35:91–5.

    Google Scholar 

  44. 44.

    The eBird/Clements Checklist of Birds of the World: v2019. 2019. Accessed 01 Dec 2020.

  45. 45.

    Coetzee J. African flora since the terminal Jurassic. In: Goldblatt P, editor. Biological relationships between Africa and South America. 1st ed. New Haven: Yale University Press; 1993. p. 37–61.

    Google Scholar 

  46. 46.

    Hartl DL, Clark AG. Principles of population genetics. Sunderland: Sinauer associates; 1998.

    Google Scholar 

  47. 47.

    Slatkin M. Gene flow and the geographic structure of natural populations. Science. 1987;236(4803):787–92.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Krüger S, Piper S, Rushworth I, Botha A, Daly B, Allan D, et al. Bearded Vulture (Gypaetus barbatus meridionalis) population and habitat viability assessment workshop report. Conservation Breeding Specialist Group (SSC/IUCN)/CBSG Southern Africa Endangered Wildlife Trust: Johannesburg; 2006.

    Google Scholar 

  49. 49.

    Negro JJ, Torres MJ. Genetic variability and differentiation of two bearded vulture Gypaetus barbatus populations and implications for reintroduction projects. Biol Conserv. 1999;87(2):249–54.

    Article  Google Scholar 

  50. 50.

    Kleinhans C, Willows-Munro S. Low genetic diversity and shallow population structure in the endangered vulture Gyps coprotheres. Sci Rep. 2019;9(1):1–11.

    CAS  Article  Google Scholar 

  51. 51.

    Mira S, Billot C, Guillemaud T, Palma L, Cancela M. Isolation and characterization of polymorphic microsatellite markers in Eurasian vulture Gyps fulvus. Mol Ecol Notes. 2002;2(4):557–8.

    CAS  Article  Google Scholar 

  52. 52.

    Chapuis M-P, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24(3):621–31.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Resour. 2006;6(1):288–95.

  54. 54.

    Goudet JF. FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995;86(6):485–6.

  55. 55.

    Raymond M, Rousset F. GENEPOP Version 1.2: Population genetics software for exat tests and ecumenicism. J Hered. 1995;86:248–9.

  56. 56.

    Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.

    CAS  PubMed  Google Scholar 

  57. 57.

    Weir BS. Genetic data analysis. Methods for discrete population genetic data. Sunderland: Sinauer Associates; 1990. p. 377.

    Google Scholar 

  58. 58.

    Raymond M, Rousset F. An exact test for population differentiation. Evolution. 1995;49(6):1280–3.

    PubMed  Article  Google Scholar 

  59. 59.

    Falush D, Stephens M. Pritchard JK. Inference of Population Structure Using Multilocus Genotype Data: Linked Loci and Correlated Allele Frequencies. Genetics; 2003. p. 164.

    Google Scholar 

  60. 60.

    Pritchard JK, Stephens M, Donelly P. Inference of population structure using multilocus genotype data. Genetics Soc Am. 2000;155:945–59.

    CAS  Google Scholar 

  61. 61.

    Li YL, Liu JX. StructureSelector: a web based software to select and visualize the optimal number of clusters by using multiple methods. Mol Ecol Resour. 2017;18(1):176–7.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Puechmaille SJ. The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem. Mol Ecol Resour. 2016;16(3):608–27.

    PubMed  Article  Google Scholar 

  63. 63.

    Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genetics Resour. 2012;4(2):359–61.

    Article  Google Scholar 

  64. 64.

    Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.

    CAS  Article  Google Scholar 

  65. 65.

    Francis RM. Pophelper: an r package and web app to analyse and visualize population structure. Mol Ecol Resour. 2016;17(1):27–32.

    PubMed  Article  CAS  Google Scholar 

  66. 66.

    Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27(2 Part 1):209–20.

  67. 67.

    Kivelä M, Arnaud-Haond S, Saramäki J. EDEN etworks: A user-friendly software to build and analyse networks in biogeography, ecology and population genetics. Mol Ecol Resour. 2015;15(1):117–22.

    PubMed  Article  Google Scholar 

  68. 68.

    Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163(3):1177–91.

    PubMed  PubMed Central  Google Scholar 

  69. 69.

    Hastings A. Complex interactions between dispersal and dynamics: lessons from coupled logistic equations. Ecology. 1993;74(5):1362–72.

    Article  Google Scholar 

Download references


The authors wish to thank David Allan from the Durban Natural Sciences Museum (DNSM, Durban, South Africa), Hein Grouw from the Natural History Museum, Tring (United Kingdom), Hans Pohlman from the Naturalis Biodiversity Centre (Leiden, Netherlands) and Paul Sweet from the American Museum of Natural History (AMNH, United States) for providing access to archival samples. Field-collected samples were donated by Sonja Krüger.


Running expenses for this study were funded by a Rated Researcher grant awarded to S. Willows-Munro by the National Research Foundation (South Africa). The first author received a free-standing MSc bursary from the NRF (SFH170613240025,) and a GC Weightman Scholarship (South Africa) from the University of KwaZulu-Natal. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




MS and SWM conceived and designed the study. FL and SK provided samples for the study. MS and SWM generated, processed, and analyzed molecular data. MS drafted the initial version of the manuscript. All authors read, edited, enhanced, and approved the final version of the manuscript.

Corresponding author

Correspondence to Sandi Willows-Munro.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained for this study from the University of KwaZulu-Natal Animal Ethics subcommittee (Reference number: AREC/054/016 M).

Consent for publication

Not applicable.

Competing interests

There are no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary data for Bearded Vulture analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Streicher, M., Krüger, S., Loercher, F. et al. Evidence of genetic structure in the wide-ranging bearded vulture (Gypaetus barbatus (Linnaeus, 1758)). BMC Ecol Evo 21, 42 (2021).

Download citation


  • Bearded vulture
  • Gene flow
  • Habitat fragmentation
  • Gypaetus barbatus
  • Population decline
  • Genetic diversity