Nemertean and phoronid genomes reveal lophotrochozoan evolution and the origin of bilaterian heads

Exceptionally Preserved Setae: A Possible Morphological Synapomorphy of Cambrian Lophotrochozoans

Cambrian Lagerstätten yield exceptionally preserved fossils that have greatly improved our understanding of the origin and evolution of animal groups. Brachiopoda, a phylum of bivalved marine invertebrates nested firmly within the lophotrochozoan protostomes, are widely recovered in such Lagerstätten. The marginal chitinous setae (or chaetae) of brachiopods are the most commonly described soft tissue and have been interpreted as performing a defensive and/or sensory role. Despite their relatively common appearance in Cambrian Lagerstätten, the origin, function, and evolution of setae in the Brachiopoda is poorly known. Here, we document exquisitely preserved setal structures from South China and Laurentia paleocontinents giving new insights into their formation, microstructure and preservation mode. New setae typically make their appearance within the follicle of a neighbouring older seta and then branches off laterally forming its own follicle. Setal microstructure is likely to be composed of many canals, highly comparable to setae of their recent counterparts. Moreover, setae recovered from these palaeo-continents present different preservation: aside from the normal preservation of iron oxides and carbonaceous ingredients, some compositions of calcium are also detected in this originally chitinous organization. Investigating the evolutionary origins of chitinous setae, a specialized type found notably in lophotrochozoans such as brachiopods and annelids, reveals its presence in early Cambrian stem groups. This character likely serves as a morphological synapomorphy in lophotrochozoan evolution. However, the dearth of morpho-ultrastructure and comparative studies in Cambrian fossils presents a challenge in fully understanding this evolutionary development.

Evidence for Multiple Independent Expansions of Fox Gene Families Within Flatworms

Expansion and losses of gene families are important drivers of molecular evolution. A recent survey of Fox genes in flatworms revealed that this superfamily of multifunctional transcription factors, present in all animals, underwent extensive losses and expansions during platyhelminth evolution. In this paper, I analyzed Fox gene complement in four additional species of platyhelminths, that represent early-branching lineages in the flatworm phylogeny: catenulids (Stenostomum brevipharyngium and Stenostomum leucops) and macrostomorphs (Macrostomum hystrix and Macrostomum cliftonense). Phylogenetic analysis of Fox genes from this expanded set of species provided evidence for multiple independent expansions of Fox gene families within flatworms. Notably, FoxG, a panbilaterian brain-patterning gene, appears to be the least susceptible to duplication, while FoxJ1, a conserved ciliogenesis factor, has undergone extensive expansion in various flatworm lineages. Analysis of the single-cell atlas of S. brevipharyngium, combined with RNA in situ hybridization, elucidated the tissue-specific expression of the selected Fox genes: FoxG is expressed in the brain, three of the Fox genes (FoxN2/3-2, FoxO4 and FoxP1) are expressed in the pharyngeal cells of likely glandular function, while one of the FoxQD paralogs is specifically expressed in the protonephridium. Overall, the evolution of Fox genes in flatworms appears to be characterized by an early contraction of the gene complement, followed by lineage-specific expansions that have enabled the co-option of newly evolved paralogs into novel physiological and developmental functions.

Digestive System Development and Posterior Hox/Parahox Gene Expression During Larval Life and Metamorphosis of the Phoronid Phoronopsis harmeri

Phoronida is a small group of marine animals, most of which are characterized by a long larval period and complex metamorphosis. As a result of metamorphosis, their body changes so much that their true anterior and posterior ends are very close to each other, and the intestine becomes long and U-shaped. Using histology and electron microscopy, we have shown that the elongation and change in shape of the digestive tract that occurs during metamorphosis in Phoronopsis harmeri larvae is accompanied by the formation of new parts and changes in ultrastructure. At the same time, our in situ hybridization data suggest that the posterior markers Cdx and Post2 are expressed in posterior tissues at larval stages, during metamorphosis, and in juveniles, and that changes in their expression correlate with remodeling of the posterior parts of the digestive tract. Our data may shed light on the evolution of body patterning in animals undergoing complex metamorphosis.

Fusion, fission, and scrambling of the bilaterian genome in Bryozoa

Groups of orthologous genes are commonly found together on the same chromosome over vast evolutionary distances. This extensive physical gene linkage, known as macrosynteny, is seen between bilaterian phyla as divergent as Chordata, Echinodermata, Mollusca, and Nemertea. Here, we report a unique pattern of genome evolution in Bryozoa, an understudied phylum of colonial invertebrates. Using comparative genomics, we reconstruct the chromosomal evolutionary history of five bryozoans. Multiple ancient chromosome fusions followed by gene mixing led to the near-complete loss of bilaterian linkage groups in the ancestor of extant bryozoans. A second wave of rearrangements, including chromosome fission, then occurred independently in two bryozoan classes, further scrambling bryozoan genomes. We also discover at least five derived chromosomal fusion events shared between bryozoans and brachiopods, supporting the traditional but highly debated Lophophorata hypothesis and suggesting macrosynteny to be a potentially powerful source of phylogenetic information. Finally, we show that genome rearrangements led to the dispersion of genes from bryozoan Hox clusters onto multiple chromosomes. Our findings demonstrate that the canonical bilaterian genome structure has been lost across all studied representatives of an entire phylum, and reveal that linkage group fission can occur very frequently in specific lineages.

Identification of Six Novel Proteins Containing a ZP Module from Nemertean Species

During fertilization, a series of reactions between the eggs and spermatozoa proceed predominantly in a species-specific manner. The molecules mediating these species-specific reactions remain unknown except in a few organisms. In this study, we focused on two species belonging to the phylum Nemertea, Kulikovia alborostrata and K. fulva, and explored molecules involved in species-specific interactions between gametes. Orthologs of molecules known to be involved in species-specific reactions were not expressed in the ovaries of these two species. In contrast, we identified six novel proteins, named NeZPL1-NeZPL6, containing a ZP module. Among these, we found that NeZPL6 is located on the surface of an unfertilized egg and is suggested to be involved in its interaction with spermatozoa. Furthermore, we found an indel of three amino acids in the EGF-like domain of NeZPL6, which possibly confers species specificity to this interaction. Our results suggested the existence of a novel system for species recognition in animal gametes.

Evolution of the Cdk4/6-Cdkn2 system in invertebrates

The cell cycle is driven by cyclin-dependent kinases (Cdks). The decision whether the cell cycle proceeds is made during G1 phase, when Cdk4/6 functions. Cyclin-dependent kinase inhibitor 2 (Cdkn2) is a specific inhibitor of Cdk4/6, and their interaction depends on D84 in Cdkn2 and R24/31 in Cdk4/6. This knowledge is based mainly on studies in mammalian cells. Here, we comprehensively analyzed Cdk4/6 and Cdkn2 in invertebrates and found that Cdk4/6 was present in most of the investigated phyla, but the distribution of Cdkn2 was rather uneven among and within the phyla. The positive charge of R24/R31 in Cdk4/6 was conserved in all analyzed species in phyla with Cdkn2. The presence of Cdkn2 and the conservation of the positive charge were statistically correlated. We also found that Cdkn2 has been tightly linked to Fas associated factor 1 (Faf1) during evolution. We discuss potential interactions between Cdkn2 and Cdk4/6 in evolution and the possible cause of the strong conservation of the microsynteny.

The Evolution and Biological Activity of Metazoan Mixed Lineage Kinase Domain-Like Protein (MLKL)

In mammals, mixed lineage kinase domain-like protein (MLKL) is the executor of necroptosis. MLKL comprises an N-terminal domain (NTD), which alone suffices to trigger necroptosis by forming pores in the plasma membrane, and a C-terminal domain that inhibits the NTD activity. Evolutionarily, MLKL is poorly conserved in animals and not found in Protostomia. Although MLKL orthologs exist in invertebrate Deuterostomia, the biological activity of invertebrate MLKL is unknown. Herein, we examined 34 metazoan phyla and detected MLKL not only in Deuterostomia but also in Protostomia (Rotifera). The Rotifera MLKL exhibited low identities with non-Rotifera MLKL but shared relatively high identities with non-metazoan MLKL. In invertebrates, MLKL formed two phylogenetic clades, one of which was represented by Rotifera. In vertebrates, MLKL expression was tissue-specific and generally rich in immune organs. When expressed in human cells, the MLKL-NTD of Rotifera, Echinodermata, Urochordata, and Cephalochordata induced strong necroptosis. The necroptotic activity of Rotifera MLKL depended on a number of conserved residues. Together these findings provided new insights into the evolution of MLKL in Metazoa and revealed the biological activity of invertebrate MLKL.

Combinatorial Wnt signaling landscape during brachiopod anteroposterior patterning

BACKGROUND

Wnt signaling pathways play crucial roles in animal development. They establish embryonic axes, specify cell fates, and regulate tissue morphogenesis from the early embryo to organogenesis. It is becoming increasingly recognized that these distinct developmental outcomes depend upon dynamic interactions between multiple ligands, receptors, antagonists, and other pathway modulators, consolidating the view that a combinatorial "code" controls the output of Wnt signaling. However, due to the lack of comprehensive analyses of Wnt components in several animal groups, it remains unclear if specific combinations always give rise to specific outcomes, and if these combinatorial patterns are conserved throughout evolution.

RESULTS

In this work, we investigate the combinatorial expression of Wnt signaling components during the axial patterning of the brachiopod Terebratalia transversa. We find that T. transversa has a conserved repertoire of ligands, receptors, and antagonists. These genes are expressed throughout embryogenesis but undergo significant upregulation during axial elongation. At this stage, Frizzled domains occupy broad regions across the body while Wnt domains are narrower and distributed in partially overlapping patches; antagonists are mostly restricted to the anterior end. Based on their combinatorial expression, we identify a series of unique transcriptional subregions along the anteroposterior axis that coincide with the different morphological subdivisions of the brachiopod larval body. When comparing these data across the animal phylogeny, we find that the expression of Frizzled genes is relatively conserved, whereas the expression of Wnt genes is more variable.

CONCLUSIONS

Our results suggest that the differential activation of Wnt signaling pathways may play a role in regionalizing the anteroposterior axis of brachiopod larvae. More generally, our analyses suggest that changes in the receptor context of Wnt ligands may act as a mechanism for the evolution and diversification of the metazoan body axis.

Diversification of von Willebrand Factor A and Chitin-Binding Domains in Pif/BMSPs Among Mollusks

Pif is a shell matrix protein (SMP) identified in the nacreous layer of Pinctada fucata (Pfu) comprised two proteins, Pif97 and Pif 80. Pif97 contains a von Willebrand factor A (VWA) and chitin-binding domains, whereas Pif80 can bind calcium carbonate crystals. The VWA domain is conserved in the SMPs of various mollusk species; however, their phylogenetic relationship remains obscure. Furthermore, although the VWA domain participates in protein-protein interactions, its role in shell formation has not been established. Accordingly, in the current study, we investigate the phylogenetic relationship between PfuPif and other VWA domain-containing proteins in major mollusk species. The shell-related proteins containing VWA domains formed a large clade (the Pif/BMSP family) and were classified into eight subfamilies with unique sequential features, expression patterns, and taxa diversity. Furthermore, a pull-down assay using recombinant proteins containing the VWA domain of PfuPif 97 revealed that the VWA domain interacts with five nacreous layer-related SMPs of P. fucata, including Pif 80 and nacrein. Collectively, these results suggest that the VWA domain is important in the formation of organic complexes and participates in shell mineralisation.

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.

Emerging trends in the study of spiralian larvae

Many animals undergo indirect development, where their embryogenesis produces an intermediate life stage, or larva, that is often free-living and later metamorphoses into an adult. As their adult counterparts, larvae can have unique and diverse morphologies and occupy various ecological niches. Given their broad phylogenetic distribution, larvae have been central to hypotheses about animal evolution. However, the evolution of these intermediate forms and the developmental mechanisms diversifying animal life cycles are still debated. This review focuses on Spiralia, a large and diverse clade of bilaterally symmetrical animals with a fascinating array of larval forms, most notably the archetypical trochophore larva. We explore how classic research and modern advances have improved our understanding of spiralian larvae, their development, and evolution. Specifically, we examine three morphological features of spiralian larvae: the anterior neural system, the ciliary bands, and the posterior hyposphere. The combination of molecular and developmental evidence with modern high-throughput techniques, such as comparative genomics, single-cell transcriptomics, and epigenomics, is a promising strategy that will lead to new testable hypotheses about the mechanisms behind the evolution of larvae and life cycles in Spiralia and animals in general. We predict that the increasing number of available genomes for Spiralia and the optimization of genome-wide and single-cell approaches will unlock the study of many emerging spiralian taxa, transforming our views of the evolution of this animal group and their larvae.

Phoronida-A small clade with a big role in understanding the evolution of lophophorates

Phoronids, together with brachiopods and bryozoans, form the animal clade Lophophorata. Modern lophophorates are quite diverse-some can biomineralize while others are soft-bodied, they could be either solitary or colonial, and they develop through various eccentric larval stages that undergo different types of metamorphoses. The diversity of this clade is further enriched by numerous extinct fossil lineages with their own distinct body plans and life histories. In this review, I discuss how data on phoronid development, genetics, and morphology can inform our understanding of lophophorate evolution. The actinotrocha larvae of phoronids is a well documented example of intercalation of the new larval body plan, which can be used to study how new life stages emerge in animals with biphasic life cycle. The genomic and embryonic data from phoronids, in concert with studies of the fossil lophophorates, allow the more precise reconstruction of the evolution of lophophorate biomineralization. Finally, the regenerative and asexual abilities of phoronids can shed new light on the evolution of coloniality in lophophorates. As evident from those examples, Phoronida occupies a central role in the discussion of the evolution of lophophorate body plans and life histories.

Integrating Complex Life Cycles in Comparative Developmental Biology

The goal of comparative developmental biology is identifying mechanistic differences in embryonic development between different taxa and how these evolutionary changes have led to morphological and organizational differences in adult body plans. Much of this work has focused on direct-developing species in which the adult forms straight from the embryo and embryonic modifications have direct effects on the adult. However, most animal lineages are defined by indirect development, in which the embryo gives rise to a larval body plan and the adult forms by transformation of the larva. Historically, much of our understanding of complex life cycles is viewed through the lenses of ecology and zoology. In this review, we discuss the importance of establishing developmental rather than morphological or ecological criteria for defining developmental mode and explicitly considering the evolutionary implications of incorporating complex life cycles into broad developmental comparisons of embryos across metazoans.

Genomic and transcriptomic survey of bryozoan Hox and ParaHox genes with emphasis on phylactolaemate bryozoans

BACKGROUND

Bryozoans are mostly sessile aquatic colonial invertebrates belonging to the clade Lophotrochozoa, which unites many protostome bilaterian phyla such as molluscs, annelids and brachiopods. While Hox and ParaHox genes have been extensively studied in various lophotrochozoan lineages, investigations on Hox and ParaHox gene complements in bryozoans are scarce.

RESULTS

Herein, we present the most comprehensive survey of Hox and ParaHox gene complements in bryozoans using four genomes and 35 transcriptomes representing all bryozoan clades: Cheilostomata, Ctenostomata, Cyclostomata and Phylactolaemata. Using similarity searches, phylogenetic analyses and detailed manual curation, we have identified five Hox genes in bryozoans (pb, Dfd, Lox5, Lox4 and Post2) and one ParaHox gene (Cdx). Interestingly, we observed lineage-specific duplication of certain Hox and ParaHox genes (Dfd, Lox5 and Cdx) in some bryozoan lineages.

CONCLUSIONS

The bryozoan Hox cluster does not retain the ancestral lophotrochozoan condition but appears relatively simple (includes only five genes) and broken into two genomic regions, characterized by the loss and duplication of serval genes. Importantly, bryozoans share the lack of two Hox genes (Post1 and Scr) with their proposed sister-taxon, Phoronida, which suggests that those genes were missing in the most common ancestor of bryozoans and phoronids.

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.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

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.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Spatial compartmentalisation of bacteria in phoronid microbiomes

The phylum Phoronida comprises filter-feeding invertebrates that live in a protective tube sometimes reinforced with particulate material from the surrounding environments. Animals with these characteristics make promising candidate hosts for symbiotic bacteria, given the constant interactions with various bacterial colonizers, yet phoronids are one of the very few animal phyla with no available microbiome data whatsoever. Here, by sequencing the V4 region of the 16S rRNA gene, we compare bacterial microbiomes in whole phoronids, including both tube and living tissues, with those associated exclusively to the isolated tube and/or the naked animal inside. We also compare these communities with those from the surrounding water. Phoronid microbiomes from specimens belonging to the same colony but collected a month apart were significantly different, and bacterial taxa previously reported in association with invertebrates and sediment were found to drive this difference. The microbiomes associated with the tubes are very similar in composition to those isolated from whole animals. However, just over half of bacteria found in whole specimens are also found both in tubes and naked specimens. In conclusion, phoronids harbour bacterial microbiomes that differ from those in the surrounding water, but the composition of those microbiomes is not stable and appears to change in the same colony over a relatively short time frame. Considering individual spatial/anatomical compartments, the phoronid tube contributes most to the whole-animal microbiome.

Updated knowledge and a proposed nomenclature for nuclear receptors with two DNA binding domains (2DBD-NRs)

Nuclear receptors (NRs) are important transcriptional modulators in metazoans. Typical NRs possess a conserved DNA binding domain (DBD) and a ligand binding domain (LBD). Since we discovered a type of novel NRs each of them has two DBDs and single LBD (2DBD-NRs) more than decade ago, there has been very few studies about 2DBD-NRs. Recently, 2DBD-NRs have been only reported in Platyhelminths and Mollusca and are thought to be specific NRs to lophotrochozoan. In this study, we searched different databases and identified 2DBD-NRs in different animals from both protostomes and deuterostomes. Phylogenetic analysis shows that at least two ancient 2DBD-NR genes were present in the urbilaterian, a common ancestor of protostomes and deuterostomes. 2DBD-NRs underwent gene duplication and loss after the split of different animal phyla, most of them in a certain animal phylum are paralogues, rather than orthologues, like in other animal phyla. Amino acid sequence analysis shows that the conserved motifs in typical NRs are also present in 2DBD-NRs and they are gene specific. From our phylogenetic analysis of 2DBD-NRs and following the rule of Nomenclature System for the Nuclear Receptors, a nomenclature for 2DBD-NRs is proposed.

Phylogenetic tracing of midbrain-specific regulatory sequences suggests single origin of eubilaterian brains

Conserved cis-regulatory elements (CREs) control Engrailed-, Pax2-, and dachshund-related gene expression networks directing the formation and function of corresponding midbrain circuits in arthropods and vertebrates. Polarized outgroup analyses of 31 sequenced metazoan genomes representing all animal clades reveal the emergence of Pax2- and dachshund-related CRE-like sequences in anthozoan Cnidaria. The full complement, including Engrailed-related CRE-like sequences, is only detectable in spiralians, ecdysozoans, and chordates that have a brain; they exhibit comparable genomic locations and extensive nucleotide identities that reveal the presence of a conserved core domain, all of which are absent in non-neural genes and, together, distinguish them from randomly assembled sequences. Their presence concurs with a genetic boundary separating the rostral from caudal nervous systems, demonstrated for the metameric brains of annelids, arthropods, and chordates and the asegmental cycloneuralian and urochordate brain. These findings suggest that gene regulatory networks for midbrain circuit formation evolved within the lineage that led to the common ancestor of protostomes and deuterostomes.

Duplication and Losses of Opsin Genes in Lophotrochozoan Evolution

Opsins are G-coupled receptors playing a key role in metazoan visual processes. While many studies enriched our understanding of opsin diversity in several animal clades, the opsin evolution in Lophotrochozoa, one of the major metazoan groups, remains poorly understood. Using recently developed phylogenetic approaches, we investigated the opsin evolution in 74 lophotrochozoan genomes. We found that the common ancestor of Lophotrochozoa possessed at least seven opsin paralog groups that underwent divergent evolutionary history in the different phyla. Furthermore, we showed for the first time opsin-related molecules in Bilateria that we named pseudopsins, which may prove critical in uncovering opsin evolution.

The localization of Toll and Imd pathway and complement system components and their response to Vibrio infection in the nemertean Lineus ruber

BACKGROUND

Innate immunity is the first line of defense against pathogens. In animals, the Toll pathway, the Imd pathway, the complement system, and lectins are well-known mechanisms involved in innate immunity. Although these pathways and systems are well understood in vertebrates and arthropods, they are understudied in other invertebrates.

RESULTS

To shed light on immunity in the nemertean Lineus ruber, we performed a transcriptomic survey and identified the main components of the Toll pathway (e.g., myD88, dorsal/dif/NFκB-p65), the Imd pathway (e.g., imd, relish/NFκB-p105/100), the complement system (e.g., C3, cfb), and some lectins (FreD-Cs and C-lectins). In situ hybridization showed that TLRβ1, TLRβ2, and imd are expressed in the nervous system; the complement gene C3-1 is expressed in the gut; and the lectins are expressed in the nervous system, the blood, and the gut. To reveal their potential role in defense mechanisms, we performed immune challenge experiments, in which Lineus ruber specimens were exposed to the gram-negative bacteria Vibrio diazotrophicus. Our results show the upregulation of specific components of the Toll pathway (TLRα3, TLRβ1, and TLRβ2), the complement system (C3-1), and lectins (c-lectin2 and fred-c5).

CONCLUSIONS

Therefore, similarly to what occurs in other invertebrates, our study shows that components of the Toll pathway, the complement system, and lectins are involved in the immune response in the nemertean Lineus ruber. The presence of these pathways and systems in Lineus ruber, but also in other spiralians; in ecdysozoans; and in deuterostomes suggests that these pathways and systems were involved in the immune response in the stem species of Bilateria.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

The gold-ringed octopus (Amphioctopus fangsiao) genome and cerebral single-nucleus transcriptomes provide insights into the evolution of karyotype and neural novelties

BACKGROUND

Coleoid cephalopods have distinctive neural and morphological characteristics compared to other invertebrates. Early studies reported massive genomic rearrangements occurred before the split of octopus and squid lineages (Proc Natl Acad Sci U S A 116:3030-5, 2019), which might be related to the neural innovations of their brain, yet the details remain elusive. Here we combine genomic and single-nucleus transcriptome analyses to investigate the octopod chromosome evolution and cerebral characteristics.

RESULTS

We present a chromosome-level genome assembly of a gold-ringed octopus, Amphioctopus fangsiao, and a single-nucleus transcriptome of its supra-esophageal brain. Chromosome-level synteny analyses estimate that the chromosomes of the ancestral octopods experienced multiple chromosome fission/fusion and loss/gain events by comparing with the nautilus genome as outgroup, and that a conserved genome organization was detected during the evolutionary process from the last common octopod ancestor to their descendants. Besides, protocadherin, GPCR, and C2H2 ZNF genes are thought to be highly related to the neural innovations in cephalopods (Nature 524:220-4, 2015), and the chromosome analyses pinpointed several collinear modes of these genes on the octopod chromosomes, such as the collinearity between PCDH and C2H2 ZNF, as well as between GPCR and C2H2 ZNF. Phylogenetic analyses show that the expansion of the octopod protocadherin genes is driven by a tandem-duplication mechanism on one single chromosome, including two separate expansions at 65 million years ago (Ma) and 8-14 Ma, respectively. Furthermore, we identify eight cell types (i.e., cholinergic and glutamatergic neurons) in the supra-esophageal brain of A. fangsiao, and the single-cell expression analyses reveal the co-expression of protocadherin and GPCR in specific neural cells, which may contribute to the neural development and signal transductions in the octopod brain.

CONCLUSIONS

The octopod genome analyses reveal the dynamic evolutionary history of octopod chromosomes and neural-related gene families. The single-nucleus transcriptomes of the supra-esophageal brain indicate their cellular heterogeneities and functional interactions with other tissues (i.e., gill), which provides a foundation for further octopod cerebral studies.

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.

Investigation of Peptide Toxin Diversity in Ribbon Worms (Nemertea) Using a Transcriptomic Approach

Nemertea is a phylum of nonsegmented worms (supraphylum: Spiralia), also known as ribbon worms. The members of this phylum contain various toxins, including peptide toxins. Here, we provide a transcriptomic analysis of peptide toxins in 14 nemertean species, including Cephalothrix cf. simula, which was sequenced in the current study. The summarized data show that the number of toxin transcripts in the studied nemerteans varied from 12 to 82. The most represented groups of toxins were enzymes and ion channel inhibitors, which, in total, reached a proportion of 72% in some species, and the least represented were pore-forming toxins and neurotoxins, the total proportion of which did not exceed 18%. The study revealed that nemerteans possess a much greater variety of toxins than previously thought and showed that these animals are a promising object for the investigation of venom diversity and evolution, and in the search for new peptide toxins.

New insights into the evolutionary dynamic and lineage divergence of gasdermin E in metazoa

Gasdermin (GSDM) is a family of pore-forming proteins that induce pyroptosis. To date, the origin and evolution of GSDM in Metazoa remain elusive. Here, we found that GSDM emerged early in Placozoa but is absent in a large number of invertebrates. In the lower vertebrate, fish, three types of GSDME, i.e., GSDMEa, GSDMEb, and a previously unreported type (designated GSDMEc), were idenitied. Evolutionarily, the three GSDMEs are distinctly separated: GSDMEa is closely related to tetrapod GSDME; GSDMEb exists exclusively in fish; GSDMEc forms the lineage root of tetrapod GSDMA/B/C/D. GSDMEc shares conserved genomic features with and is probably the prototype of GSDMA, which we found existing in all tetrapod classes. GSDMEc displays fast evolutionary dynamics, likely as a result of genomic transposition. A cross-metazoan analysis of GSDME revealed that GSDMEa shares a conserved caspase recognition motif with the GSDME of tetrapods and cnidarians, whereas GSDMEb has a unique caspase recognition motif similar to that of mammalian GSDMD, and GSDMEc exhibits no apparent caspase recognition motif. Through functional test, four highly conserved residues in vertebrate GSDME proved to be essential to auto-inhibition. Together our results provide new insights into the origin, evolution, and function of metazoan GSDMs.

Evolution of EGF-like and Zona pellucida domains containing shell matrix proteins in mollusks

Several types of shell matrix proteins (SMPs) have been identified in molluskan shells. Their diversity is the consequence of various molecular processes, including domain shuffling and gene duplication. However, the evolutionary origin of most SMPs remains unclear. In this study, we investigated the evolutionary process EGF-like and zona pellucida (ZP) domains containing SMPs. Two types of the proteins (EGF-like protein (EGFL) and EGF-like and ZP domains containing protein (EGFZP)) were found in the pearl oyster, Pinctada fucata. In contrast, only EGFZP was identified in the gastropods. Phylogenetic analysis and genomic arrangement studies showed that EGFL and EGFZP formed a clade in bivalves, and their encoding genes were localized in tandem repeats on the same scaffold. In P. fucata, EGFL genes were expressed in the outer part of mantle epithelial cells are related to the calcitic shell formation. However, in both P. fucata and the limpet Nipponacmea fuscoviridis, EGFZP genes were expressed in the inner part of the mantle epithelial cells are related to aragonitic shell formation. Furthermore, our analysis showed that in P. fucata, the ZP domain interacts with eight SMPs that have various functions in the nacreous shell mineralization. The data suggest that the ZP domain can interact with other SMPs, and EGFL evolution in pterimorph bivalves represents an example of neo-functionalization that involves the acquisition of a novel protein through gene duplication.

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.

Head transcriptome profiling of glossiphoniid leech (Helobdella austinensis) reveals clues about proboscis development

Cephalization refers to the evolutionary trend towards the concentration of neural tissues, sensory organs, mouth and associated structures at the front end of bilaterian animals. Comprehensive studies on gene expression related to the anterior formation in invertebrate models are currently lacking. In this study, we performed de novo transcriptional profiling on a proboscis-bearing leech (Helobdella austinensis) to identify differentially expressed genes (DEGs) in the anterior versus other parts of the body, in particular to find clues as to the development of the proboscis. Between the head and the body, 132 head-specific DEGs were identified, of which we chose 11 to investigate their developmental function during embryogenesis. Analysis of the spatial expression of these genes using in situ hybridization showed that they were characteristically expressed in the anterior region of the developing embryo, including the proboscis. Our results provide information on the genes related to head formation and insights into the function of proboscis-related genes during organogenesis with the potential roles of genes not yet characterized.

Lineid Heteronemerteans (Nemertea: Pilidiophora) from Sagami Bay, Japan, with Some Proposals for the Family-Level Classification System

We performed molecular phylogenetic analyses using nucleotide sequence data from five genes (16S rRNA, cytochrome c oxidase subunit I, 18S rRNA, 28S rRNA, histone H3) determined from 36 specimens representing 25 (10 named, 15 unnamed) species of lineid heteronemerteans collected in Sagami Bay, Japan, along with other sequences obtained from public databases. External features of the 25 species are briefly described and illustrated. Lineus fuscoviridis Takakura, 1898 is transferred to Notospermus Huschke, 1830. Our molecular analysis indicated that one of our Notospermus specimens from Sagami Bay and material previously collected from the Seto Inland Sea are conspecific with Notospermus geniculatus (Delle Chiaje, 1822) s. str. [the year of publication is not 1828 as previously regarded], originally described from Naples, Italy. The new species Siphonenteron nakanoi is established; our tree shows it as the sister taxon to Siphonenteron bilineatum Renier in Meneghini, 1847; it differs from the latter in having more diffuse, scattered yellow pigmentation at the anterolateral margin of the head on each side. Pros and cons of lumping and splitting Lineidae in future studies are discussed.

The genome of the poecilogonous annelid Streblospio benedicti

Streblospio benedicti is a common marine annelid that has become an important model for developmental evolution. It is the only known example of poecilogony (where two distinct developmental modes occur within a single species) that is due to a heritable difference in egg size. The dimorphic developmental programs and life-histories exhibited in this species depend on differences within the genome, making it an optimal model for understanding the genomic basis of developmental divergence. Studies using S. benedicti have begun to uncover the genetic and genomic principles that underlie developmental uncoupling, but until now they have been limited by the lack of availability of genomic tools. Here we present an annotated chromosomal-level genome assembly of S. benedicti generated from a combination of Illumina reads, Nanopore long reads, Chicago and Hi-C chromatin interaction sequencing, and a genetic map from experimental crosses. At 701.4 Mb, the S. benedicti genome is the largest annelid genome to date that has been assembled to chromosomal scaffolds. The complete genome of S. benedicti is valuable for functional genomic analyses of development and evolution, as well as phylogenetic comparison within the annelida and the Lophotrochozoa. Despite having two developmental modes, there is no evidence of genome duplication or substantial gene number expansions. Instead, lineage specific repeats account for much of the expansion of this genome compared to other annelids.

Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome

The mutualism between the giant tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has been extensively researched over the past 40 years. However, the lack of the host whole genome information has impeded the full comprehension of the genotype/phenotype interface in Riftia. Here we described the high-quality draft genome of Riftia, its complete mitogenome, and tissue-specific transcriptomic data. The Riftia genome presents signs of reductive evolution, with gene family contractions exceeding expansions. Expanded gene families are related to sulphur metabolism, detoxification, anti-oxidative stress, oxygen transport, immune system, and lysosomal digestion, reflecting evolutionary adaptations to the vent environment and endosymbiosis. Despite the derived body plan, the developmental gene repertoire in the gutless tubeworm is extremely conserved with the presence of a near intact and complete Hox cluster. Gene expression analyses establishes that the trophosome is a multi-functional organ marked by intracellular digestion of endosymbionts, storage of excretory products and haematopoietic functions. Overall, the plume and gonad tissues both in contact to the environment harbour highly expressed genes involved with cell cycle, programmed cell death, and immunity indicating a high cell turnover and defence mechanisms against pathogens. We posit that the innate immune system plays a more prominent role into the establishment of the symbiosis during the infection in the larval stage, rather than maintaining the symbiostasis in the trophosome. This genome bridges four decades of physiological research in Riftia, whilst simultaneously provides new insights into the development, whole organism functions and evolution in the giant tubeworm.

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.

The genome sequence of the bootlace worm, Lineus longissimus (Gunnerus, 1770)

We present a genome assembly from an individual Lineus longissimus (the bootlace worm; Nemertea; Pilidiophora; Heteronemertea; Lineidae). The genome sequence is 391 megabases in span. The majority of the assembly is scaffolded into 19 chromosomal pseudomolecules.

Tetrodotoxins Secretion and Voltage-Gated Sodium Channel Adaptation in the Ribbon Worm Kulikovia alborostrata (Takakura, 1898) (Nemertea)

Nemertea is a phylum of marine worms whose members bear various toxins, including tetrodotoxin (TTX) and its analogues. Despite the more than 30 years of studying TTXs in nemerteans, many questions regarding their functions and the mechanisms ensuring their accumulation and usage remain unclear. In the nemertean Kulikovia alborostrata, we studied TTX and 5,6,11-trideoxyTTX concentrations in body extracts and in released mucus, as well as various aspects of the TTX-positive-cell excretion system and voltage-gated sodium (Nav1) channel subtype 1 mutations contributing to the toxins' accumulation. For TTX detection, an immunohistological study with an anti-TTX antibody and HPLC-MS/MS were conducted. For Nav1 mutation searching, PCR amplification with specific primers, followed by Sanger sequencing, was used. The investigation revealed that, in response to an external stimulus, subepidermal TTX-positive cells released secretions actively to the body surface. The post-release toxin recovery in these cells was low for TTX and high for 5,6,11-trideoxyTTX in captivity. According to the data obtained, there is low probability of the targeted usage of TTX as a repellent, and targeted 5,6,11-trideoxyTTX secretion by TTX-bearing nemerteans was suggested as a possibility. The Sanger sequencing revealed identical sequences of the P-loop regions of Nav1 domains I-IV in all 17 studied individuals. Mutations comprising amino acid substitutions, probably contributing to nemertean channel resistance to TTX, were shown.

Genome Size Dynamics in Marine Ribbon Worms (Nemertea, Spiralia)

Nemertea is a phylum consisting of 1300 mostly marine species. Nemertea is distinguished by an eversible muscular proboscis, and most of the species are venomous. Genomic resources for this phylum are scarce despite their value in understanding biodiversity. Here, we present genome size estimates of Nemertea based on flow cytometry and their relationship to different morphological and developmental traits. Ancestral genome size estimations were done across the nemertean phylogeny. The results increase the available genome size estimates for Nemertea three-fold. Our analyses show that Nemertea has a narrow genome size range (0.43-3.89 pg) compared to other phyla in Lophotrochozoa. A relationship between genome size and evolutionary rate, developmental modes, and habitat was found. Trait analyses show that the highest evolutionary rate of genome size is found in upper intertidal, viviparous species with direct development. Despite previous findings, body size in nemerteans was not correlated with genome size. A relatively small genome (1.18 pg) is assumed for the most recent common ancestor of all extant nemerteans. The results provide an important basis for future studies in nemertean genomics, which will be instrumental to understanding the evolution of this enigmatic and often neglected phylum.

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.

Functional Characterization of the Nemertide α Family of Peptide Toxins

Peptide toxins find use in medicine, biotechnology, and agriculture. They are exploited as pharmaceutical tools, particularly for the investigation of ion channels. Here, we report the synthesis and activity of a novel family of peptide toxins: the cystine-knotted α nemertides. Following the prototypic α-1 and -2 (1 and 2), six more nemertides were discovered by mining of available nemertean transcriptomes. Here, we describe their synthesis using solid phase peptide chemistry and their oxidative folding by using an improved protocol. Nemertides α-2 to α-7 (2-7) were produced to characterize their effect on voltage-gated sodium channels (Blatella germanica BgNa_V1 and mammalian Na_Vs1.1-1.8). In addition, ion channel activities were matched to in vivo tests using an Artemia microwell assay. Although nemertides demonstrate high sequence similarity, they display variability in activity on the tested Na_Vs. The nemertides are all highly toxic to Artemia, with EC₅₀ values in the sub-low micromolar range, and all manifest preference for the insect BgNa_V1 channel. Structure-activity relationship analysis revealed key residues for Na_V-subtype selectivity. Combined with low EC₅₀ values (e.g., Na_V1.1: 7.9 nM (α-6); Na_V1.3: 9.4 nM (α-5); Na_V1.4: 14.6 nM (α-4)) this underscores the potential utility of α-nemertides for rational optimization to improve selectivity.

Nemertean, brachiopod and phoronid neuropeptidomics reveals ancestral spiralian signalling systems

Neuropeptides are diverse signaling molecules in animals commonly acting through G-protein coupled receptors (GPCRs). Neuropeptides and their receptors underwent extensive diversification in bilaterians and the relationships of many peptide–receptor systems have been clarified. However, we lack a detailed picture of neuropeptide evolution in lophotrochozoans as in-depth studies only exist for mollusks and annelids. Here, we analyze peptidergic systems in Nemertea, Brachiopoda, and Phoronida. We screened transcriptomes from 13 nemertean, 6 brachiopod, and 4 phoronid species for proneuropeptides and neuropeptide GPCRs. With mass spectrometry from the nemertean Lineus longissimus, we validated several predicted peptides and identified novel ones. Molecular phylogeny combined with peptide-sequence and gene-structure comparisons allowed us to comprehensively map spiralian neuropeptide evolution. We found most mollusk and annelid peptidergic systems also in nemerteans, brachiopods, and phoronids. We uncovered previously hidden relationships including the orthologies of spiralian CCWamides to arthropod agatoxin-like peptides and of mollusk APGWamides to RGWamides from annelids, with ortholog systems in nemerteans, brachiopods, and phoronids. We found that pleurin neuropeptides previously only found in mollusks are also present in nemerteans and brachiopods. We also identified cases of gene family duplications and losses. These include a protostome-specific expansion of RFamide/Wamide signaling, a spiralian expansion of GnRH-related peptides, and duplications of vasopressin/oxytocin before the divergence of brachiopods, phoronids, and nemerteans. This analysis expands our knowledge of peptidergic signaling in spiralians and other protostomes. Our annotated data set of nearly 1,300 proneuropeptide sequences and 600 GPCRs presents a useful resource for further studies of neuropeptide signaling.

Eighteen Coral Genomes Reveal the Evolutionary Origin of Acropora Strategies to Accommodate Environmental Changes

The genus Acropora comprises the most diverse and abundant scleractinian corals (Anthozoa, Cnidaria) in coral reefs, the most diverse marine ecosystems on Earth. However, the genetic basis for the success and wide distribution of Acropora are unknown. Here, we sequenced complete genomes of 15 Acropora species and 3 other acroporid taxa belonging to the genera Montipora and Astreopora to examine genomic novelties that explain their evolutionary success. We successfully obtained reasonable draft genomes of all 18 species. Molecular dating indicates that the Acropora ancestor survived warm periods without sea ice from the mid or late Cretaceous to the Early Eocene and that diversification of Acropora may have been enhanced by subsequent cooling periods. In general, the scleractinian gene repertoire is highly conserved; however, coral- or cnidarian-specific possible stress response genes are tandemly duplicated in Acropora. Enzymes that cleave dimethlysulfonioproprionate into dimethyl sulfide, which promotes cloud formation and combats greenhouse gasses, are the most duplicated genes in the Acropora ancestor. These may have been acquired by horizontal gene transfer from algal symbionts belonging to the family Symbiodiniaceae, or from coccolithophores, suggesting that although functions of this enzyme in Acropora are unclear, Acropora may have survived warmer marine environments in the past by enhancing cloud formation. In addition, possible antimicrobial peptides and symbiosis-related genes are under positive selection in Acropora, perhaps enabling adaptation to diverse environments. Our results suggest unique Acropora adaptations to ancient, warm marine environments and provide insights into its capacity to adjust to rising seawater temperatures.

Evidence from oyster suggests an ancient role for Pdx in regulating insulin gene expression in animals

Hox and ParaHox genes encode transcription factors with similar expression patterns in divergent animals. The Pdx (Xlox) homeobox gene, for example, is expressed in a sharp spatial domain in the endodermal cell layer of the gut in chordates, echinoderms, annelids and molluscs. The significance of comparable gene expression patterns is unclear because it is not known if downstream transcriptional targets are also conserved. Here, we report evidence indicating that a classic transcriptional target of Pdx1 in vertebrates, the insulin gene, is a likely direct target of Pdx in Pacific oyster adults. We show that one insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue. Transcriptomic comparison suggests that this tissue plays a similar role to the vertebrate pancreas. Using ATAC-seq and ChIP, we identify an upstream regulatory element of the cgILP gene which shows binding interaction with cgPdx protein in oyster hepatopancreas and demonstrate, using a cell culture assay, that the oyster Pdx can act as a transcriptional activator through this site, possibly in synergy with NeuroD. These data argue that a classic homeodomain-target gene interaction dates back to the origin of Bilateria.

In vertebrates insulin is a direct transcriptional target of Pdx: the same is true in Pacific oysters and the authors show insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue, showing this gene interaction dates back to the origin of Bilateria.

Mobilizing molluscan models and genomes in biology

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

The genome of the euryhaline rotifer Brachionus paranguensis: Potential use in molecular ecotoxicology

Brachionus spp. rotifers have been proposed as model organisms for ecotoxicological studies. We analyzed the whole-genome sequence of B. paranguensis through NextDenovo, resulting in a total length of 106.2 Mb and 71 contigs. The N50 and the GC content were 4.13 Mb and 28%, respectively. A total of 18,501 genes were predicted within the genome of B. paranguensis. Prominent detoxification-related gene families of phase I and II detoxifications have been investigated. In parallel with other Brachionus rotifers, high gene expansion was observed in CYP clan 3 and GST sigma class in B. paranguensis. Moreover, species-specific expansion of sulfotransferase (SULTs) and gain of UDP-glucuronosyltransferases (UGTs) through horizontal gene transfer has been specifically found within B. plicatilis complex. This whole-genome analysis of B. paranguensis provides a basis for molecular ecotoxicological studies and provides useful information for comparative studies of the evolution of detoxification mechanisms in Brachionus spp.

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.

Developmental processes in Ediacara macrofossils

The Ediacara Biota preserves the oldest fossil evidence of abundant, complex metazoans. Despite their significance, assigning individual taxa to specific phylogenetic groups has proved problematic. To better understand these forms, we identify developmentally controlled characters in representative taxa from the Ediacaran White Sea assemblage and compare them with the regulatory tools underlying similar traits in modern organisms. This analysis demonstrates that the genetic pathways for multicellularity, axial polarity, musculature, and a nervous system were likely present in some of these early animals. Equally meaningful is the absence of evidence for major differentiation of macroscopic body units, including distinct organs, localized sensory machinery or appendages. Together these traits help to better constrain the phylogenetic position of several key Ediacara taxa and inform our views of early metazoan evolution. An apparent lack of heads with concentrated sensory machinery or ventral nerve cords in such taxa supports the hypothesis that these evolved independently in disparate bilaterian clades.