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

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.

Small RNA analysis of Perna viridis after exposure to Prorocentrum lima, a DSP toxins-producing dinoflagellate

Diarrheic shellfish poisoning toxins (DSP toxins) are a set of the most important phycotoxins produced by some dinoflagellates. Studies have shown that DSP toxins have various toxicities such as genotoxicity, cytotoxicity, and immunotoxicity to bivalve mollusks. However, these toxicities appear decreasing with exposure time and concentration of DSP toxins. The underlying mechanism involved remains unclear. In this study, small RNA sequencing was performed in the digestive gland of the mussel Perna viridis after exposure to DSP toxins-producing dinoflagellate Prorocentrum lima for different time periods. The potential roles of miRNAs in response and detoxification to DSP toxins in the mussel were analyzed. Small RNA sequencing of 12 samples from 72 individuals was conducted by BGISEQ-500. A total of 123 mature miRNAs were identified, including 90 conserved miRNAs and 33 potential novel miRNAs. After exposure to P. lima, multiple important miRNAs displayed some alterations. Further miRNA target prediction revealed some important genes involved in cytoskeleton, apoptosis, complement system and immune stress. qPCR demonstrated that miR-71_5, miR-750_1 and novel_mir4 were significantly up-regulated at 6 h after exposure to P. lima, while miR-100_2 was significantly down-regulated after 96 h of exposure. Accordingly, putative target genes of these differentially expressed miRNAs experienced some changes. After 6 h of DSP toxins exposure, NHLRC2 and C1q-like were significantly down-regulated. After 96 h of DSP toxins exposure, NHLRC2 was significantly up-regulated. It is reasonable to speculate that the mussel P. viridis might respond to DSP toxins through miR-750_1, novel_mir4 and miR-71_5 regulating the expression of relevant target genes involved in apoptosis, cytoskeleton, and immune response, etc. This study might provide new clues to uncover the toxic response of bivalve to DSP toxins and lay a foundation for revealing the roles of miRNAs in the environmental adaptation in shellfish.

Cloning and expression of the ChGstα and ChGstκ genes in the gills of Crassostrea hongkongensis under nanoparticulate and ionic Zn stress

Nanoparticulate and ionic Zn have potential impacts on the detoxification systems of organisms, and Gst genes play key roles in the detoxification of xenobiotics. In this study, we cloned the ChGstα and ChGstκ genes of C. hongkongensis, and studied their expression in gills under nanoparticulate and ionic Zn stress. The results showed that the coding sequences of the ChGstα and ChGstκ genes were 684 and 675 bp, respectively, and had no signal peptide; ChGstα was cytoplasmic, while ChGstκ was mitochondrial. The two genes were expressed in all 8 tested samples, with the most abundant expression observed in hemocytes for ChGstα and digestive glands for ChGstκ. After ZnCl₂ or ZnoNP challenge, the expression of ChGstα decreased significantly in the ZnCl₂ groups, and its expression was higher in the ZnoNP groups than in the ZnCl₂ groups. The expression of ChGstκ was significantly decreased in the ZnCl₂ and ZnoNP groups, and its expression was higher in the ZnoNP groups than in the ZnCl₂ groups except at 3 h post metal Zn stress, which suggested that ChGstα and ChGstκ were more sensitive to ZnoNP than ZnCl₂.

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.

Toxin accumulation, detoxification and oxidative stress in bivalve (Anomalocardia flexuosa) exposed to the dinoflagellate Prorocentrum lima

Prorocentrum lima is a cosmopolitan benthic dinoflagellate capable of producing the diarrhetic shellfish toxins (DSTs) okadaic acid (OA) and dinophysistoxin (DTX). These compounds may cause oxidative stress and accumulate in bivalve tissues, which become vectors of intoxication to human consumers. We investigated DST accumulation, detoxification and oxidative stress biomarkers in clams (Anomalocardia flexuosa) experimentally exposed to P. lima cells or their compounds. Experimental diets consisted of 6000 cells mL-1 of the non-toxic chlorophyte Tetraselmis sp. (C; control condition), and combinations of C with 10 P. lima cells mL-1 (T10), 100 P. lima cells mL-1 (T100), or to a toxin concentration of ∼4 μg OA L-1 and ∼0.65 μg DTX-1 L-1 (T100d). Clams were exposed to these diets for 7 days (uptake phase), followed by a 7-day depuration period. No DSTs were detected in clams exposed to treatments C (control) nor to T100d (dissolved compounds) during either uptake or detoxification phase. Conversely, clams exposed to T10 or T100 accumulated, on average, up to 2.5 and 35 μg DST kg-1 in their whole bodies at the end of the uptake phase. These concentrations are ∼64 and ∼4.5 times lower than the regulatory level of 160 μg OA kg-1, respectively. Accumulated OA quotas were 12-22 times higher in the digestive gland (DG) than in remaining tissues over the uptake phase. Quick toxin transformation was indicated by the early detection of conjugated compounds - DTX-1 and OA esters - in the DG after 6 h of exposure, with OA-ester representing the main compound (30 - 100 %) in that tissue over the experiment. During the depuration period, detoxification rates represented 0.024 h-1, 0.04 h-1 and 0.052 h-1 for OA, DTX-1 and OA-ester, respectively. The activities of catalase, glutathione S-transferase, glutathione peroxidase and the levels of oxidative stress by lipoperoxidation varied similarly in the DG of A. flexuosa individuals subjected to T100, T100d and the control condition. However, contrasting antioxidant responses were measured in those exposed to T10. These findings indicate that no oxidative stress was primarily induced by DST-producing dinoflagellates in this clam species under laboratory conditions representative of toxic bloom situations. Even though, possible interactions should be considered under multistressor scenarios.

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.

Differential responses in the biotransformation systems of the oyster Crassostrea gigas (Thunberg, 1789) elicited by pyrene and fluorene: Molecular, biochemical and histological approach - Part II

Pyrene (PYR) and fluorene (FLU) are among the sixteen priority Polycyclic Aromatic Hydrocarbons (PAH) of the United States Environmental Protection Agency and are both frequently detected in contaminated sites. Due to the importance of bivalve mollusks in biomonitoring programs and the scarce information on the biotransformation system in these organisms, the aim of this study was to investigate the effect of PYR and FLU at the transcriptional level and the enzymatic activities of some biotransformation systems in the Pacific oyster Crassostrea gigas, and to evaluate the histological effects in their soft tissues. Oysters C. gigas were exposed for 24 h and 96 h to PYR (0.25 and 0.5 μM) and FLU (0.6 and 1.2 μM). After exposure, transcript levels of cytochrome P450 coding genes (CYP1-like, CYP2-like, CYP2AU2, CYP356A1, CYP17α-like), glutathione S tranferase genes (omega GSTO-like and microsomal, MGST-like) and sulfotransferase gene (SULT-like), and the activity of ethoxyresorufin O-deethylase (EROD), Glutathione S-transferase (GST) and microssomal GST (MGST) were evaluated in gills. Histologic changes were also evaluated after the exposure period. PYR and FLU bioconcentrated in oyster soft tissues. The half-life time of PYR in water was lower than fluorene, which is in accordance to the higher lipophilicity and bioconcentration of the former. EROD activity was below the limit of detection in all oysters exposed for 96 h to PYR and FLU. The reproductive stage of the oysters exposed to PYR was post-spawn. Exposure to PYR caused tubular atrophy in digestive diverticula, but had no effect on transcript levels of biotransformation genes. However, the organisms exposed for 96 h to PYR 0.5 μM showed higher MGST activity, suggesting a protective role against oxidative stress in gills of oysters under higher levels of PYR in the tissues. Increased number of mucous cells in mantle were observed in oysters exposed to the higher FLU concentration, suggesting a defense mechanisms. Oysters exposed for 24 h to FLU 1.2 μM were in the ripe stage of gonadal development and showed higher transcript levels of CYP2AU2, GSTO-like and SULT-like genes, suggesting a role in the FLU biotransformation. In addition, after 96 h of exposure to FLU there was a significant increase of mucous cells in the mantle of oysters but no effect was observed on the EROD, total GST and MGST activities. These results suggest that PAH have different effects on transcript levels of biotransformation genes and enzyme activities, however these differences could also be related to the reproductive stage.

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.

Differential responses in the biotransformation systems of the oyster Crassostrea gasar (Adanson, 1757) elicited by pyrene and fluorene: molecular, biochemical and histological approach - Part I

Polycyclic aromatic hydrocarbons (PAHs) are among the main contaminants in aquatic environments. PAHs can affect organisms due to their carcinogenic, mutagenic and/or teratogenic characteristics. Depending on the PAHs, concentration, and period of exposure, biological damage can occur leading to histopathologic alterations. This study aimed to evaluate the molecular, biochemical and histological responses of the oyster Crassostrea gasar exposed to pyrene (0.25 and 0.5 μM) and fluorene (0.6 and 1.2 μM), after exposure for 24 and 96 h. Concentrations of both PAHs were quantified in the water and in oyster tissues. Transcript levels of phase I (CYP3475C1, CYP2-like, CYP2AU1 and CYP356A) and phase II (GSTO-like, MGST-like and SULT-like) biotransformation-related genes and the activities of ethoxyresorufin-O-deethylase (EROD), total and microsomal glutathione S-transferase (GST and MGST) were evaluated in the gills. Also, histological changes and localization of mRNA transcripts CYP2AU1 in gills, mantle, and digestive diverticula were evaluated. Both PAHs accumulated in oyster tissues. Pyrene half-life in water was significantly lower than fluorene. Transcript levels of all genes were higher in oysters exposed to of pyrene 0.5 μM (24 h). Only CYP2AU1 gene was up-regulated by fluorene exposure. EROD and MGST activities were higher in oysters exposed to pyrene. Tubular atrophy in the digestive diverticula and an increased number of mucous cells in the mantle were observed in oysters exposed to pyrene. CYP2AU1 transcripts were observed in different tissues of pyrene-exposed oysters. A significant correlation was observed between tubular atrophy and the CYP2AU1 hybridization signal in oysters exposed to pyrene, suggesting the sensibility of the species to this PAH. These results suggest an important role of biotransformation-related genes and enzymes and tissue alterations associated to pyrene metabolism but not fluorene. In addition, it reinforces the role of CYP2AU1 gene in the biotransformation process of PAHs in the gills of C. gasar.

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.

The oyster immunity

Oysters, the common name for a number of different bivalve molluscs, are the worldwide aquaculture species and also play vital roles in the function of ecosystem. As invertebrate, oysters have evolved an integrated, highly complex innate immune system to recognize and eliminate various invaders via an array of orchestrated immune reactions, such as immune recognition, signal transduction, synthesis of antimicrobial peptides, as well as encapsulation and phagocytosis of the circulating haemocytes. The hematopoietic tissue, hematopoiesis, and the circulating haemocytes have been preliminary characterized, and the detailed annotation of the Pacific oyster Crassostrea gigas genome has revealed massive expansion and functional divergence of innate immune genes in this animal. Moreover, immune priming and maternal immune transfer are reported in oysters, suggesting the adaptability of invertebrate immunity. Apoptosis and autophagy are proved to be important immune mechanisms in oysters. This review will summarize the research progresses of immune system and the immunomodulation mechanisms of the primitive catecholaminergic, cholinergic, neuropeptides, GABAergic and nitric oxidase system, which possibly make oysters ideal model for studying the origin and evolution of immune system and the neuroendocrine-immune regulatory network in lower invertebrates.

Upregulating Nrf2-dependent antioxidant defenses in Pacific oysters Crassostrea gigas: Investigating the Nrf2/Keap1 pathway in bivalves

Analysis of the Pacific oyster Crassostrea gigas annotated genome revealed genes with conserved sequences belonging to typical cap 'n' collar Nrf2 domain, a major player in antioxidant protection, and domains belonging to Nrf2 cytoplasmic repressor (Keap1), but little is known about Nrf2/Keap1 induction in bivalves. C. gigas were exposed to waterborne 10 and 30μM curcumin, a known inducer of the mammalian Nrf2. Curcumin disappeared from the seawater after 10h, and accumulated in the gills (10h) and digestive gland (10-96h). A clear induction of glutathione (GSH)-related antioxidant defenses was observed at 96h in the gills of curcumin exposed animals (10 and 30μM), including GSH levels, and the activity of glutathione reductase (GR), glutathione peroxidase (GPx), and glutathione S-transferase (GST). This response was completely absent in the digestive gland, in line with the idea that bivalve gills act as a major site for antioxidant protection under acute exposure. The relative mRNA levels coding glutamate-cysteine ligase, GR, GPx2 and GSTpi were clearly induced by curcumin treatment (30μM, 24h). Curcumin pre-treatment for 96h increased oyster resistance to cumene hydroperoxide, but neither Nrf2 nor Keap1 genes were modulated by curcumin. However, the conserved sequences belonging to typical Nrf2 and Keap1 domains, and the notorious induction of antioxidant defense-related genes known to be controlled by Nrf2 in mammals, indicates a functional Nrf2/Keap1 pathway in bivalves, and curcumin seems to be a new tool to investigate the antioxidant response in bivalves.