Mega-phylogeny approach for comparative biology: an alternative to supertree and supermatrix approaches

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

The use of Anchored Hybrid Enrichment data to resolve higher-level phylogenetic relationships: A proof-of-concept applied to Asterales (Eudicotyledoneae; Angiosperms)

Anchored Hybrid Enrichment (AHE) is a tool for capturing orthologous regions of the nuclear genome shared in low or single copy across lineages. Despite the increasing number of studies using this method, its usefulness to estimate relationships at deeper taxonomic levels in plants has not been fully explored. Here we present a proof of concept about the performance of nuclear loci obtained with AHE to infer phylogenetic relationships and explore the use of gene sampling schemes to estimate divergence times in Asterales. We recovered low-copy nuclear loci using the AHE method from herbarium material and silica-preserved samples. Maximum likelihood, Bayesian inference, and coalescence approaches were used to reconstruct phylogenomic relationships. Dating analyses were conducted under a multispecies coalescent approach by jointly inferring species tree and divergence times with random gene sampling schemes and multiple calibrations. We recovered 403 low-copy nuclear loci for 63 species representing nine out of eleven families of Asterales. Phylogenetic hypotheses were congruent among the applied methods and previously published results. Analyses with concatenated datasets were strongly supported, but coalescence-based analyses showed low support for the phylogenetic position of families Argophyllaceae and Alseuosmiaceae. Estimated family ages were congruent among gene sampling schemes, with the mean age for Asterales around 130 Myr. Our study documents the usefulness of AHE for resolving phylogenetic relationships at deep phylogenetic levels in Asterales. Observed phylogenetic inconsistencies were possibly due to the non-inclusion of families Phellinceae and Pentaphragmataceae. Random gene sampling schemes produced consistent age estimates with coalescence and species tree relaxed clock approaches.

An Automated Bioinformatics Pipeline Informing Near-Real-Time Public Health Responses to New HIV Diagnoses in a Statewide HIV Epidemic

Molecular HIV cluster data can guide public health responses towards ending the HIV epidemic. Currently, real-time data integration, analysis, and interpretation are challenging, leading to a delayed public health response. We present a comprehensive methodology for addressing these challenges through data integration, analysis, and reporting. We integrated heterogeneous data sources across systems and developed an open-source, automatic bioinformatics pipeline that provides molecular HIV cluster data to inform public health responses to new statewide HIV-1 diagnoses, overcoming data management, computational, and analytical challenges. We demonstrate implementation of this pipeline in a statewide HIV epidemic and use it to compare the impact of specific phylogenetic and distance-only methods and datasets on molecular HIV cluster analyses. The pipeline was applied to 18 monthly datasets generated between January 2020 and June 2022 in Rhode Island, USA, that provide statewide molecular HIV data to support routine public health case management by a multi-disciplinary team. The resulting cluster analyses and near-real-time reporting guided public health actions in 37 phylogenetically clustered cases out of 57 new HIV-1 diagnoses. Of the 37, only 21 (57%) clustered by distance-only methods. Through a unique academic-public health partnership, an automated open-source pipeline was developed and applied to prospective, routine analysis of statewide molecular HIV data in near-real-time. This collaboration informed public health actions to optimize disruption of HIV transmission.

Recent progress on methods for estimating and updating large phylogenies

With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

V.PhyloMaker2: An updated and enlarged R package that can generate very large phylogenies for vascular plants

An earlier version of V.PhyloMaker has been broadly used to generate phylogenetic trees of vascular plants for botanical, biogeographical and ecological studies. Here, we update and enlarge this package, which is now called 'V.PhyloMaker2'. With V.PhyloMaker2, one can generate a phylogenetic tree for vascular plants based on one of three different botanical nomenclature systems. V.PhyloMaker2 can generate phylogenies for very large species lists (the largest species list that we tested included 365,198 species). V.PhyloMaker2 generates phylogenies at a fast speed. We provide an example (including a sample species list and an R script to run it) in this paper to show how to use V.PhyloMaker2 to generate phylogenetic trees.

V.PhyloMaker2: An updated and enlarged R package that can generate very large phylogenies for vascular plants

An earlier version of V.PhyloMaker has been broadly used to generate phylogenetic trees of vascular plants for botanical, biogeographical and ecological studies. Here, we update and enlarge this package, which is now called 'V.PhyloMaker2'. With V.PhyloMaker2, one can generate a phylogenetic tree for vascular plants based on one of three different botanical nomenclature systems. V.PhyloMaker2 can generate phylogenies for very large species lists (the largest species list that we tested included 365,198 species). V.PhyloMaker2 generates phylogenies at a fast speed. We provide an example (including a sample species list and an R script to run it) in this paper to show how to use V.PhyloMaker2 to generate phylogenetic trees.

OpenTree: A Python Package for Accessing and Analyzing Data from the Open Tree of Life

The Open Tree of Life project constructs a comprehensive, dynamic, and digitally available tree of life by synthesizing published phylogenetic trees along with taxonomic data. Open Tree of Life provides web-service application programming interfaces (APIs) to make the tree estimate, unified taxonomy, and input phylogenetic data available to anyone. Here, we describe the Python package opentree, which provides a user friendly Python wrapper for these APIs and a set of scripts and tutorials for straightforward downstream data analyses. We demonstrate the utility of these tools by generating an estimate of the phylogenetic relationships of all bird families, and by capturing a phylogenetic estimate for all taxa observed at the University of California Merced Vernal Pools and Grassland Reserve.[Evolution; open science; phylogenetics; Python; taxonomy.].

Combining molecular and geographical data to infer the phylogeny of Lamiales and its dispersal patterns in and out of the tropics

Lamiales is one of the most intractable orders of flowering plants, with several changes in family composition, and circumscription throughout history. The order is worldwide distributed, occurring in tropical forests and frozen habitats. In this study, a comprehensive phylogeny of Lamiales was reconstructed using DNA sequences. The tree was used to infer dispersal patterns, focusing on the tropics and extratropics. Molecular and species geographic data available from public repositories were combined to address both objectives. A total of 6,910 species, and 842 genera of Lamiales were sampled using the Python tool PyPHLAWD. The tree was inferred using RAxML, and recovered a monophyletic Lamiales. All 26 families were recovered as monophyletic with high support. The families Bignoniaceae, and Plantaginaceae are remarkable examples. The first emerged as monophyletic and included tribe Jacarandeae, while the later emerged as monophyletic in its sensu lato and included both the tribes Angelonieae, and Gratioleae. Distribution points for all species were retrieved from GBIF. After filtering, 1,136,425 records were retained. Species were coded as present in extratropical or tropical environments. The in and out of the tropics dispersal patterns were inferred using a maximum likelihood approach that identifies hidden rate changes. The model recovered higher rates of transition from extratropics to tropics, estimating two rates of state transitions. When ancestral states are considered, more discrete transitions from extratropics to tropics were observed. The extratropical state was also inferred for the crown node of Lamiales and old nested nodes, revealing a rare pattern of transitions to the tropics throughout the upper Cretaceous and Tertiary. A significant phylogenetic signal was recovered for the in and out of the tropics dispersal patterns, showing that state transitions are not frequent enough to erase the effect of tree structure on the data.

Physcraper: a Python package for continually updated phylogenetic trees using the Open Tree of Life

BACKGROUND

Phylogenies are a key part of research in many areas of biology. Tools that automate some parts of the process of phylogenetic reconstruction, mainly molecular character matrix assembly, have been developed for the advantage of both specialists in the field of phylogenetics and non-specialists. However, interpretation of results, comparison with previously available phylogenetic hypotheses, and selection of one phylogeny for downstream analyses and discussion still impose difficulties to one that is not a specialist either on phylogenetic methods or on a particular group of study.

RESULTS

Physcraper is a command-line Python program that automates the update of published phylogenies by adding public DNA sequences to underlying alignments of previously published phylogenies. It also provides a framework for straightforward comparison of published phylogenies with their updated versions, by leveraging upon tools from the Open Tree of Life project to link taxonomic information across databases. The program can be used by the nonspecialist, as a tool to generate phylogenetic hypotheses based on publicly available expert phylogenetic knowledge. Phylogeneticists and taxonomic group specialists will find it useful as a tool to facilitate molecular dataset gathering and comparison of alternative phylogenetic hypotheses (topologies).

CONCLUSION

The Physcraper workflow showcases the benefits of doing open science for phylogenetics, encouraging researchers to strive for better scientific sharing practices. Physcraper can be used with any OS and is released under an open-source license. Detailed instructions for installation and usage are available at https://physcraper.readthedocs.io.

Biodiversity at the global scale: the synthesis continues

Traditionally, the generation and use of biodiversity data and their associated specimen objects have been primarily the purview of individuals and small research groups. While deposition of data and specimens in herbaria and other repositories has long been the norm, throughout most of their history, these resources have been accessible only to a small community of specialists. Through recent concerted efforts, primarily at the level of national and international governmental agencies over the last two decades, the pace of biodiversity data accumulation has accelerated, and a wider array of biodiversity scientists has gained access to this massive accumulation of resources, applying them to an ever-widening compass of research pursuits. We review how these new resources and increasing access to them are affecting the landscape of biodiversity research in plants today, focusing on new applications across evolution, ecology, and other fields that have been enabled specifically by the availability of these data and the global scope that was previously beyond the reach of individual investigators. We give an overview of recent advances organized along three lines: broad-scale analyses of distributional data and spatial information, phylogenetic research circumscribing large clades with comprehensive taxon sampling, and data sets derived from improved accessibility of biodiversity literature. We also review synergies between large data resources and more traditional data collection paradigms, describe shortfalls and how to overcome them, and reflect on the future of plant biodiversity analyses in light of increasing linkages between data types and scientists in our field.

A broadly resolved molecular phylogeny of New Zealand cheilostome bryozoans as a framework for hypotheses of morphological evolution

Larger molecular phylogenies based on ever more genes are becoming commonplace with the advent of cheaper and more streamlined sequencing and bioinformatics pipelines. However, many groups of inconspicuous but no less evolutionarily or ecologically important marine invertebrates are still neglected in the quest for understanding species- and higher-level phylogenetic relationships. Here, we alleviate this issue by presenting the molecular sequences of 165 cheilostome bryozoan species from New Zealand waters. New Zealand is our geographic region of choice as its cheilostome fauna is taxonomically, functionally and ecologically diverse, and better characterized than many other such faunas in the world. Using this most taxonomically broadly-sampled and statistically-supported cheilostome phylogeny comprising 214 species, when including previously published sequences, and 17 genes (2 nuclear and 15 mitochondrial) we tested several existing systematic hypotheses based solely on morphological observations. We find that lower taxonomic level hypotheses (species and genera) are robust while our inferred trees did not reflect current higher-level systematics (family and above), illustrating a general need for the rethinking of current hypotheses. To illustrate the utility of our new phylogeny, we reconstruct the evolutionary history of frontal shields (i.e., a calcified body-wall layer in ascus-bearing cheilostomes) and ask if its presence has any bearing on the diversification rates of cheilostomes.

mtDNAcombine: tools to combine sequences from multiple studies

BACKGROUND

Today an unprecedented amount of genetic sequence data is stored in publicly available repositories. For decades now, mitochondrial DNA (mtDNA) has been the workhorse of genetic studies, and as a result, there is a large volume of mtDNA data available in these repositories for a wide range of species. Indeed, whilst whole genome sequencing is an exciting prospect for the future, for most non-model organisms' classical markers such as mtDNA remain widely used. By compiling existing data from multiple original studies, it is possible to build powerful new datasets capable of exploring many questions in ecology, evolution and conservation biology. One key question that these data can help inform is what happened in a species' demographic past. However, compiling data in this manner is not trivial, there are many complexities associated with data extraction, data quality and data handling.

RESULTS

Here we present the mtDNAcombine package, a collection of tools developed to manage some of the major decisions associated with handling multi-study sequence data with a particular focus on preparing sequence data for Bayesian skyline plot demographic reconstructions.

CONCLUSIONS

There is now more genetic information available than ever before and large meta-data sets offer great opportunities to explore new and exciting avenues of research. However, compiling multi-study datasets still remains a technically challenging prospect. The mtDNAcombine package provides a pipeline to streamline the process of downloading, curating, and analysing sequence data, guiding the process of compiling data sets from the online database GenBank.

Total Ortholog Median Matrix as an alternative unsupervised approach for phylogenomics based on evolutionary distance between protein coding genes

The increasing number of available genomic data allowed the development of phylogenomic analytical tools. Current methods compile information from single gene phylogenies, whether based on topologies or multiple sequence alignments. Generally, phylogenomic analyses elect gene families or genomic regions to construct phylogenomic trees. Here, we presented an alternative approach for Phylogenomics, named TOMM (Total Ortholog Median Matrix), to construct a representative phylogram composed by amino acid distance measures of all pairwise ortholog protein sequence pairs from desired species inside a group of organisms. The procedure is divided two main steps, (1) ortholog detection and (2) creation of a matrix with the median amino acid distance measures of all pairwise orthologous sequences. We tested this approach within three different group of organisms: Kinetoplastida protozoa, hematophagous Diptera vectors and Primates. Our approach was robust and efficacious to reconstruct the phylogenetic relationships for the three groups. Moreover, novel branch topologies could be achieved, providing insights about some phylogenetic relationships between some taxa.

A Total-Group Phylogenetic Metatree for Cetacea and the Importance of Fossil Data in Diversification Analyses

Phylogenetic trees provide a powerful framework for testing macroevolutionary hypotheses, but it is becoming increasingly apparent that inferences derived from extant species alone can be highly misleading. Trees incorporating living and extinct taxa are are needed to address fundamental questions about the origins of diversity and disparity but it has proved challenging to generate robust, species-rich phylogenies that include large numbers of fossil taxa. As a result, most studies of diversification dynamics continue to rely on molecular phylogenies. Here, we extend and apply a recently developed meta-analytic approach for synthesizing previously published phylogenetic studies to infer a well-resolved set of species level, time-scaled phylogenetic hypotheses for extinct and extant cetaceans (whales, dolphins and allies). Our trees extend sampling from the ∼ 90 extant species to over 500 living and extinct species, and therefore allow for more robust inference of macroevolutionary dynamics. While the diversification scenarios we recover are broadly concordant with those inferred from molecular phylogenies they differ in critical ways, notably in the relative contributions of extinction and speciation rate shifts in driving rapid radiations. The metatree approach provides the most immediate route for generating higher level phylogenies of extinct taxa, and opens the door to re-evaluation of macroevolutionary hypotheses derived only from extant taxa.

Machine learning predicts large scale declines in native plant phylogenetic diversity

Though substantial effort has gone into predicting how global climate change will impact biodiversity patterns, the scarcity of taxon-specific information has hampered the efficacy of these endeavors. Further, most studies analyzing spatiotemporal patterns of biodiversity focus narrowly on species richness. We apply machine learning approaches to a comprehensive vascular plant database for the United States and generate predictive models of regional plant taxonomic and phylogenetic diversity in response to a wide range of environmental variables. We demonstrate differences in predicted patterns and potential drivers of native vs nonnative biodiversity. In particular, native phylogenetic diversity is likely to decrease over the next half century despite increases in species richness. We also identify that patterns of taxonomic diversity can be incongruent with those of phylogenetic diversity. The combination of macro-environmental factors that determine diversity likely varies at continental scales; thus, as climate change alters the combinations of these factors across the landscape, the collective effect on regional diversity will also vary. Our study represents one of the most comprehensive examinations of plant diversity patterns to date and demonstrates that our ability to predict future diversity may benefit tremendously from the application of machine learning.

The macroevolutionary dynamics of symbiotic and phenotypic diversification in lichens

Symbioses are evolutionarily pervasive and play fundamental roles in structuring ecosystems, yet our understanding of their macroevolutionary origins, persistence, and consequences is incomplete. We traced the macroevolutionary history of symbiotic and phenotypic diversification in an iconic symbiosis, lichens. By inferring the most comprehensive time-scaled phylogeny of lichen-forming fungi (LFF) to date (over 3,300 species), we identified shifts among symbiont classes that broadly coincided with the convergent evolution of phylogenetically or functionally similar associations in diverse lineages (plants, fungi, bacteria). While a relatively recent loss of lichenization in Lecanoromycetes was previously identified, our work instead suggests lichenization was abandoned far earlier, interrupting what had previously been considered a direct switch between trebouxiophycean and trentepohlialean algal symbionts. Consequently, some of the most diverse clades of LFF are instead derived from nonlichenized ancestors and re-evolved lichenization with Trentepohliales algae, a clade that also facilitated lichenization in unrelated lineages of LFF. Furthermore, while symbiont identity and symbiotic phenotype influence the ecology and physiology of lichens, they are not correlated with rates of lineage birth and death, suggesting more complex dynamics underly lichen diversification. Finally, diversification patterns of LFF differed from those of wood-rotting and ectomycorrhizal taxa, likely reflecting contrasts in their fundamental biological properties. Together, our work provides a timeline for the ecological contributions of lichens, and reshapes our understanding of symbiotic persistence in a classic model of symbiosis.

Recent accelerated diversification in rosids occurred outside the tropics

Conflicting relationships have been found between diversification rate and temperature across disparate clades of life. Here, we use a supermatrix comprising nearly 20,000 species of rosids-a clade of ~25% of all angiosperm species-to understand global patterns of diversification and its climatic association. Our approach incorporates historical global temperature, assessment of species' temperature niche, and two broad-scale characterizations of tropical versus non-tropical niche occupancy. We find the diversification rates of most subclades dramatically increased over the last 15 million years (Myr) during cooling associated with global expansion of temperate habitats. Climatic niche is negatively associated with diversification rates, with tropical rosids forming older communities and experiencing speciation rates ~2-fold below rosids in cooler climates. Our results suggest long-term cooling had a disproportionate effect on non-tropical diversification rates, leading to dynamic young communities outside of the tropics, while relative stability in tropical climes led to older, slower-evolving but still species-rich communities.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution

PREMISE

Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution.

METHODS

We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements.

RESULTS

Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses.

CONCLUSIONS

A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

Geophytism in monocots leads to higher rates of diversification

Geophytes, plants with buds on underground structures, are found throughout the plant tree of life. These below ground structures allow plants to inhabit highly seasonal and disturbance-prone environments across ecosystems. Past researchers have hypothesised that the bulbous, cormous and tuberous habits promote diversification, but this had yet to be tested. Using a comprehensive monocot data set of almost 13 000 taxa, we investigated the effects of the geophytic habit on diversification using both state-dependent and state-independent models. We found that geophytes exhibit increased rates of diversification relative to nongeophytes. State-dependent analyses recovered higher yet similar rates of diversification for bulbous, cormous and tuberous taxa compared with rhizomatous and nongeophytic taxa. However, the state-independent model returned no difference in rates among the different traits. Geophytism shows higher rates of diversification relative to nongeophytes but we found little support for the hypothesis that the evolution of the bulb, corm or tuber appears to provide a diversification increase relative to rhizomatous and nongeophytic taxa. Our broad-scale analysis highlights the overall evolutionary importance of the geophytic habit (i.e. belowground bud placement). However, our results also suggest that belowground morphological diversity alone cannot explain this rate increase. In order to further test the evolutionary significance of these underground structures, future studies should consider these in combination with other biotic and abiotic factors.

Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation

Big, time-scaled phylogenies are fundamental to connecting evolutionary processes to modern biodiversity patterns. Yet inferring reliable phylogenetic trees for thousands of species involves numerous trade-offs that have limited their utility to comparative biologists. To establish a robust evolutionary timescale for all approximately 6,000 living species of mammals, we developed credible sets of trees that capture root-to-tip uncertainty in topology and divergence times. Our "backbone-and-patch" approach to tree building applies a newly assembled 31-gene supermatrix to two levels of Bayesian inference: (1) backbone relationships and ages among major lineages, using fossil node or tip dating, and (2) species-level "patch" phylogenies with nonoverlapping in-groups that each correspond to one representative lineage in the backbone. Species unsampled for DNA are either excluded ("DNA-only" trees) or imputed within taxonomic constraints using branch lengths drawn from local birth-death models ("completed" trees). Joining time-scaled patches to backbones results in species-level trees of extant Mammalia with all branches estimated under the same modeling framework, thereby facilitating rate comparisons among lineages as disparate as marsupials and placentals. We compare our phylogenetic trees to previous estimates of mammal-wide phylogeny and divergence times, finding that (1) node ages are broadly concordant among studies, and (2) recent (tip-level) rates of speciation are estimated more accurately in our study than in previous "supertree" approaches, in which unresolved nodes led to branch-length artifacts. Credible sets of mammalian phylogenetic history are now available for download at http://vertlife.org/phylosubsets, enabling investigations of long-standing questions in comparative biology.

Leaf longevity in temperate evergreen species is related to phylogeny and leaf size

Leaf longevity (LL), the amount of time a photosynthetically active leaf remains on a plant, is an important trait of evergreen species, affecting physiological ecology and ecosystem processes. A long LL gives leaves more time to fix carbon but carries higher construction costs, while a short LL allows plants to respond more rapidly to changing environmental conditions. For many evergreen taxa, LL data are not readily available, and it is not known if LL is phylogenetically conserved. To address this gap, we measured LL for 169 temperate and boreal evergreen woody species at the Arnold Arboretum, a botanical garden in Boston, Massachusetts, along with metrics of leaf size and number known to be related to LL. We hypothesized that LL is phylogenetically conserved, and that longer LL is associated with a greater numbers of leaves, smaller leaves, and a colder hardiness zone of the species' native range. We found that average LL ranged from 1.4 years in Rhododendron tomentosum to 10.5 years in Abies cilicia. LL was phylogenetically conserved, with some genera, such as Abies and Picea, exhibiting long LL (> 3 years) and others, such as Ilex and Rhododendron, exhibiting short LL (< 3 years). Leaf length was negatively correlated with LL in conifers, due to differences between Pinus and other genera; however, there was no correlation between LL and number of leaves. This study highlights the considerable variation and phylogenetic pattern in LL among temperate evergreen species, which has implications for carbon budgets and ecosystem models.

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.

For common community phylogenetic analyses, go ahead and use synthesis phylogenies

Should we build our own phylogenetic trees based on gene sequence data, or can we simply use available synthesis phylogenies? This is a fundamental question that any study involving a phylogenetic framework must face at the beginning of the project. Building a phylogeny from gene sequence data (purpose-built phylogeny) requires more effort, expertise, and cost than subsetting an already available phylogeny (synthesis-based phylogeny). However, we still lack a comparison of how these two approaches to building phylogenetic trees influence common community phylogenetic analyses such as comparing community phylogenetic diversity and estimating trait phylogenetic signal. Here, we generated three purpose-built phylogenies and their corresponding synthesis-based trees (two from Phylomatic and one from the Open Tree of Life, OTL). We simulated 1,000 communities and 12,000 continuous traits along each purpose-built phylogeny. We then compared the effects of different trees on estimates of phylogenetic diversity (alpha and beta) and phylogenetic signal (Pagel's λ and Blomberg's K). Synthesis-based phylogenies generally yielded higher estimates of phylogenetic diversity when compared to purpose-built phylogenies. However, resulting measures of phylogenetic diversity from both types of phylogenies were highly correlated (Spearman's ρ  > 0.8 in most cases). Mean pairwise distance (both alpha and beta) is the index that is most robust to the differences in tree construction that we tested. Measures of phylogenetic diversity based on the OTL showed the highest correlation with measures based on the purpose-built phylogenies. Trait phylogenetic signal estimated with synthesis-based phylogenies, especially from the OTL, was also highly correlated with estimates of Blomberg's K or close to Pagel's λ from purpose-built phylogenies when traits were simulated under Brownian motion. For commonly employed community phylogenetic analyses, our results justify taking advantage of recently developed and continuously improving synthesis trees, especially the Open Tree of Life.

Increasing phylogenetic stochasticity at high elevations on summits across a remote North American wilderness

PREMISE

At the intersection of ecology and evolutionary biology, community phylogenetics can provide insights into overarching biodiversity patterns, particularly in remote and understudied ecosystems. To understand community assembly of the high alpine flora in the Sawtooth National Forest, USA, we analyzed phylogenetic structure within and between nine summit communities.

METHODS

We used high-throughput sequencing to supplement existing data and infer a nearly completely sampled community phylogeny of the alpine vascular flora. We calculated mean nearest taxon distance (MNTD) and mean pairwise distance (MPD) to quantify phylogenetic divergence within summits, and assessed whether maximum elevation explains phylogenetic structure. To evaluate similarities between summits, we quantified phylogenetic turnover, taking into consideration microhabitats (talus vs. meadows).

RESULTS

We found different patterns of community phylogenetic structure within the six most species-rich orders, but across all vascular plants phylogenetic structure was largely not different from random. There was a significant negative correlation between elevation and tree-wide phylogenetic diversity (MPD) within summits: overdispersion degraded as elevation increased. Between summits, we found high phylogenetic turnover driven by greater niche heterogeneity on summits with alpine meadows.

CONCLUSIONS

Our results provide further evidence that stochastic processes may also play an important role in the assembly of vascular plant communities in high alpine habitats at regional scales. However, order-specific patterns suggest that adaptations are still important for assembly of specific sectors of the plant tree of life. Further studies quantifying functional diversity will be important in disentangling the interplay of eco-evolutionary processes that likely shape broad community phylogenetic patterns in extreme environments.

European mushroom assemblages are darker in cold climates

Thermal melanism theory states that dark-colored ectotherm organisms are at an advantage at low temperature due to increased warming. This theory is generally supported for ectotherm animals, however, the function of colors in the fungal kingdom is largely unknown. Here, we test whether the color lightness of mushroom assemblages is related to climate using a dataset of 3.2 million observations of 3,054 species across Europe. Consistent with the thermal melanism theory, mushroom assemblages are significantly darker in areas with cold climates. We further show differences in color phenotype between fungal lifestyles and a lifestyle differentiated response to seasonality. These results indicate a more complex ecological role of mushroom colors and suggest functions beyond thermal adaption. Because fungi play a crucial role in terrestrial carbon and nutrient cycles, understanding the links between the thermal environment, functional coloration and species' geographical distributions will be critical in predicting ecosystem responses to global warming.

The monocotyledonous underground: global climatic and phylogenetic patterns of geophyte diversity

PREMISE

Geophytes-plants that typically possess a bulb, corm, tuber, and/or rhizome-have long captured the attention of hobbyists and researchers. However, despite the economic and evolutionary importance of these traits, the potential drivers of their morphological diversity remain unknown. Using a comprehensive phylogeny of monocots, we test for correlations between climate and geophyte growth form to better understand why we observe such a diversity of underground traits in geophytes. Understanding the evolutionary factors promoting independent origins of these potentially adaptive organs will lend insights into how plants adapt to environmental hardships.

METHODS

Using a comprehensive phylogeny incorporated with global occurrence and climate data for the monocots, we investigated whether climatic patterns could explain differences between geophytes and non-geophytes, as well as differences among bulbous, cormous, tuberous, rhizomatous, and non-geophytic taxa. We used phylogenetically-informed ANOVAs, MANOVAs, and PCAs to test differences in climatic variables between the different growth forms.

RESULTS

Geophytes inhabit cooler, drier, and more thermally variable climates compared to non-geophytes. Although some underground traits (i.e., bulb, corm, and tuber) appear to inhabit particular niches, a result supported by strong phylogenetic signal, our data has limited evidence for an overall role of climate in the evolution of these traits. However, temperature may be a driving force in rhizome evolution, as well as the evolution of taxa which we considered here as non-geophytic (e.g., trees, epiphytes, etc.).

CONCLUSIONS

While precipitation patterns have played a role in the evolution of geophytism, our results suggest that temperature should be more strongly considered as a contributing factor promoting the evolution of belowground bud placement, specifically in rhizomatous and non-geophytic taxa. Bulbous, cormous, and tuberous taxa need closer examination of other mechanisms, such as anatomical constraints or genetic controls, in order to begin to understand the causes behind the evolution of their underground morphology.

Rates of niche and phenotype evolution lag behind diversification in a temperate radiation

Alternative models of evolutionary processes suggest different associations between species diversification and trait evolution, but limited empirical evidence is available to test these models across large clades at global extents. Here we investigate the relative timing of species diversification and niche and phenotypic evolution across a global plant radiation (Saxifragales) with enormous phenotypic and habitat variation. We demonstrate strong temporal lags among rates, with increased diversification occurring first, followed by niche and phenotype. Accelerated diversification rates are coincident with mid-Miocene expansion of temperate biomes. Later increases in niche and phenotypic evolutionary rates argue against density-dependent diversification alone, indicating a major role for ecological opportunity. These results have broad implications for understanding diversification processes and the origin of present-day temperate biotas.

Environmental change can create opportunities for increased rates of lineage diversification, but continued species accumulation has been hypothesized to lead to slowdowns via competitive exclusion and niche partitioning. Such density-dependent models imply tight linkages between diversification and trait evolution, but there are plausible alternative models. Little is known about the association between diversification and key ecological and phenotypic traits at broad phylogenetic and spatial scales. Do trait evolutionary rates coincide with rates of diversification, are there lags among these rates, or is diversification niche-neutral? To address these questions, we combine a deeply sampled phylogeny for a major flowering plant clade—Saxifragales—with phenotype and niche data to examine temporal patterns of evolutionary rates. The considerable phenotypic and habitat diversity of Saxifragales is greatest in temperate biomes. Global expansion of these habitats since the mid-Miocene provided ecological opportunities that, with density-dependent adaptive radiation, should result in simultaneous rate increases for diversification, niche, and phenotype, followed by decreases with habitat saturation. Instead, we find that these rates have significantly different timings, with increases in diversification occurring at the mid-Miocene Climatic Optimum (∼15 Mya), followed by increases in niche and phenotypic evolutionary rates by ∼5 Mya; all rates increase exponentially to the present. We attribute this surprising lack of temporal coincidence to initial niche-neutral diversification followed by ecological and phenotypic divergence coincident with more extreme cold and dry habitats that proliferated into the Pleistocene. A lack of density-dependence contrasts with investigations of other cosmopolitan lineages, suggesting alternative patterns may be common in the diversification of temperate lineages.

Temperature shapes opposing latitudinal gradients of plant taxonomic and phylogenetic β diversity

Latitudinal and elevational richness gradients have received much attention from ecologists but there is little consensus on underlying causes. One possible proximate cause is increased levels of species turnover, or β diversity, in the tropics compared to temperate regions. Here, we leverage a large botanical dataset to map taxonomic and phylogenetic β diversity, as mean turnover between neighboring 100 × 100 km cells, across the Americas and determine key climatic drivers. We find taxonomic and tip-weighted phylogenetic β diversity is higher in the tropics, but that basal-weighted phylogenetic β diversity is highest in temperate regions. Supporting Janzen's 'mountain passes' hypothesis, tropical mountainous regions had higher β diversity than temperate regions for taxonomic and tip-weighted metrics. The strongest climatic predictors of turnover were average temperature and temperature seasonality. Taken together, these results suggest β diversity is coupled to latitudinal richness gradients and that temperature is a major driver of plant community composition and change.

Gecko Adhesion in Space and Time: A Phylogenetic Perspective on the Scansorial Success Story

An evolutionary perspective on gecko adhesion was previously hampered by a lack of an explicit phylogeny for the group and of robust comparative methods to study trait evolution, an underappreciation for the taxonomic and structural diversity of geckos, and a dearth of fossil evidence bearing directly on the origin of the scansorial apparatus. With a multigene dataset as the basis for a comprehensive gekkotan phylogeny, model-based methods have recently been employed to estimate the number of unique derivations of the adhesive system and its role in lineage diversification. Evidence points to a single basal origin of the spinulate oberhautchen layer of the epidermis, which is a necessary precursor for the subsequent elaboration of a functional adhesive mechanism in geckos. However, multiple gains and losses are implicated for the elaborated setae that are necessary for adhesion via van der Waals forces. The well-supported phylogeny of gekkotans has demonstrated that convergence and parallelism in digital design are even more prevalent than previously believed. It also permits the reexamination of previously collected morphological data in an explicitly evolutionary context. Both time-calibrated trees and recently discovered amber fossils that preserve gecko toepads suggest that a fully-functional adhesive apparatus was not only present, but also represented by diverse architectures, by the mid-Cretaceous. Further characterization and phylogenetically-informed analyses of the other components of the adhesive system (muscles, tendons, blood sinuses, etc.) will permit a more comprehensive reconstruction of the evolutionary pathway(s) by which geckos have achieved their structural and taxonomic diversity. A phylogenetic perspective can meaningfully inform functional and performance studies of gecko adhesion and locomotion and can contribute to advances in bioinspired materials.

Darwin review: angiosperm phylogeny and evolutionary radiations

Darwin's dual interests in evolution and plants formed the basis of evolutionary botany, a field that developed following his publications on both topics. Here, we review his many contributions to plant biology—from the evolutionary origins of angiosperms to plant reproduction, carnivory, and movement—and note that he expected one day there would be a ‘true’ genealogical tree for plants. This view fuelled the field of plant phylogenetics. With perhaps nearly 400 000 species, the angiosperms have diversified rapidly since their origin in the Early Cretaceous, often through what appear to be rapid radiations. We describe these evolutionary patterns, evaluate possible drivers of radiations, consider how new approaches to studies of diversification can contribute to our understanding of angiosperm diversity, and suggest new directions for further insight into plant evolution.

Phylogenetic imprint of woody plants on the soil mycobiome in natural mountain forests of eastern China

Recent studies have detected strong phylogenetic signals in tree-fungus associations for diseased leaves and mycorrhizal symbioses. However, the extent of plant phylogenetic constraints on the free-living soil mycobiome remains unknown, especially at broad geographic scales. Here, 343 soil samples were collected adjacent to individual tree trunks, representing 58 woody plant species located in five mountain forests of eastern China. Integrating plant species identity and phylogenetic information, we aimed to unravel the relative contributions of phylogenetic relationships among tree species, abiotic environmental filtering, and geographic isolation to the geographic distribution of soil mycobiome. We found that the community dissimilarities of total fungi and each dominant guild (viz. saprotrophs, plant pathogens, and ectomycorrhizal fungi) significantly increased with increasing plant phylogenetic distance. Plant phylogenetic eigenvectors explained 11.4% of the variation in community composition, whereas environmental and spatial factors explained 24.1% and 7.2% of the variation, respectively. The communities of ectomycorrhizal fungi and plant pathogens were relatively more strongly affected by plant phylogeny than those of saprotrophs (13.7% and 10.4% vs. 8.5%). Overall, our results demonstrate how plant phylogeny, environment, and geographic space contribute to forest soil fungal distributions and suggest that the influence of plant phylogeny on fungal association may differ by guilds.

An overview of extant conifer evolution from the perspective of the fossil record

PREMISE OF THE STUDY

Conifers are an important living seed plant lineage with an extensive fossil record spanning more than 300 million years. The group therefore provides an excellent opportunity to explore congruence and conflict between dated molecular phylogenies and the fossil record.

METHODS

We surveyed the current state of knowledge in conifer phylogenetics to present a new time-calibrated molecular tree that samples ~90% of extant species diversity. We compared phylogenetic relationships and estimated divergence ages in this new phylogeny with the paleobotanical record, focusing on clades that are species-rich and well known from fossils.

KEY RESULTS

Molecular topologies and estimated divergence ages largely agree with the fossil record in Cupressaceae, conflict with it in Araucariaceae, and are ambiguous in Pinaceae and Podocarpaceae. Molecular phylogenies provide insights into some fundamental questions in conifer evolution, such as the origin of their seed cones, but using them to reconstruct the evolutionary history of specific traits can be challenging.

CONCLUSIONS

Molecular phylogenies are useful for answering deep questions in conifer evolution if they depend on understanding relationships among extant lineages. Because of extinction, however, molecular datasets poorly sample diversity from periods much earlier than the Late Cretaceous. This fundamentally limits their utility for understanding deep patterns of character evolution and resolving the overall pattern of conifer phylogeny.

OneTwoTree: An online tool for phylogeny reconstruction

Phylogeny reconstruction is a key instrument in numerous biological analyses, ranging from evolutionary and ecology research, to conservation and systems biology. The increasing accumulation of genomic data makes it possible to reconstruct phylogenies with both high accuracy and at increasingly finer resolution. Yet, taking advantage of the enormous amount of sequence data available requires the use of computational tools for efficient data retrieval and processing, or else the process could quickly become an error-prone endeavour. Here, we present OneTwoTree (http://onetwotree.tau.ac.il/), a Web-based tool for tree reconstruction based on the supermatrix paradigm. Given a list of taxa names of interest as the sole input requirement, OneTwoTree retrieves all available sequence data from NCBI GenBank, clusters these into orthology groups, identifies the most informative set of markers, searches for an appropriate outgroup, and assembles a partitioned sequence matrix that is then used for the final phylogeny reconstruction step. OneTwoTree further allows users to control various steps of the process, such as the merging of sequences from similar clusters, or phylogeny reconstruction based on markers from a specific genome type. By comparing the performance of OneTwoTree to a manually reconstructed phylogeny of the Antirrhineae tribe, we show that the use of OneTwoTree resulted in substantially higher data coverage in terms of both taxon sampling and the number of informative markers assembled. OneTwoTree provides a flexible online tool for species-tree reconstruction, aimed to assist researchers ranging in their level of prior expertise in the task of phylogeny reconstruction.

Evolutionary dynamics of host specialization in wood-decay fungi

Background

The majority of wood decomposing fungi are mushroom-forming Agaricomycetes, which exhibit two main modes of plant cell wall decomposition: white rot, in which all plant cell wall components are degraded, including lignin, and brown rot, in which lignin is modified but not appreciably removed. Previous studies suggested that brown rot fungi tend to be specialists of gymnosperm hosts and that brown rot promotes gymnosperm specialization. However, these hypotheses were based on analyses of limited datasets of Agaricomycetes. Overcoming this limitation, we used a phylogeny with 1157 species integrating available sequences, assembled decay mode characters from the literature, and coded host specialization using the newly developed R package, rusda.

Results

We found that most brown rot fungi are generalists or gymnosperm specialists, whereas most white rot fungi are angiosperm specialists. A six-state model of the evolution of host specialization revealed high transition rates between generalism and specialization in both decay modes. However, while white rot lineages switched most frequently to angiosperm specialists, brown rot lineages switched most frequently to generalism. A time-calibrated phylogeny revealed that Agaricomycetes is older than the flowering plants but many of the large clades originated after the diversification of the angiosperms in the Cretaceous.

Conclusions

Our results challenge the current view that brown rot fungi are primarily gymnosperm specialists and reveal intensive white rot specialization to angiosperm hosts. We thus suggest that brown rot associated convergent loss of lignocellulose degrading enzymes was correlated with host generalism, rather than gymnosperm specialism. A likelihood model of host specialization evolution together with a time-calibrated phylogeny further suggests that the rise of the angiosperms opened a new mega-niche for wood-decay fungi, which was exploited particularly well by white rot lineages.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1229-7) contains supplementary material, which is available to authorized users.

Predicting Extinction Risk for Data Deficient Bats

Conservation biology aims to identify species most at risk of extinction and to understand factors that forecast species vulnerability. The International Union for Conservation of Nature (IUCN) Red List is a leading source for extinction risk data of species globally, however, many potentially at risk species are not assessed by the IUCN owing to inadequate data. Of the approximately 1150 bat species (Chiroptera) recognized by the IUCN, 17 percent are categorized as Data Deficient. Here, we show that large trait databases in combination with a comprehensive phylogeny can identify which traits are important for assessing extinction risk in bats. Using phylogenetic logistic regressions, we show that geographic range and island endemism are the strongest correlates of binary extinction risk. We also show that simulations using two models that trade-off between data complexity and data coverage provide similar estimates of extinction risk for species that have received a Red List assessment. We then use our model parameters to provide quantitative predictions of extinction risk for 60 species that have not received risk assessments by the IUCN. Our model suggests that at least 20 bat species should be treated as threatened by extinction. In combination with expert knowledge, our results can be used as a quick, first-pass prioritization for conservation action.

An inverse latitudinal gradient in speciation rate for marine fishes

Far more species of organisms are found in the tropics than in temperate and polar regions, but the evolutionary and ecological causes of this pattern remain controversial^1,2. Tropical marine fish communities are much more diverse than cold-water fish communities found at higher latitudes^3,4, and several explanations for this latitudinal diversity gradient propose that warm reef environments serve as evolutionary 'hotspots' for species formation^5-8. Here we test the relationship between latitude, species richness and speciation rate across marine fishes. We assembled a time-calibrated phylogeny of all ray-finned fishes (31,526 tips, of which 11,638 had genetic data) and used this framework to describe the spatial dynamics of speciation in the marine realm. We show that the fastest rates of speciation occur in species-poor regions outside the tropics, and that high-latitude fish lineages form new species at much faster rates than their tropical counterparts. High rates of speciation occur in geographical regions that are characterized by low surface temperatures and high endemism. Our results reject a broad class of mechanisms under which the tropics serve as an evolutionary cradle for marine fish diversity and raise new questions about why the coldest oceans on Earth are present-day hotspots of species formation.

AnnotationBustR: an R package to extract subsequences from GenBank annotations

BACKGROUND

DNA sequences are pivotal for a wide array of research in biology. Large sequence databases, like GenBank, provide an amazing resource to utilize DNA sequences for large scale analyses. However, many sequence records on GenBank contain more than one gene or are portions of genomes. Inconsistencies in the way genes are annotated and the numerous synonyms a single gene may be listed under provide major challenges for extracting large numbers of subsequences for comparative analysis across taxa. At present, there is no easy way to extract portions from many GenBank accessions based on annotations where gene names may vary extensively.

RESULTS

The R package AnnotationBustR allows users to extract sequences based on GenBank annotations through the ACNUC retrieval system given search terms of gene synonyms and accession numbers. AnnotationBustR extracts subsequences of interest and then writes them to a FASTA file for users to employ in their research endeavors.

CONCLUSION

FASTA files of extracted subsequences and accession tables generated by AnnotationBustR allow users to quickly find and extract subsequences from GenBank accessions. These sequences can then be incorporated in various analyses, like the construction of phylogenies to test a wide range of ecological and evolutionary hypotheses.

ViCTree: an automated framework for taxonomic classification from protein sequences

Motivation

The increasing rate of submission of genetic sequences into public databases is providing a growing resource for classifying the organisms that these sequences represent. To aid viral classification, we have developed ViCTree, which automatically integrates the relevant sets of sequences in NCBI GenBank and transforms them into an interactive maximum likelihood phylogenetic tree that can be updated automatically. ViCTree incorporates ViCTreeView, which is a JavaScript-based visualization tool that enables the tree to be explored interactively in the context of pairwise distance data.

Results

To demonstrate utility, ViCTree was applied to subfamily Densovirinae of family Parvoviridae. This led to the identification of six new species of insect virus.

Availability and implementation

ViCTree is open-source and can be run on any Linux- or Unix-based computer or cluster. A tutorial, the documentation and the source code are available under a GPL3 license, and can be accessed at http://bioinformatics.cvr.ac.uk/victree_web/.

Supplementary information

Supplementary data are available at Bioinformatics online.