American Museum of Natural History

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies. Results for the final phase of the 1000 Genomes Project are presented including whole-genome sequencing, targeted exome sequencing, and genotyping on high-density SNP arrays for 2,504 individuals across 26 populations, providing a global reference data set to support biomedical genetics. The 1000 Genomes Project has sought to comprehensively catalogue human genetic variation across populations, providing a valuable public genomic resource. The data obtained so far have found applications ranging from association studies and fine mapping studies to the filtering of likely neutral variants in rare-disease cohorts. The authors now report on the final phase of the project, phase 3, which covers previously uncharacterized areas of human genetic diversity in terms of the populations sampled and categories of characterized variation. The sample now includes more than 2,500 individuals from 26 global populations, with low coverage whole-genome and deep exome sequencing, as well as dense microarray genotyping. They find that while most common variants are shared across populations, rarer variants are often restricted to closely related populations. The authors also demonstrate the use of the phase 3 dataset as a reference panel for imputation to improve the resolution in genetic association studies.

Freshwater biodiversity is the over-riding conservation priority during the International Decade for Action - 'Water for Life' - 2005 to 2015. Fresh water makes up only 0.01% of the World's water and approximately 0.8% of the Earth's surface, yet this tiny fraction of global water supports at least 100000 species out of approximately 1.8 million - almost 6% of all described species. Inland waters and freshwater biodiversity constitute a valuable natural resource, in economic, cultural, aesthetic, scientific and educational terms. Their conservation and management are critical to the interests of all humans, nations and governments. Yet this precious heritage is in crisis. Fresh waters are experiencing declines in biodiversity far greater than those in the most affected terrestrial ecosystems, and if trends in human demands for water remain unaltered and species losses continue at current rates, the opportunity to conserve much of the remaining biodiversity in fresh water will vanish before the 'Water for Life' decade ends in 2015. Why is this so, and what is being done about it? This article explores the special features of freshwater habitats and the biodiversity they support that makes them especially vulnerable to human activities. We document threats to global freshwater biodiversity under five headings: overexploitation; water pollution; flow modification; destruction or degradation of habitat; and invasion by exotic species. Their combined and interacting influences have resulted in population declines and range reduction of freshwater biodiversity worldwide. Conservation of biodiversity is complicated by the landscape position of rivers and wetlands as 'receivers' of land-use effluents, and the problems posed by endemism and thus non-substitutability. In addition, in many parts of the world, fresh water is subject to severe competition among multiple human stakeholders. Protection of freshwater biodiversity is perhaps the ultimate conservation challenge because it is influenced by the upstream drainage network, the surrounding land, the riparian zone, and - in the case of migrating aquatic fauna - downstream reaches. Such prerequisites are hardly ever met. Immediate action is needed where opportunities exist to set aside intact lake and river ecosystems within large protected areas. For most of the global land surface, trade-offs between conservation of freshwater biodiversity and human use of ecosystem goods and services are necessary. We advocate continuing attempts to check species loss but, in many situations, urge adoption of a compromise position of management for biodiversity conservation, ecosystem functioning and resilience, and human livelihoods in order to provide a viable long-term basis for freshwater conservation. Recognition of this need will require adoption of a new paradigm for biodiversity protection and freshwater ecosystem management - one that has been appropriately termed 'reconciliation ecology'.

By characterizing the geographic and functional spectrum of human genetic variation, the 1000 Genomes Project aims to build a resource to help to understand the genetic contribution to disease. Here we describe the genomes of 1,092 individuals from 14 populations, constructed using a combination of low-coverage whole-genome and exome sequencing. By developing methods to integrate information across several algorithms and diverse data sources, we provide a validated haplotype map of 38 million single nucleotide polymorphisms, 1.4 million short insertions and deletions, and more than 14,000 larger deletions. We show that individuals from different populations carry different profiles of rare and common variants, and that low-frequency variants show substantial geographic differentiation, which is further increased by the action of purifying selection. We show that evolutionary conservation and coding consequence are key determinants of the strength of purifying selection, that rare-variant load varies substantially across biological pathways, and that each individual contains hundreds of rare non-coding variants at conserved sites, such as motif-disrupting changes in transcription-factor-binding sites. This resource, which captures up to 98% of accessible single nucleotide polymorphisms at a frequency of 1% in related populations, enables analysis of common and low-frequency variants in individuals from diverse, including admixed, populations. This report from the 1000 Genomes Project describes the genomes of 1,092 individuals from 14 human populations, providing a resource for common and low-frequency variant analysis in individuals from diverse populations; hundreds of rare non-coding variants at conserved sites, such as motif-disrupting changes in transcription-factor-binding sites, can be found in each individual. This report by the 1000 Genomes Project describes the genomes of 1,092 individuals from 14 human populations, providing a resource for common and low-frequency variant analysis in individuals from diverse populations. Integrative analyses reveal profiles of rare and common variants in different populations. The frequencies of rare variants vary across biological pathways, and hundreds of rare, non-coding variants at conserved sites — such as changes disrupting transcription-factor motifs — can be established for each individual.

Many studies have shown how pigments and internal nanostructures generate color in nature. External surface structures can also influence appearance, such as by causing multiple scattering of light (structural absorption) to produce a velvety, super black appearance. Here we show that feathers from five species of birds of paradise (Aves: Paradisaeidae) structurally absorb incident light to produce extremely low-reflectance, super black plumages. Directional reflectance of these feathers (0.05-0.31%) approaches that of man-made ultra-absorbent materials. SEM, nano-CT, and ray-tracing simulations show that super black feathers have titled arrays of highly modified barbules, which cause more multiple scattering, resulting in more structural absorption, than normal black feathers. Super black feathers have an extreme directional reflectance bias and appear darkest when viewed from the distal direction. We hypothesize that structurally absorbing, super black plumage evolved through sensory bias to enhance the perceived brilliance of adjacent color patches during courtship display.

Abstract The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy , which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

This paper describes the Seventh Data Release of the Sloan Digital Sky Survey (SDSS), marking the completion of the original goals of the SDSS and the end of the phase known as SDSS-II. It includes 11,663 deg2 of imaging data, with most of the ~2000 deg2 increment over the previous data release lying in regions of low Galactic latitude. The catalog contains five-band photometry for 357 million distinct objects. The survey also includes repeat photometry on a 120° long, 2fdg5 wide stripe along the celestial equator in the Southern Galactic Cap, with some regions covered by as many as 90 individual imaging runs. We include a co-addition of the best of these data, going roughly 2 mag fainter than the main survey over 250 deg2. The survey has completed spectroscopy over 9380 deg2; the spectroscopy is now complete over a large contiguous area of the Northern Galactic Cap, closing the gap that was present in previous data releases. There are over 1.6 million spectra in total, including 930,000 galaxies, 120,000 quasars, and 460,000 stars. The data release includes improved stellar photometry at low Galactic latitude. The astrometry has all been recalibrated with the second version of the USNO CCD Astrograph Catalog, reducing the rms statistical errors at the bright end to 45 milliarcseconds per coordinate. We further quantify a systematic error in bright galaxy photometry due to poor sky determination; this problem is less severe than previously reported for the majority of galaxies. Finally, we describe a series of improvements to the spectroscopic reductions, including better flat fielding and improved wavelength calibration at the blue end, better processing of objects with extremely strong narrow emission lines, and an improved determination of stellar metallicities.

Full list of authors: Price-Whelan, Adrian M.; Lim, Pey Lian; Earl, Nicholas; Starkman, Nathaniel; Bradley, Larry; Shupe, David L.; Patil, Aarya A.; Corrales, Lia; Brasseur, C. E.; Noethe, Maximilian; Donath, Axel; Tollerud, Erik; Morris, Brett M.; Ginsburg, Adam; Vaher, Eero; Weaver, Benjamin A.; Tocknell, James; Jamieson, William; van Kerkwijk, Marten H.; Robitaille, Thomas P.; Merry, Bruce; Bachetti, Matteo; Gunther, H. Moritz; Aldcroft, Thomas L.; Alvarado-Montes, Jaime A.; Archibald, Anne M.; Bodi, Attila; Bapat, Shreyas; Barentsen, Geert; Bazan, Juanjo; Biswas, Manish; Boquien, Mederic; Burke, D. J.; Cara, Daria; Cara, Mihai; Conroy, Kyle E.; Conseil, Simon; Craig, Matthew W.; Cross, Robert M.; Cruz, Kelle L.; D'Eugenio, Francesco; Dencheva, Nadia; Devillepoix, Hadrien A. R.; Dietrich, Jorg P.; Eigenbrot, Arthur Davis; Erben, Thomas; Ferreira, Leonardo; Foreman-Mackey, Daniel; Fox, Ryan; Freij, Nabil; Garg, Suyog; Geda, Robel; Glattly, Lauren; Gondhalekar, Yash; Gordon, Karl D.; Grant, David; Greenfield, Perry; Groener, Austen M.; Guest, Steve; Gurovich, Sebastian; Handberg, Rasmus; Hart, Akeem; Hatfield-Dodds, Zac; Homeier, Derek; Hosseinzadeh, Griffin; Jenness, Tim; Jones, Craig K.; Joseph, Prajwel; Kalmbach, J. Bryce; Karamehmetoglu, Emir; Kaluszynski, Mikolaj; Kelley, Michael S. P.; Kern, Nicholas; Kerzendorf, Wolfgang E.; Koch, Eric W.; Kulumani, Shankar; Lee, Antony; Ly, Chun; Ma, Zhiyuan; MacBride, Conor; Maljaars, Jakob M.; Muna, Demitri; Murphy, N. A.; Norman, Henrik; O'Steen, Richard; Oman, Kyle A.; Pacifici, Camilla; Pascual, Sergio; Pascual-Granado, J.; Patil, Rohit R.; Perren, Gabriel, I; Pickering, Timothy E.; Rastogi, Tanuj; Roulston, Benjamin R.; Ryan, Daniel F.; Rykoff, Eli S.; Sabater, Jose; Sakurikar, Parikshit; Salgado, Jesus; Sanghi, Aniket; Saunders, Nicholas; Savchenko, Volodymyr; Schwardt, Ludwig; Seifert-Eckert, Michael; Shih, Albert Y.; Jain, Anany Shrey; Shukla, Gyanendra; Sick, Jonathan; Simpson, Chris; Singanamalla, Sudheesh; Singer, Leo P.; Singhal, Jaladh; Sinha, Manodeep; Sipocz, Brigitta M.; Spitler, Lee R.; Stansby, David; Streicher, Ole; Sumak, Jani; Swinbank, John D.; Taranu, Dan S.; Tewary, Nikita; Tremblay, Grant R.; De Val-Borro, Miguel; Vasovic, Zlatan; Van Kooten, Samuel J.; Verma, Shresth; Cardoso, Jose Vinicius de Miranda; Williams, Peter K. G.; Wilson, Tom J.; Winkel, Benjamin; Wood-Vasey, W. M.; Xue, Rui; Yoachim, Peter; Zhang, Chen; Zonca, Andrea; Astropy Project Contributors; TARDIS Collaboration; Astropy Coordination Comm.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Peer Reviewed

Abstract Aim Techniques that predict species potential distributions by combining observed occurrence records with environmental variables show much potential for application across a range of biogeographical analyses. Some of the most promising applications relate to species for which occurrence records are scarce, due to cryptic habits, locally restricted distributions or low sampling effort. However, the minimum sample sizes required to yield useful predictions remain difficult to determine. Here we developed and tested a novel jackknife validation approach to assess the ability to predict species occurrence when fewer than 25 occurrence records are available. Location Madagascar. Methods Models were developed and evaluated for 13 species of secretive leaf‐tailed geckos ( Uroplatus spp.) that are endemic to Madagascar, for which available sample sizes range from 4 to 23 occurrence localities (at 1 km 2 grid resolution). Predictions were based on 20 environmental data layers and were generated using two modelling approaches: a method based on the principle of maximum entropy (Maxent) and a genetic algorithm (GARP). Results We found high success rates and statistical significance in jackknife tests with sample sizes as low as five when the Maxent model was applied. Results for GARP at very low sample sizes (less than c. 10) were less good. When sample sizes were experimentally reduced for those species with the most records, variability among predictions using different combinations of localities demonstrated that models were greatly influenced by exactly which observations were included. Main conclusions We emphasize that models developed using this approach with small sample sizes should be interpreted as identifying regions that have similar environmental conditions to where the species is known to occur, and not as predicting actual limits to the range of a species. The jackknife validation approach proposed here enables assessment of the predictive ability of models built using very small sample sizes, although use of this test with larger sample sizes may lead to overoptimistic estimates of predictive power. Our analyses demonstrate that geographical predictions developed from small numbers of occurrence records may be of great value, for example in targeting field surveys to accelerate the discovery of unknown populations and species.

This software note announces a new open‐source release of the Maxent software for modeling species distributions from occurrence records and environmental data, and describes a new R package for fitting such models. The new release (ver. 3.4.0) will be hosted online by the American Museum of Natural History, along with future versions. It contains small functional changes, most notably use of a complementary log‐log (cloglog) transform to produce an estimate of occurrence probability. The cloglog transform derives from the recently‐published interpretation of Maxent as an inhomogeneous Poisson process (IPP), giving it a stronger theoretical justification than the logistic transform which it replaces by default. In addition, the new R package, maxnet, fits Maxent models using the glmnet package for regularized generalized linear models. We discuss the implications of the IPP formulation in terms of model inputs and outputs, treating occurrence records as points rather than grid cells and interpreting the exponential Maxent model (raw output) as as an estimate of relative abundance. With these two open‐source developments, we invite others to freely use and contribute to the software.

We believe that punctuational change dominates the history of life: evolution is concentrated in very rapid events of speciation (geologically instantaneous, even if tolerably continuous in ecological time). Most species, during their geological history, either do not change in any appreciable way, or else they fluctuate mildly in morphology, with no apparent direction. Phyletic gradualism is very rare and too slow, in any case, to produce the major events of evolution. Evolutionary trends are not the product of slow, directional transformation within lineages; they represent the differential success of certain species within a clade—speciation may be random with respect to the direction of a trend (Wright's rule). As an a priori bias, phyletic gradualism has precluded any fair assessment of evolutionary tempos and modes. It could not be refuted by empirical catalogues constructed in its light because it excluded contrary information as the artificial result of an imperfect fossil record. With the model of punctuated equilibria, an unbiased distribution of evolutionary tempos can be established by treating stasis as data and by recording the pattern of change for all species in an assemblage. This distribution of tempos can lead to strong inferences about modes. If, as we predict, the punctuational tempo is prevalent, then speciation—not phyletic evolution—must be the dominant mode of evolution. We argue that virtually none of the examples brought forward to refute our model can stand as support for phyletic gradualism; many are so weak and ambiguous that they only reflect the persistent bias for gradualism still deeply embedded in paleontological thought. Of the few stronger cases, we concentrate on Gingerich's data for Hyopsodus and argue that it provides an excellent example of species selection under our model. We then review the data of several studies that have supported our model since we published it five years ago. The record of human evolution seems to provide a particularly good example: no gradualism has been detected within any hominid taxon, and many are long-ranging; the trend to larger brains arises from differential success of essentially static taxa. The data of molecular genetics support our assumption that large genetic changes often accompany the process of speciation. Phyletic gradualism was an a priori assertion from the start—it was never “seen” in the rocks; it expressed the cultural and political biases of 19th century liberalism. Huxley advised Darwin to eschew it as an “unnecessary difficulty.” We think that it has now become an empirical fallacy. A punctuational view of change may have wide validity at all levels of evolutionary processes. At the very least, it deserves consideration as an alternate way of interpreting the history of life.

ABSTRACT A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence–absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS‐INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM‐GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size ( n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

We focus on the implications of movement, landscape variables, and spatial heterogeneity for food web dynamics. Movements of nutrients, detritus, prey, and consumers among habitats are ubiquitous in diverse biomes and can strongly influence population, consumer-resource, food web, and community dynamics. Nutrient and detrital subsidies usually increase primary and secondary productivity, both directly and indirectly. Prey subsidies, by movement of either prey or predators, usually enhance predator abundance beyond what local resources can support. Top-down effects occur when spatially subsidized consumers affect local resources by suppressing key resources and occasionally by initiating trophic cascades. Effects on community dynamics vary with the relative amount of input, the trophic roles of the mobile and recipient entities, and the local food web structure. Landscape variables such as the perimeter/area ratio of the focal habitat, permeability of habitat boundaries, and relative productivity of trophically connected habitats affect the degree and importance of spatial subsidization.

Spatial thinning of species occurrence records can help address problems associated with spatial sampling biases. Ideally, thinning removes the fewest records necessary to substantially reduce the effects of sampling bias, while simultaneously retaining the greatest amount of useful information. Spatial thinning can be done manually; however, this is prohibitively time consuming for large datasets. Using a randomization approach, the ‘thin’ function in the spThin R package returns a dataset with the maximum number of records for a given thinning distance, when run for sufficient iterations. We here provide a worked example for the Caribbean spiny pocket mouse, where the results obtained match those of manual thinning.

Many mountain ranges have been strongly glaciated during the Quaternary ice ages, and the locations of glacial refugia of mountain plants have been debated for a long time. A series of detailed molecular studies, investigating intraspecific genetic variation of mountain plants in the European Alps, now allows for a first synopsis. A comparison of the phylogeographic patterns with geological and palaeoenvironmental data demonstrates that glacial refugia were located along the southwestern, southern, eastern and northern border of the Alps. Additional glacial refugia were present in central Alpine areas, where high-elevation plants survived the last glaciation on ice-free mountain tops. The observed intraspecific phylogeographies suggest general patterns of glacial survival, which conform to well-known centres of Alpine species diversity and endemism. This implies that evolutionary or biogeographic processes induced by climatic fluctuations act on gene and species diversity in a similar way.

Summary Recent studies have demonstrated a need for increased rigour in building and evaluating ecological niche models ( ENM s) based on presence‐only occurrence data. Two major goals are to balance goodness‐of‐fit with model complexity (e.g. by ‘tuning’ model settings) and to evaluate models with spatially independent data. These issues are especially critical for data sets suffering from sampling bias, and for studies that require transferring models across space or time (e.g. responses to climate change or spread of invasive species). Efficient implementation of procedures to accomplish these goals, however, requires automation. We developed ENM eval , an R package that: (i) creates data sets for k ‐fold cross‐validation using one of several methods for partitioning occurrence data (including options for spatially independent partitions), (ii) builds a series of candidate models using Maxent with a variety of user‐defined settings and (iii) provides multiple evaluation metrics to aid in selecting optimal model settings. The six methods for partitioning data are n −1 jackknife, random k ‐folds ( = bins), user‐specified folds and three methods of masked geographically structured folds. ENM eval quantifies six evaluation metrics: the area under the curve of the receiver‐operating characteristic plot for test localities ( AUC TEST ), the difference between training and testing AUC ( AUC DIFF ), two different threshold‐based omission rates for test localities and the Akaike information criterion corrected for small sample sizes ( AIC c). We demonstrate ENM eval by tuning model settings for eight tree species of the genus Coccoloba in Puerto Rico based on AIC c. Evaluation metrics varied substantially across model settings, and models selected with AIC c differed from default ones. In summary, ENMeval facilitates the production of better ENM s and should promote future methodological research on many outstanding issues.

To better determine the history of modern birds, we performed a genome-scale phylogenetic analysis of 48 species representing all orders of Neoaves using phylogenomic methods created to handle genome-scale data. We recovered a highly resolved tree that confirms previously controversial sister or close relationships. We identified the first divergence in Neoaves, two groups we named Passerea and Columbea, representing independent lineages of diverse and convergently evolved land and water bird species. Among Passerea, we infer the common ancestor of core landbirds to have been an apex predator and confirm independent gains of vocal learning. Among Columbea, we identify pigeons and flamingoes as belonging to sister clades. Even with whole genomes, some of the earliest branches in Neoaves proved challenging to resolve, which was best explained by massive protein-coding sequence convergence and high levels of incomplete lineage sorting that occurred during a rapid radiation after the Cretaceous-Paleogene mass extinction event about 66 million years ago.

The Cosmic Evolution Survey (COSMOS) is designed to probe the correlated evolution of galaxies, star formation, active galactic nuclei (AGNs), and dark matter (DM) with large-scale structure (LSS) over the redshift range z>0.5-6. The survey includes multiwavelength imaging and spectroscopy from X-ray-to-radio wavelengths covering a 2 deg2 area, including HST imaging. Given the very high sensitivity and resolution of these data sets, COSMOS also provides unprecedented samples of objects at high redshift with greatly reduced cosmic variance, compared to earlier surveys. Here we provide a brief overview of the survey strategy, the characteristics of the major COSMOS data sets, and a summary the science goals. Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555 also based on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan; the XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA; the European Southern Observatory under Large Program 175.A-0839, Chile; Kitt Peak National Observatory, Cerro Tololo Inter-American Observatory, and the National Optical Astronomy Observatory, which are operated by the AURA, Inc., under cooperative agreement with the National Science Foundation; the National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.; and the Canada-France-Hawaii Telescope with MegaPrime/MegaCam operated as a joint project by the CFHT Corporation, CEA/DAPNIA, the NRC and CADC of Canada, the CNRS of France, TERAPIX, and the University of Hawaii.

We present an in silico approach for the reconstruction of complete mitochondrial genomes of non-model organisms directly from next-generation sequencing (NGS) data-mitochondrial baiting and iterative mapping (MITObim). The method is straightforward even if only (i) distantly related mitochondrial genomes or (ii) mitochondrial barcode sequences are available as starting-reference sequences or seeds, respectively. We demonstrate the efficiency of the approach in case studies using real NGS data sets of the two monogenean ectoparasites species Gyrodactylus thymalli and Gyrodactylus derjavinoides including their respective teleost hosts European grayling (Thymallus thymallus) and Rainbow trout (Oncorhynchus mykiss). MITObim appeared superior to existing tools in terms of accuracy, runtime and memory requirements and fully automatically recovered mitochondrial genomes exceeding 99.5% accuracy from total genomic DNA derived NGS data sets in <24 h using a standard desktop computer. The approach overcomes the limitations of traditional strategies for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information at hand and represents a fast and highly efficient in silico alternative to laborious conventional strategies relying on initial long-range PCR. We furthermore demonstrate the applicability of MITObim for metagenomic/pooled data sets using simulated data. MITObim is an easy to use tool even for biologists with modest bioinformatics experience. The software is made available as open source pipeline under the MIT license at https://github.com/chrishah/MITObim.

Abstract The evidentiary basis of the currently accepted classification of living amphibians is discussed and shown not to warrant the degree of authority conferred on it by use and tradition. A new taxonomy of living amphibians is proposed to correct the deficiencies of the old one. This new taxonomy is based on the largest phylogenetic analysis of living Amphibia so far accomplished. We combined the comparative anatomical character evidence of Haas (2003) with DNA sequences from the mitochondrial transcription unit H1 (12S and 16S ribosomal RNA and tRNAValine genes, ≈ 2,400 bp of mitochondrial sequences) and the nuclear genes histone H3, rhodopsin, tyrosinase, and seven in absentia, and the large ribosomal subunit 28S (≈ 2,300 bp of nuclear sequences; ca. 1.8 million base pairs; x̄ = 3.7 kb/terminal). The dataset includes 532 terminals sampled from 522 species representative of the global diversity of amphibians as well as seven of the closest living relatives of amphibians for outgroup comparisons. The primary purpose of our taxon sampling strategy was to provide strong tests of the monophyly of all “family-group” taxa. All currently recognized nominal families and subfamilies were sampled, with the exception of Protohynobiinae (Hynobiidae). Many of the currently recognized genera were also sampled. Although we discuss the monophyly of genera, and provide remedies for nonmonophyly where possible, we also make recommendations for future research. A parsimony analysis was performed under Direct Optimization, which simultaneously optimizes nucleotide homology (alignment) and tree costs, using the same set of assumptions throughout the analysis. Multiple search algorithms were run in the program POY over a period of seven months of computing time on the AMNH Parallel Computing Cluster. Results demonstrate that the following major taxonomic groups, as currently recognized, are nonmonophyletic: Ichthyophiidae (paraphyletic with respect to Uraeotyphlidae), Caeciliidae (paraphyletic with respect to Typhlonectidae and Scolecomorphidae), Salamandroidea (paraphyletic with respect to Sirenidae), Leiopelmatanura (paraphyletic with respect to Ascaphidae), Discoglossanura (paraphyletic with respect to Bombinatoridae), Mesobatrachia (paraphyletic with respect to Neobatrachia), Pipanura (paraphyletic with respect to Bombinatoridae and Discoglossidae/Alytidae), Hyloidea (in the sense of containing Heleophrynidae; paraphyletic with respect to Ranoidea), Leptodactylidae (polyphyletic, with Batrachophrynidae forming the sister taxon of Myobatrachidae + Limnodynastidae, and broadly paraphyletic with respect to Hemiphractinae, Rhinodermatidae, Hylidae, Allophrynidae, Centrolenidae, Brachycephalidae, Dendrobatidae, and Bufonidae), Microhylidae (polyphyletic, with Brevicipitinae being the sister taxon of Hemisotidae), Microhylinae (poly/paraphyletic with respect to the remaining non-brevicipitine microhylids), Hyperoliidae (para/polyphyletic, with Leptopelinae forming the sister taxon of Arthroleptidae + Astylosternidae), Astylosternidae (paraphyletic with respect to Arthroleptinae), Ranidae (paraphyletic with respect to Rhacophoridae and Mantellidae). In addition, many subsidiary taxa are demonstrated to be nonmonophyletic, such as (1) Eleutherodactylus with respect to Brachycephalus; (2) Rana (sensu Dubois, 1992), which is polyphyletic, with various elements falling far from each other on the tree; and (3) Bufo, with respect to several nominal bufonid genera. A new taxonomy of living amphibians is proposed, and the evidence for this is presented to promote further investigation and data acquisition bearing on the evolutionary history of amphibians. The taxonomy provided is consistent with the International Code of Zoological Nomenclature (ICZN, 1999). Salient features of the new taxonomy are (1) the three major groups of living amphibians, caecilians/Gymnophiona, salamanders/Caudata, and frogs/Anura, form a monophyletic group, to which we restrict the name Amphibia; (2) Gymnophiona forms the sister taxon of Batrachia (salamanders + frogs) and is composed of two groups, Rhinatrematidae and Stegokrotaphia; (3) Stegokrotaphia is composed of two families, Ichthyophiidae (including Uraeotyphlidae) and Caeciliidae (including Scolecomorphidae and Typhlonectidae, which are regarded as subfamilies); (4) Batrachia is a highly corroborated monophyletic group, composed of two taxa, Caudata (salamanders) and Anura (frogs); (5) Caudata is composed of two taxa, Cryptobranchoidei (Cryptobranchidae and Hynobiidae) and Diadectosalamandroidei new taxon (all other salamanders); (6) Diadectosalamandroidei is composed of two taxa, Hydatinosalamandroidei new taxon (composed of Perennibranchia and Treptobranchia new taxon) and Plethosalamandroidei new taxon; (7) Perennibranchia is composed of Proteidae and Sirenidae; (8) Treptobranchia new taxon is composed of two taxa, Ambystomatidae (including Dicamptodontidae) and Salamandridae; (9) Plethosalamandroidei new taxon is composed of Rhyacotritonidae and Xenosalamandroidei new taxon; (10) Xenosalamandroidei is composed of Plethodontidae and Amphiumidae; (11) Anura is monophyletic and composed of two clades, Leiopelmatidae (including Ascaphidae) and Lalagobatrachia new taxon (all other frogs); (12) Lalagobatrachia is composed of two clades, Xenoanura (Pipidae and Rhinophrynidae) and Sokolanura new taxon (all other lalagobatrachians); (13) Bombinatoridae and Alytidae (former Discoglossidae) are each others' closest relatives and in a clade called Costata, which, excluding Leiopelmatidae and Xenoanura, forms the sister taxon of all other frogs, Acosmanura; (14) Acosmanura is composed of two clades, Anomocoela (= Pelobatoidea of other authors) and Neobatrachia; (15) Anomocoela contains Pelobatoidea (Pelobatidae and Megophryidae) and Pelodytoidea (Pelodytidae and Scaphiopodidae), and forms the sister taxon of Neobatrachia, together forming Acosmanura; (16) Neobatrachia is composed of two clades, Heleophrynidae, and all other neobatrachians, Phthanobatrachia new taxon; (17) Phthanobatrachia is composed of two major units, Hyloides and Ranoides; (18) Hyloides comprises Sooglossidae (including Nasikabatrachidae) and Notogaeanura new taxon (the remaining hyloids); (19) Notogaeanura contains two taxa, Australobatrachia new taxon and Nobleobatrachia new taxon; (20) Australobatrachia is a clade composed of Batrachophrynidae and its sister taxon, Myobatrachoidea (Myobatrachidae and Limnodynastidae), which forms the sister taxon of all other hyloids, excluding sooglossids; (21) Nobleobatrachia new taxon, is dominated at its base by frogs of a treefrog morphotype, several with intercalary phalangeal cartilages—Hemiphractus (Hemiphractidae) forms the sister taxon of the remaining members of this group, here termed Meridianura new taxon; (22) Meridianura comprises Brachycephalidae (former Eleutherodactylinae + Brachycephalus) and Cladophrynia new taxon; (23) Cladophrynia is composed of two groups, Cryptobatrachidae (composed of Cryptobatrachus and Stefania, previously a fragment of the polyphyletic Hemiphractinae) and Tinctanura new taxon; (24) Tinctanura is composed of Amphignathodontidae (Gastrotheca and Flectonotus, another fragment of the polyphyletic Hemiphractinae) and Athesphatanura new taxon; (25) Athesphatanura is composed of Hylidae (Hylinae, Pelodryadinae, and Phyllomedusinae, and excluding former Hemiphractinae, whose inclusion would have rendered this taxon polyphyletic) and Leptodactyliformes new taxon; (26) Leptodactyliformes is composed of Diphyabatrachia new taxon (composed of Centrolenidae [including Allophryne] and Leptodactylidae, sensu stricto, including Leptodactylus and relatives) and Chthonobatrachia new taxon; (27) Chthonobatrachia is composed of a reformulated Ceratophryidae (which excludes such genera as Odontophrynus and Proceratophrys and includes other taxa, such as Telmatobius) and Hesticobatrachia new taxon; (28) Hesticobatrachia is composed of a reformulated Cycloramphidae (which includes Rhinoderma) and Agastorophrynia new taxon; (29) Agastorophrynia is composed of Bufonidae (which is partially revised) and Dendrobatoidea (Dendrobatidae and Thoropidae); (30) Ranoides new taxon forms the sister taxon of Hyloides and is composed of two major monophyletic components, Allodapanura new taxon (microhylids, hyperoliids, and allies) and Natatanura new taxon (ranids and allies); (31) Allodapanura is composed of Microhylidae (which is partially revised) and Afrobatrachia new taxon; (32) Afrobatrachia is composed of Xenosyneunitanura new taxon (the “strange-bedfellows” Brevicipitidae [formerly in Microhylidae] and Hemisotidae) and a more normal-looking group of frogs, Laurentobatrachia new taxon (Hyperoliidae and Arthroleptidae, which includes Leptopelinae and former Astylosternidae); (33) Natatanura new taxon is composed of two taxa, the African Ptychadenidae and the worldwide Victoranura new taxon; (34) Victoranura is composed of Ceratobatrachidae and Telmatobatrachia new taxon; (35) Telmatobatrachia is composed of Micrixalidae and a worldwide group of ranoids, Ametrobatrachia new taxon; (36) Ametrobatrachia is composed of Africanura new taxon and Saukrobatrachia new taxon; (37) Africanura is composed of two taxa: Phrynobatrachidae (Phrynobatrachus, including Dimorphognathus and Phrynodon as synonyms) and Pyxicephaloidea; (38) Pyxicephaloidea is composed of Petropedetidae (Conraua, Indirana, Arthroleptides, and Petropedetes), and Pyxicephalidae (including a number of African genera, e.g. Amietia [including Afrana], Arthroleptella, Pyxicephalus, Strongylopus, and Tomopterna); and (39) Saukrobatrachia new taxon is the sister taxon of Africanura and is composed of Dicroglossidae and Aglaioanura new taxon, which is, in turn, composed of Rhacophoroidea (Mantellidae and Rhacophoridae) and Ranoidea (Nyctibatrachidae and Ranidae, sensu stricto). Many generic revisions are made either to render a monophyle