nRCFV: a new, dataset-size-independent metric to quantify compositional heterogeneity in nucleotide and amino acid datasets.

MOTIVATION: Compositional heterogeneity-when the proportions of nucleotides and amino acids are not broadly similar across the dataset-is a cause of a great number of phylogenetic artefacts. Whilst a variety of methods can identify it post-hoc, few metrics exist to quantify compositional heterogeneity prior to the computationally intensive task of phylogenetic tree reconstruction. Here we assess the efficacy of one such existing, widely used, metric: Relative Composition Frequency Variability (RCFV), using both real and simulated data. RESULTS: Our results show that RCFV can be biased by sequence length, the number of taxa, and the number of possible character states within the dataset. However, we also find that missing data does not appear to have an appreciable effect on RCFV. We discuss the theory behind this, the consequences of this for the future of the usage of the RCFV value and propose a new metric, nRCFV, which accounts for these biases. Alongside this, we present a new software that calculates both RCFV and nRCFV, called nRCFV_Reader. AVAILABILITY AND IMPLEMENTATION: nRCFV has been implemented in RCFV_Reader, available at: https://github.com/JFFleming/RCFV_Reader . Both our simulation and real data are available at Datadryad: https://doi.org/10.5061/dryad.wpzgmsbpn .

Molecular palaeontology illuminates the evolution of ecdysozoan vision.

Colour vision is known to have arisen only twice-once in Vertebrata and once within the Ecdysozoa, in Arthropoda. However, the evolutionary history of ecdysozoan vision is unclear. At the molecular level, visual pigments, composed of a chromophore and a protein belonging to the opsin family, have different spectral sensitivities and these mediate colour vision. At the morphological level, ecdysozoan vision is conveyed by eyes of variable levels of complexity; from the simple ocelli observed in the velvet worms (phylum Onychophora) to the marvellously complex eyes of insects, spiders, and crustaceans. Here, we explore the evolution of ecdysozoan vision at both the molecular and morphological level; combining analysis of a large-scale opsin dataset that includes previously unknown ecdysozoan opsins with morphological analyses of key Cambrian fossils with preserved eye structures. We found that while several non-arthropod ecdysozoan lineages have multiple opsins, arthropod multi-opsin vision evolved through a series of gene duplications that were fixed in a period of 35-71 million years (Ma) along the stem arthropod lineage. Our integrative study of the fossil and molecular record of vision indicates that fossils with more complex eyes were likely to have possessed a larger complement of opsin genes.

Uncertain-tree: discriminating among competing approaches to the phylogenetic analysis of phenotype data.

Morphological data provide the only means of classifying the majority of life's history, but the choice between competing phylogenetic methods for the analysis of morphology is unclear. Traditionally, parsimony methods have been favoured but recent studies have shown that these approaches are less accurate than the Bayesian implementation of the Mk model. Here we expand on these findings in several ways: we assess the impact of tree shape and maximum-likelihood estimation using the Mk model, as well as analysing data composed of both binary and multistate characters. We find that all methods struggle to correctly resolve deep clades within asymmetric trees, and when analysing small character matrices. The Bayesian Mk model is the most accurate method for estimating topology, but with lower resolution than other methods. Equal weights parsimony is more accurate than implied weights parsimony, and maximum-likelihood estimation using the Mk model is the least accurate method. We conclude that the Bayesian implementation of the Mk model should be the default method for phylogenetic estimation from phenotype datasets, and we explore the implications of our simulations in reanalysing several empirical morphological character matrices. A consequence of our finding is that high levels of resolution or the ability to classify species or groups with much confidence should not be expected when using small datasets. It is now necessary to depart from the traditional parsimony paradigms of constructing character matrices, towards datasets constructed explicitly for Bayesian methods.

The Evolution of Temperature and Desiccation-Related Protein Families in Tardigrada Reveals a Complex Acquisition of Extremotolerance.

Tardigrada is an ecdysozoan lineage famed for its resilience. Tardigrades can tolerate high doses of radiation, low-oxygen environments, desiccation, and both high and low temperatures under a dormant state called "anhydrobiosis", which is a reversible halt of metabolism upon almost complete desiccation. A large amount of research has focused on the genetic pathways related to these capabilities, and a number of genes have been identified and linked to the extremotolerant response of tardigrades. However, the history of these genes is unclear, and the origins and history of extremotolerant genes within Tardigrada remain a mystery. Here, we generate the first phylogenies of six separate protein families linked with desiccation and radiation tolerance in Tardigrada: cytosolic abundant heat-soluble protein, mitochondrial abundant heat-soluble protein, secretory abundant heat-soluble protein, meiotic recombination 11 homolog, and the newly discovered Echiniscus testudo abundant heat-soluble proteins (alpha and beta). The high number of independent gene duplications found amongst the six gene families studied suggests that tardigrades have a complex history with numerous independent adaptations to cope with aridity within the limnoterrestrial environment. Our results suggest that tardigrades likely transitioned from a marine environment to a limnoterrestrial environment only twice, once in stem Eutardigrada and once in Heterotardigrada, which explains the unique adaptations to anhydrobiosis present in both classes.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

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.

New Tardigrade Opsins and Differential Expression Analyses Show Ontogenic Variation in Light Perception.

Opsins are light-sensitive proteins involved in many photoreceptive processes, including, but not limited to, vision and regulation of circadian rhythms. Arthropod (e.g., insects, spiders, centipedes, and crabs) opsins have been extensively researched, but the relationships and function of opsins found in lineages that are evolutionarily closely related to the arthropods remains unclear. Multiple, independent, opsin duplications are known in Tardigrada (the water bears), evidencing that protostome opsin duplications are not limited to the Arthropoda. However, the relationships, function, and expression of these new opsins are still unknown. Here, we use two tardigrade transcriptomes with deep coverage to greatly expand our knowledge of the diversity of tardigrade opsins. We reconstruct the phylogenetic relationships of the tardigrade opsins and investigate their ontogenetic expression. We found that while tardigrades have multiple opsins that evolved from lineage-specific duplications of well-understood arthropod opsins, their expression levels change during ontogeny such that most of these opsins are not co-temporally expressed. Co-temporal expression of multiple opsins underpins color vision in Arthropoda and Vertebrata. Our results clearly show duplications of both rhabdomeric and ciliary opsins within Tardigrada, forming clades specific to both the Heterotardigrada and Eutardigrada in addition to multiple independent duplications within genera. However, lack of co-temporal, ontogenetic, expression suggests that while tardigrades possess multiple opsins, they are unlikely to be able to distinguish color.

Comprehensive Transcriptome Sequencing of Tanaidacea with Proteomic Evidences for Their Silk.

Tanaidaceans are small benthic crustaceans that mainly inhabit diverse marine environments, and they comprise one of the most diverse and abundant macrofaunal groups in the deep sea. Tanaidacea is one of the most thread-dependent taxa in the Crustacea, constructing tubes, spun with their silk, for shelter. In this work, we sequenced and assembled the comprehensive transcriptome of 23 tanaidaceans encompassing 14 families and 4 superfamilies of Tanaidacea, and performed silk proteomics of Zeuxo ezoensis to search for its silk genes. As a result, we identified two families of silk proteins that are conserved across the four superfamilies. The long and repetitive nature of these silk genes resembles that of other silk-producing organisms, and the two families of proteins are similar in composition to silkworm and caddisworm fibroins, respectively. Moreover, the amino acid composition of the repetitive motifs of tanaidacean silk tends to be more hydrophilic, and therefore could be a useful resource in studying their unique adaptation of silk use in a marine environment. The availability of comprehensive transcriptome data in these taxa, coupled with proteomic evidence of their silk genes, will facilitate evolutionary and ecological studies.

Eye-Transcriptome and Genome-Wide Sequencing for Scolecophidia: Implications for Inferring the Visual System of the Ancestral Snake.

Molecular genetic data have recently been incorporated in attempts to reconstruct the ecology of the ancestral snake, though this has been limited by a paucity of data for one of the two main extant snake taxa, the highly fossorial Scolecophidia. Here we present and analyze vision genes from the first eye-transcriptomic and genome-wide data for Scolecophidia, for Anilios bicolor, and A. bituberculatus, respectively. We also present immunohistochemistry data for retinal anatomy and visual opsin-gene expression in Anilios. Analyzed in the context of 19 lepidosaurian genomes and 12 eye transcriptomes, the new genome-wide and transcriptomic data provide evidence for a much more reduced visual system in Anilios than in non-scolecophidian (=alethinophidian) snakes and in lizards. In Anilios, there is no evidence of the presence of 7 of the 12 genes associated with alethinophidian photopic (cone) phototransduction. This indicates extensive gene loss and many of these candidate gene losses occur also in highly fossorial mammals with reduced vision. Although recent phylogenetic studies have found evidence for scolecophidian paraphyly, the loss in Anilios of visual genes that are present in alethinophidians implies that the ancestral snake had a better-developed visual system than is known for any extant scolecophidian.

A Novel Approach to Investigate the Effect of Tree Reconstruction Artifacts in Single-Gene Analysis Clarifies Opsin Evolution in Nonbilaterian Metazoans.

Our ability to correctly reconstruct a phylogenetic tree is strongly affected by both systematic errors and the amount of phylogenetic signal in the data. Current approaches to tackle tree reconstruction artifacts, such as the use of parameter-rich models, do not translate readily to single-gene alignments. This, coupled with the limited amount of phylogenetic information contained in single-gene alignments, makes gene trees particularly difficult to reconstruct. Opsin phylogeny illustrates this problem clearly. Opsins are G-protein coupled receptors utilized in photoreceptive processes across Metazoa and their protein sequences are roughly 300 amino acids long. A number of incongruent opsin phylogenies have been published and opsin evolution remains poorly understood. Here, we present a novel approach, the canary sequence approach, to investigate and potentially circumvent errors in single-gene phylogenies. First, we demonstrate our approach using two well-understood cases of long-branch attraction in single-gene data sets, and simulations. After that, we apply our approach to a large collection of well-characterized opsins to clarify the relationships of the three main opsin subfamilies.

Pancrustacean Evolution Illuminated by Taxon-Rich Genomic-Scale Data Sets with an Expanded Remipede Sampling.

The relationships of crustaceans and hexapods (Pancrustacea) have been much discussed and partially elucidated following the emergence of phylogenomic data sets. However, major uncertainties still remain regarding the position of iconic taxa such as Branchiopoda, Copepoda, Remipedia, and Cephalocarida, and the sister group relationship of hexapods. We assembled the most taxon-rich phylogenomic pancrustacean data set to date and analyzed it using a variety of methodological approaches. We prioritized low levels of missing data and found that some clades were consistently recovered independently of the analytical approach used. These include, for example, Oligostraca and Altocrustacea. Substantial support was also found for Allotriocarida, with Remipedia as the sister of Hexapoda (i.e., Labiocarida), and Branchiopoda as the sister of Labiocarida, a clade that we name Athalassocarida (="nonmarine shrimps"). Within Allotriocarida, Cephalocarida was found as the sister of Athalassocarida. Finally, moderate support was found for Hexanauplia (Copepoda as sister to Thecostraca) in alliance with Malacostraca. Mapping key crustacean tagmosis patterns and developmental characters across the revised phylogeny suggests that the ancestral pancrustacean was relatively short-bodied, with extreme body elongation and anamorphic development emerging later in pancrustacean evolution.