Consensus and confusion in molluscan trees: evaluating morphological and molecular phylogenies

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum has long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species. The final assembled and filtered Berghia genome is comparable to other high quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes), and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low. Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Revealing morphological characteristics of Goniodorididae genera (Mollusca: Nudibranchia)

Detailed knowledge of the anatomy of the species is an essential element in taxonomic studies, since it allows the comparison and differentiation of separate groups of taxa. It becomes especially important when considering type species, as the subsequent identification of the species that compose the taxa is based on its characteristics, considered common in the group. However, despite its relevance, there are still numerous species without detailed descriptions, being especially significant among invertebrates. The family Goniodorididae is a little-known group of nudibranchs that includes eight recognized genera: Okenia, Goniodoris, Ancula, Lophodoris, Spahria, Trapania, Goniodoridella and Murphydoris. Several of their species are not completely described, including type species, and the systematics of the family is still unclear. Here we study in detail the external morphology and internal anatomy of the type species of five of the eight Goniodorididae genera using microcomputed tomography and scanning electron microscopy. We include the species Okenia elegans, Goniodoris nodosa, Ancula gibbosa, Goniodoridella savignyi and Murphydoris singaporensis as well as one species of Trapania, T. graeffei. We describe for the first time the detailed internal anatomy of the type species Goniodoridella savignyi. The diagnostic features of each genus are compared, and a preliminary framework is shown to clarify their systematics and identifications.

Host-Endosymbiont Genome Integration in a Deep-Sea Chemosymbiotic Clam

Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate coevolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 Gb, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene families, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 Mb but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host’s high dependence on the symbiont for nutrition. Overall, the host–symbiont genomes show not only tight metabolic complementarity but also distinct signatures of coevolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

Molluscan phylogenomics requires strategically selected genomes

The extraordinary diversity in molluscan body plans, and the genomic mechanisms that enable it, remains one of the great questions of evolution. The eight distinct living taxonomic classes of molluscs are each unambiguously monophyletic; however, significant controversy remains about the phylogenetic relationships among those eight branches. Molluscs are the second-largest animal phylum, with over 100 000 living species with broad biological, economic and medical interest. To date, only around 53 genome assemblies have been accessioned to NCBI GenBank covering only four of the eight living molluscan classes. Furthermore, the molluscan taxa where partial or whole-genome assemblies are available are often aberrantly fast evolving or recently derived lineages. Characteristic adaptations provide interesting targets for whole-genome projects, in animals like the scaly-foot snail or octopus, but without basal-branching lineages for comparison, the context of recently derived features cannot be assessed. The currently available genomes also create a non-optimal set of taxa for resolving deeper phylogenetic branches: they are a small sample representing a large group, and those that are available come primarily from a rarefied pool. Thoughtful selection of taxa for future projects should focus on the blank areas of the molluscan tree, which are ripe with opportunities to delve into peculiarities of genome evolution, and reveal the biology and evolutionary history of molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

The mitochondrial genome of an Endangered freshwater snail Koreoleptoxis nodifila (Caenogastropoda: Semisulcospiridae) from South Korea

The mitochondrial genome of the Endangered freshwater snail Koreoleptoxis nodifila (Caenogastropoda: Semisulcospiridae) from South Korea is determined and characterized in detail. It is 15,737 bp in length being composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one control region. It has a base composition of 31.23% for A, 16.29% for G, 17.84% for C, and 34.64% for T. The phylogenetic trees reconstructed based on the maximum-likelihood (ML) method and Bayesian inference (BI) confirmed that K. nodifila belongs to the Semisulcospiridae clade in the monophyletic caeonogastropod superfamily Cerithioidea.

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.

The complete mitochondrial genome of Harpago chiragra and Lambis lambis (Gastropoda: Stromboidea): implications on the Littorinimorpha phylogeny

The complete mitochondrial genomes of Harpago chiragra and Lambis lambis (Strombidae) were determined with the size of 15,460 bp and 15,481 bp, respectively, and both sequences contained 13 protein-coding genes, 22 tRNAs, and two rRNAs. H. chiragra and L. lambis have similar mitochondrial features, corresponding to typical gastropod mitochondrial genomes, such as the conserved gene order, a high A + T content (66.22% for H. chiragra and 66.10% for L. lambis), and preference for A + T-rich codons. The start or termination codon of same protein-coding gene in H. chiragra was consistent with that in L. lambis, except for the termination codon of cox1 gene (TAG for H. chiragra and TAA for L. lambis) and the start codon of nad4 (GTG for H. chiragra and ATG for L. lambis). Pairwise sequence alignments detected different degrees of variations in H. chiragra and L. lambis mitochondrial genomes; and the two species had lower levels of genetic distance (0.202 for nucleotide sequence) and closest relationships as compared to Strombus gigas and Oncomelania hupensis. The 13 partitioned nucleotide sequences of protein coding genes of H. chiragra and L. lambis were aligned with representatives of the main lineages of gastropods and their phylogenetic relationships were inferred. H. chiragra and L. lambis share the same gene order as Littorinimorpha species, except Vermetoidea, which demonstrate a gene rearrangement in species. The reconstructed phylogeny supports three major clades within Littorinimorpha: 1) Stromboidea, Tonnoidea, Littorinoidea, and Naticoidea, 2) Rissooidea and Truncatelloidea, and 3) Vermetoidea. In addition, a relaxed molecular clock calibrated with fossils dated the diversification of Strombidae near 112 (44–206) Mya and a possible radiation is detected to occur between 45–75 Mya, providing implications to understand the Cenozoic replacement event (65–135 Mya) of Aporrhaidae by Strombidae.

Evolutionary lineages of marine snails identified using molecular phylogenetics and geometric morphometric analysis of shells

The relationship between morphology and inheritance is of perennial interest in evolutionary biology and palaeontology. Using three marine snail genera Penion, Antarctoneptunea and Kelletia, we investigate whether systematics based on shell morphology accurately reflect evolutionary lineages indicated by molecular phylogenetics. Members of these gastropod genera have been a taxonomic challenge due to substantial variation in shell morphology, conservative radular and soft tissue morphology, few known ecological differences, and geographical overlap between numerous species. Sampling all sixteen putative taxa identified across the three genera, we infer mitochondrial and nuclear ribosomal DNA phylogenetic relationships within the group, and compare this to variation in adult shell shape and size. Results of phylogenetic analysis indicate that each genus is monophyletic, although the status of some phylogenetically derived and likely more recently evolved taxa within Penion is uncertain. The recently described species P. lineatus is supported by genetic evidence. Morphology, captured using geometric morphometric analysis, distinguishes the genera and matches the molecular phylogeny, although using the same dataset, species and phylogenetic subclades are not identified with high accuracy. Overall, despite abundant variation, we find that shell morphology accurately reflects genus-level classification and the corresponding deep phylogenetic splits identified in this group of marine snails.

Do chitons have a brain? New evidence for diversity and complexity in the polyplacophoran central nervous system

Molluscs demonstrate astonishing morphological diversity, and the relationships among clades have been debated for more than a century. Molluscan nervous systems range from simple 'ladder-like' cords to the complex brains of cephalopods. Chitons (Polyplacophora) are assumed to retain many molluscan plesiomorphies, lacking neural condensation and ganglionic structure, and therefore a brain. We reconstructed three-dimensional anatomical models of the nervous system in eight species of chitons in an attempt to clarify chiton neuroarchitecture and its variability. We combined new data with digitised historic slide material originally used by malacologist Johannes Thiele (1860-1935). Reconstructions of whole nervous systems in Acanthochitona fascicularis, Callochiton septemvalvis, Chiton olivaceus, Hemiarthrum setulosum, Lepidochitona cinerea, Lepidopleurus cajetanus and Leptochiton asellus, and the anterior nervous system of Schizoplax brandtii, demonstrated consistent and substantial anterior neural concentration in the circumoesophageal nerve ring. This is further organised into three concentric tracts, corresponding to the lateral, ventral and cerebral nerve cords. These represent homologues to the three main pairs of ganglia in other molluscs. Their relative size, shape and organisation are highly variable among the examined taxa, but consistent with previous studies of select species, and we formulated a set of neuroanatomical characters for chitons. These support anatomical transitions at the ordinal and subordinal levels; the identification of robust homologies in neural architecture will be central to future comparisons across Mollusca and, more broadly, Lophotrochozoa. Contrary to almost all previous descriptions, the size and structure of the chiton anterior nerve ring unambiguously qualify it as a true brain with cordal substructure.

The Performance of Two Supertree Schemes Compared Using Synthetic and Real Data Quartet Input

Despite impressive advancements in technological and theoretical tools, construction of phylogenetic (evolutionary) trees is still a challenging task. The availability of enormous quantities of molecular data has made large-scale phylogenetic reconstruction involving thousands of species, a more viable goal. For this goal, separate trees over different, overlapping subsets of species, representing histories of various markers of these species, are collected. These trees, typically with conflicting signals, are subsequently combined into a single tree over the full set, an operation denoted as supertree construction. The amalgamation of such trees into a single tree lies at the heart of many tasks in phylogenetics, yet remains a daunting endeavor, especially in light of conflicting signals. In this work, we study the performance of matrix representation with parsimony (MRP), the most widely used supertree method to date, when confronted with quartet trees. Quartet trees are the most basic informational unit when amalgamation of unrooted trees is attempted, and they remain relevant in more general settings even though standard supertree methods are not necessarily confined to quartets. This study involves both real and simulated data, and the effects of several parameters on the results are evaluated, revealing a number of anomalies associated with MRP. We show that these anomalies are surmountable when using a recently introduced supertree method, weighted quartet MaxCut (wQMC).

Salivary Glands in Predatory Mollusks: Evolutionary Considerations

Many marine mollusks attain or increase their predatory efficiency using complex chemical secretions, which are often produced and delivered through specialized anatomical structures of the foregut. The secretions produced in venom glands of Conus snails and allies have been extensively studied, revealing an amazing chemical diversity of small, highly constrained neuropeptides, whose characterization led to significant pharmacological developments. Conversely, salivary glands, the other main secretory structures of molluscan foregut, have been neglected despite their shared occurrence in the two lineages including predatory members: Gastropoda and Cephalopoda. Over the last few years, the interest for the chemistry of salivary mixtures increased based on their potential biomedical applications. Recent investigation with -omics technologies are complementing the classical biochemical descriptions, that date back to the 1950s, highlighting the high level of diversification of salivary secretions in predatory mollusks, and suggesting they can be regarded as a pharmaceutical cornucopia. As with other animal venoms, some of the salivary toxins are reported to target, for example, sodium and/or potassium ion channels or receptors and transporters for neurotransmitters such as, glutamate, serotonin, neurotensin, and noradrenaline, thus manipulating the neuromuscular system of the preys. Other bioactive components possess anticoagulant, anesthetic and hypotensive activities. Here, we overview available knowledge on the salivary glands of key predatory molluscan taxa, gastropods, and cephalopods, summarizing their anatomical, physiological and biochemical complexity in order to facilitate future comparative studies on main evolutionary trends and functional convergence in the acquisition of successful predatory strategies.

Molluscan models: Crepidula fornicata

Gastropod snails in the genus Crepidula have emerged as model systems for studying a metazoan super clade, the Spiralia. Recent work on one species in particular, Crepidula fornicata, has produced high-resolution cell lineage fate maps, details of morphogenetic events during gastrulation, key insights into the molecular underpinnings of early development, and the first demonstration of CRISPR/Cas9 genome editing in the Spiralia. Furthermore, invasive species of Crepidula are a significant ecological threat, while one of these, C. fornicata, is also being harvested for food. This review highlights progress towards developing these animals as models for evolutionary, developmental, and ecological studies. Such studies have contributed greatly to our understanding of biology in a major clade of bilaterians. This information may also help us to control and cultivate these snails.

From complex to simple: myogenesis in an aplacophoran mollusk reveals key traits in aculiferan evolution

Background

Recent studies suggest a bifurcation at the base of Mollusca, resulting in the primarily single-shelled Conchifera (Bivalvia, Gastropoda, Scaphopoda, Monoplacophora, Cephalopoda) and the spicule-bearing Aculifera (Polyplacophora, Neomeniomorpha, Chaetodermomorpha). A recent study revealed a complex larval musculature exclusively shared by Neomeniomorpha and Polyplacophora, supporting a close relationship of both taxa. However, the ontogenetic transition from the complex larval to the simple adult neomeniomorph musculature, which mainly consists of a three-layered body-wall musculature and serially iterated dorsoventral muscles, remains unknown. To close this gap in knowledge, we studied remodeling of the larval musculature during metamorphosis in the neomeniomorph Wirenia argentea. A comparative analysis with a novel data set of a polyplacophoran, Leptochiton asellus, allows us to infer the morphology of the last common ancestor of Aculifera and the evolution of its subclades therefrom.

Results

The complex larval musculature of Wirenia argentea persists through metamorphosis and becomes modified to form two of the three muscle layers of the adult body wall. The innermost longitudinal layer of the three-layered body wall musculature is generated by transformation and expansion of distinct larval longitudinal muscle bundles. The larval ventrolateral muscle strands are remodeled and eventually become the most ventral part of the adult longitudinal layer of the body wall musculature. The paired larval enrolling muscle forms the lateral parts and the former rectus muscle is destined to become the most dorsal part of the longitudinal layer of the body wall musculature. The transient ventromedian muscle is lost during postmetamorphic development.

Conclusions

Postmetamorphic remodeling in W. argentea supports the hypothesis of a complex myoanatomy rather than a three-layered body wall musculature at the base of Aculifera, and thus argues against homology of the body wall musculature of adult Neomeniomorpha and other potential molluscan sister groups. Our data show that the neomeniomorph body wall musculature is a derived condition and not an aculiferan or molluscan plesiomorphy.