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

Rat tight junction proteins are disrupted after subchronic exposure to okadaic acid

Okadaic acid (OA), a lipophilic phycotoxin distributed worldwide, causes diarrheic shellfish poisoning and even leads to tumor formation. Currently, the consumption of contaminated seafood is the most likely cause of chronic OA exposure, but there is a serious lack of relevant data. Here, the Sprague-Dawley rats were exposure to OA by oral administration at 100 µg/kg body weight, and the tissues were collected and analyzed to assess the effect of subchronic OA exposure. The results showed that subchronic OA administration disturbed colonic mucosal integrity and induced colitis. The colonic tight junction proteins were disrupted and the cell cycle of colonic epithelial cells was accelerated. It is inferred that disruption of the colonic tight junction proteins might be related to the development of chronic diarrhea by affecting water and ion transport. Moreover, the accelerated proliferation of colonic epithelial cells indicated that subchronic OA exposure might promote the restitution process of gut barrier or induce tumor promoter activity in rat colon.

Responses of ABCB and ABCC transporters to the toxic dinoflagellate Prorocentrum lima in the mussel Perna viridis

Bivalve mollusks can accumulate diarrheic shellfish poisoning (DSP) toxins through filter-feeding, but they exhibit some resistance to the toxins. Previous studies have suggested that the ABC transporters may have an important role in the resistance to DSP toxins, but comprehensive studies are lacking. In this study, we comprehensively analyzed the distribution of ABC transporters in the mussel Perna viridis, and observed responses of ABCB and ABCC transporters to the DSP toxins-producing dinoflagellate Prorocentrum lima. Total 39 members of ABC transporters were identified in P. viridis, including 3 full PvABCBs, 3 half PvABCBs, and 7 PvABCCs transporters. We found that PvABCBs and PvABCCs subfamilies were expressed in hemocytes, gills and digestive gland with some difference, especially in hemocytes. After exposure to P. lima, PvABCBs and PvABCCs displayed different expression changes in different tissues. The short-term (3 h) exposure to P. lima induced the transcription of PvABCB1_like1, PvABCB6, PvABCC1, PvABCC1_like and PvABCC1/3, and the longer-term (96 h) exposure increased the transcription of PvABCB1, PvABCB1_like, PvABCB10, PvABCC1 and PvABCC1_like1 in gills and PvABCC10 in digestive gland. These results suggest that different types of PvABCBs and PvABCCs in P. viridis may contribute to the detoxification of DSP toxins in different tissues at different time after exposure to DSP toxins. Our finding provides new evidence for further understanding the role of ABC transporters in the tolerance of mussel to DSP toxins.

Toxic Responses of Different Shellfish Species after Exposure to Prorocentrum lima, a DSP Toxins Producing Dinoflagellate

Prorocentrum lima is a global benthic dinoflagellate that produces diarrhetic shellfish poisoning (DSP) toxins, which can be ingested by filter-feeding bivalves, and eventually pose a great threat to human health through food chain. After being exposed to P. lima, different bivalves may accumulate various levels of DSP toxins and display different toxic responses. However, the underlying mechanism remains unclear. Here, we found that the content of okadaic acid-equivalents (OA-eq) varied in the digestive glands of the three bivalves including Crassostrea gigas, Mytilus coruscus and Tegillarca granosa after P. lima exposure. The degree of esterification of OA-eq in the three bivalves were opposite to the accumulation of OA-eq. The digestive gland tissues of the three bivalve species were damaged to different degrees. The transcriptional induction of Nrf2 targeted genes such as ABCB1 and GPx indicates the functionality of Nrf2 pathway against DSP toxins in bivalves. The oyster could protect against DSP toxins mainly through ABC transporters and esterification, while the mussel and clam reduce the damage induced by DSP toxins mainly by regulating the expression of antioxidant genes. Our findings may provide some explanations for the difference in toxic response to DSP toxins in different shellfish.

Transcriptomic Profile of the Cockle Cerastoderma edule Exposed to Seasonal Diarrhetic Shellfish Toxin Contamination

Bivalves constitute an important source of proteins for human consumption, but some accumulate biotoxins such as diarrhetic shellfish toxins (DSTs), constituting a risk to human health. The cockle Cerastoderma edule is one of the most important species harvested in the Portuguese coast but also one of the most affected species due to recurrent DSTs exposure. However, little is known regarding the effects of the toxins produced by blooming dinoflagellates on C. edule. Herein, we explore the Differentially Expressed Genes (DEGs) of two tissues (gills and digestive gland) from wild cockles sampled in Portugal, through their whole transcriptomic response in two different seasons (exposed and not exposed to DSTs). The de novo transcriptome assembly returned 684,723 contigs, N50 of 1049, and 98.53% completeness. Altogether, 1098 DEGs were identified, of which 353 DEGs were exclusive for the digestive gland, 536 unique for the gills and 209 DEGs were common. Among DEGs were identified known DSTs-biomarkers including glutathione peroxidase, glutathione S-transferase, superoxide dismutase, cytochrome P450, ABC transporters, actin and tubulin-related proteins, Heat shock proteins and complement C1Q-like proteins. This study provides the first transcriptomic profile of C. edule, giving new insights about its molecular responses under different environmental conditions of DSTs exposure.

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.

Effects of the toxic dinoflagellates Prorocentrum lima and Ostreopsis cf. ovata on immune responses of cultured oysters Crassostrea gasar

Oyster production in Brazil has been highlighted as an important economic activity and is directly impacted by the quality of the environment, which is largely the result of human interference and climate change. Harmful algal blooms occur in aquatic ecosystems worldwide, including coastal marine environments which have been increasing over the last decades as a result of global change and anthropogenic activities. In this study, the native oysters Crassostrea gasar from Northeast of Brazil were exposed to two toxic benthic dinoflagellate species, Prorocentrum lima and Ostreopsis cf. ovata. Their respective effects on C. gasar physiology and defense mechanisms were investigated. Oyster hemocytes were first exposed in vitro to different concentrations of both dinoflagellate species to assess their effects on hemocyte functions, such as phagocytosis, production of reactive oxygen species, as well as mortality. Results highlighted an alteration of hemocyte phagocytosis and viability in presence of O. cf. ovata, whereas P. lima did not affect the measured hemocyte functions. In a second experiment, oysters were exposed for 4 days in vivo to toxic culture of O. cf. ovata to assess its effects on hemocyte parameters, tissues damages and pathogenic Perkinsus spp. infection. An increase in hemocyte mortality was also observed in vivo, associated with a decrease of ROS production. Histopathological analyses demonstrated a thinning of the epithelium of the digestive tubules of the digestive gland, inflammatory reaction and a significant increase in the level of infection by Perkinsus spp. in oysters exposed to O. cf. ovata. These results indicate that oysters C. gasar seem to be pretty resilient to an exposure to P. lima and may be more susceptible to O. cf. ovata. Furthermore, the latter clearly impaired oyster physiology and defense mechanisms, thus highlighting that harmful algal blooms of O. cf. ovata could potentially lead to increased susceptibility of C. gasar oysters to parasite infections.

Cinnamaldehyde Could Reduce the Accumulation of Diarrhetic Shellfish Toxins in the Digestive Gland of the Mussel Perna viridis under Laboratory Conditions

Diarrhetic shellfish toxins (DSTs), some of the most important phycotoxins, are distributed almost all over the world, posing a great threat to human health through the food chain. Therefore, it is of great significance to find effective methods to reduce toxin accumulation in shellfish. In this paper, we observed the effects of four phytochemicals including cinnamaldehyde (CA), quercetin, oridonin and allicin on the accumulation of DSTs in the digestive gland of Perna viridis after exposure to the DSTs-producing Prorocentrum lima. We found that, among the four phytochemicals, CA could effectively decrease the accumulation of DSTs (okadaic acid-eq) in the digestive gland of P. viridis. Further evidence demonstrated that CA could reduce the histological alterations of the digestive gland of a mussel caused by DSTs. RT-qPCR showed that CA could suppress the CYP3A4 induction by DSTs, suggesting that the DSTs’ decrease induced by CA might be related to the inhibition of CYP3A4 transcription induction. However, further studies on the underlying mechanism, optimal treatment time, ecological safety and cost should be addressed before cinnamaldehyde is used to decrease the accumulation of DSTs in field.

Responses of cytochrome P450, GST and MXR in the mussel Perna viridis to the exposure of Aureococcus anophagefferens

The brown tide formed by a microscopic alga called Aureococcus anophagefferens has a devastating effect on filter-feeding bivalves, however, the related toxic principle remains an open question. In this study, we found that A. anophagefferens cells could motivate detoxification associated genes including CYP450, GST, P-gp and MVP, and induce SOD activity in the mussel Perna viridis. D1-like and D2-like receptors were expressed at high level in the gills of P. viridis, however, D2-like receptor transcript was too low to detect in digestive gland. The exposure of A. anophagefferens did not lead to any significant alterations in the expression of D1-like and D2-like receptors in both gills and digestive gland. These findings suggested that A. anophagefferens exhibited cytotoxicity toward bivalves, but did not obviously disrupt the dopamine system at transcriptional level in the acute exposure. Further studies are warranted to explore the nature of toxic compounds in A. anophagefferens affected bivalves.

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.

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.

Detoxification and Immune Transcriptomic Response of the Gill Tissue of Bay Scallop (Argopecten irradians) Following Exposure to the Algicide Palmitoleic Acid

Palmitoleic acid (PA) is an effective algicide against Alexandrium tamarense. However, the toxicological mechanism of PA exposure is unclear. The transcript abundance and differentially expressed genes (DEGs) in gills of bay scallop were investigated following 80 mg/L PA exposure up to 48 h using the Illumina HiSeq 4000 deep-sequencing platform with the recommended read length of 100 bp. De novo assembly of paired-end reads yielded 62,099 unigenes; 5414 genes were identified as being significantly increased, and 4452 were decreased. Based on gene ontology classification and enrichment analysis, the ‘cellular process’, ‘metabolic process’, ‘response to stimulus’, and ‘catalytic process’ with particularly high functional enrichment were revealed. The DEGs, which are related to detoxification and immune responses, revealed that acid phosphatase, fibrinogen C domain-containing protein, cyclic AMP-responsive element-binding protein, glutathione reductase, ATP-binding cassette, nuclear factor erythroid 2-related factor, NADPH2:quinone reductase, and cytochrome P450 4F22, 4B1, and 2C8-related gene expression decreased. In contrast, some genes related to glutathione S-transferase, C-type lectin, superoxide dismutase, toll-like receptors, and cytochrome P450 2C14, 2U1, 3A24 and 4A2 increased. The results of current research will be a valuable resource for the investigation of gene expression stimulated by PA, and will help understanding of the molecular mechanisms underlying the scallops’ response to PA exposure.

Accumulation of Dinophysis Toxins in Bivalve Molluscs

Several species of the dinoflagellate genus Dinophysis produce toxins that accumulate in bivalves when they feed on populations of these organisms. The accumulated toxins can lead to intoxication in consumers of the affected bivalves. The risk of intoxication depends on the amount and toxic power of accumulated toxins. In this review, current knowledge on the main processes involved in toxin accumulation were compiled, including the mechanisms and regulation of toxin acquisition, digestion, biotransformation, compartmentalization, and toxin depuration. Finally, accumulation kinetics, some models to describe it, and some implications were also considered.

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.

Multixenobiotic resistance in Mytilus edulis: Molecular and functional characterization of an ABCG2- type transporter in hemocytes and gills

Among the cellular protection arsenal, ABC transporters play an important role in xenobiotic efflux in marine organisms. Two pumps belonging to B and C subfamily has been identified in Mytilus edulis. In this study, we investigated the presence of the third major subtype ABCG2/BCRP protein in mussel tissues. Transcript was expressed in hemocytes and with higher level in gills. Molecular characterization revealed that mussel ABCG2 transporter shares the sequence and organizational structure with mammalian and molluscan orthologs. Overall identity of the predicted amino acid sequence with corresponding homologs from other organisms was between 49% and 98%. Moreover, protein efflux activity was demonstrated using a combination of fluorescent allocrites and specific inhibitors. The accumulation of bodipy prazosin and pheophorbide A was heterogeneous in gills and hemocytes. Most of the used blockers enhanced probe accumulation at different levels, most significantly for bodipy prazosin. Moreover, Mrp classical blocker MK571 showed a polyspecificity. In conclusion, our data demonstrate that several ABC transporters contribute to MXR phenotype in the blue mussel including ABCG2 that forms an active pump in hemocytes and gills. Efforts are needed to distinguish between the different members and to explore their single function and specificity towards allocrites and chemosensitizers.

Transcriptional response after exposure to domoic acid-producing Pseudo-nitzschia in the digestive gland of the mussel Mytilus galloprovincialis

Bivalve molluscs are filter feeding species that can accumulate biotoxins in their body tissues during harmful algal blooms. Amnesic Shellfish Poisoning (ASP) is caused by species of the diatom genus Pseudo-nitzschia, which produces the toxin domoic acid. The Mytilus galloprovincialis digestive gland transcriptome was de novo assembled based on the sequencing of 12 cDNA libraries, six obtained from control mussels and six from mussels naturally exposed to domoic acid-producing diatom Pseudo-nitzschia australis. After de novo assembly 94,727 transcripts were obtained, with an average length of 1015 bp and a N50 length of 761 bp. The assembled transcripts were clustered (homology > 90%) into 69,294 unigenes. Differential gene expression analysis was performed (DESeq2 algorithm) in the digestive gland following exposure to the toxic algae. A total of 1158 differentially expressed unigenes (absolute fold change > 1.5 and p-value < 0.05) were detected: 686 up-regulated and 472 down-regulated. Several membrane transporters belonging to the family of the SLC (solute carriers) were over-expressed in exposed mussels. Functional enrichment was performed using Pfam annotations obtained from the genes differentially expressed, 37 Pfam families were found to be significantly (FDR adjusted p-value < 0.1) enriched. Some of these families (sulfotransferases, aldo/keto reductases, carboxylesterases, C1q domain and fibrinogen C-terminal globular domain) could be putatively involved in detoxification processes, in the response against of the oxidative stress and in immunological processes. Protein network analysis with STRING algorithm found alteration of the Notch signaling pathway under the action of domoic acid-producing Pseudo-nitzschia. In conclusion, this study provides a high quality reference transcriptome of M. galloprovincialis digestive gland and identifies potential genes involved in the response to domoic acid.

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.

P-glycoprotein and its inducible expression in three bivalve species after exposure to Prorocentrum lima

P-glycoprotein (P-gp or ABCB1) belongs to the family of ATP-binding cassette (ABC) transporters responsible for multixenobiotic resistance (MXR) in aquatic organisms. To provide more information of P-gp in shellfish, in this study, complete cDNA of P-gp in three bivalve species including Ruditapes philippinarum, Scapharca subcrenata and Tegillarca granosa were cloned and its expressions in gill, digestive gland, adductor muscle and mantle of the three bivalves were detected after exposure to Prorocentrum lima, a toxogenic dinoflagellate. The complete sequences of R. philippinarum, S. subcrenata and T. granosa P-gp showed high homology with MDR/P-gp/ABCB proteins from other species, having a typical sequence organization as full transporters from the ABCB family. Phylogenetic analyses revealed that the amino acid sequences of P-gp from S. subcrenata and T. granosa had a closest relationship, forming an independent branch, then grouping into the other branch with Mytilus californianus, Mytilus galloprovincialis and Crassostrea gigas. However, P-gp sequences from R. philippinarum were more similar to the homologs from the more distantly related Aplysia californica than to homologs from S. subcrenata and T. granosa, suggesting that bivalves P-gp might have different paralogs. P-glycoprotein expressed in all detected tissues but there were large differences between them. After exposure to P. lima, the expression of P-gp changed in the four tissues in varying degrees within the same species and between different species, but the changes in mRNA and protein level were not always synchronous.