Resolving the relationships of clams and cockles: dense transcriptome sampling drastically improves the bivalve tree of life

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Comparative anatomy of siphons in tellinoidean clams (Bivalvia, Tellinoidea)

Siphons are tubular organs formed by fusion and posterior extension of the marginal mantle folds. They are supposed to have performed key roles in the evolution of bivalves by enabling these animals to occupy several ecological niches. However, anatomical details of these organs are scarce for one of the most diverse lineages of tropical bivalves, the superfamily Tellinoidea. We investigated the siphonal morphology of 15 species, sampling five tellinoidean families, by integrating scanning electron microscopy, confocal microscopy, and histology. The siphons revealed variations in length, pigmentation, tentacles, papillae, and number of nerve cords. Due to the presence of sensorial structures, such as papillae and tentacles, we reclassify the siphons of Tellinoidea from type A to A+. Additional anatomical patterns were identified at family and genus levels. For example, the incurrent siphon shorter than the excurrent and 24 tentacles are putative synapomorphies of Donacidae. We also highlight shared siphonal traits between Donacidae and Solecurtidae as well as between Semelidae and Tellinidae. In addition, our data support the idea of Psammobiidae as a paraphyletic lineage. Overall, we provide an extensive comparative data set on siphonal traits with significant relevance for bivalve taxonomy, functional anatomy, and evolutionary investigations.

A new species of the genus Yoldiella (Bivalvia, Protobranchia, Yoldiidae) from Haima Cold Seep, South China Sea, China

In present study, a previously unidentified but frequently encountered species of deep-sea protobranch, Yoldiellahaimaensis sp. nov., is described new to science from the Haima Cold Seep on the northwestern slope of the South China Sea. A morphological analysis confirmed that this species belongs to a previously undescribed species of the genus Yoldiella A.E. Verrill & K.J. Bush, 1897. It differs morphologically from other known species within the genus in its shell shape, degree of inflation, beaks, and number of hinge teeth. Furthermore, we sequenced three gene segments of Y.haimaensis sp. nov., comprising a nuclear ribosomal gene (18S rRNA), a nuclear protein-coding gene (histone H3), and a mitochondrial gene (cytochrome c oxidase subunit I, COI). Our phylogenetic analysis performed on the superfamily Nuculanoidea and family Yoldiidae indicates that the genus Yoldiella is non-monophyletic, and the widely recognized families within the superfamily Nuculanoidea are also not monophyletic. Our results provide molecular insights into the Protobranchia and highlight the necessity for further samples and data to revise the classification of families and genera within the superfamily using an integrative approach that combines morphological analysis and molecular data.

The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae)

Razor clams, belonging to the Pharidae and Solenidae families, are ecologically and economically important; however, very little research has been conducted on the Pharidae family. The genus Novaculina is a marine-derived freshwater lineage, and Novaculina chinensis is a rare freshwater species of the Pharidae family. In order to understand the phylogenetic relationships of N. chinensis, we sequenced the mitochondrial genome of the genus Novaculina, which is 16,262 bp in length and consists of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The phylogenetic relationships of 69 Imparidentian mitochondrial genomes (mitogenomes) indicated that N. chineisis is closely related to Sinonovacula constricta of the order Adapedonta. Our study also found that the Ka/Ks ratios of 12 protein-coding genes in the Pharidae family are lower than one, indicating the occurrence of negative purification selection. Morphological observations of the siphons of N. chinensis, Novaculina myanmarensis, and Novaculina gangetica indicate that N. chinensis may be the ancestral clade of the genus Novaculina, which has not been proposed in previous studies. Our study provides useful molecular information on the phylogenetic and evolutionary relationships of Pharidae and also contributes to the conservation and management of the germplasm resources of N. chinensis.

Invasive mussels fashion silk-like byssus via mechanical processing of massive horizontally acquired coiled coils

Zebra and quagga mussels (Dreissena spp.) are invasive freshwater biofoulers that perpetrate devastating economic and ecological impact. Their success depends on their ability to anchor onto substrates with protein-based fibers known as byssal threads. Yet, compared to other mussel lineages, little is understood about the proteins comprising their fibers or their evolutionary history. Here, we investigated the hierarchical protein structure of Dreissenid byssal threads and the process by which they are fabricated. Unique among bivalves, we found that threads possess a predominantly β-sheet crystalline structure reminiscent of spider silk. Further analysis revealed unexpectedly that the Dreissenid thread protein precursors are mechanoresponsive α-helical proteins that are mechanically processed into β-crystallites during thread formation. Proteomic analysis of the byssus secretory organ and byssus fibers revealed a family of ultrahigh molecular weight (354 to 467 kDa) asparagine-rich (19 to 20%) protein precursors predicted to form α-helical coiled coils. Moreover, several independent lines of evidence indicate that the ancestral predecessor of these proteins was likely acquired via horizontal gene transfer. This chance evolutionary event that transpired at least 12 Mya has endowed Dreissenids with a distinctive and effective fiber formation mechanism, contributing significantly to their success as invasive species and possibly, inspiring new materials design.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

Reduced chemosymbiont genome in the methane seep thyasirid and the cooperated metabolisms in the holobiont under anaerobic sediment

Previous studies have deciphered the genomic basis of host-symbiont metabolic complementarity in vestimentiferans, bathymodioline mussels, vesicomyid clams and Alviniconcha snails, yet little is known about the chemosynthetic symbiosis in Thyasiridae-a family of Bivalvia regarded as an excellent model in chemosymbiosis research due to their wide distribution in both deep-sea and shallow-water habitats. We report the first circular thyasirid symbiont genome, named Candidatus Ruthturnera sp. Tsphm01, with a size of 1.53 Mb, 1521 coding genes and 100% completeness. Compared to its free-living relatives, Ca. Ruthturnera sp. Tsphm01 genome is reduced, lacking components for chemotaxis, citric acid cycle and de novo biosynthesis of small molecules (e.g. amino acids and cofactors), indicating it is likely an obligate intracellular symbiont. Nevertheless, the symbiont retains complete genomic components of sulphur oxidation and assimilation of inorganic carbon, and these systems were highly and actively expressed. Moreover, the symbiont appears well-adapted to anoxic environment, including capable of anaerobic respiration (i.e. reductions of DMSO and nitrate) and possession of a low oxygen-adapted type of cytochrome c oxidase. Analysis of the host transcriptome revealed its metabolic complementarity to the incomplete metabolic pathways of the symbiont and the acquisition of nutrients from the symbiont via phagocytosis and exosome. By providing the first complete genome of reduced size in a thyasirid symbiont, this study enhances our understanding of the diversity of symbiosis that has enabled bivalves to thrive in chemosynthetic habitats. The resources will be widely used in phylogenetic, geographic and evolutionary studies of chemosynthetic bacteria and bivalves.

Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes

Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00178-x.

Convergence and contingency in the evolution of a specialized mode of life: multiple origins and high disparity of rock-boring bivalves

Evolutionary adaptation to novel, specialized modes of life is often associated with a close mapping of form to the new function, resulting in narrow morphological disparity. For bivalve molluscs, endolithy (rock-boring) has biomechanical requirements thought to diverge strongly from those of ancestral functions. However, endolithy in bivalves has originated at least eight times. Three-dimensional morphometric data representing 75 species from approximately 94% of extant endolithic genera and families, along with 310 non-endolithic species in those families, show that endolithy is evolutionarily accessible from many different morphological starting points. Although some endoliths appear to converge on certain shell morphologies, the range of endolith shell form is as broad as that belonging to any other bivalve substrate use. Nevertheless, endolithy is a taxon-poor function in Bivalvia today. This limited richness does not derive from origination within source clades having significantly low origination or high extinction rates, and today's endoliths are not confined to low-diversity biogeographic regions. Instead, endolithy may be limited by habitat availability. Both determinism (as reflected by convergence among distantly related taxa) and contingency (as reflected by the endoliths that remain close to the disparate morphologies of their source clades) underlie the occupation of endolith morphospace.

Form and function of the mantle edge in Protobranchia (Mollusca: Bivalvia)

We analyzed, by optical and transmission electron microscopy, the morphology and function of the mantle edge, including the formation of the periostracum, of ten species of protobranchs. Five species from the order Nuculida, four species from the order Nuculanida and one species from the order Solemyida were studied. A second outer fold, which seems to function as a template for the internal marginal crenulations of the valves, is present in the crenulated species of Nucula. The minute non-crenulated Ennucula aegeensis shows the glandular basal cells displaced toward the periostracal groove, resembling a minute additional fold between the outer and middle folds. Intense secretion of glycocalyx, together with active uptake of particles, have been observed in the inner epithelium of the middle mantle fold and the whole epithelium of the inner mantle fold in all the studied species. Contrary to the rest of the bivalves, all the protobranchs analyzed have two basal cells involved in the formation of the external nanometric pellicle of the periostracum, a character that would support the monophyly of protobranchs. A three-layered pattern is the general rule for the periostracum in protobranchs, like for other bivalves. The presence of pouches of translucent layer inside the tanned dark layer under periostracal folds is characteristic of the species with a folded periostracum; its function is unclear but could give flexibility to the periostracum. The non-nacreous internal shell layer and the presence of translucent pouches under periostracal folds in Sarepta speciosa resemble those found in nuculanids. However, the free periostracum is rather similar to those of N. hanleyi and E. aegeensis, with a continuous vesicular layer. All the latter supports the inclusion of Sarepta in the order Nuculanida but could indicate either a basal lineage or that the translucent vesicular layer is an adaptive trait.

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.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

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

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

Calibrating phylogenies assuming bifurcation or budding alters inferred macroevolutionary dynamics in a densely sampled phylogeny of bivalve families

Analyses of evolutionary dynamics depend on how phylogenetic data are time-scaled. Most analyses of extant taxa assume a purely bifurcating model, where nodes are calibrated using the daughter lineage with the older first occurrence in the fossil record. This contrasts with budding, where nodes are calibrated using the younger first occurrence. Here, we use the extensive fossil record of bivalve molluscs for a large-scale evaluation of how branching models affect macroevolutionary analyses. We time-calibrated 91% of nodes, ranging in age from 2.59 to 485 Ma, in a phylogeny of 97 extant bivalve families. Allowing budding-based calibrations minimizes conflict between the tree and observed fossil record, and reduces the summed duration of inferred 'ghost lineages' from 6.76 billion years (Gyr; bifurcating model) to 1.00 Gyr (budding). Adding 31 extinct paraphyletic families raises ghost lineage totals to 7.86 Gyr (bifurcating) and 1.92 Gyr (budding), but incorporates more information to date divergences between lineages. Macroevolutionary analyses under a bifurcating model conflict with other palaeontological evidence on the magnitude of the end-Palaeozoic extinction, and strongly reduce Cenozoic diversification. Consideration of different branching models is essential when node-calibrating phylogenies, and for a major clade with a robust fossil record, a budding model appears more appropriate.

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.

Comparative embryo/larval sensitivity of Australian marine bivalves to ten metals: A disjunct between physiology and phylogeny

Metal contamination within the urbanized coastal zon is one threat linked to a decline in the abundance, distribution and/or species diversity of wild marine bivalve populations. This study determined the 48-h embryo/larval sensitivity (no-effect concentration (NEC) and median-effect concentration (EC50)) of ten marine bivalve species (nine endemic to Australia) to aluminium (Al), cadmium (Cd), chromium (Cr), cobalt (Co), copper (Cu), iron (Fe), lead (Pb), manganese (Mn), nickel (Ni) and zinc (Zn), key metal contaminants impacting urbanized coastal zones in south-eastern Australia, in natural seawater (20-22 °C, 30‰ salinity, pH 7.8-7.9, 1.2 mg/L dissolved organic carbon). For all metals, except Fe, the order of sensitivity was oysters > mussels ≥ scallops ≥ cockles ≥ clams, where the economically-important oysters, Magallana gigas and Saccostrea glomerata, were 2.6 (Al) to 4.2 (Cd) times more sensitive than the least sensitive clam species. For all bivalve species, the order of metal sensitivity was Cu > Pb > Zn = Ni > Co > Cd > Al > Cr(VI) > Mn ≥ Fe(III), where Cu was eight times more toxic than Zn or Ni, 28 times more toxic than Cd, 220 times more toxic than Cr(VI) and 570 times more toxic than Fe(III). Iron, unlike the other nine soluble metals, occurred as particulate Fe(III) oxyhydroxide, where EC50 values decreased with increasing exposure time as the larval (D-veliger) stage. There was no significant (p > 0.05) effect of embryo/larval mass, or surface area/volume, on metal sensitivity. Further, there was no significant (p > 0.05) relationship between metal sensitivity and phylogeny (genetic distance). Divalent metal sensitivity was positively related (r2 = 0.87) to cell surface metal-binding affinity. The current Australian marine water quality guideline for Ni is not protective of the ten bivalve species (NECs were 2-6-fold below the guideline), while the guideline for Zn is not protective of oysters.

When size matters: the first comprehensive anatomical study of a species of "Condylocardiidae", an extremely miniaturized bivalve

'Miniaturization' is a widespread phenomenon among the Metazoa. In the molluscan class Bivalvia, records of miniaturization are numerous. Among the Archiheterodonta, Warrana besnardi (Klappenbach, 1963) has attracted attention for its tiny size, which does not exceed 1.5 mm in shell length, and because it belongs to a group with limited anatomical information and often-debated status, the "Condylocardiidae" (which recent molecular studies place deeply nested within the family Carditidae). All species of Warrana Laseron, 1953 are small-bodied, and so miniaturization presumably occurred from a large-bodied ancestor within the Carditidae sensu lato. South American W. besnardi is here studied in detail. Its small size and the enlargement of the anterodorsal region during growth, reflects (and likely led) to infaunal habit, living as a burrowing bivalve that passively feeds on deposit particles entering the pallial cavity anteriorly. Mantle glands, previously reported as a common feature of other archiheterodonts, are missing in W. besnardi, but spongiform tissue in the antero-ventral portion of the mantle lobes presumably represents a blood sinus that might compensate for the great reduction of the ctenidia. Lecithotrophy is reported, with yolky oocytes bearing a thick non-cellular capsule layer; brooding was not observed, and it is here hypothesized that the extreme miniaturization, with the great reduction of ctenidia, is responsible for a shift in the reproductive mode of condylocardiids, contrasting with the commonly reported ovoviviparity of the carditids.

Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks

Doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) in bivalve mollusks is one of the most notable departures from the paradigm of strict maternal inheritance of mtDNA among metazoans. Recently, work on the Mediterranean mussel Mytilus galloprovincialis suggested that a nucleotide motif in the control region of this species, known as the sperm transmission element (STE), helps protect male-transmitted mitochondria from destruction during spermatogenesis. Subsequent studies found similar, yet divergent, STE motifs in other marine mussels. Here, we extend the in silico search for mtDNA signatures resembling known STEs. This search is carried out for the large unassigned regions of 157 complete mitochondrial genomes from within the Mytiloida, Veneroida, Unionoida, and Ostreoida bivalve orders. Based on a sliding window approach, we present evidence that there are additional putative STE signatures in the large unassigned regions of several marine clams and freshwater mussels with DUI. We discuss the implications of this finding for interpreting the origin of doubly uniparental inheritance in ancestral bivalve mollusks, as well as potential future in vitro and in silico studies that could further refine our understanding of the early evolution of this unusual system of mtDNA inheritance.

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.

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

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

Novel gene rearrangement in the mitochondrial genome of Siliqua minima (Bivalvia, Adapedonta) and phylogenetic implications for Imparidentia

Siliqua minima (Gmelin, 1791) is an important economic shellfish species belonging to the family Pharidae. To date, the complete mitochondrial genome of only one species in this family (Sinonovacula constricta) has been sequenced. Research on the Pharidae family is very limited; to improve the evolution of this bivalve family, we sequenced the complete mitochondrial genome of S. minima by next-generation sequencing. The genome is 17,064 bp in length, consisting of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), and two ribosomal RNA genes (rRNA). From the rearrangement analysis of bivalves, we found that the gene sequences of bivalves greatly variable among species, and with closer genetic relationship, the more consistent of the gene arrangement is higher among the species. Moreover, according to the gene arrangement of seven species from Adapedonta, we found that gene rearrangement among families is particularly obvious, while the gene order within families is relatively conservative. The phylogenetic analysis between species of the superorder Imparidentia using 12 conserved PCGs. The S. minima mitogenome was provided and will improve the phylogenetic resolution of Pharidae species.

Structure, function and parallel evolution of the bivalve byssus, with insights from proteomes and the zebra mussel genome

The byssus is a structure unique to bivalves. Byssal threads composed of many proteins extend like tendons from muscle cells, ending in adhesive pads that attach underwater. Crucial to settlement and metamorphosis, larvae of virtually all species are byssate. By contrast, in adults, the byssus is scattered throughout bivalves, where it has had profound effects on morphological evolution and been key to adaptive radiations of epifaunal species. I compare byssus structure and proteins in blue mussels (Mytilus), by far the best characterized, to zebra mussels (Dreissena polymorpha), in which several byssal proteins have been isolated and sequenced. By mapping the adult byssus onto a recent phylogenomic tree, I confirm its independent evolution in these and other lineages, likely parallelisms with common origins in development. While the byssus is superficially similar in Dreissena and Mytilus, in finer detail it is not, and byssal proteins are dramatically different. I used the chromosome-scale D. polymorpha genome we recently assembled to search for byssal genes and found 37 byssal loci on 10 of the 16 chromosomes. Most byssal genes are in small families, with several amino acid substitutions between paralogs. Byssal proteins of zebra mussels and related quagga mussels (D. rostriformis) are divergent, suggesting rapid evolution typical of proteins with repetitive low complexity domains. Opportunities abound for proteomic and genomic work to further our understanding of this textbook example of a marine natural material. A priority should be invasive bivalves, given the role of byssal attachment in the spread of, and ecological and economic damage caused by zebra mussels, quagga mussels and others. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

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'.

A polyvalent and universal tool for genomic studies in gastropod molluscs (Heterobranchia)

Molluscs are the second most diverse animal phylum and heterobranch gastropods present ~ 44,000 species. These comprise fascinating creatures with huge morphological and ecological disparity. Such great diversity comes with even larger phylogenetic uncertainty and many taxa have been largely neglected in molecular assessments. Genomic tools have provided resolution to deep cladogenic events but generating large numbers of transcriptomes/genomes is expensive and usually requires fresh material. Here we leverage a target enrichment approach to design and synthesize a probe set based on available genomes and transcriptomes across Heterobranchia. Our probe set contains 57,606 70mer baits and targets a total of 2,259 ultra-conserved elements (UCEs). Post-sequencing capture efficiency was tested against 31 marine heterobranchs from major groups, including Acochlidia, Acteonoidea, Aplysiida, Cephalaspidea, Pleurobranchida, Pteropoda, Runcinida, Sacoglossa, and Umbraculida. The combined Trinity and Velvet assemblies recovered up to 2,211 UCEs in Tectipleura, up to 1,978 in Nudipleura, and up to 1,927 in Acteonoidea, the latter two being the most distantly related taxa to our core study group. Total alignment length was 525,599 bp and contained 52% informative sites and 21% missing data. Maximum-likelihood and Bayesian inference approaches recovered the monophyly of all orders tested as well as the larger clades Nudipleura, Panpulmonata, and Euopisthobranchia. The successful enrichment of diversely preserved material and DNA concentrations demonstrate the polyvalent nature of UCEs, and the universality of the probe set designed. We believe this probe set will enable multiple, interesting lines of research, that will benefit from an inexpensive and largely informative tool that will, additionally, benefit from the access to museum collections to gather genomic data.

Krüppel-like factor/specificity protein evolution in the Spiralia and the implications for cephalopod visual system novelties

The cephalopod visual system is an exquisite example of convergence in biological complexity. However, we have little understanding of the genetic and molecular mechanisms underpinning its elaboration. The generation of new genetic material is considered a significant contributor to the evolution of biological novelty. We sought to understand if this mechanism may be contributing to cephalopod-specific visual system novelties. Specifically, we identified duplications in the Krüppel-like factor/specificity protein (KLF/SP) sub-family of C2H2 zinc-finger transcription factors in the squid Doryteuthis pealeii. We cloned and analysed gene expression of the KLF/SP family, including two paralogs of the DpSP6-9 gene. These duplicates showed overlapping expression domains but one paralog showed unique expression in the developing squid lens, suggesting a neofunctionalization of DpSP6-9a. To better understand this neofunctionalization, we performed a thorough phylogenetic analysis of SP6-9 orthologues in the Spiralia. We find multiple duplications and losses of the SP6-9 gene throughout spiralian lineages and at least one cephalopod-specific duplication. This work supports the hypothesis that gene duplication and neofunctionalization contribute to novel traits like the cephalopod image-forming eye and to the diversity found within Spiralia.

First transcriptome assembly of a newly discovered vent mussel, Gigantidas vrijenhoeki, at Onnuri Vent Field on the northern Central Indian Ridge

This is the first report of a transcriptome assembly of a newly discovered hydrothermal vent mussel, Gigantidas vrijenhoeki (Bivalvia: Mytilidae), on the Central Indian Ridge. Gigantidas vrijenhoeki was identified from material collected at the newly discovered Onnuri Vent Field (OVF) on the Central Indian Ridge in 2018, and was reported as a new species, distinct from another dominant hydrothermal vent mussel, Bathymodiolus marisindicus, in 2020. We sequenced the transcriptome of G. vrijenhoeki using the Illumina HiSeq X System. De novo assembly and analysis of the coding regions predicted 25,405 genes, 84.76% of which was annotated by public databases. The transcriptome of G. vrijenhoeki will be a valuable resource in studying the ecological and biological characteristics of this new species, which is distinct from other deep-sea mussels. These data should also support the investigation of the relationship between the environmental conditions of hydrothermal vents and the unique distribution of G. vrijenhoeki in the OVF of the Central Indian Ridge.

Shedding light: a phylotranscriptomic perspective illuminates the origin of photosymbiosis in marine bivalves

BACKGROUND

Photosymbiotic associations between metazoan hosts and photosynthetic dinoflagellates are crucial to the trophic and structural integrity of many marine ecosystems, including coral reefs. Although extensive efforts have been devoted to study the short-term ecological interactions between coral hosts and their symbionts, long-term evolutionary dynamics of photosymbiosis in many marine animals are not well understood. Within Bivalvia, the second largest class of mollusks, obligate photosymbiosis is found in two marine lineages: the giant clams (subfamily Tridacninae) and the heart cockles (subfamily Fraginae), both in the family Cardiidae. Morphologically, giant clams show relatively conservative shell forms whereas photosymbiotic fragines exhibit a diverse suite of anatomical adaptations including flattened shells, leafy mantle extensions, and lens-like microstructural structures. To date, the phylogenetic relationships between these two subfamilies remain poorly resolved, and it is unclear whether photosymbiosis in cardiids originated once or twice.

RESULTS

In this study, we establish a backbone phylogeny for Cardiidae utilizing RNASeq-based transcriptomic data from Tridacninae, Fraginae and other cardiids. A variety of phylogenomic approaches were used to infer the relationship between the two groups. Our analyses found conflicting gene signals and potential rapid divergence among the lineages. Overall, results support a sister group relationship between Tridacninae and Fraginae, which diverged during the Cretaceous. Although a sister group relationship is recovered, ancestral state reconstruction using maximum likelihood-based methods reveals two independent origins of photosymbiosis, one at the base of Tridacninae and the other within a symbiotic Fraginae clade.

CONCLUSIONS

The newly revealed common ancestry between Tridacninae and Fraginae brings a possibility that certain genetic, metabolic, and/or anatomical exaptations existed in their last common ancestor, which promoted both lineages to independently establish photosymbiosis, possibly in response to the modern expansion of reef habitats.

Filling a blank in bivalve taxonomy: an integrative analysis of Cyamioidea (Mollusca: Bivalvia)

This contribution provides the first integrative study of Cyamioidea, a group of bivalves largely neglected in the past. By combining information obtained from shell morphology, gross anatomy, histology, reproductive biology and DNA sequencing, the distinctive characters of Cyamioidea are determined, and the taxa belonging to this superfamily are revised. Gaimardioidea, previously considered as a distinct superfamily, is reduced to family level (Gaimardiidae), which, along with Cyamiidae, are assigned to Cyamioidea. All studied cyamioideans share a common reproductive character: the presence of a true follicular epithelium surrounding each developing oocyte, which persists after spawning, surrounding the developing embryos and participating in their anchorage to the gill filaments. Several morphological and anatomical characters support the distinction of Gaimardiidae and Cyamiidae. Based on that information, the placement of Gaimardia and Kidderia in Gaimardiidae, and that of Cyamiocardium, Cyamiomactra, Cyamium, Heteromactra, Pseudokelly, Ptychocardia and Reloncavia in Cyamiidae, is confirmed. Jukesena, historically placed in Veneroidea, is here allocated to Cyamiidae. Neoleptonidae, previously also attributed to Cyamioidea, are reallocated to Veneroidea.

Spatial filters of function and phylogeny determine morphological disparity with latitude

The drivers of latitudinal differences in the phylogenetic and ecological composition of communities are increasingly studied and understood, but still little is known about the factors underlying morphological differences. High-resolution, three-dimensional morphological data collected using computerized micro-tomography (micro-CT) allows comprehensive comparisons of morphological diversity across latitude. Using marine bivalves as a model system, this study combines 3D shape analysis (based on a new semi-automated procedure for placing landmarks and semilandmarks on shell surfaces) with non-shape traits: centroid size, proportion of shell to soft-tissue volume, and magnitude of shell ornamentation. Analyses conducted on the morphology of 95% of all marine bivalve species from two faunas along the Atlantic coast of North America, the tropical Florida Keys and the boreal Gulf of Maine, show that morphological shifts between these two faunas, and in phylogenetic and ecological subgroups shared between them, occur as changes in total variance with a bounded minimum rather than directional shifts. The dispersion of species in shell-shape morphospace is greater in the Gulf of Maine, which also shows a lower variance in ornamentation and size than the Florida Keys, but the faunas do not differ significantly in the ratio of shell to internal volume. Thus, regional differences conform to hypothesized effects of resource seasonality and predation intensity, but not to carbonate saturation or calcification costs. The overall morphological differences between the regional faunas is largely driven by the loss of ecological functional groups and family-level clades at high latitudes, rather than directional shifts in morphology within the shared groups with latitude. Latitudinal differences in morphology thus represent a complex integration of phylogenetic and ecological factors that are best captured in multivariate analyses across several hierarchical levels.