A real data-driven simulation strategy to select an imputation method for mixed-type trait data.

Missing observations in trait datasets pose an obstacle for analyses in myriad biological disciplines. Considering the mixed results of imputation, the wide variety of available methods, and the varied structure of real trait datasets, a framework for selecting a suitable imputation method is advantageous. We invoked a real data-driven simulation strategy to select an imputation method for a given mixed-type (categorical, count, continuous) target dataset. Candidate methods included mean/mode imputation, k-nearest neighbour, random forests, and multivariate imputation by chained equations (MICE). Using a trait dataset of squamates (lizards and amphisbaenians; order: Squamata) as a target dataset, a complete-case dataset consisting of species with nearly complete information was formed for the imputation method selection. Missing data were induced by removing values from this dataset under different missingness mechanisms: missing completely at random (MCAR), missing at random (MAR), and missing not at random (MNAR). For each method, combinations with and without phylogenetic information from single gene (nuclear and mitochondrial) or multigene trees were used to impute the missing values for five numerical and two categorical traits. The performances of the methods were evaluated under each missing mechanism by determining the mean squared error and proportion falsely classified rates for numerical and categorical traits, respectively. A random forest method supplemented with a nuclear-derived phylogeny resulted in the lowest error rates for the majority of traits, and this method was used to impute missing values in the original dataset. Data with imputed values better reflected the characteristics and distributions of the original data compared to complete-case data. However, caution should be taken when imputing trait data as phylogeny did not always improve performance for every trait and in every scenario. Ultimately, these results support the use of a real data-driven simulation strategy for selecting a suitable imputation method for a given mixed-type trait dataset.

Metagenomics versus total RNA sequencing: most accurate data-processing tools, microbial identification accuracy and perspectives for ecological assessments.

Metagenomics and total RNA sequencing (total RNA-Seq) have the potential to improve the taxonomic identification of diverse microbial communities, which could allow for the incorporation of microbes into routine ecological assessments. However, these target-PCR-free techniques require more testing and optimization. In this study, we processed metagenomics and total RNA-Seq data from a commercially available microbial mock community using 672 data-processing workflows, identified the most accurate data-processing tools, and compared their microbial identification accuracy at equal and increasing sequencing depths. The accuracy of data-processing tools substantially varied among replicates. Total RNA-Seq was more accurate than metagenomics at equal sequencing depths and even at sequencing depths almost one order of magnitude lower than those of metagenomics. We show that while data-processing tools require further exploration, total RNA-Seq might be a favorable alternative to metagenomics for target-PCR-free taxonomic identifications of microbial communities and might enable a substantial reduction in sequencing costs while maintaining accuracy. This could be particularly an advantage for routine ecological assessments, which require cost-effective yet accurate methods, and might allow for the incorporation of microbes into ecological assessments.

Phylogenetic signal of sub-arctic beetle communities.

Postglacial dispersal and colonization processes have shaped community patterns in sub-Arctic regions such as Churchill, Manitoba, and Canada. This study investigates evolutionary community structure within the beetle (Coleoptera) families of Churchill and tests whether biological traits have played a role in governing colonization patterns from refugial and southerly geographic regions. This study quantifies sub-Arctic beetle phylogenetic community structure for each family using the net relatedness index (NRI) and nearest taxon index (NTI), calculated using publicly available data from the Barcode of Life Data Systems (BOLD); compares patterns across families with different traits (habitat, diet) using standard statistical analysis (ANOVA) as well as phylogenetic generalized least squares (PGLS) using a family-level beetle phylogeny obtained from the literature; and compares community structure in Churchill with a region in southern Canada (Guelph, Ontario). These analyses were also repeated at a genus level. The dominant pattern detected in our study was that aquatic families were much better represented in Churchill compared to terrestrial families, when compared against richness sampled from across Canada and Alaska. Individually, most families showed significant phylogenetic clustering in Churchill, likely due to the strong environmental filtering present in Arctic environments. There was no significant difference in phylogenetic structure between Churchill and Guelph but with a trend toward stronger clustering in the North. Fungivores were significantly more overdispersed than other feeding modes, predators were significantly more clustered, and aquatic families showed significantly stronger clustering compared to terrestrial. This study contributes to our understanding of the traits and processes structuring insect biodiversity and macroecological trends in the sub-Arctic.

Debar: A sequence-by-sequence denoiser for COI-5P DNA barcode data.

DNA barcoding and metabarcoding are now widely used to advance species discovery and biodiversity assessments. High-throughput sequencing (HTS) has expanded the volume and scope of these analyses, but elevated error rates introduce noise into sequence records that can inflate estimates of biodiversity. Denoising -the separation of biological signal from instrument (technical) noise-of barcode and metabarcode data currently employs abundance-based methods which do not capitalize on the highly conserved structure of the cytochrome c oxidase subunit I (COI) region employed as the animal barcode. This manuscript introduces debar, an R package that utilizes a profile hidden Markov model to denoise indel errors in COI sequences introduced by instrument error. In silico studies demonstrated that debar recognized 95% of artificially introduced indels in COI sequences. When applied to real-world data, debar reduced indel errors in circular consensus sequences obtained with the Sequel platform by 75%, and those generated on the Ion Torrent S5 by 94%. The false correction rate was less than 0.1%, indicating that debar is receptive to the majority of true COI variation in the animal kingdom. In conclusion, the debar package improves DNA barcode and metabarcode workflows by aiding the generation of more accurate sequences aiding the characterization of species diversity.

The Effects of Ecological Traits on the Rate of Molecular Evolution in Ray-Finned Fishes: A Multivariable Approach.

Myriad environmental and biological traits have been investigated for their roles in influencing the rate of molecular evolution across various taxonomic groups. However, most studies have focused on a single trait, while controlling for additional factors in an informal way, generally by excluding taxa. This study utilized a dataset of cytochrome c oxidase subunit I (COI) barcode sequences from over 7000 ray-finned fish species to test the effects of 27 traits on molecular evolutionary rates. Environmental traits such as temperature were considered, as were traits associated with effective population size including body size and age at maturity. It was hypothesized that these traits would demonstrate significant correlations with substitution rate in a multivariable analysis due to their associations with mutation and fixation rates, respectively. A bioinformatics pipeline was developed to assemble and analyze sequence data retrieved from the Barcode of Life Data System (BOLD) and trait data obtained from FishBase. For use in phylogenetic regression analyses, a maximum likelihood tree was constructed from the COI sequence data using a multi-gene backbone constraint tree covering 71% of the species. A variable selection method that included both single- and multivariable analyses was used to identify traits that contribute to rate heterogeneity estimated from different codon positions. Our analyses revealed that molecular rates associated most significantly with latitude, body size, and habitat type. Overall, this study presents a novel and systematic approach for integrative data assembly and variable selection methodology in a phylogenetic framework.

coil: an R package for cytochrome c oxidase I (COI) DNA barcode data cleaning, translation, and error evaluation.

Biological conclusions based on DNA barcoding and metabarcoding analyses can be strongly influenced by the methods utilized for data generation and curation, leading to varying levels of success in the separation of biological variation from experimental error. The 5' region of cytochrome c oxidase subunit I (COI-5P) is the most common barcode gene for animals, with conserved structure and function that allows for biologically informed error identification. Here, we present coil ( https://CRAN.R-project.org/package=coil ), an R package for the pre-processing and frameshift error assessment of COI-5P animal barcode and metabarcode sequence data. The package contains functions for placement of barcodes into a common reading frame, accurate translation of sequences to amino acids, and highlighting insertion and deletion errors. The analysis of 10 000 barcode sequences of varying quality demonstrated how coil can place barcode sequences in reading frame and distinguish sequences containing indel errors from error-free sequences with greater than 97.5% accuracy. Package limitations were tested through the analysis of COI-5P sequences from the plant and fungal kingdoms as well as the analysis of potential contaminants: nuclear mitochondrial pseudogenes and Wolbachia COI-5P sequences. Results demonstrated that coil is a strong technical error identification method but is not reliable for detecting all biological contaminants.

A new species of Archaebranchinecta (Anostraca: Branchinectidae) from the South American Altiplano.

A detailed morphological comparison was carried out among specimens of several samples of Archaebranchinecta Rogers Coronel, 2011 from the Altiplano of Peru, Bolivia, and Argentina. Surprisingly, striking differences were found between Peruvian samples collected near the western shore of Lake Titicaca, and those from Bolivia taken southwardly, near the east coasts of River Desaguadero and Lake Poopó. Accordingly, the new species Archaebranchinecta aimara sp. nov. is described, representing the second specific entity of a genus that so far included only A. pollicifera (Harding, 1940). The main differential features between both species include: (a) size and shape of the three processes of basal segment of the male second antenna; (b) shape and protrusion degree of the pair of medioventral bulges in the male genital segments; (c) presence or absence of strong ventrolateral spine on the second genital segment of female; and (d) presence or absence of a pair of ventrolateral outgrowths in the brood pouch. This research contributes to our understanding of the biodiversity and endemism of the unique Altiplano region of South America.

Climatic and vegetational drivers of insect beta diversity at the continental scale.

AIM: We construct a framework for mapping pattern and drivers of insect diversity at the continental scale and use it to test whether and which environmental gradients drive insect beta diversity. LOCATION: Global; North and Central America; Western Europe. TIME PERIOD: 21st century. MAJOR TAXA STUDIED: Insects. METHODS: An informatics system was developed to integrate terrestrial data on insects with environmental parameters. We mined repositories of data for distribution, climatic data were retrieved (WorldClim), and vegetation parameters inferred from remote sensing analysis (MODIS Vegetation Continuous Fields). Beta diversity between sites was calculated and then modeled with two methods, Mantel test with multiple regression and generalized dissimilarity modeling. RESULTS: Geographic distance was the main driver of insect beta diversity. Independent of geographic distance, bioclimate variables explained more variance in dissimilarity than vegetation variables, although the particular variables found to be significant were more consistent in the latter, particularly, tree cover. Tree cover gradients drove compositional dissimilarity at denser coverages, in both continental case studies. For climate, gradients in temperature parameters were significant in driving beta diversity more so than gradients in precipitation parameters. MAIN CONCLUSIONS: Although environmental gradients drive insect beta diversity independently of geography, the relative contribution of different climatic and vegetational parameters is not expected to be consistent in different study systems. With further incorporation of additional temporal information and variables, this approach will enable the development of a predictive framework for conserving insect biodiversity at the global scale.

Is molecular evolution faster in the tropics?

The evolutionary speed hypothesis (ESH) suggests that molecular evolutionary rates are higher among species inhabiting warmer environments. Previously, the ESH has been investigated using small numbers of latitudinally-separated sister lineages; in animals, these studies typically focused on subsets of Chordata and yielded mixed support for the ESH. This study analyzed public DNA barcode sequences from the cytochrome c oxidase subunit I (COI) gene for six of the largest animal phyla (Arthropoda, Chordata, Mollusca, Annelida, Echinodermata, and Cnidaria) and paired latitudinally-separated taxa together informatically. Of 8037 lineage pairs, just over half (51.6%) displayed a higher molecular rate in the lineage inhabiting latitudes closer to the equator, while the remainder (48.4%) displayed a higher rate in the higher-latitude lineage. To date, this study represents the most comprehensive analysis of latitude-related molecular rate differences across animals. While a statistically-significant pattern was detected from our large sample size, our findings suggest that the EHS may not serve as a strong universal mechanism underlying the latitudinal diversity gradient and that COI molecular clocks may generally be applied across latitudes. This study also highlights the merits of using automation to analyze large DNA barcode datasets.

Recalibrating the molecular clock for Arctic marine invertebrates based on DNA barcodes 1.

Divergence times for species assemblages of Arctic marine invertebrates have often been estimated using a standard rate (1.4%/MY) of molecular evolution calibrated using a single sister pair of tropical crustaceans. Because rates of molecular evolution vary among taxa and environments, it is essential to obtain clock calibrations from northern lineages. The recurrent opening and closure of the Bering Strait provide an exceptional opportunity for clock calibration. Here, we apply the iterative calibration approach to investigate patterns of molecular divergence among lineages of northern marine molluscs and arthropods using publicly available sequences of the cytochrome c oxidase subunit I (COI) gene and compare these results with previous estimates of trans-Bering divergences for echinoderms and polychaetes. The wide range of Kimura two-parameter (K2P) divergences among 73 trans-Bering sister pairs (0.12%-16.89%) supports multiple pulses of migration through the Strait. Overall, the results indicate a rate of K2P divergence of 3.2%/MY in molluscs, 5%-5.2%/MY in arthropods, and 3.5%-4.7%/MY in polychaetes. While these rates are considerably higher than the often-adopted 1.4%/MY rate, they are similar to calibrations (3%-5%/MY) in several other studies of marine invertebrates. This upward revision in rates means there is a need both to reevaluate the evolutionary history of marine lineages and to reexamine the impact of prior climatic changes upon the diversification of marine life.

DNA barcode data reveal biogeographic trends in Arctic non-biting midges.

Chironomid flies (non-biting midges) are among the most abundant and diverse animals in Arctic regions, but detailed analyses of species distributions and biogeographical patterns are hampered by challenging taxonomy and reliance on morphology for species-level identification. Here we take advantage of available DNA barcode data of Arctic Chironomidae in BOLD to analyse similarities in species distributions across a northern Nearctic - West Palearctic gradient. Using more than 260 000 barcodes representing 4666 BINs (Barcode Index Numbers) and 826 named species (some with interim names) from a combination of public and novel data, we show that the Greenland chironomid fauna shows affinities to both the Nearctic and the West Palearctic regions. While raw taxon counts indicate a strong Greenland - North American affinity, comparisons using Chao's dissimilarity metric support a slightly higher similarity between Greenland and West Palearctic chironomid communities. Results were relatively consistent across different definitions of species taxonomic units, including morphologically determined species, BINs, and superBINs based on a ∼4.5% threshold. While most taxa found in Greenland are shared with at least one other region, reflecting circum-Arctic dispersal, our results also reveal that Greenland harbours a small endemic biodiversity. Our exploratory study showcases how DNA barcoding efforts using standardized gene regions contribute to an understanding of broad-scale patterns in biogeography by enabling joint analysis of public DNA sequence data derived from diverse prior studies.

The Hyalella (Crustacea: Amphipoda) species cloud of the ancient Lake Titicaca originated from multiple colonizations.

Ancient lakes are renowned for their exceptional diversity of endemic species. As model systems for the study of sympatric speciation, it is necessary to understand whether a given hypothesized species flock is of monophyletic or polyphyletic origin. Here, we present the first molecular characterization of the Hyalella (Crustacea: Amphipoda) species complex of Lake Titicaca, using COI and 28S DNA sequences, including samples from the connected Small and Large Lakes that comprise Lake Titicaca as well as from a broader survey of southern South American sites. At least five evolutionarily distant lineages are present within Lake Titicaca, which were estimated to have diverged from one another 12-20 MYA. These major lineages are dispersed throughout the broader South American Hyalella phylogeny, with each lineage representing at least one independent colonization of the lake. Moreover, complex genetic relationships are revealed between Lake Titicaca individuals and those from surrounding water bodies, which may be explained by repeated dispersal into and out of the lake, combined with parallel intralacustrine diversification within two separate clades. Although further work in deeper waters will be required to determine the number of species present and modes of diversification, our results strongly indicate that this amphipod species cloud is polyphyletic with a complex geographic history.

Iterative Calibration: A Novel Approach for Calibrating the Molecular Clock Using Complex Geological Events.

During the past 50 years, the molecular clock has become one of the main tools for providing a time scale for the history of life. In the era of robust molecular evolutionary analysis, clock calibration is still one of the most basic steps needing attention. When fossil records are limited, well-dated geological events are the main resource for calibration. However, biogeographic calibrations have often been used in a simplistic manner, for example assuming simultaneous vicariant divergence of multiple sister lineages. Here, we propose a novel iterative calibration approach to define the most appropriate calibration date by seeking congruence between the dates assigned to multiple allopatric divergences and the geological history. Exploring patterns of molecular divergence in 16 trans-Bering sister clades of echinoderms, we demonstrate that the iterative calibration is predominantly advantageous when using complex geological or climatological events-such as the opening/reclosure of the Bering Strait-providing a powerful tool for clock dating that can be applied to other biogeographic calibration systems and further taxa. Using Bayesian analysis, we observed that evolutionary rate variability in the COI-5P gene is generally distributed in a clock-like fashion for Northern echinoderms. The results reveal a large range of genetic divergences, consistent with multiple pulses of trans-Bering migrations. A resulting rate of 2.8% pairwise Kimura-2-parameter sequence divergence per million years is suggested for the COI-5P gene in Northern echinoderms. Given that molecular rates may vary across latitudes and taxa, this study provides a new context for dating the evolutionary history of Arctic marine life.

Host Specificity in Subarctic Aphids.

Plants and herbivorous (or parasitic) insects form the majority of macroscopic life. The specificity of interaction between host plant and parasitic insect depends on the adaptations of both the host and the parasite. Over time, these interactions evolve and change as a result of an 'arms race' between host and parasite, and the resulting species-specific adaptations may be maintained, perpetuating these interactions across speciation events. This can lead to specialisation between species or clades. With speciation and species sorting over time, complex interactions evolve. Here, we elucidate a three-tier method to test these interactions using the aphids (Hemiptera: Aphididae) and plants of Churchill (Manitoba, Canada) as a model system. We analyzed these interactions by testing for three patterns in host specificity: monophagy, phylogenetic clustering, and cophylogeny. We defined monophagy strictly as one species feeding exclusively upon a single host plant species (an association likely driven by arms races in morphology, chemical resistance/tolerance, and visual appearance) and observed this in 7 of 22 aphid species. In all the remaining 'polyphagous' cases, there was a strong trend toward monophagy (80% of individuals were found on a single host plant species). Second, we observed two separate examples of phylogenetic clustering where groups of closely related aphid species fed upon individual plant species. Finally, we found no support for cophylogenetic relationships where both aphids and plants cospeciate to form congruent phylogenetic trees (evidence of coadaptation through an ongoing arms race). One explanation for uncovering species-specific interactions in a recently deglaciated, subarctic locality is that the species involved in the associations moved north together. Testing different levels of specificity in the most predominant species-species interactions on the planet will allow us to elucidate these patterns accurately and gives us insight into where to direct future research.

International Barcode of Life: Focus on big biodiversity in South Africa.

Participants in the 7th International Barcode of Life Conference (Kruger National Park, South Africa, 20-24 November 2017) share the latest findings in DNA barcoding research and its increasingly diversified applications. Here, we review prevailing trends synthesized from among 429 invited and contributed abstracts, which are collated in this open-access special issue of Genome. Hosted for the first time on the African continent, the 7th Conference places special emphasis on the evolutionary origins, biogeography, and conservation of African flora and fauna. Within Africa and elsewhere, DNA barcoding and related techniques are being increasingly used for wildlife forensics and for the validation of commercial products, such as medicinal plants and seafood species. A striking trend of the conference is the dramatic rise of studies on environmental DNA (eDNA) and on diverse uses of high-throughput sequencing techniques. Emerging techniques in these areas are opening new avenues for environmental biomonitoring, managing species-at-risk and invasive species, and revealing species interaction networks in unprecedented detail. Contributors call for the development of validated community standards for high-throughput sequence data generation and analysis, to enable the full potential of these methods to be realized for understanding and managing biodiversity on a global scale.

Positive and relaxed selection associated with flight evolution and loss in insect transcriptomes.

The evolution of powered flight is a major innovation that has facilitated the success of insects. Previously, studies of birds, bats, and insects have detected molecular signatures of differing selection regimes in energy-related genes associated with flight evolution and/or loss. Here, using DNA sequences from more than 1000 nuclear and mitochondrial protein-coding genes obtained from insect transcriptomes, we conduct a broader exploration of which gene categories display positive and relaxed selection at the origin of flight as well as with multiple independent losses of flight. We detected a number of categories of nuclear genes more often under positive selection in the lineage leading to the winged insects (Pterygota), related to catabolic processes such as proteases, as well as splicing-related genes. Flight loss was associated with relaxed selection signatures in splicing genes, mirroring the results for flight evolution. Similar to previous studies of flight loss in various animal taxa, we observed consistently higher nonsynonymous-to-synonymous substitution ratios in mitochondrial genes of flightless lineages, indicative of relaxed selection in energy-related genes. While oxidative phosphorylation genes were not detected as being under selection with the origin of flight specifically, they were most often detected as being under positive selection in holometabolous (complete metamorphosis) insects as compared with other insect lineages. This study supports some convergence in gene-specific selection pressures associated with flight ability, and the exploratory analysis provided some new insights into gene categories potentially associated with the gain and loss of flight in insects.

Barcode-based species delimitation in the marine realm: a test using Hexanauplia (Multicrustacea: Thecostraca and Copepoda).

DNA barcoding has been used successfully for identifying specimens belonging to marine planktonic groups. However, the ability to delineate species within taxonomically diverse and widely distributed marine groups, such as the Copepoda and Thecostraca, remains largely untested. We investigate whether a cytochrome c oxidase subunit I (COI-5P) global pairwise sequence divergence threshold exists between intraspecific and interspecific divergences in the copepods plus the thecostracans (barnacles and allies). Using publicly accessible sequence data, we applied a graphical method to determine an optimal threshold value. With these thresholds, and using a newly generated planktonic marine data set, we quantify the degree of concordance using a bidirectional analysis and discuss different analytical methods for sequence-based species delimitation (e.g., BIN, ABGD, jMOTU, UPARSE, Mothur, PTP, and GMYC). Our results support a COI-5P threshold between 2.1% and 2.6% p-distance across methods for these crustacean taxa, yielding molecular groupings largely concordant with traditional, morphologically defined species. The adoption of internal methods for clustering verification enables rapid biodiversity studies and the exploration of unknown faunas using DNA barcoding. The approaches taken here for concordance assessment also provide a more quantitative comparison of clustering results (as contrasted with "success/failure" of barcoding), and we recommend their further consideration for barcoding studies.

Forecasting pollination declines through DNA barcoding: the potential contributions of macroecological and macroevolutionary scales of inquiry.

While pollinators are widely acknowledged as important contributors to seed production in plant communities, we do not yet have a good understanding of the importance of pollinator specialists for this ecosystem service. Determination of the prevalence of pollinator specialists is often hindered by the occurrence of cryptic species and the limitations of observational data on pollinator visitation rates, two areas where DNA barcoding of pollinators and pollen can be useful. Further, the demonstrated adequacy of pollen DNA barcoding from historical records offers opportunities to observe the effects of pollinator loss over longer timescales, and phylogenetic approaches can elucidate the historical rates of extinction of specialist lineages. In this Viewpoint article, we review how advances in DNA barcoding and metabarcoding of plants and pollinators have brought important developments to our understanding of specialization in plant-pollinator interactions. We then put forth several lines of inquiry that we feel are especially promising for providing insight on changes in plant-pollinator interactions over space and time. Obtaining estimates of the effects of reductions in specialists will contribute to forecasting the loss of ecosystem services that will accompany the erosion of plant and pollinator diversity.

The importance of taxonomic resolution for additive beta diversity as revealed through DNA barcoding.

Additive diversity partitioning (α, ß, and γ) is commonly used to study the distribution of species-level diversity across spatial scales. Here, we first investigate whether published studies of additive diversity partitioning show signs of difficulty attaining species-level resolution due to inherent limitations with morphological identifications. Second, we present a DNA barcoding approach to delineate specimens of stream caddisfly larvae (order Trichoptera) and consider the importance of taxonomic resolution on classical (additive) measures of beta (ß) diversity. Caddisfly larvae were sampled using a hierarchical spatial design in two regions (subarctic Churchill, Manitoba, Canada; temperate Pennsylvania, USA) and then additively partitioned according to Barcode Index Numbers (molecular clusters that serve as a proxy for species), genus, and family levels; diversity components were expressed as proportional species turnover. We screened 114 articles of additive diversity partitioning and found that a third reported difficulties with achieving species-level identifications, with a clear taxonomic tendency towards challenges identifying invertebrate taxa. Regarding our own study, caddisfly BINs appeared to show greater subregional turnover (e.g., proportional additive ß) compared to genus or family levels. Diversity component studies failing to achieve species resolution due to morphological identifications may therefore be underestimating diversity turnover at larger spatial scales.