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

Integrative Inducer Intervention and Transcriptomic Analyses Reveal the Metabolism of Paralytic Shellfish Toxins in Azumapecten farreri

Paralytic shellfish toxins (PSTs) are widely distributed neurotoxins, and the PST metabolic detoxification mechanism in bivalves has received increasing attention. To reveal the effect of phase I (cytochrome P450)-II (GST)-III (ABC transport) metabolic systems on the PST metabolism in Azumapecten farreri, this study amplified stress on the target systems using rifampicin, dl-α-tocopherol, and colchicine; measured PST levels; and conducted transcriptomic analyses. The highest toxin content reached 1623.48 μg STX eq/kg in the hepatopancreas and only 8.8% of that in the gills. Inducer intervention significantly decreased hepatopancreatic PST accumulation. The proportional reductions in the rifampicin-, dl-α-tocopherol-, and colchicine-induced groups were 55.3%, 50.4%, and 36.1%, respectively. Transcriptome analysis showed that 11 modules were significantly correlated with PST metabolism (six positive/five negative), with phase I CYP450 and phase II glutathione metabolism significantly enriched in negatively correlated pathways. Twenty-three phase I-II-III core genes were further validated using qRT-PCR and correlated with PST metabolism, revealing that CYP46A1, CYP4F6, GSTM1, and ABCF2 were significantly correlated, while CYP4F11 and ABCB1 were indirectly correlated. In conclusion, phase I-II-III detoxification enzyme systems jointly participate in the metabolic detoxification of PSTs in A. farreri. This study provides key data support to profoundly elucidate the PST metabolic detoxification mechanism in bivalves.

A review on aquatic toxins - Do we really know it all regarding the environmental risk posed by phytoplankton neurotoxins?

Aquatic toxins are potent natural toxins produced by certain cyanobacteria and marine algae species during harmful cyanobacterial and algal blooms (CyanoHABs and HABs, respectively). These harmful bloom events and the toxins produced during these events are a human and environmental health concern worldwide, with occurrence, frequency and severity of CyanoHABs and HABs being predicted to keep increasing due to ongoing climate change scenarios. These contexts, as well as human health consequences of some toxins produced during bloom events have been thoroughly reviewed before. Conversely, the wider picture that includes the non-human biota in the assessment of noxious effects of toxins is much less covered in the literature and barely covered by review works. Despite direct human exposure to aquatic toxins and related deleterious effects being responsible for the majority of the public attention to the blooms' problematic, it constitutes a very limited fraction of the real environmental risk posed by these toxins. The disruption of ecological and trophic interactions caused by these toxins in the aquatic biota building on deleterious effects they may induce in different species is paramount as a modulator of the overall magnitude of the environmental risk potentially involved, thus necessarily constraining the quality and efficiency of the management strategies that should be placed. In this way, this review aims at updating and consolidating current knowledge regarding the adverse effects of aquatic toxins, attempting to going beyond their main toxicity pathways in human and related models' health, i.e., also focusing on ecologically relevant model organisms. For conciseness and considering the severity in terms of documented human health risks as a reference, we restricted the detailed revision work to neurotoxic cyanotoxins and marine toxins. This comprehensive revision of the systemic effects of aquatic neurotoxins provides a broad overview of the exposure and the hazard that these compounds pose to human and environmental health. Regulatory approaches they are given worldwide, as well as (eco)toxicity data available were hence thoroughly reviewed. Critical research gaps were identified particularly regarding (i) the toxic effects other than those typical of the recognized disease/disorder each toxin causes following acute exposure in humans and also in other biota; and (ii) alternative detection tools capable of being early-warning signals for aquatic toxins occurrence and therefore provide better human and environmental safety insurance. Future directions on aquatic toxins research are discussed in face of the existent knowledge, with particular emphasis on the much-needed development and implementation of effective alternative (eco)toxicological biomarkers for these toxins. The wide-spanning approach followed herein will hopefully stimulate future research more broadly addressing the environmental hazardous potential of aquatic toxins.

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.

Self-assembled monolayer-based immunoassays for okadaic acid detection in seawater as monitoring tools

Rapid and cost-effective methods to monitor the presence of diarrhetic shellfish poisoning (DSP) toxins in seawater samples in an easy and reliable manner are required to protect human health and avoid economic losses to shellfish industry. Immunoassays for the detection of okadaic acid (OA) and dinophysistoxin-1 and dinophysistoxin-2 are developed by immobilising OA on self-assembled monothiols or dithiols in an ordered and oriented way, providing an effective limit of detection of ∼1 ng OA equiv./mL seawater. The immunoassays are applied to the analysis of the particulate fraction of seawater samples from two Catalan harbours (NW Mediterranean) and samples collected periodically from the Galician Rias (E Atlantic), as well as a reference mussel sample. Results are in agreement with LC-MS/MS and the certified values. OA concentration in seawater correlates with Dinophysis cell abundance, with a 1-2 weeks lag. The immunoassays provide powerful high-throughput analytical methods potentially applicable as alternative monitoring tools.

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.

Effects of Florida Red Tides on histone variant expression and DNA methylation in the Eastern oyster Crassostrea virginica

Massive algal proliferations known as Harmful Algal Blooms (HABs) represent one of the most important threats to coastal areas. Among them, the so-called Florida Red Tides (FRTs, caused by blooms of the dinoflagellate Karenia brevis and associated brevetoxins) are particularly detrimental in the southeastern U.S., causing high mortality rates and annual losses in excess of $40 million. The ability of marine organisms to cope with environmental stressors (including those produced during HABs) is influenced by genetic and epigenetic mechanisms, the latter resulting in phenotypic changes caused by heritable modifications in gene expression, without involving changes in the genetic (DNA) sequence. Yet, studies examining cause-effect relationships between environmental stressors, specific epigenetic mechanisms and subsequent responses are still lacking. The present work contributes to increase this knowledge by investigating the effects of Florida Red Tides on two types of mechanisms participating in the epigenetic memory of Eastern oysters: histone variants and DNA methylation. For that purpose, a HAB simulation was conducted in laboratory conditions, exposing oysters to increasing concentrations of K. brevis. The obtained results revealed, for the first time, the existence of H2A.X, H2A.Z and macroH2A genes in this organism, encoding histone variants potentially involved in the maintenance of genome integrity during responses to the genotoxic effect of brevetoxins. Additionally, an increase in H2A.X phosphorylation (γH2A.X, a marker of DNA damage) and a decrease in global DNA methylation were observed as the HAB simulation progressed. Overall, the present work provides a basis to better understand how epigenetic mechanisms participate in responses to environmental stress in marine invertebrates, opening new avenues to incorporate environmental epigenetics approaches into management and conservation programs.

Characterization of mussel H2A.Z.2: a new H2A.Z variant preferentially expressed in germinal tissues from Mytilus

Histones are the fundamental constituents of the eukaryotic chromatin, facilitating the physical organization of DNA in chromosomes and participating in the regulation of its metabolism. The H2A family displays the largest number of variants among core histones, including the renowned H2A.X, macroH2A, H2A.B (Bbd), and H2A.Z. This latter variant is especially interesting because of its regulatory role and its differentiation into 2 functionally divergent variants (H2A.Z.1 and H2A.Z.2), further specializing the structure and function of vertebrate chromatin. In the present work we describe, for the first time, the presence of a second H2A.Z variant (H2A.Z.2) in the genome of a non-vertebrate animal, the mussel Mytilus. The molecular and evolutionary characterization of mussel H2A.Z.1 and H2A.Z.2 histones is consistent with their functional specialization, supported on sequence divergence at promoter and coding regions as well as on varying gene expression patterns. More precisely, the expression of H2A.Z.2 transcripts in gonadal tissue and its potential upregulation in response to genotoxic stress might be mirroring the specialization of this variant in DNA repair. Overall, the findings presented in this work complement recent reports describing the widespread presence of other histone variants across eukaryotes, supporting an ancestral origin and conserved role for histone variants in chromatin.

The characterization of macroH2A beyond vertebrates supports an ancestral origin and conserved role for histone variants in chromatin

Histone variants play a critical role in chromatin structure and epigenetic regulation. These "deviant" proteins have been historically considered as the evolutionary descendants of ancestral canonical histones, helping specialize the nucleosome structure during eukaryotic evolution. Such view is now challenged by 2 major observations: first, canonical histones present extremely unique features not shared with any other genes; second, histone variants are widespread across many eukaryotic groups. The present work further supports the ancestral nature of histone variants by providing the first in vivo characterization of a functional macroH2A histone (a variant long defined as a specific refinement of vertebrate chromatin) in a non-vertebrate organism (the mussel Mytilus) revealing its recruitment into heterochromatic fractions of actively proliferating tissues. Combined with in silico analyses of genomic data, these results provide evidence for the widespread presence of macroH2A in metazoan animals, as well as in the holozoan Capsaspora, supporting an evolutionary origin for this histone variant lineage before the radiation of Filozoans (including Filasterea, Choanoflagellata and Metazoa). Overall, the results presented in this work help configure a new evolutionary scenario in which histone variants, rather than modern "deviants" of canonical histones, would constitute ancient components of eukaryotic chromatin.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

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.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.

Identification and Characterization of a Novel Family of Cysteine-Rich Peptides (MgCRP-I) from Mytilus galloprovincialis

We report the identification of a novel gene family (named MgCRP-I) encoding short secreted cysteine-rich peptides in the Mediterranean mussel Mytilus galloprovincialis. These peptides display a highly conserved pre-pro region and a hypervariable mature peptide comprising six invariant cysteine residues arranged in three intramolecular disulfide bridges. Although their cysteine pattern is similar to cysteines-rich neurotoxic peptides of distantly related protostomes such as cone snails and arachnids, the different organization of the disulfide bridges observed in synthetic peptides and phylogenetic analyses revealed MgCRP-I as a novel protein family. Genome- and transcriptome-wide searches for orthologous sequences in other bivalve species indicated the unique presence of this gene family in Mytilus spp. Like many antimicrobial peptides and neurotoxins, MgCRP-I peptides are produced as pre-propeptides, usually have a net positive charge and likely derive from similar evolutionary mechanisms, that is, gene duplication and positive selection within the mature peptide region; however, synthetic MgCRP-I peptides did not display significant toxicity in cultured mammalian cells, insecticidal, antimicrobial, or antifungal activities. The functional role of MgCRP-I peptides in mussel physiology still remains puzzling.

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.

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.