Evolution of the tyrosinase gene family in bivalve molluscs: independent expansion of the mantle gene repertoire

Dual sgRNA-directed tyrosinases knockout using CRISPR/Cas9 technology in Pacific oyster (Crassostrea gigas) reveals their roles in early shell calcification

Biomineralization processes in bivalves, particularly the initial production of molecular components (such as matrix deposition and calcification) in the early stages of shell development are highly complex and well-organized. This study investigated the temporal dynamics of organic matrix and calcium carbonate (CaCO3) deposition in Pacific oysters (Crassostrea gigas) across various development stages. The shell-field initiated matrix secretion during the gastrula stage. Subsequent larval development triggered central shell-field calcification, accompanied by expansion of the calcium ring from its interior to the periphery. Notably, the expression patterns of CgTyrp-2 and CgTyr closely correlated with matrix deposition and calcification during early developmental stages, with peak expression occurring in oyster's gastrula and D-veliger stages. Subsequently, the CRISPR/Cas9 system was utilized to knock out CgTyrp-2 and CgTyr with more distinct phenotypic alterations observed when both genes were concurrently knocked out. The relative gene expression was analyzed post-knockout, indicating that the knockout of CgTyr or CgTyrp-2 led to reduced expression of CgChs1, along with increased expression of CgChit4. Furthermore, when dual-sgRNAs were employed to knockout CgTyrp-2, a large deletion (2 kb) within the CgTyrp-2 gene was identified. In summary, early shell formation in C. gigas is the result of a complex interplay of multiple molecular components with CgTyrp-2 and CgTyr playing key roles in regulating CaCO3 deposition.

Deciphering the molecular toolkit: regulatory elements governing shell biomineralization in marine molluscs

The intricate process of shell biomineralization in marine molluscs is governed by a complex interplay of regulatory elements, encompassing secretomes, transporters, and noncoding RNA. This review delves into recent advancements in understanding these regulatory mechanisms, emphasizing their significance in elucidating the functions and evolutionary dynamics of the molluscan shell biomineralization process. Central to this intricate orchestration are secretomes with diverse functional domains, selectively exported to the extrapallial space, which directly regulate crystal growth and morphology. Transporters are crucial for substrate transportation in the calcification and maintenance of cellular homeostasis. Beyond proteins and transporters, noncoding RNA molecules are integral components influencing shell biomineralization. This review underscores the nonnegligible roles played by these genetic elements at the molecular level. To comprehend the complexity of biomineralization in mollusc, we explore the origin and evolutionary history of regulatory elements, primarily secretomes. While some elements have recently evolved, others are ancient genes that have been co-opted into the biomineralization toolkit. These elements undergo structural and functional evolution through rapidly evolving repetitive low-complexity domains and domain gain/loss/rearrangements, ultimately shaping a distinctive set of secretomes characterized by both conserved features and evolutionary innovations. This comprehensive review enhances our understanding of molluscan biomineralization at the molecular and genetic levels.

A metallothionein from disk abalone (Haliotis discus discus): Insights into its functional roles in immune response, metal tolerance, and oxidative stress

Metallothioneins (MTs) are cysteine-rich metal-binding proteins whose expression is induced by exposure to essential and non-essential metals, making them potential biological markers for assessing metal pollution in various biomonitoring programs. However, the functional properties of these proteins are yet to be comprehensively characterized in most marine invertebrates. In this study, we identified and characterized an MT homolog from the disk abalone (Haliotis discus discus), referred to as disk abalone MT (AbMT). AbMT exhibited the same primary structural features as MTs from other mollusks containing two β-domains (β2β1-form). AbMT protein demonstrated metal-binding and detoxification abilities against Zn, Cu, and Cd, as evidenced by Escherichia coli growth kinetics, metal tolerance analysis, and UV absorption spectrum. Transcriptional analysis revealed that AbMT was ubiquitously expressed in all analyzed tissues and upregulated in gill tissue following challenge with Vibrio parahaemolyticus, Listeria monocytogenes, and viral hemorrhagic septicemia virus (VHSV). Additionally, overexpression of AbMT suppressed LPS-induced NO production in RAW264.7 macrophages, protected cells against H2O2-induced oxidative stress, and promoted macrophage polarization toward the M1 phase. Conclusively, these findings suggest an important role for AbMT in environmental stress protection and immune regulation in disk abalone.

Transcriptomic and Physiological Analysis Reveal Melanin Synthesis-Related Genes and Pathways in Pacific Oysters (Crassostrea gigas)

Shell color is one of the shell traits of molluscs, which has been regarded as an economic trait in some bivalves. Pacific oysters (Crassostrea gigas) are important aquaculture shellfish worldwide. In the past decade, several shell color strains of C. gigas were developed through selective breeding, which provides valuable materials for research on the inheritance pattern and regulation mechanisms of shell color. The inheritance patterns of different shell colors in C. gigas have been identified in certain research; however, the regulation mechanism of oyster pigmentation and shell color formation remains unclear. In this study, we performed transcriptomic and physiological analyses using black and white shell oysters to investigate the molecular mechanism of melanin synthesis in C. gigas. Several pigmentation-related pathways, such as cytochrome P450, melanogenesis, tyrosine metabolism, and the cAMP signaling pathway were found. The majority of differentially expressed genes and some signaling molecules from these pathways exhibited a higher level in the black shell oysters than in the white, especially after L-tyrosine feeding, suggesting that those differences may cause a variation of tyrosine metabolism and melanin synthesis. In addition, the in vitro assay using primary cells from mantle tissue showed that L-tyrosine incubation increased cAMP level, gene and protein expression, and melanin content. This study reveals the difference in tyrosine metabolism and melanin synthesis in black and white shell oysters and provides evidence for the potential regulatory mechanism of shell color in oysters.

Oxidation and cross-link of tyrosine-rich proteins are involved in the periostracum formation of the green mussel Perna viridis (Linnaeus)

Ocean acidification causes severe shell dissolution and threats the survival of marine molluscs. The periostracum in molluscs consists of macromolecules such as proteins and polysaccharides, and protects the inner shell layers from dissolution and microbial erosion. Moreover, it serves as the primary template for shell deposition. However, the chemical composition and formation mechanism of the periostracum is largely unknown. In this study, we applied transcriptomic, proteomics, physical, and chemical analysis to unravel the mysteries of the periostracum formation in the green mussel Perna viridis Linnaeus. FTIR analysis showed that the periostracum layer was an organic membrane mainly composed of polysaccharides, lipids, and proteins, similar to that of the shell matrix. Interestingly, the proteomic study identified components enriched in tyrosine and some enzymes that evolved in tyrosine oxidation, indicating that tyrosine oxidation might play an essential role in the periostracum formation. Moreover, comparative transcriptomics suggested that tyrosine-rich proteins were intensively synthesized in the periostracum groove. After being secreted, the periostracum proteins were gradually tanned by oxidation in the seawater, and the level of crosslink increased significantly as revealed by the ATR-FTIR. Our present study sheds light on the chemical composition and putative tanning mechanism of the periostracum layer in bivalve molluscs. SIGNIFICANCE: The periostracum layer, plays an essential role in the initiation of shell biomineralization, the protection of minerals from dissolution for molluscs and especially ocean acidification conditions in the changing global climate. However, the molecular mechanism underlying the periostracum formation is not fully understood. In this study, we revealed that the oxidation and cross-link of tyrosine-rich proteins by tyrosinase are involved in periostracum formation in the green mussel Perna viridis. This study provides some insights into the first step of mussel shell formation and the robust adaptation of P. viridis to diverse habitats. These findings also help to reveal the potential acclimation of bivalves to the projected acidifying seawater.

Single cell atlas of Xenoturbella bocki highlights limited cell-type complexity

Phylogenetic analyses over the last two decades have united a few small, and previously orphan clades, the nematodermatids, acoels and xenoturbelids, into the phylum Xenacoelomorpha. Some phylogenetic analyses support a sister relationship between Xenacoelomorpha and Ambulacraria (Xenambulacraria), while others suggest that Xenacoelomorpha may be sister to the rest of the Bilateria (Nephrozoa). An understanding of the cell type complements of Xenacoelomorphs is essential to assessing these alternatives as well as to our broader understanding of bilaterian cell type evolution. Employing whole organism single-cell RNA-seq in the marine xenacoelomorph worm Xenoturbella bocki, we show that Xenambulacrarian nerve nets share regulatory features and a peptidergic identity with those found in cnidarians and protostomes and more broadly share muscle and gland cell similarities with other metazoans. Taken together, these data are consistent with broad homologies of animal gland, muscle, and neurons as well as more specific affinities between Xenoturbella and acoel gut and epidermal tissues, consistent with the monophyly of Xenacoelomorpha.

Mitf involved in shell pigmentation by activating tyrosinase-mediated melanin synthesis in Pacific oyster (Crassostrea gigas)

Pigmentation is frequently observed in the molluscan shells, whereas the molecular regulation about these shell pigments formation is not clear. The microphthalmia-associated transcription factor (Mitf) is an important transactivator in melanin synthesis in vertebrates. Here, the Mitf containing a highly conserved basic helix-loop-helixleucine zipper (bHLH-LZ) domain was identified in an economically important marine bivalve Pacific oyster Crassostrea gigas. The Mitf was found to widespread tissue distribution and the expression was higher in the marginal mantle than in the central mantle. Particularly, the expression level of Mitf was high in black shell color oysters compared with white shell oysters. After injecting siRNA, the expression of Mitf decreased significantly, and the efficiency of RNA interference reached 53%. Besides, knockdown Mitf obviously decreased expression of tyrosinase family genes and tyrosinase activity of mantles, indicating a potential regulatory relationship between Mitf and Tyr or Typs. Simultaneously, there was a sharply reduce in the number of the melanosomes in the outer fold of mantle by silencing of Mitf. Luciferase assays in cell culture further verified that Mitf was involved in transcriptional regulation of Typ-2 and Typ-3 genes through binding to their specific promoter regions. These data argue that Mitf is involved in shell pigmentation through activating tyrosinase-mediated melanin synthesis in C. gigas.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Understanding Snail Mucus Biosynthesis and Shell Biomineralisation through Genomic Data Mining of the Reconstructed Carbohydrate and Glycan Metabolic Pathways of the Giant African Snail (Achatina fulica)

The giant African snail (Order Stylommatophora: Family Achatinidae), Achatina fulica (Bowdich, 1822), is the most significant and invasive land snail pest. The ecological adaptability of this snail involves high growth rate, reproductive capacity, and shell and mucus production, driven by several biochemical processes and metabolism. The available genomic information for A. fulica provides excellent opportunities to hinder the underlying processes of adaptation, mainly carbohydrate and glycan metabolic pathways toward the shell and mucus formation. The authors analysed the 1.78 Gb draft genomic contigs of A. fulica to identify enzyme-coding genes and reconstruct biochemical pathways related to the carbohydrate and glycan metabolism using a designed bioinformatic workflow. Three hundred and seventy-seven enzymes involved in the carbohydrate and glycan metabolic pathways were identified based on the KEGG pathway reference in combination with protein sequence comparison, structural analysis, and manual curation. Fourteen complete pathways of carbohydrate metabolism and seven complete pathways of glycan metabolism supported the nutrient acquisition and production of the mucus proteoglycans. Increased copy numbers of amylases, cellulases, and chitinases highlighted the snail advantage in food consumption and fast growth rate. The ascorbate biosynthesis pathway identified from the carbohydrate metabolic pathways of A. fulica was involved in the shell biomineralisation process in association with the collagen protein network, carbonic anhydrases, tyrosinases, and several ion transporters. Thus, our bioinformatic workflow was able to reconstruct carbohydrate metabolism, mucus biosynthesis, and shell biomineralisation pathways from the A. fulica genome and transcriptome data. These findings could reveal several evolutionary advantages of the A. fulica snail, and will benefit the discovery of valuable enzymes for industrial and medical applications.

Identifying genes associated with genetic control of color polymorphism in the pearl oyster Pinctada margaritifera var. cumingii (Linnaeus 1758) using a comparative whole genome pool-sequencing approach

For hundreds of years, the color diversity of Mollusca shells has been a topic of interest for humanity. However, the genetic control underlying color expression is still poorly understood in mollusks. The pearl oyster Pinctada margaritifera is increasingly becoming a biological model to study this process due to its ability to produce a large range of colors. Previous breeding experiments demonstrated that color phenotypes were partly under genetic control, and while a few genes were found in comparative transcriptomics and epigenetic experiments, genetic variants associated with the phenotypes have not yet been investigated. Here, we used a pooled-sequencing approach on 172 individuals to investigate color-associated variants on three color phenotypes of economic interest for pearl farming, in three wild and one hatchery populations. While our results uncovered SNPs targeting pigment-related genes already identified in previous studies, such as PBGD, tyrosinases, GST, or FECH, we also identified new color-related genes occurring in the same pathways, like CYP4F8, CYP3A4, and CYP2R1. Moreover, we identified new genes involved in novel pathways unknown to be involved in shell coloration for P. margaritifera, like the carotenoid pathway, BCO1. These findings are essential to possibly implement future breeding programs focused on individual selection for specific color production in pearl oysters and improve the footprint of perliculture on the Polynesian lagoon by producing less but with a better quality.

Evolutionary conservation and divergence of the transcriptional regulation of bivalve shell secretion across life-history stages

Adult molluscs produce shells with diverse morphologies and ornamentations, different colour patterns and microstructures. The larval shell, however, is a phenotypically more conserved structure. How do developmental and evolutionary processes generate varying diversity at different life-history stages within a species? Using live imaging, histology, scanning electron microscopy and transcriptomic profiling, we have described shell development in a heteroconchian bivalve, the Antarctic clam, Laternula elliptica, and compared it to adult shell secretion processes in the same species. Adult downstream shell genes, such as those encoding extracellular matrix proteins and biomineralization enzymes, were largely not expressed during shell development. Instead, a development-specific downstream gene repertoire was expressed. Upstream regulatory genes such as transcription factors and signalling molecules were largely conserved between developmental and adult shell secretion. Comparing heteroconchian data with recently reported pteriomorphian larval shell development data suggests that, despite being phenotypically more conserved, the downstream effectors constituting the larval shell 'tool-kit' may be as diverse as that of adults. Overall, our new data suggest that a larval shell formed using development-specific downstream effector genes is a conserved and ancestral feature of the bivalve lineage, and possibly more broadly across the molluscs.

Regulation of Tyrosinase Gene Expression by Retinoic Acid Pathway in the Pacific Oyster Crassostrea gigas

Retinoic acid (RA) plays important roles in various biological processes in animals. RA signaling is mediated by two types of nuclear receptors, namely retinoic acid receptor (RAR) and retinoid x receptor (RXR), which regulate gene expression by binding to retinoic acid response elements (RAREs) in the promoters of target genes. Here, we explored the effect of all-trans retinoic acid (ATRA) on the Pacific oyster Crassostera gigas at the transcriptome level. A total of 586 differentially expressed genes (DEGs) were identified in C. gigas upon ATRA treatment, with 309 upregulated and 277 downregulated genes. Bioinformatic analysis revealed that ATRA affects the development, metabolism, reproduction, and immunity of C. gigas. Four tyrosinase genes, including Tyr-6 (LOC105331209), Tyr-9 (LOC105346503), Tyr-20 (LOC105330910), and Tyr-12 (LOC105320007), were upregulated by ATRA according to the transcriptome data and these results were verified by real-time quantitative polymerase chain reaction (RT-qPCR) analysis. In addition, increased expression of Tyr (a melanin-related TYR gene in C. gigas) and Tyr-2 were detected after ATRA treatment. The yeast one-hybrid assay revealed the DNA-binding activity of the RA receptors CgRAR and CgRXR, and the interaction of CgRAR with RARE present in the Tyr-2 promoter. These results provide evidence for the further studies on the role of ATRA and the mechanism of RA receptors in mollusks.

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.

Expression of tyrosinase-like protein genes and their functional analysis in melanin synthesis of Pacific oyster (Crassostrea gigas)

Color polymorphism in Mollusca is of great interest for consumer preference. Although the heritability of shell color variation has been conducted by experimental crossing, little is known about molecular basis involved in these patterns. Tyrosinase-like proteins are important enzymes which are members of the type-3 copper protein superfamily. In this research, two tyrosinase-like protein genes including CgTyp-1 and CgTyp-3 were identified in the Pacific oyster Crassostrea gigas. Tissue expression analysis showed that CgTyp-1 and CgTyp-3 were dominantly expressed in the mantle. Particularly, they were expressed significantly higher in the edge mantle than that in the central mantle whether on the left or right mantles. Additionally, expressions of CgTyp-1 and CgTyp-3 were mainly found in the black shell color oysters, with relative lower levels in the white shell color oysters. In situ hybridization showed that positive signals for CgTyp-1 and CgTyp-3 were both detected within the outer epithelium of the outer fold either in the black or white shell color oysters. After interference, the expression levels of CgTyp-1 and CgTyp-3 mRNA were significantly attenuated, and the efficiency of RNAi reached 84.72% and 71.58%, respectively. Besides, knockdown CgTyp-1 or CgTyp-3, obviously decreased the tyrosinase activity of mantles. Furthermore, the number of the melanosomes within epithelium of the outer fold was sharply reduced by silencing of each Typ. These findings argue that CgTyp-1 and CgTyp-3 may be involved in the melanin synthesis, which lends insight into regulation mechanism of shell pigmentation in C. gigas.

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 Mantle Transcriptome of Chamelea gallina (Mollusca: Bivalvia) and Shell Biomineralization

Simple Summary

Chamelea gallina is a bivalve mollusc that represents one of the most important fishery resources in the Mediterranean basin. In this species, the thickness and sturdiness of the shell valves are two important characteristics as they are decisive for protection against predators and survival of specimens rejected in the sea because caught under commercial size. The aim of this work is to investigate the ability of this species to modulate the expression of genes encoding proteins involved in shell biomineralization process in response to abiotic and biotic factors. Our findings, obtained through a multidisciplinary approach, highlighted a different shell mineralization behaviour in C. gallina clams collected in sampling sites characterized by different salinity and food availability. Moreover, this study provided the first comprehensive transcriptome from mantle, the tissue responsible for shell formation. Therefore, these results contribute to increase knowledge on this process and might help in adopting ad hoc management plans for this fishery resource.

Abstract

The striped venus Chamelea gallina is a bivalve mollusc that represents one of the most important fishery resources of the Adriatic Sea. In this work, we investigated for the first time the ability of this species to modulate the expression of genes encoding proteins involved in biomineralization process in response to biotic and abiotic factors. We provided the first comprehensive transcriptome from the mantle tissue of clams collected in two sampling sites located along the Italian Adriatic coast and characterized by different environmental features. Moreover, the assessment of environmental parameters, scanning electron microscopy (SEM), and X-ray diffraction (XRD) measurements on valves were conducted to better contextualize RNA sequencing (RNA-Seq) data. Functional annotation of differentially expressed genes (DEGs) and SEM observations highlighted a different shell mineralization behaviour in C. gallina clams collected from two selected sites characterized by diverse environmental parameters.

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.

The Skeleton and Biomineralization Mechanism as Part of the Innate Immune System of Stony Corals

Stony corals are among the most important calcifiers in the marine ecosystem as they form the coral reefs. Coral reefs have huge ecological importance as they constitute the most diverse marine ecosystem, providing a home to roughly a quarter of all marine species. In recent years, many studies have shed light on the mechanisms underlying the biomineralization processes in corals, as characterizing the calicoblast cell layer and genes involved in the formation of the calcium carbonate skeleton. In addition, considerable advancements have been made in the research field of coral immunity as characterizing genes involved in the immune response to pathogens and stressors, and the revealing of specialized immune cells, including their gene expression profile and phagocytosis capabilities. Yet, these two fields of corals research have never been integrated. Here, we discuss how the coral skeleton plays a role as the first line of defense. We integrate the knowledge from both fields and highlight genes and proteins that are related to biomineralization and might be involved in the innate immune response and help the coral deal with pathogens that penetrate its skeleton. In many organisms, the immune system has been tied to calcification. In humans, immune factors enhance ectopic calcification which causes severe diseases. Further investigation of coral immune genes which are involved in skeleton defense as well as in biomineralization might shed light on our understanding of the correlation and the interaction of both processes as well as reveal novel comprehension of how immune factors enhance calcification.

Shell Biosynthesis and Pigmentation as Revealed by the Expression of Tyrosinase and Tyrosinase-like Protein Genes in Pacific Oyster (Crassostrea gigas) with Different Shell Colors

The widely recognized color polymorphisms of molluscan shell have been appreciated for hundreds of years by collectors and scientists, while molecular mechanisms underlying shell pigmentation are still poorly understood. Tyrosinase is a key rate-limiting enzyme for the biosynthesis of melanin. Here, we performed an extensive multi-omics data mining and identified two tyrosinase genes, including tyrosinase and tyrosinase-like protein 2 (Tyr and Typ-2 respectively), in the Pacific oyster Crassostrea gigas, and investigated the expression patterns of tyrosinase during adults and embryogenesis in black and white shell color C. gigas. Tissue expression analysis showed that two tyrosinase genes were both specifically expressed in the mantle, and the expression levels of Tyr and Typ-2 in the edge mantle were significantly higher than that in the central mantle. Besides, Tyr and Typ-2 genes were black shell-specific compared with white shell oysters. In situ hybridization showed that strong signals for Tyr were detected in the inner surface of the outer fold, whereas positive signals for Typ-2 were mainly localized in the outer surface of the outer fold. In the embryos and larvae, the high expression of Tyr mRNA was detected in eyed-larvae, while Typ-2 mRNA was mainly expressed at the trochophore and early D-veliger. Furthermore, the tyrosinase activity in the edge mantle was significantly higher than that in the central mantle. These findings indicated that Tyr gene may be involved in shell pigmentation, and Typ-2 is more likely to play critical roles not only in the formation of shell prismatic layer but also in shell pigmentation. In particular, Typ-2 gene was likely to involve in the initial non-calcified shell of trochophores. The work provides valuable information for the molecular mechanism study of shell formation and pigmentation in C. gigas.

Immune properties of invertebrate phenoloxidases

Melanin production from different types of phenoloxidases (POs) confers immunity from a variety of pathogens ranging from viruses and microorganisms to parasites. The arthropod proPO expresses a variety of activities including cytokine, opsonin and microbiocidal activities independent of and even without melanin production. Proteolytic processing of proPO and its activating enzyme gives rise to several peptide fragments with a variety of separate activities in a process reminiscent of vertebrate complement system activation although proPO bears no sequence similarity to vertebrate complement factors. Pathogens influence proPO activation and thereby what types of immune effects that will be produced. An increasing number of specialised pathogens - from parasites to viruses - have been identified who can synthesise compounds specifically aimed at the proPO-system. In invertebrates outside the arthropods phylogenetically unrelated POs are participating in melanization reactions obviously aimed at intruders and/or aberrant tissues.

Oyster biomineralization under ocean acidification: From genes to shell

Biomineralization is one of the key processes that is notably affected in marine calcifiers such as oysters under ocean acidification (OA). Understanding molecular changes in the biomineralization process under OA and its heritability, therefore, is key to developing conservation strategies for protecting ecologically and economically important oyster species. To do this, in this study, we have explicitly chosen the tissue involved in biomineralization (mantle) of an estuarine commercial oyster species, Crassostrea hongkongensis. The primary aim of this study is to understand the influence of DNA methylation over gene expression of mantle tissue under decreased ~pH 7.4, a proxy of OA, and to extrapolate if these molecular changes can be observed in the product of biomineralization-the shell. We grew early juvenile C. hongkongensis, under decreased ~pH 7.4 and control ~pH 8.0 over 4.5 months and studied OA-induced DNA methylation and gene expression patterns along with shell properties such as microstructure, crystal orientation and hardness. The population of oysters used in this study was found to be moderately resilient to OA at the end of the experiment. The expression of key biomineralization-related genes such as carbonic anhydrase and alkaline phosphatase remained unaffected; thus, the mechanical properties of the shell (shell growth rate, hardness and crystal orientation) were also maintained without any significant difference between control and OA conditions with signs of severe dissolution. In addition, this study makes three major conclusions: (1) higher expression of Ca²⁺ binding/signalling-related genes in the mantle plays a key role in maintaining biomineralization under OA; (2) DNA methylation changes occur in response to OA; however, these methylation changes do not directly control gene expression; and (3) OA would be more of a 'dissolution problem' rather than a 'biomineralization problem' for resilient species that maintain calcification rate with normal shell growth and mechanical properties.

The 'Shellome' of the Crocus Clam Tridacna crocea Emphasizes Essential Components of Mollusk Shell Biomineralization

Molluscan shells are among the most fascinating research objects because of their diverse morphologies and textures. The formation of these delicate biomineralized structures is a matrix-mediated process. A question that arises is what are the essential components required to build these exoskeletons. In order to understand the molecular mechanisms of molluscan shell formation, it is crucial to identify organic macromolecules in different shells from diverse taxa. In the case of bivalves, however, taxon sampling in previous shell proteomics studies are focused predominantly on representatives of the class Pteriomorphia such as pearl oysters, edible oysters and mussels. In this study, we have characterized the shell organic matrix from the crocus clam, Tridacna crocea, (Heterodonta) using various biochemical techniques, including SDS-PAGE, FT-IR, monosaccharide analysis, and enzyme-linked lectin assay (ELLA). Furthermore, we have identified a number of shell matrix proteins (SMPs) using a comprehensive proteomics approach combined to RNA-seq. The biochemical studies confirmed the presence of proteins, polysaccharides, and sulfates in the T. crocea shell organic matrix. Proteomics analysis revealed that the majority of the T. crocea SMPs are novel and dissimilar to known SMPs identified from the other bivalve species. Meanwhile, the SMP repertoire of the crocus clam also includes proteins with conserved functional domains such as chitin-binding domain, VWA domain, and protease inhibitor domain. We also identified BMSP (Blue Mussel Shell Protein, originally reported from Mytilus), which is widely distributed among molluscan shell matrix proteins. Tridacna SMPs also include low-complexity regions (LCRs) that are absent in the other molluscan genomes, indicating that these genes may have evolved in specific lineage. These results highlight the diversity of the organic molecules – in particular proteins – that are essential for molluscan shell formation.

Integrated Analysis of Coding Genes and Non-coding RNAs Associated with Shell Color in the Pacific Oyster (Crassostrea gigas)

Molluscan shell color polymorphism is important in genetic breeding, while the molecular information mechanism for shell coloring is unclear. Here, high-throughput RNA sequencing was used to compare expression profiles of coding and non-coding RNAs (ncRNAs) from Pacific oyster Crassostrea gigas with orange and black shell, which were from an F2 family constructed by crossing an orange shell male with a black shell female. First, 458, 13, and 8 differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified, respectively. Functional analysis suggested that the DEGs were significantly enriched in 9 pathways including tyrosine metabolism and oxidative phosphorylation pathways. Several genes related to melanin synthesis and biomineralization expressed higher whereas genes associated with carotenoid pigmentation or metabolism expressed lower in orange shell oyster. Then, based on the ncRNA analysis, 163 and 20 genes were targeted by 13 and 8 differentially expressed lncRNAs (DELs) and miRNAs (DEMs), severally. Potential DELs-DEMs-DEGs interactions were also examined. Seven DEMs-DEGs pairs were detected, in which tyrosinase-like protein 1 was targeted by lgi-miR-133-3p and lgi-miR-252a and cytochrome P450 was targeted by dme-miRNA-1-3p. These results revealed that melanin synthesis-related genes and miRNAs-mRNA interactions functioned on orange shell coloration, which shed light on the molecular regulation of shell coloration in marine shellfish.

A bivalve biomineralization toolbox

Mollusc shells are a result of the deposition of crystalline and amorphous calcite catalysed by enzymes and shell matrix proteins. Developing a detailed understanding of bivalve mollusc biomineralization pathways is complicated not only by the multiplicity of shell forms and microstructures in this class, but also by the evolution of associated proteins by domain co-option and domain shuffling. In spite of this, a minimal biomineralization toolbox comprising proteins and protein domains critical for shell production across species has been identified. Using a matched pair design to reduce experimental noise from inter-individual variation, combined with damage-repair experiments and a database of biomineralization shell matrix proteins (SMP) derived from published works, proteins were identified that are likely to be involved in shell calcification. Eighteen new, shared proteins likely to be involved in the processes related to the calcification of shells were identified by analysis of genes expressed during repair in Crassostrea gigas, Mytilus edulis and Pecten maximus. Genes involved in ion transport were also identified as potentially involved in calcification either via the maintenance of cell acid-base balance or transport of critical ions to the extrapallial space, the site of shell assembly. These data expand the number of candidate biomineralization proteins in bivalve molluscs for future functional studies and define a minimal functional protein domain set required to produce solid microstructures from soluble calcium carbonate. This is important for understanding molluscan shell evolution, the likely impacts of environmental change on biomineralization processes, materials science, and biomimicry research.

MolluscDB: a genome and transcriptome database for molluscs

As sequencing becomes more accessible and affordable, the analysis of genomic and transcriptomic data has become a cornerstone of many research initiatives. Communities with a focus on particular taxa or ecosystems need solutions capable of aggregating genomic resources and serving them in a standardized and analysis-friendly manner. Taxon-focussed resources can be more flexible in addressing the needs of a research community than can universal or general databases. Here, we present MolluscDB, a genome and transcriptome database for molluscs. MolluscDB offers a rich ecosystem of tools, including an Ensembl browser, a BLAST server for homology searches and an HTTP server from which any dataset present in the database can be downloaded. To demonstrate the utility of the database and verify the quality of its data, we imported data from assembled genomes and transcriptomes of 22 species, estimated the phylogeny of Mollusca using single-copy orthologues, explored patterns of gene family size change and interrogated the data for biomineralization-associated enzymes and shell matrix proteins. MolluscDB provides an easy-to-use and openly accessible data resource for the research community. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Genome-wide identification and expression profiling of TYR gene family in Ruditapes philippinarum under the challenge of Vibrio anguillarum

Tyrosinase (EC1.14.18.1, TYR) is also called phenol oxidase, is not only involved in pigmentation but also plays an important role in modulating innate immunity in invertebrates. Tyrosinase is a copper containing metalloenzyme. The tyrosinase protein has two copper binding sites and three conserved histidines. In this study, 21 tyrosinase genes (RpTYR) were obtained from the whole genome of Ruditapes philippinarum. Their open reading frames were from 951 to 5424 aa, the range of predicted relative molecular weight from 36.72 to 203.81 kDa, and the range of isoelectric point from 4.72 to 9.88. Transcriptome analysis showed that RpTYR gene was expressed specifically in different developmental stages, adult tissues, four strains and two groups with different shell colors. Besides, the expression profiles of 21 RpTYRs were investigated against the immune response of R. philippinarum to a Vibrio challenge. The qPCR results showed that RpTYRs were involved in the immune response of R. philippinarum after Vibrio anguillarum challenge. This study provides preliminary evidence that the tyrosinases genes are involved in the immune defense and the potential immune function of R. philippinarum. Overall, these findings suggested that the expansion of TYR genes may play vital roles in larval development, the formation of shell color pattern, and immune response in R. philippinarum.

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.

Identification of methionine -rich insoluble proteins in the shell of the pearl oyster, Pinctada fucata

The molluscan shell is a biomineral that comprises calcium carbonate and organic matrices controlling the crystal growth of calcium carbonate. The main components of organic matrices are insoluble chitin and proteins. Various kinds of proteins have been identified by solubilizing them with reagents, such as acid or detergent. However, insoluble proteins remained due to the formation of a solid complex with chitin. Herein, we identified these proteins from the nacreous layer, prismatic layer, and hinge ligament of Pinctada fucata using mercaptoethanol and trypsin. Most identified proteins contained a methionine-rich region in common. We focused on one of these proteins, NU-5, to examine the function in shell formation. Gene expression analysis of NU-5 showed that NU-5 was highly expressed in the mantle, and a knockdown of NU-5 prevented the formation of aragonite tablets in the nacre, which suggested that NU-5 was required for nacre formation. Dynamic light scattering and circular dichroism revealed that recombinant NU-5 had aggregation activity and changed its secondary structure in the presence of calcium ions. These findings suggest that insoluble proteins containing methionine-rich regions may be important for scaffold formation, which is an initial stage of biomineral formation.

Tracing key genes associated with the Pinctada margaritifera albino phenotype from juvenile to cultured pearl harvest stages using multiple whole transcriptome sequencing

BACKGROUND

Albino mutations are commonly observed in the animal kingdom, including in bivalves. In the black-lipped pearl oyster Pinctada margaritifera, albino specimens are characterized by total or partial absence of colouration resulting in typical white shell phenotype expression. The relationship of shell colour with resulting cultured pearl colour is of great economic interest in P. margaritifera, on which a pearl industry is based. Hence, the albino phenotype provides a useful way to examine the molecular mechanisms underlying pigmentation.

RESULTS

Whole transcriptome RNA-sequencing analysis comparing albino and black wild-type phenotypes at three stages over the culture cycle of P. margaritifera revealed a total of 1606, 798 and 187 differentially expressed genes in whole juvenile, adult mantle and pearl sac tissue, respectively. These genes were found to be involved in five main molecular pathways, tightly linked to known pigmentation pathways: melanogenesis, calcium signalling pathway, Notch signalling pathway, pigment transport and biomineralization. Additionally, significant phenotype-associated SNPs were selected (N = 159), including two located in the Pif biomineralization gene, which codes for nacre formation. Interestingly, significantly different transcript splicing was detected between juvenile (N = 1366) and adult mantle tissue (N = 313) in, e.g., the tyrosinase Tyr-1 gene, which showed more complex regulation in mantle, and the Notch1 encoding gene, which was upregulated in albino juveniles.

CONCLUSION

This multiple RNA-seq approach provided new knowledge about genes associated with the P. margaritifera albino phenotype, highlighting: 1) new molecular pathways, such as the Notch signalling pathway in pigmentation, 2) associated SNP markers with biomineraliszation gene of interest like Pif for marker-assisted selection and prevention of inbreeding, and 3) alternative gene splicing for melanin biosynthesis implicating tyrosinase.

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.

Ocean acidification inhibits initial shell formation of oyster larvae by suppressing the biosynthesis of serotonin and dopamine

Ocean acidification has severely affected the initial shell formation of marine bivalves during their larval stages. In the present study, it was found that dopamine (DA) content in early D-shape larvae was significantly higher than that in trochophore and D-shape larvae, while the serotonin (5-HT) content in early D-shape larvae and D-shape larvae was obviously higher than that in trochophore. Incubation of trochophore with 5-HT or DA could accelerate the formation of calcified shell, and the treatments with selective antagonists of receptors for 5-HT and DA (Cg5-HTR-1 and CgD1DR-1) obviously inhibited the formation of calcified shells. When oyster larvae were subjected to an experimental acidified treatment (pH 7.4), the biosynthesis of 5-HT and DA was inhibited, while the mRNA expression levels of the components in TGF-β pathway were significantly up-regulated in D-shape larvae. Moreover, the phosphorylation of TIR and the translocation of smad4 were hindered upon acidification treatments, and the expression patterns of chitinase and tyrosinase were completely reverted. These results collectively suggested that monoamine neurotransmitters 5-HT and DA could modulate the initial shell formation in oyster larvae through TGF-β smad pathway by regulating the expression of tyrosinase and chitinase to guarantee the chitin synthesis for shell formation. CO₂-induced seawater acidification could suppress the biosynthesis of 5-HT and DA, as well as the activation of TGF-β smad pathway, which would subvert the expression patterns of chitinase and tyrosinase and cause the failure of initial shell formation in oyster early D-shape larvae.

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO₃ crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO₃ precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

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.

Molecular mechanisms of biomineralization in marine invertebrates

Much recent marine research has been directed towards understanding the effects of anthropogenic-induced environmental change on marine biodiversity, particularly for those animals with heavily calcified exoskeletons, such as corals, molluscs and urchins. This is because life in our oceans is becoming more challenging for these animals with changes in temperature, pH and salinity. In the future, it will be more energetically expensive to make marine skeletons and the increasingly corrosive conditions in seawater are expected to result in the dissolution of these external skeletons. However, initial predictions of wide-scale sensitivity are changing as we understand more about the mechanisms underpinning skeletal production (biomineralization). These studies demonstrate the complexity of calcification pathways and the cellular responses of animals to these altered conditions. Factors including parental conditioning, phenotypic plasticity and epigenetics can significantly impact the production of skeletons and thus future population success. This understanding is paralleled by an increase in our knowledge of the genes and proteins involved in biomineralization, particularly in some phyla, such as urchins, molluscs and corals. This Review will provide a broad overview of our current understanding of the factors affecting skeletal production in marine invertebrates. It will focus on the molecular mechanisms underpinning biomineralization and how knowledge of these processes affects experimental design and our ability to predict responses to climate change. Understanding marine biomineralization has many tangible benefits in our changing world, including improvements in conservation and aquaculture and exploitation of natural calcified structure design using biomimicry approaches that are aimed at producing novel biocomposites.

Molecular adaptation of molluscan biomineralisation to high-CO2 oceans - The known and the unknown

High-CO₂ induced ocean acidification (OA) reduces the calcium carbonate (CaCO₃) saturation level (Ω) and the pH of oceans. Consequently, OA is causing a serious threat to several ecologically and economically important biomineralising molluscs. Biomineralisation is a highly controlled biochemical process by which molluscs deposit their calcareous structures. In this process, shell matrix proteins aid the nucleation, growth and assemblage of the CaCO₃ crystals in the shell. These molluscan shell proteins (MSPs) are, ultimately, responsible for determination of the diverse shell microstructures and mechanical strength. Recent studies have attempted to integrate gene and protein expression data of MSPs with shell structure and mechanical properties. These advances made in understanding the molecular mechanism of biomineralisation suggest that molluscs either succumb or adapt to OA stress. In this review, we discuss the fate of biomineralisation process in future high-CO₂ oceans and its ultimate impact on the mineralised shell's structure and mechanical properties from the perspectives of limited substrate availability theory, proton flux limitation model and the omega myth theory. Furthermore, studying the interplay of energy availability and differential gene expression is an essential first step towards understanding adaptation of molluscan biomineralisation to OA, because if there is a need to change gene expression under stressors, any living system would require more energy than usual. To conclude, we have listed, four important future research directions for molecular adaptation of molluscan biomineralisation in high-CO₂ oceans: 1) Including an energy budgeting factor while understanding differential gene expression of MSPs and ion transporters under OA. 2) Unraveling the genetic or epigenetic changes related to biomineralisation under stressors to help solving a bigger picture about future evolution of molluscs, and 3) Understanding Post Translational Modifications of MSPs with and without stressors. 4) Understanding carbon uptake mechanisms across taxa with and without OA to clarify the OA theories on Ω.

A Novel Tyrosinase Gene Plays a Potential Role in Modification the Shell Organic Matrix of the Triangle Mussel Hyriopsis cumingii

Although tyrosinases have been speculated to participate in the shell formation of mollusks, there is still a lack of experimental evidence to support this assumption. In this study, a novel tyrosinase designated HcTyr2 was isolated and characterized from the freshwater mussel Hyriopsis cumingii. The change in HcTyr2 mRNA expression during the process of embryonic development was detected by real-time quantitative PCR. The result showed that the expression of HcTyr2 mRNA was significantly upregulated at the stages of gastrulae and unmatured glochidia (P < 0.05), suggesting that this gene might fundamentally participate in the biogenesis and growth of the initial shell. Meanwhile, the upregulation of HcTyr2 mRNA at the stages of shell regeneration 24 h and 9 days after shell notching in the mantle edge (P < 0.05) implied that it might play an important role in shell periostracum and nacre formation by mediating the cross-linking of quinoproteins to promote the maturity of organic matrix. Additionally, the knockdown of HcTyr2 mRNA by RNA interference resulted in not only the suppression of periostracum growth but also structural disorder of nacre aragonite tablets, as detected by scanning electron microscopy. These results suggested that HcTyr2 might regulate the growth of shell by its oxidative ability to transform soluble matrix proteins into insoluble matrix proteins, then promoting the maturity of the shell organic framework in H. cumingii. In general, our results suggested the importance of HcTyr2 in the shell formation and regeneration of H. cumingii.

Characterization and expression of tyrosinase-like genes in the anemone Exaiptasia pallida as a function of health and symbiotic state

Coral disease is a major threat to reef ecosystems and therefore, understanding the cellular pathways underlying disease progression and resistance is critical to mitigating future outbreaks. This study focused on tyrosinase-like proteins in cnidarians, which contribute to melanin synthesis, an invertebrate innate immune defense. Specifically, characterization and phylogenetic analysis of cnidarian tyrosinases were performed, and their role in symbiosis and a "mystery disease" in the anemone Exaiptasia pallida was investigated using qPCR. The results reveal a diversity of tyrosinase-like proteins in cnidarians that separate into two major clades on a phylogenetic tree, suggesting functional divergence. Two E. pallida sequences, Ep_Tyr1 and Ep_Tyr2, were further investigated, and qPCR results revealed no gene expression differences as a function of symbiotic state, but decreased expression in late disease stages. Overall this work provides evidence for the participation of tyrosinases in the cnidarian immune response.

Characterization of the main steps in first shell formation in Mytilus galloprovincialis: possible role of tyrosinase

Bivalve biomineralization is a highly complex and organized process, involving several molecular components identified in adults and larval stages. However, information is still scarce on the ontogeny of the organic matrix before calcification occurs. In this work, first shell formation was investigated in the mussel Mytilus galloprovincialis. The time course of organic matrix and CaCO₃ deposition were followed at close times post fertilization (24, 26, 29, 32, 48 h) by calcofluor and calcein staining, respectively. Both components showed an exponential trend in growth, with a delay between organic matrix and CaCO₃ deposition. mRNA levels of genes involved in matrix deposition (chitin synthase; tyrosinase- TYR) and calcification (carbonic anhydrase; extrapallial protein) were quantified by qPCR at 24 and 48 hours post fertilization (hpf) with respect to eggs. All transcripts were upregulated across early development, with TYR showing highest mRNA levels from 24 hpf. TYR transcripts were closely associated with matrix deposition as shown by in situ hybridization. The involvement of tyrosinase activity was supported by data obtained with the enzyme inhibitor N-phenylthiourea. Our results underline the pivotal role of shell matrix in driving first CaCO₃ deposition and the importance of tyrosinase in the formation of the first shell in M. galloprovincialis.

The global effects of PmRunt co-located and co-expressed with a lincRNA lncRunt in pearl oyster Pinctada fucata martensii

Runt related transcription factors as trans-acting elements play critical roles in the developmental control of cell fate, hematopoiesis, bone formation and cancers. In previous study, the homologue of runt related transcription factor PmRunt has been identified from pearl oyster Pinctada fucata martensii and considered to play an important role in nacre formation. In this study, we used the same samples to perform RNA-seq to detect the global effects after the decrease of PmRunt expression. The transcription levels of several nacre shell matrix protein (NSMP) genes were significantly changed and the potential compensatory effect could happen internal gene families. Downregulation of PmRunt could also influence the biosynthesis of NSMPs through affecting amino acid metabolism, translation, protein processing and export. The inhibition of PmRunt also possibly affected the expression of caspases, IAPs and C1qs that related to apoptosis and immune. In addition, PmRunt highly expressed at 12 h and 12 d after transplantation in hemolymph, which was corresponded to transplantation immunity immune response and the morphology of pearl sac, suggested the cross-talk of biomineralization-immune regulation in hemocytes. Furthermore, a lincRNA (LncRunt) that co-located with PmRunt was identified and showed a significantly relative expression with PmRunt, which suggested the potential regulation. Therefore, these findings provided new idea to find the regulation targets of runt-related transcription factors and offers evidence of lncRNAs in potential biomineralization-immune regulation.

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.

microRNA-374 inhibits proliferation and promotes apoptosis of mouse melanoma cells by inactivating the Wnt signalling pathway through its effect on tyrosinase

Melanoma is one of the most malignant skin tumours with constantly increasing incidence worldwide. Previous studies have demonstrated that microRNA‐374 (miR‐374) is a novel biomarker for cancer therapy. Therefore, this study explores whether miR‐374 targeting tyrosinase (TYR) affects melanoma and its underlying mechanism. We constructed subcutaneous melanoma models to carry out the following experiments. The cells were transfected with a series of miR‐374 mimics, miR‐374 inhibitors or siRNA against TYR. Dual luciferase reporter gene assay was used for the verification of the targeting relationship between miR‐374 and TYR. Reverse transcription quantitative polymerase chain reaction and western blot analysis were conducted to determine the expression of miR‐374, TYR, β‐catenin, B‐cell leukaemia 2 (Bcl‐2), Bcl‐2 associated X protein (Bax), Low‐density lipoprotein receptor‐related protein 6 (LRP6), Leucine‐rich repeat G protein‐coupled receptor 5 (LGR5) and CyclinD1. Cell proliferation, migration, invasion, cell cycle distribution and apoptosis were evaluated using cell counting kit‐8 assay, scratch test, transwell assay and flow cytometry respectively. TYR was proved as a putative target of miR‐374 as the evidenced by the result. It was observed that up‐regulated miR‐374 or down‐regulated TYR increased expression of Bax and decreased expressions of TYR, β‐catenin, LRP6, Bcl‐2, CyclinD1 and LGR5, along with diminished cell proliferation, migration, invasion and enhanced apoptosis. Meanwhile, cells with miR‐374 inhibitors showed an opposite trend. These findings indicated that up‐regulated miR‐374 could inhibit the expression of TYR to suppress cell proliferation, migration, invasion and promote cell apoptosis in melanoma cells by inhibiting the Wnt signalling pathway.

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.

High-throughput identification of heavy metal binding proteins from the byssus of chinese green mussel (Perna viridis) by combination of transcriptome and proteome sequencing

The Byssus, which is derived from the foot gland of mussels, has been proved to bind heavy metals effectively, but few studies have focused on the molecular mechanisms behind the accumulation of heavy metals by the byssus. In this study, we integrated high-throughput transcriptome and proteome sequencing to construct a comprehensive protein database for the byssus of Chinese green mussel (Perna viridis), aiming at providing novel insights into the molecular mechanisms by which the byssus binds to heavy metals. Illumina transcriptome sequencing generated a total of 55,670,668 reads. After filtration, we obtained 53,047,718 clean reads and subjected them to de novo assembly using Trinity software. Finally, we annotated 73,264 unigenes and predicted a total of 34,298 protein coding sequences. Moreover, byssal samples were analyzed by proteome sequencing, with the translated protein database from the foot transcriptome as the reference for further prediction of byssal proteins. We eventually determined 187 protein sequences in the byssus, of which 181 proteins are reported for the first time. Interestingly, we observed that many of these byssal proteins are rich in histidine or cysteine residues, which may contribute to the byssal accumulation of heavy metals. Finally, we picked one representative protein, Pvfp-5-1, for recombinant protein synthesis and experimental verification of its efficient binding to cadmium (Cd2+) ions.

Insights into the Evolution of Shells and Love Darts of Land Snails Revealed from Their Matrix Proteins

Over the past decade, many skeletal matrix proteins that are possibly related to calcification have been reported in various calcifying animals. Molluscs are among the most diverse calcifying animals and some gastropods have adapted to terrestrial ecological niches. Although many shell matrix proteins (SMPs) have already been reported in molluscs, most reports have focused on marine molluscs, and the SMPs of terrestrial snails remain unclear. In addition, some terrestrial stylommatophoran snails have evolved an additional unique calcified character, called a "love dart," used for mating behavior. We identified 54 SMPs in the terrestrial snail Euhadra quaesita, and found that they contain specific domains that are widely conserved in molluscan SMPs. However, our results also suggest that some of them possibly have evolved independently by domain shuffling, domain recruitment, or gene co-option. We then identified four dart matrix proteins, and found that two of them are the same proteins as those identified as SMPs. Our results suggest that some dart matrix proteins possibly have evolved by independent gene co-option from SMPs during dart evolution events. These results provide a new perspective on the evolution of SMPs and "love darts" in land snails.

Mantle transcriptome sequencing of Mytilus spp. and identification of putative biomineralization genes

In molluscs, the shell secreted by mantle tissue during the biomineralization process is the first barrier against predators and mechanical damage. Changing environmental conditions, such as ocean acidification, influence shell strength and thus protection of the soft body within. Mussels are marine bivalves with important commercial and ecological value worldwide. Despite this importance, the proteins involved in the biomineralization and pigmentation processes in Mytilus spp. remain unclear, as does taxonomy of Mytilus taxa, though there have been many molecular studies. To further understanding in these areas, this study aimed to characterize and compare mantle transcriptomes of four mussel taxa using next generation sequencing. Mussels representing four taxa, were collected from several localities and RNA from mantle tissue was extracted. RNA sequences obtained were assembled, annotated and potential molecular markers, including simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were identified. Candidate contigs putatively related to biomineralization and pigmentation processes were then selected and several transcripts were chosen for phylogenetic analyses from the Bivalvia class. Transcriptome comparisons between Mytilus taxa, including gene ontology (GO) enrichment analysis and orthologues identification were performed. Of assembled contigs, 46.57%, 37.28% and 17.53% were annotated using NCBI NR, GO and Kyoto Encyclopedia of Genes and Genomes databases, respectively. Potential SSRs (483) and SNPs (1,497) were identified. Results presented a total of 1,292 contigs putatively involved in biomineralization and melanogenesis. Phylogenetic analyses of α-carbonic anhydrase, chitinase and tyrosinase revealed complex evolutionary history and diversity of these genes, which may be a result of duplication events or adaptation to different environments in mussels and other bivalves. Enrichment analyses revealed GO terms associated with pH and thermal response in Mytilus edulis from the North Sea and M. galloprovincialis from the Mediterranean Sea. The phylogenetic analysis within the genus Mytilus revealed M. californianus and M. coruscus to be genetically more distant from the other taxa: M. trossulus, M. edulis, M. chilensis and M. galloprovincialis. This work represents the first mantle transcriptome comparison between Mytilus taxa and provides contigs putatively involved in biomineralization.

The Molecular Differentiation of Anatomically Paired Left and Right Mantles of the Pacific Oyster Crassostrea gigas

Left-right (L-R) asymmetry is controlled by gene regulation pathways for the L-R axis, and in vertebrates, the gene Pitx2 in TGF-β signaling pathway plays important roles in the asymmetrical formation of organs. However, less is known about the asymmetries of anatomically identical paired organs, as well as the transcriptional regulation mechanism of the gene Pitx in invertebrates. Here, we report the molecular biological differences between the left and right mantles of an invertebrate, the Pacific oyster Crassostrea gigas, and propose one possible mechanism underlying those differences. RNA sequencing (RNA-seq) analysis indicated that the paired organs showed different gene expression patterns, suggesting possible functional differences in shell formation, pheromone signaling, nerve conduction, the stress response, and other physiological processes. RNA-seq and real-time qPCR analysis indicated high right-side expression of the Pitx homolog (cgPitx) in oyster mantle, supporting a conserved role for Pitx in controlling asymmetry. Methylation-dependent restriction-site associated DNA sequencing (MethylRAD) identified a methylation site in the promoter region of cgPitx and showed significantly different methylation levels between the left and right mantles. This is the first report, to our knowledge, of such a difference in methylation in spiralians, and it was further confirmed in 18 other individuals by using a pyrosequencing assay. The miRNome analysis and the TGF-β receptor/Smad inhibition experiment further supported that several genes in TGF-β signaling pathway may be related with the L/R asymmetry of oyster mantles. These results suggested that the molecular differentiation of the oyster's paired left and right mantles is significant, TGF-β signaling pathway could be involved in establishing or maintaining the asymmetry, and the cgPitx gene as one of genes in this pathway; the different methylation levels in its promoter regions between L/R mantles was the one of possible mechanisms regulating the left-right functional differentiation.

D1 dopamine receptor is involved in shell formation in larvae of Pacific oyster Crassostrea gigas

Dopamine (DA), a significant member of catecholamines, is reported to induce biomineralization of calcium carbonate vaterite microspheres via dopamine receptor (DR) in bivalves, implying the modulation of dopaminergic system on shell formation during larval development. In this research, a homologue of D1 type DR (CgD1DR-1) was identified from oyster Crassostrea gigas, whose full length cDNA was 1197 bp. It was widely expressed in various tissues of C. gigas, with the significantly higher levels in hepatopancreas, mantle, muscle and gill. During developmental stages, the mRNA transcripts of CgD1DR-1 in D-shape larvae were obviously higher (p < 0.05) than those in trochophore and umbo larvae, and CO₂ exposure could inhibit the synthesis of DA and mRNA expression of CgD1DR-1. After cell transfection and DA treatment, intracellular cAMP in cells with the expression of CgD1DR-1 increased significantly (p < 0.05). Furthermore, the incubation with SCH 23390 for the blockage of CgD1DR-1 significantly restrained the expressions of six shell formation-related genes including CgTyrosinase-1, CgTyrosinase-3, CgChitinaseLP, CgAMC, CgBMP and CgBMPR in trochophore and D-shape larvae. These results jointly suggested that DA together with its receptor CgD1DR-1 might be involved in shell formation during oyster larval development from trochophore to D-shape larvae, and CO₂-induced ocean acidification (OA) might influence marine bivalves by inhibiting the DA-D1DR pathway to prohibit their shell formation.

In-depth proteomic analyses of Haliotis laevigata (greenlip abalone) nacre and prismatic organic shell matrix

Background

The shells of various Haliotis species have served as models of invertebrate biomineralization and physical shell properties for more than 20 years. A focus of this research has been the nacreous inner layer of the shell with its conspicuous arrangement of aragonite platelets, resembling in cross-section a brick-and-mortar wall. In comparison, the outer, less stable, calcitic prismatic layer has received much less attention. One of the first molluscan shell proteins to be characterized at the molecular level was Lustrin A, a component of the nacreous organic matrix of Haliotis rufescens. This was soon followed by the C-type lectin perlucin and the growth factor-binding perlustrin, both isolated from H. laevigata nacre, and the crystal growth-modulating AP7 and AP24, isolated from H. rufescens nacre. Mass spectrometry-based proteomics was subsequently applied to to Haliotis biomineralization research with the analysis of the H. asinina shell matrix and yielded 14 different shell-associated proteins. That study was the most comprehensive for a Haliotis species to date.

Methods

The shell proteomes of nacre and prismatic layer of the marine gastropod Haliotis laevigata were analyzed combining mass spectrometry-based proteomics and next generation sequencing.

Results

We identified 297 proteins from the nacreous shell layer and 350 proteins from the prismatic shell layer from the green lip abalone H. laevigata. Considering the overlap between the two sets we identified a total of 448 proteins. Fifty-one nacre proteins and 43 prismatic layer proteins were defined as major proteins based on their abundance at more than 0.2% of the total. The remaining proteins occurred at low abundance and may not play any significant role in shell fabrication. The overlap of major proteins between the two shell layers was 17, amounting to a total of 77 major proteins.

Conclusions

The H. laevigata shell proteome shares moderate sequence similarity at the protein level with other gastropod, bivalve and more distantly related invertebrate biomineralising proteomes. Features conserved in H. laevigata and other molluscan shell proteomes include short repetitive sequences of low complexity predicted to lack intrinsic three-dimensional structure, and domains such as tyrosinase, chitin-binding, and carbonic anhydrase. This catalogue of H. laevigata shell proteins represents the most comprehensive for a haliotid and should support future efforts to elucidate the molecular mechanisms of shell assembly.

Electronic supplementary material

The online version of this article (10.1186/s12953-018-0139-3) contains supplementary material, which is available to authorized users.

The pearl oyster Pinctada fucata martensii genome and multi-omic analyses provide insights into biomineralization

Nacre, the iridescent material found in pearls and shells of molluscs, is formed through an extraordinary process of matrix-assisted biomineralization. Despite recent advances, many aspects of the biomineralization process and its evolutionary origin remain unknown. The pearl oyster Pinctada fucata martensii is a well-known master of biomineralization, but the molecular mechanisms that underlie its production of shells and pearls are not fully understood. We sequenced the highly polymorphic genome of the pearl oyster and conducted multi-omic and biochemical studies to probe nacre formation. We identified a large set of novel proteins participating in matrix-framework formation, many in expanded families, including components similar to that found in vertebrate bones such as collagen-related VWA-containing proteins, chondroitin sulfotransferases, and regulatory elements. Considering that there are only collagen-based matrices in vertebrate bones and chitin-based matrices in most invertebrate skeletons, the presence of both chitin and elements of collagen-based matrices in nacre suggests that elements of chitin- and collagen-based matrices have deep roots and might be part of an ancient biomineralizing matrix. Our results expand the current shell matrix-framework model and provide new insights into the evolution of diverse biomineralization systems.