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 towards 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 Mbp. Hi-C scaffolding yielded chromosome level scaffolds, with a scaffold N50 of 10.0 Mbp 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.

Acceleration of genome rearrangement in clitellate annelids.

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution.

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

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.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting.

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

The embryology, metamorphosis, and muscle development of Schizocardium karankawa sp. nov. (Enteropneusta) from the Gulf of Mexico.

Schizocardium karankawa sp. nov. has been collected from subtidal muds of the Laguna Madre, Texas, and the Mississippi coast, Gulf of Mexico. The Texas population is reproductive from early February to mid-April. Gametes are liberated by a small incision in a gonad. Oocyte germinal vesicle breakdown is increased in the presence of sperm, and the highest fertilization success was in the artificial seawater Jamarin U. Manually dechorionated embryos develop normally. Development was asynchronous via a tornaria larva, metamorphosis and maintained to the juvenile worm 6 gill-pore stage. Phalloidin-labeled late-stage tornaria revealed retractor muscles that connect the pericardial sac with the apical tuft anteriorly, the oesophagus ventrally, and muscle cells of the early mesocoels. The muscle development of early juvenile worms began with dorso-lateral trunk muscles, lateral trunk bands, and sphincters around the gill pores and anus. Adult worms are characterized by a stomochord that bifurcates anteriorly into paired vermiform processes, gill bars that extend almost the entire dorsal to ventral branchial region resulting in a narrow ventral hypobranchial ridge, and an elaborate epibranchial organ with six zones of discrete cell types. The trunk has up to three rows of liver sacs, and lateral gonads. The acorn worm evo-devo model species Saccoglossus kowalevskii, Ptychodera flava, and Schizocardium californicum are phylogenetically distant with disparate life histories. S. karnakawa from S. californicum are phylogenetically close, and differences between them that become apparent as adult worms include the number of gill pores and hepatic sacs, and elaborations of the heart-kidney-stomochord complex. An important challenge for evolutionary developmental biology is to form links from phylogenetically distant and large-scale differences to phylogenetically close and small-scale differences. This description of the embryology, development, and adult morphology of S. karankawa permits investigations into how acorn worm development evolves at fine scales.

Phylogenomics reveals deep relationships and diversification within phylactolaemate bryozoans.

Bryozoans are mostly sessile colonial invertebrates that inhabit all kinds of aquatic ecosystems. Extant bryozoan species fall into two clades with one of them, Phylactolaemata, being the only exclusively freshwater clade. Phylogenetic relationships within the class Phylactolaemata have long been controversial owing to their limited distinguishable characteristics that reflect evolutionary relationships. Here, we present the first phylogenomic analysis of Phylactolaemata using transcriptomic data combined with dense taxon sampling of six families to better resolve the interrelationships and to estimate divergence time. Using maximum-likelihood and Bayesian inference approaches, we recovered a robust phylogeny for Phylactolaemata in which the interfamilial relationships are fully resolved. We show Stephanellidae is the sister taxon of all other phylactolaemates and confirm that Lophopodidae represents the second offshoot within the phylactolaemate tree. Plumatella fruticosa clearly falls outside Plumatellidae as previous investigations have suggested, and instead clusters with Pectinatellidae and Cristatellidae as the sister taxon of Fredericellidae. Our results demonstrate that cryptic speciation is very likely in F. sultana and in two species of Plumatella (P. repens and P. casmiana). Divergence time estimates show that Phylactolaemata appeared at the end of the Ediacaran and started to diverge in the Silurian, although confidence intervals were large for most nodes. The radiation of most extant phylactolaemate families occurred mainly in the Palaeogene and Neogene highlighting post-extinction diversification.

Genome of the lepidopleurid chiton Hanleya hanleyi (Mollusca, Polyplacophora).

Mollusca is the second most species-rich phylum and includes animals as disparate as octopuses, clams, and chitons. Dozens of molluscan genomes are available, but only one representative of the subphylum Aculifera, the sister taxon to all other molluscs, has been sequenced to date, hindering comparative and evolutionary studies. To facilitate evolutionary studies across Mollusca, we sequenced the genome of a second aculiferan mollusc, the lepidopleurid chiton Hanleya hanleyi (Bean 1844), using a hybrid approach combining Oxford Nanopore and Illumina reads. After purging redundant haplotigs and removing contamination from this 1.3% heterozygous genome, we produced a 2.5 Gbp haploid assembly (>4X the size of the other chiton genome sequenced to date) with an N50 of 65.0 Kbp. Despite a fragmented assembly, the genome is rather complete (92.0% of BUSCOs detected; 79.4% complete plus 12.6% fragmented). Remarkably, the genome has the highest repeat content of any molluscan genome reported to date (>66%). Our gene annotation pipeline predicted 69,284 gene models (92.9% of BUSCOs detected; 81.8% complete plus 11.1% fragmented) of which 35,362 were supported by transcriptome and/or protein evidence. Phylogenomic analysis recovered Polyplacophora sister to all other sampled molluscs with maximal support. The Hanleya genome will be a valuable resource for studies of molluscan biology with diverse potential applications ranging from evolutionary and comparative genomics to molecular ecology.

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

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

Phylogenomic resolution of the root of Panpulmonata, a hyperdiverse radiation of gastropods: new insight into the evolution of air breathing.

Transitions to terrestriality have been associated with major animal radiations including land snails and slugs in Stylommatophora (>20 000 described species), the most successful lineage of 'pulmonates' (a non-monophyletic assemblage of air-breathing gastropods). However, phylogenomic studies have failed to robustly resolve relationships among traditional pulmonates and affiliated marine lineages that comprise clade Panpulmonata (Mollusca, Gastropoda), especially two key taxa: Sacoglossa, a group including photosynthetic sea slugs, and Siphonarioidea, intertidal limpet-like snails with a non-contractile pneumostome (narrow opening to a vascularized pallial cavity). To clarify the evolutionary history of the panpulmonate radiation, we performed phylogenomic analyses on datasets of up to 1160 nuclear protein-coding genes for 110 gastropods, including 40 new transcriptomes for Sacoglossa and Siphonarioidea. All 18 analyses recovered Sacoglossa as the sister group to a clade we named Pneumopulmonata, within which Siphonarioidea was sister to the remaining lineages in most analyses. Comparative modelling indicated shifts to marginal habitat (estuarine, mangrove and intertidal zones) preceded and accelerated the evolution of a pneumostome, present in the pneumopulmonate ancestor along with a one-sided plicate gill. These findings highlight key intermediate stages in the evolution of air-breathing snails, supporting the hypothesis that adaptation to marginal zones played an important role in major sea-to-land transitions.

Benchmarking Oxford Nanopore read assemblers for high-quality molluscan genomes.

Choosing the optimum assembly approach is essential to achieving a high-quality genome assembly suitable for comparative and evolutionary genomic investigations. Significant recent progress in long-read sequencing technologies such as PacBio and Oxford Nanopore Technologies (ONT) has also brought about a large variety of assemblers. Although these have been extensively tested on model species such as Homo sapiens and Drosophila melanogaster, such benchmarking has not been done in Mollusca, which lacks widely adopted model species. Molluscan genomes are notoriously rich in repeats and are often highly heterozygous, making their assembly challenging. Here, we benchmarked 10 assemblers based on ONT raw reads from two published molluscan genomes of differing properties, the gastropod Chrysomallon squamiferum (356.6 Mb, 1.59% heterozygosity) and the bivalve Mytilus coruscus (1593 Mb, 1.94% heterozygosity). By optimizing the assembly pipeline, we greatly improved both genomes from previously published versions. Our results suggested that 40-50X of ONT reads are sufficient for high-quality genomes, with Flye being the recommended assembler for compact and less heterozygous genomes exemplified by C. squamiferum, while NextDenovo excelled for more repetitive and heterozygous molluscan genomes exemplified by M. coruscus. A phylogenomic analysis using the two updated genomes with 32 other published high-quality lophotrochozoan genomes resulted in maximum support across all nodes, and we show that improved genome quality also leads to more complete matrices for phylogenomic inferences. Our benchmarking will ensure efficiency in future assemblies for molluscs and perhaps also for other marine phyla with few genomes available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

On the diversity of abyssal Dondersiidae (Mollusca: Aplacophora) with the description of a new genus, six new species, and a review of the family.

So far, of the 292 known species of solenogasters (Mollusca, Aplacophora), 62 belong to the clade Pholidoskepia Salvini-Plawen, 1978. Of these, only two have an abyssal distribution (3500-6000 m depth). Among Pholidoskepia, Dondersiidae Simroth, 1893 is the most diverse family. This study contributes to the knowledge of this family with the description of one new genus and six new species from the abyssal South Atlantic Ocean: Dondersia ? foraminosa sp. n., Nematomenia divae sp. n., Nematomenia brasiliensis sp. n., Nematomenia ? guineana sp. n., Helluoherpia vieiralaneroi sp. n. and Inopinatamenia (gen. n.) calamitosa sp. n. Specimens were collected during DIVA (Latitudinal Gradients of Deep-Sea BioDIVersity in the Atlantic Ocean) expeditions in the Guinea (DIVA 2 Me 63/2, 2005) and Brazil (DIVA 3 Me 79/1, 2008) Basins. Specimens were characterized based primarily on the sclerites and internal anatomy, which was studied using histology. The importance of the radula and mantle sclerites for taxonomy is emphasized. Amended diagnoses for the family and some genera within this family are provided. This contribution increases the described diversity of Dondersiidae to ten genera and 38 species and highlights the need for more study of solenogasters in the deep sea.

Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics.

BACKGROUND: Heterobranchia is a diverse clade of marine, freshwater, and terrestrial gastropod molluscs. It includes such disparate taxa as nudibranchs, sea hares, bubble snails, pulmonate land snails and slugs, and a number of (mostly small-bodied) poorly known snails and slugs collectively referred to as the "lower heterobranchs". Evolutionary relationships within Heterobranchia have been challenging to resolve and the group has been subject to frequent and significant taxonomic revision. Mitochondrial (mt) genomes can be a useful molecular marker for phylogenetics but, to date, sequences have been available for only a relatively small subset of Heterobranchia. RESULTS: To assess the utility of mitochondrial genomes for resolving evolutionary relationships within this clade, eleven new mt genomes were sequenced including representatives of several groups of "lower heterobranchs". Maximum likelihood analyses of concatenated matrices of the thirteen protein coding genes found weak support for most higher-level relationships even after several taxa with extremely high rates of evolution were excluded. Bayesian inference with the CAT + GTR model resulted in a reconstruction that is much more consistent with the current understanding of heterobranch phylogeny. Notably, this analysis recovered Valvatoidea and Orbitestelloidea in a polytomy with a clade including all other heterobranchs, highlighting these taxa as important to understanding early heterobranch evolution. Also, dramatic gene rearrangements were detected within and between multiple clades. However, a single gene order is conserved across the majority of heterobranch clades. CONCLUSIONS: Analysis of mitochondrial genomes in a Bayesian framework with the site heterogeneous CAT + GTR model resulted in a topology largely consistent with the current understanding of heterobranch phylogeny. However, mitochondrial genomes appear to be too variable to serve as good phylogenetic markers for robustly resolving a number of deeper splits within this clade.

The Iron-Responsive Genome of the Chiton Acanthopleura granulata.

Molluscs biomineralize structures that vary in composition, form, and function, prompting questions about the genetic mechanisms responsible for their production and the evolution of these mechanisms. Chitons (Mollusca, Polyplacophora) are a promising system for studies of biomineralization because they build a range of calcified structures including shell plates and spine- or scale-like sclerites. Chitons also harden the calcified teeth of their rasp-like radula with a coat of iron (as magnetite). Here we present the genome of the West Indian fuzzy chiton Acanthopleura granulata, the first from any aculiferan mollusc. The A. granulata genome contains homologs of many genes associated with biomineralization in conchiferan molluscs. We expected chitons to lack genes previously identified from pathways conchiferans use to make biominerals like calcite and nacre because chitons do not use these materials in their shells. Surprisingly, the A. granulata genome has homologs of many of these genes, suggesting that the ancestral mollusc may have had a more diverse biomineralization toolkit than expected. The A. granulata genome has features that may be specialized for iron biomineralization, including a higher proportion of genes regulated directly by iron than other molluscs. A. granulata also produces two isoforms of soma-like ferritin: one is regulated by iron and similar in sequence to the soma-like ferritins of other molluscs, and the other is constitutively translated and is not found in other molluscs. The A. granulata genome is a resource for future studies of molluscan evolution and biomineralization.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size.

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

New data from Monoplacophora and a carefully-curated dataset resolve molluscan relationships.

Relationships among the major lineages of Mollusca have long been debated. Morphological studies have considered the rarely collected Monoplacophora (Tryblidia) to have several plesiomorphic molluscan traits. The phylogenetic position of this group is contentious as morphologists have generally placed this clade as the sister taxon of the rest of Conchifera whereas earlier molecular studies supported a clade of Monoplacophora + Polyplacophora (Serialia) and phylogenomic studies have generally recovered a clade of Monoplacophora + Cephalopoda. Phylogenomic studies have also strongly supported a clade including Gastropoda, Bivalvia, and Scaphopoda, but relationships among these taxa have been inconsistent. In order to resolve conchiferan relationships and improve understanding of early molluscan evolution, we carefully curated a high-quality data matrix and conducted phylogenomic analyses with broad taxon sampling including newly sequenced genomic data from the monoplacophoran Laevipilina antarctica. Whereas a partitioned maximum likelihood (ML) analysis using site-homogeneous models recovered Monoplacophora sister to Cephalopoda with moderate support, both ML and Bayesian inference (BI) analyses using mixture models recovered Monoplacophora sister to all other conchiferans with strong support. A supertree approach also recovered Monoplacophora as the sister taxon of a clade composed of the rest of Conchifera. Gastropoda was recovered as the sister taxon of Scaphopoda in most analyses, which was strongly supported when mixture models were used. A molecular clock based on our BI topology dates diversification of Mollusca to ~546 MYA (+/- 6 MYA) and Conchifera to ~540 MYA (+/- 9 MYA), generally consistent with previous work employing nuclear housekeeping genes. These results provide important resolution of conchiferan mollusc phylogeny and offer new insights into ancestral character states of major mollusc clades.