Platyzoan paraphyly based on phylogenomic data supports a noncoelomate ancestry of spiralia

A Jurassic acanthocephalan illuminates the origin of thorny-headed worms

Acanthocephala (thorny-headed worms), characterized by the presence of an eversible proboscis with hooks, are a diverse endoparasitic group that infect a wide range of vertebrates and invertebrates1. Although long regarded as a separate phylum, they have several putative sister taxa based on morphological features, including Platyhelminthes (flatworms)2, Priapulida (penis worms)3 and Rotifera (wheel animals)4. Molecular phylogenies have instead recovered them within rotifers5-10, suggesting acanthocephalans are derived from free-living worms with a jaw apparatus (Gnathifera). Their only fossil record is Late Cretaceous eggs11, contributing limited palaeontological information to deciphering their early evolution. Here we describe an acanthocephalan body fossil, Juracanthocephalus daohugouensis gen. et. sp. nov., from the Middle Jurassic Daohugou biota of China. Juracanthocephalus shows unambiguous acanthocephalan characteristics, for example a hooked proboscis, a bursa, as well as a jaw apparatus with discrete elements that is typical of other gnathiferans. Juracanthocephalus shares features with Seisonidea (an epizoic member of Rotifera) and Acanthocephala, bridging the evolutionary gap between jawed rotifers and the obligate parasitic, jawless acanthocephalans. Our results reveal previously unrecognized ecological and morphological diversity in ancient Acanthocephala and highlight the significance of transitional fossils, revealing the origins of this highly enigmatic group of living organisms.

A Phylogenomic Backbone for Acoelomorpha Inferred From Transcriptomic Data

Xenacoelomorpha are mostly microscopic, morphologically simple worms, lacking many structures typical of other bilaterians. Xenacoelomorphs-which include three main groups, namely Acoela, Nemertodermatida, and Xenoturbella-have been proposed to be an early diverging Bilateria, sister to protostomes and deuterostomes, but other phylogenomic analyses have recovered this clade nested within the deuterostomes, as sister to Ambulacraria. The position of Xenacoelomorpha within the metazoan tree has understandably attracted a lot of attention, overshadowing the study of phylogenetic relationships within this group. Given that Xenoturbella includes only six species whose relationships are well understood, we decided to focus on the most speciose Acoelomorpha (Acoela + Nemertodermatida). Here, we have sequenced 29 transcriptomes, doubling the number of sequenced species, to infer a backbone tree for Acoelomorpha based on genomic data. The recovered topology is mostly congruent with previous studies. The most important difference is the recovery of Paratomella as the first off-shoot within Acoela, dramatically changing the reconstruction of the ancestral acoel. Besides, we have detected incongruence between the gene trees and the species tree, likely linked to incomplete lineage sorting, and some signal of introgression between the families Dakuidae and Mecynostomidae, which hampers inferring the correct placement of this family and, particularly, of the genus Notocelis. We have also used this dataset to infer for the first time diversification times within Acoelomorpha, which coincide with known bilaterian diversification and extinction events. Given the importance of morphological data in acoelomorph phylogenetics, we tested several partitions and models. Although morphological data failed to recover a robust phylogeny, phylogenetic placement has proven to be a suitable alternative when a reference phylogeny is available.

Multiple Displacement Amplification Facilitates SMRT Sequencing of Microscopic Animals and the Genome of the Gastrotrich Lepidodermella squamata (Dujardin 1841)

Obtaining adequate DNA for long-read genome sequencing remains a roadblock to producing contiguous genomes from small-bodied organisms, hindering understanding of phylogenetic relationships and genome evolution. Multiple displacement amplification leverages Phi29 DNA polymerase to produce micrograms of DNA from picograms of input. However, multiple displacement amplification's inherent biases in amplification related to guanine and cytosine (GC) content, repeat content and chimera production are a problem for long-read genome assembly, which has been little investigated. We explored the utility of multiple displacement amplification for generating template DNA for High Fidelity (HiFi) sequencing directly from living cells of Caenorhabditis elegans (Nematoda) and Lepidodermella squamata (Gastrotricha) containing one order of magnitude less DNA than required for the PacBio Ultra-Low DNA Input Workflow. High Fidelity sequencing of libraries prepared from multiple displacement amplification products resulted in highly contiguous and complete genomes for both C. elegans (102 Mbp assembly; 336 contigs; N50 = 868 kbp; L50 = 39; BUSCO_nematoda_nucleotide: S:96.1%, D:2.8%) and L. squamata (122 Mbp assembly; 157 contigs; N50 = 3.9 Mbp; L50 = 13; BUSCO_metazoa_nucleotide: S:80.8%, D:2.8%). Coverage uniformity for reads from multiple displacement amplification DNA (Gini Index: 0.14, normalized mean across all 100 kbp blocks: 0.49) and reads from pooled nematode DNA (Gini Index: 0.16, normalized mean across all 100 kbp blocks: 0.49) proved similar. Using this approach, we sequenced the genome of the microscopic invertebrate L. squamata (Gastrotricha), the first of its phylum. Using the newly sequenced genome, we infer Gastrotricha's long-debated phylogenetic position as the sister taxon of Platyhelminthes and conduct a comparative analysis of the Hox cluster.

Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

First Chromosome-Level Genome Assembly of a Ribbon Worm from the Hoplonemertea Clade, Emplectonema gracile, and Its Structural Annotation

Genome-wide information has so far been unavailable for ribbon worms of the clade Hoplonemertea, the most species-rich class within the phylum Nemertea. While species within Pilidiophora, the sister clade of Hoplonemertea, possess a pilidium larval stage and lack stylets on their proboscis, Hoplonemertea species have a planuliform larva and are armed with stylets employed for the injection of toxins into their prey. To further compare these developmental, physiological, and behavioral differences from a genomic perspective, the availability of a reference genome for a Hoplonemertea species is crucial. Such data will be highly useful for future investigations toward a better understanding of molecular ecology, venom evolution, and regeneration not only in Nemertea but also in other marine invertebrate phyla. To this end, we herein present the annotated chromosome-level genome assembly for Emplectonema gracile (Nemertea; Hoplonemertea; Monostilifera; Emplectonematidae), an easily collected nemertean well suited for laboratory experimentation. The genome has an assembly size of 157.9 Mb. Hi-C scaffolding yielded chromosome-level scaffolds, with a scaffold N50 of 10.0 Mb and a score of 95.1% for complete BUSCO genes found as a single copy. Annotation predicted 20,684 protein-coding genes. The high-quality reference genome reaches an Earth BioGenome standard level of 7.C.Q50.

The Evolution and Characterization of the RNA Interference Pathways in Lophotrochozoa

In animals, three main RNA interference mechanisms have been described so far, which respectively maturate three types of small noncoding RNAs (sncRNAs): miRNAs, piRNAs, and endo-siRNAs. The diversification of these mechanisms is deeply linked with the evolution of the Argonaute gene superfamily since each type of sncRNA is typically loaded by a specific Argonaute homolog. Moreover, other protein families play pivotal roles in the maturation of sncRNAs, like the DICER ribonuclease family, whose DICER1 and DICER2 paralogs maturate respectively miRNAs and endo-siRNAs. Within Metazoa, the distribution of these families has been only studied in major groups, and there are very few data for clades like Lophotrochozoa. Thus, we here inferred the evolutionary history of the animal Argonaute and DICER families including 43 lophotrochozoan species. Phylogenetic analyses along with newly sequenced sncRNA libraries suggested that in all Trochozoa, the proteins related to the endo-siRNA pathway have been lost, a part of them in some phyla (i.e. Nemertea, Bryozoa, Entoprocta), while all of them in all the others. On the contrary, early diverging phyla, Platyhelminthes and Syndermata, showed a complete endo-siRNA pathway. On the other hand, miRNAs were revealed the most conserved and ubiquitous mechanism of the metazoan RNA interference machinery, confirming their pivotal role in animal cell regulation.

A giant stem-group chaetognath

Chaetognaths, with their characteristic grasping spines, are the oldest known pelagic predators, found in the lowest Cambrian (Terreneuvian). Here, we describe a large stem chaetognath, Timorebestia koprii gen. et sp. nov., from the lower Cambrian Sirius Passet Lagerstätte, which exhibits lateral and caudal fins, a distinct head region with long antennae and a jaw apparatus similar to Amiskwia sagittiformis. Amiskwia has previously been interpreted as a total-group chaetognathiferan, as either a stem-chaetognath or gnathostomulid. We show that T. koprii shares a ventral ganglion with chaetognaths to the exclusion of other animal groups, firmly placing these fossils on the chaetognath stem. The large size (up to 30 cm) and gut contents in T. koprii suggest that early chaetognaths occupied a higher trophic position in pelagic food chains than today.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Fast cycling culture of the annelid model Platynereis dumerilii

Platynereis dumerilii, a marine annelid, is a model animal that has gained popularity in various fields such as developmental biology, biological rhythms, nervous system organization and physiology, behaviour, reproductive biology, and epigenetic regulation. The transparency of P. dumerilii tissues at all developmental stages makes it easy to perform live microscopic imaging of all cell types. In addition, the slow-evolving genome of P. dumerilii and its phylogenetic position as a representative of the vast branch of Lophotrochozoans add to its evolutionary significance. Although P. dumerilii is amenable to transgenesis and CRISPR-Cas9 knockouts, its relatively long and indefinite life cycle, as well as its semelparous reproduction have been hindrances to its adoption as a reverse genetics model. To overcome this limitation, an adapted culturing method has been developed allowing much faster life cycling, with median reproductive age at 13-14 weeks instead of 25-35 weeks using the traditional protocol. A low worm density in boxes and a strictly controlled feeding regime are important factors for the rapid growth and health of the worms. This culture method has several advantages, such as being much more compact, not requiring air bubbling or an artificial moonlight regime for synchronized sexual maturation and necessitating only limited water change. A full protocol for worm care and handling is provided.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

Single-cell transcriptomics refuels the exploration of spiralian biology

Spiralians represent the least studied superclade of bilaterian animals, despite exhibiting the widest diversity of organisms. Although spiralians include iconic organisms, such as octopus, earthworms and clams, a lot remains to be discovered regarding their phylogeny and biology. Here, we review recent attempts to apply single-cell transcriptomics, a new pioneering technology enabling the classification of cell types and the characterisation of their gene expression profiles, to several spiralian taxa. We discuss the methodological challenges and requirements for applying this approach to marine organisms and explore the insights that can be brought by such studies, both from a biomedical and evolutionary perspective. For instance, we show that single-cell sequencing might help solve the riddle of the homology of larval forms across spiralians, but also to better characterise and compare the processes of regeneration across taxa. We highlight the capacity of single-cell to investigate the origin of evolutionary novelties, as the mollusc shell or the cephalopod visual system, but also to interrogate the conservation of the molecular fingerprint of cell types at long evolutionary distances. We hope that single-cell sequencing will open a new window in understanding the biology of spiralians, and help renew the interest for these overlooked but captivating organisms.

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.

CaaX-less lamins: Lophotrochozoa provide a glance at the playground of evolution

Nuclear lamins are the main components of the nuclear lamina in many eukaryotes. They are members of the intermediate filament (IF) protein family. Lamins differ from cytoplasmic IF proteins by the presence of a nuclear localisation sequence (NLS) and a C-terminal tetrapeptide, the CaaX motif. The CaaX motif is target of post-translational modifications including isoprenylation, proteolytic processing, and carboxyl-methylation. These modifications, in conjunction with the NLS, direct lamins to the inner nuclear membrane where they assemble into filaments. Lamins lacking a CaaX motif are unable to associate independently with nuclear membranes and remain in the nucleoplasm. So far, three species have been reported to exclusively express CaaX-less lamins. All three belong to the lophotrochozoan lineage. To find out whether they represent rare exceptions, we analysed lamins of representatives of 17 lophotrochozoan phyla. Here we report that all four clades of Rotifera as well as individual taxa of Mollusca and Annelida lack CaaX-lamins, but express lamins with alternative C-termini. Of note, the respective mollusc and annelid groups occupy very different phylogenetic ranks. Most of these alternative C-termini are rich in aromatic residues. A possible function of these residues in membrane association is discussed. Alternative splicing of terebellid lamin transcripts gives rise to two lamin variants, one with a CaaX motif and one with an alternative C-terminus. A similar situation is found in Arenicolidae, Opheliidae, Capitellidae, and Echiura. This points a way, how the switch from lamins carrying a CaaX motif to lamins with alternative C-termini may have occurred.

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 .

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.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

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.

Evolution, expression patterns and distribution of novel ribbon worm predatory and defensive toxins

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

The evolution of the metazoan Toll receptor family and its expression during protostome development

BACKGROUND

Toll-like receptors (TLRs) play a crucial role in immunity and development. They contain leucine-rich repeat domains, one transmembrane domain, and one Toll/IL-1 receptor domain. TLRs have been classified into V-type/scc and P-type/mcc TLRs, based on differences in the leucine-rich repeat domain region. Although TLRs are widespread in animals, detailed phylogenetic studies of this gene family are lacking. Here we aim to uncover TLR evolution by conducting a survey and a phylogenetic analysis in species across Bilateria. To discriminate between their role in development and immunity we furthermore analyzed stage-specific transcriptomes of the ecdysozoans Priapulus caudatus and Hypsibius exemplaris, and the spiralians Crassostrea gigas and Terebratalia transversa.

RESULTS

We detected a low number of TLRs in ecdysozoan species, and multiple independent radiations within the Spiralia. V-type/scc and P-type/mcc type-receptors are present in cnidarians, protostomes and deuterostomes, and therefore they emerged early in TLR evolution, followed by a loss in xenacoelomorphs. Our phylogenetic analysis shows that TLRs cluster into three major clades: clade α is present in cnidarians, ecdysozoans, and spiralians; clade β in deuterostomes, ecdysozoans, and spiralians; and clade γ is only found in spiralians. Our stage-specific transcriptome and in situ hybridization analyses show that TLRs are expressed during development in all species analyzed, which indicates a broad role of TLRs during animal development.

CONCLUSIONS

Our findings suggest that a clade α TLR gene (TLR-Ca) and a clade β/γ TLR gene (TLR-Cβ/γ) were already present in the cnidarian-bilaterian common ancestor. However, although TLR-Ca was conserved in cnidarians, TLR-Cβ/γ was lost during the early evolution of these taxa. Moreover, TLR-Cβ/γ duplicated to generate TLR-Cβ and TLR-Cγ in the lineage to the last common protostome-deuterostome ancestor. TLR-Ca, TLR-Cβ and TLR-Cγ further expanded generating the three major TLR clades. While all three clades radiated in several spiralian lineages, specific TLRs clades have been presumably lost in other lineages. Furthermore, the expression of the majority of these genes during protostome ontogeny suggests a likely role in development.

Sweet or salty? The origin of freshwater gastrotrichs (Gastrotricha, Chaetonotida) revealed by molecular phylogenetic analysis

Chaetonotidae is the most diverse and widely distributed family of the order Chaetonotida (Gastrotricha) and includes both marine and freshwater species. Although the family is regarded as a sister taxon to the exclusively marine Xenotrichulidae, the type of environment, marine or freshwater, where Chaetonotidae originated is still not known. Here, we reconstructed the phylogeny of the family based on molecular sequence data and mapped both morphological and ecological characters to determine the ancestral environment of the first members of the family. Our results revealed that the freshwater genus Bifidochaetus is the earliest branching lineage in the paraphyletic Chaetonotidae (encompassing Dasydytidae and Neogosseidae). Moreover, we reconstructed Lepidochaetus-Cephalionotus clade as a monophyletic sister group to the remaining chaetonotids, which supports Kisielewski's morphological based hypothesis concerning undifferentiated type of body scales as a most primary character in Chaetonotidae. We also found that reversals to marine habitats occurred independently in different Chaetonotidae lineages, thus marine species in the genera Heterolepidoderma, Halichaetonotus, Aspidiophorus and subgenera Chaetonotus (Schizochaetonotus) or Chaetonotus (Marinochaetus) should be assumed as having secondarily invaded the marine environment. Character mapping revealed a series of synapomorphies that define the clade that includes Chaetonotidae (with Dasydytidae and Neogosseidae), the most important of which may be those linked to reproduction.

The parasite Schistocephalus solidus secretes proteins with putative host manipulation functions

BACKGROUND

Manipulative parasites are thought to liberate molecules in their external environment, acting as manipulation factors with biological functions implicated in their host's physiological and behavioural alterations. These manipulation factors are part of a complex mixture called the secretome. While the secretomes of various parasites have been described, there is very little data for a putative manipulative parasite. It is necessary to study the molecular interaction between a manipulative parasite and its host to better understand how such alterations evolve.

METHODS

Here, we used proteomics to characterize the secretome of a model cestode with a complex life cycle based on trophic transmission. We studied Schistocephalus solidus during the life stage in which behavioural changes take place in its obligatory intermediate fish host, the threespine stickleback (Gasterosteus aculeatus). We produced a novel genome sequence and assembly of S. solidus to improve protein coding gene prediction and annotation for this parasite. We then described the whole worm's proteome and its secretome during fish host infection using LC-MS/MS.

RESULTS

A total of 2290 proteins were detected in the proteome of S. solidus, and 30 additional proteins were detected specifically in the secretome. We found that the secretome contains proteases, proteins with neural and immune functions, as well as proteins involved in cell communication. We detected receptor-type tyrosine-protein phosphatases, which were reported in other parasitic systems to be manipulation factors. We also detected 12 S. solidus-specific proteins in the secretome that may play important roles in host-parasite interactions.

CONCLUSIONS

Our results suggest that S. solidus liberates molecules with putative host manipulation functions in the host and that many of them are species-specific.

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.

18S rRNA variability maps reveal three highly divergent, conserved motifs within Rotifera

BACKGROUND

18S rRNA is a major component of the small subunit of the eukaryotic ribosome and an important phylogenetic marker for many groups, often to the point of being the only marker available for some. A core structure across eukaryotes exists for this molecule that can help to inform about its evolution in different groups. Using an alignment of 18S rDNA for Rotifera as traditionally recognized (=Bdelloidea, Monogononta, and Seisonacea, but not Acanthocephala), I fitted sequences for three exemplar species (Adineta vaga, Brachionus plicatilis, and Seison nebaliae, respectively) to the core structure and used these maps to reveal patterns of evolution for the remainder of this diverse group of microscopic animals.

RESULTS

The obtained variability maps of the 18S rRNA molecule revealed a pattern of high diversity among the three major rotifer clades coupled with strong conservation within each of bdelloids and monogononts. A majority of individual sites (ca. 60%) were constant even across rotifers as a whole with variable sites showing only intermediate rates of evolution. Although the three structural maps each showed good agreement with the inferred core structure for eukaryotic 18S rRNA and so were highly similar to one another at the secondary and tertiary levels, the overall pattern is of three highly distinct, but conserved motifs within the group at the primary sequence level. A novel finding was that of a variably expressed deletion at the 3' end of the V3 hypervariable region among some bdelloid species that occasionally extended into and included the pseudoknot structure following this region as well as the central "square" of the 18S rRNA molecule. Compared to other groups, levels of variation and rates of evolution for 18S rRNA in Rotifera roughly matched those for Gastropoda and Acanthocephala, despite increasing evidence for the latter being a clade within Rotifera.

CONCLUSIONS

The lack of comparative data for comparable groups makes interpretation of the results (i.e., very low variation within each of the three major rotifer clades, but high variation between them) and their potential novelty difficult. However, these findings in combination with the high morphological diversity within rotifers potentially help to explain why no clear consensus has been reached to date with regard to the phylogenetic relationships among the major groups.

Expression of Piwi Genes during the Regeneration of Lineus sanguineus (Nemertea, Pilidiophora, Heteronemertea)

The transposon silencer piwi genes play important roles in germline determination and maintenance, gametogenesis, and stem-cell self-renewal, and the expression of certain piwi genes is indispensable for regeneration. Knowledge about piwi genes is needed for phylum Nemertea, which contains members (e.g., Lineus sanguineus) with formidable regeneration capacity. By searching the L. sanguineus genome, we identified six Argonaute genes including three ago (Ls-Ago2, Ls-Ago2a, and Ls-Ago2b) and three piwi (Ls-piwi1, Ls-piwi2, and Ls-piwi3) genes. In situ hybridization revealed that, in intact females, Ls-piwi2 and Ls-piwi3 were not expressed, while Ls-piwi1 was expressed in ovaries. During regeneration, Ls-piwi1 and Ls-pcna (proliferating cell nuclear antigen) had strong and similar expressions. The expression of Ls-piwi1 became indetectable while Ls-pcna continued to be expressed when the differentiation of new organs was finished. During anterior regeneration, expression signals of Ls-piwi2 and Ls-piwi3 were weak and only detected in the blastema stage. During posterior regeneration, no expression was observed for Ls-piwi2. To date, no direct evidence has been found for the existence of congenital stem cells in adult L. sanguineus. The "pluripotent cells" in regenerating tissues are likely to be dedifferentiated from other type(s) of cells.

Discovery of two ‘chimeric’ Gastrotricha and their systematic placement based on an integrative approach

Sublittoral sand from the islands of Sardinia (Italy) and Flores (Azores) – separated by more than 3700 km linear distance and 8 years between two independent sampling campaigns – yielded conspicuous specimens of two bizarre, yet undescribed, species of the marine gastrotrich clade Macrodasyida. These gastrotrichs combine several character traits that were already known from two, non-related genera. Morphological data were carefully analysed and digitally documented, and nuclear and mitochondrial DNA sequences were used for phylogenetic inference. The results of these analyses claim for the erection of a new genus. Specimens of the new taxon have a body length of less than 400 µm and are characterized by a wide, funnel-shaped mouth opening shielded dorsally by an oral hood and possess a posterior peduncle that ends with a Y-shaped pair of appendages that carry the posterior adhesive tubes. Further tubes occur as anterior, ventrolateral and lateral series; the gonads are unpaired and there is a set of two accessory reproductive organs. Molecular phylogenetic analyses confirm the results of former studies and clearly place the new taxon in Thaumastodermatidae. We hereby propose the establishment of Chimaeradasys gen. nov. and describe C. oligotubulatus sp. nov. from the Azores and C. polytubulatus sp. nov. from Sardinia.

A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

The curious and neglected soft-bodied meiofauna: Rouphozoa (Gastrotricha and Platyhelminthes)

Gastrotricha and Platyhelminthes form a clade called Rouphozoa. Representatives of both taxa are main components of meiofaunal communities, but their role in the trophic ecology of marine and freshwater communities is not sufficiently studied. Traditional collection methods for meiofauna are optimized for Ecdysozoa, and include the use of fixatives or flotation techniques that are unsuitable for the preservation and identification of soft-bodied meiofauna. As a result, rouphozoans are usually underestimated in conventional biodiversity surveys and ecological studies. Here, we give an updated outline of their diversity and taxonomy, with some phylogenetic considerations. We describe successfully tested techniques for their recovery and study, and emphasize current knowledge on the ecology, distribution and dispersal of freshwater gastrotrichs and microturbellarians. We also discuss the opportunities and pitfalls of (meta)barcoding studies as a means of overcoming the taxonomic impediment. Finally, we discuss the importance of rouphozoans in aquatic ecosystems and provide future research directions to fill in crucial gaps in the biology of these organisms needed for understanding their basic role in the ecology of benthos and their place in the trophic networks linking micro-, meio- and macrofauna of freshwater ecosystems.

Within-family plasticity of nervous system architecture in Syllidae (Annelida, Errantia)

BACKGROUND

The ground pattern underlying the nervous system of the last common ancestor in annelids was long thought to be settled, consisting of a dorsal brain, circumoesophageal connectives and a subepithelial, ladder-like ventral nerve cord with segmental ganglia connected by paired connectives. With the advent of immunocytochemical stainings and confocal laser scanning microscopy, it becomes evident that its architecture is extremely diverse, which makes the reconstruction of a ground pattern in annelida challenging. Whereas the nervous systems of many different families has already been described, only very few studies looked at the diversity of nervous systems within such clades to give a closer estimate on how plastic the annelid nervous system really is. So far, little is known on syllid nervous system architecture, one of the largest and most diverse groups of marine annelids.

RESULTS

The position of the brain, the circumoesophageal connectives, the stomatogastric nervous system, the longitudinal nerves that traverse each segment and the innervation of appendages are relatively uniform within the clade. Both the number of connectives within the ventral nerve cord and the number of segmental nerves, which in earlier studies were used to infer phylogenetic relationships and to reconstruct an annelid ground pattern, are highly diverse and differ between genera or even within a given genus. Differences in the distribution of somata of the brain, the nuchal innervation and its associated cell bodies were found between Syllinae and Exogoninae and may be subfamily-specific.

CONCLUSIONS

The nervous system morphology of syllids very likely depends on the taxon-specific ecological requirements. Thus, it is not surprising that in a clade, which occupies such diverse niches as the Annelida, we find similar patterns in phylogenetically widely separated species in similar niches and a high degree of modularity within a family. Only standardized protocols and staining methods can lead to comparable results, but so far different approaches have been taken to describe annelid nervous systems, making homologization of certain structures difficult. This study provides the first thorough description of the nervous system in the family Syllidae, allowing more detailed comparisons between annelid families in the future.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Discovery of the oldest bilaterian from the Ediacaran of South Australia

Analysis of modern animals and Ediacaran trace fossils predicts that the oldest bilaterians were simple and small. Such organisms would be difficult to recognize in the fossil record, but should have been part of the Ediacara Biota, the earliest preserved macroscopic, complex animal communities. Here, we describe Ikaria wariootia gen. et sp. nov. from the Ediacara Member, South Australia, a small, simple organism with anterior/posterior differentiation. We find that the size and morphology of Ikaria match predictions for the progenitor of the trace fossil Helminthoidichnites-indicative of mobility and sediment displacement. In the Ediacara Member, Helminthoidichnites occurs stratigraphically below classic Ediacara body fossils. Together, these suggest that Ikaria represents one of the oldest total group bilaterians identified from South Australia, with little deviation from the characters predicted for their last common ancestor. Further, these trace fossils persist into the Phanerozoic, providing a critical link between Ediacaran and Cambrian animals.

A developmental perspective on the evolution of the nervous system

The evolution of nervous systems in animals has always fascinated biologists, and thus multiple evolutionary scenarios have been proposed to explain the appearance of neurons and complex neuronal centers. However, the absence of a robust phylogenetic framework for animal interrelationships, the lack of a mechanistic understanding of development, and a recapitulative view of animal ontogeny have traditionally limited these scenarios. Only recently, the integration of advanced molecular and morphological studies in a broad range of animals has allowed to trace the evolution of developmental and neuronal characters on a better-resolved animal phylogeny. This has falsified most traditional scenarios for nervous system evolution, paving the way for the emergence of new testable hypotheses. Here we summarize recent progress in studies of nervous system development in major animal lineages and formulate some of the arising questions. In particular, we focus on how lineage analyses of nervous system development and a comparative study of the expression of neural-related genes has influenced our understanding of the evolution of an elaborated central nervous system in Bilateria. We argue that a phylogeny-guided study of neural development combining thorough descriptive and functional analyses is key to establish more robust scenarios for the origin and evolution of animal nervous systems.

An Introduction to the Study of Gastrotricha, with a Taxonomic Key to Families and Genera of the Group

Gastrotricha is a group of meiofaunal-sized, free-living invertebrates present in all aquatic ecosystems. The phylum includes over 860 species globally, of which 505 nominal species have been recorded in marine sandy sediments; another 355 taxa inhabit the freshwater environments, where they are recurrent members of the periphyton and epibenthos, and, to a lesser degree, of the plankton and interstitial fauna. Gastrotrichs are part of the permanent meiofauna and, in general, they rank among the top five groups for abundance within meiobenthic assemblages. The diversity, abundance, and ubiquity of Gastrotricha allow us to suppose an important role for these animals in aquatic ecosystems; however, ecological studies to prove this idea have been comparatively very few. This is mainly because the small size and transparency of their bodies make gastrotrichs difficult to discover in benthic samples; moreover, their contractility and fragility make their handling and morphological survey of the specimens rather difficult. Here we offer an overview, describe the basic techniques used to study these animals, and provide a key to known genera in an attempt to promote easy identification and to increase the number of researchers who may be interested in conducting studies on this understudied ecological group of microscopic organisms.

Revisiting metazoan phylogeny with genomic sampling of all phyla

Proper biological interpretation of a phylogeny can sometimes hinge on the placement of key taxa-or fail when such key taxa are not sampled. In this light, we here present the first attempt to investigate (though not conclusively resolve) animal relationships using genome-scale data from all phyla. Results from the site-heterogeneous CAT + GTR model recapitulate many established major clades, and strongly confirm some recent discoveries, such as a monophyletic Lophophorata, and a sister group relationship between Gnathifera and Chaetognatha, raising continued questions on the nature of the spiralian ancestor. We also explore matrix construction with an eye towards testing specific relationships; this approach uniquely recovers support for Panarthropoda, and shows that Lophotrochozoa (a subclade of Spiralia) can be constructed in strongly conflicting ways using different taxon- and/or orthologue sets. Dayhoff-6 recoding sacrifices information, but can also reveal surprising outcomes, e.g. full support for a clade of Lophophorata and Entoprocta + Cycliophora, a clade of Placozoa + Cnidaria, and raising support for Ctenophora as sister group to the remaining Metazoa, in a manner dependent on the gene and/or taxon sampling of the matrix in question. Future work should test the hypothesis that the few remaining uncertainties in animal phylogeny might reflect violations of the various stationarity assumptions used in contemporary inference methods.

OMA standalone: orthology inference among public and custom genomes and transcriptomes

Genomes and transcriptomes are now typically sequenced by individual laboratories but analyzing them often remains challenging. One essential step in many analyses lies in identifying orthologs—corresponding genes across multiple species—but this is far from trivial. The Orthologous MAtrix (OMA) database is a leading resource for identifying orthologs among publicly available, complete genomes. Here, we describe the OMA pipeline available as a standalone program for Linux and Mac. When run on a cluster, it has native support for the LSF, SGE, PBS Pro, and Slurm job schedulers and can scale up to thousands of parallel processes. Another key feature of OMA standalone is that users can combine their own data with existing public data by exporting genomes and precomputed alignments from the OMA database, which currently contains over 2100 complete genomes. We compare OMA standalone to other methods in the context of phylogenetic tree inference, by inferring a phylogeny of Lophotrochozoa, a challenging clade within the protostomes. We also discuss other potential applications of OMA standalone, including identifying gene families having undergone duplications/losses in specific clades, and identifying potential drug targets in nonmodel organisms. OMA standalone is available under the permissive open source Mozilla Public License Version 2.0.

Integrated data analysis allows the establishment of a new, cosmopolitan genus of marine Macrodasyida (Gastrotricha)

Macrodasyida (phylum Gastrotricha) comprises 365 species distributed across 34 genera and 10 families. However, current classification is under revision due to the contradictory results of molecular and morphological cladistic analyses. Studies aimed at bridging the gaps took advantage of supplementary assessments of poorly known species and particularly from observations of new taxa showing original traits that could help to identify plesiomorphic character states in these anatomically diverse micrometazoa. We follow this path by describing three new interesting macrodasyidan species respectively from Italy, Brazil and Sweden. In many respects, the new species resemble most closely species of the genus Macrodasys; however, details of the external morphology, in combination with the different lay-out of the reproductive system and the tiny spermatozoa lacking a visible flagellum, suggest they belong to a new genus, possibly in the family Macrodasyidae. These hypotheses are supported by the phylogenetic relationships of 47 taxa inferred from analyses of the 18S rRNA gene, which found the new species clustering with Thaidasys tongiorgii in a subset of a larger clade containing Macrodasys. Accordingly, the establishment of the following taxa is proposed: Kryptodasys gen. nov., K. marcocurinii sp. nov., K. carlosrochai sp. nov. and K. ulfjondeliusi sp. nov.

Dicyemida and Orthonectida: Two Stories of Body Plan Simplification

Two enigmatic groups of morphologically simple parasites of invertebrates, the Dicyemida (syn. Rhombozoa) and the Orthonectida, since the 19th century have been usually considered as two classes of the phylum Mesozoa. Early molecular evidence suggested their relationship within the Spiralia (=Lophotrochozoa), however, high rates of dicyemid and orthonectid sequence evolution led to contradicting phylogeny reconstructions. Genomic data for orthonectids revealed that they are highly simplified spiralians and possess a reduced set of genes involved in metazoan development and body patterning. Acquiring genomic data for dicyemids, however, remains a challenge due to complex genome rearrangements including chromatin diminution and generation of extrachromosomal circular DNAs, which are reported to occur during the development of somatic cells. We performed genomic sequencing of one species of Dicyema, and obtained transcriptomic data for two Dicyema spp. Homeodomain (homeobox) transcription factors, G-protein-coupled receptors, and many other protein families have undergone a massive reduction in dicyemids compared to other animals. There is also apparent reduction of the bilaterian gene complements encoding components of the neuromuscular systems. We constructed and analyzed a large dataset of predicted orthologous proteins from three species of Dicyema and a set of spiralian animals including the newly sequenced genome of the orthonectid Intoshia linei. Bayesian analyses recovered the orthonectid lineage within the Annelida. In contrast, dicyemids form a separate clade with weak affinity to the Rouphozoa (Platyhelminthes plus Gastrotricha) or (Entoprocta plus Cycliophora) suggesting that the historically proposed Mesozoa is a polyphyletic taxon. Thus, dramatic simplification of body plans in dicyemids and orthonectids, as well as their intricate life cycles that combine metagenesis and heterogony, evolved independently in these two lineages.

Amiskwia is a large Cambrian gnathiferan with complex gnathostomulid-like jaws

Phylogenomic studies have greatly improved our understanding of the animal tree of life but the relationships of many clades remain ambiguous. Here we show that the rare soft-bodied animal Amiskwia from the Cambrian of Canada and China, which has variously been considered a chaetognath, a nemertine, allied to molluscs, or a problematica, is related to gnathiferans. New specimens from the Burgess Shale (British Columbia, Canada) preserve a complex pharyngeal jaw apparatus composed of a pair of elements with teeth most similar to gnathostomulids. Amiskwia demonstrates that primitive spiralians were large and unsegmented, had a coelom, and were probably active nekto-benthic scavengers or predators. Secondary simplification and miniaturisation events likely occurred in response to shifting ecologies and adaptations to specialised planktonic habitats.

Inferring Ancient Relationships with Genomic Data: A Commentary on Current Practices

Contemporary phylogeneticists enjoy an embarrassment of riches, not only in the volumes of data now available, but also in the diversity of bioinformatic tools for handling these data. Here, I discuss a subset of these tools I consider well-suited to the task of inferring ancient relationships with coding sequence data in particular, encompassing data generation, orthology assignment, alignment and gene tree inference, supermatrix construction, and analysis under the best-fitting models applicable to large-scale datasets. Throughout, I compare and critique methods, considering both their theoretical principles and the details of their implementation, and offering practical tips on usage where appropriate. I also entertain different motivations for analyzing what are almost always originally DNA sequence data as codons, amino acids, and higher-order recodings. Although presented in a linear order, I see value in using the diversity of tools available to us to assess the sensitivity of clades of biological interest to different gene and taxon sets and analytical modes, which can be an indication of the presence of systematic error, of which a few forms remain poorly controlled by even the best available inference methods.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.

Sclerite-bearing annelids from the lower Cambrian of South China

Cambrian annelids are strikingly diverse and reveal important details of annelid character acquisition. Their contribution, however, to a wider understanding of the evolution of the trochozoans (encompassing the annelids as well as such groups as the brachiopods and molluscs) remains limited. Thus the early annelids had been linked to a variety of cataphract Cambrian metazoans, notably Wiwaxia and the halkieriids, but recent work assigns such fossils to stem-group molluscs. Here we report two new annelids from the Lower Cambrian Chengjiang Lagerstätte, South China. Ipoliknus avitus n. gen., n. sp. is biramous with neurochaetae and notochaetae, but significantly also bears dorsal spinose sclerites and dorso-lateral dentate sclerites. Adelochaeta sinensis n. gen., n. sp. is unique amongst Cambrian polychaetes in possessing the rod-like supports of the parapodia known as aciculae. This supports phylogenetic placement of Adelochaeta as sister to some more derived aciculate Palaeozoic taxa, but in contrast Ipoliknus is recovered as the most basal of the stem-group annelids. Sclerites and chaetae of I. avitus are interpreted respectively as the remnants and derivatives of a once more extensive cataphract covering that was a characteristic of more primitive trochozoans. The two sets of chaetae (noto- and neurochaetae) and two sets of sclerites (spinose and dentate) suggest that in a pre-annelid an earlier and more complete scleritome may have consisted of four zones of sclerites. Other cataphract taxa from the Lower Palaeozoic show a variety of scleritome configurations but establishing direct links with such basal annelids as Ipoliknus at present must remain conjectural.

The central nervous system of Oweniidae (Annelida) and its implications for the structure of the ancestral annelid brain

Background

Recent phylogenomic analyses congruently reveal a basal clade which consists of Oweniidae and Mageloniidae as sister group to the remaining Annelida. These results indicate that the last common ancestor of Annelida was a tube-dwelling organism. They also challenge traditional evolutionary hypotheses of different organ systems, among them the nervous system. In textbooks the central nervous system is described as consisting of a ganglionic ventral nervous system and a dorsally located brain with different tracts that connect certain parts of the brain to each other. Only limited information on the fine structure, however, is available for Oweniidae, which constitute the sister group (possibly together with Magelonidae) to all remaining annelids.

Results

The brain of Oweniidae is ring- shaped and basiepidermal. Ganglia, higher brain centers or complex sensory organs do not exist; instead the central nervous system is medullary. Posterior to the brain the ventral medullary cord arises directly from the ventral region of the brain in Myriowenia sp. while in Owenia fusiformis two medullary cords arise perpendicular to the brain ring, extend caudally and fuse posterior. The central nervous system is composed of a central neuropil and surrounding somata of the neurons. According to ultrastructural and histological data only one type of neuron is present in the central nervous system.

Conclusion

The central nervous system of Oweniidae is the simplest in terms of enlargement of the dorsal part of the brain and neuron distribution found among Annelida. Our investigation suggests that neither ganglia nor commissures inside the brain neuropil or clusters of polymorphic neurons were present in the annelid stem species. These structures evolved later within Annelida, most likely in the stem lineage of Amphinomidae, Sipuncula and Pleistoannelida. Palps were supposedly present in the last common ancestor of annelids and innervated by two nerves originating in the dorsal part of the brain. A broader comparison with species of each major spiralian clade shows the medullary nervous system to be a common feature and thus possibly representing the ancestral state of the spiralian nervous system. Moreover, ganglia and clusters of polymorphic neurons seemingly evolved independently in the compared taxa of Spiralia and Annelida.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0305-1) contains supplementary material, which is available to authorized users.

Evolutionary Implications of the microRNA- and piRNA Complement of Lepidodermella squamata (Gastrotricha)

Gastrotrichs-'hairy bellies'-are microscopic free-living animals inhabiting marine and freshwater habitats. Based on morphological and early molecular analyses, gastrotrichs were placed close to nematodes, but recent phylogenomic analyses have suggested their close relationship to flatworms (Platyhelminthes) within Spiralia. Small non-coding RNA data on e.g., microRNAs (miRNAs) and PIWI-interacting RNAs (piRNA) may help to resolve this long-standing question. MiRNAs are short post-transcriptional gene regulators that together with piRNAs play key roles in development. In a 'multi-omics' approach we here used small-RNA sequencing, available transcriptome and genomic data to unravel the miRNA- and piRNA complements along with the RNAi (RNA interference) protein machinery of Lepidodermella squamata (Gastrotricha, Chaetonotida). We identified 52 miRNA genes representing 35 highly conserved miRNA families specific to Eumetazoa, Bilateria, Protostomia, and Spiralia, respectively, with overall high similarities to platyhelminth miRNA complements. In addition, we found four large piRNA clusters that also resemble flatworm piRNAs but not those earlier described for nematodes. Congruently, transcriptomic annotation revealed that the Lepidodermella protein machinery is highly similar to flatworms, too. Taken together, miRNA, piRNA, and protein data support a close relationship of gastrotrichs and flatworms.

A Recurrent Motif: Diversity and Evolution of ShKT Domain Containing Proteins in the Vampire Snail Cumia reticulata

Proteins of the ShK superfamily are characterized by a small conserved domain (ShKT), first discovered in small venom peptides produced by sea anemones, and acting as specific inhibitors of voltage-dependent and calcium-activated K⁺ channels. The ShK superfamily includes both small toxic peptides and larger multifunctional proteins with various functions. ShK toxins are often important components of animal venoms, where they perform different biological functions including neurotoxic and immunosuppressive effects. Given their high specificity and effectiveness, they are currently regarded as promising pharmacological lead compounds for the treatment of autoimmune diseases. Here, we report on the molecular analysis of ShKT domain containing proteins produced by the Mediterranean vampire snail Cumia reticulata, an ectoparasitic gastropod that feeds on benthic fishes. The high specificity of expression of most ShK transcripts in salivary glands identifies them as relevant components of C. reticulata venom. These ShK proteins display various structural architectures, being produced either as single-domain secretory peptides, or as larger proteins combining the ShKT with M12 or CAP domains. Both ShKT-containing genes and their internal ShKT domains undergo frequent duplication events in C. reticulata, ensuring a high level of variability that is likely to play a role in increasing the range of their potential molecular targets.

A New Spiralian Phylogeny Places the Enigmatic Arrow Worms among Gnathiferans

Chaetognaths (arrow worms) are an enigmatic group of marine animals whose phylogenetic position remains elusive, in part because they display a mix of developmental and morphological characters associated with other groups [1, 2]. In particular, it remains unclear whether they are a sister group to protostomes [1, 2], one of the principal animal superclades, or whether they bear a closer relationship with some spiralian phyla [3, 4]. Addressing the phylogenetic position of chaetognaths and refining our understanding of relationships among spiralians are essential to fully comprehend character changes during bilaterian evolution [5]. To tackle these questions, we generated new transcriptomes for ten chaetognath species, compiling an extensive phylogenomic dataset that maximizes data occupancy and taxonomic representation. We employed inference methods that consider rate and compositional heterogeneity across taxa to avoid limitations of earlier analyses [6]. In this way, we greatly improved the resolution of the protostome tree of life. We find that chaetognaths cluster together with rotifers, gnathostomulids, and micrognathozoans within an expanded Gnathifera clade and that this clade is the sister group to other spiralians [7, 8]. Our analysis shows that several previously proposed groupings are likely due to systematic error, and we propose a revised organization of Lophotrochozoa with three main clades: Tetraneuralia (mollusks and entoprocts), Lophophorata (brachiopods, phoronids, and ectoprocts), and a third unnamed clade gathering annelids, nemerteans, and platyhelminthes. Consideration of classical morphological, developmental, and genomic characters in light of this topology indicates secondary loss as a fundamental trend in spiralian evolution.

Hox gene expression in postmetamorphic juveniles of the brachiopod Terebratalia transversa

BACKGROUND

Hox genes encode a family of homeodomain containing transcription factors that are clustered together on chromosomes of many Bilateria. Some bilaterian lineages express these genes during embryogenesis in spatial and/or temporal order according to their arrangement in the cluster, a phenomenon referred to as collinearity. Expression of Hox genes is well studied during embryonic and larval development of numerous species; however, relatively few studies focus on the comparison of pre- and postmetamorphic expression of Hox genes in animals with biphasic life cycle. Recently, the expression of Hox genes was described for embryos and larvae of Terebratalia transversa, a rhynchonelliformean brachiopod, which possesses distinct metamorphosis from planktonic larvae to sessile juveniles. During premetamorphic development, T. transversa does not exhibit spatial collinearity and several of its Hox genes are recruited for the morphogenesis of novel structures. In our study, we determined the expression of Hox genes in postmetamorphic juveniles of T. transversa in order to examine metamorphosis-related changes of expression patterns and to test whether Hox genes are expressed in the spatially collinear way in the postmetamorphic juveniles.

RESULTS

Hox genes are expressed in a spatially non-collinear manner in juveniles, generally showing similar patterns as ones observed in competent larvae: genes labial and post1 are expressed in chaetae-related structures, sex combs reduced in the shell-forming epithelium, whereas lox5 and lox4 in dorso-posterior epidermis. After metamorphosis, expression of genes proboscipedia, hox3, deformed and antennapedia becomes restricted to, respectively, shell musculature, prospective hinge rudiments and pedicle musculature and epidermis.

CONCLUSIONS

All developmental stages of T. transversa, including postmetamorphic juveniles, exhibit a spatial non-collinear Hox genes expression with only minor changes observed between pre- and postmetamorphic stages. Our results are concordant with morphological observation that metamorphosis in rhynchonelliformean brachiopods, despite being rapid, is rather gradual. The most drastic changes in Hox gene expression patterns observed during metamorphosis could be explained by the inversion of the mantle lobe, which relocates some of the more posterior larval structures into the anterior edge of the juveniles. Co-option of Hox genes for the morphogenesis of novel structures is even more pronounced in postmetamorphic brachiopods when compared to larvae.