Phylogenomics using formalin‐fixed and 100+ year‐old intractable natural history specimens

High-throughput DNA sequencing of museum specimens sheds light on the long-missing species of the Bokermannohyla claresignata group (Anura: Hylidae: Cophomantini)

The two species of the Bokermannohyla claresignata species group (Anura: Hylidae) have not been collected for the last four decades. It is the only species group of the hyline tribe Cophomantini that has not yet been analysed genetically. Its phylogenetic position is thus uncertain, and it has a combination of adult and larval character states that make this group a crucial missing piece that hinders our understanding of Cophomantini phylogenetics and character evolution. We obtained DNA sequences from a museum larval specimen of Bok. claresignata, using specialized extraction methods and high-throughput DNA sequencing, and combined the molecular phylogenetic results with available phenotypic information to provide new insights into the taxonomy and phylogenetic relationships of its species group. Our phylogenetic results place Bok. claresignata as sister to the Boana pulchella group, supporting its inclusion in Boana, together with Bokermannohyla clepsydra. In light of this new finding, we recognize a newly defined Boana claresignata group to accommodate these species, thus resolving both the polyphyly of Bokermannohyla and the paraphyly of Boana. Considering the phylogenetic relationships of the Boana claresignata group, we also discuss the evolution of suctorial tadpoles and mature oocyte/egg pigmentation in Cophomantini.

A new mouse of the Peromyscus maniculatus species complex (Cricetidae) from the highlands of central Mexico

The Peromyscus maniculatus species complex is one of the most widespread group of small mammals in North America. However, the taxonomy and phylogenetic relationships among its constituent taxa remain unclear. As part of a revision of Peromyscus specimens from the highlands of the Trans-Mexican Volcanic Belt in central Mexico, we identified five individuals collected in 1968 that differed externally from other Peromyscus specimens, although morphologically similar to P. labecula and P. melanotis, both latter in the P. maniculatus species complex. Based on cranial measurements and mitochondrial DNA sequences, we aimed to more accurately determine the phylogenetic relationships and the taxonomic status of these individuals. Molecular phylogenetic analyses showed that the specimens formed a monophyletic clade sister to the P. maniculatus species complex. Pairwise genetic distances between those specimens and other species within the P. maniculatus species complex were greater than 7.91%. In addition, morphological analyses clearly distinguished the test specimens from P. melanotis and P. labecula. Based on the results of our molecular and morphological analyses, we conclude that these specimens represent an undescribed species of the P. maniculatus species complex, which we describe herein.

A large-scale systematic framework of Chinese snakes based on a unified multilocus marker system

Snakes are one of the most diverse groups of terrestrial vertebrates, with approximately 3500 extant species. A robust phylogeny and taxonomy of snakes is crucial for us to know, study and protect them. For a large group such as snakes, broad-scale phylogenetic reconstructions largely rely on data integration. Increasing the compatibility of the data from different researches is thus important, which can be facilitated by standardization of the loci used in systematic analyses. In this study, we proposed a unified multilocus marker system for snake systematics by conflating 5 mitochondrial markers, 19 vertebrate-universal nuclear protein coding (NPC) markers and 72 snake-specific noncoding intron markers. This marker system is an addition to the large squamate conserved locus set (SqCL) for studies preferring a medium-scale data set. We applied this marker system to over 440 snake samples and constructed the currently most comprehensive systematic framework of the snakes in China. Robust snake phylogenetic relationships were recovered at both deep and shallow evolutionary depths, demonstrating the usefulness of this multilocus marker system. Discordance was revealed by a parallel comparison between the snake tree based on the multilocus marker system and that based on only the mitochondrial loci, highlighting the necessity of using multiple types of markers to better understand the snake evolutionary histories. The divergence times of different snake groups were estimated with the nuclear data set. Our comprehensive snake tree not only confirms many important nodes inferred in previous studies but also contributes new insights into many snake phylogenetic relationships. Suggestions are made for the current Chinese snake taxonomy.

Opening the door to greater phylogeographic inference in Southeast Asia: Comparative genomic study of five codistributed rainforest bird species using target capture and historical DNA

Indochina and Sundaland are biologically diverse, interconnected regions of Southeast Asia with complex geographic histories. Few studies have examined phylogeography of bird species that span the two regions because of inadequate population sampling. To determine how geographic barriers/events and disparate dispersal potential have influenced the population structure, gene flow, and demographics of species that occupy the entire area, we studied five largely codistributed rainforest bird species: Arachnothera longirostra, Irena puella, Brachypodius atriceps, Niltava grandis, and Stachyris nigriceps. We accomplished relatively thorough sampling and data collection by sequencing ultraconserved elements (UCEs) using DNA extracted from modern and older (historical) specimens. We obtained a genome-wide set of 753-4,501 variable loci and 3,919-18,472 single nucleotide polymorphisms. The formation of major within-species lineages occurred within a similar span of time (0.5-1.5 mya). Major patterns in population genetic structure are largely consistent with the dispersal potential and habitat requirements of the study species. A population break across the Isthmus of Kra was shared only by the two hill/submontane insectivores (N. grandis and S. nigriceps). Across Sundaland, there is little structure in B. atriceps, which is a eurytopic and partially frugivorous species that often utilizes forest edges. Two other eurytopic species, A. longirostra and I. puella, possess highly divergent populations in peripheral Sunda Islands (Java and/or Palawan) and India. These species probably possess intermediate dispersal abilities that allowed them to colonize new areas, and then remained largely isolated subsequently. We also observed an east-west break in Indochina that was shared by B. atriceps and S. nigriceps, species with very different habitat requirements and dispersal potential. By analyzing high-throughput DNA data, our study provides an unprecedented comparative perspective on the process of avian population divergence across Southeast Asia, a process that is determined by geography, species characteristics, and the stochastic nature of dispersal and vicariance events.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

The Global Museum: natural history collections and the future of evolutionary science and public education

Natural history museums are unique spaces for interdisciplinary research and educational innovation. Through extensive exhibits and public programming and by hosting rich communities of amateurs, students, and researchers at all stages of their careers, they can provide a place-based window to focus on integration of science and discovery, as well as a locus for community engagement. At the same time, like a synthesis radio telescope, when joined together through emerging digital resources, the global community of museums (the ‘Global Museum’) is more than the sum of its parts, allowing insights and answers to diverse biological, environmental, and societal questions at the global scale, across eons of time, and spanning vast diversity across the Tree of Life. We argue that, whereas natural history collections and museums began with a focus on describing the diversity and peculiarities of species on Earth, they are now increasingly leveraged in new ways that significantly expand their impact and relevance. These new directions include the possibility to ask new, often interdisciplinary questions in basic and applied science, such as in biomimetic design, and by contributing to solutions to climate change, global health and food security challenges. As institutions, they have long been incubators for cutting-edge research in biology while simultaneously providing core infrastructure for research on present and future societal needs. Here we explore how the intersection between pressing issues in environmental and human health and rapid technological innovation have reinforced the relevance of museum collections. We do this by providing examples as food for thought for both the broader academic community and museum scientists on the evolving role of museums. We also identify challenges to the realization of the full potential of natural history collections and the Global Museum to science and society and discuss the critical need to grow these collections. We then focus on mapping and modelling of museum data (including place-based approaches and discovery), and explore the main projects, platforms and databases enabling this growth. Finally, we aim to improve relevant protocols for the long-term storage of specimens and tissues, ensuring proper connection with tomorrow’s technologies and hence further increasing the relevance of natural history museums.

A new climbing shrew from Sulawesi highlights the tangled taxonomy of an endemic radiation

We describe a new species of Crocidura (Soricidae) from Sulawesi Island, Indonesia, documenting its novelty with both genetic and morphological characters. The new species is widespread on the island, with vouchered records from nine general localities distributed among five of the island’s areas of endemism. Morphologically, the new species is readily distinguished from all other described Sulawesi Crocidura by its intermediate body size, gray pelage, and long, hairy tail. The new species was mainly captured in pitfalls placed in the ground, but we also obtained evidence that it readily climbs trees and may be scansorial in its locomotor habits. Populations of the new species sampled from across the island are closely related, separated by < 0.02 uncorrected mitochondrial p-distances. The new species is one member of an endemic radiation of shrews on Sulawesi now known to contain six valid species and several undescribed species, all within the genus Crocidura. Resolution of species limits and phylogenetic relationships in this radiation is hindered by habitat loss at type localities, historical designation of new species using very small sample sizes, and a lack of genetic data from type specimens.

Kami mendeskripsikan spesies baru Crocidura (Soricidae) dari Pulau Sulawesi, Indonesia, sekaligus mendokumentasikan keunikan karakter secara genetik maupun morfologi dari spesies tersebut. Spesies baru ini tersebar luas di Pulau Sulawesi, diketahui berdasarkan spesimen yang berasal dari sembilan lokasi umum yang tersebar di lima kawasan endemik di pulau tersebut. Secara morfologi, spesies baru ini dapat dibedakan dari spesies Crocidura lainnya dari Sulawesi berdasarkan ukuran tubuh yang sedang, rambut tubuh berwarna abu-abu, dan ekor yang panjang dan berambut. Spesies baru ini sebagian besar diperoleh dari perangkap sumuran yang ditanam didalam tanah, selain itu kami juga mendapatkan bukti bahwa spesies ini mampu memanjat pohon dan kemungkinan memiliki perilaku sebagai pemanjat. Beberapa populasi spesies yang dikoleksi dari Sulawesi ini mempunyai kekerabatan yang dekat, hanya dipisahkan oleh jarak proporsi DNA mitokondria (tidak terkoreksi) sebesar < 0.02. Spesies baru ini merupakan salah satu anggota dari suatu kelompok radiasi endemik cecurut di Sulawesi yang sampai saat ini diketahui terdiri atas enam spesies yang valid, dan beberapa spesies yang belum dideskripsikan, semuanya termasuk didalam genus Crocidura. Kepastian dalam menetapkan batasan jarak antar spesies dan hubungan kekerabatan genetik dari radiasi kelompok cecurut di Sulawesi terkendala oleh kerusakan habitat pada lokasi spesimen tipe, sejarah penamaan spesies yang hanya berdasarkan sampel yang sedikit, dan keterbatasan data molekuler dari spesimen tipe.

Spider-specific probe set for ultraconserved elements offers new perspectives on the evolutionary history of spiders (Arachnida, Araneae)

Phylogenomic methods have proven useful for resolving deep nodes and recalcitrant groups in the spider tree of life. Across arachnids, transcriptomic approaches may generate thousands of loci, and target-capture methods, using the previously designed arachnid-specific probe set, can target a maximum of about 1,000 loci. Here, we develop a specialized target-capture probe set for spiders that contains over 2,000 ultraconserved elements (UCEs) and then demonstrate the utility of this probe set through sequencing and phylogenetic analysis. We designed the 'spider-specific' probe set using three spider genomes (Loxosceles, Parasteatoda and Stegodyphus) and ensured that the newly designed probe set includes UCEs from the previously designed Arachnida probe set. The new 'spider-specific' probes were used to sequence UCE loci in 51 specimens. The remaining samples included five spider genomes and taxa that were enriched using Arachnida probe set. The 'spider-specific' probes were also used to gather loci from a total of 84 representative taxa across Araneae. On mapping these 84 taxa to the Arachnida probe set, we captured at most 710 UCE loci, while the spider-specific probe set captured up to 1,547 UCE loci from the same taxon sample. Phylogenetic analyses using maximum likelihood and coalescent methods corroborate most nodes resolved by recent transcriptomic analyses, but not all (e.g. UCE data suggest monophyly of 'symphytognathoids'). Our preferred hypothesis based on topology tests, suggests monophyly of the 'symphytognathoids' (the miniature orb weavers), which in previous studies has only been supported by a combination of morphological and behavioural characters.

More than skin and bones: Comparing extraction methods and alternative sources of DNA from avian museum specimens

Next-generation sequencing has greatly expanded the utility and value of museum collections by revealing specimens as genomic resources. As the field of museum genomics grows, so does the need for extraction methods that maximize DNA yields. For avian museum specimens, the established method of extracting DNA from toe pads works well for most specimens. However, for some specimens, especially those of birds that are very small or very large, toe pads can be a poor source of DNA. In this study, we apply two DNA extraction methods (phenol-chloroform and silica column) to three different sources of DNA (toe pad, skin punch and bone) from 10 historical avian museum specimens. We show that a modified phenol-chloroform protocol yielded significantly more DNA than a silica column protocol (e.g., Qiagen DNeasy Blood & Tissue Kit) across all tissue types. However, extractions using the silica column protocol contained longer fragments on average than those using the phenol-chloroform protocol, probably as a result of loss of small fragments through the silica column. While toe pads yielded more DNA than skin punches and bone fragments, skin punches proved to be a reliable alternative source of DNA and might be especially appealing when toe pad extractions are impractical. Overall, we found that historical bird museum specimens contain substantial amounts of DNA for genomic studies under most extraction scenarios, but that a phenol-chloroform protocol consistently provides the high quantities of DNA required for most current genomic protocols.

Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives

BACKGROUND

Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections.

RESULTS

The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths).

CONCLUSIONS

Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.

Species complexes and the importance of Data Deficient classification in Red List assessments: The case of Hylobatrachus frogs

Taxonomy is the cornerstone of extinction risk assessments. Currently, the IUCN Red List treats species complexes either under a single overarching species name—resulting in an unhelpfully broad circumscription and underestimated threat assessment that does not apply to any one species lineage—or omits them altogether—resulting in the omission of species that should be assessed. We argue that taxonomic uncertainty alone, as in species complexes, should be grounds for assessment as Data Deficient (DD). Yet, use of the DD category is currently discouraged, resulting in assessments based on poor data quality and dismissal of the importance of taxonomic confidence in conservation. This policy may be leading to volatile and unwarranted assessments of hundreds of species across the world, and needs to be revised. To illustrate this point, we here present a partial taxonomic revision of torrent frogs from eastern Madagascar in the Mantidactylus subgenus Hylobatrachus. Two named species, Mantidactylus (Hylobatrachus) lugubris and M. (H.) cowanii, and several undescribed candidate species are recognised, but the application of the available names has been somewhat ambiguous. In a recent re-assessment of its conservation status, M. (H.) lugubris was assessed including all complex members except M. (H.) cowanii within its distribution, giving it a status of Least Concern and distribution over most of eastern Madagascar. After describing two of the unnamed lineages as Mantidactylus (Hylobatrachus) atsimo sp. nov. (from southeastern Madagascar) and Mantidactylus (Hylobatrachus) petakorona sp. nov. (from the Marojejy Massif in northeastern Madagascar), we show that Mantidactylus (Hylobatrachus) lugubris is restricted to the central east of Madagascar, highlighting the inaccuracy of its current Red List assessment. We propose to re-assess its status under a more restrictive definition that omits well-defined candidate species, thus representing the actual species to which its assessment refers, to the best of current knowledge. We recommend that for species complexes in general, (1) nominal lineages that can be confidently restricted should be assessed under the strict definition, (2) non-nominal species-level lineages and ambiguous names should be prioritised for taxonomic research, and (3) ambiguous names should be assessed as DD to highlight the deficiency in data on their taxonomic status, which is an impediment to their conservation. This would reduce ambiguity and underestimation of threats involved in assessing species complexes, and place the appropriate emphasis on the importance of taxonomy in anchoring conservation.

Revealing paraphyly and placement of extinct species within Epischura (Copepoda: Calanoida) using molecular data and quantitative morphometrics

Epischura (Calanoida: Temoridae) is a Holarctic group of copepods serving important ecological roles, but it is difficult to study because of small range sizes of individual species and widespread distribution of the genus. This genus includes Tertiary relicts, some endemic to single, isolated lakes and can play major roles in unique ecosystems like Lakes Baikal and Tahoe. We present the first molecular and morphological analysis of Epischura that reveals their spatio-temporal evolutionary history. Morphological measurements of mandibles and genetics estimated phylogenetic relationships among all species represented in Epischura, including E. massachusettsensis, whose extinction status is of concern. Analyses used three gene regions for six previously unsequenced species to infer highly-resolved and well-supported phylogenies confirming a split between Siberian and North American species. Previously published age estimates and sequence data from broad taxonomic sampling of calanoid copepods estimated divergence times between the two Epischura groups. Divergence time estimates for Epischura were consistent with earlier molecular clock estimates and late-Miocene cooling events. Additionally, we provide the first taxonomically broad estimates of divergence times within Calanoida. The paraphyletic nature of the genus Epischura (and the family Temoridae) is apparent and requires the resurrection of the genus Epischurella (Smirnov, 1936) to describe the Siberian species.

A Phylogenomic Supertree of Birds

It has long been appreciated that analyses of genomic data (e.g., whole genome sequencing or sequence capture) have the potential to reveal the tree of life, but it remains challenging to move from sequence data to a clear understanding of evolutionary history, in part due to the computational challenges of phylogenetic estimation using genome-scale data. Supertree methods solve that challenge because they facilitate a divide-and-conquer approach for large-scale phylogeny inference by integrating smaller subtrees in a computationally efficient manner. Here, we combined information from sequence capture and whole-genome phylogenies using supertree methods. However, the available phylogenomic trees had limited overlap so we used taxon-rich (but not phylogenomic) megaphylogenies to weave them together. This allowed us to construct a phylogenomic supertree, with support values, that included 707 bird species (~7% of avian species diversity). We estimated branch lengths using mitochondrial sequence data and we used these branch lengths to estimate divergence times. Our time-calibrated supertree supports radiation of all three major avian clades (Palaeognathae, Galloanseres, and Neoaves) near the Cretaceous-Paleogene (K-Pg) boundary. The approach we used will permit the continued addition of taxa to this supertree as new phylogenomic data are published, and it could be applied to other taxa as well.

Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce. In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open‐source UCE probe set for beetles of the suborder Adephaga.

Large-scale molecular phylogeny, morphology, divergence-time estimation, and the fossil record of advanced caenophidian snakes (Squamata: Serpentes)

Caenophidian snakes include the file snake genus Acrochordus and advanced colubroidean snakes that radiated mainly during the Neogene. Although caenophidian snakes are a well-supported clade, their inferred affinities, based either on molecular or morphological data, remain poorly known or controversial. Here, we provide an expanded molecular phylogenetic analysis of Caenophidia and use three non-parametric measures of support-Shimodaira-Hasegawa-Like test (SHL), Felsentein (FBP) and transfer (TBE) bootstrap measures-to evaluate the robustness of each clade in the molecular tree. That very different alternative support values are common suggests that results based on only one support value should be viewed with caution. Using a scheme to combine support values, we find 20.9% of the 1265 clades comprising the inferred caenophidian tree are unambiguously supported by both SHL and FBP values, while almost 37% are unsupported or ambiguously supported, revealing the substantial extent of phylogenetic problems within Caenophidia. Combined FBP/TBE support values show similar results, while SHL/TBE result in slightly higher combined values. We consider key morphological attributes of colubroidean cranial, vertebral and hemipenial anatomy and provide additional morphological evidence supporting the clades Colubroides, Colubriformes, and Endoglyptodonta. We review and revise the relevant caenophidian fossil record and provide a time-calibrated tree derived from our molecular data to discuss the main cladogenetic events that resulted in present-day patterns of caenophidian diversification. Our results suggest that all extant families of Colubroidea and Elapoidea composing the present-day endoglyptodont fauna originated rapidly within the early Oligocene-between approximately 33 and 28 Mya-following the major terrestrial faunal turnover known as the "Grande Coupure" and associated with the overall climate shift at the Eocene-Oligocene boundary. Our results further suggest that the caenophidian radiation originated within the Caenozoic, with the divergence between Colubroides and Acrochordidae occurring in the early Eocene, at ~ 56 Mya.

Biodiversity seen through the perspective of insects: 10 simple rules on methodological choices and experimental design for genomic studies

Massively parallel DNA sequencing opens up opportunities for bridging multiple temporal and spatial dimensions in biodiversity research, thanks to its efficiency to recover millions of nucleotide polymorphisms. Here, we identify the current status, discuss the main challenges, and look into future perspectives on biodiversity genomics focusing on insects, which arguably constitute the most diverse and ecologically important group among all animals. We suggest 10 simple rules that provide a succinct step-by-step guide and best-practices to anyone interested in biodiversity research through the study of insect genomics. To this end, we review relevant literature on biodiversity and evolutionary research in the field of entomology. Our compilation is targeted at researchers and students who may not yet be specialists in entomology or molecular biology. We foresee that the genomic revolution and its application to the study of non-model insect lineages will represent a major leap to our understanding of insect diversity.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Mitochondrial genome of Phalantus geniculatus (Hemiptera: Reduviidae): trnT duplication and phylogenetic implications

Reduviidae is the second largest family of Heteroptera and most of them are important natural enemies of agricultural and forest pests. Most of the sequenced mitochondrial (mt) genomes in this family have the typical gene arrangement of insects and encode 37 coding genes (13 protein-coding genes, 22 tRNA genes and two rRNA genes). In the present study, we sequenced the mt genome of Phalantus geniculatus from the subfamily Peiratinae through high-throughput sequencing and encountered the duplication of tRNA genes for the first time in this subfamily. We identified 23 tRNA genes, including 22 tRNAs commonly found in insect mt genomes and an extra trnT (trnT2), which has high sequence similarity (96.9%) to trnT1. The presence of a "pseudo-trnP" in the non-coding region between trnT1 and trnT2 supports the hypothesis that the presence of an extra trnT can be explained by the tandem duplication-random loss (TDRL) model. Phylogenetic results inferred from mt genome sequences supported a sister relationship between Phymatinae and the remaining sampled subfamilies, as well as a paraphyletic Reduviinae. The present study highlights the utility of mt genomes in the phylogenetic study of Reduviidae based on the large scale taxon sampling in the future.

Museum specimens of terrestrial vertebrates are sensitive indicators of environmental change in the Anthropocene

Natural history museums and the specimen collections they curate are vital scientific infrastructure, a fact as true today as it was when biologists began collecting and preserving specimens over 200 years ago. The importance of museum specimens in studies of taxonomy, systematics, ecology and evolutionary biology is evidenced by a rich and abundant literature, yet creative and novel uses of specimens are constantly broadening the impact of natural history collections on biodiversity science and global sustainability. Excellent examples of the critical importance of specimens come from their use in documenting the consequences of environmental change, which is particularly relevant considering the alarming rate at which we now modify our planet in the Anthropocene. In this review, we highlight the important role of bird, mammal and amphibian specimens in documenting the Anthropocene and provide examples that underscore the need for continued collection of museum specimens.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Combining transcriptomes and ultraconserved elements to illuminate the phylogeny of Apidae

Two increasingly popular approaches to reconstruct the Tree of Life involve whole transcriptome sequencing and the target capture of ultraconserved elements (UCEs). Both methods can be used to generate large, multigene datasets for analysis of phylogenetic relationships in non-model organisms. While targeted exon sequencing across divergent lineages is now a standard method, it is still not clear if UCE data can be readily combined with published transcriptomes. In this study, we evaluate the combination of UCEs and transcriptomes in a single analysis using genome-, transcriptome-, and UCE data for 79 bees in the largest and most biologically diverse bee family, Apidae. Using existing tools, we first developed a workflow to assemble phylogenomic data from different sources and produced two large nucleotide matrices of combined data. We then reconstructed the phylogeny of the Apidae using concatenation- and coalescent-based methods, and critically evaluated the resulting phylogenies in the context of previously published genetic, genomic, and morphological data sets. Our estimated phylogenetic trees are robustly supported and largely congruent with previous molecular hypotheses, from deep nodes to shallow species-level phylogenies. Moreover, the combined approach allows us to resolve controversial nodes of the apid Tree of Life, by clarifying the relationships among the genera of orchid bees (Euglossini) and the monophyly of the Centridini. Additionally, we present novel phylogenetic evidence supporting the monophyly of the diverse clade of cleptoparasitic Apidae and the placement of two enigmatic, oil-collecting genera (Ctenoplectra and Tetrapedia). Lastly, we propose a revised classification of the family Apidae that reflects our improved understanding of apid higher-level relationships.

Population genomic data reveal extreme geographic subdivision and novel conservation actions for the declining foothill yellow-legged frog

Genomic data have the potential to inform high resolution landscape genetic and biological conservation studies that go far beyond recent mitochondrial and microsatellite analyses. We characterize the relationships of populations of the foothill yellow-legged frog, Rana boylii, a declining, "sentinel" species for stream ecosystems throughout its range in California and Oregon. We generated RADseq data and applied phylogenetic methods, hierarchical Bayesian clustering, PCA and population differentiation with admixture analyses to characterize spatial genetic structure across the species range. To facilitate direct comparison with previous analyses, we included many localities and individuals from our earlier work based on mitochondrial DNA. The results are striking, and emphasize the power of our landscape genomic approach. We recovered five extremely differentiated primary clades that indicate that R. boylii may be the most genetically differentiated anuran yet studied. Our results provide better resolution and more spatially consistent patterns than our earlier work, confirming the increased resolving power of genomic data compared to single-locus studies. Genomic structure is not equal across the species distribution. Approximately half the range of R. boylii consists of a single, relatively uniform population, while Sierra Nevada and coastal California clades are deeply, hierarchically substructured with biogeographic breaks observed in other codistributed taxa. Our results indicate that clades should serve as management units for R. boylii rather than previously suggested watershed boundaries, and that the near-extinct population from southwestern California is particularly diverged, exhibits the lowest genetic diversity, and is a critical conservation target for species recovery.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Sequence capture phylogenomics of eyeless Cicurina spiders from Texas caves, with emphasis on US federally-endangered species from Bexar County (Araneae, Hahniidae)

Morphological, mitochondrial, and nuclear phylogenomic data were combined to address phylogenetic and species delimitation questions in cave-limited Cicurina spiders from central Texas. Special effort was focused on specimens and cave locations in the San Antonio region (Bexar County), home to four eyeless species listed as US Federally Endangered. Sequence capture experiments resulted in the recovery of ~200–400 homologous ultra-conserved element (UCE) nuclear loci across taxa, and nearly complete COI mitochondrial DNA sequences from the same set of individuals. Some of these nuclear and mitochondrial sequences were recovered from “standard” museum specimens without special preservation of DNA material, including museum specimens preserved in the 1990s. Multiple phylogenetic analyses of the UCE data agree in the recovery of two major lineages of eyeless Cicurina in Texas. These lineages also differ in mitochondrial clade membership, female genitalic morphology, degree of troglomorphy (as measured by relative leg length), and are mostly allopatric across much of Texas. Rare sympatry was confirmed in Bexar County, where members of the two major clades sometimes co-exist in the same karst feature. Both nuclear phylogenomic and mitochondrial data indicate the existence of undescribed species from the San Antonio region, although further sampling and collection of adult specimens is needed to explicitly test these hypotheses. Our data support the two following species synonymies (Cicurinavenii Gertsch, 1992 = Cicurinamadla Gertsch, 1992; Cicurinaloftini Cokendolpher, 2004 = Cicurinavespera Gertsch, 1992), formally proposed here. Overall, our taxonomy-focused research has many important conservation implications, and again highlights the fundamental importance of robust taxonomy in conservation research.

HiMAP: Robust phylogenomics from highly multiplexed amplicon sequencing

High-throughput sequencing has fundamentally changed how molecular phylogenetic data sets are assembled, and phylogenomic data sets commonly contain 50- to 100-fold more loci than those generated using traditional Sanger sequencing-based approaches. Here, we demonstrate a new approach for building phylogenomic data sets using single-tube, highly multiplexed amplicon sequencing, which we name HiMAP (highly multiplexed amplicon-based phylogenomics) and present bioinformatic pipelines for locus selection based on genomic and transcriptomic data resources and postsequencing consensus calling and alignment. This method is inexpensive and amenable to sequencing a large number (hundreds) of taxa simultaneously and requires minimal hands-on time at the bench (<1/2 day), and data analysis can be accomplished without the need for read mapping or assembly. We demonstrate this approach by sequencing 878 amplicons in single reactions for 82 species of tephritid fruit flies across seven genera (384 individuals), including some of the most economically important agricultural insect pests. The resulting filtered data set (>150,000-bp concatenated alignment, ~20% missing character sites across all individuals and amplicons) contained >40,000 phylogenetically informative characters, and although some discordance was observed between analyses, it provided unparalleled resolution of many phylogenetic relationships in this group. Most notably, we found high support for the generic status of Zeugodacus and the sister relationship between Dacus and Zeugodacus. We discuss HiMAP, with regard to its molecular and bioinformatic strengths, and the insight the resulting data set provides into relationships of this diverse insect group.

HiMAP: Robust phylogenomics from highly multiplexed amplicon sequencing

High-throughput sequencing has fundamentally changed how molecular phylogenetic data sets are assembled, and phylogenomic data sets commonly contain 50- to 100-fold more loci than those generated using traditional Sanger sequencing-based approaches. Here, we demonstrate a new approach for building phylogenomic data sets using single-tube, highly multiplexed amplicon sequencing, which we name HiMAP (highly multiplexed amplicon-based phylogenomics) and present bioinformatic pipelines for locus selection based on genomic and transcriptomic data resources and postsequencing consensus calling and alignment. This method is inexpensive and amenable to sequencing a large number (hundreds) of taxa simultaneously and requires minimal hands-on time at the bench (<1/2 day), and data analysis can be accomplished without the need for read mapping or assembly. We demonstrate this approach by sequencing 878 amplicons in single reactions for 82 species of tephritid fruit flies across seven genera (384 individuals), including some of the most economically important agricultural insect pests. The resulting filtered data set (>150,000-bp concatenated alignment, ~20% missing character sites across all individuals and amplicons) contained >40,000 phylogenetically informative characters, and although some discordance was observed between analyses, it provided unparalleled resolution of many phylogenetic relationships in this group. Most notably, we found high support for the generic status of Zeugodacus and the sister relationship between Dacus and Zeugodacus. We discuss HiMAP, with regard to its molecular and bioinformatic strengths, and the insight the resulting data set provides into relationships of this diverse insect group.

Squeezing water from a stone: high-throughput sequencing from a 145-year old holotype resolves (barely) a cryptic species problem in flying lizards

We used Massively Parallel High-Throughput Sequencing to obtain genetic data from a 145-year old holotype specimen of the flying lizard, Draco cristatellus. Obtaining genetic data from this holotype was necessary to resolve an otherwise intractable taxonomic problem involving the status of this species relative to closely related sympatric Draco species that cannot otherwise be distinguished from one another on the basis of museum specimens. Initial analyses suggested that the DNA present in the holotype sample was so degraded as to be unusable for sequencing. However, we used a specialized extraction procedure developed for highly degraded ancient DNA samples and MiSeq shotgun sequencing to obtain just enough low-coverage mitochondrial DNA (721 base pairs) to conclusively resolve the species status of the holotype as well as a second known specimen of this species. The holotype was prepared before the advent of formalin-fixation and therefore was most likely originally fixed with ethanol and never exposed to formalin. Whereas conventional wisdom suggests that formalin-fixed samples should be the most challenging for DNA sequencing, we propose that evaporation during long-term alcohol storage and consequent water-exposure may subject older ethanol-fixed museum specimens to hydrolytic damage. If so, this may pose an even greater challenge for sequencing efforts involving historical samples.

Phylogenomic reclassification of the world's most venomous spiders (Mygalomorphae, Atracinae), with implications for venom evolution

Here we show that the most venomous spiders in the world are phylogenetically misplaced. Australian atracine spiders (family Hexathelidae), including the notorious Sydney funnel-web spider Atrax robustus, produce venom peptides that can kill people. Intriguingly, eastern Australian mouse spiders (family Actinopodidae) are also medically dangerous, possessing venom peptides strikingly similar to Atrax hexatoxins. Based on the standing morphology-based classification, mouse spiders are hypothesized distant relatives of atracines, having diverged over 200 million years ago. Using sequence-capture phylogenomics, we instead show convincingly that hexathelids are non-monophyletic, and that atracines are sister to actinopodids. Three new mygalomorph lineages are elevated to the family level, and a revised circumscription of Hexathelidae is presented. Re-writing this phylogenetic story has major implications for how we study venom evolution in these spiders, and potentially genuine consequences for antivenom development and bite treatment research. More generally, our research provides a textbook example of the applied importance of modern phylogenomic research.

Museum specimens indicate genetic erosion in an endangered lizard

Genetic variability, one of the main factors that guarantees species persistence, and species’ conservation status are generally evaluated with indices calculated at the present time. Natural history collections might help compare historical and current genetic diversity so to identify major trends. Here we analysed museum specimens of the lizard Zootoca vivipara carniolica, with a specific and stringent protocol for degraded DNA, in order to contrast its past and current genetic variability, using fragments of one mitochondrial DNA gene. Part of the distributional range of Z. v. carniolica (Po Plain, Italy), heavily impacted by human activities, was investigated. We found two previously unknown haplotypes in populations that are extinct today, suggesting the loss of these haplotypes and thus an overall shrinking of genetic variability. We argue that these results, together with the increasing threats posed by climate and land use changes, suggest that specific conservation measures for the persistence of Z. v. carniolica in Northern Italian lowlands have to be considered.

Return of a giant: DNA from archival museum samples helps to identify a unique cutthroat trout lineage formerly thought to be extinct

Currently one small, native population of the culturally and ecologically important Lahontan cutthroat trout (Oncorhynchus clarkii henshawi, LCT, Federally listed) remains in the Truckee River watershed of northwestern Nevada and northeastern California. The majority of populations in this watershed were extirpated in the 1940s due to invasive species, overharvest, anthropogenic water consumption and changing precipitation regimes. In 1977, a population of cutthroat trout discovered in the Pilot Peak Mountains in the Bonneville basin of Utah, was putatively identified as the extirpated LCT lacustrine lineage native to Pyramid Lake in the Truckee River basin based on morphological and meristic characters. Our phylogenetic and Bayesian genotype clustering analyses of museum specimens collected from the large lakes (1872-1913) and contemporary samples collected from populations throughout the extant range provide evidence in support of a genetically distinct Truckee River basin origin for this population. Analysis of museum samples alone identified three distinct genotype clusters and historical connectivity among water bodies within the Truckee River basin. Baseline data from museum collections indicate that the extant Pilot Peak strain represents a remnant of the extirpated lacustrine lineage. Given the limitations on high-quality data when working with a sparse number of preserved museum samples, we acknowledge that, in the end, this may be a more complicated story. However, the paucity of remnant populations in the Truckee River watershed, in combination with data on the distribution of morphological, meristic and genetic data for Lahontan cutthroat trout, suggests that recovery strategies, particularly in the large lacustrine habitats should consider this lineage as an important part of the genetic legacy of this species.

Novel Substrates as Sources of Ancient DNA: Prospects and Hurdles

Following the discovery in the late 1980s that hard tissues such as bones and teeth preserve genetic information, the field of ancient DNA analysis has typically concentrated upon these substrates. The onset of high-throughput sequencing, combined with optimized DNA recovery methods, has enabled the analysis of a myriad of ancient species and specimens worldwide, dating back to the Middle Pleistocene. Despite the growing sophistication of analytical techniques, the genetic analysis of substrates other than bone and dentine remain comparatively "novel". Here, we review analyses of other biological substrates which offer great potential for elucidating phylogenetic relationships, paleoenvironments, and microbial ecosystems including (1) archaeological artifacts and ecofacts; (2) calcified and/or mineralized biological deposits; and (3) biological and cultural archives. We conclude that there is a pressing need for more refined models of DNA preservation and bespoke tools for DNA extraction and analysis to authenticate and maximize the utility of the data obtained. With such tools in place the potential for neglected or underexploited substrates to provide a unique insight into phylogenetics, microbial evolution and evolutionary processes will be realized.