Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks

In silico identification of neuropeptide genes encoded by the genome of Crassostrea virginica with a special emphasis on feeding-related genes

Suspension-feeding bivalves, including the oyster Crassostrea virginica, use mucosal lectins to capture food particles. For instance, oysters can increase the transcription of these molecules to enhance food uptake. However, the regulatory processes influencing food uptake remain unclear although likely involve neuropeptides. Information on the neuropeptidome of C. virginica is limited, hindering the comprehension of its physiology, including energy homeostasis. This study explored the genome of C. virginica to identify neuropeptide precursors in silico and compared these with orthologs from other mollusks. A special focus was given to genes with potential implication in feeding processes. qPCR was used to determine the main organs of transcription of feeding-related genes. To further probe the function of target neuropeptides, visceral ganglia extracts and synthetic NPF were injected into oysters to evaluate their impact on genes associated with feeding and energy homeostasis. A total of eighty-five neuropeptides genes were identified in C. virginica genome. About 50 % of these are suggested to play a role in feeding processes. qPCR analyses showed that visceral ganglia and digestive system are the main organs for the synthesis of feeding-related neuropeptides. Further, results showed that the transcription of several neuropeptide genes in the visceral ganglia, including NPF and insulin-like peptide, increased after starvation. Finally, the injection of visceral ganglia extracts and synthetic NPF increased the transcription of a mucosal lectin and a glycogen synthase, known to be involved in food capture and glucose storage. Overall, this study identifies key genes regulating oyster physiology, enhancing the understanding of the control of basic physiological mechanisms in C. virginica.

Development and validation of a 66K SNP array for the hard clam (Mercenaria mercenaria)

BACKGROUND

The hard clam (Mercenaria mercenaria), a marine bivalve distributed along the U.S. eastern seaboard, supports a significant shellfish industry. Overharvest in the 1970s and 1980s led to a reduction in landings. While the transition of industry from wild harvest to aquaculture since that time has enhanced production, it has also exacerbated challenges such as disease outbreaks. In this study, we developed and validated a 66K SNP array designed to advance genetic studies and improve breeding programs in the hard clam, focusing particularly on the development of markers that could be useful in understanding disease resistance and environmental adaptability.

RESULTS

Whole-genome resequencing of 84 individual clam samples and 277 pooled clam libraries yielded over 305 million SNPs, which were filtered down to a set of 370,456 SNPs that were used as input for the design of a 66K SNP array. This medium-density array features 66,543 probes targeting coding and non-coding regions, including 70 mitochondrial SNPs, to capture the extensive genetic diversity within the species. The SNPs were distributed evenly throughout the clam genome, with an average interval of 25,641 bp between SNPs. The array incorporates markers for detecting the clam pathogen Mucochytrium quahogii (formerly QPX), enhancing its utility in disease management. Performance evaluation on 1,904 samples demonstrated a 72.7% pass rate with stringent quality control. Concordance testing affirmed the array's repeatability, with an average agreement of allele calls of 99.64% across multiple tissue types, highlighting its reliability. The tissue-specific analysis demonstrated that some tissue types yield better genotyping results than others. Importantly, the array, including its embedded mitochondrial markers, effectively elucidated complex genetic relationships across different clam groups, both wild populations and aquacultured stocks, showcasing its utility for detailed population genetics studies.

CONCLUSIONS

The 66K SNP array is a powerful and robust genotyping tool that offers unprecedented insights into the species' genomic architecture and population dynamics and that can greatly facilitate hard clam selective breeding. It represents an important resource that has the potential to transform clam aquaculture, thereby promoting industry sustainability and ecological and economic resilience.

Bioremediation potential of the hard clam Mercenaria mercenaria as an intensive shrimp aquaculture pond polyculture condidate

Polyculture practices are important for achieving sustainable aquaculture development. Recently, hard clams polyculture in intensive shrimp ponds has been encouraged because bivalves can consume excess nutrients in aquaculture systems and sequester carbon. To evaluate the bioremediation potential of hard clams polyculture in intensive shrimp ponds, this study built an assessment model based on individual growth models and estimated the potential for nitrogen and phosphorus removal as well as CO2 fixation by hard clams. Firstly, key parameters required for model construction were obtained through field surveys and physiological experiments. Subsequently, an individual growth model for the hard clam Mercenaria mercenaria was developed based on the Dynamic Energy Budget (DEB) theory. Fitting of the growth data indicated that the model accurately replicated the growth patterns of hard clams, with relative root mean square errors of 9.87 % for shell length and 5.02 % for dry tissue weight. Finally, the assessment model for the bioremediation potential of hard clams demonstrated that, over 110 days in the intensive shrimp mariculture pond, the net removal of nitrogen and phosphorus by hard clams were 3.68 kg ha-1 and 0.81 kg ha-1, respectively, and CO2 fixation was 507.00 kg ha-1. These findings suggested that the DEB model is an effective tool for evaluating bivalve ecological remediation potential and can aid in selecting species for sustainable polyculture.

Identification and expression analysis of two steamer-like retrotransposons in the Chilean blue mussel (Mytilus chilensis)

BACKGROUND

Disseminated neoplasia (DN) is a proliferative cell disorder of the circulatory system of bivalve mollusks. The disease is transmitted between individuals and can also be induced by external chemical agents such as bromodeoxyuridine. In Mya arenaria, we have cloned and characterized an LTR-retrotransposon named Steamer. Steamer mRNA levels and gene copy number correlates with DN and can be used as a marker of the disease. So far, the only mollusk where a retrotransposon expression relates to DN is Mya arenaria. On the other hand, it has been reported that the Chilean blue mussel Mytilus chilensis can also suffers DN. Our aim was to identify retrotransposons in Mytilus chilensis and to study their expression levels in the context of disseminated neoplasia.

RESULTS

Here we show that 7.1% of individuals collected in August 2018, from two farming areas, presents morphological characteristics described in DN. Using Steamer sequence to interrogate the transcriptome of M. chilensis we found two putative retrotransposons, named Steamer-like elements (MchSLEs). MchSLEs are present in the genome of M. chilensis and MchSLE1 is indeed an LTR-retrotransposon. Neither expression, nor copy number of the reported MchSLEs correlate with DN status but both are expressed at different levels among individual animals. We also report that in cultured M. chilensis haemocytes MchSLEs1 expression can be induced by bromodeoxyuridine.

CONCLUSIONS

We conclude that SLEs present in Mytilus chilensis are differentially expressed among individuals and do not correlate with disseminated neoplasia. Treatment of haemocytes with a stressor like bromodeoxyuridine induces expression of MchSLE1 suggesting that in Mytilus chilensis environmental stressors can induce activation of LTR-retrotransposon.

Discovery and biochemical characterization of two hexokinases from Crassostrea gigas

Hexokinases (HKs) play a vital role in glucose metabolism, which controls the first committed step catalyzing the production of glucose-6-phosphate from glucose. Two HKs (CGIHK1 and CGIHK2) from the Pacific oyster Crassostrea giga were cloned and characterized. CGIHK1 and CGIHK2 were recombinantly expressed in Escherichia coli and successfully purified by the Ni-NTA column. The optimum pH of the two enzymes was pH 8.0 and 8.5, respectively. The optimum temperature of the two enzymes was 42 °C and 50 °C, respectively. Both enzymes showed a clear requirement for divalent magnesium and were strongly inhibited by SDS. CGIHK1 exhibited highly strict substrate specificity to glucose, while CGIHK2 could also catalyze other 11 monosaccharide substrates. This is the first report on the in vitro biosynthesis of glucose-6-phosphate by the hexokinases from Crassostrea gigas. The facile expression and purification procedures combined with different substrate specificities make CGIHK1 and CGIHK2 candidates for the biosynthesis of glucose-6-phosphate and other sugar-phosphates.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Taxonomic Distribution and Molecular Evolution of Mytilectins

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

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum has long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species. The final assembled and filtered Berghia genome is comparable to other high quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes), and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low. Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Invertebrate C1q Domain-Containing Proteins: Molecular Structure, Functional Properties and Biomedical Potential

C1q domain-containing proteins (C1qDC proteins) unexpectedly turned out to be widespread molecules among a variety of invertebrates, despite their lack of an integral complement system. Despite the wide distribution in the genomes of various invertebrates, data on the structure and properties of the isolated and characterized C1qDC proteins, which belong to the C1q/TNF superfamily, are sporadic, although they hold great practical potential for the creation of new biotechnologies. This review not only summarizes the current data on the properties of already-isolated or bioengineered C1qDC proteins but also projects further strategies for their study and biomedical application. It has been shown that further broad study of the carbohydrate specificity of the proteins can provide great opportunities, since for many of them only interactions with pathogen-associated molecular patterns (PAMPs) was evaluated and their antimicrobial, antiviral, and fungicidal activities were studied. However, data on the properties of C1qDC proteins, which researchers originally discovered as lectins and therefore studied their fine carbohydrate specificity and antitumor activity, intriguingly show the great potential of this family of proteins for the creation of targeted drug delivery systems, vaccines, and clinical assays for the differential diagnosis of cancer. The ability of invertebrate C1qDC proteins to recognize patterns of aberrant glycosylation of human cell surfaces and interact with mammalian immunoglobulins indicates the great biomedical potential of these molecules.

Gene presence/absence variation in Mytilus galloprovincialis and its implications in gene expression and adaptation

Presence/absence variation (PAV) is a well-known phenomenon in prokaryotes that was described for the first time in bivalves in 2020 in Mytilus galloprovincialis. The objective of the present study was to further our understanding of the PAV phenomenon in mussel biology. The distribution of PAV was studied in a mussel chromosome-level genome assembly, revealing a widespread distribution but with hotspots of dispensability. Special attention was given to the effect of PAV in gene expression, since dispensable genes were found to be inherently subject to distortions due to their sparse distribution among individuals. Furthermore, the high expression and strong tissue specificity of some dispensable genes, such as myticins, strongly supported their biological relevance. The significant differences in the repertoire of dispensable genes associated with two geographically distinct populations suggest that PAV is involved in local adaptation. Overall, the PAV phenomenon would provide a key selective advantage at the population level.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

Effect of heat and hypoxia stress on mitochondrion and energy metabolism in the gill of hard clam

Aquatic animals suffer from heat and hypoxia stress more frequently due to global climate change and other anthropogenic activities. Heat and hypoxia stress can significantly affect mitochondrial function and energy metabolism. Here, the response and adaptation characteristics of mitochondria and energy metabolism in the gill of the hard clam Mercenaria mercenaria under heat (35 °C), hypoxia (0.2 mg/L), and heat plus hypoxia stress (35 °C, 0.2 mg/L) after 48 h exposure were investigated. Mitochondrial membrane potentials were depolarized under environmental stress. Mitochondrial fusion, fission and mitophagy played a key role in maintain mitochondrion function. The AMPK subunits showed different expression under environmental stress. Acceleration of enzyme activities (phosphofructokinase, pyruvate kinase and lactic dehydrogenase) and accumulation of anaerobic metabolites in glycolysis and TCA cycle implied that the anaerobic metabolism might play a key role in providing energy. Accumulation of amino acids might help to increase tolerance under heat and heat combined hypoxia stress. In addition, urea cycle played a key role in amino acid metabolism to prevent ammonia/nitrogen toxicity. This study improved our understanding of the mitochondrial and energy metabolism responses of marine bivalves exposed to environmental stress.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria

The increasing concentration of CO₂ in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO₂ and reared under acidified (pH ~ 7.5, pCO₂ ~ 1200 ppm) or control (pH ~ 7.9, pCO₂ ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

Proteomic and Transcriptomic Responses Enable Clams to Correct the pH of Calcifying Fluids and Sustain Biomineralization in Acidified Environments

Seawater pH and carbonate saturation are predicted to decrease dramatically by the end of the century. This process, designated ocean acidification (OA), threatens economically and ecologically important marine calcifiers, including the northern quahog (Mercenaria mercenaria). While many studies have demonstrated the adverse impacts of OA on bivalves, much less is known about mechanisms of resilience and adaptive strategies. Here, we examined clam responses to OA by evaluating cellular (hemocyte activities) and molecular (high-throughput proteomics, RNASeq) changes in hemolymph and extrapallial fluid (EPF-the site of biomineralization located between the mantle and the shell) in M. mercenaria continuously exposed to acidified (pH ~7.3; pCO₂ ~2700 ppm) and normal conditions (pH ~8.1; pCO₂ ~600 ppm) for one year. The extracellular pH of EPF and hemolymph (~7.5) was significantly higher than that of the external acidified seawater (~7.3). Under OA conditions, granulocytes (a sub-population of hemocytes important for biomineralization) were able to increase intracellular pH (by 54% in EPF and 79% in hemolymph) and calcium content (by 56% in hemolymph). The increased pH of EPF and hemolymph from clams exposed to high pCO₂ was associated with the overexpression of genes (at both the mRNA and protein levels) related to biomineralization, acid-base balance, and calcium homeostasis, suggesting that clams can use corrective mechanisms to mitigate the negative impact of OA.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

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 Ca2+ content of hemocytes in the EPF were significantly higher than those in hemolymph, supporting the idea that hemocytes serve as a source of Ca2+ 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.