Zoosyst. Evol. 97 (1) 2021, 211-221 | DOI 10.3897/zse.97.61854 Zieh i 40h Bet ah eeT WE) Bape _., D> PENSOFT. Nitze A NATURKUNDE BERLIN Demography reveals populational expansion of a recently extinct Iberian ungulate Giovanni Forcina!, Kees Woutersen?, Santiago Sanchez-Ramirez?, Samer Angelone’, Jean P. Crampe°, Jesus M. Pérez®’, Paulino Fandos’®, José Enrique Granados’:?, Michael J. Jowers!'1° 1 CIBIO/InBIO, Centro de Investigagao em Biodiversidade e Recursos Genéticos, Universidade do Porto, Campus Agrario De Vairdo, 4485-661 Vairdo, Portugal Calle Ingeniero Montaner, 4-1-C 22004 Huesca, Spain Department of Ecology and Evolutionary Biology, University of Toronto, 25 Willcocks, Toronto, Ontario M5S 3B2, Canada Institute of Evolutionary Biology and Environmental Studies (IEU), University of Zurich, Winterthurerstrasse 190, Zurich, Switzerland 19 chemin de Peyborde 65400 Lau-Balagnas, France Departamento de Biologia Animal, Biologia Vegetal y Ecologia, Universidad de Jaén, Campus Las Lagunillas, s.n., E-23071 Jaén, Spain Wildlife Ecology & Health group (WE&H), Jaén, Spain Agencia de Medio Ambiente y Agua, E-41092 Sevilla, Isla de la Cartuja, Spain O ON DO FW DY Espacio Natural Sierra Nevada, Carretera Antigua de Sierra Nevada, Km 7, E-18071 Pinos Genil, Granada, Spain 10 National Institute of Ecology, 1210, Geumgang-ro, Maseo-myeon, Seocheon-gun, Chungcheongnam-do, 33657, South Korea http://zoobank. org/38 2480DE-1CE7-45A7-9C75-21036C7A3560 Corresponding author: Michael J. Jowers (michaeljowers@hotmail.com) Academic editor: M. TR Hawkins # Received 11 December 2020 # Accepted 9 March 2021 Published | April 2021 Abstract Reconstructing the demographic history of endangered taxa is paramount to predict future fluctuations and disentangle the contribut- ing factors. Extinct taxa or populations might also provide key insights in this respect by means of the DNA extracted from museum specimens. Nevertheless, the degraded status of biological material and the limited number of records may pose some constraints. For this reason, identifying all available sources, including private and public biological collections, is a crucial step forward. In this study, we reconstructed the demographic history based on cytochrome-b sequence data of the Pyrenean ibex (Capra pyrenaica pyre- naica), a charismatic taxon of the European wildlife that became extinct in the year 2000. Moreover, we built a database of the mu- seum specimens available in public biological collections worldwide and genotyped a privately owned 140-year-old trophy from the Spanish Pyrenees to confirm its origin. We found that the population of the Pyrenean ibex underwent a recent expansion approximate- ly 20,000 years ago, after which trophy hunting and epizootics triggered a relentless population decline. Our interpretations, based on the genetic information currently available in public repositories, provide a solid basis for more exhaustive analyses relying on all the new sources identified. In particular, the adoption of a genome-wide approach appears a fundamental prerequisite to disentangle the multiple contributing factors associated with low genetic diversity, including inbreeding depression, acting as extinction drivers. Key Words ancient DNA, anthropocene biodiversity crisis, biotic impoverishment, Capra pyrenaica pyrenaica, Cytochrome-b, epizootics, extinction, museum specimens, public repositories, trophy hunting Copyright Giovanni Forcina et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 22 Introduction The Anthropocene biodiversity crisis 1s causing a dras- tic biotic impoverishment worldwide (Leakey and Lewin 1996; Barnosky et al. 2011; Pievani 2014). Characteris- ing genetic variation of endangered species is key to re- construct and predict past and future demographic fluc- tuations (Collen et al. 2011). Such knowledge may shed light on the proximate and ultimate causes underlying the collapse of extinct taxa and of their populations, eluci- dating pivotal information for highly endangered species that may undergo a similar fate (Farrington et al. 2019). The Iberian ibex (Capra pyrenaica) includes two ex- tant (C. p. hispanica and C. p. victoriae) and two extinct subspecies (C. p. /usitanica and C. p. pyrenaica) (Herrero et al. 2020) distinguished on the basis of pelage charac- teristics and horn development (Cabrera 1911). Among the latter, C. p. pyrenaica (Fig. 1a), commonly known as the Pyrenean ibex or the bucardo, inhabited the Pyrenees until the year 2000 (Pérez et al. 2002). Morphological studies indicate this taxon as being the largest and stoutest of its group (with special reference to females), with cra- nial measurements placing it in an intermediate position between the Alpine ibex (C. ibex) and the other C. pyre- naica subspecies (Garcia-Gonzalez 1991, Granados et al. 1997). The distinctiveness of C. p. pyrenaica was later confirmed by genetic studies employing mitochondrial (Villalta et al. 1997, Manceau 1999) and microsatellite (Jiménez et al. 1999) markers pointing to a divergence between this subspecies and the others which is compara- ble to that between the latter and the Alpine ibex. Now extinct due to human activities, C. p. pyrenaica was, and still 1s, one of Europe’s most charismatic an- imals. Recently, interest in this subspecies has surged enough to motivate researchers to a controversial cloning attempt (Folch et al. 2009; Kupferschmidt 2014). Initia- tives to re-establish populations of the Iberian ibex have been taken across the Iberian Peninsula, where a popula- tion of the subspecies C. p. victoriae now occurs in the former range of C. p. lusitanica (Moco et al. 2006, 2015). Plans to relocate C. p. victoriae also into the historical range of C. p. pyrenaica (Crampe 1991) have been en- acted by releasing and monitoring of animals since 2014 onwards (Crampe et al. 2015). This population now con- sists of 250 individuals in the Pyrénées National Park and over 500 in the entire mountain chain (http://www.bou- quetin-pyrenees.fr/). Medieval sources such as Gaston Phébus’ hunting chronicles (1387-1389) revealed that the Pyrenean ibex was abundant (Astre 1952), but this iconic animal became a common target for trophy hunters during the 19" and 20" century, which triggered a relentless decline reported in detail by several authors (e.g. Gourdon 1908, Cabrera 1914, Labarere 1985, Crampe et Crégut-Bonnoure 1994, Garcia-Gonzalez and Herrero 1999, Jiménez 2016). More- over, infectious diseases like the sarcoptic mange have also been listed among the causes of its disappearance (Perez et al. 2002, Acevedo and Cassinello 2009), while zse.pensoft.net FORCINA, G. et al.: Extinction of an Iberian ungulate pathogens like the bluetongue virus (BTV), which has recently caused the sharp demographic collapse of other ungulates (e.g. the pronghorn antelope Antilocapra amer- icana: Thorne et al. 1988, the goitred gazelle Gazella sub- gutturosa subgutturosa. Gir 2008), have been recurrently detected in Iberian ibex populations across Spain (Garcia et al. 2009; Lorca-Oro et al. 2014; GOmez-Guillamon et al. 2020). Considered a serious threat to wildlife at continental level (Rossi et al. 2019), it is likely that this virus (most often transmitted from livestock) also contributed to the vertiginous decline of the Pyrenean ibex. Despite the fact that hunting and epizootics from sympatric livestock are deemed to have played a major role (along with habitat fragmentation) in drastically reducing its populations over the last two centuries, the relative contribution of multiple threatening factors remains largely unknown (Hidalgo and Garcia-Gonzalez 1995; Garcia-Gonzalez et al. 1996). Nonetheless, there is strong evidence that the severe bottleneck caused by hunting amplified high lev- els of inbreeding and homozygosis (Jiménez et al. 1999) which likely reduced fertility of the fast-shrinking pop- ulation brought beyond the minimum viable population size (Acevedo and Cassinello 2009) into an extinction vortex (Soulé 1987). In line with these considerations, a recent genome-wide study of all European ibex species evidenced a pattern of low genetic diversity and high inbreeding (Grossen et al. 2018), consistent with an overall well-defined genet- ic structure following genetic drift after isolation (Ange- lone-Alasaad et al. 2017). The scenario portrayed for the Pyrenean ibex would be similar to that unveiled for the cheetah (Acinonyx jubatus), in which minimal genetic diversity was correlated with major bottlenecks occurred ca. 100,000 (Dobrynin et al. 2015) and then again 12,000 years ago (O’Brien et al. 1985, Driscoll et al. 2002), and for both the Eurasian lynx (Lynx /ynx) and Iberian lynx (Lynx pardinus), which experienced a drastic decline 700,000— 100,000 thousand years ago (Abascal et al. 2016). Among ungulates, recurrent and non-anthropogenic bottlenecks have been indicated as the cause lessening genetic diversi- ty in the muskox (Ovibos moschatus) (Prewer et al. 2020). In order to shed some light on the demographic his- tory of the Pyrenean ibex and possibly get insights into the drivers underlying its extinction, we reconstructed its demographic history relying on mitochondrial DNA sequences (cytochrome-b: cyt-b) currently accessible in public repositories. We also built a database of extant Pyrenean ibex museum specimens, most of which have not been used for genetic analyses, and we genetically confirmed the identity of an additional record from a pri- vate collection. Our results suggest that the population of the Pyrenean ibex underwent a recent expansion ap- proximately 20,000 years ago, while trophy hunting and epizootics triggered a relentless decline over the last two centuries. These inferences provide a solid basis for more exhaustive analyses relying on all the new sources iden- tified and a genome-wide approach bearing potential to Zoosyst. Evol. 97 (1) 2021, 211-221 213 Figure 1. C. p. pyrenaica: a. Female individual (visitor centre of Ordesa y Monte Perdido National Park, Torla-Ordesa, Huesca, Spain). Photo: courtesy of Jose Miguel Pintor Ortego. b. Picture of the C. p. pyrenaica trophy of the 3-year-old male genotyped in this study and preserved at the English Circle of Pau (France). Photo: Jean P. Crampe. assess the role of multiple contributing factors underlying the extinction of this animal. We finally discuss our re- sults and compare them to those of other ungulates which recently experienced sharp population declines or even species extinction, and we address the key points for fur- ther studies to consider on the Pyrenean ibex. Materials and methods Demographic inferences We used GenBank C. p. pyrenaica cyt-b sequences (alignment length 1,140 bps) and run preliminary phy- logenetic trees to test for monophyly. This locus was chosen as the one with the highest number of C. p. pyre- naica records available in public repositories. We used the samples obtained by Urefia et al. (2018), which show clear monophyly and two additional samples from the Pyrenees (Manceau et al. 1999). The final data set con- sisted of 17 records (Fig. 2; Suppl. material 1: Table S1). Sequences were translated to amino acids to check against stop codons and aligned in SEAVIEW v4.2.11 (Gouy et al. 2010) under CLUSTALW2 default settings (Larkin et al. 2007). We used JMODELTEST (Posada 2008) to choose the optimal model. We specified a relaxed lognormal clock for the cyt-b data in BEAST v1.8.2 (Drummond et al. 2012, http://beast.bio.ed.ac.uk) and the nucleotide substi- tution model HKY as suggested by JMODELTEST. To time-calibrate the population tree, we used fossil dates as tip calibrations as reported in C. p. pyrenaica (Urefia et al. 2018). These dates were between GModerrand and Ker- aval 1 at 16,294 years ago, between Valdegoba 1 and all other C. p. pyrenaica at 40,093 years ago, and the clade composed by Chaves 56, 42, 61, 64, 63, Bolinkoba 3 and Chorrugues 1 at 29,386 years ago. We used a 10% stand- ard deviation for all calibrations. Changes in effective population size through time were inferred using a Bayes- ian skyline plot (BSP) model (Drummond et al. 2005). We conducted two independent Markov Chain Monte Carlo (MCMC) runs, each with 30 million states and sampling every 3,000" state. Independent runs were eval- uated for convergence and mixing by observing and com- paring traces of each statistic and parameters in TRACER v1.6 (Rambaut et al. 2007, http://beast.bio.ed.ac.uk/trac- er). We considered effective sampling size (ESS) values > 200 to be good indicators of parameter mixing. The first 10% of each run were discarded as burn-in, and samples were merged using LogCombiner v1.8.2 (Drummond et al. 2012). All BEAST analyses were performed through the CIPRES platform (Miller et al. 2010). In order to test for past population demographics, we used Tajima’s D (Tajima 1989) and Fu’s FS (Fu 1997) statistics. Demographic changes were also examined by calculating Harpending’s raggedness index (Harpending et al. 1993) and the sum of squared deviations (SSD) (Schneider and Excoffier 1999). All neutrality tests were calculated in ARLEQUIN v3.5 (Excoffier and Lischer 2010) under a demographic expansion model and 1,000 simulations. The results of mismatch distribution analy- ses were plotted in DNAsp v5.10 (Rozas 2009). zse.pensoft.net 214 FORCINA, G. et al.: Extinction of an Iberian ungulate Figure 2. C. p. pyrenaica hypothetical distribution during the 19" (dark grey) and at the beginning of the 20" (black) century, when this taxon occurred only in Ordesa Valley and La Maladeta Massif (east to west: Woutersen 2019). Numbers refer to sampling localities in Suppl. material 1: Table S1. Identification, origin and provision of museum Specimens We made a thorough bibliographic search to retrieve information on C. p. pyrenaica and built a comprehen- sive database of historical specimens preserved in mu- seums and smaller collections, both public and private, thus localising new valuable source of information to elucidate the natural history of this iconic animal. For this purpose, we first relied on literary sources (includ- ing letters exchanged by trophy hunters and museum curators) presented in a previous review on the Pyre- nean ibex (Woutersen 2012, 2019) and then expanded our search by using the platform VertNet (http://vert- net.org/). Hence, we examined the sex ratio and age of all known specimens to assess whether the occurrence of adult males, persistently persecuted for hunting, de- creased across time. For this purpose, we plotted ages of male and female specimens versus collection date to detect the occurrence of possible sex-biased discrepan- cies. Only a subset (n = 32) of the entire sample was used as some records were data deficient and were therefore excluded from analyses. For three specimens we lacked ages, but they were catalogued as ‘young’ and hence incorporated in the graph as juveniles (ca. 3 years of age). Specimens with dates ‘1835-1836’ were included in the analyses as 1836, those with ‘1837 or later’ as 1837 and those with ‘1878 or later’ as 1878 (Suppl. material 1: Table S2). In addition, we examined the sex ratio and age of all specimens in the 1800s and zse.pensoft.net 1900s to assess for possible differences in collection re- cords, and particularly whether the occurrence of adult males decreased across time. Genetic analyses We obtained a cranium from a C. p. pyrenaica trophy of a 3-year-old male preserved in a private collection (Fig. 1b) and hunted by British hunter Arthur Post on 10" April 1881 (Leslie 1894) in the Ordesa Valley (Spanish Pyr- enees), now part of the Ordesa and Monte Perdido Na- tional Park, a protected area instituted in 1918 to rescue the Pyrenean ibex from extinction. DNA extraction was performed using the DNeasy Blood and Tissue Kit (Qia- gen, Hilden, Germany) as in Forcina et al. (2015) at the National Institute of Ecology (NIE: Gunsan, South Ko- rea) by strictly adhering to the protocol for nucleic acid isolation from archival samples. A 310 bp-long Control Region (CR) fragment of the mitochondrial DNA (mtD- NA) was amplified in a single PCR with thermal condi- tions detailed in Forcina et al. (2015) using the primers D-loop-FW (GATCCCTCTTCTCGCTCCG) = and D-loop-cabra (CCATGCCTACCATTATGGGGA) (Amills et al. 2004a). Amplicon purification was per- formed with the Genelute PCR Clean-up Kit (final vol- ume 40 uwL; Sigma Aldrich) and directly sequenced on both DNA strands (BigDye Terminator v3.1 Cycle Sequencing Kit, ABI 3730 DNA automated sequenc- er, Applied Biosystems). We then used the Basic Local Zoosyst. Evol. 97 (1) 2021, 211-221 a: 10,000 1,000 N (log we 100 215 10 c EXPANSION. - 0 5,000 10,000 15,000 20,000 25,000 30,000 35,000 40,000 Time Exp ---@- Obs 15 20 25 Pairwise differences Figure 3. Demographic analyses: a. Bayesian skyline plot (Ne values were rescaled based on a generation time of 10 years and a 1/4 conversion factor for mtDNA data); b. Mismatch distributions. Exp: expected; Obs: observed. Alignment Search Tool (BLAST) as implemented in the National Center of Bioinformatic Information (NCBI) website to confirm species identity. The newly produced CR sequence was deposited in GenBank (accession num- ber: MW659860). Results Demographic inferences No stop codons were found. The runs showed high Effec- tive Sample Size (>1,500), corresponding to an adequate sampling of the posterior distribution. BSPs indicated a sta- ble population within the last 15,000 years and a population expansion between 15,000 and 30,000 years ago (Fig. 2a). The mismatch distribution analysis shows a bell shape dis- tribution suggesting population expansions (Fig. 2b). Negative neutrality tests, Fu’s Fs (-4.38, p = 0.01) and Tajima’s D (-1.94, p = 0.012) statistical values observed in the population supports for evidence of past popula- tion expansions. Non-significant values for SSD (0.008, p = 0.71) indicate that the data do not deviate from those expected under the model of expansion. Similar- ly, non-significant raggedness values (0.0296, p = 0.71) indicate population expansion. Non-significant value in goodness-of-fit distribution for all populations suggest that this phenomenon is likely recent (Rogers 1995). The gene diversity (0.95 + 0.043) and the mean number of pairwise differences (3.92 + 2.08) point to high popula- tion differentiation. Identification, origin and provision of museum specimens We located 45 specimens of C. p. pyrenaica collected between 1818 and 2000. Of these, 43 came from two massifs in the Spanish Pyrenees, La Maladeta and Monte Perdido (which includes the Ordesa Valley), while two more records from nearby sites in France were found through VertNet at the LACM Vertebrate Collection Nat- ural (History Museum of Los Angeles County) and the Cowan Tetrapod Collection (University of British Co- lumbia Beaty Biodiversity Museum) (Suppl. material 1: Table S2). Archival sources revealed that the specimens preserved in the biological collections (which in the 19" century were still mostly private) of Toulouse, Zurich, zse.pensoft.net 216 Mainz, St. Petersburg, Strasbourg and Vienna had been specifically bought for the purpose of exhibiting animals killed on specific request (Woutersen 2019). For instance, around 1835, several collectors of zoological specimens in Germany and Switzerland asked the naturalist Alfred Moquin-Tandon from Toulouse (France) for C. p. pyrena- ica records. The latter commissioned professional hunters in Benasque, a village at the feet of Maladeta Massif, to kill some individuals for the provision of museum speci- mens. Significantly, the documents also reveal that more was paid for males than for females. In 1841, the exis- tence of the population of the Ordesa Valley (then called Val d’Arras) attracted trophy hunters who ended up en- larging incidentally the list of C. p. pyrenaica specimens preserved in natural history museums worldwide until 1913, when the Pyrenean ibex was lawfully protected. Attributes of museum specimens The origin, collection date, sex and other features of C. p. pyrenaica specimens are listed in Suppl. material 1: Ta- ble S2. Two thirds of the animals currently available in museums come from La Maladeta massif and were killed between 1818 and 1892. Those hunted on Monte Perdido were killed between 1852 and 1910. The first of the eight specimens collected by the staff of Ordesa and Monte Perdido National Park is from 1958. Males collected both only prior to (71.4%, n= 25) and also after (73.3%, n= 33) the 1913 hunting ban were more numerous than females (Fig. 4). We can hence conclude that males were probably hunted more than females for their trophies before 1913, and that museum curators were perhaps more willing to secure male rather than female specimens for their collec- tion. The average age of individuals collected from 1800— 1900 (n = 25) was 7.9 years, while that for those collected from 1900—2000 (n = 7) reached 12.1 years. Age 3S 1805 1825 1845 1865 1885 1905 1925 1945 1965 1985 2005 Year Figure 4. Scatter plot representation of museum specimen ages and year of collection over time with an emphasis on before and after the 1913 hunting ban (as indicated by the gun icon). The trend line is r? = 0.21. Blue and purple stars are used for males and females, respectively. zse.pensoft.net FORCINA, G. et al.: Extinction of an Iberian ungulate Genetic analyses The NCBI-blast revealed that the C. p. pyrenaica CR frag- ment was 100% identical to the CR isolate CPP1 obtained from the last living individual of C. p pyrenaica, which died in Ordesa y Monte Perdido National Park (Central Pyrenees, Huesca, Spain) in the year 2000 (Folch et al. 2009), only two kilometres away from where the 3-year old male had been hunted 140 years before. Moreover, the CR fragment differed by two to six nucleotides when compared to the homologous regions of GenBank re- cords AJ233741 and AJ233739—AJ233740 (Manceau et al. 1999), respectively, which were also retrieved from animals belonging to the same Spanish population. Discussion The Pyrenean ibex is an emblematic example of wildlife that went extinct in recent times in spite of its popularity among the broader public and evidently belated conser- vation efforts that started only a decade before the last in- dividual died (e.g., Garcia-Gonzalez and Herrero 1999). Among European ungulates, the last two centuries have witnessed the extinction of other taxa such as the Cauca- sian moose (Alces alces caucasicus), the Caucasian wis- ent (Bison bonasus caucasicus), the tarpan (Equus ferus ferus) and the Portuguese ibex (C. p. /usitanica). However, in more recent times captive breeding programs have pre- vented the extinction of other taxa as with the Corsican red deer (Cervus elaphus corsicanus, Hmwe et al. 2006, Kidjo et al. 2007) and, outside Europe, with the scimitar-horned oryx (Oryx dammah, Little et al. 2016) and Pere David’s deer (Elaphurus davidianus, Zeng et al. 2007). Admittedly, several studies have been carried out to deepen our knowl- edge of the Pyrenean ibex over the last two decades of the 20" century, including two to assess the minimum vi- able population size (Garcia-Gonzalez et al. 1991, 1996), which nonetheless were futile in rescuing the subspecies from the brink of extinction. However, no demographic re- constructions had been performed to understand historical population sizes of the Pyrenean ibex thousands of years ago. In order to fill this knowledge gap, we infer temporal changes in its effective population size by means of the BSP, mismatch distributions, and neutrality tests. Our analyses evidenced a population expansion during the Marine Isotope Stage 2 (MIS 2) period (14,000—29,000 years ago), a time marked by abrupt climatic oscillations in Southern Europe (Last Glacial Maximum: 20,000 years ago). The BSP shows a contrasting trend to other critical- ly endangered ungulates (which, albeit, live in different ecological contexts and latitudes) such as the hirola (Be- atragus hunteri, Jowers et al. 2020) and the saiga ante- lope (Saiga tatarica, Campos et al. 2010). Both species underwent a genetic bottleneck as expected after a severe demographic contraction. DNA data from ancient and modern samples of the saiga antelope have shown increased haplotype diversity (with similar values to C. p. pyrendica) in ancient times and a Zoosyst. Evol. 97 (1) 2021, 211-221 reduction of diversity in modern times. Concordantly, the BSP analyses of combined samples reveals a Pleistocene/ Holocene bottleneck (Campos et al. 2010), a pattern not observed in C. p. pyrenaica. On the other hand, Urefia et al (2011) detected a significant loss of genetic diversity in samples of this taxon from Cueva de Chaves. This trend started between the Paleolithic and Neolithic (14,932 and 7,504 yrs BP: Urefia et al. 2018) and persisted to the present day (Manceau et al. 1999, Amills et al. 2004b). In accord- ance to the observed BSP expansion, the genetic differentia- tion as inferred from haplotypic divergence and genetic dis- tance was still high in comparison to other ungulate species with larger distribution ranges, such as the dorcas gazelle (Gazella dorcas, Godinho et al. 2012), the scimitar-horned oryx (Iyengar et al. 2007), the roan antelope (Hippotragus equinus, Alpers et al. 2004) or the recently extinct in the wild Saudi gazelle (Gazella saudiya, Hammond et al. 2001). Although we cannot rule out that the comparatively high- er genetic diversity is at least partly due to the heterogeneity of the samples used, this outcome certainly deserves inter- est and should be readdressed in further studies relying on a larger dataset of samples from the last two centuries. By pointing to a demographic expansion of a still genetically diverse population over the last 15,000—30,0000 years, these results do not allow us to draw conclusions about the pri- mary cause underlying the decline of the Pyrenean ibex nor which anthropogenic impact (including hunting but also dis- eases from livestock) played a major role for its extinction. Indeed, previous studies acknowledge that hunting, e1- ther legal or illegal, was extremely detrimental not only in terms of the reduction in population size, but also in terms of the disruption of the social structure of the population (Garcia-Gonzalez et al. 1996; Garcia-Gonzalez and Herrero 1999; Acevedo and Cassinello 2009; Herrero et al. 2020). However, multiple contributing factors, such as epizootics, are typically invoked. In this respect, it is worth mentioning that epidemics have been known to wipe out populations of this species over extremely short periods. For example, the sarcoptic mange affecting the C. p. hispanica popula- tion inhabiting the Sierra de Cazorla (Andalucia, Southern Spain) during the 1980s caused the disappearance of 90% of the population (10,000 individuals) in just a few years (Granados et al. 2001). According to Amills et al. (2004b), the bottleneck experienced by the Pyrenean ibex and the consequent loss of genetic diversity turned into depleted variation at genes of the Major Histocompatibility Com- plex (MHC), thus enhancing susceptibility to epidemics. On the other hand, the extinction of several other un- gulates since the 1800s has ultimately been linked to in- discriminate poaching. This is surely the case of the Saudi gazelle, the South African bluebuck (Hippotragus leuco- phaeus), the Queen of Sheba’s gazelle or Yemen gazelle (Gazella bilkis), and the Schomburgk’s deer (Rucervus schomburgki), while other heavily hunted taxa such as the Barbary deer (Cervus elaphus barbarous: Kingdon et al. 2013), the Alpine ibex (Grossen et al. 2017) and the Ap- ennine chamois (Rupicapra pyrenaica ornata: Pérez et al. 2014) were rescued from the brink of extinction thanks to last resort conservation initiatives. 21g Although the ultimate culprit factor responsible for a species extinction remains challenging to ascertain, it is reasonable to argue that multiple confounding anthro- pogenic causes likely contributed to the demise of these populations. In other cases, the occurrence of such det- rimental factors is widely acknowledged. For instance, the ure (Bos primigenius) was believed to have become extinct in the 16" century by unrestricted hunting, habitat reduction due to expanding farmlands, and diseases trans- mitted by domesticated cattle (Rokosz 1995; Van Vuure 2005). Likewise, the drastic decline in bluebuck popula- tions just before the 16" century has been largely attrib- uted to competition with livestock for forage and habitat deterioration (Klein 1974), as was the fate of the Quag- ga’s (Equus quagga quagga, Weddell 2002). Illnesses (e.g. rinderpest virus, sarcopses) are known to have dec- imated almost entire ungulate populations in a relatively short period of time; however, this scenario alone does not explain why other subspecies of ungulates in the Ibe- rian Peninsula were not affected by similar pathogens. Furthermore, species like the hirola and the coastal topi (Damaliscus lunatus topi), whose populations crashed by illnesses, were able to recover fully (Jowers et al. 2020). The new database of C. p. pyrenaica museum specimens allowed us to make some inferences about the extinction pattern followed by these animals across the last 200 years. Since male trophies are more attractive for collection than female ones, a male-bias is expected, with a progressive decrease of adult males and a disproportionately high num- bers of females and young individuals over time (Loveridge et al. 2006). Likewise, older animals are favoured since they correspond to bigger trophies (Geist 1966). We can cautiously state that this overall scenario is mirrored by our data (73.3% males vs 26.7% females). However, we did not observe an increase of female individuals in the biolog- ical collections over time, while we observed an increase in the average age of the collected individuals after the 1913 hunting ban. This age-related pattern could be explained by admitting that collectors were mostly interested in securing a male specimen for their museum and that after 1913 the majority of specimens were obtained from aged individuals post mortem by natural causes. The information retrieved to build our database also al- lowed us to flag one specimen of potentially high historical value from a private collection. Indeed, the trophy we gen- otyped in this study turned out to hold the same C. p. pyr- enaica CR haplotype to the attempted cloned individual in the year 2000 (Folch et al. 2009; Kupferschmidt 2014), which also originated from the last population of this taxon located on the Spanish side of the Pyrenees, and precisely from the female whose skin is exposed at the visitor centre of Ordesa y Monte Perdido National Park (Huesca, Spain). This finding conforms to the believed pattern of low lev- els of both neutral and adaptive genetic diversity of the Pyrenean ibex (Amills et al. 2004b) and confirms the ge- ographic origin of trophy. Significantly, the Pyrenean ibex was referred to as extremely rare since its first mention in an official document dating back to 1767. Low levels of genetic diversity are expected in ungulate populations zse.pensoft.net 218 experiencing strong bottlenecks, which may even turn into fixation of unique mtDNA haplotypes, as in the case of the Cypriot mouflon Ovis orientalis ophion (Guerrini et al. 2015), and inbreeding, thus triggering an extinction vortex. Experimental data from another related taxon, the bighorn sheep (Ovis canadensis) (Berger 1990) suggest that the small population size and the associated low genetic vari- ability might have well been the main extinction driver for the Pyrenean ibex. Concordantly, we found that the young male killed during a hunting trip organised by British vil- lagers based in Pau (France) beyond the Spanish border (since probably the French populations had gone already extinct or were already extremely scarce) held the same mtDNA found in a member of the same population which survived until two decades ago in the last stronghold of this animal: Ordesa y Monte Perdido National Park. This study further evidences the importance of biologi- cal collections as an invaluable repository of historical ma- terial in studies on extinct taxa. However, achieving spatial other than temporal representativeness of past populations is challenging. Private collections and small exhibits are an invaluable yet often overlooked repository of missing infor- mation (Casas-Marce et al. 2012) and this could well be the case with the Pyrenean ibex. In this context, citizen science may represent a useful tool to locate such specimens (Eu- ropean Commission 2013). Mitochondrial DNA data have recently pointed to C. p. pyrenaica as a distinct group from living heterospecifics other than the Alpine ibex (Urefia et al. 2018), thus making it a candidate Evolutionary Significant Unit (ESU: Manceau et al. 1999). We call for the need of genome-wide sequencing on already sampled and non-sam- pled specimens from private collections and small exhibits to clarify the still debated systematic relationships between the Iberian ibex and other European relatives, as well as with the potentially closely related C. caucasica (Urefia et al. 2018). Such studies would also allow to characterise the unique traits of C. p. pyrenaica underlying adaptation to extreme environments and a likely higher vulnerability to infectious diseases, as well as to quantify the amount of genetic diversity that became lost during its decline. Conclusions This case study shows the importance of identifying historical biological material in venues other than large and well-known public biological collections and cor- roborating the taxonomic identity and origin by means of preliminary genetic analyses. We call for the creation of an online public database of private collections hosting biological material to the benefit of biodiversity studies. Acknowledgements The Pyrenean ibex picture (https://commons.wikimedia. org/wiki/File:Celia_la_ %C3%BAltima_bucardo.JPG, CC- BY-SA-4.0) of Fig. 1a was taken by Jose Miguel Pintor Or- tego, while that of Fig. 3 (https://commons. wikimedia.org/ zse.pensoft.net FORCINA, G. et al.: Extinction of an Iberian ungulate wiki/File:Capra_pyrenaica_pyrenaica MHNT_ART_ 39. jpg, CC-BY-SA-4.0) by Roger Culos at the Muséum d’ His- toire Naturelle de la ville de Toulouse (MHNT). The map of the Pyrenees used in Fig. 2 is from Eric Gaba — Wikime- dia Commons user: Sting, and is free available at https:// commons. wikimedia.org/wiki/File:Pyrenees_relief_ map _ with_rivers-fr.svg under the terms of the GNU Free Doc- umentation License. This research was partially funded by the Plan Andaluz de Investigacion, Desarrollo e Inno- vacion (PAIDI) (RNM-118 group) and supported by fel- lowships from the Portuguese Foundation for Science and Technology to MJJ (FCT, SFRH/BPD/109148/2015) and GF (FCT, PTDC/BAA-AGR/28866/2017). The authors are thankful to the museum curators who assisted the col- lection of information on Pyerenan ibex specimens: Katrin Krohman (Naturmuseum und Forschungsinstitut Senck- enberg Frankfurt); Axel Schonhofer and Bettina Henrich (Naturhistorisches Museum Mainz/Landessammlung fur Naturkunde Mainz); Nigel T. Monaghan (National Muse- um of Ireland); Henri Cap (Muséum d’ Histoire Naturelle de Toulouse); Edi Stoeckli and Raffael Winkler (Naturhis- torisches Museum Basel); Frank E. Zachos and Barba- ra Herzig (Naturhistorisches Museum Wien); Martina Schenkel (Zoologisches Museum der Universitat Ztirich); Elisabeth Ludes-Fraulob (Musée Zoologique Strasbourg); Alexei Tikhonov and Gennady Baryshnikov (Zoologich- esku Institut Rossiisko1 Akademii Nauk); Laurent Lachaud and Laurent Charles (Museum d’ Histoire Naturelle de Bor- deaux); Marie-Laure Guerin, Eric Guiho and Luc Remy (Muséum d’ Histoire Naturelle de Nantes); Ricardo Garcia Gonzalez (Instituto Pirenaico de Ecologia, IPE-CSIC). References Abascal F, Corvelo A, Cruz EF, Villanueva-Cafias JL, Vlasova A, Marcet-Houben M, Martinez-Cruz B, Yu Cheng J, Prieto P, Quesada V, Quilez J, Li G, Garcia F, Rubio-Camarillo M, Frias L, Ribeca P, Ca- pella-Gutiérrez S, Rodriguez JM, Camara F, Lowy E, Cozzuto L, Erb I, Tress ML, Rodriguez-Ales JL, Ruiz-Orera J, Reverter F, Casas-Marce M, Soriano L, Arango JR, Derdak S, Galan B, Blanc J, Gut M, Lo- rente-Galdos B, Andrés-Nieto M, Lopez-Otin C, Valencia A, Gut I, José L, Garcia JL, Guigd R, Murphy WJ, Ruiz-Herrera A, Marques- Bonet T, Roma G, Notredame C, Mailund T, Alba MM, Gabaldon T, Alioto T, Godoy JA (2016) Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx. Ge- nome Biology 17: e251. https://do1.org/10.1186/s13059-016-1090-1 Acevedo P, Cassinello J (2009) Biology, ecology and status of Iberi- an ibex Capra pyrenaica: a critical review and research prospec- tus. Mammal Review 39(1): 17-32. https://doi.org/10.1111/j.1365- 2907.2008.00138.x Alpers DL, Van Vuuren BJ, Arctander P, Robinson TJ (2004) Population genetics of the roan antelope (Hippotragus equinus) with sugges- tions for conservation. Molecular Ecology 13: 1771-1784. https:// doi.org/10.1111/).1365-294X.2004.02204.x Amills M, Capote J, Tomas A, Kelly L, Obexer-Ruff G, Angiolillo A, Sanchez A (2004a) Strong phylogeographic relationships among three goat breeds from the Canary Islands. Journal of Dairy Re- search 71: 257-262. https://doi.org/10.1017/S0022029904000342 Zoosyst. Evol. 97 (1) 2021, 211-221 Amills M, Jiménez N, Jordana J, Riccardi A, Fernandez-Arias A, Guiral J, Bouzat JL, Folch J, Sanchez A (2004b) Low diversity in the major histocompatibility complex class II DRB/ gene of the Spanish ibex, Capra pyrenaica. Heredity 93: 266-272. https://doi.org/10.1038/ sj.hdy.6800499 Angelone-Alasaad S, Biebach I, Pérez JM, Soriguer RC, Granados JE (2017) Molecular analyses reveal unexpected genetic structure in Iberian ibex populations. PLoS ONE 12: e0170827. https://doi. org/10.1371/journal.pone.0170827 Astre G (1952) Quelques etapes de la disparition du bouquetin aux Pyrénées Centrales. Revue de Comminges 54: 129-146. Barnosky A, Matzke N, Tomiya S, Wogan GOU, Swartz B, Quental TB, Marshall C, McGuire JL, Lindsey EL, Maguire KC, Mersey B, Fer- rer EA (2011) Has the Earth’s sixth mass extinction already arrived? Nature 471: 51-57. https://doi.org/10.1038/nature09678 Cabrera A (1911) The subspecies of the Spanish ibex. Proceedings of the Zoological Society of London 66: 963-977. https://doi. org/10.1111/j.1096-3642.1911.tb01967.x Cabrera A (1914) Fauna ibérica: Mamiferos. Museo Nacional de Cien- cias Naturales, Madrid. Campos PF, Kristensen T, Orlando L, Sher A, Kholodova MV, Gother- str6m A, Hofreiter M, Druckerm DG, Kosintsev P, Tikhonov A, Baryshnikov GF, Willerslev E, Gilbert MT (2010) Ancient DNA sequences point to a large loss of mitochondrial genetic diversity in the saiga antelope (Saiga tatarica) since the Pleistocene. Mo- lecular Ecology 19: 4863-4875. https://doi.org/10.1111/j.1365- 294X.2010.04826.x Casas-Marce M, Revilla R, Fernandes M, Rodriguez A, Delibes M, Godoy JA (2012) The value of hidden scientific resources: preserved animal specimens from private collections and small museums. Bio- Science 62: 1077-1082. https://doi.org/10.1525/bi0.2012.62.12.9 Collen B, McRae L, Deinet S, de Palma A, Carranza T, Cooper N, Loh J, Baillie JEM (2011) Predicting how populations decline to extinc- tion. Philosophical Transactions of the Royal Society Series B 366: 2577-2586. https://doi.org/10.1098/rstb.2011.0015 Crampe J-P (1991) Le bouquetin ibérique. Eléments pour une réintro- duction au versant nord des Pyrénées occidentales. Documents sci- entifiques du Parc national des Pyrénées 26: 1-187. Crampe J-P, Crégut-Bonnoure E (1994) Le massif des Pyrénées, habi- tat naturel du bouquetin ibérique (Capra pyrenaica, Schinz, 1838). Evolution temporo-spatiale de |’espece de la préhistoire a nos jours. Ibex Speciale GSE 2: 39-48, Crampe J-P, Sourp E, Cavailhes J (2015) Réintroduction du Bouquetin ibérique dans les Pyrénées, quand le réve devient réalité. 33°™*s ren- contres du GEEFSM, Balme, Valli di Lanzo, Province de Torino, Italy, 21 to 24 May 2015. Dobrynin P, Liu S, Tamazian G, Xiong Z, Yurchenko AA, Krashenin- nikova K, Kliver S, Schmidt-Kiintzel A, Koepfli KP, Johnson W, Kuderna LFK, Garcia-Pérez R, Manuel M, Godinez R, Komissarov A, Makunin A, Brukhin V, Qiu W, Zhou L, Li F, Yi J, Driscoll C, Antunes A, Oleksyk TK, Eizirik E, Perelman P, Roelke M, Wildt D, Diekhans M, Marques-Bonet T, Marker L, Bhak J, Wang J, Zhang G, O’Brien SJ (2015) Genomic legacy of the African cheetah, Acinonyx jJubatus. Genome Biology 16: 1—20. https://doi.org/10.1186/s13059- 015-0837-4 Driscoll CA, Menotti-Raymond M, Nelson G, Goldstein D, O’Brien SJ (2002) Genomic microsatellites as evolutionary chronometers: a test in wild cats. Genome Research 12: 414-423. https://doi. org/10.1101/gr.185702 ZS Drummond A, Suchard MA, Xie D, Rambaut A (2012) Bayesian phylo- genetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution 29: 1969-1973. https://doi.org/10.1093/molbev/mss075 Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian co- alescent inference of past population dynamics from molecular se- quences. Molecular Biology and Evolution 22: 1185-1192. https:// doi.org/10.1093/molbev/msi103 European Commission (2013) Green paper on citizen science. Citizen science for Europe: towards a better society of empowered citizens and enhanced research. SOCIENTIZE Consortium, Zaragoza. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10: 564-567. https://doi. org/10.1111/j.1755-0998 .2010.02847.x Farrington HL, Lawson LP, Petren K (2019) Predicting population ex- tinctions in Darwin’s finches. Conservation Genetics 20: 825-836. https://doi.org/10.1007/s10592-019-01175-3 Folch J, Cocero MJ, Chesne P, Alabart JL, Dominguez V, Cognie Y, Roche A, Fernandez-Arias A, Marti JI, Sanchez P, Echegoyen E, Beckers JF, Bonastre AS, Vignon X (2009) First birth of an animal from an extinct subspecies (Capra pyrenaica pyrenaica) by cloning. Theriogenology 71: 1026-1034 https://do1.org/10.1016/j.theriog- enology.2008.11.005 Forcina G, Guerrini M, van Grouw H, Gupta BK, Panayides P, Hadjiger- ou P, Al-Sheikhly OF, Awan MN, Khan AA, Zeder MA, Barbanera F (2015) Impacts of biological globalization in the Mediterranean: un- veiling the deep history of human-mediated gamebird dispersal. Pro- ceedings of the National Academy of Sciences of the United States of America 112: 3296-3301. https://doi.org/10.1073/pnas. 1500677112 Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915—925. PMCID: PMC1208208. https://doi.org/10.1093/genetics/147.2.915 Garcia I, Napp S, Casal J, Perea A, Allepuz A, Alba A, Carbonero A, Arenas A (2009) Bluetongue epidemiology in wild ruminants from Southern Spain. European Journal of Wildlife Research 55: 173— 178. https://doi.org/10.1007/s10344-008-023 1-6 Garcia-Gonzalez R (1991) Inventario de la poblacion espafiola de bu- cardo. Informe inédito. CSIC-ICONA. Garcia-Gonzalez R, Alados CL, Ameztoy JM, Escos J, Herrero J, Hidal- go R, Pascual R (1991) Inventario de la poblacion espafiola de bu- cardo. Informe inédito. CSIC-ICONA, Jaca-Madrid. http://digital. csic.es/bitstream/10261/36824/1/R3_Inf_final_bucardo_1991.pdf Garcia-Gonzalez R, Escoés J, Alados CL (1996) Una poblacion en peligro: el bucardo. In: Alados CL, Escos J (Eds) Ecologia y comportamiento de la cabra montés. Consideraciones para su gestion. Monografias del Museo de Ciencias Naturales, CSIC, Madrid, 105-120. Garcia-Gonzalez R, Herrero J (1999) El Bucardo de los Pirineos: histo- ria de una extincion. Galemys 11: 17—26. Garcia-Gonzalez R, Hidalgo R, Ameztoy JM, Herrero J (1992) Census, population structure and habitat use of a chamois population in Or- desa N.P. living in sympatry with the Pyrenean wild goat. In: Spitz F, Janeau G, Gonzalez G, Aulagnier S (Eds) Ongulés/Ungulates 91. SFEPM-IRGM, Paris-Toulouse, 321-325. Geist V (1966) The evolution of horn-like organs. Behaviour 27: 175— 214. https://doi.org/10.1163/156853966X00155 Godinho R, Abaigar T, Lopes S, Essalhi A, Ouragh L, Cano M, Ferrand N (2012) Conservation genetics of the endangered dorcas gazelle (Gazella dorcas spp.) in northwestern Africa. Conservervation Ge- netics 13: 1003-1015. https://doi.org/10.1007/s10592-012-0348-8 zse.pensoft.net 220 Gomez-Guillamon F, Caballero-Gomez J, Agtiero M, Camacho-Sillero L, Risalde MA, Zorrilla I, Villalba R, Rivero-Juarez A, Garcia-Bo- canegra I (2020) Re-emergence of bluetongue virus serotype 4 in Iberian ibex (Capra pyrenaica) and sympatric livestock in Spain, 2018-2019. Transbounding Emerging Diseases 00: 1-9. https://doi. org/10.1111/tbed. 13696 Gourdon M (1908) Note sur une Série de Cranes de Mamiferes des Pyrénées. Bulletin de la Société des sciences naturelles de |’Ouest de la France 8: 1—34. Gouy M, Guindon S, Gascuel O (2010) SeaView version 4. A multiplat- form graphical user interface for sequence alignment and phyloge- netic tree building. Molecular Biology and Evolution 27: 221-224. https://doi.org/10.1093/molbev/msp259 Granados JE, Pérez JM, Soriguer RC, Fandos P, Ruiz-Martinez I (1997) On the biometry of the Spanish ibex, Capra pyrenaica, from Sierra Nevada (Southern Spain). Folia Zoologica 46: 9-14. Granados JE, Pérez JM, Marquez FJ, Serrano E, Soriguer RC, Fandos P (2001) La cabra montés (Capra pyrenaica, Schinz 1838). Galemys 13: 3-38. Grossen C, Biebach I, Angelone-Alasaad S, Keller LF, Croll D (2017) Population genomics analyses of European ibex species show lower diversity and higher inbreeding in reintroduced populations. Evolu- tionary Applications 11: 123-139. https://doi.org/10.1111/eva. 12490 Guerrini M, Forcina G, Panayides P, Lorenzini R, Garel M, Anayiotos P, Kassinis N, Barbanera F (2015) Molecular DNA identity of the mouflon of Cyprus (Ovis orientalis ophion, Bovidae): Near Eastern origin and divergence from Western Mediterranean conspecific pop- ulations. Systematics and Biodiversity 13: 472-483. https://doi.org/ 10.1080/14772000.2015.1046409 Gur S (2008) A serologic investigation of blue tongue virus (BTV) in cattle, sheep and Gazella subgutturosa subgutturosa in southeastern turkey. Tropical Animal Health and Production 40: 217-221. https:// doi.org/10.1007/s11250-007-9083-4 Hammond RL, Macasero W, Flores B, Osama BM, Wacher T, Bruford MW (2001) Phylogenetic Reanalysis of the Saudi Gazelle and Its Implications for Conservation. Conservation Biology 15: 1123- 1133. https://doi.org/10.1046/j.1523-1739.2001.0150041123.x Harpending HC (1994) Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Human Biology 66: 591-600. PMID: 8088750. Herrero J, Acevedo P, Arnal MC, Fernandez de Luco D, Fonseca C, Garcia-Gonzalez R, Pérez JM, Sourp E (2020) Capra pyrenaica. Fe- cha de acceso: | Oct 2020. https://doi.org/10.2305/IUCN.UK.2020- 2.RLTS.T3798A170192604.en Hidalgo R, Garcia-Gonzalez R (1995) Remnant Pyrenean wild goat population in Ordesa and Monte Perdido National Park, Pyrénées (Spain). Caprinae News IUCN 8/9: 9-13. Hmwe SS, Zachos FE, Eckert I, Lorenzini R, Fico R, Hartl GB (2006) Conservation genetics of the endangered red deer from Sardinia and Mesola with further remarks on the phylogeography of Cervus elaphus corsicanus. Biological Journal of the Linnean Society 88: 691-701. https://doi.org/10.1111/j.1095-8312.2006.00653.x Le retour du bouquetin (2020) Le retour du bouquetin dans les Pyrénées. http://www. bouquetin-pyrenees.fr/ [Accessed on 06/12/2020] Iyengar A, Gilbert T, Woodfine T, Knowles JM, Diniz FM, Brenneman RA, Louis Jr EE, Maclean N (2007) Remnants of ancient genetic diversity preserved within captive groups of scimitar-horned oryx (Oryx dammah). Molecular Ecology 16: 2436-2449. https://doi. org/10.1111/j.1365-294X.2007.03291.x zse.pensoft.net FORCINA, G. et al.: Extinction of an Iberian ungulate IUCN (1996) Cervus schomburgki, World Conservation Monitoring Centre 1996, IUCN Red List of Threatened Species. Jiménez J (2016) Declive y auge de las cabras montesas: una historia diferente. Quercus 368: 14—24. Jiménez N, Folch J, Fernandez-Arias A, Guiral J, Sanchez A (1999) Estudio genético mediante marcadores microsatélites de las pobla- ciones de cabra montés. ITEA, Vol. Extra 20: 300-302. Jowers MJ, Queirés J, Resende Pinto R, Ali AH, Mutinda M, Ange- lone S, Alves PC, Godinho R (2020) Genetic diversity in natural range remnants of the critically endangered hirola antelope. Zoo- logical Journal of the Linnean Society 190: 384-395. https://doi. org/10.1093/zoolinnean/zlz174 Kidjo N, Feracci G, Bideau E, Gonzalez G, Mattéi C, Marchand B, Aulagnier S (2007) Extirpation and reintroduction of the Corsican red deer Cervus elaphus corsicanus in Corsica. Oryx 41: 488-494. https://doi.org/10.1017/S0030605307012069 Kingdon J, Happold D, Butynski T, Hoffmann M, Happold M, Kalina MJ (2013) Mammals of Africa. Bloomsbury Publishing, London, 720 pp. Klein RG (1974) On the tazonomic status, distribution and ecology of the blue antelope, Hippotragus leucophaeus. Annals of the South African Museum 65: 99-143. Kupferschmidt K (2014) Can cloning revive Spain’s extinct mountain goat? Science 344: 137-138. https://doi.org/10.1126/science.344.6180.137 Labareére J (1985) Une espece en voie de disparition: le bouquetin des Pyrénées. Pyrénées 142: 115-131. Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947-2948. https://doi.org/10.1093/bioinformatics/btm404 Leakey R, Lewin R (1996) The Sixth Extinction: Patterns of Life and the Future of Humankind. Doubleday, London, 271 pp. Leslie SO (1894) Sir Victor Brooke, Sportsman & Naturalist. A Mem- oir of His Life and Extracts from his Letters and Journals. Murray, London, 294 pp. Little HA, Gilbert TC, Athorn ML, Marshall AR (2016) Evaluating Conservation Breeding Success for an Extinct-in-the-Wild Ante- lope. PLoS ONE 11: e0166912. https://doi.org/10.1371/journal. pone.0166912 Lorca-Oré C, Lépez-Olvera JR, Ruiz-Fons F, Acevedo P, Garcia-Bo- canegra I, Oleaga A, Gortazar C, Pujols J (2014) Long-Term Dy- namics of Bluetongue Virus in Wild Ruminants: Relationship with Outbreaks in Livestock in Spain, 2006—2011. PloS ONE 9: e100027. https://doi.org/10.1371/journal.pone.0100027 Loveridge AJ, Reynolds JC, Milner-Gulland EJ (2006) Does sport hunting benefit conservation? In: Macdonald DW, Service K (Eds) Key Topics in Conservation Biology. Blackwell Publishing, Oxford, 224—240. Manceau V, Crampe JP, Boursot P, Taberlet P (1999) Identification of evolutionary significant units in the Spanish wild goat Capra pyre- naica (Mammalia, Artiodactyla). Animal Conservation 2: 33-39. https://doi.org/10.1017/S1367943099000335 Miller MA, Pfeiffer W, Schwartz T (2010) Creating the CIPRES Sci- ence Gateway for inference of large phylogenetic trees. In: Gate- way Computing Environments Workshop (GCE), 1-8. https://doi. org/10.1109/GCE.2010.5676129 Moco G, Guerreiro M, Ferreira AF, Rebelo A, Loureiro A, Petrucci-Fon- seca F, Pérez JM (2006) The ibex Capra pyrenaica returns to its for- mer Portuguese range. Oryx 40: 351-354. https://doi.org/10.1017/ S0030605306000718 Zoosyst. Evol. 97 (1) 2021, 211-221 Moco G (2015) Ecology and Conservation of the Iberian ibex (Capra pyre- naica Schinz, 1838) in the Peneda-Gerés National Park, Portugal. Ph. D. Thesis, Jaén University. https://do1.org/10.1515/mammalia-2013-0139 O’Brien SJ, Roelke ME, Marker L, Newman A, Winkler CA, Meltzer D, Colly L, Evermann JF, Bush M, Wildt DE (1985) Genetic basis for species vulnerability in the cheetah. Science 227: 1428-1434. https://doi.org/10.1126/science.2983425 Pérez JM, Granados JE, Soriguer RC, Fandos P, Marquez FJ, Crampe JP (2002) Distribution, status and conservation problems of the Span- ish ibex, Capra pyrenaica (Mammalia: Artiodactyla). Mammal Re- view 32: 26-39. https://doi.org/10.1046/}.1365-2907.2002.00097.x Pérez T, Gonzalez I, Essler SE, Fernandez M, Dominguez A (2014) The shared mitochondrial genome of Rupicapra pyrenaica ornata and Rupicapra rupicapra cartusiana: old remains of a common past. Molecular Phylogenetics and Evolution 79: 375-379. https://doi. org/10.1016/j.ympev.2014.07.004 Phébus G (1387-1389) Le Livre de chasse. Manuscript 616 (beginning of 14" century). Biblioth¢que Nationale de France, Département des manuscrits, Paris, France: folios 21 and 86. http://expositions. bnf. fr/ phebus/grands/c04_616.htm [and] http://expositions.bnf.fr/phebus/ grands/c48_616.htm Pievani T (2014) The sixth mass extinction: Anthropocene and the hu- man impact on biodiversity. Rendiconti Lincei. Scienze Fisiche e Naturali 25: 85—93. https://doi.org/10.1007/s12210-013-0258-9 Posada D (2008) JModelTest: Phylogenetic model averaging. Molecu- lar Biology and Evolution 25: 1253-1256. https://doi.org/10.1093/ molbev/msn083 Prewer E, Kutz S, Leclerc LM, Kyle CJ (2020) Already at the bottom? Demographic declines are unlikely further to undermine genetic di- versity of a large Arctic ungulate: muskox, Ovibos moschatus (Ar- tiodactyla: Bovidae). Biological Journal of the Linnean Society 129: 459-469. https://doi.org/10.1093/biolinnean/blz175 Rogers AR (1995) Genetic evidence for a Pleistocene pop- ulation expansion. Evolution 49: 608-615. _https://doi. org/10.1111/j.1558-5646.1995 tb02297.x Rokosz M (1995) History of the Aurochs (Bos Taurus Primigenius) in Poland. Animal Genetic Resources Information 16: 5—12. https:// doi.org/10.1017/S1014233900004582 Rossi L, Tizzani P, Rambozzi L, Moroni B, Meneguz PG (2019) Sanitary Emergencies at the Wild/Domestic Caprines Interface in Europe. Animals (Basel) 9(11): e922. https://do1.org/10.3390/an19 110922 Rozas J (2009) DNA Sequence Polymorphism Analysis Using DnaSP. In: Posada D (Ed.) Bioinformatics for DNA Sequence Analysis: Methods in Molecular Biology Series (Vol. 537). Humana Press, New York, 337-350. https://doi.org/10.1007/978-1-59745-251-9_17 Schneider S, Excoffier L (1999) Estimation of past demographic param- eters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics 152: 1079-1089. PMCID: PMC 1460660. Soulé ME (1987) Viable Populations for Conservation. Cambridge University Press, Cambridge, 189 pp. https://doi.org/10.1017/ CBO9780511623400 Tajima F (1989) Statistical method for testing the neutral mutation hy- pothesis by DNA polymorphism. Genetics 123: 585-595. PMCID: PMC1203831. https://doi.org/10.1093/genetics/123.3.585 Thorne ET, Williams ES, Sprakerm TR, Helms W, Segerstrom T (1988) Bluetongue in free-ranging pronghorn antelope (Antilocapra americana) in Wyoming: 1976 and 1984. Journal of Wildlife Dis- eases 24: 113-119. https://doi.org/10.7589/0090-3558-24.1.113 221 Urefia I, Arsuaga JL, Galindo-Pellicena MA, Gotherstrom A, Valdios- era C (2011) Filogenia y evolucion local de la cabra montés (Capra pyrenaica) en el yacimiento Cuaternario de Chaves (Huesca, Es- pafia). Boletin de la Real Sociedad Espafiola de Historia Natural. Seccion Geoldgica 105: 5-14. UrefiaI, Ersmark E, Samaniego JA, Galindo-Pellicena MA, Crégut-Bon- noure E, Bolivar H, Gomez-Olivencia A, Rios-Garaizar J, Garate D, Dalén L, Arsuaga JL, Valdiosera CE (2018) Unravelling the genetic history of the European wild goats. Quaternary Science Review 185: 189-198. https://doi.org/10.1016/j.quascirev.2018.01.017 van Vuure C (2005) Retracing the Aurochs — History, Morphology and Ecology of an Extinct Wild Ox. Pensoft Publishers, Sofia-Moscow, 431 pp. Villalta M, Folch J, Alabart JL (1997) Estudio genético molecular de las poblaciones de cabra montés de la Peninsula ibérica. Programa LIFE — Plan de Recuperacion del Bucardo. Informe Final. SIA- DGA, Zaragoza. Weddell BJ (2002) Conserving Living Natural Resources in the Con- text of a Changing World. Cambridge University Press, Cambridge, 442 pp. https://doi.org/10.1017/CBO97805 11804298 Woutersen K (2012) El Bucardo de los Pirineos. Kees Woutersen Publi- caciones, Huesca, 208 pp. Woutersen K (2019) El camino de los Bucardos. Scribo Editorial, Hu- esca, 128 pp. Zeng Y, Jiang Z, Li C (2007) Genetic variability in relocated Pere Da- vid’s deer (Elaphurus davidianus) populations — Implications to re- introduction program. Conservation Genetics 8: 1051—1059. https:// doi.org/10.1007/s10592-006-9256-0 Supplementary material | Tables S1, S2 Authors: Giovanni Forcina, Kees Woutersen, Santiago Sanchez-Ramirez, Samer Angelone, Jean P. Crampe, Jesus M. Pérez, Paulino Fandos, José Enrique Grana- dos, Michael J. Jowers Data type: Table Explanation note: Table S1. Details of C. p. pyrenaica samples used in this study. FR.: France; SP.: Spain. Ta- ble S2. Database of the 45 C. p. pyrenaica specimens identified in this study as a reference for further re- search. IPE: Pyrenean Institute of Ecology; PNOMP: Ordesa and Monte Perdido National Park. Copyright notice: This dataset is made available under the Open Database License (http://opendatacom- mons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Data- set while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/zse.97.61854.suppl1 zse.pensoft.net