Molecular cloning, characterization and expression analysis of F-type lectin from pearl oyster Pinctada fucata

Molecular cloning, expression, and functional analysis of a putative lectin from the pearl oyster (Pinctada fucata, Gould 1850)

Marine lectins are a group of proteins that possess specific carbohydrate recognition and binding domains. They exhibit various activities, including antimicrobial, antitumor, antiviral, and immunomodulatory effects. In this study, a novel galectin-binding lectin gene named PFL-96 (GenBank: OQ561753.1) was cloned from Pinctada fucata. The PFL-96 gene has an open reading frame of 324 base pairs (bp) and encodes a protein comprising 107 amino acids. The protein has a molecular weight of 11.95 kDa and an isoelectric point of 9.27. It contains an N-terminal signal peptide and a galactose-binding lectin domain. The sequence identity to lectin proteins from fish, echinoderms, coelenterates, and shellfish ranges from 31.90 to 40.00 %. In the phylogenetic analysis, it was found that the PFL-96 protein is closely related to the lectin from Pteria penguin. The PFL-96 recombinant protein exhibited coagulation activity on 2 % rabbit red blood cells at a concentration of ≥8 μg/mL. Additionally, it showed significant hemolytic activity at a concentration of ≥32 μg/mL. The PFL-96 recombinant protein exhibited significant antibacterial activity against Bacillus subtilis, Staphylococcus aureus, Candida albicans, and Vibrio alginolyticus, with minimum inhibitory concentrations (MIC) of 4, 8, 16, and 16 μg/mL, respectively. The minimum bactericidal concentrations (MBC) were determined to be 8, 16, 32, and 32 μg/mL, respectively. Furthermore, the PFL-96 recombinant protein exhibited inhibitory effects on the proliferation of Hela tumor cells, HepG2 tumor cells, and C666-1 tumor cells, with IC50 values of 7.962, 8.007, and 9.502 μg/mL, respectively. These findings suggest that the recombinant protein PFL-96 exhibits significant bioactivity in vitro, contributing to a better understanding of the active compounds found in P. fucata. The present study establishes a fundamental basis for further investigation into the mechanism of action and structural optimization of the recombinant protein PFL-96. The aim is to develop potential candidates for antibacterial and anti-tumor agents.

A first insight into haemocytes of Pinctada imbricata radiata: A morpho-functional characterization

The pearl oyster Pinctada imbricata radiata (Leach, 1814), from the Pacific Ocean, was one of the first species to reach via Suez the Mediterranean, colonizing the eastern basin and recently spreading to the western. The species showed to be able to adapt to a wide range of climatic, hydrological, and ecological conditions. Since 2000 it reached the Strait of Messina, where is now infesting the transitional waters of the oriented natural reserve "Laguna di Capo Peloro." Due to such resistance and adaptation ability, various assays were performed. Haemocyte morpho-functional aspects were evaluated in haemolymph samples fixed with 1% and 2% glutaraldehyde for optical and electron microscopy (TEM). The following assays were carried out: cell characterization using several dyes, detection of intra- and extracellular lipids, the capability of phagocytosis using the yeast Saccharomyces cerevisiae and to produce superoxide anion (O^2- ). Detection of several enzymes, such as acid and alkaline phosphatase, arylsulfatase, chloro-acetylesterase and β-glucuronidase was also assessed. Cell count was demonstrated to be abundant with a mean of 8.263 × 106 mm²  ± 0.935 × 106 (SD). Two main cell populations were noticed: granulocytes and hyalocytes, both competent for phagocytosis, to produce O^2- , and characterized by lipids. Based on the granule analysis, enzymatic activity was also demonstrated. The observations under TEM confirmed all the results obtained. This study supports the hypothesis that P. imbricata radiata can be usefully employed as a model organism in environmental biomonitoring. Moreover, since the species represent potential threats to native species and ecosystems, further insights into its biological adaptations in invaded ecosystems are recommended.

Effect of cell-permeable grouper Manganese Superoxide Dismutase on environmental stress in fish

Nitrite levels are generally high in high-density aquaculture. Nitrite is a potential stress-inducing factor and can cause oxidative stress because excessive reactive oxygen species (ROS) formation through nitrite induction cannot be scavenged by the endogenous antioxidant system, thus leading to cell damage or death. Manganese Superoxide Dismutase (MnSOD) is a highly efficient endogenous ROS scavenger that quenches mitochondrial ROS and protective against oxidative stress. To enhance the efficiency of MnSOD in removing ROS and reducing oxidative caused by nitrite, in this study, we cloned grouper MnSOD (gMnSOD) fused with a cell-penetrating peptide, TAT, to construct a TAT-gMnSOD fusion protein and assessed its potential to eliminate excess ROS induced by high nitrite concentrations and enhance the resistance of zebrafish to environmental stressors. Our results revealed that TAT-gMnSOD penetrated the grouper fin (GF-1) cells, scavenged nitrite-induced intracellular ROS, and enhanced cell viability on NaNO₂ treatment. Furthermore, pretreatment of zebrafish with TAT-gMnSOD fusion protein reduced the MDA content and increased the survival rate. In addition, the TAT-gMnSOD fusion protein reduced 2-phenoxyethanol toxicity and attenuated excessive anesthesia among zebrafish. In conlusion, our cell-permeable TAT-gMnSOD fusion protein effectively counters oxidative stress, prevents environmental stress-induced damage, and increases aquaculture benefits.

Characterization and expression of galectin-3 in grass carp (Ctenopharyngodon idella)

Galectins are members of evolutionary conserved lectin family and play important roles in the innate and adaptive immunity of both vertebrates and invertebrates. Galectin-3 is the only chimera galectin with one C-terminal carbohydrate recognition domain (CRD) connected to the N-terminal end. Here, a galectin-3 (named CiGal3) from grass carp was identified and characterized, which encodes polypeptides 362 amino acids with a predicted molecular mass of 36.45 kDa and theoretical isoelectric point of 4.91. The sugar binding motifs involved in carbohydrate binding activity (H-N-R, V-N and W--E-R) were detected in CRD. In comparison to other species, CiGal3 showed the highest similarity and identity to Cyprinus carpio (95.3% sequence similarity and 92.5% sequence identity). The subcellular localization of CiGal3 was distributed in the cytoplasm and nucleus of transfected cells. The CiGal3 transcripts were ubiquitously expressed in all checked tissues and highly expressed in immune tissues. In addition, the expression of CiGal3 in liver and spleen was induced post grass carp reovirus (GCRV), lipopolysaccharide (LPS), and polyinosinic:polycytidylic acid (poly I:C) challenge. These results suggest that CiGal3 plays a vital role in the immune system.

Glycan Binding Profiling of Jacalin-Related Lectins from the Pteria Penguin Pearl Shell

We determined the primary structures of jacalin-related lectins termed PPL3s (PPL3A, 3B, and 3C, which are dimers consisting of sequence variants α + α, α + β, β + β, respectively) and PPL4, which is heterodimer consisting of α + β subunits, isolated from mantle secretory fluid of Pteria penguin (Mabe) pearl shell. Their carbohydrate-binding properties were analyzed, in addition to that of PPL2A, which was previously reported as a matrix protein. PPL3s and PPL4 shared only 35–50% homology to PPL2A, respectively; they exhibited significantly different carbohydrate-binding specificities based on the multiple glycan binding profiling data sets from frontal affinity chromatography analysis. The carbohydrate-binding specificity of PPL3s was similar to that of PPL2A, except only for Man₃Fuc₁Xyl₁GlcNAc₂ oligosaccharide, while PPL4 showed different carbohydrate-binding specificity compared with PPL2A and PPL3s. PPL2A and PPL3s mainly recognize agalactosylated- and galactosylated-type glycans. On the other hand, PPL4 binds to high-mannose-and hybrid-type N-linked glycans but not agalactosylated- and galactosylated-type glycans.

A manganese superoxide dismutase (MnSOD) from red lip mullet, Liza haematocheila: Evaluation of molecular structure, immune response, and antioxidant function

Manganese superoxide dismutase (MnSOD) is a nuclear-encoded antioxidant metalloenzyme. The main function of this enzyme is to dismutase the toxic superoxide anion (O₂-) into less toxic hydrogen peroxide (H₂O₂) and oxygen (O₂). Structural analysis of mullet MnSOD (MuMnSOD) was performed using different bioinformatics tools. Pairwise alignment revealed that the protein sequence matched to that derived from Larimichthys crocea with a 95.2% sequence identity. Phylogenetic tree analysis showed that the MuMnSOD was included in the category of teleosts. Multiple sequence alignment showed that a SOD Fe-N domain, SOD Fe-C domain, and Mn/Fe SOD signature were highly conserved among the other examined MnSOD orthologs. Quantitative real-time PCR showed that the highest MuMnSOD mRNA expression level was in blood cells. The highest expression level of MuMnSOD was observed in response to treatment with both Lactococcus garvieae and lipopolysaccharide (LPS) at 6 h post treatment in the head kidney and blood. Potential ROS-scavenging ability of the purified recombinant protein (rMuMnSOD) was examined by the xanthine oxidase assay (XOD assay). The optimum temperature and pH for XOD activity were found to be 25 °C and pH 7, respectively. Relative XOD activity was significantly increased with the dose of rMuMnSOD, revealing its dose dependency. Activity of rMuMnSOD was inhibited by potassium cyanide (KCN) and N-N'-diethyl-dithiocarbamate (DDC). Moreover, expression of MuMnSOD resulted in considerable growth retardation of both gram-positive and gram-negative bacteria. Results of the current study suggest that MuMnSOD acts as an antioxidant enzyme and participates in the immune response in mullet.

Molecular cloning and characteristics analysis of Pmtgfbr1 from Pinctada fucata martensii

•

This study obtains the full length of Pmtgfbr1 of the pearl oyster P. fucata martensii.

•

Pmtgfbr1 possesses the conserved domain of Tgfbr1.

•

Pmtgfbr1 holds negatively effect on the growth of P. fucata martensii.

Pinctada fucata martensii is cultured for pearl production. Growth improvement has received considerable research interest. Transforming growth factor β type Ⅰ receptor (TβR-I), which is involved in signals transmission of transforming growth factor beta (TGF-β), participates in cell proliferation and growth. In this study, we characterized a Tgfbr1 gene which encoded TβR-I from P. fucata martensii (Pmtgfbr1). Pmtgfbr1 cDNA contains an open reading frame of 1569 bp and encodes a polypeptide of 522 amino acids (aa). Pmtgfbr1 possesses a typical TβR-I structure (extracellular receptor ligand domain, transmembrane domain, and cytoplasmic tyrosine kinase catalytic domain). Pmtgfbr1 is expressed in all the studied tissues and exhibited the highest expression level in the adductor muscle. Moreover, Pmtgfbr1 exhibited the lower expression level in the larger group (L) than that in the smaller group (S) and is negatively correlated with growth traits (P < 0.01). Our results indicated that Pmtgfbr1 is a candidate functional gene associated with growth traits.

Molecular and functional analysis of PmCHST1b in nacre formation of Pinctada fucata martensii

Keratan sulfate possesses considerable amounts of negatively charged sulfonic acid groups and participates in biomineralization. In the present study, we investigated characteristics and functions of a CHST1 gene identified from the pearl oyster Pinctada fucata martensii (PmCHST1b) which participated in the synthesis of keratan sulfate. PmCHST1b amino acid sequence carried a typical sulfotransferase-3 domain (sulfotransfer-3 domain) and belonged to membrane-associated sulfotransferases. Homologous analysis of CHST1 from different species showed the conserved motif (5' PSB motif and 3' PB motif) which interacted with 3'-phosphoadenosine-5'-phosphosulfate (PAPS). Structure analysis of sulfotransferase domain indicted that PmCHST1b showed the conserved catalytic structure character and the relationships presented in the phylogenetic tree conformed to that of traditional taxonomy. Expression pattern of PmCHST1b in different tissues and development stages showed that PmCHST1b widely expressed in all the detected tissues and development stages and showed the highest expression level in the central zone of mantle (MC). PmCHST1b expressed highly in the trochophore, D-stage larvae and spat which corresponded to prodissoconch and dissoconch shell formation, respectively. RNA interference (RNAi) successfully inhibited expression level of PmCHST1b in MC (P<0.05), and sulfate polymer content in the extrapallial fluid significantly reduced (P<0.05). Crystallization of shell nacre became irregular. Results above indicated that PmCHST1b may affect nacre formation by participating in synthesis of keratan sulfate in extrapallial fluid. This study provided fundamental materials for further research on the role of sulfotransferases and keratan sulfate in nacre formation.

A novel fucolectin from Apostichopus japonicus with broad PAMP recognition pattern

F-type lectin (also known as fucolectin) is a newly identified family of fucose binding lectins with the sequence characters of a fucose binding motif and a unique lectin fold (the "F-type" fold). In the present study, a fucolectin was identified from sea cucumber Apostichopus japonicus (designated AjFL-1). The open reading frame (ORF) of AjFL-1 was of 546 bp, encoding a polypeptide of 181 amino acids with a predicted molecular mass of about 20 kDa. The deduced amino acid sequence of AjFL-1 shared 30%-40% similarity with the fucolectins from other animals. There were a typical F-type lectin domain (FLD) (residues 39-180) and a signal peptide (residues 1-24) in AjFL-1. The mRNA transcript of AjFL-1 could be detected by qRT-PCR in various tissues, such as intestinum, coelomocytes, respiratory tree, tentacle, and body wall, while undetectable in the gonads and longitudinal muscle. The mRNA expression level of AjFL-1 in coelomocytes was significantly up-regulated (47.06-fold to that in control group, p < 0.05) at 12 h after Vibrio splendidus challenge. Immunofluorescence assay showed that AjFL-1 protein was mainly distributed on the membrane, while few in cytoplasm of coelomocytes in sea cucumber. The recombinant AjFL-1 (rAjFL-1) could bind lipopolysaccharide (LPS), peptidoglycan (PGN), mannan (MAN) and fucose (FUC), and exhibited a broader binding activities towards Gram-negative bacterium Escherichia coli, Gram-positive bacterium Micrococcus luteus, as well fungus Pichia pastoris. In addition, rAjFL-1 could strongly promote the agglutination of fungus P. pastoris. These results indicated that AjFL-1 was a novel member of fucose-binding lectin family, which functioned as a pattern recognition receptor with broad spectrum of microbial recognition, and involved in innate immune response of sea cucumber.

Identification of EGFR in pearl oyster (Pinctada fucata martensii) and correlation analysis of its expression and growth traits

Marine pearl production is directly influenced by the growth speed of Pinctada fucata martensii. However, the slow growth rate of this organism remains the main challenge in aquaculture production. Epidermal growth factor receptor (EGFR), an important receptor of tyrosine kinases in animals, plays versatile functions in development, growth and tissue regeneration. In this study, we described the characteristic and function of an EGFR gene identified from P. f. martensii (PmEGFR). PmEGFR possesses a typical EGFR structure and is expressed in all studied tissues, with the highest expression level in adductor muscle. PmEGFR expression level is significantly higher in the fast-growing group than that in the slow-growing one. Correlation analysis represents that shell height and shell weight show positive correlation with PmEGFR expression (p < 0.05), and total weight and tissue weight exhibit positive correlation with it (p < 0.01). This study indicates that PmEGFR is a valuable functional gene associated with growth traits.

Structures of jacalin-related lectin PPL3 regulating pearl shell biomineralization

The nacreous layer of pearl oysters is one of the major biominerals of commercial and industrial interest. Jacalin-related lectins, including PPL3 isoforms, are known to regulate biomineralization of the Pteria penguin pearl shell, although the molecular mechanisms are largely unknown. The PPL3 crystal structures were determined partly by utilizing microgravity environments for 3 isoforms, namely, PPL3A, PPL3B, and PPL3C. The structures revealed a tail-to-tail dimer structure established by forming a unique inter-subunit disulfide bond at C-termini. The N-terminal residues were found in pyroglutamate form, and this was partly explained by the post-translational modification of PPL3 isoforms implied from the discrepancy between amino acid and gene sequences. The complex structures with trehalose and isomaltose indicated that the novel specificity originated from the unique α-helix of PPL3 isoforms. Docking simulations of PPL3B to various calcite crystal faces suggested the edge of a β-sheet and the carbohydrate-binding site rich in charged residues were the interface to the biomineral, and implied that the isoforms differed in calcite interactions.

F-Type Lectins: A Highly Diversified Family of Fucose-Binding Proteins with a Unique Sequence Motif and Structural Fold, Involved in Self/Non-Self-Recognition

The F-type lectin (FTL) family is one of the most recent to be identified and structurally characterized. Members of the FTL family are characterized by a fucose recognition domain [F-type lectin domain (FTLD)] that displays a novel jellyroll fold ("F-type" fold) and unique carbohydrate- and calcium-binding sequence motifs. This novel lectin family comprises widely distributed proteins exhibiting single, double, or greater multiples of the FTLD, either tandemly arrayed or combined with other structurally and functionally distinct domains, yielding lectin subunits of pleiotropic properties even within a single species. Furthermore, the extraordinary variability of FTL sequences (isoforms) that are expressed in a single individual has revealed genetic mechanisms of diversification in ligand recognition that are unique to FTLs. Functions of FTLs in self/non-self-recognition include innate immunity, fertilization, microbial adhesion, and pathogenesis, among others. In addition, although the F-type fold is distinctive for FTLs, a structure-based search revealed apparently unrelated proteins with minor sequence similarity to FTLs that displayed the FTLD fold. In general, the phylogenetic analysis of FTLD sequences from viruses to mammals reveals clades that are consistent with the currently accepted taxonomy of extant species. However, the surprisingly discontinuous distribution of FTLDs within each taxonomic category suggests not only an extensive structural/functional diversification of the FTLs along evolutionary lineages but also that this intriguing lectin family has been subject to frequent gene duplication, secondary loss, lateral transfer, and functional co-option.

Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii

Immunological rejection of the pearl oysters following nucleus implantation is a major issue limiting the successful rate of cultured pearls. To date, the molecular mechanism of immune tolerance during pearl formation in the pearl oysters is still largely unknown. Through the RNA sequencing platform and comparative transcriptomic analysis, we investigated the chronic gene expression changes at seven time points (0, 5, 10, 15, 20, 30, 60 days post implantation or dpi) over a period of 60 days following nucleus implantation in the pearl oyster Pinctada martensii. A total of 81,390 unique transcripts (or unigenes) with a combined length of 96.8 million bp and a N50 value of 2227 bp were obtained. When compared with sequences in the nr, nt, Swiss-Prot, KEGG, COG and GO databases, 36,380 unigenes can find homologous genes. Pairwise comparison of gene expression among all the samples showed that the largest number (or 6846) of differentially expressed genes was observed at 10 dpi. The number then decreased to below 5000 at 15, 20 and 30 dpi and increased again to 6679 at 60 dpi. PCA analysis further showed that the seven time points can be roughly divided into four groups. Comparative transcriptomic analysis between the four groups identified a variety of genes showing differential expression at different time points, including many immune-related genes such as those encoding for toll-like receptor, lectin, scavenger receptor, and peroxidase. In addition, GO and KEGG enrichment analysis revealed that these differentially expressed genes were mainly associated with metabolism, ribosome function, immune response, signaling transduction, and cytoskeleton organization. Notably, two KEGG pathways, namely "cell adhesion molecules" and "primary immunodeficiency" were significantly enriched during the whole process. This finding indicates that genes in these pathways are likely to play critical roles in the immune tolerance of the pearl oysters. To conclude, the data obtained contribute to a better understanding of the molecular mechanisms of nucleus implantation induced immune response in the pearl oysters, and will facilitate the development of effective measures to improve the performance of pearl culture.

The purplish bifurcate mussel Mytilisepta virgata gene expression atlas reveals a remarkable tissue functional specialization

BACKGROUND

Mytilisepta virgata is a marine mussel commonly found along the coasts of Japan. Although this species has been the subject of occasional studies concerning its ecological role, growth and reproduction, it has been so far almost completely neglected from a genetic and molecular point of view. In the present study we present a high quality de novo assembled transcriptome of the Japanese purplish mussel, which represents the first publicly available collection of expressed sequences for this species.

RESULTS

The assembled transcriptome comprises almost 50,000 contigs, with a N50 statistics of ~1 kilobase and a high estimated completeness based on the rate of BUSCOs identified, standing as one of the most exhaustive sequence resources available for mytiloid bivalves to date. Overall this data, accompanied by gene expression profiles from gills, digestive gland, mantle rim, foot and posterior adductor muscle, presents an accurate snapshot of the great functional specialization of these five tissues in adult mussels.

CONCLUSIONS

We highlight that one of the most striking features of the M. virgata transcriptome is the high abundance and diversification of lectin-like transcripts, which pertain to different gene families and appear to be expressed in particular in the digestive gland and in the gills. Therefore, these two tissues might be selected as preferential targets for the isolation of molecules with interesting carbohydrate-binding properties. In addition, by molecular phylogenomics, we provide solid evidence in support of the classification of M. virgata within the Brachidontinae subfamily. This result is in agreement with the previously proposed hypothesis that the morphological features traditionally used to group Mytilisepta spp. and Septifer spp. within the same clade are inappropriate due to homoplasy.

A galectin from Hyriopsis cumingii involved in the innate immune response against to pathogenic microorganism and its expression profiling during pearl sac formation

Hyriopsis cumingii is the most important freshwater pearl mussel cultured in China. The operation for implantation is one necessary technical step for pearl culture. However, implantation-induced trauma results in a series of immune responses and can enable the invasion of pathogenic microbes. Lectin proteins are found widely in nature and play important roles in innate immunity. Galectins are members of the lectin superfamily and are characterized by one or several carbohydrate recognition domains (CRDs) that produce multiple sugar binding sites on the protein. Here we cloned and characterized the H. cumingii galectin gene HcGal1, which encodes a 312 amino acid galectin protein. The HcGal1 transcript was detected in all tested H. cumingii tissues and showed higher expression specifically in immune tissues. The significant upregulation of HcGal1 expression was observed after challenging the mussel with lipopolysaccharide or Gram-negative and Gram-positive bacteria. After implantation, significant downregulation of the HcGal1 transcript was noted in the mantle, hemocytes, and pearl sac in the acute-stress stage (0-24 h) and the stage of wound healing and pearl-sac formation (24 h-7 d). In addition, significant upregulation of HcGal1 expression was observed in the liver in the stage of wound healing and pearl-sac formation. In the pearl-secretion stage (7-35 d), the HcGal1 transcript levels returned to normal in all tested tissues. We also show that recombinantly expressed and purified HcGal1 can agglutinate some Gram-negative and Gram-positive bacteria. In addition, in vivo experiments showed that the recombinant protein HcGal1 could promote phagocytosis by hemocytes. Our data suggest that HcGal1 plays a role in innate immune responses involved in pathogen recognition and wound healing.

А new Gal/GalNAc-specific lectin from the mussel Mytilus trossulus: Structure, tissue specificity, antimicrobial and antifungal activity

In the present study, a new Gal/GalNAc specific lectin from the mussel Mytilus trossulus (designated as MTL) was identified, and its expression levels, both in tissues and toward pathogen stimulation, were then characterized. The MTL primary structure was determined via cDNA sequencing. Deduced sequence of 150 amino acid residues showed 89% similarity to lectins from the mussels Crenomytilus grayanus and Mytilus galloprovincialis that were the first members of a new family of zoolectins. The results indicated that the MTL might be involved in immune response toward pathogen infection, and it might perform different recognition specificity toward bacteria or fungi.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.