A bivalve biomineralization toolbox

cAMP-Mediated CREM-MITF-TYR Axis Regulates Melanin Synthesis in Pacific Oysters

Colorful shells in bivalves are mostly caused by the presence of biological pigments, among which melanin is a key component in the formation of shell colours. Cyclic adenosine monophosphate (cAMP) is an important messenger in the regulation of pigmentation in some species. However, the role of cAMP in bivalve melanogenesis has not yet been reported. In this study, we performed in vitro and in vivo experiments to determine the role of cAMP in regulating melanogenesis in Pacific oysters. Besides, the function of cAMP-responsive element modulator (CREM) and the interactions between CREM and melanogenic genes were investigated. Our results showed that a high level of cAMP promotes the expression of melanogenic genes in Pacific oysters. CREM controls the expression of the MITF gene under cAMP regulation. In addition, CREM can regulate melanogenic gene expression, tyrosine metabolism, and melanin synthesis. These results indicate that cAMP plays an important role in the regulation of melanogenesis in Pacific oysters. CREM is a key transcription factor in the oyster melanin synthesis pathway, which plays a crucial role in oyster melanin synthesis through a cAMP-mediated CREM-MITF-TYR axis.

An historical "wreck": A transcriptome assembly of the naval shipworm, Teredo navalis Linnaeus, 1978

Historically famous for their negative impact on human-built marine wood structures, mollusc shipworms play a central ecological role in marine ecosystems. Their association with bacterial symbionts, providing cellulolytic and nitrogen-fixing activities, underscores their exceptional wood-eating and wood-boring behaviours, improving energy transfer and the recycling of essential nutrients locked in the wood cellulose. Importantly, from a molecular standpoint, a minute of omic resources are available from this lineage of Bivalvia. Here, we produced and assembled a transcriptome from the globally distributed naval shipworm, Teredo navalis (family Teredinidae). The transcriptome was obtained by sequencing the total RNA from five equidistant segments of the whole body of a T. navalis specimen. The quality of the produced assembly was accessed with several statistics, revealing a highly contiguous (1194 N50) and complete (over 90% BUSCO scores for Eukaryote and Metazoan databases) transcriptome, with nearly 38,000 predicted ORF, more than half being functionally annotated. Our findings pave the way to investigate the unique evolutionary biology of these highly modified bivalves and lay the foundation for an adequate gene annotation of a full genome sequence of the species.

Transcriptomics, proteomics, and physiological assays reveal immunosuppression in the eastern oyster Crassostrea virginica exposed to acidification stress

Ocean acidification (OA) is recognized as a major stressor for a broad range of marine organisms, particularly shell-building invertebrates. OA can cause alterations in various physiological processes such as growth and metabolism, although its effect on host-pathogen interactions remains largely unexplored. In this study, we used transcriptomics, proteomics, and physiological assays to evaluate changes in immunity of the eastern oyster Crassostrea virginica exposed to OA conditions (pH = 7.5 vs pH = 7.9) at various life stages. The susceptibility of oyster larvae to Vibrio infection increased significantly (131 % increase in mortality) under OA conditions, and was associated with significant changes in their transcriptomes. The significantly higher mortality of larvae exposed to pathogens and acidification stress could be the outcome of an increased metabolic demand to cope with acidification stress (as seen by upregulation of metabolic genes) at the cost of immune function (downregulation of immune genes). While larvae were particularly vulnerable, juveniles appeared more robust to the stressors and there were no differences in mortality after pathogen (Aliiroseovarius crassostrea and Vibrio spp.) exposure. Proteomic investigations in adult oysters revealed that acidification stress resulted in a significant downregulation of mucosal immune proteins including those involved in pathogen recognition and microbe neutralization, suggesting weakened mucosal immunity. Hemocyte function in adults was also impaired by high pCO2, with a marked reduction in phagocytosis (67 % decrease in phagocytosis) in OA conditions. Together, results suggest that OA impairs immune function in the eastern oyster making them more susceptible to pathogen-induced mortality outbreaks. Understanding the effect of multiple stressors such as OA and disease is important for accurate predictions of how oysters will respond to future climate regimes.

Biomineralization in a cold environment: Insights from shield compositions and transcriptomics of polar sternaspids (Sternaspidae, Polychaeta)

The survival and physiological functions of polar marine organisms are impacted by global climate changes. Investigation of the adaptation mechanisms underlying biomineralization in polar organisms at low temperatures is important for understanding mineralized organismal sensitivity to climate change. Here, we performed electron probe analysis on the shields of Antarctic polychaete Sternaspis sendalli and Arctic polychaete Sternaspis buzhinskajae (Sternaspidae), and sequenced the transcriptomes of the tissues surrounding shields to examine biomineral characteristics and adaptive mechanisms in persistently cold environments. Compared to the temperate relative species, the relative abundance of iron, phosphorus, calcium, magnesium, nitrogen, sulfur and silicon in two polar sternaspid shields was similar to Sternaspis chinensis. However, the diversity and expression levels of biomineralization-related shell matrix proteins differed between the polar and temperate species, suggesting distinct molecular mechanisms underlying shield formation in cold environments. Tubulin and cyclophilin were upregulated compared to the temperate species. Furthermore, 42 positively selected genes were identified in Antarctic S. sendalli, with functions in cytoskeletal structure, DNA repair, immunity, transcription, translation, protein synthesis, and lipid metabolism. Highly expressed genes in both polar species were associated with cytoskeleton, macromolecular complexes and cellular component biosynthesis. Overall, this study reveals conserved elemental composition yet distinct biomineralization processes in the shields of polar sternaspids. The unique expression of biomineralization related genes and other cold-adaptation related genes provide molecular insights into biomineralization in cold marine environments.

Tracing the evolution of tissue inhibitor of metalloproteinases in Metazoa with the Pteria penguin genome

Tissue inhibitors of metalloproteinase (TIMPs) play a pivotal role in regulating extracellular matrix (ECM) dynamics and have been extensively studied in vertebrates. However, understanding their evolution across invertebrate phyla is limited. Utilizing the high-quality Pteria penguin genome, we conducted phylogenomic orthology analyses across metazoans, revealing the emergence and distribution of the TIMP gene family. Our findings show that TIMP repertoires originated during eumetazoan radiation, experiencing independent duplication events in different clades, resulting in varied family sizes. Particularly, Pteriomorphia bivalves within Mollusca exhibited the most significant expansion and displayed the most diverse TIMP repertoires among metazoans. These expansions were attributed to multiple gene duplication events, potentially driven by the demands for functional diversification related to multiple adaptive traits, contributing to the adaptation of Pteriomorphia bivalves as stationary filter feeders. In this context, Pteriomorphia bivalves offer a promising model for studying invertebrate TIMP evolution.

Core genes of biomineralization and cis-regulatory long non-coding RNA regulate shell growth in bivalves

INTRODUCTION

Bivalve molluscs are abundant in marine and freshwater systems and contribute essential ecosystem services. They are characterized by an exuberant diversity of biomineralized shells and typically have two symmetric valves (a.k.a shells), but oysters (Ostreidae), some clams (Anomiidae and Chamidae) and scallops (Pectinida) have two asymmetrical valves. Predicting and modelling the likely consequences of ocean acidification on bivalve survival, biodiversity and aquaculture makes understanding shell biomineralization and its regulation a priority.

OBJECTIVES

This study aimed to a) exploit the atypical asymmetric shell growth of some bivalves and through comparative analysis of the genome and transcriptome pinpoint candidate biomineralization-related genes and regulatory long non-coding RNAs (LncRNAs) and b) demonstrate their roles in regulating shell biomineralization/growth.

METHODS

Meta-analysis of genomes, de novo generated mantle transcriptomes or transcriptomes and proteomes from public databases for six asymmetric to symmetric bivalve species was used to identify biomineralization-related genes. Bioinformatics filtering uncovered genes and regulatory modules characteristic of bivalves with asymmetric shells and identified candidate biomineralization-related genes and lncRNAs with a biased expression in asymmetric valves. A shell regrowth model in oyster and gene silencing experiments, were used to characterize candidate gene function.

RESULTS

Shell matrix genes with asymmetric expression in the mantle of the two valves were identified and unique cis-regulatory lncRNA modules characterized in Ostreidae. LncRNAs that regulate the expression of the tissue inhibitor of metalloproteinases gene family (TIMPDR) and of the shell matrix protein domain family (SMPDR) were identified. In vitro and in vivo silencing experiments revealed the candidate genes and lncRNA were associated with divergent shell growth rates and modified the microstructure of calcium carbonate (CaCO3) crystals.

CONCLUSION

LncRNAs are putative regulatory factors of the bivalve biomineralization toolbox. In the Ostreidae family of bivalves biomineralization-related genes are cis-regulated by lncRNA and modify the planar growth rate and spatial orientation of crystals in the shell.

Metabolic profiling of Mytilus coruscus mantle in response of shell repairing under acute acidification

Mytilus coruscus is an economically important marine bivalve mollusk found in the Yangtze River estuary, which experiences dramatic pH fluctuations due to seasonal freshwater input and suffer from shell fracture or injury in the natural environment. In this study, we used intact-shell and damaged-shell M. coruscus and performed metabolomic analysis, free amino acids analysis, calcium-positive staining, and intracellular calcium level tests in the mantle to investigate whether the mantle-specific metabolites can be induced by acute sea-water acidification and understand how the mantle responds to acute acidification during the shell repair process. We observed that both shell damage and acute acidification induced alterations in phospholipids, amino acids, nucleotides, organic acids, benzenoids, and their analogs and derivatives. Glycylproline, spicamycin, and 2-aminoheptanoic acid (2-AHA) are explicitly induced by shell damage. Betaine, aspartate, and oxidized glutathione are specifically induced by acute acidification. Our results show different metabolic patterns in the mussel mantle in response to different stressors, which can help elucidate the shell repair process under ocean acidification. furthermore, metabolic processes related to energy supply, cell function, signal transduction, and amino acid synthesis are disturbed by shell damage and/or acute acidification, indicating that both shell damage and acute acidification increased energy consumption, and disturb phospholipid synthesis, osmotic regulation, and redox balance. Free amino acid analysis and enzymatic activity assays partially confirmed our findings, highlighting the adaptation of M. coruscus to dramatic pH fluctuations in the Yangtze River estuary.

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

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

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.

Preliminary indication of the role of AHL-dependent quorum sensing systems in calcium carbonate precipitation in Gram-negative bacteria

Numerous microbial species participate in precipitation of carbonates in various natural environments, including soils, geological formations, freshwater biofilms and oceans. Despite the geochemical interest of such a biomineralization process, its molecular mechanisms and adaptive aspects remain poorly known. Many Gram-negative bacteria use cell-to-cell communication systems relying on N-acylhomoserine lactone (AHLs) signal molecules to express certain phenotypic traits in a density-dependent manner, a phenomenon referred as to quorum-sensing (QS). In this work, bacterial isolates collected from cave and rhizosphere soil were analyzed to study the occurrence of the AHL-mediated QS in bacterial calcium carbonate (CaCO3) precipitation. To test the production of AHLs signal molecules, we cross-streaked Gram-negative calcifying strains, selected among the environmental strains studied, with the AHL-negative mutant Chromobacterium subtsugae strain CV026. Only Burkholderia ambifaria LMG 11351 was able to restore violacein production in CV026 among the tested strains. The constructed AHL-negative mutant of B. ambifaria LMG 11351 could not precipitate CaCO3 on B-4 agar. Scanning Electron Microscopy (SEM) analysis on CaCO3 crystals obtained in vitro shows crystals of different morphologies, calcified biofilms and bacteria in close contact with the precipitated crystals. In the inner layers of the bioliths deposited by B. ambifaria LMG 11351, a stream-like organization of the Burkholderia imprints was not detected by SEM. Our data provide preliminary evidence that the activation of AHL-regulated genes may be a prerequisite for in vitro bacterial carbonatogenesis, in some cases, confirming the specific role of bacteria as CaCO3 precipitating agents. We enhance the understanding of bacterial CaCO3 biomineralization and its potential biotechnology implications for QS-based strategies to enhance or decrease CaCO3 precipitation through specific bacterial processes. The AHL-negative mutant of B. ambifaria LMG 11351 (a well-known plant growth-promoting bacterium) could also be used to study plant-bacteria interactions. The adaptive role of bacterial CaCO3 biomineralization was also discussed.

In situ mapping of biomineral skeletal proteins by molecular recognition imaging with antibody-functionalized AFM tips

Spatial localizing of skeletal proteins in biogenic minerals remains a challenge in biomineralization research. To address this goal, we developed a novel in situ mapping technique based on molecular recognition measurements via atomic force microscopy (AFM), which requires three steps: (1) the development and purification of a polyclonal antibody elicited against the target protein, (2) its covalent coupling to a silicon nitride AFM tip ('functionalization'), and (3) scanning of an appropriately prepared biomineral surface. We applied this approach to a soluble shell protein - accripin11 - recently identified as a major component of the calcitic prisms of the fan mussel Pinna nobilis [1]. Multiple tests reveal that accripin11 is evenly distributed at the surface of the prisms and also present in the organic sheaths surrounding the calcitic prisms, indicating that this protein is both intra- and inter-crystalline. We observed that the adhesion force in transverse sections is about twice higher than in longitudinal sections, suggesting that accripin11 may exhibit preferred orientation in the biomineral. To our knowledge, this is the first time that a protein is localized by molecular recognition atomic force microscopy with antibody-functionalized tips in a biogenic mineral. The 'pros' and 'cons' of this methodology are discussed in comparison with more 'classical' approaches like immunogold. This technique, which leaves the surface to analyze clean, might prove useful for clinical tests on non-pathological (bone, teeth) or pathological (kidney stone) biomineralizations. Studies using implants with protein-doped calcium phosphate coating can also benefit from this technology. STATEMENT OF SIGNIFICANCE: Our paper deals with an unconventional technical approach for localizing proteins that are occluded in biominerals. This technique relies on the use of molecular recognition atomic force microscopy with antibody-functionalized tips. Although such approach has been employed in other system, this is the very first time that it is developed for biominerals. In comparison to more classical approaches (such as immunogold), AFM microscopy with antibody-functionalized tips allows higher magnification and keeps the scanned surface clean for other biophysical characterizations. Our method has a general scope as it can be applied in human health, for non-pathological (bone, teeth) and pathological (kidney stone) biomineralizations as well as for bone implants coated with protein-doped calcium phosphate.

Biologically driven isotopic fractionations in bivalves: from palaeoenvironmental problem to palaeophysiological proxy

Traditional bulk stable isotope (δ18 O and δ13 C) and clumped isotope (Δ47 ) records from bivalve shells provide invaluable histories of Earth's local and global climate change. However, biologically driven isotopic fractionations (BioDIFs) can overprint primary environmental signals in the shell. Here, we explore how conventional measurements of δ18 O, δ13 C, and Δ47 in bivalve shells can be re-interpreted to investigate these physiological processes deliberately. Using intrashell Δ47 and δ18 O alignment as a proxy for equilibrium state, we separately examine fractionations and/or disequilibrium occurring in the two major stages of the biomineralisation process: the secretion of the extrapallial fluid (EPF) and the precipitation of shell material from the EPF. We measured δ18 O, δ13 C, and Δ47 in fossil shells representing five genera (Lahillia, Dozyia, Eselaevitrigonia, Nordenskjoldia, and Cucullaea) from the Maastrichtian age [66-69 million years ago (Ma)] López de Bertodano Formation on Seymour Island, Antarctica. Material was sampled from both the outer and inner shell layers (OSL and ISL, respectively), which precipitate from separate EPF reservoirs. We find consistent δ18 O values across the five taxa, indicating that the composition of the OSL can be a reliable palaeoclimate proxy. However, relative to the OSL baseline, ISLs of all taxa show BioDIFs in one or more isotopic parameters. We discuss/hypothesise potential origins of these BioDIFs by synthesising isotope systematics with the physiological processes underlying shell biomineralisation. We propose a generalised analytical and interpretive framework that maximises the amount of palaeoenvironmental and palaeobiological information that can be derived from the isotopic composition of fossil shell material, even in the presence of previously confounding 'vital effects'. Applying this framework in deep time can expand the utility of δ18 O, δ13 C, and Δ47 measurements from proxies of past environments to proxies for certain biomineralisation strategies across space, time, and phylogeny among Bivalvia and other calcifying organisms.

Outer fold is sole effective tissue among three mantle folds with regard to oyster shell colour

Molluscs constitute the second largest phylum of animals in the world, and shell colour is one of their most important phenotypic characteristics. In this study, we found among three folds on the mantle edge of oyster, only the outer fold had the same colour as the shell. Transcriptome and mantle cutting experiment indicated that the outer fold may be mainly reflected in chitin framework formation and biomineralisation. There were obvious differences in SEM structure and protein composition between the black and white shell periostraca. The black shell periostraca had more proteins related to melanin biosynthesis and chitin binding. Additionally, we identified an uncharacterized protein gene (named as CgCBP) ultra-highly expressed only in the black outer fold and confirmed its function of chitin-binding and CaCO₃ precipitation promoting. RNAi also indicated that CgCBP knockdown could change the structure of shell periostracum and reduce shell pigmentation. All these results suggest that the mantle outer fold plays multiple key roles in shell periostraca bioprocessing, and shell periostracum structure affected by chitin-binding protein is functionally correlated with shell pigmentation. The investigation of oyster shell periostracum structure and shell colour will provide a better understanding in pigmentation during biological mineralisation in molluscs.

Brachiopod and mollusc biomineralisation is a conserved process that was lost in the phoronid-bryozoan stem lineage

BACKGROUND

Brachiopods and molluscs are lophotrochozoans with hard external shells which are often believed to have evolved convergently. While palaeontological data indicate that both groups are descended from biomineralising Cambrian ancestors, the closest relatives of brachiopods, phoronids and bryozoans, are mineralised to a much lower extent and are comparatively poorly represented in the Palaeozoic fossil record. Although brachiopod and mollusc shells are structurally analogous, genomic and proteomic evidence indicates that their formation involves a complement of conserved, orthologous genes. Here, we study a set of genes comprised of 3 homeodomain transcription factors, one signalling molecule and 6 structural proteins which are implicated in mollusc and brachiopod shell formation, search for their orthologs in transcriptomes or genomes of brachiopods, phoronids and bryozoans, and present expression patterns of 8 of the genes in postmetamorphic juveniles of the rhynchonelliform brachiopod T. transversa.

RESULTS

Transcriptome and genome searches for the 10 target genes in the brachiopods Terebratalia transversa, Lingula anatina, Novocrania anomala, the bryozoans Bugula neritina and Membranipora membranacea, and the phoronids Phoronis australis and Phoronopsis harmeri resulted in the recovery of orthologs of the majority of the genes in all taxa. While the full complement of genes was present in all brachiopods with a single exception in L. anatina, a bloc of four genes could consistently not be retrieved from bryozoans and phoronids. The genes engrailed, distal-less, ferritin, perlucin, sp1 and sp2 were shown to be expressed in the biomineralising mantle margin of T. transversa juveniles.

CONCLUSIONS

The gene expression patterns we recovered indicate that while mineralised shells in brachiopods and molluscs are structurally analogous, their formation builds on a homologous process that involves a conserved complement of orthologous genes. Losses of some of the genes related to biomineralisation in bryozoans and phoronids indicate that loss of the capacity to form mineralised structures occurred already in the phoronid-bryozoan stem group and supports the idea that mineralised skeletons evolved secondarily in some of the bryozoan subclades.

Adaptive Bird-like Genome Miniaturization During the Evolution of Scallop Swimming Lifestyle

Genome miniaturization drives key evolutionary innovations of adaptive traits in vertebrates, such as the flight evolution of birds. However, whether similar evolutionary processes exist in invertebrates remains poorly understood. Derived from the second-largest animal phylum, scallops are a special group of bivalve molluscs and acquire the evolutionary novelty of the swimming lifestyle, providing excellent models for investigating the coordinated genome and lifestyle evolution. Here, we show for the first time that genome sizes of scallops exhibit a generally negative correlation with locomotion activity. To elucidate the co-evolution of genome size and swimming lifestyle, we focus on the Asian moon scallop (Amusium pleuronectes) that possesses the smallest known scallop genome while being among scallops with the highest swimming activity. Whole-genome sequencing of A. pleuronectes reveals highly conserved chromosomal macrosynteny and microsynteny, suggestive of a highly contracted but not degenerated genome. Genome reduction of A. pleuronectes is facilitated by significant inactivation of transposable elements, leading to reduced gene length, elevated expression of genes involved in energy-producing pathways, and decreased copy numbers and expression levels of biomineralization-related genes. Similar evolutionary changes of relevant pathways are also observed for bird genome reduction with flight evolution. The striking mimicry of genome miniaturization underlying the evolution of bird flight and scallop swimming unveils the potentially common, pivotal role of genome size fluctuation in the evolution of novel lifestyles in the animal kingdom.

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

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

Transcriptomic, Proteomic, and Functional Assays Underline the Dual Role of Extrapallial Hemocytes in Immunity and Biomineralization in the Hard Clam Mercenaria mercenaria

Circulating hemocytes in the hemolymph represent the backbone of innate immunity in bivalves. Hemocytes are also found in the extrapallial fluid (EPF), the space delimited between the shell and the mantle, which is the site of shell biomineralization. This study investigated the transcriptome, proteome, and function of EPF and hemolymph in the hard clam Mercenaria mercenaria. Total and differential hemocyte counts were similar between EPF and hemolymph. Overexpressed genes in the EPF were found to have domains previously identified as being part of the "biomineralization toolkit" and involved in bivalve shell formation. Biomineralization related genes included chitin-metabolism genes, carbonic anhydrase, perlucin, and insoluble shell matrix protein genes. Overexpressed genes in the EPF encoded proteins present at higher abundances in the EPF proteome, specifically those related to shell formation such as carbonic anhydrase and insoluble shell matrix proteins. Genes coding for bicarbonate and ion transporters were also overexpressed, suggesting that EPF hemocytes are involved in regulating the availability of ions critical for biomineralization. Functional assays also showed that Ca²⁺ content of hemocytes in the EPF were significantly higher than those in hemolymph, supporting the idea that hemocytes serve as a source of Ca²⁺ during biomineralization. Overexpressed genes and proteins also contained domains such as C1q that have dual functions in biomineralization and immune response. The percent of phagocytic granulocytes was not significantly different between EPF and hemolymph. Together, these findings suggest that hemocytes in EPF play a central role in both biomineralization and immunity.

Revisiting tolerance to ocean acidification: Insights from a new framework combining physiological and molecular tipping points of Pacific oyster

Studies on the impact of ocean acidification on marine organisms involve exposing organisms to future acidification scenarios, which has limited relevance for coastal calcifiers living in a mosaic of habitats. Identification of tipping points beyond which detrimental effects are observed is a widely generalizable proxy of acidification susceptibility at the population level. This approach is limited to a handful of studies that focus on only a few macro-physiological traits, thus overlooking the whole organism response. Here we develop a framework to analyze the broad macro-physiological and molecular responses over a wide pH range in juvenile oyster. We identify low tipping points for physiological traits at pH 7.3-6.9 that coincide with a major reshuffling in membrane lipids and transcriptome. In contrast, a drop in pH affects shell parameters above tipping points, likely impacting animal fitness. These findings were made possible by the development of an innovative methodology to synthesize and identify the main patterns of variations in large -omic data sets, fitting them to pH and identifying molecular tipping points. We propose the broad application of our framework to the assessment of effects of global change on other organisms.

Mussels Repair Shell Damage despite Limitations Imposed by Ocean Acidification

Bivalves frequently withstand shell damage that must be quickly repaired to ensure survival. While the processes that underlie larval shell development have been extensively studied within the context of ocean acidification (OA), it remains unclear whether shell repair is impacted by elevated pCO2. To better understand the stereotypical shell repair process, we monitored mussels (Mytilus edulis) with sublethal shell damage that breached the mantle cavity within both field and laboratory conditions to characterize the deposition rate, composition, and integrity of repaired shell. Results were then compared with a laboratory experiment wherein mussels (Mytilus trossulus) repaired shell damage in one of seven pCO2 treatments (400–2500 µatm). Shell repair proceeded through distinct stages; an organic membrane first covered the damaged area (days 1–15), followed by the deposition of calcite crystals (days 22–43) and aragonite tablets (days 51–69). OA did not impact the ability of mussels to close drill holes, nor the microstructure, composition, or integrity of end-point repaired shell after 10 weeks, as measured by µCT and SEM imaging, energy-dispersive X-ray (EDX) analysis, and mechanical testing. However, significant interactions between pCO2, the length of exposure to treatment conditions, the strength and inorganic content of shell, and the physiological condition of mussels within OA treatments were observed. These results suggest that while OA does not prevent adult mussels from repairing or mineralizing shell, both OA and shell damage may elicit stress responses that impose energetic constraints on mussel physiology.

A new gene family diagnostic for intracellular biomineralization of amorphous Ca-carbonates by cyanobacteria

Cyanobacteria have massively contributed to carbonate deposition over the geological history. They are traditionally thought to biomineralize CaCO3 extracellularly as an indirect byproduct of photosynthesis. However, the recent discovery of freshwater cyanobacteria forming intracellular amorphous calcium carbonates (iACC) challenges this view. Despite the geochemical interest of such a biomineralization process, its molecular mechanisms and evolutionary history remain elusive. Here, using comparative genomics, we identify a new gene (ccyA) and protein family (calcyanin) possibly associated with cyanobacterial iACC biomineralization. Proteins of the calcyanin family are composed of a conserved C-terminal domain, which likely adopts an original fold, and a variable N-terminal domain whose structure allows differentiating 4 major types among the 35 known calcyanin homologs. Calcyanin lacks detectable full-length homologs with known function. The overexpression of ccyA in iACC-lacking cyanobacteria resulted in an increased intracellular Ca content. Moreover, ccyA presence was correlated and/or co-localized with genes involved in Ca or HCO3- transport and homeostasis, supporting the hypothesis of a functional role of calcyanin in iACC biomineralization. Whatever its function, ccyA appears as diagnostic of intracellular calcification in cyanobacteria. By searching for ccyA in publicly available genomes, we identified 13 additional cyanobacterial strains forming iACC, as confirmed by microscopy. This extends our knowledge about the phylogenetic and environmental distribution of cyanobacterial iACC biomineralization, especially with the detection of multicellular genera as well as a marine species. Moreover, ccyA was probably present in ancient cyanobacteria, with independent losses in various lineages that resulted in a broad but patchy distribution across modern cyanobacteria.