Mediterranean mussel gene expression profile induced by okadaic acid exposure

Subcellular effects and lipid metabolism alterations in the gilthead seabream Sparus aurata fed on ovatoxins-contaminated mussels

The marine microalgae Ostreopsis cf. ovata are a well-known producer of palytoxin (PlTXs) analogues, i.e. ovatoxins (OVTXs) among others, which arouse concern for animal and human health. Both in field and laboratory studies, presence of OVTXs, detected in species directly feeding on O. cf. ovata, was frequently correlated with impairment on organisms' physiology, development and behaviour, while similar knowledge is still lacking for animals feeding on contaminated preys. In this study, transfer and toxicity of OVTXs were evaluated in an exposure experiment, in which gilthead seabream Sparus aurata was fed with bivalve mussel Mytilus galloprovincialis, contaminated by a toxic strain of O. cf. ovata. Mussels exposed to O. cf. ovata for 21 days accumulated meanly 188 ± 13 μg/kg OVTXs in the whole tissues. Seabreams fed with OVTX-contaminated mussels started to reject the food after 6 days of contaminated diet. Although no detectable levels of OVTXs were measured in muscle, liver, gills and gastro-intestinal tracts, the OVTX-enriched diet induced alterations of lipid metabolism in seabreams livers, displaying a decreased content of total lipid and fatty acid, together with overexpression of fatty acid biosynthetic genes, downregulation of β-oxidation genes and modulation of several genes related to lipid transport and regulation. Results from this study would suggest the hypothesis that OVTXs produced by O. cf. ovata may not be subject to bioaccumulation in fish fed on contaminated preys, being however responsible of significant biological effects, with important implications for human consumption of seafood products.

Exosomal miRNA analysis provides new insights into exposure to nanoplastics and okadaic acid

As an emerging environmental pollutant, nanoplastics (NPs) have attracted wide attention in terms of their impact on the ecological environment and human health. Currently, researches on the cytotoxicity of NPs mainly focus on oxidative stress, damage to the cell membrane and organelles, induction of immune response and genotoxicity. Okadaic acid (OA) is the main component of diarrheal shellfish toxin. Based on the previous combined toxicity exploration of polystyrene (PS) NPs and (OA) to human gastric adenocarcinoma (AGS) cells, cell-derived exosomes were extracted and exosomal miRNA profiles were analyzed for the first time in this study. The results showed that the composition of miRNAs varied after the exposure of NPs and OA. Specifically, the expression of miR-1-3p in both PS-Exo and PS-OA-Exo was significantly reduced. And the expression of miR-1248 was upregulated most significantly by comparing the DE miRNAs between PS-Exo and PS-OA-Exo. MiR-1-3p and miR-1248 may be the key genes for the combined toxicity of NPs and OA. After analysis, we found that both the decreased expression of miR-1-3p and the increased expression of miR-1248 can increase the expression of FN1 and affect DNA replication, which was surprisingly consistent with the results of our previous cytotoxicity studies. Since exosomal miRNAs are selectively encapsulated by donor cell, we speculate that the changes of exosomal miRNAs may due to the synchronous changes of intracellular environment and the downregulation of intracellular FN1 may be attributed to decreased expression of miR-1-3p and increased expression of miR-1248 in donor cells. Accordingly, we come to the conclusion that the changes of miRNAs in the exosomes derived from AGS cells after environmental stimulation could reflect the biological effects of donor cells.

Expression Analyses of Genes Related to Multixenobiotic Resistance in Mytilus galloprovincialis after Exposure to Okadaic Acid-Producing Dinophysis acuminata

The mussel Mytilus galloprovincialis is one of the most important aquaculture species in Europe. Its main production problem is the accumulation of toxins during coastal blooms, which prevents mussel commercialization. P-glycoprotein (ABCB1/MDR1/P-gp) is part of the multixenobiotic resistance system in aquatic organisms, and okadaic acid, the main DSP toxin, is probably a substrate of the P-gp-mediated efflux. In this study, the presence and possible role of P-gp in the okadaic acid detoxification process was studied in M. galloprovincialis. We identified, cloned, and characterized two complete cDNAs of mdr1 and mdr2 genes. MgMDR1 and MgMDR2 predicted proteins had the structure organization of ABCB full transporters, and were identified as P-gp/MDR/ABCB proteins. Furthermore, the expression of mdr genes was monitored in gills, digestive gland, and mantle during a cycle of accumulation-elimination of okadaic acid. Mdr1 significantly increased its expression in the digestive gland and gills, supporting the idea of an important role of the MDR1 protein in okadaic acid efflux out of cells in these tissues. The expression of M. galloprovincialismrp2, a multidrug associated protein (MRP/ABCC), was also monitored. As in the case of mdr1, there was a significant induction in the expression of mrp2 in the digestive gland, as the content of okadaic acid increased. Thus, P-gp and MRP might constitute a functional defense network against xenobiotics, and might be involved in the resistance mechanisms to DSP toxins.

OMICs Approaches in Diarrhetic Shellfish Toxins Research

Diarrhetic shellfish toxins (DSTs) are among the most prevalent marine toxins in Europe's and in other temperate coastal regions. These toxins are produced by several dinoflagellate species; however, the contamination of the marine trophic chain is often attributed to species of the genus Dinophysis. This group of toxins, constituted by okadaic acid (OA) and analogous molecules (dinophysistoxins, DTXs), are highly harmful to humans, causing severe poisoning symptoms caused by the ingestion of contaminated seafood. Knowledge on the mode of action and toxicology of OA and the chemical characterization and accumulation of DSTs in seafood species (bivalves, gastropods and crustaceans) has significantly contributed to understand the impacts of these toxins in humans. Considerable information is however missing, particularly at the molecular and metabolic levels involving toxin uptake, distribution, compartmentalization and biotransformation and the interaction of DSTs with aquatic organisms. Recent contributions to the knowledge of DSTs arise from transcriptomics and proteomics research. Indeed, OMICs constitute a research field dedicated to the systematic analysis on the organisms' metabolisms. The methodologies used in OMICs are also highly effective to identify critical metabolic pathways affecting the physiology of the organisms. In this review, we analyze the main contributions provided so far by OMICs to DSTs research and discuss the prospects of OMICs with regard to the DSTs toxicology and the significance of these toxins to public health, food safety and aquaculture.

De novo transcriptome analysis of the mussel Perna viridis after exposure to the toxic dinoflagellate Prorocentrum lima

Diarrheic shellfish poisoning (DSP) toxins are produced by harmful microalgae and accumulate in bivalve mollusks, causing various toxicity. These toxic effects appear to abate with increasing DSP concentration and longer exposure time, however, the underlying mechanisms remain unclear. To explore the underlying molecular mechanisms, de novo transcriptome analysis of the digestive gland of Perna viridis was performed after Prorocentrum lima exposure. RNA-seq analysis showed that 1886 and 237 genes were up- and down-regulated, respectively after 6 h exposure to P. lima, while 265 genes were up-regulated and 217 genes were down-regulated after 96 h compared to the control. These differentially expressed genes mainly involved in Nrf2 signing pathways, immune stress, apoptosis and cytoskeleton, etc. Combined with qPCR results, we speculated that the mussel P. viridis might mainly rely on glutathione S-transferase (GST) and ABC transporters to counteract DSP toxins during short-term exposure. However, longer exposure of P. lima could activate the Nrf2 signaling pathway and inhibitors of apoptosis protein (IAP), which in turn reduced the damage of DSP toxins to the mussel. DSP toxins could induce cytoskeleton destabilization and had some negative impact on the immune system of bivalves. Collectively, our findings uncovered the crucial molecular mechanisms and the regulatory metabolic nodes that underpin the defense mechanism of bivalves against DSP toxins and also advanced our current understanding of bivalve defense mechanisms.

Effects of Marine Toxin Domoic Acid on Innate Immune Responses in Bay Scallop Argopecten irradians

Domoic acid (DA) is an amnesic shellfish poisoning toxin produced by some species of the genera Pseudo-nitzschia and Nitzschia. This toxin has harmful effects on various species, especially scallops. This study aimed to investigate the effects of DA exposure on the immune and physical responses of bay scallop, Argopecten irradians. Various immunological and physical parameters were assessed (acid phosphatase (ACP), alkaline phosphatase (ALP), lactate dehydrogenase (LDH), lipid peroxide (LPO), nitric oxide (NO), and the total protein content) in the haemolymph of scallops at 3, 6, 12, 24, and 48 h post-exposure to DA at different concentrations (10, 50, and 100 ng/mL). Moreover, the expression of immune-related genes (CLT-6, FREP, HSP90, MT, PGRP, and PrxV) was assessed. The activities of ACP, ALP, and LDH and the total protein content and LPO increased upon exposure to DA at different concentrations, while NO levels were decreased. Furthermore, immune-related genes were assessed upon DA exposure. Our results showed that exposure to DA negatively impacts immune function and disrupts physiological activities in bay scallops.

Responses of CYP450 in the mussel Perna viridis after short-term exposure to the DSP toxins-producing dinoflagellate Prorocentrum lima

Diarrhetic shellfish poisoning (DSP) toxins are key shellfish toxins that cause diarrhea, vomiting and even tumor. Interestingly, bivalves such as Perna viridis have been reported to exhibit some resistances to alleviate toxic effects of DSP toxins in a species-specific manner. Nevertheless, the molecular mechanisms underlying the resistance phenomenon to DSP toxins, particularly the mechanistic role of CYP450 is scant despite its crucial role in detoxification. Here, we exposed P. viridis to Prorocentrum lima and examined the expression pattern of the CYP450 and our comprehensive analyses revealed that P. lima exposure resulted in unique expression pattern of key CYP450 genes in bivalves. Exposure to P. lima (2 × 10⁵ cells/L) dramatically orchestrated the relative expression of CYP450 genes. CYP2D14-like mRNA was significantly down-regulated at 6 h in gill, but up-regulated at 2 h in digestive gland compared with control counterparts (p < 0.05), while CYP3A4 mRNA was increased at 12 h in gill. After exposure to P. lima at 2 × 10⁶ cells/L, the expression of CYP3A4 mRNA was significantly increased in digestive gland at 2 h and 12 h, while CYP2D14-like was up-regulated at 6 h. Besides, CYP3L3 and CYP2C8 also exhibited differential expression. These data suggested that CYP3A4, CYP2D14-like, and even CYP3L3 and CYP2C8 might be involved in DSP toxins metabolism. Besides, provision of ketoconazole resulted in significant decrement of CYP3A4 in digestive gland at 2 h and 12 h, while the OA content significantly decreased at 2 h and 6 h compared to control group without ketoconazole. These findings indicated that ketoconazole could depress CYP3A4 activity in bivalves thereby altering the metabolic activities of DSP toxins in bivalves, and also provided novel insights into the mechanistic role of CYP3A4 on DSP toxins metabolism in bivalves.

Up-regulation of Nrf2-dependent antioxidant defenses in Perna viridis after exposed to Prorocentrum lima

It is well documented that diarrhetic shellfish poisoning (DSP) toxins have strong genetic toxicity, cytotoxicity and oxidative damage to bivalve species. However, these toxic effects seem to decrease with the extension of exposure time and the increment of the toxin concentration, the mechanism involved remained unclear, though. In this paper, we found that expression of the genes related to cytoskeleton and Nrf2 signaling pathway displayed different changes over time in the gill of Perna viridis after exposure to DSP toxins-producing microalga Prorocentrum lima. During the short-term exposure (3 h and 6 h), KEAP1 gene expression was significantly up-regulated, coupled with up-regulation of MRP, ABCB1 and CAT transcriptions and down-regulation of GPx1 and NQO1 mRNA. After longer exposure to high density of P. lima, Nrf2 was significantly up-regulated, accompanied with up-regulation of Nrf2 pathway related genes such as NQO1, SOD, GST-ω and ABCB1, whereas KEAP1 was down-regulated. TUBA1C and TUBB1 transcripts were significantly down-regulated after short-term exposure of P. lima, but both of them were up-regulated at 96 h after exposure to high density of P. lima. Paraffin section demonstrated that P. lima had a strong damage on the gill of mussels during the short-term exposure. However, the negative effect to the gill decreased, and the gill restored after longer exposure (96 h). Taking together, we proposed that P. lima had a negative impact on cytoskeleton of mussel gill tissue, could cause oxidative damage to the gills. However, longer exposure of P. lima in high density could activate Nrf2 signaling pathway, thereby reducing the influence of toxin on mussel. Our study might provide a novel clue for the resistance mechanism of shellfish to DSP toxins.

RNA-Seq Transcriptome Profiling of the Queen Scallop (Aequipecten opercularis) Digestive Gland after Exposure to Domoic Acid-Producing Pseudo-nitzschia

Some species of the genus Pseudo-nitzschia produce the toxin domoic acid, which causes amnesic shellfish poisoning (ASP). Given that bivalve mollusks are filter feeders, they can accumulate these toxins in their tissues. To elucidate the transcriptional response of the queen scallop Aequipecten opercularis after exposure to domoic acid-producing Pseudo-nitzschia, the digestive gland transcriptome was de novo assembled using an Illumina HiSeq 2000 platform. Then, a differential gene expression analysis was performed. After the assembly, 142,137 unigenes were obtained, and a total of 10,144 genes were differentially expressed in the groups exposed to the toxin. Functional enrichment analysis found that 374 Pfam (protein families database) domains were significantly enriched. The C1q domain, the C-type lectin, the major facilitator superfamily, the immunoglobulin domain, and the cytochrome P450 were among the most enriched Pfam domains. Protein network analysis showed a small number of highly connected nodes involved in specific functions: proteasome components, mitochondrial ribosomal proteins, protein translocases of mitochondrial membranes, cytochromes P450, and glutathione S-transferases. The results suggest that exposure to domoic acid-producing organisms causes oxidative stress and mitochondrial dysfunction. The transcriptional response counteracts these effects with the up-regulation of genes coding for some mitochondrial proteins, proteasome components, and antioxidant enzymes (glutathione S-transferases, thioredoxins, glutaredoxins, and copper/zinc superoxide dismutases).

The insect, Galleria mellonella, is a compatible model for evaluating the toxicology of okadaic acid

The polyether toxin, okadaic acid, causes diarrhetic shellfish poisoning in humans. Despite extensive research into its cellular targets using rodent models, we know little about its putative effect(s) on innate immunity. We inoculated larvae of the greater wax moth, Galleria mellonella, with physiologically relevant doses of okadaic acid by direct injection into the haemocoel (body cavity) and/or gavage (force-feeding). We monitored larval survival and employed a range of cellular and biochemical assays to assess the potential harmful effects of okadaic acid. Okadaic acid at concentrations ≥ 75 ng/larva (≥ 242 μg/kg) led to significant reductions in larval survival (> 65%) and circulating haemocyte (blood cell) numbers (> 50%) within 24 h post-inoculation. In the haemolymph, okadaic acid reduced haemocyte viability and increased phenoloxidase activities. In the midgut, okadaic acid induced oxidative damage as determined by increases in superoxide dismutase activity and levels of malondialdehyde (i.e. lipid peroxidation). Our observations of insect larvae correspond broadly to data published using rodent models of shellfish-poisoning toxidrome, including complementary LD₅₀ values: 206–242 μg/kg in mice, ~ 239 μg/kg in G. mellonella. These data support the use of this insect as a surrogate model for the investigation of marine toxins, which offers distinct ethical and financial incentives.

RNA-Seq Analysis for Assessing the Early Response to DSP Toxins in Mytilus galloprovincialis Digestive Gland and Gill

The harmful effects of diarrhetic shellfish poisoning (DSP) toxins on mammalian cell lines have been widely assessed. Studies in bivalves suggest that mussels display a resistance to the cytogenotoxic effects of DSP toxins. Further, it seems that the bigger the exposure, the more resistant mussels become. To elucidate the early genetic response of mussels against these toxins, the digestive gland and the gill transcriptomes of Mytilus galloprovincialis after Prorocentrum lima exposure (100,000 cells/L, 48 h) were de novo assembled based on the sequencing of 8 cDNA libraries obtained using an Illumina HiSeq 2000 platform. The assembly provided 95,702 contigs. A total of 2286 and 4523 differentially expressed transcripts were obtained in the digestive gland and the gill, respectively, indicating tissue-specific transcriptome responses. These transcripts were annotated and functionally enriched, showing 44 and 60 significant Pfam families in the digestive gland and the gill, respectively. Quantitative PCR (qPCR) was performed to validate the differential expression patterns of several genes related to lipid and carbohydrate metabolism, energy production, genome integrity and defense, suggesting their participation in the protective mechanism. This work provides knowledge of the early response against DSP toxins in the mussel M. galloprovincialis and useful information for further research on the molecular mechanisms of the bivalve resistance to these toxins.

Detoxification- and Immune-Related Transcriptomic Analysis of Gills from Bay Scallops (Argopectenirradians) in Response to Algal Toxin Okadaic Acid

To reveal the molecular mechanisms triggered by okadaic acid (OA)-exposure in the detoxification and immune system of bay scallops, we studied differentially-expressed genes (DEGs) and the transcriptomic profile in bay scallop gill tissue after 48 h exposure to 500 nM of OA using the Illumina HiSeq 4000 deep-sequencing platform. De novo assembly of paired-end reads yielded 55,876 unigenes, of which 3204 and 2620 genes were found to be significantly up- or down-regulated, respectively. Gene ontology classification and enrichment analysis of the DEGs detected in bay scallops exposed to OA revealed four ontologies with particularly high functional enrichment, which were ‘cellular process’ (cellular component), ‘metabolic process’ (biological process), ‘immune system process’ (biological process), and ‘catalytic process’ (molecular function). The DEGs revealed that cyclic AMP-responsive element-binding proteins, acid phosphatase, toll-like receptors, nuclear erythroid 2-related factor, and the NADPH2 quinone reductase-related gene were upregulated. In contrast, the expression of some genes related to glutathione S-transferase 1, C-type lectin, complement C1q tumor necrosis factor-related protein, Superoxide dismutase 2 and fibrinogen C domain-containing protein, decreased. The outcomes of this study will be a valuable resource for the study of gene expression induced by marine toxins, and will help understanding of the molecular mechanisms underlying the scallops’ response to OA exposure.

Crassostrea gigas exposure to the dinoflagellate Prorocentrum lima: Histological and gene expression effects on the digestive gland

Bivalve mollusks bioaccumulate toxins via ingestion of toxic dinoflagellates. In this study, Crassostrea gigas was used to investigate the effects related to Prorocentrum lima exposure. Oysters were fed with three diets Isochrysis galbana (2 × 10(6) cell mL(-1)) control treatment; algal mix of I. galbana (2 × 10(6)) and P. lima (3 × 10(3) cell mL(-1)); and P. lima alone (3 × 10(3) cell mL(-1)). Feeding behavior changes, histopathological alterations, and expression patterns changes of genes involved in cell cycle (p21, cafp55, p53), cytoskeleton (tub, act), and inflammatory process (casp1) were evaluated. Results indicated that the presence of diarrheic shellfish poisoning by P. lima cells decreased the clearance rate (p < 0.05), induced structural loss, significantly decreased tubule area of the digestive gland (p < 0.05), and up-regulated in expression all gene (p < 0.05), suggesting that toxic cells might trigger inflammatory tissue process, disturb cell cycle and cytoskeleton representing a risk to oysters integrity.

Marine Toxin Okadaic Acid Affects the Immune Function of Bay Scallop (Argopecten irradians)

Okadaic acid (OA) is produced by dinoflagellates during harmful algal blooms and is a diarrhetic shellfish poisoning toxin. This toxin is particularly problematic for bivalves that are cultured for human consumption. This study aimed to reveal the effects of exposure to OA on the immune responses of bay scallop, Argopecten irradians. Various immunological parameters were assessed (total hemocyte counts (THC), reactive oxygen species (ROS), malondialdehyde (MDA), glutathione (GSH), lactate dehydrogenase (LDH), and nitric oxide (NO) in the hemolymph of scallops at 3, 6, 12, 24, and 48 h post-exposure (hpe) to different concentrations of OA (50, 100, and 500 nM). Moreover, the expression of immune-system-related genes (CLT-6, FREP, HSP90, MT, and Cu/ZnSOD) was also measured. Results showed that ROS, MDA, and NO levels and LDH activity were enhanced after exposure to different concentrations of OA; however, both THC and GSH decreased between 24-48 hpe. The expression of immune-system-related genes was also assessed at different time points during the exposure period. Overall, our results suggest that exposure to OA had negative effects on immune system function, increased oxygenic stress, and disrupted metabolism of bay scallops.

Vibrio diabolicus challenge in Bathymodiolus azoricus populations from Menez Gwen and Lucky Strike hydrothermal vent sites

Menez Gwen (MG) and Lucky Strike (LS) deep-sea hydrothermal vents are located at 850 m and 1730 m depths respectively and support chemosynthesis-based ecosystems partially differing in heavy metal concentration, temperature range, and faunistic composition. The successfully adapted deep-sea vent mussel Bathymodiolus azoricus is found at both vent locations. In such inhospitable environments survival strategies rely on the establishment of bacteria-vent animal symbiosis In spite of the toxic nature of deep-sea vents, the problem of microbial threat and the need for immunity exist in B. azoricus. This study aims at investigating the immune system of B. azoricus from MG and LS populations by comparing immune gene expressions profiles using the deep-sea vent-related Vibrio diabolicus. Expression of nineteen immune genes was analyzed from gill, digestive gland and mantle tissues upon 3 h, 12 h and 24 h V. diabolicus challenges. Based on quantitative-Polymerase Chain Reaction (qPCR) significant gene expression differences were found among MG and LS populations and challenge times MG mussels revealed that gill and digestive gland gene expression levels were remarkably higher than those from LS mussels. Expression of Carcinolectin, Serpin-2, SRCR, IRGs, RTK, TLR2, NF-κB, HSP70 and Ferritin genes was greater in MG than LS mussels. In contrast, mantle tissue from LS mussels revealed the highest peak of expression at 24 h for most genes analyzed. The activation of immune signaling pathways demonstrated that gene expression profiles are distinct between the two mussel populations. These differences may possibly ensue from intrinsic immune transcriptional activities upon which host responses are modulated in presence of V. diabolicus. mRNA transcript variations were assessed during 24 h acclimatization taking into account the partial depuration to which mussels were subjected to. Additionally, gene expression differences may reflect still accountable effects from the presence of vent remaining microfluidic environments within the tissues analyzed.

Unbiased high-throughput characterization of mussel transcriptomic responses to sublethal concentrations of the biotoxin okadaic acid

Background. Harmful Algal Blooms (HABs) responsible for Diarrhetic Shellfish Poisoning (DSP) represent a major threat for human consumers of shellfish. The biotoxin Okadaic Acid (OA), a well-known phosphatase inhibitor and tumor promoter, is the primary cause of acute DSP intoxications. Although several studies have described the molecular effects of high OA concentrations on sentinel organisms (e.g., bivalve molluscs), the effect of prolonged exposures to low (sublethal) OA concentrations is still unknown. In order to fill this gap, this work combines Next-Generation sequencing and custom-made microarray technologies to develop an unbiased characterization of the transcriptomic response of mussels during early stages of a DSP bloom.

Methods. Mussel specimens were exposed to a HAB episode simulating an early stage DSP bloom (200 cells/L of the dinoflagellate Prorocentrum lima for 24 h). The unbiased characterization of the transcriptomic responses triggered by OA was carried out using two complementary methods of cDNA library preparation: normalized and Suppression Subtractive Hybridization (SSH). Libraries were sequenced and read datasets were mapped to Gene Ontology and KEGG databases. A custom-made oligonucleotide microarray was developed based on these data, completing the expression analysis of digestive gland and gill tissues.

Results. Our findings show that exposure to sublethal concentrations of OA is enough to induce gene expression modifications in the mussel Mytilus. Transcriptomic analyses revealed an increase in proteasomal activity, molecular transport, cell cycle regulation, energy production and immune activity in mussels. Oppositely, a number of transcripts hypothesized to be responsive to OA (notably the Serine/Threonine phosphatases PP1 and PP2A) failed to show substantial modifications. Both digestive gland and gill tissues responded similarly to OA, although expression modifications were more dramatic in the former, supporting the choice of this tissue for future biomonitoring studies.

Discussion. Exposure to OA concentrations within legal limits for safe consumption of shellfish is enough to disrupt important cellular processes in mussels, eliciting sharp transcriptional changes as a result. By combining the study of cDNA libraries and a custom-made OA-specific microarray, our work provides a comprehensive characterization of the OA-specific transcriptome, improving the accuracy of the analysis of expresion profiles compared to single-replicated RNA-seq methods. The combination of our data with related studies helps understanding the molecular mechanisms underlying molecular responses to DSP episodes in marine organisms, providing useful information to develop a new generation of tools for the monitoring of OA pollution.

Environmental epigenetics: A promising venue for developing next-generation pollution biomonitoring tools in marine invertebrates

Environmental epigenetics investigates the cause-effect relationships between specific environmental factors and the subsequent epigenetic modifications triggering adaptive responses in the cell. Given the dynamic and potentially reversible nature of the different types of epigenetic marks, environmental epigenetics constitutes a promising venue for developing fast and sensible biomonitoring programs. Indeed, several epigenetic biomarkers have been successfully developed and applied in traditional model organisms (e.g., human and mouse). Nevertheless, the lack of epigenetic knowledge in other ecologically and environmentally relevant organisms has hampered the application of these tools in a broader range of ecosystems, most notably in the marine environment. Fortunately, that scenario is now changing thanks to the growing availability of complete reference genome sequences along with the development of high-throughput DNA sequencing and bioinformatic methods. Altogether, these resources make the epigenetic study of marine organisms (and more specifically marine invertebrates) a reality. By building on this knowledge, the present work provides a timely perspective highlighting the extraordinary potential of environmental epigenetic analyses as a promising source of rapid and sensible tools for pollution biomonitoring, using marine invertebrates as sentinel organisms. This strategy represents an innovative, groundbreaking approach, improving the conservation and management of natural resources in the oceans.

Expression profile of eight glutathione S-transferase genes in Crassostrea ariakensis after exposure to DSP toxins producing dinoflagellate Prorocentrum lima

In this study, changes in eight GSTs mRNA level including GST-α, GST-σ, GST-ω, GST-π, GST-μ, GST-ρ, GST-θ and microsomal GST (mGST) in the oyster Crassostrea ariakensis after exposure to Prorocentrum lima have been evaluated by quantitative real-time PCR. Additionally, the contents of five GST isoforms were detected by ELISA. After exposure to P. lima at density of 2 × 10(5) cells/L, mGST mRNA significantly increased in gill, while GST-σ was induced in digestive gland. After exposure to P. lima at density of 2 × 10(6) cells/L, GST-ω and mGST expressions increased in gill, whereas GST-α and GST-σ were induced in digestive gland. The GST content and activity in oysters exposed to P. lima also showed a different pattern when the different isoforms and organs were compared. After exposure to P. lima (2 × 10(6) cell/L), GST-π increased in gill but decreased in digestive gland. The total GST enzyme activity increased in gill, while remained unchanged in digestive gland. These various regulation of GST gene expressions indicated that the GSTs isoenzymes might play divergent physiological roles in the detoxification of DSP toxins in C. ariakensis.

Saxitoxins and okadaic acid group: accumulation and distribution in invertebrate marine vectors from Southern Chile

Harmful algae blooms (HABs) are the main source of marine toxins in the aquatic environment surrounding the austral fjords in Chile. Huichas Island (Aysén) has an history of HABs spanning more than 30 years, but there is limited investigation of the bioaccumulation of marine toxins in the bivalves and gastropods from the Region of Aysén. In this study, bivalves (Mytilus chilenses, Choromytilus chorus, Aulacomya ater, Gari solida, Tagelus dombeii and Venus antiqua) and carnivorous gastropods (Argobuccinum ranelliformes and Concholepas concholepas) were collected from 28 sites. Researchers analysed the accumulation of STX-group toxins using a LC with a derivatisation post column (LC-PCOX), while lipophilic toxins (OA-group, azapiracids, pectenotoxins and yessotoxins) were analysed using LC-MS/MS with electrospray ionisation (+/-) in visceral (hepatopancreas) and non-visceral tissues (mantle, adductor muscle, gills and foot). Levels of STX-group and OA-group toxins varied among individuals from the same site. Among all tissue samples, the highest concentrations of STX-group toxins were noted in the hepatopancreas in V. antiqua (95 ± 0.1 μg STX-eq 100 g(-1)), T. dombeii (148 ± 1.4 μg STX-eq 100 g(-1)) and G. solida (3232 ± 5.2 μg STX-eq 100 g(-1); p < 0.05); in the adductor muscle in M. chilensis (2495 ± 6.4 μg STX-eq 100 g(-1); p < 0.05) and in the foot in C. concholepas (81 ± 0.7 μg STX-eq 100 g(-1)) and T. dombeii (114 ± 1.2 μg STX-eq 100 g(-1)). The highest variability of toxins was detected in G. solida, where high levels of carbamate derivatives were identified (GTXs, neoSTX and STX). In addition to the detected hydrophilic toxins, OA-group toxins were detected (OA and DTX-1) with an average ratio of ≈1:1. The highest levels of OA-group toxins were in the foot of C. concholepas, with levels of 400.3 ± 3.6 μg OA eq kg(-1) (p < 0.05) and with a toxic profile composed of 90% OA. A wide range of OA-group toxins was detected in M. chilensis with a toxicity < 80 μg OA eq kg(-1), but with 74% of those toxins detected in the adductor muscle. In all evaluated species, there was no detection of lipophilic toxins associated with biotransformation in molluscs and carnivorous gastropods. In addition, the STX-group and OA-group toxin concentrations in shellfish was not associated with the presence of HAB. The ranking of toxin concentration in the tissues of most species was: digestive glands > mantle > adductor muscle for the STX-group toxins and foot > digestive gland for the OA-group toxins. These results gave a better understanding of the variability and compartmentalisation of STX-group and OA-group toxins in different bivalve and gastropod species from the south of Chile, and the analyses determined that tissues could play an important role in the biotransformation of STX-group toxins and the retention of OA-group toxins.

Proteomic profile in Perna viridis after exposed to Prorocentrum lima, a dinoflagellate producing DSP toxins

In the current study, we compared protein profiles in gills of Perna viridis after exposure to Prorocentrumlima, a dinoflagellate producing DSP toxins, and identified the differential abundances of protein spots using 2D-electrophoresis. After exposure to P. lima, the level of okadaic acid (a main component of DSP toxins) in gills of P. viridis significantly increased at 6 h, but mussels were all apparently healthy without death. Among the 28 identified protein spots by MALDI TOF/TOF-MS, 12 proteins were up-regulated and 16 were down-regulated in the P. lima-exposed mussels. These identified proteins were involved in various biological activities, such as metabolism, cytoskeleton, signal transduction, response to oxidative stress and detoxification. Taken together, our results indicated that the presence of P. lima caused DSP toxins accumulation in mussel gill, and might consequently induce cytoskeletonal disorganization,oxidative stress, a dysfunction in metabolism and ubiquitination/proteasome activity.

RNA sequencing and de novo assembly of the digestive gland transcriptome in Mytilus galloprovincialis fed with toxinogenic and non-toxic strains of Alexandrium minutum

Background

The Mediterranean mussel Mytilus galloprovincialis is marine bivalve with a relevant commercial importance as well as a key sentinel organism for the biomonitoring of environmental pollution. Here we report the RNA sequencing of the mussel digestive gland, performed with the aim: a) to produce a high quality de novo transcriptome assembly, thus improving the genetic and molecular knowledge of this organism b) to provide an initial assessment of the response to paralytic shellfish poisoning (PSP) on a molecular level, in order to identify possible molecular markers of toxin accumulation.

Results

The comprehensive de novo assembly and annotation of the transcriptome yielded a collection of 12,079 non-redundant consensus sequences with an average length of 958 bp, with a high percentage of full-length transcripts. The whole-transcriptome gene expression study indicated that the accumulation of paralytic toxins produced by the dinoflagellate Alexandrium minutum over a time span of 5 days scarcely affected gene expression, but the results need further validation with a greater number of biological samples and naturally contaminated specimens.

Conclusion

The digestive gland reference transcriptome we produced significantly improves the data collected from previous sequencing efforts and provides a basic resource for expanding functional genomics investigations in M. galloprovincialis. Although not conclusive, the results of the RNA-seq gene expression analysis support the classification of mussels as bivalves refractory to paralytic shellfish poisoning and point out that the identification molecular biomarkers of PSP in the digestive gland of this organism is problematic.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-722) contains supplementary material, which is available to authorized users.

Correlation of transcriptomic responses and metal bioaccumulation in Mytilus edulis L. reveals early indicators of stress

Marine biomonitoring programs in the U.S. and Europe have historically relied on monitoring tissue concentrations of bivalves to monitor contaminant levels and ecosystem health. By integrating 'omic methods with these tissue residue approaches we can uncover mechanistic insight to link tissue concentrations to potential toxic effects. In an effort to identify novel biomarkers and better understand the molecular toxicology of metal bioaccumulation in bivalves, we exposed the blue mussel, Mytilus edulis L., to sub-lethal concentrations (0.54 μM) of cadmium, lead, and a Cd+Pb mixture. Metal concentrations were measured in gill tissues at 1, 2, and 4 weeks, and increased linearly over the 4 week duration. In addition, there was evidence that Pb interfered with Cd uptake in the mixture treatment. Using a 3025 sequence microarray for M. edulis, we performed transcriptomic analysis, identifying 57 differentially expressed sequences. Hierarchical clustering of these sequences successfully distinguished the different treatment groups demonstrating that the expression profiles were reproducible among the treatments. Enrichment analysis of gene ontology terms identified several biological processes that were perturbed by the treatments, including nucleoside phosphate biosynthetic processes, mRNA metabolic processes, and response to stress. To identify transcripts whose expression level correlated with metal bioaccumulation, we performed Pearson correlation analysis. Several transcripts correlated with gill metal concentrations including mt10, mt20, and contig 48, an unknown transcript containing a wsc domain. In addition, three transcripts directly involved in the unfolded protein response (UPR) were induced in the metal treatments at 2 weeks and were further up-regulated at 4 weeks. Overall, correlation of tissue concentrations and gene expression responses indicates that as mussels accumulate higher concentrations of metals, initial stress responses are mobilized to protect tissues. However, given the role of UPR in apoptosis, it serves as an early indicator of stress, which once overwhelmed will result in adverse physiological effects.

P-glycoprotein expression in Perna viridis after exposure to Prorocentrum lima, a dinoflagellate producing DSP toxins

Bivalves naturally exposed to toxic algae have mechanisms to prevent from harmful effects of diarrhetic shellfish poisoning (DSP) toxins. However, quite few studies have examined the mechanisms associated, and the information currently available is still insufficient. Multixenobiotic resistance (MXR) is ubiquitous in aquatic invertebrates and plays an important role in defense against xenobiotics. Here, to explore the roles of P-glycoprotein (P-gp) in the DSP toxins resistance in shellfish, complete cDNA of P-gp gene in the mussel Perna viridis was cloned and analyzed. The accumulation of okadaic acid (OA), a main component of DSP toxins, MXR activity and expression of P-gp in gills of P. viridis were detected after exposure to Prorocentrum lima, a dinoflagellate producing DSP toxins in the presence or absence of P-gp inhibitors PGP-4008, verapamil (VER) and cyclosporin A (CsA). The mussel P. viridis P-gp closely matches MDR/P-gp/ABCB protein from various organisms, having a typical sequence organization as full transporters from the ABCB family. After exposure to P. lima, OA accumulation, MXR activity and P-gp expression significantly increased in gills of P. viridis. The addition of P-gp-specific inhibitors PGP-4008 and VER decreased MXR activity induced by P. lima, but had no effect on the OA accumulation in gills of P. viridis. However, CsA, a broad-spectrum inhibitor of ABC transporter not only decreased MXR activity, but also increased OA accumulation in gills of P. viridis. Together with the ubiquitous presence of other ABC transporters such as MRP/ABCC in bivalves and potential compensatory mechanism in P-gp and MRP-mediated resistance, we speculated that besides P-gp, other ABC transporters, especially MRP might be involved in the resistance mechanisms to DSP toxins.

Effects of in vitro exposure to diarrheic toxin producer Prorocentrum lima on gene expressions related to cell cycle regulation and immune response in Crassostrea gigas

Background

Crassostrea gigas accumulates diarrheic shellfish toxins (DSP) associated to Prorocentrum lima of which Okadaic acid (OA) causes specific inhibitions of serine and threonine phosphatases 1 and 2A. Its toxic effects have been extensively reported in bivalve mollusks at cellular and physiological levels, but genomic approaches have been scarcely studied.

Methodology/Principal Findings

Acute and sub-chronic exposure effects of P. lima were investigated on farmed juvenile C. gigas (3–5 mm). The Pacific oysters were fed with three dinoflagellate concentrations: 0.3, 3, and 30×10³ cells mL⁻¹ along with a nontoxic control diet of Isochrysis galbana. The effects of P. lima on C. gigas were followed by analyzing expression levels of a total of four genes, three involved in cell cycle regulation and one in immune response by polymerase chain reaction and real time quantitative PCR, where changes in time and cell concentration were found. The highest expression levels were found in oysters fed 3×10³ cells mL⁻¹ at 168 h for the cycle regulator p21 protein (9 fold), chromatin assembly factor 1 p55 subunit (8 fold), elongation factor 2 (2 fold), and lipopolysaccharide/β-1, 3 glucan binding protein (13 fold above base line). Additionally, the transcript level of all the genes decreased in oysters fed wich the mixed diet 30×10³ cells mL⁻¹ of dinoflagellate after 72 h and was lowest in the chromatin assembly factor 1 p55 subunit (0.9 fold below baseline).

Conclusions

On C. gigas the whole cell ingestion of P lima caused a clear mRNA modulation expression of the genes involved in cell cycle regulation and immune system. Over-expression could be related to DNA damage, disturbances in cell cycle continuity, probably a genotoxic effect, as well as an activation of its innate immune system as first line of defense.

Gene expression profiles in zebrafish (Danio rerio) liver after acute exposure to okadaic acid

Okadaic acid (OA), a main component of diarrheic shellfish poisoning (DSP) toxins, is a strong and specific inhibitor of the serine/threonine protein phosphatases PP1 and PP2A. However, not all of the OA-induced effects can be explained by this phosphatase inhibition, and controversial results on OA are increasing. To provide clues on potential mechanisms of OA other than phosphatase inhibition, here, acute toxicity of OA was evaluated in zebrafish, and changes in gene expression in zebrafish liver tissues upon exposure to OA were observed by microarray. The i.p. ED50 (6 h) of OA on zebrafish was 1.54 μg OA/g body weight (bw). Among the genes analyzed on the zebrafish array, 55 genes were significantly up-regulated and 36 down-regulated in the fish liver tissue upon exposure to 0.176 μg OA/g bw (low-dose group, LD) compared with the low ethanol control (LE). However, there were no obvious functional clusters for them. On the contrary, fish exposure to 1.760 μg OA/g bw (high-dose group, HD) yielded a great number of differential expressed genes (700 up and 285 down) compared with high ethanol control (HE), which clustered in several functional terms such as p53 signaling pathway, Wnt signaling pathway, glutathione metabolism and protein processing in endoplasmic reticulum, etc. These genes were involved in protein phosphatase activity, translation factor activity, heat shock protein binding, as well as transmembrane transporter activity. Our findings may give some useful information on the pathways of OA-induced injury in fish.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

Okadaic acid meet and greet: an insight into detection methods, response strategies and genotoxic effects in marine invertebrates

Harmful Algal Blooms (HABs) constitute one of the most important sources of contamination in the oceans, producing high concentrations of potentially harmful biotoxins that are accumulated across the food chains. One such biotoxin, Okadaic Acid (OA), is produced by marine dinoflagellates and subsequently accumulated within the tissues of filtering marine organisms feeding on HABs, rapidly spreading to their predators in the food chain and eventually reaching human consumers causing Diarrhetic Shellfish Poisoning (DSP) syndrome. While numerous studies have thoroughly evaluated the effects of OA in mammals, the attention drawn to marine organisms in this regard has been scarce, even though they constitute primary targets for this biotoxin. With this in mind, the present work aimed to provide a timely and comprehensive insight into the current literature on the effect of OA in marine invertebrates, along with the strategies developed by these organisms to respond to its toxic effect together with the most important methods and techniques used for OA detection and evaluation.

The CHROMEVALOA database: a resource for the evaluation of Okadaic Acid contamination in the marine environment based on the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis

Okadaic Acid (OA) constitutes the main active principle in Diarrhetic Shellfish Poisoning (DSP) toxins produced during Harmful Algal Blooms (HABs), representing a serious threat for human consumers of edible shellfish. Furthermore, OA conveys critical deleterious effects for marine organisms due to its genotoxic potential. Many efforts have been dedicated to OA biomonitoring during the last three decades. However, it is only now with the current availability of detailed molecular information on DNA organization and the mechanisms involved in the maintenance of genome integrity, that a new arena starts opening up for the study of OA contamination. In the present work we address the links between OA genotoxicity and chromatin by combining Next Generation Sequencing (NGS) technologies and bioinformatics. To this end, we introduce CHROMEVALOAdb, a public database containing the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis (a sentinel model organism) in response to OA exposure. This resource constitutes a leap forward for the development of chromatin-based biomarkers, paving the road towards the generation of powerful and sensitive tests for the detection and evaluation of the genotoxic effects of OA in coastal areas.

DNA damage and transcriptional changes in the gills of mytilus galloprovincialis exposed to nanomolar doses of combined metal salts (Cd, Cu, Hg)

Aiming at an integrated and mechanistic view of the early biological effects of selected metals in the marine sentinel organism Mytilus galloprovincialis, we exposed mussels for 48 hours to 50, 100 and 200 nM solutions of equimolar Cd, Cu and Hg salts and measured cytological and molecular biomarkers in parallel. Focusing on the mussel gills, first target of toxic water contaminants and actively proliferating tissue, we detected significant dose-related increases of cells with micronuclei and other nuclear abnormalities in the treated mussels, with differences in the bioconcentration of the three metals determined in the mussel flesh by atomic absorption spectrometry. Gene expression profiles, determined in the same individual gills in parallel, revealed some transcriptional changes at the 50 nM dose, and substantial increases of differentially expressed genes at the 100 and 200 nM doses, with roughly similar amounts of up- and down-regulated genes. The functional annotation of gill transcripts with consistent expression trends and significantly altered at least in one dose point disclosed the complexity of the induced cell response. The most evident transcriptional changes concerned protein synthesis and turnover, ion homeostasis, cell cycle regulation and apoptosis, and intracellular trafficking (transcript sequences denoting heat shock proteins, metal binding thioneins, sequestosome 1 and proteasome subunits, and GADD45 exemplify up-regulated genes while transcript sequences denoting actin, tubulins and the apoptosis inhibitor 1 exemplify down-regulated genes). Overall, nanomolar doses of co-occurring free metal ions have induced significant structural and functional changes in the mussel gills: the intensity of response to the stimulus measured in laboratory supports the additional validation of molecular markers of metal exposure to be used in Mussel Watch programs.

Effect of okadaic acid on cultured clam heart cells: involvement of MAPkinase pathways

Okadaic acid (OA) is one of the main diarrhetic shellfish poisoning toxins and a potent inhibitor of protein phosphatases 1 and 2A. The downstream signal transduction pathways following the protein phosphatase inhibition are still unknown and the results of most of the previous studies are often conflicting. The aim of the present study was to evaluate the effects of OA on heart clam cells and to analyse its possible mechanisms of action by investigating the signal transduction pathways involved in OA cytotoxicity. We showed that OA at 1 µM after 24 h of treatment induces disorganization of the actin cytoskeleton, rounding and detachment of fibroblastic cells. Moreover, treatment of heart cells revealed a sequential activation of MAPK proteins depending on the OA concentration. We suggest that the duration of p38 and JNK activation is a critical factor in determining cell apoptosis in clam cardiomyocytes. In the opposite, ERK activation could be involved in cell survival. The cell death induced by OA is a MAPK modulated pathway, mediated by caspase 3-dependent mechanism. OA was found to induce no significant effect on spontaneous beating rate or inward L-type calcium current in clam cardiomyocytes, suggesting that PP1 was not inhibited even by the highest dose of OA.

Massively parallel amplicon sequencing reveals isotype-specific variability of antimicrobial peptide transcripts in Mytilus galloprovincialis

Background

Effective innate responses against potential pathogens are essential in the living world and possibly contributed to the evolutionary success of invertebrates. Taken together, antimicrobial peptide (AMP) precursors of defensin, mytilin, myticin and mytimycin can represent about 40% of the hemocyte transcriptome in mussels injected with viral-like and bacterial preparations, and unique profiles of myticin C variants are expressed in single mussels. Based on amplicon pyrosequencing, we have ascertained and compared the natural and Vibrio-induced diversity of AMP transcripts in mussel hemocytes from three European regions.

Methodology/Principal Findings

Hemolymph was collected from mussels farmed in the coastal regions of Palavas (France), Vigo (Spain) and Venice (Italy). To represent the AMP families known in M. galloprovincialis, nine transcript sequences have been selected, amplified from hemocyte RNA and subjected to pyrosequencing. Hemolymph from farmed (offshore) and wild (lagoon) Venice mussels, both injected with 10⁷Vibrio cells, were similarly processed. Amplicon pyrosequencing emphasized the AMP transcript diversity, with Single Nucleotide Changes (SNC) minimal for mytilin B/C and maximal for arthropod-like defensin and myticin C. Ratio of non-synonymous vs. synonymous changes also greatly differed between AMP isotypes. Overall, each amplicon revealed similar levels of nucleotidic variation across geographical regions, with two main sequence patterns confirmed for mytimycin and no substantial changes after immunostimulation.

Conclusions/Significance

Barcoding and bidirectional pyrosequencing allowed us to map and compare the transcript diversity of known mussel AMPs. Though most of the genuine cds variation was common to the analyzed samples we could estimate from 9 to 106 peptide variants in hemolymph pools representing 100 mussels, depending on the AMP isoform and sampling site. In this study, no prevailing SNC patterns related to geographical origin or Vibrio injection emerged. Whether or not the contact with potential pathogens can increase the amount of AMP transcript variants in mussels requires additional study.