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.

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.

Conservation of Gbx genes from EHG homeobox in bivalve molluscs.

Homeobox-containing genes encode a set of transcription factors that have been shown to control spatial patterning mechanisms in bilaterian organism development. The homeobox gene Gbx, included in the EHGbox cluster, is implicated in the development of the nervous system. In this study, we surveyed five different families of Bivalvia for the presence of Gbx genes by means of PCR with degenerate primers. We were able to recover seven Gbx gene fragments from five bivalve species: Solen marginatus, Mimachlamys varia, Venerupis pullastra, Ostrea edulis and Mytilus galloprovincialis (the derived amino acid sequence were designated Sma-Gbx, Cva-Gbx, Vpu-Gbx, Oed-Gbx and Mga-Gbx, respectively). These genes are orthologous to various Gbx genes present in bilaterian genomes. The Gbx genes in four Bivalvia families, namely Solenidae, Veneridae, Ostreidae and Mytilidae, are newly reported here and we also showed additional information of the Gbx genes of Pectinidae. The phylogenetic analyses by neighbour-joining, UPGMA, maximum parsimony and Bayesian analysis clearly indicated that the Gbx sequences formed a well supported clade and assigned these Gbx genes to the Gbx family. These data permit to confirm that the homeodomain of the Gbx family is highly conserved among these five distinct families of bivalve molluscs.