The shell-forming proteome ofLottia giganteareveals both deep conservations and lineage-specific novelties

The adult shell matrix protein repertoire of the marine snail Crepidula is dominated by conserved genes that are also expressed in larvae

Mollusca is a morphologically diverse phylum, exhibiting an immense variety of calcium carbonate structures. Proteomic studies of adult shells often report high levels of rapidly-evolving, 'novel' shell matrix proteins (SMPs), which are hypothesized to drive shell diversification. However, relatively little is known about the phylogenetic distribution of SMPs, or about the function of individual SMPs in shell construction. To understand how SMPs contribute to shell diversification a thorough characterization of SMPs is required. Here, we build tools and a foundational understanding of SMPs in the marine gastropod species Crepidula fornicata and Crepidula atrasolea because they are genetically-enabled mollusc model organisms. First, we established a staging system of shell development in C. atrasolea for the first time. Next, we leveraged previous findings in C. fornicata combined with phylogenomic analyses of 95 metazoan species to determine the evolutionary lineage of its adult SMP repertoire. We found that 55% of C. fornicata's SMPs belong to molluscan orthogroups, with 27% restricted to Gastropoda, and only 5% restricted at the species level. The low percentage of species-restricted SMPs underscores the importance of broad-taxon sampling and orthology inference approaches when determining homology of SMPs. From our transcriptome analysis, we found that the majority of C. fornicata SMPs that were found conserved in C. atrasolea were expressed in both larval and adult stages. We then selected a subset of SMPs of varying evolutionary ages for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval shell development in C. atrasolea. Out of the 18 SMPs analyzed, 12 were detected in the larval shell field. These results suggest overlapping larval vs. adult SMP repertoires. Using multiplexed HCR, we observed five SMP expression patterns and three distinct cell populations within the shell field. These patterns support the idea that modular expression of SMPs could facilitate divergence of shell morphological characteristics. Collectively, these data establish an evolutionary and developmental framework in Crepidula that enables future comparisons of molluscan biomineralization to reveal mechanisms of shell diversification.

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

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

Proteomic diversification of spermatostyles among six species of whirligig beetles

Seminal fluid protein composition is complex and commonly assumed to be rapidly divergent due to functional interactions with both sperm and the female reproductive tract (FRT), both of which evolve rapidly. In addition to sperm, seminal fluid may contain structures, such as mating plugs and spermatophores. Here, we investigate the evolutionary diversification of a lesser-known ejaculate structure: the spermatostyle, which has independently arisen in several families of beetles and true bugs. We characterized the spermatostyle proteome, in addition to spermatostyle and FRT morphology, in six species of whirligig beetles (family Gyrinidae). Spermatostyles were enriched for proteolytic enzymes, and assays confirmed they possess proteolytic activity. Sperm-leucylaminopeptidases (S-LAPs) were particularly abundant, and their localization to spermatostyles was confirmed by immunohistochemistry. Although there was evidence for functional conservation of spermatostyle proteomes across species, phylogenetic regressions suggest evolutionary covariation between protein composition and the morphology of both spermatostyles and FRTs. We postulate that S-LAPs (and other proteases) have evolved a novel structural role in spermatostyles and discuss spermatostyles as adaptations for delivering male-derived materials to females.

Antibacterial Action of Protein Fraction Isolated from Rapana venosa Hemolymph against Escherichia coli NBIMCC 8785

Natural products and especially those from marine organisms are being intensively explored as an alternative to synthetic antibiotics. However, the exact mechanisms of their action are not yet well understood. The molecular masses of components in the hemolymph fraction with MW 50-100 kDa from Rapana venosa were determined using ImageQuant™ TL v8.2.0 software based on electrophoretic analysis. Mainly, three types of compounds with antibacterial potential were identified, namely proteins with MW at 50.230 kDa, 62.100 kDa and 93.088 kDa that were homologous to peroxidase-like protein, aplicyanin A and L-amino acid oxidase and functional units with MW 50 kDa from R. venous hemocyanin. Data for their antibacterial effect on Escherichia coli NBIMCC 8785 were obtained by CTC/DAPI-based fluorescent analysis (analysis based on the use of a functional fluorescence probe). The fluorescent analyses demonstrated that a 50% concentration of the fraction with MW 50-100 kDa was able to eliminate 99% of the live bacteria. The antimicrobial effect was detectable even at a 1% concentration of the active compounds. The bacteria in this case had reduced metabolic activity and a 24% decreased size. The fraction had superior action compared with another mollusc product-snail slime-which killed 60% of the E. coli NBIMCC 8785 cells at a 50% concentration and had no effect at a 1% concentration. The obtained results demonstrate the high potential of the fraction with MW 50-100 kDa from R. venosa to eliminate and suppress the development of Escherichia coli NBIMCC 8785 bacteria and could be applied as an appropriate component of therapeutics with the potential to replace antibiotics to avoid the development of antibiotic resistance.

EF-Hand-Binding Secreted Protein Hdh-SMP5 Regulates Shell Biomineralization and Responses to Stress in Pacific Abalone, Haliotis discus hannai

The development of a shell is a complex calcium metabolic process involving shell matrix proteins (SMPs). In this study, we describe the isolation, characterization, and expression of SMP5 and investigate its potential regulatory role in the shell biomineralization of Pacific abalone Haliotis discus hannai. The full-length Hdh-SMP5 cDNA contains 685 bp and encodes a polypeptide of 134 amino acids. Structurally, the Hdh-SMP5 protein belongs to the EF-hand-binding superfamily, which possesses three EF-hand Ca2+-binding regions and is rich in aspartic acid. The distinct clustering patterns in the phylogenetic tree indicate that the amino acid composition and structure of this protein may vary among different SMPs. During early development, significantly higher expression was observed in the trochophore and veliger stages. Hdh-SMP5 was also upregulated during shell biomineralization in Pacific abalone. Long periods of starvation cause Hdh-SMP5 expression to decrease. Furthermore, Hdh-SMP5 expression was observed to be significantly higher under thermal stress at temperatures of 15, 30, and 25 °C for durations of 6 h, 12 h, and 48 h, respectively. Our study is the first to characterize Hdh-SMP5 comprehensively and analyze its expression to elucidate its dynamic roles in ontogenetic development, shell biomineralization, and the response to starvation and thermal stress.

Multi-omic insights into the formation and evolution of a novel shell microstructure in oysters

BACKGROUND

Molluscan shell, composed of a diverse range of architectures and microstructures, is a classic model system to study the relationships between molecular evolution and biomineralized structure formation. The shells of oysters differ from those of other molluscs by possessing a novel microstructure, chalky calcite, which facilitates adaptation to the sessile lifestyle. However, the genetic basis and evolutionary origin of this adaptive innovation remain largely unexplored.

RESULTS

We report the first whole-genome assembly and shell proteomes of the Iwagaki oyster Crassostrea nippona. Multi-omic integrative analyses revealed that independently expanded and co-opted tyrosinase, peroxidase, TIMP genes may contribute to the chalky layer formation in oysters. Comparisons with other molluscan shell proteomes imply that von Willebrand factor type A and chitin-binding domains are basic members of molluscan biomineralization toolkit. Genome-wide identification and analyses of these two domains in 19 metazoans enabled us to propose that the well-known Pif may share a common origin in the last common ancestor of Bilateria. Furthermore, Pif and LamG3 genes acquire new genetic function for shell mineralization in bivalves and the chalky calcite formation in oysters likely through a combination of gene duplication and domain reorganization.

CONCLUSIONS

The spatial expression of SMP genes in the mantle and molecular evolution of Pif are potentially involved in regulation of the chalky calcite deposition, thereby shaping the high plasticity of the oyster shell to adapt to a sessile lifestyle. This study further highlights neo-functionalization as a crucial mechanism for the diversification of shell mineralization and microstructures in molluscs, which may be applied more widely for studies on the evolution of metazoan biomineralization.

Spatial-temporal expression analysis of lineage-restricted shell matrix proteins reveals shell field regionalization and distinct cell populations in the slipper snail Crepidula atrasolea

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea , to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata -restricted SMPs contradicts the generally-held notion that an animal’s biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea . We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

Antitumor Activity of Bioactive Compounds from Rapana venosa against Human Breast Cell Lines

This study is the first report describing the promising antitumor activity of biologically active compounds isolated from the hemolymph of marine snail Rapana venosa-a fraction with Mw between 50 and 100 kDa and two structural subunits (RvH1 and RvH2), tested on a panel of human breast cell lines-six lines of different molecular subtypes of breast cancer MDA-MB-231, MDA-MB-468, BT-474, BT-549, SK-BR-3, and MCF-7 and the non-cancerous MCF-10A. The fraction with Mw 50-100 kDa (HRv 50-100) showed good antitumor activity manifested by a significant decrease in cell viability, altered morphology, autophagy, and p53 activation in treated cancer cells. An apparent synergistic effect was observed for the combination of HRv 50-100 with cis-platin for all tested cell lines. The combination of HRv 50-100 with cisplatin and/or tamoxifen is three times more effective compared to treatment with classical chemotherapeutics alone. The main proteins in the active fraction, with Mw at ~50 kDa, ~65 kDa, ~100 kDa, were identified by MALDI-MS, MS/MS analyses, and bioinformatics. Homology was established with known proteins with antitumor potential detected in different mollusc species: peroxidase-like protein, glycoproteins Aplysianin A, L-amino acid oxidase (LAAO), and the functional unit with Mw 50 kDa of RvH. Our study reveals new perspectives for application of HRv 50-100 as an antitumor agent used alone or as a booster in combination with different chemotherapies.

Characterization of larval shell formation and CgPOU2F1, CgSox5, and CgPax6 gene expression during shell morphogenesis in Crassostrea gigas

Shell formation is a dynamic process involving organic matrix secretion and calcification. In this study, we characterized shell morphogenesis during larval development in Crassostrea gigas. Using scanning electron microscopy (SEM) and fluorescence staining, we demonstrated that shell field, the first morphologically distinguishable shell-forming tissue, became visible soon after enlargement of the blastopore at the anterior end of the trochophore. Shell organic matrix namely protein polysaccharides and calcified structure appeared as a slit at the dorsal side of the embryo. The early shell field began to extend along the dorsal side of the trochophore larvae, and became a saddle shaped shell field that gave rise to the prodissoconch I embryonic shell in the early D-shaped larvae. Subsequently, prodissoconch II shell was formed in the late D-shaped larvae with a characteristic appearance of growth lines. To identify gene expression markers for studying shell formation, we isolated three potential larval shell formation genes CgPOU2F1, CgSox5, and CgPax6 and analyzed their expression during shell morphogenesis. The three potential shell formation genes possessed a similar pattern of expression. Their expression was detected in the shell gland and shell field regions in early D-shaped larvae, hereafter, their expression was detected at the larval mantle edge in the calcified shell stages. Together, these studies provide knowledge of shell morphogenesis in pacific oyster and molecular markers for studying the molecular regulation of biomineralization and shell formation.

Evolution of nacre- and prisms-related shell matrix proteins in the pen shell, Atrina pectinata

The molluscan shell is a good model for understanding the mechanisms underlying biomineralization. It is composed of calcium carbonate crystals and many types of organic molecules, such as the matrix proteins, polysaccharides, and lipids. The pen shell Atrina pectinata (Pterioida, Pinnidae) has two shell microstructures: an outer prismatic layer and an inner nacreous layer. Similar microstructures are well known in pearl oysters (Pteriidae), such as Pinctada fucata, and many kinds of shell matrix proteins (SMPs) have been identified from their shells. However, the members of SMPs that consist of the nacreous and prismatic layers of Pinnidae bivalves remain unclear. In this study, we identified 114 SMPs in the nacreous and prismatic layers of A. pectinata, of which only seven were found in both microstructures. 54 of them were found to bind calcium carbonate. Comparative analysis of nine molluscan shell proteomes showed that 69 of 114 SMPs of A. pectinata were found to have sequential similarity with at least one or more SMPs of other molluscan species. For instance, nacrein, tyrosinase, Pif/BMSP-like, chitinase (CN), chitin-binding proteins, CD109, and Kunitz-type serine proteinase inhibitors are widely shared among bivalves and gastropods. Our results provide new insights for understanding the complex evolution of SMPs related to nacreous and prismatic layer formation in the pteriomorph bivalves.

Proteomic and Transcriptomic Analyses in the Slipper Snail Crepidula fornicata Uncover Shell Matrix Genes Expressed During Adult and Larval Biomineralization

The gastropod shell is a composite composed of minerals and shell matrix proteins (SMPs). SMPs have been identified by proteomics in many molluscs, but few have been studied in detail. Open questions include (1) what gene regulatory networks regulate SMP expression, (2) what roles individual SMPs play in biomineralization, and (3) how the complement of SMPs changes over development. These questions are best addressed in a species in which gene perturbation studies are available; one such species is the slipper snail, Crepidula fornicata. Here, SEM and pXRD analysis demonstrated that the adult shell of C. fornicata exhibits crossed lamellar microstructure and is composed of aragonite. Using high-throughput proteomics we identified 185 SMPs occluded within the adult shell. Over half of the proteins in the shell proteome have known biomineralization domains, while at least 10% have no homologs in public databases. Differential gene expression analysis identified 20 SMP genes that are up-regulated in the shell-producing mantle tissue. Over half of these 20 SMPs are expressed during development with two, CfSMP1 and CfSMP2, expressed exclusively in the shell gland. Together, the description of the shell microstructure and a list of SMPs now sets the stage for studying the consequences of SMP gene knockdowns in molluscs.

An acidic protein, Hf15, from Haliotis fulgens involved in biomineralization

Biomineralization leads to the hardening of mineralized materials, such as the shell of Mollusk, to fulfill a wide range of functions, such as (but not limited to) skeletal support, protection of the soft tissues, navigation, etc. The study of the proteins responsible for this process, shell matrix proteins (SMPs), allows addressing questions related to structure-function relationship and to the mechanism of mineral formation, which is limited in gastropod species. In this study, a low molecular weight protein was isolated from the insoluble fraction after decalcification with acetic acid of the shell of Haliotis fulgens and, named Hf15. The unglycosylated protein has a theoretical molecular weight of 15 kDa, it possesses calcium and chiting binding properties. Hf15 can precipitate calcium carbonate in vitro in presence of different salts. Analysis by LC-MS of the five peptide sequences of Hf15 generated by trypsinization revealed that two peptides displayed homology to an uncharacterized protein 3-like from Haliotis rufescens, Haliotis asinia and H. sorenseni. The results obtained indicated that Hf15 is a novel SMP involved in shell mineralization in Haliotis fulgens.

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.

The proteomics of the freshwater pearl powder: Insights from biomineralization to biomedical application

The freshwater pearl is one kind of valuable organic jewelry and traditional Chinese medicine (TCM). However, the molecular basis of matrix protein in pearl biomineralization and biomedical applications are largely unknown to date. In this study, the matrix proteins of water-soluble matrix, acid-soluble matrix and acid-insoluble matrix from the freshwater seedless pearl powder were detected using liquid chromatography-tandem mass spectrometry (LC-MS/MS) respectively, and identified against the transcriptomic database of the pearl sac. The results showed that a total of 190 proteins were identified in pearl proteomics, which was divided into eight categories by their potential biomineralization functions. The composition of pearl matrix proteins and the high frequency conserved domains like carbonic anhydrase, von Willebrand factor type A, tyrosinase and chitin binding 2 in protein sequences, implying that the "chitin-silk fibroin gel proteins-acidic macromolecules" model was suitable for description the pearl biomineralization process. Meanwhile, ninety-one of pearl matrix proteins could be classified into seven categories by their potential medical functions including wound healing, osteogenic property, antioxidant activity, neuro-regulation effects, skin lightening effect, anti-inflammatory and anti-apoptotic effects and other immunomodulatory property. In general, these results provided valuable new insights into not only the diversity of pearl matrix protein for mollusc biomineralization, but the molecular basis of pearl matrix proteins responsible for their diverse biological properties in TCM application. SIGNIFICANCE: The significance of this study included the following points.

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.

Transcriptome analysis of the bivalve Placuna placenta mantle reveals potential biomineralization-related genes

The shells of window pane oyster Placuna placenta are very thin and exhibit excellent optical transparency and mechanical robustness. However, little is known about the biomineralization-related proteins of the shells of P. placenta. In this work, we report the comprehensive transcriptome of the mantle tissue of P. placenta for the first time. The unigenes of the mantle tissue of P. placenta were annotated by using the public databases such as nr, GO, KOG, KEGG, and Pfam. 24,343 unigenes were annotated according to Pfam database, accounting for 21.48% of the total unigenes. We find that half of the annotated unigenes of the mantle tissue of P. placenta are consistent to the annotated unigenes from pacific oyster Crassostrea gigas according to nr database. The unigene sequence analysis from the mantle tissue of P. placenta indicates that 465,392 potential single nucleotide polymorphisms (SNPs) and 62,103 potential indel markers were identified from 60,371 unigenes. 178 unigenes of the mantle tissue of P. placenta are found to be homologous to those reported proteins related to the biomineralization process of molluscan shells, while 18 of them are highly expressed unigenes in the mantle tissue. It is proposed that four unigenes with the highest expression levels in the mantle tissue are very often related to the biomineralization process, while another three unigenes are potentially related to the biomineralization process according to the Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) analysis. In summary, the transcriptome analysis of the mantle tissue of P. Placenta shows the potential biomineralization-related proteins and this work may shed light for the shell formation mechanism of bivalves.

Transcriptomic response of the intertidal limpet Patella vulgata to temperature extremes

Global warming is challenging wild species in land and water. In the intertidal zone, species are already living at their thermal limits, being vulnerable even to small increases in maximum habitat temperatures. Knowledge of the mechanisms by which many intertidal zone species cope with elevated temperatures is limited. We analysed the molecular thermal stress response of the limpet Patella vulgata under slight and frequent (one-day), and extreme and rare (three-day) warming events. Using RNA-seq to assess differential gene expression among treatments, differing molecular responses were obtained in the two treatments, with more changes in gene expression after the three-day event; with one-third of the differentially expressed transcripts being down-regulated. However, across treatments we observed shifts in gene expression for common aspects of the heat stress response including intra-cellular communication, protein chaperoning, proteolysis and cell cycle arrest. Of the 71,675 transcripts obtained, only 259 were differentially expressed after both heating events. From these, 218 defined the core group (i.e. genes induced by thermal stress with similar expression patterns irrespective of the magnitude of the warming event). The core group was composed of already well-studied genes in heat stress responses in intertidal organisms (e.g. heat shock proteins), but also genes from less explored metabolic pathways, e.g. the ubiquitin system, which were also fundamental regardless of the magnitude of the imposed warming. Moreover, we have also identified 41 signaling genes (i.e. a set of genes responding to both events and with expression patterns specific to the intensity of thermal stress), principally including genes involved in the maintenance of extracellular structure that have previously not been identified as part of the response to thermal stress in intertidal zone organisms. These signaling genes will be useful heat stress molecular biomarkers for monitoring heat stress in natural populations.

The potential for using shell proteins in gastropod systematics, assessed in patellogastropod limpets

This investigation of the application of shell protein information to gastropod systematics initially utilized available Lottia gigantea sequences and a transcriptome of Patelloida mimula developed here. Levels of differentiation between predicted sequences of reciprocal best-hit potential homologues in P. mimula and L. gigantea suggested that they could be useful within families, and possibly in higher taxa using some shell-associated proteins, particularly the peroxidases. Subsequently, proteomic analyses of the acid-soluble fraction of extractions from 17 shells and five tissue samples were conducted by combined liquid chromatography/mass spectrometry with nano-electrospray ionization. All proteins with abundance more than 1.2% in the L. gigantea shell proteome were identified with 100% confidence in most extractions by SearchGui/PeptideShaker analyses. In total, 259 of 379 peptides predicted from in silico digestion of L. gigantea shell proteins were represented by validated peptide spectrum matches in one or more specimens. Systematics applications were investigated by analysing metrics such as protein coverage by peptides and phylogenetic analyses of peptide presence/absence. The investigation suggested that diagnostic profiles based on fixed presence/absence differences can be used to separate species pairs. However, further development of analytical techniques and accumulation of reference databases is required for realising fully the systematics potential of the shell proteome.

Proteomics of Shell Matrix Proteins from the Cuttlefish Bone Reveals Unique Evolution for Cephalopod Biomineralization

In contrast to the external shells in bivalves and gastropods, most cephalopods are missing this external protection. The cuttlefish, belonging to class cephalopod, has an internal biomineralized structure made of mainly calcium carbonate for controlling buoyancy. However, the macromolecules, especially proteins that control cuttlebone mineral formation, are not sufficiently understood, limiting our understanding of the evolution of this internal shell. In this study, we extracted proteins from the cuttlebone of pharaoh cuttlefish Sepia pharaonis and performed liquid chromatography-tandem mass spectrometry to identify the shell matrix proteins (SMPs). In total, 41 SMPs were identified. Among them, hemocyanin, an oxygen-carrying protein, was the most abundant SMP. By comparison with SMPs of other marine biominerals, hemocyanin, apolipophorin, soul domain proteins, transferrin, FL-rich, and enolase were found to be unique to the cuttlebone. In contrast, typical SMPs of external shells such as carbonic anhydrase complement control protein, fibronectin type III, and G/A-rich proteins were lacking from the cuttlebone. Furthermore, the cluster analysis of biomineral SMPs suggests that the SMP repertoire of the cuttlebone does not resemble that of other species with external shells. Taken together, this study implies a potential relationship of the cuttlefish internal shell with other internal biominerals, which highlights a unique shell evolutionary pathway in invertebrates.

Evolution and biomineralization of pteropod shells

Shelled pteropods, known as sea butterflies, are a group of small gastropods that spend their entire lives swimming and drifting in the open ocean. They build thin shells of aragonite, a metastable polymorph of calcium carbonate. Pteropod shells have been shown to experience dissolution and reduced thickness with a decrease in pH and therefore represent valuable bioindicators to monitor the impacts of ocean acidification. Over the past decades, several studies have highlighted the striking diversity of shell microstructures in pteropods, with exceptional mechanical properties, but their evolution and future in acidified waters remains uncertain. Here, we revisit the body-of-work on pteropod biomineralization, focusing on shell microstructures and their evolution. The evolutionary history of pteropods was recently resolved, and thus it is timely to examine their shell microstructures in such context. We analyse new images of shells from fossils and recent species providing a comprehensive overview of their structural diversity. Pteropod shells are made of the crossed lamellar and prismatic microstructures common in molluscs, but also of curved nanofibers which are proposed to form a helical three-dimensional structure. Our analyses suggest that the curved fibres emerged before the split between coiled and uncoiled pteropods and that they form incomplete to multiple helical turns. The curved fibres are seen as an important trait in the adaptation to a planktonic lifestyle, giving maximum strength and flexibility to the pteropod thin and lightweight shells. Finally, we also elucidate on the candidate biomineralization genes underpinning the shell diversity in these important indicators of ocean health.

The temptin gene of the clade Lophotrochozoa is involved in formation of the prismatic layer during biomineralization in molluscs

The biomineralization mechanism of mollusc shell has been studied for a long time, but there is a lack of understanding about the relationship between the shell formation in vitro and the signaling system in vivo. In this study, we cloned a novel shell matrix protein gene (hc-temptin), which only be characterized as a water-borne protein pheromone of molluscs in previous studies, from the freshwater mussel Hyriopsis cumingii. By bioinformatics analysis we found that temptin was a gene unique to the clade Lophotrochozoa, and it exists in all mollusc taxa except Cephalopoda. The current data supported the premise that temptin was generated in the early emergence of molluscs and that it maintained a high mutation rate to evolve relative independently. The specificity of hc-temptin expression in the mantle tissue suggests its potential to participate in biomineralization. Its sequence contained typical Ca²⁺ binding sites. Our experiments involving the pearl formation process, damaged shell repair process, and RNAi experiment showed that hc-temptin was a shell matrix protein that plays an important role in formation of the prismatic layer. The results of this study provided new insights about the origin of the temptin gene and its role in molluscs.

The genome of Nautilus pompilius illuminates eye evolution and biomineralization

Nautilus is the sole surviving externally shelled cephalopod from the Palaeozoic. It is unique within cephalopod genealogy and critical to understanding the evolutionary novelties of cephalopods. Here, we present a complete Nautilus pompilius genome as a fundamental genomic reference on cephalopod innovations, such as the pinhole eye and biomineralization. Nautilus shows a compact, minimalist genome with few encoding genes and slow evolutionary rates in both non-coding and coding regions among known cephalopods. Importantly, multiple genomic innovations including gene losses, independent contraction and expansion of specific gene families and their associated regulatory networks likely moulded the evolution of the nautilus pinhole eye. The conserved molluscan biomineralization toolkit and lineage-specific repetitive low-complexity domains are essential to the construction of the nautilus shell. The nautilus genome constitutes a valuable resource for reconstructing the evolutionary scenarios and genomic innovations that shape the extant cephalopods.

Review: Post-translational modifications of marine shell matrix proteins

Shell matrix proteins (SMPs) are key components for the Mollusk shell biomineralization. SMPs function has been hypothesized in several proteins by bioinformatics analysis, and through in vitro crystallization assays. However, studies of the post-translational modifications (PTMs) of SMPs, which contribute to their structure and the function, are limited. This review provides the current status of the SMPs with the most common PTMs described (glycosylation, phosphorylation, and disulfide bond formation) and their role in shell biomineralization. Also, recent studies based on recombinant production of SMPs are discussed. Finally, recommendations for the study of SMPs and their PTMs are provided. The review showed that PTMs are widely distributed in SMPs, and their presence on SMPs may contribute to the modulation of their activity in some SMPs, contributing to the crystal growth formation and differentiation through different mechanisms, however, in a few cases the lack of the PTMs do not alter their inherent function.

Comparative Proteomics of Octocoral and Scleractinian Skeletomes and the Evolution of Coral Calcification

Corals are the ecosystem engineers of coral reefs, one of the most biodiverse marine ecosystems. The ability of corals to form reefs depends on the precipitation of calcium carbonate (CaCO₃) under biological control. However, several mechanisms underlying coral biomineralization remain elusive, for example, whether corals employ different molecular machineries to deposit different CaCO₃ polymorphs (i.e., aragonite or calcite). Here, we used tandem mass spectrometry (MS/MS) to compare the proteins occluded in the skeleton of three octocoral and one scleractinian species: Tubipora musica and Sinularia cf. cruciata (calcite sclerites), the blue coral Heliopora coerulea (aragonitic skeleton), and the scleractinian aragonitic Montipora digitata. Reciprocal Blast analysis revealed extremely low overlap between aragonitic and calcitic species, while a core set of proteins is shared between octocorals producing calcite sclerites. However, the carbonic anhydrase CruCA4 is present in the skeletons of both polymorphs. Phylogenetic analysis highlighted several possible instances of protein co-option in octocorals. These include acidic proteins and scleritin, which appear to have been secondarily recruited for calcification and likely derive from proteins playing different functions. Similarities between octocorals and scleractinians included presence of a galaxin-related protein, carbonic anhydrases, and one hephaestin-like protein. Although the first two appear to have been independently recruited, the third appear to share a common origin. This work represents the first attempt to identify and compare proteins associated with coral skeleton polymorph diversity, providing several new research targets and enabling both future functional and evolutionary studies aimed at elucidating the origin and evolution of coral biomineralization.

Quantitative measures and 3D shell models reveal interactions between bands and their position on growing snail shells

The nature of shell growth in gastropods is useful because it preserves the ontogeny of shape, colour, and banding patterns, making them an ideal system for understanding how inherited variation develops, is established and maintained within a population. However, qualitative scoring of inherited shell characters means there is a lack of knowledge regarding the mechanisms that control fine variation. Here, we combine empirical measures of quantitative variation and 3D modeling of shells to understand how bands are placed and interact. By comparing five-banded Cepaea individuals to shells lacking individual bands, we show that individual band absence has minor but significant impacts upon the position of remaining bands, implying that the locus controlling band presence/absence mainly acts after position is established. Then, we show that the shell grows at a similar rate, except for the region below the lowermost band. This demonstrates that wider bands of Cepaea are not an artifact of greater shell growth on the lower shell; they begin wider and grow at the same rate as other bands. Finally, we show that 3D models of shell shape and banding pattern, inferred from 2D photos using ShellShaper software, are congruent with empirical measures. This work therefore establishes a method that may be used for comparative studies of quantitative banding variation in snail shells, extraction of growth parameters, and morphometrics. In the future, studies that link the banding phenotype to the network of shell matrix proteins involved in biomineralization and patterning may ultimately aid in understanding the diversity of shell forms found in molluscs.

Genotyping-by-Sequencing of the regional Pacific abalone (Haliotis discus) genomes reveals population structures and patterns of gene flow

Continuous monitoring of the present genetic status is essential to preserve the genetic resource of wild populations. In this study, we sequenced regional Pacific abalone Haliotis discus samples from three different locations around the Korean peninsula to assess population structure, utilizing Genotyping-by-Sequencing (GBS) method. Using PstI enzyme for genome reduction, we demonstrated the resultant library represented the whole genome region with even spacing, and as a result 16,603 single nucleotide variants (SNVs) were produced. Genetic diversity and population structure were investigated using several methods, and a strong genetic heterogeneity was observed in the Korean abalone populations. Additionally, by comparison of the variant sets among population groups, we were able to discover 26 Korean abalone population-specific SNVs, potentially associated with phenotype differences. This is the first study demonstrating the feasibility of GBS for population genetic study on H. discus. Our results will provide valuable data for the genetic conservation and management of wild abalone populations in Korea and help future GBS studies on the marine mollusks.

Biomineral proteomics: A tool for multiple disciplinary studies

The hard tissues of animals, such as skeletons and teeth, are constructed by a biologically controlled process called biomineralization. In invertebrate animals, biominerals are considered important for their evolutionary success. These biominerals are hieratical biocomposites with excellent mechanical properties, and their formation has intrigued researchers for decades. Although proteins account for ~5 wt% of biominerals, they are critical players in biomineralization. With the development of high-throughput analysis methods, such as proteomics, biomineral protein data are rapidly accumulating, thus necessitating a refined model for biomineralization. This review focuses on biomineral proteomics in invertebrate animals to highlight the diversity of biomineral proteins (generally 40-80 proteins), and the results indicate that biomineralization includes thermodynamic crystal growth as well as intense extracellular matrix activity and/or vesicle transport. Biominerals have multiple functions linked to biological immunity and antipathogen activity. A comparison of proteomes across species and biomineral types showed that von Willebrand factor type A and epidermal growth factor, which frequently couple with other extracellular domains, are the most common domains. Combined with species-specific repetitive low complexity domains, shell matrix proteins can be employed to predict biomineral types. Furthermore, this review discusses the applications of biomineral proteomics in diverse fields, such as tissue regeneration, developmental biology, archeology, environmental science, and material science.

Teaching an Old Shell New Tricks: Extracting DNA from Current, Historical, and Ancient Mollusk Shells

The use of unconventional DNA sources has increased because the acquisition of traditional samples can be invasive, destructive, or impossible. Mollusks are one group for which novel genetic sources are crucial, but methodology remains relatively undeveloped. Many species are important ecologically and in aquaculture production. However, mollusks have the highest number of extinctions of any taxonomic group. Traditionally, mollusk shell material was used for morphological research and only recently has been used in DNA studies. In the present article, we review the studies in which shell DNA was extracted and found that effective procedures consider taxon-specific biological characteristics, environmental conditions, laboratory methods, and the study objectives. Importantly, these factors cannot be considered in isolation because of their fundamental, sometimes reciprocal, relationships and influence in the long-term preservation and recovery of shell DNA. Successful recovery of shell DNA can facilitate research on pressing ecological and evolutionary questions and inform conservation strategies to protect molluscan diversity.

Pearl Sac Gene Expression Profiles Associated With Pearl Attributes in the Silver-Lip Pearl Oyster, Pinctada maxima

Pearls are highly prized biomineralized gemstones produced by molluscs. The appearance and mineralogy of cultured pearls can vary markedly, greatly affecting their commercial value. To begin to understand the role of pearl sacs—organs that form in host oysters from explanted mantle tissues that surround and synthesize pearls—we undertook transcriptomic analyses to identify genes that are differentially expressed in sacs producing pearls with different surface and structural characteristics. Our results indicate that gene expression profiles correlate with different pearl defects, suggesting that gene regulation in the pearl sac contributes to pearl appearance and quality. For instance, pearl sacs that produced pearls with surface non-lustrous calcification significantly down-regulate genes associated with cilia and microtubule function compared to pearl sacs giving rise to lustrous pearls. These results suggest that gene expression profiling can advance our understanding of processes that control biomineralization, which may be of direct value to the pearl industry, particularly in relation to defects that result in low value pearls.

Phylogenetic comparisons reveal mosaic histories of larval and adult shell matrix protein deployment in pteriomorph bivalves

Molluscan shells are organo-mineral composites, in which the dominant calcium carbonate is intimately associated with an organic matrix comprised mainly of proteins and polysaccharides. However, whether the various shell matrix proteins (SMPs) date to the origin of hard skeletons in the Cambrian, or whether they represent later deployment through adaptive evolution, is still debated. In order to address this issue and to better understand the origins and evolution of biomineralization, phylogenetic analyses have been performed on the three SMP families, Von Willebrand factor type A (VWA) and chitin-binding domain-containing protein (VWA-CB dcp), chitobiase, and carbonic anhydrase (CA), which exist in both larval and adult shell proteomes in the bivalves, Crassostrea gigas and Pinctada fucata. In VWA-CB dcp and chitobiase, paralogs for larval and adult SMPs evolved before the divergence of these species. CA-SMPs have been taken as evidence for ancient origins of SMPs by their presumed indispensable function in biomineralization and ubiquitous distribution in molluscs. However, our results indicate gene duplications that gave rise to separate deployments as larval and adult CA-SMPs occurred independently in each lineage after their divergence, which is considerably more recent than hitherto assumed, supporting the "recent heritage and fast evolution" scenario for SMP evolution.

Evolution of Biomineralization Genes in the Prismatic Layer of the Pen Shell Atrina pectinata

Molluscan shells are composed of calcium carbonates, with small amounts of extracellular matrices secreted from mantle epithelial cells. Many types of shell matrix proteins (SMPs) have been identified from molluscan shells or mantle cells. The pen shell Atrina pectinata (Pinnidae) has two different shell microstructures, the nacreous and prismatic layers. Nacreous and prismatic layer-specific matrix proteins have been reported in Pteriidae bivalves, but remain unclear in Pinnidae. We performed transcriptome analysis using the mantle cells of A. pectinata to screen the candidate transcripts involved in its prismatic layer formation. We found Asprich and nine highly conserved prismatic layer-specific SMPs encoding transcript in P. fucata, P. margaritifera, and P. maxima (Tyrosinase, Chitinase, EGF-like proteins, Fibronectin, valine-rich proteins, and prismatic uncharacterized shell protein 2 [PUSP2]) using molecular phylogenetic analysis or multiple alignment. We confirmed these genes were expressed in the epithelial cells of the mantle edge (outer surface of the outer fold) and the mantle pallium. Phylogenetic character mapping of these SMPs was used to infer a possible evolutionary scenario of them in Pteriomorphia. EGF-like proteins, Fibronectin, and valine-rich proteins encoding genes each evolved in the linage leading to four Pteriomorphia (Mytilidae, Pinnidae, Ostreidae, and Pteriidae), PUSP2 evolved in the linage leading to three Pteriomorphia families (Pinnidae, Ostreidae, and Pteriidae), and chitinase was independently evolved as SMPs in Mytilidae and in other Pteriomorphia (Pinnidae, Ostreidae, and Pteriidae). Our results provide a new dataset for A. pectinata SMP annotation, and a basis for understanding the evolution of prismatic layer formation in bivalves.

Molecular Characterization of a Novel Shell Matrix Protein With PDZ Domain From Mytilus coruscus

Mollusk shells are products of biomineralization and possess excellent mechanical properties, and shell matrix proteins (SMPs) have important functions in shell formation. A novel SMP with a PDZ domain (PDZ-domain-containing-protein-1, PDCP-1) was identified from the shell matrices of Mytilus coruscus. In this study, the gene expression, function, and location of PDCP-1 were analyzed. PDCP-1 was characterized as an ∼70 kDa protein with a PDZ (postsynaptic density/discs large/zonula occludes) domain and a ZM (ZASP-like motif) domain. The PDCP-1 gene has a high expression level and specific location in the foot, mantle and adductor muscle. Recombinantly expressed PDCP-1 (rPDCP-1) altered the morphology of calcite crystals, the polymorph of calcite crystals, binding with both calcite and aragonite crystals, and inhibition of the crystallization rate of calcite crystals. In addition, anti-rPDCP-1 antibody was prepared, and immunohistochemistry and immunofluorescence analyses revealed the specific location of PDCP-1 in the mantle, the adductor muscle, and the aragonite (nacre and myostracum) layer of the shell, suggesting multiple functions of PDCP-1 in biomineralization, muscle-shell attachment, and muscle attraction. Furthermore, pull-down analysis revealed 19 protein partners of PDCP-1 from the shell matrices, which accordingly provided a possible interaction network of PDCP-1 in the shell. These results expand the understanding of the functions of PDZ-domain-containing proteins (PDCPs) in biomineralization and the supramolecular chemistry that contributes to shell formation.

Characterization of the chalky layer-derived EGF-like domain-containing protein (CgELC) in the pacific oyster, Crassostrea gigas

The shells of the Pacific oyster Crassostrea gigas contain calcite crystals with three types of microstructures: prismatic, chalky, and foliated layers. Many shell matrix proteins were annotated from the shells of C. gigas; however, it is unclear which SMPs play important roles in their shell mineralization. The matrix proteins have never been reported from the chalky layer. In this study, we identified a chalky layer-derived EGF-like domain-containing protein (CgELC) from the chalky layer of C. gigas shells. The gene sequence of the CgELC was encoded under CGI_ 10,017,544 of the C. gigas genome database. Only peptide fragments in the N-terminal region of CGI_ 10,017,544 were detected by LC-MS/MS analyses, suggesting that CGI_ 10,017,544 was digested at the predicted protease digestion dibasic site by post-translational modification to become a mature CgELC protein. We produced three types of CgELC recombinant proteins, namely, the full length CgELC, as well as the N-terminal and C-terminal parts of CgELC (CgELC-N or -C, respectively), for in vitro crystallization experiments. In the presence of these recombinant proteins, the aggregation of polycrystalline calcite was observed. Some fibrous organic components seemed to be incorporated into the calcite crystals in the presence of the r-CgELC protein. We also noted different features in the crystallization between CgELC-N and CgELC-C; some crystals were inhibited crystal plane formation and contained many columnar prisms inside the crystals (CgELC-N) and formed numerous holes on their surfaces (CgELC-C). These results suggest that CgELC is involved in crystal aggregation and incorporated into calcite crystals.

Mollusk shell alterations resulting from coastal contamination and other environmental factors

Effects of contamination on aquatic organisms have been investigated and employed as biomarkers in environmental quality assessment for years. A commonly referenced aquatic organism, mollusks represent a group of major interest in toxicological studies. Both gastropods and bivalves have external mineral shells that protects their soft tissue from predation and desiccation. These structures are composed of an organic matrix and an inorganic matrix, both of which are affected by environmental changes, including exposure to hazardous chemicals. This literature review evaluates studies that propose mollusk shell alterations as biomarkers of aquatic system quality. The studies included herein show that changes to natural variables such as salinity, temperature, food availability, hydrodynamics, desiccation, predatory pressure, and substrate type may influence the form, structure, and composition of mollusk shells. However, in the spatial and temporal studies performed in coastal waters around the world, shells of organisms sampled from multi-impacted areas were found to differ in the form and composition of both organic and inorganic matrices relative to shells from less contaminated areas. Though these findings are useful, the toxicological studies were often performed in the field and were not able to attribute shell alterations to a specific molecule. It is known that the organic matrix of shells regulates the biomineralization process; proteomic analyses of shells may therefore elucidate how different contaminants affect shell biomineralization. Further research using approaches that allow a clearer distinction between shell alterations caused by natural variations and those caused by anthropogenic influence, as well as studies to identify which molecule is responsible for such alterations or to determine the ecological implications of shell alterations, are needed before any responses can be applied universally.

ProminTools: shedding light on proteins of unknown function in biomineralization with user friendly tools illustrated using mollusc shell matrix protein sequences

Biominerals are crucial to the fitness of many organism and studies of the mechanisms of biomineralization are driving research into novel materials. Biomineralization is generally controlled by a matrix of organic molecules including proteins, so proteomic studies of biominerals are important for understanding biomineralization mechanisms. Many such studies identify large numbers of proteins of unknown function, which are often of low sequence complexity and biased in their amino acid composition. A lack of user-friendly tools to find patterns in such sequences and robustly analyse their statistical properties relative to the background proteome means that they are often neglected in follow-up studies. Here we present ProminTools, a user-friendly package for comparison of two sets of protein sequences in terms of their global properties and motif content. Outputs include data tables, graphical summaries in an html file and an R-script as a starting point for data-set specific visualizations. We demonstrate the utility of ProminTools using a previously published shell matrix proteome of the giant limpet Lottia gigantea.

Computationally predicted gene regulatory networks in molluscan biomineralization identify extracellular matrix production and ion transportation pathways

MOTIVATION

The molecular processes regulating molluscan shell production remain relatively uncharacterized, despite the clear evolutionary and societal importance of biomineralization.

RESULTS

Here we built the first computationally predicted gene regulatory network (GRN) for molluscan biomineralization using Antarctic clam (Laternula elliptica) mantle gene expression data produced over an age-categorized shell damage-repair time-course. We used previously published in vivo in situ hybridization expression data to ground truth gene interactions predicted by the GRN and show that candidate biomineralization genes from different shell layers, and hence microstructures, were connected in unique modules. We characterized two biomineralization modules of the GRN and hypothesize that one module is responsible for translating the extracellular proteins required for growing, repairing or remodelling the nacreous shell layer, whereas the second module orchestrates the transport of both ions and proteins to the shell secretion site, which are required during normal shell growth, and repair. Our findings demonstrate that unbiased computational methods are particularly valuable for studying fundamental biological processes and gene interactions in non-model species where rich sources of gene expression data exist, but annotation rates are poor and the ability to carry out true functional tests are still lacking.

AVAILABILITY AND IMPLEMENTATION

The raw RNA-Seq data is freely available for download from NCBI SRA (Accession: PRJNA398984), the assembled and annotated transcriptome can be viewed and downloaded from molluscDB (ensembl.molluscdb.org) and in addition, the assembled transcripts, reconstructed GRN, modules and detailed annotations are all available as Supplementary Files.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mollusc shellomes: Past, present and future

In molluscs, the shell fabrication requires a large array of secreted macromolecules including proteins and polysaccharides. Some of them are occluded in the shell during mineralization process and constitute the shell repertoire. The protein moieties, also called shell proteomes or, more simply, 'shellomes', are nowadays analyzed via high-throughput approaches. These latter, applied so far on about thirty genera, have evidenced the huge diversity of shellomes from model to model. They also pinpoint the recurrent presence of functional domains of diverse natures. Shell proteins are not only involved in guiding the mineral deposition, but also in enzymatic and immunity-related functions, in signaling or in coping with many extracellular molecules such as saccharides. Many shell proteins exhibit low complexity domains, the function of which remains unclear. Shellomes appear as self-organizing systems that must be approached from the point of view of complex systems biology: at supramolecular level, they generate emergent properties, i.e., microstructures that cannot be simply explained by the sum of their parts. A conceptual scheme is developed here that reconciles the plasticity of the shellome, its evolvability and the constrained frame of microstructures. Other perspectives arising from the study of shellomes are briefly discussed, including the macroevolution of shell repertoires, their maturation and their transformation through time.

Shell palaeoproteomics: First application of peptide mass fingerprinting for the rapid identification of mollusc shells in archaeology

Molluscs were one of the most widely-used natural resources in the past, and their shells are abundant among archaeological findings. However, our knowledge of the variety of shells that were circulating in prehistoric times (and thus their socio-economic and cultural value) is scarce due to the difficulty of achieving taxonomic determination of fragmented and/or worked remains. This study aims to obtain molecular barcodes based on peptide mass fingerprints (PMFs) of intracrystalline proteins, in order to obtain shell identification. Palaeoproteomic applications on shells are challenging, due to low concentration of molluscan proteins and an incomplete understanding of their sequences. We explore different approaches for protein extraction from small-size samples (<20 mg), followed by MALDI-TOF-MS analysis. The SP3 (single-pot, solid-phase) sample preparation method was found to be the most successful in retrieving the intracrystalline protein fraction from seven molluscan shell taxa, which belong to different phylogenetic groups, possess distinct microstructures and are relevant for archaeology. Furthermore, all the shells analysed, including a 7000-year-old specimen of the freshwater bivalve Pseudunio, yielded good-quality distinctive spectra, demonstrating that PMFs can be used for shell taxon determination. Our work suggests good potential for large-scale screening of archaeological molluscan remains. SIGNIFICANCE: We characterise for the first time the peptide mass fingerprints of the intracrystalline shell protein fraction isolated from different molluscan taxa. We demonstrate that these proteins yield distinctive PMFs, even for shells that are phylogenetically related and/or that display similar microstructures. Furthermore, we extend the range of sample preparation approaches for "shellomics" by testing the SP3 method, which proved to be well-suited to shell protein extraction from small-size and protein-poor samples. This work thus lays the foundations for future large-scale applications for the identification of mollusc shells and other invertebrate remains from the archaeological and palaeontological records.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

The calcitonin-like system is an ancient regulatory system of biomineralization

Biomineralization is the process by which living organisms acquired the capacity to accumulate minerals in tissues. Shells are the biomineralized exoskeleton of marine molluscs produced by the mantle but factors that regulate mantle shell building are still enigmatic. This study sought to identify candidate regulatory factors of molluscan shell mineralization and targeted family B G-protein coupled receptors (GPCRs) and ligands that include calcium regulatory factors in vertebrates, such as calcitonin (CALC). In molluscs, CALC receptor (CALCR) number was variable and arose through lineage and species-specific duplications. The Mediterranean mussel (Mytilus galloprovincialis) mantle transcriptome expresses six CALCR-like and two CALC-precursors encoding four putative mature peptides. Mussel CALCR-like are activated in vitro by vertebrate CALC but only receptor CALCRIIc is activated by the mussel CALCIIa peptide (EC50 = 2.6 ×10-5 M). Ex-vivo incubations of mantle edge tissue and mantle cells with CALCIIa revealed they accumulated significantly more calcium than untreated tissue and cells. Mussel CALCIIa also significantly decreased mantle acid phosphatase activity, which is associated with shell remodelling. Our data indicate the CALC-like system as candidate regulatory factors of shell mineralization. The identification of the CALC system from molluscs to vertebrates suggests it is an ancient and conserved calcium regulatory system of mineralization.

Molecular characterization of a whirlin-like protein with biomineralization-related functions from the shell of Mytilus coruscus

Mollusc shells are produced from calcified skeletons and have excellent mechanical properties. Shell matrix proteins (SMPs) have important functions in shell formation. A 16.6 kDa whirlin-like protein (WLP) with a PDZ domain was identified in the shell of Mytilus coruscus as a novel SMP. In this study, the expression, function, and location of WLP were analysed. The WLP gene was highly expressed and specifically located in the adductor muscle and mantle. The expression of recombinant WLP (rWLP) was associated with morphological change, polymorphic change, binding ability, and crystallization rate inhibition of the calcium carbonate crystals in vitro. In addition, an anti-rWLP antibody was prepared, and the results from immunohistochemistry and immunofluorescence analyses revealed the specific location of the WLP in the mantle, adductor muscle, and myostracum layer of the shell, suggesting multiple functions for WLP in biomineralization, muscle-shell attachment, and muscle attraction. Furthermore, results from a pull-down analysis revealed 10 protein partners of WLP in the shell matrices and a possible network of interacting WLPs in the shell. In addition, in this study, one of the WLP partners, actin, was confirmed to have the ability to bind WLP. These results expand the understanding of the functions of PDZ-domain-containing proteins in biomineralization and provide clues for determining the mechanisms of myostracum formation and muscle-shell attachment.

The shell matrix of the european thorny oyster, Spondylus gaederopus: microstructural and molecular characterization

Molluscs, the largest marine phylum, display extraordinary shell diversity and sophisticated biomineral architectures. However, mineral-associated biomolecules involved in biomineralization are still poorly characterised. We report the first comprehensive structural and biomolecular study of Spondylus gaederopus, a pectinoid bivalve with a peculiar shell texture. Used since prehistoric times, this is the best-known shell of Europe's cultural heritage. We find that Spondylus microstructure is very poor in mineral-bound organics, which are mostly intercrystalline and concentrated at the interface between structural layers. Using high-resolution liquid chromatography tandem mass spectrometry (LC-MS/MS) we characterized several shell protein fractions, isolated following different bleaching treatments. Several peptides were identified as well as six shell proteins, which display features and domains typically found in biomineralized tissues, including the prevalence of intrinsically disordered regions. It is very likely that these sequences only partially represent the full proteome of Spondylus, considering the lack of genomics data for this genus and the fact that most of the reconstructed peptides do not match with any known shell proteins, representing consequently lineage-specific sequences. This work sheds light onto the shell matrix involved in the biomineralization in spondylids. Our proteomics data suggest that Spondylus has evolved a shell-forming toolkit, distinct from that of other better studied pectinoids - fine-tuned to produce shell structures with high mechanical properties, while limited in organic content. This study therefore represents an important milestone for future studies on biomineralized skeletons and provides the first reference dataset for forthcoming molecular studies of Spondylus archaeological artifacts.

Proteomic similarity of the Littorinid snails in the evolutionary context

BACKGROUND

The introduction of DNA-based molecular markers made a revolution in biological systematics. However, in cases of very recent divergence events, the neutral divergence may be too slow, and the analysis of adaptive part of the genome is more informative to reconstruct the recent evolutionary history of young species. The advantage of proteomics is its ability to reflect the biochemical machinery of life. It may help both to identify rapidly evolving genes and to interpret their functions.

METHODS

Here we applied a comparative gel-based proteomic analysis to several species from the gastropod family Littorinidae. Proteomes were clustered to assess differences related to species, geographic location, sex and body part, using data on presence/absence of proteins in samples and data on protein occurrence frequency in samples of different species. Cluster support was assessed using multiscale bootstrap resampling and the stability of clustering-using cluster-wise index of cluster stability. Taxon-specific protein markers were derived using IndVal method. Proteomic trees were compared to consensus phylogenetic tree (based on neutral genetic markers) using estimates of the Robinson-Foulds distance, the Fowlkes-Mallows index and cophenetic correlation.

RESULTS

Overall, the DNA-based phylogenetic tree and the proteomic similarity tree had consistent topologies. Further, we observed some interesting deviations of the proteomic littorinid tree from the neutral expectations. (1) There were signs of molecular parallelism in two Littoraria species that phylogenetically are quite distant, but live in similar habitats. (2) Proteome divergence was unexpectedly high between very closely related Littorina fabalis and L. obtusata, possibly reflecting their ecology-driven divergence. (3) Conservative house-keeping proteins were usually identified as markers for cryptic species groups ("saxatilis" and "obtusata" groups in the Littorina genus) and for genera (Littoraria and Echinolittorina species pairs), while metabolic enzymes and stress-related proteins (both potentially adaptively important) were often identified as markers supporting species branches. (4) In all five Littorina species British populations were separated from the European mainland populations, possibly reflecting their recent phylogeographic history. Altogether our study shows that proteomic data, when interpreted in the context of DNA-based phylogeny, can bring additional information on the evolutionary history of species.

PfmPif97-like regulated by Pfm-miR-9b-5p participates in shell formation in Pinctada fucata martensii

Mollusk shell matrix proteins are important for the formation of organic frameworks, crystal nucleation, and crystal growth in Pinctada fucata martensii (P. f. martensii). MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play important roles in many biological processes, including shell formation. In this study, we obtained the full-length sequence of Pif97-like gene in P. f. martensii (PfmPif97-like). PfmPif97-like was mainly distributed in mantle pallial and mantle edge. Correlation analysis indicated that the average shell thickness and weight showed a positive correlation with PfmPif97-like expression (P < 0.05). The inner surface of the nacreous layer and prismatic layer showed atypical growth when we knocked down the expression of PfmPif97-like by RNA interference (RNAi). We used a luciferase reporter assay to identify that miR-9b-5p of P. f. martensii (Pfm-miR-9b-5p) downregulated the expression of PfmPif97-like by interacting with the 3′-untranslated region (UTR) while we obtained the same result by injecting the Pfm-miR-9b-5p mimics in vivo. After injecting the mimics, we also observed abnormal growth in nacre layer and prismatic layer which is consistent with the result of RNAi. We proposed that PfmPif97-like regulated by Pfm-miR-9b-5p participates in shell formation of P. f. martensii. These findings provide important clues about the molecular mechanisms that regulate biomineralization in P. f. martensii.

Microstructure and in-depth proteomic analysis of Perna viridis shell

For understanding the structural characteristics and the proteome of Perna shell, the microstructure, polymorph, and protein composition of the adult Perna viridis shell were investigated. The P. viridis shell have two distinct mineral layers, myostracum and nacre, with the same calcium carbonate polymorph of aragonite, determined by scanning electron microscope, Fourier transform infrared spectroscopy, and x-ray crystalline diffraction. Using Illumina sequencing, the mantle transcriptome of P. viridis was investigated and a total of 69,859 unigenes was generated. Using a combined proteomic/transcriptomic approach, a total of 378 shell proteins from P. viridis shell were identified, in which, 132 shell proteins identified with more than two matched unique peptides. Of the 132 shell proteins, 69 are exclusive to the nacre, 12 to the myostracum, and 51 are shared by both. The Myosin-tail domain containing proteins, Filament-like proteins, and Chitin-binding domain containing proteins represent the most abundant molecules. In addition, the shell matrix proteins (SMPs) containing biomineralization-related domains, such as Kunitz, A2M, WAP, EF-hand, PDZ, VWA, Collagen domain, and low complexity regions with abundant certain amino acids, were also identified from P. viridis shell. Collagenase and chitinase degradation can significantly change the morphology of the shell, indicating the important roles of collagen and chitin in the shell formation and the muscle-shell attachment. Our results present for the first time the proteome of P. viridis shell and increase the knowledge of SMPs in this genus.

Evolution and diversity of alpha-carbonic anhydrases in the mantle of the Mediterranean mussel (Mytilus galloprovincialis)

The α-carbonic anhydrases (α-CAs) are a large and ancient group of metazoan-specific enzymes. They generate bicarbonate from metabolic carbon dioxide and through calcium carbonate crystal formation play a key role in the regulation of mineralized structures. To better understand how α-CAs contribute to shell mineralization in the marine Mediterranean mussel (Mytilus galloprovincialis) we characterized them in the mantle. Phylogenetic analysis revealed that mollusc α-CA evolution was affected by lineage and species-specific events. Ten α-CAs were found in the Mediterranean mussel mantle and the most abundant form was named, MgNACR, as it grouped with oyster nacreins (NACR). Exposure of the Mediterranean mussel to reduced water salinity (18 vs 37 ppt), caused a significant reduction (p < 0.05) in mantle esterase activity and MgNACR transcript abundance (p < 0.05). Protonograms revealed multiple proteins in the mantle with α-CA hydratase activity and mapped to a protein with a similar size to that deduced for monomeric MgNACR. Our data indicate that MgNACR is a major α-CA enzyme in mantle and that by homology with oyster nacreins likely regulates mussel shell production. We propose that species-dependent α-CA evolution may contribute to explain the diversity of bivalve shell structures and their vulnerability to environmental changes.

Dual Gene Repertoires for Larval and Adult Shells Reveal Molecules Essential for Molluscan Shell Formation

Molluscan shells, mainly composed of calcium carbonate, also contain organic components such as proteins and polysaccharides. Shell organic matrices construct frameworks of shell structures and regulate crystallization processes during shell formation. To date, a number of shell matrix proteins (SMPs) have been identified, and their functions in shell formation have been studied. However, previous studies focused only on SMPs extracted from adult shells, secreted after metamorphosis. Using proteomic analyses combined with genomic and transcriptomic analyses, we have identified 31 SMPs from larval shells of the pearl oyster, Pinctada fucata, and 111 from the Pacific oyster, Crassostrea gigas. Larval SMPs are almost entirely different from those of adults in both species. RNA-seq data also confirm that gene expression profiles for larval and adult shell formation are nearly completely different. Therefore, bivalves have two repertoires of SMP genes to construct larval and adult shells. Despite considerable differences in larval and adult SMPs, some functional domains are shared by both SMP repertoires. Conserved domains include von Willebrand factor type A (VWA), chitin-binding (CB), carbonic anhydrase (CA), and acidic domains. These conserved domains are thought to play crucial roles in shell formation. Furthermore, a comprehensive survey of animal genomes revealed that the CA and VWA-CB domain-containing protein families expanded in molluscs after their separation from other Lophotrochozoan linages such as the Brachiopoda. After gene expansion, some family members were co-opted for molluscan SMPs that may have triggered to develop mineralized shells from ancestral, nonmineralized chitinous exoskeletons.

A Kunitz proteinase inhibitor (HcKuPI) participated in antimicrobial process during pearl sac formation and induced the overgrowth of calcium carbonate in Hyriopsis cumingii

Proteinase inhibitors with the ability to inhibit specific proteinases are usually closely connected with the immune system. Interestingly, proteinase inhibitors are also a common ingredient in the organic matrix of mollusk shells. However, the molecular mechanism that underlies the role of proteinase inhibitors in immune system and shell mineralization is poorly known. In this study, a Kunitz serine proteinase inhibitor (HcKuPI) was isolated from the mussel Hyriopsis cumingii. HcKuPI was specifically expressed in dorsal epithelial cells of the mantle pallium and HcKuPI dsRNA injection caused an irregular surface and disordered deposition on the aragonite tablets of the nacreous layer. These results indicated that HcKuPI plays a vital role in shell nacreous layer biomineralization. Moreover, the expression pattern of HcKuPI during LPS challenge and pearl formation indicated its involvement in the antimicrobial process during pearl sac formation and nacre tablets accumulation during pearl formation. In the in vitro calcium carbonate crystallization assay, the addition of GST-HcKuPI increased the precipitation rate of calcium carbonate and induced the crystal overgrowth of calcium carbonate. Taken together, these results indicate that HcKuPI is involved in antimicrobial process during pearl formation, and participates in calcium carbonate deposition acceleration and morphological regulation of the crystals during nacreous layer formation. These findings extend our knowledge of the role of proteinase inhibitors in immune system and shell biomineralization.

'Palaeoshellomics' reveals the use of freshwater mother-of-pearl in prehistory

The extensive use of mollusc shell as a versatile raw material is testament to its importance in prehistoric times. The consistent choice of certain species for different purposes, including the making of ornaments, is a direct representation of how humans viewed and exploited their environment. The necessary taxonomic information, however, is often impossible to obtain from objects that are small, heavily worked or degraded. Here we propose a novel biogeochemical approach to track the biological origin of prehistoric mollusc shell. We conducted an in-depth study of archaeological ornaments using microstructural, geochemical and biomolecular analyses, including ‘palaeoshellomics’, the first application of palaeoproteomics to mollusc shells (and indeed to any invertebrate calcified tissue). We reveal the consistent use of locally-sourced freshwater mother-of-pearl for the standardized manufacture of ‘double-buttons’. This craft is found throughout Europe between 4200–3800 BCE, highlighting the ornament-makers’ profound knowledge of the biogeosphere and the existence of cross-cultural traditions.

Just like people do today, prehistoric humans liked to adorn themselves with beautiful objects. Shells, from creatures like clams and snails, were used to decorate clothing or worn as jewelry at least as far back as 100,000 years ago. Later people used shells as the raw materials to make beads or bracelets. Learning where the shells came from may help scientists understand why prehistoric people chose certain shells and not others. It may also offer clues about how they used natural resources and the cultural significance of these objects. But identifying the shells is difficult because they lose many of their original distinctive features when worked into ornaments.

New tools that use DNA or proteins to identify the raw materials used to craft ancient artifacts have emerged that may help. So far, scientists have mostly used these genomic and proteomic tools to identify the source of materials made from animal hide, ivory or bone – where collagen is the most abundant protein molecule. Yet it is more challenging to extract and characterize proteins or genetic material from mollusc shells. This is partly because the amount of proteins in shells is at least 300 times lower than in bone, and also because the makeup of proteins in shells is not as well-known as in collagen.

Sakalauskaite et al. have now overcome these issues by combining the analytical tools used to study the proteins and mineral content of modern shells with those of ancient protein research. They then used this approach, which they named palaeoshellomics, to extract proteins from seven “double-buttons” – pearl-like ornaments crafted by prehistoric people in Europe. The double-buttons were made between 4200 and 3800 BC and found at archeological sites in Denmark, Germany and Romania. Comparing the extracted proteins to those from various mollusc shells showed that the double-buttons were made from freshwater mussels belonging to a group known as the Unionoida.

The discovery helps settle a decade-long debate in archeology about the origin of the shells used to make double-buttons in prehistoric Europe. Ancient people often crafted ornaments from marine shells, because they were exotic and considered more prestigious. But the results on the double-buttons suggest instead that mother-of-pearl from fresh water shells was valued and used by groups throughout Europe, even those living in coastal areas. The palaeoshellomics technique used by Sakalauskaite et al. may now help identify the origins of shells from archeological and palaeontological sites.