Extended and Continuous Decline in Effective Population Size Results in Low Genomic Diversity in the World's Rarest Hyena Species, the Brown Hyena.

Hyenas (family Hyaenidae), as the sister group to cats (family Felidae), represent a deeply diverging branch within the cat-like carnivores (Feliformia). With an estimated population size of <10,000 individuals worldwide, the brown hyena (Parahyaena brunnea) represents the rarest of the four extant hyena species and has been listed as Near Threatened by the IUCN. Here, we report a high-coverage genome from a captive bred brown hyena and both mitochondrial and low-coverage nuclear genomes of 14 wild-caught brown hyena individuals from across southern Africa. We find that brown hyena harbor extremely low genetic diversity on both the mitochondrial and nuclear level, most likely resulting from a continuous and ongoing decline in effective population size that started ∼1 Ma and dramatically accelerated towards the end of the Pleistocene. Despite the strikingly low genetic diversity, we find no evidence of inbreeding within the captive bred individual and reveal phylogeographic structure, suggesting the existence of several potential subpopulations within the species.

Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations.

BACKGROUND: Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts? RESULTS: In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these. CONCLUSIONS: The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.

Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild.

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe's largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

Multilocus phylogeny, species age and biogeography of the Lesser Antillean anoles.

Lesser Antillean anoles provide classic examples of island radiations. A detailed knowledge of their phylogeny and biogeography, in particular how the age of species relate to the ages of their respective islands and the age of their radiation, is essential to elucidate the tempo and mechanisms of these radiations. We conduct a large-scale phylogenetic and phylogeographic investigation of the Lesser Antillean anoles using multiple genetic markers and comprehensive geographic sampling of most species. The multilocus phylogeny gives the first well-supported reconstruction of the interspecific relationships, and the densely sampled phylogeography reveals a highly dynamic system, driven by overseas dispersal, with several alternative post-dispersal colonisation trajectories. These radiations currently occupy both the outer-older (Eocene to Miocene), and the inner-younger (<8mybp), Lesser Antillean arcs. The origin of these radiations corresponds with the age of the ancient outer arc. However, the ages of extant species (compatible with the age of other small terrestrial amniotes) are much younger, about the age of the emergence of the younger arc, or less. The difference between the age of the radiation and the age of the extant species suggests substantial species turnover on older arc islands, most likely through competitive replacement. Although extant anoles are extremely speciose, this may represent only a fraction of their biodiversity over time. While paraphyly enables us to infer several recent colonization events, the absence of the younger arc islands and extant species at the earlier and middle stages of the radiation, does not allow the earlier inter-island colonization to be reliably inferred. Reproductive isolation in allopatry takes a very considerable time (in excess of 8my) and sympatry appears to occur only late in the radiation. The resolved multilocus phylogeny, and relative species age, raise difficulties for some earlier hypotheses regarding size evolution, and provide no evidence for within-island speciation.

Tropical ancient DNA reveals relationships of the extinct Bahamian giant tortoise Chelonoidis alburyorum.

Ancient DNA of extinct species from the Pleistocene and Holocene has provided valuable evolutionary insights. However, these are largely restricted to mammals and high latitudes because DNA preservation in warm climates is typically poor. In the tropics and subtropics, non-avian reptiles constitute a significant part of the fauna and little is known about the genetics of the many extinct reptiles from tropical islands. We have reconstructed the near-complete mitochondrial genome of an extinct giant tortoise from the Bahamas (Chelonoidis alburyorum) using an approximately 1 000-year-old humerus from a water-filled sinkhole (blue hole) on Great Abaco Island. Phylogenetic and molecular clock analyses place this extinct species as closely related to Galápagos (C. niger complex) and Chaco tortoises (C. chilensis), and provide evidence for repeated overseas dispersal in this tortoise group. The ancestors of extant Chelonoidis species arrived in South America from Africa only after the opening of the Atlantic Ocean and dispersed from there to the Caribbean and the Galápagos Islands. Our results also suggest that the anoxic, thermally buffered environment of blue holes may enhance DNA preservation, and thus are opening a window for better understanding evolution and population history of extinct tropical species, which would likely still exist without human impact.

Multi-locus phylogeny and species delimitation of Australo-Papuan blacksnakes (Pseudechis Wagler, 1830: Elapidae: Serpentes).

Genetic analyses of Australasian organisms have resulted in the identification of extensive cryptic diversity across the continent. The venomous elapid snakes are among the best-studied organismal groups in this region, but many knowledge gaps persist: for instance, despite their iconic status, the species-level diversity among Australo-Papuan blacksnakes (Pseudechis) has remained poorly understood due to the existence of a group of cryptic species within the P. australis species complex, collectively termed "pygmy mulga snakes". Using two mitochondrial and three nuclear loci we assess species boundaries within the genus using Bayesian species delimitation methods and reconstruct their phylogenetic history using multispecies coalescent approaches. Our analyses support the recognition of 10 species, including all of the currently described pygmy mulga snakes and one undescribed species from the Northern Territory of Australia. Phylogenetic relationships within the genus are broadly consistent with previous work, with the recognition of three major groups, the viviparous red-bellied black snake P. porphyriacus forming the sister species to two clades consisting of ovoviviparous species.

The future of ancient DNA: Technical advances and conceptual shifts.

Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7× coverage (mammoth) in 2008 to more than 50× coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.

Ancient DNA reveals differences in behaviour and sociality between brown bears and extinct cave bears.

Ancient DNA studies have revolutionized the study of extinct species and populations, providing insights on phylogeny, phylogeography, admixture and demographic history. However, inferences on behaviour and sociality have been far less frequent. Here, we investigate the complete mitochondrial genomes of extinct Late Pleistocene cave bears and middle Holocene brown bears that each inhabited multiple geographically proximate caves in northern Spain. In cave bears, we find that, although most caves were occupied simultaneously, each cave almost exclusively contains a unique lineage of closely related haplotypes. This remarkable pattern suggests extreme fidelity to their birth site in cave bears, best described as homing behaviour, and that cave bears formed stable maternal social groups at least for hibernation. In contrast, brown bears do not show any strong association of mitochondrial lineage and cave, suggesting that these two closely related species differed in aspects of their behaviour and sociality. This difference is likely to have contributed to cave bear extinction, which occurred at a time in which competition for caves between bears and humans was likely intense and the ability to rapidly colonize new hibernation sites would have been crucial for the survival of a species so dependent on caves for hibernation as cave bears. Our study demonstrates the potential of ancient DNA to uncover patterns of behaviour and sociality in ancient species and populations, even those that went extinct many tens of thousands of years ago.

Widespread parallel population adaptation to climate variation across a radiation: implications for adaptation to climate change.

Global warming will impact species in a number of ways, and it is important to know the extent to which natural populations can adapt to anthropogenic climate change by natural selection. Parallel microevolution within separate species can demonstrate natural selection, but several studies of homoplasy have not yet revealed examples of widespread parallel evolution in a generic radiation. Taking into account primary phylogeographic divisions, we investigate numerous quantitative traits (size, shape, scalation, colour pattern and hue) in anole radiations from the mountainous Lesser Antillean islands. Adaptation to climatic differences can lead to very pronounced differences between spatially close populations with all studied traits showing some evidence of parallel evolution. Traits from shape, scalation, pattern and hue (particularly the latter) show widespread evolutionary parallels within these species in response to altitudinal climate variation greater than extreme anthropogenic climate change predicted for 2080. This gives strong evidence of the ability to adapt to climate variation by natural selection throughout this radiation. As anoles can evolve very rapidly, it suggests anthropogenic climate change is likely to be less of a conservation threat than other factors, such as habitat loss and invasive species, in this, Lesser Antillean, biodiversity hot spot.

Recent increase in species-wide diversity after interspecies introgression in the highly endangered Iberian lynx.

Genetic diversity is lost in small and isolated populations, affecting many globally declining species. Interspecific admixture events can increase genetic variation in the recipient species' gene pool, but empirical examples of species-wide restoration of genetic diversity by admixture are lacking. Here we present multi-fold coverage genomic data from three ancient Iberian lynx (Lynx pardinus) approximately 2,000-4,000 years old and show a continuous or recurrent process of interspecies admixture with the Eurasian lynx (Lynx lynx) that increased modern Iberian lynx genetic diversity above that occurring millennia ago despite its recent demographic decline. Our results add to the accumulating evidence for natural admixture and introgression among closely related species and show that this can result in an increase of species-wide genetic diversity in highly genetically eroded species. The strict avoidance of interspecific sources in current genetic restoration measures needs to be carefully reconsidered, particularly in cases where no conspecific source population exists.

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa.

Evidence from numerous Pan-African savannah mammals indicates that open-habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open-habitat formations throughout sub-Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.

First mitogenome phylogeny of the sun bear Helarctos malayanus reveals a deep split between Indochinese and Sundaic lineages.

The sun bear Helarctos malayanus is one of the most endangered ursids, and to date classification of sun bear populations has been based almost exclusively on geographic distribution and morphology. The very few molecular studies focussing on this species were limited in geographic scope. Using archival and non-invasively collected sample material, we have added a substantial number of complete or near-complete mitochondrial genome sequences from sun bears of several range countries of the species' distribution. We here report 32 new mitogenome sequences representing sun bears from Cambodia, Thailand, Peninsular Malaysia, Sumatra, and Borneo. Reconstruction of phylogenetic relationships revealed two matrilines that diverged ~295 thousand years ago: one restricted to portions of mainland Indochina (China, Cambodia, Thailand; "Mainland clade"), and one comprising bears from Borneo, Sumatra, Peninsular Malaysia but also Thailand ("Sunda clade"). Generally recent coalescence times in the mitochondrial phylogeny suggest that recent or historical demographic processes have resulted in a loss of mtDNA variation. Additionally, analysis of our data in conjunction with shorter mtDNA sequences revealed that the Bornean sun bear, classified as a distinct subspecies (H. m. euryspilus), does not harbor a distinctive matriline. Further molecular studies of H. malayanus are needed, which should ideally include data from nuclear loci.

Museum DNA reveals a new, potentially extinct species of rinkhals (Serpentes: Elapidae: Hemachatus) from the Eastern Highlands of Zimbabwe.

Genetic information plays a pivotal role in species recognition and delimitation, but rare or extinct animals can be difficult to obtain genetic samples from. While natural history wet collections have proven invaluable in the description of novel species, the use of these historical samples in genetic studies has been greatly impeded by DNA degradation, especially because of formalin-fixation prior to preservation. Here, we use recently developed museum genomics approaches to determine the status of an isolated population of the elapid snake genus Hemachatus from Zimbabwe. We used multiple digestion phases followed by single strand sequencing library construction and hybridisation capture to obtain 12S and 16S rDNA sequences from a poorly preserved tissue sample of this population. Phylogenetic and morphological analyses in an integrated taxonomic framework demonstrate that the Zimbabwean rinkhals population represents an old and highly distinct lineage, which we describe as a new species, Hemachatus nyangensis sp. nov. Our phylogenetic dating analysis is compatible with venom spitting having evolved in response to the threat posed by early hominins, although more data are required for a robust test of this hypothesis. This description demonstrates the power of museum genomics in revealing rare or even extinct species: Hemachatus from Zimbabwe are only known from a small area of the Eastern Highlands known for high endemism. No living specimens have been seen since the 1980s, most likely due to dramatic land-use changes in the Eastern Highlands, suggesting that the species could be extinct. In view of its recognition as a highly distinct lineage, urgent action is required to determine whether any populations survive, and to safeguard remaining habitat.

Correction: Insight into the introduction of domestic cattle and the process of Neolithization to the Spanish region Galicia by genetic evidence.

[This corrects the article DOI: 10.1371/journal.pone.0249537.].

Bacterial Adaptation to Venom in Snakes and Arachnida.

Animal venoms are considered sterile sources of antimicrobial compounds with strong membrane-disrupting activity against multidrug-resistant bacteria. However, venomous bite wound infections are common in developing nations. Investigating the envenomation organ and venom microbiota of five snake and two spider species, we observed venom community structures that depend on the host venomous animal species and evidenced recovery of viable microorganisms from black-necked spitting cobra (Naja nigricollis) and Indian ornamental tarantula (Poecilotheria regalis) venoms. Among the bacterial isolates recovered from N. nigricollis, we identified two venom-resistant, novel sequence types of Enterococcus faecalis whose genomes feature 16 virulence genes, indicating infectious potential, and 45 additional genes, nearly half of which improve bacterial membrane integrity. Our findings challenge the dogma of venom sterility and indicate an increased primary infection risk in the clinical management of venomous animal bite wounds. IMPORTANCE Notwithstanding their 3 to 5% mortality, the 2.7 million envenomation-related injuries occurring annually-predominantly across Africa, Asia, and Latin America-are also major causes of morbidity. Venom toxin-damaged tissue will develop infections in some 75% of envenomation victims, with E. faecalis being a common culprit of disease; however, such infections are generally considered to be independent of envenomation. Here, we provide evidence on venom microbiota across snakes and arachnida and report on the convergent evolution mechanisms that can facilitate adaptation to black-necked cobra venom in two independent E. faecalis strains, easily misidentified by biochemical diagnostics. Therefore, since inoculation with viable and virulence gene-harboring bacteria can occur during envenomation, acute infection risk management following envenomation is warranted, particularly for immunocompromised and malnourished victims in resource-limited settings. These results shed light on how bacteria evolve for survival in one of the most extreme environments on Earth and how venomous bites must be also treated for infections.

Ancient DNA from the Asiatic Wild Dog (Cuon alpinus) from Europe.

The Asiatic wild dog (Cuon alpinus), restricted today largely to South and Southeast Asia, was widespread throughout Eurasia and even reached North America during the Pleistocene. Like many other species, it suffered from a huge range loss towards the end of the Pleistocene and went extinct in most of its former distribution. The fossil record of the dhole is scattered and the identification of fossils can be complicated by an overlap in size and a high morphological similarity between dholes and other canid species. We generated almost complete mitochondrial genomes for six putative dhole fossils from Europe. By using three lines of evidence, i.e., the number of reads mapping to various canid mitochondrial genomes, the evaluation and quantification of the mapping evenness along the reference genomes and phylogenetic analysis, we were able to identify two out of six samples as dhole, whereas four samples represent wolf fossils. This highlights the contribution genetic data can make when trying to identify the species affiliation of fossil specimens. The ancient dhole sequences are highly divergent when compared to modern dhole sequences, but the scarcity of dhole data for comparison impedes a more extensive analysis.

Ancient Mitogenomes Provide New Insights into the Origin and Early Introduction of Chinese Domestic Donkeys.

Both molecular data and archaeological evidence strongly support an African origin for the domestic donkey. Recent genetic studies further suggest that there were two distinct maternal lineages involved in its initial domestication. However, the exact introduction time and the dispersal process of domestic donkeys into ancient China are still unresolved. To address these questions, we retrieved three near-complete mitochondrial genomes from donkey specimens excavated from Gaoling County, Shaanxi Province, and Linxia Basin, Gansu Province, China, dated at 2,349-2,301, 469-311, and 2,160-2,004 cal. BP, respectively. Maximum-likelihood and Bayesian phylogenetic analyses reveal that the two older samples fall into the two different main lineages (i.e., clade Ⅰ and clade Ⅱ) of the domestic donkey, suggesting that the two donkey maternal lineages had been introduced into Midwestern China at least at the opening of Silk Road (approximately the first century BC). Bayesian analysis shows that the split of the two donkey maternal lineages is dated at 0.323 Ma (95% CI: 0.583-0.191 Ma) using root-tip dating calibrations based on near-complete mitogenomes, supporting the hypothesis that modern domestic donkeys go back to at least two independent domestication events. Moreover, Bayesian skyline plot analyses indicate an apparent female population increase between 5,000 and 2,500 years ago for clade I followed by a stable population size to the present day. In contrast, clade II keeps a relatively stable population size over the past 5,000 years. Overall, our study provides new insights into the early domestication history of Chinese domestic donkeys.

Middle Pleistocene genome calibrates a revised evolutionary history of extinct cave bears.

Palaeogenomes provide the potential to study evolutionary processes in real time, but this potential is limited by our ability to recover genetic data over extended timescales.1 As a consequence, most studies so far have focused on samples of Late Pleistocene or Holocene age, which covers only a small part of the history of many clades and species. Here, we report the recovery of a low coverage palaeogenome from the petrous bone of a ∼360,000 year old cave bear from Kudaro 1 cave in the Caucasus Mountains. Analysis of this genome alongside those of several Late Pleistocene cave bears reveals widespread mito-nuclear discordance in this group. Using the time interval between Middle and Late Pleistocene cave bear genomes, we directly estimate ursid nuclear and mitochondrial substitution rates to calibrate their respective phylogenies. This reveals post-divergence mitochondrial transfer as the dominant factor explaining their mito-nuclear discordance. Interestingly, these transfer events were not accompanied by large-scale nuclear introgression. However, we do detect additional instances of nuclear admixture among other cave bear lineages, and between cave bears and brown bears, which are not associated with mitochondrial exchange. Genomic data obtained from the Middle Pleistocene cave bear petrous bone has thus facilitated a revised evolutionary history of this extinct megafaunal group. Moreover, it suggests that petrous bones may provide a means of extending both the magnitude and time depth of palaeogenome retrieval over substantial portions of the evolutionary histories of many mammalian clades.

Estimating the dwarfing rate of an extinct Sicilian elephant.

Evolution on islands, together with the often extreme phenotypic changes associated with it, has attracted much interest from evolutionary biologists. However, measuring the rate of change of phenotypic traits of extinct animals can be challenging, in part due to the incompleteness of the fossil record. Here, we use combined molecular and fossil evidence to define the minimum and maximum rate of dwarfing in an extinct Mediterranean dwarf elephant from Puntali Cave (Sicily).1 Despite the challenges associated with recovering ancient DNA from warm climates,2 we successfully retrieved a mitogenome from a sample with an estimated age between 175,500 and 50,000 years. Our results suggest that this specific Sicilian elephant lineage evolved from one of the largest terrestrial mammals that ever lived3 to an island species weighing less than 20% of its original mass with an estimated mass reduction between 0.74 and 200.95 kg and height reduction between 0.15 and 41.49 mm per generation. We show that combining ancient DNA with paleontological and geochronological evidence can constrain the timing of phenotypic changes with greater accuracy than could be achieved using any source of evidence in isolation.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx.

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species' evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies' current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.