TY - JOUR A1 - Zancolli, Giulia A1 - Baker, Timothy G. A1 - Barlow, Axel A1 - Bradley, Rebecca K. A1 - Calvete, Juan J. A1 - Carter, Kimberley C. A1 - de Jager, Kaylah A1 - Owens, John Benjamin A1 - Price, Jenny Forrester A1 - Sanz, Libia A1 - Scholes-Higham, Amy A1 - Shier, Liam A1 - Wood, Liam A1 - Wüster, Catharine E. A1 - Wüster, Wolfgang T1 - Is Hybridization a Source of Adaptive Venom Variation in Rattlesnakes? A Test, Using a Crotalus scutulatus x viridis Hybrid Zone in Southwestern New Mexico JF - Toxins N2 - Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter-and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A2 such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species. KW - adaptation KW - Crotalus KW - evolution KW - hybridization KW - introgression KW - Mojave toxin KW - molecular evolution KW - venom Y1 - 2016 U6 - https://doi.org/10.3390/toxins8060188 SN - 2072-6651 VL - 8 PB - MDPI CY - Basel ER - TY - JOUR A1 - Yuan, Junxia A1 - Sheng, Guilian A1 - Preick, Michaela A1 - Sun, Boyang A1 - Hou, Xindong A1 - Chen, Shungang A1 - Taron, Ulrike Helene A1 - Barlow, Axel A1 - Wang, Linying A1 - Hu, Jiaming A1 - Deng, Tao A1 - Lai, Xulong A1 - Hofreiter, Michael T1 - Mitochondrial genomes of Late Pleistocene caballine horses from China belong to a separate clade JF - Quaternary science reviews : the international multidisciplinary research and review journal N2 - There were several species of Equus in northern China during the Late Pleistocene, including Equus przewalskii and Equus dalianensis. A number of morphological studies have been carried out on E. przewalskii and E. dalianensis, but their evolutionary history is still unresolved. In this study, we retrieved near-complete mitochondrial genomes from E. dalianensis and E. przewalskii specimens excavated from Late Pleistocene strata in northeastern China. Phylogenetic analyses revealed that caballoid horses were divided into two subclades: the New World and the Old World caballine horse subclades. The Old World caballine horses comprise of two deep phylogenetic lineages, with modern and ancient Equus caballus and modern E. przewalskii forming lineage I, and the individuals in this study together with one Yakut specimen forming lineage II. Our results indicate that Chinese Late Pleistocene caballoid horses showed a closer relationship to other Eurasian caballine horses than that to Pleistocene horses from North America. In addition, phylogenetic analyses suggested a close relationship between E. dalianensis and the Chinese fossil E. przewalskii, in agreement with previous researches based on morphological analyses. Interestingly, E. dalianensis and the fossil E. przewalskii were intermixed rather than split into distinct lineages, suggesting either that gene flow existed between these two species or that morphology-based species assignment of palaeontological specimens is not always correct. Moreover, Bayesian analysis showed that the divergence time between the New World and the Old World caballoid horses was at 1.02 Ma (95% CI: 0.86-1.24 Ma), and the two Old World lineages (I & II) split at 0.88 Ma (95% CI: 0.69-1.13 Ma), which indicates that caballoid horses seem to have evolved into different populations in the Old World soon after they migrated from North America via the Bering Land Bridge. Finally, the TMRCA of E. dalianensis was estimated at 0.20 Ma (95% CI: 0.15-0.28 Ma), and it showed a relative low genetic diversity compared with other Equus species. KW - Equus dalianensis KW - Equus przewalskii KW - Pleistocene caballine horses KW - ancient DNA KW - phylogenetic relationship KW - divergence time Y1 - 2020 U6 - https://doi.org/10.1016/j.quascirev.2020.106691 SN - 0277-3791 VL - 250 PB - Elsevier CY - Amsterdam [u.a.] ER - TY - JOUR A1 - Yuan, Jun-Xia A1 - Hou, Xin-Dong A1 - Barlow, Axel A1 - Preick, Michaela A1 - Taron, Ulrike H. A1 - Alberti, Federica A1 - Basler, Nikolas A1 - Deng, Tao A1 - Lai, Xu-Long A1 - Hofreiter, Michael A1 - Sheng, Gui-Lian T1 - Molecular identification of late and terminal Pleistocene Equus ovodovi from northeastern China JF - PLOS ONE N2 - The extant diversity of horses (family Equidae) represents a small fraction of that occurring over their evolutionary history. One such lost lineage is the subgenus Sussemionus, which is thought to have become extinct during the Middle Pleistocene. However, recent molecular studies and morphological analysis have revealed that one of their representatives, E. ovodovi, did exist in Siberia during the Late Pleistocene. Fossil materials of E. ovodovi have thus far only been found in Russia. In this study, we extracted DNA from three equid fossil specimens excavated from northeastern China dated at 12,770-12,596, 29,525-28,887 and 40,201-38,848 cal. yBP, respectively, and retrieved three near-complete mitochondrial genomes from the specimens. Phylogenetic analyses cluster the Chinese haplotypes together with previously published Russian E. ovodovi, strongly supporting the assignment of these samples to this taxon. The molecular identification of E. ovodovi in northeastern China extends the known geographical range of this fossil species by several thousand kilometers to the east. The estimated coalescence time of all E. ovodovi haplotypes is approximately 199 Kya, with the Chinese haplotypes coalescing approximately 130 Kya. With a radiocarbon age of 12,770-12,596 cal. yBP, the youngest sample in this study represents the first E. ovodovi sample dating to the terminal Pleistocene, moving the extinction date of this species forwards considerably compared to previously documented fossils. Overall, comparison of our three mitochondrial genomes with the two published ones suggests a genetic diversity similar to several extant species of the genus Equus. Y1 - 2019 U6 - https://doi.org/10.1371/journal.pone.0216883 SN - 1932-6203 VL - 14 IS - 5 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Wuster, Wolfgang A1 - Chirio, Laurent A1 - Trape, Jean-Francois A1 - Ineich, Ivan A1 - Jackson, Kate A1 - Greenbaum, Eli A1 - Barron, Cesar A1 - Kusamba, Chifundera A1 - Nagy, Zoltan T. A1 - Storey, Richard A1 - Hall, Cara A1 - Wuster, Catharine E. A1 - Barlow, Axel A1 - Broadley, Donald G. T1 - Integration of nuclear and mitochondrial gene sequences and morphology reveals unexpected diversity in the forest cobra (Naja melanoleuca) species complex in Central and West Africa (Serpentes: Elapidae) JF - Zootaxa : an international journal of zootaxonomy ; a rapid international journal for animal taxonomists N2 - Cobras are among the most widely known venomous snakes, and yet their taxonomy remains incompletely understood, particularly in Africa. Here, we use a combination of mitochondrial and nuclear gene sequences and morphological data to diagnose species limits within the African forest cobra, Naja (Boulengerina) melanoleuca. Mitochondrial DNA sequences reveal deep divergences within this taxon. Congruent patterns of variation in mtDNA, nuclear genes and morphology support the recognition of five separate species, confirming the species status of N. subfulva and N. peroescobari, and revealing two previously unnamed West African species, which are described as new: Naja (Boulengerina) guineensis sp. nov. Broadley, Trape, Chirio, Ineich & Wuster, from the Upper Guinea forest of West Africa, and Naja (Boulengerina) savannula sp. nov. Broadley, Trape, Chirio & Wuster, a banded form from the savanna-forest mosaic of the Guinea and Sudanian savannas of West Africa. The discovery of cryptic diversity in this iconic group highlights our limited understanding of tropical African biodiversity, hindering our ability to conserve it effectively. KW - Integrative taxonomy KW - Africa KW - Naja melanoleuca KW - Naja guineensis sp nov. KW - Naja savannula sp nov. KW - Elapidae KW - systematics Y1 - 2018 U6 - https://doi.org/10.11646/zootaxa.4455.1.3 SN - 1175-5326 SN - 1175-5334 VL - 4455 IS - 1 SP - 68 EP - 98 PB - Magnolia Press CY - Auckland ER - TY - JOUR A1 - Westbury, Michael V. A1 - Hartmann, Stefanie A1 - Barlow, Axel A1 - Wiesel, Ingrid A1 - Leo, Viyanna A1 - Welch, Rebecca A1 - Parker, Daniel M. A1 - Sicks, Florian A1 - Ludwig, Arne A1 - Dalen, Love A1 - Hofreiter, Michael T1 - Extended and continuous decline in effective population size results in low genomic diversity in the world's rarest hyena species, the brown hyena JF - Molecular biology and evolution N2 - 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 similar to 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. KW - evolution KW - hyena KW - genomics KW - population genomics KW - diversity Y1 - 2018 U6 - https://doi.org/10.1093/molbev/msy037 SN - 0737-4038 SN - 1537-1719 VL - 35 IS - 5 SP - 1225 EP - 1237 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Westbury, Michael V. A1 - Baleka, Sina Isabelle A1 - Barlow, Axel A1 - Hartmann, Stefanie A1 - Paijmans, Johanna L. A. A1 - Kramarz, Alejandro A1 - Forasiepi, Analia M. A1 - Bond, Mariano A1 - Gelfo, Javier N. A1 - Reguero, Marcelo A. A1 - Lopez-Mendoza, Patricio A1 - Taglioretti, Matias A1 - Scaglia, Fernando A1 - Rinderknecht, Andres A1 - Jones, Washington A1 - Mena, Francisco A1 - Billet, Guillaume A1 - de Muizon, Christian A1 - Luis Aguilar, Jose A1 - MacPhee, Ross D. E. A1 - Hofreiter, Michael T1 - A mitogenomic timetree for Darwin’s enigmatic South American mammal Macrauchenia patachonica JF - Nature Communications N2 - The unusual mix of morphological traits displayed by extinct South American native ungulates (SANUs) confounded both Charles Darwin, who first discovered them, and Richard Owen, who tried to resolve their relationships. Here we report an almost complete mitochondrial genome for the litoptern Macrauchenia. Our dated phylogenetic tree places Macrauchenia as sister to Perissodactyla, but close to the radiation of major lineages within Laurasiatheria. This position is consistent with a divergence estimate of B66Ma (95% credibility interval, 56.64-77.83 Ma) obtained for the split between Macrauchenia and other Panperissodactyla. Combined with their morphological distinctiveness, this evidence supports the positioning of Litopterna (possibly in company with other SANU groups) as a separate order within Laurasiatheria. We also show that, when using strict criteria, extinct taxa marked by deep divergence times and a lack of close living relatives may still be amenable to palaeogenomic analysis through iterative mapping against more distant relatives. Y1 - 2017 U6 - https://doi.org/10.1038/ncomms15951 SN - 2041-1723 VL - 8 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Thorpe, Roger A1 - Barlow, Axel A1 - Surget-Groba, Yann A1 - Malhotra, Anita T1 - Multilocus phylogeny, species age and biogeography of the Lesser Antillean anoles JF - Molecular phylogenetics and evolution N2 - 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. KW - Anolis KW - Multilocus phylogeny KW - Lesser antilles KW - Species age KW - Species turnover KW - Island colonization Y1 - 2018 U6 - https://doi.org/10.1016/j.ympev.2018.06.014 SN - 1055-7903 SN - 1095-9513 VL - 127 SP - 682 EP - 695 PB - Elsevier CY - San Diego ER - TY - JOUR A1 - Taron, Ulrike H. A1 - Lell, Moritz A1 - Barlow, Axel A1 - Paijmans, Johanna L. A. T1 - Testing of Alignment Parameters for Ancient Samples BT - Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool JF - Genes N2 - High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present ‘TAPAS’, (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material. KW - ancient DNA KW - short-read mapping KW - palaeogenomics KW - alignment sensitivity / specificity Y1 - 2018 U6 - https://doi.org/10.3390/genes9030157 SN - 2073-4425 VL - 9 IS - 3 SP - 1 EP - 12 PB - Molecular Diversity Preservation International CY - Basel ER - TY - JOUR A1 - Taron, Ulrike H. A1 - Lell, Moritz A1 - Barlow, Axel A1 - Paijmans, Johanna L. A. T1 - Testing of Alignment Parameters for Ancient Samples BT - Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool JF - Genese N2 - High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present ‘TAPAS’, (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material. KW - ancient DNA KW - short-read mapping KW - palaeogenomics KW - paleogenomics KW - alignment sensitivity/specificity Y1 - 2018 U6 - https://doi.org/10.3390/genes9030157 SN - 2073-4425 VL - 9 IS - 3 PB - MDPI CY - Basel ER - TY - JOUR A1 - Sheng, Gui-Lian A1 - Basler, Nikolas A1 - Ji, Xue-Ping A1 - Paijmans, Johanna L. A. A1 - Alberti, Federica A1 - Preick, Michaela A1 - Hartmann, Stefanie A1 - Westbury, Michael V. A1 - Yuan, Jun-Xia A1 - Jablonski, Nina G. A1 - Xenikoudakis, Georgios A1 - Hou, Xin-Dong A1 - Xiao, Bo A1 - Liu, Jian-Hui A1 - Hofreiter, Michael A1 - Lai, Xu-Long A1 - Barlow, Axel T1 - Paleogenome reveals genetic contribution of extinct giant panda to extant populations JF - Current biology N2 - Historically, the giant panda was widely distributed from northern China to southwestern Asia [1]. As a result of range contraction and fragmentation, extant individuals are currently restricted to fragmented mountain ranges on the eastern margin of the Qinghai-Tibet plateau, where they are distributed among three major population clusters [2]. However, little is known about the genetic consequences of this dramatic range contraction. For example, were regions where giant pandas previously existed occupied by ancestors of present-day populations, or were these regions occupied by genetically distinct populations that are now extinct? If so, is there any contribution of these extinct populations to the genomes of giant pandas living today? To investigate these questions, we sequenced the nuclear genome of an similar to 5,000-year-old giant panda from Jiangdongshan, Teng-chong County in Yunnan Province, China. We find that this individual represents a genetically distinct population that diverged prior to the diversification of modern giant panda populations. We find evidence of differential admixture with this ancient population among modern individuals originating from different populations as well as within the same population. We also find evidence for directional gene flow, which transferred alleles from the ancient population into the modern giant panda lineages. A variable proportion of the genomes of extant individuals is therefore likely derived from the ancient population represented by our sequenced individual. Although extant giant panda populations retain reasonable genetic diversity, our results suggest that this represents only part of the genetic diversity this species harbored prior to its recent range contractions. Y1 - 2019 U6 - https://doi.org/10.1016/j.cub.2019.04.021 SN - 0960-9822 SN - 1879-0445 VL - 29 IS - 10 SP - 1695 EP - 1700 PB - Cell Press CY - Cambridge ER -