Phylogenomic insights into the diversity and evolution of Palearctic vipers.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of âˆ¼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Genomic basis for skin phenotype and cold adaptation in the extinct Steller's sea cow.

Steller's sea cow, an extinct sirenian and one of the largest Quaternary mammals, was described by Georg Steller in 1741 and eradicated by humans within 27 years. Here, we complement Steller's descriptions with paleogenomic data from 12 individuals. We identified convergent evolution between Steller's sea cow and cetaceans but not extant sirenians, suggesting a role of several genes in adaptation to cold aquatic (or marine) environments. Among these are inactivations of lipoxygenase genes, which in humans and mouse models cause ichthyosis, a skin disease characterized by a thick, hyperkeratotic epidermis that recapitulates Steller's sea cows' reportedly bark-like skin. We also found that Steller's sea cows' abundance was continuously declining for tens of thousands of years before their description, implying that environmental changes also contributed to their extinction.

Natural history collections recapitulate 200 years of faunal change.

Changing species assemblages represent major challenges to ecosystems around the world. Retracing these changes is limited by our knowledge of past biodiversity. Natural history collections represent archives of biodiversity and are therefore an unparalleled source to study biodiversity changes. In the present study, we tested the value of natural history collections for reconstructing changes in the abundance and presence of species over time. In total, we scrutinized 17 080 quality-checked records for 242 epibenthic invertebrate species from the North and Baltic Seas collected throughout the last 200 years. Our approaches identified eight previously reported species introductions, 10 range expansions, six of which are new to science, as well as the long-term decline of 51 marine invertebrate species. The cross-validation of our results with published accounts of endangered species and neozoa of the area confirmed the results for two of the approaches for 49 to 55% of the identified species, and contradicted our results for 9 to 10%. The results based on relative record trends were less validated. We conclude that, with the proper approaches, natural history collections are an unmatched resource for recovering early species introductions and declines.

Tail-propelled aquatic locomotion in a theropod dinosaur.

In recent decades, intensive research on non-avian dinosaurs has strongly suggested that these animals were restricted to terrestrial environments1. Historical proposals that some groups, such as sauropods and hadrosaurs, lived in aquatic environments2,3 were abandoned decades ago4-6. It has recently been argued that at least some of the spinosaurids-an unusual group of large-bodied theropods of the Cretaceous era-were semi-aquatic7,8, but this idea has been challenged on anatomical, biomechanical and taphonomic grounds, and remains controversial9-11. Here we present unambiguous evidence for an aquatic propulsive structure in a dinosaur, the giant theropod Spinosaurus aegyptiacus7,12. This dinosaur has a tail with an unexpected and unique shape that consists of extremely tall neural spines and elongate chevrons, which forms a large, flexible fin-like organ capable of extensive lateral excursion. Using a robotic flapping apparatus to measure undulatory forces in physical models of different tail shapes, we show that the tail shape of Spinosaurus produces greater thrust and efficiency in water than the tail shapes of terrestrial dinosaurs and that these measures of performance are more comparable to those of extant aquatic vertebrates that use vertically expanded tails to generate forward propulsion while swimming. These results are consistent with the suite of adaptations for an aquatic lifestyle and piscivorous diet that have previously been documented for Spinosaurus7,13,14. Although developed to a lesser degree, aquatic adaptations are also found in other members of the spinosaurid clade15,16, which had a near-global distribution and a stratigraphic range of more than 50 million years14, pointing to a substantial invasion of aquatic environments by dinosaurs.

Phylogenetic differentiation and taxonomic consequences in the Saurodactylus brosseti species complex (Squamata: Sphaerodactylidae), with description of four new species.

We divide Saurodactylus brosseti into five species, based on molecular and morphological characteristics. Phylogenetic analysis of two mitochondrial markers (ND4 and 12S rRNA) reveals five separately evolving lineages (clades) in the S. brosseti species complex. Each can be assigned to a particular geographical range within the distribution area in Morocco: North, South, Anti-Atlas, East and Southeast. A dichotomous key for all species within the genus Saurodactylus is presented. Despite the high level of color pattern variation among and within the northern, southern, Anti-Atlas and eastern clades, the southeastern lineage shows a unique coloration which is not found in other clades. The ecological traits of the southeastern clade seem to differ from those of the other four lineages since it is basically associated with Acacia species. Like the relictual ecosystem characterized by these trees, it appears to be highly threatened by extinction.

Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards.

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.

Allopatric diversification and evolutionary melting pot in a North African Palearctic relict: The biogeographic history of Salamandra algira.

North Africa is a climatically and topographically complex region with unique biotic assemblages resulting from the combination of multiple biogeographic realms. Here, we assess the role of climate in promoting intra-specific diversification in a Palearctic relict, the North African fire salamander, Salamandra algira, using a combination of phylogenetic and population genetic analyses, paleoclimatic modelling and niche overlap tests. We used mitochondrial DNA (Cyt-b), 9838 ddRADseq loci, and 14 microsatellite loci to characterize patterns of genetic diversity and population structure. Phylogenetic analyses recover two major clades, each including several lineages with mito-nuclear discordances suggesting introgressive patterns between lineages in the Middle Atlas, associated with a melting pot of genetic diversity. Paleoclimatic modelling identified putative climatic refugia, largely matching areas of high genetic diversity, and supports the role of aridity in promoting allopatric diversification associated with ecological niche conservatism. Overall, our results highlight the role of climatic microrefugia as drivers of populations' persistence and diversification in the face of climatic oscillations in North Africa, and stress the importance of accounting for different genomic regions when reconstructing biogeographic processes from molecular markers.

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.

Optimized DNA sampling of ancient bones using Computed Tomography scans.

The prevalence of contaminant microbial DNA in ancient bone samples represents the principal limiting factor for palaeogenomic studies, as it may comprise more than 99% of DNA molecules obtained. Efforts to exclude or reduce this contaminant fraction have been numerous but also variable in their success. Here, we present a simple but highly effective method to increase the relative proportion of endogenous molecules obtained from ancient bones. Using computed tomography (CT) scanning, we identify the densest region of a bone as optimal for sampling. This approach accurately identifies the densest internal regions of petrous bones, which are known to be a source of high-purity ancient DNA. For ancient long bones, CT scans reveal a high-density outermost layer, which has been routinely removed and discarded prior to DNA extraction. For almost all long bones investigated, we find that targeted sampling of this outermost layer provides an increase in endogenous DNA content over that obtained from softer, trabecular bone. This targeted sampling can produce as much as 50-fold increase in the proportion of endogenous DNA, providing a directly proportional reduction in sequencing costs for shotgun sequencing experiments. The observed increases in endogenous DNA proportion are not associated with any reduction in absolute endogenous molecule recovery. Although sampling the outermost layer can result in higher levels of human contamination, some bones were found to have more contamination associated with the internal bone structures. Our method is highly consistent, reproducible and applicable across a wide range of bone types, ages and species. We predict that this discovery will greatly extend the potential to study ancient populations and species in the genomics era.

The evolutionary and phylogeographic history of woolly mammoths: a comprehensive mitogenomic analysis.

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

The mitochondrial genomes of Atlas Geckos (Quedenfeldtia): mitogenome assembly from transcriptomes and anchored hybrid enrichment datasets.

The nearly complete mitogenomes of the two species of North African Atlas geckos, Quedenfeldtia moerens and Q. trachyblepharus were assembled from anchored hybrid enrichment data and RNAseq data. Congruent assemblies were obtained for four samples included in both datasets. We recovered the 13 protein-coding genes, 22 tRNA genes, and two rRNA genes for both species, including partial control region. The order of genes agrees with that of other geckos.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East.

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

Hybrid origin of European Vipers (Vipera magnifica and Vipera orlovi) from the Caucasus determined using genomic scale DNA markers.

BACKGROUND: Studying patterns of introgression can illuminate the role of hybridization in speciation, and help guide decisions relevant to the conservation of rare taxa. Vipera magnifica and Vipera orlovi are small vipers that have high conservation status due to their rarity and restricted distributions in an area of the Caucasus region where two other related species are present - V. kaznakovi and V. renardi. Despite numerous observations of hybridization between different species of small vipers, and the potential of a hybrid origin for V. magnifica and V. orlovi based on their distribution with respect to V. kaznakovi and V. renardi, hypotheses of a hybrid origin have not been formally tested. Here we generate genomic-scale data by performing next generation sequencing of double digest restriction-site associated DNA libraries, and use these multilocus data to test whether these two species are of hybrid origin. RESULTS: We generated over nine hundred loci for 38 specimens of six taxa, and analysed the dataset using Bayesian clustering and multivariate methods, as well as Patterson D-statistics, which can distinguish between incomplete lineage sorting and introgression as explanations for shared polymorphism. The results demonstrate a pattern of historical admixture in the two purported hybrids that is consistent with past gene flow from V. renardi into V. kaznakovi. The average admixture proportion in individuals was low (6.39 %) in the case of V. magnifica, but was higher in V. orlovi (19.02 %). We also show that the specific individual samples used in D-statistic tests can have a significant impact on inferences regarding the magnitude of introgression, suggesting the importance of including multiple individuals in these analyses. CONCLUSIONS: Our results support the conclusion that both V. orlovi and V. magnifica had formed through a hybridization event between V. kaznakovi and V. renardi. Given a low proportion of admixture and absence of clear ecological and morphological differences V. magnifica should be treated as a marginal population of V. kaznakovi. Further studies that include analyses of ecological segregation of V. orlovi from parental taxa and search for evolutionary consequences of hybridisation would clarify if V. orlovi is a distinct hybrid species. Until this we recommend preserving the current taxonomy and protection status of V. orlovi.

Mitochondrial phylogeny shows multiple independent ecological transitions and northern dispersion despite of Pleistocene glaciations in meadow and steppe vipers (Vipera ursinii and Vipera renardi).

The phylogeny and historical demography of small Eurasian vipers of the Vipera ursinii and V. renardi complexes were studied using mitochondrial DNA sequences analysed with Bayesian inference, Maximum Likelihood and Maximum Parsimony approaches, and mismatch distributions. Diversification in the group resulted from an initial dispersion in the later Pliocene - Pleistocene in two directions: north-westwards via the Balkans (V. ursinii complex) and north-eastwards from Asia Minor via the Caucasus (V. renardi complex). An independent, comparatively recent transition occurred from montane habitats to lowland grasslands in different mitochondrial lineages during the Late Pleistocene, when representatives of the both complexes had reached lowland steppes to the north. Effective population size showed clear signs of rapid growth in eastern V. renardi, triggered by colonization of vast lowland steppes, but in western V. ursinii complex grew during the Last Glaciation and experienced stabilization in Holocene. Expansion and population growth in lowland lineages of V. renardi was not strongly affected by Pleistocene climatic oscillations, when cold, dry conditions could have favoured species living in open grasslands. The high diversity of closely related haplotypes in the Caucasus and Tien-Shan could have resulted from repetitive expansion-constriction-isolation events in montane regions during Pleistocene climate fluctuations. The mitochondrial phylogeny pattern conflicts with the current taxonomy.

New ophthalmosaurid ichthyosaurs from the European Lower Cretaceous demonstrate extensive ichthyosaur survival across the Jurassic-Cretaceous boundary.

BACKGROUND: Ichthyosauria is a diverse clade of marine amniotes that spanned most of the Mesozoic. Until recently, most authors interpreted the fossil record as showing that three major extinction events affected this group during its history: one during the latest Triassic, one at the Jurassic-Cretaceous boundary (JCB), and one (resulting in total extinction) at the Cenomanian-Turonian boundary. The JCB was believed to eradicate most of the peculiar morphotypes found in the Late Jurassic, in favor of apparently less specialized forms in the Cretaceous. However, the record of ichthyosaurs from the Berriasian-Barremian interval is extremely limited, and the effects of the end-Jurassic extinction event on ichthyosaurs remains poorly understood. METHODOLOGY/PRINCIPAL FINDINGS: Based on new material from the Hauterivian of England and Germany and on abundant material from the Cambridge Greensand Formation, we name a new ophthalmosaurid, Acamptonectes densus gen. et sp. nov. This taxon shares numerous features with Ophthalmosaurus, a genus now restricted to the Callovian-Berriasian interval. Our phylogenetic analysis indicates that Ophthalmosauridae diverged early in its history into two markedly distinct clades, Ophthalmosaurinae and Platypterygiinae, both of which cross the JCB and persist to the late Albian at least. To evaluate the effect of the JCB extinction event on ichthyosaurs, we calculated cladogenesis, extinction, and survival rates for each stage of the Oxfordian-Barremian interval, under different scenarios. The extinction rate during the JCB never surpasses the background extinction rate for the Oxfordian-Barremian interval and the JCB records one of the highest survival rates of the interval. CONCLUSIONS/SIGNIFICANCE: There is currently no evidence that ichthyosaurs were affected by the JCB extinction event, in contrast to many other marine groups. Ophthalmosaurid ichthyosaurs remained diverse from their rapid radiation in the Middle Jurassic to their total extinction at the beginning of the Late Cretaceous.

Didactyl tracks of paravian theropods (Maniraptora) from the ?Middle Jurassic of Africa.

BACKGROUND: A new dinosaur tracksite from ?Middle Jurassic sediments of the Irhazer Group on the plains of Agadez (Rep. Niger, northwest Africa) revealed extraordinarily well preserved didactyl tracks of a digitigrade bipedal trackmaker. The distinct morphology of the pes imprints indicates a theropod trackmaker from a paravian maniraptoran closely related to birds. METHODOLOGY/PRINCIPAL FINDINGS: The early age and the morphological traits of the tracks allow for description of the new ichnotaxon Paravipus didactyloides. A total of 120 tracks are assigned to 5 individual trackways. The 'medium-sized' tracks with an average footprint length of 27.5 cm and footprint width of 23.1 cm are deeply imprinted into the track bearing sandstone. CONCLUSIONS/SIGNIFICANCE: A comparison with other didactyl tracks gives new insights into the foot morphology of advanced maniraptoran theropods and contributes to knowledge of their evolutionary history. The new ichnotaxon takes an important position in the ichnological fossil record of Gondwana and the mid-Jurassic biota worldwide, because it is among the earliest known records of paravian maniraptorans and of didactyl theropod tracks from Africa.

A new basal sauropod dinosaur from the middle Jurassic of Niger and the early evolution of sauropoda.

BACKGROUND: The early evolution of sauropod dinosaurs is poorly understood because of a highly incomplete fossil record. New discoveries of Early and Middle Jurassic sauropods have a great potential to lead to a better understanding of early sauropod evolution and to reevaluate the patterns of sauropod diversification. PRINCIPAL FINDINGS: A new sauropod from the Middle Jurassic of Niger, Spinophorosaurus nigerensis n. gen. et sp., is the most complete basal sauropod currently known. The taxon shares many anatomical characters with Middle Jurassic East Asian sauropods, while it is strongly dissimilar to Lower and Middle Jurassic South American and Indian forms. A possible explanation for this pattern is a separation of Laurasian and South Gondwanan Middle Jurassic sauropod faunas by geographic barriers. Integration of phylogenetic analyses and paleogeographic data reveals congruence between early sauropod evolution and hypotheses about Jurassic paleoclimate and phytogeography. CONCLUSIONS: Spinophorosaurus demonstrates that many putatively derived characters of Middle Jurassic East Asian sauropods are plesiomorphic for eusauropods, while South Gondwanan eusauropods may represent a specialized line. The anatomy of Spinophorosaurus indicates that key innovations in Jurassic sauropod evolution might have taken place in North Africa, an area close to the equator with summer-wet climate at that time. Jurassic climatic zones and phytogeography possibly controlled early sauropod diversification.

When continents collide: phylogeny, historical biogeography and systematics of the medically important viper genus Echis (Squamata: Serpentes: Viperidae).

We analyze the phylogeny of the medically important and taxonomically unresolved viper genus Echis using four mitochondrial gene fragments. The results show that the populations of the genus fall into four main clades: the Echis carinatus, E. coloratus, E. ocellatus and E. pyramidum groups. The E. pyramidum and E. coloratus groups are sister taxa but the interrelationships of this clade and the E. ocellatus and E. carinatus groups are unresolved. The initial divergence of the genus appears to coincide with the collision between Afro-Arabia and Eurasia, and that between the E. coloratus and E. pyramidum clades appears to be associated with the opening of the Red Sea. Later land connections between Africa and Arabia may have contributed to shaping the distribution of the E. pyramidum complex. The present distribution of E. carinatus may be the result of range expansion from southern India. Taxonomically, our results provide molecular evidence for the validity of Echis omanensis, E. khosatzkii, E. borkini and E. jogeri, for the presence of unsuspected genetic diversity within the E. pyramidum complex in eastern Africa, and for the conspecificity of E. carinatus and E. multisquamatus. The status of E. leucogaster remains to be confirmed.

Multiple colonization of Madagascar and Socotra by colubrid snakes: evidence from nuclear and mitochondrial gene phylogenies.

Colubrid snakes form a speciose group of unclarified phylogeny. Their almost cosmopolitan distribution could be interpreted as a product of plate-tectonic vicariance. We used sequences of the nuclear c-mos, the mitochondrial cytochrome b and the 16S rRNA genes in 41 taxa to elucidate the relationships between the endemic colubrid genera found in Madagascar and in the Socotra archipelago. The well-resolved trees indicate multiple origins of both the Malagasy and the Socotran taxa. The Malagasy genus Mimophis was nested within the Psammophiinae, and the Socotran Hemerophis was closely related to Old World representatives of the former genus Coluber. The remaining 14 genera of Malagasy colubrids formed a monophyletic sister group of the Socotran Ditypophis (together forming the Pseudoxyrhophiinae). Molecular-clock estimates place the divergence of Malagasy and Socotran colubrids from their non-insular sister groups into a time-frame between the Eocene and Miocene. Over-seas rafting is the most likely hypothesis for the origin of at least the Malagasy taxa. The discovery of a large monophyletic clade of colubrids endemic to Madagascar indicates a need for taxonomic changes. The relationship of this radiation to the Socotran Ditypophis highlights the potential of the Indian Ocean islands to act as an evolutionary reservoir for lineages that have become extinct in Africa and Asia.