Revisiting the evolution of Family B1 GPCRs and ligands: Insights from mollusca

Family B1 G protein-coupled receptors (GPCRs) are one of the most well studied neuropeptide receptor families since they play a central role in many biological processes including endocrine, gastrointestinal, cardiovascular and reproduction in animals. The genes for these receptors emerged from a common ancestral gene in bilaterian genomes and evolved via gene/genome duplications and deletions in vertebrate and invertebrate genomes. Their existence and function have mostly been characterized in vertebrates and few studies exist in invertebrate species. Recently, an increased interest in molluscs, means a series of genomes have become available, and since they are less modified than insect and nematode genomes, they are ideal to explore the origin and evolution of neuropeptide gene families. This review provides an overview of Family B1 GPCRs and their peptide ligands and incorporates new data obtained from Mollusca genomes and taking a comparative approach challenges existing models on their origin and evolution.

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.

The genome of the rayed Mediterranean limpet Patella caerulea (Linnaeus, 1758)

Patella caerulea (Linnaeus, 1758) is a mollusc limpet species of the class Gastropoda. Endemic to the Mediterranean Sea, it is considered a keystone species due to its primary role in structuring and regulating the ecological balance of tidal and subtidal habitats. It is currently being used as a bioindicator to assess the environmental quality of coastal marine waters and as a model species to understand adaptation to ocean acidification. Here, we provide a high-quality reference genome assembly and annotation for P. caerulea. We generated ∼30 Gb of Pacific Biosciences high-fidelity data from a single individual and provide a final 749.8 Mb assembly containing 62 contigs, including the mitochondrial genome (14,938 bp). With an N50 of 48.8 Mb and 98% of the assembly contained in the 18 largest contigs, this assembly is near chromosome-scale. Benchmarking Universal Single-Copy Orthologs scores were high (Mollusca, 87.8% complete; Metazoa, 97.2% complete) and similar to metrics observed for other chromosome-level Patella genomes, highlighting a possible bias in the Mollusca database for Patellids. We generated transcriptomic Illumina data from a second individual collected at the same locality and used it together with protein evidence to annotate the genome. A total of 23,938 protein-coding gene models were found. By comparing this annotation with other published Patella annotations, we found that the distribution and median values of exon and gene lengths was comparable with other Patella species despite different annotation approaches. The present high-quality P. caerulea reference genome, available on GenBank (BioProject: PRJNA1045377; assembly: GCA_036850965.1), is an important resource for future ecological and evolutionary studies.

Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

A combined phylogenetic strategy illuminates the evolution of Goniodorididae nudibranchs (Mollusca, Gastropoda, Heterobranchia)

Goniodorididae is a family of small dorid nudibranchs distributed worldwide that feed on entoprocts, ascidians, and bryozoans. The evolutionary relationships between its taxa have been uncertain due to the limited taxa available for phylogenetic analyses; some genera being paraphyletic. The family includes a remarkable number of synonymized genera in which the species richness is unequally distributed, while some genera have dozens of species others are monospecific. Some clades are very uniform morphologically while others are considered highly variable. To increase backbone phylogenetic resolution a target enrichment approach of ultra-conserved elements was aimed at representative Goniodorididae species for the first time. Additionally, we increase species representation by including mitochondrial markers cytochrome c oxidase subunit I and ribosomal RNA 16S as well as nuclear Histone 3 and ribosomal RNA 18S from 109 Goniodorididae species, out of approximately 160 currently valid species. Maximum likelihood and Bayesian inference analyses were performed to infer the phylogeny of the family. As a result, two subfamilies and eleven genera were elucidated. The synonymized genera Bermudella, Cargoa, and Ceratodoris are here resurrected and a new genus, Naisdoris gen. nov., is described. The clades included taxa with shared prey preference, showing that trophic behavior could have driven species evolution and morphological uniqueness within the family Goniodorididae.

A chromosome-level genome assembly of a deep-sea symbiotic Aplacophora mollusc Chaetoderma sp

The worm-shaped, shell-less Caudofoveata is one of the least known groups of molluscs. As early-branching molluscs, the lack of high-quality genomes hinders our understanding of their evolution and ecology. Here, we report a high-quality chromosome-scale genome of Chaetoderma sp. combining PacBio, Illumina, and high-resolution chromosome conformation capture sequencing. The final assembly has a size of 2.45 Gb, with a scaffold N50 length of 141.46 Mb, and is anchored to 17 chromosomes. Gene annotations showed a high level of accuracy and completeness, with 23,675 predicted protein-coding genes and 94.44% of the metazoan conserved genes by BUSCO assessment. We further present 16S rRNA gene amplicon sequencing of the gut microbiota in Chaetoderma sp., which was dominated by the chemoautotrophic bacteria (phylum Gammaproteobacteria). This chromosome-level genome assembly presents the first genome for the Caudofoveata, which constitutes an important resource for studies ranging from molluscan evolution, symposium, to deep-sea adaptation.

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.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Molluscan Genomes Reveal Extensive Differences in Photopigment Evolution Across the Phylum

In animals, opsins and cryptochromes are major protein families that transduce light signals when bound to light-absorbing chromophores. Opsins are involved in various light-dependent processes, like vision, and have been co-opted for light-independent sensory modalities. Cryptochromes are important photoreceptors in animals, generally regulating circadian rhythm, they belong to a larger protein family with photolyases, which repair UV-induced DNA damage. Mollusks are great animals to explore questions about light sensing as eyes have evolved multiple times across, and within, taxonomic classes. We used molluscan genome assemblies from 80 species to predict protein sequences and examine gene family evolution using phylogenetic approaches. We found extensive opsin family expansion and contraction, particularly in bivalve xenopsins and gastropod Go-opsins, while other opsins, like retinochrome, rarely duplicate. Bivalve and gastropod lineages exhibit fluctuations in opsin repertoire, with cephalopods having the fewest number of opsins and loss of at least 2 major opsin types. Interestingly, opsin expansions are not limited to eyed species, and the highest opsin content was seen in eyeless bivalves. The dynamic nature of opsin evolution is quite contrary to the general lack of diversification in mollusk cryptochromes, though some taxa, including cephalopods and terrestrial gastropods, have reduced repertoires of both protein families. We also found complete loss of opsins and cryptochromes in multiple, but not all, deep-sea species. These results help set the stage for connecting genomic changes, including opsin family expansion and contraction, with differences in environmental, and biological features across Mollusca.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm)

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

Genomic Insights into Mollusk Terrestrialization: Parallel and Convergent Gene Family Expansions as Key Facilitators in Out-of-the-Sea Transitions

Animals abandoned their marine niche and successfully adapted to life on land multiple times throughout evolution, providing a rare opportunity to study the mechanisms driving large scale macroevolutionary convergence. However, the genomic factors underlying this process remain largely unknown. Here, we investigate the macroevolutionary dynamics of gene repertoire evolution during repeated transitions out of the sea in mollusks, a lineage that has transitioned to freshwater and terrestrial environments multiple independent times. Through phylogenomics and phylogenetic comparative methods, we examine ∼100 genomic data sets encompassing all major molluskan lineages. We introduce a conceptual framework for identifying and analyzing parallel and convergent evolution at the orthogroup level (groups of genes derived from a single ancestral gene in the species in question) and explore the extent of these mechanisms. Despite deep temporal divergences, we found that parallel expansions of ancient gene families played a major role in facilitating adaptation to nonmarine habitats, highlighting the relevance of the preexisting genomic toolkit in facilitating adaptation to new environments. The expanded functions primarily involve metabolic, osmoregulatory, and defense-related systems. We further found functionally convergent lineage-exclusive gene gains, while family contractions appear to be driven by neutral processes. Also, genomic innovations likely contributed to fuel independent habitat transitions. Overall, our study reveals that various mechanisms of gene repertoire evolution-parallelism, convergence, and innovation-can simultaneously contribute to major evolutionary transitions. Our results provide a genome-wide gene repertoire atlas of molluskan terrestrialization that paves the way toward further understanding the functional and evolutionary bases of this process.

Role of maternal spiralian-specific homeobox gene SPILE-E in the specification of blastomeres along the animal-vegetal axis during the early cleavage stages of mollusks

Spiralians, one of the major clades of bilaterians, share a unique development known as spiralian development, characterized by the formation of tiers of cells called quartets, which exhibit different developmental potentials along the animal-vegetal axis. Recently, spiralian-specific TALE-type homeobox genes (SPILE) have been identified, some of which show zygotic and staggered expression patterns along the animal-vegetal axis and function in quartet specification in mollusks. However, it is unclear which maternal molecular components control the zygotic expression of these transcription factors. In this study, we focused on SPILE-E, a maternal transcription factor, and investigated its expression and function in mollusks. We found that the maternal and ubiquitous expression of SPILE-E in the cleavage stages is conserved in molluskan species, including limpets, mussels, and chitons. We knocked down SPILE-E in limpets and revealed that the expression of transcription factors specifically expressed in the first quartet (1q2 ; foxj1b) and second quartet (2q; SPILE-B) was abolished, whereas the macromere-quartet marker (SPILE-C) was ectopically expressed in 1q2 in SPILE-E morphants. Moreover, we showed that the expression of SPILE-A, which upregulates SPILE-B but represses SPILE-C expression, decreased in SPILE-E morphants. Consistent with changes in the expression pattern of the above transcription factors, SPILE-E-morphant larvae exhibited patchy or complete loss of expression of marker genes of ciliated cells and shell fields, possibly reflecting incomplete specification of 1q2 and 2q. Our results provide a molecular framework for quartet specification and highlight the importance of maternal lineage-specific transcription factors in the development and evolution of spiralians.

Mosquito (MS), a DD37E Family of Tc1/Mariner, Displaying a Distinct Evolution Profile from DD37E/TRT and DD37E/L18

Diverse Tc1/mariner elements with the DD37E signature have been detected. However, their evolutionary relationship and profiles are largely unknown. Using bioinformatics methods, we defined the evolution profile of a Tc1/Mariner family, which harbors the catalytic domain with the DD37E signature, and renamed it DD37E/Mosquito (MS). MS transposons form a separate monophyletic clade in the phylogenetic tree, distinct from the other two groups of elements with the DD37E signature, DD37E/L18 and DD37E/TRT (transposon related to Tc1), and represent a very different taxonomic distribution from that of DD37E/TRT. MS is only detected in invertebrate and is mostly present in Arthropoda, as well as in Cnidaria, Ctenophora, Mollusca, Nematoda, and Platyhelminthes, with a total length of about 1.3 kb, containing an open reading frame (ORF) encoding about 340 amino acids transposases, with a conserved DD37E catalytic domain. The terminal inverted repeat (TIR) lengths range from 19 bp to 203 bp, and the target site duplication (TSD) is TA. We also identified few occurrences of MS horizontal transfers (HT) across lineages of diptera. In this paper, the distribution characteristics, structural characteristics, phylogenetic evolution, and horizontal transfer of the MS family are fully analyzed, which is conducive to supplementing and improving the Tc1/Mariner superfamily and excavating active transposons.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

The Expression of miRNAs Involved in Long-Term Memory Formation in the CNS of the Mollusk Helix lucorum

Mollusks are unique animals with a relatively simple central nervous system (CNS) containing giant neurons with identified functions. With such simple CNS, mollusks yet display sufficiently complex behavior, thus ideal for various studies of behavioral processes, including long-term memory (LTM) formation. For our research, we use the formation of the fear avoidance reflex in the terrestrial mollusk Helix lucorum as a learning model. We have shown previously that LTM formation in Helix requires epigenetic modifications of histones leading to both activation and inactivation of the specific genes. It is known that microRNAs (miRNAs) negatively regulate the expression of genes; however, the role of miRNAs in behavioral regulation has been poorly investigated. Currently, there is no miRNAs sequencing data being published on Helix lucorum, which makes it impossible to investigate the role of miRNAs in the memory formation of this mollusk. In this study, we have performed sequencing and comparative bioinformatics analysis of the miRNAs from the CNS of Helix lucorum. We have identified 95 different microRNAs, including microRNAs belonging to the MIR-9, MIR-10, MIR-22, MIR-124, MIR-137, and MIR-153 families, known to be involved in various CNS processes of vertebrates and other species, particularly, in the fear behavior and LTM. We have shown that in the CNS of Helix lucorum MIR-10 family (26 miRNAs) is the most representative one, including Hlu-Mir-10-S5-5p and Hlu-Mir-10-S9-5p as top hits. Moreover, we have shown the involvement of the MIR-10 family in LTM formation in Helix. The expression of 17 representatives of MIR-10 differentially changes during different periods of LTM consolidation in the CNS of Helix. In addition, using comparative analysis of microRNA expression upon learning in normal snails and snails with deficient learning abilities with dysfunction of the serotonergic system, we identified a number of microRNAs from several families, including MIR-10, which expression changes only in normal animals. The obtained data can be used for further fundamental and applied behavioral research.

Exploring the Potential of Metatranscriptomics to Describe Microbial Communities and Their Effects in Molluscs

Metatranscriptomics has emerged as a very useful technology for the study of microbiomes from RNA-seq reads. This method provides additional information compared to the sequencing of ribosomal genes because the gene expression can also be analysed. In this work, we used the metatranscriptomic approach to study the whole microbiome of mussels, including bacteria, viruses, fungi, and protozoans, by mapping the RNA-seq reads to custom assembly databases (including the genomes of microorganisms publicly available). This strategy allowed us not only to describe the diversity of microorganisms but also to relate the host transcriptome and microbiome, finding the genes more affected by the pathogen load. Although some bacteria abundant in the metatranscriptomic analysis were undetectable by 16S rRNA sequencing, a common core of the taxa was detected by both methodologies (62% of the metatranscriptomic detections were also identified by 16S rRNA sequencing, the Oceanospirillales, Flavobacteriales and Vibrionales orders being the most relevant). However, the differences in the microbiome composition were observed among different tissues of Mytilus galloprovincialis, with the fungal kingdom being especially diverse, or among molluscan species. These results confirm the potential of a meta-analysis of transcriptome data to obtain new information on the molluscs' microbiome.

Structural and evolutionary insights into the multidomain galectin from the red abalone Haliotis rufescens with specificity for sulfated glycans

Galectins are an evolutionarily ancient family of lectins characterized by their affinity for β-galactosides and a conserved binding site in the carbohydrate recognition domain (CRD). These lectins are involved in multiple physiological functions, including the recognition of glycans on the surface of viruses and bacteria. This feature supports their role in innate immune responses in marine mollusks. Here, we identified and characterized a galectin, from the mollusk Haliotis rufescens (named HrGal), with four CRDs that belong to the tandem-repeat type. HrGal was purified by affinity chromatography in a galactose-agarose resin and exhibited a molecular mass of 64.11 kDa determined by MALDI-TOF mass spectrometry. The identity of HrGal was verified by sequencing, confirming that it is a 555 amino acid protein with a mass of 63.86 kDa. This protein corresponds to a galectin reported in GenBank with accession number AHX26603. HrGal is stable in the presence of urea, reducing agents, and ions such as Cu²⁺ and Zn²⁺. The recombinant galectin (rHrGal) was purified from inclusion bodies in the presence of these ions. A theoretical model obtained with the AlphaFold server exhibits four non-identical CRDs, with a β sandwich folding and the representative motifs for binding β-galactosides. This allows us to classify HrGal within the tandem repeat galectin family. On the basis of a phylogenetic analysis, we found that the mollusk sequences form a monophyletic group of tetradomain galectins unrelated to vertebrate galectins. HrGal showed specificity for galactosides and glucosides but only the sulfated sugars heparin and ι-carrageenan inhibited its hemagglutinating activity with a minimum inhibitory concentration of 4 mM and 6.25 X 10^-5% respectively. The position of the sulfate groups seemed crucial for binding, both by carrageenans and heparin.

The Complete Mitochondrial Genome of Entemnotrochus rumphii, a Living Fossil for Vetigastropoda (Mollusca: Gastropoda)

Pleurotomarioidea represents a truly isolated and basally diverging lineage in Vetigastropoda (Mollusca: Gastropoda) whose fossil record can date back to the late Cambrian, thus providing rare insights into the evolutionary history of molluscs. Here, we sequenced and assembled the complete mitochondrial genome of one representative species from Pleurotomarioidea-Entemnotrochus rumphii (Schepman, 1879)-of which the mitogenome is 15,795 bp in length, including 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes. The nucleotide composition was biased toward AT, and A + T content reached 65.2%. E. rumphii was recovered as sister to all other living vetigastropods according to mitogenome-based phylogenetic analysis. The mitochondrial gene order was consistent with major vetigastropods and the hypothetical ancestral gastropoda, suggesting the deep conservation of mitogenome arrangement in Vetigastropoda.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Convergent evolution of barnacles and molluscs sheds lights in origin and diversification of calcareous shell and sessile lifestyle

The calcareous shell and sessile lifestyle are the representative phenotypes of many molluscs, which happen to be present in barnacles, a group of unique crustaceans. The origin of these phenotypes is unclear, but it may be embodied in the convergent genetics of such distant groups (interphylum). Herein, we perform comprehensive comparative genomics analysis in barnacles and molluscs, and reveal a genome-wide strong convergent molecular evolution between them, including coexpansion of biomineralization and organic matrix genes for shell formation, and origination of lineage-specific orphan genes for settlement. Notably, the expanded biomineralization gene encoding alkaline phosphatase evolves a novel, highly conserved motif that may trigger the origin of barnacle shell formation. Unlike molluscs, barnacles adopt novel organic matrices and cement proteins for shell formation and settlement, respectively, and their calcareous shells have potentially originated from the cuticle system of crustaceans. Therefore, our study corroborates the idea that selection pressures driving convergent evolution may strongly act in organisms inhabiting similar environments regardless of phylogenetic distance. The convergence signatures shed light on the origin of the shell and sessile lifestyle of barnacles and molluscs. In addition, notable non-convergence signatures are also present and may contribute to morphological and functional specificities.

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.

Optimized Sensory Units Integrated in the Chiton Shell

The first step for animals to interact with external environment is to sense. Unlike vertebrate animals with flexibility, it is challenging for ancient animals that are less flexible especially for mollusca with heavy shells. Chiton, as an example, has eight overlapping shells covering almost the whole body, is known to incorporate sensory units called aesthetes inside the shell. We used micro-computed tomography combined with quantitative image analysis to reveal the optimized shell geometry to resist force and the aesthetes' global distribution at the whole animal levels to facilitate sense from diverse directions both in the seawater and air. Additionally, shell proteomics combined with transcriptome reveals shell matrix proteins responsible for shell construction and potentially sensory function, highlighting unique cadherin-related proteins among mollusca. Together, this multi-level evidence of sensory units in the chiton shell may shed light on the formation of chiton shells and inspire the design of hard armor with sensory function.

Multiple tachykinins and their receptors characterized in the gastropod mollusk Pacific abalone: Expression, signaling cascades, and potential role in regulating lipid metabolism

Tachykinin (TK) families, including the first neuropeptide substance P, have been intensively explored in bilaterians. Knowledge of signaling of TK receptors (TKRs) has enabled the comprehension of diverse physiological processes. However, TK signaling systems are largely unknown in Lophotrochozoa. This study identified two TK precursors and two TKR isoforms in the Pacific abalone Haliotis discus hannai (Hdh), and characterized Hdh-TK signaling. Hdh-TK peptides harbored protostomian TK-specific FXGXRamide or unique YXGXRamide motifs at the C-termini. A phylogenetic analysis showed that lophotrochozoan TKRs, including Hdh-TKRs, form a monophyletic group distinct from arthropod TKRs and natalisin receptor groups. Although reporter assays demonstrated that all examined Hdh-TK peptides activate intracellular cAMP accumulation and Ca²⁺ mobilization in Hdh-TKR-expressing mammalian cells, Hdh-TK peptides with N-terminal aromatic residues and C-terminal FXGXRamide motifs were more active than shorter or less aromatic Hdh-TK peptides with a C-terminal YXGXRamide. In addition, we showed that ligand-stimulated Hdh-TKRs mediate ERK1/2 phosphorylation in HEK293 cells and that ERK1/2 phosphorylation is inhibited by PKA and PKC inhibitors. In three-dimensional in silico Hdh-TKR binding modeling, higher docking scores of Hdh-TK peptides were consistent with the lower EC50 values in the reporter assays. The transcripts for Hdh-TK precursors and Hdh-TKR were highly expressed in the neural ganglia, with lower expression levels in peripheral tissues. When abalone were starved for 3 weeks, Hdh-TK1 transcript levels, but not Hdh-TK2, were increased in the cerebral ganglia (CG), intestine, and hepatopancreas, contrasting with the decreased lipid content and transcript levels of sterol regulatory element-binding protein (SREBP). At 24 h post-injection in vivo, the lower dose of Hdh-TK1 mixture increased SREBP transcript levels in the CG and hepatopancreas and accumulative food consumption of abalone. Higher doses of Hdh-TK1 and Hdh-TK2 mixtures decreased the SREBP levels in the CG. When Hdh-TK2-specific siRNA was injected into abalone, intestinal SREBP levels were significantly increased, whereas administration of both Hdh-TK1 and Hdh-TK2 siRNA led to decreased SREBP expression in the CG. Collectively, our results demonstrate the first TK signaling system in gastropod mollusks and suggest a possible role for TK peptides in regulating lipid metabolism in the neural and peripheral tissues of abalone.

Identification of peptidoglycan recognition proteins in hemocytes and kidney of common periwinkle Littorinalittorea

Peptidoglycan Recognition Proteins (PGRPs) are a diverse group of proteins involved in innate immunity. In particular, PGRPs have been shown to participate in immune pattern recognition in various mollusks. However, they have not been described in Caenogastropoda, a large molluscan group comprising sea, freshwater and land snails. In this study, four short PGRPs with molecular weights ranging from 21 to 34 kDa and their isoforms were identified and structurally characterized in the kidney and hemocytic transcriptomes of a caenogastropod mollusk Littorina littorea. All of them (LlPGRP1-4) are secretory, possess a signal peptide and a characteristic N-terminal N-acetylmuramoyl-l-alanine amidase (Ami) domain with conserved Zn²⁺ binding- and amidase catalytic sites. The shortest proteins, LlPGRP1 and LlPGRP2, have no additional conserved motifs on the N-terminus. In longer and most abundantly expressed LlPGRP3 and LlPGRP4 the Ami-domain is combined with an N-terminal SH3-domain and a cysteine-rich motif, respectively. Expression analysis showed that LlPGRPs of the common periwinkle were uninvolved in the immune response to infection with trematode Himasthla elongata though they might act in antibacterial defense.

67,000 years of coastal engagement at Panga ya Saidi, eastern Africa

The antiquity and nature of coastal resource procurement is central to understanding human evolution and adaptations to complex environments. It has become increasingly apparent in global archaeological studies that the timing, characteristics, and trajectories of coastal resource use are highly variable. Within Africa, discussions of these issues have largely been based on the archaeological record from the south and northeast of the continent, with little evidence from eastern coastal areas leaving significant spatial and temporal gaps in our knowledge. Here, we present data from Panga ya Saidi, a limestone cave complex located 15 km from the modern Kenyan coast, which represents the first long-term sequence of coastal engagement from eastern Africa. Rather than attempting to distinguish between coastal resource use and coastal adaptations, we focus on coastal engagement as a means of characterising human relationships with marine environments and resources from this inland location. We use aquatic mollusc data spanning the past 67,000 years to document shifts in the acquisition, transportation, and discard of these materials, to better understand long-term trends in coastal engagement. Our results show pulses of coastal engagement beginning with low-intensity symbolism, and culminating in the consistent low-level transport of marine and freshwater food resources, emphasising a diverse relationship through time. Panga ya Saidi has the oldest stratified evidence of marine engagement in eastern Africa, and is the only site in Africa which documents coastal resources from the Late Pleistocene through the Holocene, highlighting the potential archaeological importance of peri-coastal sites to debates about marine resource relationships.

DNA barcoding of brackish and marine water fishes and shellfishes of Sundarbans, the world's largest mangrove ecosystem

The present study aims to apply a DNA barcoding tool through amplifying two mitochondrial candidate genes i.e., COI and 16S rRNA for accurate identification of fish, aquatic molluscs and crustaceans of Sundarbans mangrove wetland, to build a reference library of fish and shellfishes of this unique ecosystems. A total of 185 mitochondrial COI barcode sequences and 59 partial sequences of the 16S rRNA gene were obtained from 120 genera, 65 families and 21 orders of fish, crustaceans and molluscs. The collected samples were first identified by examining morphometric characteristics and then assessed by DNA barcoding. The COI and 16S rRNA sequences of fishes and crustaceans were clearly discriminated among genera in their phylogenies. The average Kimura two-parameter (K2P) distances of COI barcode sequences within species, genera, and families of fishes are 1.57±0.06%, 15.16±0.23%, and 17.79±0.02%, respectively, and for 16S rRNA sequences, these values are 1.74±.8%, 0.97±.8%, and 4.29±1.3%, respectively. The minimum and maximum K2P distance based divergences in COI sequences of fishes are 0.19% and 36.27%, respectively. In crustaceans, the K2P distances within genera, families, and orders are 1.4±0.03%, 17.73±0.15%, and 22.81±0.02%, respectively and the minimum and maximum divergences are 0.2% and 33.93%, respectively. Additionally, the present study resolves the misidentification of the mud crab species of the Sundarbans as Scylla olivacea which was previously stated as Scylla serrata. In case of molluscs, values of interspecific divergence ranges from 17.43% to 66.3% in the barcoded species. The present study describes the development of a molecular and morphometric cross-referenced inventory of fish and shellfish of the Sundarbans. This inventory will be useful in future biodiversity studies and in forming future conservation plan.

The Enigmatic Metallothioneins: A Case of Upward-Looking Research

In the mid-1950s, Bert Lester Vallee and his colleague Marvin Margoshes discovered a molecule referred to today as metallothionein (MT). Meanwhile, MTs have been shown to be common in many biological organisms. Despite their prevalence, however, it remains unclear to date what exactly MTs do and how they contribute to the biological function of an organism or organ. We investigate why biochemical research has not yet been able to pinpoint the function(s) of MTs. We shall systematically examine both the discovery of and recent research on Dr. Vallee's beloved family of MT proteins utilizing tools from philosophy of science. Our analysis highlights that Vallee's initial work exhibited features prototypical of a developing research tradition: it was upward-looking, exploratory, and utilized mere interactions. Since the 1960s, MT research has increasingly become intervention- and hypothesis-based while it remained largely upward-looking in character. Whilst there is no reason to think that upward-looking research cannot successfully yield structure-function mappings, it has not yet been successful in the case of MTs. Thus, we suggest it might be time to change track and consider other research strategies looking into the evolution of MTs. Recent studies in mollusks render research in this direction worthy of pursuit.

Single individual structural variant detection uncovers widespread hemizygosity in molluscs

The advent of complete genomic sequencing has opened a window into genomic phenomena obscured by fragmented assemblies. A good example of these is the existence of hemizygous regions of autosomal chromosomes, which can result in marked differences in gene content between individuals within species. While these hemizygous regions, and presence/absence variation of genes that can result, are well known in plants, firm evidence has only recently emerged for their existence in metazoans. Here, we use recently published, complete genomes from wild-caught molluscs to investigate the prevalence of hemizygosity across a well-known and ecologically important clade. We show that hemizygous regions are widespread in mollusc genomes, not clustered in individual chromosomes, and often contain genes linked to transposition, DNA repair and stress response. With targeted investigations of HSP70-12 and C1qDC, we also show how individual gene families are distributed within pan-genomes. This work suggests that extensive pan-genomes are widespread across the conchiferan Mollusca, and represent useful tools for genomic evolution, allowing the maintenance of additional genetic diversity within the population. As genomic sequencing and re-sequencing becomes more routine, the prevalence of hemizygosity, and its impact on selection and adaptation, are key targets for research across the tree of life. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Potential of genomic technologies to improve disease resistance in molluscan aquaculture

Molluscan aquaculture is a major contributor to global seafood production, but is hampered by infectious disease outbreaks that can cause serious economic losses. Selective breeding has been widely used to improve disease resistance in major agricultural and aquaculture species, and has clear potential in molluscs, albeit its commercial application remains at a formative stage. Advances in genomic technologies, especially the development of cost-efficient genomic selection, have the potential to accelerate genetic improvement. However, tailored approaches are required owing to the distinctive reproductive and life cycle characteristics of molluscan species. Transgenesis and genome editing, in particular CRISPR/Cas systems, have been successfully trialled in molluscs and may further understanding and improvement of genetic resistance to disease through targeted changes to the host genome. Whole-organism genome editing is achievable on a much greater scale compared to other farmed species, making genome-wide CRISPR screening approaches plausible. This review discusses the current state and future potential of selective breeding, genomic tools and genome editing approaches to understand and improve host resistance to infectious disease in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Mobilizing molluscan models and genomes in biology

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

The primitive interferon-like system and its antiviral function in molluscs

The phylum mollusca is a very important group in the animal kingdom for the large number and diversified species. Recently, interest in molluscan immunity has increased due to their phylogenetic position and importance in worldwide aquaculture and aquatic environment. As the main aquaculture animal, most molluscs live in the water environment and they have to cope with many pathogen challenges, in which virus is one of the primary causes for the mass mortality. In vertebrates, interferon (IFN) system is generally recognized as the first line of defence against viral infection, while the antiviral mechanisms in molluscs remain to be clearly illuminated. Recently, some IFN-like proteins and IFN-related components have been characterized from molluscs, such as pattern recognition receptors (PRRs), interferon regulatory factors (IRFs), IFN-like receptors, JAK/STAT and IFN-stimulated genes (ISGs), which reinforce the existence of IFN-like system in molluscs. This system can be activated by virus or poly (I:C) challenges and further regulate the antiviral response of haemocytes in molluscs. This review summarizes the research progresses of IFN-like system in molluscs with the emphases on the uniformity and heterogeneity of IFN-like system of molluscs compared to that of other animals, which will be helpful for elucidating the antiviral modulation in molluscs and understanding the origin and evolution of IFN system.

Expression Pattern of Nitric Oxide Synthase during Development of the Marine Gastropod Mollusc, Crepidula fornicata

Nitric Oxide (NO) plays a key role in the induction of larval metamorphosis in several invertebrate phyla. The inhibition of the NO synthase in Crepidula fornicata, a molluscan model for evolutionary, developmental, and ecological research, has been demonstrated to block the initiation of metamorphosis highlighting that endogenous NO is crucial in the control of this developmental and morphological process. Nitric Oxide Synthase contributes to the development of shell gland, digestive gland and kidney, being expressed in cells that presumably correspond to FMRF-amide, serotoninergic and catecolaminergic neurons. Here we identified a single Nos gene in embryonic and larval transcriptomes of C. fornicata and studied its localization during development, through whole-mount in situ hybridization, in order to compare its expression pattern with that of other marine invertebrate animal models.

Non-collinear Hox gene expression in bivalves and the evolution of morphological novelties in mollusks

Hox genes are key developmental regulators that are involved in establishing morphological features during animal ontogeny. They are commonly expressed along the anterior-posterior axis in a staggered, or collinear, fashion. In mollusks, the repertoire of body plans is widely diverse and current data suggest their involvement during development of landmark morphological traits in Conchifera, one of the two major lineages that comprises those taxa that originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia). For most clades, and bivalves in particular, data on Hox gene expression throughout ontogeny are scarce. We thus investigated Hox expression during development of the quagga mussel, Dreissena rostriformis, to elucidate to which degree they might contribute to specific phenotypic traits as in other conchiferans. The Hox/ParaHox complement of Mollusca typically comprises 14 genes, 13 of which are present in bivalve genomes including Dreissena. We describe here expression of 9 Hox genes and the ParaHox gene Xlox during Dreissena development. Hox expression in Dreissena is first detected in the gastrula stage with widely overlapping expression domains of most genes. In the trochophore stage, Hox gene expression shifts towards more compact, largely mesodermal domains. Only few of these domains can be assigned to specific developing morphological structures such as Hox1 in the shell field and Xlox in the hindgut. We did not find traces of spatial or temporal staggered expression of Hox genes in Dreissena. Our data support the notion that Hox gene expression has been coopted independently, and to varying degrees, into lineage-specific structures in the respective conchiferan clades. The non-collinear mode of Hox expression in Dreissena might be a result of the low degree of body plan regionalization along the bivalve anterior-posterior axis as exemplified by the lack of key morphological traits such as a distinct head, cephalic tentacles, radula apparatus, and a simplified central nervous system.

Use of a gene encoding zona pellucida 4 as a female-specific marker for early stage sexual differentiation and size dimorphism in the pacific abalone Haliotis discus hannai

Growth rates of Pacific abalone Haliotis discus hannai are an important trait affecting the economic value in the abalone aquaculture industry. A reverse-transcription polymerase chain reaction (RT-PCR) analyses of tissues from H. discus hannai was conducted for sexually mature gonads to determine male- and female-specific target gene expression, including genes encoding zona pellucida domain 4 (zp4), sperm protein (sp) and lysin (lys), respectively. Sex-specific expression patterns of these gene expression, even in sexually immature abalone, indicate these genes can be used as sensitive and robust sex-specific molecular markers. The RT-PCR procedure was also performed to analyze tissues collected at various developmental stages (50-day intervals) beginning at fertilization to determine when sex differentiation and expression of sex-specific genes was initiated. Detection of zp4 transcript in tissues collected at 200 days post-fertilization (dpf) indicated egg-specific development starts at 150-200 dpf. To evaluate possible sex-specific differences in growth rate, there was conducting of a molecular marker-based sex identification of abalone from a population selected for rapid growth rate. In a group of large H. discus hannai, females were more prevalent than males. To assess the correlation between growth and sex, there was comparison of weights of 3-year-old Pacific abalone in specimens where there had been sex determinations by visual examination and molecular methods. The results indicated females weighed more (55.92 ± 9.38 g, n = 15) than males (43.64 ± 15.55 g, n = 6, P =  0.037), indicating females had a more rapid growth rate than males.

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.

The functional roles of the non-coding RNAs in molluscs

This review focus on the current knowledge of non-coding RNAs (ncRNAs) in molluscs. In this review, we provide an overview of long ncRNAs (lncRNA), microRNAs (miRNA) and piwi-interacting RNAs (piRNA), followed by evidence for the regulation of ncRNAs in variety of biological process in molluscs, including development, biomineralization and innate immune response. This review advances our understanding on the roles of ncRNAs in molluscs and suggest the future direction to fully understand the epigenetic regulatory network of molluscs.

Characterization of Insulin-Like Growth Factor Binding Protein 7 (Igfbp7) and Its Potential Involvement in Shell Formation and Metamorphosis of Pacific Abalone, Haliotis discus hannai

Insulin-like growth factor binding proteins (IGFBPs) are secreted proteins that play an important role in IGF regulation of growth and development of vertebrate and invertebrates. In this study, the IGFBP7 gene was cloned and characterized from mantle tissues of H. discus hannai, and designated as Hdh IGFBP7. The full-length cDNA sequence transcribed from the Hdh IGFBP7 gene was 1519-bp long with an open reading frame of 720-bp corresponding to a putative polypeptide of 239 amino acids. The molecular mass of its mature protein was approximately 23.44 KDa with an estimated isoelectric point (pI) of 5.35, and it shared significant homology with IGFBP7 gene of H. madaka. Hdh IGFBP7 has a characteristic IGFBP N-terminal domain (22–89 aa), a kazal-type serine proteinase inhibitor domain (77–128), and an immunoglobulin-like C2 domain (144–223). Furthermore, twelve cysteine residues and a signature motif of IGFBPs (XCGCCXXC) were found in its N-terminal domain. Phylogenetic analysis revealed that Hdh IGFBP7 was aligned with IGFBP7 of H. madaka. Tissue distribution analysis showed that the mRNA of Hdh IGFBP7 was expressed in all examined tissues, with the highest expression level observed in the mantle and gill tissues. The expression level of Hdh IGFBP7 mRNA was relatively higher at the juvenile stage during its metamorphosis period. In situ hybridization showed that Hdh IGFBP7 transcript was expressed in epithelial cells of the dorsal mantle pallial and mucus cells of the branchial epithelium in gill. These results provide basic information for future studies on the role of IGFBP7 in IGF regulation of shell growth, development and metamorphosis of abalone.

The Scaly-foot Snail genome and implications for the origins of biomineralised armour

The Scaly-foot Snail, Chrysomallon squamiferum, presents a combination of biomineralised features, reminiscent of enigmatic early fossil taxa with complex shells and sclerites such as sachtids, but in a recently-diverged living species which even has iron-infused hard parts. Thus the Scaly-foot Snail is an ideal model to study the genomic mechanisms underlying the evolutionary diversification of biomineralised armour. Here, we present a high-quality whole-genome assembly and tissue-specific transcriptomic data, and show that scale and shell formation in the Scaly-foot Snail employ independent subsets of 25 highly-expressed transcription factors. Comparisons with other lophotrochozoan genomes imply that this biomineralisation toolkit is ancient, though expression patterns differ across major lineages. We suggest that the ability of lophotrochozoan lineages to generate a wide range of hard parts, exemplified by the remarkable morphological disparity in Mollusca, draws on a capacity for dynamic modification of the expression and positioning of toolkit elements across the genome.

Is adaptation limited by mutation? A timescale-dependent effect of genetic diversity on the adaptive substitution rate in animals

Whether adaptation is limited by the beneficial mutation supply is a long-standing question of evolutionary genetics, which is more generally related to the determination of the adaptive substitution rate and its relationship with species effective population size (Ne) and genetic diversity. Empirical evidence reported so far is equivocal, with some but not all studies supporting a higher adaptive substitution rate in large-Ne than in small-Ne species. We gathered coding sequence polymorphism data and estimated the adaptive amino-acid substitution rate ωa, in 50 species from ten distant groups of animals with markedly different population mutation rate θ. We reveal the existence of a complex, timescale dependent relationship between species adaptive substitution rate and genetic diversity. We find a positive relationship between ωa and θ among closely related species, indicating that adaptation is indeed limited by the mutation supply, but this was only true in relatively low-θ taxa. In contrast, we uncover no significant correlation between ωa and θ at a larger taxonomic scale, suggesting that the proportion of beneficial mutations scales negatively with species' long-term Ne.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Increased performance of DNA metabarcoding of macroinvertebrates by taxonomic sorting

DNA-based identification through the use of metabarcoding has been proposed as the next step in the monitoring of biological communities, such as those assessed under the Water Framework Directive (WFD). Advances have been made in the field of metabarcoding, but challenges remain when using complex samples. Uneven biomass distributions, preferential amplification and reference database deficiencies can all lead to discrepancies between morphological and DNA-based taxa lists. The effects of different taxonomic groups on these issues remain understudied. By metabarcoding WFD monitoring samples, we analyzed six different taxonomic groups of freshwater organisms, both separately and combined. Identifications based on metabarcoding data were compared directly to morphological assessments performed under the WFD. The diversity of taxa for both morphological and DNA-based assessments was similar, although large differences were observed in some samples. The overlap between the two taxon lists was 56.8% on average across all taxa, and was highest for Crustacea, Heteroptera, and Coleoptera, and lowest for Annelida and Mollusca. Taxonomic sorting in six basic groups before DNA extraction and amplification improved taxon recovery by 46.5%. The impact on ecological quality ratio (EQR) scoring was considerable when replacing morphology with DNA-based identifications, but there was a high correlation when only replacing a single taxonomic group with molecular data. Different taxonomic groups provide their own challenges and benefits. Some groups might benefit from a more consistent and robust method of identification. Others present difficulties in molecular processing, due to uneven biomass distributions, large genetic diversity or shortcomings of the reference database. Sorting samples into basic taxonomic groups that require little taxonomic knowledge greatly improves the recovery of taxa with metabarcoding. Current standards for EQR monitoring may not be easily replaced completely with molecular strategies, but the effectiveness of molecular methods opens up the way for a paradigm shift in biomonitoring.

Untargeted metabolomic analysis of the carotenoid-based orange coloration in Haliotis gigantea using GC-TOF-MS

Seafood coloration is typically considered an indicator of quality and nutritional value by consumers. One such seafood is the Xishi abalone (Haliotis gigantea), which displays muscle color polymorphism wherein a small subset of individuals display orange coloration of muscles due to carotenoid enrichment. However, the metabolic basis for carotenoid accumulation has not been thoroughly investigated in marine mollusks. Here, GC-TOF-MS-based untargeted metabolite profiling was used to identify key pathways and metabolites involved in differential carotenoid accumulation in abalones with variable carotenoid contents. Cholesterol was the most statistically significant metabolite that differentiated abalones with orange muscles against those with common white muscles. This observation is likely due to the competitive interactions between cholesterol and carotenoids during cellular absorption. In addition, the accumulation of carotenoids was also related to fatty acid contents. Overall, this study indicates that metabolomics can reflect physiological changes in organisms and provides a useful framework for exploring the mechanisms underlying carotenoid accumulation in abalone types.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

The retinoic acid receptor (RAR) in molluscs: Function, evolution and endocrine disruption insights

Retinoid acid receptor (RAR)-dependent signalling pathways are essential for the regulation and maintenance of essential biological functions and are recognized targets of disruptive anthropogenic compounds. Recent studies put forward the inability of mollusc RARs to bind and respond to the canonical vertebrate ligand, retinoic acid: a feature that seems to have been lost during evolution. Yet, these studies were carried out in a limited number of molluscs. Therefore, using an in vitro transactivation assay, the present work aimed to characterize phylogenetically relevant mollusc RARs, as monomers or as functional units with RXR, not only in the presence of vertebrate bone fine ligands but also known endocrine disruptors, described to modulate retinoid-dependent pathways. In general, none of the tested mollusc RARs were able to activate reporter gene transcription when exposed to retinoic acid isomers, suggesting that the ability to respond to retinoic acid was lost across molluscs. Similarly, the analysed mollusc RAR were unresponsive towards organochloride pesticides. In contrast, transcriptional repressions were observed with the RAR/RXR unit upon exposure to retinoids or RXR-specific ligands. Loss-of-function and gain-of-function mutations further corroborate the obtained results and suggest that the repressive behaviour, observed with mollusc and human RAR/RXR heterodimers, is possibly mediated by ligand biding to RXR.

Discovery of Novel Hemocyanin-Like Genes in Metazoans

Among animals, two major groups of oxygen-binding proteins are found: proteins that use iron to bind oxygen (hemoglobins and hemerythrins) and two non-homologous hemocyanins that use copper. Although arthropod and mollusc hemocyanins bind oxygen in the same manner, they are distinct in their molecular structures. In order to better understand the range of natural variation in hemocyanins, we searched for them in a diverse array of metazoan transcriptomes by using bioinformatics tools to examine hemocyanin evolutionary history and to consequently revive the discussion about whether all metazoan hemocyanins shared a common origin with frequent losses or whether they originated separately after the divergence of Lophotrochozoa and Ecdysozoa. We confirm that the distribution of hemocyanin-like genes is more widespread than previously reported, including five putative novel mollusc hemocyanin genes in two annelid species from Chaetopteridae. For arthropod hemocyanins, 16 putative novel genes were retained, and the presence of arthropod hemocyanins in 11 annelid species represents a novel observation. Interestingly, Annelida is the lineage that presents the greatest repertoire of oxygen transport proteins reported to date, possessing all the main superfamily proteins, which could be explained partially by the immense variability of lifestyles and habitats. Work presented here contradicts the canonical view that hemocyanins are restricted to molluscs and arthropods, suggesting that the occurrence of copper-based blood pigments in metazoans has been underestimated. Our results also support the idea of the presence of oxygen carrier hemocyanins being widespread across metazoans with an evolutionary history characterized by frequent losses.

Staggered Hox expression is more widespread among molluscs than previously appreciated

Hox genes are expressed along the anterior-posterior body axis in a colinear fashion in the majority of bilaterians. Contrary to polyplacophorans, a group of aculiferan molluscs with conserved ancestral molluscan features, gastropods and cephalopods deviate from this pattern by expressing Hox genes in distinct morphological structures and not in a staggered fashion. Among conchiferans, scaphopods exhibit many similarities with gastropods, cephalopods and bivalves, however, the molecular developmental underpinnings of these similar traits remain unknown. We investigated Hox gene expression in developmental stages of the scaphopod Antalis entalis to elucidate whether these genes are involved in patterning morphological traits shared by their kin conchiferans. Scaphopod Hox genes are predominantly expressed in the foot and mantle but also in the central nervous system. Surprisingly, the scaphopod mid-stage trochophore exhibits a near-to staggered expression of all nine Hox genes identified. Temporal colinearity was not found and early-stage and late-stage trochophores, as well as postmetamorphic individuals, do not show any apparent traces of staggered expression. In these stages, Hox genes are expressed in distinct morphological structures such as the cerebral and pedal ganglia and in the shell field of early-stage trochophores. Interestingly, a re-evaluation of previously published data on early-stage cephalopod embryos and of the gastropod pre-torsional veliger shows that these developmental stages exhibit traces of staggered Hox expression. Considering our results and all gene expression and genomic data available for molluscs as well as other bilaterians, we suggest a last common molluscan ancestor with colinear Hox expression in predominantly ectodermal tissues along the anterior-posterior axis. Subsequently, certain Hox genes have been co-opted into the patterning process of distinct structures (apical organ or prototroch) in conchiferans.

Advances in DNA Barcoding of Toxic Marine Organisms

There are more than 200,000 marine species worldwide. These include many important economic species, such as large yellow croaker, ribbonfish, tuna, and salmon, but also many potentially toxic species, such as blue-green algae, diatoms, cnidarians, ctenophores, Nassarius spp., and pufferfish. However, some edible and toxic species may look similar, and the correct identification of marine species is thus a major issue. The failure of traditional classification methods in certain species has promoted the use of DNA barcoding, which uses short, standard DNA fragments to assist with species identification. In this review, we summarize recent advances in DNA barcoding of toxic marine species such as jellyfish and pufferfish, using genes including cytochrome oxidase I gene (COI), cytochrome b gene (cytb), 16S rDNA, internal transcribed spacer (ITS), and Ribulose-1,5-bisphosphate carboxylase oxygenase gene (rbcL). We also discuss the application of this technique for improving the identification of marine species. The use of DNA barcoding can benefit the studies of biological diversity, biogeography, food safety, and the detection of both invasive and new species. However, the technique has limitations, particularly for the analysis of complex objects and the selection of standard DNA barcodes. The development of high-throughput methods may offer solutions to some of these issues.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Molecular modularity and asymmetry of the molluscan mantle revealed by a gene expression atlas

Background

Conchiferan molluscs construct a biocalcified shell that likely supported much of their evolutionary success. However, beyond broad proteomic and transcriptomic surveys of molluscan shells and the shell-forming mantle tissue, little is known of the spatial and ontogenetic regulation of shell fabrication. In addition, most efforts have been focused on species that deposit nacre, which is at odds with the majority of conchiferan species that fabricate shells using a crossed-lamellar microstructure, sensu lato.

Results

By combining proteomic and transcriptomic sequencing with in situ hybridization we have identified a suite of gene products associated with the production of the crossed-lamellar shell in Lymnaea stagnalis. With this spatial expression data we are able to generate novel hypotheses of how the adult mantle tissue coordinates the deposition of the calcified shell. These hypotheses include functional roles for unusual and otherwise difficult-to-study proteins such as those containing repetitive low-complexity domains. The spatial expression readouts of shell-forming genes also reveal cryptic patterns of asymmetry and modularity in the shell-forming cells of larvae and adult mantle tissue.

Conclusions

This molecular modularity of the shell-forming mantle tissue hints at intimate associations between structure, function, and evolvability and may provide an elegant explanation for the evolutionary success of the second largest phylum among the Metazoa.

DNA barcoding reveals seasonal shifts in diet and consumption of deep-sea fishes in wedge-tailed shearwaters

The foraging ecology of pelagic seabirds is difficult to characterize because of their large foraging areas. In the face of this difficulty, DNA metabarcoding may be a useful approach to analyze diet compositions and foraging behaviors. Using this approach, we investigated the diet composition and its seasonal variation of a common seabird species on the Ogasawara Islands, Japan: the wedge-tailed shearwater Ardenna pacifica. We collected fecal samples during the prebreeding (N = 73) and rearing (N = 96) periods. The diet composition of wedge-tailed shearwater was analyzed by Ion Torrent sequencing using two universal polymerase chain reaction primers for the 12S and 16S mitochondrial DNA regions that targeted vertebrates and mollusks, respectively. The results of a BLAST search of obtained sequences detected 31 and 1 vertebrate and mollusk taxa, respectively. The results of the diet composition analysis showed that wedge-tailed shearwaters frequently consumed deep-sea fishes throughout the sampling season, indicating the importance of these fishes as a stable food resource. However, there was a marked seasonal shift in diet, which may reflect seasonal changes in food resource availability and wedge-tailed shearwater foraging behavior. The collected data regarding the shearwater diet may be useful for in situ conservation efforts. Future research that combines DNA metabarcoding with other tools, such as data logging, may provide further insight into the foraging ecology of pelagic seabirds.