Revisiting metazoan phylogeny with genomic sampling of all phyla

Brachiopod and mollusc biomineralisation is a conserved process that was lost in the phoronid-bryozoan stem lineage

BACKGROUND

Brachiopods and molluscs are lophotrochozoans with hard external shells which are often believed to have evolved convergently. While palaeontological data indicate that both groups are descended from biomineralising Cambrian ancestors, the closest relatives of brachiopods, phoronids and bryozoans, are mineralised to a much lower extent and are comparatively poorly represented in the Palaeozoic fossil record. Although brachiopod and mollusc shells are structurally analogous, genomic and proteomic evidence indicates that their formation involves a complement of conserved, orthologous genes. Here, we study a set of genes comprised of 3 homeodomain transcription factors, one signalling molecule and 6 structural proteins which are implicated in mollusc and brachiopod shell formation, search for their orthologs in transcriptomes or genomes of brachiopods, phoronids and bryozoans, and present expression patterns of 8 of the genes in postmetamorphic juveniles of the rhynchonelliform brachiopod T. transversa.

RESULTS

Transcriptome and genome searches for the 10 target genes in the brachiopods Terebratalia transversa, Lingula anatina, Novocrania anomala, the bryozoans Bugula neritina and Membranipora membranacea, and the phoronids Phoronis australis and Phoronopsis harmeri resulted in the recovery of orthologs of the majority of the genes in all taxa. While the full complement of genes was present in all brachiopods with a single exception in L. anatina, a bloc of four genes could consistently not be retrieved from bryozoans and phoronids. The genes engrailed, distal-less, ferritin, perlucin, sp1 and sp2 were shown to be expressed in the biomineralising mantle margin of T. transversa juveniles.

CONCLUSIONS

The gene expression patterns we recovered indicate that while mineralised shells in brachiopods and molluscs are structurally analogous, their formation builds on a homologous process that involves a conserved complement of orthologous genes. Losses of some of the genes related to biomineralisation in bryozoans and phoronids indicate that loss of the capacity to form mineralised structures occurred already in the phoronid-bryozoan stem group and supports the idea that mineralised skeletons evolved secondarily in some of the bryozoan subclades.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

Various evaluations of the last common bilaterian ancestor (lcba) currently suggest that it resembled either a microscopic, non-segmented motile adult; or, on the contrary, a complex segmented adult motile urbilaterian. These fundamental inconsistencies remain largely unexplained. A majority of multidisciplinary data regarding sedentary adult ancestral bilaterian organization is overlooked. The sedentary-pelagic model is supported now by a number of novel developmental, paleontological and molecular phylogenetic data: (1) data in support of sedentary sponges, in the adult stage, as sister to all other Metazoa; (2) a similarity of molecular developmental pathways in both adults and larvae across sedentary sponges, cnidarians, and bilaterians; (3) a cnidarian-bilaterian relationship, including a unique sharing of a bona fide Hox-gene cluster, of which the evolutionary appearance does not connect directly to a bilaterian motile organization; (4) the presence of sedentary and tube-dwelling representatives of the main bilaterian clades in the early Cambrian; (5) an absence of definite taxonomic attribution of Ediacaran taxa reconstructed as motile to any true bilaterian phyla; (6) a similarity of tube morphology (and the clear presence of a protoconch-like apical structure of the Ediacaran sedentary Cloudinidae) among shells of the early Cambrian, and later true bilaterians, such as semi-sedentary hyoliths and motile molluscs; (7) recent data that provide growing evidence for a complex urbilaterian, despite a continuous molecular phylogenetic controversy. The present review compares the main existing models and reconciles the sedentary model of an urbilaterian and the model of a larva-like lcba with a unified sedentary(adult)-pelagic(larva) model of the lcba.

Different phylogenomic methods support monophyly of enigmatic 'Mesozoa' (Dicyemida + Orthonectida, Lophotrochozoa)

Dicyemids and orthonectids were traditionally classified in a group called Mesozoa, but their placement in a single clade has been contested and their position(s) within Metazoa is uncertain. Here, we assembled a comprehensive matrix of Lophotrochozoa (Metazoa) and investigated the position of Dicyemida (= Rhombozoa) and Orthonectida, employing multiple phylogenomic approaches. We sequenced seven new transcriptomes and one draft genome from dicyemids (Dicyema, Dicyemennea) and two transcriptomes from orthonectids (Rhopalura). Using these and published data, we assembled and analysed contamination-filtered datasets with up to 987 genes. Our results recover Mesozoa monophyletic and as a close relative of Platyhelminthes or Gnathifera. Because of the tendency of the long-branch mesozoans to group with other long-branch taxa in our analyses, we explored the impact of approaches purported to help alleviate long-branch attraction (e.g. taxon removal, coalescent inference, gene targeting). None of these were able to break the association of Orthonectida with Dicyemida in the maximum-likelihood trees. Contrastingly, the Bayesian analysis and site-specific frequency model in maximum-likelihood did not recover a monophyletic Mesozoa (but only when using a specific 50 gene matrix). The classic hypothesis on monophyletic Mesozoa is possibly reborn and should be further tested.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

Polyzoa is back: The effect of complete gene sets on the placement of Ectoprocta and Entoprocta

The phylogenomic approach has largely resolved metazoan phylogeny and improved our knowledge of animal evolution based on morphology, paleontology, and embryology. Nevertheless, the placement of two major lophotrochozoan phyla, Entoprocta (Kamptozoa) and Ectoprocta (Bryozoa), remains highly controversial: Originally considered as a single group named Polyzoa (Bryozoa), they were separated on the basis of morphology. So far, each new study of lophotrochozoan evolution has still consistently proposed different phylogenetic positions for these groups. Here, we reinvestigated the placement of Entoprocta and Ectoprocta using highly complete datasets with rigorous contamination removal. Our results from maximum likelihood, Bayesian, and coalescent analyses strongly support the topology in which Entoprocta and Bryozoa form a distinct clade, placed as a sister group to all other lophotrochozoan clades: Annelida, Mollusca, Brachiopoda, Phoronida, and Nemertea. Our study favors the evolutionary scenario where Entoprocta, Cycliophora, and Bryozoa constitute one of the earliest branches among Lophotrochozoa and thus supports the Polyzoa hypothesis.

Comparative Single-Cell Transcriptomics Reveals Novel Genes Involved in Bivalve Embryonic Shell Formation and Questions Ontogenetic Homology of Molluscan Shell Types

Mollusks are known for their highly diverse repertoire of body plans that often includes external armor in form of mineralized hardparts. Representatives of the Conchifera, one of the two major lineages that comprises taxa which originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia), are particularly relevant regarding the evolution of mollusk shells. Previous studies have found that the shell matrix of the adult shell (teleoconch) is rapidly evolving and that the gene set involved in shell formation is highly taxon-specific. However, detailed annotation of genes expressed in tissues involved in the formation of the embryonic shell (protoconch I) or the larval shell (protoconch II) are currently lacking. Here, we analyzed the genetic toolbox involved in embryonic and larval shell formation in the quagga mussel Dreissena rostriformis using single cell RNA sequencing. We found significant differences in genes expressed during embryonic and larval shell secretion, calling into question ontogenetic homology of these transitory bivalve shell types. Further ortholog comparisons throughout Metazoa indicates that a common genetic biomineralization toolbox, that was secondarily co-opted into molluscan shell formation, was already present in the last common metazoan ancestor. Genes included are engrailed, carbonic anhydrase, and tyrosinase homologs. However, we found that 25% of the genes expressed in the embryonic shell field of D. rostriformis lack an ortholog match with any other metazoan. This indicates that not only adult but also embryonic mollusk shells may be fast-evolving structures. We raise the question as to what degree, and on which taxonomic level, the gene complement involved in conchiferan protoconch formation may be lineage-specific or conserved across taxa.

Microbiomes of microscopic marine invertebrates do not reveal signatures of phylosymbiosis

Animals and microorganisms often establish close ecological relationships. However, much of our knowledge about animal microbiomes comes from two deeply studied groups: vertebrates and arthropods. To understand interactions on a broader scale of diversity, we characterized the bacterial microbiomes of close to 1,000 microscopic marine invertebrates from 21 phyla, spanning most of the remaining tree of metazoans. Samples were collected from five temperate and tropical locations covering three marine habitats (sediment, water column and intertidal macroalgae) and bacterial microbiomes were characterized using 16S ribosomal RNA gene sequencing. Our data show that, despite their size, these animals harbour bacterial communities that differ from those in the surrounding environment. Distantly related but coexisting invertebrates tend to share many of the same bacteria, suggesting that guilds of microorganisms preferentially associated with animals, but not tied to any specific host lineage, are the main drivers of the ecological relationship. Host identity is a minor factor shaping these microbiomes, which do not show the same correlation with host phylogeny, or 'phylosymbiosis', observed in many large animals. Hence, the current debate on the varying strength of phylosymbiosis within selected lineages should be reframed to account for the possibility that such a pattern might be the exception rather than the rule.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

Background

Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.

Results

P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.

Conclusions

Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08607-4.

Ontogeny, Phylotypic Periods, Paedomorphosis, and Ontogenetic Systematics

The key terms linking ontogeny and evolution are briefly reviewed. It is shown that their application and usage in the modern biology are often inconsistent and incorrectly understood even within the “evo-devo” field. For instance, the core modern reformulation that ontogeny not merely recapitulates, but produces phylogeny implies that ontogeny and phylogeny are closely interconnected. However, the vast modern phylogenetic and taxonomic fields largely omit ontogeny as a central concept. Instead, the common “clade-” and “tree-thinking” prevail, despite on the all achievements of the evo-devo. This is because the main conceptual basis of the modern biology is fundamentally ontogeny-free. In another words, in the Haeckel’s pair of “ontogeny and phylogeny,” ontogeny is still just a subsidiary for the evolutionary process (and hence, phylogeny), instead as in reality, its main driving force. The phylotypic periods is another important term of the evo-devo and represent a modern reformulation of Haeckel’s recapitulations and biogenetic law. However, surprisingly, this one of the most important biological evidence, based on the natural ontogenetic grounds, in the phylogenetic field that can be alleged as a “non-evolutionary concept.” All these observations clearly imply that a major revision of the main terms which are associated with the “ontogeny and phylogeny/evolution” field is urgently necessarily. Thus, “ontogenetic” is not just an endless addition to the term “systematics,” but instead a crucial term, without it neither systematics, nor biology have sense. To consistently employ the modern ontogenetic and epigenetic achievements, the concept of ontogenetic systematics is hereby refined. Ontogenetic systematics is not merely a “research program” but a key biological discipline which consistently links the enormous biological diversity with underlying fundamental process of ontogeny at both molecular and morphological levels. The paedomorphosis is another widespread ontogenetic-and-evolutionary process that is significantly underestimated or misinterpreted by the current phylogenetics and taxonomy. The term paedomorphosis is refined, as initially proposed to link ontogeny with evolution, whereas “neoteny” and “progenesis” are originally specific, narrow terms without evolutionary context, and should not be used as synonyms of paedomorphosis. Examples of application of the principles of ontogenetic systematics represented by such disparate animal groups as nudibranch molluscs and ophiuroid echinoderms clearly demonstrate that perseverance of the phylotypic periods is based not only on the classic examples in vertebrates, but it is a universal phenomenon in all organisms, including disparate animal phyla.

Neurotransmitters of sleep and wakefulness in flatworms

STUDY OBJECTIVES

Sleep is a prominent behavioral and biochemical state observed in all animals studied, including platyhelminth flatworms. Investigations into the biochemical mechanisms associated with sleep-and wakefulness-are important for understanding how these states are regulated and how that regulation changed with the evolution of new types of animals. Unfortunately, beyond a handful of vertebrates, such studies on invertebrates are rare.

METHODS

We investigated the effect of seven neurotransmitters, and one pharmacological compound, that modulate either sleep or wakefulness in mammals, on flatworms (Girardia tigrina). Flatworms were exposed via ingestion and diffusion to four neurotransmitters that promote wakefulness in vertebrates (acetylcholine, dopamine, glutamate, histamine), and three that induce sleep (adenosine, GABA, serotonin) along with the H1 histamine receptor antagonist pyrilamine. Compounds were administered over concentrations spanning three to five orders of magnitude. Flatworms were then transferred to fresh water and video recorded for analysis.

RESULTS

Dopamine and histamine decreased the time spent inactive and increased distance traveled, consistent with their wake-promoting effect in vertebrates and fruit flies; pyrilamine increased restfulness and GABA showed a nonsignificant trend towards promoting restfulness in a dose-dependent manner, in agreement with their sleep-inducing effect in vertebrates, fruit flies, and Hydra. Similar to Hydra, acetylcholine, glutamate, and serotonin, but also adenosine, had no apparent effect on flatworm behavior.

CONCLUSIONS

These data demonstrate the potential of neurotransmitters to regulate sleep and wakefulness in flatworms and highlight the conserved action of some neurotransmitters across species.

Paleozoic origins of cheilostome bryozoans and their parental care inferred by a new genome-skimmed phylogeny

Phylogenetic relationships and the timing of evolutionary events are essential for understanding evolution on longer time scales. Cheilostome bryozoans are a group of ubiquitous, species-rich, marine colonial organisms with an excellent fossil record but lack phylogenetic relationships inferred from molecular data. We present genome-skimmed data for 395 cheilostomes and combine these with 315 published sequences to infer relationships and the timing of key events among c. 500 cheilostome species. We find that named cheilostome genera and species are phylogenetically coherent, rendering fossil or contemporary specimens readily delimited using only skeletal morphology. Our phylogeny shows that parental care in the form of brooding evolved several times independently but was never lost in cheilostomes. Our fossil calibration, robust to varied assumptions, indicates that the cheilostome lineage and parental care therein could have Paleozoic origins, much older than the first known fossil record of cheilostomes in the Late Jurassic.

First Modern Data on the Lophophore Nervous System in Adult Novocrania anomala and a Current Assessment of Brachiopod Phylogeny

Simple Summary

The nervous system of Novocrania anomala adults is described for the first time. A table containing data on the lophophore innervation in species from three brachiopod subphyla is presented. A comparative analysis suggests a close relationship between the Craniiformea and the Rhynchonelliformea, and thereby supports the “Calciata” hypothesis of brachiopod phylogeny.

Abstract

Although the lophophore is regarded as the main synapomorphy of all lophophorates, the evolution of the lophophore in certain groups of lophophorates remains unclear. To date, the innervation of the lophophore has been studied with modern methods only for three brachiopod species belonging to two subphyla: Linguliformea and Rhynchonelliformea. In the third subphylum, the Craniiformea, there are data for juveniles but not for adults. In the current research, the innervation of the lophophore in Novocrania anomala adults was studied by immunocytochemistry and confocal laser scanning microscopy. In the spiral lophophore of adults of the craniiform N. anomala, each arm is innervated by six brachial nerves: main, additional main, accessory, second accessory, additional lower, and lower brachial nerves. Compared with other brachiopod species, this complex innervation of the lophophore correlates with the presence of many lophophoral muscles. The general anatomy of the lophophore nervous system and the peculiarities of the organization of the subenteric ganglion of the craniiform N. anomala have a lot in common with those of rhynchonelliforms but not with those of linguliforms. These findings are consistent with the “Calciata” hypothesis of the brachiopod phylogeny and are inconsistent with the inference that the Craniiformea and Linguliformea are closely related.

Examples of Extreme Survival: Tardigrade Genomics and Molecular Anhydrobiology

Tardigrades are ubiquitous meiofauna that are especially renowned for their exceptional extremotolerance to various adverse environments, including pressure, temperature, and even ionizing radiation. This is achieved through a reversible halt of metabolism triggered by desiccation, a phenomenon called anhydrobiosis. Recent establishment of genome resources for two tardigrades, Hypsibius exemplaris and Ramazzottius varieornatus, accelerated research to uncover the molecular mechanisms behind anhydrobiosis, leading to the discovery of many tardigrade-unique proteins. This review focuses on the history, methods, discoveries, and current state and challenges regarding tardigrade genomics, with an emphasis on molecular anhydrobiology. Remaining questions and future perspectives regarding prospective approaches to fully elucidate the molecular machinery of this complex phenomenon are discussed.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

Evolution of DNA Methylation Across Ecdysozoa

DNA methylation is a crucial, abundant mechanism of gene regulation in vertebrates. It is less prevalent in many other metazoan organisms and completely absent in some key model species, such as Drosophila melanogaster and Caenorhabditis elegans. We report here a comprehensive study of the presence and absence of DNA methyltransferases (DNMTs) in 138 Ecdysozoa, covering Arthropoda, Nematoda, Priapulida, Onychophora, and Tardigrada. Three of these phyla have not been investigated for the presence of DNA methylation before. We observe that the loss of individual DNMTs independently occurred multiple times across ecdysozoan phyla. We computationally predict the presence of DNA methylation based on CpG rates in coding sequences using an implementation of Gaussian Mixture Modeling, MethMod. Integrating both analysis we predict two previously unknown losses of DNA methylation in Ecdysozoa, one within Chelicerata (Mesostigmata) and one in Tardigrada. In the early-branching Ecdysozoa Priapulus caudatus, we predict the presence of a full set of DNMTs and the presence of DNA methylation. We are therefore showing a very diverse and independent evolution of DNA methylation in different ecdysozoan phyla spanning a phylogenetic range of more than 700 million years.

Functional analysis of ADARs in planarians supports a bilaterian ancestral role in suppressing double-stranded RNA-response

ADARs (adenosine deaminases acting on RNA) are known for their adenosine-to-inosine RNA editing activity, and most recently, for their role in preventing aberrant dsRNA-response by activation of dsRNA sensors (i.e., RIG-I-like receptor homologs). However, it is still unclear whether suppressing spurious dsRNA-response represents the ancestral role of ADARs in bilaterians. As a first step to address this question, we identified ADAR1 and ADAR2 homologs in the planarian Schmidtea mediterranea, which is evolutionarily distant from canonical lab models (e.g., flies and nematodes). Our results indicate that knockdown of either planarian adar1 or adar2 by RNA interference (RNAi) resulted in upregulation of dsRNA-response genes, including three planarian rig-I-like receptor (prlr) homologs. Furthermore, independent knockdown of adar1 and adar2 reduced the number of infected cells with a dsRNA virus, suggesting they suppress a bona fide anti-viral dsRNA-response activity. Knockdown of adar1 also resulted in lesion formation and animal lethality, thus attesting to its essentiality. Simultaneous knockdown of adar1 and prlr1 rescued adar1(RNAi)-dependent animal lethality and rescued the dsRNA-response, suggesting that it contributes to the deleterious effect of adar1 knockdown. Finally, we found that ADAR2, but not ADAR1, mediates mRNA editing in planarians, suggesting at least in part non-redundant activities for planarians ADARs. Our results underline the essential role of ADARs in suppressing activation of harmful dsRNA-response in planarians, thus supporting it as their ancestral role in bilaterians. Our work also set the stage to study further and better understand the regulatory mechanisms governing anti-viral dsRNA-responses from an evolutionary standpoint using planarians as a model.

Today, more than ever, it is crucial to gain a deep understating of our anti-viral defenses. One of the ways to accomplish it is to study the principles governing anti-viral responses across various organisms. ADARs are a group of proteins that act on RNA molecules and alter their sequence compared to the genes that encode them (a process termed RNA editing). In recent years, ADARs have been shown to suppress abnormal anti-viral responses triggered by self-components of the cell (RNA encoded by the cell). Here, we show that the involvement of ADARs in anti-viral response regulation is conserved in planarians (free-living flatworms). We identified two ADAR proteins in planarians and showed that eliminating one (ADAR1) results in animal death and that an anti-viral response commenced in the absence of either ADAR1 or ADAR2. We further identified one of the proteins (PRLR1) that participate in initiating this anti-viral response in planarians, which its mammalian homolog (MDA5) serves a similar role. Thus, our work suggests that ADARs involvement in suppressing aberrant anti-viral response is an ancient evolutionary invention and is likely shared across multicellular organisms with bilateral symmetry.

Diversity and evolution of the animal virome

The COVID-19 pandemic has given the study of virus evolution and ecology new relevance. Although viruses were first identified more than a century ago, we likely know less about their diversity than that of any other biological entity. Most documented animal viruses have been sampled from just two phyla - the Chordata and the Arthropoda - with a strong bias towards viruses that infect humans or animals of economic and social importance, often in association with strong disease phenotypes. Fortunately, the recent development of unbiased metagenomic next-generation sequencing is providing a richer view of the animal virome and shedding new light on virus evolution. In this Review, we explore our changing understanding of the diversity, composition and evolution of the animal virome. We outline the factors that determine the phylogenetic diversity and genomic structure of animal viruses on evolutionary timescales and show how this impacts assessment of the risk of disease emergence in the short term. We also describe the ongoing challenges in metagenomic analysis and outline key themes for future research. A central question is how major events in the evolutionary history of animals, such as the origin of the vertebrates and periodic mass extinction events, have shaped the diversity and evolution of the viruses they carry.

Toward a genome sequence for every animal: Where are we now?

In less than 25 y, the field of animal genome science has transformed from a discipline seeking its first glimpses into genome sequences across the Tree of Life to a global enterprise with ambitions to sequence genomes for all of Earth's eukaryotic diversity [H. A. Lewin et al., Proc. Natl. Acad. Sci. U.S.A. 115, 4325-4333 (2018)]. As the field rapidly moves forward, it is important to take stock of the progress that has been made to best inform the discipline's future. In this Perspective, we provide a contemporary, quantitative overview of animal genome sequencing. We identified the best available genome assemblies in GenBank, the world's most extensive genetic database, for 3,278 unique animal species across 24 phyla. We assessed taxonomic representation, assembly quality, and annotation status for major clades. We show that while tremendous taxonomic progress has occurred, stark disparities in genomic representation exist, highlighted by a systemic overrepresentation of vertebrates and underrepresentation of arthropods. In terms of assembly quality, long-read sequencing has dramatically improved contiguity, whereas gene annotations are available for just 34.3% of taxa. Furthermore, we show that animal genome science has diversified in recent years with an ever-expanding pool of researchers participating. However, the field still appears to be dominated by institutions in the Global North, which have been listed as the submitting institution for 77% of all assemblies. We conclude by offering recommendations for improving genomic resource availability and research value while also broadening global representation.

How many long branch orders occur in Chelicerata? Opposing effects of Palpigradi and Opilioacariformes on phylogenetic stability

Excepting a handful of nodes, phylogenetic relationships between chelicerate orders remains poorly resolved, due to both the incidence of long branch attraction artifacts and the limited sampling of key lineages. It has recently been shown that increasing representation of basal nodes plays an outsized role in resolving the higher-level placement of long branch chelicerate orders. Two lineages have been consistently undersampled in chelicerate phylogeny. First, sampling of the miniaturized order Palpigradi has been restricted to a fragmentary transcriptome of a single species. Second, sampling of Opilioacariformes, a rarely encountered and key group of Parasitiformes, has been restricted to a single exemplar. These two lineages exhibit dissimilar properties with respect to branch length; Opilioacariformes shows relatively low evolutionary rate compared to other Parasitiformes, whereas Palpigradi possibly acts as another long branch order (an effect that may be conflated with the degree of missing data). To assess these properties and their effects on tree stability, we constructed a phylogenomic dataset of Chelicerata wherein both lineages were sampled with three terminals, increasing the representation of these lineages per locus. We examined the effect of subsampling phylogenomic matrices using (1) taxon occupancy, (2) evolutionary rate, and (3) a principal components-based approach. We further explored the impact of taxon deletion experiments that mitigate the effect of long branches. Here, we show that Palpigradi constitutes a fourth long branch chelicerate order (together with Acariformes, Parasitiformes, and Pseudoscorpiones), which further destabilizes the chelicerate backbone topology. By contrast, the slow-evolving Opilioacariformes were consistently recovered within Parasitiformes, with certain subsampling practices recovering their placement as the sister group to the remaining Parasitiformes. Whereas the inclusion of Opilioacariformes always resulted in the non-monophyly of Acari with support, deletion of Opilioacariformes from datasets consistently incurred the monophyly of Acari except in matrices constructed on the basis of evolutionary rate. Our results strongly suggest that Acari is an artifact of long-branch attraction.

Phylogenomic analyses of mud dragons (Kinorhyncha)

Mud dragons (Kinorhyncha) are microscopic invertebrates, inhabiting marine sediments across the globe from intertidal to hadal depths. They are segmented, moulting animals like arthropods, but grouping with the unsegmented priapulans and loriciferans within Ecdysozoa. There are more than 300 species of kinorhynchs described within 31 genera and 11 families, however, their evolutionary relationships have so far only been investigated using morphology and a few molecular markers. Here we aim to resolve the relationships and classification of major clades within Kinorhyncha using transcriptomic data. In addition, we wish to revisit the position of three indistinctly segmented, aberrant genera in order to reconstruct the evolution of distinct segmentation within the group. We conducted a phylogenomic analysis of Kinorhyncha including 21 kinorhynch transcriptomes (of which 18 are new) representing 15 genera, and seven outgroups including priapulan, loriciferan, nematode and nematomorph transcriptomes. Results show a congruent and robust tree that supports the division of Kinorhyncha into two major clades: Cyclorhagida and Allomalorhagida. Cyclorhagida is composed of three subclades: Xenosomata, Kentrorhagata comb. nov. (including the aberrant Zelinkaderes) and Echinorhagata. Allomalorhagida is composed of two subclades: Pycnophyidae and Anomoirhaga nom. nov. Anomoirhaga nom. nov. accommodates the aberrant genera Cateria (previously nested within Cyclorhagida) and Franciscideres together with five additional genera. The distant and derived positions of the aberrant Zelinkaderes, Cateria and Franciscideres species suggest that their less distinct trunk segmentation evolved convergently, and that segmentation evolved among kinorhynch stem groups.

Phylogenomic Analysis of Velvet Worms (Onychophora) Uncovers an Evolutionary Radiation in the Neotropics

Onychophora ("velvet worms") are charismatic soil invertebrates known for their status as a "living fossil," their phylogenetic affiliation to arthropods, and their distinctive biogeographic patterns. However, several aspects of their internal phylogenetic relationships remain unresolved, limiting our understanding of the group's evolutionary history, particularly with regard to changes in reproductive mode and dispersal ability. To address these gaps, we used RNA sequencing and phylogenomic analysis of transcriptomes to reconstruct the evolutionary relationships and infer divergence times within the phylum. We recovered a fully resolved and well-supported phylogeny for the circum-Antarctic family Peripatopsidae, which retains signals of Gondwanan vicariance and showcases the evolutionary lability of reproductive mode in the family. Within the Neotropical clade of Peripatidae, though, we found that amino acid-translated sequence data masked nearly all phylogenetic signal, resulting in highly unstable and poorly supported relationships. Analyses using nucleotide sequence data were able to resolve many more relationships, though we still saw discordant phylogenetic signal between genes, probably indicative of a rapid, mid-Cretaceous radiation in the group. Finally, we hypothesize that the unique reproductive mode of placentotrophic viviparity found in all Neotropical peripatids may have facilitated the multiple inferred instances of over-water dispersal and establishment on oceanic islands.

Form and Function of the Vertebrate and Invertebrate Blood-Brain Barriers

The need to protect neural tissue from toxins or other substances is as old as neural tissue itself. Early recognition of this need has led to more than a century of investigation of the blood-brain barrier (BBB). Many aspects of this important neuroprotective barrier have now been well established, including its cellular architecture and barrier and transport functions. Unsurprisingly, most research has had a human orientation, using mammalian and other animal models to develop translational research findings. However, cell layers forming a barrier between vascular spaces and neural tissues are found broadly throughout the invertebrates as well as in all vertebrates. Unfortunately, previous scenarios for the evolution of the BBB typically adopt a classic, now discredited 'scala naturae' approach, which inaccurately describes a putative evolutionary progression of the mammalian BBB from simple invertebrates to mammals. In fact, BBB-like structures have evolved independently numerous times, complicating simplistic views of the evolution of the BBB as a linear process. Here, we review BBBs in their various forms in both invertebrates and vertebrates, with an emphasis on the function, evolution, and conditional relevance of popular animal models such as the fruit fly and the zebrafish to mammalian BBB research.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Origin and significance of two pairs of head tentacles in the radiation of euthyneuran sea slugs and land snails

The gastropod infraclass Euthyneura comprises at least 30,000 species of snails and slugs, including nudibranch sea slugs, sea hares and garden snails, that flourish in various environments on earth. A unique morphological feature of Euthyneura is the presence of two pairs of sensory head tentacles with different shapes and functions: the anterior labial tentacles and the posterior rhinophores or eyestalks. Here we combine molecular phylogenetic and microanatomical evidence that suggests the two pairs of head tentacles have originated by splitting of the original single tentacle pair (with two parallel nerve cords in each tentacle) as seen in many other gastropods. Minute deep-sea snails of Tjaernoeia and Parvaplustrum, which in our phylogeny belonged to the euthyneurans’ sister group (new infraclass Mesoneura), have tentacles that are split along much of their lengths but associated nerves and epidermal sense organs are not as specialized as in Euthyneura. We suggest that further elaboration of cephalic sense organs in Euthyneura closely coincided with their ecological radiation and drastic modification of body plans. The monotypic family Parvaplustridae nov., superfamily Tjaernoeioidea nov. (Tjaernoeiidae + Parvaplustridae), and new major clade Tetratentaculata nov. (Mesoneura nov. + Euthyneura) are also proposed based on their phylogenetic relationships and shared morphological traits.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Geographical sampling bias on the assessment of endemism areas for marine meiobenthic fauna

Species distribution patterns are constrained by historical and ecological processes in space and time, but very often the species range sizes are geographical sampling biases resulting from unequal sampling effort. One of the most common definitions of endemism is based on the "congruence of distributional areas" criterion, when two or more species have the same distributional limits. By acknowledging that available data of marine meiobenthic species are prone to geographical sampling bias and that can affect the accuracy of the biogeographical signals, the present study combines analyses of inventory incompleteness and recognition of spatial congruence of Gastrotricha, Kinorhyncha, meiobenthic Annelida and Tardigrada in order to better understand the large-scale distribution of these organisms in coastal and shelf areas of the world. We used the marine bioregionalization framework for geographical operative units to quantify the inventory incompleteness effect (by modelling spatial predictions of species richness) and to recognize areas of endemism. Our models showed that the difference between observed and expected species richness in the Southern Hemisphere is much higher than in the Northern Hemisphere. Parsimony Analysis of Endemicity delimited 20 areas of endemism, most found in the Northern Hemisphere. Distribution patterns of meiobenthic species are shown to respond to events of geographical barriers and abiotic features, and their distribution is far from homogeneous throughout the world. Also, our data show that ecoregions with distinct biotas have at least some cohesion over evolutionary time. However, we found that inventory incompleteness may significantly affect the explanatory power of areas of endemism delimitation in both hemispheres. Yet, whereas future increases in sampling efforts are likely to change the spatial congruence ranges in the Southern Hemisphere, patterns for the Northern Hemisphere may prove to be relatively more resilient.

A Likelihood-Ratio Test for Lumpability of Phylogenetic Data: Is the Markovian Property of an Evolutionary Process retained in Recoded DNA?

In molecular phylogenetics, it is typically assumed that the evolutionary process for DNA can be approximated by independent and identically distributed Markovian processes at the variable sites and that these processes diverge over the edges of a rooted bifurcating tree. Sometimes the nucleotides are transformed from a 4-state alphabet to a 3- or 2-state alphabet by a procedure that is called recoding, lumping, or grouping of states. Here, we introduce a likelihood-ratio test for lumpability for DNA that has diverged under different Markovian conditions, which assesses the assumption that the Markovian property of the evolutionary process over each edge is retained after recoding of the nucleotides. The test is derived and validated numerically on simulated data. To demonstrate the insights that can be gained by using the test, we assessed two published data sets, one of mitochondrial DNA from a phylogenetic study of the ratites (Syst. Biol. 59:90-107 [2010]) and the other of nuclear DNA from a phylogenetic study of yeast (Mol. Biol. Evol. 21:1455-1458 [2004]). Our analysis of these data sets revealed that recoding of the DNA eliminated some of the compositional heterogeneity detected over the sequences. However, the Markovian property of the original evolutionary process was not retained by the recoding, leading to some significant distortions of edge lengths in reconstructed trees.

Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web

High throughput sequencing and phylogenomic analyses focusing on relationships among spiders have both reinforced and upturned long-standing hypotheses. Likewise, the evolution of spider webs-perhaps their most emblematic attribute-is being understood in new ways. With a matrix including 272 spider species and close arachnid relatives, we analyze and evaluate the relationships among these lineages using a variety of orthology assessment methods, occupancy thresholds, tree inference methods and support metrics. Our analyses include families not previously sampled in transcriptomic analyses, such as Symphytognathidae, the only araneoid family absent in such prior works. We find support for the major established spider lineages, including Mygalomorphae, Araneomorphae, Synspermiata, Palpimanoidea, Araneoidea and the Retrolateral Tibial Apophysis Clade, as well as the uloborids, deinopids, oecobiids and hersiliids Grade. Resulting trees are evaluated using bootstrapping, Shimodaira-Hasegawa approximate likelihood ratio test, local posterior probabilities and concordance factors. Using structured Markov models to assess the evolution of spider webs while accounting for hierarchically nested traits, we find multiple convergent occurrences of the orb web across the spider tree-of-life. Overall, we provide the most comprehensive spider tree-of-life to date using transcriptomic data and use new methods to explore controversial issues of web evolution, including the origins and multiple losses of the orb web.

The enigmatic Placozoa part 1: Exploring evolutionary controversies and poor ecological knowledge

The placozoan Trichoplax adhaerens is a tiny hairy plate and more simply organized than any other living metazoan. After its original description by F.E. Schulze in 1883, it attracted attention as a potential model for the ancestral state of metazoan organization, the "Urmetazoon". Trichoplax lacks any kind of symmetry, organs, nerve cells, muscle cells, basal lamina, and extracellular matrix. Furthermore, the placozoan genome is the smallest (not secondarily reduced) genome of all metazoan genomes. It harbors a remarkably rich diversity of genes and has been considered the best living surrogate for a metazoan ancestor genome. The phylum Placozoa presently harbors three formally described species, while several dozen "cryptic" species are yet awaiting their description. The phylogenetic position of placozoans has recently become a contested arena for modern phylogenetic analyses and view-driven claims. Trichoplax offers unique prospects for understanding the minimal requirements of metazoan animal organization and their corresponding malfunctions.

Molecular and morphological analysis of the developing nemertean brain indicates convergent evolution of complex brains in Spiralia

BACKGROUND

The brain anatomy in the clade Spiralia can vary from simple, commissural brains (e.g., gastrotrichs, rotifers) to rather complex, partitioned structures (e.g., in cephalopods and annelids). How often and in which lineages complex brains evolved still remains unclear. Nemerteans are a clade of worm-like spiralians, which possess a complex central nervous system (CNS) with a prominent brain, and elaborated chemosensory and neuroglandular cerebral organs, which have been previously suggested as homologs to the annelid mushroom bodies. To understand the developmental and evolutionary origins of the complex brain in nemerteans and spiralians in general, we investigated details of the neuroanatomy and gene expression in the brain and cerebral organs of the juveniles of nemertean Lineus ruber.

RESULTS

In the juveniles, the CNS is already composed of all major elements present in the adults, including the brain, paired longitudinal lateral nerve cords, and an unpaired dorsal nerve cord, which suggests that further neural development is mostly related with increase in the size but not in complexity. The ultrastructure of the juvenile cerebral organ revealed that it is composed of several distinct cell types present also in the adults. The 12 transcription factors commonly used as brain cell type markers in bilaterians show region-specific expression in the nemertean brain and divide the entire organ into several molecularly distinct areas, partially overlapping with the morphological compartments. Additionally, several of the mushroom body-specific genes are expressed in the developing cerebral organs.

CONCLUSIONS

The dissimilar expression of molecular brain markers between L. ruber and the annelid Platynereis dumerilii indicates that the complex brains present in those two species evolved convergently by independent expansions of non-homologous regions of a simpler brain present in their last common ancestor. Although the same genes are expressed in mushroom bodies and cerebral organs, their spatial expression within organs shows apparent differences between annelids and nemerteans, indicating convergent recruitment of the same genes into patterning of non-homologous organs or hint toward a more complicated evolutionary process, in which conserved and novel cell types contribute to the non-homologous structures.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host

Background

Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensulato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea.

Results

Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria).

Conclusion

Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07857-y.

Cool your jets: biological jet propulsion in marine invertebrates

Pulsatile jet propulsion is a common swimming mode used by a diverse array of aquatic taxa from chordates to cnidarians. This mode of locomotion has interested both biologists and engineers for over a century. A central issue to understanding the important features of jet-propelling animals is to determine how the animal interacts with the surrounding fluid. Much of our knowledge of aquatic jet propulsion has come from simple theoretical approximations of both propulsive and resistive forces. Although these models and basic kinematic measurements have contributed greatly, they alone cannot provide the detailed information needed for a comprehensive, mechanistic overview of how jet propulsion functions across multiple taxa, size scales and through development. However, more recently, novel experimental tools such as high-speed 2D and 3D particle image velocimetry have permitted detailed quantification of the fluid dynamics of aquatic jet propulsion. Here, we provide a comparative analysis of a variety of parameters such as efficiency, kinematics and jet parameters, and review how they can aid our understanding of the principles of aquatic jet propulsion. Research on disparate taxa allows comparison of the similarities and differences between them and contributes to a more robust understanding of aquatic jet propulsion.

Diversity in the Development of the Neuromuscular System of Nemertean Larvae (Nemertea, Spiralia)

In studies on the development of nervous systems and musculature, fluorescent labeling of neuroactive substances and filamentous actin (f-actin) of muscle cells and the subsequent analysis with confocal laser scanning microscopy (CLSM), has led to a broad comparative data set for the majority of the clades of the superphylum Spiralia. However, a number of clades remain understudied, which results in gaps in our knowledge that drastically hamper the formulation of broad-scale hypotheses on the evolutionary developmental biology (EvoDevo) of the structures in question. Regarding comparative data on the development of the peptidergic nervous system and the musculature of species belonging to the spiralian clade Nemertea (ribbon worms), such considerable knowledge gaps are manifest. This paper presents first findings on fluorescent labeling of the FMRFamide-like component of the nervous system and contributes additional data on the muscle development in the presently still underrepresented larvae of palaeo- and hoplonemertean species. Whereas the architecture of the FMRFamide-like nervous system is comparably uniform between the studied representatives, the formation of the musculature differs considerably, exhibiting developmental modes yet undescribed for any spiralian species. The presented results fill a significant gap in the spiralian EvoDevo data set and thus allow for further elaboration of hypotheses on the ancestral pattern of the musculature and a prominent component of the nervous system in Nemertea. However, with respect to the variety observed, it is expected that the true diversity of the developmental pathways is still to be discovered when more detailed data on other nemertean species will be available.

Taxonomic Sampling and Rare Genomic Changes Overcome Long-Branch Attraction in the Phylogenetic Placement of Pseudoscorpions

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

FGF signaling acts on different levels of mesoderm development within Spiralia

FGF signaling is involved in mesoderm induction in members of deuterostomes (e.g. tunicates, hemichordates), but not in flies and nematodes, in which it has a role in mesoderm patterning and migration. However, we need comparable studies in other protostome taxa in order to decipher whether this mesoderm-inducing function of FGF extends beyond the lineage of deuterostomes. Here, we investigated the role of FGF signaling in mesoderm development in three species of lophophorates, a clade within the protostome group Spiralia. Our gene expression analyses show that the mesodermal molecular patterning is conserved between brachiopods and phoronids, but the spatial and temporal recruitment of transcription factors differs significantly. Moreover, the use of the inhibitor SU5402 demonstrates that FGF signaling is involved in different steps of mesoderm development, as well as in morphogenetic movements of gastrulation and axial elongation. Our findings suggest that the mesoderm-inducing role of FGF extends beyond the group of deuterostomes.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Neural versus alternative integrative systems: molecular insights into origins of neurotransmitters

Transmitter signalling is the universal chemical language of any nervous system, but little is known about its early evolution. Here, we summarize data about the distribution and functions of neurotransmitter systems in basal metazoans as well as outline hypotheses of their origins. We explore the scenario that neurons arose from genetically different populations of secretory cells capable of volume chemical transmission and integration of behaviours without canonical synapses. The closest representation of this primordial organization is currently found in Placozoa, disk-like animals with the simplest known cell composition but complex behaviours. We propose that injury-related signalling was the evolutionary predecessor for integrative functions of early transmitters such as nitric oxide, ATP, protons, glutamate and small peptides. By contrast, acetylcholine, dopamine, noradrenaline, octopamine, serotonin and histamine were recruited as canonical neurotransmitters relatively later in animal evolution, only in bilaterians. Ligand-gated ion channels often preceded the establishment of novel neurotransmitter systems. Moreover, lineage-specific diversification of neurotransmitter receptors occurred in parallel within Cnidaria and several bilaterian lineages, including acoels. In summary, ancestral diversification of secretory signal molecules provides unique chemical microenvironments for behaviour-driven innovations that pave the way to complex brain functions and elementary cognition. This article is part of the theme issue 'Basal cognition: multicellularity, neurons and the cognitive lens'.

The Roles of Protein Structure, Taxon Sampling, and Model Complexity in Phylogenomics: A Case Study Focused on Early Animal Divergences

Despite the long history of using protein sequences to infer the tree of life, the potential for different parts of protein structures to retain historical signal remains unclear. We propose that it might be possible to improve analyses of phylogenomic datasets by incorporating information about protein structure. We test this idea using the position of the root of Metazoa (animals) as a model system. We examined the distribution of “strongly decisive” sites (alignment positions that support a specific tree topology) in a dataset comprising >1500 proteins and almost 100 taxa. The proportion of each class of strongly decisive sites in different structural environments was very sensitive to the model used to analyze the data when a limited number of taxa were used but they were stable when taxa were added. As long as enough taxa were analyzed, sites in all structural environments supported the same topology regardless of whether standard tree searches or decisive sites were used to select the optimal tree. However, the use of decisive sites revealed a difference between the support for minority topologies for sites in different structural environments: buried sites and sites in sheet and coil environments exhibited equal support for the minority topologies, whereas solvent-exposed and helix sites had unequal numbers of sites, supporting the minority topologies. This suggests that the relatively slowly evolving buried, sheet, and coil sites are giving an accurate picture of the true species tree and the amount of conflict among gene trees. Taken as a whole, this study indicates that phylogenetic analyses using sites in different structural environments can yield different topologies for the deepest branches in the animal tree of life and that analyzing larger numbers of taxa eliminates this conflict. More broadly, our results highlight the desirability of incorporating information about protein structure into phylogenomic analyses.

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Lack of support for Deuterostomia prompts reinterpretation of the first Bilateria

The bilaterally symmetric animals (Bilateria) are considered to comprise two monophyletic groups, Protostomia (Ecdysozoa and the Lophotrochozoa) and Deuterostomia (Chordata and the Xenambulacraria). Recent molecular phylogenetic studies have not consistently supported deuterostome monophyly. Here, we compare support for Protostomia and Deuterostomia using multiple, independent phylogenomic datasets. As expected, Protostomia is always strongly supported, especially by longer and higher-quality genes. Support for Deuterostomia, however, is always equivocal and barely higher than support for paraphyletic alternatives. Conditions that cause tree reconstruction errors-inadequate models, short internal branches, faster evolving genes, and unequal branch lengths-coincide with support for monophyletic deuterostomes. Simulation experiments show that support for Deuterostomia could be explained by systematic error. The branch between bilaterian and deuterostome common ancestors is, at best, very short, supporting the idea that the bilaterian ancestor may have been deuterostome-like. Our findings have important implications for the understanding of early animal evolution.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.

Systematic errors in orthology inference and their effects on evolutionary analyses

The availability of complete sets of genes from many organisms makes it possible to identify genes unique to (or lost from) certain clades. This information is used to reconstruct phylogenetic trees; identify genes involved in the evolution of clade specific novelties; and for phylostratigraphy—identifying ages of genes in a given species. These investigations rely on accurately predicted orthologs. Here we use simulation to produce sets of orthologs that experience no gains or losses. We show that errors in identifying orthologs increase with higher rates of evolution. We use the predicted sets of orthologs, with errors, to reconstruct phylogenetic trees; to count gains and losses; and for phylostratigraphy. Our simulated data, containing information only from errors in orthology prediction, closely recapitulate findings from empirical data. We suggest published downstream analyses must be informed to a large extent by errors in orthology prediction that mimic expected patterns of gene evolution.

•

Presence of shared orthologs across species is used for evolutionary analyses

•

We simulated realistic sets of orthologs with no gains or losses

•

Errors predicting shared orthologs correlate with phylogenetic relationships

•

Presence/absence datasets based on errors recapitulate findings from empirical data

PhyKIT: a broadly applicable UNIX shell toolkit for processing and analyzing phylogenomic data

MOTIVATION

Diverse disciplines in biology process and analyze multiple sequence alignments (MSAs) and phylogenetic trees to evaluate their information content, infer evolutionary events and processes, and predict gene function. However, automated processing of MSAs and trees remains a challenge due to the lack of a unified toolkit. To fill this gap, we introduce PhyKIT, a toolkit for the UNIX shell environment with 30 functions that process MSAs and trees, including but not limited to estimation of mutation rate, evaluation of sequence composition biases, calculation of the degree of violation of a molecular clock, and collapsing bipartitions (internal branches) with low support.

RESULTS

To demonstrate the utility of PhyKIT, we detail three use cases: (1) summarizing information content in MSAs and phylogenetic trees for diagnosing potential biases in sequence or tree data; (2) evaluating gene-gene covariation of evolutionary rates to identify functional relationships, including novel ones, among genes; and (3) identify lack of resolution events or polytomies in phylogenetic trees, which are suggestive of rapid radiation events or lack of data. We anticipate PhyKIT will be useful for processing, examining, and deriving biological meaning from increasingly large phylogenomic datasets.

AVAILABILITY

PhyKIT is freely available on GitHub (https://github.com/JLSteenwyk/PhyKIT), PyPi (https://pypi.org/project/phykit/), and the Anaconda Cloud (https://anaconda.org/JLSteenwyk/phykit) under the MIT license with extensive documentation and user tutorials (https://jlsteenwyk.com/PhyKIT).

SUPPLEMENTARY INFORMATION

Supplementary data are available on figshare (doi: 10.6084/m9.figshare.13118600) and are available at Bioinformatics online.

Evolutionary dynamics of transposable elements in bdelloid rotifers

Transposable elements (TEs) are selfish genomic parasites whose ability to spread autonomously is facilitated by sexual reproduction in their hosts. If hosts become obligately asexual, TE frequencies and dynamics are predicted to change dramatically, but the long-term outcome is unclear. Here, we test current theory using whole-genome sequence data from eight species of bdelloid rotifers, a class of invertebrates in which males are thus far unknown. Contrary to expectations, we find a variety of active TEs in bdelloid genomes, at an overall frequency within the range seen in sexual species. We find no evidence that TEs are spread by cryptic recombination or restrained by unusual DNA repair mechanisms. Instead, we find that that TE content evolves relatively slowly in bdelloids and that gene families involved in RNAi-mediated TE suppression have undergone significant expansion, which might mitigate the deleterious effects of active TEs and compensate for the consequences of long-term asexuality.

Evolution of Protein-Mediated Biomineralization in Scleractinian Corals

While recent strides have been made in understanding the biological process by which stony corals calcify, much remains to be revealed, including the ubiquity across taxa of specific biomolecules involved. Several proteins associated with this process have been identified through proteomic profiling of the skeletal organic matrix (SOM) extracted from three scleractinian species. However, the evolutionary history of this putative “biomineralization toolkit,” including the appearance of these proteins’ throughout metazoan evolution, remains to be resolved. Here we used a phylogenetic approach to examine the evolution of the known scleractinians’ SOM proteins across the Metazoa. Our analysis reveals an evolutionary process dominated by the co-option of genes that originated before the cnidarian diversification. Each one of the three species appears to express a unique set of the more ancient genes, representing the independent co-option of SOM proteins, as well as a substantial proportion of proteins that evolved independently. In addition, in some instances, the different species expressed multiple orthologous proteins sharing the same evolutionary history. Furthermore, the non-random clustering of multiple SOM proteins within scleractinian-specific branches suggests the conservation of protein function between distinct species for what we posit is part of the scleractinian “core biomineralization toolkit.” This “core set” contains proteins that are likely fundamental to the scleractinian biomineralization mechanism. From this analysis, we infer that the scleractinians’ ability to calcify was achieved primarily through multiple lineage-specific protein expansions, which resulted in a new functional role that was not present in the parent gene.

Nuclear Receptors and Development of Marine Invertebrates

Nuclear Receptors (NRs) are a superfamily of transcription factors specific to metazoans that have the unique ability to directly translate the message of a signaling molecule into a transcriptional response. In vertebrates, NRs are pivotal players in countless processes of both embryonic and adult physiology, with embryonic development being one of the most dynamic periods of NR activity. Accumulating evidence suggests that NR signaling is also a major regulator of development in marine invertebrates, although ligands and transactivation dynamics are not necessarily conserved with respect to vertebrates. The explosion of genome sequencing projects and the interpretation of the resulting data in a phylogenetic context allowed significant progress toward an understanding of NR superfamily evolution, both in terms of molecular activities and developmental functions. In this context, marine invertebrates have been crucial for characterizing the ancestral states of NR-ligand interactions, further strengthening the importance of these organisms in the field of evolutionary developmental biology.

A molecular phylogeny of historical and contemporary specimens of an under-studied micro-invertebrate group

Resolution of relationships at lower taxonomic levels is crucial for answering many evolutionary questions, and as such, sufficiently varied species representation is vital. This latter goal is not always achievable with relatively fresh samples. To alleviate the difficulties in procuring rarer taxa, we have seen increasing utilization of historical specimens in building molecular phylogenies using high throughput sequencing. This effort, however, has mainly focused on large-bodied or well-studied groups, with small-bodied and under-studied taxa under-prioritized. Here, we utilize both historical and contemporary specimens, to increase the resolution of phylogenetic relationships among a group of under-studied and small-bodied metazoans, namely, cheilostome bryozoans. In this study, we pioneer the sequencing of air-dried cheilostomes, utilizing a recently developed library preparation method for low DNA input. We evaluate a de novo mitogenome assembly and two iterative methods, using the sequenced target specimen as a reference for mapping, for our sequences. In doing so, we present mitochondrial and ribosomal RNA sequences of 43 cheilostomes representing 37 species, including 14 from historical samples ranging from 50 to 149 years old. The inferred phylogenetic relationships of these samples, analyzed together with publicly available sequence data, are shown in a statistically well-supported 65 taxa and 17 genes cheilostome tree, which is also the most broadly sampled and largest to date. The robust phylogenetic placement of historical samples whose contemporary conspecifics and/or congenerics have been sequenced verifies the appropriateness of our workflow and gives confidence in the phylogenetic placement of those historical samples for which there are no close relatives sequenced. The success of our workflow is highlighted by the circularization of a total of 27 mitogenomes, seven from historical cheilostome samples. Our study highlights the potential of utilizing DNA from micro-invertebrate specimens stored in natural history collections for resolving phylogenetic relationships among species.