Multiple ferritin subunit genes of the Pacific oyster Crassostrea gigas and their distinct expression patterns during early development

Integrative analysis of single-nucleus RNA-seq and bulk RNA-seq reveals germline cells development dynamics and niches in the Pacific oyster gonad

Gametogenesis drives the maturation of germ cell precursors into functional gametes, facilitated by interactions with the niche environment. However, the molecular mechanisms, especially in invertebrates, remain incompletely understood. In this study, the gonadal microenvironment and gametogenic processes in the Pacific oyster, a model for diffuse gonadal organization and periodic gametogenesis, are investigated. We combine single-nucleus RNA-seq and bulk RNA-seq to analyze gonadal microenvironments in oysters. Twenty-three male and nineteen female gonadal cell clusters are identified, revealing four male and three female germ cell types, alongside follicular cells in females and Sertoli/Leydig cells in males. The NOTCH and BMP (bone morphogenetic protein) signaling pathways play a significant role in the male germline niche, suggesting similarities with mammalian germ cell microenvironment. This study offers valuable insights into germ cell developmental transitions and microenvironmental characteristics.

The gastropod Lottia peitaihoensis as a model to study the body patterning of trochophore larvae

The body patterning of trochophore larvae is important for understanding spiralian evolution and the origin of the bilateral body plan. However, considerable variations are observed among spiralian lineages, which have adopted varied strategies to develop trochophore larvae or even omit a trochophore stage. Some spiralians, such as patellogastropod mollusks, are suggested to exhibit ancestral traits by producing equal-cleaving fertilized eggs and possessing "typical" trochophore larvae. In recent years, we developed a potential model system using the patellogastropod Lottia peitaihoensis (= Lottia goshimai). Here, we introduce how the species were selected and establish sources and techniques, including gene knockdown, ectopic gene expression, and genome editing. Investigations on this species reveal essential aspects of trochophore body patterning, including organizer signaling, molecular and cellular processes connecting the various developmental functions of the organizer, the specification and behaviors of the endomesoderm and ectomesoderm, and the characteristic dorsoventral decoupling of Hox expression. These findings enrich the knowledge of trochophore body patterning and have important implications regarding the evolution of spiralians as well as bilateral body plans.

Search for new biomarkers of tolerance to Perkinsus olseni parasite infection in Ruditapes decussatus clams

The grooved carpet shell (Ruditapes decussatus) is a clam species with high economic and social importance in several European and Mediterranean countries. Production of this species suffered a decline caused by biotic (parasite infection) and abiotic factors (environmental factors, stress, poor management methods and intensive culture of the introduced species Ruditapes philippinarum). The protozoan parasite Perkinsus olseni is also responsible for the decline of production, being nowadays one of the major issues for clam culture. Molecular biomarkers that might represent tolerance of R. decussatus to P. olseni have already been uncovered, shedding light in a possible production improvement by selecting those clams with a strongest immune response. In the present study, new tolerance biomarkers to P. olseni infection in R. decussatus were identified. The haemolymph proteomic profiles of naturally non/low-infected (tolerant) and highly-infected (susceptible) clams by the parasite across several heavy affected areas of Europe were characterized through a shotgun proteomics approach. Also, the mechanisms that might be involved in the responses against the disease in chronic infections were explored. Proteins related to energy restoration and balance, metabolic regulation, energy accumulation, ROS production, lysosomal activity, amino acid synthesis, proteolytic activity, iron regulation, iron withholding, and immune response modulation were significantly regulated in susceptible clams. In the tolerant group, proteins related to phagocytosis regulation, control of cell growth and proliferation, gonadal maturation, regulation of apoptosis, growth modulation, response to oxidative stress, iron regulation, shell development and metabolic regulation were significantly expressed. In summary, the protein expression profile of tolerant individuals suggests that an efficient pathogen elimination mechanism coupled to a better metabolic regulation leads to a tolerance to the parasite infection by limiting the spread through the tissues.

Conservative and Atypical Ferritins of Sponges

Ferritins comprise a conservative family of proteins found in all species and play an essential role in resistance to redox stress, immune response, and cell differentiation. Sponges (Porifera) are the oldest Metazoa that show unique plasticity and regenerative potential. Here, we characterize the ferritins of two cold-water sponges using proteomics, spectral microscopy, and bioinformatic analysis. The recently duplicated conservative HdF1a/b and atypical HdF2 genes were found in the Halisarca dujardini genome. Multiple related transcripts of HpF1 were identified in the Halichondria panicea transcriptome. Expression of HdF1a/b was much higher than that of HdF2 in all annual seasons and regulated differently during the sponge dissociation/reaggregation. The presence of the MRE and HRE motifs in the HdF1 and HdF2 promotor regions and the IRE motif in mRNAs of HdF1 and HpF indicates that sponge ferritins expression depends on the cellular iron and oxygen levels. The gel electrophoresis combined with specific staining and mass spectrometry confirmed the presence of ferric ions and ferritins in multi-subunit complexes. The 3D modeling predicts the iron-binding capacity of HdF1 and HpF1 at the ferroxidase center and the absence of iron-binding in atypical HdF2. Interestingly, atypical ferritins lacking iron-binding capacity were found in genomes of many invertebrate species. Their function deserves further research.

Expression of four new ferritins from grooved carpet shell clam Ruditapes decussatus challenged with Perkinsus olseni and metals (Cd, Cu and Zn)

Iron has a fundamental role in life and in its biochemical reactions but, when in excess, it can promote the formation of free radicals which can lead to cell death. Therefore, managing the levels of iron is essential to regulate the production of oxidative stress related to iron, and ferritins are one of the main protein families involved in this process. Ferritins are ≈480 kDa multimeric proteins composed by 24 subunits, each with 19-26 kDa, which can accumulate up to 4500 iron atoms. Besides their role in managing iron bioavailability, they have also developed a role in organism immunity and defence present throughout evolution. In this work, we identified and characterized, for the first time, four different ferritin subunits in the clam Ruditapes decussatus, a bivalve commercially and ecologically important along the south Atlantic coast and in the Mediterranean basin, which is a major target of the parasitic protozoa Perkinsus olseni, considered one of the main causes of high levels of clam mortality. Following phylogenetic annotation, the four ferritins subunits identified were subdivided into two cytosolic and two secreted forms. All four subunits maintain the canonical ferritin structure with four main helices α (A-D) and a small helix (E), but the secreted ferritins present an additional helix in their N-terminal region (F), located after the signal peptide and with possible antimicrobial properties. Additionally, we identified in ferritin 4 an extra helix α (G) located between helices B and C. These alpha helix domains revealed high degree of similarity with antimicrobial peptides associated with antibacterial and antifungal activities. Analysis of the expression of these subunits showed that ferritins 1 and 2 are ubiquitously expressed while ferritins 3 and 4 are present mainly in visceral mass. Ferritin 1 lacked a putative functional iron response element (IRE) and appeared to be under a tight regulation. Ferritins 2 and 3 showed a strong response to infection by parasite Perkinsus olseni in contrast to ferritin 4, whose main response was related to exposure to a combination of metals. The synergistic effect between metals and infection promoted a general upregulation of the four ferritins. In conclusion, our results suggest that ferritins, besides their function in iron and metals detoxification, may play a determinant role in clam immune response.

Identification of multiple ferritin genes in Macrobrachium nipponense and their involvement in redox homeostasis and innate immunity

Based on the transcriptome database, we screened out four ferritin subunit genes (MnFer2-5) from the oriental river prawn Macrobrachium nipponense, which encode two non-secretory and two secretory peptides. MnFer2 and 4 possess a strictly conserved ferroxidase site, and MnFer3 has a non-typical ferroxidase site. MnFer5 seems to be a number of ferritin families, which has a distinct dinuclear metal binding motif, but lacks an iron ion channel, a ferroxidase site and a nucleation site. Diverse tissue-specific transcriptions of the four genes indicate their functional diversity in the prawn. Among them, MnFer2 is mainly expressed in hepatopancreas and intestines, MnFer3 and 4 are predominantly expressed in gills, and MnFer5 is widely expressed in various tissues with high presence in intestines, hepatopancreas and haemocytes. The transcription of all the four MnFer genes can be strongly induced by doxorubicin, indicating the involvement of these ferritin subunits in protection from oxidative stress. Upon Aeromonas hydrophila infection, only MnFer5 is persistently up-regulated, while other subunits including MnFer2-4 are down-regulated during the early stage, followed by recovery and even a slight increase at 48 h post bacterial challenge. Moreover, the iron binding capacity of recombinant MnFer2 is also demonstrated in vitro. The E. coli expressing MnFer2 displays increased resistance to hydrogen peroxidase cytotoxicity. These results suggest a protective role of ferritins from M. nipponense in iron homeostasis, redox biology and antibacterial immunity and shed light on the molecule evolution of crustacean ferritin subunits.

Cloning and molecular characterization of two ferritins from red abalone Haliotis rufescens and their expressions in response to bacterial challenge at juvenile and adult life stages

Ferritins are ubiquitous proteins with a pivotal role in iron storage and homeostasis, and in host defense responses during infection by pathogens in several organisms, including mollusks. In this study, we characterized two ferritin homologues in the red abalone Haliotis rufescens, a species of economic importance for Chile, USA and Mexico. Two ferritin subunits (Hrfer1 and Hrfer2) were cloned. Hrfer1 cDNA is an 807 bp clone containing a 516 bp open reading frame (ORF) that corresponds to a novel ferritin subunit in H. rufescens. Hrfer2 cDNA is an 868 bp clone containing a 516 bp ORF that corresponds to a previously reported ferritin subunit, but in this study 5'- and 3'-UTR sequences were additionally found. We detected a putative Iron Responsive Element (IRE) in the 5'-UTR sequence, suggesting a posttranscriptional regulation of Hrfer2 translation by iron. The deduced protein sequences of both cDNAs possessed the motifs and domains required in functional ferritin subunits. Expression patterns of both ferritins in different tissues, during different developmental stages, and in response to bacterial (Vibrio splendidus) exposure were examined. Both Hrfer1 and Hrfer2 are most expressed in digestive gland and gonad. Hrfer1 mRNA levels increased about 34-fold along with larval developmental process, attaining the highest level in the creeping post-larvae. Exogenous feeding is initiated at the creeping larva stage; thus, the increase of Hrfer1 may suggest and immunity-related role upon exposure to bacteria. Highest Hrfer2 expression levels were detected at trochophore stage; which may be related with early shell formation. Upon challenge with, the bacteria an early mild induction of Hrfer2 (2 h post-challenge), followed by a stronger induction of Hrfer1 at 15 h post-challenge, was observed in haemocytes from adult abalones. While maximal upregulation of both genes in the whole individual occurred at 24 h post-challenge, in juveniles. A significant increase in ferritin protein levels from 6 h to 24 h post-challenge was also detected. Our results suggest an involvement of Hrfer1 and Hrfer2, and of ferritin proteins in the immune response of H. rufescens to bacterial infection.

A novel caspase-associated recruitment domain (CARD) containing protein (CgCARDCP-1) involved in LPS recognition and NF-κB activation in oyster (Crassostrea gigas)

Caspase-associated recruitment domain (CARD) containing proteins play critical roles in molecular interaction and regulation of various signaling pathways, such as the activation of caspase and NF-κB singling pathway in the process of apoptosis or inflammation. In the present study, a novel CARD containing protein (designed CgCARDCP-1) was identified and characterized from oyster Crassostrea gigas. Molecular feature analysis revealed that, the open reading frame (ORF) of CgCARDCP-1 gene was 759 bp encoding a polypeptide of 253 amino acids with a conserved N-terminal CARD domain and two transcriptional coactivator p15 (PC4) domains in C-terminus. Homologous alignment showed that the amino acid sequence of CgCARDCP-1 shared 30%-46% identity with that of caspase-2. By RT-PCR detection, the mRNA transcripts of CgCARDCP-1 were found to be widely distributed in various tissues of oyster with the highest expression level in hemocytes and mantle. And CgCARDCP-1 protein was mostly distributed in the cytoplasm of oyster hemocytes as shown by immunohistochemistry. Moreover, the CgCARDCP-1 mRNA expression level in hemocytes was significantly up-regulated after lipopolysaccharide (LPS) and Vibrio splendidus stimulations. The recombinant CgCARDCP-1 displayed strong binding activity with LPS in vitro. In addition, after transfected into the HEK-293T cell with luciferase reporter system, CgCARDCP-1 could significantly promote the NF-κB activation (1.29-fold, p < 0.05) compared to that in the control group. These results collectively demonstrated that the CgCARDCP-1 might serve as a recognition molecule for LPS and a regulator of NF-κB activation in the immune response of oyster.

Description of a Second Ferritin Light Chain Homologue From the Yellow Fever Mosquito (Diptera: Culicidae)

Ferritin is required for iron storage in vertebrates and for iron transport and storage in invertebrates, specifically insects. Classical ferritins consist of 24 subunits configured as a polyhedron wherein iron is held. The 24 subunits include light and heavy chains, each with specific functions. Several homologues of the light and heavy chains have been sequenced and studied in insects. In addition to iron transport and storage, ferritin has a role in dietary iron absorption, and functions as a protective agent preventing iron overload, decreasing oxidative stress, and reducing infection in these animals. The expression profile and regulation of a second ferritin light chain homologue (LCH2) in Aedes aegypti [Linnaeus (Diptera: Culicidae), yellow fever mosquito] was characterized in cells, animal developmental stages, and tissues post bloodmeal (PBM) by real-time PCR and immunoblot. Two previously studied ferritin subunits from Ae. aegypti, HCH and LCH1, along with LCH2 were immunoprecipitated and analyzed by mass spectrometry. The three Ae. aegypti ferritin subunits, HCH, LCH1, and LCH2, have different expression profiles and regulation with iron exposure, developmental stage, and tissue response PBM. Ae. aegypti expresses multiple and unique ferritin light chain subunits. Ae. aegypti, the vector for Zika, Dengue, and yellow fever, requires iron for oogenesis that is transported and stored in ferritin; this vector expresses a second light chain ferritin subunit homologue unlike any other species in which ferritin has been studied to date.

Effects of phenanthrene on early development of the Pacific oyster Crassostrea gigas (Thunberg, 1789)

Phenanthnere (PHE) is a polycyclic aromatic hydrocarbon continuously discarded in the marine environment and bioavailable to many aquatic species. Although studies about PHE toxicity have been documented for adult oysters, the effects on early developmental stages are poorly characterized in bivalves. In this study, the effects of PHE (0.02 and 2.0μg.L^-1) were evaluated on the embryogenesis and larval development of Crassostrea gigas. Toxicity bioassays, growth and deformities assessment, analysis of shell calcium abundance and transcript levels of genes related to xenobiotic biotransformation (CYP2AU2, CYP30C1), immune system (Cg-Tal) and tissue growth and shell formation (Ferritin, Insulin-like, Cg-Try, Calmodulin and Nacrein) were assayed in D-shape larvae after 24h of PHE exposure. At the highest concentration (2.0μg.L^-1), PHE decreased the frequency of normal development (19.7±2.9%) and shell size (53.5±2.8mm). Developmental deformities were mostly related to abnormal mantle and shell formation. Lower calcium levels in oyster shells exposed to PHE 2.0μg.L^-1 were observed, suggesting effects on shell structure. At this same PHE concentration, CYP30C1, Cg-Tal, Cg-Tyr, Calmodulin were upregulated and CYP2AU2, Ferritin, Nacrein, and Insulin-Like were downregulated compared to control larvae. At the lowest PHE concentration (0.02μg.L^-1), it was observed a minor decrease in normal larval development (89,6±6%) and the remaining parameters were not affected. This is the first study to provide evidences that exposure to PHE can affect early oyster development at the molecular and morphological levels, possibly threatening this bivalve species.

Assessment of housekeeping genes as internal references in quantitative expression analysis during early development of oyster

The early development of mollusks exhibits important characteristics from the developmental and evolutionary perspective. With the increasing number of genome-wide studies, accurate analyses of quantitative gene expression during development are impeded by the lack of validated reference genes. To improve the situation, in this study, we analyzed the expression stability of seven candidate housekeeping genes during early development of the Pacific oyster Crassostrea gigas: actin, glyceraldehyde-3-phosphate dehydrogenase (gapdh), α subunit of elongation factor 1 (elf1α), adp-ribosylation factor 1 (arf1), heterogeneous nuclear ribonucleoprotein q, ubiquitin-conjugating enzyme e2d2 and ribosomal protein s18. We focused on 11 stages from oocyte to D-veliger, which include crucial developmental processes such as axis determination, gastrulation and shell formation. Gene expression stabilities were assessed with the three commonly used programs geNorm, NormFinder and BestKeeper. Although the results obtained with the three programs varied to some extent, in general, arf1, elf1α and gapdh were highly ranked and actin was poorly ranked. This analysis also indicated that multiple genes should be used for normalization, and we concluded that arf1-elf1α-gapdh should be used as internal references. The findings of this study will help researchers to obtain accurate results in future quantitative gene expression analysis of development in bivalve mollusks.

Identification and characterization of a ferritin gene involved in the immune defense response of scallop Chlamys farreri

Scallop Chlamys farreri is an important aquaculture species in northern China. However, its mass mortality caused by several pathogens can result in great economic loss and negative impacts to the sustainable development of the scallop industry. Thus, improving the overall understanding of immune response mechanisms involved in host-pathogen interactions is necessary. Ferritins are conserved molecules in organisms that are involved in diverse biological processes, such as mediating host-pathogen responses. In this study, we report a novel ferritin gene from C. farreri (denoted as CfFER). The full length of CfFER is 848 bp and contains a 5'-UTR of 113 bp, a 3'-UTR of 219 bp, and a complete open reading frame (ORF) of 516 bp. The ORF encodes a polypeptide of 171 amino acid residues with a molecular weight of approximately 19.95 kDa and an isoelectric point of 5.07. The CfFER protein exhibited typical ferritin structures, namely, a ferroxidase diiron center, a ferrihydrite nucleation center, and an iron-binding response signature. Phylogenetic analysis revealed that CfFER was closely related to other mollusk ferritin proteins. Expression of CfFER in different tissues was analyzed by quantitative real-time PCR, and results showed that CfFER was ubiquitously expressed in all examined tissues. The highest and lowest expression levels of CfFER were measured in the muscle and hemocyte, respectively. The relative mRNA expression of CfFER in response to bacterial (Vibrio anguillarum) and viral (acute viral necrobiotic virus) challenges sharply increased by ca. 5-fold about12 h post-infection (hpi) and then normalized at 48 hpi. Western blot analysis with polyclonal antibodies generated from the recombinant product of CfFER also demonstrated the presence of ferritin protein in hemocytes. These findings strongly suggest that CfFER is involved in the immune response of C. farreri and protection against pathogen challenge.

Characterization and expression of a novel caspase gene: Evidence of the expansion of caspases in Crassostrea gigas

Caspases are a group of cysteine-aspartate proteases involved in apoptosis and a variety of non-apoptotic processes. In this study, a novel caspase gene was cloned and its potential role in apoptosis was investigated. The caspase gene (CgCasp 3/7) has an open reading frame of 1626bp encoding 541 amino acids containing the conserved functional domains and motifs of effector caspases. Its amino acid sequence shows low identity with the other effector caspases of Crassostrea gigas and contains a unique long intersubunit linker (IL). The CgCasp 3/7 mRNA was expressed highly in oocytes and then decreased gradually after fertilization, indicating CgCasp 3/7 could function in oocyte apoptosis. In adult tissues, it is located primarily in the gills and hepatopancreas. We examined the mRNA expression of CgCasp 3/7 in gills of oysters immersed in ambient (17°C) or heated (27°C) seawater. The thermal stress stimulated mRNA expression of CgCasp 3/7 by 2.5- and 4.1-fold at 2h and 6h post-treatment, respectively, indicating CgCasp3/7 was involved in the early response to thermal stress. To examine the function of the IL, CgCasp 3/7 and CgCasp 3/7-T (with a truncated IL) were expressed using an in vitro translation system and their DEVDase activity was measured. Both proteins showed a significantly higher level of DEVDase activity than control, but CgCasp3/7-T had lower DEVDase activity than CgCasp3/7, indicating CgCasp3/7 had DEVDase activity and the IL was required for maximal DEVDase activity. Our study adds to the complexity of caspases in C. gigas.

The importance of eukaryotic ferritins in iron handling and cytoprotection

Ferritins, the main intracellular iron storage proteins, have been studied for over 60 years, mainly focusing on the mammalian ones. This allowed the elucidation of the structure of these proteins and the mechanisms regulating their iron incorporation and mineralization. However, ferritin is present in most, although not all, eukaryotic cells, comprising monocellular and multicellular invertebrates and vertebrates. The aim of this review is to provide an update on the general properties of ferritins that are common to various eukaryotic phyla (except plants), and to give an overview on the structure, function and regulation of ferritins. An update on the animal models that were used to characterize H, L and mitochondrial ferritins is also provided. The data show that ferritin structure is highly conserved among different phyla. It exerts an important cytoprotective function against oxidative damage and plays a role in innate immunity, where it also contributes to prevent parenchymal tissue from the cytotoxicity of pro-inflammatory agonists released by the activation of the immune response activation. Less clear are the properties of the secretory ferritins expressed by insects and molluscs, which may be important for understanding the role played by serum ferritin in mammals.

Ferritin from the Pacific abalone Haliotis discus hannai: Analysis of cDNA sequence, expression, and activity

Ferritin plays an important role in iron homeostasis due to its ability to bind and sequester large amounts of iron. In this study, the gene encoding a ferritin (HdhFer2) was cloned from Pacific abalone (Haliotis discus hannai). The full-length cDNA of HdhFer2 contains a 5'-UTR of 121 bp, an ORF of 516 bp, and a 3'-UTR of 252 bp with a polyadenylation signal sequence of AATAAA and a poly(A) tail. It also contains a 31 bp iron-responsive element (IRE) in the 5'-UTR position, which is conserved in many ferritins. HdhFer2 consists of 171 amino acid residues with a predicted molecular weight (MW) ∼19.8 kDa and a theoretical isoelectric point (PI) of 4.84. The deduced amino acid sequence of HdhFer2 contains two ferritin iron-binding region signatures (IBRSs). HdhFer2 mRNA was detected in a wide range of tissues and was dominantly expressed in the gill. Infection with the bacterial pathogen Vibrio anguillarum significantly upregulated HdhFer2 expression in a time-dependent manner. Recombinant HdhFer2 (rHdhFer2) purified from Escherichia coli was able to bind ferrous iron in a concentration-dependent manner. In summary, these results suggest that HdhFer2 is a crucial protein in the iron-withholding defense system, and plays an important role in the innate immune response of abalone.

The venomous cocktail of the vampire snail Colubraria reticulata (Mollusca, Gastropoda)

Background

Hematophagy arose independently multiple times during metazoan evolution, with several lineages of vampire animals particularly diversified in invertebrates. However, the biochemistry of hematophagy has been studied in a few species of direct medical interest and is still underdeveloped in most invertebrates, as in general is the study of venom toxins. In cone snails, leeches, arthropods and snakes, the strong target specificity of venom toxins uniquely aligns them to industrial and academic pursuits (pharmacological applications, pest control etc.) and provides a biochemical tool for studying biological activities including cell signalling and immunological response. Neogastropod snails (cones, oyster drills etc.) are carnivorous and include active predators, scavengers, grazers on sessile invertebrates and hematophagous parasites; most of them use venoms to efficiently feed. It has been hypothesized that trophic innovations were the main drivers of rapid radiation of Neogastropoda in the late Cretaceous.

We present here the first molecular characterization of the alimentary secretion of a non-conoidean neogastropod, Colubraria reticulata. Colubrariids successfully feed on the blood of fishes, throughout the secretion into the host of a complex mixture of anaesthetics and anticoagulants. We used a NGS RNA-Seq approach, integrated with differential expression analyses and custom searches for putative secreted feeding-related proteins, to describe in detail the salivary and mid-oesophageal transcriptomes of this Mediterranean vampire snail, with functional and evolutionary insights on major families of bioactive molecules.

Results

A remarkably low level of overlap was observed between the gene expression in the two target tissues, which also contained a high percentage of putatively secreted proteins when compared to the whole body. At least 12 families of feeding-related proteins were identified, including: 1) anaesthetics, such as ShK Toxin-containing proteins and turripeptides (ion-channel blockers), Cysteine-rich secretory proteins (CRISPs), Adenosine Deaminase (ADA); 2) inhibitors of primary haemostasis, such as novel vWFA domain-containing proteins, the Ectonucleotide pyrophosphatase/phosphodiesterase family member 5 (ENPP5) and the wasp Antigen-5; 3) anticoagulants, such as TFPI-like multiple Kunitz-type protease inhibitors, Peptidases S1 (PS1), CAP/ShKT domain-containing proteins, Astacin metalloproteases and Astacin/ShKT domain-containing proteins; 4) additional proteins, such the Angiotensin-Converting Enzyme (ACE: vasopressive) and the cytolytic Porins.

Conclusions

Colubraria feeding physiology seems to involve inhibitors of both primary and secondary haemostasis, anaesthetics, a vasoconstrictive enzyme to reduce feeding time and tissue-degrading proteins such as Porins and Astacins. The complexity of Colubraria venomous cocktail and the divergence from the arsenal of the few neogastropods studied to date (mostly conoideans) suggest that biochemical diversification of neogastropods might be largely underestimated and worth of extensive investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1648-4) contains supplementary material, which is available to authorized users.